基于多维数据驱动的大型回转支承剩余寿命在线预测方法技术

本发明专利技术公开了一种基于多维数据驱动的大型回转支承剩余寿命在线预测方法，该方法首先对回转支承进行全寿命疲劳试验，获取其整个试验周期的多组振动加速度信号监测数据；然后利用聚类经验模态分解(EEMD)将振动信号分解成多个本征模态分量(IMF)，并利用主成分分析(PCA)选择能反应出回转支承性能退化的若干IMF进行信号重构，以达到滤波消噪的目的；接着，将不同时期的多维加速度信号与试验初期的信号进行PCA，得到一维的连续SPE(C-SPE)，并以此反应出多维加速度振动信号的变化情况，再计算出C-SPE的多个时域特征，建立起回转支承的性能退化模型。本方法人工干预过程少，预测结果更接近工程实际。

【技术实现步骤摘要】
基于多维数据驱动的大型回转支承剩余寿命在线预测方法
本专利技术属于振动信号处理及设备健康监测领域，涉及一种基于多维数据驱动的大型回转支承剩余使用寿命在线预测方法。
技术介绍
大型回转支承作为核心回转连接件常用于风力发电机、挖掘机、起重机等大重型设备中，其严重的故障会导致巨大的损失。因此，实现回转支承在线剩余寿命预测，对制定有效的主动维护或更换策略、提高生产效率有重要的意义。大型回转支承一般工作在背景噪声大，转速低且载荷重的工况下，各部件产生的非平稳随机信号中的故障通过频率可低至1Hz且能量很低，这导致其振动信号的有效信息被完全淹没在背景噪声下，一般的降噪方法很难达到理想效果。此外，大型回转支承直径可达800-5000mm，如此庞大的部件必须采用多个加速度传感器获取不同方向、不同位置的振动信息，才能完整反应出回转支承的运行状况，这就产生了多维的原始信号。如果直接计算出这些信号各自的时域、频域、时频域特征，将带来几十维至上百维的特征向量，使得寿命预测模型建模时耗过长且精度很差。因此，常见的设备性能退化模型建立方法并不适用于大型回转支承。综上所述，基于数据驱动的大型回转支承寿命预测的关键在于振动信号降噪、数据降维以及性能退化模型的建立。聚类经验模态分解(EEMD)非常适用于分析非平稳信号，首先向原始信号中添加均匀白噪声，然后通过包络算法，将一组信号分解成多个本征模态分量(IMF)，不同的IMF代表此信号在同一时间序列上不同的频域范围。主成分分析(PCA)则通过降维处理，将多维数据矩阵投影到主元子空间和残差子空间中，从而用互不相关的几个主元代表所有变量的变化情况。此外，PCA引入的统计量T2和平方预测误差(SPE)可以用于判断不同的多维数据之间的差异，其阈值报警功能可以用于处理多变量多样本的监测数据。最小二乘支持向量机(LS-SVM)作为一种高效的分类、回归工具，近年来被广泛应用在故障分类、寿命预测中。然而，常见的回转支承振动信号故障处理方法中使用EEMD只是将信号进行分解，然后分析各个IMF的能量，很少对IMF进行滤波、选择以重构信号，国外学者提出的EEMD-MSPCA能够有效地进行滤波消噪，但是其效果严重依赖于多个经验参数，技术人员通常无法获得这些参数的最优选择。另一方面，PCA也是主要用于故障诊断，以信号的T2或者SPE是否超过其阈值来判断设备是否出现异常，并没有更深入的研究和应用。因此，实现大型回转支承在线剩余寿命预测，需要充分研究并利用EEMD和PCA的优势，对多维加速度信号进行自适应的降噪、重构、降维、特征提取等，最终获得回转支承性能退化模型，建立起准确的在线剩余寿命预测模型。
技术实现思路
为实现上述目标，本专利技术充分融合利用了EEMD良好的非平衡信号处理能力和PCA对高维数据的降维、差异判断的优良性能，提出了一种基于多维数据驱动的大型回转支承剩余寿命在线预测方法。本专利技术的技术方案是：提出了一种基于多维数据驱动的大型回转支承剩余寿命在线预测方法，该方法首先对回转支承进行全寿命疲劳试验，获取其整个试验周期的多组振动加速度信号监测数据；然后利用聚类经验模态分解(EEMD)将振动信号分解成多个本征模态分量(IMF)，并利用主成分分析(PCA)选择能反应出回转支承性能退化的若干IMF进行信号重构，以达到滤波消噪的目的。接着，将不同时期的多维加速度信号与试验初期的信号进行PCA，得到一维的连续SPE(C-SPE)，并以此反应出多维加速度振动信号的变化情况，再计算出C-SPE的多个时域特征，建立起回转支承的性能退化模型。最后利用LSSVM建立起回转支承剩余寿命与性能退化模型的关系模型，从而建立起回转支承的剩余寿命预测模型。回转支承加装在线监测系统后，获取到的多维振动信号便可按照此方法进行降噪、重构、特征提取，最终实现在线剩余寿命预测，其具体实现步骤如下：1)获取多维加速度信号。对全新的回转支承进行满载荷的疲劳寿命试验，直至其完全失效卡死，使用均布在回转支承定圈内侧的K(K≥4)组加速度传感器采集回转支承内外圈滚道的相对轴向、径向振动信号，保存多维加速度传感器的全寿命试验数据a1,a2,…,ak,…,aK；2)EEMD分解。选择步骤1)中保存的某一组加速度数据ak(k∈[1,K])，按照试验的天数N将其分成N段，然后对每段进行EEMD分解，得到N个IMF矩阵，每个矩阵的列数J等于每段信号分解的IMF的个数，行数I等于此段信号的采样点数；3)IMF初步选择。将第1天的IMF矩阵中的第j(j∈[1,J])阶IMF拆分成多维矩阵M1j，同时将第n(n∈[1,N])天的IMF矩阵中的第j阶IMF拆分成多维矩阵Mnj，将Mnj与M1j进行PCA，以M1j为基础样本得到Mnj的平方预测误差向量SPEnj，求其均值并减去SPE的控制线Scl，得到其差值ERRnj，ERRnj即反应出加速度ak的第n天的第j阶IMF与第一天的同阶IMF之间的差异。需要指出的是，由于白噪声的均匀性，包含大量噪声的高频段IMF在不同时段不会有太大变化，而包含低频有效信息的IMF在性能退化的不同阶段会有较大变化。因此ERRnj越大，说明此阶IMF越能更有效地反应出回转支承性能退化的趋势；重复以上步骤，将j从1取到J，获取第n天各阶IMF与第1天对应IMF的差异值ERRnj,如前所述，ERRnj越大的IMF越能反应出回转支承性能退化趋势，因此更应该作为当天加速度信号重构的IMF。为量化这一过程，设权值为α(α∈(0,1))，当ERRnj最大的h(h∈[1,J])个IMF累加之和占所有ERRnj累加之和的比例超过α时，此h个IMF即作为第n天加速度ak的重构元素；4)IMF最终确定。由于回转支承性能退化过程中的随机性，如：不同的阶段不同的部件(内圈、外圈或滚珠)产生故障时会产生不同频率的冲击，这使得不同天数里面重构加速度ak所需要的IMF可能是不同的。为统一加速度ak重构时使用的IMF，将n从1取到N，重复步骤3)，得到所有N天内重构加速度ak所需要的IMF序列，这些IMF序列中会有若干重复出现的，也会有个别偶尔出现的，将出现次数之和占所有IMF序列个数的比例超过α的H(H∈[1,J])个IMF作为回转支承全寿命试验过程中加速度ak信号重构的依据；5)信号降噪与重构。从每天的IMF矩阵中选择步骤4)中确定的若干个IMF进行累加，便可得到当天的ak重构信号。需要指出的是，此处重构不仅是选择了能反应回转支承性能退化过程的IMF，更是将背景噪声所在的高频IMF舍去，间接达到了消噪的效果；6)多维加速度信号重构。将k从1取到K，重复步骤2)-5)，即可获得多维加速度全寿命周期信号的重构信号；7)数据降维与特征提取。将每天重构后的a1,a2,…,aK信号组成矩阵，即可得到N个K维的加速度信号矩阵，然后将第n(n∈[1,N])天的矩阵与第1天的矩阵进行PCA，便可以以一维向量SPE反应出K维加速度矩阵相对第1天的变化趋势，既保留了多维数据的特征信息，又达到了降维的目的。将n从1到N计算出的SPE连接起来，便得到回转支承整个寿命周期中的性能退化特征：连续SPE(C-SPE)；8)建立回转支承性能退化模型。步骤7)中建立的C-SPE是不具有频域特性的，因此计算其1本文档来自技高网...

