CN110457731A - 一种应用于直升机动部件振动趋势分析的方法 - Google Patents

Abstract

本发明公开了一种应用于直升机动部件振动趋势分析的方法。本发明通过对各飞行状态下的振动特征参数建立统计分布，通过计算时间历程中分布之间得到的偏离值来建立退化趋势的状态指标，偏离值在[0,1]区间上，越靠近1表示部件越接近退化失效。本发明实现了对直升机动部件的性能退化趋势和未来的退化状态的准确分析和评估，从而有效、及时的开展维修工作，防止事故的发生。

Description

一种应用于直升机动部件振动趋势分析的方法

技术领域

本发明涉及趋势分析及故障预测领域，具体是一种应用于直升机动部件振动趋势分析的方法。

背景技术

直升机动部件是直升机上的关键系统，是由大量相互联系但工作过程又彼此不同的部件组成，结构复杂，长期处于高速、交变重载的恶劣环境下，容易产生故障和损坏，由于其无冗余设计，一旦发生故障，将导致灾难性事故的发生，直接影响飞行员的生命安全。

直升机动部件发生失效故障前会经历逐步的退化过程，随着直升机运行时间的积累，动部件处于不断的性能退化过程中直至失效，为了尽可能减少动部件故障导致的经济损失和重大事故，需要在其运行过程中建立起能够指示动部件各组成部件性能退化的状态指标，对直升机动部件的性能退化趋势和未来的退化状态进行准确分析和评估，就可有效、及时地开展维修工作，防止事故的发生。

发明内容

本发明的目的在于提供一种基于直升机动部件振动趋势分析的方法。

实现本发明目的的技术解决方案为：一种应用于直升机动部件振动趋势分析的方法，包括如下步骤：

1)基于提取的直升机飞行振动特征参数和对应的飞行状态参数，采用关联分析方法，选取与飞行振动特征参数相关的飞行状态参数，并建立对应的飞行状态参数集Φ；

2)根据选取的飞行状态参数，基于直升机初始飞行架次数据建立飞行振动特征参数的基准分布数据库；

3)计算当前飞行架次下的直升机飞行振动特征参数的统计分布，同时根据当前飞行架次飞行状态参数集Φ中飞行状态参数的分布选取基准分布数据建立当前架次的基准分布，计算当前架次飞行振动特征参数统计分布与基准分布的重叠面积；

4)建立对数归一化函数，设定失效状态下直升机飞行振动特征参数的统计分布与基准分布的重叠面积为S，分布偏离值为DI，对S和DI赋值计算确定对数归一化函数参数的取值；

5)基于上述对数归一化函数得到飞行振动特征参数分布偏离值，作为趋势分析的指标。

6)步骤1)中采用对飞行状态参数和飞行振动特征参数进行关联分析，剔除非相关的飞行状态参数；具体方法为：

1.1，利用标准化后的样本数据，计算飞行振动特征参数变化Y与每个飞行状态参数的变化X_i之间的相关性系数，并根据其大小选择飞行状态参数，相关系数在0-1之间，0表示不相关，1表示相关，这里对选取飞行状态参数的相关性系数阈值不做要求，根据实际情况确定，若要少选飞行状态参数高阈值可设高，反之可设低；相关性系数计算公式为其中为飞行状态参数的方差，σ_Y为飞行振动特征参数的方差；x′_ij为第i类参数的第j个样本标准化后数据值，y′_j为飞行振动特征参数的第j个样本标准化后数据值，μ_i′为第i类参数的标准化样本的均值，μ_y′为参数的标准化样本的均值，n为样本量；

1.2，设飞行振动特征参数变化Y为随机变量，经相关性分析后选取的飞行状态参数变化为自变量X_i，且飞行状态参数与飞行振动特征参数的关系模型满足多变量线性回归模型，采用基于回归分析的方法建立飞行状态参数与飞行振动特征参数的关系模型b为该回归模型函数的相关系数，设显著性水平为α，计算各飞行状态参数的P值，若P值小于或等于显著性水平α，则以对应的飞行状态参数建立飞行状态参数集Φ

步骤3)中，计算当前架次飞行振动特征参数统计分布与基准分布的重叠面积，具体方法为：

3.1针对当前飞行架次采用核密度估计的方法对其分布进行估计；

设X₁，X₂，…，X_n是取自一元连续总体的样本，在任意点x处的统计分布为n为样本点数，其中选择高斯函数作为核函数，u为核函数自变量，窗宽h的估计值通过公式计算，其中σ为样本标准差；

3.2根据当前飞行架次中飞行状态参数集Φ中参数变化范围，选取基准分布数据建立基准状态的统计分布g(x)；在基准分布数据中已确立了飞行状态参数和飞行振动特征参数的对应关系，仅需根据当前飞行中的飞行参数到基准分布数据中去选择对应的飞行振动特征参数；

