Explosion severity parameters of hybrid mixture of methane and lycopodium dust
-
摘要: 基于改进的20 L球形粉尘爆炸装置,在相同初始条件下分别测量了甲烷、石松子粉尘和甲烷/石松子两相混合体系的爆炸压力、爆炸压力上升速率和爆炸指数等参数,系统研究了甲烷/石松子粉尘两相混合体系爆炸特性变化规律。结果表明:甲烷的添加能显著提高低质量浓度石松子粉尘爆炸压力而降低高质量浓度石松子粉尘爆炸压力;甲烷对石松子粉尘最大爆炸压力没有显著影响,但能显著提高石松子粉尘最大爆炸压力上升速率。甲烷/石松子粉尘混合体系爆炸指数高于单相石松子粉尘爆炸指数,但甲烷/石松子粉尘混合体系和单相石松子粉尘爆炸指数均低于单相甲烷爆炸指数。工业生产过程中应避免粉尘混入可燃气体以降低粉尘爆炸危险性。Abstract: A modified 20 L standard spherical dust explosion vessel was used to systematically study the explosion characteristics of the methane/lycopodium hybrid mixtures. The explosion pressure (pex), the explosion pressure rise rate (dp/dt)ex and the explosive deflagration index (Kst) of the single-phase methane, the single-phase lycopodium dust, and the methane/lycopodium hybrid mixtures were measured under the same initial conditions. The results showed that methane could obviously enhance the explosion pressure pex of low-concentration lycopodium dust but reduce the pex of high-concentration lycopodium dust. It was proved that methane had no significant effects on the maximum explosion pressure pmax of lycopodium dust. But it could significantly increase the maximum explosion pressure rise rate (dp/dt)max. By evaluating the Kst, it was found that the explosive deflagration index of the methane/lycopodium hybrid mixtures was higher than that of single-phase lycopodium dust. However, the explosive deflagration indices of both methane/lycopodium hybrid mixtures and single-phase lycopodium dusts were lower than that of single-phase methane. Therefore, the coexistence of combustible gas and combustible dust in industrial production process should be avoided to reduce the risk of dust explosions.
-
Key words:
- dust explosion /
- hybrid mixtures /
- explosion severity /
- methane /
- lycopodium
-
E.H.Lee等[1]首次利用理想刚塑性模型,对自由梁中点受三角形分布脉冲载荷作用下的动力响应进行了分析,得出初始输入能量在响应各阶段的分配。为模拟航天飞机的着陆过程,N.Jones等[2]研究了均匀截面或阶梯截面自由梁受到三角形分布冲击载荷作用下的变形及失效,指出空中自由梁在载荷相同条件下比受约束的梁更难发生失效。J.L.Yang等[3]进行了自由梁任意截面受冲击的实验和理论研究, 用高速摄影测量自由梁任意截面受弹丸冲击的初始变形,试件材料为2ALY12铝,冲击速度为20~200 m/s,并推导了自由梁求解真实加速度场的运动方程和一般离散表达式。
目前,对实心自由梁的研究成果多局限于弹速为200~300 m/s的低速碰撞,且重点集中在实心截面梁整体变形效应的动态力学响应研究方面[4],该问题在航空、航天工程中具有重要的应用背景[5-9]。而对于弹体以高速碰撞圆柱壳自由梁产生的局部破坏模式、穿孔效应等情况,研究报导甚少。
本文中,基于二级轻气炮加载弹丸高速碰撞圆柱壳自由梁实验,并利用ANSYS/LS-DYNA软件中拉格朗日法,对高速撞击圆柱壳自由梁过程进行数值模拟,通过实验和数值模拟结果的对比和量纲分析,得出弹丸高速碰撞圆柱壳自由梁产生的穿孔规律以及影响穿孔直径的主要因素。
1. 实验
1.1 实验参数
实验采用二级轻气炮加载弹丸,二级轻气炮一级驱动气体为氮气,充气压力为2~40 MPa;二级驱动气体为氢气,充气压力为0.1~1 MPa;最高发射速度可达7 km/s。速度测量采用磁感应测速法,靶舱内压力小于100 Pa。
实验时,调整好靶舱内薄壁圆柱壳的位置和高度,使铝弹丸垂直撞击圆柱壳长度方向的中点。实验布置如图 1所示。
圆柱壳自由梁的材料为304无缝不锈钢管,壁厚均为3 mm,圆柱壳的直径为30~60 mm;弹丸为2A12球形铝合金,直径为4.6 mm,质量为0.14 g;实验1~9的碰撞速度分别为2.10、2.06、2.10、2.06、1.80、2.30、2.58、2.71和3.30 km/s。
1.2 实验分组
实验分为2组,考察穿孔直径的大小:(1)相近碰撞速度,分别碰撞4种不同直径(分别为30、40、48和60 mm)的圆柱壳自由梁;(2)不同碰撞速度,分别碰撞直径为60 mm的圆柱壳自由梁。图 2为4种不同直径圆柱壳自由梁试件。
1.3 实验结果
弹丸高速碰撞圆柱壳自由梁过程中,产生的压缩波和反射稀疏波的共同作用,使弹丸和圆柱壳材料被部分抛出,从而在圆柱壳自由梁上形成穿孔。圆柱壳自由梁上均出现了近似椭圆形的穿孔,穿孔边缘整齐,无裂纹和撕裂现象。穿孔的轴向直径比径向直径稍大些,相差约1 mm。穿孔形状大致为圆锥台,外壁孔径较大,内壁孔径较小,相差3 mm左右。弹丸未对圆柱壳造成完全贯穿,在远离撞击端的一侧出现一个很小的弹坑。图 3分别给出了4种不同直径圆柱壳自由梁的穿孔效果。
2. 数值模拟
利用ANSYS/LS-DYNA软件,对高速碰撞圆柱壳自由梁进行数值模拟。弹丸和圆柱壳均采用黏塑性理论的Johnson-Cook材料本构模型和Mie-Grüneison状态方程。弹丸材料为2A12铝,圆柱壳材料为304钢,采用Lagrange法建模。Lagrange法优势在于,变形后材料的自由表面能自动被网格边界捕捉到,在网格中不存在材料的流动,纯拉格朗日法可很好地分析各种中等变形程度的问题, 具有计算精度高、效率高,界面清晰等优势[6]。圆柱壳中间受撞击部分采用细化网格的方法,以0.25 mm为单位划分网格,能较好地模拟撞击点的穿孔情况;圆柱壳自由梁两端大部分未变形或小变形区网格进行粗化,以2 mm为单位划分网格。
2.1 材料模型和结果
高速撞击过程中,弹、靶接触区域发生剧烈的塑性变形和局部穿孔破坏。数值计算中,弹丸和自由圆柱壳均采用经典的基于黏塑性理论的Johnson-Cook材料模型,模拟材料的应变硬化、应变率硬化和热软化效应。流动屈服应力可表示为:
σy=(A+Bˉεnp)(1+Cln˙ε∗)(1−T∗m) (1) 式中:σy为材料流动屈服强度,A、B、n、C、m为材料相关常数,εp为有效塑性应变,˙ε∗=˙ˉεp/˙ε0, 表示为量纲一有效塑性应变率,取参考应变率˙ε0=1 s−1。T*=(T-Tr)/(T-Tm),其中Tr和Tm分别对应为室温和材料熔化温度。破坏应变定义如下:
εf=(D1+D2exp(D3σ∗))(1+D4ln˙ε∗)(1+D5T∗) (2) 式中:σ*为压力与有效压力之比。当破坏参数D=∑Δˉεp/εf=1时,就会发生破裂。
另外,高速撞击压力远超过材料的动态屈服强度,接触区域材料的塑性变形如流体流动,可采用Mie-Grüneisen状态方程描述材料静水压力与内能关系:
p−pH=Γρ(E−EH) (3) pH=ρ0c0μ(1+μ)(1−(s−1)μ)2 (4) EH=pH2ρ0(μ1+μ) (5) Grüneisen状态方程基于Hugoniot冲击曲线,可根据冲击波和波后粒子速度的简单线性关系描述材料的冲击特性,u=c0+sup。假设Γρ=Γ0ρ0为材料常数,Γ0为Grüneisen系数,ρ、ρ0分别为材料当前和初始状态密度。设μ=ρ/ρ0-1,c0为材料体积声速,s为材料参数,由冲击实验确定。
计算模型中,铝球弹丸和钢管的部分材料参数见表 1。数值模拟结果见表 2,表中,D为圆柱壳直径,H为壁厚,v为碰撞速度,da、dr分别为轴向、径向穿孔直径。
表 1 材料参数Table 1. Material parameters材料 ρ/(kg·m-3) G/GPa A/MPa B/MPa n C m s c0/(km·s-1) γ0 2A12铝 2 780 27.6 265 426 0.34 0.015 1.0 1.34 5.330 2.0 304钢 7 830 77.0 332 245 0.36 0.022 1.0 1.92 3.570 1.67 表 2 数值模拟结果Table 2. Simulation resultsNo. D/mm H/mm v/(km·s-1) da/mm dr/mm 1 30 3 2.2 6.5 6.0 3 30 3 3.2 9.4 8.4 4 40 3 2.2 6.4 6.0 5 40 3 3.2 8.5 8.4 6 40 2 2.7 7.7 7.3 7 48 3 2.2 6.3 5.8 9 48 3 3.2 8.5 8.0 10 60 3 2.2 6.0 5.8 12 60 3 3.2 8.1 7.5 2.2 结果对比
实验和数值模拟采用的均为内壁孔径。图 4为典型工况下实验和数值模拟穿孔直径对比结果。表 3给出了穿孔直径大小,可知实验结果比数值模拟结果高出量最大为21.5%,最小2.53%,平均高出量为15.22%。因使用了单元侵蚀,数值模拟材料选取的是模拟软件材料库中铝和钢的标准性能参数,与实验所用2A12铝和304不锈钢的实际性能参数略有差别。因此实验结果和数值模拟结果也存在微小差别,由表 3可看出两者有差异是合理的。如图 5所示,数值模拟和实验结果拟合的曲线变化规律基本相同,因此可以利用数值模拟对圆柱壳自由梁高速撞击产生的穿孔规律进行分析和预测。
表 3 实验与数值模拟结果的比较Table 3. Comparison of simulation and experiment resultsD/mm v/(km·s-1) da/mm dr/mm 实验 模拟 实验 模拟 30 2.10 7.9 6.2 7.4 6.0 40 2.06 7.5 6.1 7.4 5.9 48 2.10 7.3 6.1 6.7 5.9 60 2.06 7.2 6.0 6.5 5.8 60 1.80 7.2 6.0 6.5 5.0 60 2.30 7.2 6.2 7.0 5.9 60 2.58 8.4 6.8 7.5 6.5 60 2.71 8.4 7.0 7.8 6.5 60 3.30 8.9 8.6 7.9 7.7 2.3 穿孔直径的变化规律
图 6为高速撞击圆柱壳自由梁产生的穿孔直径随碰撞速度变化的规律,当圆柱壳直径不变时,穿孔直径随碰撞速度的增大而增大,但穿孔直径增大的幅度随着碰撞速度增大而减小。
图 7为高速撞击圆柱壳自由梁产生的轴向穿孔直径随圆柱壳直径变化的规律,当碰撞速度和圆柱壳壁厚不变时,穿孔直径随圆柱壳直径的增大而减小,但减小幅度随圆柱壳自由梁直径的增大而减小。
图 8为高速撞击圆柱壳自由梁产生的轴向和径向穿孔直径随圆柱壳自由梁壁厚变化的规律,对于直径分别为40和60 mm的圆柱壳自由梁,当碰撞速度为2.7 km/s时,高速撞击圆柱壳自由梁产生的穿孔直径随圆柱壳自由梁壁厚的增加而减小,在圆柱壳壁厚大于4.5 mm后,圆柱壳将不再被穿透。
3. 经验关系式
为确定穿孔直径与圆柱壳自由梁的几何、物理参数及碰撞速度之间的关系,考察弹丸高速撞击圆柱壳自由梁的情况。穿孔直径可表示为:
d=f(dp,ρp,Ep,σyp,μp,εfp,lt,D,Di,ρt,Et,σyt,μt,εft,v) (6) 式中:dp为弹丸直径,ρp为弹丸密度,Ep为弹丸的杨氏模量,σyp为弹丸材料的屈服极限,μp为弹丸材料的泊松比,εfp为弹丸材料的屈服应变;D为圆柱壳自由梁的外径,Di为圆柱壳自由梁的内径,lt为圆柱壳自由梁的长度,ρt为圆柱壳自由梁的材料密度,Et为圆柱壳自由梁材料的弹性模量,σyt为圆柱壳自由梁材料的屈服极限,μt为圆柱壳自由梁材料的泊松比,εft为圆柱壳自由梁材料的屈服应变。
根据量纲分析,从所列的16个参数中(其中包括穿孔直径d)可以组成12个彼此独立的量纲一量。取其中ρp、σyp和dp作为基本量,可得到量纲一函数关系:
ddp=f(v√σyp/ρp,Epσyp,ltdp,Ddp,Didp,ρtρp,∂tρp,σytσyp,Etσyp,δtdp,μp,εfp,μt,εft) (7) 如果模型实验采用与原型相同种类的弹丸和圆柱壳自由梁材料,圆柱壳自由梁的长度和厚度以及弹丸的直径都不变。那么,上式可简化为:
ddp=f(v√σyp/ρp,Ddp,Didp) (8) 将这个隐式方程式写成显示方程式[7], 则:
ddp=x0(v√σyp/ρp)x1(Ddp)x2(Didp)x3 (9) 式中:xi(i=0, 1, 2, 3)均是待定系数。
利用最小二乘法,对数值模拟数据进行曲线拟合,确定公式中的待定系数,可得到轴向和径向穿孔直径的经验关系式:
dadp=2.5229(v√σyp/ρp)0.6741(Ddp)−1.4309(Didp)1.5079 (10) drdp=2.5454(v√σyp/ρp)0.8076(Ddp)−1.8079(Didp)1.8385 (11) 通过对经验关系式计算得到的穿孔直径与数值模拟得到的穿孔直径的比较,表明计算结果和数值模拟结果吻合得很好。
4. 结论
通过高速撞击圆柱壳自由梁的穿孔实验和数值模拟,并利用数值模拟和实验结果的拟合曲线,可以得出以下结论。
(1) 高速撞击圆柱壳自由梁产生的穿孔直径在径向和轴向都随碰撞速度增大而增大,二者近似成线性关系,但穿孔直径增大的幅度随着碰撞速度的增大而减小。对于直径分别为40、48和60 mm的圆柱壳自由梁:当撞击速度为2~2.5 km/s时,撞击速度对撞击产生的穿孔直径影响较大;当碰撞速度超过2.5 km/s后,弹丸撞击速度对穿孔直径的影响不再明显,穿孔直径随碰撞速度的增加逐渐接近定值。
(2) 当弹丸碰撞速度和圆柱壳壁厚不变时,撞击产生的穿孔直径随圆柱壳自由梁直径增大而减小,但减小的幅度随圆柱壳直径的增大而减小。当圆柱壳直径为30~40 mm时,圆柱壳直径对穿孔直径影响较大,当圆柱壳直径大于40 mm时,圆柱壳直径对穿孔直径的影响将不再明显,穿孔直径随圆柱壳自由梁直径的增加将接近定值。
(3) 圆柱壳自由梁的厚度对高速撞击产生的穿孔直径有很大的影响,当碰撞速度为2.7 km/s时,高速撞击圆柱壳自由梁产生的穿孔直径随圆柱壳自由梁壁厚的增加而减小,在圆柱壳壁厚大于4.5 mm后,圆柱壳不会被穿透。
-
表 1 甲烷、石松子粉尘及甲烷/石松子混合体系爆燃指数
Table 1. Explosive deflagration indices of methane, lycopodium and methane/lycopodium mixtures
φ/% ρ/(g·m-3) (dp/dt)max/(MPa·s-1) Kst(KG)/(MPa·m·s-1) 0(pure dust) 750 18.81 5.11 2 750 19.98 5.42 4 500 21.16 5.74 6 250 28.21 7.66 10(pure CH4) 0 126.94 34.46 -
[1] SANCHIRICO R, RUSSO P, SALIVA A, et al. Explosion of lycopodium-nicotinic acid-methane complex hybrid mixtures[J]. Journal of Loss Prevention in the Process Industries, 2014, 36:505-508. https://www.sciencedirect.com/science/article/pii/S0950423014002241 [2] ADDAI E K, GABEL D, KRAUSE U. Experimental investigation on the minimum ignition temperature of hybrid mixtures of dusts and gases or solvents[J]. Journal of Hazardous Materials, 2016, 301:314-326. doi: 10.1016/j.jhazmat.2015.09.006 [3] ABUSWER M, AMYOTTE P, KHAN F. A quantitative risk management framework for dust and hybrid mixture explosions[J]. Journal of Loss Prevention in the Process Industries, 2013, 26(2):283-289. doi: 10.1016/j.jlp.2011.08.010 [4] CASHDOLLAR K L. Coal dust explosibility[J]. Journal of Loss Prevention in the Process Industries, 1996, 9(1):65-76. doi: 10.1016/0950-4230(95)00050-X [5] DUFAUD O, PERRIN L, TRAORE M, et al. Explosions of vapour/dust hybrid mixtures: A particular class[J]. Powder Technology, 2009, 190(1/2):269-273. https://www.researchgate.net/publication/232374600_Explosions_of_vapourdust_hybrid_mixtures_A_particular_class [6] SANCHIRICO R, BENEDETTO A D, GARCIA-AGREDA A, et al. Study of the severity of hybrid mixture explosions and comparison to pure dust-air and vapour-air explosions[J]. Journal of Loss Prevention in the Process Industries, 2011, 24(5):648-655. doi: 10.1016/j.jlp.2011.05.005 [7] BARTKNECHT W. Dust explosions: Course, prevention, protection[M]. Springer Science & Business Media, 2012:25-56. [8] DENKEVITS A, HOESS B. Hybrid H 2 /Al dust explosions in Siwek sphere[J]. Journal of Loss Prevention in the Process Industries, 2015, 36:509-521. doi: 10.1016/j.jlp.2015.03.024 [9] AMYOTTE P, ABUSWER M, DI BENEDETTO A, et al. Determination of hybrid mixture explosion severity[C]//Proceedings of the 13th International Symposium on Loss Prevention and Safety Promotion in the Process Industries. Brugge, Be, 2010, 2: 149-156. [10] AMYOTTE P, LINDSAY M, DOMARATZKI R, et al. Prevention and mitigation of dust and hybrid mixture explosions[J]. Process Safety Progress, 2010, 29(1):17-21. https://www.researchgate.net/publication/230036461_Prevention_and_Mitigation_of_Dust_and_Hybrid_Mixture_Explosions [11] LIU Yi, SUN Jinhua, CHEN Dongliang. Flame propagation in hybrid mixture of coal and methane[J]. Journal of Loss Prevention in the Process Industries, 2007, 20(4/5/6):691-697. https://www.sciencedirect.com/science/article/pii/S0950423007000629 [12] 喻健良, 闫兴清.高静态动作压力下粉尘爆炸泄放标准的可靠性[J].东北大学学报(自然科学版), 2015, 36(9):1316-1320. doi: 10.3969/j.issn.1005-3026.2015.09.022YU Jianliang, YAN Xingqing. Reliability of dust explosion venting standards under elevated static activation overpressures[J]. Journal of Northeastern University (Natural Science), 2015, 36(9):1316-1320. doi: 10.3969/j.issn.1005-3026.2015.09.022 [13] YAN Xingqing, YU Jianliang, GAO Wei. Duct-venting of dust explosions in a 20 L sphere at elevated static activation overpressures[J]. Journal of Loss Prevention in the Process Industries, 2014, 32(1):63-69. https://www.sciencedirect.com/science/article/pii/S0950423014001259 [14] GOING J E, CHATRATHI K, CASHDOLLAR K L. Flammability limit measurements for dusts in 20-L and 1-m3 vessels[J]. Journal of Loss Prevention in the Process Industries, 2000, 13(3/4/5):209-219. https://www.cdc.gov/niosh/mining/UserFiles/works/pdfs/flmfd.pdf [15] DASTIDAR A G, AMYOTTE P R. Explosibility boundaries for fly ash/pulverized fuel mixtures[J]. Journal of Hazardous Materials, 2002, 92(2):115-26. doi: 10.1016/S0304-3894(01)00390-9 [16] GARCIA-AGREDA A, DI BENEDETTO A, RUSSO P, et al. Dust/gas mixtures explosion regimes[J]. Powder Technology, 2011, 205(1):81-86. https://www.sciencedirect.com/science/article/pii/S0032591010004705 [17] National Fire Protection Association, NFPA68. Standard on explosion protection by deflagration venting[S]. Quincy, MA: Batterymarch Parck, 2007: 10-11. [18] E 1226. Standard tests method for explosibility of dust clouds[S]. American Society for Testing Materials, 2010: 1-2. 期刊类型引用(7)
1. 苏彬,董浩伟,罗振敏,邓军,王涛,程方明. 气粉两相体系爆炸动力学特性及机理研究进展. 化工学报. 2024(06): 2109-2122 . 百度学术
2. 武林湲,于立富,王天枢,孙威,徐建航,李航. 密闭空间油页岩粉尘爆炸特性研究. 爆炸与冲击. 2022(01): 157-166 . 本站查看
3. 徐伟巍,熊静文,刘柏清,侯照文. 酒精蒸气/烟草粉尘两相混合体系最小点火能试验研究. 火灾科学. 2022(03): 137-142 . 百度学术
4. 徐伟巍,王家祎,覃欣欣,雷成文. 酒精蒸气-烟草粉尘耦合体系燃爆猛度特性研究. 工业安全与环保. 2021(05): 9-12 . 百度学术
5. 徐伟巍,覃欣欣. 浓度与粒径对烟草粉尘爆炸压力的影响. 山东化工. 2020(12): 224-225+227 . 百度学术
6. 甘波,高伟,张新燕,姜海鹏,毕明树. 甲烷浓度对PMMA/甲烷混合爆炸下限及预热区厚度的影响. 爆炸与冲击. 2019(02): 188-195 . 本站查看
7. 叶亚明,梁峻,封昌盛,李金清,张金进. 基于几种可燃粉尘自燃温度与燃烧速率的研究. 山东化工. 2019(11): 169-172 . 百度学术
其他类型引用(4)
-