基于Gamma退化过程的剩余寿命预测及维修决策优化模型研究

本资料包含pdf文件1个，下载需要1积分

Study on Residual Life Prediction and M aintenance Decision--M akingOptimal M odel Based on Gamma Deterioration ProcessGu Yu-bo.Jia Yun～xian.Zhang Ying-bo(Department of Equipment Command and Management，Ordnance Engineering Colege，Shijiazhuang 050003，China)Abstract：By studying the failure mechanism of the product，the state space model is established based on Gamma de-terioration process.The model parameters are solved with the EM -PF parameter estimation method，and the concreteform of model is confirm ed.The distribution function and probability density function of remaining life are obtained。

The prediction method combines the deterioration state with failure threshold，which reduces the eror.The wear char-acteristics of bearings are modeled；the correlations between state inform ation and remaining life are completed.Ac-cording to the probability density function of remaining life，a maintenance decision-making model is founded up，whose target is to get the most proper replacement time，optimizing the maintenance decision-making.At last，the da-ta of the bearing life experiment is taken for an example to verify the model，and the result shows that the model is prac·tical and effective。

Key words：roling bearing；Gamma deterioration process；residual life；maintenance decision-making机械设备在使用过程中，某些部件的运行状态会逐渐退化。状态的退化可以表现为机械部件的磨损、裂纹的增长、腐蚀程度的加深等，这些都是-系列理化效应 的结果。产 品状态退化 的收稿日期：2012-07-23；修回日期：2012-09-28基金 项 目：十 二 五 ”武 器 装 备 预 先 研 究 项 目(51327020101)作者简介：谷玉波(1987-)，男，河北新乐人 ，硕士研究生，研究方向为维修工程理论与应用。

最终结果会导致功能故障的发生，而且在产品的退化过程中能够直接表征退化状态的状态指标(磨损、裂纹、腐蚀等)经常是不可测量或难于直接进行测量的。为获得设备状态退化程度的数据，可以运用 CBM状态监控技术对设备的状态数据进行提取分析，间接地反映出设备的降状态，实现对设备剩余寿命的预测，最终实现维修决策的优化，使维修效果达到最佳。在此，运用 Gamma退化过程状态空间模型来描述设备状态性能的退化，通过该模型实现了剩余寿命的预测，并运用剩谷玉波，等：基于Gamma退化过程的剩余寿命预测及维修决策优化模型研究 ·45·余寿命预测结果对维修决策模型进行优化。 2 Gamma-SSM 退化模 型1 Gamma退化过程的状态空间模型 2. 1 模型的建立在常用的分布中，Gamma分布可以用来描述连续累积磨损过程，并且具有退化过程建模所需的所有属性，即非 负、平稳增 长和独立增量过程。因此，运用 Gamma分布建立状态空间模型，可以描述设备的退化过程 ，并实现剩余寿命的预测。

1.1 状态空间模型状态空间模型(State Space Model，SSM)是-种用来描述系统演化过程的模型，可以用来建立产品的状态退化模型。状态空间模型假设系统状态随时间的演化可以由-个不可观测或难以观测的状态序列确定，与状态序列伴随的是-个可以观测的测量序列，两者之间的关系可以由状态空间模型来描述，其中应用较为普遍的状态空间模型是典型相关(Canonical Corelation)方法 J。通过建立状态方程和观测方程，状态空间模型为描述系统的动态演化特性建立了-个有效的模型框架，-些复杂的问题就可以用比较简单的形式表示。状态空间模型的基本方程为1F ( ， ， ， 1)， (1)Y H ( ，u ， ， 2)， (2)式中： 为状态转移函数， 为测量函数，-般而言 和 假设是已知的； 为系统的输入变量或控制变量；0 ， 分别为状态方程和观测方程的静态参数； 和 分别为系统噪声和观测噪声，-般假设 和 8相互独立，且服从某-特定分布。

1.2 Gamma过程Gamma过程是-个具有独立、非负增量的随机过程。运用 Gamma过程对-些退化过程建立模型，能够使建模过程 中的数学计算更加明确易懂 J，适合描述在时域上 单调累积微小增量的退化过程，比如磨损、腐蚀、裂纹增长等。

