HMX基含AP浇注炸药烤燃实验与数值模拟

刘润泽; 王昕捷; 刘瑞峰; 段卓平; 黄风雷

doi:10.11858/gywlxb.20220538

HMX基含AP浇注炸药烤燃实验与数值模拟

doi: 10.11858/gywlxb.20220538

北京理工大学机电学院, 北京　100081

基金项目: 国家自然科学基金（12172051）

详细信息

作者简介:
刘润泽（1997—），男，硕士研究生，主要从事含能材料的热安全性研究.E-mail：1517603458@qq.com

通讯作者:
王昕捷（1991—），男，副教授，博士生导师，主要从事爆炸毁伤技术研究.E-mail：wangxinjie@bit.edu.cn

中图分类号: O359; V211
计量
- 文章访问数: 273
- HTML全文浏览量: 226
- PDF下载量: 54
出版历程
- 收稿日期: 2022-03-18
- 修回日期: 2022-04-08
- 网络出版日期: 2022-09-17
- 刊出日期: 2022-10-11

Cook-off Test and Numerical Simulation of HMX-Based Cast Explosive Containing AP

School of Mechatronical Engineering, Beijing Institute of Technology, Beijing 100081, China

摘要

摘要: 为研究新型奥克托今（HMX）基浇注炸药GOLA-1（HMX-Al-高氯酸铵（AP）-黏结剂（Kel-F））在热刺激下的响应特性，对GOLA-1炸药进行了温升速率为1.0和1.5 K/min的烤燃实验，获得了炸药中心点温升和点火时间等信息。在此基础上，结合烤燃数值模拟，预测了炸药的点火位置及温度。数值模拟与实验结果吻合较好，在1.0和1.5 K/min温升速率下GOLA-1炸药点火时间的相对偏差分别为1.3%和1.7%，表明所建立的数值模型较为合理。采用该模型进行了不同温升速率下的数值模拟，结果表明：当温升速率下降至0.4 K/min时，点火位置由装药底面边缘的环形区域移动至装药中心轴线上靠近下部的位置；点火位置随温升速率的下降而逐渐向药柱上部移动，温升速率对点火温度基本没有影响。
- HMX基PBX /
- 高氯酸铵 /
- 热安全性 /
- 烤燃实验
Abstract: In order to study the response characteristics of a new HMX-based cast explosive GOLA-1 (HMX-aluminum-ammonium perchlorate (AP)-binder (Kel-F)) under thermal stimulation, the cook-off tests with heating rates of 1.0 and 1.5 K/min were conducted, and the information such as the temperature rise of the explosive center point and the ignition time was obtained. On that basis, and combined with numerical simulation of cook-off, the ignition position and temperature of explosives were predicted. The numerical simulation is in good agreement with the experimental results. The ignition time deviations of GOLA-1 explosive at heating rates of 1.0 and 1.5 K/min are 1.3% and 1.7%, respectively, indicating that the established simulation model is reasonable. On above basis, numerical simulations at different heating rates were carried out. The results show that when the heating rate decreased to 0.4 K/min, the ignition position moved from the annular area at the bottom edge of the charge to the lower part of the central axis of the charge column. As the heating rate continued to decrease, the ignition position gradually moved to the upper part of the charge, while the heating rate had little effect on the ignition temperature.
- HMX-based PBX /
- ammonium perchlorate /
- thermal safety /
- cook-off test

HTML全文

弹药在储存、运输以及使用过程中，一旦遭受意外热刺激，则易发生燃烧、爆炸甚至爆轰，因此有必要对炸药的热安全性进行合理评估。在确保弹药安全的前提下提高炸药威力是目前亟待解决的问题。在混合炸药配方中添加氧化剂可改善其氧平衡，较大幅度地提升爆热。高氯酸铵（NH₄ClO₄，ammonium perchlorate，AP）具有密度较高、含氧量高、爆轰产物为气体、机械感度低等优点，是混合炸药配方的重要候选组分^[1]。研究含有AP的浇注炸药的热响应特性对其推广应用以及新型钝感高聚物黏结炸药（plastic bonded explosives，PBX）研发具有重要意义。

国内外对含AP的含能材料热分解特性以及热安全性开展了较多研究。刘子如等^[2]研究了AP-奥克托今（HMX）混合物的热分解特性，发现AP的存在使HMX的微商热重（derivative thermogravimetry，DTG）峰温大幅下降，证明AP与HMX在热分解中存在强烈的相互作用。Li等^[3]研究了铝粉的加入对AP-HMX混合物热分解特性的影响，发现足量的铝粉可以消除AP与HMX之间的相互催化作用。陈中娥等^[4]认为，AP热分解形成的孔隙是AP-HTPB推进剂慢速烤燃响应剧烈的原因。李苗苗等^[5]对不同AP含量的HTPE推进剂进行了差示扫描量热法测试及慢速烤燃实验，发现随着AP含量的降低，推进剂的热稳定性和热安全性略有提升。Kim等^[6]对含有AP的PBX炸药进行了快速和慢速烤燃实验，结果表明，在慢速烤燃中AP不参与点火，而在快速烤燃中点火由AP与硝胺类炸药共同参与。由此可见，AP在热刺激下的自热反应受到添加剂和温升速率等因素的影响，但是目前仍缺乏含AP的PBX炸药的热刺激响应特性相关研究^[7]。

