一种弹上寿命件的多变量加速贮存试验优化设计方法
【专利摘要】本发明涉及一种弹上寿命件的多变量加速贮存试验优化设计方法，该设计方法包括如下步骤：1)确定弹上寿命件的性能退化模型和加速模型；2)确定试验应力施加方式和优化变量；3)建立综合优化目标函数；4)确定约束条件；5)采用遗传算法求解，得到最优方案。该设计方法能解决弹上寿命件加速贮存试验多个变量包括试验样本量、试验时间、试验应力、各应力水平下样本分配、时间分配和检测间隔等需要综合优化设计的问题，同时建立了综合优化目标函数，使优化结果更具有工程适用性。
【专利说明】一种弹上寿命件的多变量加速贮存试验优化设计方法
【技术领域】
[0001]本发明涉及一种加速贮存试验优化设计方法，特别涉及一种弹上寿命件的多变量加速贮存试验优化设计方法，属于可靠性【技术领域】。
【背景技术】
[0002]导弹作为“长期贮存，一次使用”的产品，其贮存寿命(贮存可靠性、可靠贮存寿命)是重要的战术技术指标。目前贮存试验主要有自然(现场)贮存试验、实验室模拟贮存试验和加速贮存试验。根据导弹贮存指标的要求，可知弹上寿命件多数属于高可靠、长寿命产品，为了在有限时间内给出这些弹上寿命件的贮存寿命与可靠性评估结果，加速贮存试验成为解决这一问题的首选试验技术。加速贮存试验在型号应用中首先必然面临着试验方案的设计问题，即如何在有限的时间和费用下科学合理的安排试验应力水平、受试样本、试验时间、检测间隔等试验变量，以获得最有效的性能退化信息，使产品寿命与可靠性评估最准确。采用优化的试验方案，是得到产品寿命与可靠性准确评估结果的基础和前提，不仅为产品研制和使用方提供正确的决策依据；还可以大大提高试验效率，使试验资源得到充分利用，降低产品的研制成本。
[0003]目前国内外加速试验优化设计的研究，针对不同退化模型、不同加速模型、不同应力施加方式、不同优化目标或约束条件，已有较多研究成果。但是多数研究以各加速应力水平下的受试样本量和试验时间为重点优化变量，而不考虑应力水平或检测间隔。专利201010033998.7对应力水平进行了优化，但各应力水平下的检测间隔相等，没有进行优化设置。而有些研究对各应力水平下检测间隔进行了优化，但没有优化应力水平。
[0004]对弹上寿命件的加速贮存试验而言，应力水平和检测间隔同样是非常重要的试验变量。只有在合适的应力水平下才能更好地激发产品性能退化，获得更为有效的退化数据。而检测间隔不仅关系到试验样本量，而且影响到试验的有效性。比如，对于检测为破坏性测量的产品，如密封件、减震垫、胶粘剂等高分子材料，检测间隔小则需要大量样本，增加试验负担；检测间隔大则获得的试验信息少。对于需要通电测试的电子、机电、光电类产品或者需要恢复到正常应力下进行测试的产品，也即测试时引入了其他应力的情况下，检测间隔小则实际试验不符合贮存试验条件，会影响评估结果；检测间隔大，得到的性能退化数据少，同样会影响评估的准确性。因此，为了更好地激发产品性能退化、减少测试带来的误差，应力水平和各应力水平下的检测间隔也应作为重点变量进行优化。
[0005]针对弹上寿命件的加速贮存试验特点，提出一种有效的适用于工程应用的多变量优化方法是当前导弹贮存寿命评估与延寿领域的当务之急。
【发明内容】
[0006]本发明针对现有加速贮存试验优化设计方法研究中不对应力水平或检测间隔进行优化的问题，提供一种弹上寿命件的多变量加速贮存试验优化设计方法，以解决弹上寿命件加速贮存试验多个变量包括试验样本量、试验时间、试验应力、各应力水平下样本分配、时间分配和检测间隔等需要综合优化设计的问题。针对现有加速贮存试验优化设计方法仅依据一种优化目标时，存在试验变量不符合工程应用的问题，建立了综合优化目标函数，使优化结果更具有工程适用性。
[0007]本发明的目的通过以下技术方案来实现:
[0008]一种弹上寿命件的多变量加速贮存试验优化设计方法的具体步骤为:
[0009]第一步确定弹上寿命件的性能退化模型和加速模型
[0010]所述退化模型包括漂移布朗运动、混合效应模型、伽马(ga_a)过程等，此处确定漂移布朗运动模型为产品性能退化模型，描述产品的性能退化过程，如式(I)所示:
[0011]Y (t) = σ B (t)+d ⑶.t+yQ(I)
[0012]式中:y(l为漂移布朗运动的起始点，即产品性能初始值；Y(t)为产品性能退化过程，是一个漂移布朗运动；B(t)是均值为O,方差为时间t的标准布朗运动，B(t)~N(0，t)；σ是扩散系数，σ>0，不随应力和时间而改变，是常数；d(S)是漂移系数，也称为性能退化率；
[0013]所述加速模 型有阿伦尼斯(Arrhenius)模型、逆幂率模型、艾琳(Egring)模型等，
其形式都可表示为:
[0014]
ln(t/{S)) = Α + Β.φ(Ξ)(21
[0015]式中，(HS)为性能退化率，A、B为常数，9(s)是应力S的已知函数。例如，对阿伦尼斯模型，9>(S) = IfSi S是绝对温度，A为常数，B=-Ea/k, k是波尔兹曼常数，Ea是激活能。
[0016]以引起产品性能退化的主要应力为试验应力，根据试验应力确定相应的加速模型。
[0017]第二步确定试验应力施加方式和优化变量
[0018]工程上常用的试验应力施加方式为恒定应力和步进应力两种，随着加速试验技术的发展，步进应力施加方式以其节约试验资源的优点越来越受到青睐，步进应力加速试验剖面如图2所示，步进应力加速试验的变量有:
[0019](I)总样本量η;
[0020](2)各应力水平Sk(k=l，...，K);工程上一般选取3~5个加速应力水平。
[0021](3)各应力水平下的试验时间tk (总试验时间为,=乙二4 )；
[0022](4)各应力水平下的检测间隔Λ tk (检测次数为Mk=tk/ Δ tk)；
[0023]上述优化变量中，样本量η和各应力水平下试验时间tk是以往研究中的重点优化变量，有些研究不考虑应力水平优化，有些研究则不考虑检测间隔优化，本发明将试验剖面中的四类变量进行综合优化。
[0024]第三步建立综合优化目标函数
[0025]针对试验设计人员关注点不同，优化目标可分为两类:(1)以“预测”为目标:关注产品正常应力下可靠性或寿命相关指标的预测精度，一般以正常应力下该指标参数的渐近方差(或均方误差等)最小为目标；(2)以“评估”为目标:关注模型未知参数的评估精度，以模型参数估计误差最小为目标，主要包括D-优化(最小化参数方差-协方差矩阵行列式的值，也即最大化信息矩阵行列式的值)和A-优化(最小化方差-协方差矩阵的迹)。
[0026]根据上述两类优化目标，结合导弹弹上寿命件加速贮存试验目的，制定两个优化目标:
[0027](I)以弹上寿命件在正常贮存应力下贮存可靠度的渐进方差最小为目标；
[0028]min AsVar (R( ξ , S0, θ ))(3)
[0029]式中，ξ是寿命件在正常贮存条件下贮存寿命的估计值，Stl为正常应力水平，向量Q=Le1，…，θ5]τ为退化模型和加速模型的参数，θ可采用极大似然估计得到。R( ξ , S0, θ )是在正常应力下的忙存可靠度，AsVar是渐近方差Asymptotic Variance的缩与。
[0030]AsVar (R(l,S0, Θ )) =hTF_1 ( θ ) h(4)
[0031]式中，h =潘(dQ)为可靠度函数在参数fl=暴处的偏导，Fii(B) = -E ° 气(:’6)
* δθ-,— V νθ'?θ> J
为对数似然函数lnL(x; Θ)的各二阶偏导在0=0处的期望的负值，X为退化性能增量。
[0032](2)以模型参数方差-协方差矩阵(Σ ( Θ ))行列式的值最小为目标。
[0033]min det Σ ( Θ )(5)
[0034]由于信息矩阵(F·(0))与方差-协方差矩阵(Σ (Θ))互逆，因此也即最大化信息矩阵的行列式值。即
[0035]max detF ( Θ )(6)
[0036]根据试验所知，只考虑式(4)所示目标函数，优化结果会出现高应力下试验资源过少的情况；只考虑式(5)所示目标函数，优化结果则会出现各应力下试验资源均分的情况。因此本发明提出同时考虑两种优化目标，采用权重系数法，将两种优化目标通过线性加权和表示，如式(7)所示，通过加权后的优化目标可得到更符合工程实际的优化结果。
[0037]min u
[0038]U= ω j.AsVar (R( ξ , S0, θ )) + ω2.det Σ ( θ )(7)
[0039]其中，CO1和ω2为两种优化目标的加权系数。
[0040]第四步确定约束条件
[0041]约束条件包括试验费用约束和试验变量取值范围约束两部分。
[0042]( I)试验总费用Ct ;试验总费用主要包括三部分:
[0043]a)单位时间内的试验费用C。，包括试验设备折旧、试验人员的工资以及资源(水、电等)消耗等。单位:元/小时；
[0044]b)测试费用单价Cni,单位:元/次；
[0045]c)试验样本单价Cd，单位:元/件。
[0046]因此，试验总费用可表示为:
KK
[0047]C1= n-Cd+hC。+ H /Atk).Cm(8)
|-_J^-j
[0048]( 2 )试验变量取值范围约束
[0049]根据统计评估的需要，试验样本量不小于3 ;为使检测次数满足信息搜集需要又能将检测引入的误差控制在较小范围内，试验设计人员可设置最小和最大检测次数Hlmin和Hlmax ;为保证低应力水平下能获得足够的信息量，低应力下试验时间应高于高应力下试验时间各应力；为减少应力外推的误差，最低应力水平应接近正常应力水平，最高应力水平应在不改变失效机理的前提下尽可能提高。
[0050]约束条件表述为:
【权利要求】
1.一种弹上寿命件的多变量加速贮存试验优化设计方法，其特征在于，该设计方法包括如下步骤: 1)确定弹上寿命件的性能退化模型和加速模型； 2)对弹上寿命件的加速贮存试验确定试验应力施加方式和优化变量； 3)对弹上寿命件的加速贮存试验建立综合优化目标函数； 4)对弹上寿命件的加速贮存试验确定约束条件； 5)根据步骤3得到的优化目标函数和步骤4得到的约束条件，得到弹上寿命件的加速贮存试验优化模型，采用遗传算法求解，得到加速贮存试验的最优方案。
2.根据权利要求1所述的一种弹上寿命件的多变量加速贮存试验优化设计方法，其特征在于，所述步骤I中的性能退化模型为漂移布朗运动模型，如(I)所示:
3.根据权利要求1所述的一种弹上寿命件的多变量加速贮存试验优化设计方法，其特征在于，所述步骤I中的加速模型包括阿伦尼斯(Arrhenius)模型、逆幂率模型或艾琳(Egring)模型,如(2)所示:
4.根据权利要求1所述的一种弹上寿命件的多变量加速贮存试验优化设计方法，其特征在于，所述步骤2中优化变量包括: (1)总样本量η； (2)各应力水平
5.根据权利要求1所述的一种弹上寿命件的多变量加速贮存试验优化设计方法，其特征在于，所述步骤3中建立综合优化目标函数如式(7)所示，
6.根据权利要求1所述的一种弹上寿命件的多变量加速贮存试验优化设计方法，其特征在于，所述步骤4中的约束条件包括试验费用约束和试验变量取值范围约束，其中 (I)试验总费用Ct ;试验总费用主要包括三部分:a)单位时间内的试验费用C。，包括试验设备折旧、试验人员的工资以及水电资源消耗，单位:元/小时； b)测试费用单价Cm,单位:元/次； c)试验样本单价Cd,单位:元/件； 因此，试验总费用可表示为:
7.根据权利要求1所述的一种弹上寿命件的多变量加速贮存试验优化设计方法，其特征在于，根据综合优化目标函数和约束条件得到弹上寿命件的加速贮存试验优化模型，如式(6)所示，
【文档编号】G06F17/50GK103530449SQ201310450501
【公开日】2014年1月22日 申请日期:2013年9月27日 优先权日:2013年9月27日
【发明者】葛蒸蒸, 赵文晖, 李玉伟, 王岩, 马巍 申请人:北京电子工程总体研究所