一种滚动轴承故障智能诊断方法

1.本发明属于机械故障诊断技术领域，涉及一种滚动轴承故障智能诊断方法。


背景技术：

2.滚动轴承是现代旋转机械中的重要零部件，其主要功能为支撑旋转部件，减少旋转体之间的摩擦损失。然而，轴承也是最容易损坏的部件，据统计，旋转机械失效案例中45％-55％是由滚动轴承失效引起的。滚动轴承的故障可能会导致整个机械系统的损坏，造成大量的经济损失。因此，为了保证设备运行时的安全性和可靠性，对滚动轴承进行故障诊断和在线监测是很有必要的。传统的机械系统健康监测方法需要手动提取故障特征，在面对故障信息丰富多变、实际工况交替复杂时往往会导致诊断效果不理想和泛化性能欠佳的问题。因此，研究既能避免人工提取特征方法的不足又具有强大适应能力的智能故障诊断方法具有重要意义。
3.传统智能故障诊断方法通过对原始振动信号的时频域分析来完成对信号特征的提取和故障识别。不少学者作出了大量的研究。乔志城等采用改进经验小波变换(iewt)获取信号的模态分量并完成重要分量的重构，再通过最小熵解卷积(med)对重构信号进行滤波，最后完成轴承的故障诊断。何勇等通过遗传算法(ga)优化变分模态分解(vmd)，之后通过模态分量的包络谱判断轴承故障类型。saidi等利用经验模态分解(emd)将信号分解为本征模态函数(imfs)，再对imf信号进行双谱分析，最后提取轴承故障特征频率实现故障诊断。这些方法虽然能有效诊断轴承故障，但需根据经验人工提取故障特征，往往存在鲁棒性差，泛化性能不强，技术人员难以胜任故障诊断任务等问题。
4.传统智能故障诊断方法适应能力较弱，在复杂多变的工况下表现较差，且过分依赖于诊断专家的经验，因此，需要研究出一种适应能力强、受人工提取特征影响较小的故障诊断方法来解决传统方法的不足。2006年，hinton等正式提出了深度学习的概念，克服了人工提取特征的困难，在计算机视觉、语音识别、自然语言处理、图像识别等领域已取得了丰硕的成果。随着深度学习的发展，越来越多的学者将其引用到故障诊断领域中。李巍华等将深度信念网络(dbn)直接应用于轴承振动原始信号的处理，实现轴承故障的分类识别。verma等利用堆叠式稀疏自编码器(sae)提取稀疏特征，对轴承以及阀门进行故障诊断。
5.卷积神经网络(convolutional neural network，cnn)是深度学习的代表算法之一，因其具有强大的自动特征提取功能而被广泛应用于图像识别。近年来，cnn被故障诊断领域的部分学者应用于故障分类任务中。2016年，janssens等首次利用cnn对齿轮箱中的轴承、齿轮等进行故障诊断，将常规算法故障诊断率提高了6％。宫文峰等通过实验全面地分析了cnn中比较有效的深度学习技巧、模型超参数对滚动轴承故障诊断结果的影响，但并未考虑负载改变时的诊断效果。udmale等将包含时频能量的谱峭度图输入到cnn中进行故障诊断，提高了轴承在变工况下的识别精度，但并未考虑数据集规模大小对模型的影响。liang等提出了一种并行卷积神经网络(p-cnn)，融合了时域特征和频域特征，在不同负载和少量故障数据时均获得了较高的识别准确率，但该模型结构复杂，不利于故障的快速诊
断。wen等将原始振动信号转化为二维灰度图像，输入cnn中进行故障诊断取得了较好的诊断效果，但该方法并未考虑信号的时间相关性且存在故障信号信息丢失的问题。


技术实现要素：