在炸药的热刺激响应特性实验中，温升速率对炸药的点火时间、点火位置以及装药的温度梯度等产生很大的影响。王沛等^[8]对GHL炸药进行了慢速烤燃实验以及4种不同温升速率下的数值模拟，发现点火位置随着温升速率的增大从炸药内部移向炸药边缘。代晓淦等^[9]进行了不同温升速率下PBX-2炸药的烤燃实验和数值模拟，发现温升速率越大，炸药中心温度在烤燃过程中的变化越小。刘静等^[10]对模块装药进行了慢速烤燃实验和数值模拟，发现随着温升速率的提高，烤燃响应区域由1个中心环形区域变成关于中心对称的2个环形区域。由此可见，温升速率对炸药烤燃中的点火位置以及点火区域形状的影响较大，对获取热刺激下炸药的响应特性具有重要意义。

本研究将开展新型浇注炸药GOLA-1的慢速烤燃实验和数值模拟，分析GOLA-1炸药在不同条件下的热刺激响应特性及其在慢速烤燃中点火的主导因素，探讨温升速率对小型烤燃弹点火位置的影响。

1. 实验与数值模拟方法

1.1 小型烤燃弹实验

基于小型烤燃弹实验装置对GOLA-1炸药进行慢速烤燃实验。图1为小型烤燃弹实验测试系统示意图。该装置主要包括壳体、端盖、炸药、环形加热器、测温热电偶以及控温系统^[11]。烤燃弹的外径为58 mm，长度为106 mm，壳体厚度为4 mm，药柱的直径为50 mm，长度为88 mm。药柱与端盖之间有高12 mm的空气域，用以缓解炸药受热膨胀对壳体的压力。端盖中心以及炸药上表面至炸药中心均开有直径为1 mm的小孔，用于放置K型热电偶，监测温升过程中炸药中心点的温度变化。壳体和端盖材质均为45钢，两者采用螺纹连接。装药壳体外侧包裹环形加热器，在慢速烤燃实验中温升速率为1.0 和1.5 K/min。加热器外部使用岩棉包裹进行保温，防止热量散失影响实验结果。在壳体与加热器之间放置两根K型热电偶，用来控制加热速率以及测量壳体外壁面的温度变化。

图 1 小型烤燃弹实验装置示意图

Figure 1. Schematic diagram of small-scale cook-off bomb experiment device

下载: 全尺寸图片幻灯片

1.2 烤燃数值模拟

烤燃实验只能获得温度监测点的温度-时间关系，无法获取炸药内部的温度分布、点火位置、点火温度以及各组分的质量分数等信息。数值模拟可以较方便地获取以上信息，从而对炸药的热刺激响应特性进行较全面的评估。基于小型烤燃弹实验装置，建立如图2所示的数值模拟模型，其中端盖及壳体、空气域以及炸药分别用粉色、蓝色和红色表示。模型用六面体网格划分，网格尺寸为0.5 mm。

图 2 小型烤燃弹的计算网格模型

Figure 2. Computational grid model of small-scale cook-off bomb

下载: 全尺寸图片幻灯片

在直角坐标系下，热传导方程的形式为^[12]

$\rho {c_p}\frac{{\partial T}}{{\partial t}} = \lambda {\nabla ^2}T+S$

(1)

式中：ρ为密度，kg·m⁻³；c_p为比定压热容，J·kg⁻¹·K⁻¹；T为温度，K；λ为热导率，W·m⁻¹·K⁻¹；S为炸药的自热反应源项，W·m⁻³。

空气域与炸药之间的辐射热流q_r采用P1辐射模型描述^[13]

${q_{\rm r}} = \frac{1}{{3\left( {a+{\sigma _{\rm s}}} \right) - C{\sigma _{\rm s}}}}\nabla G$

(2)

式中：a为吸收系数，m⁻¹；σ_s为散射系数，m⁻¹；C为线性各向异性相位函数系数；G 为入射辐射。

GOLA-1为HMX基浇注炸药，主要成分为HMX，其他成分包括Al、 AP、黏结剂（Kel-F）。组分HMX、Al、AP、Kel-F的质量比为44∶25∶20∶11，密度为1.837 g/cm³。对于炸药中的主要成分HMX，在慢速烤燃过程中，随着温度上升，HMX首先由β相转变为δ相，而后热分解生成最终的气态产物。基于过渡态理论，HMX在烤燃过程中的4步化学反应动力学模型^[14]表示为：(1) β-HMX ↔ δ-HMX（反应1，单分子吸热反应）；(2) β-HMX+δ-HMX ↔ δ-HMX（反应2，双分子吸热反应）；(3) δ-HMX→产物（反应3，单分子吸热反应）；(4) δ-HMX+产物→产物（反应4，双分子放热反应）。反应1和反应2描述β相HMX转变为δ相HMX的晶型转变过程，反应3和反应4描述最终产物生成的点火过程。各步反应速率方程为

${r_1} = {w_\beta }\frac{{kT}}{h}\exp \left( {\frac{{T\Delta {S_{1{\text{f}}}} - {E_{1{\text{f}}}}}}{{RT}}} \right) - {w_\delta }\frac{{kT}}{h}\exp \left( {\frac{{T\Delta {S_{1{\text{r}}}} - {E_{1{\text{r}}}}}}{{RT}}} \right)$

(3)

${r_2} = {w_\beta }{w_\delta }\frac{{kT}}{h}\left[ {\exp \left( {\frac{{T\Delta {S_{2{\text{f}}}} - {E_{2{\text{f}}}}}}{{RT}}} \right) - \exp \left( {\frac{{T\Delta {S_{2{\text{r}}}} - {E_{2{\text{r}}}}}}{{RT}}} \right)} \right]$

