TATB基非均质炸药预冲击减敏的数值模拟*

黄奎邦，刘益儒，洪 滔，于 鑫，彭文扬，舒俊翔

（1.北京应用物理与计算数学研究所，北京100094；2.中国工程物理研究院研究生院，北京100088；3.中国工程物理研究院流体物理研究所，四川 绵阳621999）

高能炸药安全性在民用和军事领域受到广泛关注，其中多次冲击转爆轰行为是炸药安全性研究的重要方向之一。在一些非理想起爆条件下，炸药可能会发生预冲击压缩行为，如隔爆、绕爆或炸药外部包裹高声速材料的复杂装置等。在弱冲击的预压缩条件下，炸药感度会下降，后续的强冲击压缩时，炸药的爆轰距离变长或无法起爆形成爆轰。研究预冲击减敏问题，对掌握多次冲击条件下装置起爆、爆轰过程或爆轰驱动的可靠性，具有重要应用价值和意义。Campbell等[1]通过聚苯乙烯和铜组合层反射波产生二次冲击的楔形实验，研究了HMX 基黏塑性炸药在3.9 GPa 预冲击下后续10 GPa 的二次冲击无法实现炸药点火的现象，Campbell 等[2]还研究了PBX9404炸药预冲击压力与需要的减敏时间的关系为p2.2τ=1 150（p 的单位为100 MPa，τ 的单位为µs）。对HMX 基、RDX 基、TNT基和TATB基等PBX 炸药的减敏现象，已开展了许多实验研究[3-6]。通过实验及理论分析[7-9]，认为弱冲击压缩清除了潜在热点是非均质炸药预冲击减敏的主要物理机制。因此，非均质炸药的预冲击压缩减敏过程的合理描述，可能需要考虑孔隙率、空穴塌陷、热点等物理机制的细观反应速率模型。但细观反应速率模型计算量大、模型中的物理参数较难确定，难推广到工程的实际应用。在唯象反应速率模型上，开展非均质炸药的预冲击压缩减敏研究，仍是主要手段之一。

需要指出，目前被广泛使用的大部分唯象反应速率模型（如Lee-Tarver 点火增长模型[10]、JTF模型[11]和WSD 模型[12]等），尚不足描述非均质炸药一次冲击与多次冲击（到相同压力）时炸药状态的差异，无法应用于减敏模拟。因此，DeOliveira 等[13]在Lee-Tarver 点火增长模型中引入减敏模块，减敏进程作用于反应率模型中的点火项和增长项，较好地模拟了LX-17炸药的拐角效应中的死区现象。Hussain 等[14]、郝鹏程等[15]也做过类似的研究，在减敏模块方法中加入了一些经验参数来模拟炸药的减敏现象，不改变点火增长反应率形式及参数。但是，对于较弱冲击条件下，点火项和增长项都被抑制，爆轰不能正常进行。Starkenberg[16]则在HVRB模型[17]中扩展了依赖冲击压力的热点密度项和依赖压力的衰退率项，并引入表征冲击波产生的伪熵，进一步将热点密度项修改为伪熵函数即XHVRB模型，伪熵函数的引入使该模型具备模拟炸药多次冲击起爆的能力。近来逐渐发展了依赖熵或温度的唯象反应速率模型，如CREST 模型[18]、WSD(T)模型[19]及AWSD模型[20]等。这类模型无需增加减敏模块，用熵或温度表征一次冲击与多次冲击时炸药状态的差异，从而实现减敏的预测。

本文中，将AWSD 模型耦合进二维结构网格拉氏弹塑性流体力学程序，对主要成分为TATB的一种TATB基非均质炸药的预冲击起爆现象开展数值模拟研究。利用炸药及其产物的冲击雨贡纽实验数据，校验未反应炸药及产物的状态方程参数，通过一维冲击起爆的模拟，标定反应速率模型参数。然后，模拟该炸药的二次冲击减敏问题及其他一些预冲击减敏问题，以期为非均质炸药预冲击减敏现象的数值模拟提供参考。

