0 引言

含铝炸药作为一种高性能非理想炸药,其非理想组分能够在理想成分爆轰之后释放大量的能量,使得爆轰产物的温度和压力维持较长的时间而不至于过快衰减[1]。非理想含铝炸药水中爆炸能量释放形式直接影响了冲击波能与气泡能的分配,产生了特定结构的爆炸冲击波与气泡载荷,从而成为水雷、鱼雷等水中兵器常用的主装药[2-3]。RS211炸药(配方体系为RDX/TNT/Al/WAX)是水下非理想含铝炸药典型代表,其具有比TNT炸药更高的水下爆炸能量[4]。水中兵器装填RS211炸药,在国内外已有多年历史[5]。

近年来,随着计算机软硬件技术的发展,数值模拟已经成为含铝炸药爆炸性能预估的重要方法[6-7],数值模拟技术在含铝炸药水下爆炸载荷研究中发挥着重要作用。基于LS-DYNA软件平台,辛春亮等[8]利用一维ALE算法,开展了TNT炸药水中爆炸过程模拟,获得了冲击波与气泡半径时程曲线。Jing Ping Lu等[9]拟合了PBXW-115炸药爆炸产物的JWL状态方程和Lee-Tarver点火生长模型参数,通过二维拉格朗日算法计算了不同药量PBXW-115水中爆炸冲击波衰减演化过程。安丰江等[10]利用二维与三维ALE映射方法,计算分析了TNT炸药水爆冲击波、气泡脉动及水射流等载荷的时空特性。

上述水下爆炸载荷模拟多采用的是LS-DYNA软件中传统算法,对于边界与环境条件简化较多,无法高效可靠的实现全时空域水中爆炸过程连续一体化仿真;描述含铝非理想炸药P-V关系多使用基于理想C-J爆轰理论的JWL状态方程,不能准确反映铝粉的后燃效应;针对RS211炸药二次能量释放过程对水中爆炸载荷参数影响的数值模拟研究相对较少,缺乏可工程借鉴的非理想炸药模型参数。

本文中利用LS-DYNA软件的S-ALE算法,开展了基于JWL-Miller状态方程描述的非理想炸药水中爆炸载荷数值模拟,分析了炸药后燃效应对水中爆炸载荷的影响;利用正交试验法获得了Miller后燃模型参数对水中爆炸载荷特性影响的敏感性规律;通过典型算例结果与试验结果对比研究,匹配设置了RS211含铝炸药的Miller模型参数,并完成了模型的进一步验证与推广应用。

1 模型建立

1.1 数值算法

从宏观时间尺度上看,爆炸冲击波传播与气泡演化过程存在较大差异。非理想炸药全时空域水中爆炸过程连续一体化模拟的物理求解时间较长,数值计算资源的绝对消耗较大。水中爆炸载荷具备强间断与强非线性属性[11],通过数值计算捕捉冲击波压力与气泡脉动周期参数应构建较为精细化的计算模型。因此,高效高精度数值算法一直是水中爆炸载荷数值计算的现实需要。

LS-DYNA程序新发展的S-ALE算法,模型网格由关键字创建,省去了用户创建网格到程序读入的麻烦,同时也节省了大量读写操作带来的运算时间和内存需要[12]。极大地提高了计算模型的前处理及求解效率。与传统ALE算法相比,S-ALE算法计算效率提升约40%。因此,将S-ALE算法应用于水中爆炸过程仿真,是平衡计算精度与计算效率的有效技术途径。S-ALE算法实现水下爆炸过程模拟的核心关键字有[13]：

1) ALE_STRUCTURED_MESH:用于计算域的正交网格定义;

2) INITIAL_HYDROSTATIC_ALE:用于水域静水压力初始化;

3) ALE_AMBIENT_HYDROSTATIC:用于水域环境边界条件定义;

4) INITIAL_VOLUME_FRACTION_GEOMETRY:用于计算域中水、炸药及空气等材料的初始化填充。

1.2 有限元模型

