CN110908396A - 可重复使用运载器全阶段再入返回制导方法 - Google Patents
可重复使用运载器全阶段再入返回制导方法 Download PDFInfo
- Publication number
- CN110908396A CN110908396A CN201911187178.0A CN201911187178A CN110908396A CN 110908396 A CN110908396 A CN 110908396A CN 201911187178 A CN201911187178 A CN 201911187178A CN 110908396 A CN110908396 A CN 110908396A
- Authority
- CN
- China
- Prior art keywords
- angle
- section
- guidance
- rlv
- formula
- 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
- 238000000034 method Methods 0.000 title claims abstract description 89
- 230000008569 process Effects 0.000 claims description 24
- 238000012937 correction Methods 0.000 claims description 21
- 206010048669 Terminal state Diseases 0.000 claims description 18
- 230000008859 change Effects 0.000 claims description 17
- 238000013461 design Methods 0.000 claims description 14
- 230000003247 decreasing effect Effects 0.000 claims description 8
- 230000009471 action Effects 0.000 claims description 7
- 230000010354 integration Effects 0.000 claims description 4
- 230000007704 transition Effects 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000007423 decrease Effects 0.000 claims description 3
- 230000004907 flux Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- RCZJVHXVCSKDKB-OYKKKHCWSA-N tert-butyl 2-[(z)-[1-(2-amino-1,3-thiazol-4-yl)-2-(1,3-benzothiazol-2-ylsulfanyl)-2-oxoethylidene]amino]oxy-2-methylpropanoate Chemical compound N=1C2=CC=CC=C2SC=1SC(=O)\C(=N/OC(C)(C)C(=O)OC(C)(C)C)C1=CSC(N)=N1 RCZJVHXVCSKDKB-OYKKKHCWSA-N 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 description 2
- 241001643392 Cyclea Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000021715 photosynthesis, light harvesting Effects 0.000 description 1
- 238000005381 potential energy Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明属于航空航天飞行器制导领域,为使得RLV在受到外界扰动和不确定因素影响的情况下,仍能够获得使轨迹严格满足热流、过载及动压约束的制导指令,最终以恰当的速度在指定的跑道位置着陆。为此,本发明采取的技术方案是,可重复使用运载器全阶段再入返回制导方法,步骤如下:第一部分,初始下降段制导;第二部分:下滑段制导;第三部分:终端能量管理段制导;第四部分:自动着陆段制导:在横向上,保持倾侧角为0;在纵向上,自动着陆段设计为由拟平衡下滑段和指数下滑段组合而成,通过设计航迹角函数确定攻角指令,并通过割线法设计两段切换点处的高度,使RLV能够以预定的速度及航迹角触地。本发明主要应用于航空航天飞行器制导场合。
Description
技术领域
本发明属于航空航天飞行器制导领域,特别是涉及一种可重复使用运载器全阶段再入返回制导方法。
背景技术
可重复使用运载器(RLV,Reusable Launch Vehicle)是指能够自由往返于地球与空间轨道之间执行载荷运送任务的多用途可重复使用飞行器。RLV全阶段再入返回过程从120千米开始,到地面结束,整个再入返回过程飞行包线大,高度跨越太空域,亚轨道和航空域,速度跨越高超声速,超声速和亚声速。单一的再入返回制导方法无法在如此复杂的飞行环境下完成安全返回着陆任务,因此一般将再入返回过程分为三个阶段:再入段(120km-25km)、终端能量管理段(25km-3km)和自动着陆段(3km-地面)。
在再入段,RLV需要在具有极大初始再入动能和势能的条件下,经历高超声速向亚声速的过渡,需要严格满足再入过程的热流约束、动压约束和过载约束,是保证RLV实现安全稳定再入飞行较为关键的一段。根据飞行特点,再入段又可以分为初始下降段(120km-80km)和下滑段(80Km-25km)。在初始下降段,大气密度稀薄,RLV不满足拟平衡滑翔条件,需要设计倾侧角指令,使得RLV能进入拟平衡滑翔条件;在下滑段,RLV跨越大气层,在严格满足各种约束的同时,实现能量耗散,最终精确地进入终端能量管理段初始窗口。
在终端能量管理段,RLV需要耗散再入段结束后RLV所具有的能量,同时调整航向使其对准跑道,保证RLV能够以合适的高度、速度、位置、航向等状态约束条件进入预定的着陆窗口,终端能量管理段区别于另外两个飞行阶段的主要特征是存在航向校正锥来进行航向校正,通过将航天器与航向校正锥相切,设计倾侧角指令,实现能量的耗散并调整航向使RLV对准跑道。
自动着陆段是再入返回飞行过程的最后阶段,它的完成标志着整个再入返回飞行过程的圆满结束,该段主要追求的目标是RLV能够以合理的速度、位置和下沉率在指定的跑道位置安全可靠地着陆。
目前,针对RLV全阶段再入返回制导问题,国内外的学者都只是针对再入返回阶段中的某一阶段进行制导律设计而缺乏全阶段再入返回制导方法,导致存在模型及气动参数不统一,多段之间切换点处存在状态量跳变等问题。因此,针对上述问题,本发明对于全阶段再入返回制导律进行了设计,在统一的运动学模型和气动参数下,分别设计了再入段、终端能量管理段及自动着陆段的制导方法并设置了各个段之间的切换条件,使RLV能克服外界扰动不确定性从而实现实现120千米到地面的RLV全阶段再入返回。
发明内容
为克服现有技术的不足,本发明旨在提出一种RLV全阶段再入返回制导方法,使得RLV从120km开始,在受到外界扰动和不确定因素影响的情况下,仍能够获得使轨迹严格满足热流、过载及动压约束的制导指令,最终以恰当的速度在指定的跑道位置着陆。为此,本发明采取的技术方案是,可重复使用运载器全阶段再入返回制导方法,步骤如下:
第一部分,初始下降段制导:设计固定的倾侧角指令,使可重复使用运载器RLV能够进入拟平衡滑翔条件,并以RLV进入拟平衡滑翔条件状态点作为下滑段初始状态点;
第二部分:下滑段制导:考虑下滑段过程约束以及终端约束,通过设计攻角-速度剖面及基于终点预测偏差的倾侧角校正方法,完成约束条件下的三维预测校正制导律设计,以RLV积分至下滑段终端能量处实际状态点作为终端能量管理段初始状态点;
第三部分:终端能量管理段制导:分开考虑纵向及横向制导律,在纵向上,设计动压-高度剖面以保持动压较小的变化率,并基于动力学方程确定动压与攻角的对应关系;在横向上,基于地面航线轨迹,设计航向角进而确定倾侧角指令,使RLV能够对准跑道,以RLV积分至终端能量管理段末端高度处实际状态点作为自动着陆段初始状态点;
第四部分:自动着陆段制导:在横向上,保持倾侧角为0;在纵向上,自动着陆段设计为由拟平衡下滑段和指数下滑段组合而成,通过设计航迹角函数确定攻角指令,并通过割线法设计两段切换点处的高度,使RLV能够以预定的速度及航迹角触地。
具体实现过程如下:
第一步:在RLV再入返回过程,假设飞行器为无动力飞行的质点,考虑地球为旋转椭球时,忽略再入过程中侧力以及地球自转的影响,并取侧滑角为零,建立统一的RLV再入三自由度运动方程为:
式中中飞行状态r,θ,φ,ν,γ,χ分别表示地心距、经度、纬度、飞行速度、航迹角和航向角,σ表示倾侧角,m表示飞行器质量。重力加速度g=μg/r2,升力L,阻力D表示式如下:
综合考虑到热防护系统以及可操作性需求,在初始下降段及下滑段中,攻角设计为关于速度的剖面
第二步:在初始下降段,通过积分预测,对RLV是否进入拟平衡滑翔条件QEGC进行判定,采用计算量较小的递增策略设计制导指令,由动力学方程式(33)得到
进而得到满足拟平衡滑翔条件的
由式(37)和式(39),建立是否满足拟平衡滑翔条件判定如式(40),从而确定是否由初始下降段过渡到再入段
式中δ>0为常数。初始下降段的制导方法则通过设计倾侧角幅值,取初次迭代σini=0°,倾侧角符号则由航向角和视线角之间的关系得到
sign(σini)=-sign(χ-χsight) (9)
式中χsight为视线角;
与此同时,攻角指令则可由式(36)得到。在上述倾侧角和攻角的作用下,对式(33)进行积分预测,并对式(40)进行判定,若存在满足条件的状态点则结束初始下降段;否则令σini=σini+5°,由式(41)确定符号后重新进行积分预测以及判定,直到满足拟平衡下降条件;
第三步:在下滑段,考虑初始下降段的终端状态作为下滑段的初始状态约束,路径约束包括热流密度约束、动压约束和过载约束如式(42),以及下滑段终端状态约束通过对下滑段终点的积分预测以及基于落点偏差的校正策略,设计制导指令:
路径约束式:
式中:kQ为与飞行器相关的常数,qmax,nmax分别是最大允许的热流、动压和过载,由于时间不是下滑段决定性参数,无法通过约束终端时间从而确定积分预测范围,需要采用一个新的微分量,由于高度在下滑段不是一个单调的变量,无法选取为微分量;若选取速度为微分量,则在运动方程的分母中引入了一项可能等于0的项-D/m-gsinγ,因此,引入下述形式的能量变量作为独立变量:
能量变量e是一个单调递增的变量;此外,它体现了高度与速度的关系,可在运动方程中省去关于速度的表达式,增加了积分预测的速度;
考虑运动方程式(44),设计制导方法时,攻角指令由式(36)给出,因此只需要设计侧倾角的制导指令σgliding,采用预测校正的方法对其进行设计;
在预测环节中,在每个制导周期内,给定初始倾侧角σgliding (0),其符号由航向角和视线角之间的关系得到,并考虑式(44)对能量变量e进行积分,得到下滑段终点处的预测终点经纬度(θgf,φgf);
在校正环节中,设置性能指标函数为
J=J0(σ0)+kc(σ0) (13)
式中:表示当前倾侧角指令作用下终端预测经纬度与指定终端经纬度点的偏差;k为惩罚因子,表示在再入过程中对路径约束的违反,为基于着陆点反向求解的下滑段终端经纬度状态点,反向求解过程的制导指令在后续步骤中给出;
基于所给出的性能指标函数式(45),采用如式(46)的高斯牛顿法对倾侧角进行迭代求解:
式中λi为当次迭代使性能指标函数值减小的最小自然数;
上述迭代求解目的为寻求倾侧角指令σgliding (i),使满足式(47)
J(σgliding)≤Δ (15)
式中Δ为反向求解过程中,考虑着陆时速度和沿航迹方向的地面轨迹存在一定的余量所求得的可行范围。在每个制导周期内,重复上述对于倾侧角的迭代求解步骤,得到满足条件的倾侧角并对式(44)进行积分,最终得到下滑段终点处的终端状态
第四步:在终端能量管理段,考虑下滑段的终端状态作为终端能量管理段的初始状态,分别从纵向上和横向上对制导指令进行设计;在终端能量管理段中,高度h=r-Re为一个单调递减的变量,因此选取高度代替时间作为微分量,由于需要考虑地面投影轨迹变化情况,在运动方程中以RLV在地面坐标系中的位置(y,x)代替经纬度(θ,φ),得到以高度为微分量的运动方程:
在纵向上,为了避免动压以较大的速率变化,设计动压关于高度的多项式
qd(h)=Q0+Q1h+Q2h2+Q3h3 (17)
式中h为海拔高度,h=r-Re,多项式系数Q0,Q1,Q2,Q3由下列条件得到:
基于高度-动压剖面qd(h),在制导过程中,根据相应高度即可得到对应的动压值。由qd=0.5ρν2,及动力学方程式(33),得到
由式(49)可得到
联立式(50)及式(51),可得到阻力系数CD关于高度的函数,
联立式(35)与式(52),可求得相应的攻角指令,完成纵向制导;
在横向上,通过设计RLV地面投影轨迹来设计航向角的变化规律,进而确定倾侧角指令,根据地面投影轨迹的不同,侧向制导可分为4个部分,捕获段、直线飞行段、航向校正段和预着陆段:
1)在捕获段中,RLV调整航向角以对准校正圆柱相切点,倾侧角指令
式中,RAC为捕获段设计转弯半径。
2)在直线飞行段中,RLV保持航向角不变,倾侧角指令
σ=0 (22)
3)在航向校正段中,RLV绕校正圆柱飞行以消耗能量以及对准跑道,倾侧角指令
式中,RHAC为校正圆柱半径。
4)在最终段中,RLV保持航向角不变,切换至自动着陆段,倾侧角指令
σ=0 (24)
第五步:在自动着陆段,考虑终端能量管理段的终端状态作为自动着陆段的初始状态,设计倾侧角指令σlanding保持为0,即不考虑RLV的横向飞行,仅在纵向平面内设计攻角指令使得RLV在着陆场处满足安全降落条件;
自动着陆段制导方法设计可概述为:设计攻角指令αlanding使得RLV在既定地面处满足安全降落条件。具体实现可分为两段:拟平衡下降段及指数下降段。在拟平衡下降段中,RLV保持拟平衡滑翔条件通过设计攻角指令使航迹角γ保持不变;在指数下降段中,按指数变化规律设计航迹角,航迹角逐渐减小到降落航迹角γTD,与此同时,RLV高度逐渐减小直到接触降落跑道。
式中,γ0表示自动着陆段起始点处的航迹角,由式(34)和式(35),可以得到拟平衡下降段的攻角指令
式中L由式(57)给出,在指数下降段中中,设计航迹角的变化趋势
tanγ=tanγTDe-ks (27)
式中:k在下面给出,s表示在x方向上以着陆点处的s=0为基准,与着陆点之间的距离,在自动着陆过程中是一个负数;γTD为预设的触地航迹角;
由式(48)和式(60)得到
则可推导出
将式(64)代入到式(58)中,得到指数下降段的攻角指令αgliding2;
如此自动着陆段的制导问题转变为寻求拟平衡下降段与指数下降段的切换点hF,使两段可以平滑的切换,采用高斯牛顿法对hF进行求解,将性能指标设置为两段积分之间s的误差Δs,最终得到使Δs=0的hF,进而确定两段的攻角指令,完成自动着陆段的制导,实现RLV的安全着陆。
本发明的特点及有益效果是:
本发明设计了一种RLV全阶段再入返回制导方法,使得RLV能够在制导指令作用下,生成严格满足约束的飞行轨迹,最终以恰当的速度在指定的跑道位置着陆。由于本发明所涉及的制导方法不依赖于标称轨迹而是根据RLV当前状态设计制导指令,所以对外界干扰具备良好的鲁棒性。此外,不同于传统方法对再入返回进行分段并对每一段赋予初值及约束条件,本发明实现了对再入返回各个段之间的整合以及模型统一,消除了传统方法中由于各个段分开求解所产生状态量的跳变。最终本发明通过构建RLV再入返回虚拟仿真系统对所提出的制导方法的有效性及鲁棒性进行了验证。
附图说明:
附图1RLV全阶段再入返回制导方法流程图。
附图2正常飞行状态下高度变化曲线。
附图3正常飞行状态下经度变化曲线。
附图4正常飞行状态下纬度变化曲线。
附图5正常飞行状态下速度变化曲线。
附图6正常飞行状态下航迹角变化曲线。
附图7正常飞行状态下航向角变化曲线。
附图8两种情况下RLV全阶段再入返回三维轨迹。
附图9两种情况下高度-速度变化曲线。
具体实施方式
本发明以预测校正制导为主要研究手段,发明了一种可重复使用运载器全阶段再入返回制导方法并通过Matlab软件进行虚拟仿真实验,验证了本方法的有效性。本发明步骤如下:
第一部分,初始下降段制导方法设计:考虑到RLV在再入返回初期120km处空气较为稀薄,且RLV不满足拟平衡滑翔条件,设计固定的倾侧角指令,使RLV能够进入拟平衡滑翔条件,并以RLV进入拟平衡滑翔条件状态点作为下滑段初始状态点。
第二部分:下滑段制导方法设计:考虑下滑段过程约束以及终端约束,通过设计攻角-速度剖面及基于终点预测偏差的倾侧角校正方法,完成约束条件下的三维预测校正制导律设计。以RLV积分至下滑段终端能量处实际状态点作为终端能量管理段初始状态点。
第三部分:终端能量管理段制导方法设计:分开考虑纵向及横向制导律,在纵向上,设计动压-高度剖面以保持动压较小的变化率,并基于动力学方程确定动压与攻角的对应关系;在横向上,基于地面航线轨迹,设计航向角进而确定倾侧角指令,使RLV能够对准跑道。以RLV积分至终端能量管理段末端高度处实际状态点作为自动着陆段初始状态点。
第四部分:自动着陆段制导方法设计:在横向上,保持倾侧角为0;在纵向上,自动着陆段设计为由拟平衡下滑段和指数下滑段组合而成,通过设计航迹角函数确定攻角指令,并通过割线法设计两段切换点处的高度。使RLV能够以预定的速度及航迹角触地。
具体实现过程如下:
第一步:在RLV再入返回过程,假设飞行器为无动力飞行的质点,考虑地球为旋转椭球时,忽略再入过程中侧力以及地球自转的影响,并取侧滑角为零。建立统一的RLV再入三自由度运动方程为
式中中飞行状态r,θ,φ,ν,γ,χ分别表示地心距、经度、纬度、飞行速度、航迹角和航向角。σ表示倾侧角,m表示飞行器质量。重力加速度g=μg/r2,升力L,阻力D表示式如下:
综合考虑到热防护系统以及可操作性需求,在初始下降段及下滑段中,攻角设计为关于速度的剖面
第二步:在初始下降段,通过积分预测,对RLV是否进入拟平衡滑翔条件(QEGC)进行判定,采用计算量较小的递增策略设计制导指令。由动力学方程式(33)可得到
进而得到满足拟平衡滑翔条件的
由式(37)和式(39),建立是否满足拟平衡滑翔条件判定如式(40),从而确定是否由初始下降段过渡到再入段
式中δ>0为常数。初始下降段的制导方法则通过设计倾侧角幅值,取初次迭代σini=0°,倾侧角符号则由航向角和视线角之间的关系得到
sign(σini)=-sign(χ-χsight) (41)
式中χsight为视线角。
与此同时,攻角指令则可由式(36)得到。在上述倾侧角和攻角的作用下,对式(33)进行积分预测,并对式(40)进行判定,若存在满足条件的状态点则结束初始下降段;否则令σini=σini+5°,由式(41)确定符号后重新进行积分预测以及判定,直到满足拟平衡下降条件。
第三步:在下滑段,考虑初始下降段的终端状态作为下滑段的初始状态约束,路径约束包括热流密度约束、动压约束和过载约束如式(42),以及下滑段终端状态约束通过对下滑段终点的积分预测以及基于落点偏差的校正策略,设计制导指令。
路径约束式:
式中:kQ为与飞行器相关的常数,qmax,nmax分别是最大允许的热流、动压和过载。由于时间不是下滑段决定性参数,无法通过约束终端时间从而确定积分预测范围,需要采用一个新的微分量。由于高度在下滑段不是一个单调的变量,无法选取为微分量;若选取速度为微分量,则在运动方程的分母中引入了一项可能等于0的项-D/m-g sin γ。因此,引入下述形式的能量变量作为独立变量
能量变量e是一个单调递增的变量;此外,它体现了高度与速度的关系,可以在运动方程中省去关于速度的表达式,增加了积分预测的速度。
考虑运动方程式(44),设计制导方法时,攻角指令由式(36)给出,因此只需要设计侧倾角的制导指令σgliding,采用预测校正的方法对其进行设计。
在预测环节中,在每个制导周期内,给定初始倾侧角σgliding (0),其符号由航向角和视线角之间的关系得到,并考虑式(44)对能量变量e进行积分,得到下滑段终点处的预测终点经纬度(θgf,φgf)。
在校正环节中,设置性能指标函数为
J=J0(σ0)+kc(σ0) (45)
式中:表示当前倾侧角指令作用下终端预测经纬度与指定终端经纬度点的偏差;k为惩罚因子,表示在再入过程中对路径约束的违反。为基于着陆点反向求解的下滑段终端经纬度状态点,反向求解过程的制导指令在后续步骤中给出。
基于所给出的性能指标函数式(45),采用如式(46)的高斯牛顿法对倾侧角进行迭代求解:
式中λi为当次迭代使性能指标函数值减小的最小自然数。
上述迭代求解目的为寻求倾侧角指令σgliding (i),使满足式(47)
J(σgliding)≤Δ (47)
式中Δ为反向求解过程中,考虑着陆时速度和沿航迹方向的地面轨迹存在一定的余量所求得的可行范围。在每个制导周期内,重复上述对于倾侧角的迭代求解步骤,得到满足条件的倾侧角并对式(44)进行积分,最终得到下滑段终点处的终端状态
在终端能量管理段中,高度h=r-Re为一个单调递减的变量,因此可行选取高度代替时间作为微分量。此外,由于需要考虑地面投影轨迹变化情况,在运动方程中以RLV在地面坐标系中的位置(y,x)代替经纬度(θ,φ),得到以高度为微分量的运动方程
在纵向上,为了避免动压以较大的速率变化,设计动压关于高度的多项式
qd(h)=Q0+Q1h+Q2h2+Q3h3 (49)
式中h为海拔高度,h=r-Re。多项式系数Q0,Q1,Q2,Q3由下列条件得到:
基于高度-动压剖面qd(h),在制导过程中,根据相应高度即可得到对应的动压值。由qd=0.5ρν2,及动力学方程式(33),可得到
由式(49)可得到
联立式(50)及式(51),可得到阻力系数CD关于高度的函数,
联立式(35)与式(52),可求得相应的攻角指令,完成纵向制导。
在横向上,通过设计RLV地面投影轨迹来设计航向角的变化规律,进而确定倾侧角指令,根据地面投影轨迹的不同,侧向制导可分为4个部分,捕获段、直线飞行段、航向校正段和预着陆段:
5)在捕获段中,RLV调整航向角以对准校正圆柱相切点,倾侧角指令
式中,RAC为捕获段设计转弯半径。
6)在直线飞行段中,RLV保持航向角不变,倾侧角指令
σ=0 (54)
7)在航向校正段中,RLV绕校正圆柱飞行以消耗能量以及对准跑道,倾侧角指令
式中,RHAC为校正圆柱半径。
8)在最终段中,RLV保持航向角不变,切换至自动着陆段,倾侧角指令
σ=0 (56)
上述4个部分给出了完整的横向制导指令。结合纵向的攻角指令与横向的倾侧角制导,得到RLV终端能量管理段制导指令,对式(48)进行积分,最终得到终端能量管理终点处的终端状态hALI为终端能量管理段关于高度的终点。
第五步:在自动着陆段,考虑终端能量管理段的终端状态作为自动着陆段的初始状态,设计倾侧角指令σlanding保持为0,即不考虑RLV的横向飞行,仅在纵向平面内设计攻角指令使得RLV在着陆场处满足安全降落条件。
自动着陆段制导方法设计可概述为:设计攻角指令αlanding使得RLV在既定地面处满足安全降落条件。具体实现可分为两段:拟平衡下降段及指数下降段。在拟平衡下降段中,RLV保持拟平衡滑翔条件通过设计攻角指令使航迹角γ保持不变;在指数下降段中,按指数变化规律设计航迹角,航迹角逐渐减小到降落航迹角γTD,与此同时,RLV高度逐渐减小直到接触降落跑道。
式中,γ0表示自动着陆段起始点处的航迹角,由式(34)和式(35),可以得到拟平衡下降段的攻角指令
式中L由式(57)给出,在指数下降段中中,设计航迹角的变化趋势
tanγ=tanγTDe-ks (59)
式中:k在下面给出,s表示在x方向上以着陆点处的s=0为基准,与着陆点之间的距离,在自动着陆过程中是一个负数;γTD为预设的触地航迹角。
由式(48)和式(60)可以得到
则可推导出
将式(64)代入到式(58)中,可得到指数下降段的攻角指令αgliding2。
如此自动着陆段的制导问题转变为寻求拟平衡下降段与指数下降段的切换点hF,使两段可以平滑的切换。采用高斯牛顿法对hF进行求解,将性能指标设置为两段积分之间s的误差Δs,最终得到使Δs=0的hF,进而确定两段的攻角指令,完成自动着陆段的制导,实现RLV的安全着陆。
结合附图对本发明作进一步详述。
在虚拟仿真环境下,仿真参数设定如下:RLV的基础参数与再入约束参数如表1所示,升力系数CL=-0.2070+1.6756α,阻力系数CD=0.0785-0.3529α+2.0400α2。
表1飞行器再入参数
分两种情况进行仿真测试验证,第一种情况考虑正常飞行状态下的RLV再入返回制导,验证本发明所提出制导方法的有效性;第二种情况考虑外界扰动下的RLV再入返回制导,验证本发明所提出制导方法的鲁棒性。
第一种情况的给定指令为:高度从120000km变化到0km、速度从7600m/s变化到200m/s、经度从0rad变化到1.15rad、纬度从0rad变化到0.47rad。图2-图7给出了正常飞行状态下,各状态量相对于能量的变化曲线。由于时间不是全阶段返回的指标之一,因此采取能量作为自变量能更好的体现各状态量的变化情况。从图中可以看出,在正常飞行状态下,本发明提出的制导方法可以有效可靠的引导RLV飞向预期落点。
第二种情况给定指令未变,但是在下滑段初期给RLV加入外部扰动:高度偏差dr=1000m,速度偏差dν=100m/s。图8给出了正常飞行状态下和外界扰动情况下,RLV全阶段再入返回三维轨迹图;图9给出了正常飞行状态下和外界扰动情况下,RLV的高度-速度变化曲线。从图8和图9中可以看出,本发明提出的制导方法在RLV受到外界扰动的情况下仍能引导RLV飞向预期落点,具备良好的鲁棒性。
Claims (4)
1.一种可重复使用运载器全阶段再入返回制导方法,其特征是,步骤如下:
第一部分,初始下降段制导:设计固定的倾侧角指令,使可重复使用运载器RLV能够进入拟平衡滑翔条件,并以RLV进入拟平衡滑翔条件状态点作为下滑段初始状态点;
第二部分:下滑段制导:考虑下滑段过程约束以及终端约束,通过设计攻角-速度剖面及基于终点预测偏差的倾侧角校正方法,完成约束条件下的三维预测校正制导律设计,以RLV积分至下滑段终端能量处实际状态点作为终端能量管理段初始状态点;
第三部分:终端能量管理段制导:分开考虑纵向及横向制导律,在纵向上,设计动压-高度剖面以保持动压较小的变化率,并基于动力学方程确定动压与攻角的对应关系;在横向上,基于地面航线轨迹,设计航向角进而确定倾侧角指令,使RLV能够对准跑道,以RLV积分至终端能量管理段末端高度处实际状态点作为自动着陆段初始状态点;
第四部分:自动着陆段制导:在横向上,保持倾侧角为0;在纵向上,自动着陆段设计为由拟平衡下滑段和指数下滑段组合而成,通过设计航迹角函数确定攻角指令,并通过割线法设计两段切换点处的高度,使RLV能够以预定的速度及航迹角触地。
2.如权利要求1所述的可重复使用运载器全阶段再入返回制导方法,其特征是,具体实现过程如下:
第一步:在RLV再入返回过程,假设飞行器为无动力飞行的质点,考虑地球为旋转椭球时,忽略再入过程中侧力以及地球自转的影响,并取侧滑角为零,建立统一的RLV再入三自由度运动方程为:
式中中飞行状态r,θ,φ,ν,γ,χ分别表示地心距、经度、纬度、飞行速度、航迹角和航向角,σ表示倾侧角,m表示飞行器质量。重力加速度g=μg/r2,升力L,阻力D表示式如下:
综合考虑到热防护系统以及可操作性需求,在初始下降段及下滑段中,攻角设计为关于速度的剖面
第二步:在初始下降段,通过积分预测,对RLV是否进入拟平衡滑翔条件QEGC进行判定,采用计算量较小的递增策略设计制导指令,由动力学方程式(33)得到
进而得到满足拟平衡滑翔条件的
由式(37)和式(39),建立是否满足拟平衡滑翔条件判定如式(40),从而确定是否由初始下降段过渡到再入段
式中δ>0为常数。初始下降段的制导方法则通过设计倾侧角幅值,取初次迭代σini=0°,倾侧角符号则由航向角和视线角之间的关系得到
sign(σini)=-sign(χ-χsight) (9)
式中χsight为视线角;
与此同时,攻角指令则可由式(36)得到。在上述倾侧角和攻角的作用下,对式(33)进行积分预测,并对式(40)进行判定,若存在满足条件的状态点则结束初始下降段;否则令σini=σini+5°,由式(41)确定符号后重新进行积分预测以及判定,直到满足拟平衡下降条件;
第三步:在下滑段,考虑初始下降段的终端状态作为下滑段的初始状态约束,路径约束包括热流密度约束、动压约束和过载约束如式(42),以及下滑段终端状态约束通过对下滑段终点的积分预测以及基于落点偏差的校正策略,设计制导指令:
路径约束式:
式中:kQ为与飞行器相关的常数,qmax,nmax分别是最大允许的热流、动压和过载,由于时间不是下滑段决定性参数,无法通过约束终端时间从而确定积分预测范围,需要采用一个新的微分量,由于高度在下滑段不是一个单调的变量,无法选取为微分量;若选取速度为微分量,则在运动方程的分母中引入了一项可能等于0的项-D/m-gsinγ,因此,引入下述形式的能量变量作为独立变量:
能量变量e是一个单调递增的变量;此外,它体现了高度与速度的关系,可在运动方程中省去关于速度的表达式,增加了积分预测的速度;
考虑运动方程式(44),设计制导方法时,攻角指令由式(36)给出,因此只需要设计侧倾角的制导指令σgliding,采用预测校正的方法对其进行设计;
在预测环节中,在每个制导周期内,给定初始倾侧角σgliding (0),其符号由航向角和视线角之间的关系得到,并考虑式(44)对能量变量e进行积分,得到下滑段终点处的预测终点经纬度(θgf,φgf);
在校正环节中,设置性能指标函数为
J=J0(σ0)+kc(σ0) (13)
式中:表示当前倾侧角指令作用下终端预测经纬度与指定终端经纬度点的偏差;k为惩罚因子,表示在再入过程中对路径约束的违反,为基于着陆点反向求解的下滑段终端经纬度状态点,反向求解过程的制导指令在后续步骤中给出;
基于所给出的性能指标函数式(45),采用如式(46)的高斯牛顿法对倾侧角进行迭代求解:
式中λi为当次迭代使性能指标函数值减小的最小自然数;
上述迭代求解目的为寻求倾侧角指令σgliding (i),使满足式(47)
J(σgliding)≤Δ (15)
式中Δ为反向求解过程中,考虑着陆时速度和沿航迹方向的地面轨迹存在一定的余量所求得的可行范围。在每个制导周期内,重复上述对于倾侧角的迭代求解步骤,得到满足条件的倾侧角并对式(44)进行积分,最终得到下滑段终点处的终端状态
第四步:在终端能量管理段,考虑下滑段的终端状态作为终端能量管理段的初始状态,分别从纵向上和横向上对制导指令进行设计;在终端能量管理段中,高度h=r-Re为一个单调递减的变量,因此选取高度代替时间作为微分量,由于需要考虑地面投影轨迹变化情况,在运动方程中以RLV在地面坐标系中的位置(y,x)代替经纬度(θ,φ),得到以高度为微分量的运动方程:
在纵向上,为了避免动压以较大的速率变化,设计动压关于高度的多项式
qd(h)=Q0+Q1h+Q2h2+Q3h3 (17)
式中h为海拔高度,h=r-Re,多项式系数Q0,Q1,Q2,Q3由下列条件得到:
基于高度-动压剖面qd(h),在制导过程中,根据相应高度即可得到对应的动压值。由qd=0.5ρν2,及动力学方程式(33),得到
由式(49)可得到
联立式(50)及式(51),可得到阻力系数CD关于高度的函数,
联立式(35)与式(52),可求得相应的攻角指令,完成纵向制导;
在横向上,通过设计RLV地面投影轨迹来设计航向角的变化规律,进而确定倾侧角指令,根据地面投影轨迹的不同,侧向制导可分为4个部分,捕获段、直线飞行段、航向校正段和预着陆段:
1)在捕获段中,RLV调整航向角以对准校正圆柱相切点,倾侧角指令
式中,RAC为捕获段设计转弯半径。
2)在直线飞行段中,RLV保持航向角不变,倾侧角指令
σ=0 (22)
3)在航向校正段中,RLV绕校正圆柱飞行以消耗能量以及对准跑道,倾侧角指令
式中,RHAC为校正圆柱半径。
4)在最终段中,RLV保持航向角不变,切换至自动着陆段,倾侧角指令
σ=0 (24)
式中,γ0表示自动着陆段起始点处的航迹角,由式(34)和式(35),可以得到拟平衡下降段的攻角指令:
式中L由式(57)给出,在指数下降段中中,设计航迹角的变化趋势:
tanγ=tanγTDe-ks (27)
式中:k在下面给出,s表示在x方向上以着陆点处的s=0为基准,与着陆点之间的距离,在自动着陆过程中是一个负数;γTD为预设的触地航迹角;
由式(48)和式(60)得到
则可推导出
将式(64)代入到式(58)中,得到指数下降段的攻角指令αgliding2;
如此自动着陆段的制导问题转变为寻求拟平衡下降段与指数下降段的切换点hF,使两段可以平滑的切换,采用高斯牛顿法对hF进行求解,将性能指标设置为两段积分之间s的误差Δs,最终得到使Δs=0的hF,进而确定两段的攻角指令,完成自动着陆段的制导,实现RLV的安全着陆。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911187178.0A CN110908396B (zh) | 2019-11-28 | 2019-11-28 | 可重复使用运载器全阶段再入返回制导方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911187178.0A CN110908396B (zh) | 2019-11-28 | 2019-11-28 | 可重复使用运载器全阶段再入返回制导方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110908396A true CN110908396A (zh) | 2020-03-24 |
CN110908396B CN110908396B (zh) | 2023-02-21 |
Family
ID=69819777
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911187178.0A Active CN110908396B (zh) | 2019-11-28 | 2019-11-28 | 可重复使用运载器全阶段再入返回制导方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110908396B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111221254A (zh) * | 2020-03-26 | 2020-06-02 | 大连理工大学 | 一种基于割线法确定终端时间的航天器姿态调整能量-时间混合最优控制方法 |
CN111813146A (zh) * | 2020-07-01 | 2020-10-23 | 大连理工大学 | 基于bp神经网络预测航程的再入预测-校正制导方法 |
CN112256061A (zh) * | 2020-10-30 | 2021-01-22 | 北京航空航天大学 | 复杂环境及任务约束下的高超声速飞行器再入制导方法 |
CN112498744A (zh) * | 2020-11-12 | 2021-03-16 | 中国航天空气动力技术研究院 | 纵横向松耦合在线轨迹规划方法及电子设备 |
CN112783184A (zh) * | 2020-12-29 | 2021-05-11 | 中国科学院力学研究所广东空天科技研究院 | 一种临近空间垂直投放发射中阶段起始点控制方法及系统 |
CN112817334A (zh) * | 2021-01-18 | 2021-05-18 | 北京临近空间飞行器系统工程研究所 | 一种滑翔飞行器的弹道设计方法、装置及存储介质 |
CN113342044A (zh) * | 2021-06-18 | 2021-09-03 | 北京理工大学 | 一种可重复使用运载器末端能量管理段地面轨迹设计方法 |
CN113467497A (zh) * | 2021-07-08 | 2021-10-01 | 北京星途探索科技有限公司 | 一种基于满足载荷投放点多约束条件下的能量管理制导方法 |
CN113671974A (zh) * | 2021-07-18 | 2021-11-19 | 北京理工大学 | 一种跨域飞行器返回段转弯进场精确制导方法 |
CN113960926A (zh) * | 2021-10-18 | 2022-01-21 | 北京理工大学 | 一种气动捕获制导参数边界的自适应调节方法 |
CN114384935A (zh) * | 2022-01-17 | 2022-04-22 | 北京理工大学 | 一种无人航空载运飞行器多约束气动减速控制方法 |
CN117228007A (zh) * | 2023-11-13 | 2023-12-15 | 北京控制工程研究所 | Sso轨道航天器主动离轨策略的快速确定方法及装置 |
CN117311390A (zh) * | 2023-11-06 | 2023-12-29 | 天津大学 | 一种空天往返飞行器闭环跟踪智能组合制导方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030192984A1 (en) * | 2002-03-18 | 2003-10-16 | Smith Norman L. | System and method for return and landing of launch vehicle booster stage |
CN105159308A (zh) * | 2015-08-28 | 2015-12-16 | 北京控制工程研究所 | 一种可重复使用飞行器着陆段制导与控制律一体化耦合设计方法 |
CN108984907A (zh) * | 2018-07-18 | 2018-12-11 | 哈尔滨工业大学 | 一种基于偏航角条件的迭代制导方法 |
CN109253730A (zh) * | 2018-11-07 | 2019-01-22 | 广州大学 | 可重复使用运载器末端能量管理段的三维轨迹在线规划方法及系统 |
CN110059339A (zh) * | 2019-02-27 | 2019-07-26 | 天津大学 | 基于em-ekf算法的rlv再入段气动参数辨识方法 |
-
2019
- 2019-11-28 CN CN201911187178.0A patent/CN110908396B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030192984A1 (en) * | 2002-03-18 | 2003-10-16 | Smith Norman L. | System and method for return and landing of launch vehicle booster stage |
CN105159308A (zh) * | 2015-08-28 | 2015-12-16 | 北京控制工程研究所 | 一种可重复使用飞行器着陆段制导与控制律一体化耦合设计方法 |
CN108984907A (zh) * | 2018-07-18 | 2018-12-11 | 哈尔滨工业大学 | 一种基于偏航角条件的迭代制导方法 |
CN109253730A (zh) * | 2018-11-07 | 2019-01-22 | 广州大学 | 可重复使用运载器末端能量管理段的三维轨迹在线规划方法及系统 |
CN110059339A (zh) * | 2019-02-27 | 2019-07-26 | 天津大学 | 基于em-ekf算法的rlv再入段气动参数辨识方法 |
Non-Patent Citations (3)
Title |
---|
DONG, Q ET AL.: "Adaptive disturbance observer-based finite-time continuous fault-tolerant control for reentry RLV", 《INTERNATIONAL JOURNAL OF ROBUST AND NONLINEAR CONTROL》 * |
宗群等: "基于 Gauss 伪谱法的临近空间飞行器上升段轨迹优化", 《宇航学报》 * |
汤一华等: "可重复使用运载器再入在线制导方法研究", 《导弹与航天运载技术》 * |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111221254A (zh) * | 2020-03-26 | 2020-06-02 | 大连理工大学 | 一种基于割线法确定终端时间的航天器姿态调整能量-时间混合最优控制方法 |
CN111813146A (zh) * | 2020-07-01 | 2020-10-23 | 大连理工大学 | 基于bp神经网络预测航程的再入预测-校正制导方法 |
CN112256061A (zh) * | 2020-10-30 | 2021-01-22 | 北京航空航天大学 | 复杂环境及任务约束下的高超声速飞行器再入制导方法 |
CN112498744A (zh) * | 2020-11-12 | 2021-03-16 | 中国航天空气动力技术研究院 | 纵横向松耦合在线轨迹规划方法及电子设备 |
CN112783184A (zh) * | 2020-12-29 | 2021-05-11 | 中国科学院力学研究所广东空天科技研究院 | 一种临近空间垂直投放发射中阶段起始点控制方法及系统 |
CN112783184B (zh) * | 2020-12-29 | 2022-11-15 | 广东空天科技研究院 | 一种临近空间垂直投放发射中阶段起始点控制方法及系统 |
CN112817334A (zh) * | 2021-01-18 | 2021-05-18 | 北京临近空间飞行器系统工程研究所 | 一种滑翔飞行器的弹道设计方法、装置及存储介质 |
CN113342044A (zh) * | 2021-06-18 | 2021-09-03 | 北京理工大学 | 一种可重复使用运载器末端能量管理段地面轨迹设计方法 |
CN113342044B (zh) * | 2021-06-18 | 2022-07-05 | 北京理工大学 | 一种可重复使用运载器末端能量管理段地面轨迹设计方法 |
CN113467497A (zh) * | 2021-07-08 | 2021-10-01 | 北京星途探索科技有限公司 | 一种基于满足载荷投放点多约束条件下的能量管理制导方法 |
CN113467497B (zh) * | 2021-07-08 | 2023-09-19 | 北京星途探索科技有限公司 | 一种满足载荷投放点多约束条件下的能量管理制导方法 |
CN113671974A (zh) * | 2021-07-18 | 2021-11-19 | 北京理工大学 | 一种跨域飞行器返回段转弯进场精确制导方法 |
CN113960926A (zh) * | 2021-10-18 | 2022-01-21 | 北京理工大学 | 一种气动捕获制导参数边界的自适应调节方法 |
CN113960926B (zh) * | 2021-10-18 | 2024-04-16 | 北京理工大学 | 一种气动捕获制导参数边界的自适应调节方法 |
CN114384935A (zh) * | 2022-01-17 | 2022-04-22 | 北京理工大学 | 一种无人航空载运飞行器多约束气动减速控制方法 |
CN114384935B (zh) * | 2022-01-17 | 2023-12-08 | 北京理工大学 | 一种无人航空载运飞行器多约束气动减速控制方法 |
CN117311390A (zh) * | 2023-11-06 | 2023-12-29 | 天津大学 | 一种空天往返飞行器闭环跟踪智能组合制导方法 |
CN117311390B (zh) * | 2023-11-06 | 2024-03-19 | 天津大学 | 一种空天往返飞行器闭环跟踪智能组合制导方法 |
CN117228007A (zh) * | 2023-11-13 | 2023-12-15 | 北京控制工程研究所 | Sso轨道航天器主动离轨策略的快速确定方法及装置 |
CN117228007B (zh) * | 2023-11-13 | 2024-01-19 | 北京控制工程研究所 | Sso轨道航天器主动离轨策略的快速确定方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN110908396B (zh) | 2023-02-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110908396B (zh) | 可重复使用运载器全阶段再入返回制导方法 | |
CN110471450B (zh) | 在高度速度剖面内直接规划再入轨迹的方法 | |
CN111306989B (zh) | 一种基于平稳滑翔弹道解析解的高超声速再入制导方法 | |
Lu | Entry guidance: a unified method | |
CN110015446B (zh) | 一种半解析的火星进入制导方法 | |
CN106021628B (zh) | 一种运载火箭垂直返回弹道设计方法 | |
CN104714553B (zh) | 基于几何规划的滑翔飞行器末端能量管理轨迹规划方法 | |
CN109491246B (zh) | 一种基于数值优化算法的自适应救援轨迹规划方法 | |
CN107992074A (zh) | 一种基于飞行路径角规划的再入轨迹设计方法 | |
CN106227972A (zh) | 一种高超声速飞行器平稳滑翔弹道的优化方法 | |
CN104035335A (zh) | 基于高精度纵、横程解析预测方法的平稳滑翔再入制导律 | |
US10656650B2 (en) | Method for guiding and controlling drone using information for controlling camera of drone | |
JP5822676B2 (ja) | 多段式ロケット誘導装置、多段式ロケット誘導プログラム、多段式ロケット誘導方法および多段式ロケット誘導システム | |
CN107861517A (zh) | 基于线性伪谱的跳跃式再入飞行器在线弹道规划制导方法 | |
CN106371312B (zh) | 基于模糊控制器的升力式再入预测-校正制导方法 | |
CN107121929B (zh) | 基于线性协方差模型预测控制的鲁棒再入制导方法 | |
Chuang et al. | Periodic optimal cruise for a hypersonic vehicle with constraints | |
Langelaan | Long distance/duration trajectory optimization for small UAVs | |
CN109240323A (zh) | 一种实时解析构造的空天飞行器再入制导方法 | |
CN111595210A (zh) | 一种大空域高动态火箭子级落区精确垂直回收控制方法 | |
CN108298110A (zh) | 一种两级入轨空天飞行器上升段轨迹及设计方法 | |
JP5822675B2 (ja) | 多段式ロケット誘導装置、多段式ロケット誘導プログラム、多段式ロケット誘導方法および多段式ロケット誘導システム | |
CN112506226B (zh) | 一种基于温度约束条件的长航时无人机飞行航线规划方法 | |
CN107102547B (zh) | 一种基于滑模控制理论的rlv着陆段制导律获取方法 | |
CN115841051A (zh) | 一种基于太空飞行器运行轨道估测的方法及存储介质 |
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 |