基于偏最小二乘法极限学习机的滚动轴承故障预测方法技术

本发明专利技术涉及偏最小二乘法极限学习机的滚动轴承故障预测方法。首先分析时域、频域、时频域等多个特征指标，提出基于半正态分布和经验小波降噪相结合的特征提取方法，对滚动轴承进行故障诊断，贴近原始信号，达到更好的降噪效果。又针对多特征参数对滚动轴承故障衰退特征综合评价，提出残差改进的ISOMAP非线性特征降维方法与模糊C均值结合的方法，提高滚动轴承在不同衰退阶段的变化趋势和分类精度。以极限学习机理论为基础提出基于偏最小二乘法极限学习机的数据预测模型，优化ELM中参数，选取最优隐含层节点数和连接层权值，并选取Softmax激活函数。这样，预测精度高，计算用时短，聚类后的特征值预测效果好。通过以上步骤进行滚动轴承的故障阶段精确预测。

Fault Prediction Method of Rolling Bearing Based on Partial Least Squares Limit Learning Machine

The invention relates to a rolling bearing fault prediction method of a partial least squares extreme learning machine. Firstly, several characteristic indexes such as time domain, frequency domain and time-frequency domain are analyzed, and a feature extraction method based on the combination of semi-normal distribution and empirical wavelet denoising is proposed. The fault diagnosis of rolling bearing is carried out, which is close to the original signal and achieves better denoising effect. Aiming at the comprehensive evaluation of rolling bearing failure decay characteristics by multi-feature parameters, a method combining residual improved ISOMAP non-linear feature dimension reduction method with fuzzy C-means is proposed to improve the change trend and classification accuracy of rolling bearings in different decay stages. Based on the theory of extreme learning mechanism, a data prediction model of extreme learning machine based on partial least squares method is proposed. The parameters in ELM are optimized, the number of nodes in the optimal hidden layer and the weight of the connecting layer are selected, and the activation function of Softmax is selected. In this way, the prediction accuracy is high, the calculation time is short, and the prediction effect of clustered eigenvalues is good. Through the above steps, the rolling bearing fault stage is accurately predicted.

