返修台上那颗刚被热风枪取下来的BGA放大镜下面真的很难看出裂纹在哪。功能测试偶尔过、偶尔不过X-Ray扫一圈也没发现明显连锡或分层可一到高温环境就掉电。这种问题十有八九出在焊点内部的微裂纹上位置还特别固定——不是阵列的角上就是靠近约束边的外排焊球。与其等整板报废了再切片找原因不如在设计阶段就用Ansys Workbench把BGA焊点最容易坏的位置提前找出来。这篇教程我会按完整流程走一遍先讲清楚焊点失效的物理机理再从几何建模、材料参数、网格划分讲起然后分别跑热循环分析和振动分析最后教你怎么从云图里判断哪些焊球最危险、怎么把结果转化成改板建议。全程基于Ansys Workbench自带功能不需要额外装第三方插件适合刚接触电子可靠性仿真的结构工程师、封装工艺工程师也适合想给BGA失效分析加一道仿真手段的硬件工程师。1. 仿真前必须想清楚的物理问题焊点为什么总裂在角上1.1 热循环失效CTE失配制造的剪切疲劳BGA焊点是芯片和PCB之间唯一的电气和机械连接通道但它的工作环境相当恶劣。芯片是硅热膨胀系数只有2.6ppm/°C左右PCB是玻璃纤维加铜箔面内热膨胀系数大概在14到17ppm/°C焊球本身是锡银铜合金大约22ppm/°C。这几个数字看起来都不大可一旦温度从-40°C拉到125°C温差165°C乘以封装半宽、再乘以失配系数在阵列边缘的焊球上就能造成不小的剪切位移。把焊球想象成两片膨胀速度完全不同的板材之间的胶水温度一变上下两片板往不同方向跑胶水就被反复撕扯。BGA阵列里有个很关键的概念叫DNP也就是距中性点的距离。离封装中心越远的焊球上下表面之间的相对位移越大受力也越大。所以热循环工况下危险焊球几乎永远是阵列对角线最远端那几颗尤其是角上的那一颗或者两颗。这个规律在各类文献和实验里反复出现做仿真之前心里先有这根弦后面看结果就不会慌。1.2 振动失效惯性力作用下的应力集中振动工况的失效逻辑和热循环完全不一样。热循环是温度变化引起的位移失配振动则是板子弯曲导致焊球承受惯性力。PCB在振动激励下会像跳板一样上下弯曲封装体相当于跳板上压着的一块重物。板子弯得越厉害封装质量越大焊球上分到的惯性力就越大。这时候焊球阵列里危险的位置取决于板子的振型如果是一阶弯曲板中央曲率最大中间的焊球可能受力很大如果是高阶扭转或者局部模态危险点又会跑到封装边缘或者靠近固定螺栓的位置。焊球和焊盘交界处是材料刚度突变的地方应力集中天然严重所以振动分析看到的裂纹通常沿焊球根部或IMC界面扩展。1.3 两个危险区不重合必须双线作战这里有个特别容易踩的误区很多人觉得反正焊点最容易坏的就是角上那几颗做一个热循环就算完事了。但实际上振动分析发现的危险焊球可能完全不在封装角上而是出现在靠近PCB固定点那一侧的外排焊球。所以这个项目必须把热循环和振动当作两条独立的分析线来跑各自找各自的最大值位置最后再放在同一张图里对照。只用一种载荷工况下结论大概率会漏掉真正的风险点。2. 从零搭几何模型焊球数量多简化策略比耐心更重要2.1 先定模型边界整板、局部板还是单封装BGA焊点仿真的第一步不是急着画焊球而是先决定模型范围。整板全阵列最真实但一个36×36引脚的BGA就有一千多颗焊球哪怕每颗焊球只分几百个单元整体规模也轻松到几百万节点普通工作站跑热循环的非线性计算会等到怀疑人生。我常用的做法是局部板模型只建BGA封装所在的那一小块PCB板的尺寸取封装外形往外扩20到30毫米四角按真实螺栓孔位置做约束。对找出阵列里哪个焊球最危险这个目标来说局部板模型已经完全够用。为什么因为焊球之间的相对应力分布主要由封装和局部板的刚度匹配决定板边界稍远一点对焊球应力梯度的影响就迅速衰减了。如果机器实在跑不动还可以考虑1/4对称模型。这里有个前提热循环分析可以用对称模型因为温度载荷是均匀的结构也基本对称但振动分析我建议保留全模型因为激励方向和约束点的组合会破坏对称性。你要是为了省事把振动也做成1/4对称很可能把本该出现的反对称模态给硬生生抹掉。2.2 焊球别用理想球体按这个形状去切很多初学者会在三维软件里画一个完美的球然后放在焊盘上。这样建模特别好看但做出来的结果往往偏保守而且网格也不好画。实际回流焊之后BGA焊球是上下被焊盘压扁、腰部微鼓的截锥形状上下表面各有一块被焊盘润湿形成的平面。精确还原这个形貌需要扫描电镜切片工程上完全没必要。我通常直接在SpaceClaim里用旋转体拉一个截锥上下面直径0.33毫米腰部最大直径0.42毫米高度0.35毫米上下边缘各倒一个0.03~0.05毫米的圆角。圆角一定要保留后面热循环能不能收敛很多时候就看这个圆角。你要是不倒圆角焊球根部就是个90°尖角应力在那里会随着网格加密无限增大也就是所谓的应力奇异然后非线性求解器各种不收敛最后你只能在论坛上找为什么我的热循环算到一半就飞了的答案。2.3 一套可以直接用的BGA结构参数为了方便说后面的步骤我直接给一套我常用的参数你可以照抄进模型里项目取值封装外形15mm × 15mm × 1.2mm芯片尺寸8mm × 8mm × 0.3mm基板厚度0.36mm局部PCB尺寸50mm × 50mm × 1.6mm焊球阵列19 × 19共361颗球间距0.8mm焊球高度0.35mm焊球腰部直径0.42mm焊盘直径0.33mm焊盘厚度0.03mm这套参数对应的是比较典型的0.8mm pitch BGA。焊球数量看着多但每颗焊球的几何几乎一样在SpaceClaim里画好一颗剩下的用阵列复制建模工作量并不大。注意阵列复制之后一定要用布尔操作把相互重叠的面清理掉否则后面网格划分会出现一堆莫名其妙的“退化单元”。3. 材料参数与网格划分两个分析对模型的需求完全不同3.1 焊料热循环用Anand粘塑性模型振动用线弹性就够了BGA焊料在热循环工况下绝不是简单的线弹性材料。温度变化慢、保温时间长焊料会发生显著的蠕变和应力松弛这是焊点低周疲劳的主要损伤来源。这时候如果你只给一个弹性模量和泊松比算出来的应力会比实际高很多危险位置的判断也会失真。所以热循环分析里SAC305焊料要用Anand粘塑性本构模型。Anand模型能同时描述塑性变形和蠕变参数在Ansys Workbench的材料库或者Engineering Data里可以直接输入。下面是我常用的一组SAC305 Anand参数来自公开文献的标定结果用于方案对比足够参数符号取值初始变形阻力s045.9 MPa激活能/玻尔兹曼常数Q/R7460 K指前因子A4.61e6 1/s应力乘子ξ1.5应变率敏感指数m0.0942硬化常数h05840.49 MPa饱和应力常数ŝ58.3 MPa饱和应力应变率指数n0.015硬化指数a1.86这里必须提醒一句Anand参数不是随便拼凑的s0、A、Q/R这些参数必须来自同一组标定结果。我见过有人拿论文A的s0配上论文B的A算出来的蠕变应变直接翻倍。如果你对材料批次没底就留出安全系数按偏保守的取值来。振动工况就宽松多了激励是高频小应变焊料整体还处于弹性范围用线弹性本构完全可以。我习惯在同一个模型里给焊料定义两套材料一套弹性参数加Anand用于热循环一套纯弹性参数用于模态和随机振动。切换的时候只要在模型里替换材料赋值不需要重建几何。3.2 其他材料的参数怎么给除了焊料模型里还有硅芯片、塑封料、基板、铜焊盘和PCB。这些材料的参数直接影响载荷如何传递到焊球不能随便填。下表是我常用的常温参数材料弹性模量(GPa) / 其他泊松比CTE(ppm/°C)硅芯片1310.282.8塑封料EMC20Tg以下0.3015Tg以下Tg以上约40~60铜焊盘1100.3417SAC305焊料5125°C125°C降到约350.3622.5PCB(FR-4含铜层)EXEY22EZ10GXY8GXZGYZ4NUXY0.11NUXZNUYZ0.28XY向16Z向60PCB我用正交各向异性这一点经常被忽略。FR-4板的面内刚度和厚度方向刚度差很多如果你全都用一个各向同性的弹性模量振动模态频率很容易算出偏差热循环中板的翘曲趋势也会不对。Z向CTE给60ppm/°C和面内16ppm/°C差别很大这是为了模拟玻纤布和树脂在厚度方向上的真实膨胀特性。还有塑封料的Tg问题。很多环氧塑封料的玻璃化转变温度在150°C左右超过Tg之后模量会断崖式下跌CTE也会明显上升。如果你的温度循环上限在125°C用一组Tg以下参数没问题如果循环上限拉到150°C以上就必须给双段CTE和分段模量否则高温段的应力会严重失真。3.3 网格怎么切才能让热循环收敛、让模态不失真焊球是模型的灵魂网格策略必须围着它转。我的做法是每颗焊球沿高度分4层周向分8份腰部圆角处局部网格尺寸0.03~0.05毫米。焊球和上下铜焊盘之间不用接触对直接把网格节点共享这样热循环非线性计算最稳。用接触的话每个接触面上的状态变化都会拖慢收敛而且接触参数没设好还会出现抖动发散。PCB和封装体本身可以用六面体主导网格尺寸控制在0.5~1毫米。焊球附近可以做局部尺寸控制让网格从0.05毫米平滑过渡到1毫米。按19×19阵列算下来整模型大概会有一百五十万到两百万个节点。热循环跑一轮可能需要三五个小时甚至更久这取决于你的CPU核心数和内存。如果机器不够可以退一步只对角落和靠近约束边的焊球加密到0.05毫米内部焊球放宽到0.1毫米危险位置判断的精度不会损失太多。振动分析对网格没有热循环那么敏感因为线弹性计算收敛稳定。但要注意模态分析的结果精度取决于整体质量分布和刚度分布PCB的网格不能太粗我一般控制在1毫米。这样前20阶模态的频率误差可以压到5%以内对工程判断足够了。4. 第一道工况热循环分析——重点在应变范围不是温度场4.1 温度载荷曲线三个循环起步热循环分析在Ansys Workbench里用Transient Structural做。很多人上来就想做热-结构耦合先算温度场再映射到结构上。对于BGA焊点这种小尺度结构焊球和铜的热导率高温度均匀性极好整板同时升降温的假设完全合理所以直接对所有体施加一个随时间变化的Body Temperature就行。温度曲线我推荐按这个来从25°C参考温度出发以14°C/min降温到-40°C保温15分钟再以14°C/min升温到125°C保温15分钟再回到25°C。这是一个完整循环。然后把这个循环重复三次。为什么要三个循环因为第一个循环里焊料的初始蠕变状态还没稳定应力应变曲线会带一点棘轮效应从第二个循环开始才会进入稳定的滞回环。我们关心的是稳定循环里的应变范围所以取第三个循环的数据。分析设置里子步控制一定要打开自动时间步。升温段和降温段的初始步长给30到60秒最大步长不超过120秒保温段步长可以放到300到600秒。这样既保证温度转折处捕捉准确又不会因为时间步太小导致总步数爆炸。4.2 边界条件只约束刚体位移热循环分析里的边界条件是很多人翻车的地方。如果你的目的是看焊球相对应力分布就不要把PCB四角螺栓孔全部固定死。板子受热要膨胀你把它死死摁住算出来的热应力全是约束反力造成的假象焊球本体的应力分布会被严重干扰。我的做法是在四角螺栓孔处分别约束Z向位移然后选其中一个螺栓孔同时约束X和Y再选对角方向的一个螺栓孔约束另一个水平方向。这样整块板既能自由热胀冷缩又不会产生刚体漂移。没有螺栓孔就把PCB底面的四角做同样处理。关键原则很简单约束住刚体运动剩下的全交还给物理。关于大变形开关我个人习惯先关掉NLGEOM跑通一遍。热循环虽然涉及蠕变但焊球几何变形没有大到需要几何非线性的程度关掉之后收敛容易很多。等模型稳定了再试着打开NLGEOM对比一下结果如果焊球应变范围差异小于5%以后就可以放心关着用。4.3 后处理不能只看应力要提取循环应变范围热循环结束后云图里最显眼的是等效应力。但我要说判断焊点危险位置主要别看等效应力要看应变。焊料是弹塑性材料应力到达屈服后基本不再增加可塑性应变和蠕变应变却在不断累积这才是损伤的根源。具体操作是这样求解完成后插入Equivalent Total Strain结果找到第三个循环的低温端和高温端两个时间点分别读取应变值两者相减得到该循环的等效应变范围Δεeq。如果你想看整个阵列的分布那就分别导出两个时刻的云图在结果里做减法运算。或者偷懒一点直接看Equivalent Plastic Strain和Equivalent Creep Strain的累积云图看哪个焊球位置数值最大。累积量虽然不能直接当循环应变范围用但定位危险焊球的位置是完全靠谱的。实际算出来的结果里最大Δεeq的位置毫无悬念会在阵列角落。你还会发现同一颗焊球上靠近芯片一侧的根部应变往往比靠近PCB一侧更大这是因为基板和塑封料的约束比PCB更强。这个细节后面看切片验证的时候会对得上。4.4 一个粗略的寿命估算公式如果你不满足于只看位置还想大概知道焊点能活多久可以用Coffin-Manson低周疲劳公式粗算一下。[ N_f 0.5 \left( \frac{\Delta \varepsilon_{eq}}{2 \varepsilon_f} \right)^{-\frac{1}{c}} ]对SAC305粗取εf0.5c-0.6到-0.7。这里面不同文献取值差异很大所以我建议你只拿它做相对寿命评估别直接输出到测试报告里当绝对寿命。更讲究一点的做法是用Darveaux方法提取焊球界面单元的平均蠕变应变能密度增量再代入起裂和扩展寿命公式。那个精度更高但后处理复杂度也上了一个台阶初学阶段用Coffin-Manson足够。5. 第二道工况振动分析——模态开路随机振动确认5.1 为什么不能直接拿静力分析当振动分析很多没做过振动仿真的朋友会问我直接给封装施加一个加速度比如30g算一次静力不行吗答案是不行而且会错得离谱。振动环境的破坏力来自共振放大。如果激励频率刚好接近结构某阶固有频率实际响应加速度可以是输入的好几倍如果离得远响应又可能小到可以忽略。静态30g完全无法体现这种频率相关性。所以正确路线是两步走先做模态分析搞清楚结构的固有频率和振型再做随机振动分析把真实的功率谱密度PSD作为激励得到统计学意义上的应力响应。在Ansys Workbench里这两个分析模块是串联的随机振动分析必须挂在模态分析下面作为子分析。5.2 约束、PSD谱、阻尼比怎么设模态分析的约束和热循环保持一致在PCB四角螺栓孔加Fixed Support。提取阶数我习惯给前20阶或者覆盖到2000Hz。对于50mm×50mm的局部板前几阶模态一般出现在200到800Hz之间20阶足够看到我们关心的板弯曲和扭转模态了。随机振动的激励方向选择Z向也就是PCB法向这是对BGA焊球最危险的方向。面内激励会让PCB产生面内变形但面内刚度大、响应小初判时可以先不管。PSD谱如果你手上没有客户标准先用平直谱0.01g²/Hz、20到2000Hz。算一下总均方根加速度大概是sqrt(0.01×1980)≈4.45g这个量级符合很多车载电子产品的实测环境。注意单位是g²/Hz不是g/Hz填错之后整个应力结果都会对不上。阻尼比我取0.03这是PCB组件常见的模态阻尼范围取值太保守会放大共振峰太激进又会把应力压得太低。5.3 用RMS应力云图定位危险焊球随机振动求解完成后Workbench会输出1σ等效应力RMS应力这是统计意义上的标准差。云图里你会看到应力主要集中在PCB固定点附近和封装边缘。这时一定要把PCB和封装体隐藏掉只看焊球层否则你很容易被PCB上那个大大的红色区域带偏。从经验看1σ应力最大的焊球通常出现在靠近固定螺栓的那一侧最外排或者在封装角部。为什么会这样因为固定螺栓附近PCB弯曲曲率大封装边缘又是刚度突变点两重因素叠加让应力集中特别明显。工程上评估高周疲劳用3σ准则峰值应力约等于3倍RMS应力出现概率99.73%。把3σ应力和SAC305的高周疲劳极限做对比疲劳极限大约在20到30MPa。如果3σ应力超过这个范围说明焊点在高周振动下存在明显的疲劳风险如果1σ应力就超过了疲劳极限那这个设计基本可以判定不合格。5.4 热-振联合分析预应力模态是不是必须做温度会改变材料的弹性模量也会在结构里留下热应力而热应力会通过应力刚化效应改变结构的固有频率。原理上最严格的做法是先把热循环最后一个时刻的温度状态作为预应力做一次Static Structural再把它作为预应力载荷链接到Modal里去算模态和随机振动。在Workbench里这个流程叫Pre-Stress Modal操作上就是把Modal分析的Pre-Stress设置为Static Structural分析。但我的实际经验是初版本趋势评估真的没必要上预应力模态。热应力对焊球危险位置的排序影响很小主要是让固有频率移动几个百分点。你花大半天把预应力链路搭好最终得到的危险焊球位置和分开算几乎没有区别。我建议你先分开算两轮等需要输出精确频率报告的时候再做预应力耦合。6. 结果对照与工程判据怎么判断要坏和还能扛6.1 一个参考算例两种工况的危险焊球不重合下面我用一个参考算例收束前面的步骤。假设是一颗19×19、0.8mm pitch的BGA热循环取三个循环的Δεeq云图最大位置在左下角DNP最大的那颗焊球靠近芯片一侧根部Δεeq大约1.2%。振动分析采用平直PSD谱1σ应力最大位置则出现在靠近PCB固定螺栓那一侧的最外排焊球不在封装角上1σ应力大约8MPa3σ应力大约24MPa刚好压在高周疲劳极限附近。工况主导失效模式危险焊球位置损伤指标热循环蠕变疲劳低周角部DNP最大球芯片侧根部Δεeq≈1.2%振动高周疲劳螺栓侧最外排球焊盘界面附近1σ≈8MPa3σ≈24MPa这个例子完美说明了为什么要两条线同时跑如果你只做了热循环就会得出加强封装角部的结论完全漏掉螺栓侧外排焊球的风险如果你只做了振动又会忽略掉热循环带来的低周蠕变损伤。把两张云图叠在一起看你才会发现这个设计真正的软肋是角部焊球抗热循环、螺栓侧焊球抗振动两边要分别处理。6.2 我常用的几个工程阈值热循环Δεeq超过1%这个焊球就要重点跟踪超过2%基本属于设计明显不合理要不换材料、要不改封装、要不加底部填充。振动1σ应力低于疲劳极限的1/3可以认为3σ应力仍低于疲劳极限设计相对安全1σ应力超过疲劳极限必须整改。危险位置看趋势焊球和铜焊盘交界处因为材料突变和几何圆角应力会表现出网格敏感性直接用界面节点应力卡阈值不科学要看离界面一个单元距离的焊球本体应力或应变。大家注意仿真给的是相对危险程度排序而不是绝对寿命预言尤其是热循环寿命公式里常数项差异很大。做研发评审时我一般说右上角这颗焊球损伤指标是阵列平均值的三倍建议优先优化工程师们马上就能理解。6.3 怎么把仿真结论转成改板建议热循环危险核心思路是降低焊球承受的热失配位移。常见手段包括减小DNP也就是缩小焊球阵列的覆盖范围、选用低CTE的PCB材料、给BGA底部填充underfill、或者调整焊球的合金成分和直径。底部填充对热循环的改善效果非常显著因为underfill会把一部分剪切载荷转移到封装体上焊球本体损伤能下降30%到50%。振动危险核心思路是降低焊球处的动态响应。最常见的手段是增加PCB支撑点或调整螺栓位置让危险焊球远离固定点的高应力区其次是增加PCB厚度提高弯曲刚度也可以给封装和PCB之间点胶或填充缓冲材料提高阻尼。如果两个工况的危险焊球刚好都在同一个位置说明这个焊球正在承受热-振叠加的双重损伤。这种情况下只优化一个维度不行优先考虑underfill加局部支撑的组合方案同时把PCB的约束点调整一下尽可能让角部焊球的振动应力降下来。7. 这个流程里我踩过的几个坑以及当时的排查思路7.1 热循环算着算着就发散热循环不收敛太常见了尤其是第一次跑Anand模型的时候。我的排查顺序是这样的先看约束是不是哪个螺栓孔被固定死了再看网格焊球圆角处有没有出现锐角单元最后看时间步转折处步长是不是太大。有一次我折腾了两天最后发现问题出在一颗焊球复制时没做布尔合并焊球和焊盘之间隔了一层面共节点根本没建立起来整颗焊球在那里自由震动。所以建议大家建模完成后先单独选中一颗焊球检查一下共享拓扑确认它和上下焊盘是一体的。7.2 模态分析出现负频率模态分析前6阶频率如果接近0那就是刚体模态说明结构约束不够。BGA模型常见的问题是焊球和焊盘之间用了接触对而不是共节点接触面在模态分析里没有锁紧整个封装在振动方向上处于半自由状态。解决方法是把接触全部替换成Bonded或者直接共节点。如果前6阶频率是负值这个比较少见通常和求解设置有关检查一下求解器是否用了默认的迭代法必要时切到直接法。7.3 Anand参数看起来对但结果整体偏大有一段时间我图省事从不同的论文里各取了一组Anand参数凑在一起算出来的蠕变应变高到离谱所有焊球都快废了。后来对回文献才发现Anand模型的9个参数是内在耦合的同一组标定数据下的参数才有意义。从那以后我就在本地存了一组完整的SAC305 Anand参数每次新建项目都直接调用同一组除非拿到焊料供应商的实测标定数据否则不用第二套。这一点你们一定要记住。7.4 随机振动的PSD谱单位和方向搞错随机振动模块里加速度谱的单位有g²/Hz和(m/s²)²/Hz两种填错之后RMS应力会差将近100倍。我在一个项目里把0.01g²/Hz误填成0.01m/s²/Hz算出来应力小得可怜差点把一个实际有风险的设计判成合格。自查方法是看总Grms填完谱之后软件会显示总均方根加速度如果3g到10g范围还算合理如果显示0.1g以下大概率是单位错了。还有PSD激励的方向要确认你加在Z向不少人是默认加在X向的算完一看云图特别均匀那是因为面内刚度大、响应小完全没反映出真实危险。7.5 后处理时只看了总体变形云图被板中心的红色区域带偏随机振动结果里总变形云图的最大值几乎都在PCB中心因为板的一阶弯曲变形肯定中间最大。但焊球应力最大的位置未必在中心而是在封装边缘和固定点附近。新手特别容易盯着总变形云图发半天呆得出板中心底下那颗焊球最危险的结论。正确操作是在结果树里插入一个Named Selection把所有焊球都选进去后处理时只显示这个集合。这样你看到的永远是焊球层自己的应力PCB中间那些大红色块就不会干扰判断了。还有一个老生常谈的问题焊球和铜焊盘的交界处应力有很强的网格敏感性。网格加密两倍界面应力可能上涨30%甚至更多。所以做危险位置判断时不要拿界面第一层节点的应力值直接卡阈值要看焊球本体第二层单元的平均应力。虽然整个项目到最后也做不到绝对精确的寿命预测但至少能把你的注意力从361颗焊球缩小到角部那五六颗和螺栓侧那七八颗这就值回票价了。我后来做方案评审必带热循环和振动两张焊球应力云图方便直观地说明这版设计改哪里、为什么改——比干巴巴地念测试数据有说服力得多。