单位文秘网 2021-07-23 08:17:07 点击: 次
摘要: 研究战场环境下爆炸冲击过程及其对装备部件的毁伤机理,对于武器系统改进和战斗损伤评估具有重要意义。本文利用流固耦合算法进行爆炸冲击模拟,探讨了冲击波形成机理,确定了将爆炸冲击波直接毁伤靶板的方法,以及冲击波的数值定义和加载方法,并利用流固耦合算法进行了爆炸冲击仿真,其方法可以为爆炸冲击与毁伤研究提供一定参考。
Abstract: Research on explosion impact and its damage mechanism for equipment components under battlefield environment, are important for improving weapons systems and assess battle damage. In this paper, the fluid-structure interaction algorithm is used for explosive shock simulation, shock wave formation mechanism is discussed, the method of blast directly damage target, as well as numerical definition and loading method of the shock waves are defined, the fluid-structure interaction algorithm is used in simulating explosion and shock. The method can provide reference for researching blast shock and its damage.
关键词: 爆炸冲击;毁伤机理;流固耦合算法;冲击载荷
Key words: explosion and shock;damage mechanism;fluid-structure interaction algorithm;impact load
中图分类号:TH16 文献标识码:A 文章编号:1006-4311(2013)11-0023-02
0 引言
由于爆炸冲击波产生的威力巨大,并且可以通过计算炸药的各项参数得到产生爆炸冲击波的各项数据,能够较好的得到控制和利用,所以爆炸冲击问题一直是工程和军事领域的热点问题,越来越受到科研工作者的重视。
北京理工大学的赵文杰[1]通过实验和理论分析,对圆抛物面薄壳型雷达天线在爆炸冲击波作用下的变形特点及有关的变形特征参量进行了分析讨论,并给出了工程计算式。中国矿业大学力学与建筑工程学院的李清[2]进行了爆炸载荷作用下动态裂纹扩展实验研究。洛阳水利工程技术研究所[3]做了冲击波对工程结构及装备的动载实验,另外,有些理论已在工程实际中得到了验证与运用,如地下隧道试验[4]、地形对冲击波影响试验[5]及人防工事试验[6]。对装备部件进行爆炸冲击与毁伤模拟,属于结构分析的范畴,本文采用流固耦合算法对爆炸冲击问题进行仿真分析。
1 爆炸冲击波形成机理分析
炸药爆炸完成后,产生的高温、高压爆轰产物像一个超音速推进的活塞,把空气从原来的位置上迅速地排挤出去,形成一个以超音速运动的、状态参数有突跃的压缩空气层,即空气冲击波。在此过程中,爆轰产物的压力在未降到空气初始状态压力前继续膨胀,不断把能量传输给空气冲击波。当爆轰产物的压力下降到空气初始状态压力时,空气冲击波形成,此时爆轰产物的体积称为极限体积。由于惯性作用,爆轰产物体积继续膨胀,使压力低于周围空气压力,所以在空气冲击波尾部会形成一个稀疏区。
爆炸冲击波的形成和传播如图1所示,T1、T2、T3、T4分别表示爆炸后不同瞬间冲击波的状态,其中T3时刻冲击波开始与爆轰产物分离,T4时刻出现了稀疏区。
2 冲击波毁伤靶板的有限元分析
利用有限元方法对靶板和装备部件进行爆炸冲击毁伤仿真,采用如图2所示的步骤。
2.1 冲击波模型 近距离爆炸时,爆炸产物和爆炸冲击波共同对目标进行毁伤,如果建立炸药模型使之产生冲击波,则无法将爆炸产物的作用分离出来;另外,目标与炸药距离扩大,会使建模空间扩大,给计算带来不便,容易造成较大误差。为此,本文中采用了对靶板直接施加冲击波载荷的方法来研究冲击波对靶板的毁伤效果。
实际情况下的冲击波载荷相当复杂,一般和结构变形相耦合且不易精确测定,为研究方便,略去荷载中的次要因素,对实际载荷进行理想化处理,冲击载荷简化后的主要形式有:正弦脉冲、三角脉冲和矩形脉冲。针对爆炸冲击波的特点,有些学者将爆炸冲击波简化为三角脉冲波的形式。Jones在讨论船舶板的抨击损坏时采用Ochi提出的三角形压力脉冲形式的抨击时间历程曲线[7]。
为使问题的分析过程得以简化,本文将爆炸冲击波取为三角形脉冲载荷的形式,载荷-时间历程曲线如图3所示。
如图3所示,爆炸冲击波作用时间为t,在t/2时,冲击波峰值超压达到最大,由于空气中爆炸冲击波的强度一般在80~130MPa,所以峰值超压在此范围内取值,本文选用峰值超压为80MPa和100MPa的冲击波对靶板进行冲击,以研究其毁伤机理。
2.2 靶板材料模型 由于在战损仿真分析中,一般将材料等效为LY-12铝合金,为使仿真分析结果的一致性,本文中的靶板材料也定为LY-12铝合金,其材料参数如表1所示[8]:
LY-12铝合金为应变率无关材料,因此,可选择*MAT_PLASTIC_KINEMATIC材料模型进行有限元分析。
不同应变率下,LY-12铝合金的损伤应力介于0.55至0.56 GPa之间,因此将损伤应力取为0.555GPa,由弹性模量E=72GPa,可计算出该材料的失效应变为0.77%。其他材料参数为:屈服极限γ=0.31Gpa,硬化系数β=0.7。
2.3 冲击载荷定义和施加 冲击载荷的定义包括爆炸冲击波作用时间的定义和强度的定义,且强度随时间变化而变化。为了对靶板施加爆炸冲击载荷,需定义载荷—时间变量数组,赋值情况见表2。冲击载荷定义完成后,对靶板进行冲击载荷施加。
2.4 网格划分与网格密度 为确定合适的网格密度,在不同的靶板网格密度下,用强度100Mpa、作用时间0.00500s的爆炸冲击波对LY-12铝合金靶板(靶板尺寸300×300×8 mm)进行冲击实验,结果如表3所示:
由表3可以看出,随着靶板网格密度的增大,计算时间显著增长,但靶板的最大位移基本趋于稳定。为此,在后续研究中选用了靶板长和宽均50等分,厚度4等分的划分方法。
3 基于流固耦合算法的爆炸毁伤仿真
对装备部件进行毁伤仿真分析,实质上就是利用有限元分析中的流固耦合方法,分析爆炸复合物和冲击波对部件的损伤情况,研究其损伤机理。多物质流固耦合算法的原理是,对炸药及其他流体材料(如空气、水、土壤等)采用Euler算法,对其他的结构采用Lagrange算法,然后通过流固耦合方式来处理相互作用(*CONSTRAINED_LAGRANGE_IN_SOLID),该方法的优点是炸药和流体材料在Euler单元中流动,不存在单元的畸变问题,并且通过流固耦合方式来处理相互作用,能方便地建立爆炸模型。
本文利用TNT炸药近距离爆炸进行毁伤分析,将所有物体都选择为三维实体单元,选定使用LS-DYNA 中提供的SOLID164实体单元。
LS-DYNA专门为各种炸药提供了材料模型(*MAT_HIGN_EXPLOSIVE_BURN),再结合JWL状态方程对它进行描述。TNT炸药的材料模型*MAT_HIGN_EXPLOSIVE_BURN中的相关参数为:密度R0=1.7g/cm3,爆速D=0.753cm/μs,Chapman-Jouget压力PCJ=0.255×1011Pa。
JWL状态方程用于描述压力与体积应变之间的关系,其形式为
P=A1-■e■+B1-■e■+■
对于TNT炸药,在g-cm-μs单位制中,以上方程的输入参数为:A=5.4094,B=0.093726,R1=4.5,R2=1.1,ω=0.35。
4 爆轰产物的破坏作用分析
大量实验表明,距离爆炸中心10~15r0(r0为装药半径)的范围内,目标受到爆炸产物和爆炸冲击波的复合作用,而超过10~15r0范围时,目标只受到爆炸冲击波的作用[9]。
基于上述流固耦合方法进行爆炸冲击仿真,以装药密度为1.6g/cm3的TNT炸药为例,爆炸完成瞬间爆轰产物的各项参数如下:压强=196MPa,密度=2.13g/cm3,温度=3350℃,速度=1750m/s。
从这些数据可以看出,爆炸产物的各项参数都变得很大,爆炸发生瞬间,目标在10~15 r0范围内时,会受到很大的冲量,导致目标严重毁伤。
从以上分析可以看出,爆炸产物的破坏威力大,但作用范围小,只能对直接接触或近距离接触的目标发生作用。不过随着现代高精武器的出现,能够使炸药在距目标距离很小时产生爆炸,所以利用爆炸产物的巨大毁伤作用对目标进行毁伤打击,已成为现实。所以研究炸药近距离爆炸时对目标的毁伤情况,具有重要意义。
参考文献:
[1]赵文杰,蒋浩征,王秀兰.爆炸冲击波对圆抛物面薄壳雷达天线毁伤效应研究[J].弹箭与制导学报,2000,(1):8~13.
[2]李清,杨仁树等.爆炸载荷作用下动态裂纹扩展试验研究[J].岩土力学与工程学报,2005(8),24(16):2913~2916.
[3]张正宇,赵根等.三峡三期碾压混凝土围堰拆除爆破研究[J].工程爆破,2003(3),9(1):1~8.
[4]唐海,李海波等.地形地貌对爆破振动波传播的影响实验研究[J].岩石力学与工程学报,2007(9),26(9):1818~1823.
[5]郭涛,徐全军等.钢筋混凝土地下人防工事爆破拆除[J].工程爆破,2004(3),10(1):35~37.
[6]姚熊亮.舰船结构振动冲击与噪声[M].国防工业出版社. 2007(2).
[7]Jones H.Plastic Behavior on Ship Structures.Trans.SNAME.1976,84:115~145.
[8]Finite Element Modeling of the RAH-66 Comanche Helicopter Tailcone Section Using PATRAN and DYTRAN. ADA392096,2001.
[9]徐安桃,刘祥凯,袁洪涛,甘茂治.装备战斗损伤研究现状及发展趋势[J].车辆与动力技术,2001,(3):57~62.
(责任编辑:单位文秘网) )地址:https://www.kgf8887.com/show-239-70624-1.html
上一篇:铜与钢的焊接
下一篇:侵彻过程中的墙角效应
版权声明:
本站由单位文秘网原创策划制作,欢迎订阅或转载,但请注明出处。违者必究。单位文秘网独家运营 版权所有 未经许可不得转载使用