中国科技核心期刊      中国指挥与控制学会会刊     军事装备类重点期刊
Theory & Research

Research on a New Maneuver Penetration Strategy of Ballistic Missile

  • XIE Ru-heng ,
  • DING Yang ,
  • YANG Yi
Expand
  • The 28th Research Institute of China Electronics Technology Group Corporation, Nanjing 210001, China

Received date: 2021-01-26

  Revised date: 2021-02-22

  Online published: 2021-06-10

Abstract

At present, many countries are studying ballistic missiles and ballistic missile interception technologies. In particular, the United States recently launched an interceptor in the Pacific to successfully intercept intercontinental missiles, which has greatly attacked the deterrence of intercontinental ballistic missiles. This paper aims to propose a new ballistic missile penetration strategy. Firstly, the offensive and defensive confrontation model of the ballistic missile and the interceptor is constructed on the pitch plane, and the models of the offensive and defensive confrontation system of the two parties are respectively constructed based on the optimal control theory. The heading error is defined as the performance index, and the above model is solved to obtain the optimal maneuverable penetration and interception strategy of the ballistic missile and interceptor in the pitch plane. The theoretical results are simulated and verified, and the effects of the interceptor’s initial lead angle, ballistic missile maneuverability, maneuver timing and maneuver duration on the final miss are analyzed. The simulation results have successfully verified the correctness and effectiveness of the maneuver strategy in this paper, which has certain reference value.

Cite this article

XIE Ru-heng , DING Yang , YANG Yi . Research on a New Maneuver Penetration Strategy of Ballistic Missile[J]. Command Control and Simulation, 2021 , 43(3) : 12 -17 . DOI: 10.3969/j.issn.1673-3819.2021.03.003

近年来,随着全球军事强国导弹防御系统的日趋成熟,战略核力量面临的巨大挑战就是大气层外中段防御系统,为有效躲避敌方反导防御系统中外大气层杀伤拦截器(Exoatmospheric Kill Vehicle,EKV)的拦截[1-3],弹道导弹中段机动变轨技术已成为重点研究方向。
美国国家导弹防御系统的核心技术主要体现在外大气层杀伤拦截器(EKV)上。外大气层杀伤拦截器是一种小型、非核、自寻的、靠直接碰撞杀伤的拦截器,它还装有红外导引头、变轨推进系统以及数据处理制导系统[4-6]。美国多年来通过拦截试验对EKV进行了多次改进和完善,2020年在太平洋发射的拦截器成功拦截洲际弹道导弹,使其具备实战应用的能力,因此研究分析外大气层杀伤拦截器的作战特点和突防措施,具有一定的现实意义。
目前的弹道导弹中段突防措施,主要分为以下两个方面。
第一种是程序式机动突防[7-8],所谓机动突防,就是事先设定好弹道导弹的机动程序,即不管敌方是否拦截,如何拦截,弹道导弹都会在发射后的固定时间,进行固定的机动,以此来降低拦截器对弹道导弹轨迹的预测精度,提高弹道导弹的突防成功率。
第二种是自主机动突防[9-13],突防弹道导弹在探测到拦截器时,通过获取拦截器的飞行参数,由弹载计算模块实时计算突防机动指令进行突防机动[13-17]。这种机动突防方式具有主动智能化的特点,能够针对不同的拦截方式给出最优的突防策略,大大提高突防的成功率。但是由于整个突防过程时间短,对弹道导弹探测模块、计算模块以及发动机要求较高,所以具备自主机动突防能力的弹道导弹非常少。
本文基于最优控制和弹道导弹与拦截器在俯仰平面的对抗模型,构建了一种新的机动突防策略,并用Matlab仿真验证了机动突防策略的正确性。

1 俯仰平面对抗模型建立

本文构建弹道导弹和拦截器在俯仰平面的攻防对抗模型,如图1所示。
图1 俯仰平面攻防对抗示意图
图1中,aE表示拦截器机动加速度,垂直于速度方向,左侧为正,右侧为负;θE表示拦截器的弹道倾角,以x轴正方向起始,逆时针为正,顺时针为负;σE表示拦截器的前置角,弹体x轴与弹目视线之间的夹角,向左为正,向右为负;VE表示拦截器的飞行速度;gE表示拦截器的重力加速度;aM表示弹道导弹的机动加速度,垂直于速度方向,右侧为正,左侧为负;θM表示弹道导弹的弹道倾角,以x轴负方向起始,顺时针为正,逆时针为负;σM表示弹道导弹的前置角,弹体x轴与弹目视线之间的夹角,向右为正,向左为负;VM表示弹道导弹的飞行速度,gM表示弹道导弹的重力加速度;r表示拦截器与弹道导弹的距离;μT表示拦截器相对于弹道导弹的弹目视线高低角,以x轴正方向起始,逆时针为正,顺时针为负。
拦截器动力学方程建立如下:
x ¨ E = - a E s i n θ E y ¨ E = a E c o s θ E - g E θ ˙ E = ( a E - g E c o s θ E ) / V E V ˙ E = - g E s i n θ E
弹道导弹的动力学方程如下:
x ¨ M = a M s i n θ M y ¨ M = a M c o s θ M - g M θ ˙ M = ( a M - g M c o s θ M ) / V M V ˙ M = - g M s i n θ M
拦截器与弹道导弹相对距离r和接近速度VXD计算如下:
r= ( x M - x E ) 2 + ( y M - y E ) 2
VXD=VEcos σ E+VMcos σ M
攻防对抗的时间tf采用下式进行估计计算:
tf= r 0 V X D 0
式中,r0表示攻防对抗初始时刻的相对距离,VXD0表示攻防对抗初始时刻的接近速度。
剩余飞行时间计算如下:
tgo=tf-t
拦截器和弹道导弹的前置角、弹道倾角和弹目视线角的关系如下:
σ E = θ E - μ T σ M = θ M + μ T
弹目视线高低角计算如下:
r c o s μ T = x M - x E r s i n μ T = y M - y E

2 模型的求解

定义航向误差τ如下:
τ=VMsin σ M-VEsin σ E=VE V M s i n σ M V E - s i n σ E
sin ζ= V M s i n σ M V E
则式(9)变为
τ=VMsin σ M-VEsin σ E=VE(sin ζ-sin σ E)
航向误差τ表示垂直于弹目视线方向上双方的速度之差;τtgo表示在没有控制输入情况下,垂直弹目视线方向上的瞬时脱靶量。
弹目视线角速度可表示如下:
μ ˙ T= τ r
τ求导数可以得到:
τ ˙= U - M- U - E
式中 U - M U - E是中间变量,表示如下:
U - M = ( a M - g M c o s θ M ) c o s σ M - g M s i n θ M s i n σ M +   τ r V M c o s σ M U - E = ( a E - g E c o s θ E ) c o s σ E - g E s i n θ E s i n σ E -   τ r V E c o s σ E
式(13)中的aMaE分别表示弹道导弹和拦截器的机动加速度。由于
μTEEMM
c o s μ T = c o s θ E c o s σ E + s i n θ E s i n σ E c o s μ T = c o s σ M c o s θ M + s i n σ M s i n θ M
U - M = ( a M - g M c o s θ M ) c o s σ M - g M s i n θ M s i n σ M + τ r V M c o s σ M =   a M c o s σ M - g M c o s μ T + τ r V M c o s σ M U - E = ( a E - g E c o s θ E ) c o s σ E - g E s i n θ E s i n σ E - τ r V E c o s σ E =   a E c o s σ E - g E c o s μ T - τ r V E c o s σ E
选取系统的性能指标为
J=|τ|
弹道导弹采取策略使得J尽可能大以便成功突防,而动能拦截器则相反,采取策略使得J尽可能小以便成功拦截弹道导弹。
根据最优控制理论中极大值极小值原理并结合 τ ˙= U - M- U - E,以航向误差τ作为控制变量,很容易得到双方的最优策略:
U - E= U - E m a x ,   τ ( t f ) > 0 U - E m i n ,   τ ( t f ) < 0
U - M= U - M m a x ,   τ ( t f ) > 0 U - M m i n ,   τ ( t f ) < 0
上式中, U - E m a x表示 U - E的最大值; U - E m i n表示 U - E的最小值; U - M m a x表示 U - M的最大值; U - M m i n表示 U - M的最小值;详细表达如下:
UEmax= a E m a x c o s σ E - g E c o s μ T - τ r V E c o s σ E ,   c o s σ E > 0 - a E m a x c o s σ E - g E c o s μ T - τ r V E c o s σ E ,   c o s σ E < 0
UMmax= a M m a x c o s σ M - g M c o s μ T + τ r V M c o s σ M ,   c o s σ M > 0 - a M m a x c o s σ M - g M c o s μ T + τ r V M c o s σ M ,   c o s σ M < 0
UEmin= - a E m a x c o s σ E - g E c o s μ T - τ r V E c o s σ E ,   c o s σ E > 0 a E m a x c o s σ E - g E c o s μ T - τ r V E c o s σ E ,   c o s σ E < 0
UMmin= - a M m a x c o s σ M - g M c o s μ T + τ r V M c o s σ M ,   c o s σ M > 0 a M m a x c o s σ M - g M c o s μ T + τ r V M c o s σ M ,   c o s σ M < 0
其中,aEmax表示动能拦截器的最大指令加速度;aMmax表示弹道导弹的最大指令加速度。
可令
U M *= a M * cos σ M-gMcos μ T+ τ rVMcos σ M
U E *= a E * cos σ E-gEcos μ T- τ rVEcos σ E
其中 a M * a E *是最优加速度指令,结合式(20)、(21)、(22)和(23),可得如下最优加速度指令:
a M *=aMmaxsgn(τ(tf)cos σ M)
a E *=aEmaxsgn(τ(tf)cos σ E)
结合式(12)、(24)、(25),假设动能拦截器和弹道导弹的机动加速度分别为零,可得
aE=aM=0
所以
τ ˙= U - M- U - E=(gE-gM)cos μ T+ τ r(VMcos σ M+VEcos σ E)≈ τ r(VMcos σ M+VEcos σ E)
对上式进一步化简可得
τ ˙= τ r(VMcos σ M+VEcos σ E)=τ/ r V M c o s σ M + V E c o s σ E= τ t g o τ t f - t
将上述微分方程对时间t进行积分可得
τ=-C1 e l n ( t f - t )
式中,C1为常数,且C1不为零;从上式可以看出τ的符号取决于C1,且C1为常数,所以
sgn(τ(tf))=sgn(τ(t))
因此式(26)和(27)变为:
a M *=aMmaxsgn(τ(t)cos σM)
a E *=aEmaxsgn(τ(t)cos σE)
从上式中可以看出,最终弹道导弹和拦截器的攻防对抗结果只和各自的机动能力和探测范围有关,机动能力表示导弹的最大机动加速度和机动时长,探测范围则决定了导弹的机动时机。

