RDX基PBX的做功能力及JWL状态方程参数确定*

王新颖，王树山，徐豫新，胡 赛

(1.北京理工大学爆炸科学与技术国家重点实验室，北京 100081;2.沈阳理工大学装备工程学院,辽宁 沈阳 110159)

RDX基PBX的做功能力及JWL状态方程参数确定*

王新颖1,2，王树山1，徐豫新1，胡 赛1

(1.北京理工大学爆炸科学与技术国家重点实验室，北京 100081;2.沈阳理工大学装备工程学院,辽宁 沈阳 110159)

为研究RDX基PBX炸药的做功能力并确定其爆轰产物的JWL状态方程参数，对RDX基PBX炸药和TNT炸药进行∅50 mm标准圆筒实验，获得了圆筒膨胀位移和速度的时程曲线，对比得出RDX基PBX炸药的做功能力明显高于TNT炸药；基于能量守恒对实验数据进行非线性拟合，得到2种炸药爆轰产物的JWL状态方程参数。TNT炸药的拟合参数和通过AUTODYN软件计算得到的结果符合较好；将采用上述方法得到的RDX基PBX炸药爆轰产物JWL状态方程参数用于数值模拟，计算结果与实验值吻合较好，符合数值模拟标定JWL状态方程参数的要求。

爆炸力学；做功能力；圆筒实验；RDX基PBX；JWL状态方程

RDX基高聚物粘结炸药充分利用高能炸药的爆轰性能和高分子添加剂的力学性能，较传统的铸装炸药具有较高的能量密度和爆轰性能，在军事上广泛应用于导弹战斗部装药、推进剂、水雷、鱼雷、反坦克导弹等，工业上用于爆炸成型、石油射孔弹等[1-3]。炸药作为武器毁伤的起始能源，其做功能力是炸药的研制和战斗部装药设计以及威力评估的基础问题。炸药的状态方程是描述炸药爆轰做功的载体爆轰产物物理参数之间的关系。

目前已经有多种较为成熟的爆轰产物状态方程形式，如γ律状态方程、JWL、BKW、KHT及VLW状态方程等，其中JWL状态方程能够比较精确地描述爆轰产物的膨胀驱动过程，JWL状态方程是由E.L.Lee等[4]在前人工作的基础上提出的，该方程的未知参数可以通过J.W.Kury等[5]提出的圆筒实验来确定。圆筒实验在评估炸药做功能力中得到了广泛应用，文献[6-8]通过圆筒实验的数值模拟方法确定JWL状态方程参数。

本文中通过∅50 mm标准圆筒实验，以TNT炸药为基础，考察RDX基PBX炸药的做功能力，并基于实验数据和能量守恒，运用Origin软件进行非线性拟合得出JWL状态方程参数，通过数值模拟对比实验数据验证参数的准确性，以期为后续研究提供依据。

1 圆筒实验及结果分析

1.1 实验模型

实验中选用压装TNT药柱和RDX基PBX药柱，尺寸均为∅50 mm×495 mm。TNT药柱为压装，密度为1.58 g/cm3，爆速为6 836 m/s，RDX基PBX药柱为质量分数96%的RDX和4%的粘结剂压装，密度为1.67 g/cm3，爆速为8 279 m/s。采用的铜管为无氧铜，密度为8.9 g/cm3，屈服强度为200 MPa，拉伸强度为380 MPa，铜管内径为50.10 mm，外径为60.34 mm。实验装置如图1所示。狭缝位置距离起爆端295 mm，采用GSJ高速转镜相机记录圆筒壁在狭缝两端的膨胀过程，扫描速度是1.5 mm/μs，并通过固定在圆筒两端的电探针测定炸药的实际爆速。

图1 圆筒实验装置图Fig.1 Equipment diagram of the cylinder test

1.2 实验数据处理

实验获得的圆筒壁膨胀过程的扫描底片如图2所示，对底片进行判读，得到圆筒壁膨胀距离时程曲线，如图3所示。

图2 圆筒壁膨胀过程原始扫描底片Fig.2 Scanning films showing expansion process of cylindrical wall

图3 不同炸药驱动下圆筒壁膨胀距离时程曲线 Fig.3 Radial expansion displacement of cylindrical wall generated by different explosives

实验数据处理采用文献[9]中的处理方法，认为圆筒壁在爆炸冲击波和气体膨胀产物的共同作用下发生膨胀运动。冲击波使得圆筒壁在初期膨胀速度迅速增加，但持续时间较短；而气体爆轰产物膨胀时的推动力峰值虽然偏低，但持续时间较长。圆筒壁膨胀距离R-R0与膨胀时间t有如下关系：

(1)

式中：R和R0分别为圆筒在t时刻和初始时刻的半径，vs为冲击波驱动圆筒壁加速达到的渐进速度，τs为其加速时间段；vg为气体爆轰产物膨胀驱动圆筒壁加速达到的渐进速度；τ1和τ2分别为爆轰产物膨胀力波动时上升和下降的时间常数。将式(1)对时间求导，可得到圆筒壁的膨胀速度u的计算公式：

(2)

将实验获得的2种炸药的圆筒壁膨胀距离和时间数据按式(1)进行拟合，其拟合系数见表1。将表1中拟合系数代入式(2)中，可得到圆筒壁膨胀速度时程曲线，如图4所示。

表1 圆筒壁膨胀位移曲线拟合系数Table 1 Fitting parameters for radial expansion displacement curve of cylindrical wall

由图4可以看出：在圆筒膨胀初期，圆筒壁在冲击波力的作用下，RDX基PBX炸药的膨胀速度迅速增大到约840 m/s，TNT炸药的膨胀速度则增大到约630 m/s；随后圆筒壁在爆轰产物作用下，膨胀速度缓慢增大，直至膨胀破裂时，速度达到最大值即破片初速v0，广泛适用的计算破片初速的方法为G.W.Gurney等[10]提出的格尼公式：

(3)

表2 破片初速及炸药格尼系数Table 2 Initial velocity of fragment and Gurney coefficient of explosive

图4 不同炸药驱动下圆筒壁膨胀速度时程曲线 Fig.4 Radial expansion velocity of cylindrical wall generated by different explosives

2 JWL状态方程参数

2.1 实验数据拟合法

基于能量守恒，炸药的化学能全部转化为爆轰产物气体的内能，随着爆轰产物的膨胀做功逐渐转换为爆轰产物动能和金属圆筒的动能，故有

E0=Ei+Ek,g+Ek,m=Ei+Ek

(4)

式中:E0为炸药的质量化学能，Ei爆轰产物膨胀做功后自身剩余内能，Ek,g和Ek,m分别为爆轰产物动能和金属圆筒动能。

JWL状态方程不显含化学反应，能精确描述爆轰产物的等熵膨胀过程，其形式为:

(5)

式中：ps为爆轰产物压力，V为爆轰产物相对比容；A、B、C、R1、R2和ω为JWL状态方程的待定参数。

(6)

式中：Es(V)为爆轰产物等熵内能，是比容V的函数。炸药瞬时定容爆轰时，V=1，此时为驱动的初态：

(7)

则膨胀过程中爆轰气体剩余的内能为：

(8)

将式(6)和式(7)代入式(4)，可得驱动的有效能量Ek为：

(9)

根据格尼公式可有驱动金属圆筒的膨胀速度u为[11]：

(10)

根据式(9)对∅50 mm标准圆筒的TNT炸药和RDX基PBX炸药得到Ek(V)-V的关系，根据式(10)即可得到u-V关系曲线，如图5所示。

图5 不同炸药驱动下圆筒壁u-V曲线和Ek(V)-V曲线 Fig.5 u-V and Ek(V)-V curves of cylindrical wall generated by different explosives

根据式(9)应用Origin软件非线性拟合[12]即可得到2种炸药的JWL状态方程的参数，见表3，TNT炸药的拟合参数和AUTODYN软件计算的结果相比，符合较好。

表3 爆轰产物的JWL状态方程参数Table 3 Parameters in JWL state equation of detonation product

图6 数值模拟的几何模型 Fig.6 Geometry model of simulation

2.2 圆筒实验数值模拟

为验证RDX基PBX 炸药爆轰产物的JWL状态方程的参数，选用AUTODYN软件对RDX基PBX圆筒实验进行数值模拟，圆筒实验为二维轴对称模型，RDX基PBX密度为1.67 g/cm3，爆速为8 279 m/s，采用高能炸药燃烧模型和JWL状态方程，JWL状态方程参数见表3；无氧铜采用各向同性弹塑性流体动力学模型[13]。几何模型如图6所示，为避免端部效应的影响，观测点选取在圆筒壁上距爆点200～300 mm处，并等间隔设置。

图7所示为数值模拟结果与实验值的对比，可以看出，基于上述能量守恒方法拟合出的JWL状态方程参数用于数值模拟，计算结果与实验值吻合较好，符合数值模拟标定JWL状态方程参数的要求。

图7 RDX基PBX炸药圆筒实验的数值模拟与实验结果比较 Fig.7 Cylinder test results of RDX-based PBX explosive compared with numerical simulation

3 结 论

通过∅50 mm圆筒实验对RDX基PBX炸药的做功能力进行评价，并基于能量守恒，拟合JWL状态方程参数，得到以下结论：

(1)RDX基PBX爆轰驱动的破片初速比TNT炸药高24.3%，格尼系数高21.1%，其做功能力明显强于TNT炸药。

