Detonation Propagation in Hydrogen/Methane-Air Mixtures in a Round Tube Filled with Orifice Plates
-
摘要: 在内径48 mm、长度5 800 mm的含环形障碍物圆管内,进行了氢气-空气及氢气-甲烷-空气的爆轰波传播试验研究,确定了爆燃转爆轰(Deflagration-to-Detonation Transition,DDT)极限。环形障碍物阻塞比为0.56,间距分为两种,即S=D和S=2D,其中S为障碍物间距,D为管道内径。火焰的速度由安装在管道壁面上的光电二极管采集得到。试验测量得到的火焰为准爆轰或阻塞火焰。在S=2D情况下得到的火焰速度均比S=D情况下的火焰速度高,并且靠近DDT极限时速度波动更明显,表明在间距较大的情况下爆轰的重起爆循环周期更长,类似于"弛振爆轰"。对于氢气-空气,障碍物间距为D时在DDT极限处有d/λ>1(富氧条件下d/λ=1.6,贫氧条件下d/λ=1.4),间距为2D时更容易形成爆轰的重起爆,在DDT极限处与准则d/λ≈1一致;对于氢气-甲烷-空气,甲烷的添加使爆轰更不稳定,对于两种间距的障碍物得到的DDT极限均有d/λ≈1(d和λ分别为障碍物内径和爆轰胞格尺寸)。说明障碍物间距对爆轰波传播有显著的影响,即间距的增大更有利于爆轰波的传播。为形成准爆轰,障碍物内径必须至少可以容纳一个爆轰胞格,同时障碍物间距足够大从而引起爆轰的重起爆。Abstract: In this study experiments were carried out in a round tube, 5 800 mm in length and 48 mm in inner-diameter, filled with orifice plates, to investigate the detonation of hydrogen-air mixtures and stoichiometric hydrogen-methane-air mixtures, and the DDT (Deflagration-to-Detonation Transition) limits were determined.The blockage ratio of the orifice plates was 0.56, and the spacing was divided into two, i.e., S=D and S=2D, in which S and D are the obstacle spacing and the tube diameter.The flame velocity was obtained using photodiodes mounted on the tube wall.The results show that the flame regime observed is the quasi-detonation or the choked flame.The flame velocity measured for S=2D is larger than that for S=D, and the velocity fluctuation is more significant.This indicates that the cycle of the detonation failure and re-initiation is longer for S=2D, which is similar to the "galloping detonation".For hydrogen-air mixtures, detonation re-initiation occurs more aptly at S=2D, and the limits correlate well with d/λ≈1.In the case of hydrogen-methane-air mixtures, the DDT limits for S=D and S=2D are both consistent with d/λ≈1, where d and λ are the inner diameter of the orifice plate and the detonation cell size.The results indicate that the obstacle spacing has a significant effect on the propagation of detonation, i.e., detonation propagates more aptly for increased spacing.To generate the quasi-detonation, the opening of the orifice plate has to be large enough to contain at least one cell size while the spacing has to be large enough to form detonation re-initiation.
-
Key words:
- gas explosion /
- orifice plates /
- DDT limits /
- hydrogen-methane
-
结构在冲击载荷作用下的动态响应一直是工程研究人员关注的热点问题。Nurick等[1]利用金属薄板开展了系列爆炸试验,建立了薄板中心挠度与无量纲冲量之间的关系,得到了方板的典型破坏模式。Ramajeyathilagam等[2]开展了水下爆炸作用下低碳钢金属板的破坏试验,通过改变爆距,依次获取了方板的塑性大变形、边界拉伸撕裂、中心拉伸撕裂等典型破坏模式。Nurick等[3]研究了水下近距离爆炸作用下圆形板的失效模式,发现随着局部爆炸强度的增加,圆形板的破坏模式由塑性大变形、冲塞逐渐过渡到花瓣开裂。Rajendran等[4]开展了水下爆炸作用下高强合金钢板的动态响应实验研究,通过理论分析建立了钢板变形计算公式。Jacob等[5]开展了空中局部爆炸作用下矩形板的实验研究,对比了钢板厚度、长宽比以及炸药形状等对钢板最终变形的影响,并得到了修正的无量纲参数。Longère等[6]开展了空中局部爆炸作用下船用钢板动态响应的试验和仿真研究。李典等[7]给出了局部载荷作用下金属板的破坏模式。在前人开展研究工作所取得成果的基础上,Jacob等[8]全面分析了板在冲击载荷作用下的失效模式。
迄今为止,结构在冲击载荷作用下的研究成果较多,但基本忽视了边界的影响。Teeling-Smith等[9]、Shen等[10]指出,边界条件对板在边界处的撕裂起到重要作用。Bonorchis等[11]通过试验分析了空爆下无固支、单边固支、两边固支条件对固支方板传递冲量的影响规律,发现支撑边界越强,作用在方板上的冲量越大。Nurick等[12]利用圆板开展了不同边界倒角下的空爆试验,通过观察圆板的减薄率,考虑边界倒角的影响后,建议将Menkes等[13]提出的毁伤模式Model 1进一步细化为:Model 1,表示边界未出现颈缩的塑性大变形;Model 1a,表示边界出现部分颈缩的塑性大变形;Model 1b,表示边界颈缩的塑性大变形。现有的研究成果虽然指出边界条件会对板的破坏模式产生较大影响,但是边界条件对方板深层次影响机理方面的研究几乎未见报道。
本研究利用落锤试验装置,通过自行设计的冲击锤头,开展不同边界支撑下(边界倒角半径分别为6、9和12 mm)固支方板的冲击试验,依据试验结果分析支撑边界对固支方板毁伤破坏模式、整体变形轮廓、边界撕裂长度等的影响,深层次分析边界对方板破坏的影响机理。
1. 试 验
1.1 试验装置
本试验利用落锤冲击机开展,试验装置如图1所示,主要包括上导轨板、质量块(1 680 kg)、下导轨板、过渡板、冲击锤头、试验板、工装框架和试验平台。试验过程中上、下导轨板与质量块通过螺栓连接固定,冲击锤头通过M20螺栓与过渡板连接,过渡板再通过4个M22螺栓与下导轨板连接,试验板通过四周满焊固定于工装框架,同时为保证冲击过程中工装框架保持平稳状态,工装框架通过螺栓连接固定在落锤工作平台区域。试验过程中通过改变落锤的高度来调整冲击强度。
本试验主要探讨不同边界支撑下固支方板的毁伤破坏,方板的安装框架设计是一大难题,本试验设计的固支方板工装框架模型如图2所示。为研究不同边界倒角支撑效果,试验共设计了3种倒角的安装条,倒角半径R分别为 6 、9和12 mm,倒角的安装位置如图2所示。安装条通过内六角螺栓与安装框架上端面连接,为避免冲击过程中发生损伤破坏,倒角安装条采用45高强钢加工而成。试验方板材料为Q345B钢,尺寸为450 mm × 350 mm,厚度为2 mm,方板四周与工装框架模型焊接固定,焊接时应保证焊接表面紧密接触,方板在长边、短边方向与框架的接触长度均为50 mm,因此,试验过程中固支方板的有效试验区域尺寸为350 mm × 250 mm。
1.2 冲击锤头设计
冲击锤头的特征尺寸对固支方板的动态响应起到重要作用。当锤头尺寸较小时,方板通常在中心区域撕裂毁伤,边界支撑对固支方板的影响难以体现,难以满足试验要求,因此冲击锤头尺寸应与方板有效试验区域相当。本试验中冲击锤头长度方向尺寸为300 mm,宽度方向尺寸为210 mm,边界处倒角半径为30 mm,锤头表面为半径300 mm的球面,材质为45钢,球头表面采用淬火处理,硬度为HRC43。冲击锤头主要通过螺栓与过渡板连接,再通过螺栓孔与下导轨板连接。
1.3 试验工况
针对倒角半径分别为6、9和12 mm的边界支撑条件,共计开展12次冲击试验,即每类边界支撑条件下分别开展4次不同强度的冲击试验,具体工况如表1所示,其中:h为落锤高度,是指锤头平面与试验方板之间的距离;v为碰撞速度;M为落锤质量,包括质量块、上下导轨、过渡板和冲击锤头的质量;E为落锤动能;D为板中心位移。
表 1 试验工况及试验结果Table 1. Experiment cases and experiment resultsExperiment
caseSpecimen No. R/mm h/mm v/(m·s−1) M/kg E/kJ Experiment results D/mm Failure mode 1 S6-1 6 400 2.80 2 024 7.934 41.2 One-side tearing 2 S6-2 6 500 3.13 2 024 9.917 45.7 Double-side tearing 3 S6-3 6 600 3.42 2 024 11.901 52.0 One-side tearing 4 S6-4 6 700 3.70 2 024 13.884 56.4 Double-side tearing 5 S9-1 9 400 2.80 2 024 7.934 39.3 Plastic deformation 6 S9-2 9 500 3.13 2 024 9.917 45.9 One-side tearing 7 S9-3 9 600 3.42 2 024 11.901 48.6 Double-side tearing 8 S9-4 9 700 3.70 2 024 13.884 55.8 Double-side tearing 9 S12-1 12 400 2.80 2 024 7.934 39.8 Plastic deformation 10 S12-2 12 500 3.13 2 024 9.917 44.6 One-side tearing 11 S12-3 12 600 3.42 2 024 11.901 51.1 One-side tearing 12 S12-4 12 700 3.70 2 024 13.884 56.9 One-side tearing 2. 结果与讨论
2.1 毁伤破坏模式分析
强冲击载荷作用下不同边界支撑下固支方板的试验结果见表1,试验中固支方板出现的破坏模式主要有3种:塑性大变形、单边撕裂和双边撕裂,如图3和图4所示。
对于试验中出现的如图3所示的塑性大变形破坏模式,方板整体向内凹陷,长边方向为典型的抛物线形状变形模式,与短边方向的变形模式差异较大,近边界处凹陷坡度较大,中间较为平缓,这主要是由于受到冲击锤头的挤压形成的。此外,方板出现较明显的4条塑性铰线,塑性铰1、塑性铰2分别是由于试验方板与锤头边界和框架倒角边界相互作用而形成的,4条塑性铰中间部分均产生一定程度的颈缩,类似塑性铰在另两类破坏模式中同样存在。
对于试验中出现的如图4所示的单边撕裂和双边撕裂破坏模式,由于受到试验方板制造工艺、模型焊接、落锤的实际操作等的影响,试验方板的破坏并未表现出对称性,如图4(a)中所示的S6-1试件,一边裂缝长度为207 mm,一边未出现裂缝;而图4(b)中所示的S6-4试件,双边撕裂裂缝长度分别为294和184 mm。固支方板的撕裂区域位于短边方向,裂缝基本处于倒角边缘位置,通过对裂缝局部放大后观察可以发现:在长度方向,裂缝中间区域破口基本呈45°,为典型的拉伸撕裂破坏,但随着裂缝向两端扩展,裂缝在厚度方向却呈现出两种不同失效模式混合的现象,在背离锤头面为拉伸撕裂,而在迎锤头面却表现为一定的剪切失效破坏,即图中的白色区域,此区域破口基本呈90°。
图5、图6、图7分别显示同一落锤高度、不同边界倒角支撑作用下固支方板的毁伤破坏模式。综合对比可知:(1)边界倒角半径越小,固支方板在冲击载荷作用下越容易出现撕裂破坏,当边界支撑倒角半径为6 mm时,落锤高度为400 mm时方板即出现撕裂,而当倒角半径为9和12 mm时,需要落锤高度达到500 mm方板才能出现撕裂;(2)边界倒角半径对固支方板的撕裂位置具有一定影响,边界倒角半径越大,撕裂位置越向内部移动,说明边界倒角半径可改变方板边界处的应力集中位置。
2.2 变形情况分析
不同边界倒角半径和落锤高度对固支方板中心点位移的影响如图8所示。由图8可知,边界支撑对方板中心点位移的影响无明显规律,如落锤高度为400 mm时,虽然试件中心点位移随着边界倒角的增大先减小后增大,但总体而言,不同倒角下方板中心点的位移相差不大,最大偏差仅为3.4 mm,边界倒角的大小对固支方板中心点位移的影响基本可忽略,其他落锤高度同样适用。这主要是由于固支方板主要依托膜应力、塑性铰来消耗冲击能量,边界倒角虽然能改变塑性变形集中位置,但整体变形模式几乎不受影响。由图8(b)可以看出,试件中心点位移随着落锤高度的增加基本呈线性增长。
依据图8的结果,进一步分析边界倒角对固支方板整体变形模式的影响,图9和图10分别表示落锤高度为500和700 mm时,边界倒角对固支方板整体变形的影响。由图9和图10可以看出,对于同一落锤高度,无论在长边方向还是在短边方向,不同倒角支撑下固支方板的变形轮廓基本保持一致,可见边界倒角大小对方板整体变形模式的影响较小。
2.3 边界撕裂长度分析
强冲击载荷作用下边界倒角对固支方板的撕裂长度有较大影响。表2给出了不同工况下固支方板试件的裂缝长度统计结果,落锤冲击作用下,固支方板短边发生一边撕裂或两边撕裂具有一定的随机性,且即使发生两边撕裂,裂缝也呈现明显的不对称性。图11显示了边界倒角半径对试件裂缝总长度(包含所有裂缝长度)的影响规律。由图11可知,对于同一落锤高度,试件裂缝总长度随着边界倒角半径的增大基本呈减小趋势,出现的裂缝中,最长裂缝长度为294 mm,最短为125 mm。图12显示了边界倒角半径对长裂缝长度的影响规律。由图12可知,对于同一落锤高度,固支方板试件出现的长裂缝长度随着边界倒角半径的增大基本呈减小趋势。
表 2 不同工况下固支方板的裂缝长度Table 2. Crack length of square plates under different experiment casesExperiment
caseSpecimen
No.R/mm h/mm Failure mode Crack length/mm Crack 1 Crack 2 1 S6-1 6 400 One-side tearing 207 2 S6-2 6 500 Double-side tearing 216 125 3 S6-3 6 600 One-side tearing 290 4 S6-4 6 700 Double-side tearing 294 184 5 S9-1 9 400 Plastic deformation 6 S9-2 9 500 One-side tearing 195 7 S9-3 9 600 Double-side tearing 250 213 8 S9-4 9 700 Double-side tearing 184 132 9 S12-1 12 400 Plastic deformation 10 S12-2 12 500 One-side tearing 170 11 S12-3 12 600 One-side tearing 230 12 S12-4 12 700 One-side tearing 280 2.4 应力状态分析
边界支撑倒角对固支方板的撕裂破坏起着重要作用,边界支撑倒角越小,方板越容易撕裂。本节主要借助数值仿真手段分析不同边界支撑条件对固支方板的深层次影响机理。
图13显示了落锤高度为400 mm,边界倒角半径为6、9和12 mm时固支方板的仿真计算结果(等效塑性应变PEEQ云图),为便于比较不同边界倒角半径对固支方板动态响应的影响,仿真计算中并未采用失效模式,因此边界倒角半径为6 mm时的仿真计算结果与试验结果存在一定差别。考虑到应力状态对单元失效存在较大影响[14-16],将方板最大等效塑性应变所处单元(中心对称处)的应力三轴度提取出来进行分析,如图14所示。由图14可以看出,边界倒角越小,方板边界处的等效塑性应变越大。锤头与试验方板开始接触后,3种不同倒角下单元所处应力三轴度基本处于上升阶段;之后锤头与试验方板之间的作用相对平稳,各单元应力三轴度趋于稳定,波动基本处于[0.5,0.6]区间;随后由于锤头首次弹起,单元所处应力三轴度先小幅上升再下降。结合本试验研究结果来看,在锤头与试验方板相对平稳作用期间,即最大等效应变单元处的应力三轴度处于[0.5,0.6]区间时,试验方板发生撕裂,3种倒角下的等效塑性应变单元所处应力三轴度基本一致。
表3为边界倒角半径不同时最大塑性应变单元的详细结果对比。
、 分别表示方板厚度方向(1方向)、短边方向(3方向)的拉伸应变, 表示剪切应变,1、2、3所指方向如图13所示。由对比结果可以看出:边界倒角半径改变时, 、 几乎未改变,基本处于[0.130,0.141]区间,而剪切应变 由0.210上升至0.268,提高了27.6%,倒角半径越小,边界对方板的剪切作用越明显。考虑到落锤高度为400 mm、倒角半径为9和12 mm时固支方板并未发生撕裂,可推断方板边界撕裂的临界应变位于[0.191,0.241]区间,且撕裂时刻所处应力三轴度位于[0.5,0.6]区间。表 3 最大塑性应变单元结果对比Table 3. Comparison of the max PEEQ resultsExperiment
caseR/mm h/mm Plastic deformation Stress triaxility PEEQ 1 6 400 0.241 –0.137 0.139 0.268 [0.5, 0.6] 5 9 400 0.191 –0.133 0.130 0.219 [0.5, 0.6] 9 12 400 0.185 –0.141 0.133 0.210 [0.5, 0.6] 3. 结 论
通过设计新型锤头,利用落锤试验机开展了倒角半径为6、9和12 mm固支方板的冲击试验,探讨了边界倒角对固支方板毁伤破坏模式、变形轮廓、撕裂长度、临界撕裂等的影响,得出以下结论。
(1)落锤冲击作用下固支方板会出现塑性大变形、单边撕裂、两边撕裂3种典型的破坏模式,裂缝中间区域破口基本呈45°,表现出明显的拉伸撕裂破坏,随着裂缝向两端扩展,裂缝在厚度方向表现出两种不同的失效模式,在背离锤头面为拉伸撕裂,而在迎锤头面表现为一定的剪切失效破坏,这与爆炸冲击载荷方板撕裂破口区域的毁伤破坏模式存在一定的不同。
(2)边界倒角对固支方板的中心点位移、整体变形模式影响较小,主要对边界处的裂缝产生较大影响,对于同一落锤高度,试件裂缝总长度随着边界倒角半径的增大呈减小趋势,长裂缝长度随着边界倒角半径的增大呈减小趋势。
(3)结合数值仿真计算发现,边界倒角主要改变单元的剪切应变,即倒角半径越小,产生的剪切效果越明显。依据本研究的分析结果,方板边界撕裂时刻所处应力三轴度位于[0.5,0.6]区间,临界撕裂应变位于[0.191,0.241]区间。
-
表 1 氢气-空气的DDT极限
Table 1. DDT limits for hydrogen-air mixtures
Obstacle spacing Lean limit/% φ d/λ L/λ Rich limit/% φ d/λ L/λ D 22 0.67 1.6 7.0 48 2.20 1.4 6.4 2D 21 0.63 1.1 7.2 49 2.29 1.0 6.6 表 2 化学计量比下氢气-甲烷-空气的DDT极限
Table 2. DDT limits for stoichiometric hydrogen-methane-air mixtures
Obstacle spacing Limit (X) d/λ L/λ D 0.75 0.9 4.0 2D 0.75 0.9 6.1 -
[1] CICCARELLI G, DOROFEEV S.Flame acceleration and transition to detonation in ducts[J].Progress in Energy and Combustion Science, 2008, 34(4):499-550. doi: 10.1016/j.pecs.2007.11.002 [2] CROSS M, CICCARELLI G.DDT and detonation propagation limits in an obstacle filled tube[J].Journal of Loss Prevention in the Process Industries, 2015, 36:380-386. doi: 10.1016/j.jlp.2014.11.020 [3] PERALDI O, KNYSTAUTAS R, LEE J H.Criteria for transition to detonation in tubes[J].Symposium (International) on Combustion, 1988, 21(1):1629-1637. doi: 10.1016/S0082-0784(88)80396-5 [4] DOROFEEV S B, SIDOROV V P, KUZNETSOV M S, et al.Effect of scale on the onset of detonations[J].Shock Waves, 2000, 10(2):137-149. doi: 10.1007/s001930050187 [5] KARIM G A, WIERZBA I, AL-ALOUSI Y.Methane-hydrogen mixtures as fuels[J].International Journal of Hydrogen Energy, 1996, 21(7):625-631. doi: 10.1016/0360-3199(95)00134-4 [6] YU M, ZHENG K, ZHENG L, et al.Effects of hydrogen addition on propagation characteristics of premixed methane/air flames[J].Journal of Loss Prevention in the Process Industries, 2015, 34:1-9. doi: 10.1016/j.jlp.2015.01.017 [7] DI SARLI V, DI BENEDETTO A.Laminar burning velocity of hydrogen-methane/air premixed flames[J].International Journal of Hydrogen Energy, 2007, 32(5):637-646. doi: 10.1016/j.ijhydene.2006.05.016 [8] BOZIER O, SORIN R, VIROT F, et al. Detonability of binary H2/CH4-air mixtures[C]//Proceedings of Third ICHS, ID, 2009: 188. https: //www. h2tools. org/content/detonability-binary-h2-ch4-air-mixtures [9] HU E, HUANG Z, LIU B, et al.Experimental study on combustion characteristics of a spark-ignition engine fueled with natural gas-hydrogen blends combining with EGR[J].International Journal of Hydrogen Energy, 2009, 34(2):1035-1044. doi: 10.1016/j.ijhydene.2008.11.030 [10] WU L, KOBAYASHI N, LI Z, et al.Experimental study on the effects of hydrogen addition on the emission and heat transfer characteristics of laminar methane diffusion flames with oxygen-enriched air[J].International Journal of Hydrogen Energy, 2016, 41(3):2023-2036. doi: 10.1016/j.ijhydene.2015.10.132 [11] WANG L, MA H, SHEN Z, et al.Experimental investigation of methane-oxygen detonation propagation in tubes[J].Applied Thermal Engineering, 2017, 123:1300-1307. doi: 10.1016/j.applthermaleng.2017.05.045 [12] WANG L, MA H, SHEN Z, et al.Detonation characteristics of stoichiometric H2-O2 diluted with Ar/N2 in smooth and porous tubes[J].Experimental Thermal and Fluid Science, 2018, 91:345-353. doi: 10.1016/j.expthermflusci.2017.08.021 [13] GAO Y, NG H D, LEE J H S.Minimum tube diameters for steady propagation of gaseous detonations[J].Shock Waves, 2014, 24(4):447-454. doi: 10.1007/s00193-014-0505-8 [14] ZHANG B, SHEN X, PANG L, et al.Methane-oxygen detonation characteristics near their propagation limits in ducts[J].Fuel, 2016, 177:1-7. doi: 10.1016/j.fuel.2016.02.089 [15] KEE R J, RUPLEY F M, MILLER J A. Chemkin-Ⅱ: a fortran chemical kinetics package for the analysis of gas-phase chemical kinetics: SAND 89-8009[R]. Albuquerque, NM: Sandia National Laboratories, 1989. https: //www. osti. gov/biblio/5681118 [16] MORLEY C. Gaseq: a chemical equilibrium program for Windows[Z]. 2005. [17] CICCARELLI G, WANG Z, LU J, et al.Effect of orifice plate spacing on detonation propagation[J].Journal of Loss Prevention in the Process Industries, 2017, 49(B):739-744. https://www.sciencedirect.com/science/article/pii/S0950423017302619 [18] GU L S, KNYSTAUTAS R, LEE J H S.Influence of obstacle spacing on the propagation of quasi-detonation[J].Dynamics of Explosions of Progress in Astronautics and Aeronautics, 1988, 114:232-247. doi: 10.2514/5.9781600865886.0232.0247 [19] ZHANG B, KAMENSKIHS V, NG H D, et al.Direct blast initiation of spherical gaseous detonations in highly argon diluted mixtures[J].Proceedings of the Combustion Institute, 2011, 33(2):2265-2271. doi: 10.1016/j.proci.2010.06.165 [20] CICCARELLI G, CROSS M.On the propagation mechanism of a detonation wave in a round tube with orifice plates[J].Shock Waves, 2016, 26(5):587-597. doi: 10.1007/s00193-016-0676-6 [21] GAO Y, LEE J H S, NG H D.Velocity fluctuation near the detonation limits[J].Combustion and Flame, 2014, 161(11):2982-2990. doi: 10.1016/j.combustflame.2014.04.020 -