临近空间柔性飞行器输出反馈控制方法研究

王民泰, 曹云峰, 周嘉麒

PDF(2192 KB)
中国科技核心期刊      中国指挥与控制学会会刊     军事装备类重点期刊
指挥控制与仿真 ›› 2022, Vol. 44 ›› Issue (4) : 105-111. DOI: 10.3969/j.issn.1673-3819.2022.04.018
工程实践

临近空间柔性飞行器输出反馈控制方法研究

作者信息 +

Research on Output Feedback Control Method of Near Space Flexible Aircraft

Author information +
History +

摘要

柔性飞行器(Flexible Flying Wing Aircraft,FFWA)因具有高升力系数、低阻力系数等优点而受到广泛关注。受到结构动力学和刚体动力学的耦合作用,飞行器容易发生弹性形变,给柔性飞行器的建模造成困难,也给控制系统的设计带来挑战。为此,针对飞翼布局柔性飞行器模型,研究了一种自适应输出反馈控制方法,并在飞翼布局柔性飞行器上进行仿真验证。首先,搭建飞翼布局柔性飞行器纵向动力学模型,并进行配平;然后,对非线性纵向动力学模型线性化得到线性飞行器模型;然后,设计包括基于观测器的基线设计和基于修正闭环参考模型的参数自适应在内的自适应输出反馈控制器,总体上保证了存在不确定参数时的稳定性和跟踪性能;最后,验证了输出反馈控制器在加入修正闭环参考模型时的优越性能。仿真结果表明,运用输出反馈控制能够很好地处理飞行器模型存在不确定参数的问题,并且加入修正闭环参考模型能够提高控制器的参数自适应性能。

Abstract

Flexible flying wing aircraft has attracted extensive attention because of its high lift coefficient and low drag coefficient. Due to the coupling of aircraft structural dynamics and rigid body dynamics, aircraft is prone to elastic deformation, which not only makes it difficult to model flexible aircraft, but also brings challenges to the design of control system. Therefore, an adaptive output feedback control method is designed for the flexible aircraft model of the flying wing configuration, and the simulation application is carried out on the flexible flying wing aircraft. Firstly, the longitudinal dynamic model of flexible flying wing aircraft is constructed and the trim analysis is carried out. Then, the nonlinear longitudinal dynamic model is linearized to obtain the linear aircraft model. Then, an adaptive output feedback controller is designed, including observer based baseline design and parameter adaptation based on modified closed-loop reference model. The overall design ensures the stability and tracking performance in the presence of large parameter uncertainty. Finally, the superior performance of the output feedback controller under the condition of parameter adaptation with modified closed-loop reference model is verified. The simulation results show that the output feedback control can well deal with the parameter uncertainty of the aircraft model, and adding the modified closed-loop reference model can improve the parameter adaptive performance of the controller.

关键词

临近空间 飞翼布局 柔性飞行器 纵向动力学模型 输出反馈控制器

Key words

near space; flying wing layout; FFWA; Longitudinal dynamics; output feedback control

引用本文

导出引用
王民泰, 曹云峰, 周嘉麒. 临近空间柔性飞行器输出反馈控制方法研究[J]. 指挥控制与仿真. 2022, 44(4): 105-111 https://doi.org/10.3969/j.issn.1673-3819.2022.04.018
WANG Min-tai, CAO Yun-feng, ZHOU Jia-qi. Research on Output Feedback Control Method of Near Space Flexible Aircraft[J]. Command Control and Simulation. 2022, 44(4): 105-111 https://doi.org/10.3969/j.issn.1673-3819.2022.04.018

1 柔性飞行器的发展

1.1 临近空间和柔性飞行器

