CN111783242B - 一种基于rvm-kf的滚动轴承剩余寿命预测方法及装置 - Google Patents
一种基于rvm-kf的滚动轴承剩余寿命预测方法及装置 Download PDFInfo
- Publication number
- CN111783242B CN111783242B CN202010554878.5A CN202010554878A CN111783242B CN 111783242 B CN111783242 B CN 111783242B CN 202010554878 A CN202010554878 A CN 202010554878A CN 111783242 B CN111783242 B CN 111783242B
- Authority
- CN
- China
- Prior art keywords
- data
- rolling bearing
- rvm
- prediction
- predicted
- 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.)
- Active
Links
- 238000005096 rolling process Methods 0.000 title claims abstract description 114
- 238000000034 method Methods 0.000 title claims abstract description 64
- 238000006731 degradation reaction Methods 0.000 claims abstract description 67
- 230000015556 catabolic process Effects 0.000 claims abstract description 61
- 238000012544 monitoring process Methods 0.000 claims abstract description 21
- 238000012549 training Methods 0.000 claims abstract description 20
- 230000015654 memory Effects 0.000 claims description 13
- 238000004590 computer program Methods 0.000 claims description 10
- 230000000694 effects Effects 0.000 claims description 8
- 230000005653 Brownian motion process Effects 0.000 claims description 4
- 238000005537 brownian motion Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 description 14
- 239000011159 matrix material Substances 0.000 description 7
- 238000001914 filtration Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 238000012360 testing method Methods 0.000 description 6
- 238000012706 support-vector machine Methods 0.000 description 5
- 230000007774 longterm Effects 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 239000013598 vector Substances 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 229910021389 graphene Inorganic materials 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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
-
- 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)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明涉及一种基于RVM‑KF的滚动轴承剩余寿命预测方法及装置。其中方法包括:获取待预测滚动轴承的实时监测数据,以及与待预测滚动轴承同类的滚动轴承整个生命周期的历史数据;利用历史数据训练RVM模型;根据实时监测数据和SKF确定预测起始时刻,并根据实时监测数据确定卡尔曼滤波器的状态方程;将预测起始时刻输入训练好的RVM模型,得到估计的振动数据,将其作为卡尔曼滤波器的观测值,进行单步预测,将单步预测数据加入历史数据更新RVM模型;利用更新后的RVM模型进行迭代预测,直至单步预测数据中的振动数据超过失效阈值完成剩余寿命的预测。本发明结合RVM模型、卡尔曼滤波器、SKF实现了进入加速退化阶段滚动轴承剩余寿命的准确预测。
Description
技术领域
本发明涉及一种基于RVM-KF的滚动轴承剩余寿命预测方法及装置,属于故障预测与健康管理技术领域。
背景技术
滚动轴承作为一种广泛应用于机械设备的基础零部件,是决定机械设备安全可靠运行的关键部件之一。根据相关的调查研究发现旋转机械设备中大概30%的机械故障是滚动轴承失效的结果。机械故障不仅会造成经济损失,还可能导致人员伤亡,因此准确预测滚动轴承的剩余使用寿命(Remaining useful life,RUL),可以有效地避免因滚动轴承失效所导致的主体机组失效、停机维修、生产力丧失、人员伤亡等严重问题。
目前滚动轴承的剩余寿命预测主要分为基于模型的剩余寿命预测方法和基于数据的剩余寿命预测方法。
基于模型的剩余寿命预测方法通过分析轴承的失效机理或总结经验利用Paris裂纹扩展模型、Forman裂纹扩展定律、指数模型、随机模型等描述轴承的退化过程,继而进行轴承的剩余寿命预测,但滚动轴承的退化机理很难准确得到并且模型的选择对预测结果也有很大的影响。
基于数据的方法不需要建立物理模型或统计模型,而是从现有的监测数据中通过神经网络、SVM(支持向量机)和RVM(相关向量机)等模型得到轴承退化状态与剩余寿命之间的关系,从而预测滚动轴承的剩余寿命,基于数据的方法充分利用了全寿命试验数据,需要大量高质量的数据准确建立轴承退化与剩余寿命之间的关系,但是该方法无法体现个体的差异。
为此,有人提出考虑个体差异和全寿命整体特征的寿命预测方法,例如:期刊为《兵工学报》、期刊号为第39卷第5期2018年5月,名称为“基于支持向量机和卡尔曼滤波的机械零件剩余寿命预测模型研究”,该方法根据现有全寿命试验数据训练所得的SVM回归模型,建立非线性卡尔曼滤波状态更新方程,依据机械零件退化特征构造时间更新方程,设定初始剩余寿命值及其方差,通过逐步迭代计算各时刻剩余寿命估计值及一定置信水平的置信区间。该计算模型能够充分利用现有零件与同类零件全寿命试验数据和被预测零件的实时状态退化数据,实现剩余寿命预测,从而达到所得剩余寿命预测模型既能够利用全寿命数据,又考虑到研究对象个体差异(卡尔曼滤波模型)的目的。
然而上述方法为了得到准确的非线性卡尔曼滤波状态方程,对训练样本质量和SVM模型精度的要求较高,也即对模型和数据的依赖性较大,并且在长期性预测方面的准确性并不高,不能够有效满足工程的需求。
发明内容
本申请的目的在于提供一种基于RVM-KF的滚动轴承剩余寿命预测方法,用以解决现有技术剩余寿命预测不准确的问题;同时还提供一种基于RVM-KF的滚动轴承剩余寿命预测装置,用以解决现有技术剩余寿命预测不准确的问题。
为实现上述目的,本发明提出了一种基于RVM-KF的滚动轴承剩余寿命预测方法的技术方案,包括以下步骤:
1)获取待预测滚动轴承的实时监测数据,以及与待预测滚动轴承同类的滚动轴承整个生命周期的历史数据;所述实时监测数据、历史数据是指振动数据和对应的时间;
2)利用历史数据训练RVM模型;所述RVM模型的输入为时间,输出为振动数据;
3)根据实时监测数据和SKF确定预测起始时刻;并根据实时监测数据确定卡尔曼滤波器的状态方程;
4)将所述预测起始时刻输入训练好的RVM模型,得到估计的振动数据,将估计的振动数据作为卡尔曼滤波器的观测值进行单步预测,得到单步预测数据,若单步预测数据中的振动数据超出失效阈值,则完成剩余寿命的预测,若单步预测数据中的振动数据未超出失效阈值,则将单步预测数据加入步骤2)的历史数据中,重新训练RVM模型,更新RVM模型;迭代的执行本步骤直至单步预测数据中的振动数据超出失效阈值。
另外,本发明还提出了一种基于RVM-KF的滚动轴承剩余寿命预测装置的技术方案,包括处理器、存储器以及存储在所述存储器中并可在处理器上运行的计算机程序,所述处理器在执行所述计算机程序时实现上述基于RVM-KF的滚动轴承剩余寿命预测方法的技术方案。
本发明的基于RVM-KF的滚动轴承剩余寿命预测方法及装置的技术方案的有益效果是:本发明通过同类滚动轴承的历史数据建立了RVM模型,根据待预测滚动轴承的实时监测数据和SKF确定了待预测滚动轴承进入加速退化的时刻,进而确定了预测起始时刻,同时还根据实时监测数据确定了卡尔曼滤波器的状态方程,根据预测起始时刻结合建立的RVM模型对振动数据进行估计,并将估计值作为卡尔曼滤波器的观测值进行单步预测,单步预测的振动数据未超出失效阈值时,将预测数据加入RVM模型中,重新训练RVM模型,进行迭代预测,直至单步预测的振动数据超出失效阈值,表明迭代结束,得到剩余寿命。本发明通过不断更新RVM模型,结合RVM模型、卡尔曼滤波器、SKF实现了进入加速退化阶段滚动轴承剩余寿命的准确预测,减少了滚动轴承剩余寿命长期预测的不确定性,提高了滚动轴承剩余寿命长期预测的精度,而且将卡尔曼滤波与状态方程结合有效减少了预测过程对数据和模型的依赖。
进一步的,上述基于RVM-KF的滚动轴承剩余寿命预测方法及装置中,为了更好的拟合加速退化过程,卡尔曼滤波器的状态方程为随机效应指数模型。
进一步的,上述基于RVM-KF的滚动轴承剩余寿命预测方法及装置中,所述随机效应指数模型为:
其中,Y(t)为t时刻的振动数据;α为常数;θ和β为参数;σ2为方差;ε(t)为误差项,是均值为0、方差为σ2的布朗运动。
进一步的,上述基于RVM-KF的滚动轴承剩余寿命预测方法及装置中,所述振动数据包括从振动信号中提取的至少一个特征值。
进一步的,上述基于RVM-KF的滚动轴承剩余寿命预测方法及装置中,为了提高预测的准确性,并且降低模型训练的复杂程度,训练RVM模型时利用的历史数据为加速退化阶段的历史数据。
附图说明
图1是本发明基于RVM-KF的滚动轴承剩余寿命预测方法的流程图;
图2是本发明滚动轴承的退化曲线;
图3是本发明预测的滚动轴承2-3的退化轨迹和现有技术预测的滚动轴承2-3的退化轨迹的比较图;
图4是本发明预测的滚动轴承2-5的退化轨迹和现有技术预测的滚动轴承2-5的退化轨迹的比较图;
图5是本发明预测的滚动轴承2-3的剩余寿命和现有技术预测的滚动轴承2-3的剩余寿命的比较图;
图6是本发明预测的滚动轴承2-5的剩余寿命和现有技术预测的滚动轴承2-5的剩余寿命的比较图;
图7是本发明基于RVM-KF的滚动轴承剩余寿命预测装置的结构示意图。
具体实施方式
基于RVM-KF的滚动轴承剩余寿命预测方法实施例:
基于RVM-KF的滚动轴承剩余寿命预测方法的主要构思在于,本发明将RVM模型、卡尔曼滤波器以及SKF结合,对滚动轴承进入加速退化阶段的剩余寿命进行迭代预测,提高了寿命预测的准确性。
具体的,基于RVM-KF的滚动轴承剩余寿命预测方法如图1所示,包括以下步骤:
1)收集与待预测滚动轴承同类的滚动轴承整个生命周期的历史数据,历史数据包括振动信号和对应的时间(这里的时间为使用时间)。
为了减小个体差异,对于不同型号的滚动轴承运用不同的模型参数进行寿命预测,因此在建模时选用的历史使用数据对应的滚动轴承,应当与待预测的滚动轴承同类(也即同型号)。
对多个同型号的滚动轴承在相同的试验条件下进行全寿命周期试验,试验过程中周期性的采集滚动轴承工作过程中的振动信号(每1min采集1.28s的振动信号),直至滚动轴承完全失效,获得多个滚动轴承的全寿命振动信号,每个振动信号都有其对应的使用时间。
2)从步骤1)的振动信号中提取能够有效表征滚动轴承退化的特征值。
为了实现剩余寿命的预测,从振动信号中提取表征滚动轴承退化的特征值作为振动数据,然而振动信号中包含多种振动数据,在表征退化的振动数据可以用其他振动数据进行计算得到的情况下,也可以提取其他振动数据,本发明对此不做限制。
3)将步骤2)得到特征值以及每个特征值在步骤1)中对应的时间输入预先建立的SKF模型中,找出同类滚动轴承在生命周期中的加速退化时刻(即早期故障点)。
滚动轴承的退化过程如图2所示包括平稳退化阶段和加速退化阶段,SKF模型为一个动态贝叶斯网络,包括两个卡尔曼滤波器,两个卡尔曼滤波器包括用于预测某个时刻属于平稳退化概率的第一卡尔曼滤波器,适用于平稳退化阶段;以及预测某个时刻属于加速退化概率的第二卡尔曼滤波器,适用于加速退化阶段。SKF所输入的数据通过第一卡尔曼滤波器可以得到每个时刻属于平稳退化阶段的第一概率,通过第二卡尔曼滤波器可以得到每个时刻属于加速退化阶段的第二概率。
具体的卡尔曼滤波器的概率计算过程如下:
第一卡尔曼滤波器和第二卡尔曼滤波器都是基本的尔曼滤波器,其预测和更新过程如下:
预测过程为,状态预测:
协方差预测:
更新过程为,测量残差:
残差协方差:
卡尔曼增益:
状态更新:
协方差更新:
其中,yt为t时刻监测到的特征值,为t时刻估计的特征值,φt为t时刻的状态转移矩阵,/>为φt的转置矩阵,Ht为t时刻的观测矩阵,/>为Ht的转置矩阵,δt为t时刻的测量残差,Ct为t时刻的残差协方差,Ct -1为Ct的逆矩阵,Kt为t时刻的卡尔曼增益,Qt为t时刻的过程噪声的方差,Rt为t时刻测量噪声的方差,/>为协方差的估计值,I为单位矩阵,xt为t时刻更新的估计的特征值,Pt为t时刻更新的协方差的估计值。
接着,第i个卡尔曼滤波器的测量残差服从均值为0的正态分布的概率密度函数为:
得到各个卡尔曼滤波器的概率为:
状态加权
协方差加权
其中,为t时刻下特征值属于第i个卡尔曼滤波器的概率,Zij为第i个卡尔曼滤波器转移到第j个卡尔曼滤波器的概率,n为卡尔曼滤波器的数量,/>为t时刻下第i个卡尔曼滤波器的测量残差,/>为t时刻下第i个卡尔曼滤波器的残差协方差,/>为t时刻第i个卡尔曼滤波器更新的状态,Pt i为t时刻第i个卡尔曼滤波器更新的协方差。将计算的状态加权和协方差加权带入上述卡尔曼滤波过程中,得到各时刻最可能的退化状态。
得到每个时刻的第一概率和第二概率后,若连续若干个时刻(若干个时刻一般选取5个时刻,当然若干个的数量可以根据需要进行设定)的第二概率大于第一概率,则将这若干个时刻的最前端的时刻作为加速退化时刻,当然作为其他实施方式,加速退化时刻也可以根据这若干个时刻进行计算得出,或者将这若干个时刻的最后端的时刻作为加速退化时刻,本发明对此不做限制。
作为其他实施方式,也可以采用其他方法判断待预测滚动轴承所处的阶段,例如:将特征值和经验值比较,也可以判断所处的阶段。本发明对如何判断待预测滚动轴承所处的阶段的方式并不做限制。
4)将步骤3)中得到的加速退化时刻后的历史数据作为训练数据集,训练RVM模型。
RVM模型是贝叶斯框架下的一种广义线性模型,训练数据集为所建立的RVM模型为:
ln=y(tn;ω)+εn;
训练样本的似然函数为:
其中,ln为第n个目标数据(即输出数据),这里为特征值;tn为第n个输入数据,这里为时间(即使用时长);N为训练数据集中数据总数;ωn为第n个数据对应的权重,K(t,tn)为第n个数据对应的核函数;εn为均值为0、方差为σ2的高斯白噪声;Φ是一个基于核函数的N×N+1矩阵,Φ=[φ(t1),φ(t2),…,φ(tn),…,φ(tN)]T,φ(tn)为第n数据对应的核向量,φ(tn)=[1,K(tn,t1),K(tn,t2),…K(tn,tN)]K(tn,t1),K(tn,t2),…,K(tn,tN)均为K(t,tn)中的向量元素,t为输入集;l为输出集;ω为权重集,ω=(ω1,ω2,…,ωn,…ωN)T。
直接使用极大似然法对ω与σ2的求解,结果通常使ω中的元素大部分都不是0,从而导致过拟合。为避免在RVM中出现这个现象,为ω加上先决条件:它们的概率分布是落在0周围的正态分布:
其中,α为一个N+1的矢量;αi为第i个超参数;N()为正态分布。
因此,对于新给定的输入t*,其预测值l*:
p(l*|l)=∫p(l*|ω,α,σ2)p(ω,α,σ2|l)dωdαdσ2。
本步骤中,对于RVM模型的训练优选加速退化时刻后的历史数据,这样不仅可以避免数据量过多,还可以使得RVM模型针对加速退化阶段进行预估,预估的更加准确,作为其他实施方式,也可以将全部历史数据对RVM模型进行训练,本发明对此不做限制。
5)获取待预测滚动轴承的实时监测数据,实时监测数据包括振动信号以及对应的时间(这里为使用时间);从振动信号中提取能够有效表征滚动轴承退化的实时特征值,并根据实时特征值和SKF判断待预测滚动轴承所处的阶段,若判断出待预测滚动轴承进入加速退化时刻,则根据加速退化时刻确定预测起始时刻。
一般情况下预测起始时刻比加速退化时刻稍退后些,例如:加速退化时刻为使用200min对应的时刻,那么预测起始时刻为使用220min对应的时刻,当然,预测起始时刻可以根据需要进行设定,本发明不做限制。
判断待预测滚动轴承所处的阶段的过程与步骤3)相同,只不过这里输入SKF的数据为实时特征值,因此这里不对其具体判定过程进行赘述。
6)根据实时监测数据确定卡尔曼滤波器(KF)中状态方程的参数,确定状态方程,状态方程为随机效应指数模型。
确定状态方程参数的实时监测数据可以为加速退化时刻附近的监测数据,可以为预测起始时刻附近的监测数据,本发明对此不做限制。在未确定预测起始时刻前,状态方程可以根据监测数据的更新而更新状态方程,但是预测起始时刻确定后,状态方程是不变的。
本实施例中,滚动轴承的退化曲线如图2所示,包括平稳的平稳退化状态和加速退化状态,从图中可以看出,加速退化曲线与指数函数曲线类似,因此利用随机效应指数模型表示卡尔曼滤波器的状态方程,随机系数指数模型为:
其中,Y(t)为t时刻的振动数据(这里为滚动轴承的特征值);α为常数;θ和β为参数(也即随机参数);ε(t)为误差项,是均值为0、方差为σ2的布朗运动。
为了方便计算,将随机系数指数模型进行取对数,变形为:
上式中:θ'=lnθ,θ'和β'为变形后的参数,ε(t)表示误差,是一个均值为0,方差为σ2的中心布朗运动。
作为其他实施方式,第三卡尔曼滤波器的状态方程可以是现有的任意一种模型,可以符合加速退化的过程即可,本发明对此不做限制。
7)将步骤5)中确定的预测起始时刻输入步骤4)中训练好的RVM模型,得到估计的待预测滚动轴承的特征值,将估计的待预测滚动轴承的特征值作为卡尔曼滤波器的测量值,进行单步预测,得到单步预测数据,将单步预测数据加入步骤4)中的训练数据集中重新训练RVM模型,更新RVM模型。
单步预测数据中包括单步预测的时间:预测起始时刻+步长,以及单步预测的特征值;
8)利用更新后的RVM模型进行下一时刻的特征值的估计,进行迭代预测,直至单步预测数据中的特征值(即预测特征值)超过失效阈值,完成剩余寿命的预测。
步骤8)中判断单步预测数据中的振动数据是否超过失效阈值,若没超过,则继续迭代预测,若超过则根据迭代次数确定待预测滚动轴承的剩余寿命。
步骤8)中每更新一次RVM模型均是对下一使用时刻的特征值进行估计,例如:确定预测起始时刻为使用220min,根据步骤4)训练好的RVM模型,得到使用220min的估计特征值,通过KF对使用220min的估计特征值进行单步预测,得到单步预测数据;单步预测数据中包括使用220min以及对应的预测特征值;如果使用220min对应的预测特征值未超出失效阈值,则将该单步预测数据加入步骤4)的训练数据集中,重新训练RVM模型,更新RVM模型;
然后根据更新后的RVM模型进行使用220min+1min的特征值的估计,得到使用221min时的估计特征值,通过KF对使用221min的估计特征值进行单步预测,得到单步预测数据;本次迭代的单步预测数据中包括使用221min以及对应的预测特征值;
以此类推,进行迭代预测,若对使用230min的特征值估计后,单步预测数据中的预测特征值超出失效阈值,那么剩余寿命为230min-220min=10min。
本发明中为了得到不同预测起始时刻对应的剩余寿命,在第一个预测起始时刻进行寿命预测完毕后,更新预测起始时刻重新进行寿命预测,此时还需要对状态方程进行更新,也即重复步骤6)-步骤8),完成新的预测起始时刻的剩余寿命的预测。进而得到不同的预测起始时刻对应的剩余寿命。
下面结合实例来验证本发明的有效性和正确性,数据集XJTU-SY是由西安交通大学(XJTU)设计科学与基础构件研究所和浙江长兴昇阳科技股份有限公司(SY)提供。
本发明将滚动轴承的退化过程分为平稳退化和加速退化,提取工况2下的5个滚动轴承的小波包系数均方根(特征值)进行退化状态识别,5个滚动轴承的序号分别为2-1、2-2、2-3、2-4、2-5,各个滚动轴承进入加速退化阶段的时刻如表一。本发明设置滚动轴承的失效阈值为13,并且将序号2-3和序号2-5的滚动轴承作为待预测的滚动轴承,将序号2-3作为待预测的滚动轴承时,序号2-1、2-2、2-4、2-5的滚动轴承的数据作为训练集;将序号2-5作为待预测的滚动轴承时,序号2-1、2-2、2-3、2-4的滚动轴承的数据作为训练集;利用其余4个滚动轴承加速退化阶段小波包系数均方根作为RVM训练集,根据训练好的RVM模型估计待预测滚动轴承的退化特征值,将RVM的估计值作为KF的测量值进行单步预测处理,将KF滤波值加入RVM训练集重新训练RVM模型进行下一次迭代预测,直至达到滚动轴承的失效阈值完成轴承剩余寿命的预测。将利用RVM、KF和本发明迭代方法估计轴承的退化轨迹、剩余寿命进行了对比,对比结果如图3、图4、图5、图6所示。并如表二所示,利用绝对误差MAE和均方根误差RMSE给出了3种方法进行退化轨迹预测的精度。
绝对误差:
均方根误差:
式中,xi为轴承的真实退化特征量,yi为预测的退化特征量数据,n为数据的总个数。
表一5个滚动轴承的加速退化时刻
轴承序号 | 2-1 | 2-2 | 2-3 | 2-4 | 2-5 |
进入加速退化阶段时刻(min) | 456 | 55 | 305 | 31 | 124 |
表二 序号为2-3、2-5的滚动轴承采用不同的的预测方法的绝对误差和均方根误差
从图3和图4的比较结果可以看出,本发明的方法(RVM-KF)的预测结果与真实值更加接近。本发明基于SKF模型,利用分段线性逼近滚动轴承的非线性退化,跟踪滚动轴承退化过程的动态变化,不需要设置阈值和大量数据训练模型,自适应地获得各时刻最可能的退化模型;RVM是贝叶斯框架下的一种广义线性模型,相较于神经网络RVM更适合处理小样本问题,可以提供完整的预测分布;同时将KF与随机效应指数模型相结合可以更好的拟合轴承的退化过程,从而使寿命预测更准确。
本发明的方法不需要精确的模型和大量高质量的数据就能实时准确预测轴承进入加速退化阶段较早期的剩余使用寿命,同时减少了长期预测的不确定性。
基于RVM-KF的滚动轴承剩余寿命预测装置实施例:
基于RVM-KF的滚动轴承剩余寿命预测装置,如图7所示,包括处理器、存储器以及存储在所述存储器中并可在处理器上运行的计算机程序,所述处理器在执行所述计算机程序时实现基于RVM-KF的滚动轴承剩余寿命预测方法。
基于RVM-KF的滚动轴承剩余寿命预测方法的具体实施过程以及效果在上述基于RVM-KF的滚动轴承剩余寿命预测方法实施例中介绍,这里不做赘述。
也就是说,以上基于RVM-KF的滚动轴承剩余寿命预测方法实施例中的方法应理解可由计算机程序指令实现基于RVM-KF的滚动轴承剩余寿命预测方法的流程。可提供这些计算机程序指令到处理器(如通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备等),使得通过处理器执行这些指令产生用于实现上述方法流程所指定的功能。
本实施例所指的处理器是指微处理器MCU或可编程逻辑器件FPGA等的处理装置;
本实施例所指的存储器用于存储实现基于RVM-KF的滚动轴承剩余寿命预测方法而形成的计算机程序指令,包括用于存储信息的物理装置,通常是将信息数字化后再以利用电、磁或者光学等方式的媒体加以存储。例如:利用电能方式存储信息的各式存储器,RAM、ROM等;利用磁能方式存储信息的的各式存储器,硬盘、软盘、磁带、磁芯存储器、磁泡存储器、U盘;利用光学方式存储信息的各式存储器,CD或DVD。当然,还有其他方式的存储器,例如量子存储器、石墨烯存储器等等。
通过上述存储有实现基于RVM-KF的滚动轴承剩余寿命预测方法而形成的计算机程序指令的存储器、处理器构成的基于RVM-KF的滚动轴承剩余寿命预测装置,在计算机中由处理器执行相应的程序指令来实现,计算机可使用windows操作***、linux***、或其他,例如使用android、iOS***程序设计语言在智能终端实现,以及基于量子计算机的处理逻辑实现等。
作为其他实施方式,基于RVM-KF的滚动轴承剩余寿命预测装置还可以包括其他的处理硬件,如数据库或多级缓存、GPU等,本发明并不对基于RVM-KF的滚动轴承剩余寿命预测装置的结构做具体的限定。
Claims (4)
1.一种基于RVM-KF的滚动轴承剩余寿命预测方法,其特征在于,包括以下步骤:
1)获取待预测滚动轴承的实时监测数据,以及与待预测滚动轴承同类的滚动轴承整个生命周期的历史数据;所述实时监测数据、历史数据是指振动数据和对应的时间;
2)利用历史数据训练RVM模型;所述RVM模型的输入为时间,输出为振动数据;
3)根据实时监测数据和SKF确定预测起始时刻;并根据实时监测数据确定卡尔曼滤波器的状态方程,所述卡尔曼滤波器的状态方程为随机效应指数模型,所述随机效应指数模型为:
其中,Y(t)为t时刻的振动数据;α为常数;θ和β为参数;σ2为方差;ε(t)为误差项,是均值为0、方差为σ2的布朗运动;
4)将所述预测起始时刻输入训练好的RVM模型,得到估计的振动数据,将估计的振动数据作为卡尔曼滤波器的观测值,进行单步预测,得到单步预测数据,若单步预测数据中的振动数据超出失效阈值,则完成剩余寿命的预测,若单步预测数据中的振动数据未超出失效阈值,则将单步预测数据加入步骤2)的历史数据中,重新训练RVM模型,更新RVM模型;迭代的执行本步骤直至单步预测数据中的振动数据超出失效阈值。
2.根据权利要求1所述的基于RVM-KF的滚动轴承剩余寿命预测方法,其特征在于,所述振动数据包括从振动信号中提取的至少一个特征值。
3.根据权利要求1所述的基于RVM-KF的滚动轴承剩余寿命预测方法,其特征在于,训练RVM模型时利用的历史数据为加速退化阶段的历史数据。
4.一种基于RVM-KF的滚动轴承剩余寿命预测装置,其特征在于,包括处理器、存储器以及存储在所述存储器中并可在处理器上运行的计算机程序,所述处理器在执行所述计算机程序时实现如权利要求1-3中任一项所述的基于RVM-KF的滚动轴承剩余寿命预测方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010554878.5A CN111783242B (zh) | 2020-06-17 | 2020-06-17 | 一种基于rvm-kf的滚动轴承剩余寿命预测方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010554878.5A CN111783242B (zh) | 2020-06-17 | 2020-06-17 | 一种基于rvm-kf的滚动轴承剩余寿命预测方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111783242A CN111783242A (zh) | 2020-10-16 |
CN111783242B true CN111783242B (zh) | 2024-05-28 |
Family
ID=72757291
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010554878.5A Active CN111783242B (zh) | 2020-06-17 | 2020-06-17 | 一种基于rvm-kf的滚动轴承剩余寿命预测方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111783242B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112949204B (zh) * | 2021-03-22 | 2022-12-09 | 西安交通大学 | 一种数据与模型自适应匹配的滚动轴承剩余寿命预测方法 |
CN113609685B (zh) * | 2021-08-09 | 2023-04-07 | 电子科技大学 | 基于优化rvm和混合退化模型的轴承剩余寿命预测方法 |
CN114662226B (zh) * | 2022-02-24 | 2024-03-29 | 北京工业大学 | 一种基于时变卡尔曼滤波的滚动轴承剩余使用寿命预测方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1705542A1 (en) * | 2005-03-24 | 2006-09-27 | Abb Research Ltd. | Estimating health parameters or symptoms of a degrading system |
CN104318043A (zh) * | 2014-02-20 | 2015-01-28 | 河南科技大学 | 滚动轴承振动性能可靠性变异过程检测方法与装置 |
CN104615866A (zh) * | 2015-01-21 | 2015-05-13 | 北京航空航天大学 | 一种基于物理统计模型的寿命预测方法 |
CN105300692A (zh) * | 2015-08-07 | 2016-02-03 | 浙江工业大学 | 一种基于扩展卡尔曼滤波算法的轴承故障诊断及预测方法 |
CN105956236A (zh) * | 2016-04-22 | 2016-09-21 | 西安交通大学 | 一种双重更新的四因素随机退化模型齿轮寿命预测方法 |
KR20170067292A (ko) * | 2015-12-08 | 2017-06-16 | 한양대학교 산학협력단 | 기계 시스템의 잔여 수명 예측 장치 및 방법 |
CN106934125A (zh) * | 2017-02-28 | 2017-07-07 | 西安交通大学 | 一种梯形噪声分布的指数模型机械设备剩余寿命预测方法 |
CN107145720A (zh) * | 2017-04-19 | 2017-09-08 | 浙江大学 | 连续退化和未知冲击共同作用下的设备剩余寿命预测方法 |
CN108846241A (zh) * | 2018-07-04 | 2018-11-20 | 西安交通大学 | 基于弗雷歇距离动态数模联动的滚动轴承寿命预测方法 |
CN109543317A (zh) * | 2018-04-28 | 2019-03-29 | 北京航空航天大学 | 一种pemfc剩余使用寿命预测的方法及装置 |
CN109726524A (zh) * | 2019-03-01 | 2019-05-07 | 哈尔滨理工大学 | 一种基于cnn和lstm的滚动轴承剩余使用寿命预测方法 |
CN110095731A (zh) * | 2019-05-06 | 2019-08-06 | 哈尔滨工业大学 | 应用于长寿命空间锂离子电池的剩余寿命直接预测方法 |
KR20200049295A (ko) * | 2018-10-31 | 2020-05-08 | 한국전력공사 | 회전기기 고장 예지를 위한 건전성 지표 추이 및 잔존수명 예측 기법 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11262272B2 (en) * | 2018-10-10 | 2022-03-01 | Palo Alto Research Center Incorporated | Adaptive remaining useful life estimation method using constraint convex regression from degradation measurement |
TWI669617B (zh) * | 2018-10-12 | 2019-08-21 | 財團法人工業技術研究院 | 設備健康狀態監控方法及其系統 |
-
2020
- 2020-06-17 CN CN202010554878.5A patent/CN111783242B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1705542A1 (en) * | 2005-03-24 | 2006-09-27 | Abb Research Ltd. | Estimating health parameters or symptoms of a degrading system |
CN104318043A (zh) * | 2014-02-20 | 2015-01-28 | 河南科技大学 | 滚动轴承振动性能可靠性变异过程检测方法与装置 |
CN104615866A (zh) * | 2015-01-21 | 2015-05-13 | 北京航空航天大学 | 一种基于物理统计模型的寿命预测方法 |
CN105300692A (zh) * | 2015-08-07 | 2016-02-03 | 浙江工业大学 | 一种基于扩展卡尔曼滤波算法的轴承故障诊断及预测方法 |
KR20170067292A (ko) * | 2015-12-08 | 2017-06-16 | 한양대학교 산학협력단 | 기계 시스템의 잔여 수명 예측 장치 및 방법 |
CN105956236A (zh) * | 2016-04-22 | 2016-09-21 | 西安交通大学 | 一种双重更新的四因素随机退化模型齿轮寿命预测方法 |
CN106934125A (zh) * | 2017-02-28 | 2017-07-07 | 西安交通大学 | 一种梯形噪声分布的指数模型机械设备剩余寿命预测方法 |
CN107145720A (zh) * | 2017-04-19 | 2017-09-08 | 浙江大学 | 连续退化和未知冲击共同作用下的设备剩余寿命预测方法 |
CN109543317A (zh) * | 2018-04-28 | 2019-03-29 | 北京航空航天大学 | 一种pemfc剩余使用寿命预测的方法及装置 |
CN108846241A (zh) * | 2018-07-04 | 2018-11-20 | 西安交通大学 | 基于弗雷歇距离动态数模联动的滚动轴承寿命预测方法 |
KR20200049295A (ko) * | 2018-10-31 | 2020-05-08 | 한국전력공사 | 회전기기 고장 예지를 위한 건전성 지표 추이 및 잔존수명 예측 기법 |
CN109726524A (zh) * | 2019-03-01 | 2019-05-07 | 哈尔滨理工大学 | 一种基于cnn和lstm的滚动轴承剩余使用寿命预测方法 |
CN110095731A (zh) * | 2019-05-06 | 2019-08-06 | 哈尔滨工业大学 | 应用于长寿命空间锂离子电池的剩余寿命直接预测方法 |
Non-Patent Citations (6)
Title |
---|
New Particle Filter Based on GA for Equipment Remaining Useful Life Prediction;ke li et al.;sensors;20170328;第17卷(第4期);全文 * |
Satellite lithium-ion battery remaining useful life estimation with an iterative updated RVM fused with the KF algorithm;Yuchen SONG et al.;Chinese Journal of Aeronautics;第31卷(第1期);全文 * |
基于UKF的轴承剩余寿命预测方法研究;阙子俊;金晓航;孙毅;;仪器仪表学报(09);全文 * |
基于支持向量机和卡尔曼滤波的机械零件剩余寿命预测模型研究;于震梁;孙志礼;曹汝男;王鹏;;兵工学报;20180531(05);全文 * |
滚动轴承性能保持可靠性预测;夏新涛等;轴承(第6期);全文 * |
航空发动机的健康指标构建与剩余寿命预测;彭开香;皮彦婷;焦瑞华;唐鹏;;控制理论与应用(04);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111783242A (zh) | 2020-10-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111783242B (zh) | 一种基于rvm-kf的滚动轴承剩余寿命预测方法及装置 | |
CN112629863B (zh) | 变工况下动态联合分布对齐网络的轴承故障诊断方法 | |
CN110378052B (zh) | 基于循环神经网络考虑未来工况的设备剩余寿命预测方法 | |
CN108984893B (zh) | 一种基于梯度提升方法的趋势预测方法 | |
Fan et al. | A novel machine learning method based approach for Li-ion battery prognostic and health management | |
CN111414703B (zh) | 一种滚动轴承剩余寿命预测方法及装置 | |
CN113722985B (zh) | 航空发动机健康状态评估和剩余寿命预测方法及*** | |
CN111581892B (zh) | 基于gdau神经网络的轴承剩余寿命的预测方法 | |
Zhao et al. | Bearing health condition prediction using deep belief network | |
CN108961460B (zh) | 基于稀疏esgp与多目标优化的故障预测方法及装置 | |
CN114048688A (zh) | 一种风力发电机轴承寿命预测方法 | |
CN111523081A (zh) | 一种基于增强门控循环神经网络的航空发动机故障诊断方法 | |
CN110610226A (zh) | 一种发电机故障预测方法及装置 | |
CN109145373B (zh) | 基于改进esgp与预测区间的剩余寿命预测方法及装置 | |
CN112347571A (zh) | 考虑模型和数据不确定性的滚动轴承剩余寿命预测方法 | |
CN112069666B (zh) | 一种基于概率潮流法的电网短期可靠性评估方法 | |
CN111523727B (zh) | 基于不确定过程的考虑恢复效应的电池剩余寿命预测方法 | |
CN114166509A (zh) | 一种电机轴承故障预测方法 | |
CN114740388A (zh) | 一种基于改进tcn的锂电池剩余寿命状态评估方法 | |
CN114662386A (zh) | 一种轴承故障诊断方法及*** | |
CN114660497A (zh) | 一种针对容量再生现象的锂离子电池寿命预测方法 | |
CN117521512A (zh) | 一种基于多尺度贝叶斯卷积Transformer模型的轴承剩余使用寿命预测方法 | |
CN114861349A (zh) | 一种基于模型迁移和维纳过程的滚动轴承rul预测方法 | |
CN111400964B (zh) | 一种故障发生时间预测方法及装置 | |
CN117406100A (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 |