Predicting through-thickness cone cracking of reinforced concrete slabs struck normally by flat-nosed projectiles
-
摘要: 主要针对钢筋混凝土靶板在受到平头弹丸撞击下发生的厚度方向开裂的问题进行研究,并提出了一个弹丸低速撞击有限厚度板的二阶段模型。模型中第一阶段为侵彻阶段,弹丸受到混凝土介质的侵彻阻力由静阻力和速度效应引起的动阻力组成;模型中第二阶段为开裂阶段,钢筋混凝土靶板发生动态剪切破坏的最大承载力可以通过静态剪切破坏最大承载力乘以一个动态增强因子得到。该模型可以用来预测钢筋混凝土靶板发生厚度方向开裂破坏的临界能量。模型预测与实验结果吻合较好。Abstract: A two-stage model is proposed herein to predict the through-thickness cone cracking of a reinforced concrete slab struck transversely by a flat-ended projectile at relatively low velocities. The first stage is penetration, in which the resistive pressure from the concrete medium during the penetration process is composed of two parts: quasi-static resistive pressure and dynamic resistive pressure arising from velocity effect. The second stage is cone cracking, in which the dynamic resistive load is obtained through the punch shear strength for the concrete slab loaded quasi-statically by introducing a dynamic enhancement factor. A semi-analytical equation is derived for predicting the critical impact energy that causes the through-thickness cone cracking of the reinforced concrete slab subjected to impact by the flat-nosed projectile. It transpires that the present model predictions are in good agreement with available test data.
-
在爆炸,切削,高速撞击等过程中,由于变形在极短的时间内完成,大部分塑性功转化成的热能来不及散失,就会形成绝热剪切带(adiabatic shear band, ASB)。虽然对ASB内因高速形变造成绝热温升带来的连续式动态再结晶已有不少研究,但大部分都是基于变形前后的对比。在极短时间内,其变形过程尚不十分清楚,此时变形过程中的应力应变和绝热温升在剪切带中的分布扮演着重要角色。李继承等[1-2]利用ANSYS/LS-DYNA软件模拟921A钢帽形样经霍普金森杆高速冲击过程,发现绝热剪切带从剪切区两端向其中心扩展, 裂纹尖端存在一个高温高应变区, 该区域随着断裂的发展而逐渐向剪切区中心扩展。冀建平等[3]用同样方法模拟了45钢帽形样霍普金森杆高速冲击过程,并根据特征单元应力应变数据,计算特征单元温度随时间变化,从而确定45钢退火态材料所产生的绝热剪切带为形变。两者虽然目的不一样,材料不一样,但都通过ANSYS/LS-DYNA软件成功模拟了霍普金森杆高速冲击过程,不同模拟结果与实际都吻合。说明只要选择合适的本构模型,ANSYS/LS-DYNA软件模拟各种材料高速冲击过程是适合的。A.Marchand等[4]通过实验发现等效应变超过40%,流动应力开始迅速下降,这种应力坍塌现象随后也被T.W.Wright等[5]通过数值模拟发现,后来被众多学者作为判定ASB是否生成的标志。M.Zhou等[6]还得到了应力坍塌的产生所需临界应变的经验公式,并指出等效应变到达临界应变时发生应力坍塌。D.J.Rattazzi等[7]分别采用等效塑性应变等于0.5,等效应力达到最大值和等效应力下降到最大值的90%作为形成剪切带的临界条件,发现剪切带开始时间与选择的临界条件有关。高锰TRIP钢在形变过程中会发生相变,需要考虑模拟过程中不涉及相变是否对模拟结果影响较大。A.Khosravifard等[8]利用ANSYS软件包对高锰TRIP/TWIP钢高应变速率扭转进行有限元模拟,模拟过程中不考虑组织变化,仅输入力学和热学参数,并采用简易的时间力矩曲线进行加载,主要观察剪切应力应变分布及应变率软化行为和高速形变过程中绝热温升,发现模拟结果与实际吻合。说明不考虑组织变化,依然可以模拟剪切带内外应力应变场、温度场分布。
本文中选择Johnson-Cook本构模型和Grüneisen状态方程,利用ANSYS/LS-DYNA软件模拟高锰TRIP钢高速冲击过程,可对实验结果进行补充,有利于更加全面了解高锰TRIP钢高速冲击过程。
1. 热粘塑性本构模型
Johnson-Cook模型形式简单、参数较少,是目前应用最广泛且与应变率相关的材料本构模型之一,并同时考虑到了金属材料的应变强化效应、应变率强化效应和热软化效应。Johnson-Cook本构模型:
σ=(A+Bεnp)(1+Cln˙ε˙ε0)(1−T∗m) (1) 式中:εp为等效塑性应变,A为准静态下的屈服应力,B为应变硬化系数,n为应变硬化指数,˙ε0为参考应变率(可取准静态应变率),C为应变率敏感系数,T*=(T-Tr)/(Tm-Tr), Tr和Tm分别为参考温度和熔化温度,一般取Tr为300 K,m为热软化系数。
由于实验在室温进行,即T*=0,故可将式(1)转化为:
σ=(A+Bεnp)(1+Cln˙ε˙ε0) (2) 当实验为准静态状态时,可将式(2)变为转化为:
σ=A+Bεnp (3) 将式(3)两端取对数,可转化为
ln[(σ−A)/MPa]=nlnε+ln(B/MPa) (4) 最后利用准静态实验数据进行曲线拟合,即可得到A、B和n。
实验采用高锰钢各组分的质量分数为:Mn,17.22%;Si,2.87%;Al,0.48%;C,0.022%。将高锰钢准静态下的实验数据进行曲线拟合,最终得到材料屈服强度A=236.58 MPa,曲线斜率n=0.76,截距ln(B/MPa)=7.525 8, 即B=1 948 MPa。
当ε=0时,室温下Johnson-Cook本构模型可以转化为屈服强度与应变率的函数:
σs=A(1+Cln˙ε˙ε0) (5) 根据实验得到的不同应变率下的屈服强度值,利用最小二乘法拟合得出C=0.052。同理m可由准静态下不同温度下的应力应变数据拟合得到,取m=1[9]。状态方程中的参量和损伤失效参量都需要大量的实验才能确定,受实验条件限制,取值参考文献[10]。
2. 有限元模型及验证
模拟的帽形样品尺寸如图 1所示。鉴于系统的轴对称特性,数值模型简化为二维轴对称模型,如图 2所示。
为了更精确地反映数据变化规律,在变化梯度较大的强制剪切区应细化网格。利用实验得到的应力时程曲线,如图 3所示。在试样的上端施加载荷,在对实验进行初步模拟时,将载荷曲线简化为折线,将屈服之前的应力简化为一条直线,起点是原点,屈服之后的应力简化为一条直线,终点和实验得到的曲线终点重合,加载时间和实验一致,为80 μs。虚线是实验得到的应力时程曲线,实线是简化的加载曲线。试样下端固定。
通过上面的数值模拟,可以得到加载80 μs后试样的位移量为0.139 mm,然后与实验结果(实验测得变形量为0.124 mm)作对比,发现高速冲击后试样的长度有0.015 mm的误差,有一定的吻合性。
3. 模拟结果与讨论
3.1 剪切带形成判据
定义直角坐标系下特征单元位置,如图 4所示。将剪切带靠近帽顶的一端定义为剪切带的上端,靠近帽沿一端定义为剪切带的下端。将特征单元404和1941等效塑性应变达到临界应变时的平均压下量作为形成剪切带和剪切带贯穿整个剪切区域所需的平均压下量,等效塑性应变随压下率(压缩前后的长度差与原始长度的比值)变化曲线如图 5所示。根据绝热剪切带的定义,高速变形时,塑性变形功基本完全转化为热量,而热量来不及散失,引起热软化现象,变形过程中,加工硬化和热软化相互竞争,当热软化大于加工硬化时,便会出现绝热剪切带。所以当某一个单元在剪切带内,由于热软化大于加工硬化,在卸载前,该单元就会出现应力下降,而剪切带外的单元加工硬化大于热软化,只有在卸载时(当t=80 μs时)应力才会下降。
图 6所示为2个不同单元的剪切应力时程曲线,从图中可知,单元2167在卸载前剪切应力就下降了,而单元1680在卸载时剪切应力才下降。所以,单元2167在剪切带内,而单元1680不在剪切区内,且卸载前应力下降时的等效塑性应变为形成剪切带的临界等效塑性应变。图 7所示为单元2167的剪切应力等效塑性应变关系曲线,可以看出,剪切应力下降时对应的等效塑性应变为0.500。结合变形量下强制剪切区域的扫描电镜形貌图,如图 8所示,可以看出在形变量9.65%下未产生剪切带只是发生了相变,而在形变量12.40%、23.20%、25.80%下均能看到明显的剪切带,且贯穿整个强制剪切区,说明形成剪切带且贯穿整个剪切区的临界形变量处于9.65%和12.40%之间。从图 5中特征单元1941随形变量逐渐增加的等效塑性应变曲线图可以看出等效塑性应变0.500对应的形变量为12.10%,正好处于9.65%和12.4%之间,和实验观察到的现象吻合。因此,本文中把等效塑性应变为0.500作为判定ASB是否生成的标志。
3.2 剪切区内外应变分布
过单元1941,横穿剪切带取一系列单元,以x轴为横坐标(按图 4定义的坐标系),剪切应变εxy为纵坐标作曲线,如图 9所示。从图中可以看出横跨剪切带(沿x轴)的剪切应变分布类似于高斯分布,剪切区中心剪切应变最大,约为0.500,然后沿两边基体剪切应变逐渐下降到接近为零(10-3量级),且从剪切带中心到距离其563 μm处(假设此处为基体的应变)下降较快;继续向两边扩展时,下降较平缓。实际上,剪切区中心剪切应变远大于0.500,这是因为ANSYS模拟是具有网格依赖性的,要求在变化梯度较大的区域要细化网格,而上述模拟中剪切区网格细化程度远远达不到实际剪切区晶粒细化程度。模拟不同细化程度下剪切区中心剪切应变大小,发现当细化程度从1/27减小到1/9、1/3时剪切区中心剪切应变大小相应从0.500减小到0.222、0.195,说明细化程度直接影响实验结果。
图 9(b)所示为高速冲击时在不同时刻沿y轴方向(平行于剪切带)剪切带中心剪切应变的分布图,因破坏而删除的单元在图中不再显示。从图 9(b)中可以看出,剪切应变沿y轴呈两端高中间低的特点。当t=28.8 μs时,剪切区上端剪切应变高于下端,此时剪切带刚开始在强制剪切区上端形成。因此剪切带首先在施加载荷的一端形成,接着另一端也逐渐形成,然后剪切带由两端向中间扩展。
3.3 剪切区内外温度分布
通过编辑软件的k文件可以直接输出系统的温度数据。编辑k文件过程中假设材料变形能的90%转化为热能,材料的热容为450 J/(kg·K), 导热系数为12.979 W/(m·K)。发现温度的分布和等效塑性应变的分布非常类似,温升区域主要集中在剪切带区及周围形变影响区。还是过单元1941横穿剪切带和平行于剪切带取一系列单元,作出沿x轴和y轴的温升曲线,如图 10所示。因破坏而删除的单元在图中不再显示。从图中可以看出曲线的形状和应力应变基本上是一样的,横穿剪切带方向剪切带中心温度最高,约402 K,然后向两边逐渐下降。温升主要集中在x=3~3.7 mm较窄的区域内,基体的温度变化不大。这和文献[11]中模拟的绝热剪切带内温度分布特征基本吻合。平行于剪切带方向t=19.8 μs时剪切带两端温度开始上升,然后随着剪切带逐渐向中间扩展,温升区域也逐渐向中间扩展。整个试样,剪切带两端温度最高,约为590 K,这也是最容易发生断裂的地方。
通过预压缩变形30%获得的初始组织几乎没有残余奥氏体,但是在高速冲击之ASB内仍以细小的等轴的奥氏体为主,这是由于温升加上剪切力作用,发生了马氏体逆转变,从而形成奥氏体。由于冲击时间极短,很难原位检测剪切带内部的温升。文献[12]曾报道ε-马氏体向γ-奥氏体逆转变开始温度约400 K,且随着变形量的增加而逐渐降低。根据这个依据,在本文模拟条件下,剪切带内部的温升是可以诱导马氏体逆转变发生的。
4. 讨论
对样品进行高速冲击是一个很快的过程,通常实验只能对变形前和变形后试样形貌和组织进行检测,而很难捕捉到中间的变化过程。这个时候需要借助模拟这个辅助工具来加强对试样在变形过程变化的了解。绝热剪切带是材料在高速变形时一种典型的破坏形式,本文中选择形成剪切带的临界等效塑性应变为0.500,此值并不影响应变和温度的分布。临界塑性应变依赖于材料和加载条件。不同的材料要根据材料本身特性选择不同的临界条件;若高应变率下加载时能形成剪切带,那么即使应变达到临界塑性应变0.500,在低应变率下加载时,强制剪切区也不一定能形成剪切带。通过帽型样形成强制剪切区,着重观察剪切区内应变和温度分布。通过模拟结果,发现绝热剪切带和温升区域都是从剪切区两端开始形成,然后向中间扩展。横穿剪切带方向,应变和温度分布都是剪切带中心最高,然后向两边逐渐降低,类似于高斯分布。平行于剪切带方向,应变和温度分布都是呈两端高中间低的特点。
5. 结论
基于ANSYS/LS-DYNA软件模拟高锰TRIP钢在高速冲击过程中绝热剪切带的形成,得到以下规律:
(1) 随着压缩变形的增加,剪切带逐渐从帽形样的强制剪切区的两端开始形成,并向中间扩展。
(2) 横穿剪切带的方向,从剪切带中心到周围基体剪切应变的分布类似于高斯分布。剪切带中心处剪切应变最大,从剪切带中心到距离其563 μm处下降较快,然后继续向两边扩展时下降较平缓。剪切带中心处剪切应变约为0.500,刚到基体时下降为0.440,距剪切带中心约563 μm处靠近试样对称轴侧基体的剪切应变是0.058,到试样中心或两边边界时剪切应变下降到10-3量级。
(3) 平行剪切带方向,从上端到下端剪切应变的分布大体上都是呈两端高、中间低的特点。
(4) 温升区域主要集中在剪切带区域和其周围形变影响区。在高速冲击的过程中,剪切区域温升首先出现在两端,然后逐渐向中间扩展。整个试样,剪切带两端温度最高,约为590 K,这也是最容易发生断裂的位置。
-
表 1 破坏模式转换临界值
Table 1. Critical values for the transition of different failure modes
编号 m/kg H/d f′c/MPa V0/(m·s-1) Hc/d 破坏模式 模型预测 实验观察[9] 55 401 1.0 25.1 6.32 3.00 CC CC 71 404 1.8 33.2 6.80 3.30 CC CC 66 396 2.7 36.3 5.07 3.53 CC CC 61 407 4.0 32.7 6.83 3.90 P+CC P+CC 77 399 6.0 32.0 14.35 4.20 P+CC P+CC 注:P为侵彻;CC为厚度方向开裂破坏。 -
[1] Kennedy R P. A review of procedures for the analysis and design of concrete structures to resist missile impact effects[J]. Nuclear Engineering and Design, 1976, 37(2):183-203. doi: 10.1016-0029-5493(76)90015-7/ [2] Li Q M, Reid S R, Wen H M, et al. Local impact effects of hard missiles on concrete targets[J]. International Journal of Impact Engineering, 2005, 32(1/2/3/4):224-284. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=f836952d7e21bc3f5f0eb712e6968392 [3] Reid S R, Wen H M. Predicting penetration, cone cracking, scabbing and perforation of reinforced concrete targets struck by flatted faced projectiles[R]. UMIST Report ME/AM/02.01/TE/G/018507/Z, 2001. [4] Wen H M, Xian Y X. A unified approach for concrete impact[J]. International Journal of Impact Engineering, 2015, 77:84-96. doi: 10.1016/j.ijimpeng.2014.11.015 [5] 咸玉席, 文鹤鸣.平头弹侵彻半无限混凝土靶的工程模型[J].防护工程, 2012, 34(2):35-38. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK201203952455Xian YuXi, Wen Heming. An engineering model for the penetration of flat-nosed projectiles into semi-infinite concrete targets[J]. Protective Engineering, 2012, 34(2):35-38. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK201203952455 [6] Wen H M, Yang Y. A note on the deep penetration of projectiles into concrete[J]. International Journal of Impact Engineering, 2014, 66:1-4. doi: 10.1016/j.ijimpeng.2013.11.008 [7] Yankelevsky D Z, Leibowitz O. Punching shear in concrete slabs[J]. International Journal of Mechanical Sciences, 1999, 41:1-15. doi: 10.1016/S0020-7403(97)00086-6 [8] 咸玉席, 文鹤鸣.钢筋混凝土平板冲剪强度的无量纲公式[J].高压物理学报, 2016, 30(4):291-300. http://www.cqvip.com/QK/96553X/201604/669715156.htmlXian Yuxi, Wen Heming. Dimensionless formulae for punching shear strength of reinforced concrete slabs[J]. Chinese Journal of High Pressure Physics, 2016, 30(4):291-300. http://www.cqvip.com/QK/96553X/201604/669715156.html [9] Sinclair A C E. Drop tests on a reinforced concrete floor at Rogerstone Power Station[R]. TD/SEB/REP/41555/93, 1993. 期刊类型引用(12)
1. 常白雪,张元瑞,王少华,彭克锋,虞吉林,郑志军. 梯度多胞材料的动态力学性能分析与设计研究综述. 爆炸与冲击. 2024(08): 5-33 . 本站查看
2. 袁良柱,陈美多,谢雨珊,陆建华,王鹏飞,徐松林. 细观非连续介质的应力波传播研究. 爆炸与冲击. 2024(09): 50-61 . 本站查看
3. 周辉,任辉启,吴祥云,易治,黄魁,穆朝民,王海露. 成层式防护结构中分散层研究综述. 爆炸与冲击. 2022(11): 3-28 . 本站查看
4. 余同希,朱凌,许骏. 结构冲击动力学进展(2010-2020). 爆炸与冲击. 2021(12): 4-64 . 本站查看
5. 刘冕,王根伟,宋辉,王彬. 负梯度泡沫金属中的局部密实化现象. 高压物理学报. 2020(04): 117-127 . 百度学术
6. 常白雪,郑志军,赵凯,何思渊. 具有恒定冲击载荷的梯度泡沫金属材料设计. 爆炸与冲击. 2019(04): 3-11 . 本站查看
7. 胡俊,任建伟,马巍,刘建华,王爱国. 冲击荷载下含随机缺陷的梯度蜂窝材料的力学性能. 材料导报. 2019(16): 2777-2784 . 百度学术
8. 胡俊,任建伟,王爱国,吴德义. 非线性梯度胞元分布蜂窝材料的冲击力学响应. 材料导报. 2019(24): 4066-4071 . 百度学术
9. 罗伟铭,石少卿,廖瑜,孙建虎. 成层式铝蜂窝夹芯结构冲击响应试验研究. 材料导报. 2018(08): 1328-1332 . 百度学术
10. 常白雪,郑志军,赵凯,虞吉林. 梯度多胞材料耐撞性设计的简化模型和渐近解. 中国科学:物理学 力学 天文学. 2018(09): 233-241 . 百度学术
11. Xu-ke Lan,Shun-shan Feng,Qi Huang,Tong Zhou. Blast response of continuous-density graded cellular material based on the 3D Voronoi model. Defence Technology. 2018(05): 433-440 . 必应学术
12. 李志斌,卢芳云. 梯度温度场中多胞材料牺牲层的抗冲击分析. 兵工学报. 2017(12): 2463-2471 . 百度学术
其他类型引用(4)
-