CN112660426A - 一种火箭软着陆制导方法 - Google Patents
一种火箭软着陆制导方法 Download PDFInfo
- Publication number
- CN112660426A CN112660426A CN202011480607.6A CN202011480607A CN112660426A CN 112660426 A CN112660426 A CN 112660426A CN 202011480607 A CN202011480607 A CN 202011480607A CN 112660426 A CN112660426 A CN 112660426A
- Authority
- CN
- China
- Prior art keywords
- rocket
- soft landing
- equation
- engine
- constraint equation
- 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
Images
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明提供了一种火箭软着陆制导方法,属于运载火箭控制技术领域。本发明通过自适应发动机开机方法获取发动机开机指令,为发动机开机后的着陆过程创造良好初始接入条件,而且根据火箭着陆飞行状态建立的火箭动力软着陆段在线轨迹规划方程,可以实时在线滚动规划后续飞行过程偏差适应能力最强的标称轨迹;通过制导跟踪方法,可以将火箭导引向实时生成的标称轨迹,大大降低火箭着陆过程中干扰的影响;此外,本发明还通过小步长预测关机方法确定发动机关机条件,进一步提升火箭的着陆精度。
Description
技术领域
本发明属于运载火箭控制技术领域,更具体地说,是涉及一种火箭软着陆制导方法。
背景技术
目前能够实现火箭垂直着陆的Falcon系列火箭,其一级或助推级均配置有9台Merlin 1D发动机,通过在着陆过程中仅开启其中几台而不是全部发动机既能够大幅降低推重比,还可以保持推力对称,维持姿态控制稳定。火箭通过调节发动机推力,在着陆前使推重比等于1,以安全落速实现匀速着陆。然而,目前很多现役火箭一级或助推级均只并联了少量发动机,如中国CZ-5、CZ-7火箭,欧洲Ariane 6,日本H3火箭等,一级均只配置了两台发动机。若着陆过程中关闭一台,则很难保持姿态稳定;若不关闭则要求发动机具备深度节流能力,而大推力发动机深度节流一直是一个难题。如果制导方法在最小推力仍大于火箭重力的大推重比条件下仍能够实现火箭定点软着陆,则可以缓解发动机深度节流的压力,使得在役火箭通过发动机节流能力的适当调整,就可以经济的升级为可重复使用运载火箭,避免为了重复使用而设计多台发动机并联的新火箭。但是大推重比条件下,火箭动力着陆段的物理可行域被大幅度压缩,在干扰状态下易飞出可行域,从而难以实现定点软着陆。
发明内容
本发明的目的在于提供一种火箭软着陆制导方法,旨在解决现有火箭垂直着陆时,在干扰状态下易飞出可行域,导致软着陆位置不精确的问题。
为实现上述目的,本发明采用的技术方案是:
一种火箭软着陆制导方法,包括:
采用自适应发动机开机方法获取发动机开机指令;
依照所述发动机开机指令打开火箭发动机,火箭进入软着陆阶段;
根据火箭软着陆运动过程建立火箭动力软着陆段第一运动方程和火箭动力软着陆段第一状态约束方程;
根据所述火箭动力软着陆段第一运动方程和所述火箭动力软着陆段第一状态约束方程建立火箭动力软着陆段在线轨迹规划方程;
获取火箭的当前状态量和当前控制量;所述当前状态量包括速度和位置;所述当前控制量包括推力幅值和程序角;
根据所述当前状态量和所述火箭动力软着陆段在线轨迹规划方程得到火箭的理想控制量;所述理想控制量包括理想推力幅值和理想程序角;
采用制导跟踪方法修正理想控制量与当前控制量的偏差;
获取火箭动力软着陆段的标准速度关机量和标准高度关机量;
根据所述标准速度关机量、所述标准高度关机量以及小步长预测关机方法得到火箭发动机关机指令;
依据所述发动机关机指令关闭发动机,火箭软着陆完成。
优选的,所述采用自适应发动机开机方法获取发动机开机指令包括:
根据火箭软着陆运动过程建立火箭动力软着陆段第二运动方程和火箭动力软着陆段第二状态约束方程;
根据所述火箭动力软着陆段第二运动方程和所述火箭动力软着陆段第二状态约束方程建立火箭理想开机速度轨迹规划方程;
获取火箭当前位置和当前速度;
根据所述火箭当前位置和所述火箭理想开机速度轨迹规划方程得到火箭的当前位置理想着陆时间和当前位置理想开机速度;
根据所述当前位置理想着陆时间、所述当前位置理想开机速度和所述当前速度确定所述发动机开机指令。
优选的,所述火箭动力软着陆段第一运动方程包括:
D=0.5ρSrefCD||V||V
其中,r=[x,y,z]T为位置矢量,x为火箭在x轴的位置,y为火箭在y轴的位置,z为火箭在z轴的位置,V=[Vx,Vy,Vz]T为速度矢量,Vx为火箭在x轴的速度,Vy为火箭在y轴的速度,Vz为火箭在z轴的速度,m为火箭质量,ψ为偏航角,为俯仰角,ωψ为偏航角速度,为俯仰角速度,T为发动机推力幅值,δT为发动机推力幅值变化率,D为气动阻力矢量,ρ为大气密度,Sref为参考面积,CD为气动阻力系数,g为重力加速度的矢量,Isp为发动机比冲,g0为海平面重力加速度。
优选的,所述火箭动力软着陆段第一状态约束方程包括:第一初始状态约束方程、第一终端状态约束方程和第一过程约束方程;
所述第一初始状态约束方程包括:
[r0,V0,m0]=[r,V,m](t0)
其中,t0表示当前时刻,r0表示当前时刻火箭的位置,V0表示当前时刻火箭的速度,m0表示当前时刻火箭的质量;
所述第一终端状态约束方程包括:
Vyf≤Vy(tf)≤0,|Vx(tf)|≤Vxf,|Vz(tf)|≤Vzf,m(tf)≥mmin
其中,tf表示终端时间,yf表示火箭在y轴期望的终端位置,xf表示火箭在x轴期望的终端位置,zf表示火箭在z轴期望的终端位置,Vyf为火箭在y轴的终端速度下界,Vxf为火箭在x轴期望的终端速度,Vzf为火箭在z轴期望的终端速度,Ty表示火箭在y轴方向上的发动机推力幅值,表示火箭在终端位置期望的俯仰角,mmin表示火箭最小质量;
所述第一过程约束方程包括:
优选的,所述根据所述火箭动力软着陆段第一运动方程和所述火箭动力软着陆段第一状态约束方程建立火箭动力软着陆段在线轨迹规划方程,包括:
联立所述火箭动力软着陆段第一运动方程和所述火箭动力软着陆段第一状态约束方程,同时以动力软着陆段燃料最省为优化目标建立所述火箭动力软着陆段在线轨迹规划方程;
所述火箭动力软着陆段在线轨迹规划方程包括:
min J=-m(tf)
[r0,V0,m0]=[r,V,m](t0),y(tf)=yf,|x(tf)|≤xf,|z(tf)|≤zf,
Vyf≤Vy(tf)≤0,|Vx(tf)|≤Vxf,|Vz(tf)|≤Vzf,m(tf)≥mmin,
Tmin≤T≤Tmax,|δT(t)|≤δTmax,y(t)≥0,Vy(t)≤0。
优选的,所述火箭动力软着陆段第二运动方程包括:
D=0.5ρSrefCD||V||V
优选的,火箭动力软着陆段第二状态约束方程包括:第二初始状态约束方程、第二终端状态约束方程和第二过程约束方程:
所述第二初始状态约束方程包括:
r0=r(t0),m0=m(t0).
其中,t0表示初始点时间,r0表示火箭初始位置,m0表示火箭初始质量;
所述第二终端状态约束方程包括:
rf表示火箭期望的终端位置,Vf表示火箭在期望的终端位置处的速度;
所述第二过程约束方程包括:
其中,Tmid为发动机推力调节范围的中值。
优选的,所述根据所述火箭动力软着陆段第二运动方程和所述火箭动力软着陆段第二状态约束方程建立火箭理想开机速度轨迹规划方程,包括:
联立火箭动力软着陆段第二运动方程和所述火箭动力软着陆段第二状态约束方程,同时以动力软着陆段燃料最省为优化目标建立所述火箭理想开机速度轨迹规划方程;
所述火箭理想开机速度轨迹规划方程包括:
minJ=-m(tf)
本发明提供的一种火箭软着陆制导方法的有益效果在于:与现有技术相比,本发明提出了通过自适应发动机开机方法获取发动机开机指令,为发动机开机后的着陆过程创造良好初始接入条件,而且根据火箭着陆飞行状态建立的火箭动力软着陆段在线轨迹规划方程,可以实时在线滚动规划后续飞行过程偏差适应能力最强的标称轨迹;通过制导跟踪方法,可以将火箭导引向实时生成的标称轨迹,大大降低火箭着陆过程中干扰的影响;此外,本发明还通过小步长预测关机方法确定发动机关机条件,进一步提升火箭的着陆精度。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种火箭软着陆制导方法流程图;
图2为本发明实施例提供的一种火箭软着陆制导方法原理图。
具体实施方式
为了使本发明所要解决的技术问题、技术方案及有益效果更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
请一并参阅图1、2,现对本发明提供的一种火箭软着陆制导方法进行说明。一种火箭软着陆制导方法,包括:
步骤101:采用自适应发动机开机方法获取发动机开机指令;
为得到火箭进入动力软着陆段的发动机开机条件,本发明提出了一种可重复使用运载火箭动力软着陆段自适应发动机开机方法,使得火箭能够在可行域最大的条件下进入动力软着陆段,提升动力软着陆段的偏差适应能力。
步骤101具体包括:
S1:根据火箭软着陆运动过程建立火箭动力软着陆段第二运动方程和火箭动力软着陆段第二状态约束方程。
下面对着一过程进行进一步描述,定义目标坐标系,在目标坐标系下描述火箭动力软着陆段的运动。目标坐标系原点O在着陆点,OY轴垂直目标点当地水平面指向天,OX轴在目标点当地水平面内指向发射点,OZ轴满足右手定则。在描述动力软着陆段三自由度质心运动方程时,将火箭视为质点,考虑发动机推力、气动力以及质量变化对火箭运动过程的影响,忽略绕质心姿态运动的动态过程,认为火箭姿态能够实时响应程序角指令。定义发动机推力始终沿火箭轴向,推力矢量与OXY平面之间的夹角为偏航角ψ,推力矢量在OXY平面内的投影与OX轴之间的夹角为俯仰角
建立的火箭动力软着陆段第二运动方程如下:
其中,r=[x,y,z]T为位置矢量,x为火箭在x轴的位置,y为火箭在y轴的位置,z为火箭在z轴的位置,V=[Vx,Vy,Vz]T为速度矢量,Vx为火箭在x轴的速度,Vy为火箭在y轴的速度,Vz为火箭在z轴的速度,m为火箭质量,ψ为偏航角,为俯仰角,ωψ为偏航角速度,为俯仰角速度,T为发动机推力幅值,D为气动阻力矢量,ρ为大气密度,Sref为参考面积,CD为气动阻力系数,g为重力加速度的矢量(重力加速度在目标系下的投影矢量),Isp为发动机比冲,g0为海平面重力加速度。
火箭动力软着陆段第二状态约束方程包括:第二初始状态约束方程、第二终端状态约束方程和第二过程约束方程;
第二初始状态约束方程包括:
r0=r(t0),m0=m(t0). (2)
其中,t0表示初始点时间,r0表示火箭初始位置,m0表示火箭初始质量;
第二终端状态约束方程包括:
第二过程约束方程包括:
其中,Tmid为发动机推力调节范围的中值,其中下标max表示着陆过程相应状态量和控制量允许的最大值,高度和速度不等式约束表示着陆过程中高度不能低于目标点,且纵向速度始终向下。
S2:根据火箭动力软着陆段第二运动方程和火箭动力软着陆段第二状态约束方程建立火箭理想开机速度轨迹规划方程;
进一步的,根据(1)-(4)中描述的约束条件,以动力软着陆段燃料最省为优化目标,构建可行域最大的理想开机速度轨迹规划问题,具体形式如下:
其中,minJ=-m(tf)为火箭动力软着陆段燃料优化目标函数。
S3:获取火箭当前位置和当前速度;
S4:根据火箭当前位置和火箭理想开机速度轨迹规划方程得到火箭的当前位置理想着陆时间和当前位置理想开机速度;
利用数值优化算法(例如:内点法或序列二次规划算法)求解可行域最大的理想开机速度轨迹规划问题,可得到最优解。则最优解的终端时间为tf opt,即为理想着陆时间;最优解的初始速度矢量V0 opt,即为对应当前位置的理想开机速度。
S5:根据当前位置理想着陆时间、当前位置理想开机速度和当前速度确定发动机开机指令。
具体的,定义动力软着陆段可允许最短飞行时间为tfmin,若tfmin≥tf opt,则发动机开机,否则评估火箭当前时刻速度与理想开机速度之间的偏差。定义开机特征量ΔVI,用于量化火箭当前时刻速度与理想开机速度之间的偏差,具体计算方式如下,
其中px、py、pz分别表示火箭在x、y、z轴三个方向上速度偏差的权重参数,可根据任务需求调整不同方向速度偏差占开机特征量的比重,进一步的,px表示火箭在x轴方向上速度偏差的权重参数,py表示火箭在y轴方向上速度偏差的权重参数,pz表示火箭在z轴方向上速度偏差的权重参数,Vx0 opt表示火箭的当前位置理想开机速度在x轴方向上的分量,Vy0 opt表示火箭的当前位置理想开机速度在y轴方向上的分量,Vz0 opt表示火箭的当前位置理想开机速度在z轴方向上的分量。
定义上一控制周期对应是时间为told,若ΔVI(t)≥ΔVI(told),则发动机开机,否则继续在下一控制周期,重复执行S3-S5,直至发动机开机为止。
本发明提出的自适应发动机开机方法,在线规划中值推力着陆的标准轨迹,得到当前位置对应的理想初始速度,通过计算当前速度与理想速度的偏差,在偏差最小时发出发动机开机指令,同时在线生成的推力中值标准轨迹可作为后续火箭动力软着陆段在线轨迹规划方程的初始猜测值,提升在线规划算法计算效率。
步骤102:依照发动机开机指令打开火箭发动机,火箭进入软着陆阶段;
步骤103:根据火箭软着陆运动过程建立火箭动力软着陆段第一运动方程和火箭动力软着陆段第一状态约束方程;
火箭动力软着陆段第一运动方程是在目标坐标系中建立的包括:
其中,r=[x,y,z]T为位置矢量,x为火箭在x轴的位置,y为火箭在y轴的位置,z为火箭在z轴的位置,V=[Vx,Vy,Vz]T为速度矢量,Vx为火箭在x轴的速度,Vy为火箭在y轴的速度,Vz为火箭在z轴的速度,m为火箭质量,ψ为偏航角,为俯仰角,ωψ为偏航角速度,为俯仰角速度,T为发动机推力幅值,δT为发动机推力幅值变化率,D为气动阻力矢量,ρ为大气密度,Sref为参考面积,CD为气动阻力系数,g为重力加速度的矢量(重力加速度在目标系下的投影矢量),Isp为发动机比冲,g0为海平面重力加速度。
火箭动力软着陆段第一状态约束方程包括:第一初始状态约束方程、第一终端状态约束方程和第一过程约束方程;
第一初始状态约束方程包括:
[r0,V0,m0]=[r,V,m](t0) (8)
其中,t0表示当前时刻,r0表示当前时刻火箭的位置,V0表示当前时刻火箭的速度,m0表示当前时刻火箭的质量;
第一终端状态约束方程包括:
Vyf≤Vy(tf)≤0,|Vx(tf)|≤Vxf,|Vz(tf)|≤Vzf,m(tf)≥mmin (9)
其中,tf表示终端时间,yf表示火箭在y轴期望的终端位置,xf表示火箭在x轴期望的终端位置,zf表示火箭在z轴期望的终端位置,Vyf为火箭在y轴的终端速度下界,Vxf为火箭在x轴期望的终端速度,Vzf为火箭在z轴期望的终端速度,Ty表示火箭在y轴方向上的发动机推力幅值,表示火箭在终端位置期望的俯仰角,mmin表示火箭最小质量;
第一过程约束方程包括:
其中,表示俯仰角的最大值,Tmin表示发动机推力幅值的最小值,Tmax表示发动机推力幅值的最大值,δTmax表示发动机推力幅值变化率的最大值,为俯仰角速度的最大值,ωψmax为偏航角速度的最大值,其中下标max表示着陆过程相应状态量和控制量允许的最大值,高度和速度不等式约束表示着陆过程中高度不能低于目标点,且纵向速度始终向下。
步骤104:根据火箭动力软着陆段第一运动方程和火箭动力软着陆段第一状态约束方程建立火箭动力软着陆段在线轨迹规划方程;
根据(7)-(10)中描述的约束条件,以动力软着陆段燃料最省为优化目标,构建可行域最大的动力软着陆段轨迹规划命题,具体形式如下,
min J=-m(tf)
[r0,V0,m0]=[r,V,m](t0),y(tf)=yf,|x(tf)|≤xf,|z(tf)|≤zf,
Vyf≤Vy(tf)≤0,|Vx(tf)|≤Vxf,|Vz(tf)|≤Vzf,m(tf)≥mmin,
Tmin≤T≤Tmax,|δT(t)|≤δTmax,y(t)≥0,Vy(t)≤0
火箭在发动机开机后进入动力软着陆段,本发明在描述规划问题时引入推力变化率的影响,通过限制推力变化率的幅值,构建可行域最大的火箭动力软着陆段在线轨迹规划方程。在每个规划周期,根据火箭当前时刻的飞行状态,采用非线性规划算法(如牛顿法、内点法、序列二次规划等)滚动求解满足约束条件的离散化最优飞行轨迹状态量和推力控制矢量,使得最优推力幅值控制量变化尽可能平缓,并与最大值和最小值之间留出余量,使得推力始终具备正向和负向的调节能力,应对后续飞行过程中可能出现的干扰或偏差,从而在物理可行域内尽可能提升火箭对后续飞行过程不确定性的适应能力。
步骤105:获取火箭的当前状态量和当前控制量;当前状态量包括速度和位置;当前控制量包括推力幅值和程序角;
步骤106:根据当前状态量和火箭动力软着陆段在线轨迹规划方程得到火箭的理想控制量;理想控制量包括理想推力幅值和理想程序角;
步骤107:采用制导跟踪方法修正理想控制量与当前控制量的偏差;
进一步的,由于制导周期比规划周期短,因此在火箭当前飞行时刻,根据任务特点选取合适的插值自变量(如时间、速度、位置等),在火箭动力软着陆段在线轨迹规划方程最新生成的标准轨迹状态序列上,插值(如线性插值、二次曲线插值、三次样条等)得到火箭当前时刻的理想状态量和理想控制量。根据火箭当前实际状态与理想状态的偏差,利用制导跟踪方法(如PID控制、LQR控制等)实时修正理想控制量,可以大大提高火箭的着陆精度。
步骤108:获取火箭动力软着陆段的标准速度关机量和标准高度关机量;
步骤109:根据标准速度关机量、标准高度关机量以及小步长预测关机方法得到火箭发动机关机指令;
结合火箭发动机关机后效影响和火箭着陆时的理想质量,计算能够以理想速度(例如-1m/s)实现着陆对应的关机时刻速度、高度值,作为火箭动力软着陆段的标准速度关机量和标准高度关机量。在火箭进入动力软着陆段后,根据当前速度和高度,以及每个制导周期的速度变化率和高度变化率,分别预测到达标准速度关机量和标准高度关机量的剩余时间,并取两者最小值作为预测待飞时间。为降低制导周期对关机精度的影响,采用小步长预测关机方法,当预测待飞时间小于N个制导周期(如N可取1、2、3)后,启动1ms计时,使得关机指令的发出精度小于1ms,提升着陆精度。
步骤1010:依据发动机关机指令关闭发动机,火箭软着陆完成。
本发明提供的一种火箭软着陆制导方法的有益效果在于:与现有技术相比,本发明提出了通过自适应发动机开机方法获取发动机开机指令,为发动机开机后的着陆过程创造良好初始接入条件,而且根据火箭着陆飞行状态建立的火箭动力软着陆段在线轨迹规划方程,可以实时在线滚动规划后续飞行过程偏差适应能力最强的标称轨迹;通过制导跟踪方法,可以将火箭导引向实时生成的标称轨迹,大大降低火箭着陆过程中干扰的影响;此外,本发明还通过小步长预测关机方法确定发动机关机条件,进一步提升火箭的着陆精度,且本发明适用于在大推重比条件下,火箭动力反推软着陆过程。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种火箭软着陆制导方法,其特征在于,包括:
采用自适应发动机开机方法获取发动机开机指令;
依照所述发动机开机指令打开火箭发动机,火箭进入软着陆阶段;
根据火箭软着陆运动过程建立火箭动力软着陆段第一运动方程和火箭动力软着陆段第一状态约束方程;
根据所述火箭动力软着陆段第一运动方程和所述火箭动力软着陆段第一状态约束方程建立火箭动力软着陆段在线轨迹规划方程;
获取火箭的当前状态量和当前控制量;所述当前状态量包括速度和位置;所述当前控制量包括推力幅值和程序角;
根据所述当前状态量和所述火箭动力软着陆段在线轨迹规划方程得到火箭的理想控制量;所述理想控制量包括理想推力幅值和理想程序角;
采用制导跟踪方法修正理想控制量与当前控制量的偏差;
获取火箭动力软着陆段的标准速度关机量和标准高度关机量;
根据所述标准速度关机量、所述标准高度关机量以及小步长预测关机方法得到火箭发动机关机指令;
依据所述发动机关机指令关闭发动机,火箭软着陆完成。
2.如权利要求1所述的一种火箭软着陆制导方法,其特征在于,所述采用自适应发动机开机方法获取发动机开机指令包括:
根据火箭软着陆运动过程建立火箭动力软着陆段第二运动方程和火箭动力软着陆段第二状态约束方程;
根据所述火箭动力软着陆段第二运动方程和所述火箭动力软着陆段第二状态约束方程建立火箭理想开机速度轨迹规划方程;
获取火箭当前位置和当前速度;
根据所述火箭当前位置和所述火箭理想开机速度轨迹规划方程得到火箭的当前位置理想着陆时间和当前位置理想开机速度;
根据所述当前位置理想着陆时间、所述当前位置理想开机速度和所述当前速度确定所述发动机开机指令。
3.如权利要求1所述的一种火箭软着陆制导方法,其特征在于,所述火箭动力软着陆段第一运动方程包括:
D=0.5ρSrefCD||V||V
4.如权利要求3所述的一种火箭软着陆制导方法,其特征在于,所述火箭动力软着陆段第一状态约束方程包括:第一初始状态约束方程、第一终端状态约束方程和第一过程约束方程;
所述第一初始状态约束方程包括:
[r0,V0,m0]=[r,V,m](t0)
其中,t0表示当前时刻,r0表示当前时刻火箭的位置,V0表示当前时刻火箭的速度,m0表示当前时刻火箭的质量;
所述第一终端状态约束方程包括:
Vyf≤Vy(tf)≤0,|Vx(tf)|≤Vxf,|Vz(tf)|≤Vzf,m(tf)≥mmin
其中,tf表示终端时间,yf表示火箭在y轴期望的终端位置,xf表示火箭在x轴期望的终端位置,zf表示火箭在z轴期望的终端位置,Vyf为火箭在y轴的终端速度下界,Vxf为火箭在x轴期望的终端速度,Vzf为火箭在z轴期望的终端速度,Ty表示火箭在y轴方向上的发动机推力幅值,表示火箭在终端位置期望的俯仰角,mmin表示火箭最小质量;
所述第一过程约束方程包括:
5.如权利要求4所述的一种火箭软着陆制导方法,其特征在于,所述根据所述火箭动力软着陆段第一运动方程和所述火箭动力软着陆段第一状态约束方程建立火箭动力软着陆段在线轨迹规划方程,包括:
联立所述火箭动力软着陆段第一运动方程和所述火箭动力软着陆段第一状态约束方程,同时以动力软着陆段燃料最省为优化目标建立所述火箭动力软着陆段在线轨迹规划方程;
所述火箭动力软着陆段在线轨迹规划方程包括:
min J=-m(tf)
[r0,V0,m0]=[r,V,m](t0),y(tf)=yf,|x(tf)|≤xf,|z(tf)|≤zf,
Vyf≤Vy(tf)≤0,|Vx(tf)|≤Vxf,|Vz(tf)|≤Vzf,m(tf)≥mmin,
Tmin≤T≤Tmax,|δT(t)|≤δTmax,y(t)≥0,Vy(t)≤0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011480607.6A CN112660426B (zh) | 2020-12-15 | 2020-12-15 | 一种火箭软着陆制导方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011480607.6A CN112660426B (zh) | 2020-12-15 | 2020-12-15 | 一种火箭软着陆制导方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112660426A true CN112660426A (zh) | 2021-04-16 |
CN112660426B CN112660426B (zh) | 2021-09-14 |
Family
ID=75404923
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011480607.6A Active CN112660426B (zh) | 2020-12-15 | 2020-12-15 | 一种火箭软着陆制导方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112660426B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114132530A (zh) * | 2021-09-17 | 2022-03-04 | 北京空间飞行器总体设计部 | 一种基于触地关机的地外天体安全软着陆方法 |
CN114690793A (zh) * | 2022-03-28 | 2022-07-01 | 北京理工大学 | 基于滑模控制的可重复使用运载火箭垂直软着陆制导方法 |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2237583A1 (en) * | 1994-11-21 | 1996-05-30 | Michael S. Kelly | Space launch vehicles configured as gliders and towed to launch altitude by conventional aircraft |
JPH11301598A (ja) * | 1998-04-20 | 1999-11-02 | Mitsubishi Electric Corp | 宇宙航行体の誘導制御装置 |
US5987363A (en) * | 1996-03-26 | 1999-11-16 | California Institute Of Technology | Three-dimensional representation of a spacecraft's trajectory |
CN102035648A (zh) * | 2010-09-29 | 2011-04-27 | 北京航天自动控制研究所 | 一种飞行器状态信息的安全实时软回收系统及方法 |
CN102139768A (zh) * | 2010-10-28 | 2011-08-03 | 中国科学院力学研究所 | 一种亚轨道飞行器再入飞行的攻角制导方法 |
US20110297793A1 (en) * | 2008-12-22 | 2011-12-08 | Astrium Sas | Reusable module for launcher |
CN104590589A (zh) * | 2014-12-22 | 2015-05-06 | 哈尔滨工业大学 | 基于燃料最优的火星探测器着陆制导方法 |
CN105043171A (zh) * | 2015-06-30 | 2015-11-11 | 北京航天长征飞行器研究所 | 一种带倾角约束的火箭弹纵向导引方法 |
CN106227968A (zh) * | 2016-07-29 | 2016-12-14 | 北京航天自动控制研究所 | 一种航天器主发动机关机点优化方法 |
CN108100307A (zh) * | 2017-12-05 | 2018-06-01 | 北京理工大学 | 一种用于复杂约束下低能量小天体精确探测轨道转移方法 |
CN109229423A (zh) * | 2018-11-06 | 2019-01-18 | 中国科学院国家空间科学中心 | 一种返回式探空火箭的减速控制方法 |
CN109491246A (zh) * | 2018-11-19 | 2019-03-19 | 北京航天自动控制研究所 | 一种基于数值优化算法的自适应救援轨迹规划方法 |
CN110015446A (zh) * | 2019-03-05 | 2019-07-16 | 上海卫星工程研究所 | 一种半解析的火星进入制导方法 |
CN110989650A (zh) * | 2019-12-30 | 2020-04-10 | 航天科工火箭技术有限公司 | 一种基于火箭返回的着陆控制方法和装置 |
CN111753398A (zh) * | 2020-05-25 | 2020-10-09 | 北京理工大学 | 一种基于状态转移矩阵的平衡点转移轨道重构方法 |
-
2020
- 2020-12-15 CN CN202011480607.6A patent/CN112660426B/zh active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2237583A1 (en) * | 1994-11-21 | 1996-05-30 | Michael S. Kelly | Space launch vehicles configured as gliders and towed to launch altitude by conventional aircraft |
US5987363A (en) * | 1996-03-26 | 1999-11-16 | California Institute Of Technology | Three-dimensional representation of a spacecraft's trajectory |
JPH11301598A (ja) * | 1998-04-20 | 1999-11-02 | Mitsubishi Electric Corp | 宇宙航行体の誘導制御装置 |
US20110297793A1 (en) * | 2008-12-22 | 2011-12-08 | Astrium Sas | Reusable module for launcher |
CN102035648A (zh) * | 2010-09-29 | 2011-04-27 | 北京航天自动控制研究所 | 一种飞行器状态信息的安全实时软回收系统及方法 |
CN102139768A (zh) * | 2010-10-28 | 2011-08-03 | 中国科学院力学研究所 | 一种亚轨道飞行器再入飞行的攻角制导方法 |
CN104590589A (zh) * | 2014-12-22 | 2015-05-06 | 哈尔滨工业大学 | 基于燃料最优的火星探测器着陆制导方法 |
CN105043171A (zh) * | 2015-06-30 | 2015-11-11 | 北京航天长征飞行器研究所 | 一种带倾角约束的火箭弹纵向导引方法 |
CN106227968A (zh) * | 2016-07-29 | 2016-12-14 | 北京航天自动控制研究所 | 一种航天器主发动机关机点优化方法 |
CN108100307A (zh) * | 2017-12-05 | 2018-06-01 | 北京理工大学 | 一种用于复杂约束下低能量小天体精确探测轨道转移方法 |
CN109229423A (zh) * | 2018-11-06 | 2019-01-18 | 中国科学院国家空间科学中心 | 一种返回式探空火箭的减速控制方法 |
CN109491246A (zh) * | 2018-11-19 | 2019-03-19 | 北京航天自动控制研究所 | 一种基于数值优化算法的自适应救援轨迹规划方法 |
CN110015446A (zh) * | 2019-03-05 | 2019-07-16 | 上海卫星工程研究所 | 一种半解析的火星进入制导方法 |
CN110989650A (zh) * | 2019-12-30 | 2020-04-10 | 航天科工火箭技术有限公司 | 一种基于火箭返回的着陆控制方法和装置 |
CN111753398A (zh) * | 2020-05-25 | 2020-10-09 | 北京理工大学 | 一种基于状态转移矩阵的平衡点转移轨道重构方法 |
Non-Patent Citations (3)
Title |
---|
CONG WANG, ZHENGYU SONG: "Trajectory Optimization for Reusable Rocket Landing", 《2019 CHINESE AUTOMATION CONGRESS (CAC) 》 * |
MA LIN,WANG KEXIN, XU ZUHUA: "Optimal Feedback Guidance for Powered Landing of Reusable Rockets", 《THE 31TH CHINESE CONTROL AND DECISION CONFERENCE》 * |
王劲博: "火箭返回着陆问题高精度快速轨迹优化算法", 《控制理论与应用》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114132530A (zh) * | 2021-09-17 | 2022-03-04 | 北京空间飞行器总体设计部 | 一种基于触地关机的地外天体安全软着陆方法 |
CN114132530B (zh) * | 2021-09-17 | 2023-06-06 | 北京空间飞行器总体设计部 | 一种基于触地关机的地外天体安全软着陆方法 |
CN114690793A (zh) * | 2022-03-28 | 2022-07-01 | 北京理工大学 | 基于滑模控制的可重复使用运载火箭垂直软着陆制导方法 |
CN114690793B (zh) * | 2022-03-28 | 2023-01-06 | 北京理工大学 | 基于滑模控制的可重复使用运载火箭垂直软着陆制导方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112660426B (zh) | 2021-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112660426B (zh) | 一种火箭软着陆制导方法 | |
Shin et al. | Autonomous flight of the rotorcraft-based UAV using RISE feedback and NN feedforward terms | |
CN112550770B (zh) | 一种基于凸优化的火箭软着陆轨迹规划方法 | |
CN110989669A (zh) | 一种多级助推滑翔飞行器主动段在线自适应制导算法 | |
Kim et al. | Trajectory tracking controller design using neural networks for a tiltrotor unmanned aerial vehicle | |
Snyder et al. | Quad-segment polynomial trajectory guidance for impact-time control of precision-munition strike | |
CN112507461B (zh) | 一种运载火箭动力软着陆段发动机开机方法 | |
Bialy et al. | An adaptive backstepping controller for a hypersonic air-breathing missile | |
Ansari et al. | Retrospective cost adaptive control of generic transport model under uncertainty and failure | |
CN113885543A (zh) | 一种复杂海况干扰下的航行体自抗扰横滚控制器的设计方法 | |
CN109703769A (zh) | 一种基于预瞄策略的空中加油对接控制方法 | |
CN113504723B (zh) | 一种基于逆强化学习的运载火箭减载控制方法 | |
CN108255193A (zh) | 一种垂直/短距起降飞机飞行控制方法 | |
CN114637312A (zh) | 一种基于智能变形决策的无人机节能飞行控制方法及系统 | |
CN116301028B (zh) | 基于吸气式高超声速平台的多约束在线飞行轨迹规划中段导引方法 | |
Ansari et al. | Hybrid genetic algorithm fuzzy rule based guidance and control for launch vehicle | |
CN116697829A (zh) | 一种基于深度强化学习的火箭着陆制导方法及系统 | |
CN112629339B (zh) | 一种基于直接法的火箭软着陆轨迹规划方法 | |
Durham et al. | Minimum drag control allocation | |
CN114690793B (zh) | 基于滑模控制的可重复使用运载火箭垂直软着陆制导方法 | |
Anderson et al. | Design of an air to air interceptor using genetic algorithms | |
Li et al. | Actuators coupled design based adaptive backstepping control of air-breathing hypersonic vehicle | |
Fiorentini et al. | Nonlinear control of a hypersonic vehicle with structural flexibility | |
Duan et al. | Flight-path tracking control of a transportation aircraft: Comparison of two nonlinear design approaches | |
Wang et al. | Trajectory optimization for reusable rocket landing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |