CN112270108A - 一种运载火箭椭圆轨道在线规划方法 - Google Patents
一种运载火箭椭圆轨道在线规划方法 Download PDFInfo
- Publication number
- CN112270108A CN112270108A CN202011272311.5A CN202011272311A CN112270108A CN 112270108 A CN112270108 A CN 112270108A CN 202011272311 A CN202011272311 A CN 202011272311A CN 112270108 A CN112270108 A CN 112270108A
- Authority
- CN
- China
- Prior art keywords
- orbit
- constraint
- elliptical orbit
- rocket
- expressed
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Marketing (AREA)
- Geometry (AREA)
- Game Theory and Decision Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明涉及运载火箭椭圆轨道的在线规划方法,是一种火箭在线轨迹规划方法,属航天制导控制领域。能够将椭圆轨道入轨约束这一非线性约束转化为凸约束,从而能够将椭圆轨道规划问题构建为凸规划问题,便于在线求解与实现。通过序列更新入轨点猜想与建立新的椭圆轨道凸规划模型,能够保证入轨精度,精确满足椭圆轨道规划的终端约束。
Description
技术领域
本发明涉及运载火箭椭圆轨道的在线规划方法,是一种火箭在线轨迹规划方法,属航天制导控制领域。
背景技术
对于多级运载火箭,当前一级火箭关机时没有到达预定状态,造成后续火箭无法进入预定轨道时,或当火箭发生非致命故障时,需要火箭能尽量利用剩余燃料,在满足预定轨道倾角、升交点精度、近地点幅角以及近地点高度的情况下,进入轨道的远地点高度最大,从而使得火箭能够到达高度最高的远地点,便于后续任务的进一步挽救。
近年来发展的轨迹规划方法,如伪谱法、间接法等,在求解速度、初始猜想依赖性等方面的缺点,使得其难以在线求解,这样当火箭在飞行过程中遇到非正常状态时,只能进行离线求解,然后通过上传方式将轨迹参数传递到箭上,这种方式需要地面进行大量的人力物力支持,且火箭的飞行区域存在测控盲区和通讯受限情况,因此可能会错过火箭任务的最佳挽救时间,使得火箭能量损失更多,甚至导致挽救任务失败。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提出一种运载火箭椭圆轨道在线规划方法,能够在线自主规划远地点高度最大的椭圆轨道,并同时满足椭圆轨道的其它约束。通过该技术能快速可靠地求解出运载火箭非正常状态下能够到达的椭圆轨道和飞行轨迹,促进火箭或有效载荷后续任务的完成,从而尽最大努力挽救火箭任务,提高火箭的任务适应能力。
本发明的技术解决方案是:
一种运载火箭椭圆轨道在线规划方法,首先根据任务要求,确定椭圆轨道的轨道约束(包含轨道倾角、升交点经度、近地点幅角和近地点高度);其次设计一种凸化方法,将该椭圆轨道约束转化为凸约束;然后根据火箭当前本体参数和状态信息,结合动力学方程,将火箭椭圆轨道规划问题建模为一个序列凸规划问题;最后利用原始-对偶内点法,实时求解该问题,得到远地点高度最大的椭圆轨道和从当前状态到进入椭圆轨道的飞行轨迹和程序角指令。最终通过仿真实验,对该方法进行仿真验证;
详细步骤为:
(1)构建椭圆轨道入轨约束;
将动量矩矢量H=r×v表示为如下形式:
H=κh (1)
其中h表示轨道法向量,κ等价于火箭的终端能量。
其终端轨道动量矩约束可以表示为
r(tf)×v(tf)=κh (2)
采用线性化的方法将其凸化,对动量矩常矢量在原入轨点状态rs与vs处进行泰勒展开,轨道的动量矩约束可以表示为:
根据轨道方程,椭圆轨道的近地点高度约束为:
其中rp表示近地点高度,p表示轨道半通径,e表示轨道偏心率。轨道偏心率与半通径可以通过动量矩矢量与拉普拉斯矢量表示:
e=||L||/μ (5)
p=||H||2/μ (6)
其中,L为拉普拉斯常矢量
由于拉普拉斯矢量的方向即为轨道近地点方向,因此椭圆轨道的近地点高度约束和指向约束可以表示为:
其中rset表示任务指定的近地点。
考虑终端轨道近地点指向约束,令λ=||L||,则其可以表示为
在原入轨点状态rs与vs处进行泰勒展开,终端轨道近地点指向约束表示为:
近地点高度约束表示为:
其中,
(2)构建椭圆轨道凸规划问题;
建立动力学模型
推力加速度的幅值约束可以表示为
||u(t)||≤υ(t)≤Tmaxe-z(t) (13)
Tmax为发动机的最大推力值,对Tmaxe-z(t)关于z0(t)进行泰勒展开并保留线性项,得
以远地点高度最大为性能指标,可以表示为:
maximizeκ (15)
通过以上约束的凸化,得到凸化后的椭圆轨道轨迹规划问题为:
maximizeκ
0≤||u(t)||≤σ(t)≤Tmaxe-z(t)[1-(z(t)-z0(t))]
DHrr(tf)+DHvv(tf)-κh=-Hs+DHrrs+DHvvs
r(t0)=r0,v(t0)=v0,m(t0)=m0
z(tf)≥ln(mf)
利用原始-对偶内点法进行求解得到最优控制指令u。
(3)利用得到的最优控制指令u计算制导指令,用于控制火箭飞行。
T=um=[Tx,Ty,Tz]T
有益效果
能够将椭圆轨道入轨约束这一非线性约束转化为凸约束,从而能够将椭圆轨道规划问题构建为凸规划问题,便于在线求解与实现。通过序列更新入轨点猜想与建立新的椭圆轨道凸规划模型,能够保证入轨精度,精确满足椭圆轨道规划的终端约束。
附图说明
图1为进入最大远地点高度椭圆轨道的火箭在线轨迹规划示意图;
图2为仿真结果示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
一种运载火箭椭圆轨道在线规划方法,首先根据任务要求,确定椭圆轨道的轨道约束(包含轨道倾角、升交点经度、近地点幅角和近地点高度);其次设计一种凸化方法,将该椭圆轨道约束转化为凸约束;然后根据火箭当前本体参数和状态信息,结合动力学方程,将火箭椭圆轨道规划问题建模为一个序列凸规划问题;最后利用原始-对偶内点法,实时求解该问题,得到远地点高度最大的椭圆轨道和从当前状态到进入椭圆轨道的飞行轨迹和程序角指令。最终通过仿真实验,对该方法进行仿真验证;
详细步骤为:
(1)构建椭圆轨道入轨约束;
将动量矩矢量H=r×v表示为如下形式:
H=κh (1)
其中h表示轨道法向量,κ等价于火箭的终端能量。
其终端轨道动量矩约束可以表示为
r(tf)×v(tf)=κh (2)
采用线性化的方法将其凸化,对动量矩常矢量在原入轨点状态rs与vs处进行泰勒展开,轨道的动量矩约束可以表示为:
根据轨道方程,椭圆轨道的近地点高度约束为:
其中rp表示近地点高度,p表示轨道半通径,e表示轨道偏心率。轨道偏心率与半通径可以通过动量矩矢量与拉普拉斯矢量表示:
e=||L||/μ (5)
p=||H||2/μ (6)
其中,L为拉普拉斯常矢量
由于拉普拉斯矢量的方向即为轨道近地点方向,因此椭圆轨道的近地点高度约束和指向约束可以表示为:
其中rset表示任务指定的近地点。
考虑终端轨道近地点指向约束,令λ=||L||,则其可以表示为
在原入轨点状态rs与vs处进行泰勒展开,终端轨道近地点指向约束表示为:
近地点高度约束表示为:
其中,
(2)构建椭圆轨道凸规划问题;
建立动力学模型
推力加速度的幅值约束可以表示为
||u(t)||≤υ(t)≤Tmaxe-z(t) (13)
Tmax为发动机的最大推力值,对Tmaxe-z(t)关于z0(t)进行泰勒展开并保留线性项,得
以远地点高度最大为性能指标,可以表示为:
maximizeκ (15)
通过以上约束的凸化,得到凸化后的椭圆轨道轨迹规划问题为:
maximizeκ
0≤||u(t)||≤σ(t)≤Tmaxe-z(t)[1-(z(t)-z0(t))]
DHrr(tf)+DHvv(tf)-κh=-Hs+DHrrs+DHvvs
r(t0)=r0,v(t0)=v0,m(t0)=m0
z(tf)≥ln(mf)
利用原始-对偶内点法进行求解得到最优控制指令u。
(3)利用得到的最优控制指令u计算制导指令,用于控制火箭飞行。
T=um=[Tx,Ty,Tz]T
实施例
以某型火箭为对象,在如图1中的L1点需要进行椭圆轨道在线规划,规划目标为最大化椭圆的远地点高度,如图1中的rp所示,最终使得火箭在如图1中的L2点进入规划得到的椭圆轨道。
仿真结果:利用本发明的椭圆轨道在线规划方法,得到仿真结果,使得火箭进入<200km×32800km>的椭圆轨道,其中近地点高度为200km,远地点高度为32800km。仿真结果图如图2所示,其中:
1)虚线代表轨迹规划开始时刻所处的停泊轨道;
2)“○”表示的曲线,代表规划得到的当前状态到进入椭圆轨道时刻的飞行轨迹;
3)实线代表规划得到的最终椭圆入轨轨道。
Claims (5)
1.一种运载火箭椭圆轨道在线规划方法,其特征在于:首先根据任务要求,确定椭圆轨道的轨道约束,然后将该椭圆轨道约束转化为凸约束并根据火箭当前本体参数和状态信息,结合动力学方程,将火箭椭圆轨道规划问题建模为序列凸规划问题,最后利用原始-对偶内点法,实时求解,得到远地点高度最大的椭圆轨道和从当前状态到进入椭圆轨道的飞行轨迹和程序角指令。
2.根据权利要求1所述的一种运载火箭椭圆轨道在线规划方法,其特征在于:椭圆轨道的轨道约束包括轨道倾角、升交点经度、近地点幅角和近地点高度。
3.根据权利要求1所述的一种运载火箭椭圆轨道在线规划方法,其特征在于:详细步骤包括:构建椭圆轨道入轨约束;
将动量矩矢量H=r×v表示为如下形式:
H=κh (1)
其中h表示轨道法向量,κ等价于火箭的终端能量;
其终端轨道动量矩约束可以表示为
r(tf)×v(tf)=κh (2)
采用线性化的方法将其凸化,对动量矩常矢量在原入轨点状态rs与vs处进行泰勒展开,轨道的动量矩约束可以表示为:
根据轨道方程,椭圆轨道的近地点高度约束为:
其中rp表示近地点高度,p表示轨道半通径,e表示轨道偏心率;轨道偏心率与半通径可以通过动量矩矢量与拉普拉斯矢量表示:
e=||L||/μ (5)
p=||H||2/μ (6)
其中,L为拉普拉斯常矢量
由于拉普拉斯矢量的方向即为轨道近地点方向,因此椭圆轨道的近地点高度约束和指向约束可以表示为:
其中rset表示任务指定的近地点;
考虑终端轨道近地点指向约束,令λ=||L||,则其可以表示为
在原入轨点状态rs与vs处进行泰勒展开,终端轨道近地点指向约束表示为:
近地点高度约束表示为:
其中,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011272311.5A CN112270108B (zh) | 2020-11-13 | 2020-11-13 | 一种运载火箭椭圆轨道在线规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011272311.5A CN112270108B (zh) | 2020-11-13 | 2020-11-13 | 一种运载火箭椭圆轨道在线规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112270108A true CN112270108A (zh) | 2021-01-26 |
CN112270108B CN112270108B (zh) | 2023-07-18 |
Family
ID=74339939
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011272311.5A Active CN112270108B (zh) | 2020-11-13 | 2020-11-13 | 一种运载火箭椭圆轨道在线规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112270108B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113238582A (zh) * | 2021-07-09 | 2021-08-10 | 北京信息科技大学 | 飞行器的飞行轨迹的规划方法、装置及系统 |
CN115268276A (zh) * | 2022-09-26 | 2022-11-01 | 北京航天自动控制研究所 | 一种基于梯度下降的制导参数在线修正方法及系统 |
CN115469552A (zh) * | 2022-11-01 | 2022-12-13 | 北京航天自动控制研究所 | 一种基于初始协态训练的在线弹道规划方法及装置 |
CN117687306A (zh) * | 2024-02-01 | 2024-03-12 | 北京航空航天大学 | 一种基于模式选择参数的五合一火箭轨迹优化方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030222179A1 (en) * | 2002-06-03 | 2003-12-04 | Peterson Glenn E. | Method and system for controlling the eccentricity of a near-circular orbit |
JP2013107585A (ja) * | 2011-11-24 | 2013-06-06 | Mitsubishi Space Software Kk | 多段式ロケット誘導装置、多段式ロケット誘導プログラム、多段式ロケット誘導方法および多段式ロケット誘導システム |
CN109398762A (zh) * | 2018-10-17 | 2019-03-01 | 湖北航天技术研究院总体设计所 | 一种固体火箭基于椭圆转移轨道的入轨式弹道设计方法 |
CN111428372A (zh) * | 2020-03-29 | 2020-07-17 | 西北工业大学 | 基于凸规划和自适应迭代的火箭动力故障降级入轨制导方法 |
-
2020
- 2020-11-13 CN CN202011272311.5A patent/CN112270108B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030222179A1 (en) * | 2002-06-03 | 2003-12-04 | Peterson Glenn E. | Method and system for controlling the eccentricity of a near-circular orbit |
JP2013107585A (ja) * | 2011-11-24 | 2013-06-06 | Mitsubishi Space Software Kk | 多段式ロケット誘導装置、多段式ロケット誘導プログラム、多段式ロケット誘導方法および多段式ロケット誘導システム |
CN109398762A (zh) * | 2018-10-17 | 2019-03-01 | 湖北航天技术研究院总体设计所 | 一种固体火箭基于椭圆转移轨道的入轨式弹道设计方法 |
CN111428372A (zh) * | 2020-03-29 | 2020-07-17 | 西北工业大学 | 基于凸规划和自适应迭代的火箭动力故障降级入轨制导方法 |
Non-Patent Citations (4)
Title |
---|
XIAOMING CHENG: "Autonomous trajectory planning for space vehicles with a Newton–Kantorovich/convex programming approach", NONLINEAR DYN * |
YUANLI 等: "Online trajectory optimization for power system fault of launch vehicles via convex programming", AEROSPACE SCIENCE AND TECHNOLOGY * |
宋征宇 等: "运载火箭上升段推力下降故障的自主轨迹规划方法", 中国科学 * |
马昊磊 等: "一种基于任务分区的轨道转移飞行器轨迹重规划方法", 航天控制 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113238582A (zh) * | 2021-07-09 | 2021-08-10 | 北京信息科技大学 | 飞行器的飞行轨迹的规划方法、装置及系统 |
CN115268276A (zh) * | 2022-09-26 | 2022-11-01 | 北京航天自动控制研究所 | 一种基于梯度下降的制导参数在线修正方法及系统 |
CN115268276B (zh) * | 2022-09-26 | 2023-02-03 | 北京航天自动控制研究所 | 一种基于梯度下降的制导参数在线修正方法及系统 |
CN115469552A (zh) * | 2022-11-01 | 2022-12-13 | 北京航天自动控制研究所 | 一种基于初始协态训练的在线弹道规划方法及装置 |
CN117687306A (zh) * | 2024-02-01 | 2024-03-12 | 北京航空航天大学 | 一种基于模式选择参数的五合一火箭轨迹优化方法及系统 |
CN117687306B (zh) * | 2024-02-01 | 2024-04-26 | 北京航空航天大学 | 一种基于模式选择参数的五合一火箭轨迹优化方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN112270108B (zh) | 2023-07-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112270108A (zh) | 一种运载火箭椭圆轨道在线规划方法 | |
US20220147820A1 (en) | Optimal rescue orbital elements online decision-making method based on rbfnn for launch vehicles under thrust drop fault | |
Hanson | Advanced guidance and control project for reusable launch vehicles | |
Dukeman | Atmospheric ascent guidance for rocket-powered launch vehicles | |
Flaherty et al. | Operability benefits of airbreathing hypersonic propulsion for flexible access to space | |
CN110329544B (zh) | 一种用于自主快速交会对接的单脉冲制导方法、可读介质 | |
Hanson et al. | Advanced guidance and control methods for reusable launch vehicles: test results | |
CN110053788B (zh) | 一种考虑复杂摄动的星座长期保持控制频次估计方法 | |
Hanson et al. | Ascent guidance comparisons | |
Schierman et al. | In-flight entry trajectory optimization for reusable launch vehicles | |
Ryan et al. | Fundamentals and issues in launch vehicle design | |
CN112580188B (zh) | 一种运载火箭圆轨道在线规划方法 | |
Shen et al. | On-board entry trajectory planning expanded to sub-orbital flight | |
CN109976362B (zh) | 一种使用差分进化算法的火箭上面级弹道重规划方法 | |
CN113741551A (zh) | 一种基于代理模型的全过程轨迹优化方法及装置 | |
Kazmerchuk et al. | Method for optimization of trajectories including gravitational maneuvers of a spacecraft with a solar sail | |
Hanson et al. | Test results for entry guidance methods for reusable launch vehicles | |
Weigel et al. | Dual payload ascent trajectory optimization with a splash-down constraint | |
Yue et al. | Research on agile reentry method for maneuvering vehicle | |
Long et al. | Predictor-corrector strategy based energy suboptimal obstacle avoidance for landing on small bodies | |
Lisy et al. | Post-mission disposal (PMD) validation of FCC's Ka-band NGSO processing round participants | |
De Zaiacomoa et al. | Retalt: Development of key flight dynamics and gnc technologies for reusable launchers | |
CN105631229A (zh) | 一种运载器滑行段测量盲区运行轨迹填补方法 | |
Hamed et al. | Multiple Debris Orbital Collision Avoidance | |
CN117699054A (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 |