CN109974538A - 一种垂直起降可重复使用运载器多终端约束上升段制导方法 - Google Patents
一种垂直起降可重复使用运载器多终端约束上升段制导方法 Download PDFInfo
- Publication number
- CN109974538A CN109974538A CN201910227588.7A CN201910227588A CN109974538A CN 109974538 A CN109974538 A CN 109974538A CN 201910227588 A CN201910227588 A CN 201910227588A CN 109974538 A CN109974538 A CN 109974538A
- Authority
- CN
- China
- Prior art keywords
- control system
- vtol
- launch vehicle
- reusable launch
- guidance
- 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
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F42—AMMUNITION; BLASTING
- F42B—EXPLOSIVE CHARGES, e.g. FOR BLASTING, FIREWORKS, AMMUNITION
- F42B15/00—Self-propelled projectiles or missiles, e.g. rockets; Guided missiles
- F42B15/01—Arrangements thereon for guidance or control
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Aviation & Aerospace Engineering (AREA)
- Combustion & Propulsion (AREA)
- Motorcycle And Bicycle Frame (AREA)
- Feedback Control In General (AREA)
Abstract
本发明提出了一种垂直起降可重复使用运载器多终端约束上升段制导方法,属于运载火箭弹道制导控制技术领域。所述方法包括:步骤一、建立垂直起降可重复使用运载器动力学方程,并根据所述垂直起降可重复使用运载器动力学方程确定过程约束要求;步骤二、根据最优控制原理推导并获取最优控制条件;步骤三、根据制导任务需求,给定终端位置、速度和姿态角要求,并确定其满足的终端约束要求;步骤四、获取终端状态变量和协态变量;步骤五、根据所述满足终端约束的状态变量和协态变量初值,结合最优控制条件即可获得指导指令。本发明有效提高的垂直起降可重复使用运载器上升段制导方法的收敛性、工程实用性和精确性。
Description
技术领域
本发明涉及一种垂直起降可重复使用运载器多终端约束上升段制导方法,属于运载火箭弹道制导控制技术领域。
背景技术
为了降低运载器的发射成本,提升空间快速响应能力,垂直起降可重复使用运载器相关技术的发展倍受关注。可重复使用运载器最关键的环节即为对其的精确返回控制以及回收,要求运载器子级分离后能够精确返回预定着陆回收场。与传统运载火箭不同的是,为了节省推进剂,保证火箭子级返回能力和精度,垂直起降可重复使用运载器上升段终端约束不仅包括主动段结束点的位置和速度状态约束,还包括终端姿态角约束。然而传统的运载火箭上升段一般采用迭代制导方法,其终端姿态角偏差可能达到10~20°,无法适用于垂直起降重复使用运载器上升段制导任务。
发明内容
本发明目的是为了解决现有运载器上升段制导无法约束终端姿态角的问题,提供了一种基于最优控制原理的垂直起降可重复使用运载器多终端约束上升段制导方法。本发明的应用对象为小型垂直起降重复使用运载器。其基本思路为:基于最优控制原理,将垂直起降重复使用运载器的位置、速度、姿态信息视为状态量,根据垂直起降重复使用运载器受力模型建立状态方程,并推导最优控制条件,最后根据终端约束推导得到Hamilton两点边值问题,通过牛顿迭代求解即可得到姿态角指令。本发明所采取的具体技术方案为:
一种垂直起降可重复使用运载器多终端约束上升段制导方法,所述垂直起降可重复使用运载器多终端约束上升段制导方法包括:
步骤一、建立垂直起降可重复使用运载器动力学方程,并根据所述垂直起降可重复使用运载器动力学方程确定过程约束要求;
步骤二、根据最优控制原理推导并获取最优控制条件;
步骤三、根据垂直起降可重复使用运载器多终端约束上升段的制导任务,给定终端位置、速度和姿态角要求并确定所述终端位置、速度和姿态角要求满足的终端约束要求;
步骤四、在给定初值的条件下对所述垂直起降可重复使用运载器动力学方程和最优控制条件进行积分,获得终端状态变量和协态变量;其中,根据工程经验给定初值;利用牛顿迭代法,对所述终端状态变量和协态变量进行求解能够使终端状态变量和协态变量满足终端约束的状态变量和协态变量初值;
步骤五、根据所述满足终端约束的状态变量和协态变量初值,结合最优控制条件即可获得指导指令。
进一步地,步骤一所述垂直起降可重复使用运载器动力学方程为:
其中,r,V,u,a分别为垂直起降可重复使用运载器的位置、速度、姿态和姿态变化率矢量,和分别表示位置、速度和姿态的导数;T为发动机推力,m为发动机质量;
所述确定过程约束要求:保证|u|≡1,则有约束要求为:
ua=0 (2)。
进一步地,步骤二所述的根据最优控制原理推导并获取最优控制条件的过程包括:
第一步、根据最优控制原理,确定性能函数J为:
J=min tf (3)
其中,t为发射时间,下角标“f”表示终端值;
第二步、确定哈密尔顿函数H为:
其中,pr,pv,pu,μ为协态变量;
第三步、根据最优控制原理,确定最优控制条件为:
第四步、根据第三步获得的最优控制条件,则确定:
且有
则
a=-pu-μu (11)
根据约束
ua=0 (12)
有
进一步地,步骤三所述所述终端位置、速度和姿态角要求满足的终端约束要求为:
Hf+1=0 (14)。
本发明有益效果:
本发明提出的一种基于最优控制理论的垂直起降可重复使用运载器多终端约束上升段制导方法。该方法通过对垂直起降可重复使用运载器最优控制模型的推导,可实现对垂直起降可重复使用运载器上升段位置、速度、姿态角多终端约束精确制导,同时保证在初值不精确、终端约束精度要求高的情况下,对制导指令的高精度、高效率求解。本发明有效提高的垂直起降可重复使用运载器上升段制导方法的收敛性、工程实用性和精确性,在垂直起降可重复使用运载器制导领域具有广阔的应用前景。
具体实施方式
下面结合具体实施例对本发明做进一步说明,但本发明不受实施例的限制。
实施例1:
一种垂直起降可重复使用运载器多终端约束上升段制导方法,包括以下步骤:
步骤一:建立垂直起降可重复使用运载器动力学方程,并给出过程约束要求:
式中,r,V,u,a分别为垂直起降可重复使用运载器的位置、速度、姿态,姿态变化率矢量,本专利中上标“·”表示导数,T为发动机推力,m为发动机质量。
为保证|u|≡1,则有约束
ua=0 (2)
步骤二:推导最优控制条件:
根据最优控制原理,性能指标函数J取
J=min tf (3)
t为发射时间,本专利中下标“f”表示终端值。
哈密尔顿函数H
式中,pr,pv,pu,μ为协态变量。
根据最优控制原理,最优控制条件为
根据最优控制条件,有:
且有
则
a=-pu-μu (11)
根据约束
ua=0 (12)
有
步骤三:给出终端约束要求;
根据任务需求,给定终端位置、速度和姿态角要求且同时还需满足
Hf+1=0 (14)
步骤四:在给定初值的条件下,积分公式(1)、(5)~(7),即可得到终端状态变量、协态变量。应用牛顿迭代法,求解能够使终端状态变量、协态变量满足终端约束的状态变量、协态变量初值即可;
步骤五:得到状态变量、协态变量初值,并根据最优控制条件公式(11),即可得到制导指令。
本发明提出的一种基于最优控制理论的垂直起降可重复使用运载器多终端约束上升段制导方法。该方法通过对垂直起降可重复使用运载器最优控制模型的推导,可实现对垂直起降可重复使用运载器上升段位置、速度、姿态角多终端约束精确制导,同时保证在初值不精确、终端约束精度要求高的情况下,对制导指令的高精度、高效率求解。本发明有效提高的垂直起降可重复使用运载器上升段制导方法的收敛性、工程实用性和精确性,在垂直起降可重复使用运载器制导领域具有广阔的应用前景。
虽然本发明已以较佳的实施例公开如上,但其并非用以限定本发明,任何熟悉此技术的人,在不脱离本发明的精神和范围内,都可以做各种改动和修饰,因此本发明的保护范围应该以权利要求书所界定的为准。
Claims (4)
1.一种垂直起降可重复使用运载器多终端约束上升段制导方法,其特征在于,所述垂直起降可重复使用运载器多终端约束上升段制导方法包括:
步骤一、建立垂直起降可重复使用运载器动力学方程,并根据所述垂直起降可重复使用运载器动力学方程确定过程约束要求;
步骤二、根据最优控制原理推导并获取最优控制条件;
步骤三、根据垂直起降可重复使用运载器多终端约束上升段的制导任务,给定终端位置、速度和姿态角要求并确定所述终端位置、速度和姿态角要求满足的终端约束要求;
步骤四、在给定初值的条件下对所述垂直起降可重复使用运载器动力学方程和最优控制条件进行积分,获得终端状态变量和协态变量;其中,给定初值;利用牛顿迭代法,对所述终端状态变量和协态变量进行求解能够使终端状态变量和协态变量满足终端约束的状态变量和协态变量初值;
步骤五、根据所述满足终端约束的状态变量和协态变量初值,结合最优控制条件即可获得指导指令。
2.根据权利要求1所述垂直起降可重复使用运载器多终端约束上升段制导方法,其特征在于,步骤一所述垂直起降可重复使用运载器动力学方程为:
其中,r,V,u,a分别为垂直起降可重复使用运载器的位置、速度、姿态和姿态变化率矢量,和分别表示位置、速度和姿态的导数;T为发动机推力,m为发动机质量;
所述确定过程约束要求:保证|u|≡1,则有约束要求为:
ua=0 (2)。
3.根据权利要求1所述垂直起降可重复使用运载器多终端约束上升段制导方法,其特征在于,步骤二所述的根据最优控制原理推导并获取最优控制条件的过程包括:
第一步、根据最优控制原理,确定性能函数J为:
J=min tf (3)
其中,t为发射时间,下角标“f”表示终端值;
第二步、确定哈密尔顿函数H为:
其中,pr,pv,pu,μ为协态变量;
第三步、根据最优控制原理,确定最优控制条件为:
第四步、根据第三步获得的最优控制条件,则确定:
且有
则
a=-pu-μu (11)
根据约束
ua=0 (12)
有
4.根据权利要求1所述垂直起降可重复使用运载器多终端约束上升段制导方法,其特征在于,步骤三所述所述终端位置、速度和姿态角要求满足的终端约束要求为:
Hf+1=0 (14)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910227588.7A CN109974538B (zh) | 2019-03-25 | 2019-03-25 | 一种垂直起降可重复使用运载器多终端约束上升段制导方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910227588.7A CN109974538B (zh) | 2019-03-25 | 2019-03-25 | 一种垂直起降可重复使用运载器多终端约束上升段制导方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109974538A true CN109974538A (zh) | 2019-07-05 |
CN109974538B CN109974538B (zh) | 2020-09-04 |
Family
ID=67080412
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910227588.7A Active CN109974538B (zh) | 2019-03-25 | 2019-03-25 | 一种垂直起降可重复使用运载器多终端约束上升段制导方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109974538B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110989644A (zh) * | 2019-11-29 | 2020-04-10 | 上海宇航系统工程研究所 | 一种考虑目标点多终端约束的飞行器轨迹规划方法 |
CN111580555A (zh) * | 2020-05-13 | 2020-08-25 | 北京控制工程研究所 | 一种高超声速飞行器上升段分段自适应预测校正制导方法 |
CN113758383A (zh) * | 2021-08-17 | 2021-12-07 | 航天科工火箭技术有限公司 | 一种用于验证垂直起降技术的可重复使用火箭及验证方法 |
CN114459297A (zh) * | 2022-03-04 | 2022-05-10 | 航天科工火箭技术有限公司 | 可重复使用运载火箭的地面回收平台 |
CN115421388A (zh) * | 2022-09-23 | 2022-12-02 | 南京理工大学 | 一种基于凸优化的远程导弹末级多姿态约束在线轨迹规划方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE2012793A1 (de) * | 1970-03-18 | 1971-10-07 | Industrieanlagen Betriebsges | Verfahren und Anordnung zum senkrechten Start von Flugkörpern aus bewaldetem oder be bautem Gelände heraus mit anschließender schneller Umlenkung in andere Flugnchtun gen |
JPH0278761A (ja) * | 1988-09-14 | 1990-03-19 | Mitsubishi Heavy Ind Ltd | 飛しよう体ロケットエンジンの制御装置 |
US20180222580A1 (en) * | 2012-12-07 | 2018-08-09 | Delorean Aerospace, Llc | Vertical takeoff and landing aircraft |
US20190061936A1 (en) * | 2017-08-22 | 2019-02-28 | United States Of America As Represented By The Administration Of Nasa | Vertical take off and landing (vtol) aircraft with vectored thrust having continuously variable pitch attitude in hover |
CN109407688A (zh) * | 2018-12-13 | 2019-03-01 | 北京航天自动控制研究所 | 一种垂直起降火箭在线轨迹规划的质心运动解耦方法 |
CN109470252A (zh) * | 2018-10-23 | 2019-03-15 | 哈尔滨工业大学 | 一种基于凸优化的垂直起降重复使用运载器快速轨迹优化方法 |
CN109506517A (zh) * | 2018-11-21 | 2019-03-22 | 中国人民解放军空军工程大学 | 一种带约束的中制导弹道优化方法 |
-
2019
- 2019-03-25 CN CN201910227588.7A patent/CN109974538B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE2012793A1 (de) * | 1970-03-18 | 1971-10-07 | Industrieanlagen Betriebsges | Verfahren und Anordnung zum senkrechten Start von Flugkörpern aus bewaldetem oder be bautem Gelände heraus mit anschließender schneller Umlenkung in andere Flugnchtun gen |
JPH0278761A (ja) * | 1988-09-14 | 1990-03-19 | Mitsubishi Heavy Ind Ltd | 飛しよう体ロケットエンジンの制御装置 |
US20180222580A1 (en) * | 2012-12-07 | 2018-08-09 | Delorean Aerospace, Llc | Vertical takeoff and landing aircraft |
US20190061936A1 (en) * | 2017-08-22 | 2019-02-28 | United States Of America As Represented By The Administration Of Nasa | Vertical take off and landing (vtol) aircraft with vectored thrust having continuously variable pitch attitude in hover |
CN109470252A (zh) * | 2018-10-23 | 2019-03-15 | 哈尔滨工业大学 | 一种基于凸优化的垂直起降重复使用运载器快速轨迹优化方法 |
CN109506517A (zh) * | 2018-11-21 | 2019-03-22 | 中国人民解放军空军工程大学 | 一种带约束的中制导弹道优化方法 |
CN109407688A (zh) * | 2018-12-13 | 2019-03-01 | 北京航天自动控制研究所 | 一种垂直起降火箭在线轨迹规划的质心运动解耦方法 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110989644A (zh) * | 2019-11-29 | 2020-04-10 | 上海宇航系统工程研究所 | 一种考虑目标点多终端约束的飞行器轨迹规划方法 |
CN110989644B (zh) * | 2019-11-29 | 2021-04-23 | 上海宇航系统工程研究所 | 一种考虑目标点多终端约束的飞行器轨迹规划方法 |
CN111580555A (zh) * | 2020-05-13 | 2020-08-25 | 北京控制工程研究所 | 一种高超声速飞行器上升段分段自适应预测校正制导方法 |
CN113758383A (zh) * | 2021-08-17 | 2021-12-07 | 航天科工火箭技术有限公司 | 一种用于验证垂直起降技术的可重复使用火箭及验证方法 |
CN113758383B (zh) * | 2021-08-17 | 2023-03-03 | 航天科工火箭技术有限公司 | 一种用于验证垂直起降技术的可重复使用火箭及验证方法 |
CN114459297A (zh) * | 2022-03-04 | 2022-05-10 | 航天科工火箭技术有限公司 | 可重复使用运载火箭的地面回收平台 |
CN114459297B (zh) * | 2022-03-04 | 2024-06-04 | 航天科工火箭技术有限公司 | 可重复使用运载火箭的地面回收平台 |
CN115421388A (zh) * | 2022-09-23 | 2022-12-02 | 南京理工大学 | 一种基于凸优化的远程导弹末级多姿态约束在线轨迹规划方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109974538B (zh) | 2020-09-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109974538A (zh) | 一种垂直起降可重复使用运载器多终端约束上升段制导方法 | |
CN107544262B (zh) | 一种运载火箭自适应精确回收控制方法 | |
CN108984907A (zh) | 一种基于偏航角条件的迭代制导方法 | |
CN106250625B (zh) | 一种航天器迭代制导的优化方法 | |
CN108803649A (zh) | 一种垂直起降重复使用运载器自抗扰滑模控制方法 | |
CN110254753B (zh) | 一种地球静止轨道卫星电推力器及其布局优化方法 | |
CN111591470B (zh) | 一种适应推力可调模式的飞行器精确软着陆闭环制导方法 | |
CN110716566A (zh) | 一种欠驱动无人艇的有限时间轨迹跟踪控制方法 | |
CN110471456A (zh) | 高超声速飞行器俯冲段制导、姿控、变形一体化控制方法 | |
CN107450323B (zh) | 高超声速飞行器再入段神经网络复合学习控制方法 | |
CN102880055A (zh) | 基于预测模型的高超声速飞行器神经网络控制方法 | |
CN106371312A (zh) | 基于模糊控制器的升力式再入预测‑校正制导方法 | |
CN107092725B (zh) | 一种基于闭环仿真的运载器分布载荷优化设计方法 | |
CN102288177B (zh) | 一种基于角速率输出的捷联系统速度解算方法 | |
Leutenegger | Unmanned solar airplanes: Design and algorithms for efficient and robust autonomous operation | |
CN115952384A (zh) | 一种运载火箭掉头过程坐标系转换方法及其控制仿真应用 | |
CN112362289A (zh) | 飞行器分体试验方法、装置、计算机设备及可读存储介质 | |
CN113138608B (zh) | 一种使用扰动观测器和非线性速度观测器的四旋翼无人机视觉伺服控制方法 | |
Wei | Autonomous control system for the quadrotor unmanned aerial vehicle | |
CN111781820A (zh) | 基于重心动力学的航模pid算法控制方法 | |
Zheng et al. | Smooth sub-optimal trajectory generation for transition maneuvers | |
Chen et al. | A unified aerodynamic-influence-coefficient approach for aeroelastic/aeroservoelastic and MDO applications | |
CN104634183A (zh) | 一种基于升阻比实时估计的自适应制导方法 | |
CN117104542A (zh) | 一种用于垂直回收和重复使用技术的验证方法及平台 | |
CN114715437A (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 |