1 计算方法

采用二维结构网格拉氏弹塑性流体力学程序，它具有空间一阶离散精度和时间二阶离散精度，动量方程和大多数状态方程与能量方程的联立求解采用显示格式，具有较高的计算效率。以下重点介绍耦合的反应速率模型方程、未反应炸药和产物状态方程及混合法则。

1.1 反应速率模型方程

在WSD模型基础上，AWSD 模型中引入冲击波阵面后未反应炸药的温度TSH，使模型同时基于冲击压力和温度，适合模拟不同初始温度条件下的炸药冲击到爆轰及直径效应实验[20]。总反应速率方程为：

式中：Fp为压力p 的函数，表征当地压力对反应速率的影响；F1和F2为冲击温度TSH的函数，分别控制冲击起爆过程反应率以及爆轰阶段反应率；Fλ为控制反应快结束时的慢反应项。分别为：

1.2 状态方程及参数

1.2.1未反应物状态方程及参数

炸药未反应物采用Mie-Grüneisen 形式的Davis状态方程[20]，以等熵线为参考线，未反应炸药状态方程为：

式中：E0为初始比内能，A、B、C 为根据未反应炸药冲击Hugoniot 数据标定的常数。

未反应炸药的温度为：

本文的炸药冲击雨贡纽实验数据来自文献[21-23]，实验样品中炸药初始密度范围为1.883～1.898 g/cm3，炸药的D-u关系为D=2.417 km/s+2.141 u。由文献[24]，按实验样品炸药初始密度中间值ρ0=1.890 g/cm3标定参数，标定的未反应炸药状态方程参数见表1。

表1 未反应TATB 基炸药的Davis 状态方程参数Table1 Davis EOSparametersof unreacted TATB-based explosive

参数a、k、vc、n和b通过拟合爆轰产物冲击雨贡纽实验数据来标定。标定的原则为：（1）等熵线必须过CJ 点；（2）等熵线在CJ 点处必须与Rayleigh 线相切；（3）截断体积需简单地用Gurney 能量矫正；（4）气体产物完全膨胀时对外做的功等于化学能；（5）在vg较大时，等效绝热指数等于理想气体的；（6）遵循前5个原则确定了a、k、vc、pc和n后，通过调整参数b使超压爆轰条件下产物冲击雨贡纽曲线位于炸药冲击雨贡纽曲线上方。参考实验数据[25]，标定得到炸药产物状态方程参数见表2。标定的未反应炸药及产物状态方程p-v 雨贡纽曲线与实验结果的比较如图1所示，在高压段，未反应炸药p-v 曲线低于爆轰产物p-v 曲线。

表2 TATB基炸药产物的Davis 状态方程参数Table 2 Davis EOSparametersof reaction products of TATB-based explosive

图1 未反应炸药和产物的雨贡纽曲线Fig.1 Hugoniot curves of unreacted explosives and products

1.2.3混合法则

采用温度压力平衡的混合法则，假设固体未反应炸药与气体产物之间压力平衡且温度平衡。λ 为反应度，有：

2 一次冲击起爆数值模拟及模型参数标定

参考AWSD模型模拟PBX9502炸药的模型参数[10]，根据本文的TATB基炸药冲击起爆实验数据的POP图[28-29]，对个别参数进行微调。模型参数见表3，计算模型为长2 cm 铝飞片以不同初始速度撞击长4 cm 炸药，炸药初始密度为1.890 g/cm，一维冲击起爆实验及模拟的POP图见图2。拟合得到实验的冲击压力与到爆轰距离的关系为lg(p0/GPa)=(1.44±0.04)−(0.34±0.04)lg(L*/mm)，数值模拟的冲击压力与到爆轰距离的关系为lg(p0/GPa)=(1.45±0.004)−(0.35±0.004)lg(L*/mm)，图3为入射压力p0=12.514 GPa条件下，不同位置炸药粒子速度历程曲线的模拟结果与实验结果的对比。数值模拟的爆轰增长行为与实验结果符合较好，但是波速比实验结果偏高，分析认为与未反应炸药冲击雨贡纽关系有关。张琪敏等[28]统计了锰铜法、反向撞击法、电磁法及PDV 方法测D-u关系数据，认为不同方法测量的实验结果存在一定差异，采用不同类型实验数据来标定未反应炸药状态方程对模拟冲击波在炸药中的传播波速会存在偏差。

表3 TATB基炸药AWSD模型参数Table 3 Parametersof the AWSDmodel for the TATB-based explosive

图2 冲击压力与到爆轰距离的关系Fig.2 Relationships between impact pressure and run-to-detonation distance

图3 冲击起爆粒子速度曲线Fig.3 Particle velocity curves of shock to detonation

3 预冲击减敏现象的数值模拟

3.1 TATB基炸药二次冲击起爆实验的数值模拟

二次压缩实验在口径57 mm 火炮上实施，其原理与另一种炸药的二次压缩实验[30]相同。飞片由蓝宝石和约1 mm 的Kel-F粘合而成，由于蓝宝石波阻抗远高于Kel-F，撞击炸药后向Kel-F传入的冲击波在蓝宝石界面形成高压力反射波，反射波再入射炸药形成二次冲击压缩。典型的二次冲击压缩炸药粒子速度曲线如图4所示。由图4(a)可见：前驱波为不会形成爆轰的弱波，波后压力为6.136 GPa；前驱波入射后约0.45µs，第二次入射波进入炸药，其波后压力为12.650 GPa。第二次压缩波未追赶上前驱波前，炸药无明显反应增长行为。第二次压缩波在距离炸药端面约4 mm 位置，追上第一次入射波。从第二次压缩波进入炸药开始，到爆轰距离为9.85 mm，到爆轰时间为1.72µs。图3为入射压力12.514 GPa的一次冲击起爆的粒子速度曲线，入射压力比二次加载的主波略小，一次冲击起爆到爆轰距离为8.54 mm，到爆轰时间为1.60µs，可以判断预冲击压缩对炸药的减敏作用。

图5为二次冲击起爆数值模拟结果，当主波进入炸药后且未追上前驱波前，粒子速度无增长行为，而追上前驱波后，反应才明显发生并发展为爆轰。表4为二次冲击起爆实验和数值模拟的到爆轰距离L*和到爆轰时间T*，L*和T*均以主波进入炸药开始计。由表4可见，数值模拟结果与实验结果较好符合。图6为一次冲击起爆、二次冲击起爆不同位置计算单元化学反应份额的变化曲线，二次冲击时：距离炸药端面2 mm 处，受前驱波的作用发生减敏，反应份额的增长为一条缓慢曲线，可以理解为受压力和温度影响的缓慢燃烧过程；距离炸药端面5 mm 处，单元不受前驱波影响，反应份额迅速增长，理解为受冲击作用的剧烈燃烧过程。一次冲击时，距离炸药端面2、5 mm 的计算单元反应份额增长与二次冲击距离端面5 mm 的单元类似，都为迅速增长。这说明，AWSD 模型在不影响弱冲击起爆时又可以较好模拟预冲击减敏过程。图7为不同位置的温度变化，多次冲击时，炸药单元（2 mm 位置处）温升比只经过一次冲击（5 mm 位置）的温升低，因此，依赖于温度的AWSD模型可以区分多次冲击炸药状态的差异。

图4 二次冲击压缩炸药粒子速度曲线的实验结果Fig.4 Experimental particle velocity curves of explosive by double-shock compression

图5 二次冲击压缩炸药粒子速度曲线的数值模拟结果Fig.5 Simulated particle velocity curves of explosive by double-shock compression

表4 二次冲击起爆实验的到爆轰距离和到爆轰时间的实验和数值模拟结果Table 4 Experimental and simulated results of distance and timeof run to detonation by double-shock initiation

图6 不同位置计算单元的反应份额Fig.6 Reaction fractions of numerical elements at different positions

图7 不同位置计算单元的温度Fig.7 Temperatures of numerical elements at different positions

3.2 TATB基炸药拐角效应的数值模拟

钝感炸药爆轰波经过拐角后，由于侧向稀疏作用，爆轰波被稀疏为5～9 GPa 范围的无法起爆炸药的弱冲击波，使拐角附近炸药经历一次预冲击压缩；随着后续爆轰波的发展，也逐渐扩展到拐角区域炸药，但受到弱波预冲击区域的炸药始终无法形成爆轰，即为死区现象[31]。本文中采用文献[13,15]中的模拟LX-17炸药爆轰波过拐角的凸井实验[32]的相似模型。图8为实验图像[32]，由于炸药比LX-17略微钝感，且采用拉氏方法模拟，完全按照凸井实验建模会出现大变形问题，所以将模拟模型作了适当调整，但是不改变爆轰波过拐角的爆轰过程。图9为计算模型，起爆器区域给一个高压区并采用比体积起爆模拟雷管组件引爆起爆器，其他区域炸药分别采用AWSD模型和WSD模型[12]模拟，炸药起爆后爆轰波沿着金属钨滑移传播，到达金属钨右上角后发生绕爆。网格尺度为0.1 mm，可以保证一个反应区内有15 个网格。

图8 LX-17拐角效应[32]Fig.8 LX-17 corner turning effect[32]

图9 拐角效应计算模型Fig.9 Corner-turning numerical model

图10为AWSD模型模拟结果，到4µs时爆轰波到达拐角，由于左边侧向稀疏作用，爆轰波被稀疏使拐角附近炸药经历一次预冲击压缩，随着向上的爆轰波的发展，也逐渐扩展到左边区域炸药，但是受到弱波预冲击区域的炸药始终无法形成爆轰，形成了稳定的死区。死区特征与文献[13,31-33]的实验和数值模拟结果相似。

图11为WSD模型模拟结果，爆轰波过拐角后，拐角附近炸药也经历一次预冲击压缩，但上方的爆轰波迅速向预压缩区域发展，约2µs再次形成稳定爆轰，预压缩区域的炸药也逐渐反应完成，无法形成死区。对比图10～11可见，依赖温度的AWSD 模型可以较好模拟拐角效应的死区的形成过程，而依赖压力的WSD 模型无法模拟死区现象。

图10 AWSD 模型的拐角效应特征时刻密度和反应份额分布Fig.10 Density and reaction fraction distributions at character times for corner-turning by the AWSD model

图11 WSD模型的拐角效应特征时刻密度和反应份额分布Fig.11 Density and reaction fraction distributionsat character timesfor corner-turning by the WSDmodel

4 结 论

在二维结构网格拉氏弹塑性流体力学程序中，引入基于温度和压力的AWSD反应速率模型及未反应炸药、产物的Davis状态方程。通过冲击雨贡纽实验、一维冲击起爆实验标定了状态方程参数及反应率模型参数后，针对一种TATB基钝感炸药的二次冲击起爆及拐角效应等预冲击减敏问题开展数值模拟研究，得到以下结论：

（1）二次冲击的模拟，受预压缩区域的炸药反应变慢，到爆轰距离增长约为1 mm，与该炸药二次冲击实验减敏现象相符；

（2）模拟拐角效应时，爆轰波经过拐角后，在拐角附近形成稳定的不起爆区域，与炸药主要成分相同的LX-17炸药的拐角效应实验的死区特征相符；

（3）无需引入冲击减敏模块，依赖温度的AWSD模型适用于模拟预冲击减敏现象，可以为工程中炸药预冲击减敏问题的模拟提供参考，由于模型包含冲击温度，后续还可以开展不同初始温度条件下的预冲击减敏的模拟研究。