D型战斗部破片飞散性及威力场快速计算*

2019-06-05李翔宇李振铎梁民族

爆炸与冲击 2019年4期

李翔宇，李振铎，梁民族

（国防科技大学文理学院，湖南长沙 410073）

定向战斗部[1-3]一般通过特殊的结构设计和起爆控制，在破片飞散之前运用一些机构适时调整破片攻击方向，使破片在一定范围内相对集中，提高战斗部对目标的毁伤效率。目前研究较多的定向战斗部包括偏心起爆[4-5]、可变形、动能杆、破片芯和可控旋转战斗部等，其中可变形战斗部具有良好的军事应用前景。可变形战斗部是在主装药起爆前通过起爆目标方向的辅助装药使壳体变形，从原来的对称结构变成非对称的D型结构，将战斗部的环向破片更多地集中在目标方向，提高目标方向的破片密度，实现对目标的高效毁伤。因此，D型结构是可变形战斗部弹体变形型面的重要参考，研究D型战斗部的破片飞散特性对于可变形战斗部设计具有重要的指导意义。龚柏林等[6]、王马法等[7-8]、李振铎等[9]采用实验和数值模拟方法初步研究了D型结构的破片飞散特性，获得了D型结构破片能量分布特性规律。

目前，国内外多家单位采用射击线[10-12]方法建立了典型目标的毁伤评估模型及程序。采用射击线描述破片飞散区参数的方法主要有两种，一是利用经验公式进行估算，二是利用有限元软件对破片飞散区的形成过程进行数值模拟。对于传统周向均匀战斗部，经验公式能够很好地估算各种装填结构下的破片飞散分布。然而对于D型定向战斗部，能够合理描述破片飞散分布的经验公式较少。因此，本文中分别对型面宽度为90°、120° 和150° 三种D型战斗部的破片飞散规律进行实验和数值模拟研究，考察型面宽度和起爆模式对破片威力场的影响规律，拟合出快速计算D型结构破片威力场的计算公式。

1 实验

1.1 实验装置及弹靶布置

图1(a)为D型战斗部实验装置示意图，包括壳体、装药和端盖。壳体包括形成自然破片的圆弧段和形成预制破片的直线段，壳体厚度0.48 cm、高度8 cm、内径10 cm，材料为45钢；直线段壳体采用横向和纵向刻槽控制破片形状，刻槽深度为0.38 cm，相邻刻槽间距为0.48 cm；圆弧段对应的半径为r，直线段两端与圆心夹角定义为方位角α，直线段中间对应方位角为0°，图1(b)和(c)圆心角分别为120° 和90°，记为D-120° 和D-90°；直线段和圆弧段壳体边界采用焊接连接。端盖为厚度1.0 cm的LY12硬铝，用螺钉与壳体铆接。装药形状与壳体一致，外径10 cm、高度8 cm，材料为TNT。

弹靶相对位置布置如图2所示，其中弹靶间距为350 cm，直线段正方向布置3块宽100 cm、高250 cm靶板，用网格线将靶板分成多个边长25 cm的正方形，每块靶板对应方位角为8.18°。靶板从左向右依次编号为1#、2#、3#，通过靶板记录破片的空间分布。靶板上布置测速靶网，采用通断靶网和高速摄影记录破片着靶信号，进而获得破片的空间分布和速度分布。

1.2 实验结果及分析

通过靶板上破片着靶位置及数量统计，得到破片密度分布和破片数占比随方位角变化结果，如图3所示。可以看出，三种结构在8.18°方位角范围内的破片密度基本持平，而后破片密度随方位角迅速下降。统计表明，D-90°、D-120°、D-150° 三种结构在24.54° 方位角内破片数分别为破片总数的95.89%、91.92%、83.23%。利用Boltzmann分布对破片比例随方位角变化曲线进行拟合，得到D-90°、D-120° 和D-150°结构中包含 90% 破片的方位角范围分别为 21.16°、23.88° 和 30.08°。

图3 破片密度分布和破片数比例随方位角变化图Fig.3 Variation of fragment density distribution and fragment ratio with azimuthal angle