3 仿真分析

3.1 弹道导弹不机动

首先在弹道导弹不进行机动的情况下进行突防仿真验证,因为拦截器的拦截速度很大,在机动加速度限制在4g时,无法使飞行弹道快速改变,所以必然要求具有一定的前置角才能成功拦截弹道导弹。表1给出了攻防对抗最终脱靶量和拦截器初始前置角的关系。
表1 拦截器初始前置角对突防效果影响关系表
前置角/(°) 脱靶量/m 前置角/(°) 脱靶量/m
0 5 454.99 8 1.34
1 3 252.45 9 1.10
2 1 031.19 10 3.62
3 2.12 11 1.77
4 2.19 12 0.7
5 0.067 13 1.14
6 1.98 14 305.41
7 3.69 15 2632.45
表1图2中的仿真结果可以看出,弹道导弹在无机动情况下的突防效果与拦截器攻防对抗初始时刻的初始前置角有很大关系。只有当拦截器攻防对抗初始时刻的初始前置角在一定的范围之内时才可成功拦截弹道导弹;当拦截器的初始前置角满足3°≤ψM0≤13°情况下才可成功拦截弹道导弹,称该前置角范围为拦截器的最佳拦截窗口。
图2 拦截器初始前置角对突防效果影响关系图

3.2 弹道导弹机动

