CN108228977A - 一种基于飞行状态参数的直升机振动特征折算方法 - Google Patents

Abstract

本发明涉及一种基于飞行状态参数的直升机振动特征折算方法，采用RBF神经网络的方法分别在标准飞行状态和目标飞行状态下建立直升机飞行状态参数及发动机状态参数与振动特征参数之间的关系，分别计算在当前飞行状态参数及发动机状态参数条件下和标准状态下的振动特征参数理论值，通过比较两者之间的差别得出相应的振动参数折算因子，并利用折算因子对实际测得的振动特征进行折算，从而达到将振动特征参数折算到标准工况下的目的。

Description

一种基于飞行状态参数的直升机振动特征折算方法

技术领域

本发明属于直升机状态监测与健康管理领域，具体涉及一种应用于变飞行状态条件下直升机机上振动特征参数的折算方法。

背景技术

对直升机的振动特性进行分析，可以认为直升机的振动主要来自于四个方面：直升机本身的振动，飞行状态变化引起的振动，故障引发的振动和噪声。根据直升机飞行状态的不同，直升机各相关振动传感器测得的振动信号及其相应的振动特征参数等也都会不同。在采用基于特征参数趋势分析的方法进行故障预测之时，若不消除飞行状态变化引起的振动对振动特征参数的影响，将会极大的影响预测的精度和准确度。

通常，直升机在飞行过程中状态多变，扭矩、速度等参数变化频繁，导致直升机的振动特性变化较大从而掩盖了因其本身健康状态性能退化而产生的振动特性的变化。因此，为提高趋势分析和故障预测处理结果的可信度，尽可能减小飞行状态变化导致的影响，需要开展基于飞行状态参数的振动特征参数折算方法研究。

发明内容

本发明的目的是：提供一种可根据直升机飞行状态参数进行振动特征参数折算的方法，采用了基于RBF神经网络的方法，可将不同飞行状态下的振动特征参数折算至标准状态下。

考虑到现有技术的上述问题，根据本发明公开的一个方面，本发明采用以下技术方案：

一种基于飞行状态参数的直升机振动特征折算方法，包括以下步骤：

步骤一：针对直升机机上监测的各种飞行状态参数和发动机状态参数，选取与其飞行状态及振动特征相关的参数，利用标准飞行状态下的样本数据，将样本数据中的飞行状态参数和发动机状态参数作为输入，振动特征参数作为输出，对RBF神经网络进行训练，构建标准飞行状态下的RBF模型；

步骤二：基于已构建的标准飞行状态下的RBF模型，根据标准飞行状态下样本数据中的飞行状态参数和发动机状态参数进行计算，得到振动特征参数的理论值，计算其均值，作为标准飞行状态下振动特征参数的基准值；

步骤三：选取目标飞行状态下的样本数据，将样本数据中的飞行状态参数和发动机状态参数作为输入，振动特征参数作为输出，对RBF神经网络进行训练，构建目标飞行状态下的RBF模型；

步骤四：将该目标飞行状态下的实时监测数据中的飞行状态参数和发动机状态参数作为已构建的目标飞行状态下RBF模型的输入，计算得到其对应的振动特征参数的理论值，与标准飞行状态下振动特征参数的基准值作差，得到当前飞行状态参数和发动机状态参数下的折算因子，将该时刻对应的振动特征参数实测值减去折算因子，得到该时刻的振动特征参数实测值折算至标准状态下的振动特征参数。

有益效果：

根据本发明的方法，可以将在不同飞行状态下测得的直升机振动特征参数快速折算至同一标准飞行状态下，从而可以保证在进行趋势分析和故障预测过程中其输入振动特征参数处于同一标准飞行状态工况下，提升故障预测的准确性。

附图说明

图1为RBF神经网络结构图

图2为聚类中心计算示意图

图3为振动特征折算流程图

具体实施方式

下面结合实施例对本发明作进一步地详细说明，但本发明的实施方式不限于此。

本发明的基于飞行状态参数的直升机振动特征折算方法，采用RBF神经网络的方法建立直升机飞行状态参数及发动机状态参数与振动特征参数之间的关系，分别计算在当前飞行状态参数及发动机状态参数条件下和标准状态下的振动特征参数理论值，通过比较两者之间的差别得出相应的振动参数折算因子，并利用折算因子对实际测得的振动特征进行折算，从而达到将振动特征参数折算到标准工况下的目的，为直升机HUMS系统的研制提供技术支撑。具体包括以下步骤：

