发明名称 基于经验加速模型的退化量分布参数建模外推方法
摘要 本发明公开了一种基于经验加速模型的退化量分布参数建模外推方法,包括以下几个步骤:步骤一、选择对应加速寿命试验的经验加速模型;步骤二、建立退化量分布参数与应力关系的退化量分布参数加速模型;步骤三、利用所得退化量分布参数加速模型外推目标应力下的退化量分布参数值。本发明采用已经在加速寿命试验领域被大量研究和应用证明的经验加速模型进行建模,可以有效利用产品研究过程中的相关试验信息和相关产品的试验信息,提高了试验结果的可信性。本发明以模型推导为基础,避免了小样本情况下,单纯采用回归方法描述分布参数与应力关系时,准确性难以保证的问题。
申请公布号 CN103488826A 申请公布日期 2014.01.01
申请号 CN201310415072.8 申请日期 2013.09.12
申请人 北京航空航天大学 发明人 党香俊;姜同敏;封雷;孙富强
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 官汉增
主权项 1.一种基于经验加速模型的退化量分布参数建模外推方法,其特征在于,包括以下几个步骤:步骤一、选择对应加速寿命试验的经验加速模型;如果对产品施加的加速应力为电压应力,选择逆幂率模型;此时用于加速寿命试验的逆幂率模型表示为:ξ=Av<sup>-C</sup>     (1)式中:ξ是某寿命特征;A是一个正常数;C是一个与激活能有关的正常数;v是应力;步骤二、建立退化量分布参数与应力关系的退化量分布参数加速模型;用D<sub>f</sub>表示产品的退化失效阈值,则产品的寿命特征为:ξ=f(D<sub>f</sub>)<sup>-1</sup>     (2)式中f(.)为性能退化函数的通用表达式,此处采用隐函数的形式;理论情况下,退化失效阈值取任意正值,相应地有:t=f(d<sub>f</sub>)<sup>-1</sup>,d<sub>f</sub>≥0     (3)式中,变量d<sub>f</sub>表示退化失效阈值D<sub>f</sub>在退化量值为正的范围内连续变化,变量t表示相应的退化寿命ξ也是连续变化的时间量值;逆幂率模型描述了产品寿命与所受应力的关系,并不以产品寿命的变化而改变,所以在同一寿命时刻,即t=ξ时,公式(1)和公式(2)结合为:f(d<sub>f</sub>)<sup>-1</sup>=Av<sup>-C</sup>     (4)用通用退化量M代替d<sub>f</sub>,有:f(M)<sup>-1</sup>=Av<sup>-C</sup>,M=d<sub>f</sub>     (5)对公式(5)取对数得:ln(f(M)<sup>-1</sup>)=lnA-Clnv     (6)公式(5)和公式(6)描述了退化量分布参数的与应力的关系;假设寿命加速模型具有如下通用形式:<img file="FDA0000381098740000011.GIF" wi="1189" he="63" />其中<img file="FDA0000381098740000012.GIF" wi="113" he="67" />为应力S的函数,在t=ξ时,结合公式(7)和公式(3),有公式(4)和公式(5)所对应的通用表达形式:<img file="FDA0000381098740000022.GIF" wi="1188" he="78" /><img file="FDA0000381098740000023.GIF" wi="1269" he="76" />步骤三、利用所得退化量分布参数加速模型外推目标应力下的退化量分布参数值;根据退化轨迹获得退化函数f(.),应用公式(6)采用加速应力下不同应力下退化量分布参数的数据,回归得到第i时刻加速模型公式中的参数,对任意第i时刻的正常应力下退化量<img file="FDA0000381098740000024.GIF" wi="87" he="72" />的外推,将正常应力值v<sub>0</sub>代入此时的外推公式:<maths num="0001"><![CDATA[<math><mrow><mi>ln</mi><mrow><mo>(</mo><mi>f</mi><msup><mrow><mo>(</mo><msubsup><mi>M</mi><mi>i</mi><mn>0</mn></msubsup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mo>=</mo><mi>ln</mi><msub><mi>A</mi><mi>i</mi></msub><mo>-</mo><msub><mi>C</mi><mi>i</mi></msub><mi>ln</mi><msub><mi>v</mi><mn>0</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths>求解得到。
地址 100191 北京市海淀区学院路37号