CN110597068A - 考虑攻角非对称约束的高超声速飞行器鲁棒控制方法 - Google Patents
考虑攻角非对称约束的高超声速飞行器鲁棒控制方法 Download PDFInfo
- Publication number
- CN110597068A CN110597068A CN201910981928.5A CN201910981928A CN110597068A CN 110597068 A CN110597068 A CN 110597068A CN 201910981928 A CN201910981928 A CN 201910981928A CN 110597068 A CN110597068 A CN 110597068A
- Authority
- CN
- China
- Prior art keywords
- given
- designer
- angle
- attack angle
- follows
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明涉及一种考虑攻角非对称约束的高超声速飞行器鲁棒控制方法,本方法考虑高度跟踪情形将攻角作为***状态,基于非对称障碍李雅普诺夫函数设计控制律,将攻角跟踪误差非对称约束信息引入控制律,保证误差被限制在预设的非对称区间内,同时对攻角虚拟控制进行限幅,结合二者实现攻角非对称约束,保证超燃冲压发动机正常工作。考虑***存在不确定性,本方法在控制器设计中,加入鲁棒项估计未知非线性函数上界,对不确定带来的影响进行补偿。针对速度子***设计PID控制器实现速度跟踪。
Description
技术领域
本发明涉及一种飞行器控制方法,特别是涉及一种考虑攻角非对称约束的高超声速飞行器鲁棒控制方法,属于飞行控制领域。
背景技术
高超声速飞行器具有飞行速度快、突防能力强的特点,在军事方面使得全球快速打击成为可能,在民用方面可用于洲际快速运输,因而得到世界各国的广泛关注。吸气式高超声速飞行器多以超燃冲压发动机提供动力,而超燃冲压发动机的正常进气对飞行器攻角要求苛刻,且攻角为负值时发动机更易出现故障,因此需保证其幅值不超出给定范围。已有研究针对高超声速飞行器控制多未考虑攻角约束,或仅考虑对称约束进行设计,未能保证飞行器的攻角非对称约束,可能导致超燃冲压发动机无法正常工作。
《Barrier Lyapunov function-based adaptive control for hypersonicflight vehicles》(Hao An,Hongwei Xia,Changhong Wang,Nonlinear Dynamics,2017,88(3):1833–1853)一文针对高超声速飞行器攻角约束问题,基于障碍李雅普诺夫函数设计控制器,限制了俯仰角跟踪误差,但其并未对攻角本身进行限制,且仅考虑了对称形式的障碍李雅普诺夫函数,在攻角为负值与正值情形下未分别引入不同的约束,缺乏精细调整能力。
发明内容
要解决的技术问题
为解决由于超燃冲压发动机进气需求导致的高超声速飞行器攻角非对称约束问题,本发明提出了一种考虑攻角非对称约束的高超声速飞行器鲁棒控制方法。
技术方案
一种考虑攻角非对称约束的高超声速飞行器鲁棒控制方法,其特征在于步骤如下:
步骤1:考虑高超声速飞行器纵向通道动力学模型:
其中,表示推力,表示升力,表示阻力,表示俯仰力矩;V表示速度,γ表示航迹倾角,h表示高度,α表示攻角,q表示俯仰角速度;δe表示升降舵偏角,Φ表示节流阀开度; 表示动压,ρ表示空气密度,表示平均气动弦长,S表示气动参考面积;m、Iyy和g代表质量、俯仰轴的转动惯量以及重力加速度, 均为气动参数;
步骤2:定义高度跟踪误差eh=h-hd,设计航迹角指令γd:
式中,hd为高度指令,由设计者给出,为高度指令的一阶微分,kh>0,ki>0由设计者给出;考虑巡航段航迹角变化小,航迹角指令的一阶微分取为零;
取x1=γ,x2=α,x3=q;可将式(3)-(5)写成如下严格反馈形式:
其中,fi=fi0+Δfi,i=1,2,3,fi0为已知标称值,Δfi为***不确定性且满足|Δfi|<Δfim;gi,i=1,2,3为已知非线性函数;
步骤3:第1步:定义航迹角跟踪误差为:
e1=x1-γd (8)
设计攻角虚拟控制量为:
其中,z1为待设计的补偿误差项, 为不确定性上界的估计值,k1>0、λ1>0、σ>0均由设计者给定;
为了使攻角虚拟控制满足给定的约束条件,令x2c通过如下的饱和函数获得x2l:
其中,与x2c分别为x2c上下界,由设计者给出;
定义b2=x2l-x2c,z1=e1-z10,其中z10由以下自适应律得到:
设计自适应律如下:
其中,η1>0及δ1>0由设计者给出;
设计一阶微分器如下:
其中,p20>0、p21>0均由设计者给出;
第2步:定义攻角跟踪误差为:
e2=x2-x2c (14)
构造非对称障碍李雅普诺夫函数:
其中为误差e2的上界,-ka<0为误差e2的下界,均由设计者给出;
定义设计俯仰角速率虚拟控制量为:
其中, 为不确定性上界的估计值,k2>0、λ2>0均由设计者给定;
设计自适应律如下:
其中,η2>0及δ2>0由设计者给出;
设计一阶微分器如下:
其中,p30>0、p31>0均由设计者给出;
步骤4:定义俯仰角速率跟踪误差为:
e3=x3-x3c (19)
设计升降舵偏角为:
其中, 为不确定性上界的估计值,k3>0、λ3>0均由设计者给定;
设计自适应律如下:
其中,η3>0及δ3>0由设计者给出;
步骤5:定义速度跟踪误差为:
其中,Vd为速度参考指令;
设计节流阀开度Φ为:
其中,kpV>0、kiV>0和kdV>0由设计者给定;
步骤6:根据得到的升降舵偏角δe和速度子***的节流阀开度Φ,返回到飞行器动力学模型(1)-(5),对高度和速度进行跟踪控制。
有益效果
本发明提出的一种考虑攻角非对称约束的高超声速飞行器鲁棒控制方法,本方法考虑高度跟踪情形将攻角作为***状态,基于非对称障碍李雅普诺夫函数设计控制律,将攻角跟踪误差非对称约束信息引入控制律,保证误差被限制在预设的非对称区间内,同时对攻角虚拟控制进行限幅,结合二者实现攻角非对称约束,保证超燃冲压发动机正常工作。考虑***存在不确定性,本方法在控制器设计中,加入鲁棒项估计未知非线性函数上界,对不确定带来的影响进行补偿。针对速度子***设计PID控制器实现速度跟踪。与现有技术相比有益效果为:
(1)本发明通过构造非对称障碍李雅普诺夫函数,将攻角跟踪误差的不同上下界约束直接引入控制律设计中,确保了攻角跟踪误差被限制在给定的非对称区间内,保证了跟踪精度。
(2)本发明在限制攻角跟踪误差的同时对攻角虚拟控制采取非对称限幅设计,结合二者实现飞行过程中攻角能够被约束在人为给定的非对称范围内,可保证超燃冲压发动机在不同飞行情形下的进气需求。
附图说明
图1是本发明考虑攻角非对称约束的高超声速飞行器鲁棒控制的流程图。
具体实施方式
现结合实施例、附图对本发明作进一步描述:
参照图1,本发明考虑攻角非对称约束的高超声速飞行器鲁棒控制应用于一类高超声速飞行器动力学模型,通过以下步骤实现:
(a)考虑高超声速飞行器纵向通道动力学模型:
其中,表示推力,表示升力,表示阻力,表示俯仰力矩;V表示速度,γ表示航迹倾角,h表示高度,α表示攻角,q表示俯仰角速度。δe表示升降舵偏角,Φ表示节流阀开度。 表示动压,ρ=6.7429×10-5表示空气密度,表示平均气动弦长,S=17表示气动参考面积;m=300、Iyy=5×105、g=32代表质量、俯仰轴的转动惯量以及重力加速度,均为气动参数,具体定义如下:
(b)定义高度跟踪误差eh=h-hd,设计航迹角指令γd:
式中,hd为高度指令,由设计者给出,为高度指令的一阶微分,kh=0.5,ki=0.05。考虑巡航段航迹角变化小,航迹角指令的一阶微分取为零。
取x1=γ,x2=α,x3=q;可将式(3)-(5)写成如下严格反馈形式:
其中,fi=fi0+Δfi,i=1,2,3,fi0=0.9fi为已知标称值,Δfi为***不确定性且满足|Δfi|<Δfim;gi,i=1,2,3为已知非线性函数。
(c)第1步:定义航迹角跟踪误差为:
e1=x1-γd (8)
设计攻角虚拟控制量为:
其中,z1为待设计的补偿误差项, 为不确定性上界的估计值,k1=0.8,λ1=1.5,σ=10均由设计者给定。
为了使攻角虚拟控制满足给定的约束条件,令x2c通过如下的饱和函数获得x2l:
其中,与x2c=-0.02分别为x2c上下界。
定义b2=x2l-x2c,z1=e1-z10,其中z10由以下自适应律得到:
设计自适应律如下:
其中,η1=0.1,δ1=0.001。
设计一阶微分器如下:
其中,p20>0.9,p21>1.2。
第2步:定义攻角跟踪误差为:
e2=x2-x2c (14)
构造非对称障碍李雅普诺夫函数:
其中为误差e2的上界,-ka=-0.07为误差e2的下界。
定义设计俯仰角速率虚拟控制量为:
其中, 为不确定性上界的估计值,k2=1,λ2=30。
设计自适应律如下:
其中,η2=0.1,δ2=0.001。
设计一阶微分器如下:
其中,p30=0.05,p31=2。
第3步:定义俯仰角速率跟踪误差为:
e3=x3-x3c (19)
设计升降舵偏角为:
其中, 为不确定性上界的估计值,k3=35,λ3=10。
设计自适应律如下:
其中,η3=0.1及δ3=0.001由设计者给出。
(d)定义速度跟踪误差为:
其中,Vd为速度参考指令。
设计节流阀开度Φ为:
其中,kpV=5,kiV=1,kdV=0.01。
(e)根据得到的升降舵偏角δe和速度子***的节流阀开度Φ,返回到飞行器动力学模型(1)-(5),对高度和速度进行跟踪控制。
Claims (1)
1.一种考虑攻角非对称约束的高超声速飞行器鲁棒控制方法,其特征在于步骤如下:
步骤1:考虑高超声速飞行器纵向通道动力学模型:
其中,表示推力,表示升力,表示阻力,表示俯仰力矩;V表示速度,γ表示航迹倾角,h表示高度,α表示攻角,q表示俯仰角速度;δe表示升降舵偏角,Φ表示节流阀开度; 表示动压,ρ表示空气密度,表示平均气动弦长,S表示气动参考面积;m、Iyy和g代表质量、俯仰轴的转动惯量以及重力加速度, 均为气动参数;
步骤2:定义高度跟踪误差eh=h-hd,设计航迹角指令γd:
式中,hd为高度指令,由设计者给出,为高度指令的一阶微分,kh>0,ki>0由设计者给出;考虑巡航段航迹角变化小,航迹角指令的一阶微分取为零;
取x1=γ,x2=α,x3=q;可将式(3)-(5)写成如下严格反馈形式:
其中,fi=fi0+Δfi,i=1,2,3,fi0为已知标称值,Δfi为***不确定性且满足|Δfi|<Δfim;gi,i=1,2,3为已知非线性函数;
步骤3:第1步:定义航迹角跟踪误差为:
e1=x1-γd (8)
设计攻角虚拟控制量为:
其中,z1为待设计的补偿误差项, 为不确定性上界的估计值,k1>0、λ1>0、σ>0均由设计者给定;
为了使攻角虚拟控制满足给定的约束条件,令x2c通过如下的饱和函数获得x2l:
其中,与x 2c分别为x2c上下界,由设计者给出;
定义b2=x2l-x2c,z1=e1-z10,其中z10由以下自适应律得到:
设计自适应律如下:
其中,η1>0及δ1>0由设计者给出;
设计一阶微分器如下:
其中,p20>0、p21>0均由设计者给出;
第2步:定义攻角跟踪误差为:
e2=x2-x2c (14)
构造非对称障碍李雅普诺夫函数:
其中kb>0为误差e2的上界,-ka<0为误差e2的下界,均由设计者给出;
定义设计俯仰角速率虚拟控制量为:
其中, 为不确定性上界的估计值,k2>0、λ2>0均由设计者给定;
设计自适应律如下:
其中,η2>0及δ2>0由设计者给出;
设计一阶微分器如下:
其中,p30>0、p31>0均由设计者给出;
步骤4:定义俯仰角速率跟踪误差为:
e3=x3-x3c (19)
设计升降舵偏角为:
其中, 为不确定性上界的估计值,k3>0、λ3>0均由设计者给定;
设计自适应律如下:
其中,η3>0及δ3>0由设计者给出;
步骤5:定义速度跟踪误差为:
其中,Vd为速度参考指令;
设计节流阀开度Φ为:
其中,kpV>0、kiV>0和kdV>0由设计者给定;
步骤6:根据得到的升降舵偏角δe和速度子***的节流阀开度Φ,返回到飞行器动力学模型(1)-(5),对高度和速度进行跟踪控制。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910981928.5A CN110597068A (zh) | 2019-10-16 | 2019-10-16 | 考虑攻角非对称约束的高超声速飞行器鲁棒控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910981928.5A CN110597068A (zh) | 2019-10-16 | 2019-10-16 | 考虑攻角非对称约束的高超声速飞行器鲁棒控制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110597068A true CN110597068A (zh) | 2019-12-20 |
Family
ID=68867502
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910981928.5A Pending CN110597068A (zh) | 2019-10-16 | 2019-10-16 | 考虑攻角非对称约束的高超声速飞行器鲁棒控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110597068A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111158398A (zh) * | 2020-01-15 | 2020-05-15 | 哈尔滨工业大学 | 考虑攻角约束的高超音速飞行器的自适应控制方法 |
CN111273681A (zh) * | 2020-04-09 | 2020-06-12 | 中北大学 | 一种考虑攻角受限的高超声速飞行器高安全抗干扰控制方法 |
CN112068444A (zh) * | 2020-09-22 | 2020-12-11 | 中国人民解放军海军航空大学 | 一种采用非线性自适应滑模的飞行器攻角控制方法 |
CN113110543A (zh) * | 2021-04-19 | 2021-07-13 | 西北工业大学 | 一种非线性非最小相位飞行器的鲁棒飞行控制方法 |
CN113448339A (zh) * | 2020-03-25 | 2021-09-28 | 中国人民解放军海军工程大学 | 一种基于虚拟反演的飞行器攻角跟踪控制方法 |
CN114200827A (zh) * | 2021-11-09 | 2022-03-18 | 西北工业大学 | 一种超音速大机动靶标的多约束双通道控制方法 |
CN114779636A (zh) * | 2022-04-17 | 2022-07-22 | 西北工业大学 | 考虑气动伺服弹性的飞行器鲁棒自适应控制方法 |
CN116594414A (zh) * | 2023-03-28 | 2023-08-15 | 西北工业大学 | 一种高超声速飞行器纵向控制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104155990A (zh) * | 2014-08-15 | 2014-11-19 | 哈尔滨工业大学 | 考虑攻角约束的高超声速飞行器俯仰通道姿态控制方法 |
CN107450324A (zh) * | 2017-09-05 | 2017-12-08 | 西北工业大学 | 考虑攻角约束的高超声速飞行器自适应容错控制方法 |
CN107831671A (zh) * | 2017-12-06 | 2018-03-23 | 浙江工业大学 | 一种基于非对称时变障碍李雅普诺夫函数的四旋翼飞行器输出受限反步控制方法 |
CN109375639A (zh) * | 2018-11-27 | 2019-02-22 | 浙江工业大学 | 一种基于非对称改进型障碍李雅普诺夫函数的刚性飞行器姿态约束跟踪控制方法 |
-
2019
- 2019-10-16 CN CN201910981928.5A patent/CN110597068A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104155990A (zh) * | 2014-08-15 | 2014-11-19 | 哈尔滨工业大学 | 考虑攻角约束的高超声速飞行器俯仰通道姿态控制方法 |
CN107450324A (zh) * | 2017-09-05 | 2017-12-08 | 西北工业大学 | 考虑攻角约束的高超声速飞行器自适应容错控制方法 |
CN107831671A (zh) * | 2017-12-06 | 2018-03-23 | 浙江工业大学 | 一种基于非对称时变障碍李雅普诺夫函数的四旋翼飞行器输出受限反步控制方法 |
CN109375639A (zh) * | 2018-11-27 | 2019-02-22 | 浙江工业大学 | 一种基于非对称改进型障碍李雅普诺夫函数的刚性飞行器姿态约束跟踪控制方法 |
Non-Patent Citations (2)
Title |
---|
BIN XU 等: "Barrier Lyapunov Function Based Learning Control of Hypersonic Flight Vehicle With AOA Constraint and Actuator Faults", 《IEEE TRANSACTIONS ON CYBERNETICS》 * |
YUYAN GUO 等: "Adaptive Hypersonic Flight Control under Asymmetric AOA Constraint", 《第三十八届中国控制会议论文集(6)》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111158398A (zh) * | 2020-01-15 | 2020-05-15 | 哈尔滨工业大学 | 考虑攻角约束的高超音速飞行器的自适应控制方法 |
CN113448339A (zh) * | 2020-03-25 | 2021-09-28 | 中国人民解放军海军工程大学 | 一种基于虚拟反演的飞行器攻角跟踪控制方法 |
CN113448339B (zh) * | 2020-03-25 | 2022-07-01 | 中国人民解放军海军工程大学 | 一种基于虚拟反演的飞行器攻角跟踪控制方法 |
CN111273681A (zh) * | 2020-04-09 | 2020-06-12 | 中北大学 | 一种考虑攻角受限的高超声速飞行器高安全抗干扰控制方法 |
CN112068444A (zh) * | 2020-09-22 | 2020-12-11 | 中国人民解放军海军航空大学 | 一种采用非线性自适应滑模的飞行器攻角控制方法 |
CN112068444B (zh) * | 2020-09-22 | 2022-02-15 | 中国人民解放军海军航空大学 | 一种采用非线性自适应滑模的飞行器攻角控制方法 |
CN113110543A (zh) * | 2021-04-19 | 2021-07-13 | 西北工业大学 | 一种非线性非最小相位飞行器的鲁棒飞行控制方法 |
CN113110543B (zh) * | 2021-04-19 | 2022-08-09 | 西北工业大学 | 一种非线性非最小相位飞行器的鲁棒飞行控制方法 |
CN114200827A (zh) * | 2021-11-09 | 2022-03-18 | 西北工业大学 | 一种超音速大机动靶标的多约束双通道控制方法 |
CN114779636A (zh) * | 2022-04-17 | 2022-07-22 | 西北工业大学 | 考虑气动伺服弹性的飞行器鲁棒自适应控制方法 |
CN116594414A (zh) * | 2023-03-28 | 2023-08-15 | 西北工业大学 | 一种高超声速飞行器纵向控制方法 |
CN116594414B (zh) * | 2023-03-28 | 2024-03-15 | 西北工业大学 | 一种高超声速飞行器纵向控制方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110597068A (zh) | 考虑攻角非对称约束的高超声速飞行器鲁棒控制方法 | |
CN107450324B (zh) | 考虑攻角约束的高超声速飞行器自适应容错控制方法 | |
CN110568765A (zh) | 面向攻角跟踪的高超声速飞行器非对称输出受限控制方法 | |
CN108828957B (zh) | 基于切换机制的飞行器全局有限时间神经网络控制方法 | |
CN111665857B (zh) | 基于复合智能学习的变体飞行器控制方法 | |
CN111679583B (zh) | 基于气动参数估计的变体飞行器自适应控制方法 | |
CN107390531B (zh) | 参数学习有限时间收敛的高超声速飞行器控制方法 | |
CN107942651B (zh) | 一种近空间飞行器控制*** | |
CN110109361B (zh) | 一种高超声速飞行器快速平滑自适应滑模控制方法 | |
CN114564045A (zh) | 一种考虑降雨和阵风条件的无人机飞行控制律设计方法 | |
CN107065554B (zh) | 近空间可变翼飞行器自适应小翼切换控制***及工作方法 | |
CN102880052A (zh) | 基于时标功能分解的高超声速飞行器执行器饱和控制方法 | |
CN113778129B (zh) | 一种干扰补偿的高超声速变后掠翼飞行器跟踪控制方法 | |
CN108762098B (zh) | 基于复合学习的非最小相位飞行器神经网络控制方法 | |
CN107632518B (zh) | 基于高增益观测器的高超声速飞行器神经网络复合学习控制方法 | |
CN103592847A (zh) | 一种基于高增益观测器的高超声速飞行器非线性控制方法 | |
CN109062234B (zh) | 一种非最小相位飞行器复合学习滑模控制方法 | |
CN110908288A (zh) | 一种基于障碍Lyapunov函数的无人机地速约束抗扰控制方法 | |
CN114859712B (zh) | 面向油门约束的飞行器制导控制一体化方法 | |
CN113126495B (zh) | 基于地效干扰补偿的低空飞行鲁棒智能控制方法 | |
CN115079574A (zh) | 一种挠性高超声速飞行器分布式故障补偿方法 | |
CN111596686B (zh) | 高超声速飞行器纵向***的预设性能控制方法 | |
CN112149234B (zh) | 一种基于俯仰角速率输入的飞行器质点运动模型设计方法 | |
CN110989397B (zh) | 一种航空器失事搜寻仿真方法及*** | |
CN114690793B (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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20191220 |
|
WD01 | Invention patent application deemed withdrawn after publication |