Terminal hypersonic trajectory modeling and optimization for maneuvering penetration and precision strike
-
摘要:
针对高超声速飞行器再入末段机动突防、精确打击问题，从最优控制角度出发，提出了一种考虑拦截弹动力学特性的最优机动突防弹道优化方法，获得了高超声速飞行器的最大机动能力。该方法将拦截弹运动模型引入突防弹道优化的模型中，通过施加约束限制拦截弹按照比例导引律飞行。根据飞行任务和交战双方的弹道特点分段，结合各段的任务和特性，分别提出了突防性能指标和精确打击性能指标等，并通过加权函数将各个独立、矛盾的性能指标统一，建立了多对象、多段、多约束机动突防弹道优化模型，采用Radau多段伪谱法（MRPM）进行求解。针对该问题求解的初值敏感、可行域窄等问题，提出了一系列弹道优化策略，提高了收敛速度和求解精度，最终获得了最优机动弹道，并通过协态映射原理对其最优性进行了验证。结果表明，该方法能充分发挥高超声速飞行器的机动能力，获得满足落点精度要求的突防弹道，相对已有方法，将脱靶量提高了1~2个量级。灵敏度分析表明，该弹道对拦截弹的发射时间不敏感。
Abstract:Aimed at the maneuvering penetration and precision strike problem of hypersonic vehicle terminal trajectory, an optimal maneuvering trajectory optimization method considering the dynamic characteristics of intercepting was proposed from the viewpoint of optimal control, so as to obtain the maximum maneuverability of hypersonic vehicles. In this paper, the intercepting missile model was introduced into the model of penetration trajectory optimization, and a constraint was imposed to restrict the intercepting missile to fly according to the proportional guidance law. The trajectories were divided into phases according to the different missions and trajectory characteristics of the belligerents. The penetration performance index and the precision strike performance index are put forward according to the task and characteristics of each phase, and by the weighting function the independent and contradictory performance indicators are unified. Thus a multi-object, multi-phase and multi-constrained maneuvering penetration trajectory optimization model was established. And multiphase Radau pseudospectral method (MRPM) was used to solve the problem. Due to the initial sensitivity and narrow feasible region of the problem, a series of trajectory optimization strategies were proposed to improve the convergence rate and the precision of the solution. Finally, the optimal maneuvering trajectory was obtained, and the optimality of the solution was verified based on the principle of costate mapping. The results show that the method can give full play to the maneuverability of the hypersonic vehicle, and obtain the penetration trajectory which satisfies the terminal accuracy. Compared with the existing method, the miss distance is increased by 1-2 orders of magnitude. Sensitivity analysis shows that the trajectory is insensitive to the launch time of the interceptor.
-
Key words:
- maneuvering /
- penetration /
- interception /
- optimization strategy /
- optimality condition /
- sensitivity analysis
-
图 1 高超末段机动突防-对地打击示意图
Figure 1. Diagram of maneuvering penetration and ground attack in terminal phase of hypersonic trajectory
图 3 平滑性能指标对控制曲线的影响
Figure 3. Influence of smoothness performance index on control curves
图 12 突防弹攻角-时间曲线和倾侧角-时间曲线
Figure 12. Penetration missile attack angle-time curve and bank angle-time curves
图 13 不同机动形式突防弹攻角-时间曲线和倾侧角-时间曲线
Figure 13. Penetration missile attack angle-time curve and bank angle-time curves of different motivations
图 15 热流密度-时间，动压-时间，过载-时间曲线
Figure 15. Heat flux-time, dynamic pressure-time and overload-time curves
表 1 突防弹初始及终端状态约束
Table 1. Initial and terminal state constraints of penetration missile
时间 xA/m yA/m zA/m VA/(m·s-1) θA/(°) ψA/(°) αA/(°) σA/(°) t0 10 188.8 21 299.3 10 188.8 1 878.8 -18.9 135.0 10 -180 tf -9 346 0 -9 346 (800, 1 100) (-90, -75) (0, 7) 
下载: 导出CSV
表 2 拦截弹初始状态约束
Table 2. Initial state constraints of interception missile
状态约束 xD0/
myD0/
mzD0/
mVD0/
(m·s-1)θD0/
(°)ψD0/
(°)mD0/
kg数值 -8 000 0 -8 000 1 000 38 -45 315
下载: 导出CSV
表 3 突防弹常规路径约束
Table 3. Path constraints of penetration missile
路径约束 
q/(N·m-2) ny 数值 400 6×105 20
下载: 导出CSV
表 4 优化结果终端状态
Table 4. Terminal states of optimal results
性能指标 xAf/m yAf/m zAf/m VAf/(m·s-1) θAf/(°) ψAf/(°) αAf/(°) σAf/(°) J1 -9 346.99 -1.00 -9 346.99 1 450.67 -84.37 132.49 6.00 -178.64 J2 -9 344.99 -1.00 -9 344.99 1 515.68 -75.01 135.00 6.00 -180.01 J3 -9 344.99 -1.00 -9 346.99 1 377.61 -75.01 180.01 6.00 -210.65 J4 -9 344.99 -1.00 -9 346.99 1 018.95 -90.01 180.01 6.00 -180.20 常值 -9 345.99 -1.81×10-12 -9 345.99 1 462.12 -84.62 135.01 10.00 -180.00
下载: 导出CSV
表 5 突防脱靶量对比
Table 5. Comparison of penetration miss distance
性能指标 脱靶量/m J1 19.00 J2 5.18 J3 216.65 J4 234.59 常值 3.69
下载: 导出CSV
表 6 不同机动形式突防脱靶量对比
Table 6. Comparison of penetration miss distance of different motivations
机动形式 最优机动 程序机动 脱靶量/m 234.6 119.7
下载: 导出CSV
表 7 求解精度对比
Table 7. Comparison of solution accuracy
求解精度 GPOPS初值生成器 分段嵌套优化 节点数 41 105 节点误差 0.293 5 8.313 9×10-5 可行性 2.0×10-2 3.2×10-8 最优性 9.4×10-3 7.2×10-6
下载: 导出CSV
-
[1] 陈万春, 聂蓉梅, 刘佳琪, 等.PAC-3爱国者拦截弹末制导精度仿真研究[J].飞航导弹, 1999, 19(7):57-62. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=fhdd907.013&dbname=CJFD&dbcode=CJFQCHEN W C, NIE R M, LIU J Q, et al.Simulation and research on the terminal guidance precision of PAC-3 patriot missile[J].Winged Missiles Journal, 1999, 19(7):57-62(in Chinese). http://kns.cnki.net/KCMS/detail/detail.aspx?filename=fhdd907.013&dbname=CJFD&dbcode=CJFQ [2] 崔静, 姜玉宪.拦截导弹动力学特性对摆动式机动策略突防效果的影响[J].宇航学报, 2001, 22(5):33-38. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=yhxb200105005&dbname=CJFD&dbcode=CJFQCUI J, JIANG Y X.The effect of interceptor's dynamic system order on the penetration efficiency of weaving maneuver strategy[J].Journal of Astronautics, 2001, 22(5):33-38(in Chinese). http://kns.cnki.net/KCMS/detail/detail.aspx?filename=yhxb200105005&dbname=CJFD&dbcode=CJFQ [3] 姜玉宪, 崔静.导弹摆动式突防策略的有效性[J].北京航空航天大学学报, 2002, 28(2):133-136. http://bhxb.buaa.edu.cn/CN/abstract/abstract10821.shtmlJIANG Y X, CUI J.Effectiveness of weaving maneuver strategy of a missile[J].Journal of Beijing University of Aeronautics and Astronautics, 2002, 28(2):133-136(in Chinese). http://bhxb.buaa.edu.cn/CN/abstract/abstract10821.shtml [4] ZARCHAN P.Proportional navigation and weaving targets[J].Journal of Guidance, Control, and Dynamics, 1995, 18(5):969-974. doi: 10.2514/3.21492 [5] SHINAR J, STEIBERG D.Analysis of optimal evasive maneuvers based on a linearized two-dimensional kinematic model[J]. Journal of Aircraft, 1977, 14(8):795-802. doi: 10.2514/3.58855 [6] SHINAR J, TABAK R.New results in optimal missile avoidance[J]. Journal of Guidance, Control, and Dynamics, 1994, 17(5):897-902. doi: 10.2514/3.21287 [7] YANUSHEVSKY R.Analysis of optimal weaving frequency of maneuvering targets[J].Journal of Spacecraft and Rockets, 2004, 41(3):477-479. doi: 10.2514/1.6459 [8] IMADO F.Some aspects of a realistic three-dimensional pursuit-evasion game[J].Journal of Guidance, Control, and Dynamics, 1993, 16(2):125-139. doi: 10.2514/3.21002 [9] IMADO F, UEHARA S.High-g barrel roll maneuvers against proportional navigation from optimal control, viewpoint[J]. Journal of Guidance, Control, and Dynamics, 1998, 21(6):876-881. doi: 10.2514/2.4351 [10] IMADO F, KURODA T. Engagement tactics for two missiles against an optimally maneuvering aircraft[J].Journal of Guidance, Control, and Dynamics, 2011, 34(2):574-582. doi: 10.2514/1.49079 [11] BENSON D A, THORVALDSEN G T, RAO A V.Direct trajectory optimization and costate estimation via an orthogonal collocation method[J].Journal of Guidance, Control, and Dynamics, 2006, 29(6):1435-1440. doi: 10.2514/1.20478 [12] HUNTINGTON G T, RAO A V.Optimal reconfiguration of spacecraft formation via a Gauss pseudospectral method[J].Journal of Guidance, Control, and Dynamics, 2008, 31(3):689-698. doi: 10.2514/1.31083 [13] ZHANGK N, CHEN W C. Reentry vehicle constrained trajectory optimization[C]//Proceedings of AIAA International Space Planes and Hypersonic Systems and Technologies Conference. Reston: AIAA, 2011, 1: 1-16. [14] DARBY C L, RAO A V.Minimum-fuel low earth orbit aeroassisted orbital transfer of small spacecraft[J].Journal of Spacecraft and Rocket, 2011, 48(4):618-628. doi: 10.2514/1.A32011 [15] WEISS M, SHIMA T, CASTANEDA D, et al.Minimum effort intercept and evasion guidance algorithms for active aircraft defense[J]. Journal of Guidance, Control, and Dynamics, 2016, 39(10):2297-2311. doi: 10.2514/1.G000558 [16] DIVIA G. Advances in global pseudospectral methods for optimal control[D]. Gainesville: University of Floroda, 2011: 92-96. [17] ZARCHAN P.Tactical and strategic missile guidance[M]. 6th ed. Reston:AIAA, 2012:117-122. -
下载:
点击查看大图
计量
- 文章访问数: 812
- HTML全文浏览量: 108
- PDF下载量: 477
- 被引次数: 0
