CN113609770A - 基于分段线性拟合hi及lstm的滚动轴承rul预测方法 - Google Patents
基于分段线性拟合hi及lstm的滚动轴承rul预测方法 Download PDFInfo
- Publication number
- CN113609770A CN113609770A CN202110900150.8A CN202110900150A CN113609770A CN 113609770 A CN113609770 A CN 113609770A CN 202110900150 A CN202110900150 A CN 202110900150A CN 113609770 A CN113609770 A CN 113609770A
- Authority
- CN
- China
- Prior art keywords
- curve
- rolling bearing
- degradation period
- sequence
- segments
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/044—Recurrent networks, e.g. Hopfield networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/04—Ageing analysis or optimisation against ageing
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Software Systems (AREA)
- Computing Systems (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- Mathematical Physics (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Medical Informatics (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种基于分段线性拟合HI及LSTM的滚动轴承RUL预测方法,获取预设工况下若干滚动轴承的HI曲线并获取剩余寿命序列,基于聚类评价指标和拟合评价指标筛选出每个滚动轴承采用BUP时间序列分割算法进行分段的最优分段数,对HI曲线进行分段并提取出退化期HI曲线以及退化期剩余寿命序列,进行归一化得到归一化退化期HI曲线和归一化退化期剩余寿命序列,采用长短时记忆网络作为RUL预测模型,将归一化退化期HI曲线作为输入,将归一化退化期剩余寿命序列作为标签,对RUL预测模型进行训练,对于预设工况下的某个滚动轴承,采用相同方法得到归一化退化期HI曲线,将其输入RUL预测模型得到预测的剩余寿命序列。本发明可以有效提高RUL的预测准确率和稳健性。
Description
技术领域
本发明属于滚动轴承寿命预测技术领域,更为具体地讲,涉及一种基于分段线性拟合HI及LSTM的滚动轴承RUL预测方法。
背景技术
在目前旋转机械的所有关键组件中,滚动轴承一直是研究的重要对象之一。一方面由于滚动轴承使用广泛,在旋转机械设备应用中有着不可替代的作用;另一方面是由于滚动轴承相较于其他的零部件更易失效。由于滚动轴承的服役环境和工作条件的复杂性和不确定性、疲劳损伤发展的随机性和失效模式的多样性,导致滚动轴承的最大使用寿命有较大的离散性,并且滚动轴承的定期维护经常导致“维护不足”和“过度维护”的问题。而滚动轴承的RUL(Remaining Useful Life,剩余使用寿命)预测可以最大限度的延长轴承的使用寿命并降低维护成本,因此对滚动轴承的RUL预测的研究一直是该领域研究的重点。
目前已有的关于旋转机械的RUL预测的研究方法大致可以分为基于模型和基于数据驱动的预测方法两大类。基于模型的预测方法通过建立物理或数学模型来描述系统或设备的退化趋势,此方法要求大量的物理机理或经验知识,但是面对现代结构复杂的机械设备系统和多变的服役环境很难建立有效的失效模型,因此基于模型的预测方法的应用有限。基于数据驱动的RUL预测方法没有上述方面的要求,对设备的退化过程建模一般只需要设备运行过程中的采集和监测的数据。在数据驱动的方法主要分为统计模型和人工智能两大类。神经网络等人工智能方法是近年来针对各种序列分类和预测问题的最先进的模型之一,神经网络方法以其建模优势和极强的非线性处理能力,逐渐成为主流研究方法。同时,在当今大数据时代下,在工程实际应用中,现有的基于神经网络的寿命预测技术存在模型输入样本量太大、对退化期的样本定位模糊以及在多工况和多失效模式下的寿命预测的精度较低等问题,需要研究解决。
发明内容
本发明的目的在于克服现有技术的不足,提供一种基于分段线性拟合HI及LSTM的滚动轴承RUL预测方法,对HI曲线采用BUP时间序列分割算法进行分段得到退化期HI曲线,采用LSTM作为RUL预测模型,将退化期HI曲线作为预测模型的输入,实现对滚动轴承的RUL预测,提高RUL的预测准确率和稳健性。
为实现上述发明目的,本发明基于分段线性拟合HI及LSTM的滚动轴承RUL预测方法包括以下步骤:
S1:根据实际需要获取预设工况下若干滚动轴承的特征数据序列,根据特征数据序列绘制各个滚动轴承对应的HI曲线,同时获取该HI曲线对应的滚动轴承剩余寿命序列,其中每个值即为对应时刻该滚动轴承的RUL值;
S2:对于每个滚动轴承的HI曲线,分别确定该HI曲线使用BUP时间序列分割算法进行健康状态划分的最优分段数,具体方法为:
S2.1:根据实际情况设置K个备选分段数Mk,其中k=1,2,…,K,分别以每个备选分段数作为BUP时间序列分割算法的循环迭代停止条件,使用BUP时间序列分割算法对HI曲线对应的特征数据序列进行健康状态划分,得到Mk个分段特征数据序列;
S2.2:将备选分段数Mk对应的每个分段特征数据序列分别作为一个数据类别,计算得到聚类评价指标CHk;对于Mk个分段特征数据序列,计算每个分段特征数据序列的均方根RMSk,m,m=1,2,…,Mk,然后采用以下公式计算得到拟合评价指标Costk:
S2.5:从K个备选分段数Mk中筛选出聚类评价指标CHk次最大的备选分段数Mk′;
S2.8:将分段数Mk′作为最优分段数;
S3:将步骤S2得到的最优分段数作为BUP时间序列分割算法的循环迭代停止条件,使用BUP分割算法对HI曲线进行健康状态划分,得到分段线性拟合后的HI曲线;
S4:对于每个滚动轴承分段线性拟合后的HI曲线,将每两个健康状态之间的分界点作为健康状态转折点,并将第一个健康状态转折点作为初始故障发生点,以初始故障发生点区分滚动轴承的健康期与退化期,初始故障发生点之前为健康期,之后为退化期;
S5:对于每个滚动轴承分段线性拟合后的HI曲线分别提取出退化期HI曲线,并从该HI曲线对应的滚动轴承剩余寿命序列提取出退化期HI曲线对应的退化期剩余寿命序列,分别对退化期HI曲线和退化期剩余寿命序列进行归一化处理,得到归一化退化期HI曲线和归一化退化期剩余寿命序列;
S6:采用长短时记忆网络作为滚动轴承的RUL预测模型,将步骤S5中的归一化退化期HI曲线作为输入,将归一化退化期剩余寿命序列作为标签,对RUL预测模型进行训练;
S7:对于预设工况下的某个滚动轴承,获取其到当前时刻的特征数据序列,绘制对应的HI曲线,采用步骤S2中的相同方法确定其最优分段数,采用步骤S3中的相同方法得到分段线性拟合后的HI曲线,采用步骤S4中的相同方法确定退化期,然后从该滚动轴承分段线性拟合后的HI曲线分别提取出退化期HI曲线并进行归一化,得到归一化退化期HI曲线,将该归一化退化期HI曲线输入步骤S6训练好的RUL预测模型中,得到预测的剩余寿命序列。
本发明基于分段线性拟合HI及LSTM的滚动轴承RUL预测方法,获取预设工况下若干滚动轴承的特征数据序列,绘制对应的HI曲线并获取剩余寿命序列,基于聚类评价指标和拟合评价指标筛选出每个滚动轴承采用BUP时间序列分割算法进行分段的最优分段数,对HI曲线进行分段并提取出退化期HI曲线以及退化期剩余寿命序列,进行归一化得到归一化退化期HI曲线和归一化退化期剩余寿命序列,采用长短时记忆网络作为RUL预测模型,将归一化退化期HI曲线作为输入,将归一化退化期剩余寿命序列作为标签,对RUL预测模型进行训练,对于预设工况下的某个滚动轴承绘制HI曲线,采用相同方法得到归一化退化期HI曲线,将其输入RUL预测模型得到预测的剩余寿命序列。
本发明具有以下有益效果:
1)本发明利用BUP时间序列分割算法对全寿命周期的HI曲线进行了健康状态划分,识别出了初始故障发生点,仅将初始故障发生点之后的退化期HI曲线用来预测RUL,一方面避免了RUL预测模型的过大输入样本量带来的较高的硬件要求和训练负担,另一方面还优化了RUL预测模型的输入,提高了预测结果的准确度和稳健性
2)本发明提出的基于分段线性拟合HI及LSTM的滚动轴承RUL预测方法具有预测模型的结构简单、对硬件资源的要求低,并且经实验证明本发明可以有效提高在多工况和多失效模型下的滚动轴承的RUL预测的准确度和稳健度。
附图说明
图1是本发明基于分段线性拟合HI及LSTM的滚动轴承RUL预测方法的具体实施方式流程图;
图2是本发明中确定最优分段数的流程图;
图3是本实施例中滚动轴承的HI曲线图;
图4是本实施例中每个滚动轴承的健康状态的划分结果图;
图5是本实施例中测试集轴承Bearing1_5的预测标签和真实标签之间的对比图;
图6是本实施例中测试集轴承Bearing2_5的预测标签和真实标签之间的对比图;
图7是本实施例中测试集轴承Bearing3_5的预测标签和真实标签之间的对比图;
图8是本实施例中本发明与对比方法的20次预测结果的RMSE对比图;
图9是本实施例中本发明与对比方法的20次预测结果的Score对比图。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述,以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是,在以下的描述中,当已知功能和设计的详细描述也许会淡化本发明的主要内容时,这些描述在这里将被忽略。
实施例
图1是本发明基于分段线性拟合HI及LSTM的滚动轴承RUL预测方法的具体实施方式流程图。如图1所示,本发明基于分段线性拟合HI及LSTM的滚动轴承RUL预测方法的具体步骤包括:
S101:获取HI曲线:
根据实际需要获取预设工况下若干滚动轴承的特征数据序列,根据特征数据序列绘制各个滚动轴承对应的HI(健康指标)曲线,同时获取该HI曲线对应的滚动轴承剩余寿命序列,其中每个值即为对应时刻该滚动轴承的RUL值。
S102:确定最优分段数:
对于每个滚动轴承的HI曲线,分别确定该HI曲线使用BUP时间序列分割算法进行健康状态划分的最优分段数。图2是本发明中确定最优分段数的流程图。如图2所示,本发明中确定最优分段数的具体步骤包括:
S201:获取备选分段数的分段结果:
根据实际情况设置K个备选分段数Mk,其中k=1,2,…,K,分别以每个备选分段数作为BUP时间序列分割算法的循环迭代停止条件,使用BUP时间序列分割算法对HI曲线对应的特征数据序列进行健康状态划分,得到Mk个分段特征数据序列。BUP时间序列分割算法的具体原理的步骤可以参见文献“E.Keogh,S.Chu,D.Hart,et al.Segmenting Time Series:A Survey and Novel Approach[M].Singapore:Data Mining in Time SeriesDatabases,2004,1-21”,本实施例中对其具体算法简述如下:
记用于绘制HI曲线的特征数据序列为S,序列长度为N,首先创建序列S的N/2个分段(Seg0,Seg1,…SegN/2-1)来近似该长度为N的特征数据序列S,将此N/2个分段作为时间序列S的初始最佳近似,对每个相邻分段对(Segi,Segi+1)使用最小二乘方法进行拟合合并,i=0,1,…,N/2-2,并计算对应的合并代价,合并代价为合并前的相邻段对与合并后的拟合曲线之间的均方误差(Mean Square Error,MSE);
在迭代过程中,当一对合并代价最小的相邻分段(SegminInd,SegminInd+1)合并时,需要进行以下几步操作:
步骤一:记录保存临时合并代价cost(minInd);
步骤二:将合并代价最小的相邻分段对(SegminInd,SegminInd+1)合并到新的段Snew中;
步骤三:从分段集合中删除相邻合并段(SegminInd-1,SegminInd)、(SegminInd+1,SegminInd+2)的合并代价merge_cost(minInd)、merge_cost(minInd+1),将步骤一中保存的临时合并代价cost(minInd)添加至合并代价merge_cost(minInd)中;
步骤四:重新计算Snew与其左邻居SminInd-1的合并成本;
步骤五:重新计算Snew与其右邻居SminInd+2的合并成本;
当分段数达到预设的备选分段数时,循环迭代停止之后,得到分段结果。
结合一些学者的研究和大多数滚动轴承的实际退化过程,本实施例中将滚动轴承的退化过程中经历的健康状态的阶段个数限定在2-4范围内,即每个滚动轴承的备选分段数为2、3、4;。
S202:计算分段结果评价指标:
对于备选分段数Mk的分段特征数据序列,计算其聚类评价指标CHk和拟合评价指标Costk。
聚类评价指标CH的全称为Calinski-Harabaz,是一种无监督聚类评价指标。类别内部数据的协方差越小且类别之间的协方差越大CH越高,即CH分数越高说明聚类质量越好。本发明中,将备选分段数Mk对应的每个分段特征数据序列作为一个数据类别,计算得到聚类评价指标CHk,本实施例中聚类评价指标的计算公式如下:
其中,tr()表示求取矩阵的迹,Bk表示类别之间的协方差矩阵,Wk表示类别内部数据的协方差矩阵。
拟合评价指标Costk的计算方法如下:对于Mk个分段特征数据序列,计算每个分段特征数据序列的均方根RMSk,m,m=1,2,…,Mk,然后采用以下公式计算得到拟合评价指标Costk:
S203:筛选最大聚类评价指标的备选分段数:
S205:筛选次最大聚类评价指标的备选分段数:
从K个备选分段数Mk中筛选出聚类评价指标CHk次最大的备选分段数Mk′。
S208:将分段数Mk′作为最优分段数。
S103:健康状态划分:
将步骤S102得到的最优分段数作为BUP时间序列分割算法的循环迭代停止条件,使用BUP分割算法对HI曲线进行健康状态划分,得到分段线性拟合后的HI曲线。
S104:确定初始故障发生点及退化期:
对于每个滚动轴承分段线性拟合后的HI曲线,将每两个健康状态之间的分界点作为健康状态转折点,并将第一个健康状态转折点作为初始故障发生点,以初始故障发生点区分滚动轴承的健康期与退化期,初始故障发生点之前为健康期,之后为退化期。
S105:确定训练样本:
对于每个滚动轴承分段线性拟合后的HI曲线分别提取出退化期HI曲线,并从该HI曲线对应的滚动轴承剩余寿命序列提取出退化期HI曲线对应的退化期剩余寿命序列,分别对退化期HI曲线和退化期剩余寿命序列进行归一化处理,得到归一化退化期HI曲线和归一化退化期剩余寿命序列。
退化期HI曲线的归一化采用最大最小归一化,退化期剩余寿命序列归一化采用每个时刻的RUL值除以该滚动轴承的最大寿命。
S106:搭建并训练RUL预测模型:
采用长短时记忆网络(LSTM)作为滚动轴承的RUL预测模型,将步骤S105中的归一化退化期HI曲线作为输入,将归一化退化期剩余寿命序列作为标签,对RUL预测模型进行训练。
S107:RUL预测:
对于预设工况下的某个滚动轴承,获取其到当前时刻的特征数据序列,绘制对应的HI曲线,采用步骤S102中的相同方法确定其最优分段数,采用步骤S103中的相同方法得到分段线性拟合后的HI曲线,采用步骤S104中的相同方法确定退化期,然后从该滚动轴承分段线性拟合后的HI曲线分别提取出退化期HI曲线并进行归一化,得到归一化退化期HI曲线,将该归一化退化期HI曲线输入步骤S106训练好的RUL预测模型中,得到预测的剩余寿命序列。
为了更好地说明本发明的技术方案和技术效果,采用一个具体实例对本发明的工作流程和技术效果进行分析说明。轴承故障是旋转机械中的一种典型故障,因此本实施例采用西安交通大学轴承数据集XJTU-SY的轴承加速全寿命振动信号数据进行实验测试。
本实施例中使用的西安交通大学的轴承加速寿命试验共设计了3类工况,在每一种工况下对5个滚动轴承进行加速寿命测试,表1是本实施列中基于西安交通大学的轴承数据的具体试验工况。
工况编号 | 工况一 | 工况二 | 工况三 |
转速(r/min) | 2100 | 2250 | 2400 |
径向力(kN) | 12 | 11 | 10 |
表1
本实施例中振动加速度信号数据集通过单向加速度传感器进行采集,采样频率为25.6kHz,采样间隔为1min,每次的采样时长为1.28s。包含了3种工况下的15个滚动轴承的全周期寿命振动信号,且明确标注了各失效轴承的样本总数、基本额定寿命、实际寿命以及故障位置等信息。
将15个滚动轴承的水平振动加速度数据进行特征提取、特征优选后得到的优选特征进行特征融合,根据融合后的特征数据序列得到各自的HI曲线。图3是本实施例中滚动轴承的HI曲线图。
对每个滚动轴承的HI曲线经过最优分段数的选择算法确定健康状态的最佳划分个数。表2是每个滚动轴承的健康状态最佳的划分个数。
表2
将表2中每个滚动轴承的健康状态的最佳划分个数作为BUP时间序列分割算法的停止条件,对各个滚动轴承的HI曲线进行健康状态划分,并保存每条HI曲线的分段线性拟合数据。图4是本实施例中每个滚动轴承的健康状态的划分结果图。
将第一个健康状态转折点作为初始故障发生点,则可以确定本实施例中每个滚动轴承的初始故障发生点。表3是每个滚动轴承的初始故障发生点。初始故障发生点之后的HI曲线即为退化期HI曲线。
表3
将处于退化期所有HI曲线的分段线性拟合数据作为RUL预测模型的输入样本,分为训练样本集和测试样本集。本实施例选取每种工况下的前4个滚动轴承作为RUL预测模型的训练集样本,选取每种工况下的最后一个滚动轴承作为RUL预测模型的测试集样本。表4是本实施例划分后的训练集和测试集的具体情况。
表4
由于LSTM在处理长时间序列方面的优越性,利用LSTM搭建滚动轴承的RUL预测模型。考虑到本实施例中处于退化期的HI曲线的训练样本量较少,为了避免过拟合,本实施例搭建了一个结构较为简单的LSTM预测模型,即隐藏层的层数和各隐藏层的节点数都较少的LSTM预测模型。表5是LSTM预测模型的超参数。本实施例采用的LSTM预测模型结构共三层,其中包括一个LSTM层,两个Dense层,每一层的节点数分别为32-8-1,迭代次数取300,损失函数为MSE,优化器为Adam。
参数编号 | 参数名称 | 参数值 |
1 | 输入节点数 | 3 |
2 | 隐藏单元节点数 | 32-8 |
3 | 激活函数 | ‘Relu’-‘Relu’-‘Linear’ |
4 | 输出节点数 | 1 |
5 | 批尺寸 | 64 |
6 | 迭代次数 | 300 |
7 | 优化器 | Adam |
8 | 损失函数 | MSE |
表5
利用构建好LSTM预测模型对表4中训练集的12个滚动轴承的样本特征及标签完成迭代训练。然后将表4中测试集的3个滚动轴承的HI曲线的分段线性拟合数据送入训练好的LSTM预测模型进行前向传播,得到各个测试集轴承的RUL预测结果。最后与各测试集轴承真实的RUL标签进行对比以评估搭建的LSTM预测模型的能力的强弱。图5是本实施例中测试集轴承Bearing1_5的预测标签和真实标签之间的对比图。图6是本实施例中测试集轴承Bearing2_5的预测标签和真实标签之间的对比图。图7是本实施例中测试集轴承Bearing3_5的预测标签和真实标签之间的对比图。从图5、图6和图7可以看到,用来测试的3个滚动轴承的预测趋势均能描述并跟随轴承真实变化的退化趋势,预测的RUL曲线与真实RUL曲线具有较高的吻合程度,说明LSTM预测模型对滚动轴承的RUL预测效果较好。由此说明本发明基于LSTM和分段线性拟合HI的RUL预测模型能够适用于不同工况、不同退化趋势和不同失效模式下的滚动轴承的RUL预测。
为了验证本发明提出的将HI曲线的分段线性拟合数据作为LSTM预测模型的输入的优越性,将其预测结果与原始的HI曲线、振动信号提取的均方根(RMS)特征曲线作为输入的预测结果进行对比。为了避免单次实验结果的偶然性和探究不同输入在多次实验表现出的准确度和稳定度,本实施例共重复了20次相同的实验。然后分别计算了三种不同输入下的20次预测结果的评估指标:均方根误差RMSE和打分函数Score。
图8是本实施例中本发明与对比方法的20次预测结果的RMSE对比图。图9是本实施例中本发明与对比方法的20次预测结果的Score对比图。RMSE值和Score值均为3个滚动轴承的RMSE和Score之和。图8和图9的结果表明,就预测稳定度而言,HI曲线的分段线性拟合数据在20次重复实验中表现出来的预测稳定度较好,同时RMSE和Score评估指标也比其他两种输入更小,说明预测误差更小,即预测准确度更好。因此,综合预测准确度和预测稳健性考虑,本发明采用的HI曲线的分段线性拟合数据作为输入时预测得到的RUL的准确度和稳定度均更高。
表6是三种不同输入下的20次平均预测评估指标对比表。分别对比了3个滚动轴承的20次实验的平均RMSE指标和平均Score指标的得分,并计算了3个轴承RMSE指标和Score得分的总和。
表6
表6结果表明,采用原始HI数据和分段线性拟合数据作为预测模型的输入,均比直接使用单个特征RMS作为输入所获得的预测误差更低,说明通过特征融合构建的HI相对于单个特征可以提升RUL的预测效果。另外也可以看到,采用分段线性拟合的HI数据作为输入相比于原始HI数据作为输入,所获得的RMSE和Score评估指标的总计得分都更小,说明本发明采用的分段线性拟合数据作为输入获得的RUL预测误差更低,验证了本发明采用分段线性拟合数据提升了RUL预测模型的预测精度。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (1)
1.一种基于分段线性拟合HI及LSTM的滚动轴承RUL预测方法,其特征在于,包括以下步骤:
S1:根据实际需要获取预设工况下若干滚动轴承的特征数据序列,根据特征数据序列绘制各个滚动轴承对应的HI曲线,同时获取该HI曲线对应的滚动轴承剩余寿命序列,其中每个值即为对应时刻该滚动轴承的RUL值;
S2:对于每个滚动轴承的HI曲线,分别确定该HI曲线使用BUP时间序列分割算法进行健康状态划分的最优分段数,具体方法为:
S2.1:根据实际情况设置K个备选分段数Mk,其中k=1,2,…,K,分别以每个备选分段数作为BUP时间序列分割算法的循环迭代停止条件,使用BUP时间序列分割算法对HI曲线对应的特征数据序列进行健康状态划分,得到Mk个分段特征数据序列;
S2.2:将备选分段数Mk对应的每个分段特征数据序列作为一个数据类别,计算得到聚类评价指标CHk;对于Mk个分段特征数据序列,计算每个分段特征数据序列的均方根RMSk,m,m=1,2,…,Mk,然后采用以下公式计算得到拟合评价指标Costk:
S2.5:从K个备选分段数Mk中筛选出聚类评价指标CHk次最大的备选分段数Mk′;
S2.8:将分段数Mk′作为最优分段数;
S3:将步骤S2得到的最优分段数作为BUP时间序列分割算法的循环迭代停止条件,使用BUP分割算法对HI曲线进行健康状态划分,得到分段线性拟合后的HI曲线;
S4:对于每个滚动轴承分段线性拟合后的HI曲线,将每两个健康状态之间的分界点作为健康状态转折点,并将第一个健康状态转折点作为初始故障发生点,以初始故障发生点区分滚动轴承的健康期与退化期,初始故障发生点之前为健康期,之后为退化期;
S5:对于每个滚动轴承分段线性拟合后的HI曲线分别提取出退化期HI曲线,并从该HI曲线对应的滚动轴承剩余寿命序列提取出退化期HI曲线对应的退化期剩余寿命序列,分别对退化期HI曲线和退化期剩余寿命序列进行归一化处理,得到归一化退化期HI曲线和归一化退化期剩余寿命序列;
S6:采用长短时记忆网络作为滚动轴承的RUL预测模型,将步骤S5中的归一化退化期HI曲线作为输入,将归一化退化期剩余寿命序列作为标签,对RUL预测模型进行训练;
S7:对于预设工况下的某个滚动轴承,获取其到当前时刻的特征数据序列,绘制对应的HI曲线,采用步骤S2中的相同方法确定其最优分段数,采用步骤S3中的相同方法得到分段线性拟合后的HI曲线,采用步骤S4中的相同方法确定退化期,然后从该滚动轴承分段线性拟合后的HI曲线分别提取出退化期HI曲线并进行归一化,得到归一化退化期HI曲线,将该归一化退化期HI曲线输入步骤S6训练好的RUL预测模型中,得到预测的剩余寿命序列。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110900150.8A CN113609770B (zh) | 2021-08-06 | 2021-08-06 | 基于分段线性拟合hi及lstm的滚动轴承rul预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110900150.8A CN113609770B (zh) | 2021-08-06 | 2021-08-06 | 基于分段线性拟合hi及lstm的滚动轴承rul预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113609770A true CN113609770A (zh) | 2021-11-05 |
CN113609770B CN113609770B (zh) | 2023-06-02 |
Family
ID=78307359
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110900150.8A Active CN113609770B (zh) | 2021-08-06 | 2021-08-06 | 基于分段线性拟合hi及lstm的滚动轴承rul预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113609770B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116522467A (zh) * | 2023-01-16 | 2023-08-01 | 北京控制工程研究所 | 一种航天器部件寿命的预测方法和装置 |
CN116908707A (zh) * | 2023-09-14 | 2023-10-20 | 中国华能集团清洁能源技术研究院有限公司 | 一种基于sa-lstm神经网络的锂电池剩余寿命预测方法及设备 |
CN117252083A (zh) * | 2023-07-12 | 2023-12-19 | 中国科学院空间应用工程与技术中心 | 一种结合退化阶段划分和子域自适应的轴承剩余寿命预测方法及系统 |
CN118364418A (zh) * | 2024-06-20 | 2024-07-19 | 无锡中基电机制造有限公司 | 用于轴承座的耐腐性智能检测方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190384257A1 (en) * | 2018-06-13 | 2019-12-19 | Hitachi, Ltd. | Automatic health indicator learning using reinforcement learning for predictive maintenance |
CN111832759A (zh) * | 2020-07-13 | 2020-10-27 | 福建省港航勘察设计院有限公司 | 基于bim的水运工程分类编码方法 |
CN112101220A (zh) * | 2020-09-15 | 2020-12-18 | 哈尔滨理工大学 | 一种无监督模型参数迁移的滚动轴承寿命预测方法 |
-
2021
- 2021-08-06 CN CN202110900150.8A patent/CN113609770B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190384257A1 (en) * | 2018-06-13 | 2019-12-19 | Hitachi, Ltd. | Automatic health indicator learning using reinforcement learning for predictive maintenance |
CN111832759A (zh) * | 2020-07-13 | 2020-10-27 | 福建省港航勘察设计院有限公司 | 基于bim的水运工程分类编码方法 |
CN112101220A (zh) * | 2020-09-15 | 2020-12-18 | 哈尔滨理工大学 | 一种无监督模型参数迁移的滚动轴承寿命预测方法 |
Non-Patent Citations (1)
Title |
---|
黄金苑;李少波;张安思;杨万里;刘慧斌;胡建军;: "基于深度卷积网络的多工况寿命预测方法研究" * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116522467A (zh) * | 2023-01-16 | 2023-08-01 | 北京控制工程研究所 | 一种航天器部件寿命的预测方法和装置 |
CN116522467B (zh) * | 2023-01-16 | 2023-10-31 | 北京控制工程研究所 | 一种航天器部件寿命的预测方法和装置 |
CN117252083A (zh) * | 2023-07-12 | 2023-12-19 | 中国科学院空间应用工程与技术中心 | 一种结合退化阶段划分和子域自适应的轴承剩余寿命预测方法及系统 |
CN116908707A (zh) * | 2023-09-14 | 2023-10-20 | 中国华能集团清洁能源技术研究院有限公司 | 一种基于sa-lstm神经网络的锂电池剩余寿命预测方法及设备 |
CN116908707B (zh) * | 2023-09-14 | 2023-12-26 | 中国华能集团清洁能源技术研究院有限公司 | 一种基于sa-lstm神经网络的锂电池剩余寿命预测方法及设备 |
CN118364418A (zh) * | 2024-06-20 | 2024-07-19 | 无锡中基电机制造有限公司 | 用于轴承座的耐腐性智能检测方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN113609770B (zh) | 2023-06-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113609770B (zh) | 基于分段线性拟合hi及lstm的滚动轴承rul预测方法 | |
CN110361176B (zh) | 一种基于多任务特征共享神经网络的智能故障诊断方法 | |
CN107294993B (zh) | 一种基于集成学习的web异常流量监测方法 | |
CN111459700B (zh) | 设备故障的诊断方法、诊断装置、诊断设备及存储介质 | |
CN109492026B (zh) | 一种基于改进的主动学习技术的电信欺诈分类检测方法 | |
CN110082136B (zh) | 基于云遗传算法优化支持向量机的旋转机械故障诊断方法 | |
CN113326244B (zh) | 一种基于日志事件图和关联关系挖掘的异常检测方法 | |
CN110210169B (zh) | 一种基于lstm的盾构机故障预测方法 | |
CN110599459A (zh) | 基于深度学习的地下管网风险评估云系统 | |
CN112560327B (zh) | 一种基于深度梯度下降森林的轴承剩余寿命预测方法 | |
CN107679135A (zh) | 面向网络文本大数据的话题检测与跟踪方法、装置 | |
CN108683658B (zh) | 基于多rbm网络构建基准模型的工控网络流量异常识别方法 | |
CN113076738A (zh) | Gnn编码器及基于图上下文学习的异常点检测方法 | |
CN115587543A (zh) | 基于联邦学习和lstm的刀具剩余寿命预测方法及系统 | |
CN109543693A (zh) | 基于正则化标签传播的弱标注数据降噪方法 | |
CN114819260A (zh) | 一种水文时间序列预测模型动态生成方法 | |
CN115982655A (zh) | 一种基于决策树的缺失数据流程异常预测方法 | |
Fan | Data mining model for predicting the quality level and classification of construction projects | |
CN112749905A (zh) | 一种基于大数据挖掘的洪涝灾害评估方法 | |
CN113448840A (zh) | 基于预测缺陷率和模糊综合评价模型的软件质量评价方法 | |
CN117172381A (zh) | 基于大数据的风险预测方法 | |
CN116881640A (zh) | 岩心采出程度的预测方法、系统及计算机可读存储介质 | |
CN114547796B (zh) | 一种基于优化bn网络的球磨机特征融合故障诊断方法 | |
CN114298413B (zh) | 一种水电机组振摆趋势预测方法 | |
CN110348005A (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 |