2. SiP封装中的电-热-力耦合机理

SiP封装把裸芯片、无源器件、基板、焊球、塑封料和散热结构压缩在很小的体积内，多个物理场不再彼此独立：电流产生焦耳热，温度升高改变材料电阻与力学属性，热膨胀差异又反过来改变接触压力与电流路径。理解这些耦合关系，是判断“哪些场必须联合仿真、哪些可以解耦简化”的前提。

一、焦耳热：电→热的第一级耦合

焦耳热来自导体电阻对电流的耗散，单位体积发热率可写为 q = J²ρ，其中 J 为电流密度，ρ 为电阻率。对SiP而言，最需要关注的是高电流密度的RDL走线、微凸点（μBump）、铜柱和焊球。以一条宽20 μm、厚5 μm、长2 mm的RDL铜线为例，若通过0.5 A电流，电流密度约5×10⁹ A/m²，铜电阻率取1.7×10⁻⁸ Ω·m，则单位长度电阻约0.017 Ω/mm，整段功耗约8.5 mW。单条线看似很小，但SiP中往往有数百条电源/地线并行，局部叠加后温升可达十几到几十摄氏度。

焦耳热的关键特征是正反馈风险：温度升高使铜电阻率按约0.39%/℃增大，若散热不足，发热进一步增加。因此电源网络（PDN）设计不能只看压降，还要看电流密度分布是否在凸点、过孔处形成“热点”。

二、热膨胀与热应力：热→力的第二级耦合

SiP由多种热膨胀系数（CTE）差异明显的材料组成：硅约2.6 ppm/℃，铜约17 ppm/℃，FR-4基板约14~17 ppm/℃，塑封料约8~12 ppm/℃，焊料约21~24 ppm/℃。温度变化ΔT时，自由伸长量ΔL = α·L·ΔT，但界面约束使各层不能自由伸缩，于是产生热应变与热应力，典型表达式为σ = E·α·ΔT/(1−ν)（单轴约束近似）。

在回流焊冷却阶段，焊点从约250℃降到室温，ΔT约220℃，焊球与芯片、基板之间的CTE失配会在焊点根部形成剪切塑性应变，这是热疲劳裂纹的主要来源。对底部填充（underfill）结构，胶体可分担部分应变，但若胶体与焊料界面粘附不足，反而会在角部形成应力集中。

三、电迁移：电与力的耦合退化

电迁移是电子风力与静电力共同作用下金属原子沿电流方向迁移的现象，其驱动力与电流密度J和温度T强相关，常用Black方程描述失效时间：MTTF = A·J⁻ⁿ·exp(Ea/kT)，其中n通常取1~2，激活能Ea对铜约0.8~1.0 eV。当电流密度超过约10⁵ A/cm²（即10⁹ A/m²）量级，且结温较高时，电迁移会显著加速，在阴极侧形成空洞、阳极侧形成小丘，最终导致开路或短路。

电迁移与热耦合紧密：焦耳热抬高局部温度，温度又指数级缩短MTTF；同时热应力产生的应力梯度也会驱动原子迁移（应力迁移），两者叠加使失效位置常出现在电流拥挤的过孔拐角或凸点颈部。

四、耦合强弱判定与简化原则

判断耦合强弱，可从三个维度入手：时间尺度（电瞬态纳秒级、热秒级、力小时级）、空间尺度（电流拥挤区与热点是否重合）、敏感度（温度对电阻、应力对迁移率的影响是否超过阈值）。若焦耳热引起的温升小于5℃且电流密度低于10⁵ A/cm²，通常可先做电-热弱耦合，再单独校核热应力；若温升超过20℃或电流密度接近临界值，则必须做电-热-力双向耦合。

五、工程判断的落地方法

在实际SiP项目中，可先建立电-热顺序耦合模型：用有限元或专用工具求解电流密度分布，得到焦耳热源，再求解温度场。若最高温升小于5℃且最大电流密度低于10⁵ A/cm²，可认为电-热耦合较弱，热应力单独用温度载荷计算即可。若温升超过20℃，或热点恰好位于凸点、过孔等应力集中区，则应升级为双向耦合，把温度场反馈给电学模型更新电阻率，同时把热应变反馈给力学模型。

一个实用的检查清单是：① 电流密度是否超过10⁵ A/cm²；② 结温是否超过125℃；③ 焊点塑性应变是否超过0.2%；④ 电迁移MTTF是否低于设计寿命的10倍裕量。四项中任一项越界，就应提高耦合等级，而不是继续用单场近似。

综上，SiP中的电-热-力耦合不是三个独立问题的简单叠加，而是通过焦耳热、热膨胀、电迁移形成闭环反馈。掌握“先判断强弱、再选择耦合等级”的原则，才能在精度与效率之间取得平衡，避免过度仿真或漏判关键失效模式。