基于靶板冲塞式破坏的战斗部破片优化设计

余丽山; 李彦彬; 翟夕阳; 金学科; 何小九

doi:10.15892/j.cnki.djzdxb.2019.05.034

弹箭与制导学报 >

2019 , Vol. 39 >Issue 5: 143 - 146

DOI: https://doi.org/10.15892/j.cnki.djzdxb.2019.05.034

基于靶板冲塞式破坏的战斗部破片优化设计

余丽山 ,
李彦彬 ,
翟夕阳 ,
金学科 ,
何小九

展开

空军工程大学防空反导学院, 西安 710051

余丽山(1992-),男,湖北孝昌人,硕士研究生,研究方向:国防采办与项目管理。

收稿日期: 2017-03-03

网络出版日期: 2025-05-12

收起

Warhead Fragment Optimization Design Based on Overload Damage of the Target Board

YU Lishan ,
LI Yanbin ,
ZHAI Xiyang ,
JIN Xueke ,
HE Xiaojiu

Expand

Air and Missile Defense College, Air Force Engineering University, Xi’an 710051, China

Received date: 2017-03-03

Online published: 2025-05-12

Fold

摘要

针对破片杀伤战斗部破片参数的优化设计问题,对破片飞散特性以及对薄靶板的撞击和侵彻行为进行了分析,基于目标靶板毁伤的整体和局部分析提出了破片冲塞式破坏目标靶板能量消耗和比动能的双目标毁伤准则,并建立了战斗部破片参数优化分析模型。最后通过引入遗传算法对破片打击薄靶目标消耗的能量和毁伤概率进行了估算,表明该模型可有效分析破片长度、直径以及初速对打击薄靶目标的影响,研究结果可为战斗部破片的优化设计提供参考。

关键词： 破片; 优化设计; 能量消耗; 比动能; 遗传算法

本文引用格式

余丽山 , 李彦彬 , 翟夕阳 , 金学科 , 何小九 . 基于靶板冲塞式破坏的战斗部破片优化设计[J]. 弹箭与制导学报, 2019 , 39(5) : 143 -146 . DOI: 10.15892/j.cnki.djzdxb.2019.05.034

Abstract

Aiming at fragment parameters optimization design problem of fragmenting warhead, in this paper, the fragment flying characteristics as well as the impact and penetration behavior were analyzed. Based on the overall and local analysis of the overload damage of the target board, it put forward the double target damage criterion of the energy consumption of the fragment and the kinetic energy per unit area to the target board, and established parameter optimization analysis model of warhead fragment. At last it estimated the energy consumption to the thin plate target and the damage probability by introducing the genetic algorithm. It showed the model could effectively analyze fragment length, diameter and velocity influence on the thin plate target. The research results can provide reference for optimization design of warhead fragment.

Key words： fragment; optimization design; energy consumption; kinetic energy per unit area; genetic algorithm

0 引言

目前地空导弹主要采用半预制或预制破片战斗部。它主要依靠大量高速飞行的破片来命中并击穿来袭目标的关键部件,或者引燃油箱、引爆弹药进而实现对敌武器的毁伤。破片杀伤战斗部对目标的杀伤效果,除了与目标易损性和交会条件相关,还显著与战斗部装药结构尤其是破片的初速和几何尺寸相关。选择合理几何尺寸的破片,并确定其抛撒速度,最大程度的发挥战斗部破片对目标的杀伤效果,同时保证导弹的机动性是地空导弹破片杀伤战斗部设计面临的难题之一。

针对导弹破片杀伤战斗部破片的设计问题,国内外主要采用理论分析、数值模拟和试验手段进行研究。Valerio Cozzani^[1]对实际破片形状和数量进行等效优化分析,通过对舰船毁伤情况的评估来指导战斗部破片的设计;李旭峰^[2]从战斗部的装药结构和材料的角度研究具有高效毁伤作用的含能破片设计方案,并通过实验验证设计的正确性;魏峰^[3]采用数值模拟手段对球形破片在炸药爆炸驱动作用下的飞散特性进行了分析,基于仿真结果对战斗部破片进行了优化。文中针对地空导弹破片杀伤战斗部的优化设计问题,对破片飞散特性以及破片与薄靶作用行为进行了分析,基于破片能量消耗和比动能双重原则,建立了战斗部破片参数优化设计的双目标规划模型,并引入了遗传算法对模型进行求解,最后有效地得到了战斗部破片优化设计的各参数。

1 破片的运动分析