图4(a)给出了每块靶板上首次着靶时刻的高速摄影照片，图4(b)给出了基于高速摄影和靶网得到的破片速度。对比发现，D-90° 结构在0° 方位角处的破片速度最大，D-150° 结构的破片速度最小，且破片速度近似呈高斯正态分布。主要原因是型面宽度较小的结构，爆轰波驱动破片的有效装药多，因而破片获得较高速度；随着方位角的增加，爆轰波从垂直入射逐渐变成斜入射甚至发生马赫反射，致使爆轰波对破片驱动加速能力降低，造成不同型面宽度下相同方位角处破片速度的差距逐步缩小。

图4 破片着靶高速摄影照片及破片速度结果Fig.4 High-speed photography of the target and the experimental results of fragment velocity

2 数值模拟

2.1 计算模型及材料参数

针对实验装置建立相同结构的有限元模型，如图5(a)所示。采用LS-DYNA软件中流固耦合算法计算炸药爆轰驱动破片飞散过程，空气和炸药为Euler网格，端盖、壳体和破片为Lagrange网格。采用高能炸药材料和JWL状态方程描述炸药爆轰过程。采用线性多项式状态方程来描述空气介质特性，并施加压力外流边界条件。采用Grüneisen状态方程和流体弹塑性材料模型描述壳体和端盖的动态响应行为。采用随动硬化材料模型描述破片动态响应。图5(b)为直线段预制破片在某个时刻的空间分布图，可以看出破片变形很小，基本保持形状。

图5 D型战斗部有限元模型Fig.5 Finite element model of D-shaped warhead

2.2 与实验结果对比

图6(a)给出了D-90° 模型在350 cm处靶板上的破片分布的数值模拟结果，图中每一点代表一枚破片，纵坐标为破片高度位置，横坐标为破片位置对应的方位角，破片着靶速度在图上以不同颜色区分。可以看出，方位角0° 位置处的破片比较集中，速度最大；上下和左右边缘处的破片密度和速度都较小。图6(b)～(d) 分别给出了三种D型战斗部破片飞散的实验与数值模拟结果，结果表明D-90°、D-120° 和D-150° 三种结构中包含90%破片的方位角分别为18.4°、25.2° 和30°，与实验结果吻合较好。

图6 三种D型战斗部破片分布结果Fig.6 Fragment distributions of three D-shaped warheads

图7 给出了 D-90°、D-120°和 D-150°三种 D 型战斗部破片速度的实验和数值模拟结果，破片速度呈现中间高、两侧低的趋势，破片速度分布规律基本一致，近似符合高斯正态分布。

2.3 破片飞散规律及分析

为了考察起爆方式对破片分布的影响，对D-120° 结构在几何中心起爆1、上端面中心起爆2、双端面中心起爆3、中心线起爆4、偏心点起爆5、上端面偏心起爆6、双端面偏心起爆7和偏心线起爆8等8种起爆方式下破片飞散特性进行分析，八种起爆模式如图8(a)所示。以周向均匀战斗部中心起爆的能量分布为参考基准，得到破片能量密度随方位角变化的曲线，如图8(b)所示。可以看出，D型战斗部破片在目标方向上能量聚焦明显，通过选择合理的起爆模式能够使杀伤能量在目标方向上相对集中。偏心起爆下目标方向上的破片至起爆点的距离大，对应有效装药量多，作用在破片上的抛射能量相对较大，所以偏心起爆下能量聚焦效果更佳。偏心起爆的四种模式中，偏心线起爆和双端面偏心起爆在目标方向上的能量聚焦能力强且总能量输出较高，20° 方位角内破片总能量分别是圆柱形中心起爆能量的3.4倍和3.3倍。

图7 三种D型战斗部破片速度的实验和数值模拟结果Fig.7 Fragment velocities of two D-shaped warheads

图8 不同起爆方式下破片能量分布Fig.8 Fragment energy distribution in different denotation modes

