CN113642407A - 一种适用于轴承剩余使用寿命预测的特征提取优化方法 - Google Patents
一种适用于轴承剩余使用寿命预测的特征提取优化方法 Download PDFInfo
- Publication number
- CN113642407A CN113642407A CN202110799462.4A CN202110799462A CN113642407A CN 113642407 A CN113642407 A CN 113642407A CN 202110799462 A CN202110799462 A CN 202110799462A CN 113642407 A CN113642407 A CN 113642407A
- Authority
- CN
- China
- Prior art keywords
- bearing
- formula
- prediction
- monotonicity
- trend
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/04—Bearings
- G01M13/045—Acoustic or vibration analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
- G06F2218/06—Denoising by applying a scale-space analysis, e.g. using wavelet analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开一种适用于轴承剩余使用寿命(Remaining Useful Life,RUL)预测的特征提取优化方法,包括以下步骤:1)对滚动轴承原始振动信号采用离散小波变换方法进行信号去噪;2)计算出常用的轴承RUL预测统计特征;3)利用Savitzky‑Golay滤波器对统计特征进行平滑处理;4)对平滑处理后的统计特征进行累积变换;5)利用累积特征的单调性与趋势性来筛选用于轴承RUL预测的特征。本发明显著改善了特征的单调性和趋势性,优化后的特征有利于提高轴承剩余使用寿命的预测精度。
Description
技术领域
本发明涉及机械信号处理领域,尤其涉及一种适用于轴承剩余使用寿命预测的特征提取优化方法。
背景技术
随着当今社会迈入以智能制造为主导的工业4.0时代,机械设备日益朝着大型化、自动化、集成化、智能化的方向发展。为了保证设备的正常运转及便利维修,开展针对设备的在线监测和健康管理显得尤为重要。由于机械设备内部结构复杂,各子部件之间存在强耦合性,导致难以通过建立精确的物理模型用于监控和预测。因此,在传感器、存储、大数据等新技术快速发展的背景下,基于数据驱动的剩余使用寿命预测(remaining usefullife,RUL)技术近些年受到诸多关注。
滚动轴承被称为“工业的关节”,是旋转机械设备的关键零件,其性能与质量直接影响了机械设备的安全性与可靠性。因此,准确、及时地预测轴承的剩余使用寿命对工业生产安全而言具有重大研究意义。通过安装在机械设备内部的传感器实时采集数据,并提取反映轴承退化信息的特征,对于实现精确的RUL预测至关重要。特征提取是实现基于数据驱动的轴承RUL预测的前提,目的是提取出符合轴承退化趋势的特征用于表征退化过程。传统的特征提取方法是计算原始振动信号的统计特征,包括均方根、峭度、峰峰值等。近年来,深度学习凭借其出色的非线性映射能力在滚动轴承特征提取领域得到广泛应用。但是时间卷积网络、深度信念网络等深度学习方法往往需要大量标签数据进行有监督微调,标签数据的缺乏与难以获取严重制约了运用深度学习方法提取特征,而传统统计特征在单调性与趋势性上相比于轴承退化趋势往往不尽人意。因此,本发明提出一种适用于轴承剩余使用寿命预测的特征提取优化方法,通过对传统特征采用累积变换的方式优化其单调性与趋势性,使其在表征轴承退化过程中有更好的表现。
发明内容
本发明的目的在于采用一种适用于轴承剩余使用寿命预测的特征提取优化方法,通过将提取的特征变换为相应的累积形式,达到优化特征单调性与趋势性的效果,具有提高特征表征轴承退化过程能力的意义。
本发明的目的可通过以下的技术步骤来实现:
一种适用于轴承剩余使用寿命预测的特征提取优化方法,分以下步骤:
1)首先对滚动轴承原始振动信号采用离散小波变换方法进行信号去噪,包在消除异常的噪声信号的同时保留有用的退化信息;
2)在经过步骤1)信号去噪处理后得到的重构信号上,计算出常用的轴承RUL预测统计特征,包括峰峰值、标准差、均方根、三角特征等;
3)利用Savitzky-Golay滤波器对步骤2)所得的统计特征进行平滑处理,以减少提取特征的波动和进一步滤除不需要的噪声;
4)对步骤3)平滑处理后的统计特征进行累积变换,具体计算公式如下:
式中,cmn表示第m个特征ym (i)在n个样本中的累积变换结果。
5)利用累积特征的单调性与趋势性来筛选用于轴承RUL预测的特征。单调性公式计算方法如下:
式中:n表示总观测值。趋势性公式计算方法如下:
式中:xi和yi分别代表时间和特征的值,N为时间序列的长度。
本发明采用了一种有效的特征变换方法,即将提取的特征变换为相应的累积形式,通过累积函数对时间序列逐点执行累加操作和缩放操作,使用累积特征来表征退化的趋势。试验验证部分采用了2012年IEEE PHM预测挑战赛中的滚动轴承监测数据,通过对比各特征在累积变换前后的单调性、趋势性变化,得出本发明提出的特征提取优化方法可以提高特征的单调性和趋势性,有利于轴承剩余使用寿命的准确预测。
附图说明
图1为本发明提供的一种适用于轴承剩余使用寿命预测的特征提取优化方法的流程图;
图2为本发明提供的三角特征与传统特征对比图;
图3(a)为本发明提供的轴承偏度特征对比图;
图3(b)为本发明提供的轴承偏度特征累积值对比图;
图3(c)为本发明提供的轴承上边界特征对比图;
图3(d)为本发明提供的轴承上边界累积值特征对比图。
具体实施方式
以下结合附图和实施例对本发明技术方案展开完整详细的说明。
如图1所示,本发明提供一种适用于轴承剩余使用寿命预测的特征提取优化方法,包括以下步骤:
1)首先对滚动轴承原始振动信号采用离散小波变换方法进行信号去噪,包在消除异常的噪声信号的同时保留有用的退化信息;
2)在经过步骤1)信号去噪处理后得到的重构信号上,计算出常用的轴承RUL预测统计特征,包括峰峰值、标准差、均方根、三角特征等;
3)利用Savitzky-Golay滤波器对步骤2)所得的统计特征进行平滑处理,以减少提取特征的波动和进一步滤除不需要的噪声;
4)对步骤3)平滑处理后的统计特征进行累积变换,具体计算公式如下:
式中,cmn表示第m个特征ym (i)在n个样本中的累积变换结果。
5)利用累积特征的单调性与趋势性来筛选用于轴承RUL预测的特征。单调性公式计算方法如下:
式中:n表示总观测值。趋势性公式计算方法如下:
式中:xi和yi分别代表时间和特征的值,N为时间序列的长度。
进一步的优化方案,所述步骤1)的离散小波变换采用多贝西四阶小波(db4)为母小波,分解层数为4层,细节分量阈值处理中采用软阈值函数和最大最小阈值,最后通过小波重构得到去噪后的信号。
进一步的优化方案,所述步骤2)中的统计特征详见表1所示。
表1特征和相应的计算公式
进一步的优化方案,步骤3)中的Savitzky-Golay滤波器基于线性最小二乘法,通过采用低次多项式拟合相邻数据点的连续子集的方式实现平滑过程。
为了方便观察本发明提供的三角特征与传统特征之间的差异,选取反双曲正弦标准差、反正切值标准差两种三角特征和上边界、峰值两种传统特征为例进行对比。从图2可以看出,传统的统计特征表现出较低的单调性与趋势性,这不利于RUL预测模型学习轴承的退化趋势,同时较高的尺度也不利于预测模型的训练过程。而结合了标准差与三角函数的三角特征(提取方案见表1)具有更高的单调性与趋势性以及更低的尺度,有利于预测模型的训练过程。
传统的统计特征往往无法很好地表征轴承的退化趋势,如图3的(a)和(c)所示,仅经过步骤3)平滑处理后的特征单调性与趋势性各不相同,单调性较差。而从图3的(b)和(d)不难看出,累积变换后的特征表现出了更好的单调性与趋势性,说明采用的累积变换对于传统的统计特征的趋势性与单调性具有很好的优化作用。
进一步的优化方案,为了全面地说明累积变换的优势与必要性,采用步骤5)中的衡量特征单调性和趋势性的公式量化计算每种统计特征累积变化前后的单调性与趋势性,表2显示了统计特征累积变化前后的单调性与趋势性平均值,从表中可以看出,累积变化对本发明选取的统计特征的单调性与趋势性具有非常明显的改进作用。由于在滚动轴承RUL预测中,特征的单调性与趋势性越好,越能表征轴承的退化趋势,RUL预测的精度就越高。基于以上要求,舍弃经过累积变换后单调性和趋势性仍然较差的平均绝对值,保留剩余12个累积变换后的统计特征。
表2累积变换前后比较
综上所述,本发明提出一种适用于轴承剩余使用寿命预测的特征提取优化方法。首先,采用离散小波变换对滚动轴承原始振动信号进行信号去噪;其次,从去噪后的信号中提取统计特征;再次,利用Savitzky-Golay滤波器对统计特征进行平滑处理;然后,采用累积变换对统计特征进行改进,从而优化特征的单调性与趋势性;最后,筛选出单调性与趋势性效果显著的统计特征。本发明提出的适用于轴承剩余使用寿命预测的特征提取优化方法可解释性更强,优化后的特征有利于提高轴承剩余使用寿命的预测精度。该特征提取优化方法不仅能处理轴承的振动信号,还可以应用到其他机械设备关键零部件,具有良好的可扩展性。
Claims (1)
1.一种适用于轴承剩余寿命预测的特征提取优化方法,其特征在于,包括以下步骤:
1)首先对滚动轴承原始振动信号采用离散小波变换方法进行信号去噪,包在消除异常的噪声信号的同时保留有用的退化信息;
2)在经过步骤1)信号去噪处理后得到的重构信号上,计算出常用的轴承RUL预测统计特征,包括峰峰值、标准差、均方根、三角特征等;
3)利用Savitzky-Golay滤波器对步骤2)所得的统计特征进行平滑处理,以减少提取特征的波动和进一步滤除不需要的噪声;
4)对步骤3)平滑处理后的统计特征进行累积变换,具体计算公式如下:
式中,cmn表示第m个特征ym(i)在n个样本中的累积变换结果。
5)利用累积特征的单调性与趋势性来筛选用于轴承RUL预测的特征。单调性公式计算方法如下:
式中:n表示总观测值。趋势性公式计算方法如下:
式中:xi和yi分别代表时间和特征的值,N为时间序列的长度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110799462.4A CN113642407B (zh) | 2021-07-15 | 2021-07-15 | 一种适用于轴承剩余使用寿命预测的特征提取优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110799462.4A CN113642407B (zh) | 2021-07-15 | 2021-07-15 | 一种适用于轴承剩余使用寿命预测的特征提取优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113642407A true CN113642407A (zh) | 2021-11-12 |
CN113642407B CN113642407B (zh) | 2023-07-07 |
Family
ID=78417449
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110799462.4A Active CN113642407B (zh) | 2021-07-15 | 2021-07-15 | 一种适用于轴承剩余使用寿命预测的特征提取优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113642407B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150122025A1 (en) * | 2012-04-24 | 2015-05-07 | Aktiebolaget Skf | Method, computer program product & system |
CN109726524A (zh) * | 2019-03-01 | 2019-05-07 | 哈尔滨理工大学 | 一种基于cnn和lstm的滚动轴承剩余使用寿命预测方法 |
CN110610035A (zh) * | 2019-08-28 | 2019-12-24 | 武汉科技大学 | 一种基于gru神经网络的滚动轴承剩余寿命预测方法 |
-
2021
- 2021-07-15 CN CN202110799462.4A patent/CN113642407B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150122025A1 (en) * | 2012-04-24 | 2015-05-07 | Aktiebolaget Skf | Method, computer program product & system |
CN109726524A (zh) * | 2019-03-01 | 2019-05-07 | 哈尔滨理工大学 | 一种基于cnn和lstm的滚动轴承剩余使用寿命预测方法 |
CN110610035A (zh) * | 2019-08-28 | 2019-12-24 | 武汉科技大学 | 一种基于gru神经网络的滚动轴承剩余寿命预测方法 |
Non-Patent Citations (1)
Title |
---|
王正;: "基于DAE-HTPF的新能源汽车电池剩余寿命预测", 机械与电子, no. 03 * |
Also Published As
Publication number | Publication date |
---|---|
CN113642407B (zh) | 2023-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shen et al. | An automatic and robust features learning method for rotating machinery fault diagnosis based on contractive autoencoder | |
Gao et al. | Total variation on horizontal visibility graph and its application to rolling bearing fault diagnosis | |
CN111562108A (zh) | 一种基于cnn和fcmc的滚动轴承智能故障诊断方法 | |
CN111060337B (zh) | 一种基于专家系统的运转设备实时故障诊断方法 | |
CN113670616B (zh) | 一种轴承性能退化状态检测方法及系统 | |
CN116593157A (zh) | 少样本下基于匹配元学习的复杂工况齿轮故障诊断方法 | |
CN114034481A (zh) | 一种轧机齿轮箱故障诊断系统及方法 | |
CN115687973A (zh) | 一种基于多信号融合的机械设备故障诊断方法 | |
CN114819315A (zh) | 基于多参数融合健康因子和时间卷积神经网络的轴承退化趋势预测方法 | |
CN114755010A (zh) | 一种旋转机械振动故障诊断方法及其系统 | |
CN114626625B (zh) | 基于多模块U-BiLSTM网络的轴承寿命预测方法 | |
CN114462446A (zh) | 基于振动信号的滚动轴承故障诊断方法及电子设备 | |
CN117171657A (zh) | 风力发电设备故障诊断方法、装置、电子设备及存储介质 | |
CN111291918A (zh) | 平稳子空间外源矢量自回归的旋转机械退化趋势预测方法 | |
CN112648220A (zh) | 一种基于小波-近似熵的风机故障诊断方法 | |
Zhao et al. | A feature extraction method based on LMD and MSE and its application for fault diagnosis of reciprocating compressor | |
CN113642407A (zh) | 一种适用于轴承剩余使用寿命预测的特征提取优化方法 | |
CN116907838A (zh) | 基于增强关系网络的齿轮故障诊断方法 | |
CN114077850B (zh) | 一种变工况下的基于图数据的旋转机械设备状态监测方法 | |
CN112069621B (zh) | 基于线性可靠度指标的滚动轴承剩余使用寿命的预测方法 | |
CN113836476B (zh) | 一种基于智能手机的路面平整度测量方法 | |
CN114218980A (zh) | 滚动轴承剩余寿命预测方法 | |
CN110595811B (zh) | 机械设备健康状态特征图构建方法 | |
CN107704825A (zh) | 基于自适应集成多小波的机械设备故障特征提取方法 | |
CN112683533A (zh) | 一种用于轴承故障诊断的信号增强方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |