发明名称 一种基于故障物理的电子产品可靠性加速试验方法
摘要 本发明属于电子产品可靠性试验技术,涉及一种基于故障物理的电子产品可靠性加速试验方法。其特征在于:进行电子产品可靠性加速试验的步骤如下:构建电子产品的数字样机模型;热应力状态和振动应力状态分析;确定电子产品的故障机理;确定电子产品的故障物理模型;给定可靠性加速试验的温度条件;计算加速试验时间;计算可靠性加速试验的振动条件;试验实施及试验结果评估。本发明提出了一种基于故障物理的电子产品可靠性加速试验方法,能考虑应力在产品结构中的传递作用,提高了模型的准确度,进而提高了试验结果精度。
申请公布号 CN104392073A 申请公布日期 2015.03.04
申请号 CN201410771569.8 申请日期 2014.12.12
申请人 中国航空综合技术研究所 发明人 任占勇;王昕;陈新;孙瑞锋;李想;王欣;王礼沅;杜鑫
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 中国航空专利中心 11008 代理人 梁瑞林
主权项 一种基于故障物理的电子产品可靠性加速试验方法,其特征在于:进行电子产品可靠性加速试验的步骤如下:1.1、构建电子产品的数字样机模型:数字样机模型是指二维数字样机模型或三维数字样机模型,电子产品包括机箱、支架、模块、电路板和元器件;1.2、热应力状态和振动应力状态分析:采用有限元仿真分析软件分析电子产品在载荷条件下的热应力状态和振动应力状态,载荷条件是指环境载荷和工作载荷;1.3、确定电子产品的故障机理:根据电子产品所承受热应力状态和振动应力状态,结合电子产品的故障模式、影响及危害性分析报告、外场和实验室故障数据,确定电子产品的主要失效位置、失效模块、失效电路板和潜在故障元器件,及其故障机理。1.4、确定电子产品的故障物理模型:根据潜在故障元器件的故障机理,确定电子产品的故障物理模型,设置潜在故障元器件的几何结构参数、材料属性、应力参数、模型修正因子等模型参数;1.5、给定可靠性加速试验的温度条件:高温:被试验电子产品给出的极限温度范围的上限减去10℃~20℃,持续时间为温度稳定时间加上被试验电子产品测试时间;低温:被试验电子产品给出的极限温度范围的下限加上5℃~10℃,持续时间为温度稳定时间加上被试验电子产品测试时间;温度变化率:5℃/min~15℃/min;:1.6、计算加速试验时间T<sub>1</sub>:以被试验电子产品基于环境模拟的可靠性试验剖面中的温度条件和试验时间T<sub>0</sub>为输入参数,计算在给定可靠性加速试验中给出的温度条件下的加速试验时间T<sub>1</sub>,具体方法如下:1.6.1、计算温度条件加速因子τ<sub>V</sub>:根据电子产品的故障物理模型,计算潜在故障元器件在环境模拟的可靠性试验剖面中的温度条件下的损伤,得到环境模拟的可靠性试验剖面中的温度条件持续作用下热疲劳失效的n个潜在薄弱点N<sub>1</sub>,N<sub>2</sub>,…,N<sub>n</sub>,n个潜在薄弱点的首发故障循环数分别记为N<sub>T1</sub>,N<sub>T2</sub>,…N<sub>Ti</sub>,…N<sub>Tn</sub>,i=1,2,…n;计算在给定可靠性加速试验的温度条件的损伤,得到n个潜在薄弱点的首发故障循环数分别记为N'<sub>T1</sub>,N'<sub>T2</sub>,…N'<sub>Ti</sub>,…N'<sub>Tn</sub>;第i个故障点的加速因子τ<sub>Vi</sub>:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&tau;</mi><mi>Vi</mi></msub><mo>=</mo><mfrac><msub><mi>N</mi><mi>Ti</mi></msub><msub><msup><mi>N</mi><mo>&prime;</mo></msup><mi>Ti</mi></msub></mfrac><mo>.</mo><mo>.</mo><mo>.</mo><mo>[</mo><mn>1</mn><mo>]</mo></mrow>]]></math><img file="FDA0000631974080000021.GIF" wi="1643" he="133" /></maths>将n个潜在故障点的加速因子进行算术平均,得到产品温度加速因子τ<sub>V</sub>:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&tau;</mi><mi>V</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mi>&Sigma;</mi><mn>1</mn><mi>n</mi></munderover><msub><mi>&tau;</mi><mi>Vi</mi></msub><mo>.</mo><mo>.</mo><mo>.</mo><mo>[</mo><mn>2</mn><mo>]</mo></mrow>]]></math><img file="FDA0000631974080000022.GIF" wi="1637" he="134" /></maths>1.6.2、计算加速试验时间T<sub>1</sub>:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>T</mi><mn>1</mn></msub><mo>=</mo><mfrac><msub><mi>T</mi><mn>0</mn></msub><mrow><msub><mi>&tau;</mi><mi>V</mi></msub><mo>&times;</mo><msub><mi>&tau;</mi><mn>1</mn></msub></mrow></mfrac><mo>.</mo><mo>.</mo><mo>.</mo><mo>[</mo><mn>3</mn><mo>]</mo></mrow>]]></math><img file="FDA0000631974080000023.GIF" wi="1577" he="142" /></maths>其中,t<sub>1</sub>为给定可靠性加速试验的温度条件时间;1.7、计算可靠性加速试验的振动条件:以基于环境模拟的可靠性试验剖面的振动条件和对应的振动应力施加时间T′<sub>0</sub>为输入参数,计算加速试验时间T<sub>1</sub>内振动应力施加时间T′<sub>1</sub>的振动量值,具体方法如下:1.7.1、计算振动加速因子b<sub>V</sub>:根据电子产品的故障物理模型,计算潜在故障元器件在环境模拟的可靠性试验剖面中的振动条件下的损伤,得到产品在基于环境模拟的可靠性试验剖面中振动条件W<sub>a</sub>下振动疲劳失效的m个潜在薄弱点,分别在振动条件W<sub>01</sub>和振动条件W<sub>02</sub>下,得出的振动疲劳损伤模型分别计算产品在振动耐久故障点损伤率为1时的失效时间,将m个潜在薄弱点在振动条件W<sub>01</sub>下的首发故障时间记为t<sub>V1</sub>,t<sub>V2</sub>,…,t<sub>Vj</sub>…t<sub>Vm</sub>,j=1,2,…m;将m个潜在薄弱点在振动条件W<sub>02</sub>下的首发故障时间记为t'<sub>V1</sub>,t'<sub>V2</sub>,…,t'<sub>Vj,</sub>…t'<sub>Vm</sub>,将第j个薄弱环节点在振动条件W<sub>01</sub>和振动条件W<sub>02</sub>条件下的首发故障时间代入公式:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfrac><msup><msub><mi>t</mi><mi>vi</mi></msub><mo>&prime;</mo></msup><msub><mi>t</mi><mi>vi</mi></msub></mfrac><mo>=</mo><msup><mrow><mo>(</mo><mfrac><msub><mi>W</mi><mn>01</mn></msub><msub><mi>W</mi><mn>02</mn></msub></mfrac><mo>)</mo></mrow><msub><mi>b</mi><mi>vi</mi></msub></msup><mo>.</mo><mo>.</mo><mo>.</mo><mo>[</mo><mn>4</mn><mo>]</mo></mrow>]]></math><img file="FDA0000631974080000028.GIF" wi="1790" he="173" /></maths>得到第j个故障点的常数因子b<sub>vi</sub>;将m个潜在故障点的加速因子进行算术平均,得到振动加速因子b<sub>V</sub>:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>b</mi><mi>V</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>m</mi></mfrac><munderover><mi>&Sigma;</mi><mn>1</mn><mi>m</mi></munderover><msub><mi>b</mi><mi>vj</mi></msub><mo>.</mo><mo>.</mo><mo>.</mo><mo>[</mo><mn>5</mn><mo>]</mo></mrow>]]></math><img file="FDA0000631974080000025.GIF" wi="1652" he="142" /></maths>1.7.2、计算在振动应力施加时间T′<sub>1</sub>下的振动量值W<sub>b</sub>:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>W</mi><mi>b</mi></msub><mo>=</mo><msub><mi>W</mi><mi>a</mi></msub><mo>&times;</mo><msup><mrow><mo>(</mo><mfrac><msup><msub><mi>T</mi><mn>0</mn></msub><mo>&prime;</mo></msup><msup><msub><mi>T</mi><mn>1</mn></msub><mo>&prime;</mo></msup></mfrac><mo>)</mo></mrow><mfrac><mn>1</mn><msub><mi>b</mi><mi>V</mi></msub></mfrac></msup><mo>.</mo><mo>.</mo><mo>.</mo><mo>[</mo><mn>6</mn><mo>]</mo></mrow>]]></math><img file="FDA0000631974080000026.GIF" wi="1715" he="162" /></maths>1.8、试验实施及试验结果评估:综合温度和振动应力条件,并按照综合后的加速试验条件开展试验,试验结束后,按照下式对产品的平均故障间隔时间MTBF在给定置信度下的单侧置信下限θ<sub>L</sub>进行评估:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>&theta;</mi><mi>L</mi></msub><mo>&GreaterEqual;</mo><mfrac><msub><mrow><mn>2</mn><mi>T</mi></mrow><mn>0</mn></msub><msubsup><mi>&chi;</mi><mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>c</mi><mo>)</mo></mrow><mo>,</mo><mrow><mo>(</mo><mn>2</mn><mi>r</mi><mo>+</mo><mn>2</mn><mo>)</mo></mrow></mrow><mn>2</mn></msubsup></mfrac><mo>.</mo><mo>.</mo><mo>.</mo><mo>[</mo><mn>7</mn><mo>]</mo></mrow>]]></math><img file="FDA0000631974080000027.GIF" wi="1715" he="151" /></maths>式中:r是在可靠性加速试验时间T<sub>1</sub>内出现的责任故障数;C是置信度,C=0~1。
地址 100028 北京市朝阳区京顺路7号