3.3将基准状态分布的统计分布g(x)和当前架次分布的统计分布f(x)在同一坐标轴上进行表示，将f(x)和g(x)的重叠面积即重叠度作为用于评价当前飞行架次状态变化趋势的评价指标；针对每一特征参数，f(x)和g(x)重叠度S的计算公式为：S即为用于评价当前飞行架次状态变化趋势的评价指标，其值在(0，1)区间范围内呈递减的趋势。

步骤4)中建立对数归一化函数，通过设定失效所对应的分布重叠度值和分布偏离值，确定对数归一化函数参数，具体为：

建立对数归一化函数：其中S为计算得出的分布重叠度值，DI＝f(S)为分布偏离值，a为归一化参数，计算公式为：设定动部件失效对应的分布重叠度值为10^-4，对应的分布偏离值为0.99。

步骤5)中，采用对数归一化函数将重叠度映射为分布偏离值DI。

本发明与现有技术相比，其显著优点：

(1)根据直升机的飞行特点以每个飞行架次之间的状态变化建立趋势曲线，可在分析样本中包含多状态的信息并统一化各个分析样本的外部影响因素来提升趋势分析的正确性；

(2)根据当前架次的飞行状态自适应选择基准数据建立基准统计分布，减小了飞行状态变化对趋势分析结果的影响；

(3)建立了对数归一化函数，解决了状态退化过程中趋势曲线上重叠面积数值变化不明显的问题。

附图说明

图1为振动特征与飞行参数的关系。

图2为各飞行状态下振动特征的分布。

图3为时间历程上分布的重叠度。

图4为重叠度的归一化。

具体实施方式

本发明是一种可定量评价动部件性能退化状态的方法，考虑了直升机飞行过程中飞行状态的变化对振动的影响，采用了基于时间历程上各振动特征参数分布的偏离度来定量描述动部件的性能退化过程及性能退化状态。

本发明通过对各飞行状态下的振动特征参数建立统计分布，通过计算时间历程中分布之间得到的偏离值来建立退化趋势的状态指标，偏离值在[0,1]区间上，越靠近1表示部件越接近退化失效。

下面结合附图和实施例对本发明作进一步说明。

步骤一：首先利用标准化后的样本数据(直升机飞行中记录的飞行状态参数和提取的飞行振动特征参数数据，如图1所示)，计算每个飞参(旋翼转速、发动机扭矩、大气压力、飞行高度、指示空速、机外大气温度、偏航角、俯仰角、横滚角、三方向的加速度)的变化x与振动特征参数变化y之间的相关性系数，计算公式为其中x′_ij为第i类参数的第j个样本标准化后数据值；y′_j为振动特征参数的第j个样本标准化后数据值；μ_i′为第i 类参数的标准化样本的均值；μ_y′为参数的标准化样本的均值；n—样本量。其次假定振动特征参数Y为随机变量，经相关性后的飞行状态参数为自变量，且飞行状态参数与振动特征参数的关系模型满足多变量线性回归模型，采用基于回归分析的方法建立飞行状态参数与振动特征参数的关系模型经过计算得到该回归模型函数的相关系数，设显著性水平为α＝0.05，计算各因素的相关系数为作为飞行参数的选择依据，集合每个振动特征参数对应的飞行参数分布统计的基准。

步骤二：针对多数情况下振动特征分布形式未知的情况采用核密度估计的方法对振动特征分布进行估计。设X₁,X₂,…,X_n是取自一元连续总体的样本，在任意点x处的统计分布为其中选择高斯函数作为核函数，u为函数自变量，窗宽h的估计值可通过公式计算，其中σ为样本标准差，n为样本点数，得到的分布如图2 所示。

步骤三：将基准状态分布的统计分布f(x)和当前架次分布的统计分布g(x) 在同一坐标轴上进行表示，如图3所示。将f(x)和g(x)的重叠面积作为用于评价当前飞行架次状态变化趋势的评价指标。针对每一特征参数，f(x)和g(x)重叠度S的计算公式为：S即为用于评价当前飞行架次状态变化趋势的评价指标，其值在(0，1)区间范围内呈递减的趋势。

步骤四：建立归一化函数：其中S为计算得出的分布重叠度值，DI＝f(S)为分布偏离值，a为归一化参数，设定动部件失效对应的分布重叠度值为10^-4，对应的分布偏离值为0.99，则归一化参数a的计算公式为：如图4所示。

