Effects of Metal Interlayer and Air Gap on the Shock Initiation of Insensitive Explosives
-
摘要: 为了确定空气间隙和金属隔层对冲击起爆的影响,采用火炮加载蓝宝石飞片冲击起爆
∅ 50 mm × 30 mm的A型炸药,产生的冲击波通过空气间隙和金属隔层起爆∅ 50 mm的台阶型B型炸药。在B型炸药的后界面粘贴镀膜氟化锂(LiF)窗口,使用光子多普勒测速仪(PDV)测量金属和B型炸药的后界面速度,进而计算得到金属和B型炸药的冲击波透射压力,再利用阻抗匹配计算得到金属和B型炸药的入射压力。结果表明:传爆药和金属隔层间的空气间隙使冲击压缩过程转变为准等熵压缩和冲击压缩两个过程,同时使冲击波的幅值减小;确定了金属隔层厚度为5 mm时冲击波压力的衰减范围;当使用A型炸药作为传爆药,空气间隙为0.3 mm,金属隔层厚度为5 mm时,B型炸药在7~10 mm之间开始反应。-
关键词:
- 金属隔层 /
- 空气间隙 /
- 冲击起爆 /
- 光子多普勒测速仪(PDV) /
- 台阶型炸药
Abstract: The influences of air gap and metal interlayer on the shock initiation of the∅ 50 mm stepped explosive B have been investigated by means of the sapphire flyer planar impact experiment and photonic doppler velocimetry (PDV) technique. In the experiment, a lithium fluoride (LiF) window was stuck to the rear interface of explosive sample to allow the rear interface velocity between the metal and the sample explosive be measured by PDV technique. Both the transmission shock pressure and the incident shock pressure can be obtained by the impedance-match method. The experiment results have shown that the air gap divided the impact compression process into quasi-isentropic compression and impact compression, respectively, and at the same time, the amplitude of the shock pressure was declined. Specifically, the shock wave attenuation range caused by 5 mm thick metal was achieved. When the explosive A was used as the booster, and with a 0.3 mm thick air gap and 5 mm thick metal experimental setup, the reaction of sample explosive B started in a range from 7 mm to 10 mm.-
Key words:
- metal interlayer /
- air gap /
- shock initiation /
- photonic doppler velocimetry /
- stepped explosive
-
钝感高能炸药(Insensitive high explosive, IHE)具有突出的安全性能[1],是目前最具影响力的高能材料,在国防工业和工程爆破中有着广泛应用,其冲击起爆特征受到广泛关注。深入了解炸药冲击起爆过程,对研究可靠起爆炸药、避免炸药在意外冲击下发生爆炸具有重要意义。
隔层对冲击波具有明显的衰减效果,隔板起爆实验是研究炸药冲击起爆特性的有效方法之一[2],其中金属材料隔层广泛应用于冲击起爆实验。陈朗等[3]利用隔板实验对固体炸药的冲击起爆过程进行研究,得到了JO-9159炸药的起爆压力阈值和爆轰成长距离。胡湘渝等[4]运用DYNA2D程序对隔板实验过程进行了数值模拟,并与实验结果进行对比,较好地拟合了实验结果。郑波等[5]利用隔板实验对比探究了验证板对固体推进剂的冲击起爆影响。向梅等[6]通过数值模拟比较了铝、有机玻璃、钢隔层和未反应JB-9014炸药的爆轰反应边界值,发现有机玻璃对冲击波的衰减和铝近似,钢隔层对冲击波的衰减较快。以上研究均把隔层作为一种实验工具,没有针对隔层对冲击起爆的具体影响开展研究。
在冲击起爆过程中,空气间隙是一个重要的干扰因素。空气间隙与隔层有较大差异,传爆药和隔层间的空气间隙不仅会使冲击波衰减,而且会使本来的冲击压缩过程转变成准等熵压缩[7]和冲击压缩两个过程。因此,探究空气间隙和金属隔层对冲击起爆的影响,对于研究炸药的冲击起爆特性具有重要意义。
本研究以炸药驱动金属隔层和金属隔层冲击起爆B型炸药两个实验为基础,在实验中加入了空气间隙,通过对测得的金属隔层后界面速度和B型炸药的粒子速度进行分析,得到金属隔层和空气间隙对冲击起爆的影响。
1. 实验原理和方法
在火炮加载平台上对金属隔层和B型炸药进行冲击实验,采用光子多普勒测速仪[8-10](Photonic Doppler velocimetry, PDV)测量金属隔层和B型炸药样品的粒子速度。样品和PDV光纤探头粘贴于样品架,位置关系如图1所示。
实验分为3个过程,装置如图2所示。首先,测量金属隔层自由界面速度,如图2(a)所示,用飞片撞击尺寸为
∅ 50 mm × 30 mm的 A型炸药,使其达到爆轰状态,爆轰波通过厚度为0.3 mm的空气间隙到达金属隔层,用PDV测量金属隔层自由界面速度,金属隔层自由界面速度的变化可以反映传爆药和金属隔层间的空气间隙对冲击波的影响,通过界面速度计算出透射压力,再根据入射压力得到金属隔层的冲击波透射率。其次,如图2(b)所示,在金属隔层后贴合如图3所示的台阶型B型炸药,台阶厚度分别为2、3、4、5、7和10 mm,每个台阶面粘贴尺寸为∅ 15 mm × 11 mm 的LiF窗口,通过PDV测速系统测量每个台阶样品和LiF窗口的界面速度,根据速度的变化,判断炸药的反应情况,得到传爆药和金属隔层间的空气间隙对冲击起爆的影响。最后,在金属隔层和B型炸药样品间增加一个0.22 mm厚的空气间隙,与上一过程中PDV测速系统测量得到的界面速度进行对比,得出金属隔层和B型炸药样品间的空气间隙对冲击起爆的影响。实验条件如表1所示。表 1 实验条件Table 1. Experimental conditionsShot No. Dimensions/(mm × mm) Sample Booster Air gap Metal compartment Air gap 01 ∅50 × 30 ∅50 × 30 ∅50 × 5 02 ∅50 × 30 ∅50 × 30 ∅50 × 5 03 ∅50 × 30 ∅50 × 30 ∅50 × 5 ∅50 × 0.22 Explosive B 04 ∅50 × 30 ∅50 × 30 ∅50 × 5 Explosive B 2. 结果与分析
2.1 金属隔层和空气间隙对透射压力的影响
通过解析解预先估算金属隔层的速度,根据声学近似理论[11],用等熵线方程描述爆轰产物,金属视为刚体,得到金属自由界面速度u
u=Dj(1+θ−1ηθ−lθDjt) (1) 式中:Dj为A型炸药的爆轰速度;θ和η为过程量,
θ=[1+2η(1−l/Djt)]−1/2 ,η=m(k+1)k−1/ Mkk ;l 为装药长度;m为A型炸药的质量;M为金属片的质量;k为爆轰产物的多方气体指数;t为炸药起爆的零点时间。本研究中取Dj = 8.87 km/s,k = 3,l = 31.09 mm,m = 113.525 g,M = 179.47 g。把Dj、k、m、M和l的值代入式(1),得到u关于t的函数,根据该函数得出的金属隔层速度曲线解析解和由Shot 01、Shot 02测量得到的金属隔层自由界面中心点速度曲线如图4所示。
由图4可知,测量得到的金属自由界面速度与理论计算得到的速度解析解相比有一定偏差,这是因为理论计算需要满足3个前提条件:(1)金属隔层为刚体;(2)爆轰波冲击绝热刚壁,反射冲击波的熵增极小,冲击波可视为波阵面上有间断的简单压缩波,爆轰产物的熵视为常数;(3)气体产物和被抛射物体都处于管内运动。上述条件与实际情况是不符合的。首先,金属并不是刚体,反应过程中有形变,其熵是减小的;其次,炸药的爆轰波存在尖峰,不是简单的压缩波;最后,本实验并非处于管内运动,实验空间是开放的,炸药质量应为有效部分质量ma,ma小于炸药质量m。因此,解析解并不能反映爆炸的实际反应情况,只能用来估算最终的理论速度,且实际速度相较于理论值偏小。金属自由界面速度没有持续加速,只有一次起跳,是因为实验中使用的金属材料阻抗较大,且相对较厚,由后界面反射的冲击波在到达前界面的过程中,在金属隔层中衰减,冲击波到达自由界面时,比自由界面压力低,不能使自由界面速度加速。在爆轰产物的推动作用下,金属隔层保持匀速飞行。
PDV测速系统测得的是金属后界面速度,根据应力波理论[12],冲击波后的粒子速度为后界面速度的1/2。由冲击波后粒子速度可以得到后界面的透射压力p2out
p2out=ρ2Dt12u2=12(c2+12λ2u2)u2 (2) 式中:ρ2为金属隔层密度,Dt为金属隔层冲击波速度,c2和λ2为线性Hugoniot关系式中的常系数,u2为金属后界面速度。
为了得到金属隔层的冲击波透射率,需要首先得到金属隔层的入射压强。可由A型炸药与金属隔层阻抗匹配[13]得出
pj=12ut(ρ2Dt+ρ1Dj)=12ut[ρ2(c2+λ2ut)+ρ1Dj] (3) 式中:pj为A型炸药的爆轰压力,ρ1为A型炸药密度,ut为金属隔层前界面粒子速度。
pj数值见表2,其中,uj为炸药爆轰时的粒子速度,cj为炸药爆轰时的内部声速。把pj代入式(3),求得金属隔层前界面粒子速度ut。金属隔层的入射压力p2in可描述为
表 2 A型炸药和JWL状态方程参数[14]Table 2. Parameters of explosive A and JWL equation of stateA/GPa B/GPa R1 R2 ¯ω ρ1/(g·cm–3) Dj/(km·s–1) uj/(km·s–1) cj/(km·s–1) pj/GPa 934.770 12.723 4.6 1.1 0.37 1.863 8.87 2.22 6.65 36.8 p2in=ρ2Dtut=ρ2(c2+λ2ut)ut (4) 根据图4,u2取图中速度起跳的最高点,把u2代入式(2),计算得出透射压力p2out,由式(4)计算得出入射压力p2in,结果列于表3。
表 3 第1次和第2次金属自由界面速度测试计算结果Table 3. Computation results of the first and second test for free surface velocities of metalShot No. u2/(km·s–1) D2/(km·s–1) p2in/ GPa p2out/ GPa Shot No. u2/(km·s–1) D2/(km·s–1) p2in/ GPa p2out/ GPa 01 1.061 4.703 68.828 45.675 02 0.982 4.654 68.828 41.817 1.113 4.734 68.828 48.193 1.058 4.700 68.828 45.490 0.935 4.625 68.828 39.554 1.125 4.742 68.828 48.804 1.118 4.738 68.828 48.475 1.182 4.777 68.828 51.644 1.133 4.746 68.828 49.181 1.095 4.723 68.828 47.304 图5为图4中第3个点的速度开始增长时前面部分的放大图像。由图5可知,金属自由界面粒子速度刚开始增长的时候,并不是直接跳跃至最高点,而是经过一段平滑的曲线增长。这是由于在金属隔层和传爆药之间存在厚度为0.3 mm的空气间隙,当冲击波通过空气间隙时,对空气进行压缩,金属隔层压力的增加是连续的,可以把这个过程看作是对金属隔层进行准等熵压缩。
同时,空气间隙也会使冲击波的幅值衰减,从冲击波传入空气间隙到传出这段时间,爆轰产物在空气间隙中膨胀,将其近似地视为等熵膨胀,在此期间压强的改变量
Δp 为Δp=∂p∂¯vΔ¯v (5) 式中:
Δ¯v 为空气间隙的厚度和装药厚度的比值。p由A型炸药的爆轰产物JWL状态方程表示p=A(1−¯ωR1¯v)e−R1¯v+B(1−¯ωR2¯v)e−R2¯v+¯ω¯ve0 (6) 式中:A、B、R1、R2、
¯ω 为经验系数,¯v 为爆轰产物的相对比容,e0为单位初始体积的比内能。对压强求相对体积的偏导可得∂p∂¯v=−AR1(1−¯ωR1¯v)e−R1¯v−BR2(1−¯ωR2¯v)e−R2¯v+A¯ωR1¯v2e−R1¯v+B¯ωR2¯v2e−R2¯v−¯ω¯v2e0 (7) A型炸药的JWL状态方程参数如表2所示。
把
Δ¯v 和JWL状态方程参数代入式(5)和式(7),得到爆轰产物在空气间隙中做等熵膨胀时减少的压强Δp :Δp = –119.819 1 × 0.013 38 GPa = –1.598 2 GPa。将炸药和金属阻抗匹配后的入射压力近似地处理为p2in与Δp 之和,得到金属隔层厚度为5 mm时各个测量点的透射率η1,结果如表4所示。表 4 金属自由界面速度测试Table 4. Test of free surface velocities of metalShot No. p2in/GPa p2out/GPa η1/% 01 67.230 45.675 67.94 67.230 48.193 71.68 67.230 39.554 58.83 67.230 48.475 72.10 67.230 49.181 73.15 67.230 47.304 70.04 02 67.230 41.817 62.20 67.230 45.490 67.66 67.230 48.804 72.59 67.230 51.644 76.82 由表4可知,金属隔层厚度为5 mm时,冲击波透射率范围为58.83%~76.82%,平均透射率为70.46%。考虑到侧向稀疏波影响,选择两次测试中心点的透射率,分别为72.10%和72.59%,取其平均值72.35%。
2.2 金属隔层对钝感炸药冲击起爆的影响
分别测得金属隔层与台阶型样品炸药之间存在和不存在空气间隙情况下样品的后界面速度,如图6和图7所示。
当金属隔层和B型炸药紧密贴合时,对样品炸药产生作用的是A型炸药爆轰后通过金属隔层传播的冲击波,图像显示为三角波。当金属隔层和B型炸药之间存在空气间隙时,金属隔层撞击B型炸药会产生一个持续的压力,输入的应是方波,由于本实验中空气间隙过小,仅为0.22 mm,因此未能形成方波。图7为不同厚度台阶面测试位置的速度对比图像。由图7可知,B型炸药厚度为2、3、5和10 mm时,当金属隔层和B型炸药之间有空气间隙时,第一次起跳的速度高于没有空气间隙,这是由于金属隔层在通过空气间隙时产生了加速。
由图6可知,厚度小于7 mm的样品后界面速度都出现了二次起跳,波系图如图8所示。因为LiF窗口和金属隔层的阻抗比A型炸药的阻抗高,均会对冲击波产生反射;B型炸药处于高温高压的爆轰状态时,在金属前界面也会反射冲击波。利用式(8)、式(9)、式(10)可以估算金属隔层后界面反射的冲击波二次起跳的间隔时间
Δt 。取没有空气间隙、样品厚度为3 mm的后界面粒子速度692.48 m/s,代入式(8)、式(9)、式(10),得到冲击二次起跳间隔时间为0.94μs ,如果是金属隔层前界面反射的冲击波,二次起跳间隔时间还会增加1μs 左右。表5列出了从图像得到的二次起跳的间隔时间,和理论计算基本一致,所以二次起跳是由金属隔层后界面反射的冲击波引起的。表 5 起跳间隔时间Table 5. Interval between two accelerationsStep thickness/mm Interval/μs 1st test 2nd test 2 0.760 0.559 3 0.900 0.882 4 1.335 1.300 5 1.625 1.629 t3=hD3 (8) Δt1=h−utt3D3−u3+ut+h−utt3D3 (9) Δt2=h−ut(t3+Δt1)D3 (10) 式中:t3为冲击波从金属隔层后界面到达LiF窗口的时间,h为样品台阶厚度,D3为样品冲击波速度,
Δt1 为冲击波从LiF窗口到达金属后界面的时间,Δt2 为冲击波从金属后界面反射后追到样品粒子的时间。由表5可知,台阶样品厚度从2 mm增加到5 mm时,二次起跳的时间间隔逐渐增加,这是由于随着台阶厚度的增加,冲击波反射时间变长。在台阶样品厚度为7 mm和10 mm时,没有观察到二次起跳现象,这是因为相较于前面的台阶,7 mm和10 mm台阶相对较厚,冲击波追赶粒子所需的时间更长,所以在光纤探头能够测量的时间范围内冲击波未能追到后界面的粒子。
由于这4次实验都是利用飞片起爆A型炸药,然后通过空气间隙和金属隔层,所以Shot 01、Shot 02和Shot 03、Shot 04实验中金属隔层的速度相同。根据金属后界面自由速度对台阶型B型炸药进行阻抗匹配,将B型炸药参数(ρ = 1.893 g/cm3, c = 2.418 77 km/s, λ = 2.139 61)[15]代入式(3),得到B型炸药的前界面粒子速度u3为0.980 km/s,冲击波速度D3为4.515 km/s,入射压强p3为8.389 GPa。 如图6所示,仅台阶厚度为10 mm的B型炸药样品的后界面速度有上升趋势,说明当入射压力为8.389 GPa时,反应发生在样品厚度为7~10 mm之间的某个位置。
3. 结 论
在火炮平台上进行了一系列一维平面冲击实验。采用PDV测速系统对金属隔层自由界面速度和B型炸药后界面速度进行了测量,得到了速度随时间变化曲线,由此分析得出以下结论。
(1)运用阻抗匹配方法,根据金属隔层后界面速度和A型炸药爆轰粒子速度,得到了金属隔层的透射压力和入射压力,进而计算得到金属隔层厚度为5 mm时冲击波的透射率为72.35%;
(2)传爆药和金属隔层之间的空气间隙使原本的冲击压缩过程转变为准等熵压缩和冲击压缩两个过程,同时使冲击波的幅值减小;
(3)金属隔层和B型炸药间存在较小空气间隙时,会使金属隔层产生加速,炸药更加容易反应;
(4)当金属隔层厚度为5 mm,并且存在空气间隙时,A型炸药不能使厚度为10 mm的B型炸药达到爆轰状态。
-
表 1 实验条件
Table 1. Experimental conditions
Shot No. Dimensions/(mm × mm) Sample Booster Air gap Metal compartment Air gap 01 ∅50 × 30 ∅50 × 30 ∅50 × 5 02 ∅50 × 30 ∅50 × 30 ∅50 × 5 03 ∅50 × 30 ∅50 × 30 ∅50 × 5 ∅50 × 0.22 Explosive B 04 ∅50 × 30 ∅50 × 30 ∅50 × 5 Explosive B 表 2 A型炸药和JWL状态方程参数[14]
Table 2. Parameters of explosive A and JWL equation of state
A/GPa B/GPa R1 R2 ¯ω ρ1/(g·cm–3) Dj/(km·s–1) uj/(km·s–1) cj/(km·s–1) pj/GPa 934.770 12.723 4.6 1.1 0.37 1.863 8.87 2.22 6.65 36.8 表 3 第1次和第2次金属自由界面速度测试计算结果
Table 3. Computation results of the first and second test for free surface velocities of metal
Shot No. u2/(km·s–1) D2/(km·s–1) p2in/ GPa p2out/ GPa Shot No. u2/(km·s–1) D2/(km·s–1) p2in/ GPa p2out/ GPa 01 1.061 4.703 68.828 45.675 02 0.982 4.654 68.828 41.817 1.113 4.734 68.828 48.193 1.058 4.700 68.828 45.490 0.935 4.625 68.828 39.554 1.125 4.742 68.828 48.804 1.118 4.738 68.828 48.475 1.182 4.777 68.828 51.644 1.133 4.746 68.828 49.181 1.095 4.723 68.828 47.304 表 4 金属自由界面速度测试
Table 4. Test of free surface velocities of metal
Shot No. p2in/GPa p2out/GPa η1/% 01 67.230 45.675 67.94 67.230 48.193 71.68 67.230 39.554 58.83 67.230 48.475 72.10 67.230 49.181 73.15 67.230 47.304 70.04 02 67.230 41.817 62.20 67.230 45.490 67.66 67.230 48.804 72.59 67.230 51.644 76.82 表 5 起跳间隔时间
Table 5. Interval between two accelerations
Step thickness/mm Interval/μs 1st test 2nd test 2 0.760 0.559 3 0.900 0.882 4 1.335 1.300 5 1.625 1.629 -
[1] 姬广富. 高能钝感炸药分子和晶体的结构和性能的理论研究 [D]. 南京: 南京理工大学, 2002.JI G F. Theoretical study on the structure and properties of molecules and crystalline of insensitive high explosives [D]. Nanjing: Nanjing University of Science and Technology, 2002. [2] 章冠人. 凝聚炸药起爆动力学 [M]. 北京: 国防工业出版社, 1991.ZHANG G R. Initiation kinetics of condensed explosives [M]. Beijing: National Defense Industry Press, 1991. [3] 陈朗, 伍俊英, 方青, 等. 固体炸药冲击起爆研究 [J]. 火炸药学报, 2004, 27(1): 1–4. doi: 10.3969/j.issn.1007-7812.2004.01.001CHEN L, WU J Y, FANG Q, et al. Investigation of solid explosives initiation under shock waves [J]. Chinese Journal of Explosives & Propellants, 2004, 27(1): 1–4. doi: 10.3969/j.issn.1007-7812.2004.01.001 [4] 胡湘渝, 黄风雷, 张德良. 铸装B炸药的隔板实验数值模拟 [C]//全国计算流体力学会议, 2000: 115–119.HU X Y, HUANG F L, ZHANG D L. Numerical simulation of separator experiment for cast B explosives [C]//National Conference on Computational Fluid Dynamics, 2000: 115–119. [5] 郑波, 袁鹏刚. 验证板反射冲击波对固体推进剂冲击起爆影响的实验研究 [J]. 固体火箭技术, 2014, 37(5): 734–736.ZHENG B, YUAN P G. Experimental study of influence of reflected shock from witness plate on propellant SDT [J]. Journal of Solid Rocket Technology, 2014, 37(5): 734–736. [6] 向梅, 黄毅民, 饶国宁, 等. 复合装药结构隔板实验与数值模拟 [J]. 兵工学报, 2013, 34(2): 246–250.XIANG M, HUANG Y M, RAO G N, et al. Experimental and numerical simulation study of the shockwave sensitivity of composite charge explosive [J]. Acta Armamentarii, 2013, 34(2): 246–250. [7] 蔡进涛, 王桂吉, 赵剑衡, 等. 固体炸药的磁驱动准等熵压缩实验研究 [J]. 高压物理学报, 2010, 24(6): 455–460. doi: 10.11858/gywlxb.2010.06.009CAI J T, WANG G J, ZHAO J H, et al. Magnetically driven quasi-isentropic compression experiments of solid explosives [J]. Chinese Journal of High Pressure Physics, 2010, 24(6): 455–460. doi: 10.11858/gywlxb.2010.06.009 [8] GUSTAVSEN R L, DATTELBAUM D M, SHEFFIELD S A, et al. Application of photonic Doppler velocimetry (PDV) to shock and detonation physics experiments: LA-UR-11-00886 [R]. Los Alamos, NM (United States): Los Alamos National Laboratory, 2011. [9] SARGIS P D, MOLAU N E, SWEIDER D, et al. Photonic Doppler velocimetry: UCRL-ID-133075 [R]. Livermore: Lawrence Livermore National Laboratory, 1999. [10] 李建中, 王德田, 刘俊, 等. 多点光子多普勒测速仪及其在爆轰物理领域的应用 [J]. 红外与激光工程, 2016, 45(4): 0422001.LI J Z, WANG D T, LIU J, et al. Multi-channel photonic Doppler velocimetry and its application in the field of explosion physics [J]. Infrared and Laser Engineering, 2016, 45(4): 0422001. [11] 奥尔连科 Л П. 爆炸物理学 [M]. 孙承纬, 译. 北京: 科学出版社, 2011.OРЛЕНКО Л П. Explosive physics [M]. Translated by SUN C W. Beijing: Science Press, 2011. [12] 王礼立. 应力波基础 [M]. 第2版. 北京: 国防工业出版社, 2005.WANG L L. Foundation of stress waves [M]. 2nd ed. Beijing: Defense Industry Press, 2005. [13] 张震宇, 田占东, 陈军, 等. 爆轰物理 [M]. 长沙: 国防科技大学出版社, 2016.ZHANG Z Y, TIAN Z D, CHEN J, et al. Detonation physics [M]. Changsha: National University of Defense Technology Press, 2016. [14] 向梅, 黄毅民, 韩勇, 等. JO-9159与JB-9014复合药柱爆轰驱动平面飞片实验与数值模拟 [J]. 高压物理学报, 2014, 28(3): 379–384. doi: 10.11858/gywlxb.2014.03.018XIANG M, HUANG Y M, HAN Y, et al. Experimental study and numerical simulation of plane flyer driven by detonation of JO-9159 and JB-9014 composite charge [J]. Chinese Journal of High Pressure Physics, 2014, 28(3): 379–384. doi: 10.11858/gywlxb.2014.03.018 [15] 刘俊明, 张旭, 赵康, 等. 用PVDF压力计研究未反应JB-9014钝感炸药的Grüneisen参数 [J]. 高压物理学报, 2018, 32(5): 051301.LIU J M, ZHANG X, ZHAO K, et al. Using PVDF gauge to study Grüneisen parameter of unreacted JB-9014 insensitive explosive [J]. Chinese Journal of High Pressure Physics, 2018, 32(5): 051301. -