根据 Gamma过程，定义非负随机过程 (t)；t≥0，其概率密度函数为) e ， ≥0；A， >0。

IX (3)其中，F(IX)J Xa-1e-dx，则称 服从参数为J 0( ，A)的 Gamma分布， 为形状参数，A为尺度参数。简记为X(t)～Gamma(t l ，A)。其期望和方差分别为：E( )a/A，D( )IX/A 。

假设设备状态 的退化符合 Gamma过程，Y表示与 相对应的观测变量， 表示设备 t 时刻所对应状态，则建立的状态空间模型的状态方程和观测方程为- - 1～Gamma(Ol(t )-O/(t -1)，A)，(4)Y H ( ) 。 (5)假设当系统状态 达到故障阈值 墨时系统发生故障，如图 1所示。则系统状态从 0首次达到故障阈值所需的时间为Tfinft： Xf，t>0。 (6)时间图 1 系统退化过程为简化问题，假设观测噪声是可以加的。假设设备状态退化过程为平稳 Gamma过程，形状参数 (t)为线性函数，即Ol(t)0·t，并且假设观测量 Y与状态量 之间的关系为 YC· s，则在获得观测序列Y 的条件下 ，退化模型可以表示为1 - ～ Gamma(Ol·(t 1- )，A)， (7)YiC· 。 (8)由假设可知，观测噪声 服从均值为 0的正态分布。所以，Y～N(C· ， )。因此只要获得参数 。，A，C， 的值，就可以确定该退化模型的形式。

2.2 模型参数估计在对上述模型中参数进行求解的过程中，会出现新的需要估计的参数，导致运算量增大。为解决此问题，运用文献[5]提到的经验最大化(Ex。

cle Filtering，PF)算法相结合的方法求解模型中出现的参数。

EM算法是进行极大似然估计的-种有效方法，简单来说，这-算法的每-次迭代主要包括2· 46· 《轴承)2013.No.4个步骤：(1)求期望(Expectation Step)，称为 E步；(2)求极大值(Maximization)，称为 M步。EM算法通过假设潜在变量的存在，来简化似然函数极大估计 。

PF算法是-种基于 Monte Cado方法和递推Bayes估计的统计滤波方法 J，建立在序贯重要性采样(Sequence Importance Sampling，SIS)和 Bayes理论的基础上，对于解决非线性、非 Gauss问题有很好的效果 。PF算法的基本思想是：首先根据系统状态向量的经验分布在状态空间产生-组随机样本集合(即粒子)，然后根据观测数据不断调整粒子的权重和位置，通过调整后的粒子信息修正最初的经验分布。其实质是用由粒子及权重组成的离散随机测度近似相关的概率分布，并根据算法递推更新离散随机测度。当粒子容量足够大时，就近似于状态变量真实的后验概率分布。

由于EM-PF算法兼有两种算法各 自的优点，可以更加有效地对模型中的参数进行求解。

EM算法 PF算法-获得初始参数8 I初始化粒子群及权重I。 。 。 。 。 。 - - - -- t' ! lE步 -1获得完全似然函数期望酒 蓟 剩劐凄采 lQ(口l0 )lM步l计算新参数 ” - I N 计算 )， ) )， )1/ N / 是否满定、<-/ 收敛条件 / / Y 获得最优参数 IY l00(㈩) I - 输出E(x。 )，以 )2)，E(1nv图2 模型参数求解过程3 剩余寿命预测剩余寿命是指产品从被检测的某-时刻起到该产品发生故障的时间长度 ，设备剩余寿命的预测是维修管理中的-个重要环节，由于剩余寿命预测对产品使用过程中的安全性、经济性和任务性具有重要的影响，因此它是制定维修决策的重要依据。