破片杀伤战斗部对目标进行毁伤作用主要是依靠产生的高速运动破片群。当战斗部内的炸药爆炸时,在一瞬间可产生数十万大气压以上的压力,由于战斗部壳体材料的屈服强度远远小于此压力,于是壳体破裂并产生高速运动的破片群,击穿、引燃或引爆目标^[4]。

根据格尼公式,圆柱形壳体战斗部装药爆炸后,静态破片初速度为^[5]:

(1)

v 0 = 2 E β 1 + 0.5 β

式中:

2 E

为格尼常数;

β

为装药质量比。

破片在飞行过程中,其运动方程为:

(2)

m f d v f d t = − 12 c x ρ a A f v f 2

式中:

m f

和

v f

分别表示破片的质量和速度;

c x

表示阻力系数,取决于破片的形状以及飞行速度,文中假设

c x

为常数;

ρ a

表示遭遇点的空气密度;

A f

表示破片的平均有效面积(迎风面积)。令衰减系数

k = c x ρ a A f / (2 m f

),则解得破片飞行到不同距离的存速:

(3)

v f = v 0 e − k x

当战斗部爆炸后形成破片静态初速度为

v 0

时,破片在飞行过程中会受到空气阻力而导致速度衰减,速度衰减系数

k

与破片的形状、破片飞行的速度、空气密度、破片的迎风面积和质量有关。当破片的飞行距离为

x

时,破片的存速为

v f

。

2 破片对薄靶板目标的撞击和侵彻

由图1破片对薄靶板的侵彻模型可知,破片在穿孔时,在破片和目标靶板接触的地方,会形成蘑菇头效应,导致穿孔时孔径变大。

显示原图|下载原图ZIP|生成PPT

图1 破片对薄靶板的侵彻过程

破片高速飞向目标靶板时,在破片和靶板接触的同时挤压目标靶板。由于高速运动破片的作用,使得在和靶板接触的截面上,产生剪应变和剪应力,并伴随着热量生成。瞬间产生的大量热量在短时间内来不及散去,导致接触部位温度迅速升高,从而降低了目标靶板的抗剪能力,产生了冲塞式破坏^[6-7]。破片的动能主要转化为4部分:①剩余破片和击穿靶板后形成塞块的动能;②薄靶板被剪切时消耗的能量;③击穿薄靶板过程中破片损失质量的动能;④破片撞击薄靶板时由于变形所消耗的能量。根据文献[8]中的动量方程和能量方程解得破片的剩余速度为:

(4)

v r = m r m + m r v f 2 − v 502

式中:

m

、

m r

分别表示为塞块的质量、破片的剩余质量;

v 50

表示破片的弹道极限穿透速度,即破片以

v 50

的速度撞击靶板刚好可以穿透。

塞块的质量为:

(5)

m = ρ m A f T s e c θ

式中:

ρ m

为靶板的密度,

k g / m 3

;

A f

为破片的横截面积,

m 2

;

T

为靶板的厚度,

m

;

θ

为破片的入射角度,(°)。

破片贯穿靶板的极限速度为:

(6)

v 50 = C T / d b + K L / d 12

式中:

d

为侵彻体直径,

m

;

L

为侵彻体长度,

m

;C、K、b为经验常数,对于钢质破片侵彻钢靶,其值分别为1 297,-164,0.61,单位都是

m / s

。

破片的剩余质量为:

(7)

m r = m f e − ρ U v f cos θ σ L

式中:

L

为破片的尺寸系数;

σ

为破片的动态屈服强度;

ρ

为破片的密度;

U

为破片轴向塑性波波速。

U = 1.61 σ g / ρ

σ = 3.92 B, N 2 / m m 2

B

为目标靶板的布氏硬度。

3 毁伤准则

破片杀伤战斗部爆炸后,形成大量高速运动、不同形状的破片。每个破片都相当于一个能量体,击穿并摧毁目标的致命性部件。从目标的角度考虑,可将致命性部件等效为目标靶板。确定毁伤准则时,既要从目标靶板的整体(破片能量消耗准则)考虑,又要从目标毁伤的局部(比动能准则)考虑。

3.1 破片能量消耗准则

从目标整体考虑,破片的能量消耗形成对目标靶板的毁伤。破片消耗的能量越多,对目标的毁伤程度越大。

破片消耗的能量为^[9]:

(8)

Δ E = 12 m f v f 2 - 12 m r v r 2

3.2 破片比动能准则

战斗部破片在空气中飞行时除受到自身的重力外,还受到空气阻力。因此,在不断减速的同时也会旋转。从目标毁伤的局部考虑,破片与目标遭遇时接触面积对毁伤的效果也有较大的影响。