临近空间(Near Space,简称NS)是介于航空领域与航天领域之间且尚未被充分开发利用的空间领域,又被称为近空间、亚轨道空间、空天过渡带等[1]。该空间距地面20 ~100 km,包含大气平流层、中间大气层、部分电离层三个区域[2],该高度高于传统航空设备等飞行器的最高飞行高度,但是低于卫星、空间设备等航天器的轨道运行高度[3]。临近空间具有平均风速低,气象条件稳定,太阳光辐照强度大,湿度小,阻力小,电磁特性优良等优点[4]
临近空间大气稀薄的环境致使飞行器在该区域很难获得较大的升力,因此,航空飞行器无法飞行到该空间,而无动力轨道航天器在该空间又会因为大气阻力的黏滞作用而迅速坠落。柔性飞行器因具有柔性机翼,能够获得较大的升力,因此能在该空间长期飞行,进行侦察、监视、跟踪等活动。
随着世界各国军事科技的飞速发展,低空飞行的飞行器安全性不再得到保障。为能继续执行情报收集、目标侦察、网络通信等军民活动,各国迫切需要能够避免地面炮火攻击,进行高空长时间飞行的飞行器,于是,高空长航时飞行器(High Altitude Long Endurance,简称HALE)应运而生。为了完成高难度的飞行任务,需要HALE的机翼具有大展弦比,进而获得较大的升阻比,并且,机身的轻薄特点使得机翼的柔性较大,所以,这类具有大展弦比机翼的高空长航时飞行器又被称为柔性飞行器(Flexible Aircraft,简称FA)。

1.2 飞翼布局柔性飞行器

柔性飞行器具有大展弦比特点,与普通飞行器相比,其优势主要表现为结构简单,重量轻,可高空长时间航行及多任务执行等。大柔性飞行器有以下三种布局:常规翼布局、联接翼布局[5]和飞翼布局。
常规翼布局柔性飞行器与固定翼飞行器结构类似,但其机翼相比较于固定翼飞行器机翼更加细长[6],如图1所示。其优点是展弦比大,质量轻;缺点是机身阻力大,不能完成灵活性较高的任务。
图1 柔性常规翼飞行器

Full size|PPT slide

联接翼布局柔性飞行器通过后掠前翼与前掠后翼两部分联接组成[5],如图2所示。其优点是结构强度大,飞行稳定,飞行器框架结构小,灵活性强,巡航速度快,具有良好的自然姿态恢复能力和良好的气动静稳定特性,升阻比也较高;其缺点是在结构、气动弹性、控制方面存在耦合,情况复杂,在控制律设计方面难度较大。
图2 柔性联接翼飞行器

Full size|PPT slide

飞翼布局柔性飞行器是柔性飞行器设计采用最常见、历史最悠久的布局[6],如图3所示。其优点是重量轻,飞行阻力小,升力面大,飞行获得的升力大。此外,在隐蔽性方面,相比较于其他两种布局,也有一定的优势;其缺点是飞行器纵向稳定性较差,容易导致飞行器出现静不稳定的情况。此外,其存在刚体结构和气动弹性强耦合的特点,导致其控制系统设计难度较大。
图3 飞翼型柔性飞行器

Full size|PPT slide

本文针对一类飞翼布局柔性飞行器的纵向动力学模型,设计了一种自适应输出反馈控制器。借助控制器中的观测器和修正闭环参考模型对柔性飞行器的不确定性进行估计,总体上保证了存在较大参数不确定时的稳定性和跟踪性能。

2 柔性飞行器纵向动力学模型

2.1 纵向动力学方程

飞翼布局柔性飞行器的纵向运动方程表示为[7]
mV˙=(Tcosα-D)-mgsinμmVα˙=-(Tsinα+L)+mVq+mgcosμθ˙=qIyq˙=M-Iyqh˙=Vsinμη˙=ξd3ξ˙=Mj,x-κcξ-κkη+d1-d2
(1)
其中,h为高度,航迹倾斜角μ=θ-α,并且
d1=s2m[(V˙sin(α)+Vcos(α)α˙)cos(η)Vsin(α)sin(η)η˙2s3cos(η)sin(η)η˙2]d2=(IyyIzzms212)sin(η)cos(η)q2s2mcos(η)Vcos(α)qd3=Ixx+m(s24s26cos2(η))
Iy=Iyy+2Iyycosη+ms2sin2η/6T=T1+T2+T3
(2)
式中, Ti*是每块刚体翼上螺旋桨推力,κc为扭簧阻尼系数,κk为扭簧弹性常数。

