CN115328185B - 一种飞行器非线性非定常气动载荷修正系统 - Google Patents

一种飞行器非线性非定常气动载荷修正系统 Download PDF

Info

Publication number
CN115328185B
CN115328185B CN202211053328.0A CN202211053328A CN115328185B CN 115328185 B CN115328185 B CN 115328185B CN 202211053328 A CN202211053328 A CN 202211053328A CN 115328185 B CN115328185 B CN 115328185B
Authority
CN
China
Prior art keywords
angular velocity
aircraft
attitude angle
nonlinear
filter
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
Application number
CN202211053328.0A
Other languages
English (en)
Other versions
CN115328185A (zh
Inventor
刘义
王东豪
高晓琼
柴振达
闫贯博
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Jinghang Computing Communication Research Institute
Original Assignee
Beijing Jinghang Computing Communication Research Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Jinghang Computing Communication Research Institute filed Critical Beijing Jinghang Computing Communication Research Institute
Priority to CN202211053328.0A priority Critical patent/CN115328185B/zh
Publication of CN115328185A publication Critical patent/CN115328185A/zh
Application granted granted Critical
Publication of CN115328185B publication Critical patent/CN115328185B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/08Control of attitude, i.e. control of roll, pitch, or yaw
    • G05D1/0808Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
    • G05D1/0816Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability
    • G05D1/0825Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability using mathematical models
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/10Simultaneous control of position or course in three dimensions
    • G05D1/101Simultaneous control of position or course in three dimensions specially adapted for aircraft
    • G05D1/106Change initiated in response to external conditions, e.g. avoidance of elevated terrain or of no-fly zones
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling 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样条逼近突变饱和的式子,再用方向控制推进装置即喷工质舵调节器建立硬件结构舵操纵回路里的虚飞控算子,使机械舵在舵偏角很小的情况下依然能实现高机动姿控;针对干扰条件下的鲁棒稳定问题,科学家设计了估计算法有源自动修复的想法,制造额外的估计回路和算法,对大幅值低频连续特征或小幅值高频离散特征的多干扰采取平行估计策略,继而让多个干扰观测值导进综控机重构操控软件,抵消干扰危害,实时自愈合飞行器。然而,上述系统都没考虑参数摄动情况,更没有针对飞行器非线性非定常气动载荷设计专门的补偿系统。现有摄动补偿控制算法则更多聚焦于通用动态系统。例如,采用博弈论思想解决互联系统摄动补偿问题;采用分段李雅普诺夫系统证明带摄动的分段线性系统的全局稳定性;用摄动与输入函数的位置关系将摄动划分为乘性摄动和加性摄动,并用经典鲁棒系统直接被动修复加性摄动;针对多输入多输出系统中的单通道摄动,用另一无摄动通道错位补偿单摄动对整体系统的影响,这些都是理论上常见的补偿系统。而现有针对工业工程现场偏应用的补偿系统也很多,例如,针对磨矿机执行机构摄动,为保障矿物颗粒粒径分布吻合期望随机分布函数,须首先建立非高斯随机分布系统和摄动模型,进而设计随机容错控制器修复执行器故障。
综上,目前的现有技术中没有针对飞行器非线性非定常气动载荷的补偿系统,存在无法自动修复飞行器的非线性非定常气动载荷、飞行稳定性和控制精度差等缺陷。
发明内容
鉴于上述的分析,本发明实施例旨在提供一种飞行器非线性非定常气动载荷修正系统,用以解决现有技术中没有针对飞行器非线性非定常气动载荷的补偿系统,存在无法自动修复飞行器的非线性非定常气动载荷、飞行稳定性和控制精度差的问题。
本发明的目的主要是通过以下技术方案实现的:
本发明实施例提供了一种飞行器非线性非定常气动载荷修正系统,包括:
飞行器模型和内、外双反馈回路;其中,
所述飞行器模型,用于在非线性非定常气动载荷下,输出实际姿态角和实际角速度;
所述外反馈回路,用于对期望姿态角和所述实际姿态角进行处理,获得参考角速度;
所述内反馈回路,用于对所述实际角速度和参考角速度进行处理,获得非线性非定常气动载荷的估计值,并输出补偿控制指令,对非线性非定常气动载荷进行修正。
基于上述方法的进一步改进,所述外反馈回路包括外回路指令滤波器和外回路反步控制器;
所述外回路指令滤波器,用于对期望姿态角进行滤波;
所述外回路反步控制器,用于对所述实际姿态角和滤波后的期望姿态角进行处理,获得参考角速度。
基于上述方法的进一步改进,所述内反馈回路包括内回路指令滤波器、自适应控制器和内回路反步控制器;
所述内回路指令滤波器,用于对所述参考角速度进行滤波;
所述自适应控制器,用于对所述实际角速度和滤波后的参考角速度进行处理,获得非线性非定常气动载荷的估计值;
所述内回路反步控制器,用于基于所述实际角速度、滤波后的参考角速度和所述非线性非定常气动载荷的估计值,将自身重构为自适应摄动补偿控制器,并输出补偿控制指令,对非线性非定常气动载荷进行修正。
基于上述方法的进一步改进,所述飞行器模型,描述为:
Figure BDA0003824615460000041
其中,输入变量u=[δaeryz]T,δa表示副翼控制输入,δe表示升降舵控制输入,δr表示方向舵控制输入,δy表示横侧向推力矢量输入,δz表示纵向推理矢量输入;外回路状态变量x1=[α β φ]T,α表示迎角,β表示偏航角,φ表示滚转角;内回路状态变量x2=[pq r]T,p表示俯仰角速度,q表示偏航角速度,r表示滚转角速度;
Figure BDA0003824615460000042
中的各元素均为气动参数,δC是横滚静稳定性导数,δCm是纵向静稳定性导数,/>
Figure BDA0003824615460000043
是纵向过失速导数,δC是航向静稳定性导数;
f2的表达式如下所示:
Figure BDA0003824615460000044
其中,fp(·)、fq(·)和fr(·)分别表示俯仰通道、偏航通道和滚转通道的非线性状态函数,记x0,y0,z0分别为飞机质心在地轴系三个坐标轴的投影,则I11,I22,I33分别表示机身绕x0轴的转动惯量、机身绕y0轴的转动惯量和机身绕z0轴的转动惯量,I13表示机身在x0-z0轴上的惯性积,
Figure BDA0003824615460000045
分别为俯仰气动力矩、偏航气动力矩和滚转气动力矩;
G1和G2为控制效率分配矩阵;
Figure BDA0003824615460000046
为气动和结构复合参数矩阵;
角速率通道气动参数摄动模型,表示为:
θ2=θ2,0+Δθ2
上式中各项均为参数向量,各参数向量中的各元素表达式如下所示:
Figure BDA0003824615460000051
其中,带下标r的气动系数为飞机的静态气动参数,方程组等号左侧为整体气动参数,κi,i=1,2,3,4,5是气动参数摄动系数,当摄动为0时,气动参数摄动系数值为0;当摄动不为0时,满足1≤|κi|≤10,整体气动参数为静态气动参数与摄动系数的和。
基于上述方法的进一步改进,所述外回路指令滤波器和内回路指令滤波器,包括:
所述外回路指令滤波器为第1级滤波器,所述内回路指令滤波器为第2级滤波器,即i=1,2;第i级滤波器的形式如下:
Figure BDA0003824615460000052
其中,ωic为第i级滤波器频率,
Figure BDA0003824615460000053
为第i级滤波器阻尼比,二者均为率变频器待选定的参数,xid为第i级滤波前的期望状态变量,zi1和zi2的初值选取为zi1(0)=xid(0);
滤波后的期望状态变量,满足:
Figure BDA0003824615460000054
其中,xic为滤波后的期望状态变量。
基于上述方法的进一步改进,所述外回路指令滤波器,用于对期望姿态角进行滤波,包括:
所述外回路指令滤波器,用于对滤波前的期望姿态角x1d进行滤波处理,获得滤波后的飞行器的期望姿态角x1c
基于上述方法的进一步改进,所述外回路反步控制器,用于对所述实际姿态角和滤波后的期望姿态角进行处理,获得参考角速度,包括:
将经外回路指令滤波后的期望姿态角x1c,与实际姿态角x1做差,得到姿态角误差e1=x1-x1c
对姿态角误差和滤波后的期望姿态角的导数,按照如下方式进行反步动态逆控制处理,获得参考角速度x2d,即第2级滤波前的期望状态变量:
Figure BDA0003824615460000061
其中,k1>0为外回路比例控制参数。
基于上述方法的进一步改进,所述内回路指令滤波器,用于对所述参考角速度进行滤波,包括:
所述内回路指令滤波器,用于对参考角速度x2d,即第2级滤波前的期望状态变量,进行第2级指令滤波处理,即内回路指令滤波处理,处理后,获得滤波后的参考角速度x2c,即内回路状态变量滤波后的期望值。
基于上述方法的进一步改进,所述自适应控制器,用于对所述实际角速度和滤波后的参考角速度进行处理,获得非线性非定常气动载荷的估计值,包括:
将实际角速度与滤波后的参考角速度做差,得到角速度误差e2=x2-x2c
通过以下方式,获得非线性非定常气动载荷估计值的导数:
Figure BDA0003824615460000071
其中,Γ2是自适应学习率,υ2为跟踪补偿误差,且满足υ2=e2
对非线性非定常气动载荷估计值的导数进行积分,获得非线性非定常气动载荷估计值
Figure BDA0003824615460000072
基于上述方法的进一步改进,所述内回路反步控制器,用于基于所述实际角速度、滤波后的参考角速度和所述非线性非定常气动载荷的估计值,将自身重构为自适应摄动补偿控制器,并输出补偿控制指令,对非线性非定常气动载荷进行修正,包括:
所述内回路反步控制器,对非线性非定常气动载荷的估计值、角速度误差和滤波后的参考角速度的导数进行处理,将自身重构为自适应摄动补偿控制器:
Figure BDA0003824615460000073
其中,k2为比例控制参数;u为补偿控制指令,即控制输出,表示控制器期望控制输入产生的舵面偏转角;
所述自适应摄动补偿控制器,发出补偿控制指令,对飞行器进行自适应摄动补偿控制,完成对所述非线性非定常气动载荷的修复。
与现有技术相比,本发明至少可实现如下有益效果之一:
1、本发明提供了一种飞行器非线性非定常气动载荷修正系统,使得飞行器在飞行中即使出现非线性非定常气动载荷也能保持稳定。
2、本发明在公开一个自适应摄动补偿控制器的同时,也提供了一个自适应控制器,该自适应控制器估计得到的准确摄动信息不仅有助于修复摄动,也有助于驾驶员基于复杂多元任务做综合决策。
3、本发明是在标称动态逆控制的基础上设计的自适应补偿控制算法,具有动态逆控制的优点,即能解决非线性问题、结构简单易实现、可靠性高等,为剩余通道设计其它复杂控制器预留了设计空间。
本发明中,上述各技术方案之间还可以相互组合,以实现更多的优选组合方案。本发明的其他特征和优点将在随后的说明书中阐述,并且,部分优点可从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过说明书以及附图中所特别指出的内容中来实现和获得。
附图说明
附图仅用于示出具体实施例的目的,而并不认为是对本发明的限制,在整个附图中,相同的参考符号表示相同的部件。
图1为为本发明实施例中控制系统的架构图。
具体实施方式
下面结合附图来具体描述本发明的优选实施例,其中,附图构成本申请一部分,并与本发明的实施例一起用于阐释本发明的原理,并非用于限定本发明的范围。
本实施例公开了一种飞行器非线性非定常气动载荷修正系统,如图1所示,包括飞行器模型和内、外双反馈回路;其中,
所述飞行器模型,用于在非线性非定常气动载荷下,输出飞行器的实际姿态角和实际角速度;
所述外反馈回路,用于基于期望姿态角和所述实际姿态角进行处理,获得参考角速度;
所述内反馈回路,用于对所述实际角速度和参考角速度进行处理,获得非线性非定常气动载荷的估计值,并输出相应的补偿控制指令至飞行器,对非线性非定常气动载荷进行修正。
具体来说,所述飞行器模型,表示为:
Figure BDA0003824615460000091
其中,输入变量u=[δaeryz]T,δa表示副翼控制输入,δe表示升降舵控制输入,δr表示方向舵控制输入,δy表示横侧向推力矢量输入,δz表示纵向推理矢量输入;外回路状态变量x1=[α β φ]T,α表示迎角,β表示偏航角,φ表示滚转角;内回路状态变量x2=[pq r]T,p表示俯仰角速度,q表示偏航角速度,r表示滚转角速度;
Figure BDA0003824615460000092
中的各元素均为气动参数,δC是横滚静稳定性导数,δCm是纵向静稳定性导数,/>
Figure BDA0003824615460000093
是纵向过失速导数,δC是航向静稳定性导数。
f2的表达式如下所示:
Figure BDA0003824615460000094
其中,fp(·)、fq(·)和fr(·)分别表示俯仰通道、偏航通道和滚转通道的非线性状态函数,设x0,y0,z0分别为飞行器质心在地轴系三个坐标轴的投影,则I11,I22,I33分别表示机身绕x0轴的转动惯量、机身绕y0轴的转动惯量和机身绕z0轴的转动惯量,I13表示机身在x0-z0轴上的惯性积,
Figure BDA0003824615460000101
分别为俯仰气动力矩、偏航气动力矩和滚转气动力矩。
G1、G2
Figure BDA0003824615460000102
的表达式如下所示:
Figure BDA0003824615460000103
Figure BDA0003824615460000104
Figure BDA0003824615460000105
其中,G1和G2为控制效率分配矩阵,
Figure BDA0003824615460000106
表示副翼在俯仰通道的控制分配系数,
Figure BDA0003824615460000107
表示方向舵在俯仰通道的控制分配系数,/>
Figure BDA0003824615460000108
表示推力矢量横侧向舵在俯仰通道的控制分配系数,/>
Figure BDA0003824615460000109
表示升降舵在偏航通道的控制分配系数,/>
Figure BDA00038246154600001010
表示推力矢量纵向舵在偏航通道的控制分配系数,/>
Figure BDA00038246154600001011
表示副翼在滚转通道的控制分配系数,/>
Figure BDA00038246154600001012
表示方向舵在滚转通道的控制分配系数,/>
Figure BDA00038246154600001013
表示推力矢量横侧向舵在滚转通道的控制分配系数;/>
Figure BDA00038246154600001014
为气动和结构复合参数矩阵;其中,c3,c4,c7和c9是气动常数,qd是动压,b是翼展,S是参考机翼面积,c是平均气动弦长,/>
Figure BDA00038246154600001015
中的物理量均通过风洞实验测定。
角速率通道气动参数摄动模型,表示为:
Figure BDA00038246154600001016
上式中各项均为参数向量,各参数向量中的各元素表达式如下所示:
Figure BDA0003824615460000111
/>
其中,带下标r的气动系数为飞行器的静态气动参数,方程组等号左侧为整体气动参数,κi,i=1,2,3,4,5是气动参数摄动系数,当摄动为0时,气动参数摄动系数值为0,则静态气动参数与整体气动参数相等;当摄动不为0时,满足1≤|κi|≤10,整体气动参数为静态气动参数与摄动系数的和。
采用这种形式的气动参数摄动,主要出于两方面的考虑。其一是为了模仿大迎角下的气动参数的迟滞效应。飞行器气动参数的迟滞效应要是受到迎角的角速度
Figure BDA0003824615460000112
的影响,
Figure BDA0003824615460000113
的值越大,流场的迟滞效应也就越明显。而飞行器的迟滞效应主要在大迎角状态下发生,迎角α的大小对于迟滞环的大小也有显著的影响。所以在考虑迟滞环的大小时,用/>
Figure BDA0003824615460000114
来对气动参数进行修正。
另外一个考虑就是检验控制系统的鲁棒性。由于飞行器在大迎角飞行时,特别是在进行过失速机动时,流场非常复杂,要非常准确的获得飞行器的气动参数几乎不可能,所以气动参数误差的存在是难免的,而且有时候还非常大。而飞行器的控制律是基于试飞(或风洞实验)得到的气动参数而设计的,所以要使控制系统可用的条件就是飞行器控制系统具有较好的鲁棒性。所以在控制律设计时使用静态的气动参数,而在计算飞行器机动过程中的气动力和力矩时,利用
Figure BDA0003824615460000115
对其进行修正。
所述外反馈回路包括外回路指令滤波器和外回路反步控制器;
所述外回路指令滤波器,用于对期望姿态角进行滤波;
所述外回路反步控制器,用于对所述实际姿态角和滤波后的期望姿态角进行处理,获得参考角速度。
所述内反馈回路包括内回路指令滤波器、自适应控制器和内回路反步控制器;
所述内回路指令滤波器,用于对所述参考角速度进行滤波;
所述自适应控制器,用于对所述实际角速度和滤波后的参考角速度进行处理,获得非线性非定常气动载荷的估计值;
所述内回路反步控制器,用于基于所述实际角速度、滤波后的参考角速度和所述非线性非定常气动载荷的估计值,将自身重构为自适应摄动补偿控制器,并输出相应的补偿控制指令,对非线性非定常气动载荷进行修正。
所述飞行器模型输出外回路反馈,得到飞行器的实际姿态角,即外回路状态变量x1=[α β φ]T
所述外回路指令滤波器,以滤波前的期望姿态角为输入,输出滤波后的期望姿态角。
如图1所示,所述外回路指令滤波器为第1级滤波器,所述内回路指令滤波器为第2级滤波器;第i级滤波器的形式如下,i=1,2:
Figure BDA0003824615460000121
其中,ωic为第i级滤波器频率,
Figure BDA0003824615460000122
为第i级滤波器阻尼比,二者均为率变频器待选定的参数,xid为第i级滤波前的期望状态变量,zi1和zi2表示第i级滤波器状态变量,初值选取为zi1(0)=xid(0)。
滤波后的期望状态变量,满足:
Figure BDA0003824615460000123
其中,xic为第i级滤波器滤波后的期望状态变量。
此处的二阶指令滤波器,将经典的自适应反步控制中对中间信号的求导替换为求解二阶指令滤波器的状态变量,避免了对复杂运动学方程进行求导和微分爆炸,减小了计算量。
具体来说,所述外回路指令滤波器,用于对滤波前的期望姿态角x1d进行滤波处理,获得滤波后的飞行器的期望姿态角x1c
x1c是经外回路指令滤波后的期望姿态角,x1是实际姿态角,则姿态角误差e1=x1-x1c
具体来说,所述外回路反步控制器,用于对姿态角误差和滤波后的期望姿态角的导数,按照如下方式进行反步动态逆控制处理,获得参考角速度x2d,即第2级滤波前的期望状态变量:
Figure BDA0003824615460000131
其中,k1>0为外回路比例控制参数。
所述飞行器模型输出内回路反馈,得到飞行器的实际角速度,即内回路状态变量x2=[p q r]T
具体来说,所述内回路指令滤波器,用于对参考角速度x2d,即第2级滤波前的期望状态变量,进行第2级指令滤波处理,即内回路指令滤波处理,处理后,获得滤波后的参考角速度x2c,即内回路状态变量滤波后的期望值。
将实际角速度与滤波后的参考角速度做差,得到角速度误差e2=x2-x2c
需要说明的是,对非线性非定常气动载荷进行自适应估计,需要利用跟踪补偿误差,而跟踪补偿误差即角速度误差:
υ2=e2
其中,υ2为跟踪补偿误差。
具体来说,所述自适应控制器,用于通过以下方式,获得非线性非定常气动载荷估计值的导数:
Figure BDA0003824615460000141
其中,Γ2是自适应学习率。
对非线性非定常气动载荷估计值的导数进行积分,获得非线性非定常气动载荷估计值
Figure BDA0003824615460000142
具体来说,所述内回路反步控制器,用于对非线性非定常气动载荷的估计值、角速度误差和滤波后的参考角速度的导数进行处理,将自身重构为自适应摄动补偿控制器:
Figure BDA0003824615460000143
其中,k2为比例控制参数;u为补偿控制指令,即控制输出,表示控制器期望控制输入产生的舵面偏转角。
所述自适应摄动补偿控制器,用于发出相应的补偿控制指令,对飞行器进行自适应摄动补偿控制,完成对所述非线性非定常气动载荷的修复。
本领域技术人员可以理解,实现上述实施例方法的全部或部分流程,可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于计算机可读存储介质中。其中,所述计算机可读存储介质为磁盘、光盘、只读存储记忆体或随机存储记忆体等。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。

Claims (9)

1.一种飞行器非线性非定常气动载荷修正系统,其特征在于,包括飞行器模型和内反馈回路、外反馈回路;其中,
所述飞行器模型,用于在非线性非定常气动载荷下,输出飞行器的实际姿态角和实际角速度;
所述外反馈回路,用于基于期望姿态角和所述实际姿态角进行处理,获得参考角速度;
所述内反馈回路,用于对所述实际角速度和参考角速度进行处理,获得非线性非定常气动载荷的估计值,并输出相应的补偿控制指令至飞行器,对非线性非定常气动载荷进行修正,所述内反馈回路包括内回路指令滤波器、自适应控制器和内回路反步控制器;
所述内回路指令滤波器,用于对所述参考角速度进行滤波;
所述自适应控制器,用于对所述实际角速度和滤波后的参考角速度进行处理,获得非线性非定常气动载荷的估计值;
所述内回路反步控制器,用于基于所述实际角速度、滤波后的参考角速度和所述非线性非定常气动载荷的估计值,将自身重构为自适应摄动补偿控制器,输出相应的补偿控制指令。
2.根据权利要求1所述的飞行器非线性非定常气动载荷修正系统,其特征在于,所述外反馈回路包括外回路指令滤波器和外回路反步控制器;
所述外回路指令滤波器,用于对期望姿态角进行滤波;
所述外回路反步控制器,用于对所述实际姿态角和滤波后的期望姿态角进行处理,获得参考角速度。
3.根据权利要求1所述的飞行器非线性非定常气动载荷修正系统,其特征在于,所述飞行器模型,描述为:
Figure FDA0004191784190000021
其中,输入变量u=[δaeryz]T,δa表示副翼控制输入,δe表示升降舵控制输入,δr表示方向舵控制输入,δy表示横侧向推力矢量输入,δz表示纵向推理矢量输入;外回路状态变量x1=[α β φ]T,α表示迎角,β表示偏航角,φ表示滚转角;内回路状态变量x2=[p q r]T,p表示俯仰角速度,q表示偏航角速度,r表示滚转角速度;
Figure FDA0004191784190000022
中的各元素均为气动参数,δC是横滚静稳定性导数,δCm是纵向静稳定性导数,/>
Figure FDA0004191784190000023
是纵向过失速导数,δC是航向静稳定性导数;
f2的表达式如下所示:
Figure FDA0004191784190000024
其中,fp(·)、fq(·)和fr(·)分别表示俯仰通道、偏航通道和滚转通道的非线性状态函数,记x0,y0,z0分别为飞机质心在地轴系三个坐标轴的投影,则I11,I22,I33分别表示机身绕x0轴的转动惯量、机身绕y0轴的转动惯量和机身绕z0轴的转动惯量,I13表示机身在x0-z0轴上的惯性积,
Figure FDA0004191784190000025
分别为俯仰气动力矩、偏航气动力矩和滚转气动力矩;
G1和G2为控制效率分配矩阵;
Figure FDA0004191784190000026
为气动和结构复合参数矩阵;
角速率通道气动参数摄动模型,表示为:
θ2=θ2,0+Δθ2
上式中各项均为参数向量,各参数向量中的各元素表达式如下所示:
Figure FDA0004191784190000031
其中,带下标r的气动系数为飞机的静态气动参数,方程组等号左侧为整体气动参数,κi,i=1,2,3,4,5是气动参数摄动系数,当摄动为0时,气动参数摄动系数值为0;当摄动不为0时,满足1≤|κi|≤10,整体气动参数为静态气动参数与摄动系数的和。
4.根据权利要求3所述的飞行器非线性非定常气动载荷修正系统,其特征在于,所述外回路指令滤波器为第1级滤波器,所述内回路指令滤波器为第2级滤波器,第i级滤波器的形式如下,i=1,2:
Figure FDA0004191784190000032
其中,ωic为第i级滤波器频率,
Figure FDA0004191784190000033
为第i级滤波器阻尼比,二者均为率变频器待选定的参数,xid为第i级滤波前的期望状态变量,zi1和zi2表示第i级滤波器状态变量,初值选取为zi1(0)=xid(0);
滤波后的期望状态变量,满足:
Figure FDA0004191784190000034
其中,xic为第i级滤波器滤波后的期望状态变量。
5.根据权利要求4所述的飞行器非线性非定常气动载荷修正系统,其特征在于,所述外回路指令滤波器,用于对期望姿态角进行滤波,包括:
所述外回路指令滤波器,用于对滤波前的期望姿态角x1d进行滤波处理,获得滤波后的飞行器的期望姿态角x1c
6.根据权利要求5所述的飞行器非线性非定常气动载荷修正系统,其特征在于,所述外回路反步控制器,用于对所述实际姿态角和滤波后的期望姿态角进行处理,获得参考角速度,包括:
将经外回路指令滤波后的期望姿态角x1c,与实际姿态角x1做差,得到姿态角误差e1=x1-x1c
对姿态角误差和滤波后的期望姿态角的导数,按照如下方式进行反步动态逆控制处理,获得参考角速度x2d,即第2级滤波前的期望状态变量:
Figure FDA0004191784190000041
其中,k1>0为外回路比例控制参数。
7.根据权利要求6所述的飞行器非线性非定常气动载荷修正系统,其特征在于,所述内回路指令滤波器,用于对所述参考角速度进行滤波,包括:
所述内回路指令滤波器,用于对参考角速度x2d,即第2级滤波前的期望状态变量,进行第2级指令滤波处理,即内回路指令滤波处理,处理后,获得滤波后的参考角速度x2c,即内回路状态变量滤波后的期望值。
8.根据权利要求7所述的飞行器非线性非定常气动载荷修正系统,其特征在于,所述自适应控制器,用于对所述实际角速度和滤波后的参考角速度进行处理,获得非线性非定常气动载荷的估计值,包括:
将实际角速度与滤波后的参考角速度做差,得到角速度误差e2=x2-x2c
通过以下方式,获得非线性非定常气动载荷估计值的导数:
Figure FDA0004191784190000042
其中,Γ2是自适应学习率,υ2为跟踪补偿误差,且满足υ2=e2
对非线性非定常气动载荷估计值的导数进行积分,获得非线性非定常气动载荷估计值
Figure FDA0004191784190000051
9.根据权利要求8所述的飞行器非线性非定常气动载荷修正系统,其特征在于,所述内回路反步控制器,用于基于所述实际角速度、滤波后的参考角速度和所述非线性非定常气动载荷的估计值,将自身重构为自适应摄动补偿控制器,并输出补偿控制指令,对非线性非定常气动载荷进行修正,包括:
所述内回路反步控制器,对非线性非定常气动载荷的估计值、角速度误差和滤波后的参考角速度的导数进行处理,将自身重构为自适应摄动补偿控制器:
Figure FDA0004191784190000052
其中,k2为比例控制参数;u为补偿控制指令,即控制输出,表示控制器期望控制输入产生的舵面偏转角;
所述自适应摄动补偿控制器,发出补偿控制指令,对飞行器进行自适应摄动补偿控制,完成对所述非线性非定常气动载荷的修正。
CN202211053328.0A 2022-08-30 2022-08-30 一种飞行器非线性非定常气动载荷修正系统 Active CN115328185B (zh)

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 CN115328185A (zh) 2022-11-11
CN115328185B true 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)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116643501B (zh) * 2023-07-18 2023-10-24 湖南大学 稳定性约束下的空中作业机器人可变阻抗控制方法及系统

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
BRPI0517266A (pt) * 2004-11-08 2008-10-07 Bell Helicopter Textron Inc sistema de controle de vÈo que tem um projeto de laço de controle triplo
US20130311013A1 (en) * 2012-05-16 2013-11-21 Optical Air Data Systems, Llc Measurement Assisted Aerodynamic State Estimator
CN103488814B (zh) * 2013-08-16 2017-02-15 北京航空航天大学 一种适用于再入飞行器姿态控制的闭环仿真系统
WO2015180171A1 (en) * 2014-05-30 2015-12-03 SZ DJI Technology Co., Ltd. Aircraft attitude control methods
US10737770B2 (en) * 2015-02-23 2020-08-11 Arif Mir Jalal ogly PASHAYEV Method and device for increasing the stability and maneuverability of unmanned aerial vehicles (UAV) using a gyroscopic effect
CN107065539B (zh) * 2017-03-14 2020-01-31 南京航空航天大学 一种飞翼布局飞行器的操纵面故障自适应容错控制方法
CN107844123B (zh) * 2017-10-11 2021-03-02 中国空气动力研究与发展中心计算空气动力研究所 一种非线性飞行器航迹控制方法
CN109460596B (zh) * 2018-10-29 2023-12-12 成都飞机工业(集团)有限责任公司 一种飞翼无人机非线性载荷计算方法
CN109597303B (zh) * 2018-11-29 2020-09-18 南京航空航天大学 一种复合式旋翼飞行器全模式飞行控制方法
CN111367182A (zh) * 2020-04-23 2020-07-03 湖南云顶智能科技有限公司 考虑输入受限的高超声速飞行器抗干扰反步控制方法
CN113778129B (zh) * 2021-09-23 2023-09-19 北京理工大学 一种干扰补偿的高超声速变后掠翼飞行器跟踪控制方法
CN115220467A (zh) * 2022-08-24 2022-10-21 南京航空航天大学 一种基于神经网络增量动态逆的飞翼飞行器姿态控制方法

