基于状态监测的寿命预测与预防维护规划研究综述

基于状态监测的寿命预测与预防维护规划研究综述

该设备是否能够持续有效地开展工作，以及是否能够及时有效地实施，与当局、公司乃至国家的竞争力密切相关。如今,越来越复杂精巧的设备使其有效运行的可靠性成为突出问题;另一方面,竞争激烈的市场又要求设备能更可靠高效地工作。由此,高效、经济的维护策略的重要性就显而易见了。1981年,美国国内工厂针对重要设备的维护费用为6000亿美元。20年内,这一数字翻了20倍。一个更惊人的事实是,由于无效、低效的维护,有1/3至1/2的维护费用被浪费了。因此,就需要不断发展与改进现有的维护策略。如今,设备维护策略已由事后维护(correctivemaintenance),基于时间的预防维护(time-basedpreventivemaintenance),逐步发展为基于状态的维护(condition-basedmaintenance(CBM))。对于维护费用的分析显示,设备失效后维护的费用往往数倍于采用CBM策略的维护费用。另一份调查显示,对状态监测投入1万美元-2万美元,每年有望节省50万美元的维护费用。CBM的益处不仅是维护费用的降低,还包括:设备失效时间的减少(设备可用度的提高),库存的削减,和后勤与供应链的改善。例如:美国“智能维护系统”中心(centerforintelligentmaintenancesystems)估计,仅在美国国内,若采用CBM策略,仅因设备可用度提高这一项,每年就可节省50亿美元。CBM对维护决策的支持主要有两方面:诊断(diagnostics)与预测(prognostics)。Diagnostics致力于失效发生后的发现、隔离与鉴定;Prognostics致力于预测、防止失效的发生。虽然无法完全代替diagnostics,prognostics在预防失效发生方面更为有效。然而,相对diagnostics而言,prognostics方面的研究较少。Prognostics主要有两大方面的研究:寿命预测与预防维护规划。文献总结了至2008年寿命预测方面的一些工作,文献总结了至2005年寿命预测和预防维护规划两大方面的主要研究。然而,一方面prognostics的研究日新月异,新的方法、模型不断涌现,需要及时总结;另一方面需要对prognostics的各种方法、模型进行分类和归纳,将有助于掌握各类方法的特点及其适用性,并发现新的研究方向。本文对目前为止主要的寿命预测方法和预防维护规划模型进行总结、分类和简单的比较,进而探讨一些必要的、有价值的研究方向。1状态监测数据基于状态监测的设备寿命预测,具体而言,是在设备A运行的某一时刻t,根据监测的(至时刻t的)设备A的运行状态和/或同类设备的历史数据,预测设备A由当前至失效的剩余寿命(residuallife)(或剩余有效寿命(remainingusefullife))。这里,历史数据可以是同类设备从运行到失效过程中的状态监测数据(conditionmonitoringdata),可以是失效时间数据、维护时间数据等事件数据(eventdata),也可以是两者的综合。Jardine等将剩余寿命表述为:Τ-t|Τ＞t,Ζ(t)(1)式中,T为表示失效时刻的随机变量,t为当前时刻,向量Z(t)为至当前的状态数据。更全面地,可将式(1)拓展为:Τ-t|Τ＞t,Ζ(u),(Ζhj(v),Ehj)(2)式中,向量Zhj(v)为同类设备j在运行时刻ν的状态监测数据,向量Ehj为同类设备j的事件数据(尤其是失效时间,预防维护时间等数据),u∈[0,t],ν∈[0,Thj],j=1,2,…,L,L表示可获得历史数据的同类设备的个数,Thj为同类设备j的最终失效或不再服役的时刻。寿命预测,既指预测变量T-t的分布,有时也指预测T-t的期望,即:E[Τ-t|Τ＞t,Ζ(u),(Ζhj(v),Ehj)](3)文献中,设备寿命预测的方法主要可以分为:基于物理(第一原理)的方法和基于经验的方法,下文1.1节和1.2节总结近年来关于这两类方法的主要工作。1.1基于状态约束的转轴裂纹关联研究基于物理的方法试图结合:①设备特定的机械动力学知识(如:损伤增长模型);②设备在线状态监测数据,以预测设备的剩余寿命。文献中常见的基于物理的方法主要使用裂纹扩展模型,碎裂增长模型等。Ray与Tangirala提出了裂纹扩展的非线性随机模型,结合实测传感器信号估计裂纹长度、裂纹扩展速率及材料的剩余寿命。Li等提出了一个滚动轴承寿命预测的框架,该框架包括一个确定性损伤扩展模型,一个诊断模型和一个自适应算法。其中,诊断模型用于从传感性信号中估计缺陷尺寸;由此估计,损伤扩展模型可计算得缺陷扩展率,并预测下一时刻的损伤尺寸;自适应算法根据预测值和实际值的偏差在线更新损伤扩展模型的参数,这就抓住了损伤扩展的动态特性。可是,诊断模型并未予以详细讨论,在文献的实验中,直接使用了实测的损伤尺寸。在后续的工作中,Li等将确定性损伤扩展模型推广到随机损伤扩展模型。Glodez等提出了一个齿轮寿命的确定性计算模型,该模型以应变-寿命方法计算裂纹萌生所需的应力循环数,以Paris方程模拟裂纹扩展过程,以位移相关法计算应力强度因子与裂纹长度的函数关系,最后估计齿轮的寿命。基于振动响应分析和损伤力学,Qiu等提出了一个基于刚度的轴承系统预测模型。在文献中,轴承系统作为单自由度振动系统,因此其固有频率和加速度幅值可与系统刚度相关联。另一方面,通过损伤力学,可将失效时间,服役时间和刚度变化联系起来。这样,一个轴承系统的固有频率和加速度幅值便可与失效时间和服役时间联系起来。因此,可以通过振动测量在线预测轴承系统的失效时间。Oppenheimer和Loparo提出了一个基于物理的转轴裂纹诊断与预测方法,该方法包括监测器和寿命模型。假设可由实测信号计算轴端速度和轴端力,根据转轴的动力学模型,监测器可估计裂纹的长度。而后,寿命模型基于线弹性断裂力学Forman裂纹增长原理计算转轴的剩余寿命。基于简化的齿轮动力学方程,Wang提出了在线估计、更新齿轮啮合刚度、啮合系数的方法。而后,根据齿轮啮合刚度和啮合系数的变化就可判断齿轮的健康状态,并预测其剩余寿命。不过,文献指出仍需要一个裂纹扩展模型来预测齿轮的剩余寿命。针对飞机发动机的油湿部件,Orsagh等讨论了信息融合框架,融合方法等一系列有关故障诊断和预测的概念。文献中以轴承为例讨论了轴承接触失效情况下的寿命预测方法,在碎裂(spall)萌生前使用Yu-Harris模型,在碎裂萌生后则使用碎裂扩展模型。同时,文献也讨论了将Yu-Harris模型和碎裂扩展模型随机化的方法,该方法考虑了载荷、不对中、润滑和制造过程中的不确定因素。Chelidze和Cusumano等将失效过程视为分级动态系统,其中损伤扩展的过程远慢于可观测的系统动态特性,并假设快速的、可观测的系统动态特性由某常微分方程控制,但并不需要知道该方程的具体模型。另一方面,若需进行寿命预测,则需要知道损伤扩展的数学模型。这方面类似的研究还有文献。文献可视为对文献的进一步扩展,将状态监测数据与具体的物理模型结合以提高寿命预测的精确度。以H-60直升机中级齿轮箱的主动齿轮为例,文献将疲劳失效分为裂纹形成与裂纹扩展阶段。在裂纹形成阶段,用Neuber原理估计裂纹萌生的循环数;在裂纹扩展阶段,则使用Paris公式计算裂纹扩展至失效的循环数。(这样的思路与文献中相似。)文献中建议,当基于状态监测数据的裂纹长度估计模型建立后,可用该模型取代Paris公式以估计裂纹扩展至失效长度的概率分布。Li和Lee提出了一个直齿轮寿命预测方法,该方法包括一个基于有限元法的齿轮裂纹长度估计模型,一个齿轮动力学仿真器(以估计齿轮载荷),一个计算应力强度因子的仿真器,和一个基于Paris公式的寿命预测模型。文献中方法考虑了在裂纹扩展过程中,齿轮刚度变化造成的动态载荷。试验证明,该方法具有较高的预测精度。若事先将齿轮载荷与应力强度因子在各种情况下的值制成查询表,文献中的方法则适用于现场应用。Marble和Morton就轴承碎裂扩展进行了一系列实验研究,他们指出碎裂形成的过程是多个小裂纹萌生和聚结的结果,而非单一主裂纹扩展的过程。因此,传统的裂纹扩展模型(例如Paris公式)并不适用于轴承。文献中提出了碎裂扩展公式,结合有限元方法计算公式中需要的碎裂周围的应力,进而预测轴承寿命。同时,文献也讨论了根据在线诊断结果更新模型的方法。Ramakrishana和Pecht根据线性损伤理论(Palmgren-Miner理论)估计印刷电路板上焊点在振动和温度载荷下的损伤累计,并以同样的方法估计焊点仅在热循环(thermalcycling)下的剩余寿命。试验结果表明,线性损伤理论能较好地估计一批样本的剩余寿命均值,但文献中未讨论如何对单个样本的剩余寿命做出较好地预测。类似地,Gu等用应变计检测振动载荷下印刷电路板的响应(弯曲曲率),以解析模型进一步计算焊点的应变,然后以Miner原理估计焊点损伤的累积并预测焊点的剩余寿命。类似的研究还有文献。基于物理的寿命预测方法,由于具备设备特定的物理学模型,往往不需要大量同类设备的历史数据即可获得较精确的预测结果。但是,该方法有时需要进行停机检查,这往往是不经济的甚至为生产所不允许。更重要的是,对于复杂设备,建立完备的物理学模型往往非常复杂。这就催生了直接基于设备的状态监测数据,而不使用具体的物理学模型的方法:基于经验的方法。1.2a在线预测设备剩余寿命n基于经验的方法,又称数据驱动的(data-driven)方法,试图直接从状态监测数据(及同类设备的历史数据)在线预测设备的剩余寿命。对基于经验的方法,在一些文献中有进一步的分类,但往往各小类之间有所重叠,如:基于状态外推的方法和基于人工智能的方法。我们认为,按照“监测状态-失效/剩余寿命”之间的不同联系,可将目前主要的基于经验的方法进一步分为:基于状态预测/外推的方法,基于统计回归的方法,和基于相似性的方法。(1)剩余寿命vd基于状态预测/外推的方法,认为设备失效可直接定义于状态空间,即可用确定性的失效阀值(单变量情况)或失效面(多变量情况)定义设备失效。在设备A运行的当前时刻t,基于状态预测/外推的方法试图预测其状态的进一步发展,并将预测的未来状态与失效阀值(或失效面)进行比较以预测设备的剩余寿命分布,并估计剩余寿命。以一个附有张力的钢带为例,Swanson等以Kalman滤波器跟踪钢带的固有频率。给定一个固有频率的失效阀值,Swanson等进而给出了预测失效时间及计算预测不确定性的模型。Engel等以多项式模型外推特征变量以预测直升机齿轮箱的剩余寿命。文献中提出了以寿命预测为目标,评估多个特征变量的一种方法,而目前文献中大多数评估的方法则以故障的早期发现或故障诊断为目标。文献中也指出,一个有效的寿命预测方法,不仅需要考察预期剩余寿命,也需考虑该预测的不确定性(以置信区间表达)。Cheng和Pecht讨论了将基于物理的方法和基于经验的方法融合的预测框架。其中基于物理的方法主要用于确定关键参数,确定潜在的失效机制并为其排序,确定失效模型,定义失效;基于经验的方法主要用于从监测变量中提取特征,建立健康基准,异常监测,预测参数的变化趋势等。从寿命预测的角度讲,文献中的方法主要还是基于经验的。Vachtsevanos和Wang介绍了基于动态小波神经网络(dynamicwaveletneuralnetwork)的预测方法,并以一个滚动轴承为例展示了该方法的可行性。示例中,以轴承振动信号的功率谱密度(PSD)为特征变量。在预测剩余寿命时,需假设基于PSD的失效阀值。Chinnam和Baruach指出在一些情况下无法在特征变量空间定义失效(即找到确定性的失效阀值或失效面),却可以根据专家知识(expertknowledge)或操作经验来定义失效。这样的情况下,他们建议采用focusedtime-laggedfeed-forward神经网络预测特征变量,采用Sugeno模糊推理模型来定义失效,进而预测设备的性能可靠性。文献中以钻头为例描述了其方法的可行性,但是限于其神经网络的预测步长较短而无法估计剩余寿命。但是,文献中对失效定义的处理方法仍有借鉴意义,结合能进行长程预测的预测模型,可以继续探讨在此失效定义下,剩余寿命预测的精确性。给定失效阀值,Orchard和Vachtsevanos使用粒子滤波方法预测行星载板的剩余寿命,并给出该预测的置信区间。试验结果显示文献中的方法具有对系统加载条件和特征变量信噪比的“鲁棒性”。文献中,Gebraeel等以Bayes方法在线更新指数衰退模型的参数,结合失效阀值得到设备剩余寿命分布的封闭形式,并将该方法应用于轴承的加速寿命试验中监测的振动幅值数据。文献中方法的推广、改进有文献。Gebraeel和Pan考虑了设备运行环境是时变的情形。Gebraeel等考虑了同类设备历史数据中失效时间数据可得,但没有失效过程中的状态监测数据的情形。Chakraborty等考虑了设备性能衰退模型的随机参数不呈正态分布的情形,发展了性能衰退模型的随机参数符合更普遍分布形式时设备剩余寿命分布的预测方法。由以上文献可见,基于状态预测/外推的方法往往采用各类时间序列的预测方法。因此,文献中大量其他的时间序列预测方法(例如,在线支持向量回归方法,指数平滑算法等)也可以在类似的框架下应用,并可能取得更好的预测效果。(2)基于hmm的热负荷寿命预测基于统计回归的方法,认为监测的状态是影响设备失效概率的因素。该方法往往需根据同类设备的历史数据建立起“设备失效概率,设备监测状态及设备运行时间”之间的函数(简称函数B)。在设备A运行的当前时刻t,基于统计的方法一般也需预测设备状态的进一步发展,并将该预测代入函数B,进而以设备失效概率分布的期望值预测设备A的剩余寿命。相比于基于状态预测/外推的方法,基于统计回归的方法往往无需定义确定性的失效阀值。Goode等指出许多设备(如:热轧钢厂的水泵)的失效过程可分为稳态过程和非稳态过程。其中,稳态过程对应于设备正常运行的阶段,而非稳态过程对应于发现故障到设备最终失效的阶段。Goode等建议用统计过程控制的方法区分这两个过程。在稳态过程,文献中使用可靠性数据(即稳态过程和非稳态过程长度的分布)进行寿命预测。进入非稳态过程,文献假设设备的状态监测数据呈指数型增长,并同时利用状态监测数据和可靠性数据进行寿命预测。当然,文献中方法的有效性很大程度上取决于模型的假设是否成立,以及状态检测数据的质量。Wang等将专家建议(关于水泵是否需要维护)分为四等,同时假设水泵轴承与造纸机轴承的寿命分布具有相同的形状函数(shapeparameter)和不同的尺度函数(scaleparameter),并以此计算水泵中关键轴承的风险函数,进而对单个水泵进行寿命预测。Kwan等提出了一个基于隐马尔科夫模型(HMM)的诊断与预测框架,并以轴承实验数据验证该框架的有效性。但是,文献中认为预测主要是失效后对设备状态变化过程的估计,而并未讨论如何在设备失效前预测其剩余寿命的方法。在文献的基础上,Zhang等提出了轴承失效诊断与预测的集成方法,该方法包括主成分分析,HMM,和一个自适应随机预测模型。其中,用PCA提取监测信号的主要特征,用HMM诊断轴承当前的状态(即:失效类型和当前处于HMM中最终失效状态的概率(健康指标)),用自适应随机预测模型预测并更新健康指标的发展趋势,在此基础上进行寿命预测。Baruach和Chinnam也建议使用HMM进行诊断与预测,对设备失效过程中的每个健康状态都单独建立一个HMM,再估计每个样本的状态转变点,进而估计状态转变点的条件分布,但文献中并未具体讨论如何预测剩余寿命。Camci和Chinnam将分级HMM(hierarchicalHMM)视为动态贝叶斯网络进行诊断和预测。相比于文献的方法,文献中的方法可以直接估计设备健康状态之间的转移概率。文献中假设每个健康状态的区间为指数分布,采用MonteCarlo仿真的方法估计设备的剩余寿命及其置信区间,并以钻头为例展示了其方法的有效性。类似的研究还有文献。Dong和He指出传统的HMM的不足在于隐含地以指数分布刻画各健康状态的区间长度,而这往往并不符合实际。因此,就有必要明确刻画各状态的区间长度。文献使用隐半马尔科夫模型(hiddensemi-Markovmodel(HSMM))以克服此不足之处。在水泵试验中,HSMM较传统的HMM取得更优的状态识别率。Peng和Dong在的基础上进一步考虑了状态区间长度与设备老化相关的情况,在HSMM中引入老化因子以刻画不同时刻各状态区间长度的分布。在水泵试验中,在接近水泵失效时,文献中的方法较文献中的方法能更精确地预测水泵的剩余寿命。Yan等建议使用logisticregression来建立设备的特征变量与其失效概率之间的联系,同时建议采用ARMA模型对特征变量进行预测,将预测值输入经训练的logistic模型进行剩余寿命预测。但是,文献中仅给出剩余寿命预测的示意图,并未给出具体的方法。Volk等建议使用比例强度模型(proportionalintensitymodel,PIM)预测设备的寿命,并以此评价预防维护的效果。文献中同时介绍了可修系统与不可修系统寿命预测的PIM。Banjevic和Jardine使用比例风险模型(proportionalhazardsmodel,PHM)预测设备的可靠性函数与剩余寿命。文献中将内部和(或)外部协变量的发展过程刻画为Markov过程,以此在某一时刻对协变量未来的发展趋势进行预测,并结合PHM预测设备的剩余寿命。文献中以运输机的转动机构为例展示了方法的主要概念。但是,将协变量发展过程刻画为Markov过程往往需要对一“连续”的协变量进行区间划分,进而以有限个数的“状态”描述协变量的发展过程。区间划分的得当与否将很大程度上决定文献中模型的预测效果。因此,在实际应用中,必须仔细寻找适当的区间划分方案(区间个数,各区间范围)。以轴承为例,Liao等比较了Logisticregression方法和PHM模型在预测单个设备剩余时的表现。文献中以监测振动信号的均方根值和翘度值作为PHM中的协变量进行寿命预测,试验结果表明PHM的预测精度好于Logisticregression模型。Ghasemi等认为PHM的协变量应为系统的真实状态,他们假设在一CBM计划中系统的状态不直接可见,但可以观察到一些特征信号,这些特征信号与系统状态随机相关。于是,隐马尔科夫模型(HMM)便可用于描述特征信号-系统状态的映射。在此基础上,将系统状态作为PHM的协变量用于估计系统的条件可靠性函数和预测剩余寿命。Wang假设设备的状态监测信号是其剩余寿命的函数(而与之前的状态监测信号无关),在此基础上利用随机滤波方法(stochasticfiltering)建立设备剩余寿命的概率密度函数。文献中以滚动轴承为例比较了真实的剩余寿命与预测的剩余寿命分布,结果较佳。根据文献的方法,Carr和Wang继而考察了系统有多个失效模式时的情况。文献中指出,在系统存在多种失效模式时,识别系统的失效模式,进而根据具体失效模式建立相应的寿命预测模型相对于不考虑失效模式的广义模型能去的更高的预测精度。Tian等建议同时利用失效和删失样本的数据来训练神经网络,进而预测设备的剩余寿命。文献中的神经网路,以设备的寿命(运行时间)和状态监测信号为输入,以设备的寿命百分比为输出。文献中以水泵轴承的振动信号数据验证了其方法较传统的不使用删失样本的神经网络能取得更精确、可靠的预测结果。(3)基于神经网络的寿命预测方法在正伺服系统剩余寿命计算上的应用基于相似性的方法,认为某正在服役的设备A的剩余寿命可预测为同类设备(又称参考设备)在某一时刻剩余寿命的“加权平均”。其中,权值根据设备A与各参考设备之间的相似性计算,而相似性则需进一步根据各设备在失效过程中的状态监测数据确定。相比于基于状态预测/外推的方法和基于统计回归的方法,基于相似性的方法有时无需进行状态预测。Zio和Maio提出了一个基于相似性的寿命预测方法。文献中以基于模糊逻辑的“距离值”刻画①正服役系统的特征变量与②参照集中系统的特征变量之间的相似程度,进而计算权重值,并以参照集中系统的剩余寿命的加权平均预测正服役系统的剩余寿命。文献以一核电系统的应用展示了改方法的有效性。试验显示,该方法具有较高的预测精度。类似的研究还有文献,但文献中使用了不同的刻画相似程度和加权平均剩余寿命的方法。以上两篇文献根据监测状态直接计算相似性,也有一些文献以监测状态估计设备当前运行时间等参数,进而根据估计值与实际值之间的误差计算相似性和权重值,总结如下。Gebraeel等采用轴承的损伤频率及其六次谐波处的振动幅值为特征变量,提出两种基于神经网络的寿命预测方法,取得较佳的预测精度。两种方法中,均以历史数据(训练集)为基础建立起多个前馈反向传播神经网络,估计某正服役轴承的服役时间,进而根据各神经网络对服役时间的预测误差确定权重值。不同的是,第一种方法将正服役轴承的寿命计算为训练集中各轴承寿命的加权平均;第二种方法则以加权平均的方法估计服役轴承指数性能衰减模型的参数,以进一步计算轴承寿命。文献中以一振动幅值的阀值来定义轴承失效。Huang等提出了一个基于自组织映射(self-organizingmap)的振动信号特征提取方法,在此基础上以文献中的神经网络方法预测球轴承的服役时间,进而预测其剩余寿命。Gebraeel和Lawley使用前馈反向转播神经网络估计轴承特征变量指数衰退模型的参数,结合失效阀值预测轴承的失效时间。同时,文献中以Bayes方法依据某轴承的在线特征变量更新剩余寿命的分布。与其他两类不使用更新的基准方法比较,文献中方法的平均预测误差较小。2维护规划模型回顾Prognostics的目的就是要为维护方案的制定提供决策支持(decisionsupport)。因此,很自然地,需要在设备状态预测的过程中考虑下一步的维护方案。上一部分中总结的设备寿命预测的方法提供了设备剩余寿命的估计,而维护规划主要是在此基础上进一步确定恰当的执行维护的时间。可见,比起寿命预测,维护规划更直接地决定了维护策略的有效与否,更直接地影响了一个企业的产能。但是,维护规划问题往往较寿命预测更为复杂。总结目前相关文献,一个预防维护规划模型一般包含五大因素,即:①目标,如:维护费用最小化,维护费用率最小化,系统平均可用度最大化等;②维护方案,刻画了维护执行的特点(框架),如:周期性(periodic)维护,顺序型(sequential)维护,控制限度(control-limit)维护等;③维护质量(或维护效果),刻画了维护行动将设备性能提升的程度,有:修复如新(perfectmaintenance),修复非新(imperfectmaintenance),修复如旧(minimalmaintenance)等;④退化模型,刻画了设备退化的过称,有:传统的寿命分布(lifetimedistribution),Gamma过程,马尔科夫过程等;⑤维护限制,如:部件间的维护冲突,维护零件有限等。下面总结近年来主要的相关模型。基于状态监测的维护规划模型分类方法多样,例如,Jardine等将之分为完全可见系统和部分可见系统。其中,对于完全可见系统,状态监测变量完全体现了系统的真实状态;对于部分可见系统,状态监测变量不能完全体现系统的真实状态,但与之相关。我们将各种维护规划模型分类为离线预防维护规划模型和在线预防维护规划模型。2.1实际运行设备状态监测数据离线预防维护规划模型指依据反映系统群体特征的衰退信息建立的维护模型。维护策略一旦确定(如:确定最优维护状态阀值),并不随某一实际运行设备的实时监测数据改变。在应用离线预防维护规划模型时,某一实际运行设备的状态监测数据往往用于观测其状态值是否达到了先前确定的最优维护状态阀值。若未达到,则允许设备继续运行;一旦达到,则执行预防维护。目前文献中的离线预防维护规划模型,主要使用的估计设备失效概率的衰退模型有:PHM,几类随机衰退/损伤累计模型,马尔科夫衰退过程等。(1)状态监测和周期状态下的换重要Makis和Jardine假设系统的失效概率不仅依赖于其工作时间(age)也依赖于一诊断随机过程(diagnosticstochasticprocess),并使用PHM描述系统的失效概率,在此基础上建立了一个替换策略以最小化长期期望费用率。文献中,Makis和Jardine假设只在离散时刻可以观察到系统诊断随机过程的值,并以阶跃函数近似计算离散观测点间的随机过程。Kumar和Westberg提出了一个基于状态监测的替换策略。其中,用PHM确定系统失效是否与监测变量有关,进一步使用PHM估计系统的可靠性,在此基础上建立系统的预防维护模型以求取最优替换时刻。值得指出的是,文献中的监测变量指的是压力等常值环境变量。Banjevic等考察了一个基于状态监测的控制限度替换策略,以最小化长期期望费用率。文献中以PHM估计系统的失效率,同时假设监测的状态符合非齐次Markov过程,进而估计系统的失效概率以估计长期期望费用率,最后求得最优维护限度。文献中也讨论了收集、预处理和使用实际油样、振动信号的一些经验。类似的基于PHM的替换策略及其应用还有文献。Ghasemi等提出了一个基于周期状态监测的替换模型以最小化长期期望费用率。文献中考察了监测变量不能直接反映系统真实状态,但与系统状态随机相关的情况,以PHM估计系统的失效概率,进而建立系统费用率模型。Ghasemi等在文献的基础上进一步考察了非周期监测的情况。(2)状态状态下的维护策略Barbera等提出了一个基于状态检测的预防维护模型,以最小化长期期望费用。在该模型中,预防维护在系统状态值达到某一阀值时进行,并且预防维护可完全恢复系统的状态。文献中假设系统的状态为非减函数,在任一时刻,系统损伤的增长为一与累计损伤无关,且概率密度函数已知的随机变量,同时系统的可靠性为系统状态的指数函数。在文献中,使用与文献中类似的建模方法,Barbera等又考查了一个由两个串联部件组成的系统级维护模型。Shahanaghi等在文献的基础上考察了修复非新的情况,认为修复后系统的状态随机地处于全新与修复前的状态之间。Aven以计数过程方法建立了一个基于状态监测的替换策略,以最小化长期期望费用率。在文献的模型中,预防维护在系统失效率达到某一阀值时进行,替换之后系统恢复全新的状态。文献中假设在某一时刻系统的失效率是其状态的确定性函数,但并未做具体的讨论。Hontelez等考察了漂移衰退过程的基于状态监测的维护策略以最小化一段时间内的平均维护费用。文献中将漂移衰退过程分成(N+1)段,并建模成马尔科夫过程。文献中使用控制控制限度维护策略,假设维修之后系统恢复全新的状态,确定了最优监测时间间隔与最优维修时间。类似的研究有文献。Wang提出了一个基于状态监测的替换模型。文献中基于随机参数增长模型估计系统在某一时刻的失效概率,在此基础上进行建模,进而确定最优监测区间和最优维护限度,以优化系统的长期期望费用率,失效时间率或可靠性。Grall等提出了一个针对随机衰减系统的基于状态监测的替换模型,求取最优维护阀值和最优监测时间间隔,以最小化系统的长期期望费用率。文献中假设系统的状态衰退在任一时刻非负,且状态衰退增长程度的概率密度函数已知,同时系统的确定性失效阀值已知,在此基础上估计系统在某一时刻的失效概率。值得指出的是,文献中考虑了监测时间间隔非常数的情况,根据系统的状态实时地确定下一个最优监测时间间隔。类似的研究还有文献。文献的基础上的进一步的研究有文献。Deloux等考察了单部件系统存在两种失效形式的情况,其中一种为累计损伤失效,另一种为应力环境下的冲击失效,进而优化系统长期期望费用率。Zhao等考察了系统状态非单调衰退且受环境因素影响的情况下的最优监测、替换策略,以最小化长期期望费用率。文献中以类似PHM的方式刻画环境因素对系统状态衰退的影响。Nakagawa和Ito考察了多部件系统的情况,他们假设系统受到的外界的随机冲击(各冲击间隔长度符合同一概率函数,各冲击造成的损伤大小而造成损伤累积符合同一概率函数),当累积的损伤超过一定值时系统失效。在次基础上,他们建立了系统预防维护模型,计算最优维护阀值(累积损伤值)以最小化长期期望费用率。Wang等提出一个基于状态监测的订货-替换模型,考虑了维修部件需提前预定的情况。文献中假设系统衰退符合随机衰退模型,衰退增量单调递增,各时刻增量互相独立且符合同样的统计分布。文献中考察了系统的长期期望费用率,预防替换率和系统平均可用度指标,根据这些指标以优化订货限度,预防替换限度和监测周期长度。Wang等考察了基于马尔科夫衰退过程的预防替换模型及另一种库存机制。Wang等在文献的基础上考察了系统失效阀值不确定的情况。Liao等提出了一个连续监测的单部件系统的预防维护模型,寻找最优控制限度(定义于系统的状态空间)以最大化系统平均可用度。文献中基于随机衰退模型(如Gamma过程)计算系统的失效概率,认为系统在预防维护后的状态分布于全新与修复前的状态之间,且其密度符合某概率函数。(3)状态阀值的确定Marseguerra等考察了一个多部件串并联系统的基于状态监测的维护模型。该模型以同时优化系统利润与可用度为目标,在系统的马尔科夫衰退模型中考虑了非完全维护因素,以遗传算法确定各部件的最优维护状态阀值。为了描述更复杂的情况,如系统衰退与外界载荷相关,维护人员有限等情况,文献又建议使用MonteCarlo仿真的方法。在系统马尔科夫衰退模型已建立的情况下,文献中模型针对的主要是系统设计阶段,而未考虑实际运行中对某一个系统的实时维护。类似的研究还有文献,但文献中采用最小化一段时间内的期望费用为目标,并采用随机衰退模型。基于一个连续-离散随机衰退模型,Makis和Jiang提出了一个基于状态监测的优化框架以最小化长期期望费用率。文献中,隐藏的设备状态发展由连续马尔科夫过程描述,监测变量与设备状态有关且由离散随机过程描述。(4)失效概率的动态特性Zhou等提出了一个基于状态监测的单部件系统可靠性限度维护策略,以最小化长期期望费用率。文献通过混合概率模型刻画修复非新的预防维护效果。但文献中未详述如何刻画系统的衰退过程以求取其在某一时刻的失效概率,未考虑到实时运行系统失效概率的动态特性。Camci提出了一个多部件系统的基于状态监测的维护模型,考虑了一些维护限制(如部件之间的维护冲突,维护资源有限等),以最小化一段时间内的系统期望维护费用。文献指出,由于系统各部件之间的联系,优化各部件的维护策略不一定导致最优系统维护策略,因此就有必要进行系统级的维护时序规划。但是,文献中假设系统的失效概率已知,并未考虑实际运行系统不断变化的状态。类似的研究还有文献。2.2基于状态空间模型的仿真研究与离线预防维护模型不同,在线预防维护模型依据某一实际运行系统的在线状态监测数据不断更新最优预防维护策略,并在恰当的时机执行预防维护。在线预防维护模型针对每个设备制定不同的预防维护策略,而在离线预防维护模型中所有设备均使用同一预防维护策略。Christer等发展了一个感应电炉的感应器的替换策略。文献中,由于感应电炉的结构限制,感应器裂纹的尺寸无法直接测量,文献使用可以测量的“传导率”为特征变量,通过状态空间模型估计实际裂纹尺寸。文献中采用卡尔曼滤波方法预测和更新感应器的失效概率,在此基础上计算最优替换时间,以最小化期望费用率。Pedregal和Carnero描述了一个案例研究,对一个石