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

黏土中爆炸成坑地冲击耦合效应实验研究

时本军 李杰 郭纬 徐天涵 徐小辉 李干 蒋海明

罗棕木, 李克, 陈浩, 张玉武, 梁民族, 林玉亮. 爆炸冲击波作用下假人头部加速度响应测试与损伤分析[J]. 爆炸与冲击, 2024, 44(12): 121435. doi: 10.11883/bzycj-2024-0242
引用本文: 时本军, 李杰, 郭纬, 徐天涵, 徐小辉, 李干, 蒋海明. 黏土中爆炸成坑地冲击耦合效应实验研究[J]. 爆炸与冲击, 2023, 43(6): 065101. doi: 10.11883/bzycj-2022-0445
LUO Zongmu, LI Ke, CHEN Hao, ZHANG Yuwu, LIANG Minzu, LIN Yuliang. Acceleration response test and damage analysis of dummy head under explosion shock wave[J]. Explosion And Shock Waves, 2024, 44(12): 121435. doi: 10.11883/bzycj-2024-0242
Citation: SHI Benjun, LI Jie, GUO Wei, XU Tianhan, XU Xiaohui, LI Gan, JIANG Haiming. Experimental study on explosion cratering and coupled ground shock in clay[J]. Explosion And Shock Waves, 2023, 43(6): 065101. doi: 10.11883/bzycj-2022-0445

黏土中爆炸成坑地冲击耦合效应实验研究

doi: 10.11883/bzycj-2022-0445
基金项目: 国家自然科学基金(51679249, 52279120, 12072371, BK20221528)
详细信息
    作者简介:

    时本军(1993- ),男,博士研究生, benjunshi@163.com

    通讯作者:

    李 杰(1981- ),男,博士,教授, lijierf@163.com

  • 中图分类号: O389