【技术实现步骤摘要】
基于偏最小二乘法极限学习机的滚动轴承故障预测方法
本专利技术属于旋转机械故障预测方法的研究领域，具体涉及基于半正态分布与EWT相结合的数据降噪，并采用模糊C均值结合等距离映射的技术对非线性信号进行降维处理的方法。
技术介绍
现代工业设备由于其大规模、高复杂性、多变量、运行于闭环控制的特点，会存在导致设备故障的许多不可预防的影响因素。因此，设备功能丧失引起的一系列灾难性事故不在少数。人们也渐渐意识到设备的状态检测与故障预测在保证机组安全可靠地运行，预防事故中起到至关重要的作用。滚动轴承作为旋转机械重要的零部件之一，其运行状态的稳定性直接关系到整体设备的工作性能。有效预测滚动轴承的故障位置，正确识别滚动轴承性能衰退阶段，对降低危险事故的发生率和设备的健康管理都具有重要意义。国外方面，20世纪60年代，滚动轴承状态监测技术开始迅速发展。1962年，Gustafsson和Tallian研究发现可以通过滚动轴承振动加速度信号的峰值变化来振动轴承的早期故障。1974年，美国波音公司的D.R.Harting开启了共振解调
的先锋。1998年，NordenE.Huang等人提出了一种新型信号处理方法Hilbert-Huang变换，将信号分解成多个IMF分量并进行Hilbert变换。国内方面，20世纪80年代开始出现滚动轴承故障检测的相关研究。2014年，胥永刚等提出将双树复小波变换和独立分量分析(ICA)结合的故障诊断方法，该方法有效地分离和提取了滚动轴承复合故障的特征信息。隋文涛等又提出基于EMD与最大峭度解卷积(MKD)相结合的方法，对重构信号进行最大峭度解卷积处理以增强故障信息。现代机械设备故障预测的发展基础是基于信号处理与分析并结合机器学习的理论方法来实现。当机械设备运行时，由于现场环境及工况变化等因素的影响会导致信号中会掺杂大量的噪声，同时也会引起动态信号出现非平稳性。这些非平稳信号的统计量是时变函数，需要通过时频分析方法处理信号，傅里叶分析局限于线性，缓慢变换的信号；EMD(经验模态分解)的问题则是模态混叠，包络过多或过少都可能发生。EWT(经验小波变换)自适应划分初始信号的频谱，同时通过相对应的带通滤波器在各自的划分区间内构造正交带通滤波器组，但存在着区间定位模糊的问题。针对振动信号的故障特征集存在信息冗余，计算量较大，进而造成算法准确率极大降低的问题，我们需要对多特征参数进行维数约减。PCA(主成分分析)属于非监督算法，没有考虑样本类别信息，各故障阶段混杂在一起；KPCA(核主成分分析)通过与核函数相结合把数据投影至高维的空间，但核函数的选取严重影响降维效果。ISOMAP(等距离映射)将高维特征数据投影至相同的全局低维坐标中，降维后，虽可以明显区分，但还存在部分混叠现象。而采用改进的ISOMAP和模糊C均值结合的方法对滚动轴承衰退阶段的变化有很好的聚类效果，聚类中间比较集中且无混叠现象，聚类集中程度较高。机械设备健康维修主要实现手段之一是提取对机械设备的故障变化进行预测与评估，主要有基于失效物理模型的预测方法，利用产品的生命周期载荷与失效机理知识评估其可靠性，可在有些情况下并非需要建立动态模型，如此便会耗时耗力；基于统计的预测方法，所需信息来自各种概率密度函数，置信区间可以客观描述预测的准确度，但其合适的回归方程类型有时很难找到。基于数据驱动的故障预测方法是指通过分析输入、输出和状态参数之间的关系获得模型，在历史数据中学习输入与输出的映射关系，于其内部构建非线性非透明与非针对特定对象的模型，计算未来值，实现故障预测。综上，首先要提出合适的信号降噪的方法，贴近原始信号，剔除数据噪声信息。接着，提出合适的特征降维方法，提高滚动轴承衰退阶段的分类精度，且具有很好的损伤识别效果。其次，确保故障预测算法的精确度，更好地对滚动轴承的故障退化趋势进行精确预测。
技术实现思路
鉴上所述，本专利技术提出基于半正态分布和EWT(经验小波变换)相结合的降噪方法，真实地对滚动轴承故障信号进行诊断，又提出了基于C-ISOMAP的滚动轴承特征降维方法，提高滚动轴承衰退阶段的分类精度，接着是基于PLS-ELM(偏最小二乘法极限学习机)的故障预测模型，预测精度高，计算用时短，聚类后的特征值预测效果好。为达到上述目的，本专利技术提供如下技术方案：基于偏最小二乘法极限学习机的滚动轴承故障预测方法，所述的方法包括如下步骤：步骤一；基于半正态分布和的EWT信号降噪；分析传统EWT进行信号降噪时存在的局限性，提出新的EWT与半正态分布相结合的理论算法；步骤二；基于C-ISOMAP的滚动轴承数据降维；针对多特征参数对滚动轴承故障衰退特征进行评价，并对ISOMAP应用残差改进邻域大小，提出改进后的ISOMAP非线性特征降维方法与模糊C均值结合的方法，在滚动轴承整个监测过程中，区分滚动轴承不同损伤程度，对滚动轴承退化趋势进行预测和分类，将滚动轴承的正常振动信号和最终失效故障信号作为训练数据，建立模糊C聚类模型；步骤三；基于PLS-ELM的滚动轴承故障预测；在极限学习机基本理论的基础上，提出基于偏最小二乘法极限学习机的故障预测方法，采用偏最小二乘法法对ELM中的参数进行优化，选取最优隐含层神经元个数和各连接层之间的权重。本专利技术相对于现有技术的有益效果是：提出半正态分布与EWT相结合的理论算法，原始信号得到了精确重构。在保持EMD方法的完备性的基础上，相比EEMD筛选迭代次数降低，提高信噪比，从而有效地抑制噪声，更贴近原始信号，解决模态混叠问题。而改进后的ISOMAP非线性特征降维方法与模糊C均值结合的方法，更为精确诊断出滚动轴承在不同衰退阶段的变化趋势。比较经过模糊C均值聚类后的PCA、KPCA和改进后的ISOMAP效果图以及可分性指标可知，改进后的有更好的聚类效果，聚类中间比较即集中且无混叠现象，聚类集中程度较高。同时PLS-ELM的故障预测模型数据相关度更高，均方根误差更小，通过小波能谱熵可以准确预测各个故障阶段。附图说明图1为滚动轴承振动信号特征提取与故障预测研究工作的路线图。图2为EWT和半正态分布相结合的故障特征提取流程图。图3为小波信号分解图。图4为CEEMD信号分解图。图5为EWT信号分解图。图6为半正态分布-EWT降噪时频域图。图7为ISOMAP流形学习的模糊C均值流程图。图8为不同故障程度的残差分析图。图9为模糊C均值聚类后不同降维方法对比图。图10为不同降维方法的可分性指标对比图。图11为基于偏最小二乘法极限学习机的回归机预测数据。图12为基于PLS-ELM的小波能谱熵故障趋势预测图。具体实施方式具体实施方式一：如图1所示，本实施方式的基于偏最小二乘法极限学习机的滚动轴承故障预测方法，所述的方法包括如下步骤：步骤一；基于半正态分布和EWT的信号降噪；分析传统EWT进行信号降噪时存在的局限性，提出新的半正态分布(受大量、独立均匀小效应影响的随机变量都服从半正态分布。正因如此，机械故障领域，机械振动被视为半正态分布。)与EWT相结合的理论算法；步骤二；基于C-ISOMAP(建立模糊C聚类模型，用ISOMAP流行算法对非线性数据降维，处理后的低维数据均能够保持原有的拓扑关系)的滚动轴承数据降维；针对多特征参数对滚动轴承故障衰退特征进本文档来自技高网...