2.2 具体研究对象

沈华勋[7]用三个刚性翼和弹性铰链构造了一个柔性飞行器模型,以此来逼近柔性飞行器的复杂非线性特性。本文的研究对象是飞翼布局柔性飞行器[7],如图4所示,柔性飞行器共有3个刚性翼,每个刚性翼的长度和重量都一样,总质量为1 459 kg,每个刚性翼都有一个产生推力的螺旋桨、一个沿着主翼尾部运行的副翼和一个连接在吊杆末端的升降舵。
图4 柔性飞翼型飞行器

Full size|PPT slide

2.3 配平分析

本文对飞翼布局柔性飞行器进行配平分析[8-10]。公式(1)中的非线性动力学表示为
X˙=f(X,U)
(3)
式中,X是状态向量,U是控制输入。
X= Vαhθηη˙T
U= δacδaoδecδeoδTT
(4)
式中,acao分别表示副翼中心和副翼外侧, δac= δa2 δao= δai,其中i=1,3。下标ec,eo分别表示升降舵中心和升降舵外部。控制输入δT=Ti,i=1,2,3。线性动力学定义为
x˙=Ax+Bu+ε(X,U)
(5)
式中,x=X-X0,u=U-U0,调整状态X0和调整输入U0满足0=f(X0,U0),ε是线性化误差,AB分别表示为:
A=f(X,U)X|X=X0U=U0B=f(X,U)U|X=X0U=U0
(6)
公式(3)中的非线性柔性飞行器动力学在风速为9 m/s、高度为12 000 m时线性化,二面角在0至45°之间变化,配平结果如图5所示。
图5 柔性飞行器配平结果

Full size|PPT slide

3 柔性飞行器输出反馈控制

柔性飞行器运行在临近空间区域内,在执行飞行任务的过程中,飞行环境十分复杂,使得其气动和环境参数会发生巨大的改变;极高的飞行速度和快速的燃料消耗会导致飞行器的惯性参数和结构参数发生变化;对柔性飞行器飞行状态的建模研究不充分,以及实验数据不完备等因素,使得现有的柔性飞行器的数学模型中不可避免地存在诸多不确定性因素。这些不确定性因素会给飞行器动态特性分析和控制系统设计带来不利的影响。所以,有必要研究模型参数不确定时的飞行控制问题。

3.1 输出反馈控制发展背景

自适应输出反馈控制器因其能够控制具有不完全状态测量对象的优点而被广泛研究。通过观测器生成状态估计,并使用该估计来执行类似状态反馈的控制[11]。基于观测器的控制器已广泛应用于飞行器控制研究,对标称设备模型[12]的控制性能优异。当线性控制器设计有足够的稳定裕度[11]时,可以容忍一定量的模型不确定性。自适应控制在参数不确定时具有良好的稳定性能[13],因此可作为线性控制器的改进方法被研究。
自适应控制器是基于飞行器传递函数矩阵的经典方法,但通常需要知道飞行器的Hermite矩阵[14],并使用非最小观测器和参考模型,包含大量积分器,致使控制性能缓慢。与经典方法不同,有文献提出了一种基于状态空间表示的新方法,使用最小观测器来生成底层状态估计[15],并将状态估计用来实现反馈控制,达到参数自适应。并引入闭环参考模型(Closed-loop Reference Model,简称CRM),将最小观测器作为参考模型。
本文在具有二阶致动器的柔性飞行器动力学模型上进行仿真验证,该模型可以充分反映柔性飞行器中的参数不确定问题。通过采用带有CRM的输出反馈控制器,选用高阶滤波器代替标准的一阶微分方程,进行参数自适应的自适应律设计,设计出针对具有不确定参数的柔性飞行器的自适应输出反馈控制器。

