• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

楔形罩线型聚能装药侵彻钢锭

武双章 顾文彬 李旭锋

张文超, 王舒, 梁增友, 覃彬, 卢海涛, 陈新元, 卢文杰. 基于空气流场压力分析的头盔冲击波防护效能研究[J]. 爆炸与冲击, 2022, 42(11): 113201. doi: 10.11883/bzycj-2021-0411
引用本文: 武双章, 顾文彬, 李旭锋. 楔形罩线型聚能装药侵彻钢锭[J]. 爆炸与冲击, 2016, 36(3): 353-358. doi: 10.11883/1001-1455(2016)03-0353-06
ZHANG Wenchao, WANG Shu, LIANG Zengyou, QIN Bin, LU Haitao, CHEN Xinyuan, LU Wenjie. A study of blast wave protection efficiency of helmet based on air flow field pressure analysis[J]. Explosion And Shock Waves, 2022, 42(11): 113201. doi: 10.11883/bzycj-2021-0411
Citation: Wu Shuangzhang, Gu Wenbin, Li Xufeng. Penetration of steel ingot by linear shaped charge with cuniform cove[J]. Explosion And Shock Waves, 2016, 36(3): 353-358. doi: 10.11883/1001-1455(2016)03-0353-06

楔形罩线型聚能装药侵彻钢锭

doi: 10.11883/1001-1455(2016)03-0353-06
详细信息
    作者简介:

    武双章(1977—),男,硕士,讲师

    通讯作者:

    顾文彬, guwenbin1@aliyun.com

  • 中图分类号: O385;TD235.3