(4)

${r_3} = {w_\delta }{Z_3}\exp \left( { - \frac{{{E_3}}}{{RT}}} \right)$

(5)

${r_4} = {w_\delta }{w_{\rm p} }{Z_4}\exp \left( { - \frac{{{E_4}}}{{RT}}} \right)$

(6)

式中：r₁～r₄为反应1～反应4的反应速率；w_β和w_δ分别为β相和δ相HMX的质量分数，w_p为产物的质量分数；k为玻尔兹曼常数；h为普朗克常数；R为普适气体常数；ΔS为活化熵，J∙mol⁻¹∙K⁻¹；E为活化能，J∙mol⁻¹；下标f代表正反应，r代表逆反应；Z₁～Z₄为反应1～反应4的指前因子，s⁻¹。

各组分质量分数的微分方程为

$\frac{{{\text{d}}{w_\beta }}}{{{\text{d}}t}} = - {r_1} - {r_2}$

(7)

$\frac{{{\text{d}}{w_\delta }}}{{{\text{d}}t}} = {r_1}+{r_2} - {r_3} - {r_4}$

(8)

$\frac{{{\text{d}}{w_{\text{p}}}}}{{{\text{d}}t}} = {r_3}+{r_4}$

(9)

单位时间内HMX热分解反应的反应源项S_HMX为

${S_{\rm {HMX}}} = {\rho _{\rm {HMX}}}{Q_1}{r_1}+{\rho _{\rm {HMX}}}{Q_2}{r_2}+{\rho _{\rm {HMX}}}{Q_3}{r_3}+{\rho _{\rm {HMX}}}{Q_4}{r_4}$

(10)

式中：ρ_HMX为HMX的密度，kg∙m⁻³；Q₁～Q₄为反应1～反应4的反应热，J∙kg⁻¹。

Kel-F为含氟含氯黏结剂，分解时吸收热量，使用单步化学反应动力学模型描述其热分解过程。其热分解反应速率方程为^[15]

${r_5} = {Z_5}\exp \left( { - \frac{{{E_{{\rm a}5}}}}{{RT}}} \right)$

(11)

式中：r₅、Z₅、E_a5分别为Kel-F热分解反应的反应速率、指前因子和活化能。单位时间内Kel-F热分解反应的反应源项S_Kel-F为

${S_{\text {Kel-F}}} = {\rho _{\text {Kel-F}}}{Q_5}{r_5}$

(12)

式中：ρ_Kel-F为Kel-F的密度，Q₅为Kel-F热分解反应的反应热。

AP为氧化剂，其分解反应为放热反应，采用一步反应描述其热分解过程^[16]

${\text{N}}{{\text{H}}_{\text{4}}}{\text{Cl}}{{\text{O}}_{\text{4}}} \to {{\text{O}}_{\text{2}}}{\text+}{{\text{H}}_{\text{2}}}{\text{O+HCl+}}{{\text{N}}_{\text{2}}}$

AP在240 ℃左右由斜方晶型（α-AP）转变为立方晶型（β-AP）^[17]。在慢速烤燃的温升过程中，AP的温度几乎均在240 ℃以下，因此采用α-AP的热分解反应速率方程^[18]，即

${r_6} = \frac{{{\text{d}}\alpha }}{{{\text{d}}t}} = {Z_6}\exp \left( { - \frac{{{E_{{\rm a}6}}}}{{RT}}} \right){\alpha ^n}{\left( {1 - \alpha } \right)^m}$

(13)

式中：r₆、Z₆、E_a6分别为α-AP热分解反应速率、指前因子和活化能；α为AP的热分解转化率，初始阶段为零，完全分解为1；n和m为反应动力学参数，对于α-AP，n=1，m=1.38。根据文献[19]，AP的热分解开始于150 ℃左右，因此在UDF用户子程序中设α在150 ℃以下为零，150 ℃以上为10⁻⁴，以启动对AP分解反应计算。设AP的密度为ρ_AP，则单位时间内AP热分解反应的反应源项S_AP为

${S_{\rm {AP}}} = {\rho _{\rm {AP}}}{Q_6}{r_6}$

(14)

式中：Q₆为AP热分解反应的反应热。

因为Al不参与点火前的热分解反应，所以不考虑铝粉的热分解反应源项，则炸药的总自热反应源项S为

$S = 0.44{S_{{\rm{HMX}}}}+0.11{S_{{\text{Kel-F}}}}+0.20{S_{{\rm{AP}}}}$

(15)

根据叠加原理，炸药的导热系数λ和比热容C分别表示为

$\lambda = 0.44{\lambda _{{\text{HMX}}}}+0.11{\lambda _{{\text{Kel-F}}}}+0.20{\lambda _{\rm {AP}}}$

(16)

$C = 0.44{C_{{\text{HMX}}}}+0.11{C_{{\text{Kel-F}}}}+0.20{C_{\rm {AP}}}$

(17)

表1和表2列出了炸药组分HMX、Kel-F和AP的反应动力学参数，表3列出了材料的物性参数。采用Fluent软件进行计算，设置温升边界条件为1.0和1.5 K/min，通过编写用户自定义函数UDF，将热导率、边界条件以及反应动力学参数加载到主程序中进行计算。

表 1 HMX的反应动力学参数

Table 1. Reaction kinetic parameters of HMX