6.本发明的目的在于针对现有技术中的缺陷，提供了一种滚动轴承故障诊断方法，能够顾在数据集较小、负载改变的环境下对滚动轴承轻微故障诊断具有更好的泛化性能和识别效果。
7.为了实现上述目的，本发明采用的技术方案是：
8.一种滚动轴承故障智能诊断方法，包括如下步骤：
9.s1：获取滚动轴承原始振动信号，用于cnn模型的训练；
10.s2：对采集到的信号按设置的样本长度进行随机分割；
11.s3：将分割后的每段信号按马尔科夫转移场(mtf)编码方式转换为二维特征图像；
12.s4：将二维特征图像数据集按比例划分为训练集和测试集；
13.s5：搭建cnn网络模型，初始化参数；
14.s6：把训练集输入cnn模型进行预训练，利用优化算法反复对模型参数进行优化，若达到最优值则进行步骤s7，否则跳转到步骤s5，修改模型参数，直到获得最优参数为止，并保存最佳模型；
15.s7：将测试集输入到训练好的cnn模型中进行滚动轴承故障诊断，最终获得轴承故障分类结果及测试准确率。
16.所述步骤s1中通过轴承实验数据中心获取滚动轴承原始振动信号，所述轴承数据实验中心包括依次连接的三相异步电机、扭矩传感器和负载，所述三相异步电机的底座外壳上安装有用于采集故障轴承振动型号的加速度传感器。
17.所述步骤s1中驱动端轴承在载荷为1-3hp三种工况下工作时采集的故障振动信号，振动信号的采样频率为12khz，采样时间为10s。
18.所述步骤s1中轴承的最低转速为1772r/min，通过公式
19.n＝fs×
60/ns20.计算出一个周期内的采样点数。
21.所述步骤s2中采用重叠采样的方式对数据进行随机分割。
22.所述步骤s3中通过马尔科夫转移概率来实现动态转移信息的编码，具体操作步骤为：
23.s301：给定一个时间序列x＝{x1,x2,xi,
…
xn}，首先将x离散化为q个分位数单元，用分位数qj(j∈[1,q])量化时间序列的每一个值，通过识别分位数q，将每个数值xi映射到相应的分位数qi上，并构造一个q
×
q的马尔科夫转移矩阵w，通过沿时间轴以一阶马尔可夫链的方式将分位数转换成矩阵w，其表达式如下式所示：
[0024][0025]
式中：w
ij
表示分位数qi位于点分位数qj后的概率，即
[0026]wij
＝p(x
t
∈qi|x
t-1
∈qj)；
[0027]
s302：由于马尔科夫转移矩阵忽略了时间序列x的分布与时间步长ti之间的时间依赖关系，故构造mtf来弥补这一不足，mtf通过沿时间顺序排列每个概率来扩展矩阵w，即为矩阵m，表达式如下式所示，矩阵m包含了每个分位数与时间步长之间的时间相关性：
[0028][0029]
式中：m
ij
表示分位数qi转移到分位数qj的转移概率，即m
ij
＝p(qi→
qj)。
[0030]
所述步骤s4中训练集和测试集的比例为3:1。
[0031]
所述步骤s5中cnn网络模型包括卷积层、池化层、全连接层和softmax分类器，具体是通过多个能提取输入数据特征的滤波器对输入的图像进行卷积和池化运算，逐层提取出具有显著特征的特征图，将获得的特征图通过全连接层，再对模型进行误差计算和参数优化，最后通过分类器输出相应的结果。
[0032]
所述步骤s6中利用损失函数通过对预测样本和真实样本标签产生的误差进行反向传播从而优化网络参数，在反向传播的过程中，使用adam自适应优化器对网络参数进行优化，利用参数的一阶矩估计和二阶矩估计动态调整每一个参数的学习率从而更新权值与偏差，计算求得参数的最优解。
[0033]
与现有技术相比，本发明具有如下有益效果：
[0034]
1)本发明首先通过mtf将原始信号可视化，即转化为包含故障特征的二维图像，然后将特征图像作为cnn的输入进行特征提取，最后实现对不同滚动轴承故障类型的识别分类，充分结合了mtf编码方式可以保留信号时间相关性和cnn强大的自动特征提取功能的优点，解决了信号故障信息容易丢失的问题，故障识别准确率高；
[0035]
2)本发明创建的mtf-cnn模型在负载改变时仍然保持了较高的诊断精度，具有更好的泛化性能，并且该模型在数据集规模减小时依然获得了较高的识别准确率，显著提高了滚动轴承故障识别效果的准确性和稳定性。
附图说明
[0036]
图1是本发明实施例的流程图。
[0037]
图2是本发明实施例的重叠采样数据增强示意图。
[0038]
图3是本发明实施例中7种轴承故障的mtf编码特征图。
[0039]
图4是本发明实施例模型训练后的准确率曲线图。
[0040]
图5是本发明实施例模型训练后的损失曲线图。
[0041]
图6是本发明实施例的故障诊断分类结果混淆矩阵。
[0042]
图7是本发明中各算法在不同负载时的识别准确率。
[0043]
图8是本发明中各算法在不同数据集规模下的识别准确率箱型图。
具体实施方式
[0044]
下面将结合具体实施例来详细说明本发明，在此本发明的示意性实施例以及说明来解释本发明，但并不作为对本发明的限定。
[0045]
如图1所示，本发明公开了一种滚动轴承故障智能诊断方法，包括如下步骤：
[0046]
s1：获取滚动轴承原始振动信号，用于cnn模型的训练；
[0047]
s2：对采集到的信号按设置的样本长度进行随机分割；
[0048]
s3：将分割后的每段信号按mtf编码方式转换为二维特征图像；
[0049]
s4：将二维特征图像数据集按比例划分为训练集和测试集；
[0050]
s5：搭建cnn网络模型，初始化参数；
[0051]
s6：把训练集输入cnn模型进行预训练，利用优化算法反复对模型参数进行优化，若达到最优值则进行步骤s7，否则跳转到步骤s5，修改模型参数，直到获得最优参数为止，并保存最佳模型；
[0052]
s7：将测试集输入到训练好的cnn模型中进行滚动轴承故障诊断，最终获得轴承故障分类结果及测试准确率。
[0053]
为更好地理解上述方案，特做以下说明：
[0054]
1.马尔科夫转移场
[0055]
马尔科夫转移场(markov transition field，mtf)可将原始一维振动信号转化为新的时间序列二维图像，该方法主要是通过马尔科夫转移概率来实现动态转移信息的编码。
[0056]
给定一个时间序列x＝{x1,x2,xi,
…
xn}，首先将x离散化为q个分位数单元，用分位数qj(j∈[1,q])量化时间序列的每一个值，通过识别分位数q，将每个数值xi映射到相应的分位数qi上，并构造一个q
×
q的邻近加权矩阵w(马尔科夫转移矩阵)，通过沿时间轴以一阶马尔可夫链的方式将分位数转换成矩阵w，其表达式如式(1)所示
[0057][0058]
式中：w
ij
表示分位数qi位于点分位数qj后的概率，即w
ij
＝p(x
t
∈qi|x
t-1
∈qj)。
[0059]
由于马尔科夫转移矩阵忽略了时间序列x的分布与时间步长ti之间的时间依赖关系，故构造mtf来弥补这一不足，mtf通过沿时间顺序排列每个概率来扩展矩阵w，即为矩阵m，表达式如式(2)所示，矩阵m包含了每个分位数与时间步长之间的时间相关性。
[0060][0061]
式中：m
ij
表示分位数qi转移到分位数qj的转移概率，即m
ij
＝p(qi→
qj)。
[0062]
mtf通过上述编码方式将时间序列可视化，转换为具有时间相关性的二维特征图像，mtf图像编码方式主要由以下几个优点：(1)通过考虑每个分位数与时间步长之间的依赖关系，保留了原始信号在不同时间间隔内的时间相关性；(2)由于一维信号与mtf二维图像编码方式为映射关系，避免的一维信号信息的丢失；(3)用不同深浅的颜色体现出分位数之间的转移概率大小，有利于充分利用cnn在图像分类中的优势。
[0063]
2.卷积神经网络
[0064]
卷积神经网络(cnn)是一类包含卷积计算且具有深度结构的前馈神经网络，cnn模型由卷积层、池化层、全连接层和softmax分类器组成，其原理是通过多个能提取输入数据特征的滤波器对输入的图像进行卷积和池化运算，逐层提取出具有显著特征的特征图，将获得的特征图通过全连接层，再对模型进行误差计算和参数优化，最后通过分类器输出相应的结果。
[0065]
2.1卷积层
[0066]
卷积层也叫作特征提取层，是cnn的核心组成部分。卷积层通过不同尺寸的卷积核与上一层输入图像进行卷积运算，得到相应的激活图像，通过激活函数提取图像特征。卷积运算的公式如式(3)所示：
[0067][0068]
式中：l表示对应网络中的第l层；表示在第l层的第j个特征图；f表示激活函数；mj表示与第l层第j个特征图相连接的第l-1层中的特征图；表示第l层第j个特征图对应第l-1层第i个特征图输入的卷积核参数；表示第l层的偏置；*表示卷积操作。
[0069]
2.2池化层
[0070]
池化层又称为下采样层，一般与卷积层相连接，其主要作用是保留主要特征的基础上降低特征维度，并减少训练参数，降低计算量，防止过拟合。池化运算的数学模型可表示为：
[0071][0072]
式中：表示输入特征图和池化层之间的连接权重；down(
·
)表示下采样函数，本发明选用的下采样函数为最大采样函数。
[0073]
2.3全连接层
[0074]
全连接层通常与最后一个池化层相连，可以整合前一层中具有类别特征的局部信息，起到收集特征和分类的作用，全连接层进行的操作为：
[0075]
δi＝f(aiσi+bi)
ꢀꢀꢀ
(5)
[0076]
式中：i＝1,2,
…
,k；δi表示第i个输出；ai和bi分别表示第i个神经元的权值和阈值；σi表示前一层的输出；f表示激活函数。
[0077]
2.4 softmax分类器
[0078]
softmax分类器位于全连接层之后，以softmax函数作为激活函数，完成对目标的分类。softmax函数的数学表达式为：
[0079][0080]
式中：i＝1,2,
…
,k；δi表示上一层第i个节点输出，k为输出节点的个数，即分类的类别个数；pi表示每一个输出的概率，所有pi之和为1。
[0081]
2.5损失函数与训练优化算法
[0082]
损失函数通过对预测样本和真实样本标签产生的误差进行反向传播从而优化网络参数。在cnn多分类问题中，常用交叉熵函数作为softmax分类函数的损失函数来衡量模型输出结果与真实结果之间的差异性，交叉熵值越小，模型预测效果越好。其公式如式(7)所示：
[0083][0084]
式中：p(xi)表示样本的真实概率分布，q(xi)表示模型预测的概率分布。
[0085]
在反向传播的过程中，使用adam自适应优化器对网络参数进行优化，利用参数的一阶矩估计和二阶矩估计动态调整每一个参数的学习率从而更新权值与偏差，计算求得参数的最优解。
[0086]
3.滚动轴承故障诊断实验验证与分析
[0087]
3.1实验数据集
[0088]
步骤s1中通过轴承实验数据中心获取滚动轴承原始振动信号，该轴承实验数据中心主要由三相异步电机、扭矩传感器和负载组装而成，并将加速度传感器安装在电机磁性底座外壳上用于采集故障轴承的振动信号。其中选取试验台驱动端轴承在载荷为1-3hp(1hp＝0.7457kw)三种工况下工作时采集的故障振动信号，信号采样频率为12khz，采样时间为10s。轴承故障形式为点蚀，轴承的内圈、外圈和滚动体表面故障均由电火花技术加工而成，选取数据库中具有代表性的两种轴承轻微故障类型，故障直径为0.18mm和0.36mm，另外包括正常轴承共计7种故障类型。
[0089]
其中轴承最低转速为1772r/min，通过公式n＝fs×
60/ns计算出一个周期内的采样点数为406个，为充分保证每个数据样本故障信息的完整性和可信度，本发明设置每个数据样本长度为2048个采样点，并采用重叠采样的方式对数据进行随机分割。将7种原始一维振动信号进行样本数据标签化，每类信号包含400个样本，共计2800个样本，并按照3:1的比例划分为训练集和测试集，具体轴承故障样本分布如表1所示。
[0090]
表1滚动轴承故障样本分布
[0091][0092]
3.2特征图像的生成
[0093]
3.2.1数据增强
[0094]
卷积神经网络作为深度学习的代表算法之一，往往需要大量的数据样本才能保证计算的准确性，在现有数据量有限的情况下，采用数据增强的方法满足训练所需样本的数量，在增加样本多样性的同时，提高模型的鲁棒性，防止发生过拟合。本发明采用重叠采样的方式对原始一维振动信号进行数据增强，示意图如图2所示，其中实线框内整个信号表示原始信号，每个虚线框内的信号即为样本信号。
[0095]
3.2.2 mtf特征图像的生成
[0096]
传统的时域分析难以清晰地识别出滚动轴承的损伤程度，无法保证对故障类型的精确分类，因此，利用mtf编码方式考虑信号时间相关性的优点，将信号按设置的数据样本长度随机分割为不同的小段并生成mtf特征图像，如图3所示。
[0097]
3.3实验结果分析
[0098]
根据关键参数对模型识别率的影响分析，选取batch size＝10，kernel size＝5，并设置epoch＝30，卷积层采用relu激活函数，参数优化器选取adam自适应优化器，设置学习率为0.001，dropout的值为0.5，即每次训练随机选取一半神经元失活，其余一半用来正常训练，以防止网络结构过拟合。
[0099]
选取负载为1hp即数据集a进行测试，模型经过30轮训练得到如图4-5所示的结果。分析图4-5可知，模型在训练集上准确率达到99.86％，在测试集上达到99.71％，经过15轮训练后，训练损失率降到0.0068并趋于平稳。验证了mtf-cnn模型在滚动轴承故障诊断中具有较好的诊断效果。
[0100]
为了进一步考察mtf-cnn模型对滚动轴承故障识别的能力，引入混淆矩阵对诊断结果进行量化分析。混淆矩阵详细地反映了滚动轴承不同故障类型下的误判数量和误判类型，如图6所示。图6中横坐标代表不同故障的预测类别，纵坐标代表不同故障的真实标签，矩阵主对角线上的数字代表对每类故障正确分类的样本数量。
[0101]
由图6明显可知：700个测试集样本只有2个样本被错分，均为故障尺寸为0.36mm的外圈故障被误判为故障尺寸为0.18mm的滚动体故障。说明了mtf-cnn模型在滚动轴承故障诊断中具有优异的故障识别能力。
[0102]
3.4模型在不同负载情况下的泛化性能分析
[0103]
在实际工况中，滚动轴承的负载往往会发生改变，因此，在不同负载情况下故障诊断模型依然能具有良好的诊断效果是非常重要的。为了验证本发明所提出模型在不同负载时的优越性，选用三种不同负载下故障数据分别构建训练集和测试集，例如图7中的b
→
c表示以数据集b作为训练集、数据集c为测试集，测试结果均取5次计算的平均值。
[0104]
本发明与svm、bpnn和基于时域分析的cnn三种常用的智能故障诊断方法进行对比分析，最终结果如图7所示。从图7可以看出，各个算法在负载相同时的诊断结果均比不同负载时要好，本发明的识别准确率均比其他方法要高，以负载为3hp即表中c
→
c为例，本发明的识别准确率为99.83％，比其他算法中效果最好的cnn高0.59％，比svm高5.23％，比效果最差的bpnn高16.52％。在负载改变时，本发明诊断性能的优越性更加凸显。由图7可知，本发明在不同负载时的平均识别准确率均高于其他方法。cnn的诊断效果相较于其他常用算法更好，在前四种负载改变的工况下，均能保持97％以上的识别准确率，在a
→
b时最高，达到98.21％，但是依然低于本发明的99.68％，同样在这四种情况下，本发明的识别率均在98％以上，在b
→
a时更是达到99.74％，而在c
→
a和c
→
b情况下，cnn识别率降低到了86％以下，分别为85.26％和85.74％，本发明则保持稳定，达到了98.6％和99.69％，均远高于cnn；bpnn在不同负载时效果最差，识别准确率均不如其他三种算法；svm只有在a
→
c时的识别准确率高于80％，其余情况均属于70％与80％之间，均低于本发明。此外，从图7的误差线可以看出，本发明的误差均低于其他三种方法，主要原因是当负载改变时，三种常用算法的诊断结果的不确定性增加，误差也随之增加。经过上述对比分析证明了本发明所提出的mtf和cnn相结合的方法在不同负载情况下的滚动轴承故障诊断中具有更好的泛化性能。
[0105]
3.5模型在不同数据集规模下的性能分析
[0106]
在实际工况中，故障数据量通常是有限的，当采用小数据集时，由于深度学习网络参数过多，模型可能会出现过拟合，导致模型无法获得较高的识别准确率，因此，在数据集较小时保持良好的诊断效果变得非常重要。在数据集a中，分别考虑不同算法在总数据集的100％、60％、30％和10％四种不同规模数据下进行对比分析，用四种方法对缩减后的故障数据集进行了5次测试，结果如表2所示。
[0107]
表2不同数据集规模下的模型识别准确率
[0108]
[0109][0110]
本发明在四种不同规模数据下的平均识别准确率均高于其他方法，在数据规模为总数据集时，本发明为准确率为99.77％，高于svm的92.66％、bpnn的74.80％、cnn的99.12％；随着数据规模的减小，其他三种方法的识别精度均明显降低，bpnn在数据集大小为总数据集的10％时，准确率下降到了35.71％，是所有方法中最差的，比基于总数据集时降低了39.09％，svm和cnn分别也降低了26.09％和31.12％。相比之下，在数据集规模减小时，本发明在数据集大小为总数据集的10％时识别率依然能达到95.04％，只比100％时低4.73％。这些结果表明当数据集规模减小时，本发明所提方法均比其他方法具有更高的识别准确率。
[0111]
不同规模数据集的诊断性能也可通过箱型图来说明，以比较不同网络识别结果的稳定性，如图8所示。从图8可以看出，本发明所提出的mtf-cnn模型在所有情况下都具有较好诊断效果，并且方差较小，除了在100％和60％数据集下cnn的结果表现出较小方差，其他算法在所有情况下方差都较大，并伴有异常值的出现。结果表明，与几种常用的故障诊断方法相比，本发明能显著提高滚动轴承故障识别效果的准确性和稳定性。
[0112]
以上对本发明实施例所提供的技术方案进行了详细介绍，本发明中应用了具体个例对本发明实施例的原理以及实施方式进行了阐述，以上实施例的说明只适用于帮助理解本发明实施例的原理。