CN114676574B - 基于滑模控制的可重复使用运载火箭燃料计算方法 - Google Patents
基于滑模控制的可重复使用运载火箭燃料计算方法 Download PDFInfo
- Publication number
- CN114676574B CN114676574B CN202210314878.7A CN202210314878A CN114676574B CN 114676574 B CN114676574 B CN 114676574B CN 202210314878 A CN202210314878 A CN 202210314878A CN 114676574 B CN114676574 B CN 114676574B
- Authority
- CN
- China
- Prior art keywords
- rocket
- sliding mode
- fuel
- mode control
- height
- 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.)
- Active
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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Feedback Control In General (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种基于滑模控制的可重复使用运载火箭燃料计算方法,所示方法包括如下步骤:步骤一、建立以伪高度变量τ为自变量的数学模型;步骤二、制定基于滑模控制理论的RLV垂直软着陆的速度控制策略;步骤三、基于滑模控制理论的RLV垂直软着陆的速度控制策略计算火箭回收过程所需燃料。本发明提出了一类基于滑模控制理论的RLV垂直软着陆的速度控制策略,并针对该策略提供了一种燃料计算方法。该方法可以提前求得火箭回收过程所需的燃料消耗量,解决了燃料约束问题。该燃料计算方法还提供了火箭减速过程中动能、势能和外部扰动三个部分的燃料消耗量的影响因子。通过调节相应的影响因子,可以控制火箭回收过程所需的燃料量。
Description
技术领域
本发明涉及一种燃料计算方法,具体涉及一种基于滑模控制的可重复使用运载火箭燃料计算方法。
背景技术
可重复使用运载器(Reusable Launch Vehicle,RLV)成为当下研究的热门。火箭回收的方式大体上可以分为伞降回收、伞降加气囊回收、有翼水平回收、垂直回收以及旋翼方式回收等。其中,垂直回收方式最具有可行性和应用前景,已经被成功的应用于猎鹰9等火箭。
由于火箭返回过程空间跨度大,各飞行阶段任务不同,飞行环境复杂多变,存在较强的不确定性干扰,对着陆要求苛刻。除此之外,还需要满足动压、热流、过载、可用燃料等约束条件。制导技术一直是RLV研究的重点和难点。滑模控制的动态响应速度快、算法简单、物理实现容易、对参数摄动和外界干扰不敏感、鲁棒性和适应性好,适用于火箭垂直着陆。垂直回收的重点之一是实现速度可控,实现软着陆。现有针对RLV的滑模制导方法研究较少,且通常只针对无动力滑行段或作为跟踪制导的参考轨迹,并未考虑速度与燃料的影响,无法完整的实现火箭垂直回收过程,导致其实用性较小。
发明内容
本发明的目的是提出一类基于滑模控制理论的RLV垂直软着陆的速度控制策略,并针对该策略提供了一种基于滑模控制的可重复使用运载火箭燃料计算方法。该方法可以计算速度控制过程火箭需要消耗的燃料量,且该消耗量可以调节。
本发明的目的是通过以下技术方案实现的:
一种基于滑模控制的可重复使用运载火箭燃料计算方法,包括如下步骤:
步骤一、建立以伪高度变量τ为自变量的数学模型:
其中,v为速度,θ为弹道倾角,m为导弹质量,g为地球重力加速度,P为发动机推力,δ为发动机偏转角,α为攻角,Isg为发动机比冲,Fe=[Fex,Fey]为离心惯性力,Fk=[Fkx,Fky]为哥氏惯性力,dw1和dw2为外部扰动,D为气动阻力,L为气动升力,右上带“'”表示关于τ的一阶导数;τ=y0-y,y0为火箭的初始高度,y为火箭在当前时刻的高度,时间变量t与伪高度变量τ的关系为:
步骤二、制定基于滑模控制理论的RLV垂直软着陆的速度控制策略
火箭回收可分为调姿段、动力减速段、气动减速段和垂直着陆段四个阶段,其中,动力减速段和垂直着陆段需要利用发动机推力控制火箭速度,由于动力减速段空气稀薄,垂直着陆段速度较小,即在控制速度时,气动力较小,可将其作为干扰项,又因为攻角和发动机偏转角较小,可以近似为零,则公式(3)简化为:
火箭垂直软着陆的关键在于速度控制,设如下状态误差变量e为:
e=v-vf (7);
其中,vf为常值期望速度,由公式(6)和公式(7)可得如下一阶非线性系统:
针对上述一阶非线性系统设计如下任意一阶滑模函数S:
滑模函数S关于伪高度变量τ的一阶导数为:
结合滑模控制理论,可得发动机控制推力P为:
步骤三、基于滑模控制理论的RLV垂直软着陆的速度控制策略计算火箭回收过程所需燃料
由公式(9)可得e的解析式为:
由公式(7)和公式(12)可得速度v的解析式为:
由公式(5)、公式(11)、公式(12)和公式(13)可得如下以伪高度τ为自变量的一阶齐次线性微分方程:
定义如下三个变量因子:
其中,D1为势能因子,D2为动能因子,D3为抗扰因子,τf为任意终端伪高度,τf=y0-yf,y0为火箭的初始高度,yf为终端高度;
若火箭的初始质量为m0且已知,则可得火箭的理论终端质量mf为:
如果已知火箭不含燃料的空载质量为mf,则可得火箭理论所需初始质量m0:
由公式(18)、公式(19)可得理论所需燃料为Δm=m0-mf。
相比于现有技术,本发明具有如下优点:
1、本发明提出了一类基于滑模控制理论的RLV垂直软着陆的速度控制策略,并针对该策略提供了一种燃料计算方法。该方法可以提前求得火箭回收过程所需的燃料消耗量,解决了燃料约束问题。
2、该燃料计算方法还提供了火箭减速过程中动能、势能和外部扰动三个部分的燃料消耗量的影响因子。通过调节相应的影响因子,可以控制火箭回收过程所需的燃料量。
附图说明
图1为燃料计算流程图;
图2为应用实例1的速度变化曲线;
图3为应用实例1的质量变化曲线;
图4为应用实例1的发动机推力变化曲线;
图5为应用实例1的轨迹变化曲线;
图6为应用实例2的速度变化曲线;
图7为应用实例2的质量变化曲线;
图8为应用实例2的发动机推力变化曲线;
图9为应用实例2的轨迹变化曲线。
具体实施方式
下面结合附图对本发明的技术方案作进一步的说明,但并不局限于此,凡是对本发明技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,均应涵盖在本发明的保护范围中。
本发明提供了一种基于滑模控制的可重复使用运载火箭燃料计算方法,如图1所示,所述方法包括如下步骤:
步骤一:建立伪高度动力学模型
火箭回收的二维运动学方程和动力学方程为:
其中,v为速度,θ为弹道倾角,m为导弹质量,g为地球重力加速度,P为发动机推力,δ为发动机偏转角,α为攻角,Isg为发动机比冲,Fe=[Fex,Fey]为离心惯性力,Fk=[Fkx,Fky]为哥氏惯性力,dw1和dw2为外部扰动,D为气动阻力,L为气动升力。
由于在制导的过程中,时间难以精确的控制。本发明用伪高度变量τ替代时间变量t进行制导律设。令伪高度变量τ为:
τ=y0-y (25);
其中,y0为火箭在任意阶段的初始高度,y为火箭在当前时刻的高度。时间变量t与伪高度变量τ的关系为:
利用公式(26),可将公式(20)~公式(24)转化为如下以伪高度变量τ为自变量的数学模型:
其中,右上带“'”均表示关于新的自变量τ的一阶导数。
步骤二、制定基于滑模控制理论的RLV垂直软着陆的速度控制策略
火箭回收可分为调姿段、动力减速段、气动减速段和垂直着陆段四个阶段。其中,动力减速段和垂直着陆段需要利用发动机推力控制火箭速度。由于动力减速段空气稀薄,垂直着陆段速度较小,即在控制速度时,气动力较小,可将其作为干扰项。又因为攻角和发动机偏转角较小,可以近似为零。公式(29)可简化为:
其中,d为复合干扰。
火箭垂直软着陆的关键在于速度控制,设如下状态误差变量e为:
e=v-vf (33);
其中,vf为常值期望速度。由公式(32)和(33)可得如下一阶非线性系统:
针对公式(34)所示一阶非线性系统设计如下任意一阶滑模函数:
滑模函数S关于伪高度变量τ的一阶导数为:
结合滑模控制理论,可得发动机控制推力P为:
步骤三、基于滑模控制理论的RLV垂直软着陆的速度控制策略计算火箭回收过程所需燃料
由公式(35)可得e的解析式为:
由公式(33)和公式(38)可得速度v的解析式为:
由公式(31)、公式(37)、公式(38)和公式(39)可得如下以伪高度τ为自变量的一阶齐次线性微分方程:
定义如下三个变量因子:
其中,D1为势能因子,其用于抵消火箭的重力势能;D2为动能因子,其用于抵消火箭的动能;D3为抗扰因子,其用于克服不确定扰动;τf为任意终端伪高度,τf=y0-yf,y0为火箭的初始高度,yf为终端高度。
若火箭的初始质量为m0且已知,则可得火箭的理论终端质量mf为:
同理,如果已知火箭不含燃料的空载质量为mf,则可得火箭理论所需初始质量:
由公式(44)或公式(45)可得理论所需燃料为Δm=m0-mf,即在规定的高度内实现期望的速度控制,火箭因保证剩余燃料不小于Δm。且由公式(41)、公式(42)和公式(43)所示的势能因子D1,动能因子D2和抗扰因子D3可知,Δm可以通过设置不同的控制律和控制参数来满足实际需求。
应用实例1:
其中,a1、p1和q1为常增益正系数,τ0为该阶段初始伪高度,τf1终端伪高度,C1为由初始条件确定的常数:C1=-e1(τ0)-a1。
公式(35)所示滑模函数关于伪高度变量τ的一阶导数为:
结合滑模控制方法,可得发动机控制推力P为:
表1动力减速阶段初始参数
结合表1所示参数和公式(44)可求得火箭在终点高度yf1处的理论质量为mf=21938kg,即理论所需消耗的燃料为Δm=8062kg。由仿真结果得,火箭在终点高度yf1处的实际速度为999.98m/s,实际质量为21972kg,即实际消耗的燃料为Δm=8028kg。理论计算多出的34kg燃料为用于抵抗外部干扰预留的燃料,通常K1值越大,预留燃料越多,其用于克服外部干扰的能力越强。其他仿真结果如图2~图5所示。
应用实例2:
垂直降落段伪高度变化范围为τf2~τf3,其中τf2=0,τf3=yf2-yf3,yf2为垂直降落段起始高度,yf3为垂直降落段着陆高度。垂直降落阶段速度较小,攻角α≈0,故气动力较小。该阶段主要通过控制量发动机推力P和发动机偏转角δ实现软着陆。将气动力和惯性力当作不确定干扰,公式(29)和公式(30)可简化为:
由于垂直降落阶段不仅需要考虑速度控制,还需要考虑着陆时的位置和落角。所以设三个如下所示的误差变量:
e31=v-vf3 (53);
ζ31=x-xf+(τ-τf3)cotθf3 (54);
其中,θf3=-90°为该阶段的期望终端落角,vf3为火箭期望着陆速度,xf为火箭初始位置离降落点的水平距离。
同时,为了约束位置和落角,设计滑模函数S32:
为方便讨论,设置两个中间控制量Pv和Pθ为:
结合滑模控制方法,中间控制量为:
综上可得控制量P和δ为:
表2垂直降落阶段初始参数
火箭的空载质量mf3=20000kg,结合表2、滑模函数S31的仿真参数和公式(44)可求得火箭在垂直降落段的初始质量至少为mf2=24728kg。由仿真结果可得,火箭着陆时的位置误差为3.6470e-09m,落角误差为1.8956e-04°,实际速度为0.98m/s,实际质量为20534kg,大于火箭的空载质量,满足实际需求。其他仿真结果如图6~图9所示。
Claims (3)
1.一种基于滑模控制的可重复使用运载火箭燃料计算方法,其特征在于所述方法包括如下步骤:
步骤一、建立以伪高度变量τ为自变量的数学模型:
其中,v为速度,t为时间变量,y为火箭在当前时刻的高度,θ为弹道倾角,m为导弹质量,g为地球重力加速度,P为发动机推力,δ为发动机偏转角,α为攻角,Isg为发动机比冲,Fe=[Fex,Fey]为离心惯性力,Fk=[Fkx,Fky]为哥氏惯性力,dw1和dw2为外部扰动,D为气动阻力,L为气动升力,右上带“'”表示关于τ的一阶导数;
步骤二、制定基于滑模控制理论的RLV垂直软着陆的速度控制策略
火箭回收可分为调姿段、动力减速段、气动减速段和垂直着陆段四个阶段,其中,动力减速段和垂直着陆段需要利用发动机推力控制火箭速度,由于动力减速段空气稀薄,垂直着陆段速度较小,即在控制速度时,气动力较小,可将其作为干扰项,又因为攻角和发动机偏转角较小,可以近似为零,则公式(3)简化为:
火箭垂直软着陆的关键在于速度控制,设如下状态误差变量e为:
e=v-vf (7);
其中,vf为常值期望速度,由公式(6)和公式(7)可得如下一阶非线性系统:
针对上述一阶非线性系统设计如下任意一阶滑模函数S:
滑模函数S关于伪高度变量τ的一阶导数为:
结合滑模控制理论,可得发动机控制推力P为:
步骤三、基于滑模控制理论的RLV垂直软着陆的速度控制策略计算火箭回收过程所需燃料
由公式(9)可得e的解析式为:
由公式(7)和公式(12)可得速度v的解析式为:
由公式(5)、公式(11)、公式(12)和公式(13)可得如下以伪高度τ为自变量的一阶齐次线性微分方程:
定义如下三个变量因子:
其中,D1为势能因子,D2为动能因子,D3为抗扰因子,τf为任意终端伪高度,τf=y0-yf,y0为火箭的初始高度,yf为终端高度;
若火箭的初始质量为m0且已知,则可得火箭的理论终端质量mf为:
如果已知火箭不含燃料的空载质量为mf,则可得火箭理论所需初始质量m0:
由公式(18)、公式(19)可得理论所需燃料为Δm=m0-mf。
2.根据权利要求1所述的基于滑模控制的可重复使用运载火箭燃料计算方法,其特征在于所述τ=y0-y,y0为火箭的初始高度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210314878.7A CN114676574B (zh) | 2022-03-28 | 2022-03-28 | 基于滑模控制的可重复使用运载火箭燃料计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210314878.7A CN114676574B (zh) | 2022-03-28 | 2022-03-28 | 基于滑模控制的可重复使用运载火箭燃料计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114676574A CN114676574A (zh) | 2022-06-28 |
CN114676574B true CN114676574B (zh) | 2022-11-11 |
Family
ID=82076951
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210314878.7A Active CN114676574B (zh) | 2022-03-28 | 2022-03-28 | 基于滑模控制的可重复使用运载火箭燃料计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114676574B (zh) |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102011076197A1 (de) * | 2011-05-20 | 2012-11-22 | Ford Global Technologies, Llc | Brennkraftmaschine mit Ölkreislauf und Verfahren zum Betreiben einer derartigen Brennkraftmaschine |
KR20190020852A (ko) * | 2015-02-13 | 2019-03-04 | 베이징 디디 인피니티 테크놀로지 앤드 디벨럽먼트 컴퍼니 리미티드 | 운송 능력 스케줄링을 위한 방법들 및 시스템들 |
CN106647255B (zh) * | 2016-09-29 | 2019-10-11 | 国电科学技术研究院有限公司 | 按设定误差性能指标函数收敛的火电机组协调控制方法 |
CN107831653B (zh) * | 2017-10-16 | 2020-07-28 | 南京航空航天大学 | 一种抑制参数摄动的高超声速飞行器指令跟踪控制方法 |
CN112550770B (zh) * | 2020-12-15 | 2021-07-13 | 北京航天自动控制研究所 | 一种基于凸优化的火箭软着陆轨迹规划方法 |
-
2022
- 2022-03-28 CN CN202210314878.7A patent/CN114676574B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN114676574A (zh) | 2022-06-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109426146B (zh) | 高超声速飞行器的高阶非奇异Terminal滑模控制方法 | |
CN106021628A (zh) | 一种运载火箭垂直返回弹道设计方法 | |
CN113778129B (zh) | 一种干扰补偿的高超声速变后掠翼飞行器跟踪控制方法 | |
CN107491080B (zh) | 一种基于非线性反馈和微分跟踪的飞行器姿态控制方法 | |
CN107121929B (zh) | 基于线性协方差模型预测控制的鲁棒再入制导方法 | |
CN111399531A (zh) | 一种高超声速飞行器滑翔段制导与姿态控制一体化设计方法 | |
Kovach et al. | A bank-to-turn autopilot design for an advanced air-to-air interceptor | |
CN112000127B (zh) | 一种基于反步法的飞行器横侧向联合控制方法 | |
CN104567545A (zh) | Rlv大气层内主动段的制导方法 | |
Bhardwaj et al. | Thrust command based integrated reference model with envelope protections for tilt-rotor vtol transition uav | |
CN106527128A (zh) | 兼顾瞬态响应与鲁棒稳定性的飞行控制律设计新方法 | |
CN115933733A (zh) | 一种固定翼无人机纵向高度速度解耦非线性控制方法 | |
CN107678442B (zh) | 一种基于双模型下的四旋翼自主着船控制方法 | |
CN114676574B (zh) | 基于滑模控制的可重复使用运载火箭燃料计算方法 | |
CN113885543B (zh) | 一种复杂海况干扰下的航行体自抗扰横滚控制器的设计方法 | |
CN110780676A (zh) | 一种受扰小型无人直升机的复合主动抗干扰轨迹跟踪控制方法 | |
Yamasaki et al. | Robust path-following for UAV using pure pursuit guidance | |
CN116301028B (zh) | 基于吸气式高超声速平台的多约束在线飞行轨迹规划中段导引方法 | |
CN114265420B (zh) | 适于高动态、控制慢响应的制导控制一体化设计方法 | |
CN116795126A (zh) | 一种输入饱和与输出受限的变形飞行器控制方法 | |
CN110209179A (zh) | 一种高超声速飞行器的快速高度跟踪算法 | |
Kim et al. | Guidance and control system design for impact angle control of guided bombs | |
Yao et al. | Trajectory tracking controller based on PID-NLADRC | |
Dong et al. | Research on guidance and control law design of decelerating transition and vertical landing for a STOVL UAV | |
CN114610057A (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 |