Reaction No.	Z/s⁻¹	ΔS_f/ (J∙mol⁻¹∙K⁻¹)	E_f/(kJ∙mol⁻¹)	ΔS_r/ (J∙mol⁻¹∙K⁻¹)	E_r/(kJ∙mol⁻¹)	Q/(J∙g⁻¹)
1		123.0	204	89.0	189	−25
2		−41.7	102	−75.2	865	−25
3	3.16×10¹⁶		200			−120
4	1.50×10¹⁴		173			3200

下载: 导出CSV

| 显示表格

表 2 Kel-F和AP的反应动力学参数

Table 2. Reaction kinetic parameters of Kel-F and AP

Component	Reaction No.	Z/s⁻¹	E/(kJ∙mol⁻¹)	Q/(J∙g⁻¹)
Kel-F	5	9.93×10¹⁷	272.0	−5870
AP	6	6.85×10¹²	146.3	297

下载: 导出CSV

| 显示表格

表 3 材料物性参数

Table 3. Parameters of material properties

Component	Density/(kg∙m⁻³)	C/(J·kg⁻¹·K)	λ/(W·m⁻¹·K)
HMX	1850	1004.26	0.5358
Kel-F	2020	1000.43	0.0527
AP	1950	1255	0.2760
Al	2719	871	1.39
Steel	8030	502.48	16.27
Air	1.225	1006.43	0.0242

下载: 导出CSV

| 显示表格

2. 结果与讨论

2.1 烤燃实验

图3为不同温升速率下GOLA-1炸药点火后烤燃弹照片。图3(a)显示：点火后烤燃弹壳体膨胀变形，端盖被冲飞，烤燃弹底部断裂脱落，壳体螺纹完整；壳体表面有明显的烧灼痕迹，壳体内有少量未反应炸药。图3(b)显示了1.5 K/min温升速率下GOLA-1炸药点火前后烤燃弹的对比。结果显示，点火后烤燃弹壳体膨胀变形，端盖被冲飞，壳体底部与侧面有明显的烧灼痕迹，壳体内有少量未反应炸药。由实验现象判定，当温升速率为1.0和1.5 K/min时，GOLA-1炸药的热刺激响应等级均为燃烧，说明炸药在慢速烤燃条件下的热安全性较好。对比两种温升速率下炸药的响应程度，可以发现，温升速率为1.0 K/min时炸药的慢速烤燃响应更剧烈。

图 3 点火后烤燃弹实物

Figure 3. Cook-off bombs after ignition

下载: 全尺寸图片幻灯片

2.2 数值模拟

基于实验进行慢速烤燃数值计算，在1.0和1.5 K/min温升速率下实验与数值模拟得到的炸药中心测温点温度-时间曲线如图4和图5所示。可以看出，实验与数值模拟得到的温升曲线基本重合，点火温度和点火时间也基本吻合，说明数值模型及参数比较合理。图4显示，9900 s时实验测得的温升曲线的温升放缓，而后恢复正常温升速率。分析认为，此阶段的温升放缓是由于HMX在由β相转变为δ相的过程中吸收热量，温度约为451 K。在慢速烤燃实验中，炸药于11946 s时发生点火反应，点火时药柱中心点温度为488 K；计算得到的点火时间为11790 s，与实验结果的相对偏差为1.3%。图5显示：HMX相变造成的温升速率放缓出现在约6800 s，温度约为451 K；实验结果显示，炸药在约8207 s发生点火反应，点火时刻药柱中心温度约为480 K；计算得到的点火时间为8070 s，与实验结果的相对偏差为1.7%。由此可见，温升速率对点火时间的影响较大，随着温升速率的上升，点火时间缩短，然而温升速率对HMX的相变温度以及GOLA-1炸药的点火温度没有影响。

图 4 GOLA-1炸药中心点的温度-时间曲线（1.0 K/min）

Figure 4. Temperature-time curves at the center point of GOLA-1 explosive (1.0 K/min)

下载: 全尺寸图片幻灯片

图 5 GOLA-1炸药中心点的温度-时间曲线（1.5 K/min）

Figure 5. Temperature-time curves at the center point of GOLA-1 explosive (1.5 K/min)

下载: 全尺寸图片幻灯片

1.0和1.5 K/min温升速率下不同时刻炸药的温度分布如图6和图7所示。相比于炸药，壳体的比热容小而热导率大，因此壳体温度高于炸药温度，热量由壳体传递到装药内部。装药的部分热量散失至烤燃弹上部的空气域中，造成装药上部温度较低。图6(d)和图7(d)为点火时烤燃弹的温度分布，可见，点火区域位于装药底部的环形区域。这是由于药柱底部受到侧面壳体与底面的共同加热，使得装药底部温度相对较高。温度的升高使得炸药的自热分解反应速率上升，造成热量进一步累积，如此循环往复，致使点火最终发生在装药底部的环形区域。对比两者发现，在1.5 K/min的温升速率下装药的温度梯度较大，点火位置相对靠外。分析认为，这是由于炸药的导热系数远低于钢壳体，而温升速率相对较大，热量来不及传递到炸药内部造成的。

图 6 1.0 K/min温升速率下不同时刻GOLA-1炸药的温度分布

Figure 6. Temperature distribution of GOLA-1 explosive at different time with the heating rate of 1.0 K/min

下载: 全尺寸图片幻灯片

图 7 1.5 K/min温升速率下不同时刻GOLA-1炸药的温度分布

Figure 7. Temperature distribution of GOLA-1 explosive at different time with the heating rate of 1.5 K/min

下载: 全尺寸图片幻灯片

