CN103921957A - 一种探月飞船跳跃式再入的跃起点能量管理方法 - Google Patents
一种探月飞船跳跃式再入的跃起点能量管理方法 Download PDFInfo
- Publication number
- CN103921957A CN103921957A CN201410019094.7A CN201410019094A CN103921957A CN 103921957 A CN103921957 A CN 103921957A CN 201410019094 A CN201410019094 A CN 201410019094A CN 103921957 A CN103921957 A CN 103921957A
- Authority
- CN
- China
- Prior art keywords
- airship
- lunar exploration
- energy
- exploration airship
- std
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
Abstract
本发明公开了一种探月飞船跳跃式再入的跃起点能量管理方法,首先采用最大倾侧角再入,以利于探月飞船捕获;然后控制探月飞船的能量耗散速率与能量耗散量,并根据在线参数辨识结果实时预估跃起点处的能量值,确定能量管理段的退出时机;最后按预设的常值倾侧角飞行,使飞船在预定能量值处跃起。本发明的方法能对跃起点处的能量实施精细管理,从而提高探月飞船大气捕获的概率,为后续航程预留合适的能量。
Description
技术领域
本发明涉及制导控制技术领域,特别是一种探月飞船跳跃式再入的跃起点能量管理方法。
背景技术
探月飞船返回再入地球时,到达大气层边界的初始状态与近地航天器返回有很大不同,最明显的特征是探月飞船的再入速度更大,达到11km/s,而近地航天器的再入速度一般约为7.8km/s。为满足再入过程中的热流密度、过载等约束条件,同时增大再入航程,提高再入的安全性与灵活性,探月飞船多采用跳跃式再入方式,如美国的“阿波罗”飞船、成员探测飞行器(CEV)等。
与直接式再入相比,跳跃式再入的过程更加复杂,制导控制的难度更大。根据跳跃式再入过程中高度的变化,可以将整个过程分为三段:一次再入段、大气层外的开普勒段和二次再入段(附图1)。研究表明,跃起点处的能量对再入弹道特性有着决定性的影响。首先,由于飞船再入速度接近第二宇宙速度,因此在一次再入段必须耗散过剩的能量,确保飞船的大气捕获;其次,由于飞船的升阻比较小,控制能力较弱,跃起后的能量过剩或不足都可能导致后续航程偏差无法修正,致使落点偏差过大。因此,对探月飞船跃起点的能量进行精细管理是十分必要的。
在“阿波罗”飞船的再入制导方法中,在跃起点前采用常数阻力段实施能量管理,同时解析预测跃起点应保留的能量值,直到将能量耗散到预留值附近。由于预测过程中采用了解析公式,因此难以实现能量的精细管理,方法的航程适用范围有限。喷气推进实验室(JPL)针对CEV飞行器的任务,在“阿波罗”制导方法的基础上用数值预测代替解析预测,从而能更准确的得到跃起点处需要的能量值,提高了制导方法的精确性和航程适用性。但数值预测对船载计算机的性能要求较高,目前还未在在工程实践中得到应用。
发明内容
本发明所要解决的技术问题是,针对现有技术不足,提供一种探月飞船跳跃式再入的跃起点能量管理方法,对跃起点处的能量实施精细管理,从而提高探月飞船大气捕获的概率,为后续航程预留合适的能量。
为解决上述技术问题,本发明所采用的技术方案是:一种探月飞船跳跃式再入的跃起点能量管理方法,该方法为:
1)探月飞船再入大气层前,令一次再入段的倾侧角为常数,设计满足航程约束、过程约束和终端约束的标称轨道,并存储以下参数:一次再入段探月飞船的速度值及该速度值对应的阻力加速度值、升力加速度值及标称倾侧角值νstd,跃起点的标称能量值Estd;设定探月飞船的气动过载值nset;
2)从探月飞船的再入点起,控制探月飞船按允许的最大倾侧角νmax飞行,判断探月飞船的气动过载值n是否大于预先设定的气动过载值nset;若是,则设定期望的探月飞船能量耗散速率CE、阻力加速度比例因子初值和升阻比比例因子初值进入步骤3);否则,进入步骤11);其中,CE为常数;
3)测量得到探月飞船的视加速度根据惯性导航原理得到探月飞船当前时刻的飞行速度地心距水平高度和高度变化率由计算得到探月飞船阻力加速度的测量值和升力加速度的测量值从而得到升阻比的测量值
4)根据上述步骤1)中存储的一次再入段探月飞船的速度值及该速度值对应的阻力加速度值、升力加速度值,得到探月飞船当前时刻的飞行速度对应的标称阻力加速度值和标称升力加速度值从而得到升阻比的标称值
5)对再入过程中的大气密度偏差和气动系数偏差进行辨识,计算阻力加速度比例因子和升阻比比例因子由低通滤波器获得阻力加速度比例因子的估计值和升阻比比例因子的估计值其中, 为阻力加速度的测量值;
6)利用CE求探月飞船的标准纵向升阻比
其中,β为大气密度系数,g为地球引力加速度;m为探月飞船质量;
7)利用下式计算期望的纵向升阻比
其中,KE为能量偏差增益系数,为能量变化率偏差增益系数;
则飞船制导指令倾侧角的余弦为:
8)利用下式计算探月飞船跃起点处能量的期望值EExp:
根据下式计算能量管理段退出时刻探月飞船的速度VExit:
9)若探月飞船当前时刻的飞行速度则重复步骤3)~10),直到退出能量管理段;
11)探月飞船进入标称倾侧角值νstd飞行段,直至探月飞船高度变化率
与现有技术相比,本发明所具有的有益效果为:本发明在线辨识大气密度偏差和飞船气动系数偏差,根据辨识结果预测并控制跃起点处的能量值,从而提高探月飞船跳跃式再入的大气捕获概率,为后续航程预留合适的能量,提高制导方法的精确性和航程适用性;本发明通过常值能量耗散速率段管理跃起点处的能量,与常值阻力能量管理方法相比,对跃起点能量的估计和控制更准确,能实现能量的精细管理;与数值预测管理方法相比,大大减小了能量预测的计算量,更利于工程实现。
附图说明
图1为跳跃式再入过程中跃起点能量管理分段示意图;
图2为本发明方法流程图;
图3为采用本发明的5000km航程探月飞船开伞点位置精度。
具体实施方式
本发明的基本思路是:首先采用最大倾侧角再入,以利于探月飞船捕获;然后控制探月飞船的能量耗散速率与能量耗散量,并根据在线参数辨识结果实时预估跃起点处的能量值,确定能量管理段的退出时机;最后按预设的常值倾侧角飞行,使飞船在预定能量值处跃起。
假定某探月飞船通过质心偏置配置方式在再入过程中获得升力,升阻比为0.3~0.4,制导系统的控制量为倾侧角ν。再入过程从高度120km处开始,此时飞船速度为11km/s,当地速度倾角为-5.8。飞船再入过程的运动方程可表示为:
其中,λ和φ分别表示飞船的经度和纬度,θ和ψ分别表示飞船速度倾角和航迹偏航角,r表示地心距,g为地球引力加速度,ω为地球自转角速度,m为飞船质量。L和D分别为飞船再入过程中的升力和阻力,计算公式为
其中,CL、CD分别为飞船的升力系数和阻力系数,S为飞船的参考面积。飞船采用惯性导航系统,可以测量到飞船体坐标系中的视加速度和飞船相对于惯性坐标系的姿态角速度根据惯性导航的基本原理,可以得到飞船的速度、地心距、高度、高度变化率等信息。
如图2所示,针对本实施例中的飞船,本发明的具体实施步骤如下:
S1.在飞船再入前,令一次再入段的倾侧角为常值,设计满足航程约束、过程约束和终端约束的标称轨道(见赵汉元.飞行器再入动力学和控制[M].长沙:国防科学技术大学出版社,1997)。在船载计算机上存储如下数值:一次再入段的阻力加速度-速度(D-V)和升力加速度-速度(L-V)的变化规律,标称倾侧角值νstd,跃起点的标称能量值Estd。
S2.从再入点120km高度起,控制飞船按倾侧角ν=180°飞行,保证飞船在一次再入段的捕获概率。当飞船的气动过载n大于设定值nset=1.5g时,进入能量管理段。
S3.进入能量管理段后,首先设定以下值:期望的飞船常值能量耗散速率CE设定为标称轨道跃起点处的能量耗散速率值阻力加速度比例因子初值 升阻比比例因子初值
然后循环执行下述步骤,直至满足终止条件:
(1)利用飞船上安装的加速度计测量得到视加速度根据惯性导航原理可以得到飞船当前时刻的飞行速度地心距高度和高度变化率将投影到飞船相对于地球的速度矢量的方向,得到阻力加速度的测量值计算公式为:
利用公式(3)计算升力加速度的测量值
由此得到升阻比的测量值
(2)根据船载计算机上存储的阻力加速度-速度(D-V)和升力加速度-速度(L-V)的变化规律,以飞船速度V为自变量,分别以阻力加速度D和升力加速度L为因变量,通过线性插值的方法获得当前速度对应的标称阻力加速度值和标称升力加速度值从而得到升阻比的标称值
(3)计算阻力加速度比例因子和升阻比比例因子利用一阶低通滤波器,得到阻力加速度比例因子的估计值
其中,初值取为利用一阶低通滤波器,得到飞船升阻比比例因子的估计值
其中,初值取为KgD和KgLD为滤波增益,取为0.4。
(4)若设地球表面为零势能面,则根据能量的计算公式可求得能量的变化率为
由于在一次再入段速度倾角很小,可近似令sinθ≈0,cosθ≈1。同时,由于势能的变化为小量,可忽略不计,由此得到简化的能量耗散速率表达式
若设期望的飞船能量耗散速率为常值CE,则有
CE=-VD (9)
继续对能量耗散速率求导,忽略高阶小量可得到能量耗散加速度为
由于能量耗散速率为常值,所以能量耗散加速度由此可得高度变化率为
继续对上式求导,并考虑到可得高度的变化加速度为
由飞行器运动方程(1)的第一式和第六式,又可得
综合式(12)和式(13),可得当前时刻能使飞船的能量耗散速率为常值的标准纵向升阻比为
实际飞行过程中,为实现常值能量耗散速率,可通过反馈能量耗散率和能量耗散加速度,得到期望的纵向升阻比
KE为能量偏差增益系数,为能量变化率偏差增益系数;本实施例中,KE=1.30;
由此得到控制量倾侧角为
(5)根据当前时刻的参数辨识结果,由式(17)计算跃起点处能量的期望值
(6)假设飞船运动到跃起点处的时间为Tall,则Tall分别由能量耗散率控制段时间T1和常值倾侧角段时间T2组成,即
Tall=T1+T2 (18)
由式(13)知,当倾侧角为常值,速度变化不大时,高度加速度近似为常值。
当进入常值倾侧角段的高度变化率确定时,则常值倾侧角段的时间T2为
结合式(9)和式(10),可得到能量管理段末端处高度变化率与速度VExit的关系
根据飞船当前时刻的速度和高度并忽略高度变化对能量产生的影响,可得到能量耗散速率控制段的时间T1为
结合式(13)、(18)~(21),可以得到跃起点处能量的期望值EExp与速度VExit满足如下时间方程
迭代求解上式,即可得到能量管理段退出时刻飞船的速度VExit。若飞船当前时刻的速度小于或等于VExit,则退出能量管理段,进入常值倾侧角飞行段。
S4.飞船按标称轨道设计得到的常值倾侧角νstd飞行,直至高度变化率也即实现跃起为止。
为进一步说明本发明对探月飞船跳跃式再入制导的作用,此处给出一个仿真算例。仿真条件设置为:飞船到达再入点的高度为120km,再入速度为11km/s,再入角为-5.9°,气动参数参照美国成员探测飞行器(CEV),标准航程为5000km,制导算法采用标称轨道制导算法。
仿真过程中,考虑的各种误差的类型、分布及大小如表1所示。
表1误差参数类型及分布
随机抽取1000组误差进行MonteCarlo数值打靶仿真,是否采用能量管理的打靶参数统计如表2所示。从表中数据对比可以知道,进行能量管理能有效提高跃起点能量的控制精度和开伞点的位置精度,并能在一定程度降低再入最大过载。附图3给出了开伞点位置分布图。
表2Monte-Carlo数值仿真打靶统计结果
本发明的名词概念附注:
(1)气动过载:飞行器所受气动力与飞行器重量的比值。
(2)倾侧角:绕速度轴旋转的一个欧拉角,反映了升力对铅垂面的倾侧。
(3)地心距:飞行器质心到地心的距离。
(4)参考面积:飞行器的气动有效面积,一般由气动计算部门给出。
升力系数、阻力系数(气动系数):飞行器所受气动力与气流动压和参考面积乘积之比。
Claims (3)
1.一种探月飞船跳跃式再入的跃起点能量管理方法,其特征在于,该方法为:
1)探月飞船再入大气层前,令一次再入段的倾侧角为常数,设计满足航程约束、过程约束和终端约束的标称轨道,并存储以下参数:一次再入段探月飞船的速度值及该速度值对应的阻力加速度值、升力加速度值及标称倾侧角值νstd,跃起点的标称能量值Estd;预先设定探月飞船的气动过载值nset;
2)从探月飞船的再入点起,控制探月飞船按允许的最大倾侧角νmax飞行,判断探月飞船的气动过载值n是否大于上述预先设定的探月飞船的气动过载值nset;若是,则设定期望的探月飞船能量耗散速率CE、阻力加速度比例因子初值和升阻比比例因子初值进入步骤3);否则,进入步骤11);其中,CE为常数;
3)测量得到探月飞船的视加速度根据惯性导航原理得到探月飞船当前时刻的飞行速度地心距水平高度和高度变化率由计算得到探月飞船阻力加速度的测量值和升力加速度的测量值从而得到升阻比的测量值
4)根据上述步骤1)中存储的一次再入段探月飞船的速度值及该速度值对应的阻力加速度值、升力加速度值,得到探月飞船当前时刻的飞行速度对应的标称阻力加速度值和标称升力加速度值从而得到升阻比的标称值
5)对再入过程中的大气密度偏差和气动系数偏差进行辨识,计算阻力加速度比例因子和升阻比比例因子由低通滤波器获得阻力加速度比例因子的估计值和升阻比比例因子的估计值其中, 为阻力加速度的测量值;
6)利用CE求探月飞船的标准纵向升阻比
其中,β为大气密度系数,g为地球引力加速度;m为探月飞船质量;
7)利用下式计算期望的纵向升阻比
其中,KE为能量偏差增益系数,为能量变化率偏差增益系数;
则飞船制导指令倾侧角的余弦为:
8)利用下式计算探月飞船跃起点处能量的期望值EExp:
根据下式计算能量管理段退出时刻探月飞船的速度VExit:
9)若探月飞船当前时刻的飞行速度则重复步骤3)~10),直到≤VExit,退出能量管理段;
11)探月飞船进入标称倾侧角值νstd飞行段,直至探月飞船高度变化率
2.根据权利要求1所述的探月飞船跳跃式再入的跃起点能量管理方法,其特征在于,所述步骤5)中,利用一阶低通滤波器计算阻力加速度比例因子的估计值和升阻比比例因子的估计值
其中,KgD和KgLD为一阶低通滤波器的滤波增益。
3.根据权利要求2所述的探月飞船跳跃式再入的跃起点能量管理方法,其特征在于,KgD和KgLD均取为0.4。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410019094.7A CN103921957B (zh) | 2014-01-16 | 2014-01-16 | 一种探月飞船跳跃式再入的跃起点能量管理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410019094.7A CN103921957B (zh) | 2014-01-16 | 2014-01-16 | 一种探月飞船跳跃式再入的跃起点能量管理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103921957A true CN103921957A (zh) | 2014-07-16 |
CN103921957B CN103921957B (zh) | 2015-12-30 |
Family
ID=51140405
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410019094.7A Active CN103921957B (zh) | 2014-01-16 | 2014-01-16 | 一种探月飞船跳跃式再入的跃起点能量管理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103921957B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104597756A (zh) * | 2014-12-18 | 2015-05-06 | 北京控制工程研究所 | 一种跳跃式再入二次再入段航程预估方法 |
CN104787361A (zh) * | 2015-04-02 | 2015-07-22 | 北京航天自动控制研究所 | 升力式飞行器再入制导的阻力加速度变化率的确定方法 |
CN106292332A (zh) * | 2016-09-08 | 2017-01-04 | 上海卫星工程研究所 | 基于阻力测量的可分离式探测器气动捕获制导方法 |
CN107861517A (zh) * | 2017-11-01 | 2018-03-30 | 北京航空航天大学 | 基于线性伪谱的跳跃式再入飞行器在线弹道规划制导方法 |
CN110733670A (zh) * | 2019-11-05 | 2020-01-31 | 中国人民解放军国防科技大学 | 一种短航程低过载的再入轨迹设计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7219859B2 (en) * | 2005-04-01 | 2007-05-22 | Spacehab, Incorporated | Multipurpose modular spacecraft |
US20100250031A1 (en) * | 2009-03-27 | 2010-09-30 | The Charles Stark Draper Laboratory, Inc. | Propulsive guidance for atmospheric skip entry trajectories |
CN102139766A (zh) * | 2010-10-28 | 2011-08-03 | 中国科学院力学研究所 | 一种亚轨道飞行器再入飞行的设计攻角获取方法 |
-
2014
- 2014-01-16 CN CN201410019094.7A patent/CN103921957B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7219859B2 (en) * | 2005-04-01 | 2007-05-22 | Spacehab, Incorporated | Multipurpose modular spacecraft |
US20100250031A1 (en) * | 2009-03-27 | 2010-09-30 | The Charles Stark Draper Laboratory, Inc. | Propulsive guidance for atmospheric skip entry trajectories |
CN102139766A (zh) * | 2010-10-28 | 2011-08-03 | 中国科学院力学研究所 | 一种亚轨道飞行器再入飞行的设计攻角获取方法 |
Non-Patent Citations (2)
Title |
---|
李惠峰: "探月飞船预测-校正再入制导律设计", 《空间控制技术与应用》, vol. 35, no. 1, 15 February 2009 (2009-02-15), pages 19 - 24 * |
郭敏文: "探月返回器跳跃式再入过载抑制算法研究", 《宇航学报》, vol. 34, no. 12, 15 December 2013 (2013-12-15), pages 1569 - 1577 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104597756A (zh) * | 2014-12-18 | 2015-05-06 | 北京控制工程研究所 | 一种跳跃式再入二次再入段航程预估方法 |
CN104787361A (zh) * | 2015-04-02 | 2015-07-22 | 北京航天自动控制研究所 | 升力式飞行器再入制导的阻力加速度变化率的确定方法 |
CN104787361B (zh) * | 2015-04-02 | 2015-11-04 | 北京航天自动控制研究所 | 升力式飞行器再入制导的阻力加速度变化率的确定方法 |
CN106292332A (zh) * | 2016-09-08 | 2017-01-04 | 上海卫星工程研究所 | 基于阻力测量的可分离式探测器气动捕获制导方法 |
CN107861517A (zh) * | 2017-11-01 | 2018-03-30 | 北京航空航天大学 | 基于线性伪谱的跳跃式再入飞行器在线弹道规划制导方法 |
CN110733670A (zh) * | 2019-11-05 | 2020-01-31 | 中国人民解放军国防科技大学 | 一种短航程低过载的再入轨迹设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103921957B (zh) | 2015-12-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103708045B (zh) | 一种探月飞船跳跃式再入的在线参数辨识方法 | |
Benito et al. | Reachable and controllable sets for planetary entry and landing | |
CN104035335A (zh) | 基于高精度纵、横程解析预测方法的平稳滑翔再入制导律 | |
US20160238481A1 (en) | Aerodynamic Modeling Using Flight Data | |
Liesk et al. | Design and experimental validation of a nonlinear low-level controller for an unmanned fin-less airship | |
CN103921957B (zh) | 一种探月飞船跳跃式再入的跃起点能量管理方法 | |
CN110015446B (zh) | 一种半解析的火星进入制导方法 | |
CN105182985A (zh) | 高超声速飞行器俯冲段全量一体化制导控制方法 | |
Wenz et al. | Moving horizon estimation of air data parameters for UAVs | |
Wolf et al. | Performance trades for Mars pinpoint landing | |
CN104536457A (zh) | 基于小型无人机导航的滑模控制方法 | |
CN103453907B (zh) | 基于分层大气模型的行星进入段导航滤波方法 | |
CN103542853B (zh) | 一种估计加速度计漂移的绝对导航滤波方法 | |
Putnam et al. | PredGuid entry guidance for Orion return from low Earth orbit | |
Wenz et al. | Estimation of wind velocities and aerodynamic coefficients for UAVs using standard autopilot sensors and a moving horizon estimator | |
CN106197432A (zh) | 一种基于FastSLAM算法的无人机着陆方法 | |
US11029709B1 (en) | Adaptive wind estimation, trajectory generation, and flight control for aerial systems using motion data | |
Borobia et al. | Flight-path reconstruction and flight test of four-line power kites | |
CN104156595B (zh) | 飞行器的飞行轨迹指令的确定方法和装置 | |
CN107804487A (zh) | 一种基于自适应偏差控制的跳跃式再入返回落点预报方法 | |
Holanda et al. | Trajectory control system for the LAICAnSat-3 mission | |
CN115542746A (zh) | 高超声速飞行器的能量管控再入制导方法及装置 | |
Whitmore et al. | Development of a closed-loop strap down attitude system for an ultrahigh altitude flight experiment | |
Ahsan et al. | Verification of equilibrium point stability for linearization of an aircraft model | |
CN104634183B (zh) | 一种基于升阻比实时估计的自适应制导方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |