一种变工况大负荷下滚动轴承故障诊断方法技术

本发明专利技术公开了一种变工况大负荷下滚动轴承故障诊断方法，包括如下步骤，设定轴承加工过程中的变工况条件，获得训练样本数据；对样本数据进行划分训练集和测试集，并对振动信号进行归一化；搭建ResNet18和BiLSTM相结合的特征提取器，用Morlet小波核卷积层替代ResNet18的本身所含有的第一层卷积神经网络；将注意力机制下的软阈值化作为非线性层插入到ResNet18的残差块中；将经过特征提取器提取到的特征向量输入到多特征空间分布对齐的域自适应网络；将经过多特征空间对齐的域自适应网络后的特征向量输入到全连接层分类器进行滚动轴承故障识别。本发明专利技术无需大量标记样本，仅收集部分标记工况样本，即可实现对其他无标记工况数据的故障诊断，具有诊断结果准确的优点。点。点。

【技术实现步骤摘要】
一种变工况大负荷下滚动轴承故障诊断方法


[0001]本专利技术涉及加工过程监测
，特别是一种变工况大负荷下滚动轴承故障诊断方法。

技术介绍

[0002]滚动轴承广泛应用于旋转机械，如泵、涡轮、齿轮箱、压缩机、发动机等，主要起到支撑和保护机械的作用。在变工况大负荷工作条件下，机械的连续运行会导致运行部件不断与轴承碰撞，从而导致轴承损坏。据统计，轴承故障占所有电机故障的40％
‑
50％，轴承一旦发生故障，如果没有及时发现和修复，势必会导致更严重的机械故障，影响机械的正常运行，甚至引发事故。因此，创建一种在变工况大负荷下有效的轴承故障识别方法具有重要的意义和进一步开发的价值。
[0003]智能故障诊断已经成功应用于轴承故障的识别，无需人工干扰和专家的先验知识。在轴承故障诊断方面，深度学习已经超越了支持向量机(SVM)，人工神经网络(ANN)等传统机器学习，在故障识别领域取得了优异的成就。但是深度学习成功应用的前提是学习数据和测试数据需要遵循相同的特征分布，或者说来自相同的工程条件。但是实际轴承运行条件往往复杂而且随时间变化，并且不同传感器的测量能力不同，对不同类型的轴承故障信号的敏感程度不同，这导致用于网络训练的源域数据和用于实际诊断的目标域中的轴承故障数据之间的数据分布不同。
[0004]此外，在不同的工作条件下获取各种轴承故障标记数据是非常耗时和耗能的。也就是说，目标域工作条件下的轴承故障数据通常是无标记的。因此，在可变的工作条件下，传统的深度学习方法无法保证在标记的源工作条件数据上工作良好的网络对未标记的目标工作条件数据具有相同的效果。这就导致了深度学习在实际工程场景中难以应用。因此针对轴承不同工况下的特征分布差异，为了有效地将用标记的源域数据训练的网络应用于未标记的目标域数据是当前智能故障诊断应用的重大障碍。

技术实现思路

