大型水平管道中环氧丙烷-铝粉-空气混和物爆燃转爆轰的实验研究

陈默 白春华 刘庆明

杜明燃, 陈宇航, 陆少锋, 梁进, 李基锐, 王尹军, 王天照, 陈智凡. 基于正交试验法的气泡帷幕削波特性研究[J]. 高压物理学报, 2023, 37(6): 065302. doi: 10.11858/gywlxb.20230684
引用本文: 陈默, 白春华, 刘庆明. 大型水平管道中环氧丙烷-铝粉-空气混和物爆燃转爆轰的实验研究[J]. 高压物理学报, 2011, 25(4): 359-364 . doi: 10.11858/gywlxb.2011.04.012
DU Mingran, CHEN Yuhang, LU Shaofeng, LIANG Jin, LI Jirui, WANG Yinjun, WANG Tianzhao, CHEN Zhifan. Bubble Curtain Clipping Characteristics Based on Orthogonal Test Method[J]. Chinese Journal of High Pressure Physics, 2023, 37(6): 065302. doi: 10.11858/gywlxb.20230684
Citation: CHEN Mo, BAI Chun-Hua, LIU Qing-Ming. Experimental Research of Deflagration to Detonation Transition of Epoxypropane-Aluminum-Air Mixtures in Large-Scale Tube[J]. Chinese Journal of High Pressure Physics, 2011, 25(4): 359-364 . doi: 10.11858/gywlxb.2011.04.012

大型水平管道中环氧丙烷-铝粉-空气混和物爆燃转爆轰的实验研究

doi: 10.11858/gywlxb.2011.04.012
详细信息
    通讯作者:

    刘庆明

Experimental Research of Deflagration to Detonation Transition of Epoxypropane-Aluminum-Air Mixtures in Large-Scale Tube

