Ignition and burning mechanisms of RDX/HMX particlessubjected to drop-weight impact
-
摘要: 对单质炸药受低速撞击的力学和化学响应研究，是进行炸药敏感性评价的基础。利用配备了光学观测的落锤撞击装置，实现了频率为1.5×105 s-1的实时观测，不但可以区分样品的“爆”或“不爆”，而且可以获取RDX和HMX颗粒炸药受落锤低速撞击变形、破碎、溅射、点火和燃烧随时间演化的特征。结果表明:RDX颗粒是在液相中点火，而HMX颗粒在固相中点火; 燃烧反应前常常发生剧烈的溅射现象，溅射是由气相反应产物释放能量推动破碎的颗粒所致。对比了单个和单层颗粒炸药响应的特点，多个颗粒由于热点密集和破碎后相互作用，其燃烧反应比单个颗粒燃烧反应更剧烈。根据图像处理估算燃烧波传播速度，很好地表征了样品宏观燃烧反应的剧烈程度。Abstract: The mechanical and chemical response of explosives under low velocity impact is the basis for the evaluation of their sensitivity. The drop-weight impact apparatus equipped with the optical photography was used to achieve a frequency of 150 thousand seconds of real time observation. It was capable of distinguishing the samples of "Go" or "no Go" and observing such characteristics of deformation, crushing and breakage, jetting, ignition and combustion evolution, of RDX and HMX particles under low velocity impact. The selected photographic frames show that ignition occurred in the partially melted RDX phase. But for the HMX particles, ignition mainly occurred in the solid phase. A violent jetting phenomenon often occurs before the reaction of combustion. The occurrence of jetting primarily results from the energy released by gaseous products, which push the pulverized or melted explosives splash. The response characteristics of single and multiple granular explosives were compared. Because of the interaction of the hot spots broken, the combustion reaction of the particles is more intense than that of the single particles. The size ratio of the image to the actual length can be used to estimate combustion wave propagation velocity in each case, which is very suitable for characterizing the intensity of the macro-combustion reaction.
-
Key words:
- combustion /
- drop-weight impact /
- RDX/HMX /
- granular explosives /
- ignition /
- high-speed photography
-
振动测试过程中，由于周边环境及系统自身而导致的振动波形偏离基线中心的现象，称为信号趋势项。在深孔台阶爆破震动观测中，近区信号受趋势项的影响比远区更大。受趋势项影响，波形最大峰值将发生变化，频谱分析的精度降低，同时也会对信号的积分变换造成一定的影响。
近年来，有不少学者对趋势项的消除方法进行了研究[1-2]，其中，小波法和经验模态分解(empirical mode decomposition, EMD)法应用于非平稳信号处理的案例较多[3-4]。二者区别在于，小波法需针对信号特征事先选择小波基函数及分解层数，分解得到的低频系数可能与趋势项频带差距较大，分解过程对操作者的先验知识水平要求较高；而EMD法为自适应分解，对先验知识水平要求较低，但EMD法在分解过程中，会受到端点震荡和模态混叠效应的影响，最重要的是，目前尚无依据确定哪个固有模态函数(intrinsic mode function, IMF)分量为趋势项的主要组成部分。
基于此，本文中拟采用集合经验模态分解(ensemble empirical mode decomposition, EEMD)并结合小波阈值去噪方法，辅以人工判别及自相关分析手段，探索并建立适用于爆破近区的振动信号趋势项去除方法。
1. 趋势项去除方法
1.1 基于人工判别的EEMD趋势项去除方法
1.1.1 EEMD方法对EMD方法的改进[5]
EEMD方法较EMD方法的改进主要在于：在进行EMD分解时，加入高斯白噪声，原信号经多次添加噪声序列，并多次分解后，对同一尺度下各IMF分量求和，可得到各IMF分量的“集体”。由于加入的i次噪声是不相关随机序列，其统计均值为零，因此，对各IMF分量整体平均后，加入的i次噪声可近似相互抵消。
1.1.2 EEMD方法的分解步骤
(1) 添加标准白噪声信号wi(t)至原信号x(t)，二者等长。
xi(t)=x(t)+wi(t)i=1,2…,n (1) 式中：i为加入噪声的次数，xi(t)为原信号x(t)加入i次噪声后的信号。
(2) 对xi(t)进行EMD自适应分解，得到相应的各IMF分量，记为cij(t)，余项为ri(t)，cij(t)为第i次加入白噪声后分解得到的第j个IMF，j=1, 2…, k，其中k为分解尺度。
(3) 对各IMF分量进行总体平均，即得真实的IMF分量:
cj(t)=1nn∑i=1cij(t) (2) 1.1.3 人工判别
针对分解得到的各IMF分量，目前尚缺少明确的判别准则来识别哪一个分量为趋势项有效组成部分，趋势项的去除存在一定的盲目性。通过对深孔台阶爆破近区大量实测振动信号的分析，总结出趋势项产生的原因，结合现场振动测试系统的物理指标，可建立人工判别准则。
(1) 近区爆破振动的特点是瞬时输入能量大、持时短、振幅波动大。传感器在较大脉冲的作用下，受内部惯性原件和电路部分的输出限制，测量指标很难保持线性输出；
(2) 目前广泛使用的的爆破测振仪，其可测频率的最低值一般为5 Hz。现场观测得到的爆破信号，低于5 Hz的部分将很难保持线性输出，因此可视为低频干扰。
1.2 基于自相关分析的小波阈值去噪
利用EEMD方法对原始信号分解并消除趋势项分量后，引入的高斯白噪声将导致重构信号的信噪比明显下降。此时，可通过小波阈值法对保留的IMF分量进行去噪。为了有效甄别噪声分量，可借助自相关函数的特性，对保留IMF分量进行辅助判别，并对判别得到的具有噪声特征的分量进行去噪处理，最后，将所有IMF分量重构，整个去噪过程结束。
1.2.1 自相关分析
随机变量在时间域的相互依赖关系可以用自相关函数来表示，如下所示：
Rx(t1,t2)=E[x(t1)⋅x(t2)]=1N−1∑i[x(t1)⋅x(t2)] (3) 式中：x(t)代表某一随机过程，N为采样点个数，t1和t2为不同时刻。
对于高斯白噪声，其自相关函数在图像上通常具备如下特点：最大值集中在坐标零点，其余各点近似为零。
1.2.2 小波阈值去噪
小波阈值去噪是通过比较小波系数与设定阈值的大小，以此来判断小波系数主要是由信号还是噪声引起。去噪前首先需要确定阈值门限并选择合理的去噪方法。阈值门限的选取规则主要有固定阈值门限准则、无偏风险估计准则、混合阈值门限准则以及极大极小阈值门限准则；去噪方法主要有为软、硬阈值法。将去噪后的小波系数重构，即得到了去噪后的信号。
1.3 趋势项去除流程
根据上述分析，整理得到了深孔台阶爆破近区振动信号的趋势项去除方法，通过Matlab平台可编程对实测波形进行批量化处理，如图 1所示。
2. 实例分析
2.1 测试条件
选取现场某次垂向振动信号，利用前述方法对实测信号进行趋势项去除，并对重构信号进行评价。测试的条件见表 1。
表 1 信号测试条件Table 1. Conditions of the test signal测区 爆心距/m 最大段药量/kg 近区 65 2280 2.2 趋势项去除
2.2.1 EEMD方法分解
初始设置白噪声标准差为0.1，集成次数100。分解得到14个IMF分量及1个余量r，如图 2所示。
从图中不难看出，EEMD方法分解得到的各IMF分量，端点振荡效应及模态混叠效应已基本消除。
2.2.2 IMF分量尺度排序
利用FFT变换，可以得到各IMF分量的尺度信息。按照由大到小的顺序，各IMF分量的主频，如表 2所示。由于趋势项是低频干扰，为简化显示范围，可略去高频分量IMF1~IMF4，仅在频谱曲线图中显示IMF5~IMF15分量，如图 3所示。
表 2 各IMF分量主频Table 2. Dominant frequency of each IMF componentIMF分量 主频/Hz IMF1 33.40 IMF2 14.40 IMF3 39.20 IMF4 32.40 IMF5 36.40 IMF6 19.40 IMF7 10.40 IMF8 7.60 IMF9 1.60 IMF10 1.00 IMF11 0.80 IMF12 0.60 IMF13 0.40 IMF14 0.20 r 0.00 2.2.3 人工干预判别
由图 5可以看出，IMF5~IMF8频段较高，频带较宽；IMF9~IMF14频段较低，频带较窄。由表 2中各IMF分量的主频值可知，IMF9~IMF14分量主要集中在0~5 Hz频带，已明显超出爆破测振仪的有效监测范围，上述分量将引起信号趋势项的产生。因此，需将IMF9~IMF14分量去除。
2.2.4 趋势项去除后剩余分量重构
对保留的IMF1~IMF8分量重构，如图 4~5所示。从图中可以看出，波形重新回到基线中心，主频也由低频缓慢过渡到合理区间，主频幅值突高的情况消失。
2.2.5 自相关分析
对保留的IMF1~IMF8分量进行自相关分析，如图 6所示。由图中可以看出，8个IMF分量中，IMF1~IMF3分量的自相关函数符合高斯白噪声的特征，因此，主要对IMF1~IMF3分量进行小波阈值去噪。
2.2.6 小波阈值去噪
阈值函数确定方法较多[6]。本文中阈值门限准则采用heursure函数，去噪方法采用软阈值函数。
2.2.7 去噪后剩余分量重构
IMF1~IMF3分量去噪结束后，将其与其余IMF分量一并进行重构，即可得到最终爆破振动信号，如图 7所示。从图中可以看出，去噪后波形的信噪比得到一定提升，滤波效果显著。
3. 结论
(1) 深孔台阶爆破条件下，爆破近区振动信号产生趋势项的原因为瞬时大能量输入引起的非线性失真，在此基础上以测试仪器有效监测范围作为识别趋势项组成部分的判别准则。
(2) 借助EEMD方法及小波分解，提出一种基于各IMF分量频带指标，辅以人工判别的趋势项去除方法，以及利用自相关分析识别噪声特征的小波阈值去噪方法，并通过实例验证了方法的有效性。
-
图 1 扫描电镜观测的不同放大倍数下的RDX和HMX颗粒样品
Figure 1. Microscopic morphologies of RDX and HMX particles observed by scanning electron microscopy
图 2 落锤加载装置及光路观测示意图
Figure 2. Schematic diagram showing the drop-weight loadingand optical measurement system
图 4 单个RDX颗粒受撞击的响应过程，落锤下落高度为15 cm
Figure 4. elected photographic frames for an individual RDX particlein response to a drop-weight impact with 15 cm height
图 6 单层(10颗)RDX颗粒受撞击后的点火燃烧图像，落锤下落高度为15 cm
Figure 6. Selected photographic frames for ten RDX particles in response to a drop-weight impact of 15 cm height
图 8 单个HMX颗粒受撞击的响应过程，落锤下落高度为15 cm
Figure 8. Selected photographic frames for an individual HMX particle in response to a drop-weight impact with 15 cm
图 10 单层(10颗)HMX颗粒受撞击后的点火燃烧图像，落锤下落高度为15 cm
Figure 10. Selected photographic frames for ten HMX particles in response to a drop-weight impact with 15 cm
表 1 火焰传播速度
Table 1. Propagation velocity of combustion flame
实验方案 落锤高度/cm 火焰速度/(m·s-1) 单个RDX颗粒 15 173.9 单层多个RDX颗粒 15 235.4 单个HMX颗粒 15 93.0 单层多个HMX颗粒 15 308.9 
下载: 导出CSV
-
[1] Balzer J E, Proud W G, Walley S M.et al.High-speed photographic study of the drop-weight impact response of RDX-DOS mixtures[J].Combustion and Flame, 2003, 135(4):547-555. doi: 10.1016/j.combustflame.2003.08.009 [2] Buntain G A, Mckinney T, Rivera T, et al.Decomposition of energetic materials on the drop-weight-impact machine[C]//Proceedings of the 9th Symposium (International) on Detonation.Portland, Oregon, United States, 1989: 1037-1043. [3] Rice B M, Hare J J.A quantum mechanical investigation of the relation between impact sensitivity and the charge distribution in energetic molecules[J].Journal of Physical Chemistry A, 2002, 106(9):1770-1783. doi: 10.1021/jp012602q [4] Keshavarz M, Jaafari M.Investigation of the various structure parameters for predicting impact sensitivity of energetic molecules via artificial neural network[J].Propellants, Explosives, Pyrotechnics, 2006, 31(3):216-225. doi: 10.1002/(ISSN)1521-4087 [5] Swallowe G M, Field J E.The ignition of a thin layer of explosive by impact:The effect of polymer particles[J].Proceedings of the Royal Society of London A:Mathematical and Physical Sciences, 1982, 379(1777):389-408. doi: 10.1098/rspa.1982.0022 [6] Swallowe G M, Field J E, Hutchinson C D.Impact experiments on thin layers of polymers and intermediate explosives[C]//Gupta Y M.Proceedings of the Fourth American Physical Society Topical Conference on Shock Waves in Condensed Matter.New York: Plenum Press, 1985: 891-898. doi: 10.1007/978-1-4613-2207-8_131 [7] Heavens S N, Field J E.The ignition of a thin layer of explosive by impact[J].Proceedings of the Royal Society of London A:Mathematical and Physical Sciences, 1982, 379(1612):77-93. http://www.jstor.org/stable/78552 [8] Field J E, Swallowe G M, Heavens S N.Ignition mechanisms of explosives during mechanical deformation[J].Proceedings of the Royal Society of London A:Mathematical and Physical Sciences, 1982, 382(1782):231-244. doi: 10.1098/rspa.1982.0099 [9] Field J E, Palmer S J P, Pope P H, et al.Mechanical properties of PBXs and their behaviour during drop-weight impact[C]//Short J M.8th Symposium (International) on Detonation.White Oak, Maryland, USA: Naval Surface Weapons Center, 1985: 635-644. [10] Balzer J E, Field J E, Gifford M J, et al.High-speed photographic study of the drop-weight impact response of ultrafine and conventional PETN and RDX[J].Combustion and Flame, 2002, 130(4):298-306. doi: 10.1016/S0010-2180(02)00373-5 [11] Hamdan S, Swallowe G M.The strain-rate and temperature dependence of the mechanical properties of polyetherketone and polyetheretherketone[J].Journal of Materials Science, 1996, 31(6):1415-1423. doi: 10.1007/BF00357847 -
下载:
计量
- 文章访问数: 4482
- HTML全文浏览量: 1353
- PDF下载量: 405
- 被引次数: 0