导弹与航天运载技术（中英文）

统计类别	落地倾角/(°)
均值	-80.164
均方差	0.0837

统计类别	落地倾角/(°)
均值	-80.164
均方差	0.0837

远程机动飞行器地心角修正的制导方法

PDF下载

王煜东 , 方岳 , 郭珂 , 颜楚雄 , 李瑾

导弹与航天运载技术 | 导航、制导与控制 2024,47(6): 42-46

收起

导弹与航天运载技术 | 导航、制导与控制 2024, 47(6): 42-46

远程机动飞行器地心角修正的制导方法

全屏

王煜东, 方岳, 郭珂, 颜楚雄, 李瑾

作者信息

北京航天长征飞行器研究所,北京,100076

王煜东(1994—),男,工程师,主要研究方向为导航、制导与控制。

方岳(1978—),男,博士,研究员,主要研究方向为飞行器总体设计。

郭珂(1983—),男,博士,高级工程师,主要研究方向为飞行器总体设计,导航、制导与控制。

颜楚雄(1991-),男,博士,工程师,主要研究方向为导航、制导与控制。

李瑾(1996—),男,工程师,主要研究方向为导航、制导与控制。

Research of Guidance Method with Geocentric Angle Correction for Long Range Maneuvering Vehicle

Yudong WANG, Yue FANG, Ke GUO, Chuxiong YAN, Jin LI

Affiliations

Beijing Institute of Space Long March Vehicle,Beijing,100076

出版时间: 2024-12-25 doi: 10.7654/j.issn.2097-1974.20240607

文章导航

摘要

收起

针对远程机动飞行器长时间在大气层内飞行,传统的带落角约束的比例导引法在末制导段的期望落角受飞行航程地心角的影响较大,不能较好地适应远程机动飞行器落角约束的问题,故提出一种地心角修正的落角约束最优制导律。仿真结果表明该制导方法可以使远程机动飞行器的落角满足落角约束,并且具有一定的鲁棒性。

关键词

飞行器 / 比例导引 / 地心角 / 落角 / 最优制导律

Abstract

收起

The guidance method for long range maneuvering vehicle is focused. The traditional proportional guidance method with falling angle in the terminal guidance phase is greatly affected by the geocentric angle of the flight range and cannot well adapt to the falling angle constraint of long range maneuvering vehicle. An optimal guidance law with falling angle constraint using geocentric angle is proposed. The simulation results show that the guidance method can make the falling angle of long range maneuvering vehicle content the falling angle constraint, and has certain robustness.

Key words

vehicle / proportional guidance / geocentric angle / falling angle / optimal guidance law

引用本文

王煜东, 方岳, 郭珂, 颜楚雄, 李瑾. 远程机动飞行器地心角修正的制导方法. 导弹与航天运载技术, 2024 , 47 (6) : 42 -46 . DOI: 10.7654/j.issn.2097-1974.20240607

Yudong WANG, Yue FANG, Ke GUO, Chuxiong YAN, Jin LI. Research of Guidance Method with Geocentric Angle Correction for Long Range Maneuvering Vehicle[J]. Missiles and Space Vehicles, 2024 , 47 (6) : 42 -46 . DOI: 10.7654/j.issn.2097-1974.20240607

正文

收起

0 引言

收起

相比于常规再入机动飞行器而言, 远程机动飞行器需要在大气层中完成长时间、长距离飞行^[1]。传统的再入机动飞行器再入距离较近, 再入制导时采用带落角约束的比例导引法对该种情况适应性较好^[2],如胡正东等^[3]针对再入机动飞行器提出一种带落角约束的变结构最优导引律, 利用遗传算法离线优化制导参数, 尤浩等[4]基于非奇异快速终端滑模设计了一种带落角约束的优先时间滑膜制导律, 可以有效避免制导律设计过程中出现奇异问题; 陈琦等^[5]将传统的带落角约束最优制导律与加权函数相结合, 提出了一种间接Gauss伪普法的最优末制导; 张嘉文^[6]采用均值聚类神经网络的方法, 在飞行过程中不断调整聚类中心, 实现对落角约束制导律的优化。上述方法均对传统方法进行了一些修改, 具有一定的理论指导意义。但随着远程机动飞行器在大气层内机动飞行的距离不断增加, 飞行器航程所对应的地心角也逐渐增大,当期望落角与地心角之和大于${90}^{\circ }$时,在再入系下建立的带落角约束的最优制导律会发生奇异, 视线角与期望落角之间的反馈不理想, 这会导致末制导的效果有所降低。

