加速退化模型与可靠度评定结果准确性的定量验证方法技术

本发明专利技术涉及一种加速退化模型与可靠度评定结果准确性的定量验证方法，该方法主要包括以下步骤：(1)利用加速退化数据验证Wiener‑Arrhenius加速退化模型准确性；(2)利用常应力可靠性数据验证Wiener‑Arrhenius加速退化模型准确性；(3)利用常应力可靠性数据建立标准可靠度函数；(4)建立基于面积比的可靠度评估结果验证模型；(5)设计蒙特卡罗算法计算面积比。该方法为解决加速退化试验中的验证难题提供了一种有效手段。

Accelerated degradation model and quantitative verification method for accuracy of reliability evaluation results

The present invention relates to a quantitative verification method for the accuracy of accelerated degradation model and reliability evaluation results. The method mainly comprises the following steps: (1) verifying the accuracy of Wiener Arrhenius accelerated degradation model with accelerated degradation data; (2) verifying the accuracy of Wiener Arrhenius accelerated degradation model with constant stress reliability data. (3) Establish the standard reliability function by using the constant stress reliability data; (4) Establish the reliability evaluation result verification model based on area ratio; (5) Design the Monte Carlo algorithm to calculate the area ratio. This method provides an effective means to solve the problem of verification in accelerated degradation test.

【技术实现步骤摘要】
加速退化模型与可靠度评定结果准确性的定量验证方法
本专利技术属于可靠性工程领域，涉及一种加速退化模型与可靠度评定结果准确性的定量验证方法。
技术介绍
与常应力可靠性试验相比，加速退化试验无论是在试验操作还是统计模型方面都变得复杂，增加了可靠性评定结果的不确定性。加速应力可靠性试验实质上是牺牲部分评估精度换取试验效率，外推到常应力下的可靠度评估结果通常会与真实值存在一定的偏差，因此需要验证此偏差是否在可接受的范围内。利用加速退化数据推导产品在常应力下的可靠度函数，首先要建立产品的性能退化模型。假定某产品的性能参数具有退化特点并且其性能退化量Y(t)服从Wiener过程，则可将Y(t)表示为Y(t)＝μΛ(t)+σB(Λ(t))，(1)其中μ为漂移参数，σ(σ>0)为扩散参数，Λ(t)＝tΛ表示时间函数，Λ为形状参数，B(·)表示标准布朗运动。采用加速系数不变原则推导出μ,σ与环境应力相关Λ与环境应力无关，并且μ,σ值应该随应力水平呈比例变化。据此以绝对温度T为协变量，对μ,σ建立基于Arrhenius方程的加速模型为μ(T)＝exp(γ1-γ2/T)(2)σ(T)＝exp(γ3-0.5γ2/T)(3)式中，γ1,γ2,γ3为待估系数。建立Wiener-Arrhenius加速退化模型为Y(t；T)＝exp(γ1-γ2/T)tΛ+exp(γ3-0.5γ2/T)B(tΛ)(4)产品的可靠度函数为式中，Φ(·)表示标准正态分布函数。
技术实现思路
本专利技术的目的在于建立一套较为科学的加速退化模型与可靠度评定结果的验证方法，该方法的具体技术方案为：步骤一：利用加速退化数据验证Wiener-Arrhenius加速退化模型准确性利用全部加速退化数据验证Wiener-Arrhenius加速退化模型是否合理。步骤一1：估计出加速退化模型的参数向量θ＝(γ1,γ2,γ3,Λ)。根据Wiener过程的统计特性，独立增量△Y(t；T)服从如下形式的正态分布，△Y(t；T)～N(exp(γ1-γ2/T)△Λ(t),exp(2γ3-γ2/T)△Λ(t))，(6)式中，△Λ(t)＝(t+△t)Λ-tΛ。设yijk为在第k个温度水平Tk下第j个产品的第i次性能退化测量值，tijk为对应的测量时间，△yijk＝yijk-y(i-1)jk为性能退化增量，为时间增量，其中i＝1,2,…,Hk；j＝1,2,…,N；k＝1,2,…,M。为了估计出θ，根据式(6)建立如下对数似然函数极大化上式可获得模型参数估计值步骤一2：构建标准正态分布统计量。根据Wiener过程的统计特性，式中构建标准正态分布统计量zijk为步骤一3：检验zijk～N(0,1)是否成立。设原假设为，H0：zijk～N(0,1)成立；备选假设为，H1：zijk～N(0,1)不成立。采用Anderson-Darling法进行假设检验，设显著性水平为α＝0.05，当p≥α时不能拒绝原假设，认为建立的Wiener-Arrhenius加速退化模型是正确的，当p<α时拒绝原假设，认为建立的Wiener-Arrhenius加速退化模型是不正确的。步骤二：利用常应力可靠性数据验证Wiener-Arrhenius加速退化模型准确性随着时间的推移，能够累积一些产品在常应力下的可靠性数据，如性能退化数据与失效时间数据，将这些可靠性数据作为标准对Wiener-Arrhenius加速退化模型进行验证。设产品在常应力水平T0下的失效时间数据为ξj，j＝1,2,…,J；性能退化数据为xij,tij，i＝1,2,…,H；j＝1,2,…,N，具体验证步骤为步骤二1：建立基于失效时间数据的被检关系式。如果所建的Wiener-Arrhenius加速退化模型是准确的，应该服从如下形式的InverseGaussian分布，将上式作为被检关系式。步骤二2：采用Anderson-Darling法检验式(9)是否成立。步骤二3：建立基于性能退化数据的被检关系式。令如果所建的Wiener-Arrhenius加速退化模型是准确的，zij应该服从标准正态分布，据此建立被检关系式为：zij～N(0,1)(11)步骤二4：采用Anderson-Darling法验证式(11)是否成立。如果式(9)与式(11)同时成立，证明所建立的Wiener-Arrhenius加速退化模型是准确的，否则，证明所建立的Wiener-Arrhenius加速退化模型不准确。步骤三：利用常应力可靠性数据建立标准可靠度函数为了利用常应力T0下的失效时间数据ξj与性能退化数据xij,tij估计出可靠度函数的参数值，建立如下似然函数极大化似然函数L(μ,σ,Λ)，获得参数的极大似然估计值建立标准可靠度函数为步骤四：建立基于面积比的可靠度评估结果验证模型利用加速退化数据外推出的产品在T0下的可靠度函数记为根据此可靠度函数评估出产品的可靠度。为了验证可靠度评估结果的准确性，提出利用面积比s定量表征评估结果的相对误差，建立基于面积比的可靠度评估结果验证模型为其中，ε为判定阈值，如果s≤ε则认为可靠度评估结果是准确的，s(R(t；T0),R(*)(t))表示外推出的可靠度曲线R(t；T0)与标准可靠度曲线R(*)(t)围成的面积，s(R(*)(t))表示标准可靠度曲线R(*)(t)与坐标轴围成的面积，ε根据工程经验设为0.2。步骤五：设计蒙特卡罗算法计算面积比定积分难以通过传统的数学解析方法计算出，如何确定面积比s成为了一个难题。为了解决此难题，设计了蒙特卡罗算法计算面积比s，具体步骤为步骤五1：选取一个较大的横坐标值t*，使得R(t*；T0)＝0，R(*)(t*)＝0，由点(0,0)，(0,1)，(t*,1)，(t*,0)构成矩形区域；步骤五2：每次通过设t＝UNI(0,t*)，R＝UNI(0,1)在矩形区域内选一个随机点(t,R)；步骤五3：判断(t,R)是否落入R(*)(t)与横坐标之间或R(t；T0)与R(*)(t)之间；步骤五4：执行步骤五2至步骤五3共N次，统计(t,R)落入R(*)(t)与横坐标之间的次数K1，及落入R(t；T0)与R(*)(t)之间的次数K2；当N足够大的时候，如下关系式成立通过以上5个步骤，能够较为客观、科学得验证建立的加速退化模型与可靠度评定结果是否准确。附图及说明图1加速退化模型与可靠度评定结果准确性的定量验证方法的步骤流程示意图具体实施方式下面结合附图对本专利技术实现步骤进行进一步说明。实施例1：某惯导系统伺服电路在长期贮存过程中具有性能退化现象，表现为电路电压逐步降低，根据产品设计规范，伺服电路的电压测量值x相对于初始值x0的相对百分比变化量y＝100|x-x0|/x0到达失效阈值D＝10时，产品失效。温度的长期作用能够导致伺服电路的电路参数漂移、磁性减弱，是伺服电路电压值降低的主要环境敏感应力。伺服电路电压在常温T0＝298.16K下的退化过程较为缓慢，为了高效评估出此产品的可靠性，设计了温度应力加速退化试验用于加速产品的退化过程。随机抽取的22个试验样品被分配到3个加速温度应力T1＝323.16K，T2＝348.16K，T3＝368.16K下进行加速退化试验，表1中列出了各样品的测量时间和退化数据测本文档来自技高网...