Also Published As

Publication number Publication date
CN115328185A (zh) 2022-11-11

Similar Documents

Publication Publication Date Title
Yu et al. Safe control of trailing UAV in close formation flight against actuator fault and wake vortex effect
Nguyen et al. Aeroelastic modeling of elastically shaped aircraft concept via wing shaping control for drag reduction
CN107065539B (zh) 一种飞翼布局飞行器的操纵面故障自适应容错控制方法
Weisshaar et al. Induced drag reduction using aeroelastic tailoring with adaptive control surfaces
CN104298109B (zh) 基于多控制器融合的无尾飞行器协调转弯控制方法
Stojakovic et al. Minimal safe speed of the asymmetrically loaded combat airplane
Virgilio Pereira et al. Model predictive control architectures for maneuver load alleviation in very flexible aircraft
CN114721266B (zh) 飞机舵面结构性缺失故障情况下的自适应重构控制方法
CN115328185B (zh) 一种飞行器非线性非定常气动载荷修正系统
Nguyen et al. Multi-objective flight control for drag minimization and load alleviation of high-aspect ratio flexible wing aircraft
Ferrier et al. Active gust load alleviation of high-aspect ratio flexible wing aircraft
CN114906349A (zh) 高机动飞机升降舵效率损失故障的自适应容错控制方法
Raj et al. Iterative learning based feedforward control for transition of a biplane-quadrotor tailsitter UAS
Nguyen et al. A multi-objective flight control approach for performance adaptive aeroelastic wing
CN115327916B (zh) 一种高机动飞机气动参数摄动的自适应补偿控制方法
CN114637203A (zh) 一种针对中高速、大机动无人机的飞行控制系统
Kim et al. Development of flight control law for improvement of uncommanded lateral motion of the fighter aircraft
VanZwieten et al. Design of a prototype ocean current turbine—Part II: flight control system
Nguyen et al. Adaptive linear quadratic gaussian optimal control modification for flutter suppression of adaptive wing
Zheng et al. Adaptive fractional order non-singular terminal sliding mode anti-disturbance control for advanced layout carrier-based UAV
Lemon et al. Model reference adaptive fight control adapted for general aviation: controller gain simulation and preliminary flight testing on a bonanza fly-by-wire testbed
Wahid et al. Comparative assesment using LQR and fuzzy logic controller for a pitch control system
CN116360258A (zh) 基于固定时间收敛的高超声速变形飞行器抗干扰控制方法
Nguyen et al. Multi-objective adaptive control for load alleviation and drag minimization of flexible aircraft
Drew et al. Multi-objective optimal control of the 6-DoF aeroservoelastic common research model with aspect ratio 13.5 wing

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