More Information
    Corresponding author: LIU Qing-Ming
  • 摘要: 在长为32.4 m、内径为0.199 m的大型长直水平管道中,对环氧丙烷-空气两相流云雾及环氧丙烷-铝粉-空气三相流云雾的爆燃转爆轰(DDT)过程进行了实验研究。对弱点火条件下多相混和物DDT过程的不同阶段特征进行了分析,对比研究了不同浓度时混和物的燃爆情况。结果表明:浓度为513 g/m3的环氧丙烷-空气混和物及浓度为237和643 g/m3的环氧丙烷-铝粉-空气混和物均能在管道中完成爆燃向爆轰的转变,进入自持爆轰阶段,其胞格尺寸分别为0.28和0.50 m。

     

  • 港口码头、水利水电设施和航道疏浚等基础建设工程在我国的发展规划中占据重要地位,然而,其作业环境复杂,难以通过机械拆除手段实现工程目标。水下爆破技术具有适用范围广、价格低廉等优势,能够很好地适应复杂环境,被广泛地应用于工程建设[13]。水下爆炸与其他介质中的爆炸不同,所产生的冲击波具有峰值压力更高、影响范围更广、毁伤效应更强的特点[45]。水下爆破作业过程的主要危害因素有水中冲击波、水面涌浪、爆破飞石等[6]。其中,水中冲击波具有毁伤作用强、作用时间长的特点[7],如何减小水中冲击波的危害是目前水下爆破领域亟待解决的重要问题之一[89]

    气泡帷幕的概念最早由加拿大工程师Adolph提出,用于Oratario水电站的水下爆破,由于其对水中冲击波的削减效果优异,因而得到广泛认可和应用[10]。为进一步研究气泡帷幕的削波作用,国内学者也开展了相关研究。刘欣等[11]通过改变气泡帷幕与被保护对象之间的距离,研究了气泡帷幕对水中冲击波衰减效果的影响,发现气泡帷幕与被保护物的间距越小,水中冲击波的削减效果越明显。胡伟才等[12]利用LS-DYNA软件,分析了气泡帷幕的数量、间距和防护距离对桥梁振动速度削减作用的影响,发现气泡帷幕数量的影响最大,气泡帷幕间距次之,气泡帷幕防护距离的影响最小。李泽华等[13]研究了气泡与水中冲击波的相互作用过程,发现气泡的强力压缩和破碎是削减水中冲击波的最佳方式。谢达建等[14]利用LS-DYNA软件建立了水下钻孔爆破模型,结合长江九朝段炸礁工程,分析了气泡帷幕的距离对削波效果的影响,发现在距离被保护对象较近处设置气泡帷幕时的防护效果更好。陆少锋等[15]研究了不同供风量形成的气泡帷幕对水中冲击波的削减效应,发现气泡帷幕的削波效果随着供风量的增加而增强。

    综上所述,气泡帷幕的削波效果与众多因素有关,而气泡帷幕层数、气泡帷幕爆心距和药包深度作为工程实践中的重要因素,对气泡帷幕削波效果的影响以及影响程度尚待深入研究。为此,本研究利用AUTODYN软件建立水下爆炸模型,设计现场试验对模型进行验证,通过3因素3水平正交试验,对气泡帷幕层数、气泡帷幕爆心距和药包深度对削波效果的影响进行敏感性分析,以期获得削波效果最优的组合方案,为气泡帷幕在水下爆炸中的应用提供参考和理论依据。

    1.1.1   炸药材料模型

    炸药采用JWL状态方程[16]描述,其表达式为

    p=A(1ωR1V)eR1V+B(1ωR2V)eR2V+EV
    (1)

    式中:p为冲击波压力,A、B、R1、R2ω为JWL状态方程参数,E为炸药内能,V为当前的相对体积。工程中常用2号岩石乳化炸药,具体参数见表1,其中:ρ为密度,D为爆速,pC-J为C-J爆轰压力。

    表  1  炸药的材料参数
    Table  1.  Material parameters of explosive
    ρ/(g·cm−3)D/(km·s−1)pC-J/GPaA/GPaB/GPaR1R2ωE/(kJ·g−1)
    1.304.59.80214.40.1824.150.950.154.19
    下载: 导出CSV 
    | 显示表格
    1.1.2   水材料模型

    水采用多项式状态方程[16]描述。当水被压缩时(压缩比μ>0),其状态方程为

    pw=A1μ+A2μ+A3μ+(B0+B1μ)ρ0e
    (2)

    当水膨胀时(μ<0),其状态方程为

    pw=T1μ+T2μ+B0ρ0e
    (3)

    当水既不压缩也不膨胀时(μ=0),其状态方程可简化为

    pw=B0ρ0e
    (4)

    式中:pw为水的压力;压缩比μ=ρ/ρ0−1,ρ为水的当前密度,ρ0为水的初始密度;A1A2A3B0B1T1T2均为常数,如表2所示;e为水的比内能,e=(p0 +ρ0gH)/B0ρ0,其中p0为大气压,g为重力加速度,H为水深(即药包深度)。根据药包深度调节水的比内能e,当H=4.5 m时,水的比内能为519.38 J/kg。

    表  2  水的材料参数
    Table  2.  Material parameters of water
    ρ/(g·cm−3)A1/GPaA2/GPaA3/GPaB0B1T1/GPaT2/GPaω
    1.0002.29.5414.570.280.282.200.150
    下载: 导出CSV 
    | 显示表格
    1.1.3   气泡帷幕材料模型

    由刘欣等[11]和谢达建等[14]的研究可知,数值模拟中可以利用空气层近似替代气泡帷幕。选用AUTODYN材料库中的AIR空气模型和理想气体状态方程[17],其表达式为

    p=(γ1)ρaρa0Ea
    (5)

    式中:Ea为空气的比内能;γ为绝热指数,取1.4;ρaρa0分别为空气的当前密度和初始密度,ρa0取1.225 kg/m3

    药包在无限水域中爆炸时,球形药包在平面内关于x轴和y轴对称,为了减少运算量,利用AUTODYN建立如图1所示的1/4轴对称计算模型,模型尺寸为20 m×20 m。炸药和水均采用欧拉网格划分,网格大小设为10 mm。炸药选用300 g的2号岩石乳化炸药,起爆点设在药包球心,选用Transmit边界条件,采用mm-mg-ms单位制[1819]。通过调节空气层的厚度实现不同气泡帷幕层数的模拟:气泡帷幕层数为1(N=1)时,空气层厚度为10 cm,气泡帷幕层数为2(N=2)时,空气层厚度为20 cm,依此类推。

    图  1  数值计算模型示意图
    Figure  1.  Schematic diagram of the numerical simulation model

    为验证数值模拟的合理性,在某河道进行了气泡帷幕水下爆炸冲击波测试。试验设置3种工况:工况1为1层气泡帷幕,工况2为2层气泡帷幕,工况3为3层气泡帷幕。每种工况测试2次,共计6次试验。试验区域的水深为9~10 m,平均水流速度为0.1 m/s。传感器和药包均设置在4.5 m水深处,且位于同一水平线上。传感器选用PCB-138A01型压力传感器。冲击波测试仪为4通道Blast-PRO型测试仪,测试过程中,触发电平设置为0.2%,量程设置为10 V,记录时长为0.4 s。炸药选用300 g 2号岩石乳化炸药。气泡帷幕设置在距药包中心6 m处。选取C1C2两个测点,其中,测点C1位于药包左侧12 m处,测点C2位于药包右侧12 m处。试验布局如图2所示。现场试验与数值模拟结果列于表3,其中pmax为冲击波峰值压力,δ为冲击波峰值压力削减率。以1层气泡帷幕为例,测点C1C2的压力时程曲线如图3所示。

    图  2  气泡帷幕试验布局
    Figure  2.  Layout of the bubble curtain test
    表  3  不同气泡帷幕层数下水下爆炸冲击波峰值压力
    Table  3.  Peak pressure of underwater blast shock wave for different bubble curtain layers
    Test No. N On-site test Simulation
    pmax/MPa δ/% pmax/MPa δ/%
    C1 C2 C1 C2
    1 1 1.883 0.154 91.82 1.750 0.291 83.37
    2 1 1.789 0.228 87.26 1.750 0.291 83.37
    3 2 1.675 0.157 90.63 1.750 0.116 93.37
    4 2 1.486 0.108 92.73 1.750 0.116 93.37
    5 3 1.862 0.247 86.73 1.750 0.248 85.83
    6 3 1.764 0.196 88.89 1.750 0.248 85.83
    下载: 导出CSV 
    | 显示表格
    图  3  设置1层气泡帷幕时测点C1C2处的压力时程曲线
    Figure  3.  Shock wave pressure-time curves at measurement points C1 and C2 for one layer of bubble curtain

    表3可以看出:设置1层气泡帷幕时,2次现场试验测得C1处的峰值压力分别为1.883和1.789 MPa,C2处的峰值压力分别为0.154和0.228 MPa,削减率分别为91.82%和87.26%;数值模拟得到C1C2处的峰值压力分别为1.750和0.291 MPa,削减率为83.37%。以峰值压力的削减率为误差评判指标,数值模拟与试验结果的相对误差分别为9.20%和4.46%。

    设置2层气泡帷幕时,现场试验得到的削减率为90.63%和92.73%,数值模拟得到削减率为93.37%,相对误差为3.02%和0.69%。设置3层气泡帷幕时,现场试验得到的削减率为86.73%和88.89%,数值模拟得到的削减率为85.83%,相对误差为1.04%和3.44%。

    由上述分析可知,采用空气层替代本试验工况中的气泡帷幕时,误差较小,说明本数值模型可以很好地模拟实际工况,因此后续正交试验均采用该方法进行模拟运算。

    为减少试验次数,同时保证试验的可靠性,选用L9(34)正交表,利用AUTODYN软件设计了3因素3水平正交试验方案。选取气泡帷幕层数(N)、气泡帷幕爆心距(D)、药包深度(H)3个因素,分别设为因素A、B、C,其中:气泡帷幕层数分为3种水平,即N取1、2、3;气泡帷幕爆心距有3种水平,分别取1、3、5 m;药包深度有3种水平,分别为2.5、6.5、10.5 m。各因素之间无交互作用。此外,设置未加气泡帷幕的空白对照组。如图1所示,距药包中心6、9和12 m分别设置测点1、测点2和测点3,取冲击波峰值压力削减率δ作为评价指标,正交试验设计因素和因素水平见表4,正交试验方案见表5

    表  4  正交试验设计因素和水平
    Table  4.  Orthogonal test design factors and levels
    LevelFactor
    ND/mH/m
    1112.5
    2236.5
    33510.5
    下载: 导出CSV 
    | 显示表格
    表  5  正交试验方案
    Table  5.  Orthogonal test protocols
    No. Level Test programme No. Level Test programme
    N D H N D H
    1 1 1 1 A1B1C1 6 2 3 1 A2B3C1
    2 1 2 2 A1B2C2 7 3 1 3 A3B1C3
    3 1 3 3 A1B3C3 8 3 2 1 A3B2C1
    4 2 1 2 A2B1C2 9 3 3 2 A3B3C2
    5 2 2 3 A2B2C3
    下载: 导出CSV 
    | 显示表格

    利用水的多项式状态方程,可以计算出药包深度为2.5、6.5和10.5 m时,水的比内能分别为449.38、589.38和729.38 J/kg。由此通过修改AUTODYN材料库中水介质的比内能[2022]模拟炸药的不同水深环境。在空白对照组中,水下爆炸冲击波峰值压力(pmax)及其到达时间(tp)如表6所示,压力时程曲线如图4所示。

    表  6  空白对照组数据
    Table  6.  Data of blank control group
    No. H/m Point pmax/MPa tp/ms No. H/m Point pmax/MPa tp/ms
    1 2.5 1 4.108 3.952 3 10.5 1 4.187 3.954
    2 2.456 6.001 2 2.533 6.001
    3 1.531 8.006 3 1.708 8.002
    2 6.5 1 4.148 3.954
    2 2.494 6.002
    3 1.688 8.016
    下载: 导出CSV 
    | 显示表格
    图  4  无气泡帷幕时冲击波压力时程曲线
    Figure  4.  Time course curves of shock wave pressure without bubble curtain

    表6图4可以看出:冲击波压力时程曲线的波形一致;在浅水区域,随着水下爆炸药包深度的增加,冲击波到达测点的时间基本保持不变,而冲击波峰值压力增大。以测点1为例,当药包深度为2.5、6.5和10.5 m时,冲击波到达峰值压力的时间分别为3.952、3.954和3.954 ms,峰值压力分别为4.108、4.148和4.187 MPa。

    按照表5设计正交试验,将所得数据绘制压力时程曲线,如图5所示。可以看出,气泡帷幕层数、气泡帷幕爆心距和药包深度对气泡帷幕的削波效果均有影响。以正交试验1为例:对于空白对照试验1(未设置气泡帷幕,H=2.5 m),测点1、测点2和测点3的冲击波峰值压力分别为4.108、2.456和1.531 MPa;对于正交试验1(N=1,D=1 m,H=2.5 m),测点1、测点2和测点3的冲击波峰值压力分别为0.431、0.314和0.233 MPa,削减率分别为89.51%、87.21%和84.78%,平均削减率为87.17%。同理,求得其余8个试验的冲击波峰值压力削减率,如表7所示。

    图  5  有气泡帷幕时的冲击波压力时程曲线
    Figure  5.  Time course curves of shock wave pressure with bubble curtain
    表  7  正交试验数据
    Table  7.  Orthogonal test data
    Test No.ProgrammeND/mH/mpmax/MPaδ/%
    Point 1Point 2Point 3
    1A1B1C1112.50.3410.3140.29586.55
    2A1B2C2136.50.3200.2780.26188.56
    3A1B3C31510.50.3550.2860.26488.26
    4A2B1C2216.50.1220.1130.10695.42
    5A2B2C32310.50.1320.1180.10795.31
    6A2B3C1252.50.1460.1210.10994.80
    7A3B1C33110.50.2750.2600.25089.51
    8A3B2C1332.50.3720.3320.30885.77
    9A3B3C2356.50.4120.3050.27487.20
    下载: 导出CSV 
    | 显示表格

    表7所示,气泡帷幕防护技术能够很好地削弱水下爆炸冲击波的峰值压力,削减率可达95.42%;在所设置的9组正交试验中,正交试验4(N=2,D=1 m,H=6.5 m)的削减率最大,防护效果最好。设置1、2、3层气泡帷幕时,平均削减率分别为86.55%~88.56%、94.80%~95.42%、85.77%~89.51%。由此可见,气泡帷幕的冲击波峰值压力削减率存在最大值,到达最佳防护效果之后,继续增加气泡帷幕数量反而出现负效应。

    另外,以正交试验1、2、3为例,测点1处,D为1、3、5 m时,冲击波到达峰值压力的时间分别为3.507、2.058、0.652 ms,说明气泡帷幕爆心距越小,水下爆炸冲击波到达峰值压力的时间越长。

    Ki表示因素j取第i个水平时的评价指标之和,ki表示因素j取第i个水平时评价指标的平均值。通过极差Rj分析各个因素对评价指标的影响程度,Rj越大,则表明因素j对评价指标的影响越明显[2325]

    设评价指标为冲击波峰值压力的平均削减率δave,根据表7中的正交试验数据,得到极差分析结果,如表8所示,水平与指标的关系如图6所示。

    表  8  极差分析结果
    Table  8.  Results of the variance analysis
    FactorK1K2K3k1k2k3R
    N263.37285.53262.4887.7995.1887.497.68
    D271.48269.64270.2690.4989.8890.090.61
    H267.12271.18273.0889.0490.3991.031.99
    下载: 导出CSV 
    | 显示表格
    图  6  水平与指标的关系
    Figure  6.  Relationship between level and indicators

    表8图6可知,气泡帷幕层数、气泡帷幕爆心距和药包深度的极差分别为7.68、0.61、1.99,说明气泡帷幕层数对气泡帷幕削波效应的影响最大,药包深度次之,气泡帷幕爆心距的影响最小。图6(a)显示,削减率随气泡帷幕层数的增加先增大后减小,进一步验证了当气泡帷幕削波效果达到最大后继续增加气泡帷幕层数反而会降低削波效果的结论。另外,气泡帷幕爆心距越小,药包深度越大,则气泡帷幕削波效果越好。根据表8可以推断,当气泡帷幕层数为2,气泡帷幕爆心距为1 m,药包深度为10.5 m时,削波效果最好。

    在方差分析中,取显著性水平P为0.05,查表得F0.05为19。通过计算因素j的离差平方和Sj、自由度μj和均方ˉSj,得到因素jFj值,与显著性水平F0.05进行对照,若Fj>F0.05,则说明因素j对气泡帷幕削波效果的影响显著;反之,则影响不显著。因素j的均方¯Sj越大,其对气泡帷幕削波效果的影响越大。

    Sj=1nni=1K2ij1h(hi=1yi)2
    (6)
    S=hi=1y2i1h(hi=1yi)2
    (7)
    Se=S3j=1Sj
    (8)

    式中:Kij为第j个因素下第i个水平所对应的试验指标;yi为第i次试验的试验指标;n为水平数,n=3;h为正交试验次数,h=9;S为总离差平方和;Se为误差平方和。

    fe=ft3j=1fj
    (9)

    式中:fe为误差项自由度;ft为总自由度,ft=h−1;fj=n−1。

    ˉSj=Sjfj
    (10)
    ˉSe=Sefe
    (11)
    Fj=ˉSjˉSe
    (12)

    式中:ˉSe为误差项的均方。F统计量服从第一自由度为fj、第二自由度为feF分布,显著性水平取0.05时,对F统计量进行检验。

    利用式(9)~式(17)计算离差平方和S、自由度f、均方ˉSF[2628],所得结果见表9

    表  9  方差分析结果
    Table  9.  Results of variance analysis
    Sources of varianceSfˉSFSignificance-effect
    N113.68256.8438.81Yes
    D3.4821.741.19No
    H6.1823.092.11No
    Error2.9321.46
    Sum126.278
    下载: 导出CSV 
    | 显示表格

    根据表9中的ˉS可以看出,正交试验所选因素对气泡帷幕削波效果的影响由大到小依次为气泡帷幕层数、药包深度、气泡帷幕爆心距。通过方差分析可以看出:气泡帷幕层数的显著性水平F为38.81,大于F0.05,表明气泡帷幕层数对削波效果有显著影响;药包深度和气泡帷幕爆心距的显著性水平分别为1.19和2.11,均小于F0.05,表明二者对气泡帷幕削波效果的影响不大,与极差分析所得结论一致。

    基于正交试验极差分析,利用AUTODYN显式有限元分析程序,建立了自由水域内药包爆炸的1/4轴对称计算模型,探究了气泡帷幕层数、气泡帷幕爆心距和药包深度对气泡帷幕削波效应的影响和敏感性,得出以下主要结论。

    (1) 气泡帷幕可以有效地削减水中冲击波峰值压力,削减率可达95.42%,从而有效降低水中冲击波对爆破区域附近被保护对象的毁伤作用。

    (2) 气泡帷幕的削波效果与气泡帷幕层数、气泡帷幕爆心距和药包深度均有关。气泡帷幕爆心距越小,药包深度越大,则气泡帷幕的削波效果越好;而气泡帷幕层数与气泡帷幕削波效果并不成正比,当气泡帷幕层数为1、2、3时,冲击波峰值压力的平均削减率分别为87.79%、95.17%和87.49%,在实际工程中应合理选择气泡帷幕层数。

    (3) 正交试验分析显示:气泡帷幕层数对气泡帷幕削波效应的影响最大,药包深度次之,气泡帷幕爆心距的影响最小;当气泡帷幕层数为2,气泡帷幕爆心距为1 m,药包深度为10.5 m时,削波效果最好。

  • Borisov A A, Khasainov B A, Veyssiere B. On Detonation of Aluminum Dust in Air and Oxygen [J]. Sov J Chem Phys, 1992, 10(2): 369-402.
    Zhang F, Gronig H. Transition to Detonation in Corn Starch Dust-Oxygen and -Air Mixtures [J]. Combust Flame, 1991, 86(1-2): 21-32.
    Zhang F, Gronig H. Two-Headed Detonation in Reactive Particle-Oxidizing Gas Flow [J]. Phys Fluids A, 1992, 4(10): 2308-2315.
    Zhang F, Gronig H. Spin Detonation in Reactive Particles-Oxidizing Gas Flow [J]. Phys Fluids A, 1991, 3(8): 1983-1990.
    Hui J M, Liu R H, Ge G L. Study on Increasing the Power of FAE (Ⅰ)-Selection of High Energy Fuels [J]. Journal of Nanjing University of Science and Technology, 1995, 19(5): 472-476. (in Chinese)
    惠君明, 刘荣海, 葛桂兰. 提高FAE威力的研究(Ⅰ)--高能量燃料的选择 [J]. 南京理工大学学报, 1995, 19(5): 472-476.
    Mogi T, Horiguchia S. Explosion and Detonation Characteristics of Dimethyl Ether [J]. J Hazard Mater, 2009, 164(1): 114-119.
    Xu X F, Xie L F, Peng J H, et al. Study on the Detonation Cellular Structure of Propylene Epoxide-Air [J]. Explosion and Shock Waves, 2004, 24(2): 158-162. (in Chinese)
    徐晓峰, 解立峰, 彭金华, 等. 环氧丙烷-空气混合物爆轰波胞格结构的研究 [J]. 爆炸与冲击, 2004, 24(2): 158-162.
    Liu Q M, Bai C H, Jiang L, et al. Deflagration-to-Detonation Transition in Nitromethane Mmist/Aluminum Dust/Air Mixtures [J]. Combust Flame, 2010, 157(1): 106-117.
    Liu Q M, Li X D, Bai C H. Deflagration to Detonation Transition in Aluminum Dust-Air Mixture under Weak Ignition Condition [J]. Combust Flame, 2009, 156(4): 914-921.
    Jiang L, Bai C H, Liu Q M. Experimental Study on DDT Process in 3-Phase Suspensions of Gas/Solid Particle/Liquid Mist Mixture [J]. Explosion and Shock Waves, 2010, 30(2): 588-592. (in Chinese)
    蒋丽, 白春华, 刘庆明. 气/固/液三相混合物燃烧转爆轰过程实验研究 [J]. 爆炸与冲击, 2010, 30(2): 588-592.
  • 加载中
计量
  • 文章访问数:  8116
  • HTML全文浏览量:  386
  • PDF下载量:  671
出版历程
  • 收稿日期:  2010-03-22
  • 修回日期:  2010-05-06
  • 刊出日期:  2011-08-15

目录

/

返回文章
返回