【技术保护点】
1.一种加速退化模型与可靠度评定结果准确性的定量验证方法，其特征在于，步骤一：利用加速退化数据验证Wiener‑Arrhenius加速退化模型准确性；步骤二：利用常应力可靠性数据验证Wiener‑Arrhenius加速退化模型准确性，常应力可靠性数据分为常应力性能退化数据与常应力失效时间数据两类；步骤三：利用常应力可靠性数据建立标准可靠度函数，为

【技术特征摘要】
1.一种加速退化模型与可靠度评定结果准确性的定量验证方法，其特征在于，步骤一：利用加速退化数据验证Wiener-Arrhenius加速退化模型准确性；步骤二：利用常应力可靠性数据验证Wiener-Arrhenius加速退化模型准确性，常应力可靠性数据分为常应力性能退化数据与常应力失效时间数据两类；步骤三：利用常应力可靠性数据建立标准可靠度函数，为步骤四：建立基于面积比的可靠度评估结果验证模型，为步骤五：设计蒙特卡罗算法计算面积比。2.如权利要求1所述的一种加速退化模型与可靠度评定结果准确性的定量验证方法，其特征在于，步骤一通过以下3个分步骤验证Wiener-Arrhenius加速退化模型的准确性：步骤一1：估计出加速退化模型的参数向量θ＝(γ1,γ2,γ3,Λ)；根据Wiener过程的统计特性，独立增量△Y(t；T)服从如下形式的正态分布，△Y(t；T)～N(exp(γ1-γ2/T)△Λ(t),exp(2γ3-γ2/T)△Λ(t))，通过△Y(t；T)的概率密度函数建立如下对数似然函数极大化上式获得模型参数估计值步骤一2：构建标准正态分布统计量；由其中构建标准正态分布统计量zijk为步骤一3：检验zijk～N(0,1)是否成立；设原假设为，H0：zijk～N(0,1)成立；备选假设为，H1：zijk～N(0,1)不成立；采用Anderson-Darling法检验zijk～N(0,1)，计算出p值，设显著性水平为α＝0.05，当p≥α时不能拒绝原假设，验证结论为Wiener-Arrhenius加速退化模型准确；当p<α时拒绝原假设，验证结论为Wiener-Arrhenius加速退化模型不准确。3.如权利要求1所述的一种加速退化模型与可靠度评定结果准确性的定量验证方法，其特征在于...

【专利技术属性】
技术研发人员：王浩伟，刘萍，周源，
申请(专利权)人：刘萍，
类型：发明
国别省市：山东,37