步骤五：以直升机初始飞行前25小时的直升机飞行实测数据利用步骤一确定飞行参数集，利用步骤二确定每个振动特征参数的基准分布。Y为当前状态下直升机动部件振动特征参数实测数据，利用步骤二得到当前状态下每个振动特征参数的统计分布，利用步骤三得到当前状态下振动特征分布与基准分布之间的重叠度，采用对数归一化函数将重叠度映射为分布偏离值DI。

Claims (6)

1.一种应用于直升机动部件振动趋势分析的方法，其特征在于包括如下步骤：

1)基于提取的直升机飞行振动特征参数和对应的飞行状态参数，采用关联分析方法，选取与飞行振动特征参数相关的飞行状态参数，并建立对应的飞行状态参数集Φ；

2)根据选取的飞行状态参数，基于直升机初始飞行架次数据建立飞行振动特征参数的基准分布数据库；

3)计算当前飞行架次下的直升机飞行振动特征参数的统计分布，同时根据当前飞行架次飞行状态参数集Φ中飞行状态参数的分布选取基准分布数据建立当前架次的基准分布，计算当前架次飞行振动特征参数统计分布与基准分布的重叠面积；

4)建立对数归一化函数，设定失效状态下直升机飞行振动特征参数的统计分布与基准分布的重叠面积为S，分布偏离值为DI，对S和DI赋值计算确定对数归一化函数参数的取值；

5)基于上述对数归一化函数得到飞行振动特征参数分布偏离值，作为趋势分析的指标。

2.根据权利要求1所述的基于分布的直升机振动趋势分析方法，其特征在于：步骤1)中采用对飞行状态参数和飞行振动特征参数进行关联分析，剔除非相关的飞行状态参数；具体方法为：

1.1，利用标准化后的样本数据，计算飞行振动特征参数变化Y与每个飞行状态参数的变化X_i之间的相关性系数，并根据其大小选择飞行状态参数，相关性系数计算公式为其中为飞行状态参数的方差，σ_Y为飞行振动特征参数的方差；x′_ij为第i类参数的第j个样本标准化后数据值，y′_j为飞行振动特征参数的第j个样本标准化后数据值，μ′_i为第i类参数的标准化样本的均值，μ_y′为参数的标准化样本的均值，n为样本量；

1.2，设飞行振动特征参数变化Y为随机变量，经相关性分析后选取的飞行状态参数变化为自变量X_i，且飞行状态参数与飞行振动特征参数的关系模型满足多变量线性回归模型，采用基于回归分析的方法建立飞行状态参数与飞行振动特征参数的关系模型b为该回归模型函数的相关系数，设显著性水平为α，计算各飞行状态参数的P值，若P值小于或等于显著性水平α，则以对应的飞行状态参数建立飞行状态参数集Φ。

3.根据权利要求1所述的基于分布的直升机振动趋势分析方法，其特征在于：步骤3)中，计算当前架次飞行振动特征参数统计分布与基准分布的重叠面积，具体方法为：

3.1针对当前飞行架次采用核密度估计的方法对其分布进行估计；

设X₁,X₂,…,X_n是取自一元连续总体的样本，在任意点x处的统计分布为n为样本点数，其中选择高斯函数作为核函数，u为核函数自变量，窗宽h的估计值通过公式计算，其中σ为样本标准差；

3.2根据当前飞行架次中飞行状态参数集Φ中参数变化范围，选取基准分布数据建立基准状态的统计分布g(x)；

3.3将基准状态分布的统计分布g(x)和当前架次分布的统计分布f(x)在同一坐标轴上进行表示，将f(x)和g(x)的重叠面积即重叠度作为用于评价当前飞行架次状态变化趋势的评价指标；针对每一特征参数，f(x)和g(x)重叠度S的计算公式为：S即为用于评价当前飞行架次状态变化趋势的评价指标，其值在(0,1)区间范围内呈递减的趋势。

4.根据权利要求1所述的基于分布的直升机振动趋势分析方法，其特征在于：步骤4)中建立对数归一化函数，通过设定失效所对应的分布重叠度值和分布偏离值，确定对数归一化函数参数，具体为：

建立对数归一化函数：其中S为计算得出的分布重叠度值，DI＝f(S)为分布偏离值，a为归一化参数，计算公式为：设定动部件失效对应的分布重叠度值为10^-4，对应的分布偏离值为0.99。

5.根据权利要求1所述的基于分布的直升机振动趋势分析方法，其特征在于：步骤5)中，采用对数归一化函数将重叠度映射为分布偏离值DI。

6.根据权利要求1或2所述的基于分布的直升机振动趋势分析方法，其特征在于：飞行状态参数包括旋翼转速、发动机扭矩、大气压力、飞行高度、指示空速、机外大气温度、偏航角、俯仰角、横滚角和三方向的加速度。