CN115328185A - 一种飞行器非线性非定常气动载荷修正系统 - Google Patents
一种飞行器非线性非定常气动载荷修正系统 Download PDFInfo
- Publication number
- CN115328185A CN115328185A CN202211053328.0A CN202211053328A CN115328185A CN 115328185 A CN115328185 A CN 115328185A CN 202211053328 A CN202211053328 A CN 202211053328A CN 115328185 A CN115328185 A CN 115328185A
- Authority
- CN
- China
- Prior art keywords
- angular velocity
- aircraft
- aerodynamic load
- attitude angle
- unsteady aerodynamic
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 16
- 238000012545 processing Methods 0.000 claims abstract description 28
- 238000000034 method Methods 0.000 claims abstract description 22
- 238000001914 filtration Methods 0.000 claims description 40
- 230000003068 static effect Effects 0.000 claims description 19
- 230000003044 adaptive effect Effects 0.000 claims description 17
- 230000008569 process Effects 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000005096 rolling process Methods 0.000 claims description 5
- 230000008439 repair process Effects 0.000 claims description 4
- 239000002131 composite material Substances 0.000 claims description 3
- 238000013016 damping Methods 0.000 claims description 3
- 230000006870 function Effects 0.000 claims description 3
- 230000004048 modification Effects 0.000 claims 1
- 238000012986 modification Methods 0.000 claims 1
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 abstract description 3
- 238000009826 distribution Methods 0.000 description 10
- 230000006872 improvement Effects 0.000 description 9
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000013461 design Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000012905 input function Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000003278 mimic effect Effects 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 235000008113 selfheal Nutrition 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007921 spray Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/08—Control of attitude, i.e. control of roll, pitch, or yaw
- G05D1/0808—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
- G05D1/0816—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability
- G05D1/0825—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability using mathematical models
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
- G05D1/106—Change initiated in response to external conditions, e.g. avoidance of elevated terrain or of no-fly zones
-
- 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)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明涉及一种飞行器非线性非定常气动载荷修正系统,属于飞行控制技术领域,解决了现有技术中没有针对飞行器非线性非定常气动载荷的补偿系统,存在无法自动修复飞行器的非线性非定常气动载荷、飞行稳定性和控制精度差等问题。通过输出飞行器在非线性非定常气动载荷下的实际姿态角和实际角速度,利用外反馈回路,对期望姿态角和实际姿态角进行处理,获得参考角速度;利用内反馈回路,对实际角速度和参考角速度进行处理,获得非线性非定常气动载荷的估计值,并输出补偿控制指令,对非线性非定常气动载荷进行修正。实现了针对飞行器对非线性非定常气动载荷的修正,提升飞行器的飞行稳定性和控制精度。
Description
技术领域
本发明涉及飞行控制技术领域,尤其涉及一种飞行器非线性非定常气动载荷修正系统。
背景技术
飞行器结构损伤会引起气动参数变化,进而影响系统的静稳定性和控制精度。飞行器在大迎角区域做机动时,气动力具有明显的非线性和非定常性,对正常式布局飞行器,气流以大迎角流过飞行器头部和具有大后掠角的边条翼时会产生旋涡绕流结构,机翼上的气流则从表面发生分离,旋涡区和分离区气流的相互作用,产生了很大的非线性非定常空气动力载荷,即摄动。进而影响飞行器的静稳定性和控制精度。
现有技术主要针对飞行器无非线性非定常气动载荷条件下的标称控制技术,或一般动态系统、流程控制系统的摄动补偿技术进行研究。例如,现存前沿飞行操纵系统提出了解决经典控制论问题的方案,针对大包线、全空域、全速域、急转弯等严格需求工况,科学家们设计了自适应动态逆系统,加入被观测器估计得到的状态和输入参数以抵消实际系统的不确定状态和输入,还加入了根据期望指令实时改变比例控制参数的自适应律,使算法结构简单且适应各种高机动模式;针对偏航舵、升降舵、副翼等机械舵面输入饱和问题,工程师提出了额外调控系统,用B样条逼近突变饱和的式子,再用方向控制推进装置即喷工质舵调节器建立硬件结构舵操纵回路里的虚飞控算子,使机械舵在舵偏角很小的情况下依然能实现高机动姿控;针对干扰条件下的鲁棒稳定问题,科学家设计了估计算法有源自动修复的想法,制造额外的估计回路和算法,对大幅值低频连续特征或小幅值高频离散特征的多干扰采取平行估计策略,继而让多个干扰观测值导进综控机重构操控软件,抵消干扰危害,实时自愈合飞行器。然而,上述系统都没考虑参数摄动情况,更没有针对飞行器非线性非定常气动载荷设计专门的补偿系统。现有摄动补偿控制算法则更多聚焦于通用动态系统。例如,采用博弈论思想解决互联系统摄动补偿问题;采用分段李雅普诺夫系统证明带摄动的分段线性系统的全局稳定性;用摄动与输入函数的位置关系将摄动划分为乘性摄动和加性摄动,并用经典鲁棒系统直接被动修复加性摄动;针对多输入多输出系统中的单通道摄动,用另一无摄动通道错位补偿单摄动对整体系统的影响,这些都是理论上常见的补偿系统。而现有针对工业工程现场偏应用的补偿系统也很多,例如,针对磨矿机执行机构摄动,为保障矿物颗粒粒径分布吻合期望随机分布函数,须首先建立非高斯随机分布系统和摄动模型,进而设计随机容错控制器修复执行器故障。
综上,目前的现有技术中没有针对飞行器非线性非定常气动载荷的补偿系统,存在无法自动修复飞行器的非线性非定常气动载荷、飞行稳定性和控制精度差等缺陷。
发明内容
鉴于上述的分析,本发明实施例旨在提供一种飞行器非线性非定常气动载荷修正系统,用以解决现有技术中没有针对飞行器非线性非定常气动载荷的补偿系统,存在无法自动修复飞行器的非线性非定常气动载荷、飞行稳定性和控制精度差的问题。
本发明的目的主要是通过以下技术方案实现的:
本发明实施例提供了一种飞行器非线性非定常气动载荷修正系统,包括:
飞行器模型和内、外双反馈回路;其中,
所述飞行器模型,用于在非线性非定常气动载荷下,输出实际姿态角和实际角速度;
所述外反馈回路,用于对期望姿态角和所述实际姿态角进行处理,获得参考角速度;
所述内反馈回路,用于对所述实际角速度和参考角速度进行处理,获得非线性非定常气动载荷的估计值,并输出补偿控制指令,对非线性非定常气动载荷进行修正。
基于上述方法的进一步改进,所述外反馈回路包括外回路指令滤波器和外回路反步控制器;
所述外回路指令滤波器,用于对期望姿态角进行滤波;
所述外回路反步控制器,用于对所述实际姿态角和滤波后的期望姿态角进行处理,获得参考角速度。
基于上述方法的进一步改进,所述内反馈回路包括内回路指令滤波器、自适应控制器和内回路反步控制器;
所述内回路指令滤波器,用于对所述参考角速度进行滤波;
所述自适应控制器,用于对所述实际角速度和滤波后的参考角速度进行处理,获得非线性非定常气动载荷的估计值;
所述内回路反步控制器,用于基于所述实际角速度、滤波后的参考角速度和所述非线性非定常气动载荷的估计值,将自身重构为自适应摄动补偿控制器,并输出补偿控制指令,对非线性非定常气动载荷进行修正。
基于上述方法的进一步改进,所述飞行器模型,描述为:
其中,输入变量u=[δa,δe,δr,δy,δz]T,δa表示副翼控制输入,δe表示升降舵控制输入,δr表示方向舵控制输入,δy表示横侧向推力矢量输入,δz表示纵向推理矢量输入;外回路状态变量x1=[α β φ]T,α表示迎角,β表示偏航角,φ表示滚转角;内回路状态变量x2=[p q r]T,p表示俯仰角速度,q表示偏航角速度,r表示滚转角速度;中的各元素均为气动参数,δClβ是横滚静稳定性导数,δCm是纵向静稳定性导数,是纵向过失速导数,δCnβ是航向静稳定性导数;
f2的表达式如下所示:
其中,fp(·)、fq(·)和fr(·)分别表示俯仰通道、偏航通道和滚转通道的非线性状态函数,记x0,y0,z0分别为飞机质心在地轴系三个坐标轴的投影,则I11,I22,I33分别表示机身绕x0轴的转动惯量、机身绕y0轴的转动惯量和机身绕z0轴的转动惯量,I13表示机身在x0-z0轴上的惯性积,分别为俯仰气动力矩、偏航气动力矩和滚转气动力矩;
角速率通道气动参数摄动模型,表示为:
θ2=θ2,0+Δθ2
上式中各项均为参数向量,各参数向量中的各元素表达式如下所示:
其中,带下标r的气动系数为飞机的静态气动参数,方程组等号左侧为整体气动参数,κi,i=1,2,3,4,5是气动参数摄动系数,当摄动为0时,气动参数摄动系数值为0;当摄动不为0时,满足1≤|κi|≤10,整体气动参数为静态气动参数与摄动系数的和。
基于上述方法的进一步改进,所述外回路指令滤波器和内回路指令滤波器,包括:
所述外回路指令滤波器为第1级滤波器,所述内回路指令滤波器为第2级滤波器,即i=1,2;第i级滤波器的形式如下:
滤波后的期望状态变量,满足:
其中,xic为滤波后的期望状态变量。
基于上述方法的进一步改进,所述外回路指令滤波器,用于对期望姿态角进行滤波,包括:
所述外回路指令滤波器,用于对滤波前的期望姿态角x1d进行滤波处理,获得滤波后的飞行器的期望姿态角x1c。
基于上述方法的进一步改进,所述外回路反步控制器,用于对所述实际姿态角和滤波后的期望姿态角进行处理,获得参考角速度,包括:
将经外回路指令滤波后的期望姿态角x1c,与实际姿态角x1做差,得到姿态角误差e1=x1-x1c;
对姿态角误差和滤波后的期望姿态角的导数,按照如下方式进行反步动态逆控制处理,获得参考角速度x2d,即第2级滤波前的期望状态变量:
其中,k1>0为外回路比例控制参数。
基于上述方法的进一步改进,所述内回路指令滤波器,用于对所述参考角速度进行滤波,包括:
所述内回路指令滤波器,用于对参考角速度x2d,即第2级滤波前的期望状态变量,进行第2级指令滤波处理,即内回路指令滤波处理,处理后,获得滤波后的参考角速度x2c,即内回路状态变量滤波后的期望值。
基于上述方法的进一步改进,所述自适应控制器,用于对所述实际角速度和滤波后的参考角速度进行处理,获得非线性非定常气动载荷的估计值,包括:
将实际角速度与滤波后的参考角速度做差,得到角速度误差e2=x2-x2c;
通过以下方式,获得非线性非定常气动载荷估计值的导数:
其中,Γ2是自适应学习率,υ2为跟踪补偿误差,且满足υ2=e2;
基于上述方法的进一步改进,所述内回路反步控制器,用于基于所述实际角速度、滤波后的参考角速度和所述非线性非定常气动载荷的估计值,将自身重构为自适应摄动补偿控制器,并输出补偿控制指令,对非线性非定常气动载荷进行修正,包括:
所述内回路反步控制器,对非线性非定常气动载荷的估计值、角速度误差和滤波后的参考角速度的导数进行处理,将自身重构为自适应摄动补偿控制器:
其中,k2为比例控制参数;u为补偿控制指令,即控制输出,表示控制器期望控制输入产生的舵面偏转角;
所述自适应摄动补偿控制器,发出补偿控制指令,对飞行器进行自适应摄动补偿控制,完成对所述非线性非定常气动载荷的修复。
与现有技术相比,本发明至少可实现如下有益效果之一:
1、本发明提供了一种飞行器非线性非定常气动载荷修正系统,使得飞行器在飞行中即使出现非线性非定常气动载荷也能保持稳定。
2、本发明在公开一个自适应摄动补偿控制器的同时,也提供了一个自适应控制器,该自适应控制器估计得到的准确摄动信息不仅有助于修复摄动,也有助于驾驶员基于复杂多元任务做综合决策。
3、本发明是在标称动态逆控制的基础上设计的自适应补偿控制算法,具有动态逆控制的优点,即能解决非线性问题、结构简单易实现、可靠性高等,为剩余通道设计其它复杂控制器预留了设计空间。
本发明中,上述各技术方案之间还可以相互组合,以实现更多的优选组合方案。本发明的其他特征和优点将在随后的说明书中阐述,并且,部分优点可从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过说明书以及附图中所特别指出的内容中来实现和获得。
附图说明
附图仅用于示出具体实施例的目的,而并不认为是对本发明的限制,在整个附图中,相同的参考符号表示相同的部件。
图1为为本发明实施例中控制系统的架构图。
具体实施方式
下面结合附图来具体描述本发明的优选实施例,其中,附图构成本申请一部分,并与本发明的实施例一起用于阐释本发明的原理,并非用于限定本发明的范围。
本实施例公开了一种飞行器非线性非定常气动载荷修正系统,如图1所示,包括飞行器模型和内、外双反馈回路;其中,
所述飞行器模型,用于在非线性非定常气动载荷下,输出飞行器的实际姿态角和实际角速度;
所述外反馈回路,用于基于期望姿态角和所述实际姿态角进行处理,获得参考角速度;
所述内反馈回路,用于对所述实际角速度和参考角速度进行处理,获得非线性非定常气动载荷的估计值,并输出相应的补偿控制指令至飞行器,对非线性非定常气动载荷进行修正。
具体来说,所述飞行器模型,表示为:
其中,输入变量u=[δa,δe,δr,δy,δz]T,δa表示副翼控制输入,δe表示升降舵控制输入,δr表示方向舵控制输入,δy表示横侧向推力矢量输入,δz表示纵向推理矢量输入;外回路状态变量x1=[α β φ]T,α表示迎角,β表示偏航角,φ表示滚转角;内回路状态变量x2=[p q r]T,p表示俯仰角速度,q表示偏航角速度,r表示滚转角速度;中的各元素均为气动参数,δClβ是横滚静稳定性导数,δCm是纵向静稳定性导数,是纵向过失速导数,δCnβ是航向静稳定性导数。
f2的表达式如下所示:
其中,fp(·)、fq(·)和fr(·)分别表示俯仰通道、偏航通道和滚转通道的非线性状态函数,设x0,y0,z0分别为飞行器质心在地轴系三个坐标轴的投影,则I11,I22,I33分别表示机身绕x0轴的转动惯量、机身绕y0轴的转动惯量和机身绕z0轴的转动惯量,I13表示机身在x0-z0轴上的惯性积,分别为俯仰气动力矩、偏航气动力矩和滚转气动力矩。
其中,G1和G2为控制效率分配矩阵,表示副翼在俯仰通道的控制分配系数,表示方向舵在俯仰通道的控制分配系数,表示推力矢量横侧向舵在俯仰通道的控制分配系数,表示升降舵在偏航通道的控制分配系数,表示推力矢量纵向舵在偏航通道的控制分配系数,表示副翼在滚转通道的控制分配系数,表示方向舵在滚转通道的控制分配系数,表示推力矢量横侧向舵在滚转通道的控制分配系数;为气动和结构复合参数矩阵;其中,c3,c4,c7和c9是气动常数,qd是动压,b是翼展,S是参考机翼面积,c是平均气动弦长,中的物理量均通过风洞实验测定。
角速率通道气动参数摄动模型,表示为:
上式中各项均为参数向量,各参数向量中的各元素表达式如下所示:
其中,带下标r的气动系数为飞行器的静态气动参数,方程组等号左侧为整体气动参数,κi,i=1,2,3,4,5是气动参数摄动系数,当摄动为0时,气动参数摄动系数值为0,则静态气动参数与整体气动参数相等;当摄动不为0时,满足1≤|κi|≤10,整体气动参数为静态气动参数与摄动系数的和。
采用这种形式的气动参数摄动,主要出于两方面的考虑。其一是为了模仿大迎角下的气动参数的迟滞效应。飞行器气动参数的迟滞效应要是受到迎角的角速度的影响,的值越大,流场的迟滞效应也就越明显。而飞行器的迟滞效应主要在大迎角状态下发生,迎角α的大小对于迟滞环的大小也有显著的影响。所以在考虑迟滞环的大小时,用来对气动参数进行修正。
另外一个考虑就是检验控制系统的鲁棒性。由于飞行器在大迎角飞行时,特别是在进行过失速机动时,流场非常复杂,要非常准确的获得飞行器的气动参数几乎不可能,所以气动参数误差的存在是难免的,而且有时候还非常大。而飞行器的控制律是基于试飞(或风洞实验)得到的气动参数而设计的,所以要使控制系统可用的条件就是飞行器控制系统具有较好的鲁棒性。所以在控制律设计时使用静态的气动参数,而在计算飞行器机动过程中的气动力和力矩时,利用对其进行修正。
所述外反馈回路包括外回路指令滤波器和外回路反步控制器;
所述外回路指令滤波器,用于对期望姿态角进行滤波;
所述外回路反步控制器,用于对所述实际姿态角和滤波后的期望姿态角进行处理,获得参考角速度。
所述内反馈回路包括内回路指令滤波器、自适应控制器和内回路反步控制器;
所述内回路指令滤波器,用于对所述参考角速度进行滤波;
所述自适应控制器,用于对所述实际角速度和滤波后的参考角速度进行处理,获得非线性非定常气动载荷的估计值;
所述内回路反步控制器,用于基于所述实际角速度、滤波后的参考角速度和所述非线性非定常气动载荷的估计值,将自身重构为自适应摄动补偿控制器,并输出相应的补偿控制指令,对非线性非定常气动载荷进行修正。
所述飞行器模型输出外回路反馈,得到飞行器的实际姿态角,即外回路状态变量x1=[α β φ]T。
所述外回路指令滤波器,以滤波前的期望姿态角为输入,输出滤波后的期望姿态角。
如图1所示,所述外回路指令滤波器为第1级滤波器,所述内回路指令滤波器为第2级滤波器;第i级滤波器的形式如下,i=1,2:
滤波后的期望状态变量,满足:
其中,xic为第i级滤波器滤波后的期望状态变量。
此处的二阶指令滤波器,将经典的自适应反步控制中对中间信号的求导替换为求解二阶指令滤波器的状态变量,避免了对复杂运动学方程进行求导和微分爆炸,减小了计算量。
具体来说,所述外回路指令滤波器,用于对滤波前的期望姿态角x1d进行滤波处理,获得滤波后的飞行器的期望姿态角x1c。
x1c是经外回路指令滤波后的期望姿态角,x1是实际姿态角,则姿态角误差e1=x1-x1c。
具体来说,所述外回路反步控制器,用于对姿态角误差和滤波后的期望姿态角的导数,按照如下方式进行反步动态逆控制处理,获得参考角速度x2d,即第2级滤波前的期望状态变量:
其中,k1>0为外回路比例控制参数。
所述飞行器模型输出内回路反馈,得到飞行器的实际角速度,即内回路状态变量x2=[p q r]T。
具体来说,所述内回路指令滤波器,用于对参考角速度x2d,即第2级滤波前的期望状态变量,进行第2级指令滤波处理,即内回路指令滤波处理,处理后,获得滤波后的参考角速度x2c,即内回路状态变量滤波后的期望值。
将实际角速度与滤波后的参考角速度做差,得到角速度误差e2=x2-x2c。
需要说明的是,对非线性非定常气动载荷进行自适应估计,需要利用跟踪补偿误差,而跟踪补偿误差即角速度误差:
υ2=e2
其中,υ2为跟踪补偿误差。
具体来说,所述自适应控制器,用于通过以下方式,获得非线性非定常气动载荷估计值的导数:
其中,Γ2是自适应学习率。
具体来说,所述内回路反步控制器,用于对非线性非定常气动载荷的估计值、角速度误差和滤波后的参考角速度的导数进行处理,将自身重构为自适应摄动补偿控制器:
其中,k2为比例控制参数;u为补偿控制指令,即控制输出,表示控制器期望控制输入产生的舵面偏转角。
所述自适应摄动补偿控制器,用于发出相应的补偿控制指令,对飞行器进行自适应摄动补偿控制,完成对所述非线性非定常气动载荷的修复。
本领域技术人员可以理解,实现上述实施例方法的全部或部分流程,可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于计算机可读存储介质中。其中,所述计算机可读存储介质为磁盘、光盘、只读存储记忆体或随机存储记忆体等。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
Claims (10)
1.一种飞行器非线性非定常气动载荷修正系统,其特征在于,包括飞行器模型和内、外反馈回路;其中,
所述飞行器模型,用于在非线性非定常气动载荷下,输出飞行器的实际姿态角和实际角速度;
所述外反馈回路,用于基于期望姿态角和所述实际姿态角进行处理,获得参考角速度;
所述内反馈回路,用于对所述实际角速度和参考角速度进行处理,获得非线性非定常气动载荷的估计值,并输出相应的补偿控制指令至飞行器,对非线性非定常气动载荷进行修正。
2.根据权利要求1所述的飞行器非线性非定常气动载荷修正系统,其特征在于,所述外反馈回路包括外回路指令滤波器和外回路反步控制器;
所述外回路指令滤波器,用于对期望姿态角进行滤波;
所述外回路反步控制器,用于对所述实际姿态角和滤波后的期望姿态角进行处理,获得参考角速度。
3.根据权利要求1或2所述的飞行器非线性非定常气动载荷修正系统,其特征在于,所述内反馈回路包括内回路指令滤波器、自适应控制器和内回路反步控制器;
所述内回路指令滤波器,用于对所述参考角速度进行滤波;
所述自适应控制器,用于对所述实际角速度和滤波后的参考角速度进行处理,获得非线性非定常气动载荷的估计值;
所述内回路反步控制器,用于基于所述实际角速度、滤波后的参考角速度和所述非线性非定常气动载荷的估计值,将自身重构为自适应摄动补偿控制器,输出相应的补偿控制指令。
4.根据权利要求3所述的飞行器非线性非定常气动载荷修正系统,其特征在于,所述飞行器模型,描述为:
其中,输入变量u=[δa,δe,δr,δy,δz]T,δa表示副翼控制输入,δe表示升降舵控制输入,δr表示方向舵控制输入,δy表示横侧向推力矢量输入,δz表示纵向推理矢量输入;外回路状态变量x1=[α β φ]T,α表示迎角,β表示偏航角,φ表示滚转角;内回路状态变量x2=[p q r]T,p表示俯仰角速度,q表示偏航角速度,r表示滚转角速度;中的各元素均为气动参数,δClβ是横滚静稳定性导数,δCm是纵向静稳定性导数,是纵向过失速导数,δCnβ是航向静稳定性导数;
f2的表达式如下所示:
其中,fp(·)、fq(·)和fr(·)分别表示俯仰通道、偏航通道和滚转通道的非线性状态函数,记x0,y0,z0分别为飞机质心在地轴系三个坐标轴的投影,则I11,I22,I33分别表示机身绕x0轴的转动惯量、机身绕y0轴的转动惯量和机身绕z0轴的转动惯量,I13表示机身在x0-z0轴上的惯性积,分别为俯仰气动力矩、偏航气动力矩和滚转气动力矩;
角速率通道气动参数摄动模型,表示为:
θ2=θ2,0+Δθ2
上式中各项均为参数向量,各参数向量中的各元素表达式如下所示:
其中,带下标r的气动系数为飞机的静态气动参数,方程组等号左侧为整体气动参数,κi,i=1,2,3,4,5是气动参数摄动系数,当摄动为0时,气动参数摄动系数值为0;当摄动不为0时,满足1≤|κi|≤10,整体气动参数为静态气动参数与摄动系数的和。
6.根据权利要求5所述的飞行器非线性非定常气动载荷修正系统,其特征在于,所述外回路指令滤波器,用于对期望姿态角进行滤波,包括:
所述外回路指令滤波器,用于对滤波前的期望姿态角x1d进行滤波处理,获得滤波后的飞行器的期望姿态角x1c。
8.根据权利要求7所述的飞行器非线性非定常气动载荷修正系统,其特征在于,所述内回路指令滤波器,用于对所述参考角速度进行滤波,包括:
所述内回路指令滤波器,用于对参考角速度x2d,即第2级滤波前的期望状态变量,进行第2级指令滤波处理,即内回路指令滤波处理,处理后,获得滤波后的参考角速度x2c,即内回路状态变量滤波后的期望值。
10.根据权利要求9所述的飞行器非线性非定常气动载荷修正系统,其特征在于,所述内回路反步控制器,用于基于所述实际角速度、滤波后的参考角速度和所述非线性非定常气动载荷的估计值,将自身重构为自适应摄动补偿控制器,并输出补偿控制指令,对非线性非定常气动载荷进行修正,包括:
所述内回路反步控制器,对非线性非定常气动载荷的估计值、角速度误差和滤波后的参考角速度的导数进行处理,将自身重构为自适应摄动补偿控制器:
其中,k2为比例控制参数;u为补偿控制指令,即控制输出,表示控制器期望控制输入产生的舵面偏转角;
所述自适应摄动补偿控制器,发出补偿控制指令,对飞行器进行自适应摄动补偿控制,完成对所述非线性非定常气动载荷的修复。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211053328.0A CN115328185B (zh) | 2022-08-30 | 2022-08-30 | 一种飞行器非线性非定常气动载荷修正系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211053328.0A CN115328185B (zh) | 2022-08-30 | 2022-08-30 | 一种飞行器非线性非定常气动载荷修正系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115328185A true CN115328185A (zh) | 2022-11-11 |
CN115328185B CN115328185B (zh) | 2023-05-30 |
Family
ID=83927820
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211053328.0A Active CN115328185B (zh) | 2022-08-30 | 2022-08-30 | 一种飞行器非线性非定常气动载荷修正系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115328185B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116643501A (zh) * | 2023-07-18 | 2023-08-25 | 湖南大学 | 稳定性约束下的空中作业机器人可变阻抗控制方法及系统 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080097658A1 (en) * | 2004-11-08 | 2008-04-24 | Shue Shyhpyng J | Flight Control System Having a Three Control Loop Design |
US20130311013A1 (en) * | 2012-05-16 | 2013-11-21 | Optical Air Data Systems, Llc | Measurement Assisted Aerodynamic State Estimator |
CN103488814A (zh) * | 2013-08-16 | 2014-01-01 | 北京航空航天大学 | 一种适用于再入飞行器姿态控制的闭环仿真系统 |
US20170068252A1 (en) * | 2014-05-30 | 2017-03-09 | SZ DJI Technology Co., Ltd. | Aircraft attitude control methods |
CN107065539A (zh) * | 2017-03-14 | 2017-08-18 | 南京航空航天大学 | 一种飞翼布局飞行器的操纵面故障自适应容错控制方法 |
US20180029695A1 (en) * | 2015-02-23 | 2018-02-01 | Arif Mir Jalal ogly PASHAYEV | Method and Device for Increasing the Stability and Maneuverability of Unmanned Aerial Vehicles (UAV) using a Gyroscopic Effect |
CN107844123A (zh) * | 2017-10-11 | 2018-03-27 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种非线性飞行器航迹控制方法 |
CN109460596A (zh) * | 2018-10-29 | 2019-03-12 | 成都飞机工业(集团)有限责任公司 | 一种飞翼无人机非线性载荷计算方法 |
CN109597303A (zh) * | 2018-11-29 | 2019-04-09 | 南京航空航天大学 | 一种复合式旋翼飞行器全模式飞行控制方法 |
CN111367182A (zh) * | 2020-04-23 | 2020-07-03 | 湖南云顶智能科技有限公司 | 考虑输入受限的高超声速飞行器抗干扰反步控制方法 |
CN113778129A (zh) * | 2021-09-23 | 2021-12-10 | 北京理工大学 | 一种干扰补偿的高超声速变后掠翼飞行器跟踪控制方法 |
CN115220467A (zh) * | 2022-08-24 | 2022-10-21 | 南京航空航天大学 | 一种基于神经网络增量动态逆的飞翼飞行器姿态控制方法 |
-
2022
- 2022-08-30 CN CN202211053328.0A patent/CN115328185B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080097658A1 (en) * | 2004-11-08 | 2008-04-24 | Shue Shyhpyng J | Flight Control System Having a Three Control Loop Design |
US20130311013A1 (en) * | 2012-05-16 | 2013-11-21 | Optical Air Data Systems, Llc | Measurement Assisted Aerodynamic State Estimator |
CN103488814A (zh) * | 2013-08-16 | 2014-01-01 | 北京航空航天大学 | 一种适用于再入飞行器姿态控制的闭环仿真系统 |
US20170068252A1 (en) * | 2014-05-30 | 2017-03-09 | SZ DJI Technology Co., Ltd. | Aircraft attitude control methods |
US20180029695A1 (en) * | 2015-02-23 | 2018-02-01 | Arif Mir Jalal ogly PASHAYEV | Method and Device for Increasing the Stability and Maneuverability of Unmanned Aerial Vehicles (UAV) using a Gyroscopic Effect |
CN107065539A (zh) * | 2017-03-14 | 2017-08-18 | 南京航空航天大学 | 一种飞翼布局飞行器的操纵面故障自适应容错控制方法 |
CN107844123A (zh) * | 2017-10-11 | 2018-03-27 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种非线性飞行器航迹控制方法 |
CN109460596A (zh) * | 2018-10-29 | 2019-03-12 | 成都飞机工业(集团)有限责任公司 | 一种飞翼无人机非线性载荷计算方法 |
CN109597303A (zh) * | 2018-11-29 | 2019-04-09 | 南京航空航天大学 | 一种复合式旋翼飞行器全模式飞行控制方法 |
CN111367182A (zh) * | 2020-04-23 | 2020-07-03 | 湖南云顶智能科技有限公司 | 考虑输入受限的高超声速飞行器抗干扰反步控制方法 |
CN113778129A (zh) * | 2021-09-23 | 2021-12-10 | 北京理工大学 | 一种干扰补偿的高超声速变后掠翼飞行器跟踪控制方法 |
CN115220467A (zh) * | 2022-08-24 | 2022-10-21 | 南京航空航天大学 | 一种基于神经网络增量动态逆的飞翼飞行器姿态控制方法 |
Non-Patent Citations (2)
Title |
---|
徐倩 等: "基于扩张状态观测器和指令滤波器的导弹姿态控制", 飞行力学 * |
梅瀚桐 等: "远程导弹的自适应反步姿态控制系统设计", 空天防御 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116643501A (zh) * | 2023-07-18 | 2023-08-25 | 湖南大学 | 稳定性约束下的空中作业机器人可变阻抗控制方法及系统 |
CN116643501B (zh) * | 2023-07-18 | 2023-10-24 | 湖南大学 | 稳定性约束下的空中作业机器人可变阻抗控制方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN115328185B (zh) | 2023-05-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109541941B (zh) | 一种针对垂直起降运载器主动段飞行的自适应增广抗扰容错方法 | |
Yu et al. | Safe control of trailing UAV in close formation flight against actuator fault and wake vortex effect | |
CN109189087B (zh) | 一种垂直起降重复使用运载器的自适应容错控制方法 | |
CN105629734B (zh) | 一种近空间飞行器的轨迹跟踪控制方法 | |
CN110347170B (zh) | 可重复使用运载器再入段鲁棒容错制导控制系统及工作方法 | |
CN110850887B (zh) | 一种四旋翼无人机复合动态逆抗干扰姿态控制方法 | |
CN109144084A (zh) | 一种基于固定时间收敛观测器的垂直起降重复使用运载器姿态跟踪控制方法 | |
CN102880060A (zh) | 再入飞行器自适应指数时变滑模姿态控制方法 | |
CN106773679A (zh) | 一种基于角速度观测器的航天器容错控制方法 | |
CN103105850A (zh) | 一种近空间飞行器故障诊断与容错控制方法 | |
CN113741188B (zh) | 执行器故障下固定翼无人机反步自适应容错控制方法 | |
CN113485304A (zh) | 一种基于深度学习故障诊断的飞行器分层容错控制方法 | |
CN114906349B (zh) | 高机动飞机升降舵效率损失故障的自适应容错控制方法 | |
CN114721266B (zh) | 飞机舵面结构性缺失故障情况下的自适应重构控制方法 | |
CN115328185B (zh) | 一种飞行器非线性非定常气动载荷修正系统 | |
CN115327916B (zh) | 一种高机动飞机气动参数摄动的自适应补偿控制方法 | |
Lu et al. | Robustness and tuning of incremental backstepping approach | |
Nguyen et al. | A multi-objective flight control approach for performance adaptive aeroelastic wing | |
Emami et al. | Disturbance observer-based adaptive neural guidance and control of an aircraft using composite learning | |
Nguyen et al. | Adaptive linear quadratic gaussian optimal control modification for flutter suppression of adaptive wing | |
CN116360258A (zh) | 基于固定时间收敛的高超声速变形飞行器抗干扰控制方法 | |
Sonneveldt et al. | Constrained adaptive backstepping flight control: Application to a nonlinear F-16/MATV model | |
CN116736692A (zh) | 一种具有延迟状态约束的四旋翼无人机滑模容错控制方法 | |
Crespo et al. | Design of a model reference adaptive controller for an unmanned air vehicle | |
CN111061283A (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 |