Damage Evolution Equation of Concrete Materials at High Temperatures and High Strain Rates
-
摘要: 采用
∅ 74 mm大口径分离式霍普金森压杆(SHPB)对不同温度(20、200、400 ℃)下的C45混凝土材料进行动态力学性能实验,得到了不同温度、不同应变率下混凝土材料的应力-应变曲线。实验结果表明:在20~400 ℃温度范围内,混凝土材料具有温度硬化和应变率硬化现象。基于上述实验数据给出了损伤变量关于塑性应变的关系式,并通过相关实验数据确定了不同温度、不同应变率下损伤演化方程的材料参数。将该损伤演化方程应用于混凝土材料的本构关系中,预测结果与实验数据具有较好的一致性,证明了所提出的高温、高应变率下混凝土材料损伤演化方程的合理性。Abstract: In this paper, the dynamic mechanical properties of C45 concrete materials at different temperatures (20, 200, 400 °C) are carried out on split Hopkinson pressure bar (SHPB) equipment with a large diameter of 74 mm. The stress-strain curves of concrete materials at different temperatures and strain rates are obtained through experiments. As the expansion of microcracks inside concrete materials is inhibited by the increase of strain rate, the concrete specimens exhibit strain rate hardening effect. The experimental results show that the concrete material has temperature hardening and strain rate hardening in the temperature range of 20 °C to 400 °C. Through the relevant theoretical derivation, the SHPB experimental data of concrete materials are transformed into the relationship between damage variables and plastic strain. Then the material parameters of the damage evolution equation at different temperatures and different strain rates are determined by relevant experimental data. Finally, the damage evolution equation of concrete materials at high temperatures and high strain rates are applied to the constitutive relation of concrete materials. The prediction results are in good agreement with the experimental data. -
目前,国内武装贩毒、劫持人质、抢劫银行等严重武装犯罪以及恐怖活动均不断出现。为应对这些危险的武装犯罪事件,需要一些停止作用大、侵彻力低的武器装备,以求能形成较高的制止犯罪能力。而我国目前主要的警用武器装备均存在枪弹的穿透力太大,导致停止作用不足的缺点。国外警方在遇到这种紧急情况时都会使用特制的低侵彻手枪弹头[1]。Kneubuehl等[2]在试验中发现,将美国M193步枪弹的撞击速度提升到750 m/s以上时,弹体在软介质中会破碎,从而达到低侵彻的效果。美国温彻斯特公司生产的“黑爪”弹,其弹体为纯黑色,射入目标后被甲会均匀地向后翻开成6瓣,就像6个带有倒钩的爪子,弹头直径最终可膨胀至口径的两倍左右。近几年国内学者也注意到了低侵彻枪弹的研究意义。李晓杰等[3]对国产5.8 mm口径手枪弹头进行改进,提出了新式的软尖弹和空尖弹,并与标准弹对比,使用ANSYS/LS-DYNA进行弹丸入水的数值仿真,得到3种弹型的速度、位移曲线,以及子弹入水产生的瞬时空腔,结果表明,空尖弹具有良好的低侵彻性能。金永喜等[4]提出了子弹在软目标介质内的翻滚模型,证明弹丸翻滚时产生的离心力能使子弹在撞击目标过程中头部发生破碎,从而使其侵彻深度得到明显的降低,并通过与制式枪弹的实验对比进行了论证。尽管近几年我国开始对低侵彻枪弹进行相关研究,但都仅局限在仿真或者理论的基础上,相关的实验验证少之又少,因此,本工作通过实验与仿真相结合,对低侵彻枪弹做进一步的研究。
根据创伤弹道学[5],子弹在撞击目标介质过程中,弹头的翻转、变形、破碎都能使弹头的能量更快地向介质传递。本工作研究一种空心开花型低侵彻弹。该弹击中目标后弹体头部会向外开裂,利用子弹头部在侵彻过程中产生的大变形,实现子弹的高能量传递性,使开花弹在低速侵彻下初速低、子弹头部变形小、侵彻阻力低,而在高速侵彻下初速高、子弹头部变形大、侵彻阻力高,最终达到开花弹在低速和高速侵彻条件下,侵彻深度始终控制在某一范围之内的目标。因此,该空心开花型低侵彻弹无论近距离还是远距离击中目标,都不会将其射穿,具备良好的低侵彻性能。在应对紧急情况时,只会制止犯罪人员而不会穿透目标误伤群众。
1. 弹丸入水实验
1.1 实验设计
鉴于生物组织的复杂性,采用水介质模拟生物组织。因为大多数生物组织含水80%左右,其密度与水相近,而且水具有均匀、透明、便于直接观察的特点。以5.8 mm口径弹道枪作为发射平台进行弹体侵彻水介质实验,实验现场布置如图 1所示,通过装药量控制弹体着靶速度,用锡箔靶和双通道测试仪测量速度。
实验中空心开花型低侵彻弹的直径为5.8 mm,长径比为5,质量为5.6 g,弹丸头部开有矩形槽,边长为2 mm,凹槽边界上沿轴向设有预开槽,如图 2所示。
1.2 实验结果及分析
通过改变子弹的装药量,得到开花弹在撞击速度为408、501、610、702、814 m/s下弹丸头部的变形情况,实验结果如表 1所示。从表 1中可以看出,开花弹入射水箱后,弹体头部发生不同程度的变形。这是由于开花弹在侵彻靶体过程中会受到较大的侵彻阻力,产生沿弹体轴向的剪切力,弹丸头部发生变形,由于弹头开有矩形凹槽,在矩形槽的轴向边界处会产生相应的拉应力,拉应力使凹槽不断向外膨胀,达到弹体材料的抗拉强度后,弹头沿矩形槽边界断裂,开裂成4瓣,断裂的部分在水的阻力作用下向后翻转,最后呈“花瓣状”。
表 1 实验结果和仿真结果对比Table 1. Comparison of experiment and simulation resultsExperiment Simulation Velocity/(m·s-1) Results Velocity/(m·s-1) Results 408 400 501 500 610 600 702 700 814 800 对比表 1可知,弹丸头部的变形程度与子弹撞击速度有关。在408 m/s较低速度下入水时,弹丸头部部分开裂,矩形槽沿径向略微膨胀。子弹以501~610 m/s速度入水时,弹丸头部开裂成4瓣并伴有外翻现象,但头部凹槽未开裂完全,仍有清晰的孔洞,此时弹丸头部最大横截面处直径为弹体直径的2倍。子弹以702~814 m/s速度入水时,弹头沿矩形槽轴向边界完全开裂成“花瓣状”,且凹槽完全消失,此时弹丸头部最大横截面处直径达到弹体直径的2.4倍。因此,开花弹在低速入射时,初速低,侵彻阻力较小,弹头变形小;在高速入射时,初速高,侵彻阻力大,弹头变形大。
2. 数值模拟的有限元建模及算法
弹头入水的模型如图 3所示。弹头在空气中,空气和水的边界采用无反射边界,空气和水的接触面采用共节点方式。
弹头由紫铜构成,保证了弹体材料的一致性。在算法的选择上:铜选用Lagrange算法,空气和水选用欧拉算法,弹头和水之间采用耦合算法[6]。在材料模型的选择上:空气和水采用LS-DYNA提供的MAT_NULL流体模型;紫铜选用Johnson-Cook材料模型,这种本构方程考虑了高速下的应变率效应和温度效应,适用于高速下的流固耦合。在状态方程的选择上:紫铜、空气和水都采用Grüneisen状态方程[5-7]。
3. 仿真结果比较
3.1 不同速度入水的头部变形比较
使用LS-DYNA软件可以较理想地模拟开花弹的入水过程,图 4为开花弹撞击水介质时弹头的变形过程。从图 4中可以看出,开花弹撞击水介质后,弹头在侵彻阻力作用下发生变形,并沿着矩形槽发生开裂,最后以花瓣状弹头继续运动。
通过模拟开花弹的入水过程,得到开花弹以不同速度撞击水介质后稳定侵彻阶段的弹体变形图,并与实验结果进行对比,如表 1所示。可以看出,仿真结果与实验中弹体变形情况基本吻合。开花弹在较低速撞击下,弹头轻微变形,矩形槽口向外扩张,但未达到开裂的程度,这是因为此时开花弹受到的侵彻阻力过小,弹头矩形槽轴向边界产生的拉应力未达到弹体材料的抗拉强度,无法使弹头开裂。开花弹以600 m/s速度撞击水靶时,弹头受到的拉应力达到一定强度,使头部沿槽口开裂,但从表 1中可以看出,此时弹丸头部未完全开裂,这是由于子弹的撞击速度仍然不足,达不到子弹需要的侵彻阻力。进一步提高开花弹的侵彻速度,在700 m/s速度撞击下,侵彻阻力较大,弹头完全开裂成4瓣,矩形槽消失。
表 2是开花弹在不同入射速度下弹丸头部最大直径的实验数据和仿真数据的对比,其中:Dmax为弹丸头部最大直径,Dmax是其平均值。可以看出,实验结果和仿真结果基本吻合。开花弹速度低于500 m/s时,弹丸头部变形较小;速度为500 m/s时,弹丸头部最大直径为弹体直径的1.2倍;子弹速度为600 m/s时,弹头变形较大,弹丸头部最大直径为弹体直径的2倍;子弹速度达到700 m/s时,开花弹弹头开裂完全,变形达到极限,弹丸头部最大直径为弹体直径的2.4倍。弹头变形程度与实验结果一致,说明数值仿真可以较理想地模拟开花弹的入水过程。
表 2 不同速度入水的弹丸头部变形比较Table 2. Comparison of deformation of projectiles entering water at different velocitiesExperiment Simulation Velocity/(m·s-1) Dmax/mm Dmax/mm Velocity/(m·s-1) Dmax/mm 408 5.10,5.08 5.090 400 5.18 501 6.71,6.79 6.750 500 6.82 610 11.30,11.23 11.265 600 11.48 702 13.82,13.80 13.810 700 13.83 814 13.81,13.81 13.810 800 13.82 3.2 不同速度入水的速度衰减比较
图 5为开花弹在不同速度下撞击水介质的速度衰减曲线,开花弹从空气中射入水后,速度开始降低。从图 5中可以看出,子弹的入射速度越高,速度衰减越明显。在400和500 m/s速度下,子弹速度衰减得比较缓慢,说明低速条件下开花弹头部未发生开裂,变形较小,子弹受到的侵彻阻力低,1 ms后子弹速度分别为182.13和143.58 m/s;600 m/s速度下,子弹速度衰减的速率有所提高,此时开花弹头部已开裂,但仍未变形完全,1 ms后子弹速度为106.36 m/s;700和800 m/s侵彻速度下,开花弹速度衰减十分明显,说明弹体头部完全开裂成“花瓣状”,与水的接触面增大,侵彻阻力高,1 ms后子弹速度分别降低到87.74和75.78 m/s。可以看出,开花弹的入射速度越高,1 ms后子弹的速度反而越低。这是因为开花弹低速入射时,弹头开裂程度小,侵彻过程中弹体最大横截面积小,与水的接触面小,受到的侵彻阻力低,因此速度衰减得较缓慢;随着入射速度的提高,弹头变形增大,与水的接触面增大,开花弹受到的侵彻阻力较大,导致速度衰减得更快。说明开花弹弹头开裂变形对提高子弹的侵彻阻力、降低子弹的速度起很大作用。
3.3 不同速度入水的位移比较
图 6为开花弹在不同速度下撞击水介质的位移曲线。从图 6中可以看出:开花弹以400 m/s入射时,1 ms内在水中的位移是22.13 cm;以500 m/s入射时位移有所增加,达到23.55 cm。这是因为子弹的初速增大,动能变大,而开花弹以500 m/s速度撞击水介质时弹体头部变形很小,矩形槽口发生略微扩张,达不到增大与水的接触面,进而增大侵彻阻力的目的,所以在水中的侵彻深度有所增加。开花弹以600、700和800 m/s速度撞击水介质时,1 ms内在水中的位移分别为20.36、20.55和19.95 cm,高速撞击下子弹在水中的位移明显低于低速撞击下的位移,这是因为开花弹在较高速度下撞击水介质时,弹头发生开裂变形,呈花瓣状,弹体的最大横截面积提高1~2倍,增大了子弹撞击过程中与水的接触面积,极大地提高了子弹在水中受到的侵彻阻力,降低子弹的运动速度,使子弹在高速入射条件下的位移比低速入射时更低。不同速度入射下1 ms后子弹的位移都在20~24 cm范围内,由此可以说明开花弹具有良好的低侵彻特性。
4. 结论
对开花型低侵彻弹进行了不同速度侵彻水介质的实验研究,结果表明:开花弹以500 m/s以下低速撞击时,弹丸头部变形较小,弹形保持良好,稳定侵彻阶段弹头最大直径为弹体直径的1.0~1.2倍;以501~610 m/s中速撞击时,弹头发生开裂,但未变形完全,弹头最大直径达到弹体直径的2倍,受到的侵彻阻力较大;以702~814 m/s高速撞击时,弹头完全开裂,呈花瓣状,弹头最大直径达到弹体直径的2.4倍,受到的侵彻阻力最大。利用LS-DYNA软件模拟开花弹在不同速度下的入水过程,仿真结果与实验结果基本一致,说明数值仿真可以较理想地模拟开花弹的入水过程;得到了子弹的速度衰减曲线和位移曲线,表明开花弹弹头开裂变形对提高子弹的侵彻阻力、降低子弹速度起很大作用,证明开花弹具有良好的低侵彻特性。
-
表 1 混凝土材料的配比
Table 1. Mixture ratio of the concrete material
kg/m3 Stone Sand Water Cement Water reducing agent 950 900 173 500 12.5 表 2 测试相关信息
Table 2. Summary of the test information
No. Temperature/℃ Strain rate/s−1 Peak stress/MPa 1 20 37 52.42 2 20 56 68.41 3 20 76 78.91 4 200 33 50.83 5 200 41 62.42 6 200 49 74.95 7 400 32 66.83 8 400 47 89.29 9 400 56 97.93 表 3 不同温度、应变率条件下的材料参数A和B
Table 3. Material parameters A and B at different temperatures and strain rates
Temperature/℃ Strain rate/s−1 A B 20 37 0.698 240.9 56 0.411 441.9 76 0.364 536.0 200 33 0.597 314.4 41 0.507 388.0 49 0.434 418.1 400 32 0.418 271.7 47 0.369 314.3 56 0.311 412.0 -
[1] CHEN L, FANG Q, JIANG X, et al. Combined effects of high temperature and high strain rate on normal weight concrete [J]. International Journal of Impact Engineering, 2015, 86: 40–56. doi: 10.1016/j.ijimpeng.2015.07.002 [2] SU H, XU J, REN W. Experimental study on the dynamic compressive mechanical properties of concrete at elevated temperature [J]. Materials and Design, 2014, 56: 579–588. [3] 刘传雄, 李玉龙, 吴子燕, 等. 高温后混凝土材料的动态压缩力学性能 [J]. 土木工程学报, 2011, 44(4): 78–83.LIU C X, LI Y L, WU Z Y, et al. Dynamic compression behavior of heated concrete [J]. China Civil Engineering Journal, 2011, 44(4): 78–83. [4] YOUSSEF M A, MOFTAH M. General stress–strain relationship for concrete at elevated temperatures [J]. Steel Construction, 2007, 29(10): 2618–2634. [5] HUO J S, HE Y M, XIAO L P, et al. Experimental study on dynamic behaviours of concrete after exposure to high temperatures up to 700 ℃ [J]. Materials & Structures, 2013, 46(1-2): 255–265. [6] 黄瑞源, 李永池, 章杰, 等. 压剪耦合损伤演化方程在混凝土本构模型中的应用 [J]. 北京理工大学学报, 2013, 33(6): 551–555. doi: 10.3969/j.issn.1001-0645.2013.06.001HUANG R Y, LI Y C, ZHANG J, et al. Application of damage evolution equation under coupled compression-shear to constitutive model of concrete [J]. Journal of Beijing Institute of Technology, 2013, 33(6): 551–555. doi: 10.3969/j.issn.1001-0645.2013.06.001 [7] 王道荣, 胡时胜. 冲击载荷下混凝土材料损伤演化规律的研究 [J]. 岩石力学与工程学报, 2003, 22(2): 223–226. doi: 10.3321/j.issn:1000-6915.2003.02.011WANG D R, HU S S. Study on damage evolution of concrete under impact load [J]. Chinese Journal of Rock Mechanics and Engineering, 2003, 22(2): 223–226. doi: 10.3321/j.issn:1000-6915.2003.02.011 [8] DOUGILL J W. On stable progressively fracturing solids [J]. Journal of Applied Mathematics and Physics, 1976, 27(4): 423–437. [9] LOLAND L E. Continuum damage model for load response estimation of concrete [J]. Cement and Concrete Research, 1980, 10(3): 395–402. doi: 10.1016/0008-8846(80)90115-5 [10] XIE H. Effects of fractal crack [J]. Theoretical & Applied Fracture Mechanics, 1995, 23(3): 235–244. [11] HOLMQUIST T J, JOHNSON G R, COOK W H. A computational constitutive model for concrete subjected to large strains, high strain rates, and high pressures[C]//Proceeding of 14th International Symposium on Ballistics. Quebec, Canada, 1993: 591−600. [12] 王春来, 徐必根, 李庶林, 等. 单轴受压状态下钢纤维混凝土损伤本构模型研究 [J]. 岩土力学, 2006, 27(1): 151–154. doi: 10.3969/j.issn.1000-7598.2006.01.030WANG C L, XU B G, LI S L, et al. Study of damage constitutive model for steel fiber reinforced concrete under uniaxial compression condition [J]. Rock and Soil Mechanics, 2006, 27(1): 151–154. doi: 10.3969/j.issn.1000-7598.2006.01.030 [13] BAI Y L, BAI J, LI H L, et al. Damage evolution, localization and failure of solids subjected to impact loading [J]. International Journal of Impact Engineering, 2000, 24(6): 685–701. [14] 李志武, 许金余, 戴双田, 等. 高温下混凝土冲击加载试验研究 [J]. 高压物理学报, 2013, 27(3): 417–422. doi: 10.11858/gywlxb.2013.03.016LI Z W, XU J Y, DAI S T, et al. Experimental study on concrete exposed to high temperature under impact loading [J]. Chinese Journal of High Pressure Physics, 2013, 27(3): 417–422. doi: 10.11858/gywlxb.2013.03.016 [15] 李业学, 谢和平, 彭琪, 等. 活性粉末混凝土力学性能及基本构件设计理论研究进展 [J]. 力学进展, 2011, 41(1): 51–59.LI Y X, XIE H P, PENG Q, et al. Progress in mechanic property and design theory of elementary struture of reactive powder concrete (RPC) [J]. Advanced in Mechanics, 2011, 41(1): 51–59. [16] OLIVIA M, NIKRAZ H. Properties of fly ash geopolymer concrete designed by Taguchi method [J]. Materials & Design, 2012, 36: 191–198. [17] LI J, WU C, HAO H, et al. A study of concrete slabs with steel wire mesh reinforcement under close-in explosive loads [J]. International Journal of Impact Engineering, 2017, 110: 242–254. doi: 10.1016/j.ijimpeng.2017.01.016 [18] 康亚明, 贾延, 罗玉财, 等. 莫尔-库仑准则下高强度混凝土的临界爆裂蒸汽压力 [J]. 爆炸与冲击, 2018, 38(1): 224–232. doi: 10.11883/bzycj-2016-0305KANG Y M, JIA Y, LUO Y C, et al. Critical vapour pressure for explosive spalling of high-strength concretebased on Mohr-Coulomb criterion [J]. Explosion and Shock Waves, 2018, 38(1): 224–232. doi: 10.11883/bzycj-2016-0305 [19] 李胜林, 刘殿书, 李祥龙, 等. Ф75 mm分离式霍普金森压杆试件长度效应的试验研究 [J]. 中国矿业大学学报, 2010, 39(1): 93–97.LI S L, LIU D S, LI X L, et al. The effect of specimen length in Ф75 mm split hopkinson pressure bar experiment [J]. Journal of China University of Mining & Technology, 2010, 39(1): 93–97. [20] NEMAT-NASSER S, ISAACS J B. Direct measurement of isothermal flow stress of metals at elevated temperatures and high strain rates with application to Ta and Ta, W alloys [J]. Acta Materialia, 1997, 45(3): 907–919. doi: 10.1016/S1359-6454(96)00243-1 [21] 贾彬, 陶俊林, 李正良, 等. 高温混凝土动态力学性能的SHPB试验研究 [J]. 兵工学报, 2009, 30(Suppl 2): 208–212.JIA B, TAO J L, LI Z L, et al. SHPB test for dynamic mechanical performances of concrete at high temperatures [J]. Acta Armamcntarii, 2009, 30(Suppl 2): 208–212. [22] 李世超, 黄瑞源, 唐奎, 等. 一种基于围压和应变率效应的动态本构模型在钢纤维混凝土中的应用 [J]. 兵工学报, 2017, 10(Suppl 1): 66–72.LI S C, HUANG R Y, TANG K, et al. Application of a dynamic constitutive model with confining pressure and strain rate effects in steel fiber reinforced concrete [J]. Acta Armamcntarii, 2017, 10(Suppl 1): 66–72. [23] 胡亮亮, 黄瑞源, 高光发, 等. 混凝土类材料SHPB实验中确定应变率的方法 [J]. 爆炸与冲击, 2019, 39(6): 063102. doi: 10.11883/bzycj-2018-0142HU L L, HUANG R Y, GAO G F, et al. A novel method for determining strain rate of concrete-like materials in SHPB experiment [J]. Explosion and Shock Waves, 2019, 39(6): 063102. doi: 10.11883/bzycj-2018-0142 [24] LEMAITRE J. How to use damage mechanics [J]. Nuclear Engineering & Design, 1984, 80(2): 233–245. [25] 李永池, 曹结东, 董杰, 等. 45钢的损伤演化方程和层裂准则研究 [J]. 力学季刊, 2006, 27(2): 329–334. doi: 10.3969/j.issn.0254-0053.2006.02.023LI Y C, CAO J D, DONG J, et al. Research of damage evolution equation and spallation criterion of 45 steel [J]. Chinese Quarterly of Mechanics, 2006, 27(2): 329–334. doi: 10.3969/j.issn.0254-0053.2006.02.023 [26] 刘俊新, 张可, 刘伟, 等. 不同围压及应变速率下页岩变形及破损特性试验研究 [J]. 岩土力学, 2017(Suppl 1): 49–58.LIU J X, ZHANG K, LIU W, et al. Experimental study of mechanical behaviours of shale under different confining pressures and different strain rates [J]. Rock and Soil Mechanics, 2017(Suppl 1): 49–58. -