【技术保护点】
1.一种基于偏最小二乘法极限学习机的滚动轴承故障预测方法，其特征在于，包括以下步骤：步骤一；基于半正态分布‑EWT的信号降噪；分析时域、频域、时频域等多个特征指标，反应滚动轴承的故障诊断能力，提出新的半正态分布与经验小波变换相结合的理论算法；步骤二；基于C‑ISOMAP的特征降维；同时针对多特征参数对滚动轴承故障衰退特征进行评价，并对特征降维的多种方法进行理论分析，提出了残差改进的ISOMAP非线性特征降维方法与模糊C均值结合的方法；在滚动轴承整个监测过程中，不仅可以区分轴承不同损伤程度，而且可以对滚动轴承退化趋势进行预测和分类；将滚动轴承的正常振动信号和最终失效故障信号作为训练数据，建立模糊C聚类模型；步骤三；基于PLS‑ELM的滚动轴承故障预测；在极限学习机基本理论的基础上，提出基于偏最小二乘法极限学习机的故障阶段预测方法，采用偏最小二乘法对ELM中的参数进行优化，选取最优隐含层节点数和各连接层权值，并选取Softmax激活函数。

【技术特征摘要】
1.一种基于偏最小二乘法极限学习机的滚动轴承故障预测方法，其特征在于，包括以下步骤：步骤一；基于半正态分布-EWT的信号降噪；分析时域、频域、时频域等多个特征指标，反应滚动轴承的故障诊断能力，提出新的半正态分布与经验小波变换相结合的理论算法；步骤二；基于C-ISOMAP的特征降维；同时针对多特征参数对滚动轴承故障衰退特征进行评价，并对特征降维的多种方法进行理论分析，提出了残差改进的ISOMAP非线性特征降维方法与模糊C均值结合的方法；在滚动轴承整个监测过程中，不仅可以区分轴承不同损伤程度，而且可以对滚动轴承退化趋势进行预测和分类；将滚动轴承的正常振动信号和最终失效故障信号作为训练数据，建立模糊C聚类模型；步骤三；基于PLS-ELM的滚动轴承故障预测；在极限学习机基本理论的基础上，提出基于偏最小二乘法极限学习机的故障阶段预测方法，采用偏最小二乘法对ELM中的参数进行优化，选取最优隐含层节点数和各连接层权值，并选取Softmax激活函数。2.根据权利要求1所述的一种基于偏最小二乘法极限学习机的滚动轴承故障预测方法；首先要分析轴承振动信号，滚动轴承中当噪声信号足够大而覆盖了有用信息时，直接采用经验小波降噪去掉噪声信号，但经验小波变换应用传统的极大值在划分频率带时，紧支撑和宽支撑若是相邻的频带，其分界线会落在宽支撑内，无法辨别；其特征在于：步骤一具体步骤如下：步骤一一；由于滚动轴承原始振动信号是受大量、独立均匀小效应影响的随机变量，服从半正态分布；应用正态分布的公式确定频率带区间个数：式中，σ代表数据的标准差；步骤一二；确定区间个数后，应用EWT划分频率带区间得到含有噪声的高频AM-FM分量，并对其信...

【专利技术属性】
技术研发人员：王亚萍，崔巍，葛江华，
申请(专利权)人：哈尔滨理工大学，
类型：发明
国别省市：黑龙江,23