起爆偏心对聚能装药射流成型过程及威力参量的影响

徐恒威 梁斌 刘俊新 卢永刚 欧小红

任贤达, 刘佳琼, 唐震, 武晓刚, 陈维毅. 不同退火温度下输电导线疲劳性能的实验分析[J]. 高压物理学报, 2019, 33(4): 045902. doi: 10.11858/gywlxb.20180566
引用本文: 徐恒威, 梁斌, 刘俊新, 卢永刚, 欧小红. 起爆偏心对聚能装药射流成型过程及威力参量的影响[J]. 高压物理学报, 2023, 37(1): 015102. doi: 10.11858/gywlxb.20220635
REN Xianda, LIU Jiaqiong, TANG Zhen, WU Xiaogang, CHEN Weiyi. Experimental Analysis of Fatigue Performance in Transmission Lines at Different Annealing Temperatures[J]. Chinese Journal of High Pressure Physics, 2019, 33(4): 045902. doi: 10.11858/gywlxb.20180566
Citation: XU Hengwei, LIANG Bin, LIU Junxin, LU Yonggang, OU Xiaohong. Effect of Initiation Eccentricity on Shaped Charge Jet Forming Process and Power Parameters[J]. Chinese Journal of High Pressure Physics, 2023, 37(1): 015102. doi: 10.11858/gywlxb.20220635

起爆偏心对聚能装药射流成型过程及威力参量的影响

doi: 10.11858/gywlxb.20220635
基金项目: 国家自然科学基金(11672278);国家重点研发计划(2017YFC0804600)
详细信息
    作者简介:

    徐恒威(1998-),男,硕士研究生,主要从事高效毁伤技术研究. E-mail:xuhw@mails.swust.edu.cn

    通讯作者:

    梁 斌(1976-),男,博士,研究员,主要从事常规武器研制及毁伤研究. E-mail:lb_110119@163.com

  • 中图分类号: O383; TJ410.33

