一种航空发动机轴承故障在线诊断方法与流程

1.本技术属于航空发动机设计领域，特别涉及一种航空发动机轴承故障在线诊断方法。


背景技术：

2.轴承作为航空发动机转子系统的重要部分，长期处于高速波动剧烈、负载大且变化明显的工况，将导致性能衰退甚至引发各类故障。因此，为了使航空发动机安全高效地运行并节省维护成本，必须掌握其运行状态和变化规律，自动准确的诊断航空发动机高速轴承故障有助于提升航空发动机转子系统的运行安全性和维修经济性。
3.现有诊断方法具有具有如下缺点：
4.1、无法识别轴承早期故障的微弱信号；
5.2、难以提取微弱的故障特征；
6.3、实际轴承故障样本少。
7.以往的故障诊断方法需依托于大量的故障诊断算法以及专家工程经验，难以适应当前轴承诊断应用。面对当前海量的数据，故障信息复杂多变，包含内外激励以及多个故障的耦合，需采用基于数据驱动的故障诊断方法直接利用采集到的参数，挖掘数据识别发动机轴承的关键特征信息识别故障。


技术实现要素：

8.本技术的目的是提供了一种航空发动机轴承故障在线诊断方法，以解决轴承故障样本少、信号微弱而导致诊断准确度不高的问题。
9.本技术的技术方案是：一种航空发动机轴承故障在线诊断方法，包括：
10.建立振动特征数据集，获取轴承不同运行工况状态下的振动数据构造工况振动特征值；对轴承不同振动数据进行时域指标提取，得到时域振动特征值；对轴承不同振动数据进行频域指标提取，得到频域振动特征值；将工况振动特征值、时域振动特征值和频域振动特征值输入振动特征数据集；
11.基于振动特征数据集，再次提取案例库数据的振动特征值输入振动特征数据集，形成离线训练集，使用pca算法对离线数据集进行降维处理，将离线训练集内的数据降维至m维，形成降维数据集；
12.对降维数据集进行训练得到第一softmax分类器；
13.建立自注意力模型，将降维数据集输入带有自注意力模型的acgan算法模型中，训练得到相应的生成器和判别器，同时将训练生成的合成数据输入至第一softmax分类器内，形成轴承故障评估模型；
14.利用轴承故障评估模型评估λ
mean
值并筛选出判定有效的合成数据，与降维后的离线训练集组合形成增强数据集，λ
mean
值计算公式为：
[0015][0016]
式中，a
amc1
为预测轴承正常正确分类占总分类数目的比例；a
amc2
为预测轴承故障分类占总分类数目的比例；f
11
为故障轴承被预测为故障的数量；f
00
为正常轴承被预测为正常的数量；f
10
为故障轴承被预测为正常的数量；f
01
为正常轴承被预测为故障的数量；λ
mean
为a
amc1
和a
amc2
的几何平均；
[0017]
对增强数据集与降维数据集进行融合得到扩充数据集，通过扩充数据集训练得到进行在线故障诊断的第二softmax分类器；
[0018]
实时采集在线数据，通过提取振动特征得到数据数据集，而后对采集的数据进行pca降维，形成在线诊断集，将在线训练集输入至第二softmax分类器内，得到轴承故障诊断模型，通过轴承故障诊断模型进行轴承在线诊断，得到轴承在线诊断结果。
[0019]
优选地，所述增强数据集的生成方法包括：
[0020]
在高斯分布中随机采样生成随机噪声，对随机噪声进行分类并将随机噪声及其对应的辅助分类标签输入至生成器，得到生成数据；
[0021]
将生成数据与通过采集轴承振动特征的真实数据输入至判别器内，形成混合数据，对混合数据进行判别，得到真/假数据与不同的故障类型数据，完成一次训练；
[0022]
分别依次更新第一softmax分类器和判别器的网络参数，再次获取生成数据和真实数据进行循环判别，并将判定有效的合成数据输入至增强数据集内，直至完成设定的循环数或者增强数据集内的数据量满足设计要求。
[0023]
优选地，所述工况振动特征值包括样本输入轴转速的平均值、样本输入轴转速的方差、样本扭矩的平均值和样本扭矩的方差；所述时域振动特征值包括加速度有效值、加速度波形指标、加速度峰值指标、加速度脉冲指标、加速度裕度指标、加速度峭度指标、速度有效值、速度波形指标、速度峰值指标、速度脉冲指标、速度裕度指标和速度峭度指标；所述频域振动特征值包括轴承所在轴转频1、2、3阶幅值和频谱中幅值前100的频率值。
[0024]
优选地，所述自注意力模型的建立方法包括：
[0025]
设置三个通道节点f(x)、g(x)和h(x)，在不同的通道节点上分别设置不同的权重矩阵wk、wq和wv，将输入x分别乘上不同的权重矩阵，得到：
[0026][0027]
计算f(x)和g(x)之间的关联矩阵s：
[0028]
s＝f
t
·g[0029]
计算f(x)和g(x)之间的相关程度：
[0030][0031]
式中：s
ij
为关联矩阵s的元素；
[0032]
得到自注意力模型的输出特征：
[0033][0034]
加入权重系数γ，进行领域信息与远距离特征的整合，得到：
[0035]
yi＝γoi+xi。
[0036]
优选地，所述acgan算法模型对降维数据集进行训练的方法为：
[0037]
定义降维数据集x＝{x1,x2,
…
,xn},其中xi为第i个样本，n为样本总数；生成器g的损失函数设置为：
[0038][0039]
将辅助分类标签c加入生成器中，得到表征数据真实与否的损失函数ls为：
[0040][0041]
判别器d的损失函数ld设置为：
[0042][0043]
设置辅助分类器c并引入至判别器d内，得到表征数据分类准确度的损失函数lc为：
[0044][0045]
引入沃瑟斯坦度量，构建整体训练的目标函数，度量生成样本和真实样本的差别，构建基于自注意力的acgan模型，则基于自注意力的acgan整体训练的目标函数为：
[0046][0047]
优选地，在建立目标函数的同时，构建梯度惩罚项：
[0048][0049]
式中：||
·
||
p
表示p范数；λ为正则项系数；其中ε～u(0,1)，u为均匀分布。
[0050]
本技术的一种航空发动机轴承故障在线诊断方法，通过先采集不同类型的振动数据，建立振动特征数据集，而后通过再次提取案例库数据的振动特征值输入该数据集，形成离线训练集，使用pca算法对离线数据集进行降维处理，对降维数据集进行训练得到第一softmax分类器；将降维数据集输入带有自注意力模型的acgan算法模型中，训练得到相应
的生成器和判别器，利用轴承故障评估模型评估λ
mean
值并筛选出判定有效的合成数据，与降维后的离线训练集组合形成增强数据集；对增强数据集与降维数据集进行融合得到扩充数据集，训练得到第二softmax分类器；将在线训练集输入至第二softmax分类器内，得到轴承故障诊断模型，通过轴承故障诊断模型进行轴承在线诊断，得到轴承在线诊断结果。可以在线做出准确、快速的判断，诊断轴承故障，具有较高识别精度和稳定性，适用于现场实时诊断应用。
附图说明
[0051]
为了更清楚地说明本技术提供的技术方案，下面将对附图作简单地介绍。显而易见地，下面描述的附图仅仅是本技术的一些实施例。
[0052]
图1为本技术整体流程示意图；
[0053]
图2为本技术自注意力模型结构示意图；
[0054]
图3为本技术增强数据集生成方法示意图。
具体实施方式
[0055]
为使本技术实施的目的、技术方案和优点更加清楚，下面将结合本技术实施例中的附图，对本技术实施例中的技术方案进行更加详细的描述。
[0056]
一种航空发动机轴承故障在线诊断方法，如图1所示，包括如下步骤：
[0057]
步骤s100，选取振动特征值输入振动特征数据集
[0058]
①
由于航空发动机在不同工况下轴承的运行状态存在一定差异，通过引入运行工况状态反应轴承运行状态差异构造工况振动特征值输入振动特征数据集。因此选取：样本输入轴转速的平均值、样本输入轴转速的方差、样本扭矩的平均值、样本扭矩的方差，共计4个工况振动特征值。
[0059]
②
对轴承振动情况进行时域指标提取，通过振动时域指标反映表征轴承的不同运行状态。因此选取：加速度有效值、加速度波形指标、加速度峰值指标、加速度脉冲指标、加速度裕度指标、加速度峭度指标、速度有效值、速度波形指标、速度峰值指标、速度脉冲指标、速度裕度指标、速度峭度指标，共计12个时域振动特征值。
[0060]
③
对轴承振动情况进行频域指标提取，通过振动频域指标反映表征轴承的不同健康状态。因此选取：轴承所在轴转频1、2、3阶幅值、频谱中幅值前100的频域振动特征值(共计103个特征值)。因此共构建119维特征构成振动特征值输入振动特征数据集。
[0061]
步骤s200，pca降维
[0062]
基于振动特征数据集，再次提取案例库数据的振动特征值输入振动特征数据集，形成离线训练集，使用pca算法对离线数据集进行降维处理，将离线训练集内的数据降维至m维，m值可根据效果自行选定数值大小，形成降维数据集。
[0063]
pca降维以函数封装，并通过传递数据集、设置降维的维数来进行数据降维，pca降维能够消除原始数据之间相互影响的因素，减少计算量，pca算法采用现有算法，具体不再赘述。
[0064]
步骤s300，对降维数据集进行训练，通过对降维数据集内的振动特征数据设置相应的有效性标准或者阈值来得到第一softmax分类器，该分类器用于评估基于自注意力的
acgan在训练过程中生成的合成数据有效性，softmax分类器的训练方法为现有技术，具体不再赘述。
[0065]
步骤s400，建立自注意力模型，将降维数据集输入带有自注意力模型的acgan算法模型中，训练得到相应的生成器和判别器，同时将训练生成的合成数据输入至第一softmax分类器内，形成轴承故障评估模型；
[0066]
传统gan为无监督学习，由两个网络即生成器和判别器组成。但其模式过于自由导致训练过程不可控。
[0067]
cgan在gan的基础上，将生成器和判别器中加入辅助分类标签c，用其来指导数据生成的方向，实现监督学习。
[0068]
acgan在cgan的基础上加以改进，标签信息只作用于生成器，将辅助分类器c引入判别器用于判别样本的类别。能够通过辅助分类器在生成过程中，控制样本生成的方向，以生成高质量的结果；但由于卷积核尺寸有限，只能学习样本局部区域的关系，模型的学习效率低下，并且可能会丢失细节。
[0069]
基于自注意力的acgan在acgan的基础上，将自注意力机制加入g和d，解决了acgan只能计算特定邻域内的数据特征的局限性，帮助模型捕捉样本长距离特征间的联系。
[0070]
优选地，如图2所示，自注意力模型的建立方法包括：
[0071]
设置三个通道节点f(x)、g(x)和h(x)，在不同的通道节点上分别设置不同的权重矩阵wk、wq和wv，将输入x分别乘上不同的权重矩阵，得到：
[0072][0073]
计算f(x)和g(x)之间的关联矩阵s：
[0074]
s＝f
t
·g[0075]
计算f(x)和g(x)之间的相关程度：
[0076][0077]
式中：s
ij
为关联矩阵s的元素；
[0078]
得到自注意力模型的输出特征：
[0079][0080]
加入权重系数γ，进行领域信息与远距离特征的整合，得到：
[0081]
yi＝γoi+xi。
[0082]
优选地，基于自注意力的acgan算法模型对降维数据集进行训练的方法为：
[0083]
在文中下标“g”表示生成器相关的参数，下标“d”表示判别器相关参数。
[0084]
定义降维数据集x＝{x1,x2,
…
,xn},其中xi为第i个样本，n为样本总数；生成器g生成接近真实分布的逼真数据，为实现从数据集x的实际分布pr(x)中学习到生成分布pg(z)。生成器g由多层神经网络构成，加入注意力机制后，得到生成器g的损失函数设置为：
[0085][0086]
将辅助分类标签c加入生成器中，用其来指导数据生成的方向，得到表征数据真实与否的损失函数ls为：
[0087][0088]
判别器d判别输入是实际量还是生成量，判别器d也由多层神经网络构成，加入自注意力机制后，得到判别器d的损失函数ld设置为：
[0089][0090]
设置辅助分类器c用于判别样本的类别，将辅助分类器c引入至判别器d内，得到表征数据分类准确度的损失函数lc为：
[0091][0092]
引入沃瑟斯坦度量(wasserstein)构建整体训练的目标函数，度量生成样本和真实样本的差别，构建基于自注意力的acgan模型，生成高质量的故障样本，用作航空发动机轴承故障诊断。则基于自注意力的acgan整体训练的目标函数为：
[0093][0094]
通过建立自注意力模型，通过对三个通道的权重大小进行动态加权融合，来增强重要通道的特征，弱化非重要通道的特征，实现全局信息的关联，从而降低计算量、提高模型训练速度。
[0095]
基于博弈论的方式生成样本更接近实际故障样本，解决实际诊断中轴承故障样本少，有效增强故障样本数量。
[0096]
优选地，针对梯度爆炸或梯度消失问题，构建梯度惩罚项：
[0097][0098]
式中：||
·
||
p
表示p范数；λ为正则项系数；其中ε～u(0,1)，u为均匀分布。
[0099]
步骤s500，利用轴承故障评估模型评估λ
mean
值并筛选出判定有效的合成数据，与降维后的离线训练集组合形成增强数据集，λ
mean
值计算公式为：
[0100]
[0101]
式中，a
amc1
为预测轴承正常正确分类占总分类数目的比例；a
amc2
为预测轴承故障分类占总分类数目的比例；f
11
为故障轴承被预测为故障的数量；f
00
为正常轴承被预测为正常的数量；f
10
为故障轴承被预测为正常的数量；f
01
为正常轴承被预测为故障的数量；λ
mean
为a
amc1
和a
amc2
的几何平均；
[0102]
优选地，如图3所示，增强数据集的生成方法包括：
[0103]
1)在高斯分布中随机采样生成随机噪声，对随机噪声进行分类并将随机噪声及其对应的辅助分类标签输入至生成器，得到生成数据；
[0104]
2)将生成数据与通过采集轴承振动特征的真实数据输入至判别器内，形成混合数据，对混合数据进行判别，得到真/假数据与不同的故障类型数据，完成一次训练；
[0105]
3)分别依次更新第一softmax分类器和判别器的网络参数，再次获取生成数据和真实数据进行循环判别，并将判定有效的合成数据输入至增强数据集内，直至完成设定的循环数或者增强数据集内的数据量满足设计要求。
[0106]
步骤s600，对增强数据集与降维数据集进行融合得到扩充数据集，通过扩充数据集训练得到进行在线故障诊断的第二softmax分类器；第二softmax分类器内通过设置故障诊断标准参数来对振动数据进行故障诊断，通过对振动特征数据集内的数据进行持续采集并在线诊断，能够持续提升故障诊断的准确度。
[0107]
步骤s700，实时采集在线数据，通过提取振动特征得到数据数据集，而后对采集的数据通过步骤s200的方式进行pca降维，形成在线诊断集，将在线训练集输入至第二softmax分类器内，得到轴承故障诊断模型，通过轴承故障诊断模型进行轴承在线诊断，得到轴承在线诊断结果。
[0108]
本技术通过先采集不同类型的振动数据，建立振动特征数据集，而后通过再次提取案例库数据的振动特征值输入该数据集，形成离线训练集，使用pca算法对离线数据集进行降维处理，以减少计算量，对降维数据集进行训练得到第一softmax分类器；将降维数据集输入带有自注意力模型的acgan算法模型中，训练得到相应的生成器和判别器，同时提升样本数量，利用轴承故障评估模型评估λ
mean
值并筛选出判定有效的合成数据，与降维后的离线训练集组合形成增强数据集；对增强数据集与降维数据集进行融合得到扩充数据集，通过扩充数据集训练得到进行在线故障诊断的第二softmax分类器；将在线训练集输入至第二softmax分类器内，得到轴承故障诊断模型，通过轴承故障诊断模型进行轴承在线诊断，得到轴承在线诊断结果。通过引入自注意力机制，能够学习原始故障没有的故障特征分布特性，提高了轴承故障样本的生成质量和学习效率；建立的轴承故障在线诊断模型具有诊断准确率高的特点，能够准确识别航空发动机轴承故障的特点，可以在线做出准确、快速的判断，诊断轴承故障，具有较高识别精度和稳定性，适用于现场实时诊断应用。
[0109]
以上所述，仅为本技术的具体实施方式，但本技术的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本技术揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本技术的保护范围之内。因此，本技术的保护范围应以所述权利要求的保护范围为准。