1.0 K/min温升速率下不同位置的温升曲线和质量分数曲线如图8、图9和图10所示。从图8可以看出，HMX在约9300 s发生相变，相变温度与中心点温度相同，均为451 K。对比中心测温点与点火点的温升曲线可以发现，点火点处HMX相变所造成的温度平台较小，且持续时间较短。点火点处的热量累积速率大于热量传递速率，温度上升较快，相同时刻点火点处的温度高于其他位置，从而使HMX相变时间早且速率快。临近点火时，由于HMX热分解反应释放大量的热量，造成点火点的温度急剧升高，点火温度约为517 K。图9和图10给出了中心测温点和点火点的HMX质量分数变化曲线。在点火时刻，β-HMX已经全部转变为δ-HMX，相变过程基本完成。点火时中心测温点仅有微量的HMX热分解产物，而点火位置则有质量分数为0.2的产物。说明在点火位置，第4步强放热反应已经开始进行，释放大量热量，导致温度急剧上升，主导了点火的发生。

图 8 1.0 K/min温升速率下点火点处的温度-时间曲线

Figure 8. Temperature-time curve of the ignition position with the heating rate of 1.0 K/min

下载: 全尺寸图片幻灯片

图 9 1.0 K/min温升速率下中心测温点处HMX的质量分数变化曲线

Figure 9. Mass fraction curves of HMX at the center point with the heating rate of 1.0 K/min

下载: 全尺寸图片幻灯片

图 10 1.0 K/min温升速率下点火点处HMX的质量分数变化曲线

Figure 10. Mass fraction curves of HMX at the ignition position with the heating rate of 1.0 K/min

下载: 全尺寸图片幻灯片

1.5 K/min温升速率下不同位置的温升曲线和质量分数曲线如图11、图12和图13所示。图11中的点火点温度-时间曲线特征与图8基本相同，点火温度约为518 K。图12和图13显示：相比于温升速率为1.0 K/min的情况，点火时刻中心测温点处β-HMX的质量分数为0.04，δ-HMX的质量分数为0.96，表明部分β-HMX尚未转变为δ-HMX，第4步强放热反应尚未发生；而点火位置的HMX已完全相变，并且HMX分解的第4步强放热反应已经开始进行，有质量分数为0.19的产物生成，主导了点火的发生。

图 11 1.5 K/min温升速率下点火点处的温度-时间曲线

Figure 11. Temperature-time curve of the ignition position with the heating rate of 1.5 K/min

下载: 全尺寸图片幻灯片

图 12 1.5 K/min温升速率下中心测温点处HMX的质量分数变化曲线

Figure 12. Mass fraction curves of HMX at center point with the heating rate of 1.5 K/min

下载: 全尺寸图片幻灯片

图 13 1.5 K/min温升速率下点火点处HMX的质量分数变化曲线

Figure 13. Mass fraction curves of HMX at ignition position with the heating rate of 1.5 K/min

下载: 全尺寸图片幻灯片

图14给出了温升速率为1.5 K/min时中心测温点处AP的转化率-时间曲线。转化率α与分解速率dα/dt均随着温度的升高而增大。在点火时刻，AP的转化率为1.06×10⁻⁴，表明此时中心点的AP基本不发生分解。7000 s左右，AP的分解速率放缓，而后快速上升，这是HMX在晶型转变过程中吸收热量而导致温升放缓造成的。图15为点火点处AP的转化率-时间曲线。可以看出：转化率在临近点火时刻迅速上升；点火时刻AP的转化率为2.26×10⁻⁴，相较于中心点处的转化率高约1倍。不论是中心测温点还是点火点，发生分解的AP的量均很少，表明在慢速烤燃过程中AP基本不参与反应，点火反应由HMX主导。这与Kim等^[6]得到的含硝胺类及AP浇注炸药慢速烤燃实验结论一致。

图 14 1.5 K/min温升速率下中心测温点的AP转化率曲线

Figure 14. Conversion rate curve of AP at the center point with the heating rate of 1.5 K/min

下载: 全尺寸图片幻灯片

图 15 1.5 K/min温升速率下点火点的AP转化率曲线

Figure 15. Conversion rate curve of AP at ignition position with the heating rate of 1.5 K/min

下载: 全尺寸图片幻灯片

2.3 温升速率对点火位置的影响

将实验与数值模拟相结合，分析温升速率对炸药热刺激响应特性的影响。在1.0和1.5 K/min两种温升速率的基础上，设定不同的温升速率，对GOLA-1炸药进行计算。当温升速率为0.5和0.4 K/min时，炸药点火位置的变化如图16所示。0.5 K/min时，点火位置仍为装药底部的环形区域，相比1.0和1.5 K/min的情况，环形区域缩小，并向装药的上部及中心移动。这是由于当温升速率下降时，炸药整体受热比较均匀，装药内部的温度梯度减小，相比于温升速率为1.0 K/min的情况，相同加热温度下炸药内部高温区域的体积较大，炸药的自热反应速率上升，炸药内部积累热量，导致点火发生在装药底部的环形区域。当温升速率降至0.4 K/min时，点火位置的形状发生改变，由环形变为扁的椭球形，并且位置移动至药柱下部中心。

图 16 不同温升速率下点火位置的对比

Figure 16. Comparison of ignition positions with different heating rates

下载: 全尺寸图片幻灯片

