基于物质点法的射弹侵彻双层隔板充液结构的数值模拟

谢桂兰 宋慕清 龚曙光 侯昆 左立来 肖芳昱

雷慧茹, 张立宏. 高压下ReN2的弹性性质[J]. 高压物理学报, 2019, 33(4): 042401. doi: 10.11858/gywlxb.20180647
引用本文: 谢桂兰, 宋慕清, 龚曙光, 侯昆, 左立来, 肖芳昱. 基于物质点法的射弹侵彻双层隔板充液结构的数值模拟[J]. 高压物理学报, 2023, 37(1): 015101. doi: 10.11858/gywlxb.20220602
LEI Huiru, ZHANG Lihong. Elastic Properties of ReN2 under High Pressure[J]. Chinese Journal of High Pressure Physics, 2019, 33(4): 042401. doi: 10.11858/gywlxb.20180647
Citation: XIE Guilan, SONG Muqing, GONG Shuguang, HOU Kun, ZUO Lilai, XIAO Fangyu. Numerical Simulation of Projectile Penetrating Double-Layer Plate Liquid-Filled Structure Based on Material Point Method[J]. Chinese Journal of High Pressure Physics, 2023, 37(1): 015101. doi: 10.11858/gywlxb.20220602

基于物质点法的射弹侵彻双层隔板充液结构的数值模拟

doi: 10.11858/gywlxb.20220602
基金项目: 国家自然科学基金(51475403)
详细信息
    作者简介:

    谢桂兰(1966-),女,博士,教授,主要从事材料成形过程数值模拟研究. E-mail:xieguilan@xtu.edu.cn

    通讯作者:

    龚曙光(1964-),男,博士,教授,主要从事CAE技术的理论与应用研究. E-mail:gongsg@xtu.edu.cn

  • 中图分类号: O385; TB122