图9 不同型面宽度下破片能量密度和破片总能量随方位角变化Fig.9 Variation of energy density and total energy of fragments with azimuthal angle at different profile widths

为了考察型面宽度对破片飞散区域的影响，对60°、90°、120° 和150° 的四种型面宽度的战斗部在端面偏心起爆下破片飞散过程进行研究，得到破片能量密度和破片总能量随方位角变化结果，如图9所示。结果表明，D-60° 结构由于接近于圆柱形，其破片密度变化相对不明显，能量分布规律趋向于周向均匀战斗部。型面宽度在90°～120° 之间时目标方向能量增益较为明显，在0° 方位角附近型面宽度越小破片能量密度越大，在12° 方位角附近破片能量密度基本持平并开始以较快的速度衰减，40°～60° 方位角范围内破片的能量密度存在一个平台段。D-150° 型面宽度较大，但在350 cm的弹靶距离下破片总能量和能量聚焦能力并不能达到最佳；D-90° 和D-120° 两种结构在42° 方位角内的破片总能量与能量聚焦能力接近且为几种结构中最高。相对于周向均匀战斗部，D-90° 战斗部在43.2° 位角范围内的破片总能量增益达 241.7%；D-120° 战斗部在 47.6° 方位角范围内的破片能量增益为232.6%。

3 破片飞散区分布快速计算

为了快速得到D型战斗部在不同型面宽度和起爆模式下破片飞散区域，建立如图10所示的分析模型，图中装药半径为7 cm；O1、O2、O3为起爆点其中O1位于装药最左端，O3位于装药圆心，O2位于O1和O3之间；直线段破片编号从下向上依次为1， 2，…，11，其中A处对应破片编号 9；θ2、θ3分别为A处爆轰波传播方向和破片飞散方向与对称轴O1O3的夹角，爆距d定义为起爆点到破片所在平面的垂直距离。

图10 D型战斗部破片飞散区快速计算分析模型Fig.10 Analysis model of rapid calculation of fragment force field of D-shaped fragmentation warhead

对O1、O2、O3三种起爆模式下的破片飞散特性进行数值模拟，由于对称位置的破片在起爆后沿水平方向飞散，边界处破片受侧向稀疏波作用飞散方向随机性较大，仅分析图中编号1～11的破片。得到各破片初始位置θ2值以及对应角度θ3的数值模拟结果如表1所示。

表1 不同起爆模式下破片飞散特性Table 1 Dispersion characteristic of fragment in different detonation modes

为了快速预测不同爆距d下的破片飞散分布情况，将爆距d与装药半径r的比值与表1中的sinθ2、sinθ3进行二次曲面拟合，拟合得到D-120°结构的破片分布快速计算公式：

式中：θ3[120°]和θ2[120°]分别为D-120°结构的θ3和θ2。

为了验证拟合公式的准确性，将D-120° 结构装药半径调整为5 cm，分析在爆距6 cm和7.5 cm下的破片飞散结果。图11给出了两种爆距情况下数值模拟结果与公式(1)计算结果的对比图。结果表明，拟合的破片飞散区域公式与数值模拟结果吻合较好。

图11 不同爆距下破片分布结果Fig.11 Fragment distributions at different burst distances

对D-90° 和D-150° 结构采用同样的方式进行拟合，得到破片分布快速计算公式：

对于其他型面宽度D型战斗部，可以根据已有的三个拟合关系式，将对应工况的破片飞散状态进行计算，再通过构造二次插值函数来获得所求结构的破片空间分布数据，其他型面宽度对应的破片分布快速计算公式为：

式中：x为D型面宽度对应的圆心角。

针对型面宽度105° 的D型战斗部，在装药半径5 cm，爆距7.5 cm情况下，将x=105° 代入公式（4），即可快速计算型面宽度105° 下破片分布结果，两种方法得到的sinθ3结果如图12所示。利用构造二次插值函数得到的破片飞散角度与数值模拟结果吻合较好，因此通过构造二次插值函数来计算其他型面的破片分布是可行的。采用此方法可为不同型面宽度的D型战斗部破片分布快速计算提供一种途径。