CN113191506B - 一种考虑设备检测不确定性情况下的非周期视情维修方法 - Google Patents
一种考虑设备检测不确定性情况下的非周期视情维修方法 Download PDFInfo
- Publication number
- CN113191506B CN113191506B CN202110395374.8A CN202110395374A CN113191506B CN 113191506 B CN113191506 B CN 113191506B CN 202110395374 A CN202110395374 A CN 202110395374A CN 113191506 B CN113191506 B CN 113191506B
- Authority
- CN
- China
- Prior art keywords
- detection
- maintenance
- equipment
- uncertainty
- degradation
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 193
- 238000012423 maintenance Methods 0.000 title claims abstract description 96
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000006731 degradation reaction Methods 0.000 claims abstract description 69
- 230000015556 catabolic process Effects 0.000 claims abstract description 38
- 238000009826 distribution Methods 0.000 claims abstract description 23
- 238000005457 optimization Methods 0.000 claims abstract description 16
- 239000002245 particle Substances 0.000 claims abstract description 15
- 238000012360 testing method Methods 0.000 claims abstract description 10
- 230000003449 preventive effect Effects 0.000 claims description 19
- 230000001965 increasing effect Effects 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000008439 repair process Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 12
- 238000012544 monitoring process Methods 0.000 abstract description 4
- 230000000737 periodic effect Effects 0.000 abstract description 4
- 238000011069 regeneration method Methods 0.000 abstract description 3
- 238000007726 management method Methods 0.000 abstract 1
- 239000000463 material Substances 0.000 abstract 1
- 230000009471 action Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000007689 inspection Methods 0.000 description 5
- 230000005653 Brownian motion process Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 238000013398 bayesian method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/20—Administration of product repair or maintenance
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N7/00—Computing arrangements based on specific mathematical models
- G06N7/01—Probabilistic graphical models, e.g. probabilistic networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Artificial Intelligence (AREA)
- Data Mining & Analysis (AREA)
- Quality & Reliability (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Marketing (AREA)
- Software Systems (AREA)
- Mathematical Physics (AREA)
- Tourism & Hospitality (AREA)
- General Engineering & Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Computing Systems (AREA)
- Evolutionary Computation (AREA)
- Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computational Linguistics (AREA)
- Development Economics (AREA)
- Probability & Statistics with Applications (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Testing And Monitoring For Control Systems (AREA)
Abstract
本发明公开了一种考虑设备检测不确定性情况下的非周期视情维修方法,具体为:首先根据已有的设备退化数据集得到相应的退化分布,并求出分布类型;然后假设退化分布的各参数所服从的分布,利用假设检验方法,获得此参数的先验信息;假定维修过程阈值以及采用的检测规划方程参数为决策变量,考虑每一次检测服从高斯分布的随机冲击,再利用半再生过程建立非周期维修优化模型;最后采用基于灾变策略的粒子群算法,得到决策变量的相关最优值解,并通过贝叶斯理论更新模型的后验参数。本发明降低了传统维修策略的连续监测与周期检测所需人力,物力等成本的影响,并考虑检测带来的不确定性的影响,为视情维修检测以及相应的预测与健康管理等提供了技术参考。
Description
技术领域
本发明涉及一种非周期视情维修方法。
背景技术
随着现代工业的发展,设备(组件和系统)变得越来越可靠。通常,很难准确估计高度可靠的设备的可靠性和寿命。随机过程通常用于预测组件和系统的寿命分布。尤其是伽马过程,由于其独立,非负,严格增加的特性而备受关注。实施维护是为了控制和延长设备的剩余使用寿命(RUL)。通常,维护可以分为两大类:替换性维护(CM)和预防性维护(PM)。替换性维护可能会导致较高的生产损失和维护成本,因此PM策略通常受到更多的关注。在监视设备状况时,可以考虑基于状态的维护(CBM),它会根据设备观察到的系统状态动态地做出维护决定。基于时间的维护(TBM)是另一种PM策略,该策略根据故障时间分析的结果定期设置检测和维护:设备的预期寿命是根据故障时间或基于使用情况的数据估计的。
在定义维护决策策略时,还需要考虑由于测量设备故障引起的不完善的检查。特别是采用视情维护的方法时,由于其维护决策都是基于相应检测得到,因此检测的不确定性可以明显的影响视情维护的决策过程。成本效率是维护决策的一项重要标准。通常,寻求最优的维护策略时,其要求是以最小的成本满足相应的维护要求。而收集的信息可以进一步支持此种维护策略。而最直观的信息是系统的状态,它有助于确定是否需要维修或更换系统。系统信息的收集需要进行相应的安排,而连续状态监视可以允许收集大量信息以指导维护决策,但是这样做可能会花费巨大成本。周期性的检测可以有效的避免连续监测的高成本问题,但在检测阶段仍然花费可观。
发明内容
发明目的:为了解决上述现有技术存在的问题,本发明提供了一种考虑设备检测不确定性情况下的非周期视情维修方法。
技术方案:一种考虑设备检测不确定性情况下的非周期视情维修方法,具体包括如下步骤:
步骤1:进入第p个检测周期,针对设备单调性变化的退化数据,采用伽马函数对退化过程进行退化分析,得到伽马退化过程;
步骤2:采用假设检验的方法得到伽马退化过程中参数的先验信息;
步骤3:根据设备的状态Z和设备的早期退化速率,建立检测规划函数M(Z),该检测规划函数用于计算下一次设备检测的时刻;在检测不确定情况下建立设备的状态方程;
步骤4:进行第k次检测时,根据预设的设备失效阈值Lf,得到检测不确定情况下的设备退化过程首达时分布;
步骤5:建立考虑检测不确定情况下的非周期视情维修优化模型C;
步骤6:采用基于灾变策略的粒子群算法对非周期视情维修优化模型C进行优化计算,得到非周期视情维修优化模型C中决策变量的值,所述决策变量包括预防性维护阈值Lp;
步骤7:若第p个检测周期中第k次检测时设备的状态值<Lp,则认为设备不需要维修,计算第p个检测周期中进入第k+1次检测的时刻,同时更新伽马退化过程中的参数的信息,令k=k+1,并返回步骤4;
如果Lf>第p个检测周期中第k次检测时设备的状态值≥Lp,则对设备进行预防性维护,令p=p+1,k=1,计算第p个检测周期中进入第1次检测的时刻,然后返回步骤1;
如果第p个检测周期中第k次检测时设备状态值≥Lf,则认为设备故障,且执行替换性维护,令p=p+1,k=1,计算第p个检测周期中进入第1次检测的时刻,然后返回步骤1。
进一步的,所述步骤1中采用伽马函数对退化过程进行退化分析得到伽马退化过程f(x|αΛ(t),β):
其中,为完全伽马函数,y为积分变量,Λ(t)=tb为非负增函数,t为时间,且Λ(0)=0,x为设备的退化路径,exp(.)表示指数函数,α、β和b为伽马退化过程的参数。
进一步的,步骤3中所述检测规划函数M(Z)为:
当设备早期退化速率大于预设的阈值时,M(Z)为:
其中p1,p2为参数;
当设备早期退化速率小于等于预设的阈值时,M(Z)为:
所述设备的状态方程为:
为第p个检测周期中进行第k次检测的时刻,/>为时刻/>时设备的状态,为时刻/>时设备的真实状态,/>为第p个检测周期中第i次检测时的干扰不确定项。
进一步的,所述步骤4中设备退化过程首达时分布为:
其中为第p个检测周期中从第一次检测到第k次检测所有干扰不确定项的均值,为第p个检测周期中从第一次检测到第k次检测所有干扰不确定项的残差;Φ[·]为标准正态分布,/>
进一步的,所述步骤5中非周期视情维修优化模型C为:
其中 表示第一个检测周期中第一次检测的时间,CI表示进行设备检测的权值,Cp表示预防维修的权值,Cc表示完美维修的权值,Cd表示失效到维修之间的停机状态的权值;所述完美维修为使系统的状态恢复至初始状态的维修;
其中s为积分变量;
π(dZ)由马尔可夫链稳态分布得到,表示为:
π(dZ)=vδ0(dZ)+(1-v)r(Z)dZ
其中δ0(.)为Dirac函数,fM(0)(Z)为在检测规划函数M(0)下状态为Z的退化过程的概率密度函数,fM(Z)(Z-x)为在检测规划函数M(Z)下状态为Z-x的退化过程的概率密度函数。
进一步的,步骤7中设备不需要维修时,计算第p个检测周期中进入第k+1次检测的时刻
进行预防性维护或者替换性维护时,计算第p个检测周期中进入第1次检测的时刻为:
其中T为第p-1个检测周期中最后一次检测的时刻,Zp-1为第p-1个检测周期中最后一次检测时设备的状态;
通过如下贝叶斯公式,对参数α,β和b进行更新:
表示当前检测周期内第k次检测与第k-1次检测之间的退化增量,/>为第p个检测周期中进行第k-1次检测的时刻,L(.)为似然函数,g(α,β,b)表示参数α,β,b的先验信息。
有益效果:本方法能够将视情维修下的检测不确定性影响纳入到维修优化模型中,此模型考虑到了更贴合实际应用情况,使优化模型更精确。并且,采用非周期的检测策略,可以进一步节省维修成本。通过贝叶斯方法进行更新,可以进一步实现模型的在线优化。
附图说明
图1是本方法整体流程图;
图2是激光数据随时间的退化波形图;
图3是相应维护方案的动作;
图4是检测规划函数示意图;
图5是在不同检测规划函数下不考虑检测不确定的C值随迭代的演化图;
图6是不同预防性维护阈值下及不同检测规划函数下不考虑检测不确定的C值演化图,其中图(a)为m1p(Z)函数的演示规划图;图(b)为m2p(x)函数的演示规划图;
图7是在不同检测规划函数下考虑检测不确定的C值随迭代的演化图;
图8是不同预防性维护阈值下及不同检测规划函数下考虑检测不确定的C值演化图,其中图(a)为m1p(Z)函数的演示规划图,图(b)为m2p(Z)函数的演示规划图。
具体实施方式
构成本发明的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。
如图1所示,本实施提供了一种考虑设备检测不确定的非周期视情维修方法,包括以下步骤:
(1)根据设备退化的数据,选取合适的退化分布。通常,非严格单调性的数据可以选取维纳过程拟合退化数据,单调性时可以选取伽马过程等拟合退化数据。
(2)将步骤(1)中所述的退化分布选取为伽马退化过程,并考虑其分布参数具有随机效应,并假定分布参数的先验信息。
(3)根据视情维修的相关要求,假定维修阈值Lp与检测规划函数中的参数p1,p2为维修优化模型中的决策变量。并且维修动作是根据检测的当前系统状态而执行,即当检测的当前状态超出预防性维修阈值,则进行完美维护;所述完美维护为使系统的状态恢复至初始状态的维护。
(4)考虑检测的不确定性,并假定检测的不确定性可以表示为服从高斯分布的随机冲击。因此,对于退化过程来说,检测的不确定性可以同等的视为将此不确定性传递给失效阈值与维修阈值,使固定阈值变为动态随机阈值,从而可以将检测不确定性考虑到维修优化模型中。使未考虑检测不确定情况下的系统固定阈值变为动态随机阈值,从而可以将检测不确定性纳入到维修优化模型中。
(5)利用半再生过程以及步骤(4)所考虑的检测不确定条件下的退化模型,构建基于非周期检测的视情维修优化模型。
(6)通过步骤(5)建立的考虑检测不确定条件下的非周期维修优化模型,采用基于灾变策略的粒子群优化算法,可以得到维修优化模型中的决策变量最优值解。若当前检测时设备的状态值小于Lp,则认为设备不需要维修,利用现场数据,通过贝叶斯理论对伽马函数的参数进行更新,得到更新后的退化过程参数,并重复步骤(3)-(6),实现模型的“在线”优化。如果Lf>当前检测时设备状态值≥Lp,则对设备进行预防性维护,并结束当前检测周期进行下一个检测周期,然后返回步骤1;如果第k次检测时设备状态值≥Lf,则认为设备故障,且执行替换性维护,并结束当前检测周期进行下一个检测周期,然后返回步骤1;
本方法的整体流程如图2所示,具体实施步骤如下:
步骤1,激光器退化数据的可用分布选择:退化过程可以选择维纳过程,伽马过程等。通过对数据进行分析,伽马过程更适合于单调递增的退化过程,因此采用伽马过程描述退化数据如下:
其中,为完全伽马函数,y为积分变量,Λ(t)=tb为非负增函数,t表示时间,且Λ(0)=0,x为设备的退化路径,exp(.)表示指数函数,退化路径x的期望和方差则分别为:
E[x]=β·αΛ(t)
Var[x]=β2·αΛ(t)
步骤2,利用假设检验的方法,可以得到退化过程的参数的最优拟合分布也既先验信息,其中α最优拟合于对数正态分布,即α~LN(-2.6343,0.4505);β最优拟合于伽马分布,即,β~Ga(5.4815,0.0108);b最优拟合于正态分布,即,b~N(0.9504,0.0045)。
如图3所示,需要制定有关非定期检查的维护策略,以便在每个检查时刻根据系统的状态采取必要的措施。令表示非周期性检查的顺序,/>表示第p个检测周期中第k次检测的时刻,并假设维护时间可以忽略不计,且只能通过检查来判断系统是否故障。当进行第k次检测时的系统状态/>超过预防性阈值Lp但小于失效阈值Lf且Lp<Lf时,将采取预防性维护措施。从图3可以推断出,如果检查是可靠的,则在第p个检测周期中进行第k次检测的时刻/>处以概率1实施PM动作。如果系统在检查时的状态超过故障阈值Lf,则执行替换性维护动作。同样,在图3中,如果在第C个检测周期中进行第k次检测的时刻/>处的检查是可靠的,则认为替换性维护动作的概率为1。如果在检查时的系统条件低于Lp,则认为系统仍在运行,并且不需要维护措施。
步骤3、检测规划函数M(Z)如图4所示,其中两个函数分别为凹函数与凸函数,适用于不同的条件。m1p(Z)对应的检测时间间隔随着系统退化下降速率较快,也即系统退化下降速率大于阈值;其初始检测间隔比较小,因此适用于早期退化速率较大的系统。m2p(Z)对应的检测时间间隔随着系统退化下降速率较缓慢,也即系统退化下降速率小于等于阈值,其初始检测间隔比较大,因此适用于早期退化比较缓慢的系统。检测规划函数表示如下:
p1,p2为参数;
步骤4、考虑检测不确定条件下,假设其不确定性由服从高斯分布的随机冲击决定,因此,系统的状态可以由下式表示:
为第p个检测周期中进行第k次检测的时刻,/>为时刻/>时设备的状态,为时刻/>时设备的真实状态,/>为第p个检测周期中第i次检测时的干扰不确定项。
此外,令
表示由全部检测k次导致的系统退化水平。通常,由测量设备引起的不可靠检查被视为正态分布的外部冲击。
根据专家经验可知设备的失效阈值为Lf,其检测不确定下的系统退化过程首达时分布可以表示为:
其中为第p个检测周期中从第一次检测到第k次检测所有干扰不确定项的均值,为第p个检测周期中从第一次检测到第k次检测所有干扰不确定项的残差;Φ[·]为标准正态分布,/>
步骤5、采用半再生过程建立考虑检测不确定情况下的非周期视情维修优化模型C,其可以建模为:
其中 表示第一个检测周期中第一次检测的时间,CI表示进行设备检测的权值,Cp表示预防维修的权值,Cc表示完美维修的权值,Cd表示失效到维修之间的停机状态的权值;所述完美维修为使系统的状态恢复至初始状态的维修;
其中s为积分变量;
π(dZ)由马尔可夫链稳态分布得到,表示为:
π(dZ)=vδ0(dZ)+(1-v)r(Z)dZ
其中δ0(.)为Dirac函数,fM(0)(Z)为在检测规划函数M(0)下状态为Z的退化过程的概率密度函数,fM(Z)(Z-x)为在检测规划函数M(Z)下状态为Z-x的退化过程的概率密度函数;x为设备的退化路径。
步骤6、通过步骤5得到的优化模型,采用基于灾变策略的粒子群算法,获得步骤3中所述的决策变量的最优值解。粒子群算法相应的操作为:
(1)对基于灾变策略的粒子群算法中的加速因子,种群数量,迭代次数,决策变量的初值进行初始化,并随机选择全局最优点和个体最优点。
(2)计算适应度函数,即步骤5中的模型C,使其有最小值。
(3)经过上述计算,如果当前最优值比历史最优值更优,则历史最优值将被当前最优值替代。
(4)根据速度更新公式调整每个粒子的速度,并根据位置公式调整每个粒子的位置,从而更新速度和位置。
(5)判断灾变发生条件,即当前粒子的最佳位置。如果否则继续进行下一步。
(6)检查停止迭代的条件,当迭代次数达到最大值时,停止迭代,得到步骤3中决策变量的最优值解。否则,回到第三步。
用于灾变策略的判断条件为:
ifx(i’,:)==gbest
L_p=6*rand(1,1);
p1=randi’([1,100],1,1);
p2=randi’([2,100],1,1);
giv=[L_p p1 p2];
x(i’,:)=giv;
end
其中x(i’,:)为总群中第i’个粒子位置,gbest为当前粒子群最佳位置,L_p,p1,p2分别为预防性维修阈值,以及检测规划函数中的参数。当达到判断条件时实施灾变策略,对粒子群进行重新赋值运算。
设备不需要维修时,计算第p个检测周期中中进入第k+1次检测的时刻
进行预防性维护或者替换性维护时,令p=p+1,k=1,计算第p个检测周期中进入第1次检测的时刻为:
其中T为第p-1个检测周期中最后一次检测的时刻,Zp-1为第p-1个检测周期中最后一次检测时设备的状态;
通过如下贝叶斯公式,对参数α,β和b进行更新:
表示当前检测周期内第k次检测与第k-1次检测之间的退化增量,/>为第p个检测周期中进行第k-1次检测的时刻,L(.)为似然函数,g(α,β,b)表示参数α,β,b的先验信息。
图5,图6,图7,图8为不同条件下的C值演化图。其中图5和图7所示分别为不同检测规划函数下不考虑/考虑检测不确定时的C值随迭代的演化,从图中可以看出,不考虑检测不确定时的C值要少于考虑检测不确定时的C值。因此,实际中,检测不确定的存在会增加系统损耗。
图6和图8分别表示不同预防性维护阈值下及不同检测规划函数下不考虑/考虑检测不确定的C值演化图。从图中可以看出,不同决策变量p1,p2值会导致不同的C值,且通过粒子群算法可以得到最优的决策变量的解。
上面结合附图对本发明的实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下做出各种变化。
Claims (6)
1.一种考虑设备检测不确定性情况下的非周期视情维修方法,其特征在于,具体包括如下步骤:
步骤1:进入第p个检测周期,针对设备单调性变化的退化数据,采用伽马函数对退化过程进行退化分析,得到伽马退化过程;
步骤2:采用假设检验的方法得到伽马退化过程中参数的先验信息;
步骤3:根据设备的状态Z和设备的早期退化速率,建立检测规划函数M(Z),该检测规划函数用于计算下一次设备检测的时刻;在检测不确定情况下建立设备的状态方程;
步骤4:进行第k次检测时,根据预设的设备失效阈值Lf,得到检测不确定情况下的设备退化过程首达时分布;
步骤5:建立考虑检测不确定情况下的非周期视情维修优化模型C;
步骤6:采用基于灾变策略的粒子群算法对非周期视情维修优化模型C进行优化计算,得到非周期视情维修优化模型C中决策变量的值,所述决策变量包括预防性维护阈值Lp;
步骤7:若第p个检测周期中第k次检测时设备的状态值<Lp,则认为设备不需要维修,计算第p个检测周期中进入第k+1次检测的时刻,同时更新伽马退化过程中的参数的信息,令k=k+1,并返回步骤4;
如果Lf>第p个检测周期中第k次检测时设备的状态值≥Lp,则对设备进行预防性维护,令p=p+1,k=1,计算第p个检测周期中进入第1次检测的时刻,然后返回步骤1;
如果第p个检测周期中第k次检测时设备状态值≥Lf,则认为设备故障,且执行替换性维护,令p=p+1,k=1,计算第p个检测周期中进入第1次检测的时刻,然后返回步骤1。
2.根据权利要求1所述的一种考虑设备检测不确定性情况下的非周期视情维修方法,其特征在于,所述步骤1中采用伽马函数对退化过程进行退化分析得到伽马退化过程f(x|αΛ(t),β):
其中,为完全伽马函数,y为积分变量,Λ(t)=tb为非负增函数,t为时间,且Λ(0)=0,x为设备的退化路径,exp(.)表示指数函数,α、β和b为伽马退化过程的参数。
3.根据权利要求2所述的一种考虑设备检测不确定性情况下的非周期视情维修方法,其特征在于,步骤3中所述检测规划函数M(Z)为:
当设备早期退化速率大于预设的阈值时,M(Z)为:
其中p1,p2为参数;
当设备早期退化速率小于等于预设的阈值时,M(Z)为:
所述设备的状态方程为:
为第p个检测周期中进行第k次检测的时刻,/>为时刻/>时设备的状态,/>为时刻/>时设备的真实状态,/>为第p个检测周期中第i次检测时的干扰不确定项。
4.根据权利要求3所述的一种考虑设备检测不确定性情况下的非周期视情维修方法,其特征在于,所述步骤4中设备退化过程首达时分布为:
其中为第p个检测周期中从第一次检测到第k次检测所有干扰不确定项的均值,/>为第p个检测周期中从第一次检测到第k次检测所有干扰不确定项的残差;Φ[·]为标准正态分布,/>
5.根据权利要求4所述的一种考虑设备检测不确定性情况下的非周期视情维修方法,其特征在于,所述步骤5中非周期视情维修优化模型C为:
其中表示第一个检测周期中第一次检测的时间,CI表示进行设备检测的权值,Cp表示预防维修的权值,Cc表示完美维修的权值,Cd表示失效到维修之间的停机状态的权值;所述完美维修为使系统的状态恢复至初始状态的维修;
其中s为积分变量;
π(dZ)由马尔可夫链稳态分布得到,表示为:
π(dZ)=vδ0(dZ)+(1-v)r(Z)dZ
其中δ0(.)为Dirac函数,fM(0)(Z)为在检测规划函数M(0)下状态为Z的退化过程的概率密度函数,fM(z)(Z-x)为在检测规划函数M(Z)下状态为Z-x的退化过程的概率密度函数。
6.根据权利要求4所述的一种考虑设备检测不确定性情况下的非周期视情维修方法,其特征在于,步骤7中设备不需要维修时,计算第p个检测周期中进入第k+1次检测的时刻
进行预防性维护或者替换性维护时,计算第p个检测周期中进入第1次检测的时刻为:
其中T为第p-1个检测周期中最后一次检测的时刻,Zp-1为第p-1个检测周期中最后一次检测时设备的状态;
通过如下贝叶斯公式,对参数α,β和b进行更新:
表示当前检测周期内第k次检测与第k-1次检测之间的退化增量,/>为第p个检测周期中进行第k-1次检测的时刻,L(.)为似然函数,g(α,β,b)表示参数α,β,b的先验信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110395374.8A CN113191506B (zh) | 2021-04-13 | 2021-04-13 | 一种考虑设备检测不确定性情况下的非周期视情维修方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110395374.8A CN113191506B (zh) | 2021-04-13 | 2021-04-13 | 一种考虑设备检测不确定性情况下的非周期视情维修方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113191506A CN113191506A (zh) | 2021-07-30 |
CN113191506B true CN113191506B (zh) | 2023-11-07 |
Family
ID=76975606
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110395374.8A Active CN113191506B (zh) | 2021-04-13 | 2021-04-13 | 一种考虑设备检测不确定性情况下的非周期视情维修方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113191506B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115330001B (zh) * | 2022-10-11 | 2022-12-23 | 合肥工业大学 | 考虑多任务模式的动态视情维修策略优化方法和系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103336903A (zh) * | 2013-07-01 | 2013-10-02 | 中国石油大学(华东) | 一种基于贝叶斯理论的石化设备失效率推断方法 |
CN110320893A (zh) * | 2019-07-23 | 2019-10-11 | 浙江浙能嘉华发电有限公司 | 一种基于rcm的火电机组设备检修方法 |
CN111428360A (zh) * | 2020-03-23 | 2020-07-17 | 西安建筑科技大学 | 一种不同维修深度分布下海底管道剩余寿命预测方法 |
-
2021
- 2021-04-13 CN CN202110395374.8A patent/CN113191506B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103336903A (zh) * | 2013-07-01 | 2013-10-02 | 中国石油大学(华东) | 一种基于贝叶斯理论的石化设备失效率推断方法 |
CN110320893A (zh) * | 2019-07-23 | 2019-10-11 | 浙江浙能嘉华发电有限公司 | 一种基于rcm的火电机组设备检修方法 |
CN111428360A (zh) * | 2020-03-23 | 2020-07-17 | 西安建筑科技大学 | 一种不同维修深度分布下海底管道剩余寿命预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113191506A (zh) | 2021-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10521490B2 (en) | Equipment maintenance management system and equipment maintenance management method | |
CN104966141B (zh) | 更新用于生成工业资产健康状况简档的模型的方法和系统 | |
Coble et al. | Identifying optimal prognostic parameters from data: a genetic algorithms approach | |
Huynh et al. | A parametric predictive maintenance decision-making framework considering improved system health prognosis precision | |
EP2837984B1 (en) | Process to optimize the maintenance of technical systems | |
Compare et al. | Predictive maintenance by risk sensitive particle filtering | |
Li et al. | Preventive maintenance decision model of urban transportation system equipment based on multi-control units | |
Duan et al. | Optimal failure mode-based preventive maintenance scheduling for a complex mechanical device | |
Chen et al. | Imperfect maintenance policy considering positive and negative effects for deteriorating systems with variation of operating conditions | |
CN113191506B (zh) | 一种考虑设备检测不确定性情况下的非周期视情维修方法 | |
Cao et al. | Optimization of condition-based maintenance for multi-state deterioration systems under random shock | |
Montgomery et al. | Minor maintenance actions and their impact on diagnostic and prognostic CBM models | |
Baraldi et al. | Modelling the effects of maintenance on the degradation of a water-feeding turbo-pump of a nuclear power plant | |
CN105678089A (zh) | 模型自匹配融合健康预测方法 | |
CN109033569B (zh) | 一种用于舰载机传感器系统预防检修阈强度和次数优化的方法 | |
Escobet et al. | A fault/anomaly system prognosis using a data-driven approach considering uncertainty | |
Kirschenmann et al. | Decision dependent stochastic processes | |
CN116128145A (zh) | 电力设备状态检修策略优化方法 | |
CN112966785B (zh) | 一种智能化星座状态识别方法和系统 | |
CN109933890B (zh) | 一种产品综合维修方法和装置 | |
Wang et al. | Maintenance decision based on data fusion of aero engines | |
Eker et al. | State-based prognostics with state duration information of cracks in structures | |
Blair et al. | Hybrid fault prognostics for nuclear applications: addressing rotating plant model uncertainty | |
Belhaj Salem et al. | Prognostic and Classification of the Dynamic Degradation in a Mechanical System Using Variance Gamma Process. Mathematics 2021, 9, 254 | |
Khan et al. | Prognostics of crack propagation in structures using time delay neural network |
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 |