Weibull分布综合应力恒加试验中参数的简化MLE方法与流程

本发明属于加速寿命试验技术领域，具体是涉及weibull分布综合应力恒加试验中参数的简化mle方法。

背景技术

加速寿命试验(acceleratedlifetest，alt)是在保持失效机理不变的条件下，通过加大试验应力来加速产品失效的一种寿命试验方法；其目的是运用加速模型外推加速应力下的试验数据，对产品在正常应力水平下的各种可靠性指标进行统计推断。alt可缩短试验周期、提高试验效率、降低试验成本，alt技术的引入为解决高可靠度长寿命产品的寿命与可靠性评估提供了新途径。alt分为恒定应力alt(简称恒加试验，csalt)，步进应力alt(简称步加试验，ssalt)和序进应力alt(简称序加试验)；alt的研究主要集中在加速模型、alt的优化设计和alt数据的统计分析三个方面。alt数据的统计分析方法主要有经典统计推断方法和贝叶斯(bayesian)统计推断方法，其中经典统计推断方法中mle方法应用最为广泛。目前，不少学者在alt数据的统计分析方面方面已作了大量研究工作。

nelson(1990)最先对alt进行了全面性研究，提出了weibull分布、正态分布和对数正态分布下完全样本和截尾样本情形alt数据的mle方法和最小二乘估计方法，且给出模型参数的置信区间；watkins(1994)研究了weibull分布下恒加试验的精简mle理论，并采用所提出的精简mle方法刻画了weibull分布的尺度参数与试验应力间的对数线性关系，算例分析结果表明所提方法是有效可行的；wang等(2000)采用mle方法研究了恒加试验weibull分布对数线性模型的参数估计问题，提出了恒加试验weibull分布对数线性模型参数的mle方法，真实试验数据分析结果表明所给方法是有效可行的；王炳兴等(2002)提出了定数截尾恒加试验weibull分布参数的渐进无偏估计及渐进置信区间估计的估计理论；watkins等(2008)提出了weibull分布下定数截尾恒加试验的精简mle理论，给出了与似然函数相关的fisher信息矩阵的精确计算方法，算例分析表明所提精简mle方法较传统mle提高了计算效率且降低了参数mle值的标准偏差；汤银才等(2009)提出了三参数weibull分布bayesian估计的laplace数值积分方法和gibbs抽样模拟方法，模拟算例说明了所提方法的可行性和有效性；xu等(2012)研究了相依指数分布竞争失效情形下恒加试验的mle方法，真实试验数据分析结果表明所提mle方法是可行的；陈文华等(2012)在weibull分布失效下，利用mle理论以产品中位寿命mle值渐进方差的均值和标准离差分别作为估计精度和稳定性的考核指标，提出了航天电连接器步加试验的模拟评价理论与方法，算例分析表明基于mle理论的模拟评价方法具有可行性；武东等(2013)研究了weibull分布下步加试验的bayesian估计方法，montecarlo仿真算例说明了所提bayesian估计方法的有效性；zheng等(2013)采用em算法和最小二乘法对两参数广义指数分布逐步定数混合截尾恒加试验数据进行了统计分析研究，获得了参数的渐进无偏估计和用于构造参数置信区间的fisher信息矩阵；张详坡(2014)研究了三参数weibull分布竞争失效场合序加试验的mle方法，实例分析结果表明mle方法是正确可行的且具有很好的估计效果。这些研究基本上都集中在单应力alt方面。

然而，在单应力alt中，为使产品在试验过程中失效机理保持不变，试验应力水平不宜过高，这样势必会造成试验周期较长；另一方面，运用传统的mle法对weibull分布的alt数据进行统计分析时，参数估计往往解析困难、计算冗繁，通常需要借助数值计算方法。因此，有必要在单应力alt数据统计分析方法的基础上，研究综合应力alt数据的统计分析方法。

技术实现要素：

本发明的目的，就是针对上述问题，以航天电连接器为研究对象，针对单应力alt的试验周期较长和weibull分布的最大似然估计难以求解问题，提出weibull分布下综合应力恒加试验的优化mle方法，并通过对航天电连接器综合应力恒加试验的montecarlo仿真失效时间进行统计分析，验证所提weibull分布下综合应力恒加试验优化mle方法的可行性和有效性。

本发明的技术方案为：

weibull分布综合应力恒加试验中参数的简化mle方法，该方法用于加速寿命试验中，weibull分布参数的极大似然估计；其特征在于，包括以下步骤：

a、采用热应力(温度)和非热应力(如电压、电流、湿度、振动等)同时作为加速应力进行试验，则建立产品寿命特征与试验应力间关系的广义艾林(eyring)模型为：

其中，η为特征寿命；t为热应力(绝对温度)，s为非热应力；a、b均为待定常数，e为激活能(ev)，k＝8.617×10^-5ev/℃为波耳兹曼常数；

并设定：