以文献[14]中5 kg球形RS211炸药在水下5 m深处爆炸为参考工况条件,炸药采用中心点起爆方式。考虑一维计算模型对边界条件的简化太多,三维计算模型网格规模庞大,故优选二维轴对称方法构建计算模型。依据试验工况,设置计算水域宽度20 m,深度15 m,空气域厚度3 m。为提高计算准确性,在计算水域施加静水压力环境条件,在计算域外侧(除对称轴位置)设置环境边界条件。计算模型网格尺寸约为炸药半径的1/5,网格单元总数约为160万。距炸药中心1、2、3、4、5、6、7 m水平位置处设置压力观测点。图1为水中爆炸计算模型。

1.3 材料模型

非理想炸药采用高能炸药本构模型及带Miller后燃能量项的JWL状态方程描述,具体形式为

(1)

Miller后燃项计算式为

(2)

式(2)中:Q0为添加的后燃能量; λ为反应率; a为能量释放常数;m为能量释放指数;n为压力指数;p为压力。一般取m=1/2,n=1/6。a取值与炸药特性和非理想成分有关[15]。当Q=0或a=0时,JWL-Miller状态方程退化为JWL方程。

水材料采用空材料本构模型与Mie-Grüneisen状态方程描述。对于压缩状态

(3)

对于膨胀状态

p=ρ0C2μ+(γ0+αμ)E

(4)

式(4)中:C为曲线的截距;S1、S2、S3分别为曲线斜率的一阶、二阶和三阶系数;γ0为Gruneision常数;a为对γ0的一阶体积修正; ρ0为初始密度;E为初始内能。

空气材料采用空材料本构模型与线性多项式状态方程描述,其由内能线性关系式表示,压力计算如下:

P=C0+C1μ+C2μ2+C3μ3+

(C4+C5μ+C6μ2)Ε

(5)

式(5)中,如果μ<0,C2μ2和C6μ2设置为0,μ=ρ/ρ0-1,ρ/ρ0是当前密度比上的参考密度。

典型材料参数[14,16-17]见表1—表3。

2 结果与讨论

2.1 后燃效应对典型爆炸载荷特征参数影响

针对RS211炸药,研究非理想炸药JWL状态方程中的Miller后燃项对水中爆炸冲击波压力及气泡脉动载荷参数的影响。图2所示为有无Miller后燃项条件下典型观测点(距炸药中心水平距离3 m处观测点的)冲击波压力时程曲线对比。表4所示为炸药JWL状态方程中有无Miller后燃项条件下典型观测点水中爆炸载荷参数对比。

由图2可以看出:虽然冲击波在其峰值过后的计算值出现压力的多峰尾随振荡[18],但带Miller后燃项状态方程(Q=0.015 Mbar,a=0.006 5 Mbar6/μs,m=1/2,n=1/6)计算所得爆炸冲击波压力时程曲线始终处于上方,说明非理想炸药的冲击波衰减速度比理想炸药慢,冲击波冲量较高。后燃效应使得炸药水中爆炸冲击波压力峰值增加,气泡脉动周期增大。

由表4可知:增加Miller后燃项后,相同观测点的爆炸冲击波压力峰值增加,最大增幅约为4.7%,气泡周期增幅约为6.7%。分析认为:Miller后燃模型考虑了非理想炸药中铝粉的后燃反应对爆轰产物能量输送[19]。爆炸冲击波压力峰值主要是含铝炸药中理想成分的理想爆轰过程决定的。非理想成分在理想成分爆轰后发生后燃反应并释放能量,其对冲击波压力峰值略有贡献(峰值小幅增加)。炸药后燃效应可使气泡获得持续能量并不断膨胀,导致其脉动周期增加。若在数值计算中不考虑非理想炸药的后燃效应,则其水中爆炸气泡载荷相关的特征参数(如气泡周期、气泡半径等)可能会偏小,因此,在非理想炸药水中爆炸载荷模拟时,推荐选用JWL-Miller状态方程,以更准确的描述非理想炸药能量释放特性。

2.2 基于正交试验法的Miller后燃模型参数敏感性分析

由于JWL-Miller状态方程参数的试验测试较为复杂,工程仿真参数较为缺乏。不同的Miller模型参数取值代表了不同的后燃能量释放特性。明确判断各参数对水中爆炸载荷特征参量的影响程度及主次顺序,是更高效准确的进行Miller模型参数匹配设置的关键。

采用正交试验法设计(仿真)试验,研究Miller模型参数对典型水中爆炸载荷参数的影响。查阅相关文献[20-23]知:典型非理想含铝炸药后燃模型参数的取值(cm-g-μs单位制)范围为:Q0≤0.1,a≤0.065,m≤0.75,n≤0.254。考虑一般适用性,在正交试验中适当扩展部分参数的变化范围,表5所示为正交试验因素水平表。

设计4因素5水平共25工况的正交试验,试验评价指标为距爆心3 m处观测点冲击波压力峰值与气泡脉动周期。正交试验结果如表6所示。

通过正交试验数据计算每个因素及其水平所对应试验结果的均值及其极差,极差的数值直接表示了该因素水平对试验结果的敏感性,数值越大说明该因素对评价指标的敏感性越高。表7为冲击波压力峰值极差分析,表8为气泡脉动周期极差分析。图3为各因素对冲击波压力峰值均值影响趋势图。图4为各因素对气泡脉动周期均值影响趋势图。

分析表7与表8结果知:对于冲击波超压峰值与气泡脉动周期,影响因素的敏感程度依次为:Q0>a>n>m。由图3知:在正交试验条件范围内,冲击波压力峰值随Q0与a增大而增大,随m增大先增大后减小,随n增大而减小。由图4知:气泡脉动周期随Q0增大而增大。随a、m与n增大,气泡脉动周期未呈现明显变化规律。分析认为:在基于试验数据的Miller模型参数匹配时,可按照Q0、a、n与m的先后顺序以及各参数对目标值影响规律进行各参数匹配设置。

2.3 典型非理想炸药RS211后燃模型参数匹配

基于文献[14]中的试验工况条件,利用TNT炸药水中爆炸冲击波压力峰值与气泡脉动周期试验数据进行计算模型的边界尺寸、网格尺寸、粘性系数及时间步长等参数[24]的匹配设置。以计算结果与试验结果误差小于10%为达标判据,实现水中爆炸载荷计算模型及参数的标定。图5所示为典型时刻冲击波压力云图。图6所示为典型时刻气泡脉动形态。由图5与图6可以看出计算模型能够反映水域中爆炸冲击波与气泡演化过程。

表9所示为TNT炸药标定模型的计算与试验结果对比。由表9知:经模型标定后,TNT炸药爆炸载荷参数计算结果误差小于7.4%。综上,标定后的计算模型能够有效模拟炸药水中爆炸过程,准确预测冲击波压力峰值与气泡脉动周期等爆炸载荷参数。

基于上述标定后的计算模型,开展非理想含铝炸药RS211水中爆炸载荷计算。利用文献[14]中的试验数据及前述Miller后燃模型参数敏感性研究结论,以计算结果与试验结果误差小于10%为评价指标,匹配得到RS211炸药后燃模型参数为:Q0=0.01 Mbar,a=0.006 Mbar1/0.085/μs,m=0.5,n=0.085。表10所示为RS211炸药计算与试验结果对比。由表10可知:在文献[14]中的试验工况条件下,RS211炸药的爆炸冲击波峰值及气泡脉动周期的计算值与试验值误差均小于10%。故参数匹配后的JWL-Miller后燃模型可反映非理想炸药RS211水中爆炸能量释放特性,实现炸药水中爆炸载荷参数准确预测。

2.4 模型进一步验证与推广应用