【技术保护点】
一种基于多维数据驱动的大型回转支承剩余寿命在线预测方法，其特征在于，包括获取全寿命疲劳试验振动信号，基于EEMD‑PCA的信号分解、消噪及重构，多维振动信号的降维及C‑SPE特征的提取，基于C‑SPE及时域特征建立回转支承性能退化模型，基于LSSVM剩余寿命预测模型的建立与训练，以及回转支承在线剩余寿命预测方法，其主要步骤如下：(1)获取多维加速度信号：对全新的回转支承进行满载荷的疲劳寿命试验，直至其完全失效卡死，使用K(K≥4)组加速度传感器采集回转支承内外圈滚道的相对轴向、径向振动信号，保存多维加速度传感器的全寿命试验数据；(2)EEMD分解：选择一组加速度数据ak(k∈[1,K])，按照试验的天数N将其分成N段，然后对每段进行EEMD分解，得到N个IMF矩阵，此矩阵的列数J等于分解的IMF的个数，行数I等于此段信号的采样点数；(3)IMF初步选择：将第1天的IMF矩阵中的第j(j∈[1,J])阶IMF拆分成多维矩阵M1j，同时将第n(n∈[1,N])天的IMF矩阵中的第j阶IMF拆分成多维矩阵Mnj，将Mnj与M1j进行PCA，以M1j为基础样本得到Mnj的平方预测误差向量SPEnj，求其均值并减去SPE的控制线Scl，得到其差值ERRnj，ERRnj即反应出加速度ak的第n天的第j阶IMF与第一天的同阶IMF之间的差异；需要指出的是，由于白噪声的均匀性，包含大量噪声的高频段IMF在不同时段不会有太大变化，因此其ERR较小，而包含低频有效信息的IMF在性能退化的不同阶段会有较大变化；因此ERRnj越大，说明此阶IMF越能更好地反应出回转支承性能退化的趋势；重复以上步骤，将j从1取到J，获取第n天各阶IMF与第1天对应IMF的差异值ERRnj,如前所述，ERRnj越大的IMF越能反应出回转支承性能退化趋势，因此更应该作为当天加速度信号重构的IMF。为量化这一过程，设权值为(α∈(0,1))，当ERRnj最大的h(h∈[1,J])个IMF累加之和占所有ERRnj累加之和的比例超过α时，此h个IMF即作为第n天加速度ak的重构元素；(4)IMF最终确定：由于回转支承性能退化过程中的随机性，如：不同的阶段不同的部件(内圈、外圈或滚珠)产生故障时会产生不同频率的冲击，这使得不同天数里面重构加速度ak所需要的IMF可能是不同的。为统一加速度ak重构时使用的IMF，将n从1取到N，重复步骤(3)，得到所有N天内重构加速度ak所需要的IMF序列，这些IMF序列中会有若干重复出现的，也会有个别偶尔出现的，将出现次数之和占所有IMF序列个数的比例超过α的H(H∈[1,J])个IMF作为回转支承全寿命试验过程中加速度ak信号重构的依据；(5)信号降噪与重构：从每天的IMF矩阵中选择步骤(4)中确定的IMF进行累加，便可得到当天的ak重构信号；需要指出的是，此处重构不仅是选择了能反应回转支承性能退化过程的IMF，更是将背景噪声所在的高频IMF舍去，间接达到了消噪的效果；(6)多维加速度信号重构：将k从1取到K，重复步骤(2)‑(5)，即可获得多维加速度全寿命周期信号的重构信号；(7)数据降维与特征提取：将每天重构后的a1,a2,…,aK信号组成矩阵，即可得到N个K维的加速度信号矩阵，然后将第n(n∈[1,N])天的矩阵与第1天的矩阵进行PCA，便可以以一维向量SPE反应出K维加速度矩阵相对第1天的变化趋势，既保留了多维数据的特征信息，又达到了降维的目的；将n从1到N计算出的SPE连接起来，便得到回转支承整个寿命周期中的性能退化特征：连续SPE(C‑SPE)；(8)建立回转支承性能退化模型：步骤(7)中建立的C‑SPE是不具有频域特性的，因此计算其10个时域指标，并连同C‑SPE本身建立起回转支承的性能退化模型；(9)建立回转支承剩余寿命预测模型：若试验中回转支承的转速是v(r/min)，则当试验进行到第n(n∈[1,N])天时，回转支承已经运行了n×24×60×v转，其剩余寿命为(N‑n)×24×60×v转，而第n天的C‑SPE及其时域特征已经从步骤(1)‑(8)中获得，因此将C‑SPE及其时域特征作为输入，将剩余寿命作为输出，便可利用LS‑SVM建立起回转支承剩余寿命预测模型；(10)回转支承在线剩余寿命预测：回转支承剩余寿命预测模型建立后，同型号的回转支承在实际使用时，可将其多维加速度信号进行EEMD分解，然后按照步骤(4)确定的IMF序列进行信号重构，再按照步骤(7)得到当天的C‑SPE，计算其时域特征并代入到步骤(9)建立的模型中，便可在线预测回转支承的剩余使用寿命。...