单位面积厚度靶板上所接受的比动能为^[10]:

(9)

E b = m f v f 2 2 S h

式中:

S

为破片与靶板的期望接触面积,

S = k m 2 / 3

,圆柱形状系数

k = 3.35 × 10 - 3 m 2 · k g - 2 / 3

;

h

为靶板击穿处的侵彻厚度,

h = T / c o s θ

。

破片对目标靶板进行贯穿损伤时,靶板在单次打击下的毁伤概率经验公式为:

(10)

P = 0, E b < 4.5 1 + 2.65 e − 0.34 E b − 2.96 e − 0.14 E b, E b ≥ 4.5

4 算例验证

4.1 目标函数

战斗部破片的优化设计以考虑对目标的毁伤程度为基础。若对目标毁伤程度最大,应保证破片以消耗最大的能量撞击目标以及破片贯穿目标时对目标的毁伤概率最大,即:

(11)

max Δ E max P

假设战斗部爆炸时,破片的形状为圆柱体,破片的初速度及结构参数在一定的区间内变化,即:

(12)

s . t . v 0 min ≤ v 0 ≤ v 0 max d min ≤ d ≤ d max L min ≤ L ≤ L max

对于双目标优化的问题,可以根据目标的重要性程度不同,将双目标转化为单目标优化。

当目标靶板上接受比动能

E b ≥ 4.5

时,靶板在破片打击下的毁伤概率如图2。

显示原图|下载原图ZIP|生成PPT

图2 毁伤概率曲线

由图2可知,当目标的毁伤概率

P ≤ 95 %

时,随着比动能的增大,目标的毁伤概率迅速增加。但当

P > 95 %

时,尽管比动能大幅度增加,目标的毁伤概率增加的却并不明显。因此可以将

Δ E

作为主目标优化函数,将毁伤概率

P

经过处理后加入到主目标优化函数的约束条件中。具体处理方法:确定常数

ξ

(

ξ

的确定可根据对目标的期望毁伤概率),令P≤ξ,v₀∈

1 000 m / s, 2 000 m / s

,d∈

0.7 c m, 1.4 c m

;L∈

0.7 c m, 1.4 c m

,求主目标函数ΔE的最大值。即:

max Δ E s.t. 1 000 ≤ v 0 ≤ 2000 0.7 ≤ d ≤ 1.4 0.7 ≤ L ≤ 1.4 P < 95 %

4.2 遗传算法

遗传算法是根据生物界的自然选择以及遗传变异等机制,提出的一种具有全局性概率搜索能力的智能算法。它具有很高的智能性、很强的鲁棒性、较高的扩充性^[11-13]。使用遗传算法可有效地对本案例求解。

遗传算法步骤为:

①首先定义适应度函数

f (x) = 12 m f v f 2 - 12 m r v r 2

,在给定空间

U

上搜索。确定种群的规模

N = 10

,代数

T = 3 000

,交叉率

P c = 0.8

,变异率P_m=0.1。

②在搜索空间

U

中随机产生10个个体,

s 1, s 2, …, s 10

,置代数器

t = 1

,初始种群

S (1) = {s 1, s 2, …, s 10}

。

③计算

S

中个体

s i

的适应度函数

f i = f (s i)

。

④确定算法停止条件,若满足,则算法终止,所求结果即为

S

中适应度最大的个体。否则,计算概率:

(13)

P s i = f i ∑ j = 1 10 f j, i = 1, 2, …, 10

按照上述选择概率从

S

中随机选择1个染色体进行复制,共做10次得群体

S 1

。

⑤从

S 1

中按照交叉率

P c

确定进行交叉的染色体

c

,配对进行交叉得群体

S 2

。在群体

S 2

中按照变异率

P m

随机确定

m

个染色体,变异

m

次后得群体

S 3

。

⑥将

S 3

作为新一代种群,用

S 3

替代

S

,t=t+1,转③。

流程图如图3所示。

显示原图|下载原图ZIP|生成PPT

图3 遗传算法流程图

4.3 算例结果

假设钢质圆柱体破片在20 m处正对钢质目标靶板飞去,即破片的入射角度

θ = 0 °

,钢的密度

ρ = 7.85 g · c m - 3

,阻力系数

c x = 1.12

,薄板的厚度为2 mm。得到破片撞击目标靶板消耗的能量随破片的各参数变化迭代次数如图4。

显示原图|下载原图ZIP|生成PPT

图4 破片能量消耗随迭代次数的变化