前述研究基于文献[14]工况条件,构建了5 kg球形RS211炸药在5 m水深处爆炸载荷计算模型(含Miller模型参数)。为实现工程推广应用,验证上述计算模型在不同药量与不同水深工况条件下的适用性,利用文献[25]中的试验数据,开展3 kg球形RS211炸药在12 m水深处爆炸载荷计算。基于试验工况设置大气环境及水域边界条件。表11所示为典型工况计算结果与试验结果对比。

由表11可以看出:计算所得炸药爆炸冲击波峰值及气泡脉动周期与试验值相比,误差小于10%。表10与表11数据对比结果表明:研究所得非理想炸药水中爆炸载荷计算模型及Miller模型参数可适用于文献[14]与文献[25]中不同工况条件下RS211炸药水中爆炸载荷模拟。

3 结论

1) JWL-Miller状态方程能够描述非理想炸药水中爆炸能量释放特性,其后燃效应增大了冲击波压力峰值及气泡脉动周期。

2) 对于非理想炸药水中爆炸载荷参数(如冲击波压力峰值与气泡脉动周期)数值计算,炸药JWL-Miller状态方程的参数敏感性顺序为:Q0>a>n>m。在模型参数优化时可优先匹配后燃能量Q0及能量释放速率系数a。

3) RS211炸药Miller后燃模型参数可匹配设置为Q=0.01,a=0.006,m=0.5,n=0.85(cm-g-μs单位制)。典型工况条件下的冲击波压力峰值与气泡脉动周期的计算结果与参考试验结果误差均小于10%,后燃模型参数满足工程应用需求。研究成果可支撑相关水中兵器战斗部威力性能分析及其快速优化设计。

参考文献：

[1] 田俊宏,孙远翔,张之凡.铝氧比对含铝炸药水下爆炸载荷及能量输出结构的影响[J].高压物理学报,2019,33(6):148-156.

TIAN Junhong,SUN Yuanxiang,ZHANG Zhifan.Effect of Al/O ratio on underwater explosion load and energy output configuration of aluminized explosive[J].Chinese Journal of High Pressure Physics,2019,33(6):146-154.

[2] 伍俊,庄铁栓,杨程宇.非理想炸药水下爆炸能量释放简化模型及数值分析研究[J].防护工程,2017(3):29-32.

WU Jun,ZHUANG Tieshuan,YANG Chengyu.Research on the simplified model of energy release and numerical analysis for nonideal explosives underwater explosion[J].Protective Engineering,2017(3):29-32.

[3] 胡宏伟,严家佳,陈朗,等.铝粉含量和粒度对CL-20含铝炸药水中爆炸反应特性的影响[J].爆炸与冲击,2017,37(1):157-161.

HU Hongwei,YAN Jiajia,CHEN Lang,et al.Effect of aluminum powder content and its particle size on reaction characteristics for underwater explosion of CL-20-based explosives containing aluminum[J].Explosion and Shock Waves,2017,37(1):157-161.

[4] 鲁忠宝,黎勤,马军利,等.鱼雷战斗部装药特点与发展[J].水下无人系统学报,2018,26(1):10-15.

LU Zhongbao,LI Qin,MA Junli,et al.Research on charge in torpedo warhead[J].Journal of Unmanned Undersea Systems,2018,26(1):10-15.

[5] 黄江华,张光法.海萨尔—一种适合水中兵器装填的新型炸药[J].水雷战与舰船防护,2003(1):27-29.

HUANG Jianghua,ZHANG Guangfa.A new explosive suitable for underwater weapon loading[J].Mine Warfare &Ship Self-Defence,2003(1):27-29.

[6] 周正青,聂建新,覃剑锋,等.铝氧比对含铝炸药性能影响的数值模拟[J].爆炸与冲击,2015,35(4):514-520.

ZHOU Zhengqing,NIE Jianxin,QIN Jianfeng,et al.Numerical simulations on effects of Al/O ratio on performance of aluminized explosives[J].Explosion and Shock Waves,2015,35(4):514-520.

[7] 项大林,荣吉利,李 健,等.基于KHT程序的RDX基含铝炸药JW状态方程参数预测研究[J].北京理工大学学报,2013,33(3):239-243.

