未反应炸药JOB-9003的JWL状态方程*

王延飞，刘 杰，张 旭，张 蓉，钟 斌

(中国工程物理研究院流体物理研究所冲击波物理与爆轰物理实验室，四川绵阳 621999)

状态方程是炸药的基本性能，对研究炸药的冲击起爆、确定化学反应速率函数以及爆轰数值模拟是不可或缺的。目前，主要采用半理论半经验的方法进行状态方程研究，即先测出炸药的冲击Hugoniot线，再根据实测的Hugoniot线对理论模型中的参数进行标定。

在未反应炸药的冲击Hugoniot线研究方面，国外通常采用轻气炮加载技术，测试技术则以VISAR(Velocity Interferometer System for Any Reflector)和组合式电磁粒子速度计为主。例如：Millett等人[1]利用一级轻气炮和PVDF(Polyvinylidene Fluoride)压力计研究了高聚物粘结炸药的冲击Hugoniot线；Burns等人[2]采用电磁粒子速度计对EDC32炸药进行了一维平面冲击实验，获得了0.59～7.5 GPa压力范围内分段式线性Hugoniot关系。国内未反应炸药的冲击Hugoniot线研究则主要采用火炮加载和平面波透镜加载，测试技术从最初的锰铜压力计发展到现在的电磁粒子速度计。例如：张旭等人[3]采用锰铜压力计测量了JB9014钝感炸药和LY12铝标准样品在LY12铝飞片同时撞击下的界面压力，通过压力对比法确定了JB9014钝感炸药的冲击Hugoniot线；傅华等人[4]利用粒子速度对比法，通过电磁粒子速度计测量了PMMA(Polymethyl Methacrylate)飞片同时撞击PMMA样品和JOB-9003炸药界面的粒子速度，获得了未反应JOB-9003炸药的冲击Hugoniot关系。

目前，炸药的状态方程大多采用三项式状态方程，为此本研究利用实测的JOB-9003炸药的冲击Hugoniot线，对三项式状态方程进行标定，以确定未反应炸药JOB-9003的JWL形式的状态方程参数。

2 实验方法

采用火炮加载技术和组合式电磁粒子速度计，对JOB-9003炸药进行一维平面冲击实验。如图1所示，实验装置主要包括加载装置、真空靶室、∅57 mm的Lexon弹托、∅55 mm的蓝宝石飞片、磁通量为0.14 T的磁场装置、组合电磁粒子速度计、待测JOB-9003炸药样品。将铝基组合式电磁粒子速度计镶入具有30°倾角的炸药样品中，感应单元之间的间隔为2 mm，以保证相邻感应单元在冲击波传播方向的间距为1 mm。采用光测技术测量飞片速度，当飞片通过测速环时，先后阻断间隔为25 mm的两束信号光，通过示波器反馈的时间差计算飞片速度。组合电磁粒子速度计是基于法拉第电磁感应定律，若已知磁场强度B和切割磁力线的感应单元长度l，通过测量感生电动势ε，即可计算粒子速度up铝基组合式电磁粒子速度计(见图2)由1个8测点组合速度计、3个单点速度计和3个冲击波示踪器构成。单点速度计贴在炸药的前表面，用于记录冲击前端的界面粒子速度；锯齿状元件为冲击波示踪器，相邻感应单元的间距为0.5 mm，总长20 mm，共有40个测量位置，用于记录冲击波到达样品中各测量位置的时刻。

图1 实验装置示意图Fig.1 Schematic illustration of experimental devices

图2 铝基组合式电磁粒子速度计Fig.2 Al-based electromagnetic particle velocity gauge

3 实验结果与分析

随着冲击波深入炸药，炸药中的粒子速度峰值不断增加，直至冲击转为爆轰。由前表面的单点速度计可测得界面处的粒子速度，根据冲击波到达0、1、2 mm处的时刻，可推算出界面处的冲击波速度。图3为组合式电磁粒子速度计在示波器上显示的典型原始信号，图4为典型的冲击波示踪器信号。通过(1)式，可获得炸药的粒子速度up随时间t变化的曲线，如图5～图8所示。在图5中，由于1、2 mm处的粒子速度信号丢失，因此采用冲击波示踪器信号(见图4)计算界面处的冲击波速度。

图3 粒子速度计记录的原始信号Fig.3 Original particle velocity signals

图4 冲击波示踪器信号Fig.4 Tracer signals

图5 6.28 GPa压力下的粒子速度波形Fig.5 Particle velocity waveform at 6.28 GPa

图6 6.10 GPa压力下的粒子速度波形Fig.6 Particle velocity waveform at 6.10 GPa

图7 4.17 GPa压力下的粒子速度波形Fig.7 Particle velocity waveform at 4.17 GPa