本节在弹道导弹机动的情况下,仿真分析弹道导弹在不同机动加速度、不同机动时机以及不同机动时长情况下的最终脱靶量,分析弹道导弹在何种情况下能够成功突防。设置拦截器的初始前置角为8°,弹道导弹最大机动加速度分别为2g、3g、4g,最大机动时长7.25 s,最远探测距离100 km,不同情况下的突防效果结果如图3表2所示。
图3 弹道导弹最大机动加速度为2g情况下突防结果
表2 弹道导弹最大机动加速度为2g下突防结果
T/s R/km
10 20 30 40 50 60 70 80 90 100
0.25 0.28 0.97 2.3 3.7 2.74 1.20 0.41 2.14 3.9 2.17
1.25 1.63 0.4 2.01 3.25 3.79 3.27 3.11 3.3 3.9 3.1
2.25 1.63 2.62 2.85 0.67 3.07 1.75 0.7 2.51 3.63 3.78
3.25 1.63 3.17 2.31 0.25 1.03 2.76 2.49 1.1 2.28 3.23
4.25 1.63 3.17 5.17 4.8 0.91 0.98 1.64 1.08 0.72 3.77
5.25 1.63 3.17 5.17 9.76 4.32 0.92 2.49 0.51 1.97 1.92
6.25 1.63 3.17 5.17 9.76 14.31 4.69 2.63 3.46 1.8 2.69
7.25 1.63 3.17 5.17 9.76 14.35 19.66 3.19 3.54 1.9 2.36
由仿真结果可知,在弹道导弹机动加速度为2g的情况下,弹道导弹至少需要机动6.25 s,机动时机弹目距离在50 km左右才可以成功突防,而且此时的脱靶量还不到20 m,其余情况均被拦截器成功拦截。
相比于最大机动加速度为2g的突防结果,最大加速度为3g的突防成功率大大提升,最短机动3.25 s便可以成功突防,此时机动时刻弹目距离为30 km。由表3可知,最大机动加速度的提升能够有效提高弹道导弹的突防成功率。
表3 弹道导弹加速度为3g情况下突防结果
T/s R/km
10 20 30 40 50 60 70 80 90 100
0.25 0.25 2.12 3.81 1.72 0.48 2.79 2.71 0.12 2.51 2.65
1.25 4.82 3.80 2.73 0.8 1.67 3.23 0.32 3.51 1.18 1.57
2.25 4.82 14.85 1.29 0.13 1.9 3.28 0.43 2.82 2.88 1.58
3.25 4.82 15.95 31.54 15.93 2.36 3.88 1.16 3.11 0.89 0.1
4.25 4.82 15.95 35.82 54.13 36.64 3.47 3.22 0.31 3.49 3.15
5.25 4.82 15.95 35.82 63.48 82.46 63.12 5.57 3.88 1.08 3.92
6.25 4.82 15.95 35.82 63.48 98.95 116.71 95.44 35.34 0.94 3.94
7.25 4.82 15.95 35.82 63.48 98.98 140.89 156.83 133.7 71.56 3.3
表4中可以看出在弹道导弹最大机动加速度为4g情况下最少只需机动2.25 s左右即可成功突防,此时机动时刻弹目距离在20 km左右。在机动时长固定时,机动时机选择很重要,过早和过晚均会被拦截器拦截,过早机动,由于机动时长一定,使得拦截器有足够的时间进行机动拦截;过晚机动,弹道导弹来不及改变飞行状态而被成功拦截。定义机动时长一定情况下,能够成功规避拦截器拦截的机动时刻弹目距离为最佳逃逸窗口。从图5中可以看出,在机动时长一定的情况下,最佳逃逸窗口中存在峰值,并定义该峰值为最佳逃逸时机,从图中可以看出,最佳逃逸时机随着机动时长增加向右偏移,即机动时长越长,机动时机就越早。
表4 弹道导弹加速度为4g情况下突防结果
T/s R/km
10 20 30 40 50 60 70 80 90 100
0.25 0.770 3.26 2.01 0.760 3.68 1.16 2.05 2.42 1.09 3.17
1.25 9.92 4.95 3.72 1.28 0.92 2.4 3.16 3.22 2.47 1.07
2.25 9.93 35.32 36.26 1.60 2.45 2.14 0.54 2.34 1.34 3.78
3.25 9.93 36.81 76.89 83.16 45.3 1.36 3.89 0.37 1.44 1.34
4.25 9.93 36.81 82.55 134.12 145.68 112.96 36.05 0.23 0.73 0.43
5.25 9.93 36.81 82.55 146.60 206.99 223.85 196.25 124.31 8.3 2.46
6.25 9.93 36.81 82.55 146.60 228.94 295.54 317.69 295.25 228.23 116.9
7.25 9.93 36.81 82.55 146.60 228.88 327.81 399.83 427.28 410 347.91
图4 弹道导弹加速度为3g情况下突防结果
图5 弹道导弹加速度为4g情况下突防结果

4 结束语

