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]研究发现，球形装药的超压、比冲量、正压作用时间以90°为分界线增强或减弱。聂源等[2]得到了球形装药的动态爆炸冲击波超压场计算模型表达式。蒋海燕等[3]研究了球形装药动爆速度影响冲击波场空间和位置参数的规律，其中关于柱形装药的动爆特征研究较少，且普遍为二维平面的数值模拟研究。
在柱形装药的爆炸冲击波研究中，郗洪柱等[4]基于镜像法研究了柱装药近地空爆冲击波峰值超压空间转换模型，并进行了实验验证；郑监等[5]对不同圆锥角和不同柱形装药质量下锥形激波管内的冲击波传播过程进行了模拟。在近地爆炸研究中，有较多关于近地爆炸冲击波的反射理论，廖真等[6]对不同装药类型和形状产生的三波点轨迹进行了分析，并结合实验得到了不同装药峰值压力的差别；陈鑫等[7]研究了地面介质对近地马赫波超压的影响，并与刚性地面进行了对比。Wang等[8]基于ANSYS对近地爆炸进行了数值模拟，分析了不同地面对冲击波的反射作用；Lu等[9]研究了近地爆炸对地面建筑物结构响应的过程分析，并分析了土壤中冲击波的作用；Deng等[10]通过研究带壳装药的近地爆炸，得到了马赫波的传播与温度场变化的对比；Li等[11]通过AME（adaptive mesh enlargement）方法对网格进行了自适应调整，从而可以对大规模流体场进行计算；You等[12]研究了一种改进的有限测点爆炸超压数据插值重建方法，从而可以得到地面附近爆炸的峰值压力分布；任朋飞等[13]对炸药在刚性地面爆炸进行了实验与数值模拟研究，验证了ALE（arbitrary Lagrange-Euler）方法算法的正确；赵蓓蕾等[14]利用LS-DYNA分析了地面冲击波传播规律，得到了马赫杆、地面超压的衰减规律；成凤生等[15]比较了柱形与球形装药超压变化，研究了波阵面的传播规律；段晓瑜等[16]开展了不同材质炸药的近地反射超压的实验研究；杜红棉等[17]拟合出了近地爆炸冲击波马赫杆的三维三波点马赫波轨迹，验证了数值模拟结果的可靠性；姚成宝等[18]使用网格自适应技术，对炸药的空中强爆炸在不同爆炸高度下的冲击波地面反射过程进行了数值模拟；吴赛等[19]得出刚性地面的反射作用会显著增强爆炸效应，可对结构或构件造成更大的危害，提高了对工程中的反射效应最低炸高理论值。
综合之前的研究成果可以看出，目前关于柱形装药的冲击波传播理论研究较少，且没有考虑实际情况下柱形装药战斗部的动态参数特征。因此，本文中，基于AUTODYN-3D软件对不同末弹道参数的柱形装药近地爆炸进行数值模拟，研究动爆条件下的柱形装药各参数对近地冲击波传播的规律，并对峰值压力、反射波产生、马赫波等冲击波参数和传播过程进行分析。
1. 数值模拟方案设计
1.1 数值模拟建模
以常规榴弹战斗部装药量作为参考（依据口径不同，装药量一般在1.75～5.00 kg），装药药量选取4.00 kg，柱形装药近地爆炸示意图如图1所示。炸药选择TNT炸药，起爆方式为中心点起爆，材料参数为AUTODYN材料库中的默认参数，TNT与空气的材料模型均采用Euler网格方法计算，TNT的状态方程参数见表1，表中：A、B、R1、R2和ω为状态方程独立参数，E0、V分别为炸药比内能和比容。为避免网格尺寸影响计算精度，将0.859 kg的装药在网格尺寸分别为0.5、1.0、1.5、2.0 cm进行空爆数值模拟。将数值模拟结果与Henrych和Brode空爆公式[20]进行对比，随着网格精度的提高，近地压力值与Henrych公式计算结果逐渐逼近。结合廖真等[6]在近地爆炸数值模拟中的网格尺寸分析，验证了1.0 cm的网格尺寸与郭炜等[21]实验中的三波点高度具有较好的一致性。同时在三维模型中要兼顾空气域尺寸和计算速度，故选取空气域网格尺寸为1.0 cm。
图 1 柱形装药近地动爆示意图Figure 1. Schematic diagram of near ground dynamic explosion of cylindrical charge表 1 TNT炸药状态方程参数Table 1. Parameters of equation of state for TNT explosiveA/GPa B/GPa R1 R2 E0/GPa ω V 3.7377 3.747 4.15 0.9 0.06 0.3 1 建模分为2个方向：第1个方向为柱形装药前后方冲击波传播方向，以x轴为基准方向，在x轴的距原点5 m处放置柱形炸药；第2个方向为炸药侧方，以y轴为基准方向，在原点位置放置炸药。2个方向的建模网格数量分别为2 000 000和3 600 000。x方向与y方向的空气域宽度均可将炸药完全覆盖，前后方向建模尺寸分别为10.00、0.10和2.00 m；侧方向建模尺寸分别为5.00、0.36和2.00 m。空气的材料参数如表2所示，材料模型见表3。表中：ρ为空气密度，γ为绝热指数，e为空气的初始比内能。
表 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 无 无 TNT炸药和空气分别采用JWL状态方程和理想气体状态方程。
p=A(1−ωR1V)e−R1V+B(1−ωR2V)e−R2V+ωE0V (1) 式中：p为爆轰产物压力。
在空气域前后面、侧面和上面设置物质流出边界，剩余两面为自由反射边界，模拟对称方向和刚性地面。设置以装药为中心的辐射型测点不能满足不同姿态近地冲击波的传播特征，所以在相同高度上设置等距的测点，并沿z方向间隔相同高度设置，测点设置如图2所示。
图 2 柱形装药近地爆炸建模及测点布置Figure 2. Modeling and measuring point diagram of near ground explosion of cylindrical charge1.2 柱形装药在刚性地面近地爆炸实验对比
将数值计算模型与任朋飞等[13]所做的刚性地面爆炸冲击场实验进行对比，在混凝土路面上放置90 mm厚的钢板来模拟刚性地面，并对长径比为1∶3的0.2、0.8和2.0 kg等3种装药量的柱形装药进行了实验。在距爆心位置处的1.12和1.95 m设置了测点，建模如图3所示。测点压力的计算值与实验值的对比如表4所示。
表 4 测点压力的计算值与实验值对比Table 4. Comparison between calculated values and experimental values on measuring point pressureTNT质量/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 实际情况下的刚性地面反射强度弱于数值模拟中的绝对反射，装药质量小导致网格精度变差，故数值模拟值与实际值存在误差。可以看出，随着装药质量的增加，误差不断减小，2.0 kg装药的平均误差为9.43%，所以对4.0 kg装药的数值模拟符合实际的刚性地面反射情况。
1.3 弹前激波和地面沙土对近地爆炸的影响
1.3.1 弹前激波对近地爆炸的影响
飞行状态下的柱形装药在前端面产生压缩波，随着运动速度的不断提高，在弹前产生的波逐渐从弱的压缩波变为强的压缩波，该冲击波是由于弹前压缩波的聚集、反射与叠加[22]形成的。在数值模拟中，先让弹丸在空气中飞行1.0 ms以上，使得装药在空气中产生较稳定的激波及压力场，然后对其进行起爆。在数值模拟中对比了静爆、动爆、静爆加激波和动爆加激波4种起爆类型，图4～5给出了长径比为2∶1的柱形装药在600和1 200 m/s速度下的稳定激波压力云图和矢量图。
图 4 不同飞行速度下柱形装药的压力云图Figure 4. Pressure nephogram of cylindrical charge at different flight speeds
图 5 不同飞行速度下柱形装药的压力矢量图Figure 5. Pressure vector diagrams of cylindrical charge at different flight speeds柱形装药在距离地面1.2 m处起爆，并在刚性界面反射传播形成马赫波，观测上述4种类型产生的爆炸压力场及峰值压力，4.0 ms时1 200 m/s速度下4种起爆类型产生的马赫波如图6所示。
图 6 不同起爆条件下柱形装药的压力云图Figure 6. Pressure nephograms of cylindrical charge under different initiation conditions由于爆炸产生的入射压力较大，是激波最大压力的200倍以上，在与地面接触产生的反射波是其压力的20倍，所以激波对马赫波的传播过程、马赫波的峰值压力影响有限，仅对起爆点上方入射波的传播有一定干扰。
从表5中可以看出，激波和动爆的联合作用与动爆的峰值压力差别很小，平均差0.12%，并且对马赫杆高度的主要影响为装药动爆速度的施加，因此在数值模拟中忽略弹前激波的影响。
表 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 1.3.2 地面沙土对近地爆炸的影响
实际战场环境地面情况多样，沙土等软质地面不能视作刚体。AUTODYN中沙土选择Linear材料模型和von Mises强度模型[23]，密度为1.6 g/cm3。计算中沙土采用Lagrange模型，数值模拟采用流固耦合计算方法，在0.2、0.4、0.6、0.8和1.0 m等5种近地高度下与刚性地面进行对比。
刚性和沙土地面的对比如图7所示，可以看出沙土介质对马赫波的影响有限。
图 7 不同地面条件下柱形装药的压力云图Figure 7. Pressure nephograms of cylindrical charge under different ground conditions从表6中看出，沙土对近地爆炸冲击波的峰值压力影响较为有限，0.5 ms时，峰值压力平均差1.76%，马赫杆高度平均差6.48%；1 ms时，峰值压力平均差1.17%，马赫杆高度平均差4.79%，说明数值模拟中沙土对近地冲击波的影响较小。
表 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 2. 近地动爆的冲击波规律探究
2.1 爆心高度对冲击波传播的影响
近地爆炸时冲击波到达地面并反射形成反射波，反射波由于空气的压缩加热其速度大于正常产生的入射波，将逐渐追上入射波；之后反射波与入射波叠加形成马赫波；马赫波在水平方向传播的同时其高度不断扩大，反射波及马赫波形成的示意如图8所示。
图 8 近地爆炸马赫波的产生过程示意图Figure 8. Schematic diagram of Mach wave generation process of near ground explosion图9为炸高1.0 m，装药质量为4 kg，长径比为2 : 1的柱形装药在倾斜角为0°的马赫波产生过程。
图 9 不同时刻下柱形装药的马赫波压力云图Figure 9. Mach wave pressure nephogram of cylindrical charge at different times以0.2 m为间隔设置了0.2～1.8 m范围内的炸高，柱形装药长径比为2∶1的近地爆炸数值模拟对照。起爆高度会直接影响冲击波到达地面的时间，导致入射波到达地面的峰值压力随起爆高度增加而减少。国防工程设计草案[19]中，是否考虑刚性地面的反射增强效应的表达式为：
HW1/3＞0.35m/kg1/3 (2) 式中：H为爆心距地面的距离，m；W为TNT炸药当量，kg。
按照式(2)计算得到高度为0.56 m，实际上反射效应在此高度仍然较大，因此对柱形装药在不同起爆高度下的近地峰值压力进行对比，如图10所示。从图10可以看出，峰值压力随爆心距和起爆高度整体呈现下降趋势，在2.0 m爆心距后的峰值压力差距明显减小。当爆心距为0.5 m时，0.4 m爆高的峰值压力是1.8 m爆高的3.6倍。峰值压力的不同造成了马赫杆的高度差距，2.0 ms时0.4 m爆高的马赫杆高度与1.6 m相比相差1.93 m。图11为爆高为1.0 m时水平方向上不同高度测点柱形装药的峰值压力-时间曲线，从0 m高度峰值压力曲线中可以看到，3.5 m测点处存在明显的2个波峰，并向后传递；在1.0 m高度的峰值压力有一定程度的衰减，而在2.0 m高度处测点的峰值压力几乎与自由空气中爆炸相近。
图 10 不同爆心高度下柱形装药的峰值压力的对比Figure 10. Comparison of peak pressures of the cylindrical charge under different heights of the explosion center
图 11 不同高度测点柱形装药的压力曲线Figure 11. Pressure curves of cylindrical charge at measuring points of different heights马赫波的产生与装药的比例炸高相关，依据极限入射角与比例炸高倒数的关系曲线来确定马赫波反射区的范围。当装药质量大并且炸心高度较低时，比例炸高倒数较大。从图12中看出，当比例炸高倒数大于0.8时，入射角将趋近于40°，所以马赫波在水平方向的发生位置与炸高近似。
图 12 入射角与比例炸高的关系Figure 12. Relationship between incident angle and proportional explosion height2.2 倾斜角度对冲击波传播的影响
近地面冲击波斜反射的传播情况较复杂，将入射波质点速度u1分解为平行于地面的分量和垂直于地面的分量，垂直于地面的分量形成的反射波在被压缩和加热的空气中传播从而速度增大。反射波在入射角θ1较小时，水平方向的速度分量较小，随着θ1的增加，反射波在水平方向的速度分量增大。从某个θ1开始，反射波质点速度u2的垂直分量将大于u1的垂直分量，2个波阵面的交点A将离开地面。斜反射过程的示意图如图13所示。
图 13 近地爆炸冲击波斜反射示意图Figure 13. Schematic diagram of oblique reflection of shock wave produced by near ground explosion柱形装药径向冲击波比轴向冲击波传播速度快，在装药倾斜角度不大时，轴向入射波将率先到达地面并进行反射。随着倾斜角度的增加，轴向入射波与地面接触的位置不断向前移动，使得斜反射与地面的夹角不断减小，水平方向的分量增加而反射波的峰值压力不断减弱。所以装药前侧的马赫杆高度在开始大于后侧的马赫杆，之后被逐渐追上并超越。以15°和75°静爆条件下马赫波传播过程为例可以看出这一趋势，如图14所示。
图 14 倾斜角度不同时柱形装药的马赫波压力云图Figure 14. Mach wave pressure nephograms of cylindrical charges with different inclined angles在大于45°附近的倾斜装药，不会再产生两个波同时抵达地面的情况，此时径向入射波将率先在后侧产生反射波。装药后侧产生的马赫杆高度将始终大于前侧的马赫杆高度，图15所示为2种不同的入射波抵近地面的状态。
图 15 倾斜角度不同时柱形装药的近地状态压力云图Figure 15. Near ground pressure nephograms of cylindrical charges with different inclined angles随着倾斜角度的增加，装药两侧的峰值压力较为接近，传播速度及马赫杆高度随着角度的增加而逐渐增大，在3.5 ms时，90°比0°的马赫杆高出77.1%。通过倾斜角相同时装药前方、后方和侧方3个方向的传播过程对比，可以看出在同一时刻冲击波的传播差别，图16所示为45°倾斜角、4.0 ms时刻3个方向的冲击波传播过程。
图 16 倾斜角为45°时不同方向柱形装药的压力云图Figure 16. Pressure nephograms of cylindrical charge in different directions when the inclination angle is 45°图17给出了4.0 m处随倾斜角度变化3个方向的峰值压力变化。在装药前方的峰值压力随着倾斜角度的增加不断增大，并在75°时达到最大值；而在装药的侧方和后方的峰值压力随着倾斜角度的增加整体呈现递减的趋势，分别在60°及90°取得最小值。
图 17 柱形装药的峰值压力随倾斜角度的变化Figure 17. Peak pressure of cylindrical charge varied with inclined angle表7为倾斜角度对马赫杆高度的影响，其中3.0 ms时前方马赫杆高度最大增长幅度为24.8%，后方马赫杆高度最大增长幅度为88.1%。
表 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 由于柱形装药形成的冲击波在径向的传播速度比轴向的快，随着倾斜角度的增大，冲击波在各方向呈现出不同的波形和峰值压力。图18为不同倾斜角度下柱形装药在3个方向的峰值压力变化。
图 18 不同方向柱形装药的峰值压力随倾斜角度的变化Figure 18. Variation of peak pressure of cylindrical charge with inclined angle in different directions2.3 动爆对近地冲击波传播的影响
在实际作战中战斗部一般为动态爆炸，使整个冲击波场向速度方向产生移动，影响冲击波流场的演化过程，造成装药沿速度方向前端压力升，后端压力降。以倾斜角为30°，爆心高度为1.0 m的2∶1长径比装药进行研究，柱形装药的初速度依次设置为0、300、600、900、1 200和1 500 m/s。
图19为0和300 m/s的冲击波场生成过程对比，动爆会使整个冲击波场向速度方向移动，增大近地侧入射波的峰值压力而减小远地侧的入射波压力。动爆会加速装药前方的轴向冲击波，驱使其率先抵达地面，使得冲击波能量向前侧聚集，从而前侧马赫波峰值压力及马赫杆高度追上后侧。
图 19 不同动爆速度下柱形装药的冲击波压力云图Figure 19. Shock wave pressure nephogram of cylindrical charge under different dynamic explosion velocities当装药特性、测点位置等条件都相同时，在压力上升区域，装药速度越大，超压峰值的上升幅度越大，在压力下降区域，装药速度越大，超压峰值的下降幅度也越大。图20为4.0 m处不同方向的峰值压力随速度变化关系，从图20可以看出，随着动爆速度的提升，前方峰值压力显著增大。比较3.0 ms时不同速度的马赫杆高度发现，前方马赫杆高度最大增长幅度为156.2%，后方马赫杆高度最大幅度为76.2%，侧方马赫杆高度最大幅度为388.2%。
图 20 柱形装药的峰值压力随动爆速度变化Figure 20. Diagram of peak pressure of cylindrical charge varying with dynamic explosion velocity由于动爆速度指向地面，侧方冲击波受此影响使得压力区域快速下移并与地面接触，速度越快区域与地面的夹角越大，使得速度快的侧向马赫波高度将迅速升高。300和1 500 m/s的侧方冲击波传播过程对比见图21。
图 21 不同动爆速度下柱形装药在侧方的马赫波压力云图Figure 21. Mach wave pressure nephogram of cylindrical charge on the side under different dynamic explosion velocities对比不同动爆速度在3个方向的峰值压力变化曲线，如图22所示。随着动爆速度的提升，在装药的前方及侧方的不同距离峰值压力都呈现增加的趋势，而装药后方趋势较不明显，且在距爆心2.0 m距离外的峰值压力无较大变化。
图 22 不同方向柱形装药的峰值压力随动爆速度的变化Figure 22. Peak pressure varying with dynamic explosion velocity for cylindrical charge in different directions2.4 装药长径比对近地爆炸冲击波的影响
装药长径比会对自由空气场中冲击波波形产生影响，由于柱形装药的入射波特点为径向波传播快、轴向波传播慢，随着高度的增加其波的差异会减弱。对1∶1、1.5∶1、2 ∶1、2.5∶1的4种长径比柱形装药在1.0 m高度、0 m/s、30°倾斜角条件下对装药前部、后部、侧部3个方向进行分析，图23为2.0 ms时刻不同长径比装药在前后方向的传播对比。
图 23 不同长径比柱形装药在前、后方的马赫波压力云图Figure 23. Mach wave pressure nephogram of cylindrical charges with different aspect ratios at the front and rear随着柱形装药长径比的增加，装药前侧的马赫杆高度逐渐降低，而装药后侧的马赫杆高度逐渐增加。前侧马赫杆高度差距在3.0 ms时可达1倍，后侧差距则不到0.5倍。图24为4.0 ms时刻装药侧方的冲击波传播过程，随着装药长径比的增加，侧方的马赫杆高度逐渐降低，3.0 ms时装药比例1∶1的马赫杆平均高度是2.5∶1高度的1.43倍。
图 24 不同长径比柱形装药在侧方的马赫波压力云图Figure 24. Mach wave pressure nephogram of cylindrical charges with different aspect ratios on the side图25所示为不同比例装药在前、后、侧3个方向的峰值超压随距爆心距离的变化关系，不同的比例装药在各个方向呈现出不同的变化趋势，在前侧的变化最显著。从表8不同长径比装药的马赫杆高度对比可以看出，随着柱形装药长径比的增加，在前方和侧方的马赫杆高度均不断下降，而后方的马赫杆高度不断增加。
图 25 不同方向柱形装药的峰值压力随长径比变化Figure 25. Peak pressure varying with length-to-diameter ratio for cylindrical charge in different directions表 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 3. 正交优化结果分析
对爆心高度、倾斜角度、动爆速度、装药长径比进行16组正交优化的数值模拟，对各组选取距爆心2.0、3.0、4.0 m处测点，对装药前后方向的峰值压力进行分析，并测量在柱形装药前方、后方马赫杆高度。装药前方、后方距爆心水平距离4.0 m处的近地峰值超压的各因素影响如图26~27所示。
图 26 正交优化下对装药前方近地峰值压力的各因素影响Figure 26. Influence of various factors of near ground peak pressure in front of the charge under orthogonal optimization
图 27 正交优化下对装药后方近地峰值压力的各因素影响Figure 27. Influence of various factors of near ground peak pressure at the rear of the charge under orthogonal optimization图26中动爆速度的峰值超压极差最大，平均极差为1.21 MPa。动爆速度越高，前侧的峰值压力越大，其次影响因素为装药长径比和爆心高度。图27中各因素对压力的影响趋势基本相同，其中倾斜角度的极差最大，对峰值压力的影响占据主导。其余因素较倾斜角度的变化幅度较小，动爆速度及装药长径比对后侧峰值压力基本呈现负影响的趋势，爆心高度对后侧峰值压力的影响不明显。
由于建模尺寸的影响，对1.0 ms时刻装药前方马赫杆高度、后方马赫杆高度进行极差分析，结果如图28~29所示。
图 28 正交优化下对装药前方马赫杆高度的各因素影响Figure 28. Influence of various factors on Mach stem height in front of the charge under orthogonal optimization
图 29 正交优化下对装药后方马赫杆高度的各因素影响Figure 29. Influence of various factors on Mach stem height at the rear of the charge under orthogonal optimization从装药前方和装药后方马赫杆高度的极差可以看出，爆心越低，则马赫杆越高。倾斜角度是影响马赫杆高度的次要因素，且在30°附近取得最大值。动爆速度对前后高度呈现相反的趋势，这点与前面的结论相同。
4. 结 论
通过分析柱形装药近地爆炸中的4个参数：爆心高度、倾斜角度、动爆速度和装药长径比，得到了各参数对近地爆炸中各指标的影响规律。分析前、后、侧3个方向的近地冲击波传播规律及峰值压力和马赫杆高度，对各个方向的压力变化和马赫波传播进行了量化对比。综合4个变量的正交优化结果，各因素对不同的近地爆炸指标有较强的规律性，在柱形装药的近地爆炸研究中得到以下主要结论。
(1) 柱形装药静爆时，随装药倾斜角的增大，前方的峰值压力先减小后增大，后方呈现下降，而侧方呈上升趋势。马赫杆高度在装药前后的比值先增大后减小，侧方的马赫杆逐渐增高。倾斜角度与装药长径比是影响静爆冲击波前后侧传播差异的主要因素，3 ms时对前、后、侧3个方向的马赫杆高度的最大提升幅度依次为29.7%、88.1%和96.4%。
(2) 柱形装药动爆的趋势较一致，即动爆速度越高，装药前方和侧方的峰值压力及马赫杆高度越大，但对后方的影响较小。3.0 ms时，前方马赫杆高度最大增长幅度为156.2%，后方马赫杆高度最大幅度为76.2%，侧方马赫杆高度最大提升为388.2%。
(3) 动爆速度对装药前方马赫波的增益最大，平均极差为1.21 MPa；倾斜角度对装药后方的马赫波的增益最大，平均极差为1.38 MPa；爆心高度则对前后侧的马赫杆高度影响最大，平均极差为1 051.25 mm。
-
图 1 柱形装药近地动爆示意图
Figure 1. Schematic diagram of near ground dynamic explosion of cylindrical charge
图 2 柱形装药近地爆炸建模及测点布置
Figure 2. Modeling and measuring point diagram of near ground explosion of cylindrical charge
图 4 不同飞行速度下柱形装药的压力云图
Figure 4. Pressure nephogram of cylindrical charge at different flight speeds
图 5 不同飞行速度下柱形装药的压力矢量图
Figure 5. Pressure vector diagrams of cylindrical charge at different flight speeds
图 6 不同起爆条件下柱形装药的压力云图
Figure 6. Pressure nephograms of cylindrical charge under different initiation conditions
图 7 不同地面条件下柱形装药的压力云图
Figure 7. Pressure nephograms of cylindrical charge under different ground conditions
图 8 近地爆炸马赫波的产生过程示意图
Figure 8. Schematic diagram of Mach wave generation process of near ground explosion
图 9 不同时刻下柱形装药的马赫波压力云图
Figure 9. Mach wave pressure nephogram of cylindrical charge at different times
图 10 不同爆心高度下柱形装药的峰值压力的对比
Figure 10. Comparison of peak pressures of the cylindrical charge under different heights of the explosion center
图 11 不同高度测点柱形装药的压力曲线
Figure 11. Pressure curves of cylindrical charge at measuring points of different heights
图 12 入射角与比例炸高的关系
Figure 12. Relationship between incident angle and proportional explosion height
图 13 近地爆炸冲击波斜反射示意图
Figure 13. Schematic diagram of oblique reflection of shock wave produced by near ground explosion
图 14 倾斜角度不同时柱形装药的马赫波压力云图
Figure 14. Mach wave pressure nephograms of cylindrical charges with different inclined angles
图 15 倾斜角度不同时柱形装药的近地状态压力云图
Figure 15. Near ground pressure nephograms of cylindrical charges with different inclined angles
图 16 倾斜角为45°时不同方向柱形装药的压力云图
Figure 16. Pressure nephograms of cylindrical charge in different directions when the inclination angle is 45°
图 17 柱形装药的峰值压力随倾斜角度的变化
Figure 17. Peak pressure of cylindrical charge varied with inclined angle
图 18 不同方向柱形装药的峰值压力随倾斜角度的变化
Figure 18. Variation of peak pressure of cylindrical charge with inclined angle in different directions
图 19 不同动爆速度下柱形装药的冲击波压力云图
Figure 19. Shock wave pressure nephogram of cylindrical charge under different dynamic explosion velocities
图 20 柱形装药的峰值压力随动爆速度变化
Figure 20. Diagram of peak pressure of cylindrical charge varying with dynamic explosion velocity
图 21 不同动爆速度下柱形装药在侧方的马赫波压力云图
Figure 21. Mach wave pressure nephogram of cylindrical charge on the side under different dynamic explosion velocities
图 22 不同方向柱形装药的峰值压力随动爆速度的变化
Figure 22. Peak pressure varying with dynamic explosion velocity for cylindrical charge in different directions
图 23 不同长径比柱形装药在前、后方的马赫波压力云图
Figure 23. Mach wave pressure nephogram of cylindrical charges with different aspect ratios at the front and rear
图 24 不同长径比柱形装药在侧方的马赫波压力云图
Figure 24. Mach wave pressure nephogram of cylindrical charges with different aspect ratios on the side
图 25 不同方向柱形装药的峰值压力随长径比变化
Figure 25. Peak pressure varying with length-to-diameter ratio for cylindrical charge in different directions
图 26 正交优化下对装药前方近地峰值压力的各因素影响
Figure 26. Influence of various factors of near ground peak pressure in front of the charge under orthogonal optimization
图 27 正交优化下对装药后方近地峰值压力的各因素影响
Figure 27. Influence of various factors of near ground peak pressure at the rear of the charge under orthogonal optimization
图 28 正交优化下对装药前方马赫杆高度的各因素影响
Figure 28. Influence of various factors on Mach stem height in front of the charge under orthogonal optimization
图 29 正交优化下对装药后方马赫杆高度的各因素影响
Figure 29. Influence of various factors on Mach stem height at the rear of the charge under orthogonal optimization
表 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 
下载: 导出CSV
表 2 空气状态方程参数
Table 2. Parameters of equation of state for air
ρ/(kg·m−3) γ e/(MJ·kg−1) 1.225 1.4 206.8
下载: 导出CSV
表 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
下载: 导出CSV
表 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
下载: 导出CSV
表 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
下载: 导出CSV
表 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
下载: 导出CSV
表 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
下载: 导出CSV
-
[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. 期刊类型引用(1)
1. 王明涛，程月华，吴昊. 柱形装药空中运动爆炸冲击波荷载计算模型. 爆炸与冲击. 2024(07): 107-125 . 
本站查看其他类型引用(2)
-
下载:
计量
- 文章访问数: 593
- HTML全文浏览量: 120
- PDF下载量: 276
- 被引次数: 3