Effect of Initiation Eccentricity on Shaped Charge Jet Forming Process and Power Parameters

  • 摘要: 为了研究起爆偏心对聚能射流的影响,运用有限元软件LS-DYNA模拟了不同起爆偏心量(0.025Dk~0.125DkDk为装药直径)下射流成型及其破甲过程,探究了药型罩非对称压垮程度、射流形态以及横向速度的变化规律,建立了理论模型以分析不同偏心量下射流横向速度分布情况,并基于正交试验设计理论和方差分析法揭示了各因素对评价指标影响程度的显著差异。结果表明:药型罩非对称压垮程度及射流横向速度均与偏心量呈正相关变化趋势。偏心量为0.025Dk时,射流侵彻深度仅下降0.7%;偏心量为0.050Dk时,侵彻深度下降突跃为12.4%;随着偏心量的增加,侵彻深度继续下降。此外,适当增大壁厚、罩顶装药高度可削弱起爆偏心对射流横向速度的影响。

     

  • 长直空间具有横截面尺寸小、长度尺寸远大于截面尺寸、泄爆面积有限等特点,常见于煤矿巷道、隧道洞库、管涵暗渠、城市综合管廊等。在煤矿采挖、隧道施工以及油气管道输运过程中存在燃气泄漏的可能,泄漏的燃气在长直空间内与空气混合后遇点火源极易引发爆炸,进而造成人员伤亡和财产损失。近年来相关事故时有发生,例如:2013年青岛排水暗渠油气爆炸,2015年成都洛带古镇隧道瓦斯爆炸,2016年重庆金山沟煤矿瓦斯爆炸等。燃气爆炸产生的压力波是造成人员及结构损害的主要原因之一[1-2],有学者发现特定条件下长直空间燃气爆炸存在明显的超压振荡现象[3-8],而超压振荡将引起荷载峰值及升压速率的明显升高[9-10]、流场湍流程度增强[11],并可能造成结构高温蠕变[3]以及引发结构共振[12],从而进一步加剧爆炸对人员和结构的危害。考虑到长直空间燃气爆炸事故可能造成的严重后果,开展爆炸超压及其振荡抑制措施研究十分必要。

    Xing等[12]利用长度为2 m的圆管进行了甲烷/空气爆炸实验,发现CH4的体积分数为9.5%~12.0%时,密闭管道爆炸出现了剧烈的压力振荡;朱传杰等[13-14]通过实验和数值模拟(AutoReaGas)研究了不同尺寸密闭方管内CH4/空气爆炸超压的振荡特性,认为超压的振荡是压力波在管道两端产生反射的结果;韦世豪等[15-16]进行了小尺寸密闭管道内汽油蒸气/空气爆炸实验,认为前驱波和压缩波在容器内的反射叠加是产生压力振荡的主要原因,而火焰和压力波耦合所形成的正反馈现象进一步加剧了振荡。研究人员在长直尾端泄爆空间内同样观察到了压力振荡现象[6-8];Xing等[12]发现随着泄爆面积的增加,压力振荡相比于密闭管道逐渐衰减,但另外出现了亥姆霍兹振荡;李毅[17]、路长等[18]通过实验发现气体流入流出的交替变化对压力的振荡具有影响;王亚磊等[19]在长0.5 m的尾端管道内开展了燃气爆炸实验,结果表明管内气体的快速泄放和回流造成了压力的振荡。

    目前应用较为广泛的抑爆措施主要包括水系抑爆及惰性气体抑爆,二者分别通过在可燃气体中加入一定的水雾/水蒸气和惰性气体对爆炸产生抑制。用于抑爆的惰性气体主要为氮气(N2)和二氧化碳(CO2[20-21],实验发现随着N2和CO2浓度的增加,惰性气体对甲烷(CH4)爆炸的抑制作用增强[22-23],CO2和N2对甲烷/空气混合气体的抑爆浓度分别为22%~23%和30%~37%[24-26]。李成兵[27]、Xie等[28]认为CO2的比热容更高以及CO2作为燃烧产物能够从化学角度对燃烧产生抑制是CO2抑爆性能优于N2的主要原因。水系抑爆主要包括水雾抑爆和水蒸气抑爆,秦文茜等[29]、康泉胜等[30]分别通过实验研究了超细水雾对瓦斯和丙烷爆炸的抑制效果,发现添加超细水雾能够有效降低爆炸超压及其升压速率,并使爆炸下限大幅提高,随着超细水雾添加量的增加,抑制效果增强。高旭亮[31]利用Fluent软件研究了超细水雾对CH4爆炸的抑制效果,得到了与上述实验相同的结论,此外还发现当施加量达到一定程度时,超细水雾甚至会完全抑制CH4的爆炸。刘丹等[32]通过实验研究了水蒸气对CH4爆炸超压及爆炸极限的影响,发现水蒸气同样能够对爆炸产生抑制,且随着水蒸气浓度的增加,爆炸超压及升压速率降低,爆炸极限范围缩小。

    研究人员对长直空间燃气爆炸超压的振荡现象进行了研究,但大多关心超压的振荡特性及其产生原因,鲜有关于超压振荡抑制措施的研究报道。虽然添加N2、CO2和水蒸气能够在一定程度上对燃气爆炸超压及其升压速率产生抑制,但其能否抑制超压的振荡仍有待深入探讨。本研究拟通过FLACS数值模拟软件,建立长直密闭和长直尾端泄爆管道模型,并利用已有实验数据对模型进行验证;基于已验证的数值模型,以化学当量比甲烷/空气混合气体为研究对象,分析不同浓度CO2、N2及水蒸气对长直空间燃气爆炸超压及升压时间的抑制作用,对比3种气体对超压振荡的抑制效果,探究不同气体抑爆效果差异的原因。

    甲烷/空气混合气体的爆炸是一个快速释能的化学过程。FLACS描述可压缩流体流动的数学模型包括质量守恒方程、动量守恒方程以及焓的输运方程[33]

    质量守恒方程:

    t(βvρ)+xj(βjvρuj)=˙mV
    (1)

    动量守恒方程:

    t(βvρui)+xj(βjρuiuj)=βvpxi+pxi(βjσij)+Fo,i+Fw,i+βv(ρρ0)gi
    (2)

    焓的输运方程:

    t(βvρh)+xj(βjρujh)=xj(βjμeffσhhxj)+βvdpdt+˙QV
    (3)

    式中:t为时间;βv为体积孔隙率;ρ为密度;ρ0为初始密度;m为质量;β为火焰模型中的转换系数;V为比容;h为比焓;σh为焓对应的Prandtl数,取0.7; μeff为有效黏度;Q为热量流率;p为压强;uiuj分别为i方向和j方向的平均速度;σij为应力张量分量;gii方向的重力加速度;Fw,iFo,i分别为容器壁和亚网格阻塞引起的流动阻力。

    FLACS的湍流模型为标准k-ε 模型[33],属于涡黏模型中的两方程模型,湍动能k的运输方程为

    (ρk)t+(ρkui)xi=xj[(μ+μtσk)kxj]+Gk+GbρεYM+Sk
    (4)

    耗散率 ε 的运输方程为

    (ρε)t+(ρεui)xi=xj[(μ+μtσε)εxj]+C1εεk(Gk+C3εGb)C2ερε2k+Sε
    (5)

    式中:μ为动力黏度;μt为湍动黏度;Gk是由平均速度梯度引起的湍动能k的产生项;Gb是由浮力引起的湍动能k的产生项;YM代表可压湍流中脉动扩张的贡献;σk为湍动能k对应的Prandtl数,取为1.0;σε为耗散率ε对应的Prandtl数,取为1.3;SkSε为用户定义的源项;C1εC2εC3ε为经验常数,FLACS根据Launder等[34]的研究结果,分别取1.44、1.92、0.80。

    FLACS燃烧模型由火焰模型和燃烧速度模型组成[33]。火焰模型采用β模型,通过定义扩散系数D和反应速率RF确定燃烧标准和反应速率在数值火焰区域的空间分布。燃烧速度取准层流燃烧速度和湍流燃烧速度中的较大值。

    建模过程中对FLACS管道数值模型做以下简化:(1)模型管道壁面为刚性,不考虑其变形和损坏;(2)假设爆炸过程为绝热爆炸,不考虑由管道壁面散热引起的能量损失;(3)忽略模型管道壁上的各种可关闭预留孔口、传感器以及管道连接部位等细节,忽略管道壁对压力波的阻碍作用,认为容器内壁光滑。

    考虑到长直密闭和尾端泄爆管道内燃气爆炸均可能存在超压的振荡现象,建立如图1所示的完全密闭和尾端泄爆两种管道模型,两种管道内部尺寸均为0.8 m × 0.8 m × 30 m,尾端泄爆管道泄爆口尺寸为0.8 m × 0.8 m,泄爆压力为0 kPa。点火点设置在管道端部中央,分别在每节管道侧壁距点火点1、4、7、···、28 m每隔3 m处布置压力测点,依次编为测点1~测点10。流速尺度(CFLV)和声速尺度(CFLC)分别取为0.5和5.0,初始温度为20 ℃,初始压强100 kPa,相对湍流强度和湍流长度尺度均为零,其余设置均采用默认值。

    图  1  管道模型及其网格划分
    Figure  1.  Tube model and its meshing

    密闭模型各边界均设置为EULER边界,尾端泄爆管道气体泄放到开敞空间,故将除z轴负方向(地面)外的其他方向均设置为PLAN_WAVE边界。计算域及网格尺寸设置如下:将管道内部设置为核心区,核心区分别采用0.10、0.05、0.03 m 3种尺寸的网格。密闭管道核心区以外0.2 m范围为拉伸区,尾端泄爆管道核心区外x轴正负方向3 m,y轴正方向15 m,y轴负方向0.5 m,z轴正方向10 m,z轴负方向0.2 m为拉伸区,拉伸系数为1.2。

    研究人员对FLACS模拟惰性气体抑爆的可靠性进行了验证,陈晓坤等[35]利用FLACS软件对多组分可燃气体添加CO2后的爆炸参数进行了模拟研究,并将模拟结果与已有实验数据进行对比,发现不同配比条件下的最大爆炸超压模拟值与实验值的相对误差均小于16%,说明利用 FLACS 模拟多组分可燃性气体爆炸可行。罗振敏等[36]在内径88 mm、长1600 mm的尾端泄爆管道内进行了CO2抑制甲烷/空气爆炸实验,并通过FLACS软件对实验进行了模拟,结果表明,除CO2添加量为18%时误差较大外,模拟得到的火焰传播时间和火焰传播速度均与实验值基本吻合且变化趋势一致,认为误差主要来源于模型的理想环境假设以及FLACS对小尺寸管道模拟的偏差。张印等[37]利用FLACS软件对20 L密闭容器内添加CO2的CH4/空气预混爆炸进行了数值模拟,发现最大爆炸压力的模拟值与理论计算值的相对误差在6.5%~10.5%之间,爆炸温度的模拟值与计算值的相对误差在3%~13%之间。

    为验证FLACS模拟大尺寸长直管道燃气爆炸的准确性,对Li等[38]开展的长直空间燃气爆炸实验进行了模拟,模拟工况见图2。该实验管道尺寸为0.8 m × 0.8 m × 30 m,由10节长度为3 m的钢制管道组成。工况A和工况B均在管道端部充入3 m长的9.5%甲烷/空气混合气体,点火位置位于管道端部中央,点火能量100 mJ。工况A中管道完全密闭,工况C中第4、7、10节管道顶部泄爆口(尺寸0.6 m × 0.6 m)打开,泄爆压力0 kPa。需要指出的是,工况A的数值模型与完全密闭管道模型(图1(a))相同;工况C的数值模型除泄爆口数量与分布尾端泄爆管道模型(图1(b))不同外,其他设置均与之相同。

    图  2  模型验证工况[38]
    Figure  2.  Cases of model verification[38]

    工况A和工况C的模拟结果与实验值的对比如图3所示,以靠近点火端的测点1为例,可以看到密闭和泄爆条件下FLACS均能够正确地模拟爆炸超压随时间的变化规律。对于密闭管道(工况A),网格尺寸为0.10 m时模拟值与实验值最接近,第1个超压峰值(65 kPa)与实验值完全相同,第2个超压峰值(110 kPa)的相对误差为14.6%。对于泄爆管道(工况C),网格尺寸为0.05和0.03 m时误差较小且接近,超压峰值的误差分别为9.5%、11.9%。造成模拟值与实验值差异的主要原因:实验管道壁为散热良好的钢制粗糙壁面,而数值模型做了光滑、绝热的假设。模拟结果表明,数值模型基本可靠,模拟误差在可接受范围内。综合考虑结果的准确性和计算效率,后续密闭管道和泄爆管道分别采用0.10和0.05 m尺寸的网格进行计算。

    图  3  模型验证(测点1)
    Figure  3.  Numerical model verification (pressure sensor 1)

    为研究添加惰性气体及水蒸气对长直空间内燃气爆炸超压及其振荡的抑制效果,利用FLACS软件对添加了不同体积分数N2、CO2和水蒸气的化学当量比CH4/空气混合气体进行数值模拟,添加气体的体积分数如表1所示。充气体积如图4所示,在密闭管道内充满混合气体,在尾端泄爆管道点火端充入3 m长的混合气体。由于密闭管道尾端[14]和泄爆管道点火端[19]的爆炸超压及其振荡相较管道其他位置更明显,如未特别说明,本研究所讨论的密闭管道内的超压时程曲线均取自管道尾端测点(测点10),所讨论的尾端泄爆管道内的超压时程曲线均取自管道点火端测点(测点1)。

    表  1  添加气体的体积分数(密闭/泄爆)
    Table  1.  Volume fraction of added gas (closed/end-vented) % 
    CO2N2Water vapor (H2O) CO2N2Water vapor (H2O)
    000 202520
    454 243024
    8108 253126
    121512 27
    162016
    下载: 导出CSV 
    | 显示表格
    图  4  混合气体体积示意图
    Figure  4.  Diagram of premixed gas volume
    2.1.1   CO2的抑爆效果

    通过模拟发现,密闭管道内CO2的体积分数超过24%时不能点火,爆炸被完全抑制,与前人实验结果[24-26]非常接近。图5为添加不同体积分数CO2预混气体的爆炸超压时程曲线。可以发现,随着CO2体积分数的增加,超压的振荡现象逐渐减弱,振荡持续时间逐渐变短,当CO2添加量达到16%时,振荡完全消失。CO2对超压振荡的抑制可归因于其对混合气体燃烧速率的降低作用[27, 39]:首先,添加CO2稀释了甲烷及氧气,使反应速率降低,燃烧产生的热量减少;其次,CO2将吸收燃烧产生的热量使反应区温度降低,抑制链式反应的进行;最后,CO2作为燃烧产物之一将与燃烧区自由基反应,使自由基浓度降低,爆炸链传递速度减缓。密闭管道内超压的振荡主要是压力波来回反射传播的结果[13-16],燃烧速率的减缓将造成压力波强度(大小及升压速率)下降,使管内压力能够及时分布均匀,所以最终未表现出压力振荡现象。

    图  5  不同CO2体积分数条件下的超压时程曲线
    Figure  5.  Overpressure-time history curves of CO2 with different volume fractions

    图6所示,当CO2体积分数(φ)从零逐渐上升到24%时,升压时间(达到压力峰值的时间)从1.28 s逐渐增加到19.60 s,最终爆炸超压由785 kPa呈线性降低到590 kPa,即CO2的体积分数每增加10%,升压时间平均增加7.63 s,爆炸超压下降81 kPa。升压时间的增加是燃烧速率降低的结果,爆炸超压的降低则与添加CO2后管内燃料质量的降低有关。需要指出的是,未添加CO2时,超压峰值(897 kPa)受压力振荡的影响高于最终爆炸超压,但当CO2的体积分数达到4%时,受CO2抑爆作用的影响,超压的振荡减弱并不再使超压峰值升高,此后管内最终压力即为爆炸最大超压。

    图  6  最终超压及升压时间随CO2体积分数的变化
    Figure  6.  Change of final overpressure and pressure rising time with CO2 volume fraction

    超压振幅随CO2体积分数的变化见图7。可以看到:未添加CO2时压力振幅存在一个明显的上升阶段(0~0.96 s),该阶段点火后气体燃烧产生的压力波在管道内的往复反射形成了超压的振荡;燃烧与压力波的耦合作用使压力振幅不断升高[15-16],并在0.96 s达到最大值(337 kPa);0.96 s之后压力振幅进入下降阶段,此时燃料被大量消耗,预混气体浓度减小,燃烧速率降低,燃烧对压力波的增强逐渐小于压力波的摩擦损耗,振幅呈指数函数下降。需要指出的是,数值模型假设管道壁绝热,所以不考虑散热造成的压力振幅降低。随着CO2添加量的增加,超压振幅逐渐减小,振幅随时间的变化由升高和降低两个阶段逐渐变为仅有降低阶段,这是因为燃烧产生的压力波强度随着CO2体积分数的升高而降低,当CO2的体积分数达到一定值(约8%),燃烧对压力波的增强作用将小于气体摩擦损耗对压力波的衰减作用,振幅上升阶段消失。

    图  7  不同CO2体积分数下超压振幅随时间的变化
    Figure  7.  Variation of overpressure amplitude of different volume fractions of CO2 with time

    图8为不同CO2体积分数条件下超压振荡频率随时间的变化。可以看到,振荡频率随时间逐渐升高,并最终趋于一定值。此外,可以发现当CO2添加量增加时,超压振荡频率的上升速率和最大超压振荡频率降低。超压振荡归因于压力波的往复反射传播,在管道长度不变的情况下,超压振荡频率的大小主要取决于压力波的传播速度。在管道模型的绝热假设条件下,随着燃烧的进行,管内温度逐渐上升,并在燃烧结束后维持在一定值,而压力波传播速度与气体温度的1/2次方正相关[40],故超压振荡频率表现为先上升后维持的变化规律。另一方面,随着CO2体积分数的增加,甲烷的质量减小,燃烧产生的热量减少,造成管内温度降低,故最大超压振荡频率随CO2体积分数的增加而降低。

    图  8  不同CO2体积分数下超压振荡频率随时间的变化
    Figure  8.  Variation of overpressure oscillation frequency of different volume fractions of CO2 with time
    2.1.2   N2的抑爆效果

    图9为添加不同体积分数N2预混气体的爆炸超压时程曲线。模拟发现:N2的体积分数等于30%时,管内压力在10 s内均无明显变化;当N2的体积分数超过30%时,无法点火,可以认为模拟得到的N2完全抑爆浓度为30%,与已有实验结果[24-26]基本一致。图10显示了最终爆炸超压及升压时间随N2体积分数的变化。当N2的体积分数由零增加到25%时,升压时间从1.28 s逐渐增长到11.58 s,最终爆炸超压由785 kPa线性降低到667 kPa,即N2的体积分数每增加10%,升压时间平均增加4.12 s,爆炸超压下降47 kPa。随着N2体积分数的升高,超压的振荡逐渐减弱,振荡持续时间逐渐变短,N2的体积分数为20%时,振荡被完全抑制。

    图  9  不同N2体积分数条件下的超压时程曲线
    Figure  9.  Overpressure-time history curves of N2 with different volume fractions
    图  10  最终超压及升压时间随N2体积分数的变化
    Figure  10.  Change of final overpressure and pressure rising time with N2 volume fraction

    图11为超压振幅随N2体积分数的变化。随着N2添加量的增加,超压振幅逐渐降低,N2的体积分数达到10%时,振幅上升阶段消失,高于CO2(8%)。图12为超压振荡频率随N2体积分数的变化。随着N2添加量的增加,管内温度降低,压力波传播速度降低,超压振荡频率减小。N2主要通过稀释甲烷和氧气浓度抑制爆炸,其次是吸收燃烧产生的热量,相比CO2缺少了化学反应抑制一项,故N2对爆炸超压及其振荡的抑制作用弱于CO2[27, 39]

    图  11  不同N2体积分数下超压振幅随时间的变化
    Figure  11.  Variation of overpressure amplitude of different volume fractions of N2 with time
    图  12  不同N2体积分数下超压振荡频率随时间的变化
    Figure  12.  Variation of overpressure oscillation frequency of different volume fractions of N2 with time
    2.1.3   水蒸气的抑爆效果

    图13为不同水蒸气体积分数下甲烷/空气混合气体的爆炸超压时程曲线。可见,当水蒸气的体积分数达到16%时,超压振荡被完全抑制。此外,当水蒸气的体积分数达到26%时,管内压力在16 s内几乎未上升;当水蒸气的体积分数超过26%时,混合气体无法点火,由此可知模拟得到的水蒸气完全抑爆浓度为26%。

    图  13  不同水蒸气体积分数条件下的超压时程曲线
    Figure  13.  Overpressure-time history curves of N2 with different volume fractions

    图14所示,当水蒸气的体积分数由零增加到24%时,升压时间从1.28 s逐渐增长到13.79 s,最终爆炸超压由785 kPa线性降低到629 kPa,即水蒸气的体积分数每增加10%,升压时间平均增加5.21 s,爆炸超压下降65 kPa。

    图  14  最终超压及升压时间随水蒸气体积分数的变化
    Figure  14.  Change of final overpressure and pressure rising time with water vapor volume fraction

    图15图16分别显示了超压振幅和超压振荡频率随水蒸气体积分数的变化。与添加CO2和N2类似,当水蒸气的体积分数约为8%时,超压振幅的上升阶段消失,随着水蒸气体积分数的继续增大,最大超压振幅减小,超压振荡频率降低。

    图  15  不同水蒸气体积分数下超压振幅随时间的变化
    Figure  15.  Variation of overpressure amplitude of different volume fractions of water vapor with time
    图  16  不同水蒸气体积分数下超压振荡频率随时间的变化
    Figure  16.  Variation of overpressure oscillation frequency of different volume fractions of water vapor with time
    2.2.1   CO2的抑爆效果

    图17为不同CO2体积分数条件下混合气体爆炸超压时程曲线。可以发现,当尾端泄爆管道内CO2的体积分数达到24%时,混合气体仅发生缓慢燃烧,管内压力不发生明显变化;当CO2的体积分数为25%时,气体将不能点火,可以认为,尾端泄爆管道CO2的完全抑爆浓度为24%,与密闭管道相同。

    图  17  不同CO2体积分数条件下的超压时程曲线
    Figure  17.  Overpressure-time history curves of CO2 with different volume fractions

    图18所示,随着CO2的体积分数由零增长到24%,升压时间从0.24 s逐渐增长到0.71 s,超压峰值由57 kPa接近线性降低到0 kPa,即CO2的体积分数每增加10%,爆炸超压下降约24 kPa。

    图  18  最终超压和升压时间随CO2体积分数的变化
    Figure  18.  Change of final overpressure and pressure rising time with CO2 volume fraction

    图19为爆炸超压振幅随时间的变化。可见,不同CO2体积分数条件下混合气体的压力振幅均随时间呈指数函数下降。泄爆管道内压力的振荡类似于阻尼振动[41],振荡可归因于管内气体过度泄放引起气体往复的泄出与回流[19],因此可以认为,压力振幅的下降与振荡过程中气体动能被耗散有关。此外,随着CO2的体积分数由零增加到20%,最大压力振幅由90 kPa减小为7 kPa,这是CO2导致的混合气体燃烧速率降低、爆炸超压和管内气体流速减小的结果。

    图  19  不同CO2体积分数下超压振幅随时间的变化
    Figure  19.  Variation of overpressure amplitude of different volume fractions of CO2 with time

    图20给出了压力振荡频率随时间的变化。可以发现,振荡频率随CO2体积分数的增加而降低,可解释为:CO2的体积分数越大,燃烧速率越低,管内气体通过燃烧获取的速度越小,气体往复运动一个周期所需时间更长,最终导致超压振荡频率降低。不考虑由超压时程曲线得到压力振荡频率的误差,可以发现,各工况下的压力振荡频率几乎不随时间而明显上升或下降,如果认为尾端泄爆管道内压力的振荡为阻尼振动,则压力振荡频率的不变可能与阻尼振动的等时性[41]有关。

    图  20  不同CO2体积分数下超压振荡频率随时间的变化
    Figure  20.  Variation of overpressure oscillation frequency of different volume fractions of CO2 with time
    2.2.2   N2的抑爆效果

    图21为不同N2体积分数下混合气体爆炸超压时程曲线。可以发现,当N2的体积分数达到30%时,混合气体仅发生缓慢燃烧,管内压力无明显变化;当N2的体积分数为31%时,爆炸被完全抑制,说明尾端泄爆管道内N2的完全抑爆浓度与密闭管道相同,均为30%。

    图  21  不同N2体积分数条件下的超压时程曲线
    Figure  21.  Overpressure-time history curves of N2 with different volume fractions

    图22所示,当N2的体积分数由零增加到30%,升压时间由0.24 s增长到0.79 s,超压峰值由57 kPa接近线性降低到0 kPa,即CO2的体积分数每增加10%,爆炸超压下降20 kPa。

    图  22  最终超压及升压时间随N2体积分数的变化
    Figure  22.  Change of final overpressure and pressure rising time with N2 volume fraction

    压力振幅和压力振荡频率随时间的变化分别如图23图24所示。可以看到,随着N2体积分数的升高,超压振幅和振荡频率均逐渐减小,说明N2对尾端泄爆管道内爆炸超压的振荡效应具有一定的抑制作用。

    图  23  不同N2体积分数下超压振幅随时间的变化
    Figure  23.  Variation of overpressure amplitude of different volume fractions of N2 with time
    图  24  不同N2体积分数下超压振荡频率随时间的变化
    Figure  24.  Variation of overpressure oscillation frequency of different volume fractions of N2 with time
    2.2.3   水蒸气的抑爆效果

    图25为不同水蒸气体积分数下甲烷/空气混合气体的爆炸超压时程曲线。可以发现,尾端泄爆管道内水蒸气完全抑制爆炸的体积分数与密闭管道相同,均为26%。

    图  25  不同水蒸气体积分数条件下的超压时程曲线
    Figure  25.  Overpressure-time history curves of water vapor with different volume fractions

    图26可以看出,随着水蒸气体积分数由零增加到24%,升压时间从0.24 s逐渐增长到0.68 s,最终爆炸超压由57 kPa接近线性降低到0 kPa,即水蒸气的体积分数每增加10%,爆炸超压下降约25 kPa。

    图  26  最终超压及升压时间随水蒸气体积分数的变化
    Figure  26.  Change of final overpressure and pressure rising time with water vapor volume fraction

    图27图28分别显示了不同水蒸气体积分数下超压振幅和振荡频率随时间的变化。可以看出,振幅和频率均随着水蒸气含量的增加而减小,说明添加水蒸气同样能够有效抑制尾端泄爆管道内甲烷爆炸超压的振荡效应。

    图  27  不同水蒸气体积分数下超压振幅随时间的变化
    Figure  27.  Variation of overpressure amplitude of different volume fractions of water vapor with time
    图  28  不同水蒸气体积分数下超压振荡频率随时间的变化
    Figure  28.  Variation of overpressure oscillation frequency of different volume fractions of water vapor with time

    由2.1、2.2节可知,CO2、水蒸气和N2的完全抑爆浓度分别为25%、26%、30%,且泄爆条件对完全抑爆浓度没有影响。密闭条件下完全抑制超压振荡时CO2、水蒸气和N2的体积分数分别为16%、16%、20%,尾端泄爆条件下压力振荡被完全抑制时的体积分数则与完全抑爆浓度相同,为25%、26%、30%。图29图30分别为添加CO2、水蒸气和N2对长直空间内甲烷爆炸超压抑制作用的对比。可以发现,当添加气体浓度相同时,密闭和尾端泄爆管道均是添加CO2的爆炸超压及其振幅下降最明显,升压时间增加最长,添加水蒸气次之,添加N2最差,说明CO2对爆炸超压及其振荡的抑制作用最强,N2最弱。

    图  29  密闭管道内CO2、N2和水蒸气对爆炸的抑制作用对比
    Figure  29.  Comparison of suppression of CO2, N2 and water vapor on the explosion in the closed tube
    图  30  尾端泄爆管道内CO2、N2和水蒸气对爆炸的抑制作用对比
    Figure  30.  Comparison of suppression of CO2, N2 and water vapor on the explosion in the end-vented tube

    3种气体抑爆效果的差异可解释为[27, 39]:N2仅能从物理角度(稀释、吸热)对超压及其振荡产生抑制,而CO2与水蒸气作为甲烷爆炸的燃烧产物,除物理抑爆外,还能够与燃烧区自由基反应,降低自由基浓度,抑制爆炸链传递速度,进而从化学角度对爆炸产生抑制,故抑爆效果优于N2。水蒸气和CO2抑爆效果的差异可能与二者比热容不同及参与链式反应的种类和程度不同有关。

    利用FLACS软件建立了长直(密闭/尾端泄爆)空间管道模型,并对模型进行了验证。基于已验证的模型研究了添加CO2、N2和水蒸气的化学当量比CH4/空气混合气体的爆炸情况,讨论了添加不同浓度气体对爆炸超压及其振荡的抑制作用,并对比了3种气体的抑爆效果。

    (1)管道内气体爆炸超压随添加气体体积分数的增加而降低,CO2、水蒸气和N2的体积分数每增加10%,密闭管道气体爆炸的最终超压分别下降81、65、47 kPa,尾端泄爆管道超压峰值分别下降24、25、20 kPa,3种气体的完全抑爆浓度分别为25%、26%和30%(体积分数)。

    (2)CO2、水蒸气和N2对长直空间燃气爆炸超压的振荡具有抑制作用,压力振幅和压力振荡频率均随添加气体体积分数的增加而减小,密闭管道内3种气体完全抑制超压振荡时的体积分数分别为16%、16%、20%,尾端泄爆管道分别为25%、26%、30%。

    (3)CO2对爆炸超压及其振荡的抑制效果最好,水蒸气次之,N2最弱。这是因为N2仅能从物理角度对爆炸产生抑制,而CO2与水蒸气作为甲烷爆炸的燃烧产物,还能通过参与链式反应而抑制爆炸。水蒸气和CO2抑爆效果的差异可能与二者比热容以及参与链式反应的种类和程度不同有关。

  • 图  药型罩不对称压垮几何示意图

    Figure  1.  Geometric representation of the asymmetric collapse of shaped charge

    图  数值模型

    Figure  2.  Numerical simulation model

    图  示踪点位置

    Figure  3.  Location of static tracer points

    图  脉冲X射线摄影试验

    Figure  4.  X-ray photography test

    图  试验结果与模拟结果对比(30 μs)

    Figure  5.  Comparison of test and numerical simulaiton results (30 μs)

    图  药型罩两侧横向压垮速度差

    Figure  6.  Lateral collapse velocity difference between two sides of liner

    图  不同偏心量下的射流横向速度

    Figure  7.  Lateral velocity of jet under different initiation offsets

    图  射流速度云图

    Figure  8.  Velocity nephogram of shaped charge jet

    图  示踪点运动轨迹

    Figure  9.  Trace line of tracer points at liner surface

    图  10  罩材沿轴线方向的射流速度分布曲线

    Figure  10.  Distribution curve of the jet velocity of liner material along the axis direction

    图  11  数值模拟结果与理论计算值的对比

    Figure  11.  Comparison of numerical simulation results and theoretical calculated values

    图  12  不同偏心量下射流的破甲威力

    Figure  12.  Armor-breaking power of jet under different offsets

    图  13  不同因素水平对射流轴向速度的影响

    Figure  13.  Effects of factor levels on axial velocity of jet

    图  14  不同因素水平对射流有效长度的影响

    Figure  14.  Effects of factor levels on the effective length of jet

    图  15  不同因素水平对射流横向速度的影响

    Figure  15.  Effects of factor levels on lateral velocity of jet

    图  16  经优化后得到的射流形态

    Figure  16.  Jet shape after optimization

    表  1  聚能装药的几何尺寸

    Table  1.   Structural parameters of shaped charge

    Charge diameterExplosive height above linerCone angle/(°)Thickness of linerRadius of liner curvatureBrust height
    Dk0.5Dk600.020Dk0.037 5Dk3Dk
    下载: 导出CSV

    表  2  炸药模型及状态方程参数[15]

    Table  2.   Parameters of explosive and equation of state[15]

    ρ/(g·cm−3)D/(m·s−1)pCJ/GPaA0/GPaB0/GPaR1R2
    1.798 73033.8795.220.24.641.29
    下载: 导出CSV

    表  3  药型罩模型及状态方程参数[15]

    Table  3.   Parameters of liner and equation of state[15]

    ρ/(g·cm−3)G0/GPaσy/GPaC0/(km·s−1) s1 s2γ0
    8.9647.70.643.941.4902
    下载: 导出CSV

    表  4  射流成型参数对比(30 μs)

    Table  4.   Comparison of jet formation parameters (30 μs)

    Head jet Tail jet Length Jet deviation
    vh,sim/(m·s−1)vh,exp/(m·s−1)δ/%vt,sim/(m·s−1)vt,exp/(m·s−1)δ/%Lsim/cmLexp/cmδ/%ξsim/(°)ξexp/(°)δ/%
    6 4046 600−2.9 622745−16.5 13.4513.60−1.1 1.251.158.6
    下载: 导出CSV

    表  5  正交试验设计

    Table  5.   Orthogonal test design

    Project Factor Evaluation indicator
    123va/(m·s−1)vb/(m·s−1)L/cm
    1Ⅰ(0.012 5Dk)Ⅰ(0.020Dk)Ⅰ(0.5Dk) 6 38445.18.68
    2Ⅰ(0.012 5Dk)Ⅱ(0.025Dk)Ⅱ(0.8Dk)6 26026.29.02
    3Ⅰ(0.012 5Dk)Ⅲ(0.030Dk)Ⅲ(1.0Dk)6 13017.98.93
    4Ⅱ(0.025Dk)Ⅰ(0.020Dk)Ⅱ(0.8Dk)6 55478.29.18
    5Ⅱ(0.025Dk)Ⅱ(0.025Dk)Ⅲ(1.0Dk)6 40157.49.02
    6Ⅱ(0.025Dk)Ⅲ(0.03Dk)Ⅰ(0.5Dk)5 73875.58.76
    7Ⅲ(0.05Dk)Ⅰ(0.020Dk)Ⅲ(1.0Dk)6 710154.19.11
    8Ⅲ(0.05Dk)Ⅱ(0.025Dk)Ⅰ(0.5Dk)6 023187.08.54
    9Ⅲ(0.05Dk)Ⅲ(0.030Dk)Ⅱ(0.8Dk)5 993135.48.94
    下载: 导出CSV

    表  6  计算结果

    Table  6.   Calculation results

    Factorva/(m·s−1) L/cm vb/(m·s−1)
    Rj Rj Rj
    K118 77418 69318 726 26.6326.9626.59 89.2211.1476.5
    K219 64718 68417 862 26.9726.5826.63 277.4270.6228.8
    K318 14418 80719 242 25.9827.1427.06 307.6239.8229.4
    k1 6 258 6 231 6 242 27 8.878.988.850.13 29.770.4158.8129.1
    k2 6 549 6 228 5 954595 8.998.868.880.13 92.590.276.316.2
    k3 6 048 6 269 6 414366 8.669.049.020.38 102.579.976.526.0
    下载: 导出CSV
  • [1] 王钰婷, 黄正祥, 贾鑫, 等. 椭圆形截面聚能装药射流成型及侵彻特性 [J]. 含能材料, 2021, 29(2): 96–106. doi: 10.11943/CJEM2020276

    WANG Y T, HUANG Z X, JIA X, et al. Jet formation and penetration characteristics of shaped charge with elliptical cross-section [J]. Chinese Journal of Energetic Materials, 2021, 29(2): 96–106. doi: 10.11943/CJEM2020276
    [2] AYISIT O, CORUH M M. Investigation of off-axis initiation of long L/D hemispherical shaped charge warheads [J]. Procedia Engineering, 2013, 58: 487–495. doi: 10.1016/j.proeng.2013.05.056
    [3] 李如江, 张晋红, 王建波, 等. 偏心起爆对聚能射流形成的影响 [J]. 弹箭与制导学报, 2012, 32(1): 85–88. doi: 10.3969/j.issn.1673-9728.2012.01.025

    LI R J, ZHANG J H, WANG J B, et al. The effect of eccentric initiation on shaped charge jet formation [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2012, 32(1): 85–88. doi: 10.3969/j.issn.1673-9728.2012.01.025
    [4] AYISIT O. The influence of asymmetries in shaped charge performance [J]. International Journal of Impact Engineering, 2008, 35(12): 1399–1404. doi: 10.1016/j.ijimpeng.2008.07.027
    [5] 刘健峰, 龙源, 纪冲, 等. 含偏心起爆对EFP战斗部飞行特性的影响 [J]. 爆炸与冲击, 2015, 35(3): 335–342. doi: 10.11883/1001-1455-(2015)03-0335-08

    LIU J F, LONG Y, JI C, et al. Effect of eccentric initiation on the flight characteristics and ballistic dispersion of EFP [J]. Explosion and Shock Waves, 2015, 35(3): 335–342. doi: 10.11883/1001-1455-(2015)03-0335-08
    [6] PACK D C, CURTIS J P. On the effect of asymmetries on the jet from a linear shaped charge [J]. Journal of Applied Physics, 1990, 67(11): 6701–6704. doi: 10.1063/1.345107
    [7] 秦承森, 石艺娜, 冯其京, 等. 非对称射流形成的几何理论(一)——几何封闭条件 [J]. 爆炸与冲击, 2007, 27(6): 493–500. doi: 10.3321/j.issn:1001-1455.2007.06.003

    QIN C S, SHI Y N, FENG Q J, et al. The geometrical theory for the formation of asymmetric jet (Ⅰ): the geometrical close condition [J]. Explosion and Shock Waves, 2007, 27(6): 493–500. doi: 10.3321/j.issn:1001-1455.2007.06.003
    [8] BROWN J, SOFTLEY I D, EDWARDS P. Experimental study of shaped charges with built-in asymmetries [J]. Propellants, Explosives, Pyrotechnics, 1993, 18(5): 255–258. doi: 10.1002/prep.19930180504
    [9] BROWN J, CURTIS J P, COOK D D. The formation of jets from shaped charges in the presence of asymmetry [J]. Journal of Applied Physics, 1992, 72(6): 2136–2143. doi: 10.1063/1.352329
    [10] FEDOROV S V. Numerical simulation of the formation of shaped-charge jets from hemispherical liners of degressive thickness [J]. Combustion, Explosion, and Shock Waves, 2016, 52(5): 600–612. doi: 10.1134/S0010508216050117
    [11] 贾伟, 吴国东, 王志军, 等. 制造误差对金属射流威力影响的仿真研究 [J]. 弹箭与制导学报, 2011, 31(4): 99–101, 132. doi: 10.3969/j.issn.1673-9728.2011.04.029

    JIA W, WU G D, WANG Z J, et al. Strength simulation study of metallic jet effected by technical errors [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2011, 31(4): 99–101, 132. doi: 10.3969/j.issn.1673-9728.2011.04.029
    [12] 刘建荣, 张国伟, 徐立新, 等. 药型罩加工精度对破甲战斗部威力影响的研究 [J]. 弹箭与制导学报, 2012, 32(3): 114–117. doi: 10.3969/j.issn.1673-9728.2012.03.032

    LIU J R, ZHANG G W, XU L X, et al. The research of influence of liner’s precision on the power of armor-penetrating warhead [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2012, 32(3): 114–117. doi: 10.3969/j.issn.1673-9728.2012.03.032
    [13] PROSKURYAKOV E V, SOROKIN M V, FOMIN V M. Rebounding of a shaped-charge jet [J]. Journal of Applied Mechanics and Technical Physics, 2007, 48(5): 633–635. doi: 10.1007/s10808-007-0080-1
    [14] BROWN J, EDWARDS P J, LEE P R. Studies of shaped charges with built-in asymmetries. Part Ⅱ: modelling [J]. Propellants, Explosives, Pyrotechnics, 1996, 21(2): 59–63. doi: 10.1002/prep.19960210202
    [15] 时党勇, 李裕春, 张胜民. 基于ANSYS/LS-DYNA 8.1进行显式动力分析 [M]. 北京: 清华大学出版社, 2005: 295−296.

    SHI D Y, LI Y C, ZHANG S M. Explicit dynamic analysis based on ANSYS/LS-DYNA 8.1 [M]. Beijing: Tsinghua University Press, 2005: 295–296.
    [16] 黄正祥. 聚能装药理论与实践 [M]. 北京: 北京理工大学出版社, 2014.

    HUANG Z X. Theory and practice of shaped charge [M]. Beijing: Beijing Institute of Technology Press, 2014.
    [17] 侯秀成, 蒋建伟, 陈智刚. 有效射流结构模式的数值模拟 [J]. 爆炸与冲击, 2014, 34(1): 35–40. doi: 10.3969/j.issn.1001-1455.2014.01.007

    HOU X C, JIANG J W, CHEN Z G. Numerical simulation on structure modules of effective jet [J]. Explosion and Shock Waves, 2014, 34(1): 35–40. doi: 10.3969/j.issn.1001-1455.2014.01.007
    [18] 李磊, 马宏昊, 沈兆武. 基于正交设计方法的双锥罩结构优化设计 [J]. 爆炸与冲击, 2013, 33(6): 567–573. doi: 10.3969/j.issn.1001-1455.2013.06.002

    LI L, MA H H, SHEN Z W. Optimal design of biconical liner structure based on orthogonal design method [J]. Explosion and Shock Waves, 2013, 33(6): 567–573. doi: 10.3969/j.issn.1001-1455.2013.06.002
  • 加载中
图(16) / 表(6)
计量
  • 文章访问数:  246
  • HTML全文浏览量:  152
  • PDF下载量:  34
出版历程
  • 收稿日期:  2022-07-21
  • 修回日期:  2022-08-18
  • 录用日期:  2023-01-16
  • 网络出版日期:  2023-02-28
  • 刊出日期:  2023-02-05

目录

/

返回文章
返回