通过确定状态空间模型参数，就能得到 Gam-ma状态空间模型的具体形式，进而得到产品的剩余寿命分布函数。令 为粒子个数； 为设备故障时间；W 为第 i个粒子的权重；X 为故障阈值 ；： 为直至当前时刻所获得的观测序列； ： 为当前获得的观测值序列；t 为第 i次状态检测时间；为 t 时刻的状态值； 。： 为 n次检测所获得的状态序列，即 1. 1， l， 2，， ；Y 为t 时刻的退化量观测值；Y 为n次检测所获得的观测序歹 ，即Yo. Y。，Y。，Y：，，Y ，贝0 t 时亥0的乘0余寿命 丁 分布函数为F(丁 I Yo： )P( t ≥ I Yo： )Ip(tr t ≥Tf I )·p(x I Yo： )dx ， (9)P( t ≥Tf l )p(x( t )≥Xf l ， (10)p(x l y )f(x 0-f I Yo： )-∑ s 6( -；)。 (11)另外，预测 t 时刻的状态概率密度函数可表示为p(x l Y0： )Jp( l )P( l Yo： )dx N∑p(x )· 。 (12)由粒子滤波算法过程可知，(11)式可通过滤波粒子 ；s1，2，， 进行表示，则剩余寿命分布函数及概率密度函数为F( I Yo： )P(r t ≥ l y0： )Jp(丁 Nt ≥ I )p( l Yo： )dx ∑p(z ti≥1，(。 ， )- Ⅷs ，(13)肌 l Yo：i) 。 (14)在Y。. 条件下，平均剩余寿命 为uL I (1-F(r I Y ))dt。 (15)4 维修决策优化模型构建剩余寿命概率密度函数的最终目的是根据剩余寿命与实际需求做出维修决策，确定合理的维修或更换时间。平时对设备进行维修管理的过程中，人们总是希望设备的维修管理费用旧能达到最低。以费用最小为目标的维修决策模型可以根据设备当前的运行状态确定在下-问隔期谷玉波，等：基于Gamma退化过程的剩余寿命预测及维修决策优化模型研究内是否进行维修以及何时进行维修，可最大程度地避免故障的发生。

在维修策略中，产品有2种维修方式：(1)修复性维修，即产品在达到预定维修时刻之前发生故障而进行的维修或更换活动，也称故障维修或故障更换；(2)预防性维修，即产品在达到预定的维修时刻进行的维修或更换活动。因此按照维修决策的目标，并根据剩余寿命概率密度函数，建立以最小费用为 目标的维修决策模型，其决策过程如图3所示。

当T) - 下- 检测时刻 时刻- 图3 以费用最小为目标的维修决策示意图4.1 假设和参数(1)每隔 进行-次状态检测并获取当前时刻的状态数据。

(2)如果在检测间隔期内发生故障，则立刻进行修复性维修。

(3)设备的预防性维修费用小于故障后的修复性维修费用，即

(4)预防性维修和修复性维修可以使设备恢复初始状态，即修复如新。

(5)参数：E(c)为更新周期内的期望总费用；E(T)为期望更新周期长度；At为状态检测间隔期，Att -t ； I y0： )为t 时刻状态信息为Y。： 时产品的剩余寿命概率密度分布函数；P为部件进行预防性维修的概率；P 为部件进行修复性维修的概率；c。为预防性维修费用；c 为修复性维修费用；TR为最佳维修更换时间；C(Ta)为部件维修更换时间为 时，长期使用下的单位时间费用；t 为当前时刻的状态监控点，t。0；为预防性维修所需的平均时间；To为修复性维修所需的平均时间；c 为每次实施状态检i见0的费用。

4.2 以费用为目标的决策模型设备部件在长期使用下的单位时间费用可表示为c( ) ， (16)E(C)cfPfcpPpnC coC(r < -t I Y0： )cp[1- < -t l Yo： )]nC Cp(cf-cp) 下i

J 0(19)将(17)-(18)式代入(16)式中，则费用模型为c( ) [cp(Cf-Cp)· r < Ta-t I Yo： )nC ]/t ( -t )·[1- <- t J：)] (r To)f(丁 l y0： ) 。

0(20)其中， ci为更新周期内监控所产生的费用，由于状态检测中检测费用大多为设备投入费用， c 常成以忽略不计，因此在实际决策中，近似取c 0。如果 和 相对于 来说很小，在费用模型中也可以将其忽略。

