CN103863579B - 一种深空探测返回过程的预测校正制导方法 - Google Patents
一种深空探测返回过程的预测校正制导方法 Download PDFInfo
- Publication number
- CN103863579B CN103863579B CN201410125842.XA CN201410125842A CN103863579B CN 103863579 B CN103863579 B CN 103863579B CN 201410125842 A CN201410125842 A CN 201410125842A CN 103863579 B CN103863579 B CN 103863579B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- reentry
- section
- roll angle
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 230000008569 process Effects 0.000 claims description 21
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000001174 ascending effect Effects 0.000 claims description 3
- 230000009191 jumping Effects 0.000 claims description 2
- 239000002775 capsule Substances 0.000 abstract 5
- 238000004422 calculation algorithm Methods 0.000 description 9
- 238000004088 simulation Methods 0.000 description 8
- 230000001629 suppression Effects 0.000 description 8
- 230000005484 gravity Effects 0.000 description 3
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000004927 fusion Effects 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种深空探测返回过程的预测校正制导方法,步骤为:(1)确定所述深空探测返回过程轨迹预测采用的倾侧角剖面;(2)迭代计算当前点倾侧角的幅值|σ0|;(3)将返回器的跳跃式再入飞行轨迹分为以下几个阶段:初始再入段、首次再入下降段、首次再入上升段、开普勒段和最终再入段,并且确定返回器的最终制导律为:当返回器处于初始再入段、首次再入下降段、首次再入上升段或最终再入段时,使用步骤(2)确定的所述当前点倾侧角σ0的幅值|σ0|进行制导;当返回器处于开普勒段时,使用|σ0|=180°进行制导,从而完成所述深空探测返回过程的预测校正制导。本发明针对深空探测返回过程,在轨迹的开普勒段设计制导律实现过载抑制,以保证返回器安全、精确着陆。
Description
技术领域
本发明涉及一种深空探测返回过程的预测校正制导方法,属于航天器深空探测返回过程的制导领域。
背景技术
深空探测返回器再入地球大气层边界,距地面高度约120km处速度约为11km/s,初始动能约为近地轨道航天器再入时的2倍,这给再入制导带来一系列新的困难和挑战,如航程偏差对初始状态偏差更加敏感,再入过程热流密度峰值增加,过载峰值增加等等,这些都将严重影响返回器着陆的精度、安全性和可靠性。
返回器一般采用小升力体构型,靠重心的横移来产生攻角,而横移不可能很大,所以攻角不能很大,则升阻比较小,一般在0.3~0.5之间。返回器再入为配平攻角飞行,即忽略了攻角和侧滑角的控制,仅仅通过调整倾侧角来改变升力在纵向剖面的分量,从而控制飞行轨迹。对于深空探测返回,升阻比很小会严重限制飞行器的机动能力,影响着陆精度。要实现大航程、高精度、强约束的再入飞行,对小升阻比深空探测返回器的制导控制设计提出了全新的要求。
再入制导方法主要可分为两类:标准轨道法和预测校正法。预测校正制导具有较高的落点精度,且对再入初始条件不敏感。随着深空探测飞行器返回时航程、精度、可靠性等要求的提高,近年来该算法已成为国内外的研究热点。然而在预测校正制导律设计的过程中,过载约束问题尤其是针对高速返回时跳跃式再入轨迹仍未能得到很好地解决。在上个世纪60年代,研究再入返回问题的相关文献中提出的过载抑制方法有:减小初始再入角和利用气动升力进行控制。然而对于有限升力的返回器以第二宇宙速度再入时,减小初始再入角可能导致返回器因为首次再入不够深入而飞出大气层,可见亟需新的过载抑制算法的提出。
因此针对深空探测返回设计满足过载约束、具有高落点精度的再入制导律变得十分重要。
发明内容
本发明的技术解决问题是:克服现有技术的不足,针对深空探测返回器再入初始速度大,且需实现大航程、高精度、安全的跳跃式再入制导问题,提供了一种深空探测返回过程的预测校正制导方法,在轨迹的开普勒段设计制导律实现过载抑制,以保证返回器安全、精确着陆。
本发明的技术解决方案是:
一种深空探测返回过程的预测校正制导方法,步骤如下:
(1)确定所述深空探测返回过程轨迹预测采用的倾侧角剖面;
倾侧角幅值剖面的数学表达式为
其中为总的目标飞行航程,Stogo为当前点的剩余航程,σf为70°,Sthres=2000km,σ为倾侧角,σ0为当前点倾侧角;
(2)迭代计算当前点倾侧角的幅值|σ0|,具体为:
(2.1)确定考虑地球自转影响的三自由度无量纲再入动力学方程为:
其中,ρ为大气密度,r为地心距,r=h+Re,h为跳跃高度,hf为开伞高度,V为返回器相对地球速度;地球平均半径Re为6378.135km,海平面重力加速度g0为9.81m/s2;θ为经度,φ为纬度,γ为飞行路径角,即飞行速度矢量与当地水平方向的夹角;ψ为速度方位角,即飞行速度矢量与当地正北方向的夹角,顺时针为正;σ为倾侧角,从返回器内部来看,右侧为正;D,L分别为无量纲化的气动阻力和气动升力;CD,CL分别为阻力系数和升力系数;返回器最大横截面积为Sref,质量为m;参数上方的“.”代表一阶导数,S为飞行航程;Ω为地球自转角速度;r和V为无量纲值;
(2.2)采用步骤(1)确定的倾侧角剖面,数值积分步骤(2.1)确定的再入动力学方程,直到开伞点高度hf=10km,得到飞行航程S的积分终值:预测航程Sf,且Sf是关于当前点倾侧角σ0的函数;
(2.3)将预测航程Sf与当前点的剩余航程Stogo进行比较,若预测航程Sf的偏差Smiss=Stogo-Sf(σ0)满足预设的着陆精度要求,即以|σ0|作为当前点倾侧角幅值;若航程偏差Smiss不满足预设的着陆精度要求,则需要更新当前倾侧角幅值|σ0|;
(3)将返回器的跳跃式再入飞行轨迹分为以下几个阶段:初始再入段、首次再入下降段、首次再入上升段、开普勒段和最终再入段,并且确定返回器的最终制导律为:当返回器处于初始再入段、首次再入下降段、首次再入上升段或最终再入段时,使用步骤(2)确定的所述当前点倾侧角σ0的幅值|σ0|进行制导;当返回器处于开普勒段时,使用|σ0|=180°进行制导,从而完成所述深空探测返回过程的预测校正制导。
所述步骤(2.3)中更新当前倾侧角幅值|σ0|具体为:
更新当前倾侧角幅值|σ0|即对方程Smiss=Stogo-Sf(σ0)=f(cosσ0)=0求根;
采用割线法迭代求解方程Smiss=Stogo-Sf(σ0)=f(cosσ0)=0,迭代表达式为
其中,xi为第i步迭代的值且f(xi)为第i步所求落点偏差Smiss,即预测落点与实际落点的偏差。
所述步骤(3)中初始再入段、首次再入下降段、首次再入上升段、开普勒段和最终再入段具体为:
当过载n为0.05<n<0.2时,返回器进入初始再入段;
当过载n>0.2且飞行路径角γ<0时,返回器进入首次再入下降段;
当过载n>0.2且γ>0时,返回器进入首次再入上升段;
当过载n再次小于0.2时,返回器进入开普勒段;
当过载n再次大于0.2时,返回器进入最终再入段;
其中,
本发明与现有技术相比的有益效果是:
本发明中的预测校正制导算法通过增加开普勒的倾侧角幅值来减小开普勒段的飞行航程,从而增加了最终再入段的剩余航程,可使剩余段飞行更加平缓,过载后峰值得到抑制,满足约束条件,同时没有增加预测校正制导算法的复杂性,具有较好的着陆精度,提高了飞船着陆的安全性。
附图说明
图1为本发明适用的跳跃式再入轨迹预测倾侧角剖面;
图2为本发明适用的跳跃式再入轨迹示意图;
图3为本发明适用的跳跃式再入轨迹开普勒段示意图;
图4为本发明标称参数下不同情况时倾侧角指令比较;
图5为本发明标称参数下不同情况时最终再入段总航程比较;
图6为本发明标称参数下不同情况时最终再入段高度-速度剖面比较;
图7为本发明标称参数下不同情况时过载随时间变化曲线比较;
图8为本发明标称参数下不同情况时过载值的历经时间曲线比较;
图9(a)、图9(b)为本发明偏差参数下不同情况时承受过载的情况比较;
图10(a)、图10(b)为本发明偏差参数下不同情况时着陆精度情况比较;
图11为本发明方法流程图。
具体实施方式
如图11所示,本发明提供了一种深空探测返回过程的预测校正制导方法,设计开普勒段过载抑制环节,形成融合制导方案,包括如下步骤:
(1)确定所述深空探测返回过程轨迹预测采用的倾侧角剖面;
确定所述轨迹预测倾侧角剖面包括倾侧角幅值剖面的确定和倾侧角符号确定两部分内容。倾侧角幅值|σ|剖面为“线性+常值”的形式,如图1所示,其数学表达式为
式(1)中为总的目标飞行航程,Stogo为当前点的剩余航程,σ为倾侧角,σf为常值,可取为70°,Sthres取为2000km,σ0为当前点倾侧角,需通过迭代算法在每个制导周期进行更新,迭代初值可选为上一制导周期所确定的值,在初始再入点可选为0。
另外,倾侧角符号由横程制导逻辑决定,这里所用到的横程制导逻辑类似阿波罗再入制导方案,详见1972年由Graves,C.A.和Harpold,J.C所写NASA报告ApolloExperienceReport-MissionPlanningforApolloEntry。
(2)迭代计算当前点倾侧角的幅值|σ0|,具体为:
(2.1)确定考虑地球自转影响的三自由度无量纲再入动力学方程为:
其中,ρ为大气密度,r为地心距,r=h+Re,h为跳跃高度,hf为开伞高度,V为返回器相对地球速度;地球平均半径Re为6378.135km,海平面重力加速度g0为9.81m/s2;θ为经度,φ为纬度,γ为飞行路径角,即飞行速度矢量与当地水平方向的夹角;ψ为速度方位角,即飞行速度矢量与当地正北方向的夹角,顺时针为正;σ为倾侧角,从返回器内部来看,右侧为正;D,L分别为无量纲化的气动阻力和气动升力;CD,CL分别为阻力系数和升力系数;返回器最大横截面积为Sref,质量为m;参数上方的“.”代表一阶导数,即对时间t的无量纲值求导;S为飞行航程;Ω为地球自转角速度;r和V为无量纲值,分别通过实际值除以Re、得到;
(2.2)采用步骤(1)确定的倾侧角剖面,数值积分步骤(2.1)确定的再入动力学方程(2),直到开伞点高度hf=10km,得到飞行航程S的积分终值:预测航程Sf,且Sf是关于当前点倾侧角幅值|σ0|的函数;
(2.3)将预测航程Sf与当前点的剩余航程Stogo进行比较,若预测航程Sf的偏差Smiss=Stogo-Sf(σ0)满足预设的着陆精度要求,即以|σ0|作为当前点倾侧角幅值;若航程偏差Smiss不满足预设的着陆精度要求,则需要更新当前倾侧角幅值|σ0|;
更新当前倾侧角幅值|σ0|即对方程Smiss=Stogo-Sf(σ0)=f(cosσ0)=0求根;
采用割线法迭代求解方程Smiss=Stogo-Sf(σ0)=f(cosσ0)=0,迭代表达式为
其中,xi为第i步迭代的值且f(xi)为第i步所求落点偏差Smiss,即预测落点与实际落点的偏差。
(3)如图2所示,将返回器的跳跃式再入飞行轨迹分为以下几个阶段:初始再入段、首次再入下降段、首次再入上升段、开普勒段和最终再入段,并且确定返回器的最终制导律为:当返回器处于初始再入段、首次再入下降段、首次再入上升段或最终再入段时,使用步骤(2)确定的所述当前点倾侧角σ0的幅值|σ0|进行制导;当返回器处于开普勒段时,使用|σ0|=180°进行制导,且倾侧角σ0符号由横程制导逻辑决定,这里均采用阿波罗再入制导方案,从而完成所述深空探测返回过程的预测校正制导。
初始再入段,首次再入下降段,首次再入上升段,开普勒段,最终再入段,主要特点分别体现为:
当过载n为0.05<n<0.2时,返回器进入初始再入段;
当过载n>0.2且飞行路径角γ<0时,返回器进入首次再入下降段;
当过载n>0.2且γ>0时,返回器进入首次再入上升段;
当过载n再次小于0.2时,返回器进入开普勒段;
当过载n再次大于0.2时,返回器进入最终再入段;
其中,
以下分析当返回器处于开普勒段时,使用|σ0|=180°进行制导的原理:
(a)分析飞行航程S与过载n之间的关系;
忽略地球自转的影响,处理(Re+hf)/r≈1,只考虑纵向剖面运动,简化动力学方程(2)为:
返回器的在轨能量e无量纲化后可表示为
式(4)中,r一直约等于1,V从约变化至0,可见能量从1变化至0。又因为
由简化动力学方程(3)可得
由方程(3)和(6),且假设cosγ≈1,有
则总的飞行航程为
由式(8)可见在能量初值ei和终端值ef相同的情况下,增大飞行航程,则气动阻力D将减小,过载 将得到抑制(升阻比L/D=CL/CD为常值)。因此,若能在气动力很小的开普勒段施加制导,减小开普勒段的飞行航程,则能有效增加最终再入段(开普勒段之后)的航程,实现最终再入段的过载抑制。
(b)分析在开普勒段如何调整倾侧角,以减小开普勒段飞行航程;
返回器再次跳出大气层,进入开普勒段飞行轨迹形式如图3所示,A为航程角;γexit、Vexit、rexit分别为跳出点处的飞行路径角、相对速度和地心距;μ为万有引力常数。值得注意的是,飞行器要完成跳出并再次再入大气层的跳跃式再入轨迹,跳出点处的飞行路径角必然大于零,速度也必须小于第一宇宙速度,且有
结合文献中的航程角计算公式可知
设小升力体返回器的升阻比L/D=CL/CD约为0.3,当过载n<0.2时,根据可得升力L<0.06,沿着矢径方向升力分量Lcosσcosγ显然远远小于万有引力。因此,这里将该升力分量视为引力的扰动项ε,且扰动力取指向地心为正,则等效的万有引力常数μ′为
显然,随着倾侧角σ的增大,ε和μ′也随之增大。又因为
由式(12)可知航程角A随着μ′的增大而减小,可知开普勒段航程随着倾侧角σ的增大而减小,最终再入段航程将随着倾侧角σ的增大而增大。因此返回器处于开普勒段时,使用|σ0|=180°进行制导,以最大程度实现剩余段飞行的过载抑制。
在开普勒段采用升力完全向下的制导律,能在抑制过载的同时,更快速地使飞行路径角减小至小于零,从而在该过程耗散了返回器的能量,确保其能安全着陆。
以探月返回器为背景,通过仿真,说明本发明所述方法的有效性。再入返回器质量为9500kg,最大横截面积为23.8m2。当马赫数大于25时,配平攻角保持约为160.2°,此时的升阻比约为0.2887。随着马赫数的变化,升阻比在整个再入过程中大约在0.2229和0.4075之间变化。地球半径为6378.135km,重力常数为9.8N/kg,其它参数设置见表1。
表1仿真参数设置表
首先对如下四种情况在标称参数下进行仿真,以比较整个再入过程的过载和着陆精度:
Case1:开普勒段,当n<0.2时,|σ|=0°;
Case2:开普勒段,σ幅值继续由预测校正算法获得;
Case3:开普勒段,当n<0.2时,|σ|=90°;
Case4:开普勒段,当n<0.2时,|σ|=180°。
由图4可知上述四种情况在进入开普勒段之前制导指令相同,因此进入开普勒段时各状态变量的值均相等,进入开普勒段后***分别采用上述四种不同的制导律。图5可知最终再入段的航程随着开普勒段倾侧角σ的增大而增大,验证了原理步骤(b)中所述。图6反映了在最终再入段Case4的航程最大,该段飞行最平缓。图7可见跳跃式再入过程过载随时间变化的曲线会出现前后两个峰值,为方便后续描述,分别称为前峰值和后峰值,该图显示Case4过载后峰值最小。图8为过载值的历经时间(飞行过程中大于等于某过载值的总飞行时间)曲线,可见Case4过载抑制效果最为显著,满足了最严格的再入过载约束(再入时航天员身体条件欠佳时所需满足的过载值历经时间约束)。
为了验证该方法在各类偏差(初始条件偏差、气动特性偏差、质量特性偏差等)的扰动下的鲁棒性和制导精度,这里在仿真中引入参数偏差,其分布类型如表2所示。
表2蒙特卡洛仿真中的误差源分布情况
将case2和case4均进行100次MonteCarlo打靶仿真,比较它们在偏差情况下的过载和着陆精度。图9(a)、图9(b)体现过载比较的情况,case4采用了融合制导算法,经统计80%的曲线满足了过载约束要求,明显优于case2的情况。且从反映着陆精度对比的曲线图8可知case4中加入过载抑制算法着陆精度较高,各散布仿真均满足精度为5km的要求。图9(a)、图9(b)、图10(a)、图10(b)中出现超出约束范围的情况,主要是因为在散布仿真时,存在气动系数、质量、大气密度等偏差效果的叠加,引起了较大的过载峰值及航程偏差。
仿真验证了该方法可简单有效地实现过载抑制,散布仿真也表明了该算法的实用性和鲁棒性。
Claims (3)
1.一种深空探测返回过程的预测校正制导方法,其特征在于步骤如下:
(1)确定所述深空探测返回过程轨迹预测采用的倾侧角剖面;
所述倾侧角剖面的数学表达式为
其中为总的目标飞行航程,Stogo为当前点的剩余航程,σf为70°,Sthres=2000km,σ0为当前点倾侧角,σ为倾侧角;
(2)迭代计算当前点倾侧角σ0的幅值,具体为:
(2.1)确定考虑地球自转影响的三自由度无量纲再入动力学方程为:
其中,ρ为大气密度,r为地心距,r=h+Re,h为跳跃高度,hf为开伞高度,V为返回器相对地球速度;地球平均半径Re为6378.135km,海平面重力加速度g0为9.81m/s2;θ为经度,φ为纬度,γ为飞行路径角,即飞行速度矢量与当地水平方向的夹角;ψ为速度方位角,即飞行速度矢量与当地正北方向的夹角,顺时针为正;σ为倾侧角,从返回器内部来看,右侧为正;D,L分别为无量纲化的气动阻力和气动升力;CD,CL分别为阻力系数和升力系数;返回器最大横截面积为Sref,质量为m;参数上方的“.”代表一阶导数;S为飞行航程;Ω为地球自转角速度;r和V为无量纲值;
(2.2)采用步骤(1)确定的倾侧角剖面,数值积分步骤(2.1)确定的再入动力学方程,直到开伞点高度hf=10km,得到飞行航程S的积分终值:预测航程Sf,且Sf是关于当前点倾侧角σ0的函数;
(2.3)将预测航程Sf与当前点的剩余航程Stogo进行比较,若预测航程Sf的偏差Smiss=Stogo-Sf(σ0)满足预设的着陆精度要求,即以|σ0|作为当前点倾侧角幅值;若航程偏差Smiss不满足预设的着陆精度要求,则需要更新当前倾侧角幅值|σ0|;
(3)将返回器的跳跃式再入飞行轨迹分为以下几个阶段:初始再入段、首次再入下降段、首次再入上升段、开普勒段和最终再入段,并且确定返回器的最终制导律为:当返回器处于初始再入段、首次再入下降段、首次再入上升段或最终再入段时,使用步骤(2)确定的当前点倾侧角σ0的幅值|σ0|进行制导;当返回器处于开普勒段时,使用|σ0|=180°进行制导,从而完成所述深空探测返回过程的预测校正制导。
2.根据权利要求1所述的一种深空探测返回过程的预测校正制导方法,其特征在于在:所述步骤(2.3)中更新当前倾侧角幅值|σ0|具体为:
更新当前倾侧角幅值|σ0|即对方程Smiss=Stogo-Sf(σ0)=f(cosσ0)=0求根;
采用割线法迭代求解方程Smiss=Stogo-Sf(σ0)=f(cosσ0)=0,迭代表达式为
其中,xi为第i步迭代的值且xi=(cosσ0)i,f(xi)为第i步落点偏差Smiss,即预测落点与实际落点的偏差。
3.根据权利要求1所述的一种深空探测返回过程的预测校正制导方法,其特征在于在:所述步骤(3)中初始再入段、首次再入下降段、首次再入上升段、开普勒段和最终再入段具体为:
当过载n为0.05<n<0.2时,返回器进入初始再入段;
当过载n>0.2且飞行路径角γ<0时,返回器进入首次再入下降段;
当过载n>0.2且γ>0时,返回器进入首次再入上升段;
当过载n再次小于0.2时,返回器进入开普勒段;
当过载n再次大于0.2时,返回器进入最终再入段;
其中,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410125842.XA CN103863579B (zh) | 2014-03-31 | 2014-03-31 | 一种深空探测返回过程的预测校正制导方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410125842.XA CN103863579B (zh) | 2014-03-31 | 2014-03-31 | 一种深空探测返回过程的预测校正制导方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103863579A CN103863579A (zh) | 2014-06-18 |
CN103863579B true CN103863579B (zh) | 2015-11-25 |
Family
ID=50902707
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410125842.XA Active CN103863579B (zh) | 2014-03-31 | 2014-03-31 | 一种深空探测返回过程的预测校正制导方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103863579B (zh) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104648695B (zh) * | 2014-12-18 | 2016-08-24 | 北京控制工程研究所 | 一种基于倾侧角可用性的再入走廊最优规划方法 |
CN104597756B (zh) * | 2014-12-18 | 2016-03-02 | 北京控制工程研究所 | 一种跳跃式再入二次再入段航程预估方法 |
CN104850129B (zh) * | 2014-12-19 | 2017-10-24 | 北京控制工程研究所 | 一种跳跃式再入的射向预偏置横向制导方法 |
CN105115512B (zh) * | 2015-09-23 | 2017-10-10 | 北京理工大学 | 一种火星大气进入段侧向预测校正制导方法 |
CN105223962B (zh) * | 2015-10-27 | 2017-09-01 | 湖北航天技术研究院总体设计所 | 一种减小外形跳变时弹体过载的动态攻角限幅方法 |
CN105629731B (zh) * | 2016-01-26 | 2018-10-12 | 西安电子科技大学 | 一种低升阻比返回器混合制导卸载方法 |
CN106054604B (zh) * | 2016-06-02 | 2018-12-25 | 西北工业大学 | 基于模型预测控制理论的再入飞行器鲁棒最优制导方法 |
CN106371312B (zh) * | 2016-09-12 | 2019-09-20 | 中国人民解放军国防科学技术大学 | 基于模糊控制器的升力式再入预测-校正制导方法 |
CN107804487B (zh) * | 2017-10-16 | 2020-09-01 | 中国西安卫星测控中心 | 一种基于自适应偏差控制的跳跃式再入返回落点预报方法 |
CN107861517B (zh) * | 2017-11-01 | 2020-06-16 | 北京航空航天大学 | 基于线性伪谱的跳跃式再入飞行器在线弹道规划制导方法 |
CN108036676B (zh) * | 2017-12-04 | 2019-08-23 | 北京航空航天大学 | 一种基于三维再入弹道解析解的全射向自主再入制导方法 |
CN108674695B (zh) * | 2018-07-19 | 2021-01-05 | 北京空间技术研制试验中心 | 航天器再入返回轨道规划方法 |
CN109446705B (zh) * | 2018-11-09 | 2023-04-07 | 北京空间技术研制试验中心 | 航天器再入段航程预报方法 |
CN109941460B (zh) * | 2019-04-09 | 2020-08-07 | 北京空间技术研制试验中心 | 航天器亚轨道返回再入过载降低设计方法 |
CN110413000B (zh) * | 2019-05-28 | 2020-10-13 | 北京航空航天大学 | 一种基于深度学习的高超声速飞行器再入预测校正容错制导方法 |
CN110304278B (zh) * | 2019-06-20 | 2022-05-24 | 南京航空航天大学 | 一种跳跃式再入轨迹的高精度优化方法 |
CN111361760B (zh) * | 2020-01-13 | 2021-07-02 | 北京理工大学 | 小天体表面移动轨迹抗差优化方法 |
CN112141366B (zh) * | 2020-09-23 | 2022-03-25 | 西北工业大学 | 一种地球轨道中航天器的轨迹预测方法及*** |
CN112269390B (zh) * | 2020-10-15 | 2021-09-21 | 北京理工大学 | 考虑弹跳的小天体表面定点附着轨迹规划方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2170194C1 (ru) * | 2000-07-20 | 2001-07-10 | Государственный космический научно-производственный центр им. М.В. Хруничева | Способ старта ракеты |
CN103587723A (zh) * | 2013-11-07 | 2014-02-19 | 北京临近空间飞行器***工程研究所 | 一种再入初始段解析式纵向在线轨迹设计及跟踪方法 |
CN103662090A (zh) * | 2013-12-13 | 2014-03-26 | 北京控制工程研究所 | 一种智能动力下降轨迹在线规划方法 |
CN103662092A (zh) * | 2013-12-13 | 2014-03-26 | 北京控制工程研究所 | 一种衔接主减速和接近段的预测校正方法 |
-
2014
- 2014-03-31 CN CN201410125842.XA patent/CN103863579B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2170194C1 (ru) * | 2000-07-20 | 2001-07-10 | Государственный космический научно-производственный центр им. М.В. Хруничева | Способ старта ракеты |
CN103587723A (zh) * | 2013-11-07 | 2014-02-19 | 北京临近空间飞行器***工程研究所 | 一种再入初始段解析式纵向在线轨迹设计及跟踪方法 |
CN103662090A (zh) * | 2013-12-13 | 2014-03-26 | 北京控制工程研究所 | 一种智能动力下降轨迹在线规划方法 |
CN103662092A (zh) * | 2013-12-13 | 2014-03-26 | 北京控制工程研究所 | 一种衔接主减速和接近段的预测校正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103863579A (zh) | 2014-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103863579B (zh) | 一种深空探测返回过程的预测校正制导方法 | |
Brunner et al. | Skip entry trajectory planning and guidance | |
CN104035335B (zh) | 基于高精度纵、横程解析预测方法的平稳滑翔再入制导方法 | |
US9671789B1 (en) | Method and system for anti-disturbance composite on-line guidance for atmosphere-entering phase of a Mars lander | |
Xue et al. | Constrained predictor-corrector entry guidance | |
Li et al. | Review and prospect of guidance and control for Mars atmospheric entry | |
CN107941087A (zh) | 一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法 | |
CN106371312B (zh) | 基于模糊控制器的升力式再入预测-校正制导方法 | |
CN105005313B (zh) | 一种基于路径点规划的火星大气进入段预测制导方法 | |
CN110015446B (zh) | 一种半解析的火星进入制导方法 | |
CN109080854B (zh) | 航天器返回预定落点的大椭圆轨道变轨规划方法 | |
Li et al. | Radio beacons/IMU integrated navigation for Mars entry | |
Zhang et al. | Multi-constrained suboptimal powered descent guidance for lunar pinpoint soft landing | |
CN109927941B (zh) | 一种基于预测离轨点精度的自主允许离轨判断方法 | |
Li et al. | Time-coordination entry guidance for multi-hypersonic vehicles | |
CN109941460B (zh) | 航天器亚轨道返回再入过载降低设计方法 | |
CN104635488A (zh) | 一种火星大气进入段鲁棒预测制导方法 | |
CN109101725B (zh) | 航天器受控再入落区预示方法 | |
Figueroa-González et al. | Landing guidance strategy of space rider | |
Benito et al. | Nonlinear predictive controller for drag tracking in entry guidance | |
CN103921957B (zh) | 一种探月飞船跳跃式再入的跃起点能量管理方法 | |
Shen et al. | On-board entry trajectory planning for sub-orbital flight | |
CN107102547A (zh) | 一种基于滑模控制理论的rlv着陆段制导律获取方法 | |
CN107796401B (zh) | 跳跃式再入飞行器线性伪谱参数修正横向制导方法 | |
Bairstow et al. | Orion reentry guidance with extended range capability using PredGuid |
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 |