为进一步分析温升速率对GOLA-1点火位置的影响，以0.05 K/min为间隔，进行0.1～0.4 K/min 7个不同温升速率的烤燃数值模拟。设温升速率为β，以壳体外侧底面的中心位置为坐标零点，沿烤燃弹中心轴线从壳体底面指向端盖为正方向，点火位置x为温度最高区域的中心，单位为m。温升速率对点火位置的影响如表4所示。随着温升速率的下降，点火位置向药柱上部移动。根据表4绘制点火位置随温升速率的散点图，如图17所示。观察图17发现，温升速率β与点火位置x之间存在线性关系，拟合关系式为 $x = - 0.052\,79\,\beta +0.041\,71$ ，相关系数 R²=0.966。

表 4 温升速率与点火位置的关系

Table 4. Heating rate vs. ignition position

No.	β/(K·min⁻¹)	x/m
1	0.40	0.0212
2	0.35	0.0220
3	0.30	0.0273
4	0.25	0.0287
5	0.20	0.0298
6	0.15	0.0334
7	0.10	0.0374

下载: 导出CSV

| 显示表格

图 17 温升速率对点火位置的影响

Figure 17. Effect of heating rate on ignition position

下载: 全尺寸图片幻灯片

3. 结　论

(1) 对新型HMX基含AP浇注炸药GOLA-1进行慢速烤燃实验，结果显示，当温升速率为1.0和1.5 K/min时炸药的烤燃响应等级均为燃烧，并且1.0 K/min下壳体的破坏更严重，响应更剧烈。从热刺激响应角度来看，温升速率越低，GOLA-1炸药的响应越剧烈。

(2) 开展了GOLA-1炸药烤燃数值模拟，得到的温升曲线与实验结果吻合较好。当温升速率为1.0 K/min时，GOLA-1炸药的点火时间与实验结果的相对偏差为1.3%，1.5 K/min时为1.7%，证明建立的烤燃模型及参数较合理。数值模拟结果显示，点火反应主要由HMX的第4步强放热反应主导，AP在慢速烤燃过程中基本不参与点火反应。

(3) 对不同温升速率下GOLA-1炸药的烤燃过程进行了数值模拟，结果表明：当温升速率为0.5 K/min时，点火仍发生在装药底部的环形区域，并向装药中心及上部移动；当温升速率下降至0.4 K/min时，点火区域的形状发生变化，由环形变为椭球形；随着温升速率继续下降，点火位置沿着药柱中心轴线向上部移动；当温升速率处在0.1～0.4 K/min之间时，温升速率与点火位置之间存在线性关系。

图 1 小型烤燃弹实验装置示意图

Figure 1. Schematic diagram of small-scale cook-off bomb experiment device

下载: 全尺寸图片幻灯片

图 2 小型烤燃弹的计算网格模型

Figure 2. Computational grid model of small-scale cook-off bomb

下载: 全尺寸图片幻灯片

图 3 点火后烤燃弹实物

Figure 3. Cook-off bombs after ignition

下载: 全尺寸图片幻灯片

图 4 GOLA-1炸药中心点的温度-时间曲线（1.0 K/min）

Figure 4. Temperature-time curves at the center point of GOLA-1 explosive (1.0 K/min)

下载: 全尺寸图片幻灯片

图 5 GOLA-1炸药中心点的温度-时间曲线（1.5 K/min）

Figure 5. Temperature-time curves at the center point of GOLA-1 explosive (1.5 K/min)

下载: 全尺寸图片幻灯片

图 6 1.0 K/min温升速率下不同时刻GOLA-1炸药的温度分布

Figure 6. Temperature distribution of GOLA-1 explosive at different time with the heating rate of 1.0 K/min

下载: 全尺寸图片幻灯片

图 7 1.5 K/min温升速率下不同时刻GOLA-1炸药的温度分布

Figure 7. Temperature distribution of GOLA-1 explosive at different time with the heating rate of 1.5 K/min

下载: 全尺寸图片幻灯片

图 8 1.0 K/min温升速率下点火点处的温度-时间曲线

Figure 8. Temperature-time curve of the ignition position with the heating rate of 1.0 K/min

下载: 全尺寸图片幻灯片

图 9 1.0 K/min温升速率下中心测温点处HMX的质量分数变化曲线

Figure 9. Mass fraction curves of HMX at the center point with the heating rate of 1.0 K/min

下载: 全尺寸图片幻灯片

图 10 1.0 K/min温升速率下点火点处HMX的质量分数变化曲线

Figure 10. Mass fraction curves of HMX at the ignition position with the heating rate of 1.0 K/min

下载: 全尺寸图片幻灯片

图 11 1.5 K/min温升速率下点火点处的温度-时间曲线

Figure 11. Temperature-time curve of the ignition position with the heating rate of 1.5 K/min

下载: 全尺寸图片幻灯片

图 12 1.5 K/min温升速率下中心测温点处HMX的质量分数变化曲线

Figure 12. Mass fraction curves of HMX at center point with the heating rate of 1.5 K/min

下载: 全尺寸图片幻灯片

图 13 1.5 K/min温升速率下点火点处HMX的质量分数变化曲线

Figure 13. Mass fraction curves of HMX at ignition position with the heating rate of 1.5 K/min

下载: 全尺寸图片幻灯片

图 14 1.5 K/min温升速率下中心测温点的AP转化率曲线

Figure 14. Conversion rate curve of AP at the center point with the heating rate of 1.5 K/min

下载: 全尺寸图片幻灯片

图 15 1.5 K/min温升速率下点火点的AP转化率曲线

Figure 15. Conversion rate curve of AP at ignition position with the heating rate of 1.5 K/min

下载: 全尺寸图片幻灯片

图 16 不同温升速率下点火位置的对比

