许兴春,高欣宝,张俊坤

(军械工程学院,河北 石家庄 050003)



膨胀石墨燃爆剂JWL状态方程参数拟合

许兴春,高欣宝,张俊坤

(军械工程学院,河北 石家庄 050003)

用LS-DYNA对膨胀石墨燃爆剂爆炸过程进行数值模拟时，需要提供膨胀石墨燃爆剂的JWL状态方程的参数。把凝聚炸药等熵线物态方程作为目标方程，通过差分进化法拟合得出膨胀石墨燃爆剂JWL状态方程的参数。通过圆筒实验对所求得的参数进行验证，结果表明，理论值与实验值最大误差不超过3.3%，能够满足膨胀石墨燃爆剂爆炸模拟研究的需要。

爆炸力学;JWL状态方程;LS-DYNA;膨胀石墨

可膨胀石墨受热膨胀后成为中空状粒子，形似蠕虫，漂浮在空中能够有效干扰毫米波；不同规格的可膨胀石墨膨胀后可以得到1～8 mm甚至更长的粒子，因此有望使用一种材料遮蔽干扰不同的波长，它在发烟剂的应用中具有很大的潜力。实验表明，爆炸法能够快速并有效地在指定空域形成膨胀石墨型气溶胶云团，但爆炸实验成本高、风险大[1-3]。

若能计算出膨胀石墨燃爆剂的JWL[4-5]状态方程参数，就可以用LS-DYNA数值模拟研究代替部分实验研究，从而达到降低实验成本、提高研究效率的目的。一般来说，JWL状态方程参数需要通过圆筒实验及二维流体动力学程序确定。本文中提出一种JWL状态方程的参数近似解法，即通过凝聚炸药等熵线的物态方程推导出目标方程，通过差分进化法拟合得出膨胀石墨燃爆剂JWL状态方程的参数。

1 JWL状态方程及应用

(1)

式中：p为爆轰产物的压力,V为爆轰产物的相对比容,E为爆轰产物的比内能。A、B、R1、R2、ω为待拟合参数，也称经验(调节)常数。

在LS-DYNA中，对于炸药材料的定义需要输入JWL状态方程的A、B、R1、R2、ω等参数值。

2 膨胀石墨燃爆剂JWL状态方程关键参数的拟合

要拟合得到膨胀石墨燃爆剂JWL状态方程的参数，需要找到能够描述膨胀石墨燃爆剂爆炸时p-V关系的方程作为目标方程。那么对于凝聚炸药，忽略初始压力p0和初始内能E0，凝聚炸药的质量守恒、动量守恒、能量守恒与C-J条件(相切条件)方程，即凝聚炸药爆轰参数方程[8]为：

(2)

式中:u0、uH分别表示原始爆炸物的质点速度和爆轰波反应末端介质的质点速度；p0、T0、ρ0、v0分别表示原始爆炸物的压力、温度、密度和比容；pH、TH、ρH、vH分别表示爆轰波反应末端断面处的压力、温度、密度和比容；e0、eH分别表示原始爆炸物和爆轰波反应末端的内能；Qv表示爆轰反应释放出的化学能。

凝聚炸药爆轰产物状态方程为：

(3)

(4)

方程(4)为凝聚炸药爆轰产物的近似状态方程，由于其中没有温度项，该方程可近似为等熵方程，本文中称为凝聚炸药等熵线物态方程。其中，k是等熵常数，在CJ点的k值一般为3左右，本文中取为3。

一段还原设备为单管马弗炉，配带入口氢气加湿系统，入口氢气露点是-20 ℃以下的干燥氢气，在通过温度传感控制水温的热水水浴后实现加湿，加湿系统最高可将露点增加到略高于+20 ℃；二段还原设备是十八管还原炉。

那么方程组(2)可以简化成如下形式：

(5)

方程组(5)即凝聚态炸药爆轰波参数近似计算方程。

(6)

(7)

(8)

那么只要已知炸药的密度ρ0、爆速D，就可以根据方程(5)和方程(7)求解得到目标方程：

(9)

实验时采用的膨胀石墨燃爆剂配方为黑火药和可膨胀石墨的混合物，黑火药/可膨胀石墨质量之比为3∶2，其装药密度ρ=1.2 g/cm3；根据经验公式[9]可计算出膨胀石墨燃爆剂的爆速D≈850 m/s。当k=3时，根据公式(5)和(8)，计算得到pH=216.8 MPa，M=0.091 5，那么拟合目标函数方程可表示为：

(10)

利用LstOpt或MATLAB数据处理软件，采用差分进化法和遗传算法相结合的方法，对方程(1)参数进行拟合。参数拟合结果如表1所示，表中Emax表示最大误差。与其对应的拟合曲线如图1所示。

从表1中可以看出,Value(d)组数据最大误差为2.30%，精度比较高，选取其作为最终拟合结果。当然这5个参数值不是确定不变的，在满足一定精度条件下，只要选取其中精度较高的一组参数即可。

3 数值模拟与分析

圆筒为轴对称结构，可在柱坐标系中建立轴对称计算模型，如图2所示。图2中OCDE区为圆筒内混合炸药部分，ABCD区为铜管部分，在E点直接起爆。OE边的长度为圆筒长度，ED和EA的长度为铜管的内径和外径。JWL参数采用Value(d)组数据带入计算。

图3为圆筒初分网格图，图4为45 μs时圆筒网格变形图，此刻圆筒上端部分已经充分膨胀，炸药爆轰波继续向下传播。

利用后处理软件LS-PREPOST，得出圆筒A点的Δr-t(半径变化-时间)曲线如图5所示。

图3 网格划分图Fig.3 Diagram of meshing

图4 圆筒45 μs时网格变形图Fig.4 Gridding distortion at 45 μs

图5 数值模拟得到的Δr -t曲线Fig.5 Simulation curve of Δr -t

4 圆筒实验验证

圆筒实验[10-11]是专门用于确定炸药爆轰产物JWL状态方程参数和评定炸药做功的标准化实验，其实验原理图如图6所示，圆筒平行放置于支架上，高速扫描相机通过金属板狭缝记录燃爆剂稳定爆轰段圆筒膨胀距离。

圆筒实验数据一般被拟合成如下形式：

(11)

式中：t为圆筒壁膨胀的时间；(R-R0)为圆筒壁膨胀的距离，用Δr表示；a、b、c、d为根据实验数据得到的拟合系数；vc为不同膨胀距离(R-R0)相对应的圆筒壁的速度；Ekc为不同膨胀距离相对应的比动能。

实验时，圆筒半径R0=25 mm，膨胀石墨燃爆剂装药密度ρ0=1.2g/cm3，爆速D≈850m/s。膨胀石墨燃爆剂圆筒实验结果如表2所示。

实验所得圆筒Δr-t曲线如图7所示。

分析图5、7中数据，可以得出模拟曲线与实验曲线符合非常好，经过MATLAB对2组数据进行分析后，得出其误差最大为3.3%。

图6 圆筒实验原理示意图Fig.6 Diagram of cylinder test

图7 实验得到的Δr -t曲线Fig.7 Test curve of Δr -t

t/μs(R-R0)/mmt/μs(R-R0)/mm250．0654．0300．0704．5350．5755．0401．0805．5452．0856．0502．5906．6553．0957．2603．51008．5

5 结 束 语

把凝聚炸药等熵线物态方程作为目标方程拟合出膨胀石墨燃爆剂的JWL状态方程的关键参数，然后利用拟合的JWL参数对圆筒实验模型进行数值模拟得到Δr-t曲线，最后用圆筒实验得出圆筒的Δr-t曲线。通过MATLAB对两组数据进行分析，得出其误差最大为3.3%。实验结果表明：基于凝聚炸药等熵线物态方程拟合膨胀石墨燃爆剂JWL状态方程参数的方法可行，满足实际应用需求。

