基于极大似然估计的定点测试综合应力加速寿命试验方案优化设计.docxVIP

基于极大似然估计的定点测试综合应力加速寿命试验方案优化设计 1 验证试验结果 对于高可靠、高成本的设备，在储存和使用过程中受到各种环境因素的影响。为了快速评估实际环境的可靠性，有必要加快综合负荷试验的进程。在加速寿命试验中,由于恒定应力加速寿命试验,具有试验方法简单,试验设备要求不高,试验理论较为成熟,试验容易取得成功,试验中得到的信息最多,试验结果较为准确的优点,因此被作为标准试验方法广泛采用。为克服其试验时间较长、所需试验样本量大的缺点,使产品在正常应力水平时寿命估计值的精度最高,已有众多学者对恒定应力加速寿命试验方案的优化设计进行了研究,其中包括综合应力加速寿命试验方案的优化设计,但都是基于产品性能参数能够连续监测的加速寿命试验。然而在实际试验过程中,产品性能参数往往无法实现在线的连续监测,只能在事先设定的若干时间点上进行检测,即定时测试,因而不能获得产品确切的失效时间,得到的只是某一时段内产品的失效数。为使所设计的试验方案更加贴近工程实际,本文以航天电连接器为研究对象,以产品中位寿命极大似然估计值的方差最小为目标,讨论了基于定时测试的综合应力加速寿命试验方案优化设计的理论与方法,并对Y11X系列航天电连接器定时测试综合应力加速寿命试验方案进行了优化设计。 2 产品寿命的统计模型1 航天电连接器在贮存和使用中对其可靠性影响最大的环境因素是环境温度和振动应力。在温度和振动应力的综合作用下,电连接器的接触寿命服从两参数的Weibull分布,其概率密度函数为: 式中:m为形状参数,η为特征寿命。 在不同的温度应力水平和振动应力水平的综合作用下,若电连接器的失效机理保持不变,则其具体的统计模型可表述为: 1)产品寿命t在统计上相互独立且服从两参数Weibull分布; 2)在不同应力水平下,Weibull分布的形状参数m保持不变; 3)产品特征寿命η与两应力之间的关系满足: 式中:ΔE为激活能(e V),k为玻尔兹曼常数(0.861 7×10–4e V/K),T为环境温度(℃),S为随机振动应力(g2/Hz))A和α为待定系数。 若令θ=lnt,则式(1)的Weibull分布就转化为极值分布,其概率密度函数为: 式中:μ=lnη,σ=1/m。则上述的统计模型可转化为更具一般性的线性-极值统计模型,具体可表述为: 1)产品寿命θ(即为lnt)在统计上独立且同服从极值分布; 2)在不同应力水平下,产品寿命的分布参数σ保持不变; 3)产品寿命的分布参数μ是应力x、y(一般为转化应力)的函数,即: 式中:β0=ln A,β1=ΔE/1 000k,β2=-α,x=1 000/T+273.15,y=ln S。 3 采用定时测试综合效应加快负荷快速通明试验方案优化设计标准 3.1 应力加速寿命试验 对于综合应力加速寿命试验,若有q个试验应力(即加速因子),每个试验应力均取K个应力水平,按照所有应力水平组合做全数寿命试验,需要进行Kq次寿命试验。为在估计精度基本保持不变的前提下减少试验次数,现采用均匀正交试验理论对应力水平进行组合,所需进行的试验次数即为相应的应力水平数K。对于2个试验应力、4个应力水平的综合应力加速寿命试验,仅需进行4次试验即可,其应力组合方式如表1所示,其中T1…T4,S1…S4。很显然,对于转化应力x和y,有: 3.2 加速寿命试验 为便于对基于定时测试的加速寿命试验方案进行优化,现假定: 1)在每次试验中,所有试样同时进行试验,除非试样快速失效,否则均在截尾时间τi(i=1,,K)终止。当然,为了对产品的可靠性做进一步的分析,若产品在截尾时间尚未失效,试验可继续进行。 2)在加速寿命试验中,最高应力水平TK和SK与正常应力水平T0和S0一般是给定的,且最高应力水平TK和SK的选择原则以不改变产品的失效机理为准则。 3)在定时截尾寿命试验中,由于可靠性特征量估计值的方差依赖于试验截尾时产品的失效数,而截尾时间长短又会影响产品失效的数量;其次,在总试验时间相同的条件下,由于各组试验截尾时间相等的试验方案其估计精度低于不等截尾时间的试验方案,现将截尾时间作为试验设计的一部分,在给定截尾时间上限τS的前提下,将各组试验的截尾时间τi(i=1,,K)作为优化变量。 4)在定时测试加速寿命试验中,各组试验的测试区间数相等,均为M;第i(i=1,,K)组试验的测试时间为ti1ti2ti M=τi,时间间隔按等概率原则确定。 5)为了保证足够的估计精度,每组试验截尾前的失效数均不小于某一期望的最小失效数rmin。 3.3 响应面试验方案 为了提高产品寿命估计值的置信程度,并有效缩短试验时间,现按照产品在正常应力水平T0和S0时其寿命分布的P阶分位数估计值的方差最小为准则,对综合应力加速寿命试验方案进行优化。一般取P=0.5。 产品在正常