Figure 16. Comparison of ignition positions with different heating rates

下载: 全尺寸图片幻灯片

图 17 温升速率对点火位置的影响

Figure 17. Effect of heating rate on ignition position

下载: 全尺寸图片幻灯片

表 1 HMX的反应动力学参数

Table 1. Reaction kinetic parameters of HMX

Reaction No.	Z/s⁻¹	ΔS_f/ (J∙mol⁻¹∙K⁻¹)	E_f/(kJ∙mol⁻¹)	ΔS_r/ (J∙mol⁻¹∙K⁻¹)	E_r/(kJ∙mol⁻¹)	Q/(J∙g⁻¹)
1		123.0	204	89.0	189	−25
2		−41.7	102	−75.2	865	−25
3	3.16×10¹⁶		200			−120
4	1.50×10¹⁴		173			3200

下载: 导出CSV

表 2 Kel-F和AP的反应动力学参数

Table 2. Reaction kinetic parameters of Kel-F and AP

Component	Reaction No.	Z/s⁻¹	E/(kJ∙mol⁻¹)	Q/(J∙g⁻¹)
Kel-F	5	9.93×10¹⁷	272.0	−5870
AP	6	6.85×10¹²	146.3	297

下载: 导出CSV

表 3 材料物性参数

Table 3. Parameters of material properties

Component	Density/(kg∙m⁻³)	C/(J·kg⁻¹·K)	λ/(W·m⁻¹·K)
HMX	1850	1004.26	0.5358
Kel-F	2020	1000.43	0.0527
AP	1950	1255	0.2760
Al	2719	871	1.39
Steel	8030	502.48	16.27
Air	1.225	1006.43	0.0242

下载: 导出CSV

表 4 温升速率与点火位置的关系

Table 4. Heating rate vs. ignition position

No.	β/(K·min⁻¹)	x/m
1	0.40	0.0212
2	0.35	0.0220
3	0.30	0.0273
4	0.25	0.0287
5	0.20	0.0298
6	0.15	0.0334
7	0.10	0.0374

下载: 导出CSV

参考文献(19)