[1] 乔小晶,张同来,任慧,等.爆炸法制备膨胀石墨及其干扰性能[J].火炸药学报,2003,26(1):70-73. Qiao Xiao-jing, Zhang Tong-lai, Ren Hui, et al. Preparation by explosion method and the interference function of expanded graphite[J]. Chinese Journal of Explosives & Propellants, 2003,26(1):70-73.

[2] 潘功配,关华,朱晨光,等.可膨胀石墨用作抗红外/毫米波双模发烟剂的研究[J].含能材料,2007,15(1):70-72. Pan Gong-pei, Guan Hua, Zhu Chen-guang, et al. Study on anti-IR/millimeter wave smoke composition with expansive graphite[J]. Chinese Journal of Energetic Materials, 2007,15(1):70-72.

[3] 潘功配,杨硕.烟火学[M].北京：北京理工大学出版社，2004.

[4] 黄辉,王泽山,蒋小华,等.挤柱炸药爆轰产物JWL状态方程研究[J].弹道学报,2007,19(1):45-47. Huang Hui, Wang Ze-shan, Jiang Xiao-hua, et al. Cylinder experiments and JWL EOS of detonation product for ECX[J]. Journal of Ballistics, 2007,19(1):45-47.

[5] 于川,李良忠,黄毅民.含铝炸药爆轰产物JWL状态方程研究[J].爆炸与冲击,1999,19(3):274-278. Yu Chuan, Li Liang-zhong, Huang Yi-min. Studies on JWL equation of state of detonation product for aluminized explosive[J]. Explosion and Shock Waves, 1999,19(3):274-278.

[6] Kury J W, Hornig H C, Lee E L, et a1. Metal accelaration by chemical explosives[C]∥Proceedings of the 4th International Symposium on Detonation. White Oak, Maryland, 1966:3-13.

[7] Kury J W. Metal acceleration by chemical explosives[C]∥4th Syrup (International) on Detonation. White Oak: Office of Naval Research, 1965:3-12． [8] 奥尔连科.爆炸物理学[M].孙承纬,译.北京:科学出版社,2011.

[9] 陈智刚,赵太勇,侯秀成.爆炸及其终点效应[M].北京：兵器工业出版社,2004.

[10] 孙占峰,李庆忠,孙学林,等.标准圆筒试验技术与数据处理方法研究[J].高压物理学报,2008,22(2):160-166. Sun Zhan-feng, Li Qing-zhong, Sun Xue-lin, et al. Study on standard cylinder test technology and data processing method[J]. Chinese Journal of High Pressure Physics, 2008,22(2):160-166.

[11] 陈朗,李泽仁.激光速度干涉仪测量法在炸药圆筒试验中的应用[J].爆炸与冲击,2001,31(3):229-232. Chen Lang, Li Ze-ren. The cylinder tests measured by VISAR interferometer[J]. Explosion and Shock Waves, 2001,31(3):229-232.

(责任编辑 曾月蓉)

Parameters fitting for the JWL EOS of expanded graphite bums agent

Xu Xing-chun, Gao Xin-bao, Zhang Jun-kun

(OrdnanceEngineeringCollege,Shijiazhuang050003,Hebei,China)

The simulation of the explosion of expanded graphite bums agent requires the determination of the parameters in the JWL EOS of expanded graphite bums agent. The JWL EOS of expanded graphite bums agent are fitted through the differential evolution method taking condensed explosives isentropic line EOS as the target equation. The cylinder test results show that the maximum error between the theoretical value and experimental value is less than 3.3%, meeting the needs of the expanded graphite bums agent explosion simulation research.

mechanics of explosion; JWL EOS; LS-DYNA; expanded graphite

10.11883/1001-1455(2015)01-0124-06

2013-06-07；

2013-10-25

许兴春(1986— )，男，博士研究生,doctxu@163.com。

O389;TJ5 国标学科代码： 13035

A