(2)基于能量守恒原理，通过非线性拟合得到2种炸药爆轰产物的JWL状态方程参数，TNT炸药的拟合参数与已知结果较符合；采用上述方法得到的RDX基PBX炸药爆轰产物JWL状态方程参数用于数值模拟，计算结果与实验值吻合较好，符合数值模拟标定JWL状态方程参数的要求。

[1] 孙国祥.高分子混合炸药[M].北京:国防工业出版社,1985:3-4.

[2] 孙华,郭志军.PBX炸药技术特性及在水中兵器上的应用[J].装备指挥技术学院学报,2009,20(3):108-111. Sun Hua, Guo Zhijun. Characteristics of PBX dynamite and its application in undersea weaponry[J]. Journal of the Academy of Equipment Command & Technology, 2009,20(3):108-111.

[3] 孙业兵,曹新茂.军用混合炸药[M].北京:兵器工业出版社,1995:10-11.

[4] Lee E L, Hornig H C, Kury J W. Adiabatic expansion of high explosive detonation products: UCRL-50422[R]. San Francisco: University of California, 1968.

[5] Kury J W, Hornig H C, Lee E L. Metal acceleration by chemical explosives[C]∥Proceedings of 4th Symposium on Detonation. White Oak, Maryland, 1965:3-13.

[6] 陈清畴,蒋小华,李敏,等.RDX基高聚物黏结炸药JWL状态方程[J].含能材料,2011,19(2):213-216. Chen Qingchou, Jiang Xiaohua, Li Min, et al. JWL equation of state for RDX based PBX[J]. Chinese Journal of Energetic Materials, 2011,19(2):213-216.

[7] 谭凯元,韩勇,罗观,等.HMX基PBX的作功能力及其JWL状态方程[J].火炸药学报,2013,36(3):42-45. Tan Kaiyuan, Han Yong, Luo Guan, et al. Power ability and JWL equation of state of a HMX-based PBX[J]. Chinese Journal of Explosives & Propellant, 2013,36(3):42-45.

[8] 沈飞,王辉,袁建飞.一种确定JWL 状态方程参数的简易算法[J].振动与冲击,2013,33(9):107-110. Shen Fei, Wang Hui, Yuan Jianfei. A simple method for determining parameters of JWL EOS[J]. Journal of Vibration and Shock, 2013,33(9):107-110.

[9] Lindsay C M, Butler G C, Rumchik C G. Increasing the utility of the copper cylinder expansion test[J]. Propellants, Explosives, Pyrotechnics, 2010,35(5):433-439.

[10] Gurney G W. The initial velocities of fragments from bombs, shells and grenades[R]. Army Ballistic Research Laboratory, 1943.

[11] 王新颖,王树山,徐豫新,等.爆轰驱动金属圆筒的能量转换与破片初速模型[J].兵工学报,2015,36(8):1417-1422. Wang Xinying, Wang Shushan, Xu Yuxin, et al. The energy conversion and fragment initial velocity model of metal cylinder driven by detonation[J]. Acta Armamentarii, 2015,36(8):1417-1422.

[12] 方安平.Origin 8.0实用指南[M].北京:机械工业出版社,2009.

[13] 韩勇,黄辉,黄毅民,等.含铝炸药圆筒试验与数值模拟[J].火炸药学报,2009,32(4):14-17. Han Yong, Huang Hui, Huang Yimin, et al. Cylinder test of aluminized explosives and its numerical simulation[J]. Chinese Journal of Explosives & Propellants, 2009,32(4):14-17.

(责任编辑 王易难)

Power capability and parameters of JWL equation of state for RDX-based PBX

Wang Xinying1,2, Wang Shushan1, Xu Yuxin1, Hu Sai1

(1.StateKeyLaboratoryofExplosionScienceandTechnology,BeijingInstituteofTechnology,Beijing100081,China;2.SchoolofEquipmentEngineering,ShenyangLigongUniversity,Shenyang110159,Liaoning,China)

In order to study the power capability of RDX-based PBX and determine the parameters of the JWL equation of state, the ∅50 mm standard cylinder test of RDX-based PBX and TNT were carried out, the expansion distance-time curve and velocity-time curve of the cylinder wall were obtained. Compared with the TNT, PBX had obviously a higher power capability. Based on the theory of energy conservation and the nonlinear fitting of the experimental data, the parameters of the JWL equation of state were obtained. Compared with the known software parameters and by numerical simulation of the cylinder test by AUTODYN, a fairly good agreement was reached between the experimental value and the actual test results, proving that this method for obtaining the JWL parameters is viable.

mechanics of explosion; power capability; cylinder test; RDX-based PBX; JWL equation of state

10.11883/1001-1455(2016)02-0242-06

2014-11-18；

2015-03-10

王新颖(1980— )，女，博士研究生，讲师，wxy801003@163.com。

O381 国标学科代码： 13035

A