XIANG Dalin,RONG Jili,LI Jian,et al.JWL equation of state parameters prediction of RDX-based aluminized explosive based on KHT code[J].Transactions of Beijing Institute of Technology,2013,33(3):239-243.

[8] 辛春亮,秦健,刘科种等.基于LS-DYNA软件的水下爆炸数值模拟研究[J].弹箭与制导学报,2008,28(3):156-158.

XIN Chunliang,QIN Jian,LIU Kezhong,et al.Research on UNDEX numerical simulation based on LS-DYNA[J].Journal of Projectiles Rockets Missiles and Guidance,2008,28(3):156-158.

[9] JING Pinglu,DAVID L K.Modeling of PBXW-115 using kinetie CHEETAH and DYNA Codes[Z].1996.DSTO-TR-1496.

[10] 安丰江,施惠基,刘琼,等.水下爆炸近场压力特性及其与结构的流固耦合作用[J].兵工学报,2015,36(1):13-23.

AN Fengjiang,SHI Huiji,LIU Qiong,et al.Load characteristics of near-field of underwater explosion and its fluid-structure interaction[J].Acta Armamentarii,2015,36(1):13-23.

[11] 陈龙伟,张华,汪旭光.水中多物质爆炸场的三维数值模拟[J].兵工学报,2009,30(2):1-4.

CHEN Longwei,ZHANG Hua,WANG Xuguang.Three-dimensional numerical simulation of multi-material explosive field in water[J].Acta Armamentarii,2009(2):1-4.

[12] 汪春辉,王嘉安,王超,等.基于S-ALE方法的圆柱体垂直出水破冰研究[J].力学学报,2021(11):3110-3123.

WANG Chunhui,WANG Jiaan,WANG Chao,et al.Research on vertical movement of cylindrical structure out of water and breaking through ice layer based on S-ALE method[J].Chinese Journal of Theoretical and Applied Mechanics,2021,53(11):3110-3123.

[13] LSTC.LS-DYNA keyword user’s manual[M].Lawrence Livermore National Laboratory,2019.

[14] 秦健,艾东民,吴成,等.几种炸药水下爆炸能量损失特性分析研究[J].兵工学报,2015,36(2):38-42.

QIN Jian,AI Dongmin,WU Cheng,et al.Research on energy dissipation characteristics of underwater explosion of explosives[J].Acta Armamentarii,2015,36(2):38-42.

[15] LI Peng.Study on energy output characteristic of non-ideal explosives[D].Beijing:Beijing Institute of Technology,2006.

[16] 段超伟,宋浦,胡宏伟,等.水下爆炸气泡脉动周期的简便计算方法[J].高压物理学报,2022,36(1):1-6.

DUAN Chaowei,SONG Pu,HU Hongwei,et al.Simple method for the calculation of bubble pulsation period in underwater explosion[J].Chinese Journal of High Pressure Physics,2022,36(1):1-6.

[17] 宋浦,杨凯,梁安定,等.国内外TNT炸药的JWL状态方程及其能量释放差异分析[J].火炸药学报,2013,36(2):42-45.

SONG Pu,YANG Kai,LIANG Anding,et al Difference analysis on JWL-EOS and energy release of different TNT charge[J].Chinese Journal of Explosive and Propellants,2013,36(2):42-45.

[18] 李健,荣吉利,杨荣杰,等.水中爆炸冲击波传播与气泡脉动的实验及数值模拟[J].兵工学报,2008,29(12):1437-1443.

LI Jian,RONG Jili,YANG Rongjie,et al.Experiment and numerical simulation for the bubble impulse in underwater explosion[J].Transactions of Beijing Institute of Technology,2008,28(12):1035-1038.

[19] 杨坤,陈朗,伍俊英,等.计算网格与人工粘性系数对炸药水中爆炸数值模拟计算的影响分析[J].兵工学报,2014,35(2):237-243.

