本发明涉及可靠性工程技术领域,特别是涉及一种考虑状态检测影响的设备剩余寿命预测方法及系统。
背景技术:
为研究工业设备和军事装备的可靠性和剩余寿命,学者们和工程技术人员侧重于研究基于性能退化数据的方法。随着传感器的不断升级,性能退化数据(陀螺仪漂移、li-ion电池电容值等)大多通过状态检测技术获得。然而由于技术或经济成本的限制,实际工程中的状态检测通常并不是连续的,仅在一些离散的时间点对设备进行状态检测。对一些特定的随机退化设备,状态检测可能会改变设备的运行应力,从而影响设备的退化过程。然而,现有的方法并未在退化建模及剩余寿命预测的过程中充分考虑状态检测对设备退化过程的影响。
不同设备的状态检测技术各不相同,因此,状态检测对设备退化过程的影响也存在差异。事实上,状态检测既有可能加剧设备的运行应力也有可能释放设备的运行压力,从而改变设备原有的退化过程。以惯性导航设备中机械陀螺仪为例,在其长期贮存过程中,为保证其性能指标满足使用需要,需要周期性地对陀螺仪进行状态检测。这些状态检测需要对陀螺仪通电与加温,从而引入额外的电应力,并改变陀螺仪所处的温度场。此外,状态检测引入的额外电应力与温度应力将会影响陀螺仪的性能退化过程。再如,li-ion电池在电动汽车与混合电动汽车中应用广泛,其性能在使用过程中也存在类似的现象。作为li-ion电池的重要性能指标,电容值通常基于对充/放电过程的记录计算得到。研究表明,除了随着贮存时间累计,电池的电容会有所下降,充/放电循环也会加速电容的损耗,从而消耗电池的寿命。针对含有此类状态检测的退化设备,必须将状态检测的影响融入到退化建模过程中,从而提高退化模型的合理性与剩余寿命预测的准确性。为此,需要解决两个问题:一是如何描述每次状态检测对设备退化过程的影响;二是如何将状态检测的影响融入到设备的退化模型与剩余寿命预测中。本发明假设相较于两次状态检测之间的时间相隔,单次状态检测的持续时间可忽略不计。实际设备大多在离散的时间点进行检测,因此该假设是合理的。同时,受设备复杂性与运行过程不确定性等因素的共同影响,单次状态检测对设备退化过程的影响也具有一定的随机性。因此,可将状态检测对设备退化过程的影响描述为退化过程中的随机冲击。需要说明的是,根据每次状态检测的间隔是否一致,工程中可将状态检测分为周期性状态检测和随机状态检测,本发明重点关注周期性状态检测,如无特殊说明,状态检测均指周期性状态检测。周期性检测中,相邻两次检测的时间间隔相等且为固定值,但单次状态检测对设备退化过程的影响是随机的。
技术实现要素:
本发明的目的是提供一种考虑状态检测影响的设备剩余寿命预测方法及系统,采用幅值服从正态分布的随机冲击刻画单次状态检测对设备退化水平的影响,有效克服了现有研究难以描述检测活动与退化过程间定量关系的不足。
为实现上述目的,本发明提供了如下方案:
一种考虑状态检测影响的设备剩余寿命预测方法,包括:
采用线性wiener过程描述法建立设备自身的退化模型;
根据所述设备自身的退化模型采用幅值服从正态分布的随机冲击刻画单次状态检测法,建立状态检测影响下的退化模型;
采用转换阈值方法,建立状态检测影响下的设备剩余寿命的分布模型;
根据所述状态检测影响下的退化模型,得到状态空间模型;
采用最大期望算法对所述状态空间模型的参数进行估计,得到估计后的状态空间模型;
根据所述估计后的状态空间模型对所述设备剩余寿命的分布模型进行自适应更新,得到更新后的设备剩余寿命的分布模型;
根据所述更新后的设备剩余寿命的分布模型对设备的剩余寿命进行预测。
可选的,所述采用线性wiener过程描述法建立设备自身的退化模型,具体包括:
采用线性wiener过程描述法建立设备自身的退化模型x′(t)=ηt+σbb(t)
其中,x′(t)为t时刻设备的退化水平;b(t)为标准布朗运动;η与σb分别为wiener过程的漂移系数与扩散系数。
可选的,所述根据所述设备自身的退化模型采用幅值服从正态分布的随机冲击刻画单次状态检测,建立状态检测影响下的退化模型,具体包括:
根据所述设备自身的退化模型采用幅值服从正态分布的随机冲击刻画单次状态检测,建立状态检测影响下的退化模型
其中,b(t)为标准布朗运动;η与σb分别为wiener过程的漂移系数与扩散系数,
可选的,所述采用转换阈值方法,建立状态检测影响下的设备剩余寿命的分布模型,具体包括:
采用转换阈值方法,建立状态检测影响下的设备剩余寿命的分布模型,所述分布模型包括概率密度函数模型和累积分布函数模型;
所述概率密度函数模型通过下式表示:
所述累积分布函数模型通过下式表示:
其中,η为wiener过程的漂移系数,σb为wiener过程的扩散系数,μkl为时间区间[tk,tk+l)内状态检测引起的退化水平变化量之和的均值,
可选的,所述根据所述状态检测影响下的退化模型,得到状态空间模型,具体包括:
根据所述状态检测影响下的退化模型,得到状态空间模型
其中,ηk为设备在第k次状态检测的瞬时退化速率;υk为第k次状态检测导致的设备性能突变,γk+1为第k+1次状态检测导致的设备性能突变,δxk为第k次状态检测与k+1次状态检测之间设备退化状态的增量;εk为随机效用,
可选的,所述根据所述估计后的状态空间模型对所述设备剩余寿命的分布模型进行自适应更新,得到更新后的设备剩余寿命的分布模型,具体包括:
根据所述估计后的状态空间模型对所述设备剩余寿命的分布模型进行自适应更新,得到更新后的设备剩余寿命的分布模型,所述更新后的设备剩余寿命的分布模型包括更新后的概率密度函数模型和更新后的累积分布函数模型;
所述更新后的概率密度函数模型通过下式表示:
所述更新后的累积分布函数模型通过下式表示:
其中,
一种考虑状态检测影响的设备剩余寿命预测系统,包括:
第一退化模型建立模块,用于采用线性wiener过程描述法建立设备自身的退化模型;
第二退化模型建立模块,用于根据所述设备自身的退化模型采用幅值服从正态分布的随机冲击刻画单次状态检测法,建立状态检测影响下的退化模型;
第一分布模型建立模块,用于采用转换阈值方法,建立状态检测影响下的设备剩余寿命的分布模型;
第一状态空间模型确定模块,用于根据所述状态检测影响下的退化模型,得到状态空间模型;
第二状态空间模型确定模块,用于采用最大期望算法对所述状态空间模型的参数进行估计,得到估计后的状态空间模型;
第二分布模型建立模块,用于根据所述估计后的状态空间模型对所述设备剩余寿命的分布模型进行自适应更新,得到更新后的设备剩余寿命的分布模型;
预测模块,用于根据所述更新后的设备剩余寿命的分布模型对设备的剩余寿命进行预测。
可选的,所述第一退化模型建立模块,具体包括:
第一退化模型建立单元,用于采用线性wiener过程描述法建立设备自身的退化模型x′(t)=ηt+σbb(t),
其中,x′(t)为t时刻设备的退化水平;b(t)为标准布朗运动;η与σb分别为wiener过程的漂移系数与扩散系数。
可选的,所述第二退化模型建立模块,具体包括:
第二退化模型建立单元,用于根据所述设备自身的退化模型采用幅值服从正态分布的随机冲击刻画单次状态检测,建立状态检测影响下的退化模型
其中,b(t)为标准布朗运动;η与σb分别为wiener过程的漂移系数与扩散系数,
可选的,所述第一分布模型建立模块,具体包括:
第一分布模型建立单元,用于采用转换阈值方法,建立状态检测影响下的设备剩余寿命的分布模型,所述分布模型包括概率密度函数模型和累积分布函数模型;
所述概率密度函数模型通过下式表示:
所述累积分布函数模型通过下式表示:
其中,η为wiener过程的漂移系数,σb为wiener过程的扩散系数,μkl为时间区间[tk,tk+l)内状态检测引起的退化水平变化量之和的均值,
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明提供一种考虑状态检测影响的设备剩余寿命预测方法,通过采用幅值服从正态分布的随机冲击刻画单次状态检测对设备退化水平的影响,有效克服了现有研究难以描述检测活动与退化过程间定量关系的不足。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例考虑状态检测影响的设备剩余寿命预测方法流程图;
图2为本发明实施例陀螺仪的性能退化数据;
图3为本发明实施例参数a的估计过程;
图4为本发明实施例参数σ的估计过程;
图5为本发明实施例参数θ0|0的估计过程;
图6为本发明实施例参数p0|0的估计过程;
图7为本发明实施例参数
图8为本发明实施例滤波结果;
图9为本发明实施例观测值与估计值的对比;
图10为本发明实施例不考虑检测影响且更新参数时剩余寿命的概率密度函数;
图11为本发明实施例考虑检测影响且更新参数时剩余寿命的概率密度函数;
图12为本发明实施例不考虑检测影响且未更新参数时剩余寿命的概率密度函数;
图13为本发明实施例考虑检测影响但未更新参数时剩余寿命的概率密度函数;
图14为本发明实施例剩余寿命预测的均方误差的对比;
图15为本发明实施例考虑状态检测影响的设备剩余寿命预测系统结构图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种考虑状态检测影响的设备剩余寿命预测方法及系统,采用幅值服从正态分布的随机冲击刻画单次状态检测对设备退化水平的影响,有效克服了现有研究难以描述检测活动与退化过程间定量关系的不足。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明实施例考虑状态检测影响的设备剩余寿命预测方法流程图。如图1所示,一种考虑状态检测影响的设备剩余寿命预测方法,包括:
步骤101:采用线性wiener过程描述法建立设备自身的退化模型;
步骤102:根据所述设备自身的退化模型采用幅值服从正态分布的随机冲击刻画单次状态检测法,建立状态检测影响下的退化模型;
步骤103:采用转换阈值方法,建立状态检测影响下的设备剩余寿命的分布模型;
步骤104:根据所述状态检测影响下的退化模型,得到状态空间模型;
步骤105:采用最大期望算法对所述状态空间模型的参数进行估计,得到估计后的状态空间模型;
步骤106:根据所述估计后的状态空间模型对所述设备剩余寿命的分布模型进行自适应更新,得到更新后的设备剩余寿命的分布模型;
步骤107:根据所述更新后的设备剩余寿命的分布模型对设备的剩余寿命进行预测。
步骤101,具体包括:
采用线性wiener过程描述法建立设备自身的退化模型x′(t)=ηt+σbb(t)
其中,x′(t)为t时刻设备的退化水平;b(t)为标准布朗运动;η与σb分别为wiener过程的漂移系数与扩散系数。
步骤102,具体包括:
根据所述设备自身的退化模型采用幅值服从正态分布的随机冲击刻画单次状态检测,建立状态检测影响下的退化模型
其中,b(t)为标准布朗运动;η与σb分别为wiener过程的漂移系数与扩散系数,
步骤103,具体包括:
采用转换阈值方法,建立状态检测影响下的设备剩余寿命的分布模型,所述分布模型包括概率密度函数模型和累积分布函数模型;
所述概率密度函数模型通过下式表示:
所述累积分布函数模型通过下式表示:
其中,η为wiener过程的漂移系数,σb为wiener过程的扩散系数,μkl为时间区间[tk,tk+l)内状态检测引起的退化水平变化量之和的均值,
步骤104,具体包括:
根据所述状态检测影响下的退化模型,得到状态空间模型
其中,ηk为设备在第k次状态检测的瞬时退化速率;υk为第k次状态检测导致的设备性能突变,γk+1为第k+1次状态检测导致的设备性能突变,δxk为第k次状态检测与k+1次状态检测之间设备退化状态的增量;εk为随机效用,
步骤106,具体包括:
根据所述估计后的状态空间模型对所述设备剩余寿命的分布模型进行自适应更新,得到更新后的设备剩余寿命的分布模型,所述更新后的设备剩余寿命的分布模型包括更新后的概率密度函数模型和更新后的累积分布函数模型;
所述更新后的概率密度函数模型通过下式表示:
所述更新后的累积分布函数模型通过下式表示:
其中,
由于wiener过程即维纳过程具有良好的数学特性且在工程中应用广泛,本发明基于wiener过程构建考虑状态检测影响时设备的退化模型。为了更加直观理解考虑状态检测影响时设备的退化模型,首先研究不考虑状态检测影响的情形,然后对每次状态检测的影响进行单独建模,进而将状态检测的影响融入退化模型中。
基于wiener过程研究不考虑状态检测影响时的退化模型,基本模型为:
x′(t)=ηt+σbb(t)(1)
其中x′(t)为t时刻设备的退化水平;b(t)为标准布朗运动;η与σb分别为wiener过程的漂移系数与扩散系数;当设备退化过程呈现出非线性趋势时,可尝试通过特定的时间变换λ(t),能够将非线性退化过程转换为线性即:
x′(λ(t))=ηλ(t)+σbb(λ(t))(2)
将非线性退化过程线性化,转换后的退化过程仍旧能够利用式(1)对应的结论进行剩余寿命预测。因此,这里基于式(1)进行讨论。
为将状态检测的影响融入到设备的退化过程模型中,这里将状态检测的影响描述为导致设备退化水平突变的冲击。每次状态检测的时间通常远小于两次状态检测之间的间隔。因此,在建模过程中假设每次状态检测时间的可忽略不计。
需要说明的是,状态检测导致的退化水平突变量为随机变量。为简化模型,假设该随机变量服从均值与方差分别为μi和
基于上述描述,考虑状态检测的影响时,设备的退化过程可描述为:
其中
表示时间区间[0,t]内状态检测引起的退化水平变化量之和。根据正态分布的性质可知,ξ0(t)服从正态分布,即
首先研究不考虑状态检测影响的情形。根据wiener过程的基本特性,在首达时间的意义下,设备的寿命t服从逆高斯分布的随机变量,其概率密度函数为:
累积分布函数为:
其中,φ(·)为标准正态分布的累积分布函数,w为设备预先设定的失效阈值。
同理,设备在tk时刻的剩余寿命lk可定义为lk=inf{lk:x(tk+lk)≥wxk},其概率密度函数和累积分布函数分别为:
其中,wk为阈值与当前退化值的差值,wk=w-xk,xk为当前时刻的退化量。
引入状态检测的影响后,式(3)所定义的退化过程首达固定阈值w的时间可转化为式(1)所示随机过程首达时变随机阈值
根据全概率公式,在式(3)、(4)的基础上对转换后的失效阈值求期望,能够得到考虑状态检测影响下设备寿命t*的概率密度函数和累积分布函数,分别为:
同理,在预测设备的剩余寿命时,需要考虑未来一段时间内可能存在的状态检测对退化过程的影响。假定在时间区间[tk,tk+l)内,共有
为实现同时考虑时变不确定性、随机效用、以及状态检测等因素时设备剩余寿命的自适应在线预测,需要对退化模型进行改进。同时,根据工程经验,受到一系列冲击的作用后,设备的退化速率可能会加剧,而同样的冲击对不同退化速率设备的影响也截然不同,这里分别通过wiener过程的漂移系数及随机冲击描述设备的退化速率及状态检测对退化过程的影响,且退化模型改进为如下状态空间模型的形式:
其中ηk设备在第k次状态检测与k+1次状态检测之间的瞬时退化速率;δxk第k次状态检测与k+1次状态检测之间设备退化状态的增量;随机变量
其中,vk=σbb(δtk),yk=δxk,
这里,θk服从二元正态分布bvn(0,σ),其均值与方差分别为
令y1k=[y1,…,yk]t,接下来基于y1k采用卡尔曼滤波算法对漂移系数和状态检测影响进行在线更新。
步骤1:状态估计
pk|k-1=apk-1|k-1at+σk;
步骤2:协方差阵更新
其中,
其中,
需要说明的是,式(11)和(12)中的剩余寿命结果并未考虑了漂移系数的随机效用,而式(15)和式(16)考虑了漂移系数的期望和方差,实现了考虑周期性状态检测影响的设备剩余寿命在线预测。
基于以上分析,待估计参数包括a与σ,观测方程中的
假设截止tk时刻共对设备进行了k次状态检测,对应的k个观测值分别记为
漂移系数及检测影响的
e步:
其中
m步:最大化
为计算
其中,j≤k。
基于式(14)所构建的状态空间模型,上述待计算的条件期望值可通过卡尔曼方法和最优固定区间平滑方法计算。
步骤1:根据前述的卡尔曼滤波算法获取
步骤2:执行如下后向平滑步骤:
其中,sj为平滑的增益矩阵。
步骤3:初始化协方差阵
mk|k=(i-kkckt)ak-1pk-1|k-1;
步骤4:根据下式对协方差阵进行更新
其中,mj|k=cov(θj,θj-1|y1k)。
基于上一次迭代的参数估计值执行rts算法后,所需条件期望值可根据协方差阵的性质确定。具体地,有
其中,tr表示对矩阵求迹运算。求式(23)关于未知参数向量的偏导数,并令结果为零,解所得方程组,即能得到m步的参数估计值。最终的结果依次为:
其中,
具体实施例1:
陀螺仪是惯性导航系统的重要组成部件。惯性陀螺仪造价昂贵,其性能直接决定着导航系统的性能,且惯性器件的性能退化是导致惯性导航系统失效的主要原因之一。通常,在将陀螺仪安装在导弹等一次使用的武器装备上之前,会经历一段长期贮存期。在贮存环境与自身材料老化等因素的共同作用下,陀螺仪的性能会随着其贮存时间的增加而发生退化。为了保证贮存的陀螺仪安装后能够达到预期的性能指标,通常会周期性的对其进行状态检测,状态检测会改变陀螺仪的贮存环境与工作应力,从而加剧陀螺仪的退化过程。为提高陀螺仪寿命预测的精度,为陀螺仪的维修、替换及备件订购提供科学合理的依据,必须考虑状态检测对陀螺仪退化的影响。
与陀螺仪的贮存时间相比,每次状态检测的持续时间可以忽略不计。因此,可以采用所提方法预测陀螺仪的剩余寿命。这里所用的陀螺仪性能退化数据如图2所示。图2为本发明实施例陀螺仪的性能退化数据。当陀螺仪的一次项漂移超过阈值w=0.366(°/hour)时,将不再满足使用需求,故陀螺仪在第72个观测时刻失效。为说明与验证所提方法,下面用前67个数据估计退化模型的参数,并用剩余的数据对模型的参数及陀螺仪的剩余寿命进行在线更新。
选择相邻两次状态检测的差值作为状态空间模型的观测值,并利用所提的em算法估计未知参数。参数估计的结果分别如图3至图7所示。图3为本发明实施例参数a的估计过程。图4为本发明实施例参数σ的估计过程。图5为本发明实施例参数θ0|0的估计过程。图6为本发明实施例参数p0|0的估计过程。图7为本发明实施例参数
图8为本发明实施例滤波结果。图9为本发明实施例观测值与估计值的对比。利用估计得到的参数a、σ、
图10为本发明实施例不考虑检测影响且更新参数时剩余寿命的概率密度函数。图11为本发明实施例考虑检测影响且更新参数时剩余寿命的概率密度函数。考虑与不考虑状态检测影响时剩余寿命的概率密度函数及均值如图10和图11所示。图12为本发明实施例不考虑检测影响且未更新参数时剩余寿命的概率密度函数。图13为本发明实施例考虑检测影响但未更新参数时剩余寿命的概率密度函数。不进行在线更新时,所得的剩余寿命预测结果分别见图12和图13。从图中可以清楚地看到,考虑状态检测的影响能够提高剩余寿命预测的精度。
此外,为了对不同方法的剩余寿命预测结果进行量化比较,采用均方误差进行衡量,即:
其中
图14为本发明实施例剩余寿命预测的均方误差的对比。从图14中可以看出,考虑状态检测的影响能够降低剩余寿命预测的均方误差,而对参数及剩余寿命进行在线更新能够进一步降低估计的均方误差。因而,有必要在引入状态检测的影响后根据新的观测数据对模型参数进行更新,同时也验证了本发明所提考虑状态检测的随机系统退化建模与剩余寿命预测方法的有效性。
图15为本发明实施例考虑状态检测影响的设备剩余寿命预测系统结构图。如图15所示,一种考虑状态检测影响的设备剩余寿命预测系统,包括:
第一退化模型建立模块201,用于采用线性wiener过程描述法建立设备自身的退化模型;
第二退化模型建立模块202,用于根据所述设备自身的退化模型采用幅值服从正态分布的随机冲击刻画单次状态检测法,建立状态检测影响下的退化模型;
第一分布模型建立模块203,用于采用转换阈值方法,建立状态检测影响下的设备剩余寿命的分布模型;
第一状态空间模型确定模块204,用于根据所述状态检测影响下的退化模型,得到状态空间模型;
第二状态空间模型确定模块205,用于采用最大期望算法对所述状态空间模型的参数进行估计,得到估计后的状态空间模型;
第二分布模型建立模块206,用于根据所述估计后的状态空间模型对所述设备剩余寿命的分布模型进行自适应更新,得到更新后的设备剩余寿命的分布模型;
预测模块,用于根据所述更新后的设备剩余寿命的分布模型对设备的剩余寿命进行预测。
所述第一退化模型建立模块201,具体包括:
第一退化模型建立单元,用于采用线性wiener过程描述法建立设备自身的退化模型x′(t)=ηt+σbb(t),
其中,x′(t)为t时刻设备的退化水平;b(t)为标准布朗运动;η与σb分别为wiener过程的漂移系数与扩散系数。
所述第二退化模型建立模块202,具体包括:
第二退化模型建立单元,用于根据所述设备自身的退化模型采用幅值服从正态分布的随机冲击刻画单次状态检测,建立状态检测影响下的退化模型
其中,b(t)为标准布朗运动;η与σb分别为wiener过程的漂移系数与扩散系数,
所述第一分布模型建立模块203,具体包括:
第一分布模型建立单元,用于采用转换阈值方法,建立状态检测影响下的设备剩余寿命的分布模型,所述分布模型包括概率密度函数模型和累积分布函数模型;
所述概率密度函数模型通过下式表示:
所述累积分布函数模型通过下式表示:
其中,η为wiener过程的漂移系数,σb为wiener过程的扩散系数,μkl为时间区间[tk,tk+l)内状态检测引起的退化水平变化量之和的均值,
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。