专利名称：：步进应力加速退化试验优化设计方法
技术领域：
：本发明涉及一种步进应力加速退化试验的优化设计方法，属于加速退化试验
技术领域：
。
背景技术：
：随着现代科技的发展，高可靠、长寿命产品越来越多，为了评估这些产品的可靠性与寿命，美国罗姆航空中心在20世纪60年代就正式提出加速寿命试验(ALT)，目前ALT技术的研究已经取得了大量研究成果，但长寿命和高可靠性产品在ALT中只有少量失效出现或根本没有失效出现，这就对产品的寿命与可靠性评估带来了极大的困难。为解决这一问题，加速退化试验(ADT)应运而生。ADT克服了ALT只记录产品失效时间而不考虑产品性能退化情况的不足，通过对加速退化数据处理可以对高可靠长寿命产品的可靠性及寿命进行较好的评估。ADT需要在了解产品退化失效机理并建立加速模型之后实施，通过监控并记录产品的性能退化参数，了解其性能退化趋势，从而获得产品的退化数据，基于退化数据对产品进行可靠性与寿命的评估。因此，在不改变产品失效机理的前提下，如何对ADT方案进行科学合理的设计，以便快速有效的获得产品的性能退化数据，提高产品可靠性评估精度是当前加速退化试验研究的一个热点。恒定应力加速退化试验(CSADT)是指选择一组高于正常应力的加速应力水平，将一定数量的试件分为N组，每组在一个应力下实施加速退化试验；而步进应力加速退化试验(SSADT)是指选择一组高于正常应力的加速应力水平，在试验开始时将一定数量的样品都放置在试验箱内，所选的加速应力呈台阶式递增实施。CSADT的研究和应用相对成熟，但由于SSADT比CSADT具有节省样本、试验时间、监测仪器和试验箱的使用，大大减少试验费用等优势，因此目前SSADT受到更为广泛的关注，SSADT的相关研究也越来越多。SSADT方案主要包括试验样本量、试验应力、性能参数的监测记录频率、每一应力下的监测次数等内容。加速退化试验方案设计的研究文献大约从1994年开始，Boulanger&Escobar两位学者针对一类经典的退化模型对ADT进行了方案设计的研究，讨论了试验应力、总试验样本量、每一应力下样本量比例以及测试时间的确定问题。之后TsengST和YuHF两位学者在1997和1998年分别对退化试验和加速退化试验的截尾时间进行了探讨。2002年Shuo-JyeWu&Chun-TaoChang在费用约束条件下以p阶分位寿命的均方误差最小为目标，提出了退化率服从指数分布的退化试验优化设计方法。2004年LCTang，GYYang&MXie针对两应力水平的SSADT提出以试验费用最少为目标，以一定可靠度估计精度为约束，退化量服从随机过程的ADT优化设计方法，确定试验样本量和各应力水平的监测次数。同年，YuHF&TsengST以p阶分位寿命方差最小为目标，在试验费用约束条件下，对退化率服从倒数Weibull分布的退化试验进行了优化设计研究，得到了最优的样本量、监测频率和试验监测次数。2006年Chen-MaoLiao&TsengST以试验费用不超过总预算为目标，以产品p分位寿命估计的渐进方差最小为约束，用随机扩散过程对典型SSADT问题建模，来确定试验样本量，监测频率和试验截止时间。2009年YashunWang;ChunhuaZhang&XunChen提出基于MonteCarlo仿真的SSADT优化方法，以p阶分位寿命渐进方差的局部估计最小为目标，以试验费用为约束，利用MonteCarlo仿真方法模拟退化试验过程生成数据，通过统计分析计算P分位寿命渐进方差，选取渐进方差最小的方案为最优方案，方案给出试验样本量，监测频率和监测次数。Tseng,ST&Balakrishnan，N&Tsai，CC等学者对服从gamma过程的疲劳失效产品进行SSADT优化设计，以试验费用不超过预算为约束，以产品MTTF估计的渐进方差最小为目标，进行试验优化，最终确定试验样本量，监测频率和试验截止时间。XiaoyangLi&TongminJiang讨论了在多种退化失效模式下制订SSADT方案的方法，根据竞争失效机理，以漂移布朗运动建立多态退化系统的可靠性模型。以试验费用不超过预算为约束，以P分位寿命估计的渐进方差最小为目标，进行试验优化，最优方案给出试验样本，监测频率和监测次数。到目前，对SSADT方案优化设计的研究中，优化设计变量都没有包含试验应力水平，试验优化在试验应力水平给定的情况下进行，最优方案中只提供试验样本，监测频率，试验截尾时间或监测次数。但是试验应力水平是试验方案中很重要的一项，只有在合适的应力水平下才能更好地激发产品性能退化，获得更为有效的退化数据，从而提高产品可靠性评估精度。因此，在加速退化试验中将试验应力水平也作为一项优化变量，对试验样本量、试验应力水平、监测频率、各应力下监测次数进行系统的优化设计是目前加速退化试验优化设计领域的当务之急。
发明内容本发明的目的是为了解决上述问题，提出了将试验应力水平作为优化变量的步进应力加速退化试验(SSADT)优化方法。将试验方案中所涉及到的所有变量，包括试验样本量、试验应力水平、监测频率、各应力下监测次数都考虑在内，以P分位寿命可靠度的渐进方差最小为目标，以试验总费用不超过预算为约束，进行系统的优化设计，使SSADT优化设计方法更为全面完善，优化结果更具有工程实际意义和应用价值。本发明是步进应力加速退化试验优化设计方法，包括以下几个步骤步骤一、确定产品性能退化模型，加速模型和可靠度函数；；步骤二、建立目标函数；步骤三、确定约束条件和优化变量；步骤四、实现优化算法。本发明的优点在于(1)建立目标函数以p分位寿命可靠度的渐进方差最小为目标。计算p分位寿命可靠度的渐进方差时，采用差分方式对可靠度函数求偏导，该方法解决了复杂函数求导困难，容易出错等问题。(2)根据工程经验提出对每一应力下监测次数的约束条件，使试验优化结果更具有工程实际意义。由于低应力下产品性能退化慢，高应力下退化快，故低应力下的监测次数应大于高应力下的监测次数，而且根据工程实际情况，各应力下的监测次数应不小于10次，使每一应力下都能采集到足够的性能退化数据，即M工>M2>……>MK>IO，这一约束也大大减少了计算次数，提高了试验优化速率。(3)将试验应力作为优化变量，并根据实际情况给出应力水平的约束条件，对试验方案更系统更全面地进行优化。在实施SSADT前，应了解到产品的极限应力，即已知正常应力和工作极限应力Smax的情况下，试验应力约束条件可设为S。<S:<S2<……<SK《Smax。试验的优化模型中，目标函数minAsVar(R(lp))是一个关于样本量，监测间隔，监测次数，以及试验应力的函数，将试验应力作为优化变量，可将试验方案进行综合的全面的优化设计。图1是本发明步进应力加速退化试验优化设计方法的流程图2是本发明中步骤四的流程图。具体实施例方式下面将结合附图和实施例对本发明作进一步的详细说明。本发明是一种步进应力加速退化试验优化设计方法，流程如图1所示，包括以下几个步骤步骤一、确定产品性能退化模型、加速模型和可靠度函数；根据产品特点、敏感应力和性能参数退化趋势等，确定产品性能退化模型和产品的加速模型，进而确定退化模型的概率密度分布函数f(t，C)和可靠性函数R(t);产品的性能退化模型是指对产品性能参数随时间的退化趋势进行拟合的数学描述。目前，性能参数退化趋势拟合方法主要采用对性能参数随时间变化的确定函数中添加随机过程分量描述。性能退化模型主要有布朗(brown)运动、伽马(gamma)过程等，其中漂移布朗运动(带有漂移系数的布朗运动)的应用最多。加速模型用来描述产品性能参数退化速率与应力之间的关系，是一个仅与应力相关的确定性函数。常用的加速模型有阿伦尼斯(Arrhenius)模型、逆幂率模型、艾琳(Egring)模型等，其形式都可表示为ln(d(")-^+S^(s)，其中，伊(力是应力s的某一已知函数，d(s)为性能退化率，A、B为常数。若产品性能退化模型可采用对性能参数随时间变化的确定函数中添加随机过程分量来描述，加速模型通过变形可表示为ln(rf("^^+S，W，则可采用本发明的方法进行试验优化设计。例如当选取漂移布朗运动模型为产品性能退化模型，描述产品的性能退化过程，如式(1)所示:Y(t)=oB(t)+d(s)t+y0(1)其中y。_漂移布朗运动的起始点，即产品性能初始值；Y(t)-产品性能退化过程，是一个漂移布朗运动；B(t)-均值为O，方差为时间t的标准布朗运动，B(t)N(0，t);o-扩散系数，o>0，不随应力和时间而改变，是常数；d(S)_漂移系数，也称为性能退化率；本发明中选取Arrhenius模型为加速模型，加速应力为温度T，则d(s)为d(Tk)=exp(A+B/Tk)(2)A为常数，B=-Ea/k，k是波尔兹曼常数8.6171X10_5eV/K，Ea是激活能，单位为eV，Tk是绝对温度，单位为K。漂移布朗运动模型的概率密度分布函数为<formula>formulaseeoriginaldocumentpage12</formula>其中C为产品失效临界值；o为扩散系数；d(S)为应力S下的性能退化率；y。为性能初始值。具有形如上式的概率密度函数的分布，通常称为逆高斯分布，则产品的可靠度函数为卿=0<formula>formulaseeoriginaldocumentpage12</formula>①表示标准正态分布的累积概率分布函数，R(t)即为基于漂移布朗运动的可靠性评估模型。步骤二、建立目标函数；选取p阶分位寿命估计值的可靠度的渐进方差最小作为目标函数；艮卩minAsVar(R("))其中，Ip是产品在正常条件下的p阶分位寿命的估计值；具体为(1)建立产品性能退化模型的对数似然函数；设K个应力水平的步进应力加速退化试验，样本量为n，产品性能监测时间间隔为At，试验应力为Sk，各应力水平的性能监测次数为Mk，试验中共监测M次，M=Z:^^,k=1，2，…，K，每一应力水平的试验时间为MkAt，总试验时间为M.At;每次进行监控的时间为tkij，k=1，...，K;i=1，...，n，j=1，...，Mk，监控到的性能值为yikj，性能增量ADikj=yikj-yik(j—d。由采用性能增量表示概率密度分布函数f(At，AD)，得到极大似然函数L，进而得到对数似然函数InL;(2)求解产品在正常条件下的p阶分位寿命的估计值Ip;由步骤一确定的可靠性函数R(t)，对给定的分位点p，则^-i—乂p)是产品在正常条件下的P分位寿命的估计值；(3)求解渐进方差AsVar(R(Ip));Ip的估计值^在n—①时，服从均值为《p，方差为h11—1(9)h，记为AsVar(R(Ip))的渐进正态分布即，渐进方差为AsVar(R(《p))=h11—1(e)h;其中e为可靠度函数R(t)中的未知参数，假设存在，n个未知参数，则9=;<formula>formulaseeoriginaldocumentpage13</formula>(5)<formula>formulaseeoriginaldocumentpage13</formula>(6)狄I(e)为9的Fisher信息矩阵，是一个三阶正定方阵；①采用差分方法求h;采用差分方法对可靠度函数求偏导，得呵》-lim附,'+AH《)牟o(7)的lO—6或更小辨式中eig[epe2，，en，]，aei趋于o，取值为e②推导Fisher信息矩阵I(0);根据对数似然函数lnL，对参数9i求一阶偏导和二阶偏导，求解式(6)，得到信息矩阵I(e);③得到p阶分位寿命可靠度估计的渐进方差，建立目标函数；AsVar(R(Ip))是为关于样本量n，监测间隔At，监测次数[M!，M2，......，MK]和试验应力[Sp&，......，SJ的函数，最后将p阶分位寿命可靠度的渐进方差最小作为目标函数，即minAsVar(R(lp));例如，针对步骤一中的漂移布朗运动模型，步骤二具体为(1)建立产品性能退化模型的对数似然函数；将加速模型式(2)带入概率密度函数式(3)，用性能增量表示为1其极大似然函数为<formula>formulaseeoriginaldocumentpage13</formula>对数似然函数为<formula>formulaseeoriginaldocumentpage13</formula>(8)<formula>formulaseeoriginaldocumentpage13</formula>(9)<formula>formulaseeoriginaldocumentpage13</formula>(10)(2)求解产品在正常条件下的p阶分位寿命的估计值Ip;由式(4)，对给定的分位点p，则t=》—'(/7)是产品在正常条件下的p分位寿命的估1(3)求解渐进方差AsVar(R(Ip));^①时，服从均值为I计值。Ip的估计值《p在nAsVar(Rs(Ip))的渐进正态分布SP，渐进方差为AsVar(R("))其中朋(。朋(。a4,朋-^方差为hT"(e)h(记为9=i(e):(A，B，o2)(11)五五d4朋.朋2symmetrical五五(12)I(e)为9的Fisher信息矩阵，是一个三阶正定方阵。①采用差分方法求h;采用差分方法对可靠度函数求偏导，得朋(。二limi(^"+A^)—i(&，^a4m朋朋(。limA5(13)=lim由于可靠度函数R(t)的表达式复杂，故对其参数求偏导的解析式难以给出，推导过程极易出错，Chen-MaoLiao&TsengST对SSADT优化设计方法的研究中，附录里给出了对产品寿命分布函数求偏导的解析式，在对模型参数a求导时直接套用小(t，a)即令$=州,")，其他参数13，&，一也是如此，但①为标准正态分布，对1求导，="0，其他参数不能直接套用小(t)，故该文献中计算渐进方差有误。本发明中采用差分方式对可靠度函数求偏导，只要AA，AB，A一趋于0，取值为参数A，B，o2的10—6或更小，则该方法适用。②推导Fisher信息矩阵1(9);在At时间内的性能增量AD服从均值为d(s)At，方差为o2*At的正态分布。可知、广exp(乂+5/7;)-A/五=0(14)31n丄一朋—=1(15)利用公式(14)、公式(15)求解Fisher信息矩阵，公式推导如下「—厂exp(J+B/7;).M]S力—ZjZjZjM1=14=1产l々f^〖[AZ),&-expQ4+g/;).Af]eXp(j+g/7;)、,-l4=1产lCT々ln丄朋23ct22,=it=i户i1)五a2ln丄的求解为(=1A=l产l'exp(」+万/7;)—-2)五,5A,《二".l.^X.exp[204+5/7;)];的求解为32ln、^f》l"，p(j+B/。&4朋^五3)五AD,&—exp(J+B/7;)A/ln丄3厕.朋2Af+hexp(2^+25/7;)的求解为乂工ss=>五朋24)五五a2ln丄3扁的求解为=五「3MnZ)L朋&4JLa4朋J<formula>formulaseeoriginaldocumentpage16</formula>的求解为<formula>formulaseeoriginaldocumentpage16</formula>有由上式可见，除了方差自身的二阶偏导数，一切与方差相关的偏导数均等于零，故<formula>formulaseeoriginaldocumentpage16</formula>由以上可知，则Fisher信息矩阵为<formula>formulaseeoriginaldocumentpage16</formula>③得到p阶分位寿命可靠度估计的渐进方差，建立目标函数；<formula>formulaseeoriginaldocumentpage16</formula>将p阶分位寿命可靠度的渐进方差最小作为目标函数，即minAsVar(R(IP))。步骤三、确定约束条件和优化变量；①确定约束条件(1)试验总费用Ct;试验费用的主要因素包括两部分a)单位时间内的试验费用C。，包括试验设备折旧、产品性能测试仪器的使用、试验人员的工资以及资源(水、电等)消耗等。单位时间元/小时；b)试验样本单价Cd。因此试验总费用可以表示为<formula>formulaseeoriginaldocumentpage17</formula>(2)监测次数M！>M2>......>MK>10;根据工程经验每一应力下的性能检测次数应不小于10次，若某一应力下监测次数过少，则在评估时该应力下退化率的回归拟合精度难以确保；在失效机理不变的前提下，产品退化率随应力的升高而加快，为使每一应力下能够检测到足够的退化量来外推产品的寿命与可靠性，应使低应力水平下的监测次数应高于高应力水平。(3)应力水平S。<Si<S2<<SK《S隨；S工应高于产品使用时的正常应力S。，Smax为产品的工作应力极限，通过预试验或者产品的设计和使用说明等相关信息可以获得S^，进而根据实际情况确定最高应力S^各应力下产品的失效机理不变。②确定优化变量；(1)样本量n;(2)监测时间间隔At;(3)监测次数Mk，(k=1，2，…K);(4)应力水平Sk，(k=1，2，".K)。即步进应力加速退化试验共投入n个样本，对每一样本每隔At小时监测一次，进行K步应力水平试验，每应力水平Sk下，监测次数为Mk。步骤四、实现优化算法；由步骤一至步骤三，确定SSADT试验方案的优化问题，用数学描述为minAsVar(R("))<formula>formulaseeoriginaldocumentpage17</formula>上式为本发明中的优化模型。其中，需要优化的决策变量为n，At，Mk，Sk，它们的取值均为正整数。优化算法具体为首先根据工程经验确定退化模型式中的参数[^，92，…，9n，]、产品性能初值y。、产品退化临界值C、单位时间内的试验费用C。、试验样本单价Cd、试验总费用C。产品试验正常应力S。、最高应力Sp试验应力水平数K;然后根据步骤三所述的约束条件，确定所有满足约束条件的试验应力S2，，SK]和监测次数[M15M2，…，MJ，计算每个满足约束条件时相应的P阶分位寿命可靠度的渐进方差，最后选取方差最小的试验方案作为最优方案。优化算法的流程如图2所示，包括以下几个步骤St印l:根据历史数据、相似产品信息、工程经验、手册、标准等与试验样本相关的信息确定初始参数，包括产品性能退化模型的参数[epe2，…，en，];产品性能初值y。；产品退化临界值C;单位时间内的试验费用S。；试验样本单价Cd;试验总费用Ct;产品试验正常应力S。、最高应力Sp试验应力水平数K。St印2:计算p阶分位寿命Ip;根据试验最高应力SK和试验应力水平数K，以及应力水平的约束条件S。<S:<S2<……<SK《S^，求解[SnS2，，SK]的所有组合，记为解空间Qs。例如，当Sk表示温度吋，令从5(TC开始(正常应力S。=25°C)，并令各应力之间间隔为5的倍数，当Sk碗定后，S2，，SK—!的取值空间为{50，55，60，......，SK_5}，以Si<S2<<SK—!<SK为约束，可得到[SnS2，，SK—pSK]的解空间Qs。以K=3，最高应力S3=70为例，S"S2的取值空间为{50，55，60，65}，以S丄<S2<S3为约束，[S"S2，S3]的解空间为{[50，55，70]，[50，60，70]，[50，65，70]，[55，60，70]，[55，65，70]，[60，65，70]}共6个解。St印3:令样本量n从3开始(一般样本量不宜少于3)，由试验费用约束条件"《+ArfMrC。SC;,为求n的最大值，若令Mk取最小值10，At=1，则n的上限为nniax=(Ct-C0K10)/Cd。St印4:若n《nmax继续下一步，否则，跳到St印14。St印5:令监测间隔At=1，根据试验费用约束条件"'Q+A^f;M4'C。^C,，若令Mk取最小值，即Mk二10(k=1，…，K)则可得到At的最大值，其上限AL=^:f、:。St印6:若At《Atmax继续下一步，否则，跳到St印12。St印7:计算fM,[(C,-C,")/(C,&)],根据约束条件M!>M2>……A:=l>MK>10，则若f;A^2〖'10，继续下一步，否则跳到St印12。St印8:由监测次数的约束条件M！>M2>>MK>10和St印7计算得到的,求解[M!，M2，……，MK]的所有组合，记为解空间QM。at以K=3，Z^4^35为例，根据约束条件M！>M2>M3>10，[MnM2，M3]的解空间为{[15，10，10]，[14，11，10]，[13，12，10]，[13，11，11]，[12，12，11]}共5个解。St印9:由St印8求得的解空间Qm和St印2得到的Qs，两个解空间的所有组合QMXQs，按St印2和St印8中的例子，则监测次数[MnM2，......，MK]和试验应力[S丄，S2，…，Sk]的所有组合共6X5二30个。根据步骤二目标函数的建立过程，AsVar(R(lp))是一个关于样本量n，监测间隔At，监测次数[MpMy......，MK]和试验应力[SpS2，…，Sj的函数，将所有组合带入Fisher信息矩阵I(9)，计算渐进方差AsVar(R(IP))=h11—1(9)h，记为:g(AsVarR，n，At，M2，......，MK]，S2，，SK])(其中AsVarR为求得的渐进方差)；St印lO:取g中渐进方差最小，记为G(At，)=ming(AsVarR，n，At，[M丄，M2，......，Mj，[S"S2，…，SK]);Stepll:At=At+l，重复循环St印6St印lO;St印12:若At>At隨，取G(At，)中渐进方差最小，记为Gn(n，)=minG(At，)；St印13:令n=n+l，重复循环St印4St印12;St印14:若n>n隨，取Gn(n，:)中渐进方差最小，记为minGn=minGn(n，)，即为达到目标函数渐进方差最小的最优试验。实施例步骤一、确定产品性能退化模型，加速模型和可靠度函数；根据产品特点、敏感应力和性能参数退化趋势，选用漂移布朗运动模型描述产品性能退化过程，试验应力为温度T，对应的加速模型选用Arrhenius模型。则确定产品性能退化模型为式(1)。产品在正常应力下的可靠度函数为一步骤二、建立目标函数；建立目标函数minAsVar(R(lp))p阶分位寿命可靠度渐进方差AsVar(R(Ip))的计算过程见具体实施方式中的步骤二。步骤三、确定约束条件和优化变量；①确定约束条件(1)试验总费用Ct=100Xl(f元；试验费用的主要因素包括两部分a)单位时间内的试验费用C。，包括试验设备折旧、产品性能测试仪器的使用、试验人员的工资以及资源(水、电等)消耗等。单位时间元/小时；C。=0.6X103元/小时b)试验样本单价CdCd=3X103元/小时。试验总费用表示为<formula>formulaseeoriginaldocumentpage19</formula>(2)监测次数Mi>M2>......>MK>10。(3)应力水平T。<1\<T2<......<TK《T隨。②确定优化变量(1)样本量n;(2)监测时间间隔At;(3)监测次数Mk，(k=1，2，…K);(4)应力水平为温度Tk，(k=1，2，…K)。步骤四、实现优化算法。St印l:确定初始参数表1退化模型参数信息表模型参数A激活能Ea扩散系数o性能初始值y。性能临界值c估计值11.440.50.0710095B=_Ea/k=-5802.4;表2试验费用表<table>tableseeoriginaldocumentpage20</column></row><table>试验总费用Ct=100X103元；<table>tableseeoriginaldocumentpage20</column></row><table>选应力水平数K=4;St印2:计算p阶分位寿命《p:p=0.5，《。.5=R—1(0.5)=6349;寻找所有符合要求的[1\，T2，……，TK—"1k]的解空间Qp经计算中共165组解。St印3:令样本量n从3开始，其上限为nmax=(Ct_C。K10)/Cd=25;St印4到St印14按照优化算法的具体步骤实施，求出产品p分位寿命可靠度渐进方差最小时对应的最优试验。运算结果如下<table>tableseeoriginaldocumentpage21</column></row><table>权利要求步进应力加速退化试验优化设计方法，其特征在于，包括以下几个步骤步骤一、确定产品性能退化模型，加速模型和可靠度函数；确定产品性能退化模型和产品的加速模型，进而确定退化模型的概率密度分布函数f(t，C)和可靠性函数R(t)；步骤二、建立目标函数；选取p阶分位寿命估计值的可靠度的渐进方差最小作为目标函数；即minAsVar(R(ξp))(1)其中，ξp是产品在正常条件下的p阶分位寿命的估计值；具体为(1)建立产品性能退化模型的对数似然函数；设K个应力水平的步进应力加速退化试验，样本量为n，产品性能监测时间间隔为Δt，试验应力为Sk，各应力水平的性能监测次数为Mk，试验中共监测M次，k＝1，2，…，K，每一应力水平的试验时间为Mk·Δt，总试验时间为M·Δt，每次进行监控的时间为tkij，k＝1，...，K；i＝1，...，n，j＝1，...，Mk，监控到的性能值为yikj，性能增量ΔDikj＝yikj-yik(j-1)，由采用性能增量表示概率密度分布函数f(Δt，ΔD)，得到极大似然函数L，进而得到对数似然函数lnL；(2)求解产品在正常条件下的p阶分位寿命的估计值ξp；由步骤一确定的可靠性函数R(t)，对给定的分位点p，则是产品在正常条件下的p分位寿命的估计值；(3)求解渐进方差AsVar(R(ξp))；ξp的估计值在n→∞时，服从均值为ξp，方差为hTI-1(θ)h，记为AsVar(R(ξp))的渐进正态分布即，渐进方差为AsVar(R(ξp))＝hTI-1(θ)h；其中θ为可靠度函数R(t)中的未知参数，假设存在n个未知参数，则θ＝[θ1，θ2，…，θn，]；<mrow><msup><mi>h</mi><mi>T</mi></msup><mo>=</mo><mrow><mo>(</mo><mfrac><mrow><mo>∂</mo><mi>R</mi><mrow><mo>(</mo><msub><mi>ξ</mi><mi>p</mi></msub><mo>)</mo></mrow></mrow><mrow><mo>∂</mo><msub><mi>θ</mi><mn>1</mn></msub></mrow></mfrac><mo>,</mo><mfrac><mrow><mo>∂</mo><mi>R</mi><mrow><mo>(</mo><msub><mi>ξ</mi><mi>p</mi></msub><mo>)</mo></mrow></mrow><mrow><mo>∂</mo><msub><mi>θ</mi><mn>2</mn></msub></mrow></mfrac><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mfrac><mrow><mo>∂</mo><mi>R</mi><mrow><mo>(</mo><msub><mi>ξ</mi><mi>p</mi></msub><mo>)</mo></mrow></mrow><mrow><mo>∂</mo><msub><mi>θ</mi><mi>n</mi></msub></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>I(θ)为θ的Fisher信息矩阵，是一个三阶正定方阵；①采用差分方法求h；采用差分方法对可靠度函数求偏导，得<mrow><mfrac><mrow><mo>∂</mo><mi>R</mi><mrow><mo>(</mo><msub><mi>ξ</mi><mi>p</mi></msub><mo>)</mo></mrow></mrow><mrow><mo>∂</mo><msub><mi>θ</mi><mi>i</mi></msub></mrow></mfrac><mo>=</mo><munder><mi>lim</mi><mrow><mi>Δ</mi><msub><mi>θ</mi><mi>i</mi></msub><mo>→</mo><mn>0</mn></mrow></munder><mfrac><mrow><mi>R</mi><mrow><mo>(</mo><msub><mi>ξ</mi><mi>p</mi></msub><mo>,</mo><msub><mi>θ</mi><mi>i</mi></msub><mo>+</mo><msub><mi>Δθ</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><mi>R</mi><mrow><mo>(</mo><msub><mi>ξ</mi><mi>p</mi></msub><mo>,</mo><msub><mi>θ</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow><msub><mi>Δθ</mi><mi>i</mi></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>式中θi∈[θ1，θ2，…，θn，]，Δθi趋于0，取值为θi的10-6或更小；②推导Fisher信息矩阵I(θ)；根据对数似然函数lnL，对参数θi求一阶偏导和二阶偏导，求解式(3)，得到信息矩阵I(θ)；③得到p阶分位寿命可靠度估计的渐进方差，建立目标函数；AsVar(R(ξp))是为关于样本量n，监测间隔Δt，监测次数[M1，M2，......，MK]和试验应力[S1，S2，......，SK]的函数，最后将p阶分位寿命可靠度的渐进方差最小作为目标函数，即minAsVar(R(ξp))；步骤三、确定约束条件和优化变量；①确定约束条件(1)试验总费用Ct包括a)单位时间内的试验费用Co，单位时间元/小时；b)试验样本单价Cd；试验总费用Ct表示为<mrow><msub><mi>C</mi><mi>t</mi></msub><mo>=</mo><mi>n</mi><mo>·</mo><msub><mi>C</mi><mi>d</mi></msub><mo>+</mo><mi>Δt</mi><mo>·</mo><munderover><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>M</mi><mi>k</mi></msub><mo>·</mo><msub><mi>C</mi><mi>o</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>(2)监测次数M1≥M2≥……≥MK≥10；每一个应力下的性能检测次数大于等于10次；(3)应力水平S0＜S1＜S2＜……＜SK≤Smax；其中，S1至SK表示K个应力水平，S0为产品使用时的正常应力，Smax为产品的工作应力极限；②确定优化变量；(1)样本量n；(2)监测时间间隔Δt；(3)监测次数Mk，(k＝1，2，…K)；(4)应力水平Sk，(k＝1，2，…K)；即步进应力加速退化试验共投入n个样本，对每一样本每隔Δt小时监测一次，进行K步应力水平试验，每应力水平Sk下，监测次数为Mk；步骤四、实现优化算法；由步骤一至步骤三，确定步进应力加速退化试验方案的优化问题，用数学描述为minAsVar(R(ξp))<mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mi>n</mi><mo>·</mo><msub><mi>C</mi><mi>d</mi></msub><mo>+</mo><mi>Δt</mi><mo>·</mo><munderover><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>M</mi><mi>k</mi></msub><mo>·</mo><msub><mi>C</mi><mi>o</mi></msub><mo>≤</mo><msub><mi>C</mi><mi>t</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>M1≥M2≥……≥MK≥10S0＜S1＜S2＜……＜SK≤Smax上式为本发明中的优化模型，其中，需要优化的决策变量为n，Δt，Mk，Sk，它们的取值均为正整数；具体为首先根据工程经验确定退化模型式中的参数[θ1，θ2，…，θn]、产品性能初值y0、产品退化临界值C、单位时间内的试验费用Co、试验样本单价Cd、试验总费用Ct、产品试验正常应力S0、最高应力SK、试验应力水平数K；然后根据步骤三所述的约束条件，确定所有满足约束条件的试验应力[S1，S2，…，SK]和监测次数[M1，M2，…，MK]，计算每个满足约束条件时相应的p阶分位寿命可靠度的渐进方差，最后选取方差最小的试验方案作为最优方案。F2010100339987C00011.tif,F2010100339987C00012.tif,F2010100339987C00013.tif,F2010100339987C00021.tif2.根据权利要求1所述的步进应力加速退化试验优化设计方法，其特征在于，当选取漂移布朗运动模型为产品性能退化模型，所述的步骤一具体为产品的退化模型，如式(7)所示Y(t)=oB(t)+d(s)t+y0(7)其中y。_漂移布朗运动的起始点，即产品性能初始值；Y(t)-产品性能退化过程，是一个漂移布朗运动；B(t)-均值为0，方差为时间t的标准布朗运动，B(t)N(0，t);o-扩散系数，o>0，不随应力和时间而改变，是常数；d(s)-漂移系数，也称为性能退化率，它是一个仅与应力相关的确定性函数，是加速模假定d(s)为Arrhenius模型，加速应力为温度T，则d(s)为d(Tk)=exp(A+B/Tk)(8)A为常数，B=-Ea/k，k是波尔兹曼常数8.6171X10_5eV/K，Ea是激活能，单位为eV，Tk是绝对温度，单位为K;漂移布朗运动模型的概率密度分布函数为<formula>formulaseeoriginaldocumentpage4</formula>(9)其中C为产品失效临界值，o为扩散系数，d(s)为应力s下的性能退化率，y。为性能初始值，具有形如上式的概率密度函数的分布，称为逆高斯分布，则产品的可靠度函数为<formula>formulaseeoriginaldocumentpage4</formula>(10)①表示标准正态分布的累积概率分布函数，R(t)即为基于漂移布朗运动的可靠性评估模型。3.根据权利要求1所述的步进应力加速退化试验优化设计方法，其特征在于，当选取漂移布朗运动模型为产品性能退化模型，所述的步骤二为根据步骤一中的可靠性函数R(t)，选取p阶分位寿命估计值的可靠度的渐进方差最小作为目标函数即式(1);具体为(1)建立产品性能退化模型的对数似然函数；将加速模型式(8)带入概率密度函数式(9)，用性能增量表示为<formula>formulaseeoriginaldocumentpage4</formula>(11)其极大似然函数为<formula>formulaseeoriginaldocumentpage5</formula>对数似然函数为<formula>formulaseeoriginaldocumentpage5</formula>(13)(2)求解产品在正常条件下的p阶分位寿命的估计值Ip;由式(10)，对给定的分位点p，则^乂p)是产品在正常条件下的p分位寿命的估计值；(3)求解渐进方差AsVar(R(Ip));"的估计值^在n—①时，服从均值为《p，方差为hTl—ne)h，记为AsVar(Rs("))的渐进正态分布则渐进方差为<formula>formulaseeoriginaldocumentpage5</formula>(15)1(9)为9的Fisher信息矩阵，是一个三阶正定方阵①采用差分方法求h;采用差分方式对可靠度函数求偏导，得，)a4朋(。<formula>formulaseeoriginaldocumentpage5</formula>②推导Fisher信息矩阵I(e);在At时间内的性能增量AD服从均值为d(s)At，方差为o△t的正态分布，则<formula>formulaseeoriginaldocumentpage6</formula>利用公式(17)、公式(18)求解Fisher信息矩阵，具体过程中如下<formula>formulaseeoriginaldocumentpage6</formula>的求解为'[AZ^-exp04+B/7;).A,]<formula>formulaseeoriginaldocumentpage6</formula>2)五<formula>formulaseeoriginaldocumentpage6</formula>的求解为<formula>formulaseeoriginaldocumentpage6</formula>的求解为<formula>formulaseeoriginaldocumentpage6</formula>的求解为<formula>formulaseeoriginaldocumentpage7</formula>由上式，除了方差自身的二阶偏导数，一切与方差相关的偏导数均等于零，故有<formula>formulaseeoriginaldocumentpage7</formula>由以上，则Fisher信息矩阵为<formula>formulaseeoriginaldocumentpage7</formula>③得到P阶分位寿命可靠度估计的渐进方差，建立目标函数<formula>formulaseeoriginaldocumentpage7</formula>将P阶分位寿命可靠度的渐进方差最小作为目标函数，即minAsVar(R(Ip))。4.根据权利要求1所述的步进应力加速退化试验优化设计方法，其特征在于，所述的步骤四具体包括以下几个步骤St印l:根据工程经验确定初始参数，包括产品性能退化模型的参数92，…，en，]、产品性能初值y。、产品退化临界值c、单位时间内的试验费用c。、试验样本单价&、试验总费用Ct、产品试验正常应力S。、最高应力Sp试验应力水平数K;St印2:计算p阶分位寿命lp;根据试验最高应力SK和试验应力水平数K，以及应力水平的约束条件S。〈S工<S2<……<SK《S^，得到[SpS2，…，SJ的所有组合，记为解空间Qs;St印3:令样本量n从3开始，由试验费用约束条件"'(^ikr4-C。2C,,求n的最大值，令Mk=10，At=l，则n的上限为n隨=(C「C。K10)/Cd;St印4:若n《nmax继续下一步，否则，跳到St印14，其中nmax为总样本量；St印5:令监测间隔At=1，根据试验费用约束条件"《+A/'^MrC。SC;,若令Mk取最小值10，得到at的最大值，A"=;St印6:若At《Atmax继续下一步，否则，跳到St印12;St印7:计算f^=[(C<-C,")/(C。△,)]，根据约束条件M！>M2>……>MK>IO，若堂M^/M0，则继续下一步，否则跳到St印12;/:St印8:由监测次数的约束条件Mi>M2>……>MK>10和St印7计算得到的2]M4，求解[MpM^，MK]的所有组合，记为解空间QM;St印9:由St印8求得的解空间Qm和St印2得到的Qs，两个解空间的所有组合QMXQs，将所有组合带入Fisher信息矩阵I(e)，计算渐进方差AsVar(RUp))=h11—1(9)h，记为g(AsVarR，n，At，[M"M2，......，MK]，[S丄，S2，…，SK]);St印lO:取g中渐进方差最小，记为G(At，:)=ming(AsVarR，n，At，[M"M2，......，MJ，[S"S2，…，SK]);St印ll:At=At+l，重复循环St印6St印lO;St印12:若At>At隨，取G(At，)中渐进方差最小，记为Gn(n，)=minG(At，);St印13:令n=n+l，重复循环St印4St印12;St印14:若n>nmax，取Gn(n，:)中渐进方差最小，记为minGn=minGn(n，:)，则Gn(n，:)即为达到目标函数渐进方差最小的最优试验。全文摘要本发明公开了一种步进应力加速退化试验优化设计方法，包括以下几个步骤，步骤一、确定产品性能退化模型，加速模型和可靠度函数；步骤二、建立目标函数；步骤三、确定约束条件和优化变量；步骤四、实现优化算法；本发明以p分位寿命可靠度的渐进方差最小为目标。计算p分位寿命可靠度的渐进方差时，采用差分方式对可靠度函数求偏导，该方法解决了复杂函数求导困难，容易出错等问题；根据工程经验提出对每一应力下监测次数的约束条件，使试验优化结果更具有工程实际意义，M1≥M2≥……≥MK≥10的约束也大大减少了计算次数，提高了试验优化速率；将试验应力作为优化变量，并根据实际情况给出应力水平的约束条件，对试验方案更系统更全面地进行优化。文档编号G01R31/00GK101793927SQ201010033998公开日2010年8月4日申请日期2010年1月12日优先权日2010年1月12日发明者姜同敏,李晓阳,葛蒸蒸申请人:北京航空航天大学