CN103488881B - 一种不确定退化测量数据下的设备剩余寿命预测方法 - Google Patents
一种不确定退化测量数据下的设备剩余寿命预测方法 Download PDFInfo
- Publication number
- CN103488881B CN103488881B CN201310409735.5A CN201310409735A CN103488881B CN 103488881 B CN103488881 B CN 103488881B CN 201310409735 A CN201310409735 A CN 201310409735A CN 103488881 B CN103488881 B CN 103488881B
- Authority
- CN
- China
- Prior art keywords
- degradation
- time
- measurement data
- measurement
- state
- 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
- 238000006731 degradation reaction Methods 0.000 title claims abstract description 105
- 230000015556 catabolic process Effects 0.000 title claims abstract description 86
- 238000000034 method Methods 0.000 title claims abstract description 61
- 238000005259 measurement Methods 0.000 claims abstract description 105
- 238000012544 monitoring process Methods 0.000 claims description 12
- 238000009826 distribution Methods 0.000 claims description 10
- 230000005653 Brownian motion process Effects 0.000 claims description 8
- 238000005315 distribution function Methods 0.000 claims description 8
- 238000005537 brownian motion Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 6
- 230000001186 cumulative effect Effects 0.000 claims description 4
- 238000007476 Maximum Likelihood Methods 0.000 claims description 3
- 238000009792 diffusion process Methods 0.000 claims description 3
- 239000000654 additive Substances 0.000 claims description 2
- 230000000996 additive effect Effects 0.000 claims description 2
- 238000012423 maintenance Methods 0.000 abstract description 2
- 238000012360 testing method Methods 0.000 description 4
- 238000013179 statistical model Methods 0.000 description 2
- 230000007850 degeneration Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于可靠性工程技术领域,涉及一种不确定退化测量数据下的设备剩余寿命预测方法。该方法包括以下步骤:建立不确定退化测量数据下的设备性能退化模型;估计模型参数;估计设备潜在性能退化状态;预测剩余寿命。本发明给出了一种不确定退化测量数据下的预测设备剩余寿命的方法,同时考虑了退化状态不确定性和测量不确定性,不仅可以对设备的个体寿命和总体可靠性寿命特征量进行预测分析,还可以作为预测设备剩余寿命的一种有效分析工具,为设备基于状态的维修保障提供有力的理论依据和技术支撑,从而节约经费开支,避免不必要的经济损失,有很好的工程应用价值。
Description
技术领域
本发明属于可靠性工程技术领域,涉及一种不确定退化测量数据下的对高可靠性设备进行剩余寿命预测的方法。
背景技术
现有的剩余寿命预测方法,主要分为两大类:基于模型的方法和基于数据的方法。基于模型的方法包括基于机理模型的方法和基于统计模型的方法。其中,基于机理模型的方法由于***的复杂性和随机性,难以得到***准确的机理模型,使得该方法的应用受到限制:基于统计模型的方法综合利用了***机理知识和状态监测数据,可以充分利用已有信息对剩余寿命进行预测。基于数据的方法充分利用了设备的有效测试数据,不需要建立***性能退化过程的机理模型,然而,由于测量得到的数据会受到各种干扰以及噪声的影响,使测量数据存在一定的不确定性,这将降低对设备剩余寿命的预测的准确性。
在工程实践中,基于设备的测量退化数据进行剩余寿命预测得到了很大的关注.由于***随机退化过程和测量误差的影响,测量数据中不可避免包含退化随机性和测量不确定性.然而,现有基于观测数据的剩余寿命估计研究中,没有将退化随机性和测量不确定性对估计的剩余寿命分布的影响同时考虑,因此基于数据的剩余寿命预测方法在解决实际问题上的应用受到了很大的限制,至今没有见到相关的报道。
发明内容
针对上述现有技术状况,本发明的目的是:提供一种充分利用设备在性能退化过程中得到的测量数据,科学预测高可靠性设备个体寿命和总体寿命特征量,以解决不确定测量数据下设备的剩余寿命预测问题。
现将本发明构思及技术解决方案叙述如下:
本发明的基本构思是,对于潜在性能退化过程可以采用带漂移的维纳模型近似地描述,而不确定测量的影响可以在性能退化状态的基础上通过引入加性噪声的方式考虑;然后利用滤波技术,实现了潜在退化状态的实时估计;在退化状态估计的基础上,实现同时考虑退化状态不确定性和测量不确定性的剩余寿命预测,以提高剩余寿命预测的准确性。
根据上述发明构思,本发明一种不确定退化测量数据下的设备剩余寿命预测方法,其特征在于:包括以下步骤:
步骤1:建立不确定退化测量数据下的设备性能退化模型:
步骤1.1:潜在性能退化过程,将性能退化过程{X(t),t≥0}选取为基于维纳过程的性能退化模型:
X(t)=X(0)+θt+σB(t) (1)
其中,θ是漂移系数,σ(σ>0)是扩散系数;{B(t),t≥0}是标准布朗运动,且有σB(t)□N(0,σ2t),用来刻画退化过程的随机动态及退化过程随时间变化的不确定性.令X(0)=x0表示初始退化状态.不失一般性,令x0=0;
步骤1.2:不确定测量数据,令{Y(t),t≥0}表示对应的测量过程,通过测量过程建立潜在的退化状态与不确定测量数据之间的关系。采用以下退化领域最常用的测量模型,即t时刻的潜在退化状态与测量之间的关系为:
Y(t)=X(t)+ε (2)
其中,ε是随机测量误差,假设在任一时刻t测量误差是独立同分布的,且有ε□N(0,γ2).进一步假设ε与B(t)是相互独立的;
步骤1.3:基于以上的模型描述,一旦随机退化过程{X(t),t≥0}首次达到预先给定的失效阈值ω,则认为设备失效,根据首达时间的概念,设备的寿命T可以定义为
T=inf{t:X(t)≥ω|X(0)<ω}, (3)
其中,ω为预先给定的失效阈值;
步骤1.4:在以上的框架下,基于实时测量的退化数据实现对服役设备剩余寿命的预测,并在新的退化测量数据获取后可以实现剩余寿命分布的更新;假设获取退化测量数据的离散监测时间点为0=t0<t1<...<tk,令yk=Y(tk)表示tk时刻的退化测量.因此,到时刻tk的所有退化测量数据的集合可以表示为Y1∶k={y1,y2,...,yk},对应的退化状态的集合为X1:k={x1,x2,...,xk},其中xk=X(tk).由式(2),进一步可以将tk时刻的测量方程描述为yk=xk+εk,其中εk是ε的独立同分布实现;
步骤2:估计模型参数:
步骤2.1:令Θ=(σ2,γ2,θ)′表示未知参数向量.假设有N个测量设备,且第i个设备的测量时间为t1,...,tM,对应的测量数据为{Yi(tj)=yi,j,i=1,...,N,j=1,...,M},其中表示每个设备的测量数据的个数.由式(2)知,第i个设备的在时刻tj的测量可以表示为
Yi(tj)=θtj+σB(tj)+εij, (4)
其中,εij是测量误差,且有εij□N(0,γ2);
步骤2.2:令t=(t1,...,tM)′,yi=(yi,1,...,yi,M)′.根据独立性假设和布朗运动的独立增量性质,可知yi是多变量正态分布的,其均值和方差特征如下:
yi□N(μ,Ω),μ=θt (5)
Ω=σ2Q+γ2IM,Q=[min{ti,tj}]1≤i,j≤M (6)
其中,IM是M阶单位阵;
步骤2.3:对于第i个设备的测量数据,有
yi□N(θt,σ2Q+γ2IM). (7)
步骤2.4:令Y=(y1,y2,...,yN)′表示N个设备所有的测量数据.根据不同设备测量数据之间的独立性假设,关于对应所有测量数据Y的似然函数为
因此,σ2,γ2,θ的极大似然估计可以通过最大化式(8)得到;
步骤3:估计设备潜在性能退化状态:
步骤3.1:将性能退化状态方程和测量方程在监测时刻转换为离散时间方程.然后在离散时间点tk,k=1,2,...上可以得到变换后的退化模型:
其中vk=σ[B(tk)-B(tk-1)],εk是ε在tk时刻的实现.{vk}k≥1和{εk}k≥1是独立同分布的噪声序列,且有vk□N(0,σ2(tk-tk-1))和εk□N(0,γ2):
步骤3.2:定义和Pk|k=var(xk|Y1:k)分别为通过测量Y1:k对退化状态xk估计的期望和方差;定义和Pk|k-1=var(xk|Y1:k-1)分别为一步预测的期望和方差:
步骤3.3:在时间tk,基于卡尔曼滤波的潜在退化状态估计和更新过程如下:
步骤3.3.1:状态估计:
K(k)=Pk|k-1(Pk|k-1+γ2)-1
Pk|k-1=Pk-1|k-1+σ2(tk-tk-1)
步骤3.3.2:方差更新:
Pk|k=(1-K(k))Pk|k-1.
步骤3.4-应用以上的Kalman滤波算法,基于到tk时刻所有测量Y1:k的状态xk的后验估计是高斯分布的,且可以解析显式表示,即
步骤4:预测剩余寿命:
步骤4.1:通过首达时间,tk时刻的剩余寿命Lk可以定义为
Lk=inf{lk>0:X(lk+tk)≥ω}. (10)
步骤4.2:对应剩余寿命Lk的概率密度函数和累积分布函数分别为和且有
其中,Y1:k是到时刻tk的所有退化测量数据的集合;
步骤4.3:基于前面的模型、测量数据、参数估计、潜在性能退化状态估计,tk时刻预测的剩条寿命的概率密度函数累积分布函数分别及对应的剩余寿命的期望和方差如下:
其中,Φ(·)为标准正态分布的累积分布函数;基于以上过程,即可实现不确定测量数据下预测设备的剩余寿命;
步骤4.4:基于以上过程,即可实现不确定测量数据下预测设备的剩余寿命。
本发明给出了在设备性能退化测量数据不确定的情况下复杂设备的剩余寿命预测方法。可以对设备的个体寿命和总体可靠性寿命特征量进行预测分析,为设备基于状态的维修保障提供有力的理论依据和技术支撑,从而节约经费开支,避免不必要的经济损失,有很好的工程应用价值。
附图说明
图1:本发明实施例陀螺漂移的监测数据曲线图
图2:本发明的方法基于不确定测量数据的潜在退化状态估计曲线图
图3:本发明的方法不同测量点的剩余寿命概率密度曲线图
具体实施方式
现结合附图对本发明的具体实施方式作进一步说明
实施例
下面就通过本实施例测试得到一组陀螺仪漂移监测数据来验证发明提出的方法。该数据包括三个陀螺仪在实际使用过程中记录的在210个不同状态监测时间点上的漂移系数数据。在实验中,按照该陀螺仪的技术指标,选择失效阈值为w=0.36(°/hour)。按照通电时间来记,监测时间间隔为2.5小时。在测试中,一个运行到超过失效阈值,对应的寿命为528小时,其余两个在完成测试后仍满足技术指标的要求。对应陀螺漂移监测全寿命周期的测量数据如图1所示。本实施例不确定退化测量数据下的设备剩余寿命预测方法,包括以下步骤:
步骤1:建立不确定退化测量数据下的设备性能退化模型:
基于观测数据进行剩余寿命预测时,需要将退化随机性和测量不确定性对预测的剩余寿命分布的影响同时考虑,当潜在性能退化过程超过失效阈值ω时,即认为设备失效;
1)潜在性能退化过程:潜在性能退化过程{X(t),t≥0}可以描述为:
X(t)=X(0)+θt+σB(t) (1)
其中,θ是漂移系数:σ(σ>0)是扩散系数;{B(t),t≥0}是标准布朗运动,令X(0)=x0表示初始退化状态。不失一般性,令x0=0;
2)不确定测量数据:令{Y(t),t≥0}表示对应的测量过程,t时刻的潜在退化状态与测量之间的关系为:
Y(t)=X(t)+ε (2)
其中,ε是随机测量误差,且有ε□N(0,γ2)。进一步假设ε与B(t)是相互独立的。
设获取退化测量数据的离散监测时间点为0=t0<t1<...<tk,令yk=Y(tk)表示tk时刻的退化测量。因此,到时刻tk的所有退化测量数据的集合可以表示为Y1:k={y1,y2,...,yk},对应的退化状态的集合为X1:k={x1,x2,...,xk),其中xk=X(tk)。由式(29),进一步可以将tk时刻的测量方程描述为yk=xk+εk,其中εk是ε的独立同分布实现;
步骤2:估计模型参数
令Θ=(σ2,γ2,θ)′表示未知参数向量。假设有N个测量设备,且第i个设备的测量时间为t1,...,tM,对应的测量数据为{Yi(tj)=yi,j,i=1,...,N,j=1,...,M},其中表示每个设备的测量数据的个数。由式(29)知,第i个设备的在时刻tj的测量可以表示为
Yi(tj)=θtj+σB(tj)+εij, (4)
其中,εij是测量误差,且有εij□N(0,γ2);
令t=(t1,...,tM)′,yi=(yi,1,...,yi,M)′。根据独立性假设和布朗运动的独立增量性质,可知yi是多变量正态分布的,其均值和方差特征如下:
yi□N(μ,Ω),μ=θt (5)
Ω=σ2Q+γ2IM,Q=[min{ti,tj}]1≤i,j≤M (6)
其中,IM是M阶单位阵:
因此,对于第i个设备的测量数据,有
yt□N(θt,σ2Q+γ2IM). (7)
令Y=(y1,y2,...,yN)′表示N个设备所有的测量数据。根据不同设备测量数据之间的独立性假设,关于Θ对应所有测量数据Y的似然函数为
因此,σ2,γ2,θ的极大似然估计可以通过最大化式(34)得到。并利用AIC准则和对数似然函数值进行模型比较。表1给出了在考虑不确定测量和不考虑不确定测量时的参数估计结果。
表1.退化模型参数估计结果
通过表1可以看出,考虑不确定测量时模型对数据的拟合结果较不考虑不确定测量的结果有较大的提高,因为考虑测量不确定性的结果对应更大的似然函数值和更小的AIC值。
步骤3:估计设备潜在性能退化状态
为了估计设备的潜在性能退化状态,将性能退化状态方程和测量方程在监测时刻转换为离散时间方程.然后在离散时间点tk,k=1,2,...上可以得到变换后的退化模型:
xk=xk-1+θ(tk-tk-1)+vk (9)
yk=xk+εk
其中vk=σ[B(tk)-B(tk-1)],且有vk□N(0,σ2(tk-tk-1))和εk□N(0,γ2).
根据模型(22),可以利用卡尔曼滤波技术实现潜在退化状态的估计.;首先,定义和Pk|k=var(xk|Y1:k)分别为通过测量Y1:k对退化状态xk估计的期望和方差.此外,定义和Pk|k-1=var(xk|Y1:k-1)分别为一步预测的期望和方差。在时间tk,基于卡尔曼滤波的潜在退化状态估计和更新过程如下:
状态估计:
K(k)=Pk|k-1(Pk|k-1+γ2)-1
Pk|k-1=Pk-1|k-1+σ2(tk-tk-1)
方差更新:
Pk|k=(1-K(k))Pk|k-1.
应用以上的Kalman滤波算法,基于到tk时刻所有测量Y1:k的状态xk的后验估计是高斯分布的,且可以解析显式表示,即
基于表1给出的模型参数,利用图1中的全周期数据就可以实现本文提出剩余寿命估计方法。首先给出本文方法基于测量数据估计的潜在退化状态的结果,如图2所示。
步骤4:预测剩余寿命
通过首达时间,tk时刻的剩余寿命Lk可以定义为
Lk=inf{lk>0:X(lk+tk)≥ω}. (10)
对应剩余寿命Lk的概率密度函数和累积分布函数分别为和且有
其中,Y1:k是到时刻tk的所有退化测量数据的集合;
基于前面的模型、测量数据、参数估计、潜在性能退化状态估计,tk时刻预测的剩余寿命的概率密度函数累积分布函数分别及对应的剩余寿命的期望和方差如下:
估计出式(25)-(27)中的参数以后,就能够得到不确定测量数据下的设备剩余寿命分布。由式(25)-(27)可以看出,预测的剩余寿命同时考虑了退化过程本身的随机性和测量不确定性引起的潜在退化状态估计的不确定性,且在γ2=0时,仅考虑退化过程本身的随机性的结果是以上结果的特例。
在退化状态估计的基础上,可以通过式(25)计算在不同测量数据点时的剩余寿命的概率密度函数。图3给出了几个不同测量点对应的剩余寿命估计的概率密度曲线。
从图2和图3可以看出本发明方法可以实现基于不确定测量的潜在退化状态和剩余寿命分布的估计。为了验证本发明方法可以提高剩余寿命估计的准确性,在寿命的30%、50%、70%、90%点上分别给出了考虑不确定测量和不考虑不确定测量时估计的平均剩余寿命以及对应的估计相对误差,如表2所示。
表2.在不同寿命分位点估计的平均剩余寿命和相对误差
通过表2中的结果可以看出,本文提出的方法可以提高剩余寿命估计的准确性,而且考虑不确定测量的影响可以减少剩余寿命估计的相对误差,进而验证了本发明方法的有效性。
Claims (6)
1.一种不确定退化测量数据下的设备剩余寿命预测方法,其特征在于:对于潜在性能退化过程用带漂移的维纳模型近似地描述,不确定测量的影响在性能退化状态的基础上通过引入加性噪声的方式考虑;利用滤波技术,实现潜在退化状态的实时估计;在退化状态估计的基础上,实现退化状态不确定性和测量不确定性的剩余寿命预测,提高剩余寿命预测的准确性;包括以下步骤:
步骤1:建立不确定退化测量数据下的设备性能退化模型,
步骤2:估计模型参数,
步骤3:估计设备潜在性能退化状态,
步骤4:预测剩余寿命。
2.根据权利要求1所述的一种不确定退化测量数据下的设备剩余寿命预测方法,其特征在于:步骤1中所述的“建立不确定退化测量数据下的设备性能退化模型”的具体步骤如下:
步骤1.1:潜在性能退化过程,将性能退化过程{X(t),t≥0}选取为基于维纳过程的性能退化模型:
X(t)=X(0)+θt+σB(t) (1)
其中,θ是漂移系数,σ是扩散系数;{B(t),t≥0}是标准布朗运动,且有σB(t)~N(0,σ2t),用来刻画退化过程的随机动态及退化过程随时间变化的不确定性,令X(0)=x0表示初始退化状态,不失一般性,令x0=0;
步骤1.2:不确定测量数据,令{Y(t),t≥0}表示对应的测量过程,通过测量过程建立潜在的退化状态与不确定测量数据之间的关系,采用以下退化领域最常用的测量模型,即t时刻的潜在退化状态与测量之间的关系为:
Y(t)=X(t)+ε (2)
其中,ε是随机测量误差,假设在任一时刻t测量误差是独立同分布的,且有ε~N(0,γ2),进一步假设ε与B(t)是相互独立的;
步骤1.3:基于以上的模型描述,一旦随机退化过程{X(t),t≥0}首次达到预先给定的失效阈值ω,则认为设备失效,根据首达时间的概念,设备的寿命T定义为
T=inf{t:X(t)≥ω|X(0)<ω}, (3)
其中,ω为预先给定的失效阈值;
步骤1.4:在以上的框架下,基于实时测量的退化数据实现对服役设备剩余寿命的预测,并在新的退化测量数据获取后实现剩余寿命分布的更新;假设获取退化测量数据的离散监测时间点为0=t0<t1<...<tk,令yk=Y(tk)表示tk时刻的退化测量,因此,到时刻tk的所有退化测量数据的集合表示为Y1:k={y1,y2,...,yk},对应的退化状态的集合为X1:k={x1,x2,...,xk},其中xk=X(tk),由式(2),进一步将tk时刻的测量方程描述为yk=xk+εk,其中εk是ε的独立同分布实现。
3.根据权利要求1所述的一种不确定退化测量数据下的设备剩余寿命预测方法,其特征在于:步骤2中所述的“估计模型参数”的具体步骤如下:
步骤2.1:令Θ=(σ2,γ2,θ)′表示未知参数向量.假设有N个测量设备,且第i个设备的测量时间为t1,...,tM,对应的测量数据为{Yi(tj)=yi,j,i=1,...,N,j=1,...,M},其中表示每个设备的测量数据的个数.由式(2)知,第i个设备的在时刻tj的测量表示为
Yi(tj)=θtj+σB(tj)+εij, (4)
其中,εij是测量误差,且有εij~N(0,γ2);
步骤2.2:令t=(t1,...,tM)′,yi=(yi,1,...,yi,M)′.根据独立性假设和布朗运动的独立增量性质,可知yi是多变量正态分布的,其均值和方差特征如下:
yi~N(μ,Ω),μ=θt (5)
Ω=σ2Q+γ2IM,Q=[min{ti,tj}]1≤i,j≤M (6)
其中,IM是M阶单位阵;
步骤2.3:对于第i个设备的测量数据,有
yi~N(θt,σ2Q+γ2IM) (7)
步骤2.4:令Y=(y1,y2,...,yN)′表示N个设备所有的测量数据.根据不同设备测量数据之间的独立性假设,关于Θ对应所有测量数据Y的似然函数为
因此,σ2,γ2,θ的极大似然估计通过最大化式(8)得到。
4.根据权利要求1所述的一种不确定退化测量数据下的设备剩余寿命预测方法,其特征在于:步骤3中所述的“估计设备潜在性能退化状态”的具体步骤如下:
步骤3.1:将性能退化状态方程和测量方程在监测时刻转换为离散时间方程,然后在离散时间点tk,k=1,2,...,k为自然数,得到变换后的退化模型:
其中vk=σ[B(tk)-B(tk-1)],εk是ε在tk时刻的实现,{vk}k≥1和{εk}k≥1是独立同分布的噪声序列,且有vk~N(0,σ2(tk-tk-1))和εk~N(0,γ2);
步骤3.2:定义和Pk|k=var(xk|Y1:k)分别为通过测量Y1:k对退化状态xk估计的期望和方差;定义和Pk|k-1=var(xk|Y1:k-1)分别为一步预测的期望和方差;
步骤3.3:在时间tk,基于卡尔曼滤波的潜在退化状态估计和更新过程:
步骤3.4:应用以上的卡尔曼滤波算法,基于到tk时刻所有测量Y1:k的状态xk的后验估计是高斯分布的,且解析显式表示,即
5.根据权利要求4所述的一种不确定退化测量数据下的设备剩余寿命预测方法,其特征在于:步骤3.3:中所述的“在时间tk,基于卡尔曼滤波的潜在退化状态估计和更新过程”的具体步骤如下:
步骤3.3.1:状态估计:
K(k)=Pk|k-1(Pk|k-1+γ2)-1
Pk|k-1=Pk-1|k-1+σ2(tk-tk-1)
步骤3.3.2:方差更新:
Pk|k=(1-K(k))Pk|k-1。
6.根据权利要求1所述的一种不确定退化测量数据下的设备剩余寿命预测方法,其特征在于:步骤4中所述的“预测剩余寿命”的具体步骤如下:
步骤4.1:通过首达时间,tk时刻的剩余寿命Lk定义为
Lk=inf{lk>0:X(lk+tk)≥ω} (10)
步骤4.2:对应剩余寿命Lk的概率密度函数和累积分布函数分别为和且有
其中,Y1:k是到时刻tk的所有退化测量数据的集合;
步骤4.3:基于前面的模型、测量数据、参数估计、潜在性能退化状态估计,tk时刻预测的剩余寿命的概率密度函数累积分布函数及对应的剩余寿命的期望和方差如下:
其中,Φ(·)为标准正态分布的累积分布函数;基于以上过程,即可实现不确定测量数据下预测设备的剩余寿命;
步骤4.4:基于以上过程,即可实现不确定测量数据下预测设备的剩余寿命。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310409735.5A CN103488881B (zh) | 2013-09-06 | 2013-09-06 | 一种不确定退化测量数据下的设备剩余寿命预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310409735.5A CN103488881B (zh) | 2013-09-06 | 2013-09-06 | 一种不确定退化测量数据下的设备剩余寿命预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103488881A CN103488881A (zh) | 2014-01-01 |
CN103488881B true CN103488881B (zh) | 2017-05-17 |
Family
ID=49829097
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310409735.5A Active CN103488881B (zh) | 2013-09-06 | 2013-09-06 | 一种不确定退化测量数据下的设备剩余寿命预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103488881B (zh) |
Families Citing this family (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104200091B (zh) * | 2014-08-27 | 2017-02-01 | 北京航空航天大学 | 一种基于寿命模型分散性的柱塞泵加速因子区间计算方法 |
CN104573881B (zh) * | 2015-02-10 | 2018-01-09 | 广东石油化工学院 | 一种基于退化数据建模的服役设备剩余寿命自适应预测方法 |
CN104850736A (zh) * | 2015-04-27 | 2015-08-19 | 大连理工大学 | 一种基于状态空间模型的高速数控铣床刀具寿命预测方法 |
TWI570587B (zh) * | 2015-12-07 | 2017-02-11 | 財團法人工業技術研究院 | 半導體機台零件剩餘壽命預測系統與方法 |
CN105354615A (zh) * | 2015-12-10 | 2016-02-24 | 青岛农业大学 | 防喷阀的可变成本实时预测维护与备件订购方法及*** |
CN105701350B (zh) * | 2016-01-13 | 2018-09-18 | 安徽理工大学 | 间歇性工作开关电源剩余寿命预测方法 |
CN105868557A (zh) * | 2016-03-29 | 2016-08-17 | 浙江大学 | 一种两阶段退化情况下的机电设备剩余寿命在线预测方法 |
CN106484949B (zh) * | 2016-09-12 | 2019-08-16 | 西安理工大学 | 基于退化数据的动量轮可靠性分析与剩余寿命预测方法 |
CN106353691A (zh) * | 2016-10-31 | 2017-01-25 | 首都师范大学 | 一种计及自愈影响的锂电池退化建模及寿命预测方法 |
CN106526492A (zh) * | 2016-10-31 | 2017-03-22 | 首都师范大学 | 计及时变环境温度影响的锂电池退化建模及寿命预测方法 |
CN106291391A (zh) * | 2016-10-31 | 2017-01-04 | 首都师范大学 | 一种计及随机变电流的锂电池退化建模及寿命预测方法 |
FR3061307B1 (fr) * | 2016-12-22 | 2021-05-07 | Electricite De France | Caracterisation perfectionnee d'un dispositif electrochimique en operation pour un pronostic de fonctionnement futur du dispositif |
CN107451392B (zh) * | 2017-06-23 | 2020-01-31 | 山东科技大学 | 一种含有多个相关退化过程的剩余寿命预测方法 |
CN107480442B (zh) * | 2017-08-08 | 2020-01-14 | 山东科技大学 | 依赖于时间和状态的长程相关退化过程剩余寿命预测方法 |
CN108304685A (zh) * | 2018-04-13 | 2018-07-20 | 中国人民解放***箭军工程大学 | 一种非线性退化设备剩余寿命预测方法及*** |
CN110654949B (zh) * | 2018-06-29 | 2021-02-19 | 哈尔滨工业大学 | 一种维保情况下的电梯安全剩余使用寿命确定方法 |
CN109387779B (zh) * | 2018-09-28 | 2020-08-04 | 河北工业大学 | 一种基于统计数据驱动的万能式断路器操作附件剩余寿命预测方法 |
CN109101778B (zh) * | 2018-11-22 | 2019-04-02 | 中国人民解放军国防科技大学 | 基于性能退化数据和寿命数据融合的Wiener过程参数估计方法 |
CN109992875B (zh) * | 2019-03-28 | 2020-11-17 | 中国人民解放***箭军工程大学 | 一种切换设备剩余寿命的确定方法及*** |
EP3959571A1 (de) * | 2019-04-23 | 2022-03-02 | Volkswagen Aktiengesellschaft | Verfahren zum bestimmen von restnutzungszyklen, restnutzungszyklusbestimmungsschaltung, restnutzungszyklusbestimmungsvorrichtung |
CN110705139A (zh) * | 2019-08-27 | 2020-01-17 | 华东光电集成器件研究所 | 一种基于多应力耦合下电子设备可靠性评估方法 |
CN110569623B (zh) * | 2019-09-19 | 2023-03-24 | 西安工程大学 | 一种细纱机电机使用寿命预测方法 |
CN111967142B (zh) * | 2020-07-20 | 2023-07-07 | 北京航空航天大学 | 一种考虑认知不确定性的性能退化实验建模与分析方法 |
CN112214934B (zh) * | 2020-12-08 | 2021-03-02 | 中国人民解放军国防科技大学 | 一种基于多传感器的悬浮***寿命预测方法及相关装置 |
CN112949057B (zh) * | 2021-02-26 | 2023-09-22 | 中国人民解放***箭军工程大学 | 一种融合失效寿命数据的设备剩余寿命预测方法 |
CN113011036B (zh) * | 2021-03-26 | 2023-03-21 | 中国人民解放***箭军工程大学 | 一种针对线性维纳退化过程模型参数的无偏估计方法 |
CN113033015B (zh) * | 2021-04-09 | 2024-05-14 | 中国人民解放***箭军工程大学 | 考虑两阶段自适应Wiener过程的退化设备剩余寿命预测方法 |
CN113065675B (zh) * | 2021-04-13 | 2023-07-18 | 中国人民解放军空军工程大学 | 一种基于剩余寿命预测的设备最优维护方法 |
CN113158471B (zh) * | 2021-04-26 | 2023-01-24 | 中国人民解放***箭军工程大学 | 一种考虑测量不确定性的退化设备剩余寿命预测方法 |
CN113392507A (zh) * | 2021-05-25 | 2021-09-14 | 西安科技大学 | 一种锂离子动力电池剩余寿命预测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101870076A (zh) * | 2010-07-02 | 2010-10-27 | 西南交通大学 | 一种基于性能退化模型的数控机床导轨副寿命预测方法 |
CN102789545A (zh) * | 2012-07-12 | 2012-11-21 | 哈尔滨工业大学 | 基于退化模型匹配的涡轮发动机剩余寿命的预测方法 |
-
2013
- 2013-09-06 CN CN201310409735.5A patent/CN103488881B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101870076A (zh) * | 2010-07-02 | 2010-10-27 | 西南交通大学 | 一种基于性能退化模型的数控机床导轨副寿命预测方法 |
CN102789545A (zh) * | 2012-07-12 | 2012-11-21 | 哈尔滨工业大学 | 基于退化模型匹配的涡轮发动机剩余寿命的预测方法 |
Non-Patent Citations (4)
Title |
---|
A comparison of models for measurable deterioration:An application to coatings on steel structures;Robin P.Nicolai et al.;《reliability engineering & system safety》;20061120;第1635-1650页 * |
A Wiener-process-based degradation model with a recursive filter algorithm for remaining useful life estimation;xiaosheng si et al.;《Mechanical Systems and Signal Processing》;20120901;第219-237页 * |
一种基于EMD去噪的维纳过程剩余寿命预测方法;张建勋等;《proceedings of the 32nd chinese control conference》;20130726;摘要、第6113页第3节-第6114页左栏 * |
带测量误差的非线性退化过程建模与剩余寿命估计;司小胜等;《自动化学报》;20130531;第39卷(第5期);第530-541页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103488881A (zh) | 2014-01-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103488881B (zh) | 一种不确定退化测量数据下的设备剩余寿命预测方法 | |
Ye et al. | A Bayesian approach to condition monitoring with imperfect inspections | |
Yu | State-of-health monitoring and prediction of lithium-ion battery using probabilistic indication and state-space model | |
CN109829137B (zh) | 一种周期应力下非线性退化设备的寿命预测方法及*** | |
US9658291B1 (en) | Methods and apparatus for dynamic estimation of battery open-circuit voltage | |
Feng et al. | Remaining useful lifetime prediction based on the damage-marker bivariate degradation model: A case study on lithium-ion batteries used in electric vehicles | |
CN112949026B (zh) | 一种考虑年龄和状态依赖的退化设备剩余寿命预测方法 | |
Orchard et al. | A particle filtering-based framework for real-time fault diagnosis and failure prognosis in a turbine engine | |
JP2019537923A (ja) | 電池の健全性状態の判定および警報の発生 | |
CN107515965B (zh) | 一种基于不确定过程的加速退化建模评估方法 | |
CN110851980A (zh) | 一种设备剩余寿命预测方法及*** | |
CN109977552B (zh) | 一种考虑状态检测影响的设备剩余寿命预测方法及*** | |
Tang et al. | Novel metrics and methodologies for the verification and validation of prognostic algorithms | |
Edwards et al. | Impact of input uncertainty on failure prognostic algorithms: Extending the remaining useful life of nonlinear systems | |
CN113158471B (zh) | 一种考虑测量不确定性的退化设备剩余寿命预测方法 | |
US12001180B2 (en) | Condition-based method for malfunction prediction | |
CN111460692A (zh) | 考虑退化速率相互影响的设备剩余寿命预测方法及*** | |
Sankararaman et al. | Uncertainty quantification in remaining useful life of aerospace components using state space models and inverse form | |
Daigle et al. | Advanced methods for determining prediction uncertainty in model-based prognostics with application to planetary rovers | |
CN112100574A (zh) | 一种基于重采样的aakr模型不确定度计算方法及*** | |
Kumar et al. | Power system dynamic state estimation using kalman filtering technique | |
Xu et al. | An adaptive gamma process based model for residual useful life prediction | |
CN115884895A (zh) | 用于预报电动车辆车队的车辆电池的剩余使用寿命的方法 | |
CN112949057B (zh) | 一种融合失效寿命数据的设备剩余寿命预测方法 | |
Song et al. | A sliding sequence importance resample filtering method for rolling bearings remaining useful life prediction based on two Wiener-process models |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |