3. 焊点热疲劳物理机制

焊点热疲劳并不是单一因素造成的失效，而是“热膨胀系数失配 → 循环塑性应变 → 蠕变-塑性交互 → 裂纹萌生与扩展 → 失效判据与寿命指标”这一条物理链在温度循环下逐步累积的结果。理解这条链，是后续建立仿真模型、选择本构方程和判据的前提。

一、热膨胀系数失配：驱动力的来源

焊点连接芯片与基板，两侧材料的热膨胀系数（CTE）差异显著。以典型倒装焊为例：硅芯片 CTE 约 2.6 ppm/℃，FR-4 基板 CTE 约 16~18 ppm/℃，无铅焊料 SAC305 约 21~23 ppm/℃。当温度从 -40℃ 升到 125℃ 时，芯片与基板的自由伸长量差异可达数微米量级，而焊点高度通常只有几十微米，因此焊点被迫承担这部分不匹配变形。

工程上常用 Engelmaier 经验式估算剪切应变幅：

Δγ = (L_D · Δα · ΔT) / h
其中：
L_D  —— 器件到中性点距离（mm）
Δα  —— 芯片与基板 CTE 之差（ppm/℃）
ΔT  —— 温度循环幅值（℃）
h   —— 焊点有效高度（mm）

示例：取 L_D = 5 mm，Δα = 14 ppm/℃，ΔT = 165℃，h = 0.06 mm，则 Δγ ≈ 0.192，即约 19% 的剪切应变幅。这个量级远超焊料弹性极限，说明焊点必然进入塑性甚至蠕变主导区间。注意：该数值为教学示例估算，实际需结合具体封装尺寸与材料参数。

二、循环塑性应变与蠕变-塑性交互

在温度循环中，焊点经历“升温软化—降温硬化”的反复加载。由于同系温度 T/T_m 较高（SAC305 熔点约 217℃，125℃ 时 T/T_m ≈ 0.82），蠕变与塑性变形同时存在且相互耦合：升温阶段以蠕变为主，降温阶段以塑性屈服为主，二者共同消耗焊点的延性储备。

描述这种交互常用 Anand 本构模型，其核心参数包括：

• s₀：初始变形阻力，SAC305 约 20~30 MPa；
• Q/R：蠕变激活能项，约 9000~11000 K；
• A：指数前因子，约 1e5~1e7 s⁻¹；
• ξ、m、n：应力敏感与应变率敏感指数。

在仿真中，这些参数决定每个温度循环内非弹性应变增量的大小。若参数偏软，会高估损伤；偏硬则低估裂纹萌生时间。因此，参数标定应优先使用与目标温度区间一致的蠕变与拉伸数据。

三、裂纹萌生与扩展

裂纹通常萌生于焊点与铜焊盘之间的金属间化合物（IMC）界面，尤其是芯片侧界面，因为该处应力三轴度高、IMC 本身脆性大。萌生阶段可用非弹性应变能密度 ΔW 作为控制量：当某单元在一个循环内的累积 ΔW 超过临界值 W_c 时，认为该处萌生微裂纹。SAC305 的 W_c 在文献中常取 0.5~2 mJ/mm³ 量级，具体值需通过等温疲劳或热循环试验反标定。

裂纹扩展阶段常用 Paris 律描述：da/dN = C(ΔK)^m，其中 ΔK 为应力强度因子幅值。在焊点中，裂纹多沿界面或焊料内部扩展，扩展路径受 IMC 形貌、晶粒取向和孔洞分布影响。仿真中可通过 VCCT 或内聚力单元模拟界面开裂。

四、失效判据与寿命指标

工程上判断焊点失效通常采用三类判据：

• 电学判据：焊点电阻上升超过初始值的 10%~20%，或出现开路；
• 力学判据：裂纹面积占焊点截面积比例超过 20%；
• 损伤判据：累积非弹性应变能密度达到临界值。

寿命指标方面，常用 Nf（失效循环数）、N50（50% 失效概率寿命）和 B10（10% 失效概率寿命）。在加速试验与仿真之间，需用加速因子 AF 换算，例如 Norris-Landzberg 模型：

AF = (ΔT_test / ΔT_use)^n · (f_use / f_test)^m · exp[Ea/k · (1/T_max,use - 1/T_max,test)]
典型取值：n ≈ 1.9，m ≈ 0.33，Ea ≈ 0.12 eV

示例：若试验条件 ΔT=165℃、频率 1 cycle/h，使用条件 ΔT=80℃、频率 1 cycle/day，取 n=1.9、m=0.33，则 AF 约为 4~6 倍量级。这意味着 1000 次试验循环大致对应 4000~6000 次使用循环。该结果为教学示例，实际需根据材料与载荷谱重新标定。