本文针对上述问题提出了一种考虑地心角的落角约束最优制导律, 这种方法具有很强的工程可实现性, 能够较好地解决远程机动飞行器期望落角不理想的问题。通过仿真可以验证该方法能够满足远程机动飞行器的落角要求, 并且具有较强的鲁棒性。

1 建立飞行器数学模型

收起

以目标和飞行器质心为基准, 将空间运动分解为纵向平面和横向平面, 由于纵向平面内的运动对航程影响较大, 为简化分析, 仅分析纵向平面内的运动, 如图1所示。

其中,$M$为飞行器再入起始位置;$T$为目标点位置;${M}^{* }$;${T}^{* }$分别为虚拟飞行器、虚拟目标点位置;$v$为速度矢量;$\lambda ,{\lambda }^{* }$分别为纵向平面内起始位置与终点位置的视线角；$\gamma ,{\gamma }^{* }$分别为速度矢量在纵向平面内起始位置与终点位置的方位角,即速度倾角;$\eta$为速度矢量前置角。为简化分析, 认为地球模型是均质圆球模型,$\beta$为飞行器航程所对应的地心角,根据弦切角定理可知,再入系$x$轴与目标系${x}^{\prime }$之间的夹角也为$\beta$。因此可以得出虚拟飞行器${M}^{* }$与真实目标点$T$之间的相对位置关系:

