CN107688687A - 一种考虑长程相关性和部件不确定性的寿命预测方法 - Google Patents

一种考虑长程相关性和部件不确定性的寿命预测方法 Download PDF

Info

Publication number
CN107688687A
CN107688687A CN201710554292.7A CN201710554292A CN107688687A CN 107688687 A CN107688687 A CN 107688687A CN 201710554292 A CN201710554292 A CN 201710554292A CN 107688687 A CN107688687 A CN 107688687A
Authority
CN
China
Prior art keywords
mrow
msub
msup
msubsup
mfrac
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.)
Pending
Application number
CN201710554292.7A
Other languages
English (en)
Inventor
周东华
陈茂银
席霄鹏
卢晓
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University of Science and Technology
Original Assignee
Shandong University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shandong University of Science and Technology filed Critical Shandong University of Science and Technology
Priority to CN201710554292.7A priority Critical patent/CN107688687A/zh
Publication of CN107688687A publication Critical patent/CN107688687A/zh
Pending legal-status Critical Current

Links

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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M15/00—Testing of engines
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/392—Determining battery ageing or deterioration, e.g. state of health

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

本发明公开了一种考虑长程相关性和部件不确定性的寿命预测方法，属于预测与健康管理中剩余寿命预测领域，包括以下步骤：读入状态监测退化数据，针对指定漂移结构(如指数函数、幂函数以及线性函数)初始化模型参数；利用极大似然算法估计模型的未知参数；基于弱收敛准则对所得退化过程进行近似，并计算剩余寿命的概率密度函数；最后，利用AIC准则和均方误差判断模型的拟合效果及预测能力。本方法较以往方法具有更高的普适性和预测精度。

Description

一种考虑长程相关性和部件不确定性的寿命预测方法
技术领域
本发明属于预测与健康管理中剩余寿命预测领域，具体涉及一种考虑长程相关性和部件不确定性的寿命预测方法。
背景技术
针对实际工业系统，剩余寿命预测是有效降低维护费用的一项重要技术。基于状态监测数据，剩余寿命预测技术的核心是估计退化过程距离首次到达预定失效阈值的时间的概率密度函数，高精度的剩余寿命预测方法是预测维护领域的理论基础。
目前，预测与健康管理领域国内外学者提出的剩余寿命预测算法按照信息的获取方式可分为基于直接监测数据和基于间接监测数据两类。其中，基于直接监测数据的方法简单有效，可信度相对更高。因此，很多学者提出了基于直接监测数据的退化建模和剩余寿命预测算法，如利用多特征融合技术对马尔科夫模型进行剩余寿命预测(Yan,2011)、利用EM算法和贝叶斯平滑器对非线性维纳过程进行剩余寿命预测(Huang,2015)、结合维纳过程和主成分回归利用拉姆波传感器数据进行剩余寿命预测(Mishra,2016)。
但是，实际工业系统通常会与其所处环境发生交互作用，致使过去的退化状态能够影响未来的退化轨迹。对于这种情况，无记忆性的马尔科夫过程将无法描述此类退化过程。一个典型的例子是大型高炉的生产过程，该过程包含一系列复杂的物理和化学反应，这些反应最直观的体现就在于炉内温度的变化。调研发现，炉温曲线通常呈现某种确定性的趋势，而不是简单的随机游走，这意味着退化状态之间可能具有长程相关性。
近几年利用带漂移项的布朗运动进行退化建模的方法取得了较为广泛的应用(Si,2012；Benkedjouh,2013；Valivs,2015；Li,2015)。但是，布朗运动实际上隶属于马尔科夫过程，即完全忽略了记忆效应。仍然以高炉炉温为例，炉内复杂的热传导大大地影响了温度的测量，使得温度退化非平稳且存在相关性。
显然，已有的基于监测数据的退化建模和剩余寿命预测算法无法描述含记忆效应的非马尔科夫过程。此外，考虑一些实际因素，比如环境与操作条件的影响，各个设备之间可能会产生一定的差异性。为了解决此问题，本发明搭建了一种新的基于分数布朗运动的非线性退化模型，同时引入一个随机的漂移系数来描述不同部件之间的不确定性。由于分数布朗运动具有长程相关性和自相似性的分形特征，能够用于描述一类非平稳的非马尔科夫过程。本发明提供的剩余寿命预测算法相比传统方法适应性更强，预测精度相对更高。
发明内容
本发明的目的在于建立一个更为普适的退化模型，同时提高剩余寿命的预测精度，以便其能够处理工业退化过程中存在的长程相关性和部件不确定性，并为后续的维护工作提供合理的建议。
为了实现上述目的，本发明采用如下技术方案：
一种考虑长程相关性和部件不确定性的寿命预测方法，在处理退化数据时是按照如下步骤依次进行的：
步骤1：读入M组状态监测退化数据，记为X＝(X1,X2,…,XM)；
步骤2：对退化数据进行零初值化处理，针对如下模型初始化参数，包括漂移系数的均值μα和标准差σα、非线性系数β、扩散系数σ以及Hurst指数H；
x(t)＝αtβ+σBH(t) (1)；
其中，x(t)是t时刻下的退化状态，α是服从高斯分布的漂移系数，β是固定的非线性系数，σ是扩散系数，BH(t)表示标准分数布朗运动；
步骤3：根据零初值化处理后的退化数据列写对数似然函数，如公式(2)所示；
其中，Xm＝(xm(t1),xm(t2),…,xm(tN))T,m＝1,2,…,M，Ψ＝σ2Q，Qi,j＝cov(BH(ti),且有
分别对μα和σα求偏导并令导数等于零，得到μα和σα的解析表达式，如下所示：
将(5)式和(6)式代入(2)式，利用多维搜索方法求出β,σ,H的极大似然估计，再将β,σ,H代回μα和σα的表达式，即(5)式和(6)式，求得μα和σα的极大似然估计；
步骤4：基于弱收敛准则将退化过程近似为马尔科夫过程，计算剩余寿命lm,k的概率密度函数fm,k(lm,k)，如公式(7)所示：：
其中，
其中，
Δl表示外推的采样间隔，ωm,k＝ωm-xm(tk)；
步骤5：根据公式(12)、(13)计算模型的AIC(Akaike information criterion，赤池信息量准则)和均方误差MSE，即
其中，p为模型未知参数的个数，M为系统部件的总数，N为退化数据的长度，为第m个部件在tk时刻下真实的剩余寿命，输出fm,k(lm,k)作为最终的剩余寿命预测结果。
本发明所带来的有益技术效果：
本发明所述方法能够处理含长程相关性和部件不确定性的退化过程，剩余寿命的预测精度大大提高；本发明所述方法普适性更强，准确度有很大提高，原因在于新的退化模型将布朗运动升级为分数布朗运动，考虑了长程相关性，具体来说：当0.5＜H≤1时，未来状态更趋向于沿着当前方向发展；当0≤H＜0.5时，未来的退化更趋向于均值回归过程；当H＝0.5时，分数布朗运动退化为无记忆性的一般布朗运动；此外，对于存在部件不确定性的退化过程，本发明仍能较为准确的预测系统的剩余寿命。
本发明可以应用到以下场合：
i高速列车、航空航天设备等的性能退化过程；
ii各种化工生产过程。
附图说明
图1为本发明方法的流程图。
图2为锂电池容量原始退化轨迹示意图。
图3为修正后的锂电池容量退化轨迹示意图。
图4为锂电池剩余寿命预测结果示意图。
图5为涡扇发动机测量信号示意图。
图6为修正后的涡扇发动机退化轨迹示意图。
图7为涡扇发动机剩余寿命预测结果示意图。
具体实施方式
下面结合附图以及具体实施方式对本发明作进一步详细说明：
1、如图1所示，一种考虑长程相关性和部件不确定性的寿命预测方法，具体包括如下步骤：
步骤1：读入M组状态监测退化数据，记为X＝(X1,X2,…,XM)；
步骤2：对退化数据进行零初值化处理，针对如下模型初始化参数，包括漂移系数的均值μα和标准差σα，非线性系数β，扩散系数σ，以及Hurst指数H；
x(t)＝αtβ+σBH(t) (1)；
其中，x(t)是t时刻下的退化状态，α是服从高斯分布的漂移系数，β是固定的非线性系数，σ是扩散系数，BH(t)表示标准分数布朗运动；
步骤3：根据零初值化处理后的退化数据列写对数似然函数如下
其中，Xm＝(xm(t1),xm(t2),…,xm(tN))T,m＝1,2,…,M，Ψ＝σ2Q，Qi,j＝cov(BH(ti),且有
分别对μα和σα求偏导并令导数等于零，得到μα和σα的解析表达式，
将式(5)和式(6)代入式(2)，得到如下表达式：
利用MATLAB提供的多维搜索方法fminsearch函数求出β,σ,H的极大似然估计，再将β,σ,H代回式(5)和式(6)，求得μα和σα的极大似然估计；
步骤4：基于弱收敛准则将退化过程近似为马尔科夫过程，计算剩余寿命lm,k的概率密度函数fm,k(lm,k)：
其中，
其中，
Δl表示外推的采样间隔，ωm,k＝ωm-xm(tk)；
步骤5：计算模型的赤池信息量准则(Akaike information criterion，AIC)和均方误差MSE，即
其中，p为模型未知参数的个数，M为系统部件的总数，N为退化数据的长度，为第m个部件在tk时刻下真实的剩余寿命，输出fm,k(lm,k)作为最终的剩余寿命预测结果。
2、实例寿命预测结果
下面结合实例对所述寿命预测方法的具体实施步骤做进一步的具体说明。
下述实例的仿真环境如下：
机型：Intel Core i7-4790(CPU 3.60Ghz，8.00GB RAM)；
操作系统：Windows 10；
软件：Matlab R2016a。
a、实例一：以锂电池退化数据为例，给出剩余寿命预测的仿真步骤说明以及预测结果：
1)读入四组锂电池容量退化数据，如图2所示；
2)为简便起见，我们首先将时间数据转化为循环周期数据。注意到锂电池的退化存在一个弛豫效应，也就是说，电池的电量可能会在一个循环周期过后的存储阶段发生回升的现象，该现象会对剩余寿命预测造成巨大误差。因此，这里采用线性等分的方法对原始数据进行预处理，以消除所谓的弛豫效应，最终得到的退化数据如图3所示。对退化数据进行零初值化处理，并利用式(1)的模型结构初始化参数{μα,σα,β,σ,H}；
3)利用式(2)-式(6)求出μα,σα,β,σ,H的极大似然估计，参数估计结果为：μa＝-0.0072，σa＝0.0012，b＝0.9007，σB＝0.0042；
4)利用式(7)-式(11)求取剩余寿命lm,k的概率密度函数fm,k(lm,k)，以电池B0005为例，得到的剩余寿命预测结果如图4所示，可以看出，在各个监测时刻下，估计的剩余寿命能够很好地对应于各自的真值。同时，估计的概率密度函数的形状会随着监测时刻的增加而逐渐变得尖锐，说明预测的方差逐渐减小，预测精度趋于准确；
5)分别利用式(12)和式(13)计算模型的AIC和MSE，可得模型的拟合度为AIC＝-2.9714×103，均方误差为MSE＝1.3666×103，输出fm,k(lm,k)作为最终的剩余寿命预测结果。
b、实例二：以涡扇发动机数据退化数据为例，给出剩余寿命的预测结果：
1)利用NASA Ames预测数据库的涡扇发动机退化数据进一步验证所提模型的实用性。选取C-MAPSS模块的“FD003”数据集模拟真实商业涡扇发动机的退化，主要包括高压腔(HPC)退化和风扇(Fan)退化。涡扇发动机的全部退化数据由21个不同类型的传感器进行监测，每个传感器采样了100台同批次的涡扇发动机的退化数据。以传感器15为例，原始测量信号如图5所示；
2)注意到传感器8，传感器13以及传感器15的测量数据对应于同样类型的退化过程，因此，我们提取一个指定测试单元下这三个传感器的测量数据作为退化信号进行剩余寿命预测。考虑到复杂噪声源的存在，难以建立准确的噪声模型，因此我们首先利用平滑预处理的方法对退化数据进行去噪，同时对初值进行零位调整，修正后的涡扇发动机退化轨迹如图6所示。利用式(1)的模型结构初始化参数{μα,σα,β,σ,H}；
3)利用式(2)-式(6)求出μα,σα,β,σ,H的极大似然估计，模型参数的估计结果如下：μa＝2.8372×10-24，σa＝1.1572×10-24，b＝10.0089，σB＝0.0046；
4)利用式(7)-式(11)求取剩余寿命lm,k的概率密度函数fm,k(lm,k)，以传感器13为例，得到的剩余寿命预测结果如图7所示，显然，在各个监测时刻下，预测的概率密度函数的期望都能够较好地跟随真实的剩余寿命，验证了所提算法的有效性；
5)分别利用式(12)和式(13)计算模型的AIC和MSE，可得模型的拟合度为AIC＝-4.1261×103，均方误差为MSE＝318.1819，输出fm,k(lm,k)作为最终的剩余寿命预测结果。
当然，上述说明并非是对本发明的限制，本发明也并不仅限于上述举例，本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换，也应属于本发明的保护范围。

