一种飞行器弹性新型控制方法研究

2018-11-21郑勇斌魏明英

导航定位与授时 2018年6期

郑勇斌，魏明英

(北京电子工程总体研究所，北京 100854)

0 引言

为了扩大飞行器的杀伤范围和提高执行多种飞行任务的能力，新一代飞行器普遍朝着高速、轻质化、飞得更远、机动能力强等方面发展，导致飞行器出现了下列特征：1)装载更多燃料；2)存在多级飞行器；3)结构上采用复合材料；4)长细比增大；5)飞行速度更快；6)静不稳定增加。

由此产生了新的问题：

1)长细比增大，多级飞行器质量增加，而且飞行速度更快，飞行包络线更大，气动加热引起结构刚度显著下降，导致飞行器弹性频率越来越低；

2)多级飞行器连接方式易带来明显个体差异，很难保证燃料消耗规律预示精度，而且飞行过程中结构复合材料存在各向异性的特征，大幅增加了结构弹性变化的不确定性，导致了飞行器弹性频率变化与理论预示存在很大差异；

3)静不稳定增加、飞行器大机动能力要求和稳定裕度之间矛盾无法协调，严重影响了飞行器的总体设计。

这些新问题给飞行器弹性抑制和总体设计带来了严峻挑战，稳定裕度不足和弹性频率预示不准极易引起系统出现弹性振动现象，影响姿态控制精度要求，严重时会导致飞行试验失败，而且最重要的一点是严重影响了飞行器总体作战能力的提升。

目前飞行器传统弹性抑制控制方法是根据飞行器推进剂预示消耗规律理论计算弹性频率变化规律，设计时变宽频凹陷滤波器抑制飞行器弹性[1]。通常串连几个凹陷滤波器，兼容理论预示和飞行过程中的弹性频率差异，一般差异要求在2Hz范围内。这种抑制方法是以牺牲系统稳定性和快速性为代价，无法满足新一代飞行器大机动和高精度的作战需求，使得高精度快响应的飞行器弹性控制方法成为一个具有迫切需求而又充满挑战的难题。

国内外其他的弹性控制研究方法主要集中在利用带通滤波器辨识频率、弹性弹体建模的模型辨识、自适应滤波器辨识频率[2-4]、利用傅立叶变换计算弹性频率等方面，在线辨识出弹性频率，然后设计滤波器加以消除，从而可减小为了兼容弹性频率散布而带来的系统相位滞后。文献[5]介绍了美国土星五号通过在线辨识弹性频率，设计自适应弹性滤波器实现了弹性振动的抑制；文献[6]介绍了利用波音747飞机气动弹性模型，证明设计自适应滤波器可以有效提高飞行控制系统的带宽和气动伺服弹性系统的稳定裕度。文献[7]和文献[8]采用H∞设计方法处理弹性振动模态的参数以及结构不确定性，设计目标包括幅相裕度、弹性模态抑制以及一些时域指标。类似的设计方法也出现在处理超高声速飞行器弹性模态不确定性的文献[9-11]。文献[12-14]采用DFT方法找到弹性频率，然后在线调整滤波器参数实现弹性抑制。文献[15]采用插值FFT方法辨识出弹性振动信号，从而加以消除。这些方法建立在高精度弹性频率辨识上，一旦弹性频率辨识不准，则无法在工程上应用，而且目前研究结果是基于一定假设基础上得到的，因此在辨识精度、辨识速度、辨识灵敏度和时延等方面研究还不够深入，工程应用性不强。

本文在工程上创新性地首次提出一种飞行器弹性新型控制方法，即飞行器刚体信息高精度快速提取方法，实现了飞行器弹性抑制和飞行控制。此方法不同于参考文献中的其他使用方法，不依赖于弹性频率高精度辨识，工程应用性强，解决了传统弹性滤波方法带来的控制系统稳定裕度严重不足、机动能力弱的难题，对静不稳定大的飞行器具有更好的适应性，而且极大地提升了飞行器作战边界和控制精度。

1 传统弹性抑制控制方法局限

为了充分抑制飞行器弹性频率，传统弹性抑制控制方法往往通过设计宽频凹陷滤波器加以抑制，凹陷滤波器形式如下：

对应滤波器幅相特性如图1所示。

