Numerical Simulation of Disturbance by Double-Layer Explosive Reactive Armor with Wedged Flying-Plate on Jet
-
摘要: 为得到干扰聚能射流能力更好的爆炸反应装甲,在经典爆炸反应装甲的基础上,设计了一种双层楔形飞板爆炸反应装甲。利用ANSYS/LSDYNA-3D仿真软件对3种不同方案进行了模拟计算,分别对各方案中飞板飞行形态、逃逸射流特性、射流的动能变化以及聚能射流对靶板的侵彻深度进行了分析。结果表明:夹层炸药引爆后,楔形飞板在向外飞出的同时具有一定的旋转特征,合理的摆放结构能够增大飞板与射流的作用面积;聚能射流在穿过反应装甲后,动能急剧下降,穿深能力降低,方案二聚能射流侵彻深度最浅,方案三次之,方案一最深,表明方案二具有良好的防护效果。对楔形飞板的研究丰富了爆炸反应装甲的结构设计,为反应装甲的进一步研究提供了理论参考。Abstract: Based on the traditional explosive reactive armor (ERA), we designed a structure of double-layer ERA with wedged flying-plates, to achieve the best performance in jet disturbing.By using the simulation software ANSYS/LSDYNA-3D, we simulated 3 different schemes of ERA structures, and analyzed the flying states of flying-layers, the character of escaped jet, the kinetic energy of jet and the penetrating depth of target by jet separately.The results show that the flying-layer tended to rotate after the sandwich explosive were detonated, and a proper structure arrangement helped to enlarge the active area between the jet and the flying-layer.The kinetic energy of jet declined greatly after across the ERA, resulting in a lower penetrating depth.Of the 3 schemes, Scheme 2 has the smallest penetrating depth by jet into the target, the next is Scheme 3, while Scheme 1 has the maximum penetrating depth, indicating that the structure of Scheme 2 has the best protection effect.This study on the wedged double-layer ERA enriched the structure design of ERA, providing theoretical reference for its further research.
-
Key words:
- wedged flying-layer /
- explosive reactive armor /
- shaped charge jet /
- penetration
-
作为液体和气体长距离运输的一种重要方式,管道运输在国家经济发展和国民生活中发挥着重要作用。然而,随着城镇化进程的加速,城市管网系统密集分布,爆炸作用引起的管道安全问题受到国内外广泛关注[1-3]。都的箭等[4]通过实验研究发现,正对爆心管段背面受到很大的轴向拉应力作用,且管道受爆炸载荷的影响主要与爆心距有关。Ji等[5]研究了X70钢管在局部爆炸载荷下的动力响应,发现管道的挠度和损伤程度随炸药量和接触面积的增大而增大,且壁厚对管道损伤和失效后的运动有重要作用。数值模拟是研究爆炸问题的一种重要方法,只要方法得当,模拟效果可与实际情况相吻合[6-7]。为此,梁政等[8]利用数值模拟方法研究了管道埋深、药量和管道壁厚因素对爆炸载荷下的埋地管道动力响应的影响。房冲[9]通过模拟研究发现,在爆炸载荷下充水管道的变形量、位移和峰值压强都比内空管道小。余洋等[10]采用野外实验与数值计算相结合的方法研究了初始条件对钢质方管在侧向局部爆炸载荷作用下损伤破坏效应的影响。
迄今为止,对爆炸载荷作用下焊缝区附近埋地钢管的动力响应的相关研究鲜有报道。基于此,以两种含Y型焊缝(坡口有2 mm余高焊缝和坡口无余高焊缝)的埋地X70钢管为例,采用有限元软件ANSYS/LS-DYNA,数值模拟研究爆炸载荷作用下焊缝区附近埋地X70钢管的动力响应规律,以期为埋地管线附近的爆破施工设计和埋地管线的安全防护提供一定的理论参考。
1. 有限元模型
1.1 计算模型
采用cm-g-μs单位制,建立由TNT炸药、黄土和焊接管道组成的计算模型,如图1所示。模型纵向长38.4 cm,管道中心到模型侧面的宽度为130.0 cm,模型整体高271.6 cm,其中:TNT炸药为边长14.0 cm的正方体,采用中心起爆方式;焊接管道为外径1 016.0 mm、壁厚14.6 mm的X70钢管。焊缝选取两种尺寸[11],分别为Y型坡口有余高(H = 2.0 mm)焊缝和Y型坡口无余高(H = 0)焊缝,如图2所示。为了提高计算收敛速度,将焊缝尺寸设计图进行适当的简化,简化模型如图3所示。两种焊缝均不考虑分层焊接工艺的影响,焊缝与管道采用共节点方式连接。
考虑到计算模型的对称性,取1/2模型建模。炸药、黄土、管道及焊缝选用SOLID164六面体实体单元,用扫掠方式划分网格,并对焊缝位置进行网格细化处理。炸药和黄土采用欧拉网格,焊接管道和焊缝采用拉格朗日网格,运用任意拉格朗日-欧拉算法及管土间流固耦合算法模拟爆炸载荷作用下埋地焊接管道的动力响应。在土体外侧和底面设置透射边界条件,模型对称面施加对称约束。
1.2 计算工况
为了初步揭示爆炸载荷作用下两种焊缝形式的埋地焊接管道的动力响应规律,选取药包尺寸为14.0 cm × 14.0 cm × 14.0 cm的TNT炸药,对埋深为1.5 m的焊缝有余高(H = 2.0 mm)管道(管道A)和焊缝无余高(H = 0)管道(管道B),在炸高分别为60.0、85.0和110.0 cm的3种条件下的6种工况进行模拟计算,如表1所示,其中,hB为炸高。
表 1 计算工况Table 1. Calculation conditionsWeld type Buried depth of pipeline/m Size of TNT/(cm × cm × cm) hB/cm No weld reinforcement (H = 0) 1.5 14.0 × 14.0 × 14.0 60.0, 85.0, 110.0 Weld reinforcement (H = 2.0 mm) 1.5 14.0 × 14.0 × 14.0 60.0, 85.0, 110.0 1.3 材料参数
TNT炸药选用高能炸药模型(Mat_High_Explosive_Burn)和JWL状态方程定义。JWL状态方程表达式为
pz=A(1−ωR1ν)e−R1ν+B(1−ωR2ν)e−R2ν+ωEν (1) 式中:pz为爆炸产物的压力,A、B、R1、R2、ω为TNT材料常数,v为爆炸产物的相对比容,E为炸药初始内能。炸药密度ρz、爆速D以及JWL状态方程参数见表2[12]。
黄土选用泡沫模型(Mat_Soil_and_Foam)描述。该材料模型的应力屈服常数f为
f=Sijδij/2−(a0+a1pt+a2pt2) (2) 式中:Sij为土体材料的Cauchy偏应力张量,δij为土体材料的Kronecker系数,a0、a1、a2分别为土体摩擦角、土体黏聚力和土体爆炸动载效应的影响系数,pt为土体压力。a0、a1、a2由土工实验测得的内摩擦角和土壤黏聚力参数确定,土体密度ρt、剪切模量G、体积模量K等参数见表3[13-14]。
X70钢管道和焊缝均采用双线性随动材料模型(Mat_Plastic_Kinematic)描述,遵循von Mises屈服准则,其表达式为
σ={Esεε⩽εeσy+Et(ε−εe)ε>εe (3) 式中:σ为应力;σy为屈服应力;Es为弹性模量;Et为切线模量,0 < Et < Es;ε为应变;εe为弹性极限应变。管道和焊缝的具体材料参数见表4[15-18],其中,μ为泊松比。
2. 结果与讨论
2.1 管道应力分析
图4为边长14.0 cm的正方体TNT炸药爆炸时,炸高hB为60.0 cm,埋深为1.5 m的两种X70管道焊缝附近的von-Mises应力云图。由图4可以看出:当传播时间为1 440 μs时,爆炸应力波阵面已经接触管道;当传播时间为1 600 μs时,焊缝有2.0 mm余高的管道A和焊缝无余高的管道B的最大应力增幅分别为81.4 MPa和43.0 MPa;当传播时间为1 920 μs时,管道A和管道B的最大应力均大于焊缝与管道的材料屈服应力,且应力沿迎爆面正对爆心位置向外扩展,其中管道A的应力呈“十”字形扩展,而管道B的应力以椭圆形向四周扩展;在3 520 μs时,管道应力集中主要沿裂缝位置发展,管道A和管道B的应力最大值分别为601.2 MPa和591.0 MPa;在6 080 μs时,管道A和管道B继续变形但应力减小,应力最大值分别减小到581.8 MPa和565.8 MPa;在9 120 μs时,管道A和管道B的应力集中基本消失。
图5和图6分别为两种管道外表面上正对爆心位置的焊缝与管道分界面处焊缝单元与管道单元的应力时程曲线。在管道受爆炸载荷作用阶段,两种管道的应力在大约480 μs内呈跳跃式上升。其主要原因是管道为瞬时受力,一部分爆炸能量使管道变形并向管道四周传递,导致焊缝与管道分界面处两个典型单元的应力呈降低趋势,此现象与图4的应力云图吻合。根据应力集中系数和余高关系的经验公式[19]可得:管道A和管道B的应力集中系数分别为1.016和1.008,即随着余高增大,应力集中系数逐渐增大。对比图5和图6可知,管道A的焊缝单元应力峰值较高,应力下降趋势相对较陡。这也说明焊缝余高的存在使得焊缝与管道分界面的截面尺寸突变增大,从而导致焊缝有余高的焊接管道受应力集中的影响较大。在1 912 μs时,图5(管道A)和图6(管道B)的焊缝单元应力最大,分别约为560.0 MPa(焊缝的屈服强度为550 MPa)和545.6 MPa。同时,管道A的焊缝处首先达到管道屈服强度(480 MPa),按照von Mises屈服准则,管道A开始进入局部塑性变形阶段,此时管道B的应力尚未达到材料的屈服强度。
2.2 管道位移分析
埋地X70管道的迎爆面和背爆面的最大位移如表5所示。从表5可知,由于爆炸冲击波的一部分能量在土中被耗散,且随着爆炸冲击波在土中传播距离的增大,两种焊缝形式管道的迎爆面和背爆面的最大位移均呈现减小的趋势。当炸高hB从60.0 cm增加到85.0 cm以及从85.0 cm增加到110.0 cm时,管道A和管道B迎爆面的最大位移减小量分别为2.303 cm、0.715 cm和2.300 cm、0.572 cm,而管道A和管道B背爆面的最大位移减小量分别为0.391 cm、0.235 cm和0.373 cm、0.280 cm。两种焊缝形式管道迎爆面的最大位移减小量大于背爆面,这是由于爆炸冲击载荷在土中传播后直接作用于管道迎爆面,对管道迎爆面产生的影响较大,土体对管道背爆面具有一定的支撑作用,从而减小了管道背爆面位移。在相同炸高下管道A比管道B的最大位移大,且在炸高为60.0、85.0和110.0 cm时,两种焊缝形式管道迎爆面的最大位移差值分别为0.270、0.267和0.124 cm,即随着炸高的增大,两种焊缝形式的埋地焊接管道最大位移的差值逐渐减小,也说明当炸高较小时,管道A整体抵抗变形的能力弱于管道B。然而,随着炸高的增大,作用于管道的能量减小[20],管道塑性变形较小,使得这种现象逐渐模糊。
表 5 埋地X70管道的迎爆面和背爆面的最大位移Table 5. Maximum displacement of explosion-front and explosion-back surfaces of buried X70 pipelineTypes of weld hB/cm Maximum displacement/cm Explosion-front surface Explosion-back surface Weld reinforcement (H = 2.0 mm) 60.0 5.482 0.846 85.0 3.179 0.455 110.0 2.464 0.220 No weld reinforcement (H = 0) 60.0 5.212 0.943 85.0 2.912 0.570 110.0 2.340 0.290 2.3 管道等效应变分析
表6为两种不同类型焊缝的埋地焊接管道在不同炸高下的最大等效应变统计。从表6可知,管道A和管道B的最大等效应变均随炸高的增大而减小。当炸高从60.0 cm增大到85.0 cm时,管道A和管道B的最大等效应变分别减小约58.12%和61.13%;当炸高从85.0 cm增大到110.0 cm时,管道A和管道B的最大等效应变分别减小约45.92%和38.05%,在炸高相同时,管道A的最大等效应变大于管道B,且管道A的最大等效应变位于焊缝余高表面,而管道B的最大等效应变在焊缝与管道处一定范围内沿纵向分布。这在一定程度上说明管道B能更好地协调焊缝与管道分界处的应变,有利于保障焊缝与管道的局部协同变形性能。
表 6 不同炸高时埋地X70管道的最大等效应变Table 6. Maximum effective strain of buried X70 pipeline with different blasting heightsTypes of weld hB/cm Peak effective strain/10–3 Weld reinforcement
(H = 2.0 mm)60.0 9.937 85.0 4.162 110.0 2.251 No weld reinforcement
(H = 0)60.0 6.877 85.0 2.673 110.0 1.656 2.4 管道振速分析
表7为不同炸高下两种焊缝形式管道的迎爆面和背爆面处焊缝位置中心单元X方向的最大振动速度。从表7可以看出,两种焊缝形式管道的迎爆面和背爆面的最大振动速度均随着炸高增大而减小,且迎爆面的最大振速均大于背爆面。这说明迎爆面受爆炸地震波的影响较大。当炸高hB为60.0、85.0和110.0 cm时,管道B的迎爆面的最大振动速度较管道A大,迎爆面差值分别为1.600、0.539和0.329 m/s,而背爆面差值在0.200 m/s以内。管道峰值速度随着管壁厚度的增大而减小[12],由于管道A增加了管道在焊缝位置的径向厚度,可将其视为管道焊缝位置的壁厚增大导致管道A的峰值振速减小。这说明管道A抵抗爆炸振动的性能优于管道B,且在炸高为60.0 cm时,管道A抵抗振动性能的优势较为明显。
表 7 埋地X70管道的迎爆面和背爆面最大振速Table 7. Maximum vibration velocity of explosion-front and explosion-back surfaces of buried X70 pipelineTypes of weld hB/cm Maximum vibration velocity/(m·s–1) Explosion-front surface Explosion-back surface Weld reinforcement (H = 2.0 mm) 60.0 22.748 4.431 85.0 9.316 2.817 110.0 4.503 1.693 No weld reinforcement (H = 0) 60.0 24.348 4.294 85.0 9.855 2.867 110.0 4.832 1.746 图7为不同炸高时两种焊缝形式的管道典型单元的速度时程曲线。当炸高hB分别为60.0、85.0和110.0 cm时,管道A和管道B达到最大振速的时间分别为2 560 μs和2 560 μs、4 500 μs和4 600 μs、7 200 μs和7 200 μs,两种焊缝形式的管道达到最大振速的时间差值均在100 μs以内。这说明两种焊缝形式的管道达到最大振速的时间主要受炸高的影响,受焊缝形式的影响较小。
3. 结 论
(1)当炸高为60.0 cm时,两种焊缝形式的埋地X70焊接管道在爆炸载荷作用下焊缝位置均出现应力集中,但焊缝有余高的管道受应力集中影响较大,且会先于焊缝无余高管道进入屈服阶段。
(2)当炸高为60.0~110.0 cm时,由于爆炸载荷直接作用于迎爆面,且管土间的相互作用对管道背爆面具有一定的支撑作用,两种焊缝形式管道迎爆面的最大位移均大于背爆面的最大位移。当炸高为60.0、85.0 cm时,焊缝有余高的管道整体抵抗变形的能力明显弱于焊缝无余高的管道。
(3)焊缝无余高管道较焊缝有余高管道在焊缝与管道分界处的应变更为协调,能更好地保障焊缝与管道的局部协同变形性能。
(4)在相同的爆炸载荷下,焊缝有余高管道抵抗振动的性能优于焊缝无余高管道。药量相同条件下,相对于焊缝形式,炸高对含焊缝区管道的最大振速起主要作用。
-
表 1 B炸药材料模型及状态方程参数
Table 1. Parameters of material model and equation of state of explosive B
ρ/(g·cm-3) pCJ/GPa D/(m·s-1) AJWL/GPa BJWL/GPa RJWL1 RJWL2 E0 1.717 29.5 7 980 524.2 7.678 4.2 1.1 0.085 表 2 铜、钢材料模型及状态方程参数
Table 2. Parameters of material model and equation of state of copper and steel
Material ρ/(g·cm-3) G/GPa AJ-C/GPa BJ-C/GPa S1 S2 C/(m·s-1) V0 Copper 8.96 46 0.09 0.92 1.489 0 3 940 1 Steel 7.785 77.5 0.175 0.376 1.49 0 4 570 1 表 3 夹层炸药的材料模型及状态方程参数
Table 3. Parameters of material model and equation of state of sandwich explosive
ρ/(g·cm-3) G/GPa σY/GPa A/GPa B/GPa R1 R2 R3 1.712 3.54 0.2 524.2 7.678 778.1 -5.031×10-2 2.223×10-5 R5 R6 cP/(J·kg-1·K-1) cR/(J·kg-1·K-1) GROW2 AR2 ES1 ES2 11.3 1.13 10-3 2.487×10-3 300 1 0.222 0.333 -
[1] 王凤英.装甲防护技术的发展[J].测试技术学报, 2002, 16(2):144-147. http://www.cqvip.com/QK/98543X/200202/6495303.htmlWANG F Y.The development of amor fence[J].Journal of Test and Measurement Technology, 2002, 16(2):144-147. http://www.cqvip.com/QK/98543X/200202/6495303.html [2] 王凤英, 岳继伟, 王志远, 等.多层变角度反应装甲对聚能射流的干扰作用[J].高压物理学报, 2017, 31(5):566-572. doi: 10.11858/gywlxb.2017.05.009WANG F Y, YUE J W, WANG Z Y, et al.Interference effect of multi-layerd reactive armor with a variable angle on shaped charge jet[J].Chinese Journal of High Pressure Physics, 2017, 31(5):566-572. doi: 10.11858/gywlxb.2017.05.009 [3] 叶春辉, 刘天生.影响爆炸反应装甲两飞板运动过程的因素分析[J].科技信息, 2012(28):117. http://mall.cnki.net/magazine/Article/KJXX201228102.htmYE C H, LIU T S.Analysis of factors affecting explosive reactive armor moving process of two flying plates[J].Science and Technology Information, 2012(28):117. http://mall.cnki.net/magazine/Article/KJXX201228102.htm [4] 甄金朋, 刘天生, 张硕, 等.爆炸反应装甲驱动飞板运动的数值模拟[J].火炸药学报, 2010, 33(2):78-81. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hzyxb201002020ZHEN J P, LIU T S, ZHANG S, et al.Numerical simulation of the flyer plate propelled by explosive reactive amor[J].Chinese Journal of Explosive and Propellants, 2010, 33(2):78-81. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hzyxb201002020 [5] 李歌, 王凤英, 刘天生.爆炸反应装甲干扰穿甲弹的试验研究[J].火炸药学报, 2010, 33(2):50-52. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=bgxb201002013&dbname=CJFD&dbcode=CJFQLI G, WANG F Y, LIU T S.Experimental study on explosive reactive amor disturbing penetrators[J].Chinese Journal of Explosive and Propellants, 2010, 33(2):50-52. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=bgxb201002013&dbname=CJFD&dbcode=CJFQ [6] 黄鹤, 刘炜, 陈建申.爆炸反应装甲含能材料的发展及探索[J].兵工自动化, 2013, 32(1):55-57. doi: 10.7690/bgzdh.2013.01.017HUANG H, LIU W, CHEN J S.Explosive reactive armor of energetic materials development and exploration[J].Ordnance Industry Automation, 2013, 32(1):55-57. doi: 10.7690/bgzdh.2013.01.017 [7] 毛东方, 李向东, 宋柳丽.V型夹层炸药对射流干扰的数值模拟[J].爆炸与冲击, 2008, 28(11):86-91. http://www.cqvip.com/Main/Detail.aspx?id=26640867MAO D F, LI X D, SONG L L.Numerical simulation of disturbance by sandwich explosive on jet[J].Explosion and Shock Waves, 2008, 28(11):86-91. http://www.cqvip.com/Main/Detail.aspx?id=26640867 [8] 刘宏伟, 夏松林, 赵靖.V形反应装甲与射流作用过程分析[J].兵器材料科学与工程, 2011, 34(4):20-22. http://www.cqvip.com/QK/95120X/201104/38832989.htmlLIU H W, XIA S L, ZHAO J.Interaction process between jet and V-shaped double ERA[J].Ordnance Material Science and Engineering, 2011, 34(4):20-22. http://www.cqvip.com/QK/95120X/201104/38832989.html [9] 吴成, 蒋建伟, 冯顺山, 等.爆炸反应装甲运动规律的数值仿真及研究[J].兵工学报, 2002, 23(1):35-38. http://d.wanfangdata.com.cn/Periodical_bgxb200201009.aspxWU C, JIANG J W, FENG S S, et al.A study on the moving features of explosive reactive armor by numerical simulation and experiments[J].Acta Armamentarii, 2002, 23(1):35-38. http://d.wanfangdata.com.cn/Periodical_bgxb200201009.aspx [10] PAIK S H, KIM S J, YOO Y H, et al.Protection performance of dual flying oblique plates against a yawed long-rod penetrator[J].International Journal of Impact Engineering, 2007, 34(8):1413-1422. doi: 10.1016/j.ijimpeng.2006.06.006 [11] MAYSELESS M.Effectiveness of explosive reactive armor[J].Journal of Applied Mechanics, 2011, 78(5):1-11. http://adsabs.harvard.edu/abs/2011JAM....78e1006M [12] 黄正祥.聚能装药理论与实践[M].北京:北京理工大学出版社, 2014:376-387.HUANG Z X.Theory and practice of shaped charge[M].Beijing:Beijing Institute of Technology Press, 2014:376-387. [13] 武海军, 陈利, 王江波, 等.反应装甲对射流干扰的数值模拟研究[J].北京理工大学学报, 2006, 26(7):565-568, 605. http://industry.wanfangdata.com.cn/yj/Detail/Periodical?id=Periodical_bjlgdxxb200607001WU H J, CHEN L, WANG J B, et al.Numerical simulation on reactive armor disturbing jet[J].Transactions of Beijing Institute of Technology, 2006, 26(7):565-568, 605. http://industry.wanfangdata.com.cn/yj/Detail/Periodical?id=Periodical_bjlgdxxb200607001 [14] 姬龙, 黄正祥, 顾晓辉.双层楔形爆炸反应装甲飞板的运动规律[J].爆炸与冲击, 2013, 33(4):387-393. doi: 10.11883/1001-1455(2013)04-0387-07JI L, HUANG Z X, GU X H.Analysis and experimental study on the explosive field of double-layer explosive reactive armor[J].Explosion and Shock Waves, 2013, 33(4):387-393 doi: 10.11883/1001-1455(2013)04-0387-07 -