一种基于DBN-SVDD模型的滚动轴承性能退化评估方法与流程

一种基于dbn-svdd模型的滚动轴承性能退化评估方法
技术领域
1.本发明涉及滚动轴承性能评估领域，尤其是涉及一种基于dbn-svdd模型的滚动轴承性能退化评估方法。


背景技术：

2.滚动轴承设备是机械系统中的重要部件之一，其运行状态直接影响到机械系统的安全性，而轴承状态从正常运行到最终失效通常要经历不同的性能退化过程。因此，对轴承设备进行可靠的性能退化评估，准确追踪其运行状态，对机械设备的维护维修具有重要意义。
3.在滚动轴承性能退化评估中，最为关键的是退化特征的提取和退化特征指标的评估。常用的提取退化特征指标的方法有时、频域特征和信号分解方法。时域特征指标如有效值(rms)和峰度。rai等人采用基于经验模式分解和单值分解的模型提取轴承的初步退化趋势，然后采用基于k-medoids聚类模型构建最终退化特征指标，评价轴承的健康状态。周建民等人提出一种融合fcm-svdd模型的方法，利用自回归模型(ar)对轴承全寿命周期数据进行特征提取，将得到的退化指标作为特征矩阵输入到fcm和svdd模型,得到融合方法的性能退化曲线。但是这些退化特征指标构建模型高度依赖于手工经验，且在滚动轴承振动信号取自不同的操作系统时，要自适应地改变时频域指标和指标评估的方法，才能对轴承设备做出可靠的性能退化评估。
4.深度学习是近年来关注的热点之一，其在面对大数据时具有较强的自适应非线性特征提取能力，可以减少对手工经验的依赖，而dbn作为深度学习的基本模型之一被广泛应用于滚动轴承故障分类诊断，但用于轴承性能退化评估上的并不多见，dong等人提出了一种基于迁移学习和dbn的滚动轴承退化评估模型，将退化评价的研究转化为退化模式的分类任务，完成不同工况下的退化模式评估。peng等人提出使用不带输出层的dbn轴承信号特征得到性能退化指标做出性能退化评估，并使用改进的粒子滤波模型来预测滚动轴承的剩余寿命。
5.但是上述提到的基于dbn提取特征的滚动轴承性能退化评估方法，在建立dbn模型时，所用信号数据多为试验所得的轴承全寿命周期数据，而轴承实际运行中正常状态居多，获得的大数据也多为正常状态下的振动信号。
6.基于此，应用正常数据建模，并能评估整个滚动轴承寿命周期的性能退化情况显得尤为重要。本文提出一种基于dbn-svdd的滚动轴承性能退化评估方法，使用轴承正常数据对dbn进行建模实现特征向量的提取，提取的正常状态下的特征向量作为训练样本建立svdd评估模型，用轴承全寿命周期数据测试训练后的dbn-svdd模型，完成其性能退化评估。


技术实现要素：

