A study on the circumferential propagation law of the shock waves produced by the near ground dynamic explosion of cylindrical charge
-
摘要: 战斗部复杂的末弹道参数会影响近地爆炸冲击波的周向传播规律及对目标的毁伤程度,研究柱形装药近地爆炸冲击波传播规律对精确评估毁伤效能具有重要的工程意义。基于AUTODYN-3D软件对不同末弹道参数的柱形装药近地爆炸进行了数值模拟,通过对2个方向分别建模,获得了柱形装药近地爆炸下前、后、侧3个方向的冲击波压力数据;研究了战斗部的落速、落角、爆心高度和装药长径比4个参数对柱形装药近地爆炸冲击波传播的影响规律,分析了冲击波的演化过程、峰值压力和马赫杆高度。研究结果表明:静爆时,爆心高度是影响冲击波马赫杆高度的主要因素,落角与装药长径是影响马赫杆高度方向差异的主要因素;动爆时,能够增大周向马赫杆高度,前方最显著;另外,随着动爆速度的提升,前向冲击波峰值线性增大。正交优化的结果显示,4种变量中,动爆速度对柱形装药前方峰值压力极差最大,落角对后方峰值压力极差最大,爆心高度对马赫杆高度影响最大。通过研究柱形装药近地动爆冲击波周向传播规律,表明合理的调整装药参数和近地爆炸姿态对实现某方向的最大毁伤或减小超压伤害具有借鉴意义。Abstract: The complex terminal ballistic parameters of the warhead will affect the circumferential propagation law of the near ground explosive wave and the damage degree to the target. Studying the propagation law of the near ground explosive wave of the cylindrical charge has important engineering significance to accurately evaluate the damage efficiency. By using AUTODYN-3D software, the near ground explosion of cylindrical charge with different terminal ballistic parameters is simulated and calculated, and the data of shock wave pressure in the front, back and side directions produced by the near ground explosion of cylindrical charge are obtained by modeling in two directions respectively. Thus, the influences of four parameters, namely, the velocity of the battle group, the impact angle, the height of the explosion center and the ratio of the length to diameter of the charge, on the propagation of the shock wave produced by the near ground explosion of the cylindrical charge are studied. The evolution process of the shock wave, the peak pressure and the height of the Mach stem are analyzed. The results show that the height of the explosion center is the main factor affecting the height of the shock wave Mach stem during static explosion, and the impact angle and the length-to-diameter ratio of the charge are the main factors affecting the difference in the height direction of the Mach stem. During dynamic explosion, the height of circumferential Mach stem can be increased, especially in the front; in addition, the peak value of forward shock wave increases linearly with the increase of dynamic detonation velocity. The results of orthogonal optimization show that the dynamic detonation velocity has the largest range to the front peak pressure of the cylindrical charge among the four variables; the impact angle has the largest range to the rear peak pressure; and the height of explosion center has the greatest influence on the height of Mach stem. By studying the circumferential propagation law of the shock wave produced by near ground dynamic explosion of the cylindrical charge, the results show that reasonable adjustment of the charge parameters and the front of the near ground explosion can be used for reference to achieve the maximum damage or reduce the hyper-pressure damage in a certain direction.
-
Key words:
- cylindrical charge /
- near ground explosion /
- orthogonal optimization /
- Mach wave
-
国内学者对矿震机理[1-3]、矿震和采空区影响下围岩动力响应[4-5]、矿震和采空区诱发矿山地质灾害[6-7]等问题展开了大量研究[8-11]。贾宝新等[12]采用人工爆破方式模拟矿震进行微震定位方法研究,王帅等[13]采用数值微震试验的方式,基于到时时差开发层状介质中的震源定位方法,贾宝新等[14]利用钢球撞击块系夹层模型的方式模拟冲击动荷载,研究冲击地压下深部块系构造中摆型波的波形参数影响因素。潘一山等[15]、贾宝新等[16]从理论层面研究震波的三维传播规律及在非均匀介质条件下震波传播模型,考察微震监测实际接收到的信号与理论上的预期信号的差别,以服务于矿震震级精确计算和震源精确定位。
当前研究多集中在矿震、采空区等诱发的矿山地质灾害的机理、规律、预测预报和防治等方面,对采空区或采动影响下围岩冲击震动灾害的诱发问题、对矿震导致采空区围岩矿山地质灾害加剧问题等,均有相关研究。但在矿震与采空区耦合作用下,用直观的试验模型揭示各种矿山地质灾害的内在联系和统一致灾机理的相关研究较少。矿震与采空区耦合作用下,导致地表沉陷和地表建筑物破坏具有相关性,其作用机制具有紧密联系,这方面的研究缺乏试验支持。本文中建立力学模型开展了带孔混凝土圆板冲击试验,并将该试验模型用于揭示矿震与采空区同时存在时围岩的动力响应规律。
本文中的模型与文献[14]中的钢球撞击模型研究思路相似,不同于用真实爆破和数值计算来模拟矿震,而是采用一种折中方案——采用重锤敲击钢筋模拟矿震。该方法比真实爆破更易操作和实现、又比数值模拟计算更符合实际。本文中以文献[15-16]为模型设计的理论基础着重考察矿震影响下采空区围岩的动力响应规律和耦合致灾机理,侧重研究矿震发生后产生的结果,矿震震源机理研究的具体实现和应用。
1. 建模设计思路
通过现场调研勘测,以及在阜新某煤矿的实地考察,综合分析多方面收集的资料,验证了一些经验和共识:单一岩爆的发生,往往对事件发生点附近的人员和设备造成较大危害,但一般不引起其他较大的矿山地质灾害。废旧采空区在没有外界扰动的情况下,最常见的地质灾害就是地表沉陷,但这种沉陷的应变率很低,需要较长考察期才能观察到明显变化。但地震或者冲击地压发生后,处于准静态的采空区围岩经常会突然发生失稳,造成的后果包括采空区覆岩冒落、采空区地表突然塌陷、采空区中夹岩柱失稳、采空区地表建筑物损毁等,各种地质灾害同时发生,造成严重人员伤亡和财产损失。
工程技术人员在长期的实践和摸索中,逐渐形成了一些成熟的问题:矿震发生时,采空区覆岩冒落、地表沉陷加剧、中夹岩柱失稳、地表建筑物破坏加剧等地质灾害,是否与矿震和采空区的耦合作用有紧密联系?如果有,这些矿山地质灾害是否有内在关联和统一机理?为了回答这个问题,需要对矿区现场矿震与采空区耦合作用的真实环境、条件进行层层简化和抽象概括。排除干扰因素和次要矛盾,有针对性地合理简化和等效,提炼出只包含主要矛盾和关键因素的力学、试验模型,并参数化为变量,在变化中加以考察,即本文建模设计的基本思路。图 1形象地说明了这种抽象和简化的思路。
本文中,试验模型的本质是一种弹性动力学模型。设想一个足够大的薄板,在某一位置开设圆孔,距圆孔一定距离施加一个瞬时冲击胀缩动荷载。考察胀缩动荷载作用下,圆孔周边的动力响应规律。由于动荷载下带孔圆板的动力响应问题暂时还没有解析解,因此建立试验模型进行研究。
2. 试件制备和试验条件
混凝土等级为C40。用薄铁板制作模具,用PVC圆管预留圆孔,浇筑混凝土,并在圆板中心立带肋钢筋。海螺牌42.5水泥;设计坍落度:(100±30) mm;石子粒径5~25 mm,连续级配;中砂。设计水灰比为2.433:1。试拌并测定坍落度,如果坍落度不符合规范要求,则重新试拌,直到坍落度满足要求。需调整单位用水量、砂率时,需重新计算调整后的基准配合比并测表观密度。试件混凝土组分用量如表 1所示,调整后的混凝土基准配合比:水泥为438 kg/m3,水为180 kg/m3,砂为577 kg/m3,石子为1 227 kg/m3,外加剂为1.31 kg/m3。坍落度为115 mm,表观密度为2 423 kg/m3。
表 1 试件混凝土组分用量Table 1. Concrete mixture dosage of specimen试件类别 水泥用量/kg 水用量/kg 黄砂用量/kg 石子用量/kg 外加剂用量/kg 实心圆板 8.595 75 3.532 5 11.323 6 24.079 9 0.025 708 8 带孔圆板 8.179 72 3.361 5 10.775 6 22.914 4 0.024 464 4 本试验中爆破测振仪的传感器接收到的是加速度信号,需要将重锤的敲击动荷载换算为加速度,以方便后文分析比较。控制敲击动荷载大小的两个变量是重锤重量和下落高度。本试验重锤质量m=3 kg,连接细杆质量忽略不计,通过调整重锤下落的高度,来调节敲击动荷载量值的大小。重锤下落高度分别取h1=8 cm, h2=125 cm,撞击时间取t=0.01 s。撞击速度v=√2gh,则下落动能12mv2=mgh,由动量定理:Ft=mv-0,而F=ma,所以撞击加速度为:a=vt=√2ght。将上述条件代入可得,a1≈125 m/s2,a2≈500 m/s2。
3. 动荷载试验
3.1 带孔圆板冲击试验
试件制作过程、带孔圆板冲击试验过程、不带孔圆板对照组试验见图 2。采用固定高度落锤敲击法,敲击试验操作过程为用敲击支架上的重锤敲击钢筋,钢筋受到动荷载并传递给混凝土圆板,在混凝土圆板边缘布置4个加速度信号传感器,其中一个传感器布置在圆孔所在半径端部,与之相对的直径另一端布置一个传感器,与该直径垂直的直径两端各布置一个传感器。加速度信号监测及接收系统采用爆破测振仪。敲击动荷载分别为125和500 m/s2。每种动荷载、每个试件各敲击2个循环,每个循环敲击10次,结果取平均值。
3.2 实心圆板对比试验及相应理论模型
3.2.1 实心圆板对比试验
为了更加严谨合理,设置了不带孔圆板敲击试验。不带孔模型的大小、用料与带孔试件相同,重锤下落高度1.25 m。对比试验如图 2所示。对照组试验结果见表 2。试验中传感器的ar分量为径向加速度信号,传感器的aφ分量为切线方向加速度信号。径向方向中以从圆心向外为正,切线方向以逆时针方向为正。从表 2数据可以看出,环向响应加速度可以近似忽略不计;4个传感器接收到的径向响应加速度大小接近、方向一致。可初步认为,不带孔板的响应加速度分布与中央钢筋呈中心对称。
表 2 对比试验敲击响应加速度Table 2. Response acceleration of knocking test on contrast test传感器 第一轮 第二轮 aφ/(m·s-2) ar/(m·s-2) aφ/(m·s-2) ar/(m·s-2) 1 0 32.88 0 31.52 2 0 29.73 0 29.73 3 0 31.01 0 44.56 4 0 27.71 0 33.26 3.2.2 不带孔圆板对照组试验对应的点源球面纵波模型
如果弹性体具有圆球形的孔洞,而在孔洞内受到球对称的爆炸之类的作用,则由于对称,只可能发生径向位移uR,不可能发生切向位移,而且径向位移uR将只是径向坐标R和时间t的函数。这样,由孔洞向外传播的弹性波将是球对称的,即所谓球面波[17]。点源球面纵波微分方程的通解是[17]:
RΨ=f1(R−c1t)+f2(R+c1t) (1) 式中:Ψ是位移的势函数且uR=∂Ψ∂R;f1、f2为任意函数,表示沿径向传播的球面波,传播速度c1=√E(1−μ)(1+μ)(1−2μ)ρ。函数f1表示由内向外传播的球面波,适用于圆球形孔洞内受球对称动力作用时的情况[17]。根据理论模型,中心胀缩荷载产生的激励响应关于荷载中心对称,即与荷载中心相同距离的圆周上,响应加速度都大小相等。对照组试验结果与理论预期一致。
4. 试验数据与结果分析
4.1 单组试验结果
传感器编号如图 1中所示。表 3所示为试件敲击响应加速度。传感器1仅有径向加速度,且是负值;传感器2仅有径向加速度,且为正值。传感器1、2处无切向加速度,并非因圆孔阻隔,而是由于模型对称。传感器3、4加速度分布一致,也是因为对称。传感器3、4均存在正径向加速度,且传感器4还存在正切向加速度,传感器3还存在负切向加速度。根据表 2对比可以看出,传感器1~4位置只有径向加速度,无切向加速度,因为实心圆板关于钢筋对称。
表 3 对比试验敲击响应加速度Table 3. Response acceleration of knocking test on contrast test传感器 第一轮 第二轮 第三轮 第四轮 aφ/(m·s-2) ar/(m·s-2) aφ/(m·s-2) ar/(m·s-2) aφ/(m·s-2) ar/(m·s-2) aφ/(m·s-2) ar/(m·s-2) 1 0 -8.06 0 -5.77 0 -33.65 0 -25.47 2 0 6.55 0 2.92 0 48.48 0 41.46 3 -3.52 5.41 -3.77 3.04 -20.69 20.72 -42.03 50.86 4 3.92 7.93 12.82 10.21 36.55 16.55 23.87 34.98 4.2 单组试验结果分析
4.2.1 加剧采空区覆岩的沉陷和冒落
根据试验模型,震源到采空区连线上,采空区远离震源一侧,会出现指向震源的负加速度,这在试验角度提出了矿震动荷载加剧采空区沉陷的动力模型。考察一种理想情况,震源、采空区、地表依次自下向上竖直分布,符合浅部被采空、深部继续开采的常见开采结构。这种情况下,容易认为重力是采空区沉陷的主要原因,而矿震动荷载只是重力效应的催化剂。但根据本文的试验结果,震源、采空区、地表依次自下向上竖直分布,并不意味着沉陷仅仅由于重力效应,本试验中的对应点分别是震源、圆孔和圆孔远离震源侧,三者在水平直线排列的情况下即产生指向震源的负加速度,与重力无关。
4.2.2 对相邻采空区中夹岩柱的影响
震源与采空区呈左右分布,会对采空区中夹岩柱产生影响。当采空区侧帮发生矿震,水平临侧矿柱产生负加速度,采空区顶板产生剪切加速度。这时采空区顶板的沉陷和冒落,一方面是顶板下部岩体掏空和顶板岩体自身重力,矿震动荷载加剧了顶板围岩的应力扰动,使顶板围岩不稳定;另外采空区上部的剪切加速度也会对顶板的稳定性产生影响。当多采空区相邻,水平临侧继续开采,矿震动荷载会使相邻采空区中夹岩柱产生指向爆源的加速度,这对中夹岩柱稳定分析及加固支护有一定意义。
4.2.3 对采空区地表建筑物的影响
采空区与矿震动荷载同时存在会对地表建筑物的稳定安全产生不利影响。地震波研究中,面波比体波对地表建筑物的损伤更大;对于体波,横波比纵波对地表建筑物的损伤更大。归结原因,是因为地表的剪切震动对建筑物的危害更大。根据试验模型,当矿震震源与采空区在同一个水平面时,会在震源与采空区水平连线相垂直的上方(下方)产生水平方向的剪切震动,因此该试验模型可以解释矿山微震动荷载与采空区同时存在时,会对地表建筑物产生更不利的影响。当采空区与动载荷不在同一水平面上时,地表也会产生有所衰减的剪切加速度。本文考察的切向加速度,是震源与采空区同时存在的结果,采空区的存在改变了围岩的结构和应力状态,在矿震波等的综合作用下,会产生剪切加速度。因此可以这样认为,矿震通过体波和面波使地表产生剪切加速度,使地表建筑发生破坏,当有采空区存在时,在矿震和采空区的综合作用下,会增大地表的剪切加速度,因而会加剧地表建筑物的破坏。
4.3 试验数据的综合分析
图 3~6中的横坐标表示孔径,纵坐标表示响应加速度值。图 3~6分别表示传感器1~4位置处的响应加速度在不同动荷载下,随圆孔直径的变化趋势。图 3~4中两组点分别为敲击动荷载125和500 m/s2时,响应加速度数据。图 5~6中两组点分别为径向、环向响应加速度的变化数据。根据图 3可以看出,孔径很小时,传感器1处位置加速度在小范围出现正值,原因是圆孔较小,对加速度扰动的分流再聚合过程帮助不大,但从该部分加速度量值角度观察,圆孔对加速度扰动的阻隔作用仍然很明显。根据曲线整体,传感器1处位置的负加速度随孔径呈指数增长。根据图 4及前文分析,传感器2处位置加速度分布特征与不带孔圆板最为接近,在全部孔径范围内,加速度未出现负值,传感器2处位置加速度随孔径呈指数增长。
根据图 5所示,传感器3处的径、环向加速度随孔径均呈指数式增长。环向加速度的增长速度明显大于径向加速度。说明采空区越大,震源与采空区形成的动剪切荷载越大,对建筑物危害越大。
从图 5中可以看出,动荷载量值增加,径向、环向加速度按比例线性增加,随孔径变化规律不变。图 6与图 5反映的规律一致,更明显地,可以看出,孔径很小时径向位移大于环向位移,孔径增加,环向位移快速增加,并很快超过径向位移。原因在于,孔径较小时,圆孔边界对震动波的反射量有限,而直接传递过来的径向加速度相对较大,两者矢量合成,就使环向加速度不占主导,量值比径向加速度小。
综合上述结果:动荷载量值增加,关键位置加速度线性增加,且不同孔径下规律一致。随着孔径增加,径向、环向加速度在4个关键位置均呈指数增长,环向位移增长速率明显大于径向位移。
研究选取不同大小动荷载、不同孔径,孔径与动荷载的相对比例覆盖范围较大,基本包括实际矿震动荷载与实际采空区的相对比例。因此可以认为实际尺寸的采空区与实际大小的震源,所在围岩加速度分布规律相同,试验模型的结论能应用于工程实际。
5. 结论
为了研究矿震与采空区耦合致灾机理,进行了模型试验。共得到如下几个结论:
(1) 在采空区孔洞远离震源的一侧的半径上监测到指向圆心的与重力效应无关的负加速度。
(2) 浅部采空、深部继续开采并发生矿震,会在采空区上部产生向下的负加速度加剧沉陷和冒落。
(3) 水平方向同时存在多个采空区并发生矿震,采空区之间的岩柱会产生指向震源的加速度。
(4) 垂直震源与圆孔连线方向的圆板边缘处,不但产生径向加速度,而且产生切向加速度。当采空区与矿震动荷载处于同一水平面时,会在采空区与矿震发生区的上部地表产生水平剪切加速度。
(5) 随着圆孔直径的增大,圆板周围的径向加速度呈指数增长,当存在环向加速度时,环向加速度增长的趋势远大于径向加速度。因此采空区范围增大,对地表建筑物的影响将显著增长。
(6) 采空区尺寸效应和动荷载放缩效应,不影响上述试验结论,试验结论具有现实的工程意义。
-
表 1 TNT炸药状态方程参数
Table 1. Parameters of equation of state for TNT explosive
A/GPa B/GPa R1 R2 E0/GPa ω V 3.7377 3.747 4.15 0.9 0.06 0.3 1 表 2 空气状态方程参数
Table 2. Parameters of equation of state for air
ρ/(kg·m−3) γ e/(MJ·kg−1) 1.225 1.4 206.8 表 3 材料模型
Table 3. Material model
材料 状态方程 强度模型 失效模型 空气 理想气体 无 无 TNT JWL 无 无 表 4 测点压力的计算值与实验值对比
Table 4. Comparison between calculated values and experimental values on measuring point pressure
TNT质量/kg 测点距离/m 压力/MPa 误差/% 实验值 计算值 2.0 1.12 1.23 1.41 14.56 2.0 1.95 0.69 0.72 4.30 0.8 1.12 0.80 0.90 12.86 0.8 1.95 0.48 0.60 24.17 0.2 1.12 0.30 0.32 6.67 0.2 1.95 0.09 0.14 55.56 表 5 不同起爆条件下柱形装药的峰值压力
Table 5. Peak pressure of cylindrical charge under different initiation conditions
初始速度/(m·s−1) 2 ms时的压力/MPa 4 ms时的压力/MPa 激波加动爆 动爆 静爆 激波加动爆 动爆 静爆 600 1.04 1.04 0.98 0.40 0.40 0.39 1 200 0.98 0.98 1.00 0.39 0.39 0.40 表 6 不同地面条件下柱形装药的峰值压力和马赫杆高度
Table 6. Peak pressure and Mach stem height of cylindrical charge under different ground conditions
地面类型 起爆高度/m 峰值压力/MPa 马赫杆高度/mm 0.5 ms 1.0 ms 0.5 ms 1.0 ms 刚性地面 0.2 4.21 2.14 650 1 380 沙土地面 0.2 4.12 2.06 610 1 340 刚性地面 1.0 7.95 3.23 20 90 沙土地面 1.0 8.85 3.27 20 90 表 7 不同倾斜角度下柱形装药在前、后方的马赫杆高度
Table 7. Mach stem heights of cylindrical chargeswith different inclined angles at the front and rear
角度/(°) 马赫杆高度/mm 1.0 ms 2.0 ms 3.0 ms 前方 后方 前方 后方 前方 后方 0 90 90 480 480 1 010 1 010 30 310 170 740 1 040 1 180 1 850 60 130 390 750 1 180 1 260 1 900 90 130 130 610 610 1 250 1 250 表 8 不同长径比柱形装药的马赫杆高度
Table 8. Mach stem height of cylindrical charge with different aspect ratio
长径比 马赫杆高度/mm 1.0 ms 2.0 ms 3.0 ms 前方 后方 侧方 前方 后方 侧方 前方 后方 侧方 1∶1 480 150 90 1 250 870 520 1 970 1 520 1 060 1.5∶1 380 170 60 940 990 410 1 480 1 740 890 2∶1 310 170 50 740 1 050 350 1 180 1 850 790 2.5∶1 260 170 40 610 1 080 310 960 1 920 740 -
[1] 陈龙明, 李志斌, 陈荣. 装药动爆冲击波特性研究 [J]. 爆炸与冲击, 2020, 40(1): 013201. DOI: 10.11883/bzycj-2019-0029.CHEN L M, LI Z B, CHEN R. Characteristics of dynamic explosive shock wave of moving charge [J]. Explosion and Shock Waves, 2020, 40(1): 013201. DOI: 10.11883/bzycj-2019-0029. [2] 聂源, 蒋建伟, 李梅. 球形装药动态爆炸冲击波超压场计算模型 [J]. 爆炸与冲击, 2017, 37(5): 951–956. DOI: 10.11883/1001-1455(2017)05-0951-06.NIE Y, JIANG J W, LI M. Overpressure calculation model of sphere charge blastingwith moving velocity [J]. Explosion and Shock Waves, 2017, 37(5): 951–956. DOI: 10.11883/1001-1455(2017)05-0951-06. [3] 蒋海燕, 李芝绒, 张玉磊, 等. 运动装药空中爆炸冲击波特性研究 [J]. 高压物理学报, 2017, 31(3): 286–294. DOI: 10.11858/gywlxb.2017.03.010.JIANG H Y, LI Z R, ZHANG Y L, et al. Characteristics of air blast wave field for explosive chargemoving at different velocities [J]. Chinese Journal of High Pressure Physics, 2017, 31(3): 286–294. DOI: 10.11858/gywlxb.2017.03.010. [4] 郗洪柱, 孔德仁, 乐贵高, 等. 圆柱装药近地空爆冲击波峰值超压空间转换模型 [J]. 高压物理学报, 2021, 35(3): 032301. DOI: 10.11858/gywlxb.20200652.XI H Z, KONG D R, YUE G G, et al. Space conversion model of peak overpressure in near-earth air blast shockwave with cylindrical charge [J]. Chinese Journal of High Pressure Physics, 2021, 35(3): 032301. DOI: 10.11858/gywlxb.20200652. [5] 郑监, 卢芳云, 陈荣. 柱形装药条件下锥形水中爆炸激波管内的冲击波特性 [J]. 爆炸与冲击, 2021, 41(10): 103201. DOI: 10.11883/bzycj-2020-0316.ZHENG J, LU F Y, CHEN R. Shock wave characteristics in a conical water explosion shock tube under cylindrical charge condition [J]. Explosion and Shock Waves, 2021, 41(10): 103201. DOI: 10.11883/bzycj-2020-0316. [6] 廖真, 唐德高, 李治中, 等. 近地面空中爆炸马赫反射数值模拟研究 [J]. 振动与冲击, 2020, 39(5): 164–169, 176. DOI: 10.13465/j.cnki.jvs.2020.05.022.LIAO Z, TANG D G, LI Z Z, et al. Numerical simulation for Mach reflection in air explosion near ground [J]. Journal of Vibration and Shock, 2020, 39(5): 164–169, 176. DOI: 10.13465/j.cnki.jvs.2020.05.022. [7] 陈鑫, 高轩能. 炸药近地爆炸的数值模拟及影响参数的分析 [J]. 华侨大学学报(自然科学版), 2014, 35(5): 570–575. DOI: 10.11830/ISSN.1000-5013.2014.05.0570.CHEN X, GAO X N. Numerical simulation and analysis of influence parameters for explosions near ground [J]. Journal of Huaqiao University (Natural Science), 2014, 35(5): 570–575. DOI: 10.11830/ISSN.1000-5013.2014.05.0570. [8] WANG Y, WANG H, CUI C Y, et al. Investigating different grounds effects on shock wave propagation resulting from near-ground explosion [J]. Applied Sciences, 2019, 9(17): 3639. DOI: 10.3390/app9173639. [9] LU Y, WANG Z Q. Characterization of structural effects from above-ground explosion using coupled numerical simulation [J]. Computers & Structures, 2006, 84(28): 1729–1742. DOI: 10.1016/j.compstruc.2006.05.002. [10] DENG G Q, YU X. Numerical study on the case effect of a bomb air explosion [J]. Defence Technology, 2021, 17(4): 1461–1470. DOI: 10.1016/j.dt.2020.08.003. [11] LI T, WANG C, YANG T H, et al. A novel construction method of computational domains on large-scale near-ground explosion problems [J]. Journal of Computational Physics, 2020, 407: 109226. DOI: 10.1016/j.jcp.2019.109226. [12] YOU W B, DING Y H, YAO Y. Static explosion field reconstruction based on the improved biharmonic spline interpolation [J]. IEEE Sensors Journal, 2020, 20(13): 7235–7240. DOI: 10.1109/JSEN.2020.2978502. [13] 任朋飞, 王显会, 张进成, 等. 刚性地面爆炸冲击场仿真与试验研究 [J]. 科学技术与工程, 2017, 17(21): 30–36. DOI: 10.3969/j.issn.1671-1815.2017.21.006.REN P F, WANG X H, ZHANG J C, et al. Simulation and experimental study on explosive shock field on rigid ground [J]. Science Technology and Engineering, 2017, 17(21): 30–36. DOI: 10.3969/j.issn.1671-1815.2017.21.006. [14] 赵蓓蕾, 崔村燕, 陈景鹏, 等. 近地爆炸地面冲击波传播规律的数值研究 [J]. 四川兵工学报, 2015, 36(9): 45–48. DOI: 10.11809/scbgxb2015.09.012.ZHAO B L, CUI C Y, CHEN J P, et al. Numerical study of propagation law of ground shock wave in near surface explosion [J]. Journal of Ordnance Equipment Engineering, 2015, 36(9): 45–48. DOI: 10.11809/scbgxb2015.09.012. [15] 成凤生, 宋浦, 顾晓辉, 等. TNT装药爆炸波在刚性平面上方传播反射的数值研究 [J]. 爆破器材, 2011, 40(4): 1–4. DOI: 10.3969/j.issn.1001-8352.2011.04.001.CHENG F S, SONG P, GU X H, et al. Numerical investigation into the propagation and reflection of TNT blast wave above rigid plane [J]. Explosive Materials, 2011, 40(4): 1–4. DOI: 10.3969/j.issn.1001-8352.2011.04.001. [16] 段晓瑜, 崔庆忠, 郭学永, 等. 炸药在空气中爆炸冲击波的地面反射超压实验研究 [J]. 兵工学报, 2016, 37(12): 2277–2283. DOI: 10.3969/j.issn.1000-1093.2016.12.013.DUAN X Y, CUI Q Z, GUO X Y, et al. Experimental investigation of ground reflected overpressure of shock wave in air blast [J]. Acta Armamentarii, 2016, 37(12): 2277–2283. DOI: 10.3969/j.issn.1000-1093.2016.12.013. [17] 杜红棉, 曹学友, 何志文, 等. 近地爆炸空中和地面冲击波特性分析和验证 [J]. 弹箭与制导学报, 2014, 34(4): 65–68. DOI: 10.3969/j.issn.1673-9728.2014.04.018.DU H M, CAO X Y, HE Z W, et al. Analysisand validation for characteristics of air and ground shock wave near field explosion [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2014, 34(4): 65–68. DOI: 10.3969/j.issn.1673-9728.2014.04.018. [18] 姚成宝, 王宏亮, 浦锡锋, 等. 空中强爆炸冲击波地面反射规律数值模拟研究 [J]. 爆炸与冲击, 2019, 39(11): 112201. DOI: 10.11883/bzycj-2018-0287.YAO C B, WANG H L, PU X F, et al. Numerical simulation of intense blast wave reflected on rigid ground [J]. Explosion and Shock Waves, 2019, 39(11): 112201. DOI: 10.11883/bzycj-2018-0287. [19] 吴赛, 赵均海, 张冬芳, 等. 刚性地面对爆炸冲击波传播规律的影响 [J]. 工程爆破, 2020, 26(2): 1–8. DOI: 10.3969/j.issn.1006-7051.2020.02.001.WU S, ZHAO J H, ZHANG D F, et al. Influence of rigid ground on propagation characteristics of blast wave [J]. Engineering Blasting, 2020, 26(2): 1–8. DOI: 10.3969/j.issn.1006-7051.2020.02.001. [20] 辛春亮, 徐更光, 刘科种, 等. 考虑后燃烧效应的TNT空气中爆炸的数值模拟 [J]. 含能材料, 2008, 16(2): 160–163. DOI: 10.3969/j.issn.1006-9941.2008.02.011.XIN C L, XU G G, LIU K Z, et al. Numerical simulation of TNT explosion with post-detonation burning effect in air [J]. Chinese Journal of Energetic Materials, 2008, 16(2): 160–163. DOI: 10.3969/j.issn.1006-9941.2008.02.011. [21] 郭炜, 俞统昌, 金朋刚. 三波点的测量与实验技术研究 [J]. 火炸药学报, 2007, 30(4): 55–57, 61. DOI: 10.3969/j.issn.1007-7812.2007.04.014.GUO W, YU T C, JIN P G. Test of triple point and study on its test technology [J]. Chinese Journal of Explosives and Propellants, 2007, 30(4): 55–57, 61. DOI: 10.3969/j.issn.1007-7812.2007.04.014. [22] 畅博, 谷鸿平, 牛晨伟, 等. 运动炸药近地爆炸冲击波场特性研究 [J]. 爆破, 2018, 35(3): 49–54. DOI: 10.3963/j.issn.1001-487X.2018.03.008.CHANG B, GU H P, NIU C W, et al. Characteristics investigation on shock wave field of near ground blasting of moving explosive charge [J]. Blasting, 2018, 35(3): 49–54. DOI: 10.3963/j.issn.1001-487X.2018.03.008. [23] 曾一鑫, 白春华, 陈健, 等. 含铝炸药震源激发地震波的数值模拟 [J]. 高压物理学报, 2013, 27(6): 872–876. DOI: 10.11858/gywlxb.2013.06.012.ZENG Y X, BAI C H, CHEN J, et al. Simulations of seismic waves stimulated by aluminized explosive [J]. Chinese Journal of High Pressure Physics, 2013, 27(6): 872–876. DOI: 10.11858/gywlxb.2013.06.012. -