Studies on the state equation of the underwater detonation products for PBX-01 explosive
-
摘要: 提出通过水中实验确定炸药的水中爆轰产物JWL状态方程参数的方法；选择PBX-01高能炸药进行水中实验，利用ANSYS/LS-DYNA程序建立炸药的水中实验模型，将实验结果与数值计算结果进行对比，确定PBX-01炸药水中爆轰产物的JWL状态方程参数。研究结果显示，圆筒实验确定的JWL参数在反映炸药水中爆轰产物的膨胀状态时有所不足，水中实验确定的JWL状态方程参数能够更准确地描述PBX-01炸药水中爆轰产物的膨胀过程，因此对水中爆炸的研究需要通过水中爆炸实验建一套状态方程参数。Abstract: A method was raised by taking out the underwater explosion test to determin the JWL state equation parameters of detonation products. The PBX-01 explosive was applied and a calculating model of underwater explosion test was established by using ANSYS/LS-DYNA program. The JWL state equation parameters of detonation products for PBX-01 explosive were determined according to the comparison of the experimental and numerical results. Parameters determined by underwater explosion test and cylinder expansion test were used for detonation products driving water test. The results show that the underwater explosion test numerical simulation is accordant better in initial stage. Specially appointed patameters need to be established in the underwater explosions.
-
Key words:
- mechanics of explosion /
- equation of state /
- underwater explosion /
- detonation products
-
在常规武器系统中，聚能装药战斗部扮演着越来越重要的角色，提高炸药能量利用率是战斗部设计者追求的目标之一。带隔板装药不仅结构简单而且可以在不增加装药量以及不改变装药长径比的条件下，提高其对药型罩的压垮做功能力。典型的带隔板装药结构如图 1所示，根据惠更斯原理[1]，爆轰波经过隔板的绕射后将以O′点开始传播，以较小的入射角驱动药型罩形成聚能侵彻体。这类带隔板装药在驱动药型罩形成高速杆式射流、大长径比爆炸成型弹丸以及利用起爆控制技术形成多模毁伤元等方面应用非常广泛。
图 1 带隔板柱形装药结构及爆轰波示意图Figure 1. Schematic diagram for cylindrical charge with waveshaper and detonation wave早在20世纪70年代，汤明钧等[2]利用高速摄影测定了带隔板破甲弹中装药端面上爆轰波的爆速和爆压，得出马赫波的波速达到正常爆速的1.9倍，爆压达到正常CJ爆压的3倍。杨全中等[3]对带隔板装药爆轰波的相互作用进行了研究，并利用Whitham方法[4]计算了爆轰波传播过程中的波形和碰撞临界角等参数。B.B.Dunne等[5]指出了发生马赫反射的临界入射角取决于爆轰产物的多方指数，并给出了马赫杆压力等参数的计算方法。L.M.Hull等[6]计算了球面爆轰波的马赫反射参数，通过实验给出了PBX9501和PBXN-9等炸药的马赫反射参数。张俊秀等[7]进行了爆轰波马赫反射的实验研究，并根据修改的Whitham方法[4]计算了马赫反射的压力、爆速和增长角等参数。Zhang Xianfeng等[8]、张先锋等[9]分别利用理论计算和数值模拟分析了夹层装药爆轰波的传播过程，并基于质量守恒给出了爆轰波传播过程中马赫杆高度的计算模型。
本文中，基于三波理论[3]和Whitham方法[4]，对带隔板装药爆轰波相互作用的正规反射和马赫反射进行理论分析，并给出马赫杆高度的计算模型。在此基础上，利用AUTODYN有限元计算程序中的Lee-Tarver模型[10]对带隔板装药爆轰波传播过程进行数值模拟, 以期研究成果可为描述带隔板装药驱动典型介质提供理论基础和数据支持。
1. 带隔板装药爆轰波马赫反射理论模型
1.1 带隔板装药爆轰波相互作用
2个等强度爆轰波的正规斜碰撞，可视为其中一侧爆轰波对刚性壁面的正规斜碰撞，当入射角φⅠ较小时，流入界面I和界面R的粒子流的偏转角相等。入射角不断增大，达到马赫反射的临界入射角φⅠ, c时，偏转角θ2小于θ1。粒子流的不断积累迫使碰撞点P移动到Q点，便形成了图 2所示的马赫反射。基于F.Grasso等[11]关于爆轰波反射的描述，爆轰前反应区用带下标1的参数表示，反射波前反应区用带下标2的参数表示，超压爆轰反应区用带下标3的参数表示。
假设炸药完全反应且满足CJ条件，区域(1)中的爆轰参数可由炸药的CJ参数、初始密度ρ0、多方指数γ和入射角获得。粒子流通过区域(1)发生偏转流入区域(2)，偏转角和反射角满足下列关系:
tanθ1=[Ma21sin2(φⅡ+ε)−1]tan−1(φⅡ+ε)Ma21[(γ+1)/2−sin2(φⅡ+ε)]+1 (1) 式中Ma1=√1+(γ+1)2/(γtanφⅠ)2为马赫数; 发生正规反射时，ε=θ1；发生马赫反射时，ε=θ1-θ2。
发生正规反射和非正规反射时, 式(1)中ε的表达式是不同的，若用发生正规反射时ε=θ1求解发生非正规反射时式子(1)的表达式，在0＜φⅡ+ε＜π/2范围内必定无解。根据此条件可以获得发生马赫反射的临界入射角，不同炸药临界入射角的理论计算值和实验结果如表 1所示, 表中DCJ为炸药CJ爆轰速度。通过比较发现，两者比较接近，从而可以验证马赫反射临界入射角计算的正确性。通过计算还发现，发生马赫反射时的临界入射角与炸药多方指数成一一对应关系，随着多方指数的增大，发生马赫反射的临界入射角不断减小，如图 3所示，此结论与B.B.Dunne[5]给出的结果一致。
表 1 不同炸药装药中发生马赫反射时临界入射角的计算值与实验结果Table 1. Calculated and experimental results of critical incident angles corresponding to Mach reflection in different explosive charges
图 3 马赫反射临界入射角随多方指数的变化Figure 3. Critical incident angle for the onset of Mach reflection varying with polytropic exponent大量理论和实验研究发现[5, 11-12]，马赫杆是一个弯曲的曲面，交于对称轴P点，相切于爆轰波(I)Q点。在P点, 偏转角α=θ1-θ2=0°; 在Q点, α=θ1-θ2；马赫杆与对称轴切线的夹角从π/2变化到φⅠ。由于碰撞点P处的压力和密度等都远大于CJ值，难以准确获取，本文中借助炸药爆轰产物的状态方程p=p(E, v)来求得此处的压力, 其中p为压力，E为炸药能量，v为比容。由H.M.Sternberg等[12]的定义，产生马赫反射时爆轰产物中将从三波点产生一条滑动线。此直线通过的地方，压力是连续的，但速度和密度非连续，因此，在p3=p2处存在一个三波点。
1.2 马赫杆高度模型
基于Whitham方法[4, 13]获取马赫反射的增长角和高度等参数，并将马赫波简化为平面波。将爆轰波的每一元素沿每个微元射线管的传播当作固壁管道中的激波传播来处理，在二维情况下，爆轰波的位置和与其垂直的射线位置如图 4所示。爆轰波的传播过程可以看成是在变截面流管中的传播过程，当截面发生微小变化时，流动符合等熵流性质，得到马赫反射的A(z)关系式为：
图 4 马赫反射后爆轰波形位置与其垂直射线的关系简图Figure 4. Relationship of detonation wave and its vertical line after Mach reflectionAf(z)=A0f(1)=const (2) 式中:z=p/pCJ为爆轰波强度，且：
f(z)=z1/γ(2z−1)1/2B(z)C(z);C(z)=exp{γ−1/2sin−1[(3−γ)z−2(γ+1)z]} B(z)={[2(γ−1)(2z−1)]1/2+2[(γ−1)z+1]1/2}[2γ/(γ−1)]1/2 基于上述讨论，如图 4所示, 假设Δt时刻后，由爆轰波NS变成马赫波PQ，φⅠ为爆轰波入射角，χ为三波点增长角，N和Q点为不同时刻的三波点。设NS管道截面积为ACJ，PQ处管道截面积为AM，根据几何关系知:
tanφⅠ=(DCJDM+AMACJ){[1−(DCJDM)2][1−(AMACJ)2]}−1/2 (3) cosχsin(φⅠ−χ)=DMDCJ (4) sinχcos(φⅠ−χ)=AMACJ (5) 由式(5)知:
AMACJ=f(1)f(z) (6) 参考B.D.Lambourn等[14]的研究，将爆速DM和DCJ定义为:
DMDCJ=z(2z−1)1/2 (7) 上述计算中未考虑爆轰波后稀疏波的影响，计算结果误差较大，对Whitham方法[4, 13]中的A(z)关系式进行修正:
ACJAM=(2z−1)12k(z) (8) 式中:k(z)是慢变函数，一般取常数, 通常由实验确定。
根据张俊秀等[7]提供的入射角和爆速参数计算方法，通过方程(2)~(8)，得到马赫波参数DM、z和χ的计算结果如表 2所示。利用Whitham方法[4, 13]计算获得的马赫波速度和三波点增长角与实验结果[7]吻合较好，从而验证此方法是可行的。
φⅠ/(°) DCJ/(km·s-1) DM, exp/(km·s-1) DM, cal/(km·s-1) DM,exp−DM,calDM,exp/% z χexp/(°) χcal/(°) χexp−χcalχexp/% 46.5 4.14 5.89 5.86 0.5 3.42 1.5 1.55 -3.2 48.5 4.20 5.77 5.79 -0.3 3.22 2.27 2.13 6.6 49.5 4.48 6.07 6.11 -0.7 3.13 2.80 2.45 14.3 同一装药中炸药的CJ爆速为常数，三波点增长角只与入射角和k(z)有关。根据带隔板装药爆轰波传播的特点，将入射角离散化，网格点(n, n)上的值可由网格点(n-1, n-1)上的值求得，马赫杆的高度可以定义为:
{Hn=Rtanχn−1tanφn−1−Ln−1tanφn−1+Hn−11+tanχn−1tanφn−1Ln=Ln−1tanχn−1tanφn−1+(R−Hn−1)tanφn−11+tanχn−1tanφn−1 (9) 式中:R为隔板半径；Hn和Ln分别为马赫杆的高度和传播距离，n=1, 2, 3, …。
当n=1时，H0=1, L1=Rtanφ0。L.M.Hull[5]也认为，当入射角为临界入射角时，马赫杆高度为零或非常小。结合表 1的数据和方程(3)~(9)，可以确定马赫杆高度。
爆轰波绕射过隔板后，爆轰波不断向轴线方向汇聚，入射角不断增大。当入射角小于临界入射角时，发生正规反射，马赫杆高度为零。当入射角大于临界入射角时，发生马赫反射，马赫杆高度随着入射角的增大而增大，相同炸药中，随着k(z)的变化，马赫杆高度计算结果如图 5所示。马赫杆在不断增大的过程中，马赫杆上的压力和速度迅速降低，当马赫杆上的压力和速度降低到CJ值时，马赫杆消失，此时，马赫反射变成正常的CJ爆轰。
图 5 马赫杆高度随入射角变化的计算结果Figure 5. Calculated variation of Mach stem height with incident angle2. 带隔板装药爆轰波传播及马赫反射过程数值模拟
2.1 有限元计算方案
根据前述对带隔板装药爆轰波传播过程的分析，在保证隔板外侧爆轰波可以正常传播的条件下，隔板直径越大，调整爆轰波形状的能力越强，在装药结构中爆轰波马赫杆的形状就越明显。带隔板装药爆轰波的传播过程是一种多物质相互作用的大变形运动过程，因此本文中采用多物质耦合的Euler算法描述带隔板装药爆轰波的传播过程。图 6所示为装药的有限元模型。
图 6 带隔板装药数值模拟的有限元模型Figure 6. A finite element model for numerical simulation on charge with waveshaper2.2 材料模型及参数
凝聚炸药爆轰产物状态方程描述炸药爆轰达到CJ状态后爆轰产物中压力、温度和体积等物理量之间的关系，不仅体现炸药本身的做功能力，而且是研究爆炸力学各种问题的基础。本文中分别采用点火与增长模型的Lee-Tarver状态方程[10]和JWL状态方程[15]来描述爆轰波绕射过隔板后在炸药中的传播过程。炸药参数主要来自文献[10, 15]及AUTODYN中内置的参数。
为了获取环形爆轰波形并保证透射波不起爆主装药，隔板厚度应足够大，这里隔板厚度选20 mm。隔板直径选42 mm，药柱直径选50 mm，使得隔板外侧装药厚度大于爆轰波正常传播时装药的临界直径。隔板材料为环氧树脂，其状态方程为Puff状态方程[16]，主要参数分别为：ρ0=1.985 g/cm3，A1=17.89 MPa，A2=201.2 MPa，A3=-6.053 MPa，Grüneisen系数为1.72，膨胀系数为0.25。
2.3 数值模拟结果与分析
针对上述讨论，对主装药分别采用JWL[15]和Lee-Travel[10]2种不同的状态方程进行描述，爆轰波波形的模拟结果如图 7所示。基于JWL模型[15]所形成的爆轰波的波阵面存在多段线组合形式，如图 7(a)所示。产生该结果的主要原因是：由于JWL模型为程序起爆模式，而炸药微元的爆轰采用预先设定时间起爆，JWL模型在描述带隔板装药爆轰(特别是隔板直径较大)时存在局限性，较难精确描述炸药微元的超压爆轰行为。基于炸药冲击起爆点火与增长的Lee-Travel模型[10]形成的爆轰波形较光滑，如图 7(b)所示。通过对比发现，Lee-Travel模型[10]能更真实地描述带隔板装药爆轰波的传播过程。
图 7 采用不同的炸药状态方程模拟得到的爆轰波波形Figure 7. Detonation waveform of explosive simulated by different equations of state对带隔板装药中的主副装药分别采用JWL[15]和Lee-Travel[10]2种不同的状态方程进行描述，爆轰波传播过程中不同时刻的波形如图 8所示。装药起爆后，爆轰波开始以球面波向前传播，经过隔板绕射后变成汇聚爆轰波。汇聚爆轰波在装药轴线处发生相互碰撞反射，当t=86 μs时，爆轰波已经发生马赫反射，随着爆轰波的传播, 马赫杆的高度不断增加，如图 8(b)和(c)所示。随着爆轰波的传播，马赫杆的高度不断增加，但其压力和速度迅速降低，当马赫杆上的压力和爆轰速度等参数降低到CJ值时，马赫反射消失，马赫反射变成正常爆轰, 如图 8(d)所示。
为了进一步了解带隔板装药爆轰波传播过程中发生马赫发射后参数的变化，在装药轴线以及距轴线2.5、5.0、7.5和10.0 mm处分别设有一系列高斯点。以TNT炸药为例，通过设置高斯点可以获取爆轰产物的参数变化，当马赫反射发生后，随着爆轰波的传播，不同位置处压力和密度的变化见图 9~10。发生马赫反射时，装药轴线处(即马赫杆碰撞点处)的压力和密度均达到最大值，其中最大压力与CJ压力的比值达到4.3，最大密度与CJ密度的比为2.05，与Zhang Xianfeng等[8]对TNT装药发生马赫反射时最大压力和密度的理论计算值接近。随着爆轰波的传播，马赫杆碰撞点的压力和密度逐渐降低。通过对距轴线不同距离处设置高斯点发现，压力和密度先增大后减小，这说明马赫杆存在一定高度且是不断变化的。
图 9 马赫反射发生后，随着爆轰波的传播，在轴线不同位置处压力的变化Figure 9. Pressure at different positions of the axis after Mach reflection
图 10 马赫反射发生后，随着爆轰波的传播，在轴线不同位置处密度的变化Figure 10. Density at different positions of the axis after Mach reflection为了获取带隔板装药中马赫杆高度的变化，分别对带隔板装药中的Comp B、TNT和PBX9404炸药爆轰波的传播过程进行了数值模拟。随着爆轰波的传播，马赫杆高度的精确捕捉越来越难，其数值模拟结果和理论计算结果如图 11~13所示。
当入射角较小时，马赫杆的高度为零，随着入射角不断增大并达到临界入射角时，马赫反射发生。此时，马赫杆碰撞点的压力和密度均达到最大值，随爆轰波的传播，马赫杆高度不断增加，马赫杆上的压力和密度不断减小。当马赫杆上的压力和密度等降到CJ值时，马赫反射消失，变成正常爆轰。
3. 结论
通过对带隔板装药爆轰波马赫反射过程的理论分析以及数值模拟，得到的主要结论如下：
(1) 基于三波理论和修正的Whitham方法，得到了马赫反射临界入射角和三波点增长角等参数，最终获得了马赫杆高度的计算模型。计算结果表明：爆轰波相互作用发生马赫反射时，马赫杆的高度为零；随着爆轰波的传播，马赫杆的高度不断增加，马赫杆上的压力和速度降低，当其降低到CJ值时，马赫反射消失。
(2) 基于点火增长模型的Lee-Travel模型，带隔板装药爆轰波发生马赫反射后马赫杆高度的数值模拟结果与理论计算结果吻合较好，这说明利用Lee-Travel模型可以较准确地描述带隔板装药爆轰波的传播过程。
(3) 带隔板装药技术是一种利用爆轰波相互作用产生超压爆轰的装药技术。马赫杆高度模型的建立不仅丰富和完善了带隔板装药爆轰波马赫反射的传播规律，而且可为描述带隔板装药驱动药型罩的过程提供理论基础。
-
图 2 PBX-01炸药水中爆轰产物膨胀的高速扫描底片
Figure 2. Photograph of the underwater detonation products expansion process
图 5 PBX-01炸药爆轰产物膨胀位移曲线与水中实验值
Figure 5. Calculated and tested expansion distance-time curves in underwater explosion test of PBX-01 explosive
表 1 PBX-01炸药爆轰产物JWL状态方程参数
Table 1. The JWL state equation parameters of the PBX-01 detonation products
实验 A/GPa B/GPa R1 R2 ω E0 圆筒实验 406.4 16.3 3.90 1.45 0.5 11.48 水中实验 356.5 26.3 3.40 1.14 0.5 12.48 
下载: 导出CSV
-
[1] 穆哈迈德·苏塞斯卡.炸药测试方法[M].邵英斌, 译.绵阳: 中国工程物理研究院流体物理研究所, 1998: 151. [2] 安丰江, 吴成, 王宁飞.水下爆炸能量耗散特性分析研究[J].北京理工大学学报, 2011, 31(4): 379-382. http://www.cqvip.com/QK/95548X/20114/37717930.htmlAn Feng-jiang, Wu Cheng, Wang Ning-fei. A research on the energy dissipation of underwater explosion[J]. Transactions of Beijing Institute of Techbology, 2011, 31(4): 379-382. http://www.cqvip.com/QK/95548X/20114/37717930.html [3] 李建, 荣吉利, 雷旺.水下爆炸气泡运动的理论研究[J].应用力学学报, 2010, 27(1): 120-124. http://d.wanfangdata.com.cn/Periodical/yylxxb201001022Li Jian, Rong Ji-li, Lei Wang. Research on the bubble motion in underwater explosion[J]. Chinese Journal of Applied Mechanics, 2010, 27(1): 120-124. http://d.wanfangdata.com.cn/Periodical/yylxxb201001022 [4] 牟金磊, 朱锡, 李海涛, 等.炸药水下爆炸能量输出特性实验研究[J].高压物理学报, 2010, 24(2): 87-92.Mu Jin-lei, Zhu Xi, Li Hai-tao, et al. Experimental research on underwater explosion energy output of explosive[J]. Chinese Journal of High Pressure Physics, 2010, 24(2): 87-92. [5] 饶国宁, 陈网桦, 胡毅亭, 等.不同炸药水下能量输出特性的实验研究[J].爆破器材, 2007, 36(1): 9-11. http://d.wanfangdata.com.cn/Periodical/bpqc200701003Rao Guo-ning, Chen Wang-hua, Hu Yi-ting, et al. Experimental study on underwater energy output characteristics of different explosives[J]. Explosive Materials, 2007, 36(1): 9-11. http://d.wanfangdata.com.cn/Periodical/bpqc200701003 [6] 蒋国岩, 金辉, 李兵, 等.水下爆炸研究现状及发展方向展望[J].科技导报, 2009, 27(9): 87-91. http://www.cnki.com.cn/Article/CJFDTotal-KJDB200909033.htmJiang Guo-yan, Jin Hui, Li Bing, et al. Review of studies on underwater explosion[J]. Journal of Science and Technology, 2009, 27(9): 87-91. http://www.cnki.com.cn/Article/CJFDTotal-KJDB200909033.htm [7] Kury J W, Horing H C, Lee E L, et al. Metal acceleration by chemical explosives[C]//Fourth Symposium(International)on Detonation. Washington DC, USA. 1965: 3-13. [8] 董海山, 周芬芬.高能炸药及相关物性能[M].北京: 科学出版社, 1989. [9] 龙新平, 韩勇, 蒋志海, 等.炸药驱动水的初期过程[J].爆炸与冲击, 2010, 30(1): 12-16. doi: 10.11883/1001-1455(2010)01-0012-05Long Xin-ping, Han Yong, Jiang Zhi-hai, et al. Measurement and numerical simulation of initial stage about detonation products driving water[J]. Explosion and Shock Waves, 2010, 30(1): 12-16. doi: 10.11883/1001-1455(2010)01-0012-05 [10] 韩勇, 黄辉, 黄毅民, 等.含铝炸药圆筒实验与数值模拟[J].火炸药学报, 2009, 32(4): 14-17. http://www.cqvip.com/QK/90400B/20094/31907828.html [11] 龙新平. VLW爆轰产物状态及纳米级含铝炸药爆炸特性研究[D].北京: 北京理工大学, 1999. [12] Lee E, Breithaupt D, Mcmillan C, et al. The motion of thin metal walls and the equation of state of detonation products[R]. Albuquerque, NA, USA: Lawrence Livemore National Laboratory, 1985. [13] Walsh J M, Rice M H. Dynamic compression of liquid from measurements on strong shock waves[J]. Journal of Chemistry Physics, 1957, 26(4): 815-823. doi: 10.1063/1.1743414 [14] 彭福泉.实用非金属材料手册[M].吉林: 吉林科学技术出版社, 1991: 389. 期刊类型引用(6)
1. 李昱锦，缪玉松，李斌，黄飞飞，葛帅，杨涛，吴佳怡. 对称双线性起爆下爆轰波的轴向传播特性. 爆破器材. 2025(01): 7-13 . 
百度学术2. 张朝平，张先锋，谈梦婷，侯先苇，熊玮，刘闯，顾晓辉. 聚能杆式射流侵彻混凝土和岩石靶体试验与数值模拟. 含能材料. 2023(08): 773-785 . 百度学术3. 李瑞，李伟兵，靳洪忠，王桂林，洪晓文，朱建军，李文彬，王晓鸣. 基于Jones-Wilkins-Lee状态方程的爆轰波相互作用参数理论分析. 兵工学报. 2019(03): 516-521 . 百度学术4. 潘建，张建军，侯云辉，冯雪磊. 带隔板装药EFP成型数值模拟与实验研究. 振动与冲击. 2019(18): 89-94+121 . 百度学术5. 潘建，张建军，侯云辉，刘燕宁. 带隔板装药爆轰波驱动药型罩数值模拟研究. 弹箭与制导学报. 2019(04): 149-154 . 百度学术6. 成乐乐，赵太勇，陈智刚，印立魁，王维占，郭光全，雷文星. 分层装药对杆式射流的影响分析. 爆破器材. 2018(04): 35-39 . 百度学术其他类型引用(5)
-
下载:
百度学术
计量
- 文章访问数: 3151
- HTML全文浏览量: 351
- PDF下载量: 596
- 被引次数: 11