图8 3.11 GPa压力下的粒子速度波形Fig.8 Particle velocity waveform at 3.11 GPa

图9 JOB-9003炸药的冲击Hugoniot曲线Fig.9 Hugoniot curve of JOB-9003 explosive

由于平面冲击波是一种有物质流通过的间断面，由动量守恒定律可知

式中：pH为初始冲击压力，ρ0为炸药的初始密度，us为入射冲击波速度，up为波后界面粒子速度。在界面处，由(2)式可以得到每发实验的初始冲击压力pH，进而得到4发实验的us和up，如表1所示。

利用us-up实验数据，绘制出实验压力范围内炸药的us-up关系

式中：C0和λ为常数。对4组实验数据进行线性拟合，得到JOB-9003炸药的冲击Hugoniot关系(见图9)

表1 JOB-9003炸药平面冲击实验参数Table 1 Parameters of plane impact experiments on JOB-9003 explosive

4 三项式状态方程

一般而言，无论是固体还是液体，压力p都可通过三项式状态方程表示为比容v和温度T的函数

式中：px(v)为冷压，pTN(v,T)为晶格热振动贡献，pTe(v,T)为自由电子项贡献。对于非金属材料，可以不考虑自由电子项贡献，晶格热振动贡献可通过固体德拜模型确定，而冷压则由Grüneisen物态方程和热力学关系式求得。Grüneisen物态方程是描述晶格热振动贡献的一种特性方程

式中：E为内能，px和Ex分别为冷能和冷压，γ(v)为Grüneisen参数。将(3)式代入平面冲击波的质量守恒、动量守恒和能量守恒关系式，则(2)式可写为

式中：η=1-v/v0，其中v0为炸药的初始比容。通过px和Ex之间的关系以及冲击的能量守恒方程，可将(6)式变换为

图曲线

由(10)式可知，在本实验压力范围内γ(v)的变化很小，取定值2.76。拟合得到的冷压曲线见图10。

三项式状态方程中的晶格热振动贡献可由德拜模型得到[6]

式中：w是与温度有关的量，cV为定容比热容。当冲击压力较低时，cV和w可取定值。对于JOB-9003炸药，cV取1.02 J/(g·K)[7]；由初始条件p(v0,T0)=0，可得到w。联立(5)式和(11)式，得到本实验压力范围内简化的JWL形式状态方程

式中：A、B、C、D为拟合参数，如表2所示。

表2 JOB-9003炸药的状态方程参数Table 2 Equation of state parameters of JOB-9003 explosive

5 结 论

利用火炮加载技术及组合式电磁粒子速度计，进行了JOB-9003炸药的爆轰反应研究。实验原理简单，信噪比高，干扰小，对于研究低压下炸药的爆轰反应是可行的。对JOB-9003炸药进行三项式状态方程标定，在实验压力范围内与实测的Hugoniot曲线相近，效果较好。

[1] MILLETT J C F,BOURNE N K.The shock Hugoniot of a plastic bonded explosive and inert simulants [J].J Phys D,2004,37(18):2613-2617.

[2] BURNS M J,GUSTAVSEN R L,BARTRAM B D.One-dimensional plate impact experiments on the cyclotetramethylene tetranitramine (HMX) based explosive EDC32 [J].J Appl Phys,2012,112(6):064910.

[3] 张 旭,池家春,冯民贤.JB9014钝感炸药冲击绝热线测量 [J].高压物理学报,2001,15(4):304-308.

ZHANG X,CHI J C,FENG M X.Hugoniot relation of JB9014 insensitive high explosive [J].Chinese Journal of High Pressure Physics,2001,15(4):304-308.

[4] 傅 华,谭多望,李金河,等.未反应JOB-9003炸药冲击Hugoniot关系测试 [J].高压物理学报,2009,23(6):427-432.

FU H,TAN D W,LI J H,et al.Hugoniot relation of unreacted JOB-9003 explosive [J].Chinese Journal of High Pressure Physics,2009,23(6):427-432.

[5] 张宝坪,张庆明,黄风雷.爆轰物理学 [M].北京:兵器工业出版社,1997:375-402.

ZHANG B P,ZHANG Q M,HUANG F L.Detonation physics [M].Beijing:Weapon Industry Press,1997:375-402.

[6] 经福谦.实验物态方程导引 [M].第2版.北京:科学出版社,1999:29-38.

JING F Q.Introduction to experimental equation of state [M].2nd ed.Beijing:Science Press,1999:29-38.

[7] 董海山,周芬芬.高能炸药及相关物性能 [M].北京:科学出版社,1989:297-299.

DONG H S,ZHOU F F.High energy explosives and correlative physical properties [M].Beijing:Science Press,1989:297-299.