[发明专利]一种多应力多退化量步进加速退化试验方案优化设计方法

摘要

针对优化目标函数的解析形式难以推导的难题，本发明提供一种多应力多退化量步进加速退化试验方案优化设计方法，其基于蒙特卡罗统计仿真理论，建立多应力多退化量的步进加速退化试验方案优化设计优化模型，并提出相应的优化算法解决此优化问题，最终提出多应力多退化量步进加速退化试验方案优化设计方法。本发明所提出的方法易于流程化、便于工程应用，可为多应力多退化量场合产品寿命预测提供优化的试验方案支撑，以最小的试验代价实现最准确的寿命预测。

主权项

1.一种多应力多退化量步进加速退化试验方案优化设计方法，其特征在于：包括以下步骤：步骤1、获取产品加速退化试验相关信息1‑1)产品的退化量及失效阈值信息产品在工作或贮存过程中有m个退化量Yi(i＝1,2,...,m)随时间逐渐退化，一旦某个退化量Yi超过失效阈值Di(i＝1,2,...,m)，产品就会发生失效；1‑2)产品退化量的联合概率密度函数信息时刻t产品退化量Y＝(Y1,Y2,…,Ym)T服从多维正态分布，其联合概率密度函数表示为其中，μ＝(μ1,μ2,…,μm)T为均值向量，Σ为方差协方差矩阵σij(i＝1,...,m，j＝1,...,m)为退化量Yi与Yj的协方差；当i＝j时，σij为退化量Yi的方差；|Σ|为Σ的为行列式值；1‑3)产品退化模型及加速模型信息在不同应力水平组合α下，这里其中表示第i种加速应力的第li个水平，i＝1,...,s，li＝1,...,L，L为应力水平的个数，s为加速应力的个数；产品的m维正态分布均值向量μ的第j维元素与试验时间的关系满足如下退化模型μj＝bj+ajt j＝1,2,…,m     (3)式中，bj为截距参数，aj为退化速率参数，t为试验时间；退化速率参数aj与不同应力水平组合α之间满足如下多应力加速模型其中ηj0、ηji为加速模型的系数，Ti(·)为任意的单调函数，Si为第i种加速应力；式(4)两边取自然对数，有即式中，Aj(x)＝lnaj，γ′j0＝lnηj0，γ′ji＝ηji，xi＝ln[T(Si)]，i＝1,2,…,s；xi称为等效应力水平；将xi其进行标准化，得到取值范围为[0,1]的标准化等效应力水平ξi其中，xiL和xiH分别为应力xi的最低水平和最高水平；因此，式(6)可重写为其中，1‑4)产品性能退化的累积损伤模型信息令νi表示第i个应力水平组合αi下退化轨迹的起始时间，且此时的退化量与第i‑1个应力水平组合αi‑1结束时的退化量相等，则ν1是以下方程的解μj(ν1|α2)＝μj(t1|α1)    (9)类似地，νi满足μj(νi|αi+1)＝μj(νi‑1+ti‑ti‑1|αi)    (10)其中，i＝1,...,K‑1；因此，μj(t)能表示为其中，j＝1,...,m；因此，获取的产品加速退化试验模型参数先验信息描述为I＝(Σ,bj,γji,Dj),j＝1,…,m；i＝0,1,…,s    (12)步骤2、设计产品多应力多退化量步进加速退化试验基本方案；Yi的退化受到S1,S2,...,Ss种应力的影响，高于使用条件或贮存条件的s种应力组合能加速Yi退化过程；在进行多应力步进加速退化试验时，这s种加速应力的应力水平数均取为L；s种加速应力的最高应力水平设置应不使加速退化试验过程中产品的退化机理发生改变，即产品在这s种加速应力的加速退化试验中的退化机理与正常使用过程中的退化机理保持一致；令表示一种应力水平组合；根据均匀设计及正交设计原则，选取一系列应力水平组合α₁,α₂,…,α_K形成试验方案，其中K为应力水平组合的个数；如果试验方案是分式析因设计方案，则K＝L^s‑1；在开展多应力多退化量步进加速退化试验时，随机抽取N个样品在应力水平组合α1下进行试验，每隔F单位时间测试一次性能参数即监测频率为F，监测次数为M1；当试验进行到时间τ1时，应力水平组合由α1变为α2，继续进行试验，监测频率为F，监测次数为M2；当试验进行到时间τ2时，应力水平组合由α2变为α3，继续进行试验，监测频率为F，监测次数为M3；试验按如此方式进行，直到预定的时间结束；应力水平组合最终变为τK，监测频率为F，监测次数为MK，试验到时间τK时试验全部结束；步进加速退化试验每一应力水平组合下的试验时间为τi(i＝1,2,...,K)，且τi＝F·Mi·tu，其中tu为单位时间，为1天或1小时；因此，总试验时间τ可表示为于是步进加速退化试验的应力水平组合随时间的变化规律可表示为其中，步骤3、建立多应力多退化量步进加速退化试验方案优化模型；3‑1)确定优化模型的目标函数将产品在使用应力水平组合下的p阶分位寿命估计的均方误差平方根RMSE作为优化的目标函数：其中，E[·]表示数学期望；τ_p0为p阶分位寿命；为p阶分位寿命估计；及τ_p0的求解方法描述如下；p阶分位寿命τp0指产品在时刻τp0时的失效概率为p，此时产品的可靠度为1‑p，即R0(τp0)＝1‑p     (16)而产品在使用应力水平组合下时刻τ_p0的可靠度为通过联立方程(16)、(17)求解τ_p0；式(17)中，τ_p0为p阶分位寿命，R₀(τ_p0)为τ_p0时刻产品的可靠度，Y₀₁(t),...,Y_0m(t)为产品在使用应力水平组合下时刻t的m个退化量，D₁,...,D_m为产品m个退化量阈值，f(y₀₁,y₀₂,…,y_0m)为产品m个退化量的联合分布密度函数；f(y₀₁,y₀₂,…,y_0m|τ_p0)为τ_p0时刻产品m个退化量的联合分布密度函数；3‑2)确定优化模型的设计变量多应力多退化量步进加速退化试验的每一要素都能作为设计变量：①试验应力S1,S2,…,Ss；②试验应力水平及其组合α1,…,αK；③试验样品个数N；④监测频率F；⑤应力水平组合αj下的监测次数Mj；因此，试验方案可表示为d＝(s,Si,L,K,αj,N,F,Mj)，i＝1,…,sj＝1,…K；3‑3)确定优化模型的约束条件优化模型的约束条件如下：①试验总费用CT不超过试验预算Cb，CT≤Cb；②试验的应力数不低于2，s≥2；③每一应力的水平数不少于2，L≥2；④应力水平的组合数与应力水平数满足K＝Ls‑1；⑤试验样品个数不少于5，N≥5；⑥监测频率不小于1个时间单位，F≥1；⑦监测次数不少于3次，Mj≥3；⑧总的试验费用CT由试验运行费用、测量费用、试验样品费用组成，由下式计算：其中，Cop表示单位时间试验运行费用，CM表示单次测量的费用，Cd表示样品单价；综上所述，建立多应力多退化量步进加速退化试验方案优化模型可描述为步骤4、对优化模型进行优化；4‑1)根据约束条件构造可行试验方案集D，输入选取的试验方案数Z及蒙特卡罗仿真次数Nmc，令z＝1；4‑2)在D中选取一个方案dz＝(s,Si,L,K,αj,N,F,Mj)，i＝1,…,s，j＝1,…,K，z＝1,…,Z，令nmc＝1；4‑3)根据方案dz及前述获取的先验参数信息I＝(Σ,bj,γji,βj,Dj),j＝1,…,m；i＝0,1,…,s计算仿真参数：4‑3‑1)将α1,α2,…,αK中的标准化等效应力水平代入加速模型(8)计算aji(j＝1,...,m，i＝1,...,K)；4‑3‑2)根据式(3)、(9)、(10)、(11)可推导出其中j＝1,2,...,m；4‑3‑3)令t＝F,2F,…,MF，其中根据式(21)计算μ(t)＝[μ₁(t),μ₂(t),…,μ_m(t)]^T；4‑4)根据μ(t)及先验信息I中的Σ参数，对于t＝F,2F,…,MF生成N个m维正态分布N(μ(t),Σ)向量如下其中，Yn(tj)＝[Yn1(tj),Yn2(tj),…,Ynm(tj)]T且n＝1,…,N，j＝1,…,M，t1＝F,t2＝2F,…,tM＝MF；Yn(tj)从m维正态分布N(μ(tj),Σ)中抽样获得；4‑5)分析仿真数据(22)，按如下步骤计算4‑5‑1)对t＝F,2F,…,MF估计均值向量和方差协方差矩阵其中因此，方差协方差矩阵可用下式估计其中可令t＝kF，由式(24)计算；4‑5‑2)用模型(21)拟合4‑5‑1)得到的估计模型(21)的参数4‑5‑3)令即根据式(7)将它们转化为标准应力水平并代入式(8)计算使用应力水平组合下的4‑5‑4)将和代入式(3)确定使用应力水平组合下和t的关系，然后联立方程(16)、(17)求解令n_mc＝n_mc+1；4‑5‑5)如果n_mc≤N_mc则返回步骤4‑3)，重复步骤4‑3)～4‑5‑4)；否则，由以上步骤可得根据先验信息I中的参数b_j,γ_ji,i＝0,1,…,s,j＝1,…,m，再由步骤4‑5‑3)、4‑5‑4)求解τ_p0；由下式计算方案d_z对应的优化目标函数值U_z令z＝z+1；4‑6)如果z≤Z，则返回步骤4‑2)选取另一方案，重复步骤4‑3)～4‑5)，否则转到步骤4‑7)；4‑7)选取使U(d)最小的试验方案作为最优试验方案d*。