其中，凹陷滤波器1的中心频率为22Hz，在频率10rad/s处相位滞后6.5°，在频率20rad/s处相位滞后14°；凹陷滤波器2的中心频率为11Hz，在频率10rad/s处相位滞后17°，在频率20rad/s处相位滞后36°。弹性频率越低，系统相位滞后越严重，从而使得飞行器控制幅值裕度和相位裕度越难以协调保证，严重制约了飞行器总体设计和控制系统设计。

2 扩维滤波控制方法

2.1 模型建立

其中,az代表角加速度。

弹性变形引起的角速度可以表示为2个频率下的正弦信号。

(1)

对式(1)取时间的二阶导数，得

(2)

将式(2)写成状态方程的形式为

其中：

设x=ωz为系统状态(刚体角速度)，并设增广状态向量为

则增广系统

(3)

测量方程写作

z(t)=HX(t)+ζ(t)

(4)

其中，z(t)代表陀螺测量值，H=[1 1 0

1 0]，此处考虑了弹性变形对角速度信号量测的影响。ζ(t)是量测噪声，为零均值高斯白噪声过程。

将方程式(3)、式(4)按采样间隔时间T离散化得到

X[k+1]=ΦX[k]+Bdu[k]+Γζ[k]

(5)

z[k]=HX[k]+ζ[k]

(6)

对于增广系统式(5)和式(6)，Kalman滤波算法为

K(k+1)=P(k+1/k)HT(k+1)[H(k+1)·

P(k+1/k)HT(k+1)+R(k+1)]-1

P(k+1)=[I-K(k+1)H(k+1)]P(k+1/k)

初始状态X(0)的协方差矩阵为

噪声的统计特性为

则新的动态噪声的协方差阵为

把提取出的刚体角速度信息引入过载控制，至此就建立了扩维滤波控制方法。

2.2 可观性分析

可观性是扩维滤波渐近稳定的基本条件。

Rank(O)=5

所以系统是可观的。

2.3 可控性分析

Rank(O1)=1

所以系统不完全可控，不能任意配置极点。

3 数学仿真与分析

采用扩维滤波控制方法搭建过载控制仿真模型，与传统弹性滤波控制方法进行比较分析。

3.1 弹性抑制效果分析

假设飞行器一阶弹性频率为11Hz，按照传统弹性滤波器抑制方法设计了宽频凹陷滤波器，并比较了弹性频率预示精度2Hz和4Hz差异经弹性滤波器后输出信息，如图2～图4所示。

从上述曲线可以看出，传统弹性滤波控制方法在频率散布下弹性抑制能力较弱，扩维滤波控制方法对弹性频率散布不敏感；而且弹性频率越低，传统滤波器带来的相位滞后越严重，严重影响了过载跟踪响应，系统出现振荡，而扩维滤波控制方法下过载跟踪品质良好。

另外一点需要说明的是，采用传统弹性滤波控制方法控制参数设计具有很大的局限性，难以协调刚体裕度和弹性裕度，因此系统往往响应较慢，且系统容易因稳定裕度不足出现振荡发散；而采用扩维滤波控制方法控制参数设计具有更大的灵活性，系统响应快且平稳，尤其适合飞行器弹性频率预示不准或者弹性频率发生跳变的情况。

3.2 系统性能分析

工程应用时需分析陀螺输出噪声模型、角加速度计噪声模型、系统方程参数变化、输出方程参数变化、误差协方差初值、噪声协方差初值、增益矩阵初值、延时等对弹性抑制效果和系统性能的影响。

本文暂以陀螺输出噪声和角加速度计噪声开展性能分析。

1)陀螺输出噪声影响分析

考虑陀螺输出包括弹性分量1(°)/s以及动态环境下角运动噪声方差0.5(°)/s，得到下列仿真曲线，如图5～图7所示。从曲线可以看出，陀螺输出经传统弹性滤波器后，依然含有陀螺角运动噪声信息，表现形式为低频振荡，经过控制通道计算导致执行机构中也含有低频振荡信息，严重时会导致系统发散；而扩维滤波控制方法能够很好地提取刚体信息，对陀螺角运动噪声不敏感。

2)角加速度计输出噪声影响分析

考虑角加速度计噪声方差30(°)/s2的新型滤波方法的效果，如图8～图10所示。从曲线可以看出，角加速度计噪声对过载跟踪响应影响较小。图10表明扩维滤波控制方法实现了刚体信息的快速提取，且滞后小于15ms，有利于提高系统响应快速性。