[1]	李玉斌, 黄辉, 潘丽萍, 等. 高氯酸铵的包覆降感与应用研究 [J]. 含能材料, 2014(6): 792–797. doi: 10.11943/j.issn.1006-9941.2014.06.015 LI Y B, HUANG H, PAN L P, et al. Desensitizing technology of AP by coating and its application [J]. Chinese Journal of Energetic Materials, 2014(6): 792–797. doi: 10.11943/j.issn.1006-9941.2014.06.015
[2]	刘子如, 施震灏, 阴翠梅, 等. 热红联用研究AP与RDX和HMX混合体系的热分解 [J]. 火炸药学报, 2007, 30(5): 57–61. doi: 10.3969/j.issn.1007-7812.2007.05.015 LIU Z R, SHI Z H, YIN C M, et al. Investigation on thermal decomposition of mixed systems of AP with RDX and HMX by DSC-TG-FTIR [J]. Chinese Journal of Explosives & Propellants, 2007, 30(5): 57–61. doi: 10.3969/j.issn.1007-7812.2007.05.015
[3]	LI Y B, PAN L P, YANG Z J, et al. The effect of wax coating, aluminum and ammonium perchlorate on impact sensitivity of HMX [J]. Defence Technology, 2017, 13(6): 422–427. doi: 10.1016/j.dt.2017.05.022
[4]	陈中娥, 唐承志, 赵孝彬. 固体推进剂的慢速烤燃行为与热分解特性的关系研究 [J]. 含能材料, 2005, 13(6): 393–396. doi: 10.3969/j.issn.1006-9941.2005.06.013 CHEN Z E, TANG C Z, ZHAO X B. Relationship between slow cook-off behaviour and thermal decomposition characteristics of solid propellant [J]. Chinese Journal of Energetic Materials, 2005, 13(6): 393–396. doi: 10.3969/j.issn.1006-9941.2005.06.013
[5]	李苗苗, 郑亭亭, 陈静静, 等. HMX含量对HTPE推进剂热安全性的影响 [J]. 固体火箭技术, 2020, 43(2): 229–236. doi: 10.7673/j.issn.1006-2793.2020.02.015 LI M M, ZHENG T T, CHEN J J, et al. Influence of HMX content on thermal safety characteristics of HTPE propellant [J]. Journal of Solid Rocket Technology, 2020, 43(2): 229–236. doi: 10.7673/j.issn.1006-2793.2020.02.015
[6]	KIM Y, PARK Y, YOH J J. Slow and rapid thermal decomposition characteristics of enhanced blast explosives for burning in fuel-rich, oxygen-rich conditions [J]. Thermochimica Acta, 2019, 678: 178300. doi: 10.1016/j.tca.2019.178300
[7]	王琼, 丁黎, 张冬梅, 等. 含AP的浇铸PBX炸药的热安全性 [J]. 含能材料, 2015, 23(7): 693–696. doi: 10.11943/j.issn.1006-9941.2015.07.016 WANG Q, DING L, ZHANG D M, et al. Thermal safety of casted PBX containing AP [J]. Chinese Journal of Energetic Materials, 2015, 23(7): 693–696. doi: 10.11943/j.issn.1006-9941.2015.07.016
[8]	王沛, 陈朗, 冯长根. 不同升温速率下炸药烤燃模拟计算分析 [J]. 含能材料, 2009, 17(1): 46–49, 54. doi: 10.3969/j.issn.1006-9941.2009.01.012 WANG P, CHEN L, FENG C G. Numerical simulation of cook-off for explosive at different heating rates [J]. Chinese Journal of Energetic Materials, 2009, 17(1): 46–49, 54. doi: 10.3969/j.issn.1006-9941.2009.01.012
[9]	代晓淦, 黄毅民, 吕子剑, 等. 不同升温速率热作用下PBX-2炸药的响应规律 [J]. 含能材料, 2010, 18(3): 282–285. doi: 10.3969/j.issn.1006-9941.2010.03.010 DAI X G, HUANG Y M, LYU Z J, et al. Reaction behavior for PBX-2 explosive at different heating rate [J]. Chinese Journal of Energetic Materials, 2010, 18(3): 282–285. doi: 10.3969/j.issn.1006-9941.2010.03.010
[10]	刘静, 余永刚. 不同升温速率下模块装药慢速烤燃特性的数值模拟 [J]. 兵工学报, 2019, 40(5): 990–995. doi: 10.3969/j.issn.1000-1093.2019.05.011 LIU J, YU Y G. Simulation of slow cook-off for modular charges at different heating rates [J]. Acta Armamentarii, 2019, 40(5): 990–995. doi: 10.3969/j.issn.1000-1093.2019.05.011
[11]	吴浩, 段卓平, 白孟璟, 等. DNAN基含铝炸药烤燃实验与数值模拟 [J]. 含能材料, 2021, 29(5): 414–521. doi: 10.11943/CJEM2020298 WU H, DUAN Z P, BAI M J, et al. Small-scale cook-off experiments and simulations of DNAN-based aluminized explosives [J]. Chinese Journal of Energetic Materials, 2021, 29(5): 414–521. doi: 10.11943/CJEM2020298
[12]	冯长根. 热爆炸理论 [M]. 北京: 科学出版社, 1988.
[13]	寇永锋, 陈朗, 马欣, 等. 黑索今基含铝炸药烤燃试验和数值模拟 [J]. 兵工学报, 2019, 40(5): 978–989. doi: 10.3969/j.issn.1000-1093.2019.05.010 KOU Y F, CHEN L, MA X, et al. Cook-off experimental and numerical simulation of RDX-based aluminized explosives [J]. Acta Armamentarii, 2019, 40(5): 978–989. doi: 10.3969/j.issn.1000-1093.2019.05.010
[14]	PERRY W L, GUNDERSON J A, BALKEY M M, et al. Impact-induced friction ignition of an explosive: infrared observations and modeling [J]. Journal of Applied Physics, 2010, 108(8): 084902. doi: 10.1063/1.3487932
[15]	TARVER C M, KOERNER J G. Effects of endothermic binders on times to explosion of HMX- and TATB-based plastic bonded explosives [J]. Journal of Energetic Materials, 2007, 26(1): 1–28. doi: 10.1080/07370650701719170
[16]	KIM K H, KIM C K, YOO J C, et al. Test-based thermal decomposition simulation of AP/HTPB and AP/HTPE propellants [J]. Journal of Propulsion and Power, 2011, 27(4): 822–827. doi: 10.2514/1.B34099
[17]	张端庆. 火药用原材料性能与制备 [M]. 北京: 北京理工大学出版社, 1995.
[18]	RAJIĆ M, SUĆESKA M. Study of thermal decomposition kinetics of low-temperature reaction of ammonium perchlorate by isothermal TG [J]. Journal of Thermal Analysis and Calorimetry, 2000, 63(2): 375–386. doi: 10.1023/A:1010136308310
[19]	BOLDYREV V V. Thermal decomposition of ammonium perchlorate [J]. Thermochimica Acta, 2006, 443(1): 1–36. doi: 10.1016/j.tca.2005.11.038

施引文献

资源附件(0)

访问统计

点击查看大图

图(17) / 表(4)

计量

文章访问数: 273
HTML全文浏览量: 226
PDF下载量: 54

1. 实验与数值模拟方法
1.1 小型烤燃弹实验
1.2 烤燃数值模拟
2. 结果与讨论
2.1 烤燃实验
2.2 数值模拟
2.3 温升速率对点火位置的影响
3. 结　论

HMX基含AP浇注炸药烤燃实验与数值模拟

doi: 10.11858/gywlxb.20220538

作者简介: 刘润泽（1997—），男，硕士研究生，主要从事含能材料的热安全性研究.E-mail：1517603458@qq.com

通讯作者: 王昕捷（1991—），男，副教授，博士生导师，主要从事爆炸毁伤技术研究.E-mail：wangxinjie@bit.edu.cn

计量

出版历程

Cook-off Test and Numerical Simulation of HMX-Based Cast Explosive Containing AP

1. 实验与数值模拟方法

1.1 小型烤燃弹实验

1.2 烤燃数值模拟

2. 结果与讨论

2.1 烤燃实验

2.2 数值模拟

2.3 温升速率对点火位置的影响

3. 结 论

计量

出版历程

目录

1. 实验与数值模拟方法

1.1 小型烤燃弹实验

1.2 烤燃数值模拟

2. 结果与讨论

2.1 烤燃实验

2.2 数值模拟

2.3 温升速率对点火位置的影响

3. 结 论

作者简介:
刘润泽（1997—），男，硕士研究生，主要从事含能材料的热安全性研究.E-mail：1517603458@qq.com

通讯作者:
王昕捷（1991—），男，副教授，博士生导师，主要从事爆炸毁伤技术研究.E-mail：wangxinjie@bit.edu.cn

3. 结　论

3. 结　论