Penetration of steel ingot by linear shaped charge with cuniform cove

  • 摘要: 为了获得楔形罩线型聚能装药射流侵彻钢锭的特点和规律,利用ANSYS/LS-DYNA软件建立了三维数值计算模型,并进行了数值模拟。在此基础上,进行了实际的切割实验:对于某一炸高条件,切割深度随时间的增加先快速增加,后缓慢增加;对于不同的炸高,切割深度随炸高的增加而缓慢增加,当炸高为40~60 mm和70~120 mm时,侵彻深度表现出对炸高的不敏感性;同时,获得了不同炸高条件下不同时刻的钢锭切口断面形状。结果表明,数值模拟和切割实验结果有很好的一致性,可以用该三维数值计算模型模拟实际切割器侵彻钢靶的过程,并获得其特点和规律。
  • 近年来国外战争及恐怖袭击中,各种爆炸性武器尤其是简易爆炸装置(improvised explosive device, IED)的大量使用使得爆炸冲击伤成为战场环境中人员遭受的主要创伤,由于个人防护装备(主要保护躯干)的更新和发展,冲击伤最敏感的靶器官如肺、肠等的损伤及其所致死亡大大减少,但颅脑冲击伤却凸显出来,成为现代战争特征性创伤[1]。爆炸冲击波作用人体诱发的颅脑冲击伤(blast-induced traumatic brain injury, BTBI)可导致人体记忆缺失、注意力不集中、情绪紧张等问题,还会导致睡眠障碍、疲劳、头晕、易怒、头痛和癫痫发作,严重影响士兵作战能力,是伊拉克和阿富汗最常见的军事伤害[2-3]

    目前不同种类战斗头盔如步兵头盔、飞行员头盔、坦克乘员头盔、摩托兵头盔等被用来保护颅脑免受伤害,但是现代战斗头盔的设计主要是用来避免钝物或枪弹对头部的直接打击,防止可能对头部造成的非贯穿性伤害等,它们对爆炸冲击波的防护效果很大程度上是未知的[4]。学者们对头盔的爆炸冲击波防护能力已进行了一些研究,Zhang等[5]利用验证的头部、头盔有限元模型,对戴头盔和不戴头盔的情况下,头部对爆炸载荷的响应进行了研究,模拟结果表明,先进战斗头盔可使头部颅内压减轻10%~35%,同时导致平均峰值脑应变和应变与应变率乘积降低30%。Nyein等[6]对超压为1 MPa的冲击波冲击未戴头盔、戴头盔、戴头盔和面罩的头部这三种工况进行数值模拟研究,结果表明现有的先进作战头盔对爆炸冲击波没有显著防护效果,但是也不会产生负面效果;头盔和面罩结构组合能够阻碍冲击波直接传递到面部,使应力传递到颅内腔时出现衰减和延迟,降低了颅内应力,防护效果较好。Sarvghad-Moghaddam等[7]通过计算流体动力学模拟,研究了可压缩空气从不同方向接近头盔结构时的湍流、超声速流动,评估了不同入射超压和头盔间隙大小对脑后冲击波压力汇聚的影响。Tan等[8]利用颅脑-头盔仿真模型进行爆炸冲击波作用头部研究,得出爆炸冲击波能够通过颅骨传播,造成较高的颅内压,来自正面爆炸的压力波可能会通过头盔的侧面进入头盔间隙并对颅脑后方进行第二次冲击。Li等[9]通过试验及数值模拟对头盔的冲击波防护机理进行了研究,研究结果表明冲击波正面作用头盔防护下头部时,头部颅顶区域超压峰值减小,但颅后区域超压峰值增大。康越等[10]对头盔-头模系统前部、前额部、顶部、后部、耳部以及眼部等重点区域的冲击波超压峰值和持续作用时间进行了对比分析,得出提高单兵头面部防护装备结构密闭性(如佩戴眼镜、耳罩或者防护面罩),可有效阻止冲击波进入头盔-头模系统内部,减弱叠加汇聚效应,提高单兵头面部装备防护性能。

    基于以上可知,头盔的冲击波防护能力有一些结论得到了学者的统一认可,但是部分结论存在矛盾与冲突。国内相关研究起步较晚,冲击波作用于佩戴头盔的头部时,空气流场压力分布及冲击波传播规律尚不明确,为降低士兵在战场环境中因爆炸冲击波致颅脑创伤概率,提高战斗头盔的冲击波防护能力,亟需对战斗头盔的冲击波防护效能进行研究。本文通过爆炸试验与数值模拟方法相结合,对典型战斗头盔的冲击波防护能力进行研究,结合试验结果对数值模拟结果有效性进行验证,同时采用数值模拟方法对冲击波作用有无防护头部后流场压力演化规律进行描述,进一步通过数值仿真方法对衬垫结构对头盔的冲击波防护能力影响进行研究,以期对冲击波防护装备优化设计提供技术指导及数据支撑。

    爆炸环境下头部模型周围空气压力场测试试验在室外环境下进行,爆炸源为50 g TNT药柱。进行了两种工况的爆炸试验,每种工况下,药柱与头部之间的距离是1 m。工况1为无防护头部模型(裸头模),工况2为配戴头盔防护头部模型(盔壳头模),为防止传感器引起的测量误差,每组均设置对称裸头模型作为对照组。在头模上安装三个压力传感器(PCB M102/106B系列),测量前额、颅顶和颅后的压力。这些传感器的前端与头模表面重合。战斗头盔外壳厚度为6.5 mm,头模具有典型的人体头部外观结构,由塑胶材料制成。试验场地布置及头部模型分别如图1图2所示。

    图  1  试验场地布置
    Figure  1.  Layout of test device
    图  2  头部替代模型
    Figure  2.  The dummy head models

    冲击波作用于无防护头部时,对前额部位直接进行冲击,由于波阻抗不匹配,冲击波在前额部位反射叠加产生较高的超压峰值(277 kPa),随后冲击波继续向前传播并逐渐衰减,先后作用于颅顶(95 kPa)、颅后(52 kPa),颅后超压峰值相比前额衰减81.3%(如图3(a)所示)。冲击波作用于头盔防护下头部时,由于头盔的阻挡作用,冲击波不能直接作用前额部位,前额部位超压峰值(151 kPa)比无防护时低,随后由于绕流作用冲击波传入头-头盔间隙先后作用于颅顶(64 kPa)、颅后(114 kPa)部位,颅后超压相较前额衰减24.5%(如图3(b)所示),颅后超压大于颅顶超压(增加50 kPa)。由于头盔的存在,颅后超压出现压力增强现象。

    图  3  实验获得的压力-时间曲线
    Figure  3.  Experimental time histories of pressure

    分别取不同工况下相同位置处的超压峰值进行对比,由图4可知,头盔防护下头部前额、颅顶处超压峰值较无防护时都降低,分别为无防护时的54.5%、67.3%。颅后空气超压出现增高,增加到无防护时的2.19倍。测试结果说明头盔防护可以有效减小冲击波对头部的直接冲击,但是会对颅后冲击波超压起到增强作用,对颅后的冲击波防护起到负面效果。下文将对头盔防护下颅后冲击波超压增强现象结合数值模拟进行分析说明。

    图  4  峰值超压对比
    Figure  4.  Peak pressure contrast

    所建立的头部有限元模型中,皮肤、颅骨几何数据来源于陆军军医大学可视人体切片数据集(CVH)[11]。基于皮肤、颅骨几何模型,利用Geomagic对几何模型进行优化,利用Hypermesh对几何进行头部有限元网格划分。建立的头部有限元模型包括皮肤、颅骨、面颅、脑脊液、脑组织等组织器官及硬脑膜、软脑膜等膜结构,膜结构采用四边形壳单元划分,其余部分均采用6面体实体单元,模型单元尺寸约3 mm,有限元模型如图5所示。

    图  5  头部有限元模型
    Figure  5.  Finite element model of head

    爆炸空气环境用边长60 cm的立方体区域表征,头部放置空气域中心,冲击波与头部模型发生作用部分进行网格加密,以保证耦合界面不出现压力泄露,其他部分采用渐变网格,网格逐渐变大从而减少计算成本,空气域中心区域单元尺寸4 mm,此时相应的计算结果收敛。

    冲击波作用头部过程采用三维多物质ALE流固耦合方法[12]来实现,流固耦合算法用于定义空气和头部之间的相互作用,此算法使空气给头部施加压力载荷,而头部则相当于空气的边界条件,用于约束空气的运动。有限元模型中空气单元采用Euler网格描述,头部单元采用Lagrangian网格描述,利用*CONSTRAINED_LAGRANGE_IN_SOLID关键字进行罚耦合算法的定义,定义空气网格边界参考压力为标准大气压101.325 kPa,这样冲击波负压段压力才会回升,与实际情况相符。

    试验使用的头盔具有地面部队人员佩戴军用头盔典型结构,其内部吊带悬挂系统能够避免头部与盔壳直接接触,对钝器冲击起到缓冲效果,头盔实物图如图6所示。

    图  6  典型战斗头盔实物图
    Figure  6.  Actual picture of typical combat helmet

    冲击波作用头部属于瞬态冲击,可以忽略头部悬挂在冲击波作用时对头部力的传导,在建模过程中将内部吊带去除,对头盔几何模型进行简化。战斗头盔凯夫拉盔壳厚6.9 mm,泡沫衬垫厚1~3 cm,头盔及泡沫衬垫数值模型如图7

    图  7  头盔、泡沫衬垫数值模型
    Figure  7.  Simulation model of helmet and foam pad

    脑组织材料接近不可压缩,是一种质地很软的生物软组织,其表现为不可压缩性、非线性、各向异性及粘弹性,大量脑组织实验证明脑组织的变形只取决于其剪切模量,本文脑组织采用粘弹性本构模型[13,14](*MAT_VISCOELASTIC)对材料性能进行表征,其剪切弹性模量的计算公式为

    G(t)=G0+(G0G)eβt
    (1)

    式中:G0为短效剪切模量,G为长效剪切模量, β为衰减系数。

    脑脊液采用弹性流体本构模型[15](*MAT_ELASTIC_FLUID)表征。

    头部模型皮肤、颅骨、面颅、硬脑膜、软脑膜均采用线弹性模型进行表征,参考文献[16-20],各组织结构材料参数如表13所示。

    表  1  弹性模型材料参数
    Table  1.  Material parameters of elastic model
    结构ρ/(kg·m−3)E/MPaμ
    皮肤120016.70.42
    面颅171053700.19
    密质骨2000150000.22
    松质骨130010000.24
    硬脑膜113031.50.45
    软脑膜113011.50.45
    注:ρ为密度,E为弹性模量,μ为泊松比.
    下载: 导出CSV 
    | 显示表格
    表  2  黏弹性模型材料参数[13-14]
    Table  2.  Material parameters of viscoelastic model
    结构ρ/(kg·m−3)G0/kPaG/kPaβ/(s−1)K/MPa
    脑组织10401.660.92816.95557
    注:K为体积模量.
    下载: 导出CSV 
    | 显示表格
    表  3  弹性流体模型材料参数[15]
    Table  3.  Material parameters of elastofluid model
    结构ρ/(kg·m−3)K/MPa
    脑脊液999.81960
    下载: 导出CSV 
    | 显示表格

    头盔盔壳材料使用复合材料模型*MAT_COMPOSITE_DAMAGE进行表征。头盔泡沫衬垫是一种聚氨酯硬泡沫,具有加载率敏感性的应力-应变行为,使用MAT_LOW_DENSITY_FOAM材料模型进行表征,其材料密度为6.1×10−2 g/cm3,杨氏模量为8.4 MPa。参考文献[5,13],确定头盔及泡沫衬垫具体材料参数。头盔盔壳结构材料参数见表4,其中:ρ为材料密度,E11E22E33为各向弹性模量,ν12ν13ν32为各向泊松比,G12G23G13为各向剪切模量,S11S22 为压缩模量,Sc为剪切强度,Sn为压缩强度。泡沫衬垫应力-应变曲线如图8所示。

    表  4  头盔盔壳模型材料参数[5,13]
    Table  4.  Material parameters of helmet model
    ρ/(kg·m−3)E11/GPaE22/GPaE33/GPaν12ν13ν32
    123018.518.560.250.330.33
    G12/GPaG23/GPaG13/GPaS11/MPaS22/GPaSc/GPaSn/GPa
    0.772.502.500.5580.5550.5551.086
    下载: 导出CSV 
    | 显示表格
    图  8  硬泡沫应力-应变曲线
    Figure  8.  Hard foam stress-strain curve

    利用LS-DYNA显式动力分析软件模拟了50 g TNT距离头部1 m处爆炸时冲击波与头部相互作用过程。空气参考压力为标准大气压101.325 kPa,在欧拉域迎爆面设置压力加载边界,冲击波从此边界输入。为了避免压力波在欧拉域边界反射对头部产生不真实的二次冲击,在欧拉域周围设置非反射边界[21]

    空气冲击波超压在正压段大致按指数规律衰减,一些经验公式可以描述此衰减过程,其中Friedlander方程[8]较接近实际且又简单易于计算:

    p=pSet/t(1tt)
    (2)

    式中:pS为入射超压;t为正压作用时间。

    冲击波即将与头部接触时冲击波超压为106 kPa,正压持续时间为690 μs。冲击波输入面及作用于头部的载荷曲线如图9

    图  9  头部冲击波载荷
    Figure  9.  Acting blast wave loads on the head

    进行典型战斗头盔对爆炸冲击波的防护效果模拟研究,需对冲击波作用有无防护头部后流场压力分布进行验证,首先利用LS-DYNA显式动力分析软件对上文试验工况进行数值模拟,试验结果及仿真结果对比见图10表5

    图  10  仿真与试验超压曲线对比
    Figure  10.  Comparison of simulation and experimental results
    表  5  观测点超压峰值及误差
    Table  5.  Overpressure peak and error at observation points
    模型方法pfront/kPaηfront/%ptop/kPaηtop/%pback/kPaηback/%
    裸头试验27719.4954.2525.7
    数值模拟2239155
    盔壳试验1511.3641.51141.7
    数值模拟14965116
    注:pfrontptoppback分别为为前额、颅顶和颅后的超压峰值,ηfrontηtopηback分别为对应的误差。
    下载: 导出CSV 
    | 显示表格

    图10表5可得,有无头盔防护下头部前额、颅顶、颅后观测点仿真结果和试验结果峰值压力误差较小,且压力变化趋势与试验结果基本一致,说明数值仿真能较好反映爆炸冲击波作用头部后空气压力场分布情况。

    基于上述有限元模型及爆炸冲击波加载,利用LS-DYNA显式动力分析软件对爆炸冲击波正面作用无防护、无衬垫头盔防护、带衬垫头盔防护下头部3种工况进行数值模拟,在头部前额、颅后、耳侧、顶侧选取6个压力采集点提取压力时程曲线,结合冲击波流场压力云图分析爆炸冲击波与头部相互作用过程,得出不同工况下冲击波流场压力变化规律,对泡沫衬垫对头盔的冲击波防护效果影响进行研究。

    2.5.1   无头盔防护空气流场压力分布

    无头盔防护时头部周围特征部位压力采集点位置如图11所示,采集点超压-时程曲线如图12所示,结合空气流场压力云图对冲击波作用头部过程进行分析。

    图  11  头部周围采样点位置示意图
    Figure  11.  Schematic diagram of sampling points around the head
    图  12  采样点超压曲线
    Figure  12.  Pressure curves of sampling points

    爆炸冲击波正面作用无防护头部时,空气流场压力演化过程如图13所示。t=400 μs时冲击波即将与头部接触,此时冲击波峰值超压为106 kPa;随后冲击波与头部相互作用,由于波阻抗不匹配,冲击波在空气-皮肤介质界面反射叠加,超压出现增强;t=430 μs时,冲击波作用于A点,超压峰值达到164 kPa,为自由场相同位置(98 kPa)的1.67倍;t=600 μs时,由于头部阻挡,冲击波在头部周围除发生反射外还发生绕流作用向颅后传播;由图13可以看出,在稀疏波作用下面部空气受影响而改变了运动方向,形成围绕头部向颅后传播的的高压气流,同时又和相邻的入射波一起作用,变成绕流向颅后传播。t=1000 μs时,绕流波会在颅后部位重新汇聚,汇聚作用的合成波区会出现峰值超压增强,对头部造成二次冲击,但是相比于面部首次冲击,二次冲击超压峰值(56 kPa)显著降低。耳部凹面结构使得冲击波反射叠加增强,出现更高的超压,峰值超压增大为150 kPa,为自由场相同位置(86 kPa)的1.7倍。冲击波对头部作用顺序为A(164 kPa)→B(116 kPa)→(E(150 kPa)与F(91 kPa))→C(39 kPa)→D(56 kPa)。

    图  13  爆炸冲击波压力云图
    Figure  13.  Pressure nephogram of blast wave
    2.5.2   头盔(无衬垫)防护下空气流场压力分布

    头盔(无衬垫)防护下头部周围特征部位压力采集点位置如图14所示,采集点超压-时程曲线如图15所示,结合空气流场压力云图对冲击波作用头盔防护下头部过程进行分析。

    图  14  头部周围采样点位置示意图
    Figure  14.  Schematic diagram of sampling points around the head
    图  15  采样点超压曲线
    Figure  15.  Pressure curves of sampling points

    爆炸冲击波正面作用头盔(无衬垫)防护下头部时,空气流场压力演化过程如图16所示。t=400 μs时冲击波即将与头部接触,此时冲击波峰值超压为106 kPa与上述工况一致。t=480 μs时冲击波作用至A点,由于头盔的存在避免了冲击波对前额的直接冲击,前额(A点)部位冲击波超压峰值为92 kPa,较无防护时(164 kPa)降低43.9%,同时作用时刻较无防护出现延迟。t=600 μs时由于头盔-头部阻挡,冲击波除发生反射透射外还围绕头盔-头部发生绕流作用,冲击波锋面分为2个部分,一部分在头盔外围绕头部传播,另一部分绕流进入头-头盔间隙进行传播,同时头盔的阻挡使得冲击波在头盔与颅顶间隙传播滞后于头盔外部,冲击波绕流到颅后会从头部后方向颅顶方向传播,由前后不同方向传入头-头盔间隙的高压气流会在颅后一处发生碰撞,碰撞区的压力骤然升高,C点压力(130 kPa)上升为无防护时(39 kPa)的3.3倍。头盔外部绕流波到达头部后方会形成一个压力较高的合成波区,合成的冲击波继续向后传播。

    图  16  爆炸冲击波压力云图
    Figure  16.  Pressure nephogram of blast wave

    由于头盔的阻挡高压气流不能很快扩散,C点冲击波正压作用时间显著增长(由560 μs增加到740 μs)。此结构头盔对耳部无正面遮挡,与无防护相比,耳部E点超压变化(由150 kPa变为151 kPa)较小,对侧顶F位置超压(由91 kPa变为49 kPa)有衰减效果,但是会使正压作用时间增长(由440 μs变为970 μs)。头盔防护下冲击波对头部的作用顺序为A(92 kPa)→B(70 kPa)→E(151 kPa)→F(49 kPa)→D(85 kPa)→C(130 kPa)。由此解释了试验中头盔对颅后冲击波防护的负面效果形成机制。

    2.5.3   头盔(带衬垫)防护下空气流场压力分布

    头盔(带衬垫)防护下头部周围特征部位压力采集点位置如图17所示,采集点超压-时程曲线如图18所示,结合空气流场压力云图对冲击波作用带衬垫头盔防护下头部过程进行分析。

    图  17  头部周围采样点位置示意图
    Figure  17.  Schematic diagram of sampling points around the head
    图  18  特征部位采样点超压曲线
    Figure  18.  Pressure curves of sampling points

    爆炸冲击波正面作用头盔(带衬垫)防护下头部时,空气流场压力演化过程如图19所示。由于冲击波绕流作用即使有泡沫衬垫阻挡高压气流仍会进入头-头盔间隙,在间隙内形成冲击波,t=480 μs时冲击波作用至A点,头盔防护避免了冲击波对前额的直接冲击,相较无防护前额部位空气超压降低23 kPa (由164 kPa变为141 kPa),但由于前额处衬垫阻挡冲击波在衬垫-前额间隙反射叠加使压力增强,故前额处空气超压较无衬垫时升高53.2%(由92 kPa变为141 kPa)。与无衬垫头盔防护时同理,冲击波通过头-头盔间隙传至颅后,在颅后(C点)位置仍会出现高压气流碰撞作用使得空气超压相对裸头模升高1.4倍(由56 kPa变为81 kPa),但是衬垫对高压气流的传播有阻挡作用,故相对于无衬垫头盔防护时C点压力减小37.6%(由130 kPa变为81 kPa),同时侧顶F位置空气超压也出现减小(由49 kPa变为35 kPa)。由以上可知加入衬垫结构后颅后的压力汇聚效应仍然存在,但是相同位置空气超压都出现一定程度减小。由于此结构头盔对耳部无正面遮挡,E点超压在不同工况下基本一致。加衬垫头盔防护下冲击波对头部的作用顺序为A(141 kPa)→E(154 kPa)→B(51 kPa)→F(35 kPa)→D(86 kPa)→C(81 kPa)。

    图  19  爆炸冲击波压力云图
    Figure  19.  Pressure nephogram of blast wave
    2.5.4   脑组织压力分析

    为了对头盔防护效果有更清晰的认识,在脑组织额叶、顶叶、枕叶三个特征部位各取一个压力观测点,观察脑组织压力变化情况,观测点位置如图20。所使用头部模型脑组织压力响应在先前研究中已通过有效性验证[22]

    图  20  压力观测点
    Figure  20.  Pressure observation points

    提取每组工况下3个观测点压力数据并进行分析,组织压力变化曲线图21~图23所示。

    图  21  额叶观测点压力变化
    Figure  21.  Pressure curves of frontal lobe observation points
    图  22  顶叶观测点压力变化
    Figure  22.  Pressure curves of parietal lobe observation points
    图  23  枕叶观测点压力变化
    Figure  23.  Pressure curves of occipital lobe observation points

    脑组织所受压力过高时会出现脑出血造成严重脑损伤,过高的负瞬态压力也可引起脑空化性损伤[23],根据颅内压耐受标准,当颅内压峰值超过235 kPa时,会发生严重脑损伤,低于173 kPa时,会发生轻微损伤或不发生损伤[24-25]。由以上不同工况下脑组织压力曲线可知,冲击波作用头部后,头盔防护下(无衬垫)额叶、顶叶组织压力峰值分别较无防护时降低59%、63.3%,枕叶组织压力峰值较无防护升高10.8%。头盔结构可以有效减小额叶及顶叶脑组织压力使其降于损伤阈值以下,但在头盔防护下,枕叶组织压力峰值会升高,头盔防护结构在颅后起到负面效果。有无衬垫头盔防护下组织压力峰值相差较小,但通过空气流场压力分布可知衬垫结构可使颅后空气超压降低,因此可以通过设计泡沫衬垫结构来减弱头盔颅后冲击波防护的负面效果,提高头盔对爆炸冲击波的防护能力。

    通过试验及数值模拟方法进行了头盔防护颅脑冲击伤效能研究,研究结果如下:

    (1) 头盔防护能避免前额部位直接受到冲击波冲击,使前额、颅顶部位空气超压峰值分别衰减为无防护时的54.5%、67.3%,但是会使颅后空气超压增强为无防护的2.19倍;

    (2) 头盔的阻挡作用导致头盔与头顶间隙冲击波传播滞后于头盔外部,冲击波传至头部后方时会进行压力汇聚,高压气流会在颅后头盔间隙内发生碰撞,致使颅后冲击波超压峰值显著升高,头盔的阻挡使高压流不能很快扩散,冲击波正压作用时间显著增长;

    (3) 头盔防护下泡沫衬垫可使颅后冲击波超压减弱,消弱头盔对颅后冲击波防护的负面效果,使头盔的冲击波防护效能得到提升;冲击波正面冲击下耳廓结构会使冲击波超压放大,为自由场相同位置1.7倍,耳是冲击波作用重要靶器官。

  • 图  1  LSC和钢锭有限元模型

    Figure  1.  Finite model of LSC and steel ingot

    图  2  实验设置和切割效果

    Figure  2.  Experimental setup and cutting effect

    图  3  不同炸高时的侵彻深度

    Figure  3.  Penetration depths at different standoffs

    图  4  侵彻深度随炸高的变化

    Figure  4.  Penetration depths vs standoffs

    图  5  70 mm炸高时LSC切割钢锭过程

    Figure  5.  Processes of LSC cutting steel ingot at 70 mm standoff

    图  6  110 mm炸高时LSC切割钢锭过程

    Figure  6.  Processes of LSC cutting steel ingot at 110 mm standoff

    图  7  炸高不同时的切口断面形状

    Figure  7.  Shapes of cutting sections at different standoffs

    表  1  材料模型和状态方程

    Table  1.   Material models and state equations

    部位 材料 模型 状态方程 单元尺寸/(mm×mm×mm) 单元数
    装药 B炸药 HIGH_EXPLOSIVE_BURN JWL 2.00×2.00×2.00 38 750
    药形罩 纯铁 JOHNSON_COOK Grüneisen 1.50×0.75×2.00 6 000
    空隙 空气 NULL Grüneisen 2.00×2.00×2.00 399 400~520 900
    45钢 JOHNSON_COOK Grüneisen 2.00×2.00×2.00 120 000
    下载: 导出CSV

    表  2  实验结果

    Table  2.   Experimental results

    h/mm Pnum/mm Pexp/mm δΡ/%
    50 68.32 62.60 9.14
    80 70.20 64.80 7.69
    110 72.14 66.40 7.96
    下载: 导出CSV
  • [1] Dave N, Meryl M. Historical development of linear shaped charge[C]//43rd AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit. Cincinnati, OH: AIAA, 2007. https://www.researchgate.net/publication/268482273_Historical_Development_of_Linear_Shaped_Charge
    [2] 武双章, 顾文彬, 李旭锋, 等.切割防护门的线型聚能装药参数正交优化设计研究[J].爆破器材, 2012, 41(5):12-16. http://d.old.wanfangdata.com.cn/Periodical/bpqc201205004

    Wu Shuangzhang, Gu Wenbin, Li Xufeng, et al. Study on the orthogonal optimization design of the parameters of linear shaped charge cutting the protective door[J]. Explosive Material, 2012, 41(5):12-16. http://d.old.wanfangdata.com.cn/Periodical/bpqc201205004
    [3] 关焕文, 张絮, 封文春, 等.聚能爆破切割清理弹射通道技术[J].航空学报, 2013, 34(1):52-57. http://d.old.wanfangdata.com.cn/Periodical/hkxb201301007

    Guan Huanwen, Zhang Xu, Feng Wenchun, et al. Technology of explosive cutting to clear ejection path[J]. Acta Aeronautica et Astronautica Sinica, 2013, 34(1):52-57. http://d.old.wanfangdata.com.cn/Periodical/hkxb201301007
    [4] 夏红兵, 李磊, 马宏昊, 等.多向线性聚能切割弹的研究与应用[J].含能材料, 2013, 21(1):85-91. doi: 10.3969/j.issn.1006-9941.2013.01.019

    Xia Hongbing, Li Lei, Ma Honghao, et al. Application of the multi-directional linear cumulative cutter[J]. Chinese Journal of Energetic Materials, 2013, 21(1):85-91. doi: 10.3969/j.issn.1006-9941.2013.01.019
    [5] 李志强, 刘晓明, 赵勇刚, 等.微爆索线性切割航空有机玻璃的实验研究和有限元分析[J].爆炸与冲击, 2007, 27(5):385-389. doi: 10.3321/j.issn:1001-1455.2007.05.001

    Li Zhiqiang, Liu Xiaoming, Zhao Yonggang, et al. Experimental study and finite element analysis of linear cutting aerial PMMA using micro detonation cord[J]. Explosion and Shock Waves, 2007, 27(5):385-389. doi: 10.3321/j.issn:1001-1455.2007.05.001
    [6] Meryl M, Tom K. Performance analysis of linear shaped charge for aerospace applications[C]//41st AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit. Tucson, Arizona: AIAA, 2005. https://www.researchgate.net/publication/268481607_Performance_Analysis_of_Linear_Shaped_Charge_for_Aerospace_Applications
    [7] 王志华, 李志强, 任会兰, 等.有机玻璃的几何形状对爆炸切割影响的研究[J].高压物理学报, 2010, 24(1):43-48. doi: 10.3969/j.issn.1000-5773.2010.01.008

    Wang Zhihua, Li Zhiqiang, Ren Huilan, et al. Finite element analysis on the influence of PMMA geometry on explosion cutting[J]. Chinese Journal of High Pressure Physics, 2010, 24(1):43-48. doi: 10.3969/j.issn.1000-5773.2010.01.008
    [8] 武双章, 顾文彬.端面效应对楔形罩切割器射流成型影响的试验和数值模拟[J].含能材料, 2013, 21(6):806-811. doi: 10.3969/j.issn.1006-9941.2013.06.023

    Wu Shuangzhang, Gu Wenbin. Experimental and three-dimensional numerical simulation for influence of end constraint on jet formation of cuneiform cover cutter[J]. Chinese Journal of Energetic Materials, 2013, 21(6):806-811. doi: 10.3969/j.issn.1006-9941.2013.06.023
    [9] 李成兵, 裴明敬, 沈兆武.高速杆式弹丸三维数值模拟[J].高压物理学报, 2007, 21(2):165-172. doi: 10.3969/j.issn.1000-5773.2007.02.008

    Li Chengbing, Pei Mingjing, Shen Zhaowu. Three-dimensional numerical simulation of the high velocity rod-shaped projectile[J]. Chinese Journal of High Pressure Physics, 2007, 21(2):165-172. doi: 10.3969/j.issn.1000-5773.2007.02.008
    [10] 武双章, 顾文彬, 刘建青, 等.椭圆形罩LSC侵彻钢锭试验与三维数值模拟[J].含能材料, 2014, 22(5):600-606. http://d.old.wanfangdata.com.cn/Periodical/hncl201405006

    Wu Shuangzhang, Gu Wenbin, Liu Jianqing, et al. Test and three-dimensional numerical simulation of penetrating steel ingot by linear shaped charge with oval cover[J]. Chinese Journal of Energetic Materials, 2014, 22(5):600-606. http://d.old.wanfangdata.com.cn/Periodical/hncl201405006
    [11] 武双章, 顾文彬, 李旭锋, 等.单角度楔形罩LSC对X光底片盒影响的试验研究[J].工程爆破, 2014, 20(4):1-5. doi: 10.3969/j.issn.1006-7051.2014.04.001

    Wu Shuangzhang, Gu Wenbin, Li Xufeng, et al. Experimental study on the effect of linear shaped charge with single angle cuniform cover on X-ray cassette[J]. Engineering Blasting, 2014, 20(4):1-5. doi: 10.3969/j.issn.1006-7051.2014.04.001
  • 期刊类型引用(8)

    1. 蒋八运,程扬帆,李世周,钱家祺,韦箫,薛茹君. 环氧丙烷/空气混合物气-液两相燃爆特性. 含能材料. 2023(07): 699-706 . 百度学术
    2. 李峰,张晨雨,王悦,王博,张梦雨,荆亚东. 20L球型爆炸装置气液输送管段结构的优化设计. 高压物理学报. 2023(04): 178-189 . 百度学术
    3. 臧小为,吴峰,虞浩,吕启申,潘旭海,蒋军成. 20L近球形容器中甲醇喷雾液滴爆炸特性实验研究. 中国安全生产科学技术. 2021(11): 25-31 . 百度学术
    4. 刘雪岭,张奇. 预点火湍流对正戊烷云雾爆炸参数的影响. 爆炸与冲击. 2019(03): 4-13 . 本站查看
    5. 张成均,白春华. 基于20 L球罐的多相混合物扩散模拟. 中国安全生产科学技术. 2019(04): 52-58 . 百度学术
    6. 白春华,张成均,刘楠,姚宁. 环境温度对多相混合物爆炸特性影响的实验研究. 高压物理学报. 2019(04): 182-188 . 百度学术
    7. 吕启申,臧小为,潘旭海,马鹏,虞浩,蒋军成. 温度和浓度对甲醇喷雾爆炸特性参数的影响. 爆炸与冲击. 2019(09): 149-157 . 本站查看
    8. 郑秋雨,孙永强,王旭,蓝真亮. 可燃液体雾滴粒径分布对其引燃性的影响. 安全与环境学报. 2018(05): 1779-1783 . 百度学术

    其他类型引用(2)

  • 加载中
图(7) / 表(2)
计量
  • 文章访问数:  4101
  • HTML全文浏览量:  1196
  • PDF下载量:  444
  • 被引次数: 10
出版历程
  • 收稿日期:  2014-06-27
  • 修回日期:  2014-10-16
  • 刊出日期:  2016-05-25

目录

/

返回文章
返回