2. 界面断裂力学基础

封装结构中的分层，本质上是两种材料界面在应力与环境作用下发生裂纹萌生与扩展的过程。要判断界面裂纹是否会扩展、向哪里扩展、扩展多快，必须先把“裂纹驱动力”量化。界面断裂力学提供了四类互补的表征工具：应力强度因子 K、能量释放率 G、J 积分，以及界面裂纹闭合积分。它们分别从应力场、能量、路径守恒和界面本构四个角度描述同一物理过程。

一、界面裂纹的应力场与应力强度因子

对于均匀各向同性材料中的裂纹，裂纹尖端应力场具有 r^-1/2 奇异性，可用单一应力强度因子 K_I、K_II、K_III 描述。界面裂纹则不同：由于两侧材料弹性模量和泊松比不匹配，裂纹面在张开的同时会产生相对滑移，形成 I/II 混合型加载。更关键的是，界面裂纹尖端应力场出现振荡奇异性，其形式为 r^-1/2±iε，其中 ε 为双材料振荡指数：

ε = (1/2π)·ln[(κ₁/μ₁ + 1/μ₂) / (κ₂/μ₂ + 1/μ₁)]

其中 μ 为剪切模量，κ 为 Kolosov 常数：
平面应变 κ = 3 − 4ν，平面应力 κ = (3 − ν)/(1 + ν)

振荡指数意味着裂纹面在尖端附近出现物理上不合理的相互嵌入，因此工程上通常引入一个很小的“接触区”或采用复应力强度因子 K = K₁ + iK₂ 来表征。K₁ 对应张开分量，K₂ 对应面内剪切分量，二者之比决定混合度 ψ = arctan(K₂/K₁)。混合度是界面断裂中最关键的参数之一，因为界面韧性 G_c 强烈依赖 ψ。

二、能量释放率 G：从全局能量角度衡量驱动力

能量释放率定义为裂纹扩展单位面积所释放的弹性能：

G = −∂Π/∂A

Π 为系统总势能，A 为裂纹面积。
在线弹性条件下，界面裂纹的 G 与复应力强度因子满足：

G = (|K|²) / (16 cosh²(πε)) · (κ₁+1)/μ₁ + (κ₂+1)/μ₂

与 K 相比，G 是一个标量，不区分 I/II 型，因此更适合作为扩展判据：当 G ≥ G_c(ψ) 时裂纹扩展。在封装热机械载荷下，G 通常由热膨胀失配、固化收缩、湿气膨胀和残余应力共同贡献。例如，一个典型的塑封料–铜界面在 −55 ℃ 到 125 ℃ 温度循环中，边缘区域的 G 可达到数 J/m² 量级，而界面韧性往往只有 5～30 J/m²，因此边缘和角点是最易分层的区域。

三、J 积分：路径无关的守恒量

J 积分由 Rice 提出，定义为围绕裂纹尖端的回路积分：

J = ∫_Γ (W·n₁ − σ_ij·n_j·∂u_i/∂x₁) ds

W 为应变能密度，Γ 为逆时针包围裂纹尖端的路径，
n 为路径外法线，u 为位移场。

在线弹性条件下 J = G。J 积分的价值在于：它不依赖路径选择，可用于弹塑性材料，并且能通过有限元后处理直接计算。对于界面裂纹，J 积分仍然成立，但需要保证积分路径完全位于某一材料内部或跨越界面时正确计入界面本构。实际封装仿真中，常沿裂纹尖端取多个同心路径计算 J，若各路径结果偏差小于 5%，说明网格和路径设置合理。

四、界面裂纹闭合积分：分离混合模式

J 积分给出总能量，但无法直接分离 I 型和 II 型分量。界面裂纹闭合积分通过虚拟裂纹闭合技术实现模式分解：

G_I = lim(Δa→0) (1/(2Δa)) ∫₀^Δa σ_yy(r)·Δu_y(Δa−r) dr
G_II = lim(Δa→0) (1/(2Δa)) ∫₀^Δa τ_xy(r)·Δu_x(Δa−r) dr
G = G_I + G_II

其中 σ_yy、τ_xy 为裂纹尖端前方应力，Δu_y、Δu_x 为裂纹面相对位移。在有限元中，通常取裂纹尖端后方的节点力和节点位移，用“力×位移/（2×虚拟裂纹扩展面积）”估算 G_I 和 G_II。这样得到的混合度 ψ 可直接与界面韧性测试数据对应。

五、四类表征量的工程选用

在实际封装可靠性分析中，四类量各有适用场景。若关注裂纹尖端局部应力场和混合模式，用 K 和 ψ；若关注扩展判据和能量平衡，用 G；若涉及塑性或复杂路径，用 J；若需要把总能量拆成 I/II 分量并与界面韧性曲线对应，用闭合积分。它们之间并非互斥，而是同一驱动力的不同投影。

理解这四类表征量之后，下一步才能把它们与具体的封装结构、载荷谱和界面韧性数据结合起来，建立可预测的分层判据。本节的重点不是记住公式，而是建立“驱动力—混合模式—扩展判据”之间的对应关系，并在仿真或实验中正确提取和解释这些参数。