3.2 线性柔性飞行器模型

将非线性模型围绕平衡点进行线性化,从而得出线性飞行器模型。
x˙p=Apxp+Bpθp Txp+BpΛupyp=Cpxpz=Cpcxp+Dzθp Txp+DzΛup
(7)
式中,xp Rnp是状态,upRm是控制输入,yp Rpp是测量输出,zRd是跟踪输出。通过设计u,使得z在不确定性的情况下,能够跟踪参考模型的跟踪轨迹zmΛ*是执行机构效率的对角矩阵, θp*表示状态相关的不确定性。

3.3 设计自适应输出反馈控制器

设计自适应输出反馈控制的第一步是选择参考模型,其使用修改的闭环参考模型(CRM)[16]:
x˙m=Amxm+Bzzcmd+Ley+Fa(ey)+Fb(ey)
ym=Cxm,zm=Czxm
(8)
式中,Am=A-BKT,其中,KRn×m为不包含不确定性的系统设计的标称反馈增益矩阵,使用线性二次调节器(Linear Quadratic Regulator,简称LQR)技术。Fa(ey)、Fb(ey)是稍后要设计的函数。ymy的参考输出轨迹,ey=y-ym是跟踪误差,CRM误差反馈增益L设计为
L= B1aR-1S
(9)
其中,
S=(C B1a)T, C-=SC
(10)
R-1=(C-B1a)-1[C-A B1a+(C-A B1a)T](C-B1a)-1+εI
(11)
Fa(ey)和Fb(ey)设计为:
Fa(ey)=F2(ψ31T (t) e-ψy[1])+F3(ψ31T (t) e-ψy[1][2])+F3(ψ32T (t) e-ψy[2])
(12)
Fb(ey)=-F2T(t) e-Ψ1y10[1])-F3(ψ31T (t) e-Ψ1y10[1][2])
(13)
式中,Fi(ω(t))定义为:
F2(ω(t))=( B2as+ B3a a20)[ω(t)]
F3(ω(t))= B3a π22 (s)[ω(t)]
(14)
式中, π22(s)表示s中的一个多项式,将 Bia定义为:
B1a=A2B3 a22+3 a21+B3 a20
B2a=AB3 a22+B3 a21
B3a=B3 a22
(15)
B1a, B2a是{A, Bia,C}的相关度i输入路径。 ψ31(t)和 ψ32(t)是 ψ31*Rm×m ψ32*Rm×m的估计值,也是转换后 ψ3*的元素。Λ(t)是Λ*的估计值。CRM中的FaL一起容纳了通过形成误差模型而引入的不确定性,Fb容纳了 x¨m的不可测部分。
通过从控制对象模型公式(7)中减去CRM模型公式(8)得出误差方程:
e˙x=(ALC)exFa(ey)Fb(ey)+B1Ψ1Txm+B3(Ψ3+K)Txm+B3Λu
(16)
式中,A*=(A+B1 ψ1*T+B3 ψ3*T)。

4 Matlab/Simulink仿真验证

4.1 不确定性飞行器模型及参数