7.为了克服现有技术的不足，本发明提供了一种基于dbn-svdd模型的滚动轴承性能退化评估方法。
8.为了实现以上目的，本发明提供如下技术方案：
9.一种基于dbn-svdd的滚动轴承性能退化评估方法，包括步骤：
10.(1)提取特征：获取正常状态下的滚动轴承振动信号数据，对所获取的振动信号进行fft和归一化处理，将所得到的数据作为训练样本建立dbn特征提取模型，获得特征向量；
11.(2)建立模型：用dbn模型中dbn无监督预训练部分提取出的正常状态下的特征向量作为训练样本建立svdd评估模型，并获得滚动轴承的性能退化指标；
12.(3)实时评估：把将同一滚动轴承的未知振动信号输入上述训练建立完成的模型中，得到性能退化指标，描绘出滚动轴承的性能退化曲线。
13.本技术方案进一步的优化，所述dbn是由多个受限玻尔兹曼机(rbm)逐个堆叠而形成，每个rbm都由个可见层和个隐含层构成，每个可见层由多个可见层单元构成，每个隐含层也由多个隐含层单元构成，所述可见层为每个rbm的输入层，所述隐含层为每个rbm的输出层，每个rbm的输出层作为下一个rbm的输入层。
14.本技术方案进一步的优化，步骤(1)具体包括：
15.(a)将经过fft处理后的滚动轴承振动信号幅值范围规范到0～1，归一化方式为：其中y
min
和y
max
表示归一化的范围值，这里y
min
＝0，y
min
＝0。x
min
和x
max
表示被归一化数据中的最小值和最大值。
16.(b)将经过fft和归一化预处理后的滚动轴承振动信号的数据，输入dbn模型中，完成对dbn模型中无监督预训练部分的训练；
17.(c)将滚动轴承数据输入训练好的dbn模型中，利用dbn模型中堆叠的多个rbm对输入的滚动轴承数据进行逐层特征提取，最后一个受限玻尔兹曼机提取出dbn提取的特征向量。
18.本技术方案进一步的优化，svdd评估模型建立的具体内容为：
19.将通过dbn模型提取的特征向量作为svdd模型的训练样本，完成svdd模型训练，构建包含正常状态下特征向量的半径为r的超球体，将新样本输入训练好的svdd模型中，得出新样本与超球体球心之间的相对距离d，将新样本与超球体球心之间的相对距离d作为轴承的性能退化指标di，当d1≤r，新样本被识别为轴承在正常状态下运行的目标，当d1≥r，表明轴承处于退化状态，且di可以表示轴承的故障严重程度，即di越大故障严重程度越大。
20.与现有技术相比本发明具有以下有益效果：对振动信号进行fft变换和归一化处理，再输入dbn特征提取模型中，自动提取出特征向量，减少了对人工经验的依赖；并且利用svdd模型对特征向量的描述可以较好的表征轴承性能退化情况。以正常状态下的振动信号作为训练样本进行dbn-svdd性能退化评估模型的构建，克服实际轴承设备运行中故障样本较难获取的问题，对实际生产中轴承早期故障的发生和性能退化评估具有很好的指导意义，对于设备性能监测具有重要意义。
附图说明
21.图1为本发明滚动轴承性能退化评估流程图；
22.图2为本发明dbn结构图；
23.图3为本发明svdd原理图；
24.图4为本发明具体实施例所用数据原始信号图；
25.图5为本发明具体实施例第四层rbm提取的特征向量图；
26.图6为本发明具体实施例基于dbn-svdd的性能退化评估结果图；
27.图7为本发明具体实施例基于dbn-svdd的性能退化评估结果局部图；
28.图8为本发明实施例的对比方法有效值做性能退化指标结果图。
具体实施方式
29.为详细说明技术方案的技术内容、构造特征、所实现目的及效果，以下结合具体实施例并配合附图详予说明。
30.参阅图1所示，为本发明滚动轴承性能退化评估流程图，本发明为一种基于dbn-svdd的滚动轴承性能退化评估方法，该方法包括以下步骤：首先提取特征向量，的获取正常状态下的滚动轴承振动信号数据，对所获取的振动信号进行fft和归一化处理，将所得到的数据作为训练样本建立dbn特征提取模型，获得特征向量。然后建立dbn模型，用dbn模型中dbn无监督预训练部分提取出的正常状态下的特征向量作为训练样本建立svdd评估模型，并获得滚动轴承的性能退化指标。最后实时评估，把将同一滚动轴承的未知振动信号输入上述训练建立完成的模型中，得到性能退化指标，描绘出滚动轴承的性能退化曲线。
31.参阅图2所示，为本发明dbn结构图，dbn是由多个受限玻尔兹曼机rbm单元组成的多隐层神经网络，dbn的结构如图2所示。其中的rbm由两层组成，用于输入数据的可视化层v和提取特征的隐含层h，可视化层v和隐含层h中分别包含几个独立不相互连接的神经元节点，但层层之间是完全连接的。
32.例如，一个rbm单元中，可视化层v有m个可见单元，隐含层h有n个隐藏单元，a和b分别表示可视化层v和隐含层h的偏置项，w是可视化层v和隐含层h之间的连接权重，则该联合结构网络的层间单元节点具有的能量表示为：其中w
ij
表示第i个节点vi和第j个节点hj之间的权重值。
33.rbm中可见层v与隐含层h之间的联合概率分布可由网络的节点能量表示为：
[0034][0035]
rbm中可见层v与隐含层h之间的条件概率分布可以表示为：
[0036][0037][0038][0039]
其中，f(x)为激活函数，采用sigmoid函数。
[0040]
dbn模型预训练时，当前的rbm的输出被用作下一个rbm的输入，每个rbm通过使用对比发散(cd-k)进行无监督学习单独训练，以此更新网络参数，降低重构误差，该算法更新
准则为：
[0041][0042][0043][0044]
其中，k为迭代次数，η为学习率。
[0045]
构成dbn的多个rbm对轴承正常状态数据进行无监督贪婪学习来确定每层rbm之间的连接权重值及可见层和隐藏层的偏置项，以此完成对dbn特征提取模型的训练。
[0046]
步骤(2)建立模型：用dbn模型中dbn无监督预训练部分提取出的正常状态下的特征向量作为训练样本建立svdd评估模型，并获得滚动轴承的性能退化指标。
[0047]
参阅3所示，本发明svdd原理图，其主要思想是找到一个包含所有或大部分目标的最小体积的最优超球，针对一个目标类样本集{xi，i＝1，2，...，n}，寻找包含所有或大多数目标类样本的最优超球，该超球由中心c和半径r描述，满足以下优化函数：
[0048][0049]
其中c是一个惩罚参数，它控制超球体和误差之间的权衡，ξi是松弛变量，允许一些训练数据在超球体之外。
[0050]
一般来说，方程(9)是通过引入拉格朗日乘子来求解的，它可以引入拉格朗日乘子αi转化为以下函数：
[0051][0052]
由于输入空间中的数据并不总是线性预测的，引入核函数k把原来线性不可分的问题转换成高维空间中线性可分的问题，同时并没有增加太多的计算量，高斯核是最常用的函数，定义如下：
[0053][0054]
其中σ是宽度参数，将其用于方程(10)得到：
[0055][0056]
当样本点位于超球体内，则αi＝0；当样本点位于超球体的边界上，则0＜αi＜c；当样本点位于超球体之外时，则αi＝c。在超球体界面上的样本为支持向量(x
sv
)。通过核函数的方法将训练样本由低维空间映射到高维空间中，从而使得新的高维空间中对目标样本进行很好的描述，在映射空间中超球体半径r可由任意支持向量到中心的距离计算得到：
[0057][0058]
对于新样本z，它与球心的距离可表示为：
[0059][0060]
求出新样本与超球体球心之间的相对距离d，将其作为轴承的性能退化指标di，如果di≤r，则新样本被识别为轴承在正常状态下运行的目标，否则它是一个离群点，表明轴承处于退化状态。
[0061]
本发明利用在ablt系列轴承寿命强化试验机上实验获取的一组全寿命周期滚动轴承信号对本发明进行验证，其中轴承型号为6307，转速为3000r/min，径向载荷为11.13kn，采样频率为25.6hz，采用加速度传感器每隔1min采集一组长度为20480的数据，。数据采集从轴承正常状态一直持续到寿命结束，共采集了10620组数据，不包括磨合期数据(试验设备开机启动预热，逐步增加载荷，数据波动不稳定的初始阶段)，最终失效形式均为内圈严重点蚀。
[0062]
参阅图4所示，为本发明具体实施例所用数据原始信号图，由于现有的滚动轴承全寿命周期数据划分长度10620过长，对计算机要求过高，分别将每组的轴承数据再次分割为十组，长度为2048，因此进行实验验证轴承全寿命周期数据10620组长度为2048的数据。
[0063]
本次实验选取前3000组正常状态下的轴承数据作为训练样本经实验确定dbn特征提取模型设置为四层，即其包含四个rbm，结构为500-200-100-10，迭代次数为50，初始学习率为0.01。svdd评估模型里的惩罚参数设置为0.5，核宽度设置为1。将10620组轴承全寿命周期数据输入上述训练好的dbn-svdd模型，其中经过dbn特征提取模型的数据大小为10620*10，将其按照顺序十组合为一组变为1062*100，再输入svdd评估模型，得出该滚动轴承的性能退化评估结果。
[0064]
参阅图5、图6和图7所示，分别为本发明具体实施例第四层rbm提取的特征向量图、本发明具体实施例基于dbn-svdd的性能退化评估结果图和本发明具体实施例基于dbn-svdd的性能退化评估结果局部图。从中可以看出，当样本处于0～513min期间，由于其di值在r值之下，此阶段处于正常运行状态；在513～978min期间的样本的di值超出r值并在一定范围内保持稳定，说明此阶段处于早期微弱故障阶段；在978～1309min段的di值突然增加，且在小范围内di值保持稳定，此阶段轴承处于故障加剧第一阶段；从1309min开始的di值又有所上升，此时轴承进入故障加剧第二阶段。进一步证明本发明的优越性，选取常用性能退化指标有效值(rms)对该组轴承数据进行评估。
[0065]
参阅图8所示，为本发明实施例的对比方法有效值做性能退化指标结果图。该组轴承数据的rms在0～980min期间波动趋势基本保持稳定，即rms指标未能检测出即轴承早期微弱故障发生的时刻点；在980～1040min期间，rms值突增，且波动稳定，可认为此阶段为轴承故障加剧第一阶段；从1040min开始，轴承rms值较于前一阶段大幅增加且波动混乱，此刻轴承进入故障加剧第二阶段。与本发明发现轴承故障加剧第一阶段的时刻点978min相比晚了2min，与本发明发现故障加剧第二阶段的时刻1039min晚了1min。
[0066]
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存
在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者终端设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者终端设备所固有的要素。在没有更多限制的情况下，由语句“包括
……”
或“包含
……”
限定的要素，并不排除在包括要素的过程、方法、物品或者终端设备中还存在另外的要素。此外，在本文中，“大于”、“小于”、“超过”等理解为不包括本数；“以上”、“以下”、“以内”等理解为包括本数。
[0067]
尽管已经对上述各实施例进行了描述，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例做出另外的变更和修改，所以以上仅为本发明的实施例，并非因此限制本发明的专利保护范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围之内。