YANG Kun,CHEN Lang,WU Junying,et al.The effects of computing grid and artificial viscosity coefficient on underwater explosion numerical simulation[J].Acta Armamentarii,2014(2):237-243.

[20] LU Yong.The law and characterization of energy release in high aluminized explosives[D].Nanjing:Nanjing University of Science &Technology,2014.

[21] XIE Xinyuan.Research on the energy release law of aluminized explosive[D].Mianyang:China Academy of Engineering Physics,2018.

[22] 黄菊,王伯良,仲倩,等.温压炸药能量输出结构的初步研究[J].爆炸与冲击,2012,32(2):164-168.

HUANG Ju,WANG Boliang,ZHONG Qian,et al.A preliminary investigation on energy output structure of a thermobaric explosive[J].Explosion and ShockWaves,2012,32(2):164-168.

[23] 辛春亮,徐更光,刘科种,等.含铝炸药Miller能量释放模型的应用[J].含能材料,2008,16(4):436-440.

XIN Chunliang,XU Gengguang,LIU Kezhong,et al.Application of miller energy release model for aluminized explosive[J].Chinese Journal of Energetic Materials,2008,16(4):436-440.

[24] 荣吉利,项大林,李健,等.含铝炸药水下爆炸特性研究[J].北京理工大学学报,2012,32(3):221-225.

RONG Jili,XIANG Dalin,Li Jian,et al.Study on under water explosion character of aluminized explosive[J].Transactions of Beijing Institute of Technology,2012,32(3):221-225.

[25] 张远平,池家春,龚晏青,等.含铝炸药水下爆炸性能的实验研究[J].高压物理学报,2010,24(4):316-320.

ZHANG Yuanping,CHI Jiachun,GONG Yanqing,et al.Experimental study on underwater explosion performance of aluminized explosive[J].Chinese Journal of High Pressure Physics,2010,24(4):316-320.

Numerical simulation of underwater explosion load produced by nonideal explosives based on JWL-Miller equation of state

GU Hongping, CHEN Da, ZHANG Lijian, LYU Yongzhu, LI Baohua

(Xi’an Modern Chemistry Research Institute, Xi’ an 710065, China)

Abstract：To study on the underwater explosion load produced by nonideal explosives, the computing models of underwater explosion load produced by nonideal explosives with JWL-Miller stateequation are constructed based on S-ALE algorithm of LS-DYNA software, and the load characteristics of underwater explosion are analyzed. The influence sensitivity of the underwater explosion load produced by parameters of Miller afterburn model is studied based on orthogonal test method. Through the results comparison of simulation and experiment, the parameters of miller model are matched. The computing models are verified and applied, and the results show that the JWL-Miller state equation is appropriate for simulation of the underwater explosion load produced by nonideal explosives. The error between the numerical results and the experimental results of the shock wave pressure peak and bubble pulsation period produced by typical RS211 explosive in water is less than 10%. The methods, models and parameters obtained can provide support for the explosive power analysis and optimization design of relevant underwater weapon warhead.

Key words： explosion mechanics; underwater explosion load; nonideal explosives; numerical simulation; JWL-Miller

本文引用格式：谷鸿平,陈达,张立建,等.基于JWL-Miller状态方程非理想炸药水中爆炸载荷模拟[J].兵器装备工程学报,2023,44(10):293-299.

Citation format：GU Hongping, CHEN Da, ZHANG Lijian, et al.Numerical simulation of underwater explosion load produced by nonideal explosives based on JWL-Miller equation of state[J].Journal of Ordnance Equipment Engineering,2023,44(10):293-299.

中图分类号：TJ630

文献标识码：A

文章编号：2096-2304(2023)10-0293-07

收稿日期：2023-04-02;

修回日期：2023-05-20

作者简介：谷鸿平(1983—),男,硕士,研究员,E-mail:ghp911@sina.com。

doi：10.11809/bqzbgcxb2023.10.040

科学编辑 杨继森 博士(重庆理工大学教授)

责任编辑 唐定国