在第2节中,将三段式刚体假设和传统刚体飞机六自由度方程相结合,得到了飞翼布局柔性飞行器的纵向动力学模型。该模型的状态包括风速V、迎角α、俯仰角θ、俯仰速率q、二面角η。围绕模型的线性化得到具有不确定参数的柔性飞行器模型。测量值为飞行器的垂直加速度Az、二面角η和俯仰速率q。其他状态如α η˙是不可测量的,也不能用于控制。控制目标是使中心升降舵δe和外部副翼δa能够跟踪垂直加速度Az的指令,同时调节η
柔性飞行器模型的初始参数为V0=9.144 m/s,α0=0°,θ0=0°,q0=0°/s,线性化模型的数值在公式(15)中给出。其中,{A,B3,C}在η0=11°的线性化模型可以用下列公式近似:
Θp*= 0.6-4.5200.050.411.470.11.830-0.02-0.35-0.59
Λ*= 0.910.520.520.79
η≥11°时,柔性飞行器的俯仰模态不稳定,模型(15)包含固有频率ωn=2 deg/s和阻尼比ζ=0.7的二阶致动器动力学。ue是升降舵指令,ua是副翼指令。
x˙=Ax+B3uy=Cx
(17)
式中,AB3C分别表示为:
A= -0.2793.476-32.2-0.0150.5140.52500-2.57-6.4700-0.070-4.10401.0130.1930.10000-0.795-0.079000001000000000-54.0400.2551.84521.41005.991-6.363000000010000000.0020.04400.819-0.075-6.518000.195-0.03400000000100000000000010000000000-10-1.40000000000-10-1.4000001000000000-123.12000000-23.84-2.37600
B3= 000000001000000000000100T
C= 000100000000000000000010000000000001
式中,xyu分别表示为:
x=[Vαθqηη˙δeδaδ˙eδ˙aωηωAs]y=[qωηωAz]u=[ueua]
(18)
对于控制设计,使用非方阵飞行器的平方法推导出控制对象的配平控制参数。Ls3采用 a22=1, a21=2, a20=1,ε=30,Λmax=2,Ψmax=30确定。
本次仿真是在包含标称二阶致动器的飞翼布局柔性飞行器模型上,在充分考虑外部环境以及飞行器动力学耦合作用,致使飞行器出现参数不确定问题的基础上,对比使用LQR方法,产生基于观测器的线性控制器(没有加入自适应)和加入修正闭环参考模型的参数自适应的反馈控制器在不确定性控制方向上的结果,凸显加入闭环参考模型的自适应反馈控制器在柔性飞行器不确定性问题上优越的控制性能。

4.2 基线LQR控制器仿真

本节作为对照组,采用基线LQR控制器来控制柔性飞行器的二阶标称致动器模型(ωn=1 rad/sec和阻尼比ζ=0.7)。控制目标使中心升降舵δe和外部副翼δa能够跟踪垂直加速度Az的指令,同时调节η。

4.3 自适应输出反馈控制器仿真

本节作为实验组,采用设计的自适应输出反馈控制器来控制柔性飞行器的二阶标称致动器模型(ωn=1 rad/sec和阻尼比ζ=0.7)。控制目标使中心升降舵δe和外部副翼δa能够跟踪垂直加速度Az的指令,同时调节η。
图6~13所示,通过对比分析可知,本文设计的自适应输出反馈控制在针对具有参数不确定的柔性飞行器模型控制方面,具有更加优越的性能,相比于基线LQR控制器,自适应输出反馈控制器能够更好地控制二面角和垂直加速度跟踪指令,并且中心升降舵和外部副翼也能更好地跟踪垂直加速度的指令。相关的自适应输出反馈控制器参数如图14所示。
图6 基线LQR控制下的二面角

Full size|PPT slide

图7 基线 LQR控制器下的垂直加速度

Full size|PPT slide

图8 基线LQR控制器下的外副翼

Full size|PPT slide

图9 基线LQR控制器下的升降舵

Full size|PPT slide

图10 自适应输出反馈控制下的二面角

Full size|PPT slide

图11 自适应输出反馈控制下的垂直加速度

Full size|PPT slide

图12 自适应输出反馈控制下的外副翼

Full size|PPT slide

图13 自适应输出反馈控制下的升降舵

Full size|PPT slide

图14 自适应输出反馈控制器的参数轨迹

Full size|PPT slide

5 结束语

本文针对具有不确定参数的飞翼布局柔性飞行器模型,提出了一种新的自适应输出反馈控制方法。其中,设计了闭环参考模型,引入附加滤波器和可调参数,通过生成正实数传递函数矩阵,进而实现全局稳定和渐近跟踪,并在带有二阶致动器的飞翼布局柔性飞行器模型上进行仿真验证,结果显示其具有良好的控制性能。下一步的工作计划是将自适应输出反馈控制应用在其他布局柔性飞行器上,并将相关度扩展到更高的阶数,以及搭建柔性飞行器的自适应输出反馈控制物理实验平台,进行实物验证。

