煤尘层最低着火温度变化规律实验研究

刘天奇; 李雨成; 罗红波

doi:10.11858/gywlxb.20170649

煤尘层最低着火温度变化规律实验研究

doi: 10.11858/gywlxb.20170649

刘天奇^{1, 2, 3,},
李雨成^{1, 2, 3},
罗红波^{1, 2, 3}

1.
辽宁工程技术大学安全科学与工程学院, 辽宁阜新 123000
2.
辽宁工程技术大学矿山热动力灾害与防治教育部重点实验室, 辽宁阜新 123000
3.
辽宁工程技术大学煤炭资源安全开采与洁净利用工程研究中心, 辽宁阜新 123000

基金项目:

国家自然科学基金 51774168

国家自然科学基金 51204089

辽宁省自然科学基金 201602355

辽宁省教育厅项目 LJYL002

辽宁省教育厅项目 551610001276

详细信息

作者简介:
刘天奇(1990-), 男, 博士研究生, 主要从事矿井通风与粉尘放置理论及技术研究.E-mail:ltq613@163.com

中图分类号: X932
计量
- 文章访问数: 7061
- HTML全文浏览量: 2936
- PDF下载量: 482
出版历程
- 收稿日期: 2017-09-27
- 修回日期: 2017-10-30

Variation Law for Minimum Ignition Temperature of Coal Dust Layer

LIU Tianqi^{1, 2, 3
,},
LI Yucheng^{1, 2, 3},
LUO Hongbo^{1, 2, 3}

1.
College of Safety Science and Engineering, Liaoning Technical University, Fuxin 123000, China
2.
Key Laboratory of Mine Thermodynamic Disasters and Control of Ministry of Education, Liaoning Technical University, Fuxin 123000, China
3.
Research Center of Coal Resources Safe Mining and Clean Utilization, Liaoning Technical University, Fuxin 123000, China

摘要

摘要: 为研究煤尘层最低着火温度随煤样变质程度、煤尘粒径及煤尘层厚度的变化规律，采用煤尘层最低着火温度测定系统进行实验研究，结果表明：随着煤样由褐煤到无烟煤变质程度逐渐增大，煤尘层最低着火温度由290℃上升到400℃以上，同时，褐煤、长焰煤、不粘煤、气煤煤尘层着火时观察到明显的火焰。随着煤尘粒径不断减小，不同煤质的煤尘层最低着火温度明显减小，煤尘层厚度为15 mm时，随着煤尘粒径由0.5 mm减小至0.075 mm，不同煤质煤尘层最低着火温度分别减小了31.0%、26.7%、28.1%、25.8%、28.6%、27.8%、18.9%和15.0%，煤尘粒径影响作用十分显著。随着煤尘层厚度的增大，不同煤质在不同粒径下的煤尘层最低着火温度都减小，其中无烟煤的变化最不明显。
- 煤尘层厚度 /
- 着火温度 /
- 变化规律 /
- 煤尘粒径
Abstract: In this paper we studied the variation of the minimum ignition temperature of coal dust layer with its metamorphism, particle size and thickness using the minimum ignition temperature measurement system.The results showed that, with the gradual increase of its degree of metamorphism, the minimum ignition temperature of coal dust layer varied from 290℃ to above 400℃; and that when the coal dust layer of lignite, long flame coal, non-stick coal, gas coal were ignited, an obvious flame was observed; and that, as the coal dust particle size decreased, the minimum ignition temperature of different kinds of the coal dust layer exhibited a significant reduction trend.It was found that when the thickness of the layer is 15 mm, with the particle size decreasing from 0.5 mm to 0.075 mm, the minimum ignition temperature decreased by 31.0%, 26.7%, 28.1%, 25.8%, 28.6%, 27.8%, 18.9% and 15.0%, respectively, indicating a very obvious effect of the coal dust particle size.With the increase of the layer's thickness, the minimum ignition temperature of different kinds of the coal dust layer exhibited a decreasing trend, while the trend of anthracite was the least significant.
- thickness of coal dust layer /
- ignition temperature /
- change law /
- coal dust particle size

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. Coal dust layer MIT tester

下载: 全尺寸图片幻灯片

图 2 煤尘层最低着火温度测定装置结构

Figure 2. Structure of coal dust layer MIT tester

1.Thermocouple height adjustment knob; 2.Spring; 3.Heater base; 4.Recording thermocouple on hot plate surface; 5.Hot plate surface; 6.Heater; 7.Metal ring; 8.Control thermocouple on hot plate surface; 9.Dust layer thermocouple; 10.Heater leads

下载: 全尺寸图片幻灯片

图 3 热板表面上煤尘层着火的典型温度-时间曲线

Figure 3. T-t curve of coal dust layer on hot surface

下载: 全尺寸图片幻灯片

图 4 煤尘层温度-时间曲线

Figure 4. T-t curve of coal dust layer

下载: 全尺寸图片幻灯片

图 5 煤尘层最低着火温度与粒径关系

Figure 5. T_L-r curve of coal dust layer

下载: 全尺寸图片幻灯片

表 1 8种煤质的煤尘层最低着火温度

Table 1. Coal dust layer MIT of 8 kinds of coal quality

