本发明专利技术涉及一种基于逆高斯过程的性能退化试验贝叶斯可靠性评估方法，包含：S1、建立逆高斯过程退化统计模型，确定可靠度函数和平均失效时间；S2、收集产品在研制各阶段中试验件的性能退化数据；S3、将试验件性能退化数据折合到正样产品在正常工作应力下的性能退化数据，作为贝叶斯可靠性估计的先验信息；S4、将正样产品的性能退化数据作为后验信息，结合先验信息，得到计算可靠度函数的参数；S5、计算可靠度函数和平均失效时间，得到产品可靠性的贝叶斯估计结果。本发明专利技术采用逆高斯模型刻画产品性能的退化过程，能更好的拟合数据；将研制各阶段的试验件性能退化数据纳入贝叶斯估计先验信息范畴，使产品贝叶斯可靠性评估结果更精确、更全面。
【技术实现步骤摘要】
基于逆高斯过程的性能退化试验贝叶斯可靠性评估方法
本专利技术涉及一种针对高可靠长寿命产品的性能退化试验贝叶斯可靠性评估技术，具体是指一种基于逆高斯过程的性能退化试验贝叶斯可靠性评估方法，属于性能退化试验分析、可靠性评估
技术介绍
随着技术的发展，产品的可靠性越来越高，寿命也越来越长，因此利用传统的寿命试验来评估产品的可靠性也已经越来越困难。然而，产品的功能是由其物理或化学性能参数动态表征的，失效与性能退化之间存在着必然的联系，性能退化信息中蕴含着丰富的与产品寿命有关的关键信息，有助于对产品的可靠性做出更准确的评估。对于性能退化试验的可靠性评估问题，目前有许多学者开展了研究。其中，李晓阳、孙富强等提出了基于模糊理论的加速退化试验预测方法，该方法针对每个应力水平下的性能退化数据进行回归分析，然后利用加速模型外推产品在正常工作应力水平下的性能退化率，并采用极大似然估计法估计漂移布朗运动的扩散系数，建立基于模糊理论的加速退化试验寿命与可靠性预测模型。王立志、李晓阳等提出了基于漂移布朗运动模型的加速退化试验贝叶斯估计方法，该方法通过确定加速模型参数的初值，建立各应力水平下模型参数的先验及后验分布，然后拟合加速模型参数的最终值，最后评估产品的寿命及可靠度。该方法结合相应的试验优化方法，可以对各应力水平分阶段进行优化试验设计。高军、黄道平等人提出了一种基于改进的漂移布朗运动的加速退化试验方法，该方法通过定义退化增量，扩展了常见的漂移布朗运动模型，对加速退化试验处理方法进行改进，用异点准确位置检验替代常用的统计分析检验，并执行异点剔除步骤，由此获得退化模型的参数估计值和失效概率。高军还提出了一种基于灰色系统理论的加速退化试验可靠性评估方法，该方法利用灰色系统理论对预处理后的数据进行分析，获得灰色系统理论分析数据进行可靠性评估，该方法能够实现对系统运行行为正确认识和有效控制。以上对产品性能退化过程的各种分析方法，主要采用了维纳过程、伽玛过程等随机过程，然而有时候维纳过程或者伽玛过程的拟合效果并不是非常好。Ye和Chen指出逆高斯过程是复合Poisson过程的极限，并进行了很好的物理解释，可以应用于退化试验中，并且比伽玛过程灵活。另一方面，对于复杂的高可靠性产品往往有多个具有高度继承性的研制阶段，研制各阶段试验件的性能退化信息在很大程度能够帮助评估最终产品的可靠性，即正样产品的可靠性，然而传统的性能退化试验分析往往忽略了这些宝贵的信息。因此，本专利技术提出一种基于逆高斯过程的性能退化试验贝叶斯可靠性评估方法，将研制各阶段的试验件的性能退化数据纳入贝叶斯估计先验信息范畴，结合正样产品性能退化试验数据，使产品可靠性评估结果更加精确、全面。
技术实现思路
本专利技术的目的是提供一种基于逆高斯过程的性能退化试验贝叶斯可靠性评估方法，采用逆高斯模型刻画产品性能的退化过程，能更好的拟合数据；将研制各阶段的试验件的性能退化数据纳入贝叶斯估计先验信息范畴，使产品贝叶斯可靠性评估结果更精确、更全面。为实现上述目的，本专利技术提供一种基于逆高斯过程的性能退化试验贝叶斯可靠性评估方法，包含以下步骤：S1、建立逆高斯过程退化统计模型，确定可靠度函数和平均失效时间的表达式；S2、收集产品在研制各阶段中的试验件的性能退化数据；S3、将获取的试验件性能退化数据折合转换到正样产品在正常工作应力下的性能退化数据，作为贝叶斯可靠性估计的先验信息；S4、将正样产品的性能退化数据作为贝叶斯可靠性估计的后验信息，结合S3中得到的先验信息，得到计算可靠度函数的参数；S5、计算可靠度函数和平均失效时间，得到产品可靠性的贝叶斯估计结果。所述的S1中，具体包含以下步骤：S11、将产品性能的退化过程看作一个逆高斯过程，退化量Y(t)，t>0具有如下性质：对于Y(0)＝0，以概率1成立；对于所有t>τ>u，有Y(t)-Y(τ)≥0以及Y(τ)-Y(u)≥0成立，且Y(t)-Y(τ)与Y(τ)-Y(u)相互独立；对于所有t>τ>u，有Y(t)-Y(τ)～IG(∧(t)-∧(τ),λ(∧(t)-∧(τ))2)；其中，IG表示逆高斯分布；∧(t)表示单调递增函数，且∧(0)＝0；λ表示超参数；符号～表示Y(t)-Y(τ)服从逆高斯分布IG；S12、假设产品性能退化的逆高斯过程为一个平稳过程，则∧(t)＝ηt，其中η为超参数，计算得到退化增量的期望和方差：E[Y(t+Δt)-Y(t)]＝∧(t+Δt)-∧(t)＝ηΔt；S13、假设ω为产品性能的失效阈值，产品性能的失效时间T为性能退化量Y(t)首次达到失效阈值ω的时间，即：T＝inf{t|Y(t)≥ω,t>0}；当t值接近无限大时，Y(t)则渐近服从均值为∧(t)，方差为的正态分布；S14、可靠度函数R(t)近似为：其中，P(Y(t)≥ω)表示Y(t)≥ω的概率，Φ(·)表示标准正态分布的分布函数；产品性能的平均失效时间MTTF为：所述的S2中，具体包含以下步骤：S21、假设正样产品的技术成熟度为1，通过专家打分法确定在研制各阶段中的试验件的技术成熟度m；具体为：记表示第i个专家zi对维度pk的技术等级打分，i＝1,2,...,ni，k＝1,2,...,nk，对维度pk的技术等级进行加权平均，得到试验件的技术成熟度m为：其中，qi表示专家zi的权重，qk表示维度pk的权重；S22、记ti,j,k为应力条件k下第i个产品第j次测量的时间，Yi,j,k为对应的产品性能测量值，保持相同的性能退化量Yi,j,k-Yi,j-1,k，对试验时间进行折合：其中，表示技术成熟度折合因子。所述的S3中，具体为：根据性能退化试验选用的应力类型，选择加速模型，对试验时间再一次进行折合：其中，Ak为加速因子，由所选择的加速模型确定；将非正常工作应力下的性能退化数据转换为正常工作应力下的性能退化数据，并作为贝叶斯可靠性估计的先验信息。所述的S4中，具体包含以下步骤：S41、将在研制各阶段中的试验件性能退化数据通过技术成熟度折合因子和加速因子折合后得到的性能退化数据作为先验信息，确定先验分布；S42、将正样产品的性能退化数据通过加速因子折合后得到的性能退化数据作为后验信息，确定后验分布；S43、设定逆高斯分布中的参数η服从伽玛分布，服从正态分布；利用贝叶斯可靠性估计法得到参数η、λ的后验分布所述的S5中，根据S14中得到的可靠度函数和平均失效时间的表达式，可得：综上所述，本专利技术所提供的基于逆高斯过程的性能退化试验贝叶斯可靠性评估方法，采用逆高斯模型刻画产品性能的退化过程，能更好的拟合数据；并且利用技术成熟度折合因子和加速因子将研制各阶段中的试验件的性能退化数据转化成正样产品在正常工作应力下的性能退化数据，拓展了先验信息的获取，使产品贝叶斯可靠性评估结果更精确、更全面。附图说明图1为本专利技术中的基于逆高斯过程的性能退化试验贝叶斯可靠性评估方法的流程图。具体实施方式以下结合图1，详细说明本专利技术的一个优选实施例。如图1所示，为本专利技术所提供的基于逆高斯过程的性能退化试验贝叶斯可靠性评估方法，包含以下步骤：S1、建立逆高斯过程退化统计模型，确定可靠度函数和平均失效时间的表达式；S2、收集产品在研制各阶段中的试验件的性能退化数据；S本文档来自技高网...
【技术保护点】
一种基于逆高斯过程的性能退化试验贝叶斯可靠性评估方法，其特征在于，包含以下步骤：S1、建立逆高斯过程退化统计模型，确定可靠度函数和平均失效时间的表达式；S2、收集产品在研制各阶段中的试验件的性能退化数据；S3、将获取的试验件性能退化数据折合转换到正样产品在正常工作应力下的性能退化数据，作为贝叶斯可靠性估计的先验信息；S4、将正样产品的性能退化数据作为贝叶斯可靠性估计的后验信息，结合S3中得到的先验信息，得到计算可靠度函数的参数；S5、计算可靠度函数和平均失效时间，得到产品可靠性的贝叶斯估计结果。
【技术特征摘要】
1.一种基于逆高斯过程的性能退化试验贝叶斯可靠性评估方法，其特征在于，包含以下步骤：S1、建立逆高斯过程退化统计模型，确定可靠度函数和平均失效时间的表达式；S2、收集产品在研制各阶段中的试验件的性能退化数据；S3、将获取的试验件性能退化数据折合转换到正样产品在正常工作应力下的性能退化数据，作为贝叶斯可靠性估计的先验信息；S4、将正样产品的性能退化数据作为贝叶斯可靠性估计的后验信息，结合S3中得到的先验信息，得到计算可靠度函数的参数；S5、计算可靠度函数和平均失效时间，得到产品可靠性的贝叶斯估计结果。2.如权利要求1所述的基于逆高斯过程的性能退化试验贝叶斯可靠性评估方法，其特征在于，所述的S1中，具体包含以下步骤：S11、将产品性能的退化过程看作一个逆高斯过程，退化量Y(t)，t>0具有如下性质：对于Y(0)＝0，以概率1成立；对于所有t>τ>u，有Y(t)-Y(τ)≥0以及Y(τ)-Y(u)≥0成立，且Y(t)-Y(τ)与Y(τ)-Y(u)相互独立；对于所有t>τ>u，有Y(t)-Y(τ)～IG(∧(t)-∧(τ),λ(∧(t)-∧(τ))2)；其中，IG表示逆高斯分布；∧(t)表示单调递增函数，且∧(0)＝0；λ表示超参数；符号～表示Y(t)-Y(τ)服从逆高斯分布IG；S12、假设产品性能退化的逆高斯过程为一个平稳过程，则∧(t)＝ηt，其中η为超参数，计算得到退化增量的期望和方差：E[Y(t+Δt)-Y(t)]＝∧(t+Δt)-∧(t)＝ηΔt；S13、假设ω为产品性能的失效阈值，产品性能的失效时间T为性能退化量Y(t)首次达到失效阈值ω的时间，即：T＝inf{t|Y(t)≥ω,t>0}；当t值接近无限大时，Y(t)则渐近服从均值为∧(t)，方差为的正态分布；S14、可靠度函数R(t)近似为：其中，P(Y(t)≥ω)表示Y(t)≥ω的概率，Φ(·)表示标准正态分布的分布函数；产品性能的平均失效时间MTTF为：3.如权利要求2所述的基于逆高斯过程的性能退化试验贝叶斯可靠性评估方法，其特征在于，所述的S2中，具体包含以下步骤：S21、假设正样产品的技术成熟度为1，通过专家打分法确定在研制各阶段中的试验件的技术成熟度m；具...
【专利技术属性】
技术研发人员：李鸿志，周郁，黄勇，张衡，李雁斌，张志俊，
申请(专利权)人：上海无线电设备研究所，
类型：发明
国别省市：上海,31
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。