参考文献

[1]
孟新宇. 临近空间飞艇建模与控制关键技术研究[D]. 南京: 南京航空航天大学, 2017.
[2]
张健, 张德虎. 高空长航时太阳能无人机总体设计要点分析[J]. 航空学报, 2016, 37(S1):1-7.
高空长航时(HALE)太阳能无人机(UAV)的基本工作原理与常规动力飞机相比有显著的区别,体现在飞机总体设计方法上有其独到之处,方案设计中还需要对一些关键技术细节进行认真权衡。阐述了高空长航时太阳能无人机总体设计中遵循的重量不变和能量平衡原则的机理,同时从太阳能无人机巡航平飞功率需求、布局形式选取、飞行剖面优化和临近空间使用环境影响等方面出发,研究了太阳能无人机总体设计中的若干注意事项,主要结论可用于高空长航时太阳能无人机总体设计和方案优化。
[3]
张健, 王江三, 耿延升, 等. 高空长航时太阳能无人机的技术挑战[J]. 航空科学技术, 2020, 31(4):14-20.
[4]
杨铭超. 近空间飞行器起飞段抗风及模态切换鲁棒控制器设计[D]. 南京: 南京航空航天大学, 2017.
[5]
李军, 李占科, 宋笔锋. 联翼高空长航时无人机总体布局设计研究[J]. 飞行力学, 2009, 27(4):1-4.
[6]
吕蓉蓉. 大柔性飞行器阵风扰动影响分析及稳定控制研究[D]. 南京: 南京航空航天大学, 2018.
[7]
沈华勋, 徐亮, 陆宇平, 等. 飞翼式柔性飞机纵向动力学建模与稳定性分析[J]. 动力学与控制学报, 2016, 14(3):241-246.
[8]
Patil M J, Hodges D H. Flight Dynamics of Highly Flexible Flying Wings[J]. Journal of Aircraft, 2006, 43(6): 1790-1799.
[9]
Raghavan B, Patil M J. Flight Dynamics of High Aspect-ratio Flying Wings: Effect of Large Trim Deformation[J]. Journal of Aircraft, 2009, 46(5): 1808-1812.
[10]
Su W, Cesnik C E S. Dynamic Response of Highly Flexible Flying Wings[J]. AIAA Journal, 2011, 49(2): 324-339.
[11]
Doyle J, Stein G. Multivariable Feedback Design: Concepts for a Classical/Modern Synthesis[J]. IEEE transactions on Automatic Control, 1981, 26(1): 4-16.
[12]
Rynaski E. Flight Control Synthesis using Robust Output Observers[C]// Guidance and Control Conference, San Diego,CA,USA, 1982: 1575.
[13]
靳瑾, 张景瑞, 刘藻珍. 带有执行机构不确定性的航天器自适应姿态跟踪与参数辨识方法[C]// 第四届全国动力学与控制青年学者研讨会论文摘要集,辽宁,大连, 2010:102-103.
[14]
Morse A S. Parameterizations for Multivariable Adaptive Control[C]// 1981 20th IEEE Conference on Decision and Control Including the Symposium on Adaptive Processes, IEEE,San Diego, CA, USA, 1981: 970-972.
[15]
胡志坤, 孙岩, 姜斌, 等. 一种基于最优未知输入观测器的故障诊断方法[J]. 自动化学报, 2013, 39(8):1225-1230.
[16]
Gibson T, Annaswamy A, Lavretsky E. Improved Transient Response in Adaptive Control Using Projection Algorithms and Closed Loop Reference Models[C]// AIAA Guidance, Navigation, and Control Conference.Phoenix, USA, 2012: 4775.
PDF(2192 KB)

2966

Accesses

0

Citation

Detail

段落导航
相关文章

/