基于Darveaux模型的BGA焊球温循寿命仿真预测
在电子产品可靠性工程中，温度循环（Thermal Cycling）测试几乎是所有BGA、CSP、SiP等封装形式必须经历的“生死考验”。
为什么要做温循测试？ 因为现实世界中，芯片会经历极端的环境温度变化：汽车电子从-40℃的寒冬到125℃的发动机舱，消费电子从空调房到烈日下的户外，5G基站一年四季都要经受昼夜温差……这些温度循环会让芯片与PCB板因热膨胀系数（CTE）严重失配而产生反复的剪切应力，最终导致BGA焊球发生疲劳断裂。这正是BGA封装最常见的失效模式之一。
传统温循测试的做法是：做出真实样品，放入温循箱，按照JEDEC JESD22-A104等标准设定温度范围、循环速率和循环次数（常常需要2000~3000小时甚至更长），然后通过电测和剖面分析判断焊球是否开裂。
听起来严谨，但实际操作却暴露出明显劣势： 周期极长、成本高昂、无法在设计早期介入。一旦试验失败，往往已经进入量产阶段，改版成本巨大。更糟糕的是，它只能告诉你“坏了”还是“没坏”，却很难告诉你“为什么坏”“还能优化到什么程度”。
而有限元仿真（FEA）技术的应用，正在彻底改变这一局面。 它可以在产品设计阶段，就通过建立精准的有限元模型，模拟焊球在完整温度循环过程中的应力-应变响应，计算累积蠕变应变能密度，再结合Coffin-Manson、Engelmaier或Darveaux等疲劳寿命模型，提前定量预测BGA焊球的温循寿命。
接下来，就以一个典型的大尺寸FCBGA为例，手把手带你建立有限元模型，计算焊球的温循寿命。首先我们建立一个简化但保留主要结构的BGA焊接在PCB上的几何模型，如下图所示。
完成材料的赋予，下一步就是加载随时间变化的温度载荷。一般来说可以参考测试国标，也可以先测量温循测试中焊球附近的实际温度再施加到模型上。为了得到单个温度循环内焊球稳定的损伤增量，温循仿真一般需要计算3~5个循环，并取最后一个循环的损伤增量计算疲劳寿命。本文采用的-40~125℃的标准温度循环，如下图所示。
考虑到焊球在温循工况下的总应变较大，需要激活大变形分析开关。计算得到的BGA在温循过程中的动态变形结果如下图所示（Z向变形放大系数10）。可以看到焊球在温循过程中受到循环往复的热应力，并且边缘焊球受到的拉扯更明显。
选择角部的一个焊球做单独分析，提取其在温循过程中的蠕变应变能演化过程，如下图所示。可以看到焊球的蠕变应变能密度最大的地方出现在焊球的对角线上，特别是与焊盘接触的位置。
随着温循的持续，微裂纹最开始在蠕变应变能最大的位置产生，然后进一步扩展，直至扩展到整个与焊盘的接触面上，致使焊球失效断裂，如下面实物剖面图所示。
选择蠕变应变能最大位置绘制随时间变化的蠕变应变能变化，可以认为最后一个循环蠕变应变能达到了稳定增加，取最后一个循环的蠕变应变能增量进行焊球疲劳寿命计算。
在实际工程中，一般不会按每个单元计算寿命，而是计算焊球与焊盘接触的那一层单元的平均蠕变应变能增量，进而计算焊球的整体疲劳寿命。用第3个循环末的结果线性做差减去第2个循环末的结果，提取与上焊盘接触的焊球整层单元的蠕变应变能增量总和为：0.119，整层单元的体积为0.0669，得到平均蠕变应变能密度增量为1.779。
对于常用的 SAC305焊料，可用业界经典的 Darveaux 模型计算焊球的温循疲劳寿命，模型具体内容及典型参数如下图所示：
代入仿真得到的平均蠕变应变能密度增量1.779，焊球/焊盘直径0.8mm，得到焊球的温循疲劳寿命为：裂纹萌生寿命54.5 cycles，裂纹扩展寿命304.2 cycles，焊球总温循寿命为358.7 cycles。
如果采用温箱进行实际温循测试的话，需要进行358.7h，即半个月才能得到结果，按照带多通道在线电性能监测的温箱设备使用费用1000元/小时来算，需要花费的费用为36万元。而通过有限元仿真预测焊球疲劳寿命的成本几乎可以忽略不计。
有限元仿真不是要取代真实的温循试验，而是让试验从“盲目试错”升级为“有目标的验证”。通过仿真，我们可以提前看到焊球的应力集中位置、蠕变应变能累积规律，以及不同设计参数（焊球直径、高度、材料、Underfill、基板结构等）对寿命的影响，从而实现真正的设计优化。
在芯片封装越来越小型化、高密度、异构集成化的今天，传统靠“做出来再试”的可靠性流程已经越来越难以满足苛刻的量产需求。
希望这篇文章能帮你打开一扇从“试验驱动”转向“仿真驱动”可靠性设计的窗口。 如果你在项目中也遇到焊球寿命不达标、温循试验反复失败等问题，欢迎在评论区留言你的具体痛点，一起交流讨论。
感谢阅读！欢迎关注“先进电子封装仿真”。
【免责声明】本文仅为有限元仿真方法的一般性技术分享与讨论，所举案例均为常见行业示例，不针对任何特定公司或产品。如与贵公司实际项目存在相似之处，纯属正常技术共性，特此说明。