利用数值方法最携(2O)式，得出设备进行基于状态的维修时可达到的最小单位费用及所对应的维修时刻 。当获取新的状态信息值 ：时，需要将数据代入(20)式重新计算更新结果。

4.3 最佳维修时间确定方法通过费用决策模型可以得到各检测时刻对应的维修时间 ，使得单位时间维修费用最校当- t >At时，不对部件进行维修，直到下次正常检测时刻；当首次出现 TR-t ≤At时，则到达时对部件进行维修或更换，此时的 即为基于状态维修的最佳维修时间。

5 实例分析试验采用3套 Rexnord ZA-2155双列滚子轴承，对2套轴承进行全寿命试验，用于模型参数估计和模型验证，对第 3套轴承进行截尾寿命试验，用于维修决策，分别记为轴承 1、轴承2和轴承3。

全寿命试验过程中，轴承的转速和载荷保持在2 000 r/min和10 000 LB(约44 500 N)，并采用PCB加速度传感器对轴承的振动信号进行采集，每1 h采集1次信号，采样频率为20 kHz，长度为1 s。

试验过程中，轴承 1和轴承 2分别在 966 h，982 h出现强烈的噪声，此时认为轴承故障。用振《轴承)2013.No.4动信号的能量均值作为间接状态观测数据Y，磨损量为状态数据 ，2组全寿命数据的振动信号能量均值在整个寿命周期上的变化如图4所示。

0 l00 20o 3O0 400 500 600 700 300 900 l 000时间/h图4 试验 轴承 平均 能量 的变化过程应用轴承 1的间接状态观测数据对模型参数进行估计，并采用 EM-PF算法对剩余寿命模型中的参数进行计算，结果见表 1。

表 1 模型参数计算结果查墼 ! !望笪 ： ! ： ! 箜 !! !确定模型参数之后就能得到基于 Gamma退化过程的状态空间模型，通过计算求得轴承剩余寿命概率密度函数 l ： )。对于轴承2的观测数据，运用 Matlab 7.0对剩余寿命预测方法进行实现，不同状态检测时刻(500-980 h)的剩余寿命概率密度函数如图5所示 ，其中 表示剩余图5 不同检测时刻剩余寿命概率密度寿命的估计值， 表示剩余寿命真实值。剩余寿命预测值与实际值见表2。

表 2 剩余寿命预测值与实际值比较(轴承 2) h堕型盛 !实际值 482 402 322 242 162 82 2堡型笪 丝! !： 兰 堑由上述结果可以看出，该模型对于此类试验背景下的轴承剩余寿命预测具有-定的实用性，下面利用该模型对轴承 3进行维修决策，确定使得单位时间费用最小的最佳维修时间。假设 c 600元，c 1 200元，在得到轴承的剩余寿命概率密度函数之后，将 J Yo. )代人到(2)式中，计算各维修时刻的单位时间费用，得到每次检测时费用最低的更换时间 71R，如图6所示。图中可以看出，在某-检测时刻，随着更换时间的增加，单位时间费用先减少后增加，在极值点处可以取得该次检测的最佳维修时间。图中 表示各检测时刻使单位时间费用最低的维修时间。

暖糕厘茁瓣图 6 不 同检 测时刻所确定的维修 更换时间结合最佳维修时间确定的方法，通过 Matlab计算，可以得出首次满足 Ta-t ≤At条件时，T 701 h。所以在 701 h对轴承进行更换，可以使单位时间费用 C(T )达到最低，约为289.6元。

6 结束语为解决维修决策的实际问题，重点研究了剩余寿命预测方法和维修决策优化模型。首先根据产品状态的退化过程的特征建立了基于 Gamma退化过程的状态空间模型，并对模型参数的估计方法进行了说明，进而通过该模型实现剩余寿命的预测，得到产品的剩余寿命概率密度函数。最后以滚子轴承寿命试验得到的数据为例，验证了剩余寿命预测模型，并以费用最小为 目标建立了维修决策模型，对另-轴承进行了维修决策的优化，实现了剩余寿命预测与维修决策的有机结合，0 O O O 0 O 0 O O O 删 露 删∞l露 m m m稍讲 捕《藻