由图4可知,随着破片的各参数变化,破片撞击目标靶板消耗的能量增加,迭代到833次后趋于稳定,并稳定在全局最优处。基于战斗部破片各参数设计要求,选取3组有代表性数据如表1所示。

表1 破片参数

破片参数	破片组序号
破片参数	1	2	3
$d / c m$	1.16	1.32	1.39
$L / c m$	1.30	1.17	0.80
$v 0 / (m · s - 1)$	1 333.73	1 532.81	1 684.50
$P$	0.80	0.93	0.94
$Δ E / J$	2 734.76	4 709.41	5 553.91

通过对表1战斗部破片的参数分析可知,当破片的直径为

1.39 c m

,长度为

0.8 c m

,战斗部爆炸后的初速度为

1 684.5 m / s

时,破片对目标的毁伤概率最大,为

94 %

,破片撞击目标靶板消耗的能量最多,为

5 553.91 J

。

5 结论

1)对破片飞散特性以及破片与薄靶板的撞击与侵彻行为进行了分析,基于破片能量消耗和比动能双重毁伤准则,建立了破片参数优化分析的双目标规划模型;

2)运用遗传智能算法对本模型进行求解,可以有效控制目标函数的变化趋势,使目标函数最终稳定在全局最优处,并得到全局最优解;

3)算例与分析表明,在破片参数最优的情况下,单个破片对目标的毁伤概率为

94 %

,破片撞击目标靶板消耗的能量为

5 553.91 J

。若继续增加破片的质量和速度,毁伤概率几乎没有增加,即毁伤作用并没有增大。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序

[1]	GUBINELLI Gianfilippo, COZZANI Valerio. Assessment of missile hazards: evaluation of the fragment number and drag factors[J]. Journal of Hazardous Materials, 2009, 161(1):439-449. DOI PMID

[2]	李旭锋. 含能破片对模拟战斗部的引爆机理研究[D]. 南京: 南京理工大学, 2006.

[3]	魏峰. 某战斗部预制破片设计及破片威力数值模拟[D]. 沈阳: 东北大学, 2010.

[4]	韩晓明, 高峰. 导弹战斗部原理及应用[M]. 西安: 西北工业大学出版社, 2012:92-96.

[5]	GONG Chaoan, Chen Zhigang, YIN Likui. Analysis of movement laws of fragment and shock wave from a blast fragmentation warhead[J]. Journal of Measurement Science and Instrumentation, 2015, 6(3):218-222.

[6]	熊冉, 高欣宝, 许兴春, 等. 破片侵彻金属薄板后的剩余速度研究[J]. 爆破, 2013, 30(4):41-44.

[7]	洪建华, 刘彤. 杀伤破片侵彻击穿和引爆靶弹的分析与研究[J]. 兵工学报, 2003, 24(3):316-321.

[8]	李向东, 杜忠华. 目标易损性[M]. 北京: 北京理工大学出版社, 2013:173-180.

[9]	陈长海, 朱锡, 侯海量, 等. 舰船舷侧复合装甲结构抗动能穿甲模拟实验[J]. 爆炸与冲击, 2011, 31(1):11-18.

[10]	高波, 党东东, 韩凯, 等. 飞机易损性定量分析[J]. 海军航空工程学院学报(军事版), 2013, 11(2):43-47.

[11]	PAL Sankar K, BANDYOPADHYAY Sanghamitra, MURTHY C A. Genetic algorithms for generation of class boundaries[J]. IEEE Transactions on Systems, Man, and Cybernetics: Part B Cybernetics, 1998, 28(6):816-828.

[12]	BÄCK T, FORGEL D B, MICHALEWICZ Z. Handbooks of Evolutionary Computation[M]. New York: Oxford University Press, 1997.

[13]	FONSECA C M, FLEMING P J. Multiobjective optimization and multiple constraint handling with evolutionary algorithm:I A unified formulation[J]. IEEE Transactions on Systems, Man, and Cybernetics:Part A Systems and Humans, 1998, 28(1):26-37.

Options

文章导航

模态框（Modal）标题

摘要

本文引用格式

Abstract

0 引言

1 破片的运动分析

2 破片对薄靶板目标的撞击和侵彻

图1 破片对薄靶板的侵彻过程

3 毁伤准则

3.1 破片能量消耗准则

3.2 破片比动能准则

4 算例验证

4.1 目标函数

图2 毁伤概率曲线

4.2 遗传算法

图3 遗传算法流程图

4.3 算例结果

图4 破片能量消耗随迭代次数的变化

表1 破片参数

5 结论

参考文献