Experimental study on explosion cratering and coupled ground shock in clay

  • 摘要: 为获得黏土中爆炸成坑体积与耦合地冲击能量的关系,采用10.5 g TNT厘米级球形炸药球作为爆炸源,在1500 mm×1490 mm分层式爆炸装置中开展了变埋深条件下的爆炸实验,利用3D扫描设备记录不同埋深下弹坑的真实体积,并通过动态土压力传感器测得地冲击传播衰减规律。实验结果表明:随埋深增大,耦合至黏土中的有效地冲击能量急剧增大,装药中心下方的有效弹坑体积与耦合至黏土中的有效地冲击能量基本呈正比关系,当装药比例埋深与封闭爆炸条件下爆炸空腔半径相当时,耦合至黏土中的有效地冲击能量基本达到饱和。结合实验结果给出了黏土中爆炸耦合地冲击能量分配随装药比例埋深的变化规律,建立了地下爆炸等效封闭当量计算方法,为地下工程抗爆设计提供了理论依据。
  • 在军事行动、恐怖袭击、意外事故等情形中,爆炸伤已成为人员伤亡的主要模式[1]。除了眼眶、听器及肺脏等空腔脏器容易遭受爆炸冲击波的直接损害外,爆炸引起的创伤性脑损伤(blast-induced traumatic brain injury, bTBI)的比例显著上升,引起了广泛关注[2-5]。2000-2020年,约43万名美国军人确诊创伤性脑损伤[6];在国际恐怖袭击事件中,爆炸袭击为最主要的形式之一,约占29%[7];2019年江苏响水“3·21”化工厂特大爆炸事故造成了78人遇难,640人受伤[8]。此外,bTBI不仅导致脑组织挫伤、轴突弥漫性损伤等原发性损伤,也会引起焦虑、易怒、睡眠障碍和抑郁等继发性损害,严重影响伤者的身心健康。因此,研究爆炸冲击波的头部致伤机制,通过可测物理量预测爆炸事件中头部损伤等级具有重要意义。

    原发性脑损伤的力学机制复杂且影响因素较多,研究人员提出了爆炸冲击波直接作用[9-10]、颅骨弯曲变形[11-12]、空化效应[13-14]、胸腔压迫[15]和爆炸产生的加速度[16-18]等多种作用机制,利用模型试验和仿真计算的方法研究爆炸脑损伤的致伤机理。其中,宏观尺度的加速度损伤机制重点关注头部的线性加速度和旋转加速度,由于惯性作用,脑组织和颅骨间发生相对运动,从而导致脑组织挫伤和轴突弥漫性损伤。Gullotti等[16]通过大鼠试验发现,在相同爆炸条件下,头部加速度降低使大鼠神经元损伤减轻。Mao等[17]通过有限元仿真手段,将大鼠头部定义为刚性材料以获取加速度引起头部的力学响应,发现加速度单独作用可产生显著的颅内压。除了动物试验和计算模型外,Du等[18]和Sarvghad-Moghaddam等[19]建立了人体颅脑有限元模型研究以加速度为主导的损伤模式,发现爆炸加速度与颅内压、最大剪应力、最大主应变以及轴突应变率有密切联系。

    虽然爆炸加速度引起的bTBI的确切机制尚不清楚,但其可作为评价爆炸脑损伤严重程度的预测指标。在运动学和汽车碰撞领域,许多损伤评估指标基于头部的运动加速度,如基于线性加速度的头部损伤标准(head injury criterion, HIC)[20]、考虑脑损伤耐受力的广义加速度模型(generalized acceleration model for brain injury tolerance, GAMBIT)[21]、头部撞击功率(head impact power, HIP)[22]、基于角速度和角加速度的大脑旋转损伤标准(brain rotation injury criterion, BRIC)[23]和基于累积角加速度的旋转损伤标准(rotation injury criterion, RIC)[24]。Lockhart等[25]利用爆炸试验数据中的峰值加速度和HIC,对爆炸载荷下的人体数值模型进行了验证,结果表明爆炸载荷可产生显著的加速度响应和潜在的头部损伤。Shridharani等[26]开展了一系列激波管加载猪头试验,发现头部峰值加速度和HIC值的对数形式与峰值入射超压有良好的线性相关性。Singh等[27]通过建立爆炸-头部耦合模型对比头部矢状面和横切面的头部响应,发现横切面的运动学参数以及HIC值略高于矢状面。

    综上所述,人体头部作为极易遭受爆炸伤害的靶目标,其爆炸致伤机制及损伤预测判定备受关注。然而,以往的研究大多是采用动物或头部模型试验,忽略了人体总质量以及头-颈连接对加速度测量的影响。本研究旨在运用标准人体参数的假人模型开展静爆试验,通过测量不同比例距离下假人头部经历的加速度响应,探索头部加速度与爆炸参数的内在联系。基于峰值线性加速度、HIC和HIP损伤评估指标和冲击波对人员损伤判据,对比分析试验工况下的头部损伤风险等级,评价基于加速度参数的头部损伤预测方法的适用性和有效性,以期为爆炸加速度引起的头部损伤评估提供参考。

    静爆试验选用消防训练假人作为被测对象以模拟标准人体,假人净身高165 cm、质量60 kg,表面为防爆PU材质,内部由碎布和沙袋填充,假设为等密度体。假人头部、躯干部分和腿部分别近似为18 cm×21 cm、34 cm×60 cm、17 cm×84 cm的圆柱体,假人头部质量约占整体质量的7.2%,与成年人头部质量占比6.8%[28]接近。

    爆炸试验中,将假人模型放置于高度可调节的定位装置上,以支撑模型保持自由站立姿态且使待测部位所在平面尽可能与地面保持垂直,如图1所示,铝合金支架高170 cm、宽70 cm,横梁高度115 cm与假人模型手臂高度相当。当爆炸冲击波作用时,支架允许模型自由向后跌落,模拟假人与地面真实接触条件,使得冲击波作用下假人模型的自由运动模式不被干扰。

    图  1  假人模型
    Figure  1.  Dummy model

    假人头部加速度的测量选用扬州科动电子型号KD1020的PE型加速度传感器,测量范围为±500g,频率响应为0.5~4 kHz,工作温度为−40~120 ℃,满足试验工况下的加速度量程及环境温度。为了确保加速度传感器与假人头部刚性附着,利用喉箍标准件将固定在三向块上的加速度传感器栓接在头部冠状面背侧以监测头部的整体运动,如图2(a)所示,其中假人面向爆心为x轴方向,假人与爆心连线的水平垂直方向为y轴方向,与地面垂直的假人体长方向为z轴方向。

    图  2  测量设备及其布设情况
    Figure  2.  Measuring instrument and layout

    自由场冲击波超压的测量选用PCB型号137B21B的IEPE型压力传感器,测量范围为±6.7 MPa、分辨率为0.059 kPa、工作温度为−73~135 ℃。将压力传感器布设在假人与爆心相同距离的空旷位置,防止冲击波反射和绕射影响测量信号,如图2(b)所示。所有传感器均通过屏蔽电缆线与东华测试DH5960型数据采集仪连接,数据采集仪的采样频率设置为1 MHz,触发模式为负延时500 ms的通道触发,以保证获取的信号完整可靠。

    本次爆炸试验中采用TNT球形装药方式,其炸高设定为1.5 m,与加速度传感器和压力传感器固定安装高度一致,确保测试数据准确性。表1展示了8发静爆试验工况,TNT药量分别为1、2、3和4 kg,每种药量均重复进行2次试验,以增强试验结果的可靠性。基于爆炸冲击波对人员损伤的超压判据[29]和峰值超压经验公式[30],设置假人模型和自由场压力测点距爆心均为4.2 m(图2(c)),随着TNT当量的逐级递增,理论超压分别为0.045、0.068、0.089和0.108 MPa,对应的人员损伤程度为出现听觉器官损伤等的轻度损伤[29]

    表  1  试验工况
    Table  1.  Test conditions
    试验编号 TNT质量/
    kg
    测试距离/
    m
    比例距离/
    (m∙kg−1/3)
    理论超压/
    MPa
    1#,4# 1 4.2 4.20 0.045
    2#,7# 2 4.2 3.33 0.068
    3#,5# 3 4.2 2.91 0.089
    6#,8# 4 4.2 2.65 0.108
    下载: 导出CSV 
    | 显示表格

    TNT炸药由雷管引爆后,产生的爆炸冲击波迅速在空气中传播,掠过自由场传感器的敏感面并通过数据采集仪输出压力信号,图3给出了不同TNT当量下PCB传感器获取的爆炸自由场压力时程曲线。由图3可见,爆炸试验测得的实际超压在0.044~0.120 MPa之间,与预估超压范围接近。可以看出,在相同TNT当量下,传感器测量的峰值超压、正压持续时间接近,表明传感器测量的一致性和可靠性较好。同时,随着当量增加,冲击波的到达时间越早、峰值超压越高,峰值超压均值由0.046 MPa急剧增大到0.118 MPa。此外,可以观察到每一发试验测量的超压时程曲线中均存在双峰值的情况。这是由于,冲击波传播到地面后发生反射,反射冲击波经过压力传感器形成第2道压力峰值,Qi等[31]在试验中也观察到了这一现象,在试验中将压力传感器的敏感面朝上以尽量降低反射波的影响。

    图  3  不同TNT当量下自由场冲击波超压时间演化曲线
    Figure  3.  Overpressure-time histories of free-field shock waves under different TNT masses

    根据萨道夫斯基经验公式[30],空气中爆炸冲击波峰值压力可表示为:

    pm=0.084(3WL)+0.27(3WL)2+0.7(3WL)3
    (1)

    式中:pm为爆炸冲击波峰值压力,MPa;W为TNT当量,kg;L为被测点与爆心的距离,m。

    图4展示了爆炸试验中测得的冲击波峰值压力与理论公式的对比。可以看出,试验值与理论曲线吻合程度较高,测得的试验值总体高于理论值,与理论曲线的最大误差δ不超过13.3%,进而验证了自由场压力测试的可靠性。

    图  4  自由场冲击波超压随比例爆距的变化
    Figure  4.  Variation of free-field shock wave overpressure with scaled distance away from explosion center

    爆炸冲击波压力引起假人整体和局部运动的加速度,在毫秒时间尺度诱发原发性脑损伤;此后,由于人体惯性运动引起头部旋转和撞击地面发生的继发性脑损伤的时间尺度分别在百毫秒和秒量级[32],在本次试验研究中重点关注由加速度主导的原发性脑损伤。在汽车工业领域,头部损伤标准(HIC)由1 650 Hz以下的低频加速度决定,但爆炸载荷往往能够引起更高频率的加速度响应且对头部损伤影响更大[33]。根据加速度传感器4000 Hz的频率响应滤波加速度信号如图58所示。由图58可以看出,模型头部的加速度响应为连续的振荡信号,在相同比例距离下,xyz等3个方向的加速度振荡趋势一致,且峰值加速度值接近,表明重复试验结果的一致性较好。随着TNT当量的增加,峰值加速度有上升的趋势。此外,经过分析自由场冲击波测量结果,在本文所开展的试验工况下冲击波掠过假人头部的时间非常短(在0.5 ms以内),而从加速度曲线存在的明显周期性振荡(周期在15 ms左右)来看,试验测量的加速度反映的是假人头部的整体加速度。

    图  5  1 kg TNT的加速度时程曲线
    Figure  5.  Acceleration time histories of 1 kg TNT
    图  6  2 kg TNT的加速度时程曲线
    Figure  6.  Acceleration time histories of 2 kg TNT
    图  7  3 kg TNT的加速度时程曲线
    Figure  7.  Acceleration time histories of 3 kg TNT
    图  8  4 kg TNT的加速度时程曲线
    Figure  8.  Acceleration time histories of 4 kg TNT

    图9对比了不同TNT当量下各轴向的峰值加速度,可以看出,各轴向峰值加速度随TNT当量增大而急剧升高,x轴方向峰值加速度由最小16.29g升高到最大70.11g,除4 kg TNT当量下x轴的峰值加速度差异较大外,相同爆炸条件下同轴向加速度的标准差均在4.71g以下。为理清爆炸冲击波作用下头部加速度与入射超压之间的联系,提取各个工况下测得的峰值加速度和峰值超压并绘出关系图如图10所示。可以发现,排除异常数据后,xz轴方向的峰值加速度与峰值超压的线性相关性较强(R2>0.92),而y轴方向的相关性较弱。此外,正对爆心方向(x轴)峰值加速度关于超压的敏感性最强,表明爆炸冲击波作用在假人头部时,对迎爆面方向的加速度贡献最大。

    图  9  不同TNT质量下各轴向峰值加速度
    Figure  9.  Peak axial accelerations at different TNT masses
    图  10  各轴向峰值加速度与峰值超压关系
    Figure  10.  Relation between peak axial acceleration and peak overpressure

    在当前研究中,爆炸环境下加速度引起的头部损伤评估缺乏统一的判定标准,本文中结合爆炸试验数据,通过峰值线性加速度、头部损伤标准(HIC)、头部撞击功率(HIP)等基于平动加速度的损伤评估指标开展头部损伤分析。

    线性加速度是用于头部损伤风险评估和防护装备安全性评估的最基本运动学参数,峰值线性加速度指标来源于体育赛事中的损伤研究,该指标是忽略加速度持续影响时间、考虑所有方向的总加速度,但峰值线性加速度的评判标准给出了持续时间的限制,持续时间一般为3 ms以上:

    am={a2x+a2y+a2z}max
    (2)

    式中:axayaz为沿假人头部局部坐标系三轴的线加速度分量。头部严重损伤阈值am为200g(持续时间长于2 ms)[20]

    HIC是基于线性加速度的头部损伤评估指标,区别于仅考虑加速度峰值的判据,其考虑了加速度场的持续作用,记HIC值为γ,可由下式[20]计算:

    γn={(t2t1)[1t2t1t2t1a(t)dt]2.5}max
    (3)

    式中:n为积分的持续时间,t1t2分别为积分的起止时间,a(t)为加速度关于时间的函数。γ取加速度积分的最大值,本研究中n取15 ms,头部严重损伤阈值为700[20]

    HIP是一种综合运动头部损伤指标,其假设头部损伤严重程度与头部受到的撞击功率有关,定义为具有质量和转动惯性刚体的动能变化率,HIP值记为η,可由下式[22]计算:

    η=maxaxdt+mayaydt+mazazdt+Ixxαxαxdt+Iyyαyαydt+Izzαzαzdt
    (4)

    式中:η为HIP计算值,αxαyαz为同坐标系下三轴的角加速度分量,m为假人头部质量,IxxIyyIzz为三轴的转动惯量。由式(4)可知,η是一个随时间变化的物理量,损伤评估阈值取最大值ηm,头部严重损伤阈值为30 kW[22]

    根据加速度测试结果与假人模型头部尺寸(头部质量m=4.32 kg, Ixx=0.0246, Iyy=0.0246, Izz=0.0175),表2总结了不同工况下各损伤评估指标的特征值,amγ15ηm的范围分别为22.51g~71.41g, 4.66~81.46, 0.91~6.26 kW。本次工况下计算的各指标特征值远小于严重损伤阈值,因此根据文献[22, 34-35]列举了各损伤评估指标关于轻度创伤型脑损伤(mTBI)的损伤风险阈值如表3所示。

    表  2  不同工况下各损伤评估指标的特征值
    Table  2.  Characteristic values of injury assessment indii under different conditions
    试验编号 TNT质量/kg 比例距离/(m∙kg−1/3) 峰值超压/MPa am/g γ15 ηm/kW
    1# 1 4.20 0.0486 22.77 4.66 1.30
    2# 2 3.33 0.0693 39.70 28.30 2.60
    3# 3 2.91 0.0984 58.68 45.92 3.97
    4# 1 4.20 0.0443 22.51 6.81 0.91
    5# 3 2.91 0.0956 64.62 51.14 4.35
    6# 4 2.65 0.1161 71.41 81.46 6.26
    7# 2 3.33 0.0775 38.54 23.98 2.50
    8# 4 2.65 0.1199 61.77 57.53 5.27
    下载: 导出CSV 
    | 显示表格
    表  3  各损伤评估指标mTBI的损伤风险阈值
    Table  3.  Injury risk threshold of mild traumatic brain injury for each injury assessment index
    风险概率/% am/g 来源 风险概率/% γ15 来源 风险概率/% ηm/kW 来源
    5 39 文献[34] 5 10 文献[35] 5 4.7 文献[22]
    10 50 10 52 10 6.7
    25 66 25 136 25 9.7
    50 82 50 235 50 12.8
    80 106 75 333 75 15.9
    下载: 导出CSV 
    | 显示表格

    图11为各工况下总加速度的时程曲线,可以发现总加速度的第一峰值均发生在5 ms以内,而后加速度经历先下降后陡然上升,在10 ms左右达到第2峰值,这与自由场冲击波超压的规律一致,第2峰值是地面反射冲击波所致。此外,除比例距离4.2 m/kg1/3条件下的最大加速度为第1峰值外,其余工况下的第2峰值为最大值且持续时间更长。图12为峰值线性加速度随峰值超压的变化,从图中可以看出,总加速度峰值亦与峰值超压呈线性关系,峰值超压小于0.07 MPa时的轻度头部损伤风险低于5%,峰值超压大于0.09 MPa时的轻度头部损伤风险高于10%,本次试验工况下最高损伤风险高于25%。

    图  11  不同工况下总加速度时程曲线
    Figure  11.  Time histories of resultant acceleration at different conditions
    图  12  峰值线性加速度与峰值超压关系
    Figure  12.  Relation between peak linear acceleration and peak overpressure

    图13展示了HIC15值与比例距离的关系,HIC15与比例距离按幂函数拟合效果较好。可以看出,随着比例距离的降低快速增大,比例距离由4.20 m/kg1/3减小至2.65 m/kg1/3时,HIC15平均增大了12倍,最大值81.46远低于FMVSS目录中规定的头部严重损伤阈值700[20]。在相同比例距离下,以HIC15的平均值为参考进行误差分析,误差在±5.0%以内,表明相同工况下的一致性较好。由图14可知,HIC15随自由场峰值超压呈指数形式上升,峰值超压在0.05~0.10 MPa之间时,轻度头部损伤风险为5%,对应的超压损伤判据为出现耳膜破裂;当峰值超压超过0.103 5 MPa时,轻度头部损伤风险为10%,对应的超压损伤判据为50%耳膜破裂。

    图  13  HIC15与比例距离的关系
    Figure  13.  Relation between HIC15 and scaled distance
    图  14  HIC15与峰值超压的关系
    Figure  14.  Relation between HIC15 and peak overpressure

    根据式(4)绘制出HIP的时程曲线如图15所示,HIP曲线随时间总体呈先上升后下降的趋势,且相同工况下的吻合程度较高,HIP峰值时间与总加速度时程曲线的第二峰值时间相近。根据η最大值HIPm随峰值超压的变化关系(图16)可以发现,HIPm随峰值超压呈二次项上升的关系。当峰值超压小于0.1 MPa时,HIPm均未达到5%轻度损伤风险阈值4.7 kW,本次工况下HIPm最大为6.26 kW,亦低于10%轻度损伤风险阈值。与冲击波超压损伤判据相比,HIP评估指标在针对峰值超压低于0.12 MPa的爆炸环境下,表现出较低的灵敏度以及风险等级。

    图  15  不同TNT当量下HIP时程曲线
    Figure  15.  Time histories of HIP at different TNT masses
    图  16  HIPm与峰值超压的关系
    Figure  16.  Relation between HIPm and peak overpressure

    根据冲击波超压对人员的损伤判据[29],死亡率1%和50%的超压范围分别为0.138~0.241 MPa和0.276~0.345 MPa,结合上述3种损伤评估指标与峰值超压的拟合公式,绘制出不同指标预测此超压范围内的头部损伤风险如图17所示。可以看出,在死亡率1%的超压范围内,峰值线性加速度指标预测的轻度损伤风险由50%快速上升至80%,而HIC和HIP指标预测的轻度损伤风险经历25%后缓慢上升至75%。这表明峰值线性加速度指标预测此超压范围内的损伤风险相对较高。当3种评估指标达到头部严重损伤阈值时,对应的峰值超压依次为0.322、0.300和0.332 MPa,均分布在死亡率50%的超压范围内,其中HIC指标下的峰值超压最低,表明HIC指标预测头部严重损伤的敏感性最强,与超压损伤判据的一致性较高。

    图  17  不同损伤评估指标的预测情况
    Figure  17.  Predictions of different injury assessment indexes

    通过球形TNT装药静爆试验,研究了爆炸冲击波作用下假人模型头部的加速度响应及其与冲击波超压的内在联系,基于峰值线性加速度、HIC和HIP评估指标分析了爆炸加速度引起的头部损伤风险,得到的主要结论如下。

    (1)自由场爆炸产生的冲击波峰值超压的试验值和理论曲线吻合较好,试验值略高于理论值但误差在13.3%以内。随着TNT当量的增加,x轴方向峰值加速度由最小16.29g升高到最大70.11g。各方向的峰值加速度与峰值超压呈良好的线性关系,x轴方向峰值加速度关于超压的敏感性最强,爆炸冲击波对迎爆面方向的加速度贡献最大。

    (2)试验工况下,总加速度峰值随峰值超压线性增大,轻度损伤风险最高为25%;HIC15随峰值超压呈指数形式上升,0.05~0.10 MPa超压的轻度损伤风险为5%,超过0.1035 MPa超压的轻度损伤风险为10%;HIPm随峰值超压呈二次项上升关系,当峰值超压低于0.1 MPa时,头部轻度损伤风险低于5%,表现出较低的灵敏度及风险等级。

    (3)在更严重损伤程度的超压范围内,各损伤评估指标与峰值超压的拟合关系具有良好的适用性。在0.138~0.241 MPa的超压范围内,HIC和HIP指标的轻度损伤风险经历25%后缓慢上升至75%,峰值线性加速度指标预测的损伤风险相对较高。加速度类型指标预测头部严重损伤对应的峰值超压依次为0.322、0.300和0.332 MPa,其中HIC指标下的峰值超压更接近50%死亡率的起始阈值,表明其预测头部严重损伤的效果较好。

  • 图  1  黏土的粒径级配曲线

    Figure  1.  Particle size gradation curve of clay

    图  2  分层式爆炸容器

    Figure  2.  Layered explosion vessel

    图  3  球形爆源

    Figure  3.  Spherical explosion source

    图  4  药球位置及传感器布置

    Figure  4.  Drug package location and sensor arrangement

    图  5  测量采集记录图

    Figure  5.  Measurement acquisition record diagram

    图  6  成坑的俯视与剖面图

    Figure  6.  Overhead view and profile of the crater

    图  7  3D扫描成坑形貌

    Figure  7.  3D scanning into crater morphology

    图  8  可视弹坑实测值与ConWep计算值的对比

    Figure  8.  Comparison of the measured values of the visible burst crater with those computed by ConWep

    图  9  随装药埋深增加有效弹坑的变化规律

    Figure  9.  The change rule of effective crater with the increase of charge depth

    图  10  不同装药比例埋深时实测地冲击应力波形

    Figure  10.  Measured stress waveform of clay at different scaled buried depths

    图  11  不同比例爆心距离处地冲击应力峰值比值随装药比例埋深增加的变化情况

    Figure  11.  Peak of the ground impact stress varied with scaled butied depth charge at different scaled blast center distances

    图  12  不同装药比例埋深时地冲击应力峰值衰减曲线

    Figure  12.  Peak groud impact stress attenuation with scaled blast center distance at different scaled buried depths of charge

    图  13  地下爆炸的三种形式

    Figure  13.  Three types of underground explosion

    图  14  地冲击能量耦合系数随比例埋深的变化

    Figure  14.  Variation of ground impact energy coupling coefficient with scaled burial depth

    图  15  黏土中耦合系数随装药比例埋深变化关系

    Figure  15.  Relationship between coupling coefficient and scaled burial depth in clay

    表  1  装药埋深(h)及爆心距(R)设计

    Table  1.   Design of burial depth of charge (h) and burst core distance (R)

    工况h/(m·kg−1/3)R/(m·kg−1/3)
    1#2#3#4#5#
    1−0.0560.7991.4272.2953.2744.000
    200.7991.4272.2953.2744.000
    30.140.7991.4272.2953.2744.000
    40.370.7991.4272.2953.2744.000
    50.550.7991.4272.2953.2744.000
    61.190.7991.4272.2953.2744.000
    71.460.4791.1552.0683.0463.772
    下载: 导出CSV

    表  2  不同埋深条件下弹坑尺寸数据

    Table  2.   Size data of craters under different burial depths

    工况h/
    (m·kg−1/3)
    rv /
    (m·kg−1/3)
    dv /
    (m·kg−1/3)
    Vv /
    (m3·kg−1)
    ra /
    (m·kg−1/3)
    d/
    (m·kg−1/3)
    V/
    (m3·kg−1)
    1−0.0560.2610.2470.0350.1790.2510.017
    200.2800.2920.0480.1840.2600.018
    30.140.5800.5390.3800.2890.3790.066
    40.370.6820.6850.6670.3400.4160.101
    50.550.8960.7761.3040.4530.4590.197
    61.190.6150.5020.3970.4780.4780.228
    71.46000.0000.4790.4790.230
    注:rvdvVv分别为可视弹坑的半径、深度和体积,radV为有效弹坑的半径、深度和体积。
    下载: 导出CSV

    表  3  黏土中各比例埋深下地冲击应力峰值数据

    Table  3.   Subsurface impact stress peak data of each proportion buried depth in clay

    工况h/
    (m·kg−1/3)
    σpk/MPa工况h/
    (m·kg−1/3)
    σpk/MPa
    1#2#3#4#5#1#2#3#4#5#
    1−0.0560.0470.0150.0190.0110.00250.550.3580.1550.1150.0860.022
    200.0500.0380.0210.0150.00661.190.3340.1560.1260.1150.032
    30.140.1410.0920.0740.0540.01571.461.0800.1960.1240.1000.024
    40.370.1840.1000.0700.0520.014(0.389)(0.166)(0.047)(0.023)(0.015)
    注:(1)1~6炮次1#、2#、3#、4#、5#测点比例距离分别为0.799、1.427、2.295、3.274、4.000;(2)第7炮次1#、2#、3#、4#、5#测点比例距离分别为0.479、1.155、2.068、3.046、3.772;(3)第7炮括号内数据为将第7炮次数据在比例距离0.799、1.427、2.295、3.274、4.000处换算数据。
    下载: 导出CSV

    表  4  不同装药比例埋深条件下拟合参数

    Table  4.   Fitting parameters with different scaled buried depths of charge

    h/(m·kg−1/3)nAˉnA'
    −0.0561.160.0351.14
    0.031
    01.050.0380.041
    0.141.090.1200.123
    0.371.070.1240.148
    0.551.230.2610.258
    1.191.160.2670.258
    1.461.190.2690.260
    下载: 导出CSV
  • [1] ADUSHKIN V V, SPIVAK A. Underground explosions [M]. Lexington: Weston Geophysical Corp, 2015: 431-480.
    [2] Department of the Army. Fundamentals of protective design for conventional weapons: TM5-855-1 [R]. Washington: US Department of the Army, 1986.
    [3] BALADI G Y, NELSON I. Ground shock calculation parameter study: S-71-4 [R]. Mississippi: US Army Waterway Experimental Station, 1974.
    [4] JAMES K. Cense explosion test program: report 2, case 2, explosions in soil [R]. Mississippi: US Army Waterway Experimental Station, 1977.
    [5] 乔登江. 地下核爆炸现象学概论 [M]. 北京: 国防工业出版社, 2002.
    [6] 梁霍夫 Γ M. 岩土中爆炸动力学基础 [M]. 刘光寰, 王明洋, 译. 南京: 工程兵工程学院, 1993.
    [7] YANKELEVSKY D Z, KARINSKI Y S, FELDGUN V R. Re-examination of the shock wave’s peak pressure attenuation in soils [J]. International Journal of Impact Engineering, 2011, 38(11): 864–881. DOI: 10.1016/j.ijimpeng.2011.05.011.
    [8] 穆朝民, 任辉启, 辛凯, 等. 变埋深条件下土中爆炸成坑效应 [J]. 解放军理工大学学报(自然科学版), 2010, 11(2): 112–116. DOI: 10.3969/j.issn.1009-3443.2010.02.003.

    MU C M, REN H Q, XIN K, et al. Effects of crater formed by explosion in soils [J]. Journal of PLA University of Science and Technology (Natural Science Edition), 2010, 11(2): 112–116. DOI: 10.3969/j.issn.1009-3443.2010.02.003.
    [9] 穆朝民, 任辉启, 李永池, 等. 变埋深条件下饱和土爆炸能量耦合系数的试验研究 [J]. 岩土力学, 2010, 31(5): 1574–1578. DOI: 10.16285/j.rsm.2010.05.006.

    MU C M, REN H Q, LI Y C, et al. Experiment study of explosion energy coupling coefficient with different burial depths in saturated soils [J]. Rock and Soil Mechanics, 2010, 31(5): 1574–1578. DOI: 10.16285/j.rsm.2010.05.006.
    [10] 施鹏, 邓国强, 杨秀敏, 等. 土中爆炸地冲击能量分布研究 [J]. 爆炸与冲击, 2006, 26(3): 240–244. DOI: 10.11883/1001-1455(2006)03-0240-05.

    SHI P, DENG G Q, YANG X M, et al. Study on ground shock energy distribution of explosion in soil [J]. Explosion and Shock Waves, 2006, 26(3): 240–244. DOI: 10.11883/1001-1455(2006)03-0240-05.
    [11] 叶亚齐, 任辉启, 李永池, 等. 砂质黏土中不同深度爆炸自由场地冲击参数预计方法研究 [J]. 岩石力学与工程学报, 2011, 30(9): 1918–1923.

    YE Y Q, REN H Q, LI Y C, et al. Study of prediction of ground shock parameters in free field at different depths of burst in sandy clay [J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(9): 1918–1923.
    [12] 赵红玲, 侯爱军, 童怀峰, 等. 石灰岩中不同埋深爆炸自由场直接地冲击参数的预计方法 [J]. 爆炸与冲击, 2011, 31(3): 290–294. DOI: 10.11883/1001-1455(2011)03-0290-05.

    ZHAO H L, HOU A J, TONG H F, et al. Prediction method of the direct ground shock parameters of explosion at different buried depths in free field of limestone [J]. Explosion and Shock Waves, 2011, 31(3): 290–294. DOI: 10.11883/1001-1455(2011)03-0290-05.
    [13] 何翔, 吴祥云, 李永池, 等. 石灰岩中爆炸成坑和地冲击传播规律的试验研究 [J]. 岩石力学与工程学报, 2004, 23(5): 725–729. DOI: 10.3321/j.issn:1000-6915.2004.05.004.

    HE X, WU X Y, LI Y C, et al. Testing study on crater formed by explosion and propagation laws of ground shock in limestone [J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(5): 725–729. DOI: 10.3321/j.issn:1000-6915.2004.05.004.
    [14] JAMES K. Cense explosion test program: report 1, case 1, explosions in sandstone [R]. Mississippi: U S Army Waterway Experimental Station, 1977.
    [15] Department of the Army. Structures to resist the effects of accidental explosions: TM 5-1300 [R]. Washington: US Department of the Army, 1990.
    [16] 李晓军, 张殿臣, 李清献, 等. 常规武器破坏效应与工程防护技术 [M]. 洛阳: 总参工程兵科研三所, 2001.
    [17] 李重情, 穆朝民, 石必明. 变埋深条件下混凝土中爆炸应力传播规律的研究 [J]. 振动与冲击, 2017, 36(6): 140–145. DOI: 10.13465/j.cnki.jvs.2017.06.021.

    LI Z Q, MU C M, SHI B M. Investigate on shock stress propagation in concrete at different depths under blasting [J]. Journal of Vibration and Shock, 2017, 36(6): 140–145. DOI: 10.13465/j.cnki.jvs.2017.06.021.
    [18] LEONG E C, ANAND S, CHEONG H K, et al. Re-examination of peak stress and scaled distance due to ground shock [J]. International Journal of Impact Engineering, 2007, 34(9): 1487–1499. DOI: 10.1016/j.ijimpeng.2006.10.009.
    [19] 中华人民共和国住房和城乡建设部. GB/T 50123—2019 土工试验方法标准 [S]. 北京: 中国计划出版社, 2019: 16–59.

    Ministry of Housing and Urban-Rural Development of the People's Republic of China. GB/T 50123—2019 Standard for geotechnical testing method [S]. Beijing: China Planning Press, 2019: 16–59.
    [20] KRAUTHAMMER T. Modern protective structures [M]. Boca: CRC Press, 2008.
    [21] DAVID H W. User’s guide for microcomputer programs CONWEP and FUNPRO: applications of TMS-855-1, “fundamentals of protective design for conventional weapons”: AD-A195867 [R]. 1988.
    [22] 徐天涵. 钻地核爆地冲击效应与防护研究 [D]. 南京: 中国人民解放军陆军工程大学, 2021: 41–79.
    [23] HASKELL N A. Analytic approximation for the elastic radiation from a contained underground explosion [J]. Journal of Geophysical Research, 1967, 72(10): 2583–2587. DOI: 10.1029/JZ072i010p02583.
    [24] 王明洋, 李杰. 爆炸与冲击中的非线性岩石力学问题Ⅲ: 地下核爆炸诱发工程性地震效应的计算原理及应用 [J]. 岩石力学与工程学报, 2019, 38(4): 695–707. DOI: 10.13722/j.cnki.jrme.2018.1078.

    WANG M Y, LI J. Nonlinear mechanics problems in rock explosion and shock. Part Ⅲ: The calculation principle of engineering seismic effects induced by underground nuclear explosion and its application [J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(4): 695–707. DOI: 10.13722/j.cnki.jrme.2018.1078.
    [25] SHISHKIN N I. Seismic efficiency of a contact explosion and a high-velocity impact [J]. Journal of Applied Mechanics and Technical Physics, 2007, 48(2): 145–152. DOI: 10.1007/s10808-007-0019-6.
  • 加载中
图(15) / 表(4)
计量
  • 文章访问数:  491
  • HTML全文浏览量:  130
  • PDF下载量:  134
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-10-14
  • 修回日期:  2023-02-20
  • 网络出版日期:  2023-03-22
  • 刊出日期:  2023-06-05

目录

/

返回文章
返回