一种基于改进SSA优化VMD和CNN参数的滚动轴承故障诊断方法

一种基于改进ssa优化vmd和cnn参数的滚动轴承故障诊断方法
技术领域
1.本发明属于机械故障诊断领域，特别是涉及一种基于改进ssa优化vmd和cnn参数的滚动轴承故障诊断方法。


背景技术：

2.制造业是立国之本、强国之基。随着现代科技工业技术尤其是信息技术的迅速发展，航空、通信、工业应用等各个领域的工程系统日趋复杂，系统的综合化、智能化程度不断提高，组成环节和影响因素的不断增加，发生故障和功能失效的几率也逐渐加大，因此，故障诊断逐渐成为研究者关注的焦点。
3.关于故障诊断方法的分类，国内外不同的研究机构和组织的提法不尽一致，从目前主流的技术和应用研究工作来看，主要可以分为：基于模型的故障诊断技术；基于统计可靠性的故障诊断技术及基于数据驱动的故障诊断技术。
4.由于数据驱动诊断技术是对数据直接进行处理，对人为因素和专家知识依赖较少，具有普遍适用性且能达到较高的诊断精度，因此现有的故障诊断方法研究中，大多采用此方法。其中，dragomirestkiy提出了一种完全非递归，且具有完备数学理论支撑的变分模态分解方法(vmd)。由于vmd自身不存在端点效应，而且可以自适应的将复合信号中两个频率相似的谐波信号分隔开，因此该方法已经广泛应用于机械故障诊断中。
5.传统的vmd算法受参数人为设置的限制，具有很大的主观性，导致故障诊断率不高。采用vmd进行分解信号过程中，需要预先设定分解个数k和惩罚因子α，这两个参数直接影响到vmd的分解效果，当k值取值较小时，会出现“欠分解”现象，若α取值较小，会出现模态混叠现象，反之，当k值取值较大时，会出现“过分解”现象，若α取值较大，会出现虚假分量现象。因此对vmd参数进行优化，以确定最优的k和α是亟需解决的问题。
6.20世纪八九十年代，由于计算机水平和其他技术的限制，深度学习并没有表现出优异的性能。2006年hinton提出一种贪婪逐层预训练的策略，为深度网络的训练提供了新思路，并被证实可以用来训练许多其他类型的深度网络。自此以后，深度学习得到广泛关注，并在模式识别、图像处理等领域取得了良好的发展成果。目前常用的深度学习模型有深度信念网络(deep belief network，dbn)，卷积神经网络(convolutional neuralnetworks，cnn)，循环神经网络(recurrent neural network，rnn)，稀疏自动编码器(sparse auto-encoders，sae)等。
7.卷积神经网络是一种受生物的视知觉机制启发而构建的深层学习模型，擅长处理图像数据。在故障诊断领域，目前主要的研究方向有两大类，一类是通过一定的信号处理方法得到将机械设备振动信号的时频特征图，将cnn作为分类器使用；另一类是将未处理的原始信号直接作为cnn的输入，利用cnn强大的非线性特征自学习能力，自适应逐层提取输入数据特征，并进行分类任务。
8.cnn主要分为滤波层和分类层。滤波层由卷积层、池化层和激活层构成，而分类层
大多直接由全连接形成。cnn的层级分布，决定了cnn具有较多的超参数，所谓超参数优化或调整问题就是在网络模型中选取一个最优化超参数搭配的问题。对于给定的问题，很难了解超参数的最佳组合方案，但用以往的经验法复制在其他问题上使用的值，或者通过反复试验寻找最佳值都会耗费大量的实验成本和时间，而且可能得到的并不是最优解。因此对cnn超参数进行优化，以确定分类效果最优参数是亟需解决的问题。
9.随着现代工业技术的快速发展，设备自身的结构和运行环境逐渐复杂，传感器采集到的设备状态监测数据也呈现多样化趋势，这使得故障特征的难度不断增加。因此，如何充分利用好vmd完成对信号降噪，使cnn更好的实现故障特征的提取与分类，是今后机械设备故障诊断方法的研究重点。


技术实现要素：