(1)在各应力水平组合(ti,sj)下，产品的寿命tij服从双参数weibull分布w(mij,ηij)，即tij～w(mij,ηij)，i＝0,1,2,…,k，j＝0,1,2,…,l，其累积分布函数为：

其中，mij>0为形状参数，ηij>0为特征寿命或尺度参数；

设定寿命tij～w(mij,ηij)，则对数服从极值分布g(μij,σij)，即yij～g(μij,σij)，其累积分布函数为：

其中，μij＝lnηij为极值分布g(μij,σij)的位置参数，σij＝1/mij为极值分布g(μij,σij)的尺度参数；

(2)在各应力水平组合(ti,sj)下,产品的失效机理保持不变，即所有应力水平组合下weibull分布的形状参数mij相同：这等价于极值分布g(μij,σij)的尺度参数σij相同：i＝0,1,2,…,k，j＝0,1,2,…,l；

(3)在各应力水平组合(ti,sj)下，weibull分布的特征寿命ηij与应力水平组合(ti,sj)间满足对数线性广义艾林加速模型：

其中，γ0＝lna，γ1＝e/10³k和γ2＝-b均为待估参数；φ(sj)＝lnsj均为转化应力水平，记

b、根据步骤a设定的条件，对产品进行定时截尾综合应力恒加试验：

设定试验总样本量为n，在应力水平组合(ti,sj)下，投入nij个样品进行定时截尾寿命试验，试验截止时间为τij，在[0,τij]内观测到rij个样品发生失效，失效时间依次为：

同时设定其余nij-rij个样品将在(τij,∞)内发生失效；

c、获取样本数据的对数似然函数：

根据产品的失效机理保持不变，和weibull分布的概率密度函数：

可得应力水平组合(ti,sj)下样本数据的似然函数为：

从而应力水平组合(ti,sj)下样本数据的对数似然函数为：

可得所有应力水平组合下全部样本数据的对数似然函数为：

其中：

d、获得全部样本数据的优化对数似然函数：

定义非负整数a、b和c，即a,b,c>0，由于

利用qa,b,c可求得对数似然函数lnl的一阶偏导数令：

可得：

将(9)式代入(8)可得全部样本数据的优化对数似然函数为：

lnl^*＝d0,0lnm+(m-1)qe-m(γ1d1,0+γ2d0,1)-d0,0lnq0,0,0(10)

e、获得weibull分布尺度参数的极大似然估计：

根据步骤d，获得公式(10)的一阶偏导数分别为：

lnl^*的二阶偏导数分别为：

lnl^*的二阶混合偏导数分别为：

通过令对数似然函数lnl^*的一阶偏导数为零，然后再解方程组，即可获得m、γ1和γ2的极大似然估计值和然后再将和代入式(9)即可得γ0的极大似然估计值

根据极大似然估计的不变性，即可得公式(1)的极大似然估计为：

在应力水平组合(t0,s0)下，weibull分布尺度参数的极大似然估计为：

其中，

可靠度为r的可靠寿命tr,00的极大似然估计为

可靠度r00(t)的极大似然估计值为

失效率λ00(t)的极大似然估计值为

平均寿命te的极大似然估计值为

上述方案中，根据mle理论，参数mle值和的协方差矩阵为

式(12)中：协方差矩阵∑用fisher信息矩阵f来估计，而fisher信息矩阵f是对数似然函数lnl的负二阶偏导数的数学期望，即

当对数似然函数二阶偏导数的数学期望的精确分布难以确定时，fisher信息阵f的元素可用对数似然函数lnl的负二阶偏导数值来近似，其计算公式为

可靠性统计模型参数的渐进置信区间：

(1)特征寿命η00的置信区间

根据mle的渐进正态性，统计量的分布可以采用正态分布近似，即

其中，

方差与协方差可通过fisher信息矩阵f的逆矩阵，即协方差矩阵∑得到。

由此得到μ00的置信度为1-α的双侧置信线为

从而由μ00＝lnη00得到η00的置信度为1-α的双侧置信线为

(2)形状参数m的置信区间

根据mle理论，统计量的分布亦可用正态分布近似，即

因此，在置信水平为1-α时，形状参数m的双侧置信线为

其中，的方差是协方差矩阵∑的元素，u1-α2是标准正态分布的1-α/2分位数。

(3)加速模型系数γi的置信区间

根据mle理论，统计量的分布亦可用正态分布近似，即

给定置信水平1-α，对数线性加速模型系数γi的双侧置信线为

与

其中，的方差由协方差矩阵∑给出。

本发明的方案中，与现有技术最直接的差异为，现有技术求解参数是通过令对数似然函数lnl的4个一阶偏导数为零，然后采用newton法解对数似然方程组求得参数的mle值和这种直接求解对数似然方程组的方法，计算效率十分低下。

本发明的有益效果为，通过对对数似然函数的优化，极大的提高了参数估计的效率，因此提高了加速寿命试验的效率，有利于产品的快速测试。