技术特征：
1.一种航空发动机轴承故障在线诊断方法，其特征在于，包括：建立振动特征数据集，获取轴承不同运行工况状态下的振动数据构造工况振动特征值；对轴承不同振动数据进行时域指标提取，得到时域振动特征值；对轴承不同振动数据进行频域指标提取，得到频域振动特征值；将工况振动特征值、时域振动特征值和频域振动特征值输入振动特征数据集；基于振动特征数据集，再次提取案例库数据的振动特征值输入振动特征数据集，形成离线训练集，使用pca算法对离线数据集进行降维处理，将离线训练集内的数据降维至m维，形成降维数据集；对降维数据集进行训练得到第一softmax分类器；建立自注意力模型，将降维数据集输入带有自注意力模型的acgan算法模型中，训练得到相应的生成器和判别器，同时将训练生成的合成数据输入至第一softmax分类器内，形成轴承故障评估模型；利用轴承故障评估模型评估λ
mean
值并筛选出判定有效的合成数据，与降维后的离线训练集组合形成增强数据集，λ
mean
值计算公式为：式中，a
amc1
为预测轴承正常正确分类占总分类数目的比例；a
amc2
为预测轴承故障分类占总分类数目的比例；f
11
为故障轴承被预测为故障的数量；f
00
为正常轴承被预测为正常的数量；f
10
为故障轴承被预测为正常的数量；f
01
为正常轴承被预测为故障的数量；λ
mean
为a
amc1
和a
amc2
的几何平均；对增强数据集与降维数据集进行融合得到扩充数据集，通过扩充数据集训练得到进行在线故障诊断的第二softmax分类器；实时采集在线数据，通过提取振动特征得到数据数据集，而后对采集的数据进行pca降维，形成在线诊断集，将在线训练集输入至第二softmax分类器内，得到轴承故障诊断模型，通过轴承故障诊断模型进行轴承在线诊断，得到轴承在线诊断结果。2.如权利要求1所述的航空发动机轴承故障在线诊断方法，其特征在于，所述增强数据集的生成方法包括：在高斯分布中随机采样生成随机噪声，对随机噪声进行分类并将随机噪声及其对应的辅助分类标签输入至生成器，得到生成数据；将生成数据与通过采集轴承振动特征的真实数据输入至判别器内，形成混合数据，对混合数据进行判别，得到真/假数据与不同的故障类型数据，完成一次训练；分别依次更新第一softmax分类器和判别器的网络参数，再次获取生成数据和真实数据进行循环判别，并将判定有效的合成数据输入至增强数据集内，直至完成设定的循环数或者增强数据集内的数据量满足设计要求。3.如权利要求1所述的航空发动机轴承故障在线诊断方法，其特征在于：所述工况振动
特征值包括样本输入轴转速的平均值、样本输入轴转速的方差、样本扭矩的平均值和样本扭矩的方差；所述时域振动特征值包括加速度有效值、加速度波形指标、加速度峰值指标、加速度脉冲指标、加速度裕度指标、加速度峭度指标、速度有效值、速度波形指标、速度峰值指标、速度脉冲指标、速度裕度指标和速度峭度指标；所述频域振动特征值包括轴承所在轴转频1、2、3阶幅值和频谱中幅值前100的频率值。4.如权利要求1所述的航空发动机轴承故障在线诊断方法，其特征在于，所述自注意力模型的建立方法包括：设置三个通道节点f(x)、g(x)和h(x)，在不同的通道节点上分别设置不同的权重矩阵w
k
、w
q
和wv，将输入x分别乘上不同的权重矩阵，得到：计算f(x)和g(x)之间的关联矩阵s：s＝f
t
·
g计算f(x)和g(x)之间的相关程度：式中：s
ij
为关联矩阵s的元素；得到自注意力模型的输出特征：加入权重系数γ，进行领域信息与远距离特征的整合，得到：y
i
＝γo
i
+x
i
。5.如权利要求1所述的航空发动机轴承故障在线诊断方法，其特征在于，所述acgan算法模型对降维数据集进行训练的方法为：定义降维数据集x＝{x1,x2,
…
,x
n
},其中x
i
为第i个样本，n为样本总数；生成器g的损失函数设置为：将辅助分类标签c加入生成器中，得到表征数据真实与否的损失函数l
s
为：判别器d的损失函数l
d
设置为：设置辅助分类器c并引入至判别器d内，得到表征数据分类准确度的损失函数l
c
为：引入沃瑟斯坦度量，构建整体训练的目标函数，度量生成样本和真实样本的差别，构建
基于自注意力的acgan模型，则基于自注意力的acgan整体训练的目标函数为：6.如权利要求5所述的航空发动机轴承故障在线诊断方法，其特征在于，在建立目标函数的同时，构建梯度惩罚项：式中：||
·
||
p
表示p范数；λ为正则项系数；其中ε～u(0,1)，u为均匀分布。

技术总结
本申请属于航空发动机设计领域，为一种航空发动机轴承故障在线诊断方法，先建立振动特征数据集和离线训练集，使用PCA算法对离线数据集进行降维处理，对降维数据集进行训练得到第一Softmax分类器；将降维数据集输入带有自注意力模型的ACGAN算法模型中，训练得到相应的生成器和判别器，利用轴承故障评估模型评估λ


技术研发人员：郭小鹏 刘永泉 何俊杰 杜少辉 安中彦 武鹏 张效语 李修文 金贻涛
受保护的技术使用者：中国航发沈阳发动机研究所
技术研发日：2023.05.04
技术公布日：2023/8/24