步骤一：X₀是直升机标准飞行状态下与其飞行状态及振动特征参数相关的飞行状态参数和发动机状态参数样本数据，为n×m矩阵，其中m为相关的飞行状态参数和发动机参数的个数，n为样本量，即在标准飞行状态下在n个时刻测得了m个飞行状态参数和发动机状态参数的样本数据。Y₀是其对应的振动特征参数样本，为n维向量。将X₀作为输入参数，Y₀作为输出参数，对RBF神经网络进行训练。RBF神经网络的结构如图1所示，其中输入节点的个数为m，输出节点的个数为1，隐节点的个数为2m，ω_i(i＝1,2,…,2n)为隐层到输出层的权值，f_i(i＝1,2,…,2m)为隐节点的基函数，采用高斯函数其中，x是m维输入向量；c_i是第i个基函数的中心，与x具有相同维数的向量，σ_i决定第i个基函数围绕中心点的宽度，基函数中心c_i采用K均值算法计算得到，具体算法如下：

a)在X₀的n行行向量中，随机选取2m个行向量作为基函数中心c_i的初始值。

b)将样本X₀中所有行向量x^p(p＝1，2，…，n)按最近的聚类中心分组，即：如果意味着x^p到第j个聚类中心cj的距离最近，则将样本x^p划归到第j类。

c)将X₀中所有行向量重新归类后，计算各类的样本均值作为新的聚类中心，即：其中N_j为第j类的样本数。

d)重复步骤b)、c)直至所有聚类中心不再变化，从而得到最终的隐节点基函数中心c_i(i＝1,2,…,2m)。

半径σ_i决定了RBF单元接受域的大小，对网络的精度有极大的影响，计算方法可令每个基函数中心c_i与该类训练样本之间的平均距离作为相应的半径σ_i，即：其中N_i为第i类的样本数。权值的初始值可以取均值权值的调整采用梯度下降法，其迭代公式为：ω_i(t+1)＝ω_i(t)+η(y-u)f_i(x)，其中η为学习速率，y为网络的期望输出即振动特征参数的实际监测值，u为权值调整前网络的输出，f_i(x)为相应隐节点的输出。

根据样本数据X₀和Y₀计算完成隐层到输出层的权值ω_i(i＝1,2,…,2m)，隐节点的基函数中心c_i(i＝1,2,…,2m)和半径σ_i(i＝1,2,…,2m)后，即完成标准飞行状态下的RBF模型构建。

步骤二：基于根据步骤一已构建的标准飞行状态下的RBF模型，根据标准飞行状态下样本数据中的飞行状态参数和发动机状态参数X₀作为模型输入进行计算，得到振动特征参数的模型输出值Y₀′，计算其均值将作为标准飞行状态下振动特征参数的基准值；

步骤三：选取目标飞行状态下的样本数据，其中X₁是该目标飞行状态下与步骤一中的X₀相同类型飞行状态参数和发动机状态参数的样本数据，为s×m矩阵，其中m为相关的飞行状态参数和发动机参数的个数，s为样本量，即在该目标飞行状态下在s个时刻测得了m个飞行状态参数和发动机状态参数的样本数据。Y₁是其对应的与步骤一中同一类振动特征参数的样本数据，为s维向量。RBF神经网络的结构与步骤一中相同，如图1所示，将X₁作为输入参数，Y₁作为输出参数，对RBF神经网络进行训练。

其中输入节点的个数为m，输出节点的个数为1，隐节点的个数为2m，w_i(i＝1,2,…,2n)为隐层到输出层的权值，g_i(i＝1,2,…,2m)为隐节点的基函数，采用高斯函数其中，x是m维输入向量；c′_i是第i个基函数的中心，与x具有相同维数的向量，σ′_i决定第i个基函数围绕中心点的宽度，基函数中心c′_i采用K均值算法计算得到，具体算法如下：

a)在X₁的s行行向量中，随机选取2m个行向量作为基函数中心c′_i的初始值。

b)将样本X₁中所有行向量x₁ ^p(p＝1，2，…，s)按最近的聚类中心分组，即：如果意味着x₁ ^p到第j个聚类中心c′_j的距离最近，则将样本x₁ ^p划归到第j类。

c)将X₁中所有行向量重新归类后，计算各类的样本均值作为新的聚类中心，即：其中N_j为第j类的样本数，如图2所示。

d)重复步骤b)、c)直至所有聚类中心不再变化，从而得到最终的隐节点基函数中心c′_i(i＝1,2,…,2m)。