10.针对上述现有技术的不足，本发明提供了一种基于改进ssa优化vmd和cnn参数的滚动轴承故障诊断方法。本发明通过改进麻雀搜索算法优化vmd和cnn参数，使得vmd得到最优参数组合[k0,α0]，以此为基础进行vmd分解、重构信号，能够得到更加准确的信号特征。同样的，以验证集分类错误率为适应度函数求其最小值，目的是找到一组超参数，使得网络错误率最低，使得cnn得到最优超参数，从而完成滚动轴承的故障诊断。
[0011]
为了实现上述目的，本发明采用如下技术予以实现：
[0012]
一种优化vmd和cnn参数的故障诊断方法，包括如下过程：
[0013]
s1、获取滚动轴承振动信号；
[0014]
s2、基于改进麻雀搜索算法优化vmd算法分解个数和惩罚因子，得到vmd的分解个数和惩罚因子的最优参数组合；
[0015]
s3、利用改进参数后的vmd将振动信号分解成若干含故障信息的信号分量；
[0016]
s4、基于改进的麻雀搜索算法优化cnn的超参数，包括学习率、训练次数、各个卷积层的卷积核大小与卷积核数量和全连接层的节点数；
[0017]
s5、输入s3的信号分量，利用改进参数后的cnn实现故障特征的提取和故障分类。
[0018]
其中，在s2所述的改进麻雀搜索算法中，采用猫映射产生初始种群，其表达式为：式中:a1、b1为任意实数；mod1表示求a1小数部分。利用反向选择策略生成初始反向解加入初始种群，利用精英反向学习策略对麻雀种群进行选择，以此提高麻雀种群的质量，使种群增加多样性，为较好的全局搜索奠定基础，同时可以利用反向解引导算法跳出局部最优，并且通过动态边界不断缩小搜索空间，有利于提高算法收敛速度。利用改进的tent混沌和柯西变异扰动策略加入算法，避免算法陷入局部最优，提高了全局搜索能力和寻优精度。改进后tent混沌映射表达式为：其中n为序列内粒子数。柯西变异公式为：mutation(x)＝x(1 tan(π(u-0.5)))，式中：x为原来个体位置，u为(0,1)区间内的随机数。改进探索
者位置更新公式，令其中，α是区间(0,1]的随机数，iter
max
为迭代次数最大值。使探索者跟随者自适应调整策略，其公式为：其中，b为比例系数，k为扰动偏离因子，对非线性递减的r值进行扰动。将vmd分解完的信号的包络熵值作为改进的麻雀算法的适应度函数，计算vmd最优参数组合[k0,α0]。本发明的特点及有益效果：
[0019]
1、本发明改进了麻雀搜索算法，在寻优参数的收敛速度与精度上都能取得良好的优化效果。
[0020]
2、本发明基于改进的麻雀搜索算法优化了vmd的参数，相比于传统的vmd方法能筛选出含故障信号更多的信号分量，实现对故障特征信息的准确提取。
[0021]
3、本发明基于改进的麻雀搜索算法优化了cnn的超参数，能够使优化参数后的cnn更有效的提取轴承故障特征，更有效的判断轴承故障类别，应用于轴承信号分析中，具有较为广泛的实用性。
附图说明
[0022]
图1为本发明流程图
[0023]
图2为改进vmd参数的流程图
[0024]
图3为改进cnn超参数的流程图
[0025]
图4为改进vmd参数后分解出信号分量波形图
[0026]
图5为改进vmd参数后分解出信号分量频谱图
[0027]
图6为改进cnn超参数后的网络准确率与损失率图
[0028]
图7为未改进cnn超参数的网络准确率与损失率图
具体实施方式
[0029]
本发明是一种基于改进ssa优化vmd和cnn参数的滚动轴承故障诊断方法，其流程图如图2和图3所示。一种优化vmd和cnn参数的故障诊断方法过程包括如下步骤：
[0030]
(1)优化vmd过程包括如下步骤
[0031]
步骤1，建立适应度为包络熵和峭度指标的综合函数模型，初始化改进的麻雀算法参数，种群数量，最大迭代次数，参数上下限，参数维度，以及输入参数的取值范围；
[0032]
步骤2，猫映射初始化种群，再通过精英反向解选出优质种群作为初始种群；
[0033]
步骤3，在vmd算法参数的取值范围内初始化麻雀的位置坐标，然后计算每一个麻雀的适应度值和麻雀种群的平均适应度值；
[0034]
步骤4，如果个体适应度值小于群体平均适应度值，则表明出现“聚集”趋势，进行柯西变异，如果比变异之前的个体更优，则用变异后的个体替代变异前的个体，否则保持原个体不变。如果个体适应度值大于等于群体平均适应度值则表明出现“发散”趋势，对个体进行tent混沌扰动，如果扰动后的个体性能更优，则用扰动后的代替之前个体，否则保持不变；
[0035]
步骤5，根据麻雀种群当前的状态，更新整个种群所经历的最优位置及其适应度，最差位置及其适应度；
[0036]
步骤6，判断算法运行是否达到最大迭代次数或者求解精度，若是，循环结束，输出寻优结果并得到分解后的信号分量如图3所示，其频率如图4所示，得到vmd的分解个数和惩罚因子的最优参数组合，否则返回步骤4。
[0037]
其中，步骤1包络熵作为适应度函数，是一种很好的评价信号稀疏特性的标准，熵值大小反映了概率分布的均匀性轴承早期故障信号经vmd算法处理后，如果所得imf分量中包含的噪声较多，与故障相关的周期性冲击特征不明显，则分量信号的稀疏性较弱，包络熵值较大。如果imf分量中包含的故障特征信息较多，波形中出现规律性冲击脉冲，则信号将呈现出较强的稀疏特性，包络熵值较小。包络熵公式为：公式中：pj是a(j)的归一化形式；a(j)是信号x(j)经希尔伯特变化解调后的包络信号。峭度指标是一种无量纲指标，对故障信号微弱变化敏感，非常适用于轴承的早期故障诊断，分量的峭度定义为式中:μi为imf的均值；σi为imf的方差。当轴承无故障时，振动信号接近正态分布，峭度指标越等于3。当滚动轴承发生局部故障时，振动信号在故障冲击作用下将偏离正态分布，而轴承的故障信息也往往包含在由故障冲击成分引起的幅值调制信号中。因此，峭度值越大冲击成分与故障信息也就越丰富，故障特征更易于提取。因此，将两种指标进行融合可以更好的有利于对振动信号中的故障特征进行提取。建立综合目标函数c可表示为:c＝ei 1/abs(qd
i-3)。
[0038]
步骤2中精英反向学习策略是针对反向学习不一定比当前搜索空间更易搜索到全局最优解这一问题提出的，该策略利用优势个体来构造反向种群，以此增强种群多样性。具体是设xi＝(x
i1
,x
i2
,
…
,x
id
)是一普通粒子，对应的自身极值即为精英粒子则精英反向解可定义为：其中：μ∈(0,1)为服从正态分布的随机数；[daj,dbj]为第j维搜索空间的动态边界，可按此式得到：利用动态边界来代替搜索空间的固定边界，能够积累搜索经验，使生成的反向解位于逐步缩小的搜索空间，加速算法收敛。当生成的反向解位于边界之外时，使用随机生成的方法进行重置。如此式所示：
[0039]
(2)优化cnn超参数过程包括如下步骤
[0040]
步骤1，初始化参数，确定麻雀种群数、最大迭代次数以及边界等基本参数，确定cnn的层级结构；
[0041]
步骤2，输入已经通过vmd处理好的数据，并使训练集，验证集，测试集按照7：2：1的比例分配；
[0042]
步骤3，计算此时麻雀个体的适应度值且同时排序，适应度函数为验证集的误差率；
[0043]
步骤4，更新麻雀个体位置包括发现者、跟随者、警戒者的位置；
[0044]
步骤5，是否满足适应度条件，若是则输出最优超参数，否则返回步骤3；
[0045]
步骤6，利用优化参数后的cnn训练模型，其准确率与损失率如图5所示；
[0046]
步骤7，将测试集输入训练好的模型，完成故障诊断。
[0047]
为验证所提出的方法在滚动轴承故障诊断的可行性，本专利采用凯斯西储大学的滚动轴承数据集，通过matlab软件对数据进行实验及测试。滚动轴承故障数据标书及样本分布如表1所示。表1滚动轴承样本分布表
[0048]
实验样本选用48khz的轴承故障数据，分别选取7mil、14mil和21mil故障体量级别下，三种包含轴承内圈、外圈及滚动体的故障位置，外加一组正常数据共10组，按7∶2∶1的比例划分训练集、验证集及测试集(其中训练集和测试集无任何交集)。分别对训练集、验证集和测试集进行处理。利用优化参数后的cnn训练模型，其准确率与损失率如图5所示与其未改进前的损失率和准确率如图6对比有显著提升，因此本专利采用的改进麻雀搜索算法优化vmd和cnn参数的滚动轴承故障诊断方法有效，能够提高轴承故障诊断精度。