本文构建了弹道导弹和拦截器在俯仰平面的攻防对抗,并以航向误差τ作为攻防对抗双方性能指标。拦截器欲使航向误差最小以成功拦截弹道导弹;而弹道导弹欲使航向误差最大以躲避拦截器的拦截。基于此对上述攻防对抗模型进行求解,获取双方在二维平面内的最优机动策略。并在Matlab仿真环境下仿真分析了拦截器初始前置角、弹道导弹机动能力、可机动时间和机动时机对突防效果的影响,结合仿真分析总结如下结论。
1)拦截器对弹道导弹的拦截存在最佳拦截窗口,只有拦截器初始前置角位于最佳拦截窗口内才能成功拦截弹道导弹。
2)弹道导弹在机动时长和机动加速度一定时突防存在最佳逃逸窗口和最佳逃逸时机,在最佳逃逸窗口中突防成功率较大,机动时长越长,最佳逃逸时机越早。
3)弹道导弹的机动能力越强,越有利于成功规避拦截器的拦截。
4)在弹道导弹机动能力和机动时长一定时,并非越早机动越好,也非越晚机动越好,选择合适的机动时机仍有可能成功规避拦截器的拦截。越早机动,机动时长一定,拦截器有足够的时间进行机动拦截;越晚机动,由于速度太大无法迅速改变飞行轨迹而被拦截成功。
5)基于本文设计的突防机动策略,弹道导弹在某些情况下只需机动2.25 s便可成功突破拦截器的拦截,为突防后的弹道回归预留了足够的时间。
上述仿真结果验证了本文机动突防策略的正确性和可行性,为发展弹道导弹中段自主机动突防提供理论参考。
[1]
方喜龙, 刘新学, 张高瑜, 等. 国外弹道导弹机动突防策略浅析[J]. 飞航导弹, 2011(12):17-22.

[2]
鲜勇, 斯文辉. 对GMD系统拦截时机与空域的分析与研究[J]. 火力与指挥控制, 2011, 36(5):192-194,200.

[3]
吴启星, 张为华. 弹道导弹中段机动突防研究[J]. 宇航学报, 2006(6):1243-1247.

[4]
罗珊, 方群, 涂良辉. 弹道导弹反动能拦截器机动突防研究[J]. 弹箭与制导学报, 2006, 26(2):1121-1124.

[5]
王刚, 周凤岐, 周军. 远程弹道导弹大气层外机动突防[J]. 弹箭与制导学报, 2006, 26(1):364-366,370.

[6]
雍恩米, 唐国金, 罗亚中. 弹道导弹中段机动突防制导问题的仿真研究[J]. 导弹与航天运载技术, 2005(4):13-18.

[7]
张广明, 高普云, 唐乾刚, 等. 弹道导弹在被动段突防的平面脉冲式变轨方法[J]. 弹道学报, 2008, 20(3):83-86,95.

[8]
管千山, 谷良贤, 阎代维. 跳跃式弹道导弹突防能力[J]. 火力与指挥控制, 2008, 33(1):84-85,89.

[9]
Shima T, Weiss H. Differential Game Guidance Considering EKV’s Direct Lift and Range Dependent Estimation Delay[C]// IEEE Conference on Decision and Control, 2002.

[10]
Anderson G M. Comparison of Optimal Control and Differential Game Intercept Missile Guidance Laws[J]. Journal of Guidance Control & Dynamics, 1981, 4(2):109-115.

[11]
Shima T, Shinar J. Time-Varying Linear Pursuit-Evasion Game Models with Bounded Controls[J]. Journal of Guidance Control & Dynamics, 2002, 25(3):425-432.

[12]
Oshman Y, Shinar J, Weizman S A. Using a Multiple-Model Adaptive Estimator in a Random Evasion Missile/Aircraft Encounter[J]. Journal of Guidance Control & Dynamics, 2001, 24(6):1176-1186.

[13]
Turetsky V. Capture Zones of Linear Feedback Pursuer Strategies[J]. Automatica, 2008, 44(2):560-566.

DOI

[14]
Turetsky V. Upper Bounds of the Pursuer Control Based on a Linear-Quadratic Differential Game[J]. Journal of Optimization Theory & Applications, 2004, 121(1):163-191.

[15]
鲜勇, 田海鹏, 王剑, 等. 基于微分对策的导弹智能机动突防研究[J]. 飞行力学, 2014, 32(1):70-73+83.

[16]
梁建波. 弹道导弹实时机动突防与制导方法研究[D]. 西安: 第二炮兵工程大学, 2011.

[17]
雷刚, 田海鹏, 鲜勇, 等. 基于博弈论的导弹机动突防策略研究[J]. 固体火箭技术, 2014, 37(3):291-294.

Outlines

/