附图说明

图1为实施例在各应力水平组合下的weibull分布直线；

图2为实施例在正常应力水平组合下产品的可靠度曲线；

具体实施方式

下面结合附图和实施例进一步详细描述本发明的技术方案：

实施例

本例为仿真y11x系列航天电连接器的失效过程，以环境温度和振动为加速寿命试验的加速应力，并采用广义艾林(eyring)模型为加速模型。有关y11x系列航天电连接器加速试验的研究表明，在环境温度和振动应力综合作用下，y11x系列航天电连接器失效机理的改变点为(158℃，1.0g²/hz)，此即y11x系列航天电连接器的最高应力水平组合点。另一方面，根据国标gjb101a—1997，y11x系列航天电连接器的正常应力水平组合可选取为(85℃，0.06g²/hz)。

通过大量摸底试验数据的统计分析，以及参考现有参考文献在y11x系列航天电连接器加速试验与可靠性方面的研究，选取

作为广义艾林(eyring)—weibull可靠性统计模型参数的真值。

当加速应力水平数k＝l＝4时，采用均匀设计表u4(4²)来进行试验，当给定最高试验应力水平组合后，试验方案的设计变量为和φ2，即设计变量水平组合为若温度与振动应力的转化应力水平取为等间隔水平，即

则温度与振动应力均为4水平数时的均匀设计方案如表1所示：

表1综合应力恒加试验均匀设计方案

若各应力水平组合下试验的截尾时间均为50h，则均匀设计方案下航天电连接器的montecarlo仿真失效时间如表2所示：

表2均匀设计方案的montecarlo仿真失效时间

可靠性统计模型的检验

表2列举的仿真失效时间是否满足可靠性统计模型的基本假设，需要进行检验。首先，weibull分布的拟合优度检验和失效机理一致性的检验，可采用weibull概率纸法：

采用中位秩公式：

f(tijh)＝(h-0.3)/(nij+0.4)，h＝1,2,…,rij

估计产品的失效概率f(tijh)，并将各应力水平组合下的试验数据[tijh，f(tijh)]画在同一张weibull概率纸上，如图1所示。由图1容易发现，各应力水平组合下试验数据点的趋势均为直线，且数据点的趋势大致平行。因此可以认为在各温度与振动综合应力作用下，航天电连接器的montecarlo仿真失效时间均服从weibull分布，且形状参数m保持不变，即失效机理具有一致性。

weibull分布的拟合优度检验亦可采用范·蒙特福特检验法，其检验结果如表3所示：

表3weibull分布的拟合优度检验结果

由表3的结果可知，范·蒙特福特统计量f值均介于f分布的0.05分位点f0.05(f1,f2)与0.95分位点f0.95(f1,f2)之间，这表明表2列举的航天电连接器的仿真失效时间均服从weibull分布。

其次，失效机理一致性的检验可采用巴特利特检验法；在显著性水平α＝0.05下，巴特利特检验统计量的值b²/c＝2.5508，它小于卡方分布临界值χ²(0.95,3)＝7.8147。因此可以接受失效机理一致性的原假设，即可以认为在各温度与振动应力水平组合下weibull分布的形状参数相等。

模型参数的mle和渐进置信区间：

采用本发明的优化mle方法，并运用newton算法来估计可靠性统计模型的未知参数，γ1、γ2和m的初值为和newton方法迭代运算过程中，参数γ1、γ2和m，对数似然函数lnl^*，以及lnl^*的所有一、二阶偏导数值的收敛过程如表4所示：

表4newton方法迭代运算过程

由表4的结果可知，经过7步迭代运算后，lnl^*的一阶偏导数和均收敛到零，且参数γ1、γ2和m收敛到稳定值，由此可得γ1、γ2和m的mle值为

再将mle值和代入式(9)可得γ0的mle值为

从而可得加速模型的估计为

而采用传统mle方法，netwon算法经过39步迭代后，参数γ0、γ1、γ2和m才能收敛到稳定值和

由此可得本发明的方法相对于传统方法极大的提高了效率。

已知mle值和后，根据式(13)可获得与似然函数相关的fisher信息矩阵

再由式(12)可得与的协方差矩阵

估计出协方差矩阵∑后，根据第4部分介绍的渐进置信区间理论，可容易获得参数的渐进置信区间。在置信水平1-α＝90％下，参数γ0、γ1、γ2、m和正常应力水平组合下weibull分布特征寿命η0x的mle值，渐进置信区间如表5所示：

表5模型参数的mle和90％置信区间

正常应力水平组合下可靠性指标的mle：

利用加速模型(15)式可得正常应力水平组合下特征寿命η00的mle值为

可靠度r00(t)的mle值为

失效率λ00(t)的mle值为

可靠度为r的可靠寿命tr,00的mle为

平均寿命te的mle值为

te＝89.173h

可靠度估计的可靠度曲线如图2所示。