一种基于长短期记忆网络的设备可靠性与性能评估方法转让专利
申请号 : CN202110261185.1
文献号 : CN112966443B
文献日 : 2022-04-29
发明人 : 徐廷学 , 王瑞奇 , 李海君 , 顾钧元
申请人 : 中国人民解放军海军航空大学
摘要 :
权利要求 :
1.一种基于长短期记忆网络的设备可靠性与性能评估方法,其特征在于:该方法包括以下步骤:
步骤S10:对n台样品设备进行全寿命周期监测,得到监测数据数据,并进行数据分批与初始化,以及标注前两批与最后两批数据标注性能退化量;其中样品设备全寿命周期监测数据,分别记作Ai,i=1,2,…,n;通过滑窗将每个数据Ai分为mi个监测批次,记作Bij,其中j=1,2,…,mi;取各设备前两批数据Bi1与Bi2,并标注其性能退化量c1与c2为区间[a1,a2]的随机数,取各设备的最后二批数据Bid,di=mi‑1,标注其性能退化量c3为[a3,a4]的随机数,选取设备寿命周期监测数据的倒数第一批数据 标注其性能退化量c4为[a5,a6]之间的随机数;参数选取为a1=0.01,a2=0.05,a3=0.89,a4=0.9,a5=0.9,a6=0.91;
步骤S20:根据所述的设备批次监测数据与性能退化量,建立长短期记忆LSTM神经网络,并输入数据进行网络训练;
步骤S30:根据所述训练好的LSTM网络,输入各设备全寿命监测数据,得到不同时刻设备性能退化指标;并求解各时刻性能退化指标的均值与方差,其计算如下:其中DHIij为将各设备全寿命监测数据Ai用训练好的LSTM模型进行预测,得到i设备j时刻的预测结果,该结果即为i设备j时刻性能退化指标DHIij,i=1,2,....n,j=1,2,....mi,性能退化指标代表了设备的性能指标; 画出所有样品设备全寿命周期中性能退化指标;通过统计n台试验数据的性能退化指标,获得各设备j时刻下的性能退化均值记作μ(j),获得2
各设备j时刻下的性能退化方差记作σ(j);
步骤S40:根据所述的各时刻性能退化指标的均值与方差数据,采用最小二乘方法拟合均值与方差随时间变化的斜率数据,并计算区间中性能退化指标的最大值,作为设备性能失效阈值如下:
DHIfault=max(DHIij),i=1,2...n,j=1,2,....mlinear;
2
其中 与 为通过最小二乘拟合均值μ和方差σ随时间变化的斜率,即 和mliner为均值和方差同时满足一次函数的最大时间值; max为求取最大值的函数,DHIfault为失效阈值;
步骤S50,根据所述的性能退化指标均值和方差的斜率数据以及设备性能失效阈值,计算设备失效概率随时间变化的分布密度函数,再通过积分运算计算故障失效阈值预测时间如下:
其中 和 为所述的性能退化指标均值和方差的斜率,DHIfault为设备性能失效阈值,f(t)为离散随机损伤累积模型的失效概率密度函数; MTBF(DHIfault)为设备从性能退化指标为0的健康状态退化到故障失效阈值DHIfault的预测时间;
步骤S60,根据所述的LSTM模型的性能退化指标计算每台设备退化到当前状态的预测时间;然后计算每台设备的当前设备状态到故障状态的平均时间如下:MTBFij=MTBF(DHIfault)‑MTBF(DHIij);
其中DHIij表示根据根据所述LSTM模型得到的i设备时刻j的性能退化指标值,MTBFij为i设备时刻j的设备状态到故障状态的平均故障时间;
步骤S70,根据所述的每台设备当前时刻到故障状态的平均预测时间,计算在不同剩余寿命假设情况下的性能指标与可靠性指标均值与协方差,然后计算不同寿命长度下的故障概率先验概率;
步骤S80,根据所述的训练好的LSTM网络输入监测到的正在运行的设备信息,得到性能退化评估值,然后计算其可靠性指标,再通过贝叶斯公式,根据可靠性指标与性能退化评估值,计算不同寿命假设下的后验概率如下:其中P(L)为通过统计获得样品设备在各寿命长度L下的故障概率先验概率,Lmax为设备的可能最长剩余寿命,可以取试验数据的最大寿命值即Lmax=max(mij), 和 为在剩余寿命RUL=l条件下,性能指标和可靠性指标的均值, 为在剩余寿命RUL=l条件下,性能指标和可靠性指标的协方差,最后P(MTBF,DHI|RUL=l)为所述的寿命长度L下的故障概率先验概率;
步骤S90,根据所述的后验概率与剩余寿命,联合求解最值的基于可靠性与性能一体化的剩余寿命预测结果如下:
MTBFnew=MTBF(DHIfault)‑MTBF(DHInew);
其中DHInew为一台正在运行的设备,当前监测信息输入训练好的LSTM网络,得到性能退化评估值, MTBFnew为将性能退化评估值DHInew代入可靠性评估模型得到的可靠性评估指标, P(RUL=l|MTBFnew,DHInew)为通过贝叶斯公式,计算在MTBFnew和DHInew条件下,各剩余寿命l=1,2,....RULmax的后验概率, 为根据设备后验概率分布P(RUL=l|MTBFnew,DHInew),求剩余寿命的期望值,也就是基于可靠性与性能一体化的剩余寿命预测结果。
说明书 :
一种基于长短期记忆网络的设备可靠性与性能评估方法
技术领域
背景技术
然后使用统计方法,选择最合适的统计分析模型,最后通过系统可靠性结构模型和部件寿
命分布模型,得到产品的可靠性。然而,一般情况下,随着产品的使用时间增加,表征产品性
能的特征参数会出现退化情况,性能退化现象是自然而大量存在的,并且由于不同设备使
用条件和运行历史的不同,运行同样时间后,其性能退化情况也存在不同,因此,会出现有
些设备寿命较短,有些设备寿命较长的现象。而传统的可靠性分析方法,仅考虑了“时间”这
一个变量,缺乏对设备运行条件、环境、当前状态等细节信息的考虑。这种可靠性评估方法,
造成对于长期运行后的设备可靠性评估存在较大误差,将状态仍然良好的设备,归为报废
设备,仅仅是因为使用次数临近或超过设计可靠性极限值;或者将性能已经出现明显退化
的设备,仍然评估为满足可靠性要求。当前,一些研究中借助长短期记忆网络(LSTM)具有捕
捉序列数据变化趋势特征的能力,将其用于设备的剩余寿命预测当中,对于存在明显性能
退化的设备,其预测结果较为准确。然而,对于设备退化早期,由于设备性能退化并不明显,
且存在着随机性的波动,因此,预测结果并不准确,误差较大。另一方面,对于受单位时间或
单位冲击损伤导致的设备性能退化过程,基于Pisson假设的离散随机损伤累积模型,可以
在缺乏或没有寿命数据的情况下预测设备的平均故障时间,即可以基于设备退化的早期数
据,预测设备的故障时间,但离散随机损伤模型和传统可靠性分析方法相同,缺乏对设备性
能退化情况的具体情况信息,因此在设备性能退化末期,模型误差较大,甚至存在超出模型
预测范围的情况。
发明内容
大的问题。
性能失效阈值;
测时间;
故障概率先验概率;
评估值,计算不同寿命假设下的后验概率;
括:
2,....mi,性能退化指标代表了设备的性能指标。画出所有样品设备全寿命周期中性能退
化指标,如图3所示。通过统计n台试验数据的性能退化指标,获得各时刻j下的性能退化均
2
值记作μ(j),方差记作σ(j)。
标的最大值包括:
DHIfault为失效阈值。
计算故障失效阈值预测时间包括:
指标为0的健康状态退化到故障失效阈值DHIfault的预测时间。
包括:
计算不同寿命长度下的故障概率先验概率包括:
在剩余寿命RUL=l条件下,性能指标和可靠性指标的协方差,最后P(MTBF,DHI|RUL=l)为
所求的寿命长度L下的故障概率先验概率。
可靠性指标与性能退化评估值,计算不同寿命假设下的后验概率,并联合求解最值的基于
可靠性与性能一体化的剩余寿命预测结果包括:
指标。P(RUL=l|MTBFnew,DHInew)为通过贝叶斯公式,计算在MTBFnew和DHInew条件下,各剩余
寿命l=1,2,....RULmax的后验概率。 为根据设备后验概率分布P(RUL=l|MTBFnew,
DHInew),求剩余寿命的期望值,也就是基于可靠性与性能一体化的剩余寿命预测结果。
叶斯信息融合技术,将可靠性与性能评估结果在剩余寿命这一指标上得以融合,形成可靠
性与性能一体化评估方法。该方法融合了离散随机损伤模型和LSTM对早期和晚期剩余寿命
预测的结果,从而可以更准确的评估各阶段设备的可靠性或剩余寿命。
附图说明
的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据
这些附图获得其他的附图。
具体实施方式
Bi1与Bi2,并标注其性能退化量c1与c2为区间[a1,a2]的随机数,取各设备的最后二批数据
Bid,di=mi‑1,标注其性能退化量c3为[a3,a4]的随机数,选取设备寿命周期监测数据的倒数
第一批数据 标注其性能退化量c4为[a5,a6]之间的随机数。其中a1=0.01,a2=0.05,a3
=0.89,a4=0.9,a5=0.9,a6=0.91。
2,....mi,性能退化指标代表了设备的性能指标。画出所有样品设备全寿命周期中性能退
化指标,如图3所示。通过统计n台试验数据的性能退化指标,获得各时刻j下的性能退化均
2
值μ(j)和方差σ(j),其中性能退化指标的均值和方差按照如下公式计算:
阈值按照如下公式计算:
变化的分布密度函数。通过如下公式计算设备从性能退化指标为0的健康状态退化到故障
失效阈值DHIfault的预测时间MTBF(DHIfault):
化指标值。
在剩余寿命RUL=l条件下,性能指标和可靠性指标的协方差,相应的计算公式如下:
MTBF(DHIfault)‑MTBF(DHInew),得到可靠性评估指标MTBFnew。最后,通过贝叶斯公式,计算在
MTBFnew和DHInew条件下,各剩余寿命l=1,2,....RULmax的后验概率P(RUL=l|MTBFnew,
DHInew):
在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保
护范围之内。