CN110733670A - 一种短航程低过载的再入轨迹设计方法 - Google Patents
一种短航程低过载的再入轨迹设计方法 Download PDFInfo
- Publication number
- CN110733670A CN110733670A CN201911070440.3A CN201911070440A CN110733670A CN 110733670 A CN110733670 A CN 110733670A CN 201911070440 A CN201911070440 A CN 201911070440A CN 110733670 A CN110733670 A CN 110733670A
- Authority
- CN
- China
- Prior art keywords
- reentry
- roll angle
- overload
- section
- track
- 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 37
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000010354 integration Effects 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 241001061260 Emmelichthys struhsakeri Species 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
- B64G1/247—Advanced control concepts for autonomous, robotic spacecraft, e.g. by using artificial intelligence, neural networks or autonomous agents
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G99/00—Subject matter not provided for in other groups of this subclass
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Aviation & Aerospace Engineering (AREA)
- General Physics & Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Automation & Control Theory (AREA)
- Pure & Applied Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Data Mining & Analysis (AREA)
- Combustion & Propulsion (AREA)
- Theoretical Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Operations Research (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Robotics (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种短航程低过载的再入轨迹设计方法,考虑对于以接近第二宇宙速度再入大气的飞行器再入过程,引进大升阻比航天器滑翔式再入的概念。首先设定初始下降段常值倾侧角飞行,在滑翔段和末端引入轨迹参数对轨迹进行描述;然后根据轨迹参数与控制指令倾侧角的关系,得到需求的控制指令;为实现低过载,引入过载限制,当过载超出一定范围时,通过适当减小倾侧角来保证过载不超过最大限制的要求;最后通过校正初始下降段倾侧角,实现不同航程任务需求。
Description
技术领域
本发明涉及轨迹规划技术,属于飞行器再入制导领域,可应用于小升阻比的再入航天器再入大气层过程的制导策略中,对于设计的此类再入轨迹,因为其航程短和过载低的特点,对现有载人飞船应急返回过程有一定参考价值。
背景技术
探月返回器以第二宇宙速度再入大气层,面临剧烈的热约束问题,对于载人返回器,对应的过载约束、落点精度要求更高。现有的地球大气再入方式主要有弹道式再入、弹道-升力式再入、升力式再入。对于载人探月返回器主要以弹道- 升力式中的跳跃式再入方法实现返回过程,在升力作用下再次冲出大气层,做一段弹道式飞行后,再一次进入大气层的返回再入。通过跳跃式再入可以使得飞船再入的可设计航程范围变宽、再入最大过载和热流变小,但是对于应急需要以短航程快速返回的情况,存在再入时间长,总吸热量大的问题。现有的以弹道式再入的应急返回方式虽然比较简单,但存在再入过载较大,可能达到7~13g,返回舱着陆点散布范围大的问题。
再入轨迹规划是指按照一定的方法,规划出一条满足再入过程约束和再入终端约束的轨迹。现有轨迹规划技术主要侧重于轨迹规划的最优性和快速性上,即分为离线最优轨迹规划和在线快速轨迹规划技术。在再入最优轨迹规划方面,主要体现在优化方法的发展,以邻近极值法、多重打靶法等为代表的间接法,和以配点法、伪谱法等为代表的直接法。在快速轨迹规划上,主要研究如何提高轨迹规划速度,由离线轨迹规划向在线轨迹规划发展和如何拓展轨迹维度,由二维轨迹向三维轨迹发展的问题。
发明内容
本发明的目的是,针对现有再入方式存在的问题,通过引入大升阻比飞行器中滑翔飞行的概念,通过设定滑翔段轨迹形式,提供一种短航程低过载的再入轨迹设计方法,为载人探月飞船再入返回提供一种应急返回弹道。
为解决上述技术问题,本发明所采用的技术方案是:一种短航程低过载的再入轨迹设计方法,包括以下步骤:
Step1:在地球圆球(r=R0+h)假设,大气密度指数模型(ρ=ρ0e-βh)假设下,以高度为自变量的再入运动方程描述形式如下:
其中v为相对地球的速度大小,为当地速度倾角,ψ为速度方位角,r为地心距,h为飞行器高度,λ、分别为地心经度、地心纬度,σ为倾侧角,m为飞行器质量;ρ为大气密度;Sref为飞行器参考面积;CL、CD分别为飞行器的气动升力系数和气动阻力系数,R0为地球平均半径,ω为地球自转角速度。其余符号下标0代表再入点初始值或标准值,下标f代表再入终端点的值,下标id代表初始下降段参数,gf代表滑翔段参数;为方便进行解析推导,引入如下记号:
上式中,Ωid为满足初始再入状态的积分常数,计算公式如下:
将代入到运动方程的第三式,积分得到速度与高度的关系:
上式中,Πid为满足初始状态的积分常数,计算公式如下:
Πid=v0 2+2gh0;
Step3:在滑翔段和末端都按照相同的轨迹形式飞行,飞行轨迹形式设定如下:
对上式积分得到速度倾角与高度的关系式
其中,Ωgf为满足初始状态的积分常数,计算公式如下:
结合运动方程第三式,
因为在滑翔阶段和末端,阻力加速度的影响远大于引力的影响,因此忽略引力项的影响,并代入速度倾角θ的表达式,适当变形得到
积分上式得到速度
其中Πgf,为满足初始状态的积分常数,计算公式如下:
Step5:当飞行器的过载超过设定的上限时,通过调整倾侧角来实现对再入过载的控制,得到校正后的制导指令σc1。
Step6:若设计出的轨迹不满足任务航程要求,转到Step2,校正初始下降段倾侧角大小,实际航程偏大则增大初始下降段倾侧角大小,反之减小,重复上述步骤,直至满足要求;
Step7:在横向控制上,参考一般的倾侧角翻转策略,输出倾侧角符号;
Step8:设计出的轨迹满足任务航程,过载未超过最大限制,结束。
与现有技术相比,本发明所具有的有益效果为:
(1)本发明提出了一种短航程低过载的再入轨迹设计方法,设计原则形式简单,易于操作,通过引入滑翔的概念,在解析解的基础上,实现再入轨迹设计,运算速度较之于直接数值积分速度更快。
(2)本发明设计出的再入轨迹航程较之于跳跃式再入有更短的航程,航程范围约2000~3000km,能够实现快速返回,同时利用设计的过载控制方案,实现对再入过程最大过载的控制,能够实现再入过程过载n<6g,为飞船的应急返回方案提供参考。
附图说明
图1为一种短航程低过载的再入轨迹设计方法流程图;
图2为2100km航程任务设计轨迹的高度-航程图;
图3为2100km航程任务设计轨迹的过载-时间图。
具体实施方式
下面以成员探测飞行器CEV返回过程为例,说明本发明实施的方式。
再入任务参数设置如下表:
表1再入点初始状态及任务参数
飞行器的质量m为9500kg,参考面积S为23.8m2,飞行器的升阻比约为0.35。
S2:根据解析公式得到再入过程中的速度-高度(V-H),速度倾角-高度 (θ-H)关系,用于快速积分求解航程;
S3:根据轨迹参数μ与倾侧角σ的关系换算得到滑翔段和末段的控制指令;
S5:若设计出的轨迹不满足任务航程(2100km)要求,转到S1,校正初始下降段倾侧角大小,实际航程偏大则增大初始下降段倾侧角大小,反之降低,步长在本例中0.5度,取重复上述步骤,在本例中,最终得到航程 2098km;
S6:在横向参考一般的倾侧角翻转策略,输出倾侧角符号;
S7:设计出的轨迹高度-航程结果如图1,过载-时间如图2所示,满足任务航程,再入过程过载最大5.6g,未超过最大限制(6g),结束。
与传统的方法相比,本发明采用的方法设计出的以第二宇宙速度再入的航天器,其再入轨迹航程更短,再入返回时间在400~450s左右,对于载人飞船应急情况需要快速返回有一定的参考意义,同时在滑翔飞行和过载控制的基础上,实现整个再入过程过载均小于6g,为航天员的安全和舒适性提供保证。
Claims (5)
1.一种短航程低过载的再入轨迹设计方法,其特征在于,该方法包括以下步骤:
1)对于小升阻比的再入飞行器,将从再入点到开伞的过程分为三段,分别为初始下降段、滑翔段和末段;其中初始下降段以常值倾侧角飞行,通过速度和高度、速度倾角和高度的解析公式积分得到初始下降段内的相关状态参数;设定滑翔段和末端都按照相同的轨迹形式飞行,计算公式如下:其中μ为轨迹参数,β=1/7110为大气密度参数,两个阶段的轨迹参数取值不同,滑翔段的轨迹参数μ1小于末段的轨迹参数μ2;θ为速度倾角;
2)根据设定的轨迹形式,推导出滑翔段和末段的速度和高度、速度倾角和高度的解析公式;
4)当飞行器的过载超过设定的上限时,通过调整倾侧角大小实现对再入过载的控制,得到校正后的倾侧角σc1;
5)对步骤2)中的解析公式积分获得实际航程,与任务航程进行比较;若设计出的再入轨迹不满足航程任务需求,则当航程偏大时,增大初始下降段倾侧角幅值;航程偏小时,减小初始下降段倾侧角幅值;若设计出的再入轨迹满足航程任务需求,则输出倾侧角幅值;
6)参考横向倾侧角翻转逻辑,输出倾侧角符号,输出最终满足要求的倾侧角制导指令,结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911070440.3A CN110733670B (zh) | 2019-11-05 | 2019-11-05 | 一种短航程低过载的再入轨迹设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911070440.3A CN110733670B (zh) | 2019-11-05 | 2019-11-05 | 一种短航程低过载的再入轨迹设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110733670A true CN110733670A (zh) | 2020-01-31 |
CN110733670B CN110733670B (zh) | 2020-12-29 |
Family
ID=69272251
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911070440.3A Active CN110733670B (zh) | 2019-11-05 | 2019-11-05 | 一种短航程低过载的再入轨迹设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110733670B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112278336A (zh) * | 2020-11-03 | 2021-01-29 | 北京空间飞行器总体设计部 | 高速再入热流飞行试验验证方法 |
CN112596537A (zh) * | 2020-11-27 | 2021-04-02 | 中国人民解放军国防科技大学 | 用于在线轨迹规划的模型误差补偿方法、系统及存储介质 |
CN112651103A (zh) * | 2020-11-27 | 2021-04-13 | 中国人民解放军国防科技大学 | 用于提高飞行器在线轨迹规划成功率的方法、系统及介质 |
CN113093776A (zh) * | 2021-03-04 | 2021-07-09 | 北京航天飞行控制中心 | 一种航天器的离轨参数确定方法及装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103708045A (zh) * | 2014-01-16 | 2014-04-09 | 中国人民解放军国防科学技术大学 | 一种探月飞船跳跃式再入的在线参数辨识方法 |
CN103921957A (zh) * | 2014-01-16 | 2014-07-16 | 中国人民解放军国防科学技术大学 | 一种探月飞船跳跃式再入的跃起点能量管理方法 |
CN104156538A (zh) * | 2014-08-21 | 2014-11-19 | 天津大学 | 可重复使用运载器再入制导与控制系统性能评估方法 |
CN104648695A (zh) * | 2014-12-18 | 2015-05-27 | 北京控制工程研究所 | 一种基于倾侧角可用性的再入走廊最优规划方法 |
CN105573337A (zh) * | 2015-12-14 | 2016-05-11 | 中国人民解放军国防科学技术大学 | 一种满足再入角和航程约束的离轨制动闭路制导方法 |
CN106628263A (zh) * | 2016-11-23 | 2017-05-10 | 北京电子工程总体研究所 | 一种再入返回航天器推进系统优化配置方法 |
CN108036676A (zh) * | 2017-12-04 | 2018-05-15 | 北京航空航天大学 | 一种基于三维再入弹道解析解的全射向自主再入制导方法 |
-
2019
- 2019-11-05 CN CN201911070440.3A patent/CN110733670B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103708045A (zh) * | 2014-01-16 | 2014-04-09 | 中国人民解放军国防科学技术大学 | 一种探月飞船跳跃式再入的在线参数辨识方法 |
CN103921957A (zh) * | 2014-01-16 | 2014-07-16 | 中国人民解放军国防科学技术大学 | 一种探月飞船跳跃式再入的跃起点能量管理方法 |
CN104156538A (zh) * | 2014-08-21 | 2014-11-19 | 天津大学 | 可重复使用运载器再入制导与控制系统性能评估方法 |
CN104648695A (zh) * | 2014-12-18 | 2015-05-27 | 北京控制工程研究所 | 一种基于倾侧角可用性的再入走廊最优规划方法 |
CN105573337A (zh) * | 2015-12-14 | 2016-05-11 | 中国人民解放军国防科学技术大学 | 一种满足再入角和航程约束的离轨制动闭路制导方法 |
CN106628263A (zh) * | 2016-11-23 | 2017-05-10 | 北京电子工程总体研究所 | 一种再入返回航天器推进系统优化配置方法 |
CN108036676A (zh) * | 2017-12-04 | 2018-05-15 | 北京航空航天大学 | 一种基于三维再入弹道解析解的全射向自主再入制导方法 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112278336A (zh) * | 2020-11-03 | 2021-01-29 | 北京空间飞行器总体设计部 | 高速再入热流飞行试验验证方法 |
CN112278336B (zh) * | 2020-11-03 | 2022-03-22 | 北京空间飞行器总体设计部 | 高速再入热流飞行试验验证方法 |
CN112596537A (zh) * | 2020-11-27 | 2021-04-02 | 中国人民解放军国防科技大学 | 用于在线轨迹规划的模型误差补偿方法、系统及存储介质 |
CN112651103A (zh) * | 2020-11-27 | 2021-04-13 | 中国人民解放军国防科技大学 | 用于提高飞行器在线轨迹规划成功率的方法、系统及介质 |
CN112596537B (zh) * | 2020-11-27 | 2022-03-29 | 中国人民解放军国防科技大学 | 用于在线轨迹规划的模型误差补偿方法、系统及存储介质 |
CN112651103B (zh) * | 2020-11-27 | 2022-10-18 | 中国人民解放军国防科技大学 | 用于提高飞行器在线轨迹规划成功率的方法、系统及介质 |
CN113093776A (zh) * | 2021-03-04 | 2021-07-09 | 北京航天飞行控制中心 | 一种航天器的离轨参数确定方法及装置 |
CN113093776B (zh) * | 2021-03-04 | 2024-02-02 | 北京航天飞行控制中心 | 一种航天器的离轨参数确定方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN110733670B (zh) | 2020-12-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110733670B (zh) | 一种短航程低过载的再入轨迹设计方法 | |
CN106021628B (zh) | 一种运载火箭垂直返回弹道设计方法 | |
de Paiva et al. | A control system development environment for AURORA's semi-autonomous robotic airship | |
CN106371312B (zh) | 基于模糊控制器的升力式再入预测-校正制导方法 | |
Yang | Positioning control for stratospheric satellites subject to dynamics uncertainty and input constraints | |
CN104155990A (zh) | 考虑攻角约束的高超声速飞行器俯仰通道姿态控制方法 | |
Tsukerman et al. | Optimal rendezvous guidance laws with application to civil autonomous aerial refueling | |
CN110334439B (zh) | 一种利用气动力辅助的超轨道速度飞行技术方法 | |
CN112455720B (zh) | 一种空天飞行器气动力辅助变轨设计方法 | |
No et al. | Control and simulation of arbitrary flight trajectory-tracking | |
Scheuermann et al. | Combined lateral and longitudinal control of parafoils using upper-surface canopy spoilers | |
Jonckheere et al. | Propulsion control of crippled aircraft by H/sub/spl infin//model matching | |
Yamasaki et al. | Robust trajectory-tracking method for UAV guidance using proportional navigation | |
Mooij | Characteristic motion of re-entry vehicles | |
Movahhed et al. | Output tracking of a 6-DOF flying wing UAV in longitudinal motion using LQR optimized sliding mode control with integral action | |
Lu et al. | Multi-Phase Vertical Take-Off and Landing Trajectory Optimization with Feasible Initial Guesses | |
Manjarrez et al. | Low level control architecture for automatic takeoff and landing of fixed wing UAV | |
Comer et al. | Total energy flight control architecture optimization for a tilt-wing aircraft | |
Pravitra et al. | Adaptive control for attitude match station-keeping and landing of A fixed-wing UAV onto A maneuvering platform | |
CN114690793B (zh) | 基于滑模控制的可重复使用运载火箭垂直软着陆制导方法 | |
Lee et al. | Matching trajectory optimization and nonlinear tracking control for HALE | |
CN113741551B (zh) | 一种基于代理模型的全过程轨迹优化方法及装置 | |
Labonté | Airplanes at constant speeds on inclined circular trajectories | |
CN114077259B (zh) | 太阳能无人机无动力下滑控制方法 | |
Alatorre et al. | Nonconventional angle-of-attack control strategy for reducing the airspeed during the fixed-wing drone 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 |