半径σ′_i决定了RBF单元接受域的大小，对网络的精度有极大的影响，计算方法可令每个基函数中心c′_i与该类训练样本之间的平均距离作为相应的半径σ′_i，即：其中N_i为第i类的样本数。权值的初始值可以取均值权值的调整采用梯度下降法，其迭代公式为：w_i(t+1)＝w_i(t)+η₁(y₁-u₁)g_i(x)，其中η₁为学习速率，y₁为网络的期望输出即振动特征参数的实际监测值，u₁为权值调整前网络的输出，g_i(x)为相应隐节点的输出。根据样本数据X₁和Y₁计算完成隐层到输出层的权值w_i(i＝1,2,…,2m)，隐节点的基函数中心c′_i(i＝1,2,…,2m)和半径σ′_i(i＝1,2,…,2m)后，即完成该目标飞行状态下的RBF模型构建。

步骤四：在与步骤三相同的目标飞行状态下的某一时刻，测得与步骤三中的X₁相同类型飞行状态参数和发动机状态参数的实测数据为x₂，为一m维向量，其对应与步骤三中相同类型的振动特征参数的实测数据为y₂，将x₂作为步骤三中已构建的该目标飞行状态下的RBF模型的输入，计算得到其对应的振动特征参数的理论值u₂，与步骤二中得到的标准飞行状态下振动特征参数的基准值作差，得到当前飞行状态参数和发动机状态参数下的折算因子将该时刻对应的振动特征参数实测值y₂减去折算因子Δy，得到的值即为该时刻的振动特征参数实测值折算至标准状态下的振动特征参数y′₂，即y′₂＝y₂-Δy，此时已完成振动特征参数的折算，可以认为y′₂是当前时刻在标准飞行状态下的振动特征参数值如图3所示。按照上述方法将可将所有典型飞行状态下的振动特征参数值折算至标准飞行状态下，从而在趋势分析和故障预测中保证了直升机所有振动监测数据都在同一工况之下。

Claims (5)

1.一种基于飞行状态参数的直升机振动特征折算方法，其特征在于，包括如下步骤：

1)针对直升机机上监测的各种飞行状态参数和发动机状态参数，选取与其飞行状态及振动特征相关的参数，利用标准飞行状态下的样本数据，构建标准飞行状态下的RBF模型；

2)基于步骤1)构建的标准飞行状态下的RBF模型，根据标准飞行状态下样本数据中的飞行状态参数和发动机状态参数进行计算，得到振动特征参数的基准值；

3)采用目标飞行状态下的样本数据，构建飞行状态参数及发动机状态参数与振动特征参数在目标飞行状态下的RBF关系模型；

4)将目标飞行状态下的实时监测数据中的飞行状态参数和发动机状态参数作为步骤3)中已构建的目标飞行状态下RBF模型的输入，基于模型输出和步骤2)中的标准飞行状态下振动特征参数的基准值将该时刻对应的振动特征参数实测值折算至标准飞行状态下。

2.根据权利要求1所述的一种基于飞行状态参数的直升机振动特征折算方法，其特征在于，步骤1)具体为：采用标准飞行状态下的样本数据，构建飞行状态参数及发动机状态参数与振动特征参数在标准飞行状态下的RBF关系模型。

3.根据权利要求2所述的一种基于飞行状态参数的直升机振动特征折算方法，其特征在于，步骤2)具体为：根据标准飞行状态下样本数据中的飞行状态参数和发动机状态参数进行计算，得到振动特征参数的理论值，计算其均值，作为标准飞行状态下振动特征参数的基准值。

4.根据权利要求3所述的一种基于飞行状态参数的直升机振动特征折算方法，其特征在于，步骤3)具体为：选取目标飞行状态下的样本数据，将目标飞行状态下的样本数据中的飞行状态参数和发动机状态参数作为输入，振动特征参数作为输出，对RBF神经网络进行训练，构建目标飞行状态下的RBF模型。

5.根据权利要求4所述的一种基于飞行状态参数的直升机振动特征折算方法，其特征在于步骤4)具体为：将该目标飞行状态下的实时监测数据中的飞行状态参数和发动机状态参数作为已构建的目标飞行状态下RBF模型的输入，将模型输出与标准飞行状态下振动特征参数的基准值作差，得到当前飞行状态参数和发动机状态参数下的折算因子，将该时刻对应的振动特征参数实测值减去折算因子，得到折算后的振动特征参数。