【技术特征摘要】
1.一种基于多维数据驱动的大型回转支承剩余寿命在线预测方法，其特征在于，包括获取全寿命疲劳试验振动信号，基于EEMD-PCA的信号分解、消噪及重构，多维振动信号的降维及C-SPE特征的提取，基于C-SPE及时域特征建立回转支承性能退化模型，基于LSSVM剩余寿命预测模型的建立与训练，以及回转支承在线剩余寿命预测方法，其主要步骤如下：(1)获取多维加速度信号：对全新的回转支承进行满载荷的疲劳寿命试验，直至其完全失效卡死，使用K(K≥4)组加速度传感器采集回转支承内外圈滚道的相对轴向、径向振动信号，保存多维加速度传感器的全寿命试验数据；(2)EEMD分解：选择一组加速度数据ak(k∈[1,K])，按照试验的天数N将其分成N段，然后对每段进行EEMD分解，得到N个IMF矩阵，此矩阵的列数J等于分解的IMF的个数，行数I等于此段信号的采样点数；(3)IMF初步选择：将第1天的IMF矩阵中的第j(j∈[1,J])阶IMF拆分成多维矩阵M1j，同时将第n(n∈[1,N])天的IMF矩阵中的第j阶IMF拆分成多维矩阵Mnj，将Mnj与M1j进行PCA，以M1j为基础样本得到Mnj的平方预测误差向量SPEnj，求其均值并减去SPE的控制线Scl，得到其差值ERRnj，ERRnj即反应出加速度ak的第n天的第j阶IMF与第一天的同阶IMF之间的差异；需要指出的是，由于白噪声的均匀性，包含大量噪声的高频段IMF在不同时段不会有太大变化，因此其ERR较小，而包含低频有效信息的IMF在性能退化的不同阶段会有较大变化；因此ERRnj越大，说明此阶IMF越能更好地反应出回转支承性能退化的趋势；重复以上步骤，将j从1取到J，获取第n天各阶IMF与第1天对应IMF的差异值ERRnj,如前所述，ERRnj越大的IMF越能反应出回转支承性能退化趋势，因此更应该作为当天加速度信号重构的IMF;为量化这一过程，设权值为(α∈(0,1))，当ERRnj最大的h(h∈[1,J])个IMF累加之和占所有ERRnj累加之和的比例超过α时，此h个IMF即作为第n天加速度ak的重构元素；(4)IMF最终确定：由于回转支承性能退化过程中的随机性，如：不同的阶段不同的部件(内圈、外圈或滚珠)产生故障时会产生不同频率的冲击，这使得不同天数里面重构加速度ak所需要的IMF可能是不同的;为统一加速度ak重构时使用的IMF，将n从1取到N，重复步骤(3)，得到所有N天内重构加速度ak所需要的IMF序列，这些IMF序列中会有若干重复出现的，也会有个别...

【专利技术属性】
技术研发人员：黄筱调，封杨，陈捷，王华，洪荣晶，
申请(专利权)人：南京工业大学，南京工大数控科技有限公司，
类型：发明
国别省市：江苏;32