发明名称 一种热振联合载荷下BGA焊点疲劳寿命预测方法
摘要 一种热振联合载荷下BGA焊点疲劳寿命预测方法,其步骤如下:1.采用有限元分析软件建立单个BGA焊点的有限元模型;2.按照工作环境条件对BGA焊点进行热循环仿真;3.根据热循环仿真结果计算热循环载荷下BGA焊点的损伤率;4.在不同的温度条件下对BGA焊点进行随机振动仿真;5.计算热循环载荷影响下由于随机振动载荷导致的BGA焊点损伤率;6.计算BGA焊点在热循环载荷与随机振动载荷下共同作用下总损伤率;7.计算得到BGA焊点在热振联合载荷下的预测寿命。本发明能够提高实际生产中热振联合载荷下BGA焊点的疲劳寿命的预测准确度,具有一定的工程应用价值。
申请公布号 CN103778292A 申请公布日期 2014.05.07
申请号 CN201410031170.6 申请日期 2014.01.23
申请人 北京航空航天大学 发明人 胡薇薇;孟祥坤;孙宇锋;赵广燕;牟浩文
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京慧泉知识产权代理有限公司 11232 代理人 王顺荣;唐爱华
主权项 1.一种热振联合载荷下BGA焊点疲劳寿命预测方法,其特征在于:该方法的具体步骤如下:步骤一:在有限元分析软件中输入所要研究的BGA焊点的初始条件参数,建立出BGA焊点的有限元模型;步骤二:按照BGA焊点实际工作时的环境条件,在有限元分析软件中通过设置热循环的温度范围和每一个热循环周期的时间,包括高低温的保温时间和升温降温时间来对BGA焊点的有限元模型进行热循环仿真;步骤三:热循环仿真完成后,提取BGA焊点的有限元模型的塑性应变量,将其代入Coffin‐Mason模型中,如公式(2)所示,得到BGA焊点在热循环下的疲劳寿命,再利用公式<img file="FDA0000460723130000014.GIF" wi="309" he="128" />计算得到热循环下BGA焊点的损伤率D<sub>TC</sub><maths num="0001"><![CDATA[<math><mrow><msub><mi>N</mi><mi>r</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mrow><mo>(</mo><mfrac><mi>&Delta;&gamma;</mi><msub><mrow><mn>2</mn><mi>&epsiv;</mi></mrow><mi>f</mi></msub></mfrac><mo>)</mo></mrow><mfrac><mn>1</mn><mi>c</mi></mfrac></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,D为损伤率,TTF为在某一载荷作用下BGA焊点的疲劳寿命;N<sub>r</sub>为热循环下BGA焊点疲劳寿命,C是疲劳塑性指数,ε<sub>f</sub>是疲劳塑性系数,△γ为塑性应变量;步骤四:在热循环的温度范围中选择几个关键的温度,在这几个温度下按照实际工作时的环境条件对BGA焊点的有限元模型进行随机振动仿真,随机振动仿真具体流程如下所示:(1)在有限元分析软件中设置受温度影响的参数;(2)利用有限元分析软件对BGA焊点的有限元模型进行模态分析;(3)在步骤(2)的基础上再利用此有限元分析软件对BGA焊点的有限元模型进行随机振动分析;(4)在步骤(3)的分析结果中读取BGA焊点的有限元模型应变最大处的应变响应功率谱密度;(5)把得到的应变响应功率谱密度通过逆傅里叶变换转化为时域应变信息;(6)采用雨流计数法将时域应变信息进行简化,从中获得应变范围及其频率,得到应变范围分布函数;步骤五:根据随机振动仿真得到的结果,利用修正后Coffin‐Manson‐Basquin方程,如公式(3)所示,通过解此方程得到BGA焊点在随机振动下的疲劳寿命;<maths num="0002"><![CDATA[<math><mrow><mfrac><mi>&Delta;&epsiv;</mi><mn>2</mn></mfrac><mo>=</mo><mfrac><msub><mi>&Delta;&epsiv;</mi><mi>e</mi></msub><mn>2</mn></mfrac><mo>+</mo><mfrac><msub><mi>&Delta;&epsiv;</mi><mi>P</mi></msub><mn>2</mn></mfrac><mo>=</mo><mfrac><mrow><msub><mi>&sigma;</mi><mi>f</mi></msub><mo>-</mo><msub><mi>&sigma;</mi><mn>0</mn></msub></mrow><mi>E</mi></mfrac><msup><mrow><mo>(</mo><msub><mrow><mn>2</mn><mi>N</mi></mrow><mi>f</mi></msub><mo>)</mo></mrow><mi>b</mi></msup><mo>+</mo><msub><mi>&epsiv;</mi><mi>f</mi></msub><msup><mrow><mo>(</mo><msub><mrow><mn>2</mn><mi>N</mi></mrow><mi>f</mi></msub><mo>)</mo></mrow><mi>c</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,△ε是总的应变范围,△ε<sub>e</sub>是弹性应变范围,△ε<sub>p</sub>是塑性应变范围,E是弹性模量,σ<sub>f</sub>是应力强度系数,σ<sub>0</sub>是温度导致塑性应力,b是疲劳强度指数,c是疲劳塑性指数,ε<sub>f</sub>是疲劳塑性系数,N<sub>f</sub>为BGA焊点在随机振动下的疲劳寿命;之后再计算得到在热循环的影响下由随机振动导致的BGA焊点的损伤率D<sub>RV</sub>,具体计算过程如下:(1)选择确定的温度;(2)选择该温度下BGA焊点的疲劳损伤系数;(3)根据有限元软件仿真分析得到该温度下随机振动导致的BGA焊点上的应变;(4)根据公式(3)计算该温度下BGA焊点在随机振动载荷下的疲劳寿命;(5)计算该温度保温时间占热循环过程总时间的比值;(6)把步骤(4)计算的结果带入公式(1),得到该温度下随机振动导致的BGA焊点的损伤率;(7)对于其他温度重复步骤(1)~(6);(8)最后将所有温度下由随机振动造成的BGA焊点的损伤率相累加;步骤六:将计算得到的热循环导致的BGA焊点损伤率和热循环影响下由随机振动导致的BGA焊点损伤率相加,得到BGA焊点在热振联合载荷下的总损伤率,如公式(4)所示D<sub>total</sub>=D<sub>TC</sub>+D<sub>RV</sub>         (4)其中,D<sub>TC</sub>是热循环导致的BGA焊点损伤率;D<sub>RV</sub>是热循环影响下由于随机振动导致的BGA焊点损伤率;D<sub>total</sub>为热振联合载荷下BGA焊点的总损伤率;步骤七:将总损伤率代入公式(5),得到热振联合载荷下BGA焊点预测的疲劳寿命,完成热振联合载荷下BGA焊点疲劳寿命的预测;<maths num="0003"><![CDATA[<math><mrow><msub><mi>N</mi><mi>t</mi></msub><mo>=</mo><mfrac><mn>1</mn><msub><mi>D</mi><mi>total</mi></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,D<sub>total</sub>为热振联合载荷下BGA焊点的总损伤率,N<sub>t</sub>为热振联合载荷下BGA焊点预测的疲劳寿命。
地址 100191 北京市海淀区学院路37号