[0005]本专利技术的目的在于，提供一种变工况大负荷下滚动轴承故障诊断方法。本专利技术无需大量标记样本，仅收集部分标记工况样本，即可实现对其他无标记工况数据的故障诊断，具有诊断结果准确的优点。
[0006]本专利技术的技术方案：一种变工况大负荷下滚动轴承故障诊断方法，包括如下步骤：
[0007]步骤1、设定轴承加工过程中的变工况条件，采集设备运行过程中的轴承运作的振动信号及对应的轴承故障状态，获得少量实验训练样本数据；
[0008]步骤2、对样本数据进行划分训练集和测试集，并对振动信号进行归一化；
[0009]步骤3、搭建ResNet18和BiLSTM相结合的特征提取器，构建Morlet小波核卷积层，用Morlet小波核卷积层替代ResNet18的本身所含有的第一层卷积神经网络，作为特征提取器的第一层；
[0010]步骤4、将注意力机制下的软阈值化作为非线性层插入到ResNet18的残差块中，加强特征提取器的特征提取能力；
[0011]步骤5、将经过特征提取器提取到的特征向量输入到多特征空间分布对齐的域自适应网络，实现源域和目标域的特征分布对齐；
[0012]步骤6、将经过多特征空间对齐的域自适应网络后的特征向量输入到全连接层分类器进行滚动轴承故障识别。
[0013]上述的变工况大负荷下滚动轴承故障诊断方法，步骤2中，归一化的公式如下：
[0014][0015]式中：x
i
为输入数据；为x
i
的均值；为x
i
的标准差；为归一化数据。
[0016]前述的变工况大负荷下滚动轴承故障诊断方法，步骤3中，所述BiLSTM和ResNet18相结合的特征提取器是将ResNet18作为前端特征提取器，提取振动信号中的局部特征，然后将提取出的特征序列作为BiLSTM的输入，利用BiLSTM对时序关系进行建模和预测。
[0017]前述的变工况大负荷下滚动轴承故障诊断方法，所述Morlet小波核卷积层的构建过程是：
[0018]定义小波基函数：
[0019][0020]式中：t为时间，S为与频率成反比的尺度参数，u为小波基函数的平移参数，ψ表示预定义函数；
[0021]小波核驱动的卷积操作由以下方程定义：
[0022]h＝ψ
u，s
(t)*x；
ꢀꢀ
(3)
[0023]式中：h表示卷积后的结果；
[0024]在梯度反向传播的过程中，小波核驱动的卷积层作为神经网络的第一层，更新尺度参数S和平移参数u，如下所示：
[0025][0026][0027]其中，是求导，是长度为L的小波核卷积层第k个小波核，u
k
和S
k
表示相应的平移参数和尺度参数；
[0028]所述Morlet小波核卷积层的尺度参数S和平移参数u的更新通过将偏导数
和带入式(4)来获得。
[0029]前述的变工况大负荷下滚动轴承故障诊断方法，所述偏导数和的获取是：
[0030]计算Morlet小波的实部，如下：
[0031][0032]其中，C是信号重构的归一化常数；
[0033]再根据式(2)获取小波字典：
[0034][0035]计算偏导数，如下：
[0036][0037]前述的变工况大负荷下滚动轴承故障诊断方法，步骤4中，所述软阈值化的公式表达为：
[0038][0039]式中：x是输入特征；y是输出特征；τ是阈值；
[0040]所述残差块中设有用于估计软阈值化所需要的估计模块，在估计模块中，对特征映射的绝对值进行GAP运算，得到一维向量；然后将一维向量传播到二层FC网络中，得到一个缩放参数，在二层FC网络的末端应用Sigmod函数，将缩放参数缩放到(0.1)的范围，该范围表示为：
[0041][0042]其中，Z对应二层FC网络的输出，α对应缩放参数；
[0043]再将α乘以|x|的平均值，得到阈值：
[0044][0045]式中：i，j和c分别为特征映射的宽度、高度和通道的索引。
[0046]前述的变工况大负荷下滚动轴承故障诊断方法，步骤5中，所述域自适应网络进行源域和目标域的特征分布对齐的过程首先进行源域和目标域的特征分布的二阶统计量保持一致，二阶统计量对齐损失定义为源域特征和目标域特征的二阶统计量之间的距离可以被定义为：
[0047][0048]其中表示协方差矩阵之间的Frobenius范数，为源域特征，为目标域特征，d为每个样本的维度；C
S
为源域数据的协方差矩阵；C
T
为目标域数据间的协方差矩阵，定义如下；
[0049][0050][0051]其中，1表示所有元素都为1的列向量；源域数据个数为n
S
，目标域数据个数为n
T
；X
s
表示源域所有样本的并集；X
t
示目标域所有样本的并集；
[0052]实现初步特征分布对齐后，使用本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种变工况大负荷下滚动轴承故障诊断方法，其特征在于：包括如下步骤：步骤1、设定轴承加工过程中的变工况条件，采集设备运行过程中的轴承运作的振动信号及对应的轴承故障状态，获得少量实验训练样本数据；步骤2、对样本数据进行划分训练集和测试集，并对振动信号进行归一化；步骤3、搭建ResNet18和BiLSTM相结合的特征提取器，构建Morlet小波核卷积层，用Morlet小波核卷积层替代ResNet18的本身所含有的第一层卷积神经网络，作为特征提取器的第一层；步骤4、将注意力机制下的软阈值化作为非线性层插入到ResNet18的残差块中，加强特征提取器的特征提取能力；步骤5、将经过特征提取器提取到的特征向量输入到多特征空间分布对齐的域自适应网络，实现源域和目标域的特征分布对齐；步骤6、将经过多特征空间对齐的域自适应网络后的特征向量输入到全连接层分类器进行滚动轴承故障识别。2.根据权利要求1所述的变工况大负荷下滚动轴承故障诊断方法，其特征在于：步骤2中，归一化的公式如下：式中：x
i
为输入数据；为x
i
的均值；为x
i
的标准差；为归一化数据。3.根据权利要求1所述的变工况大负荷下滚动轴承故障诊断方法，其特征在于：步骤3中，所述BiLSTM和ResNet18相结合的特征提取器是将ResNet18作为前端特征提取器，提取振动信号中的局部特征，然后将提取出的特征序列作为BiLSTM的输入，利用BiLSTM对时序关系进行建模和预测。4.根据权利要求1所述的变工况大负荷下滚动轴承故障诊断方法，其特征在于：所述Morlet小波核卷积层的构建过程是：定义小波基函数：式中：t为时间，s为与频率成反比的尺度参数，u为小波基函数的平移参数，ψ表示预定义函数；小波核驱动的卷积操作由以下方程定义：h＝ψ
u，s
(t)*x；
ꢀꢀꢀꢀ
(3)式中：h表示卷积后的结果；在梯度反向传播的过程中，小波核驱动的卷积层作为神经网络的第一层，更新尺度参数s和平移参数u，如下所示：
其中，是求导，是长度为L的小波核卷积层第k个小波核，u
k
和s
k
表示相应的平移参数和尺度参数；所述Morlet小波核卷积层的尺度参数s和平移参数u的更新通过将偏导数和带入式(4)来获得。5.根据权利要求4所述的变工况大负荷下滚动轴承故障诊断方法，其特征在于：所述偏导数和的获取是：计算Morlet小波的实部，如下：其中，C是信号重构的归一化常数；再根据式(2)获取小波字典：计算偏导数，如下：计算偏导数，如下...

【专利技术属性】
技术研发人员：周余庆，黑振栋，娄勇坚，王公海，钟美鹏，杨海洋，
申请(专利权)人：嘉兴南湖学院，
类型：发明
国别省市：