Claims (1)

1.一种考虑长程相关性和部件不确定性的寿命预测方法，其特征在于：在处理退化数据时是按照如下步骤依次进行的：
步骤1：读入M组状态监测退化数据，记为X＝(X1,X2,…,XM)；
步骤2：对退化数据进行零初值化处理，针对如下模型初始化参数，包括漂移系数的均值μα和标准差σα、非线性系数β、扩散系数σ以及Hurst指数H；
x(t)＝αtβ+σBH(t) (1)；
其中，x(t)是t时刻下的退化状态，α是服从高斯分布的漂移系数，β是固定的非线性系数，σ是扩散系数，BH(t)表示标准分数布朗运动；
步骤3：根据零初值化处理后的退化数据列写对数似然函数，如公式(2)所示；
其中，Xm＝(xm(t1),xm(t2),…,xm(tN))T,m＝1,2,…,M，Ψ＝σ2Q，i,j＝1,2,…,N，且有
<mrow> <mo>|</mo> <mo>&amp;Sigma;</mo> <mo>|</mo> <mo>=</mo> <mo>|</mo> <mi>&amp;Psi;</mi> <mo>|</mo> <mrow> <mo>(</mo> <msubsup> <mi>&amp;sigma;</mi> <mi>&amp;alpha;</mi> <mn>2</mn> </msubsup> <msup> <mi>&amp;Phi;</mi> <mi>T</mi> </msup> <msup> <mi>&amp;Psi;</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mi>&amp;Phi;</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
<mrow> <msup> <mo>&amp;Sigma;</mo> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>=</mo> <msup> <mi>&amp;Psi;</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>-</mo> <mfrac> <msubsup> <mi>&amp;sigma;</mi> <mi>&amp;alpha;</mi> <mn>2</mn> </msubsup> <mrow> <msubsup> <mi>&amp;sigma;</mi> <mi>&amp;alpha;</mi> <mn>2</mn> </msubsup> <msup> <mi>&amp;Phi;</mi> <mi>T</mi> </msup> <msup> <mi>&amp;Psi;</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mi>&amp;Phi;</mi> <mo>+</mo> <mn>1</mn> </mrow> </mfrac> <msup> <mi>&amp;Psi;</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mi>&amp;Phi;&amp;Phi;</mi> <mi>T</mi> </msup> <msup> <mi>&amp;Psi;</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
分别对μα和σα求偏导并令导数等于零，得到μα和σα的解析表达式，如下所示：
<mrow> <msub> <mover> <mi>&amp;mu;</mi> <mo>^</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <msup> <mi>&amp;Phi;</mi> <mi>T</mi> </msup> <msup> <mi>&amp;Psi;</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msub> <mi>X</mi> <mi>m</mi> </msub> </mrow> <mrow> <msup> <mi>M&amp;Phi;</mi> <mi>T</mi> </msup> <msup> <mi>&amp;Psi;</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mi>&amp;Phi;</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
<mrow> <msub> <mover> <mi>&amp;sigma;</mi> <mo>^</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mo>=</mo> <msup> <mrow> <mo>{</mo> <mfrac> <mn>1</mn> <mrow> <mi>M</mi> <msup> <mrow> <mo>(</mo> <msup> <mi>&amp;Phi;</mi> <mi>T</mi> </msup> <msup> <mi>&amp;Psi;</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mi>&amp;Phi;</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mi>m</mi> </msub> <mo>-</mo> <msub> <mi>&amp;mu;</mi> <mi>&amp;alpha;</mi> </msub> <mi>&amp;Phi;</mi> <mo>)</mo> </mrow> <mi>T</mi> </msup> <msup> <mi>&amp;Psi;</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mi>&amp;Phi;&amp;Phi;</mi> <mi>T</mi> </msup> <msup> <mi>&amp;Psi;</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mi>m</mi> </msub> <mo>-</mo> <msub> <mi>&amp;mu;</mi> <mi>&amp;alpha;</mi> </msub> <mi>&amp;Phi;</mi> <mo>)</mo> </mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mrow> <msup> <mi>&amp;Phi;</mi> <mi>T</mi> </msup> <msup> <mi>&amp;Psi;</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mi>&amp;Phi;</mi> </mrow> </mfrac> <mo>}</mo> </mrow> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
将(5)式和(6)式代入(2)式，利用多维搜索方法求出β,σ,H的极大似然估计，再将β,σ,H代回μα和σα的表达式，即(5)式和(6)式，求得μα和σα的极大似然估计；
步骤4：基于弱收敛准则将退化过程近似为马尔科夫过程，计算剩余寿命lm,k的概率密度函数fm,k(lm,k)，如公式(7)所示：：
<mrow> <msub> <mi>f</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>l</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>g</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>l</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>&amp;infin;</mi> </msubsup> <msub> <mi>g</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>l</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <msub> <mi>dl</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
其中，
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>g</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>l</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>&amp;omega;</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <msub> <mi>u</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>v</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <msub> <mi>c</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> </mrow> <mrow> <mo>&amp;lsqb;</mo> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mi>k</mi> </msub> <mo>+</mo> <msub> <mi>l</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <msqrt> <mrow> <mn>2</mn> <msubsup> <mi>&amp;pi;u</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> <mn>3</mn> </msubsup> </mrow> </msqrt> </mrow> </mfrac> <mfrac> <mrow> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mi>k</mi> </msub> <mo>+</mo> <msub> <mi>l</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>+</mo> <mi>&amp;Delta;</mi> <mi>l</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mi>k</mi> </msub> <mo>+</mo> <msub> <mi>l</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>l</mi> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>&amp;times;</mo> <mi>exp</mi> <mo>{</mo> <mo>-</mo> <mfrac> <msup> <mrow> <mo>&amp;lsqb;</mo> <msub> <mi>&amp;omega;</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>&amp;mu;</mi> <mi>&amp;alpha;</mi> </msub> <msup> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mi>k</mi> </msub> <mo>+</mo> <msub> <mi>l</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mi>&amp;beta;</mi> </msup> <mo>+</mo> <msub> <mi>&amp;mu;</mi> <mi>&amp;alpha;</mi> </msub> <msubsup> <mi>t</mi> <mi>k</mi> <mi>&amp;beta;</mi> </msubsup> <mo>&amp;rsqb;</mo> </mrow> <mn>2</mn> </msup> <mrow> <mn>2</mn> <msub> <mi>u</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> </mrow> </mfrac> <mo>}</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
其中，
<mrow> <msub> <mi>u</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>=</mo> <msup> <mrow> <mo>&amp;lsqb;</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mi>k</mi> </msub> <mo>+</mo> <msub> <mi>l</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mi>&amp;beta;</mi> </msup> <mo>-</mo> <msubsup> <mi>t</mi> <mi>k</mi> <mi>&amp;beta;</mi> </msubsup> <mo>&amp;rsqb;</mo> </mrow> <mn>2</mn> </msup> <msubsup> <mi>&amp;sigma;</mi> <mi>&amp;alpha;</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msup> <mi>&amp;sigma;</mi> <mn>2</mn> </msup> <mo>&amp;lsqb;</mo> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mi>k</mi> </msub> <mo>+</mo> <msub> <mi>l</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
<mrow> <msub> <mi>v</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>=</mo> <mo>&amp;lsqb;</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mi>k</mi> </msub> <mo>+</mo> <msub> <mi>l</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mi>&amp;beta;</mi> </msup> <mo>-</mo> <msubsup> <mi>t</mi> <mi>k</mi> <mi>&amp;beta;</mi> </msubsup> <mo>&amp;rsqb;</mo> <msubsup> <mi>&amp;sigma;</mi> <mi>&amp;alpha;</mi> <mn>2</mn> </msubsup> <msub> <mi>&amp;omega;</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>&amp;mu;</mi> <mi>&amp;alpha;</mi> </msub> <msup> <mi>&amp;sigma;</mi> <mn>2</mn> </msup> <mo>&amp;lsqb;</mo> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mi>k</mi> </msub> <mo>+</mo> <msub> <mi>l</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
<mrow> <msub> <mi>c</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>=</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mi>k</mi> </msub> <mo>+</mo> <msub> <mi>l</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mi>&amp;beta;</mi> </msup> <mo>-</mo> <msubsup> <mi>t</mi> <mi>k</mi> <mi>&amp;beta;</mi> </msubsup> <mo>-</mo> <mfrac> <mrow> <mi>&amp;beta;</mi> <msup> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mi>k</mi> </msub> <mo>+</mo> <msub> <mi>l</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mrow> <mi>&amp;beta;</mi> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mi>&amp;Delta;</mi> <mi>l</mi> <mo>&amp;lsqb;</mo> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mi>k</mi> </msub> <mo>+</mo> <msub> <mi>l</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mi>k</mi> </msub> <mo>+</mo> <msub> <mi>l</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>+</mo> <mi>&amp;Delta;</mi> <mi>l</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mi>k</mi> </msub> <mo>+</mo> <msub> <mi>l</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
Δl表示外推的采样间隔，ωm,k＝ωm-xm(tk)；
步骤5：根据公式(12)、(13)计算模型的AIC(Akaike information criterion，赤池信息量准则)和均方误差MSE，即
<mrow> <mi>M</mi> <mi>S</mi> <mi>E</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>M</mi> <mi>N</mi> </mrow> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>&amp;infin;</mi> </msubsup> <msup> <mrow> <mo>(</mo> <msub> <mi>l</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>-</mo> <msub> <mover> <mi>l</mi> <mo>~</mo> </mover> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msub> <mi>f</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>l</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <msub> <mi>dl</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
其中，p为模型未知参数的个数，M为系统部件的总数，N为退化数据的长度，为第m个部件在tk时刻下真实的剩余寿命，输出fm,k(lm,k)作为最终的剩余寿命预测结果。
CN201710554292.7A 2017-07-10 2017-07-10 一种考虑长程相关性和部件不确定性的寿命预测方法 Pending CN107688687A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710554292.7A CN107688687A (zh) 2017-07-10 2017-07-10 一种考虑长程相关性和部件不确定性的寿命预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710554292.7A CN107688687A (zh) 2017-07-10 2017-07-10 一种考虑长程相关性和部件不确定性的寿命预测方法

Publications (1)

Publication Number Publication Date
CN107688687A true CN107688687A (zh) 2018-02-13

Family

ID=61152835

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710554292.7A Pending CN107688687A (zh) 2017-07-10 2017-07-10 一种考虑长程相关性和部件不确定性的寿命预测方法

Country Status (1)

Country Link
CN (1) CN107688687A (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108829983A (zh) * 2018-06-21 2018-11-16 四川大学 基于多隐藏状态分数布朗运动的设备剩余寿命预测方法
CN109187038A (zh) * 2018-08-06 2019-01-11 大连理工大学 一种仿真获得船舶柴油机故障数据的方法
CN110197288A (zh) * 2019-05-30 2019-09-03 重庆大学 故障影响下设备的剩余使用寿命预测方法
CN111523727A (zh) * 2020-04-22 2020-08-11 北京航空航天大学 基于不确定过程的考虑恢复效应的电池剩余寿命预测方法
CN112068003A (zh) * 2020-11-16 2020-12-11 中南大学 基于线性维纳过程的镉镍蓄电池寿命预测方法和装置
CN112327169A (zh) * 2020-11-09 2021-02-05 上海工程技术大学 一种锂电池剩余寿命预测方法
CN112949058A (zh) * 2021-02-27 2021-06-11 中国人民解放军火箭军工程大学 一种基于隐含非线性维纳过程的锂电池启发式剩余寿命预测方法
CN113158471A (zh) * 2021-04-26 2021-07-23 中国人民解放军火箭军工程大学 一种考虑测量不确定性的退化设备剩余寿命预测方法
CN113492724A (zh) * 2021-07-28 2021-10-12 合肥国轩高科动力能源有限公司 基于赫斯特指数的动力电池sop预测方法、存储介质及设备
CN114925894A (zh) * 2022-05-09 2022-08-19 重庆大学 一种基于维纳过程和人工智能的发动机剩余寿命预测方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070239407A1 (en) * 2006-01-12 2007-10-11 Goldfine Neil J Remaining life prediction for individual components from sparse data
CN101976311A (zh) * 2010-11-22 2011-02-16 北京航空航天大学 基于漂移布朗运动模型的加速退化试验贝叶斯评估方法
CN105868557A (zh) * 2016-03-29 2016-08-17 浙江大学 一种两阶段退化情况下的机电设备剩余寿命在线预测方法
CN106484949A (zh) * 2016-09-12 2017-03-08 西安理工大学 基于退化数据的动量轮可靠性分析与剩余寿命预测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070239407A1 (en) * 2006-01-12 2007-10-11 Goldfine Neil J Remaining life prediction for individual components from sparse data
CN101976311A (zh) * 2010-11-22 2011-02-16 北京航空航天大学 基于漂移布朗运动模型的加速退化试验贝叶斯评估方法
CN105868557A (zh) * 2016-03-29 2016-08-17 浙江大学 一种两阶段退化情况下的机电设备剩余寿命在线预测方法
CN106484949A (zh) * 2016-09-12 2017-03-08 西安理工大学 基于退化数据的动量轮可靠性分析与剩余寿命预测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
司小胜 等: "《数据驱动的设备剩余寿命预测理论及应用》", 30 April 2016, 国防工业出版社 *

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108829983A (zh) * 2018-06-21 2018-11-16 四川大学 基于多隐藏状态分数布朗运动的设备剩余寿命预测方法
CN109187038A (zh) * 2018-08-06 2019-01-11 大连理工大学 一种仿真获得船舶柴油机故障数据的方法
CN109187038B (zh) * 2018-08-06 2019-07-16 大连理工大学 一种仿真获得船舶柴油机故障数据的方法
CN110197288A (zh) * 2019-05-30 2019-09-03 重庆大学 故障影响下设备的剩余使用寿命预测方法
CN110197288B (zh) * 2019-05-30 2023-06-02 重庆大学 故障影响下设备的剩余使用寿命预测方法
CN111523727B (zh) * 2020-04-22 2022-05-10 北京航空航天大学 基于不确定过程的考虑恢复效应的电池剩余寿命预测方法
CN111523727A (zh) * 2020-04-22 2020-08-11 北京航空航天大学 基于不确定过程的考虑恢复效应的电池剩余寿命预测方法
CN112327169A (zh) * 2020-11-09 2021-02-05 上海工程技术大学 一种锂电池剩余寿命预测方法
CN112068003B (zh) * 2020-11-16 2021-02-19 中南大学 基于线性维纳过程的镉镍蓄电池寿命预测方法和装置
CN112068003A (zh) * 2020-11-16 2020-12-11 中南大学 基于线性维纳过程的镉镍蓄电池寿命预测方法和装置
CN112949058A (zh) * 2021-02-27 2021-06-11 中国人民解放军火箭军工程大学 一种基于隐含非线性维纳过程的锂电池启发式剩余寿命预测方法
CN113158471A (zh) * 2021-04-26 2021-07-23 中国人民解放军火箭军工程大学 一种考虑测量不确定性的退化设备剩余寿命预测方法
CN113492724A (zh) * 2021-07-28 2021-10-12 合肥国轩高科动力能源有限公司 基于赫斯特指数的动力电池sop预测方法、存储介质及设备
CN113492724B (zh) * 2021-07-28 2023-12-08 合肥国轩高科动力能源有限公司 基于赫斯特指数的动力电池sop预测方法、存储介质及设备
CN114925894A (zh) * 2022-05-09 2022-08-19 重庆大学 一种基于维纳过程和人工智能的发动机剩余寿命预测方法

Similar Documents

Publication Publication Date Title
CN107688687A (zh) 一种考虑长程相关性和部件不确定性的寿命预测方法
Chen et al. A risk-averse remaining useful life estimation for predictive maintenance
CN110909822B (zh) 一种基于改进的高斯过程回归模型的卫星异常检测方法
Wen et al. A new multi-sensor fusion with hybrid Convolutional Neural Network with Wiener model for remaining useful life estimation
CN110807257A (zh) 航空发动机剩余寿命预测方法
CN102270302B (zh) 一种基于灰色支持向量机的多应力加速寿命试验预测方法
Son et al. RUL prediction for individual units based on condition monitoring signals with a change point
Huang et al. A novel Bayesian deep dual network with unsupervised domain adaptation for transfer fault prognosis across different machines
CN109033499A (zh) 一种多阶段一致性检验的航空发动机剩余寿命预测方法
CN103197663B (zh) 一种故障预测方法及系统
CN104778354B (zh) 定期抽检型产品贮存寿命评估方法
Zhang et al. Remaining Useful Life Prediction of Rolling Bearings Using Electrostatic Monitoring Based on Two‐Stage Information Fusion Stochastic Filtering
CN105867360A (zh) 一种机电控制系统的初值预估迭代学习故障诊断算法
CN108446714A (zh) 一种多工况下的非马尔科夫退化系统剩余寿命预测方法
Wang et al. An improved fusion prognostics method for remaining useful life prediction of bearings
Wu et al. Real-time monitoring and diagnosis scheme for IoT-enabled devices using multivariate SPC techniques
Lv et al. A spatial–temporal variational graph attention autoencoder using interactive information for fault detection in complex industrial processes
CN106971058A (zh) 一种泵站运行监控数据异常点检测方法
Yan et al. Two‐Stage Degradation Assessment and Prediction Method for Aircraft Engine Based on Data Fusion
CN104156612A (zh) 基于粒子滤波的正向与逆向预测误差的故障预报方法
CN110532629A (zh) 一种多模态退化过程的剩余寿命预测方法
CN105718733B (zh) 基于模糊贴近度和粒子滤波的故障预报方法
Mishra et al. Hybrid models for rotating machinery diagnosis and prognosis: estimation of remaining useful life
CN108932197A (zh) 基于参数Bootstrap重抽样的软件失效时间预测方法
Luo et al. A novel method for remaining useful life prediction of roller bearings involving the discrepancy and similarity of degradation trajectories

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20180213

RJ01 Rejection of invention patent application after publication