4. 热应力理论基础

在温度循环可靠性测试中，TSV（硅通孔）结构由硅基体、铜柱、绝缘层（SiO₂）及阻挡层（Ta/Ti）等多种材料构成。由于各材料的热膨胀系数（CTE）存在显著差异，当温度发生变化时，不同材料层之间会产生不协调的膨胀或收缩，从而在界面处引发热应力。这种热应力是导致铜柱挤出、界面开裂、空洞萌生等失效模式的根本驱动力。

4.1 热膨胀系数（CTE）失配

热膨胀系数（CTE）定义为材料在单位温度变化下的长度相对变化率，单位为 ppm/°C（即 10⁻⁶/°C）。在 TSV 结构中，关键材料的 CTE 对比如下：

CTE 失配的程度通常用失配因子 Δα 表示：

Δα = α_Cu - α_Si ≈ 14 ppm/°C

对于直径 10 μm、深宽比 10:1 的 TSV，在 -55°C ~ 125°C 的温度循环范围内，铜柱与硅孔壁之间的最大热失配位移量可达 0.2 ~ 0.5 μm，这一量级的位移足以在界面处产生显著的塑性应变。

4.2 热应力与应变的基本公式

热应力分析建立在连续介质力学和热弹性理论基础上。以下为 TSV 结构中最核心的公式体系：

4.2.1 热应变公式

当材料受到温度变化 ΔT 时，产生的自由热应变为：

ε_th = α · ΔT

其中：

• ε_th — 热应变（无量纲）
• α — 材料的热膨胀系数（1/°C）
• ΔT — 温度变化量（°C），ΔT = T_current − T_reference

在 TSV 结构中，由于铜柱被硅孔壁约束，实际应变不等于自由热应变，而是由约束产生的机械应变与热应变共同决定：

ε_total = ε_mech + ε_th = ε_mech + α · ΔT

4.2.2 热应力本构方程（广义胡克定律）

对于各向同性线弹性材料，考虑热效应的三维应力-应变关系为：

σ_ij = 2G · ε_ij + λ · ε_kk · δ_ij − (3λ + 2G) · α · ΔT · δ_ij

简化为一维形式（沿 TSV 轴向）：

σ = E · (ε_total − α · ΔT) = E · (ε_mech)

其中：

• σ — 热应力（MPa）
• E — 弹性模量（GPa）
• ε_total — 总应变
• ε_mech — 机械应变（由约束产生）

4.2.3 界面剪切应力公式

在 TSV 铜柱与硅孔壁的界面处，剪切应力分布沿轴向不均匀，两端大、中间小。经典剪滞模型（Shear-Lag Model）给出：

τ(z) = (Δα · ΔT · E_Cu · r) / (2 · t_ox) · [sinh(β · z) / cosh(β · L/2)]

其中：

• τ(z) — 沿轴向位置 z 处的界面剪切应力
• r — 铜柱半径
• t_ox — 绝缘层厚度
• L — 铜柱长度
• β — 剪滞参数，与材料刚度比和几何尺寸相关

4.3 热应力分布特征与失效关联

基于上述公式，TSV 结构在温度循环中的热应力分布具有以下特征：

• 轴向应力（σ_zz）： 铜柱内部以拉/压应力为主，低温时铜收缩受约束产生拉应力，高温时铜膨胀受约束产生压应力。
• 径向应力（σ_rr）： 在铜-硅界面处产生径向挤压或分离应力，影响界面粘接强度。
• 剪切应力（τ_rz）： 集中在铜柱两端（距端面约 10% ~ 20% 长度范围内），是界面开裂的主要驱动力。
• von Mises 等效应力： 用于判断铜柱是否进入塑性屈服，通常在铜柱两端及拐角处达到峰值。

4.4 热应力分析流程

在 TSV 可靠性评估中，热应力分析通常遵循以下步骤：

1. 确定温度载荷： 获取温度循环曲线的上下限（T_max, T_min）、升降温速率及驻留时间。
2. 建立材料本构模型： 铜采用弹塑性模型（考虑随动强化或等向强化），硅和 SiO₂ 采用线弹性模型。
3. 计算等效温度变化： ΔT = T_current − T_stress_free（应力自由温度通常取电镀或退火温度，约 150°C ~ 200°C）。
4. 求解热应力场： 使用有限元方法（FEM）或解析模型计算各应力分量。
5. 评估失效风险： 将计算应力与材料强度准则（如 Tresca、von Mises、最大主应力）进行比较。

4.5 热应力与失效模式的关联