作为随控布局飞机设计中的一项关键技术,颤振主动抑制方法一直受到众多研究者的强烈关注。随着控制理论的发展,颤振主动抑制的控制律设计方法也逐渐从最优控制过渡到鲁棒控制。鲁棒控制通过描绘对象的不确定性,并在不确定性允许的摄动范围内综合设计控制律以使系统保持稳定性和性能鲁棒,其中的结构奇异值方法则是在H∞理论上发展起来的一种重要技术,且已经成为颤振主动抑制的一个有效工具[1, 2, 3]。因飞行器的颤振边界依赖于动压和马赫数等飞行参数,一种应用于颤振主动控制律设计的简单方法便是在μ框架中引入对飞行参数的扰动并确定使系统不稳定的最小参数变化。为进行气动伺服弹性系统的稳定性分析,Lind等[4]建立了包含动压不确定性的线性分式变换模型。Moulin[5]在气动伺服弹性系统的状态空间模型基础上,引入空速和空气密度扰动并建立线性分式变换模型,以此为基础设计了颤振抑制的μ控制器。为便于设计鲁棒控制器,Qian等[6]在Moulin的方法的基础上通过模型变换有效降低了不确定块的维数,并设计了颤振抑制效果更好的鲁棒控制器。但前述的模型均以刚性机翼或低速机翼为研究对象,很少以亚声速范围内的机翼为研究对象。其中一个重要原因是在亚声速范围内马赫数对广义非定常气动力系数的影响无法解析表达,难于构造相应的线性分式变换模型。除鲁棒控制方法外,变增益控制[7]、非线性控制[8]、自适应控制[9]等方法也在逐渐获得深入研究并应用于颤振主动抑制。
Blue等[10]对一个二元翼段建立了随马赫数和动压变化的近似线性参数变化模型,在较大的飞行范围内都取得了足够的精度,但没有考虑广义气动力系数的变化特征和弹性模态的影响。本文克服了上述缺点并以线性参数变化模型为基础构造了包含动压和马赫数参数不确定性的线性分式变换模型,随后结合鲁棒控制理论设计了用于颤振抑制和阵风减缓的μ控制器,并对结果进行了时域仿真验证。
1 气动伺服弹性系统模型开环气动伺服弹性系统模型是控制器设计的出发点。在构造系统的开环运动方程之前,需要先将频域的广义气动力系数从频域有理延拓为拉氏变量s的有理函数,这可以通过最小状态法[11]来实现,其拟合公式为
式中, =sb/V,s为拉氏变量,b为机翼参考半弦长,V为来流速度,Q0、Q1、Q2、D和E均为实系数矩阵,R为气动滞后根组成的对角矩阵。Qξξ()、Qξδ()和Qg()分别对应于升力面模态、控制面偏转和阵风激励的广义非定常气动力系数矩阵。注意在拟合过程中可约束Q2中阵风激励的对应列为零以消除阵风速度的二阶导数项。
基于控制面和阵风激励的弹性体状态空间方程为
可简写为
式中
式中,ξ为机翼弹性模态对应的n维广义坐标,xa为na维气动力状态扩充项,Mξξ广义质量矩阵,Kξξ为广义刚度矩阵,Cξξ为广义阻尼矩阵,Mξδ为结构模态和控制面的耦合质量阵,δ为控制面偏转模态坐标,wg为垂直阵风速度。
以翼尖加速度为输出时,输出方程为
式中Φ为翼尖节点的模态向量。
舵机的传递函数可表示为
式中, Aac、Bac和Cac分别表示系统矩阵、输入矩阵和输出矩阵,xac为舵机状态向量,uc为舵机输入信号。
连续阵风的状态空间方程形式为
式中各矩阵可由Dryden阵风模型[12]或Von Karman阵风模型[12]得到。其中Ag、Bg、Cg和Dg分别表示系统矩阵、输入矩阵、输出矩阵和直接转换矩阵,xg为阵风状态向量,η为白噪声输入信号。
引入舵机和阵风环节,可得到弹性机翼的开环状态空间方程
可简写为
式中
2 随马赫数和动压变化的线性参数变化模型由(2)式、(3)式和(8)式可知,矩阵A和B中部分元素受马赫数Ma、空速V和空气密度ρ的影响,其中马赫数通过广义气动力系数矩阵间接影响系数矩阵的取值。但在标准大气情况下,Ma、V和ρ三者的取值可通过2个独立变量Ma和q得到,因此(9)式可描述为
式中
式中,Δq和ΔMa分别表示马赫数和动压的最大变化量,δ1和δ2分别表示动压和马赫数的归一化不确定性,取值范围为[-1,1]。
(11)式中飞行参数对系数矩阵的影响可通过(13)式来近似[10]以构造线性参数变化模型
在给定Ma和q的取值范围后,(13)式右端的各未知系数矩阵可通过最小二乘法求解得到。
仔细观察(2)式和(8)式可发现,仅和a的对应行(即第n+1行到第2n+na行)的各矩阵系数受Ma和q影响,故可将(9)式的第n+1行到第2n+na行重写为
3 线性分式变换模型为便于鲁棒控制器的设计,需要将上节得到的线性参数变化模型转换为线性分式变换模型。首先将(14)式重组如下
再定义如下的输入和输出变量z和v
显然有
式中
现将(16)式代入(15)式可得
再将(19)式代入(9)式并结合(16)式可得
(20)式中各矩阵可简单推得,此处恕不详述。
将(20)式中对应行代入(5)式可得
此处省略(21)式右端各矩阵的详细描述。
组合(20)式与(21)式可得
令P(s)为S的拉氏变换,则有
结合(23)式和(17)式可得
式中,Fl(P,Δ)表示P(s)与Δ构成的下线性分式变换,其表达式见下式,可用图 1表示。
为便于控制器设计,需要将不确定块Δ归一化以使其所有元素绝对值不大于1,这可在(12)式的基础上通过线性分式变换将Fl(P,Δ)转换为Fl(P1,Δ1),其中为Δ1为归一化的不确定块。
至此已构造出能够计及飞行参数变化的统一模型,下面将通过算例构造模型并设计用于颤振主动抑制和阵风减缓的鲁棒控制器。
4 算 例 4.1 开环颤振分析算例机翼的有限元模型如图 2所示,机翼根部固支。采用MSC/Nastran软件计算得到前9阶结构模态,其中前4阶模态特征如表 1所示。采用ZAERO软件计算各马赫数下的非定常气动力系数,图 3给出了机翼气动力模型及控制面的位置。机翼在马赫数0.5~0.7范围内由P-K法计算得到的颤振边界如表 2所示。
Ma | 颤振速度/(m·s -1) | 颤振频率/Hz | 颤振动压/kPa |
0.5 | 206.92 | 17.07 | 139.12 |
0.55 | 223.48 | 15.24 | 138.56 |
0.6 | 239.28 | 15.12 | 135.42 |
0.65 | 254.57 | 14.53 | 131.36 |
0.7 | 269.32 | 13.99 | 126.40 |
在构造随动压和马赫数变化的线性参数变化模型时,需要在给定的飞行范围内选取部分飞行状态并用最小二乘法求解(13)式右端的各系数矩阵。本文选取的飞行状态点如表 3所示。从表 3可以看出,开环颤振点附近的动压增量被适当缩小以改善精度。在构造各飞行状态点的状态空间模型时,对非定常气动力系数的有理拟合均采用4个滞后根,采用Dryden连续阵风模型,阵风强度为1.5 m/s,阵风尺度为760 m,输出为翼尖前缘点的过载,并将其作为反馈控制器的输入信号。舵机传递函数见(26)式
Ma | 动压/kPa | ||||
0.5 | 100 | 125 | 134 | 144 | 165 |
0.55 | 100 | 125 | 133 | 143 | 165 |
0.6 | 100 | 125 | 130 | 140 | 165 |
0.65 | 100 | 126 | 136 | 141 | 165 |
0.7 | 100 | 120 | 131 | 140 | 165 |
为验证拟合结果的精度,本文比较了线性参数变化模型和原状态空间模型在各马赫数下的颤振点特征。结果表明颤振动压误差均小于1%,这说明线性参数变化模型的精度是足够的。图 4比较了部分飞行条件下,原状态空间模型与线性参数变化模型中翼尖前缘的阵风频响特性曲线。从图中可以看出两者几乎一致,这从另一个侧面验证了模型的精度。
4.2 μ控制器设计对于当前模型,反馈控制的主要目标是抑制开环气动弹性系统的颤振及减缓阵风响应。在设计鲁棒控制律之前,需要构造完整的系统模型。图 5给出了闭环鲁棒控制模型的结构,其中K即为待设计的鲁棒控制器。图中P1即为名义系统模型,Δ1为代表动压和马赫数不确定性的且维数为39×39的对角不确定块,Wy为翼尖过载的输出加权函数,表达式为
因连续阵风频谱主要集中于低频范围,因此在低频范围施加较高的权重能够充分发挥控制器阵风减缓的能力。本文将低频范围取为0~70 rad/s,计算可得突风到翼尖过载的传递函数在稳定飞行范围内的峰值为0.04,则为了将对应飞行状态下的翼尖过载降低至少50%,需将Wy的系数取为2/0.04。
图 5中wn代表传感器的量测白噪声信号,添加输出噪声有2个方面的作用:①提高鲁棒控制器的抗噪声能力;②完整化模型结构以避免控制器设计过程中的可能出现的数值问题[13]。ε为对量测噪声的加权小量,取值为0.001。
Wu为舵面控制信号的输出加权函数。为抑制翼尖的高频振荡部分和量测噪声中的高频信号引起的舵面偏转信号,需要对舵面控制信号的高频部分施加较高的权重,本文取为
式中,Wu的系数为突风白噪声输入到舵面控制信号闭环传递函数最大允许增益的倒数。
图 5中的模型可通过Matlab鲁棒控制工具箱[14]中的sysic命令转换成适于鲁棒控制器设计的结构。取Ma0=0.5,ΔMa=0.1,通过多次调整,在q0=115 kPa且Δq=15 kPa时控制器效果达到最佳。最终的鲁棒控制器通过平衡截断法[14]被降至11阶,且闭环系统的最大μ值为0.977,达到鲁棒稳定性和鲁棒性能要求。下节将通过闭环颤振及阵风响应分析进一步检验μ控制器的性能。
4.3 闭环颤振及阵风响应分析图 6给出了机翼在马赫数在0.5~0.7范围内开环系统和闭环系统的颤振边界,从图中可以看出,鲁棒控制器将颤振动压在给定范围内平均提高了约10%。
为验证μ控制器的阵风减缓效果,令飞行马赫数在20 s内从0.7线性变化到0.5,同时动压从100 kPa线性变化到130 kPa,易知在此过程中开环和闭环系统都是稳定的。图 7给出了开环和闭环系统的翼尖过载时域响应。从图 7中可以看出控制器有明显的阵风减缓效果,开环系统的翼尖过载均方根为1.07g,而闭环系统为0.52g,减小了51.4%。
图 8给出了闭环系统的舵面转角和速度响应,从图中可以看出,控制面的转角和速度均在合理范围内,这也可以为舵机的选择提供参考。
5 结 论提出了一种亚声速下弹性机翼的颤振主动抑制与阵风减缓方法。因机翼的非定常气动力同时受马赫数和动压的影响,本文以线性参数变化模型为基础构造了含马赫数和动压不确定性的线性分式变换模型。最终设计出的鲁棒控制器,不仅将马赫数0.5~0.7范围内的颤振动压平均提高10%,还有效减缓了机翼的阵风响应。
[1] | Borglund D. Robust Aeroelastic Stability Analysis Considering Frequency-Domain Aerodynamic Uncertainty[J]. Journal of Aircraft, 2003, 40(1):189-193 |
Click to display the text | |
[2] | Karpel M, Moulin B, Idan M. Robust Aeroservoelastic Design with Structural Variations and Modeling Uncertainties[J]. Journal of Aircraft, 2003, 40(5):946-954 |
Click to display the text | |
[3] | Borglund D. The μ-k Method for Robust Flutter Solutions[J]. Journal of Aircraft, 2004, 41(5):1209-1216 |
Click to display the text | |
[4] | Lind R, Brenner M. Aeroelastic and Aeroservoelastic Models, Robust Aeroservoelastic Stability Analysis:Flight Test Applications[M]. London:Springer, 1999:55-63 |
Click to display the text | |
[5] | Moulin B. Robust Controller Design for Active Flutter Suppression[R]. AIAA-2004-5115 |
Click to display the text | |
[6] | Qian W M, Huang R, Hu H Y, et al. New Method of Modeling Uncertainty for Robust Flutter Suppression[J]. Journal of Aircraft, 2013, 50(3):994-998 |
Click to display the text | |
[7] | Barker J M, Blalas G J, Blue P A. Gain-Scheduled Linear Fractional Control for Active Flutter Suppression[J]. Journal of Guidance, Control, and Dynamics, 1999, 22(4):507-512 |
Click to display the text | |
[8] | Ronch A D, Tantaroudas N D, Jiffri S, et al. A Nonlinear Controller for Flutter Suppression:from Simulation to Wind Tunnel Testing[R]. AIAA-2014-0345 |
Click to display the text | |
[9] | Cassaro M, Battipede M, Marzocca P, et al. Comparison of Adaptive Control Architectures for Flutter Suppression[J]. Journal of Guidance, Control, and Dynamics, 2015, 38(2):346-355 |
Click to display the text | |
[10] | Blue P, Balas G J. Linear Parameter-Varying Control for Active Flutter Suppression[R], AIAA-1997-3640 |
Click to display the text | |
[11] | Azoulay D, Karpel M. Characterization of Method for Computation of Aeroservoelastic Response to Gust Excitation[R]. AIAA-2006-1938 |
Click to display the text | |
[12] | Hoblit F M, Gust Loads on Aircraft:Concepts and Applications[M]. Washington D C, American Institute of Aeronautics and Astronautics, Inc, 1988 |
Click to display the text | |
[13] | Doyle J C, Glover K, Khargonekar P, et al. State-Space Solutions to Standard H2 and H∞ Control Problems[C]//IEEE Trans on Automatic Control, 1989:831-847 |
Click to display the text | |
[14] | Moore B. Principal Component Analysis in Linear Systems:Controllability, Observability, and Model Reduction[C]//IEEE Trans on Automatic Control, 1981:17-31 |
Click to display the text |