法律状态公告日
法律状态信息
法律状态
2015-02-11
未缴年费专利权终止 IPC(主分类):G06F17/50 授权公告日:20130807 终止日期:20131219 申请日:20111219
专利权的终止
2013-08-07
授权
授权
2012-09-05
实质审查的生效 IPC(主分类):G06F17/50 申请日:20111219
实质审查的生效
2012-07-04
公开
公开
技术领域
本发明涉及一种快速判别温度应力加速实验失效机理一致性的方法,属于加速退化实验 设计技术领域。
技术背景
随着科学技术的发展,电子器件的可靠性越来越高,为了评价其可靠性参数,加速实验 的应用越来越广。然而,加速条件下器件潜在的失效机理可能会被激发出来,成为主要失效 机理,使得加速条件下失效机理发生改变,从而导致由加速实验获得的器件寿命不能代表器 件的真实寿命。目前,恒定应力加速实验失效一致性判别方法需要在实验结束后,再进行失 效机理一致性的判别,不能避免无效实验。因此,加速实验经常在未确定失效机理和未保证 与正常工作条件下具有相同失效机理的情况下进行。而本方法可在实验初期判别失效机理一 致性,避免无效加速实验,节约实验时间。
发明内容
本发明的目的在于:提供一种快速判别温度应力加速实验失效机理一致性的方法,使得 加速寿命实验中避免因退化过程中失效机理改变而导致错误的可靠性数据,保证了加速实验 设计方案的科学性和结果分析的有效性,避免无效加速实验,节约实验时间。
1.一种快速判别温度应力加速实验失效机理一致性的方法,其特征在于,包括以下步骤:
步骤一、确定L个加速实验温度应力水平T1<T2<...<Tl<...TL,其中l=1,2,...L;
步骤二、加速实验初期采集T1、T2、...Tl、...TL水平下的器件参数退化数据;
退化数据设定为失效敏感参数退化至任意程度下的退化时间,此处取失效敏感参数退化 量达3%时的退化时间为退化数据;
步骤三、根据步骤二中确定的参数退化数据拟合参数退化分布模型,判断其是否符合威 布尔退化分布模型,并判断应力水平Tl下的参数退化分布模型与应力T1,T2,...Tl-1下的参数 退化分布模型是否相同;
步骤四、若退化模型相同,对温度应力Tl下的退化数据进行威布尔退化分布模型参数估 计,具体步骤如下:
①确定威布尔模型参数估计的方法;
参数估计方法选择为极大似然估计;
威布尔分布模型有如下形式:
式中ti为温度应力Tl下失效敏感参数退化量,mi为形状参数,mi>0,i=0,1,2,…n,n为 温度应力Tl下的实验样品数,ηi为尺度参数,ηi>0;
威布尔模型参数极大似然估计-MLE方程组如下:
其中m为形状参数,η为尺度参数,r为产生故障数,ts定时截尾时间;
②形状参数m的参数估计方程表达式如下:
设实验用n个样品,共有n个产生故障,即r=n,可简化为:
形状参数m的估计方程为超越方程,用迭代法求解;
记
a.形状参数m0取
b.
c.判断|m2k+1-m2k|<ε是否成立,成立则转到下述步骤e,否则则转到下述步骤d;
d.m2k+2=(m2k+1+m2k)/2,k=k+1,并转到上述步骤b;
e.迭代结束,输出结果;
③尺度参数η的参数估计方程表达式如下:
设实验用n个样品,共有n个产生故障,即r=n,则
可以直接计算得到尺度参数η的参数估计;
步骤五、基于步骤四得到的迭代结果,判断是否服从失效机理不变的条件,具体步骤如 下:
①保持失效机理不变的两个条件的确定:
温度能够加快产品内部物理化学反应,促使产品参数加速退化,因此在加速寿命实验中 常用温度作为加速应力,此处,产品的参数退化速率与温度的关系常用Arrhenius模型表示:
式中M为温度敏感参数;t为实验时间;dM/dt为参数反应速率;A为常数;EA为激活能;kB为玻尔兹曼常数;T为实验温度;
加速条件下和正常工作条件下,器件的失效判据是一样的,即达到失效时器件的损伤累 积量是相等的,所不同的仅仅是加速系数,即参数的退化速率不同,因此:
ΔM=ΔM0 (13)
式中ΔM为加速条件下的敏感参数的退化量,ΔM0为正常工作条件下的敏感参数的退化量;
根据Arrhenius方程:
以上两式相比整理可得外推常温下的工作寿命t0:
令
t=t0/AF (17)
同一失效机理对应相同失效激活能,因此如果在加速应力条件下和正常工作下的激活能 EA不发生改变,则加速因子将严格保证与不同温度应力的倒数之差呈指数关系;
在加速实验的参数退化过程中,取失效敏感参数退化量为3%时的退化时间以服从威布尔 分布模型,设FU(t0)为正常工作水平下退化时间的分布函数,服从参数为m、η的威布尔分 布,其中f(t)为该威布尔分布的概率密度函数,则有:
可以得到:
FU(t0)=FS(t0/AF)=FS(t) (18)
施加应力后退化时间仍然服从威布尔分布,记为FS(t),但是其分布参数转变为m、AFη; 若该温度范围内失效机理未发生改变,参数估计与分布函数应满足上述条件;
根据获得的参数估计可作出如下判别:若威布尔分布的参数估计同时满足以下两个关系
1)形状参数m服从:m1=m2=...=ml;
2)尺度参数η服从:
则可以确定在T1-Tl范围内失效机理没有发生改变,加速实验是合理的;反之,失效机理 在加速实验过程中发生改变;
②判断第l个应力Tl下退化模型分布函数形状参数m是否服从关系1);
若退化模型分布函数形状参数m服从关系1),则继续判断尺度参数ηl是否服从关系2), 否则其失效机理在Tl下发生改变;
③判断第l个应力Tl下退化模型分布函数形状参数ηl是否服从关系2);
若退化模型分布函数尺度参数ηl服从关系2),则判断失效机理在Tl下未发生改变,否 则失效机理在Tl下发生改变,加速实验在Tl-Tl-1内失效机理保持一致;
步骤六、若退化模型不同,则认为在应力Tl下失效机理发生改变,不能继续进行加速退 化实验,加速实验在T1-Tl-1内失效机理保持一致。
附图说明
图1是本发明加速机理一致性判断流程图
图2是本发明迭代法求极大似然估计量m流程图
图3是本发明恒定温度应力下威布尔分布的概率密度曲线图
图4是本发明样本退化量分布示意图
图5是本发明威布尔分布下的Arrhenius寿命-温度关系图
具体实施方式
下面将结合附图1和示例对本发明做进一步的详细说明。
1.一种快速判别温度应力加速实验失效机理一致性的方法,其特征在于,包括以下步骤:
步骤一、确定L个加速实验温度应力水平T1<T2<...<Tl<...TL,其中l=1,2,...L;
步骤二、加速实验初期采集T1、T2、...Tl、...TL水平下的器件参数退化数据;
退化数据设定为失效敏感参数退化至任意程度下的退化时间,此处取失效敏感参数退化 量达3%时的退化时间为退化数据;
步骤三、根据步骤二中确定的参数退化数据拟合参数退化分布模型,判断其符合威布尔 退化分布模型(本发明以参数退化量服从威布尔分布模型为例),并判断应力水平Tl下的参数 退化分布模型与应力T1,T2,...Tl-1下的参数退化分布模型是否相同;
步骤四、若退化模型相同,对温度应力Tl下的退化数据进行威布尔退化分布模型参数估 计,具体步骤如下:
①确定威布尔模型参数估计的方法;
参数估计方法选择为极大似然估计;
威布尔分布模型有如下形式:
式中ti为温度应力Tl下失效敏感参数退化量,mi为形状参数,mi>0,i=0,1,2,…n,n为 温度应力Tl下的实验样品数,ηi为尺度参数,ηi>0;
威布尔模型参数极大似然估计-MLE方程组如下:
其中m为形状参数,η为尺度参数,r为产生故障数,ts定时截尾时间;
②形状参数m的参数估计方程表达式如下:
设实验用n个样品,共有n个产生故障,即r=n,可简化为:
形状参数m的估计方程为超越方程,用迭代法求解;
记
a.形状参数m0取
b.
c.判断|m2k+1-m2k|<ε是否成立,成立则转到下述步骤e,否则则转到下述步骤d;
d.m2k+2=(m2k+1+m2k)/2,k=k+1,并转到上述步骤b;
e.迭代结束,输出结果;
③尺度参数η的参数估计方程表达式如下:
设实验用n个样品,共有n个产生故障,即r=n,则
可以直接计算得到尺度参数η的参数估计;
步骤五、基于步骤四得到的迭代结果,判断是否服从失效机理不变的条件,具体步骤如 下:
①保持失效机理不变的两个条件的确定:
温度能够加快产品内部物理化学反应,促使产品参数加速退化,因此在加速寿命实验中 常用温度作为加速应力,此处,产品的参数退化速率与温度的关系常用Arrhenius模型表示:
式中M为温度敏感参数;t为实验时间;dM/dt为参数反应速率;A为常数;EA为激活能;kB 为玻尔兹曼常数;T为实验温度;
加速条件下和正常工作条件下,器件的失效判据是一样的,即达到失效时器件的损伤累 积量是相等的,所不同的仅仅是加速系数,即参数的退化速率不同,因此:
ΔM=ΔM0 (13)
式中ΔM为加速条件下的敏感参数的退化量,ΔM0为正常工作条件下的敏感参数的退化量;
根据Arrhenius方程:
以上两式相比整理可得外推常温下的工作寿命t0:
令
t=t0/AF (17)
同一失效机理对应相同失效激活能,因此如果在加速应力条件下和正常工作下的激活能 EA不发生改变,则加速因子将严格保证与不同温度应力的倒数之差呈指数关系;
在加速实验的参数退化过程中,取失效敏感参数退化量为3%时的退化时间以服从威布尔 分布模型,设FU(t0)为正常工作水平下退化时间的分布函数,服从参数为m、η的威布尔分 布,其中f(t)为该威布尔分布的概率密度函数,则有:
可以得到:
FU(t0)=FS(t0/AF)=FS(t) (18)
施加应力后退化时间仍然服从威布尔分布,记为FS(t),但是其分布参数转变为m、AFη; 若该温度范围内失效机理未发生改变,参数估计与分布函数应满足上述条件;
根据获得的参数估计可作出如下判别:若威布尔分布的参数估计同时满足以下两个关系
1)形状参数m服从:m1=m2=...=ml;
2)尺度参数η服从:
则可以确定在T1-Tl范围内失效机理没有发生改变,加速实验是合理的;反之,失效机理 在加速实验过程中发生改变;
②判断第l个应力Tl下退化模型分布函数形状参数m是否服从关系1);
若退化模型分布函数形状参数m服从关系1),则继续判断尺度参数ηl是否服从关系2), 否则其失效机理在Tl下发生改变;
③判断第l个应力Tl下退化模型分布函数形状参数ηl是否服从关系2);
若退化模型分布函数尺度参数ηl服从关系2),则判断失效机理在Tl下未发生改变,否 则失效机理在Tl下发生改变,加速实验在T1-Tl-1内失效机理保持一致;
步骤六、若退化模型不同,则认为在应力Tl下失效机理发生改变,不能继续进行加速退 化实验,加速实验在T1-Tl-1内失效机理保持一致。
示例:
下面通过赋予威布尔分布参数与激活能等参数获得参数退化量的理论数据,结合本发明 提出的失效一致性判别方法对各温度应力下的失效机理一致性判别进行详细说明。
步骤一、取4个不同应力390、400、410、420K进行恒定加速寿命实验,实验样本量每 组10个,共40只样品。赋予300K下的理论初值:激活能FA=1eV,形状参数m0=2,η0=1×107, 390、400、410K下的理论数据均由此外推得到。同时在420K下取激活能0.9eV作为理论失 效对比数据;
步骤二、计算机辅助取得在短期内各应力下器件失效敏感参数退化量达到3%时的退化时 间如表1、图4所示。
表1、失效敏感参数退化量达到3%时的退化时间
步骤三、根据步骤二中确定的参数退化数据进行拟合退化分布模型,确定该分布为威布尔 分布;
步骤四、若退化模型相同,进行分布模型参数估计,形状参数的具体估计步骤如图2的 迭代流程图所示;
对该理论退化数据的威布尔分布的形状参数m、尺度参数η进行参数估计,并计算加速 因子AF如表1所示,尺度参数可以由式(11)计算得到:
表2各温度应力下的形状参数m、尺度参数η和加速因子AF
利用参数估计在Arrhenius寿命-温度关系图上分别绘制390、400、410K下概率密度曲线 如图5所示。
步骤五、基于步骤五得到的迭代结果,判断是否服从失效机理不变的条件, 具体步骤如下:
①保持失效机理不变两个条件的确定:
威布尔分布的参数估计需要同时满足以下两个关系
1)形状参数m服从:m1=m2=...=ml;
2)尺度参数η服从:
则可以确定在T1-Tl范围内失效机理没有发生改变,加速实验是合理的;反之,失效机理在加 速实验过程中发生改变;
②判断应力390、400、410、420K下退化模型分布函数形状参数m是否服从关系1);
390、400、410K下的形状参数满足m390=m400=m410;420K下的形状参数不满足判别条 件1),即mi≠m420,i=390、400、410K;
③判断应力390、400、410、420K下退化模型分布函数形状参数η是否服从关系2);
390、400、410K下的尺度参数满足关系
420K下的形状参数不满足判别条件2),即
因此,可以判断失效机理在420K下发生改变,加速实验在390K-410K内失效机理保持不变。
机译: 分布式共享内存系统中缓存一致性的快速失效
机译: 快速判别装置,其加速方法及其程序
机译: 判别室性和室上性心动过速的方法和装置,以及判别窦性起源的快速心律和非窦性起源的快速心律的方法和装置