(1)$\left\{\begin{array}{l}\lambda ={\lambda }^{* }+ \beta \\\gamma ={\gamma }^{* }+ \beta \\\eta =\lambda +\gamma ={\lambda }^{* }+ {\gamma }^{* }\end{array}\right.$

注意到再入系$o -{xyz}$与目标系$o -{x}^{\prime }{y}^{\prime }{z}^{\prime }$均为右手系,因此在图1中逆时针为正,飞行器从$M$位置运动到$T$位置为顺时针,其航程所对应的地心角为负值, 即$\beta < 0$。

为使远程机动飞行器终端落地角度达到一定值, 即希望落地点当地速度倾角满足约束条件, 需在目标系$o -{x}^{\prime }{y}^{\prime }{z}^{\prime }$中描述运动关系,令$r$表示飞行器与真实目标点之间的矢径, 则有:

(2)$\left\{\begin{array}{l}\dot{\mathbf{r}}= -\mathbf{v}\cos \eta \\\mathbf{r}{\dot{\lambda }}^{* }= \mathbf{v}\sin \eta \end{array}\right.$

将式(2)中的第2式对时间求导,并将式(1)带入, 可得纵向平面内的相对运动方程:

(3)${\ddot{\lambda }}^{* }= \left({\frac{\dot{\mathbf{v}}}{\mathbf{v}}- \frac{2\dot{\mathbf{r}}}{\mathbf{r}}}\right){\dot{\lambda }}^{* }- \frac{\dot{\mathbf{r}}}{\mathbf{r}}{\dot{\gamma }}^{* }$

2 最优制导律

收起

机动飞行器的最优制导律是终端有约束的最优导引规律。其终端约束包括机动飞行器的落地倾角和落地速度大小^[7]。本文研究的最优导引规律是对落地倾角有约束, 对地面固定目标进行攻击的机动飞行器导引规律。

终端约束条件取视线角与期望的当地速度倾角相等, 且视线角速度为零, 即:

(4)$\left\{\begin{array}{l}{\lambda }^{* }\left({t}_{\mathrm{f}}\right)= -{\gamma }_{\mathrm{f}}^{* }\\{\dot{\lambda }}^{* }\left({t}_{\mathrm{f}}\right)= 0 \end{array}\right.$

其中,${t}_{\mathrm{f}}$为飞行终端时间;${\gamma }_{\mathrm{f}}$为飞行终端速度倾角。令${x}_{1}= {\lambda }^{* }+ {\gamma }_{\mathrm{f}}^{* },{x}_{2}= {\lambda }^{* }$,则可得状态方程:

(5)$\left\{\begin{array}{l}{\dot{x}}_{1}= {x}_{2}\\{\dot{x}}_{2}= \left({\frac{\dot{\mathbf{v}}}{\mathbf{v}}- \frac{2\dot{\mathbf{r}}}{\mathbf{r}}}\right){x}_{2}- \frac{\dot{\mathbf{r}}}{\mathbf{r}}{\dot{y}}^{* }\end{array}\right.$

为得到显示解,假设$\dot{v}/v = 0$,并定义:

(6)${T}_{\mathrm{g}}= - r/\dot{r}$

状态方程简化为

(7)$\left\{\begin{array}{l}{\dot{x}}_{1}= {x}_{2}\\{\dot{x}}_{2}= \frac{2}{{T}_{\mathrm{g}}}{x}_{2}+ \frac{1}{{T}_{\mathrm{g}}}{\dot{\gamma }}^{* }\end{array}\right.$

令$\mathbf{A}= \left\lbrack \begin{matrix} 0 & 1 \\ 0 & 2/{T}_{\mathrm{g}}\end{matrix}\right\rbrack ,\mathbf{B}= \left\lbrack \begin{matrix} 0 \\ 1/{T}_{\mathrm{g}}\end{matrix}\right\rbrack ,\mathbf{x}= {\left({x}_{1},{x}_{2}\right)}^{\mathrm{T}},\mathbf{u}= {\dot{\gamma }}^{* }$, 则式(7)可改写为

(8)$\left\{\begin{array}{l}\dot{\mathbf{x}}= \mathbf{A}\mathbf{x}+ \mathbf{B}\mathbf{u}\\\mathbf{x}\left({t}_{\mathrm{f}}\right)= 0 \end{array}\right.$

由于远程机动飞行器终端对落角、落速均有一定约束, 在保证落角要求的同时使速度损失尽可能小, 进而可以有富余能量进行减速, 因此可以设计性能指标$J$:

(9)$ J ={\mathbf{x}}^{\mathrm{T}}\left({t}_{\mathrm{f}}\right)\mathbf{x}\left({t}_{\mathrm{f}}\right)+ \frac{1}{2}{\int }_{0}^{{t}_{\mathrm{f}}}{\left({\dot{\gamma }}^{* }\right)}^{2}\mathrm{\;d}t $

根据极大值原理,线性二次型性能指标最优控制为

(10)$\mathbf{u}= -{\mathbf{R}}^{-1}{\mathbf{B}}^{\mathrm{T}}\mathbf{{Px}}$

其中,$\mathbf{R},\mathbf{B}$和$\mathbf{P}$为线性二次型最优控制参数矩阵,当性能指标为式(9)时,$\mathbf{R}= 1,\mathbf{u}= {\dot{\gamma }}^{* }$,于是有:

(11)${\dot{\gamma }}^{* }= -{\mathbf{B}}^{\mathrm{T}}\mathbf{{Px}}$

通过逆解黎卡第矩阵,并忽略终端时刻${T}_{\mathrm{g}}$小量, 有:

(12)$\mathbf{P}= \left\lbrack \begin{matrix}\frac{4}{{T}_{\mathrm{g}}}& 2 \\ 2 & 4{T}_{\mathrm{g}}\end{matrix}\right\rbrack $

将$\mathbf{P}$带入式(11)展开,得到最优导引律:

(13)${\dot{\gamma }}^{* }= - 4{\dot{\lambda }}^{* }-2\frac{\left({\lambda }^{* }+ {\gamma }_{\mathrm{f}}^{* }\right)}{{T}_{\mathrm{g}}}$

3 地心角修正后的最优制导律

收起

注意到式(1)中${\gamma }^{* }$、${\lambda }^{* }$与$\gamma$、$\lambda$均相差一个地心角$\beta$,设$X\text{、}Y\text{、}Z$为再入系$o -{xyz}$下的飞行器与目标点之间的相对位置之差,${X}^{* }\text{、}{Y}^{* }\text{、}{Z}^{* }$为目标系$o -{x}^{\prime }{y}^{\prime }{z}^{\prime }$中的飞行器与目标点之间的相对位置之差, 根据图1中几何关系, 可得:

(14)$\sin {\lambda }^{* }= \frac{{Y}^{* }}{r}$

两边对时间$t$求导可得:

(15)${\lambda }^{* }\cos {\lambda }^{* }= \frac{{\dot{Y}}^{* }\mathbf{r}- {Y}^{* }\dot{\mathbf{r}}}{{\mathbf{r}}^{2}}$

根据余弦函数定义可以得到:

(16)${\lambda }^{* }= \frac{{\dot{Y}}^{* }\mathbf{r}- {Y}^{* }\dot{\mathbf{r}}}{{\mathbf{r}}^{2}\sqrt{1 -\frac{{\left({Y}^{* }\right)}^{2}}{{\mathbf{r}}^{2}}}}$

根据相对位置关系可得:

(17)${\lambda }^{* }= \arctan \left(\frac{{Y}^{* }}{{X}^{* }}\right)$

由于${r}^{2}= {\left({X}^{* }\right)}^{2}+ {\left({Y}^{* }\right)}^{2}+ {\left({Z}^{* }\right)}^{2}$,两边对$t$求导, 可得:

(18)${2r}\dot{r}= 2{X}^{* }{\dot{X}}^{* }+ 2{Y}^{* }{\dot{Y}}^{* }+ 2{Z}^{* }{\dot{Z}}^{* }$

两边同时除以$2{\mathbf{r}}^{2}$,并结合式(6),可以得到:

(19)${T}_{\mathrm{g}}= \frac{-{\mathbf{r}}^{2}}{{X}^{* }{\dot{X}}^{* }+ {Y}^{* }{\dot{Y}}^{* }+ {Z}^{* }{\dot{Z}}^{* }}$

根据坐标旋转矩阵, 注意到再入系与目标系之间绕$z$轴顺时针旋转地心角$\beta$,易得坐标转换矩阵为

(20)$\mathbf{M}\left(\beta \right)= \left\lbrack \begin{matrix}\cos \beta &\sin \beta & 0 \\- \sin \beta &\cos \beta & 0 \\ 0 & 0 & 1 \end{matrix}\right\rbrack $

注意$\beta$为负值,因此有:

(21)$\left\lbrack \begin{matrix}{X}^{* }\\{Y}^{* }\\{Z}^{* }\end{matrix}\right\rbrack =\left\lbrack \begin{matrix}\cos \beta &\sin \beta & 0 \\- \sin \beta &\cos \beta & 0 \\ 0 & 0 & 1 \end{matrix}\right\rbrack \left\lbrack \begin{matrix} X \\ Y \\ Z \end{matrix}\right\rbrack $

将式(1)、式(16)、式(17)、式(19)、式(21)带入式(13),得到带地心角修正的最优制导律。

(22)$\dot{\gamma }= - 4{f}_{1}\left({X, Y,\beta }\right)-2\frac{\left({f}_{2}\left( X, Y,\beta \right)+ {\gamma }_{\mathrm{f}}\right)}{{f}_{3}\left({X, Y,\beta }\right)} +\dot{\beta }$

其中,${f}_{1}= {\lambda }^{* },{f}_{2}= {\lambda }^{* },{f}_{3}= {T}_{\mathrm{g}}$。

式(22)中的地心角与地心角的导数可以根据远程机动飞行器的位置根据大地主题解算^[8]实时计算求得。

可以根据图2实现远程机动飞行器基于带地心角在线修正的落角约束最优制导策略。

4 轨迹仿真分析

收起

根据图2的制导策略进行轨迹仿真, 采用传统带落角约束的最优制导律与带地心角修正的落角约束最优制导律这两种方法进行对比, 分别考核小航程与大航程情况下的终端约束条件满足情况。期望落地倾角设置为$-{80}^{\circ }$。

4.1 典型工况对比

在小航程情况下,航程${400}\mathrm{\;{km}}$,对比曲线如图3至图4所示。

根据图3、图4可以看出, 在小航程情况下, 带地心角修正的最优制导律与传统方法相差不大, 因为此时航程对应的地心角较小, 但可以看出在设置同样期望落地倾角的情况下,带地心角修正的最优制导律在落角上比传统方法略大,采用地心角修正后的落角为${81.3}^{\circ }$,未修正的制导律落角为${77.8}^{\circ }$,大约相差${3.5}^{\circ }$,而${400}\mathrm{\;{km}}$航程所对应的地心角也大约为${3.5}^{\circ }$, 此时按传统方法可提高预设的期望落地倾角, 使传统方法与带地心角修正后的制导律具有同样效果。

大航程情况下,航程约${3000}\mathrm{\;{km}}$,对比曲线如图5、图6所示。

根据图6可以看出, 在大航程的情况下, 采用带地心角修正的制导律的飞行器落地倾角能够达到$-{80}^{\circ }$,但是采用传统方法,落角在$-{53}^{\circ }$左右,根据飞行轨迹可以看出,航程增大后,${3000}\mathrm{\;{km}}$航程对应地心角大约为${27}^{\circ }$,此时地心角与落角的绝对值之和大于${90}^{\circ }$后,传统方法出现奇异,速度倾角先减小到一定程度后即保持不变, 不能达到期望落角, 仿真效果不理想。但是采用地心角在线修正的制导律能够很好地满足约束条件。

4.2 蒙特卡洛打靶

为验证该制导律在扰动条件下的鲁棒性, 对远程机动飞行器在初始位置误差、初始速度误差、气动系数误差以及大气密度扰动情况下^[9-10],采用蒙特卡洛打靶法^[11-13]仿真1000条子样飞行轨迹,假设制导系统对误差的大小是不可知的, 仍按标称状态的参数计算制导指令, 其他仿真参数与4.1节相同。

1000条打靶仿真结果所绘制的飞行轨迹如图7所示, 落角情况如图8所示。

根据仿真可以看出, 在开始启动带地心角修正的最优制导律后, 飞行器的落地倾角能够很好地控制在期望落角上, 飞行器落地倾角的打靶仿真统计结果如表1所示。

仿真结果表明, 通过采用地心角修正后的落角约束最优制导律,飞行器落地倾角偏差不大于${10}^{\circ }$,这对远程机动飞行器制导而言精度是较高的, 故在加入近乎严苛的随机干扰下, 本文所提出的方法能够满足精度要求, 具有较强的鲁棒性。

5 结束语

收起

本文针对远程机动飞行器落地倾角不能满足理想的约束条件的问题, 提出了一种新型的、工程上易于实现的制导方法。通过将飞行器航程所对应地心角实时反馈在制导律中, 仿真结果表明, 该方法可有效针对远程机动飞行器落角不能达到期望落角的情况, 并且具有一定的鲁棒性。

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

[1]

魏毅寅, 李瑜, 段广仁. 远程滑翔导弹滑翔段制导算法研究[J]. 宇航学报, 2009,30(2): 481-485.

WEI

Yiyin

, LI

, DUAN

Guangren

. Research on guidance method of glide phase for long-range glide missile[J]. Journal of Astronautics, 2009,30(2): 481-485.

[2]

陈刚, 董龙雷, 闫桂荣, 等. 航天器再入制导方法研究现状与发展评述[J]. 飞行力学, 2008,26(1): 1-4.

CHEN

Gang

, DONG

Longlei

, YAN

Guirong

, et al. Recent development and review of spacecraft reentry guidance method[J]. Flight Dynamics, 2008,26(1): 1-4.

[3]

胡正东, 曹渊, 蔡洪. 带落角约束的再入机动弹头的变结构导引律[J]. 系统工程与电子技术, 2009,31(2): 393-398.

Zhengdong

, CAO

Yuan

, CAI

Hong

. Variable structure guidance law of reentry maneuvering warhead with terminal angular constraint[J]. System Engineering and Electronics, 2009,31(2): 393-398.

[4]

尤浩, 赵久奋, 李朋, 等. 带落角约束的非奇异快速终端二阶滑膜制导律[J]. 弹箭与制导学报, 2019,39(6): 155-159.

YOU

Hao

, ZHAO

Jiufen

, LI

Peng

, et al. Nonsingular fast terminal second-order sliding mode guidance law with impact angle constraints[J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2019,39(6): 155-159.

[5]

陈琦, 杨靖, 王中原, 等. 带有双曲正切加权函数的落角约束最优制导律[J]. 哈尔滨工业大学学报, 2020,52(4): 92-99.

CHEN

, YANG

Jing

, WANG

Zhongyuan

, et al. Impact angle constrained optimal guidance law based on hyperbolic tangent weighting functions[J]. Journal of Harbin Institute of Technology, 2020,52(4): 92-99.

[6]

张嘉文. 带落角约束的均值聚类神经网络滑模制导律研究[J]. 电光与控制, 2021,28(3): 47-55.

ZHANG

Jiawen

. Sliding mode guidance law of means clustering neural network with falling angle constraint[J]. Electronics Optics & Control, 2021,28(3): 47-55.

[7]

赵汉元. 飞行器再入动力学和制导[M]. 长沙: 国防科技大学出版社, 1997.

ZHAO

Hanyuan

. Reentry dynamics and guidance of aerial vehicles[M]. Changsha: National Defense Science and Technology University Press, 1997.

[8]

汪绍航, 边少锋, 金立新, 等. 椭球大地测量常用幂级数的第三扁率展开[J]. 测绘科学技术学报, 2021(6): 571-578.

WANG

Shaohang

, BIAN

Shaofeng

, JIN

Lixin

, et al. The third flattening expansion of power series commonly used in ellipsoid geodesy[J]. Journal of Geomatics Science and Technology, 2021(6): 571-578.

[9]

童伟, 贾平会, 张宁, 等. 远程机动飞行器能量管控制导方法[J]. 导弹与航天运载技术, 2018(4): 88-92.

TONG

Wei

, JIA

Pinghui

, ZHANG

Ning

, et al. Research of energy management and control strategy for long range maneuvering vehicle[J]. Missile and Space Vehicles, 2018(4): 88-92.

[10]

严东升, 贾平会, 郭珂, 等. 基于小样本试验数据再入偏差折合及评估方法研究[J]. 导弹与航天运载技术, 2021(1): 45-49.

YAN

Dongsheng

, JIA

Pinghui

, GUO

, et al. Research on reentry error reduction and evaluation method based on small sample test data[J]. Missile and Space Vehicles, 2021(1): 45-49.

[11]

钱杏芳, 林瑞雄. 导弹飞行力学[M]. 北京: 北京理工大学出版社, 2000.

QIAN

Xingfang

, LIN

Ruixiong

. Missile flight dynamics[M]. Beijing: Beijing Institute of Technology Press, 2000.

[12]

杨军, 杨晨, 段朝阳, 等. 现代导弹制导控制系统设计[M]. 北京: 航空工业出版社, 2005.

YANG

Jun

, YANG

Chen

, DUAN

Chaoyang

, et al. Design of modern missile guidance and control system[M]. Beijing: Aviation Industry Press, 2005.

[13]

胡寿松. 自动控制原理[M]. 北京: 国防工业出版社, 1994.

Shousong

. Principle of automatic control[M]. Beijing: National Defence Industry Press, 1994.

2024年第47卷第6期

PDF下载

418

174

引用本文

BibTeX

文章信息

doi: 10.7654/j.issn.2097-1974.20240607

接收时间：2022-07-01
首发时间：2025-07-04
出版时间：2024-12-25

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2022-07-01
修回日期：2023-06-28

基金

作者信息

北京航天长征飞行器研究所,北京,100076

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/ddyht/CN/10.7654/j.issn.2097-1974.20240607

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT