一种基于Hugoniot关系的爆轰产物等熵状态方程*

刘益儒,胡晓棉

(北京应用物理与计算数学研究所，北京 100088)

炸药爆轰产物等熵状态方程描述的是热力学平衡条件下爆轰产物在C-J点后近似等熵膨胀过程中的热力学状态参量之间的关系，是炸药做功能力的主要表征。研究爆轰产物等熵状态方程对掌握炸药爆轰性能及发展爆炸力学具有重要意义。

目前爆轰产物等熵状态方程分为两类：一类是从理论模型出发的状态方程，如BKW、JCZ和Williamsburg等熵状态方程等[1-3]，这类方程形式复杂且计算量较大，不适用于工程计算；另一类为经验或半经验状态方程，最常用的有常γ以及JWL形式的等熵状态方程[4-5]，这类方程形式简单、比较适用于工程计算，但往往需要特定的标定实验确定方程经验参数，例如常γ状态方程的经验参数通过拟合平面爆轰波冲击靶板的初始自由表面速度实验数据确定，JWL状态方程的经验参数通过拟合圆筒实验数据确定。这些标定实验往往成本高、周期长，且在获得实验数据后，还需考虑如何合理拟合这些实验数据，并且目前获得较好的标定结果还依赖于复杂的二维流体动力学数值计算，这又增加计算成本。

为此，有必要发展一种不需要实验标定，形式简单，又能满足工程计算的爆轰产物等熵状态方程。本文中基于压力-粒子速度平面上爆轰产物Hugoniot经验关系，利用Riemann积分，发展一种新的爆轰产物等熵状态方程形式，该方程的参数仅为炸药的初始比容和C-J状态量，这些参数通过简单且成本较低的炸药密度和爆轰参数测量实验即可直接获得，避免复杂的参数标定过程，和传统的经验或半经验等熵状态方程相比，可节约标定参数计算成本。

1 爆轰产物等熵状态方程建模

1.1 爆轰产物Hugoniot经验关系

Cooper[6]对136组炸药爆轰产物压力粒子速度实验数据进行收集并建立数据库，这些实验数据涉及14种理想或非理想炸药，初始密度范围为1.13～7.47 g/cm3。采用的实验方法包括自由表面速度法和阻抗匹配法：两者均采用平面爆轰波冲击不同厚度的的惰性靶板，但前者测量靶板内自由表面粒子速度，将数据外推到零厚度靶板，再利用靶板冲击Hugoniot关系求得爆轰产物压力粒子速度数据；后者(如水箱法)测量距界面不同位置处的冲击波速度或压力，并外推到界面位置，再利用靶板冲击Hugoniot关系求得p-u数据。Cooper[6]采用相应炸药的C-J爆压pCJ和粒子速度uCJ对这些数据进行无量纲化后，发现，p/pCJ和u/uCJ在高压和低压时分别较好地符合二次方函数和幂函数关系，因此拟合得到如下分段形式的爆轰产物爆轰产物压力粒子速度Hugoniot经验关系：

(1)

值得注意的是，表征函数分段位置d值的选取与人为判断有关，并没有明确物理意义，除了约束条件，d值的选取亦可直接影响拟合曲线的拟合优度。文献[6]中的拟合曲线的拟合优度值为R2|F1=0.957和R2|F2=0.806，而文献[7]中的拟合曲线的拟合优度值为R2|F1=0.947和R2|F2=0.789。可见后者的拟合优度不如前者。因此，对函数分段位置做适当调整，即令d=0.06，对实验数据进行重新分段拟合，得到新的拟合系数为a=2.264，b=-1.528，c=0.264，m=112.5，n=-8.186，从而拟合优度值达到R2|F1=0.966 9和R2|F2=0.901 5。图1为线性坐标和对数坐标下无量纲参数间关系曲线。

1.2 爆轰产物等熵状态方程

利用式(1)推出爆轰产物等熵状态方程的过程如下：

(2)

(3)

根据等熵线上的关系：

(4)

式中:下标S代表等熵状态量，c为等熵声速，ρ为产物密度。可以推出如下Riemann积分：

(5)

式中:v= 1/ρ为产物比容。将式(3)代入式(5)，并考虑约束条件a+b+c=1和b+2c=1,整理得：

(6)

式(6)即为提出的爆轰产物等熵状态方程的具体形式。式中：v0是炸药初始比容，V=v/v0是相对比容，c、m和n均为1.1节中给出的爆轰产物在p-u平面上的Hugoniot经验关系拟合系数，是与具体炸药性能无关的常数，因此不需要针对特定炸药进行实验标定。另外可以看到，方程参数仅为炸药的初始比容v0和爆轰产物C-J点的状态量即pCJ、uCJ和vCJ，其中v0和pCJ可分别通过简单的密度测量实验和爆压测量实验获得，uCJ和vCJ则可以利用下式进行估算：

(7)

式中：D为爆速，通过简单的爆速测量实验即可获得。这些实验与传统的爆轰产物状态方程标定实验相比，成本显著降低；另外，由于不需要利用实验数据拟合经验参数，从而节约了计算成本。

将本文中得到的拟合系数，即c=0.264，m=112.5，n=-8.186代入式(6)，进一步化简得：

(8)

为了在模型验证中进行对比分析，这里将Lanterman[7]的拟合系数，即c=0.258 5，m=206.5，n=-8.62，代入式(6)，得到等熵状态方程：

(9)

2 模型验证

为验证模型的合理性，分别采用式(8)、(9)给出的等熵状态方程以及JWL等熵状态方程在p-V平面上绘制了Comp-B、HMX、PETN、ANFO、TNT以及LX-14这6种爆轰性能有明显差异的炸药爆轰产物等熵膨胀曲线，如图2所示。这里采用JWL等熵状态方程进行对比，是因为该方程能够比较精确的描述爆轰产物的等熵膨胀过程。炸药组分质量比、初始密度、C-J爆速和爆压等参数可参阅文献[8] 。

描述p-V关系的JWL等熵状态方程的形式[5]为：

pS(V)=Aexp(-R1V)+Bexp(-R2V)+CV-(ω+1)

(10)

式中:A、B、C、R1、R2和ω是经验参数，需要通过拟合圆筒实验数据确定；这6种炸药的JWL状态方程经验参数值列于表1[8]。

从图2中可以看到，不同于JWL状态方程描述的连续光滑膨胀曲线，式(8)～(9)给出的状态方程所描述的膨胀曲线分别在p/pCJ=0.06和p/pCJ=0.1(或V≈2和V≈1.5)位置处出现分段，这很有可能是在拟合分段函数形式的产物Hugoniot经验关系时，未考虑分段位置处的连续性造成的。另外，与式(9)给出状态方程曲线相比，式(8)给出状态方程曲线的连续性以及与JWL等熵状态方程曲线的符合程度均更好，这说明在拟合Hugoniot关系曲线时，将分段位置从d=0.1调整为d=0.06，不仅提高了Hugoniot关系曲线的拟合优度，还使在此基础上推出的等熵状态方程更符合实际情况，因此应优先选用式(8)给出的等熵状态方程。

表1 炸药爆轰产物等熵JWL状态方程参数[8]Table 1 JWL parameter for explosives[8]

现比较分析式(8)给出的状态方程与JWL状态方程。可以看到，从总体上看，这2个状态方程描述的等熵膨胀曲线符合较好，特别是在高压阶段(p/pCJ>0.06)，两者基本重合。在低压阶段(p/pCJ≤0.06)特别是接近函数分段位置处存在一定差异，可能的原因是：(1)在拟合确定低压阶段爆轰产物Hugoniot经验关系时，采用的实验数据样本不足；(2)从在产物Hugoniot经验关系基础上推出的新状态方程形式上看，高压段时，方程为指数函数形式(见式(8))，这与JWL状态方程高压阶段起主要作用的方程形式类似(见式(10)等号右边第1项)，而在低压段时，方程为幂函数形式(见式(8))，而JWL状态方程在这个压力范围内(对应于JWL方程的中压段)起主要作用的方程形式为指数函数和幂函数之和的形式(见式(10)等号右边第2和第3项)，依此可以判断，采用的低压段Hugoniot经验关系的函数形式可能不是很合理，需要进一步改进。

3 结 论

本文中基于拟合得到的爆轰产物Hugoniot经验关系，提出一种新的爆轰产物等熵状态方程，该方程不包含需要实验标定的经验参数，方程参数仅为炸药的初始比容v0和爆轰产物C-J点状态量，这些参数的测量实验与传统的经验爆轰产物等熵状态方程参数的标定实验相比，成本显著降低；另外，由于不需要拟合方程经验参数，极大地节约了计算成本。

采用该状态方程给出了Comp-B、HMX、PETN、ANFO、TNT以及LX-14炸药的爆轰产物等熵膨胀曲线，与JWL状态方程描述的等熵膨胀曲线进行比较，发现两类曲线符合较好，尤其是在高压阶段，两者高度重合。

[1] SUCESKA M, ANG H G, CHAN S, et al. Study of the effect of covolumes in BKW equation of state on detonation properties of CHNO explosives[J]. Propellants, Explosives, Pyrotechnics, 2013,38(1):103-112.

[2] HOBBS M L, BRUNDAGE A, YARRINGTON C. JCZS2i: An improved JCZ database for EOS calculations at high temperatures and pressures[R]. SAND2014-19288PE 540872. Sandia National Laboratories, Albuquerque, New mexico, United States, 2014.

[3] BRAITHWAITE M, SHARPE G J. Reduced, chemistry implicit, equations of state for explosion and detonation products[C]∥Proceedings of 15th International Detonation Symposium. San Francisco, California, United States, 2014:869-877.

[4] 孙承纬,赵锋,文尚刚,等.一般物态方程形式下爆轰产物的一维等熵流动[J].爆炸与冲击,2003,23(6):481-487.

SUN Chengwei, ZHAO Feng, WEN Shanggang, et al. One dimensional isentropic flow of detonation products with general equation of state[J]. Explosion and Shock Waves, 2003,23(6):481-487.

[5] WANG H, GAO J, LUO Y, et al. Metal accelerating ability and JWL EOS of DNTF/HMX explosive[J]. Ordnance Industry Automation, 2014(7):64-66.

[6] COOPER P W. Shock behavior of explosives about the C-J (Chapman-Juget) point[C]∥Proceedings of 9th International Symposium on Detonation. Portland, Oregon, United States, 1989:379-387.

[7] LANTERMAN D D. An improved empirical fit of the detonation product isentrope near the CJ point[C]∥Proceedings of 15th International Detonation Symposium. San Francisco, California, United States, 2014:466-470.

[8] HORNBERG H. Determination of fume state parameters from expansion measurements of metal tubes[J]. Propellants, Explosives, Pyrotechnics, 1986,11(1):23-31.