No.	Coal quality	Particle size/mm	Thickness/mm	Minimum ignition temperature/℃	Ignition type	Ignition time/min
1	Lignite	0.2	5	290	a	18
2	Long flame coal	0.2	5	310	a	22
3	Not sticky coal	0.2	5	330	a	18
4	Gas coal	0.2	5	320	a	20
5	Coking coal	0.2	5	340	c	11
6	Lean coal	0.2	5	360	c	9
7	Meagre coal	0.2	5	370	c	13
8	Anthracite	0.2	5	>400	No fire	>30

下载: 导出CSV

参考文献(12)

[1]	金龙哲.矿井粉尘防治理论[M].北京:科学出版社, 2010:22-30. JIN L Z.Dust prevention theory in mine[M].Beijing:Science Press, 2010:22-30.
[2]	毕明树.气体和粉尘爆炸防治工程学[M].北京:化学工业出版社, 2012:39-48. BI M S.Gas and dust explosion prevention engineering sciences[M].Beijing:Chemical Industry Press, 2012:39-48.
[3]	毕明树, 李江波.密闭管内甲烷-煤粉复合爆炸火焰传播规律的实验研究[J].煤炭学报, 2010, 35(8):1298-1302. http://www.oalib.com/paper/4235181 BI M S, LI J B.Experimental study on flame propagation law of methane-pulverized coal compound explosion in closed tube[J].Journal of China Coal Society, 2010, 35(8):1298-1302. http://www.oalib.com/paper/4235181
[4]	宫广东, 刘庆明, 胡永利.管道中煤尘爆炸特性实验[J].煤炭学报, 2010, 35(4):609-612. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=mtxb201004020&dbname=CJFD&dbcode=CJFQ GONG G D, LIU Q M, HU Y L.Experimental study on explosion characteristics of coal dust in pipeline[J].Journal of China Coal Society, 2010, 35(4):609-612. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=mtxb201004020&dbname=CJFD&dbcode=CJFQ
[5]	牛芳, 刘庆明, 白春华.甲烷-煤尘爆炸物火焰传播特性[J].高压物理学报, 2012, 26(4):455-461. doi: 10.11858/gywlxb.2012.04.015 NIU F, LIU Q M, BAI C H.Flame propagation characteristics of methane-coal dust explosion[J].Chinese Journal of High Pressure Physics, 2012, 26(4):455-461. doi: 10.11858/gywlxb.2012.04.015
[6]	李雨成.矿井粉尘防治理论及技术[M].北京:煤炭工业出版社, 2015:35-56. LI Y C.Dust prevention theory and technology in mine[M].Beijing:Coal Industry Press, 2015:35-56.
[7]	陈金健. 煤尘云最低着火温度及抑制技术研究[D]. 太原: 中北大学, 2016. http://cdmd.cnki.com.cn/Article/CDMD-10110-1016179809.htm CHEN J J. Study on lowest ignition temperature and suppression technology of coal dust cloud[D]. Taiyuan: North University of China, 2016. http://cdmd.cnki.com.cn/Article/CDMD-10110-1016179809.htm
[8]	李润之.不同挥发分煤尘层最低着火温度变化规律研究[J].安全与环境学报, 2017, 17(3):954-957. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=aqhj201703029&dbname=CJFD&dbcode=CJFQ LI R Z.Study on variation law of minimum ignition temperature of different volatile coal dust[J].Journal of Safety and Environment, 2017, 17(3):954-957. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=aqhj201703029&dbname=CJFD&dbcode=CJFQ
[9]	司荣军, 王春秋.瓦斯对煤尘爆炸特性影响的实验研究[J].中国安全科学学报, 2006, 16(12):86-91. doi: 10.3969/j.issn.1003-3033.2006.12.016 SI R J, WANG C Q.Experimental study on the effect of gas on coal dust explosion characteristics[J].China Safety Science Journal, 2006, 16(12):86-91. doi: 10.3969/j.issn.1003-3033.2006.12.016
[10]	屈姣. 甲烷和煤尘爆炸特性实验研究[D]. 西安: 西安科技大学, 2015. http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=D719360 QU J. Study on explosion characteristics of gas and coal dust[D]. Xi'an: Xi'an Technology University, 2015. http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=D719360
[11]	宋广朋. 瓦斯煤尘共存的爆炸特性与传播研究[D]. 青岛: 山东科技大学, 2011. http://cdmd.cnki.com.cn/Article/CDMD-10424-1012276655.htm SONG G P. Study on ignition mechanism and experimental study of coal powder gas coupling system[D]. Qingdao: Shandong University of Science and Technology, 2011. http://cdmd.cnki.com.cn/Article/CDMD-10424-1012276655.htm
[12]	平洋. 煤粉瓦斯耦合体系着火机理和实验研究[D]. 沈阳: 东北大学, 2011. http://cdmd.cnki.com.cn/Article/CDMD-10145-1013116384.htm PING Y. Study on ignition mechanism and experimental study of coal powder gas coupling system[D]. Shenyang: Northeastern University, 2011. http://cdmd.cnki.com.cn/Article/CDMD-10145-1013116384.htm