Numerical Simulation of Projectile Penetrating Double-Layer Plate Liquid-Filled Structure Based on Material Point Method

  • 摘要: 在充液结构内设置双层隔板空气夹层可以有效减弱水锤效应带来的危害。为了研究双层隔板间距及位置对水锤效应的影响机理,通过已发表的试验数据对数值模拟方法进行验证后,基于物质点法开展了射弹侵彻双层隔板充液结构的数值模拟,分析了空穴形成过程、射弹剩余速度、固定点液体压力峰值、出入射壁及双隔板壁变形等情况。结果表明:随着双层隔板间距增大,充液结构的变形程度呈先减轻后加重的趋势;双层隔板位置越靠近入射壁,对压力冲击波传导的阻碍作用越强,充液结构的抗侵彻能力越好。

     

  • 随着社会高速发展,特别是我国“双碳目标”的提出,对清洁能源的需求与日俱增,氢能作为一种环保、高效的新能源受到广泛关注。考虑到氢气管道建设成本高,利用现有的天然气管道,向甲烷中掺入一定量的氢气进而提升混合燃料热值,已成为一种经济可行的解决方案[1]。值得关注的是,近年来天然气输运事故屡有发生,若向其掺入氢气,有可能进一步加剧爆炸的危险性。倪靖等[2]讨论了不同掺氢比对甲烷-氧气爆轰特性的影响,发现掺氢后能够提高爆轰波的传播速度和爆轰敏感性。余明高等[3]通过实验研究了障碍物对甲烷-氢气爆炸特性的影响,发现最大爆炸超压和火焰传播速度随着障碍物阻塞率以及氢气体积分数的增大而增大。Yu等[4]研究了掺氢对甲烷-空气预混火焰传播特性的影响,发现随着氢气含量的增加,火焰前沿速度和爆炸超压显著升高。

    鉴于可燃气体在运输过程中的危险性,大量学者对可燃气体的抑爆问题开展了研究。Li等[5]发现CO2对甲烷爆炸的抑制效果好于N2。Luo等[6]研究BC粉对氢气-甲烷-空气预混气的抑制后发现,BC粉对掺氢比低的预混气体有较好的抑制作用。陈鹏等[7]发现当金属丝网层数大于3时,甲烷-空气预混气爆炸火焰经过金属丝网时会发生淬熄。徐海顺等[8]通过研究发现铝镁合金泡沫对甲烷-空气预混气体爆炸有较好的抑制作用,材料对传播火焰的影响机制主要体现在湍流促进和冷却抑制两方面。袁必和等[9] 研究了新型多孔聚丙烯复合材料对瓦斯爆炸的影响,发现多孔材料的抑爆性能受填充位置、材料内径以及填充长度的影响。An等[10]通过对比球形非金属材料和铝合金材料对可燃气体的抑爆结果,发现由于球形非金属材料具有较高的结构强度,因而抑爆效果较好。邵继伟等[11]研究多孔材料对可燃气体的抑爆效果后发现,组合型多孔材料在密闭容器管道系统内的抑爆效果更为突出。综合以上研究发现,非金属多孔材料对抑制甲烷爆炸具有良好的效果,但甲烷掺氢的爆炸特性较甲烷出现明显改变,并且前人对非金属多孔材料在甲烷掺氢抑爆方面的研究较少。传统的非金属材料存在易燃、爆炸后熔融物易黏附管道内壁等缺陷,严重限制了多孔非金属材料在可燃气体抑爆中的实际应用。球形抑爆材料中空多孔结构能够有效增大比表面积,扩大热损失,并且耐高温耐火焰,具有良好的阻火隔爆性质[12-13]

    本研究通过实验探究掺氢比对甲烷掺氢爆炸特性的影响,比较单一球形多孔非金属材料和组合球形多孔非金属材料的抑爆效能,并通过改变材料填充长度,分析最佳阻火抑爆性能的搭配参数,为球形多孔非金属材料在混合燃料阻隔防爆领域的应用提供理论支撑与实验依据。

    为研究掺氢比对甲烷掺氢预混气体爆炸的影响并测试球形多孔材料的抑爆性能,自主搭建了气体管道爆炸实验平台,如图1所示。实验平台主要包括爆炸实验管道、点火系统、配气系统、同步控制、高速摄像系统(拍摄速度4 000帧每秒)以及数据采集系统。

    图  1  实验装置示意图
    Figure  1.  Schematic diagram of the experiment setup

    实验管道由3节内径60 mm、管壁厚度10 mm、外径80 mm、长度2 m的水平圆形管道以及1节内径60 mm、长0.5 m、配备可视窗的方形管道组成,各个管道之间通过法兰-螺栓连接,实验管道为封闭状态。实验管道共安装2台压力传感器(CYG508微型高压传感器),分别距离管道左端1.2和6.3 m。点火系统(距离管道右端0.5 m处)由高压点火器以及2根钨丝组成,实验中点火能采用20 J。配气系统通过气瓶连接配气仪进行配气,预混气体中可燃气体的体积分数为10%。通过高速摄像系统采集火焰图像。实验过程如下:首先,清洁管道,连接管道和仪器;然后,检查管道气密性,采用抽真空法配气,向管道中通入预混气体;接着点火触发,数据采集系统对实验数据进行采集;最后,排出废气,重复实验,每组实验重复多次。

    本实验的主要材料为球形多孔非金属材料和聚氨酯(polyurethane,PU)材料,如图2所示。球形多孔非金属材料的主要材质为聚偏氟乙烯,化学性质稳定,具有很强的抗冲击性和耐高温性,并且球形材料为多孔结构,孔隙率大,比表面积大。组合球形内部填充材料主要为聚氨酯,聚氨酯材料的可塑性较强,对爆燃火焰的传播具有一定的抑制作用[14]

    图  2  多孔非金属抑爆材料
    Figure  2.  Porous non-metallic explosion suppression material

    最大爆炸压力(pmax)是评价爆炸强度以及材料抑爆性能的重要指标。为了判断掺氢比对可燃气体爆炸的影响以及单一球形多孔非金属材料和组合球形非金属材料的抑爆性能,设计如下实验方案:首先研究不同掺氢比xpmax的影响,分别对掺氢比为0%、5%、10%的甲烷掺氢预混气体进行测试,记录管道内爆炸压力并进行对比分析;在此基础上,研究多孔抑爆材料对pmax的影响,选取pmax最大的预混气体,分别对单一球形材料和组合球形材料进行抑爆实验,填充于距管道左端1.5 m的位置,填充长度分别为20、30、40 cm,对不同填充长度的同一材料以及相同填充长度的不同材料进行对比分析,判断其抑爆性能。

    掺氢比x表示为

    x=φ(H2)φ(CH4)+φ(H2)×100%
    (1)

    式中:φ(CH4)φ(H2)分别为甲烷-氢气混合气体中甲烷和氢气的体积分数。

    管道中不同掺氢比对预混气体爆炸压力随时间的变化曲线如图3所示。从图3可以看出,1号压力传感器所采集的压力峰值随着氢气体积分数的增大而增大。不同掺氢比对火焰图像的影响如图4所示。对比发现,随着掺氢比的提高,火焰的亮度明显变强,甲烷-空气预混气体掺氢后火焰传播特性及爆炸特性发生显著改变,氢气对甲烷爆炸具有促进作用[15]

    图  3  掺氢比对预混气体爆炸压力的影响
    Figure  3.  Effect of hydrogen doping ratio on explosion pressure of premixed gas
    图  4  不同掺氢比条件下甲烷掺氢爆炸火焰传播图像
    Figure  4.  Flame propagation images of methane hydrogen-doped syngas under different hydrogen doping ratios

    预混气体爆炸压力是反映爆炸破坏效应的重要指标,爆炸压力越大,爆炸造成的后果往往越严重。掺氢比对预混气体最大爆炸压力pmax和最大爆炸压力上升速率(dp/dt)max的影响如图5所示。当掺氢比为0%、5%、10%时,1号压力传感器处的pmax分别为200、235、245 kPa,(dp/dt)max分别为2 250、2 875、3 250 kPa/s。与未掺氢(掺氢比为0%)相比,掺氢比为10%时,1号压力传感器处的pmax提高22.50%,(dp/dt)max提高44.44%,2号传感器处的(dp/dt)max提高26.67%,说明掺氢比较低时,pmax和(dp/dt)max均随着氢气含量的增加而增大,爆炸强度增大。这是由于氢气的活性高于甲烷,反应比甲烷剧烈,发生爆炸反应的时间更短,随着掺氢比的增加,单位体积内氢气的量增多,氢气比例的升高增强了能量释放的集中程度,气体燃烧速率加快,缩短了热量损失时间,导致爆炸压力上升速率提高[16]

    图  5  不同掺氢比条件下最大爆炸压力和最大爆炸压力上升速率曲线
    Figure  5.  Maximum explosion pressure and maximum pressure rise rate for different hydrogen doping ratios

    为研究球形非金属材料对甲烷掺氢后预混气体爆炸的抑制作用,采用爆炸压力较大的掺氢比为10%的预混气体进行对照。空管道与填充不同长度的抑爆材料对爆炸压力的影响如图6所示,其中0 cm代表无抑爆材料。由图6可知,填充抑爆材料后,各个工况下的pmax均小于无抑爆材料。1号、2号压力传感器测得无抑爆材料的pmax分别为245、200 kPa。当抑爆材料的填充长度分别为20、30、40 cm时,1号压力传感器处的pmax分别为205、185、175 kPa,对比无抑爆材料时的pmax分别降低了16.33%、24.49%、28.57%,爆炸压力得到有效抑制。这是由于球形抑爆材料是多孔中空塑料球形结构,有较大的比表面积,并且球形材料把管道空间切割成多个小空间,增大了反应面与球形材料表面的接触面积,促进了球形材料与火焰的热交换,增加了热损失,并且冲击波通过球形材料时,被切割的反应面与后面的球形材料发生碰撞导致能量损失,进而使爆炸压力降低[12]。2号压力传感器处各个工况下的pmax相较于无抑爆材料分别降低5.00%、5.00%、12.50%,抑制效果较弱。这是由于随着火焰传播,火焰通过球形抑爆材料后再无障碍物,燃烧速率加快,火焰到达2号压力传感器时已充分反应,致使降压效果不明显。

    图  6  不同填充长度条件下单一球形材料对预混气体爆炸压力的影响
    Figure  6.  Effect of single spherical material with different filling lengths on explosion pressure of premixed gas

    不同填充长度对(dp/dt)max的影响如图7所示。在低填充长度下,(dp/dt)max有所上升,当填充长度为20 cm时,1号压力传感器处的(dp/dt)max提高7.69%,2号压力传感器处提高36.84%。重复实验均显示(dp/dt)max提高,这是由于当球形材料的填充长度较低时,火焰在穿越过程中的湍流度增大,燃烧传质传热进程加快,燃烧强度增大,火焰传播速度进一步加快[8]。随着填充长度的增加,火焰被分离成多个离散的湍流火焰,最大爆炸压力上升速率逐渐降低,填充长度为40 cm时,相较于无抑爆材料,1号压力传感器处的(dp/dt)max降低33.85%。

    图  7  不同填充长度条件下最大爆炸压力和最大爆炸压力上升速率曲线
    Figure  7.  Maximum explosion pressure and maximum pressure rise rate for different filling lengths

    不同长度的组合球形材料对爆炸压力的影响曲线如图8所示。从图8可以看出:与无抑爆材料相比,填充组合球形材料后pmax明显降低,说明组合球形材料对管道内可燃气体爆炸有较好的抑制效果。当填充长度分别为20、30、40 cm时,1号压力传感器测得的pmax分别为160、155、120 kPa,与无抑爆材料相比,分别衰减了34.69%、36.73%、51.02%,相较于单一球形非金属材料,抑爆性能分别提升了112.43%、49.98%、78.58%;2号压力传感器测得的pmax为135、135、110 kPa,与无抑爆材料相比,分别衰减了32.50%、32.50%、45.00%。可以看到,当填充长度为20和30 cm时,抑爆性能差别不大,这是多孔材料冷却抑制与障碍物加压共同作用的结果: 当火焰经过球形多孔抑爆材料时,材料会吸收热量从而对爆炸产生抑制效果,材料的填充长度越长,抑制效果越好;球形材料在管道中相当于障碍物,火焰经过障碍物时层流会转变成湍流,增加障碍物数量能够明显增大火焰传播过程中的湍流强度,导致压力上升[3]。当填充长度较低时,材料冷却抑制效果与障碍物加压效果相互抵消,最终导致抑制作用较接近。

    图  8  不同填充长度的组合球形材料对预混气体爆炸压力的影响
    Figure  8.  Effect of combined spherical material with different filling lengths on explosion pressure of premixed gas

    图9中可以看出,(dp/dt)max也得到了较大程度的抑制,与无抑爆材料相比,填充长度为20、30、40 cm时,1号压力传感器处的(dp/dt)max分别衰减34.62%、46.15%、53.85%,2号压力传感器处的(dp/dt)max最高衰减了47.37%,抑制效果明显。与单一球形材料相比,组合球形材料的抑爆性能得到大幅增强。这是由于爆炸穿越填充材料时,被球形材料多孔结构离散后的湍流火焰和冲击波能量与球形多孔材料中的聚氨酯材料接触,聚氨酯材料粗糙的壁面消耗了链式反应中的自由基,阻碍燃烧链式反应的进行,并使部分火焰热量转移到聚氨酯材料中,在火焰与聚氨酯材料发生碰撞的过程中,聚氨酯材料与冲击波产生摩擦,导致部分冲击波能量转化为热量并消散掉[14]

    图  9  不同填充长度条件下最大爆炸压力和最大爆炸压力上升速率曲线
    Figure  9.  Maximum explosion pressure and maximum pressure rise rate for different filling lengths

    基于自主搭建的气体爆炸管道平台,研究了不同掺氢比条件下甲烷掺氢的爆炸特性以及不同类型球形抑爆材料的抑爆性能,得到以下结论。

    (1) 氢气对甲烷-空气爆炸具有一定的促进作用。在低掺氢比下,随着可燃气体中氢气体积分数的增加,火焰传播速率加快,最大爆炸压力和最大爆炸压力上升速率升高。

    (2) 单一球形多孔非金属材料对可燃气体爆炸具有两方面作用:一方面为抑制作用,体现在填充材料后,球形材料会吸收能量,导致最大爆炸压力降低,并且填充长度越长,抑制效果越好;另一方面为促进作用,体现在填充长度较低时,会导致燃烧传质传热进程加快,火焰传播加速,导致最大爆炸压力上升速率上升。

    (3) 组合球形多孔非金属材料与单一球形多孔非金属材料相比,抑爆性能大幅提升,降压效果明显,当填充长度为40 cm时,最大爆炸压力衰减51.02%。球形材料抑爆性能受填充长度影响,在冷却抑制与障碍物增压的共同作用下,填充20和30 cm的组合球形多孔非金属材料的抑爆效果相差不大。

  • 图  水锤效应过程[1]

    Figure  1.  Process of hydrodynamic ram[1]

    图  试验装置[1]、几何模型、1/2物质点离散模型和压力测点位置[1]

    Figure  2.  Test device[1], geometric profile of cross section, 1/2 discrete model of MPM and location of pressure measuring points[1]

    图  水锤效应过程对比

    Figure  3.  Comparison of hydrodynamic ram process

    图  100%充液率下测点压力时程曲线

    Figure  4.  Time-history curves of pressure at measuring points with 100% filling rate

    图  双层隔板充液结构模型和压力测点位置

    Figure  5.  Model of double-layer plate liquid-filled structure and the location of the pressure measuring points

    图  双层隔板间距分布

    Figure  6.  Setting of double-layer plate spacing

    图  不同双层隔板间距情况下空穴的形成过程

    Figure  7.  Process of cavitation with different settings of double-layer plate spacing

    图  双层隔板间距不同时射弹的剩余速度

    Figure  8.  Residual velocity of projectile under different settings of double-layer plate spacing

    图  双层隔板间距不同时不同测点的压力时程曲线

    Figure  9.  Time-history curves of pressure at different measuring points under different settings of double-layer plate spacing

    图  10  出入射壁及双隔板的变形情况

    Figure  10.  Deformation of the entry wall, exit wall and the double-layer plates

    图  11  双层隔板位置

    Figure  11.  Position setting of the double-layer plate

    图  12  不同双层隔板位置情况下空穴的形成过程

    Figure  12.  Process of cavitation with different position settings of the double-layer plate

    图  13  双层隔板位置不同时射弹的剩余速度

    Figure  13.  Residual velocity of projectile with different position settings of the double-layer plates

    图  14  双层隔板位置不同时不同测点的压力时程曲线

    Figure  14.  Time-history curve of pressure at different measuring points with different position settings of the double-layer plates

    图  15  出、入射壁及双隔板变形情况

    Figure  15.  Deformation of the entry wall, exit wall and the double-layer plates

    表  1  油箱和射弹的材料参数[7]

    Table  1.   Material parameters of fuel tank and projectile[7]

    Materialρ/(kg∙m−3)E/GPaνA/MPaB/MPanCmD1
    Al6063-T52700710.332001440.62010.2
    4340 steel7830207 0.28
    下载: 导出CSV

    表  2  水的材料参数[7]

    Table  2.   Material parameters of water[7]

    Materialρ0/(kg∙m−3)vd/(mPa∙s)c0/(m∙s−1)S1S2S3γ0a
    Water10000.8914481.979000.113.0
    下载: 导出CSV
  • [1] VARAS D, LÓPEZ-PUENTE J, ZAERA R. Experimental analysis of fluid-filled aluminium tubes subjected to high-velocity impact [J]. International Journal of Impact Engineering, 2009, 36(1): 81–91. doi: 10.1016/j.ijimpeng.2008.04.006
    [2] BALL R E . The fundamentals of aircraft combat survivability: analysis and design [M]. 2nd ed. Reston, VA: American Institute of Aeronautics and Astronautics, 2003: 799−864.
    [3] BALL R E. Aircraft fuel tank vulnerability to hydraulic ram: modification of the northrop finite element computer code BR-1 to include fluid-structure interaction; theory and user’s manual for BR-1HR: NPS-57BP74071 [R]. California: Naval Postgraduate School, 1974.
    [4] FUHS A E, BALL R E , POWER H L. FY 73 hydraulic ram studies: AD776536 [R]. California: Naval Postgraduate School, 1974.
    [5] LUNDSTROM E A. Hydraulic ram model for high explosive ammunition [C]//ASME Pressure Vessels and Piping Conference. Honolulu, HI, 1989.
    [6] 李亚智, 陈钢. 充液箱体受弹丸撞击下动态响应的数值模拟[J]. 机械强度, 2007, 29(1): 5.

    LI Y Z, CHEN G. Numerical simulation of liquid-filled tank response to projectile impact [J]. Journal of Mechanical Strength, 2007, 29(1): 5.
    [7] VARAS D. Numerical modelling of the hydrodynamic ram phenomenon [J]. International Journal of Impact Engineering, 2009, 36(3): 363–374. doi: 10.1016/j.ijimpeng.2008.07.020
    [8] VARAS D, ZAERA R, LÓPEZ-PUENTE J. Numerical modelling of partially filled aircraft fuel tanks submitted to hydrodynamic ram [J]. Aerospace Science and Technology, 2012, 16(1): 19–28. doi: 10.1016/j.ast.2011.02.003
    [9] ARTERO-GUERRERO J A, VARAS D, PERNAS-SANCHEZ J, et al. Experimental analysis of an attenuation method for hydrodynamic ram effects [J]. Materials & Design, 2018, 155: 451–462.
    [10] 张宇, 王彬文, 刘小川, 等. 充液格栅结构抗射弹冲击特性研究 [J]. 科学技术与工程, 2020, 20(18): 9.

    ZHANG Y, WANG B W, LIU X C, et al. Research on anti-projectile impact characteristics of liquid-filled grid structure [J]. Science Technology and Engineering, 2020, 20(18): 9.
    [11] 李营, 赵鹏铎, 张春辉, 等. 空气夹层对含液结构在球形弹体侵彻作用下动态响应的影响 [J]. 振动与冲击, 2018, 37(3): 186–195. doi: 10.13465/j.cnki.jvs.2018.03.030

    LI Y, ZHAO P D, ZHANG C H, et al. Influences of air-contain structure on dynamic responses of liquid-filled structures under spherical projectile penetration [J]. Journal of Vibration and Shock, 2018, 37(3): 186–195. doi: 10.13465/j.cnki.jvs.2018.03.030
    [12] 张雄. 物质点法 [M]. 北京: 清华大学出版社, 2013: 38−42.

    ZHANG X. Material point method [M]. Beijing: Tsinghua University Press, 2013: 38−42.
  • 加载中
图(15) / 表(2)
计量
  • 文章访问数:  209
  • HTML全文浏览量:  129
  • PDF下载量:  43
出版历程
  • 收稿日期:  2022-06-06
  • 修回日期:  2022-07-05
  • 网络出版日期:  2023-02-07
  • 刊出日期:  2023-02-05

目录

/

返回文章
返回