3. 热机械有限元控制方程与弱形式

封装结构的热机械分析，本质上是在求解一个温度场与位移场相互耦合的边值问题。焊点、铜柱、塑封料、基板这些材料的热膨胀系数差异很大，温度变化会在界面附近产生显著的应力集中。有限元方法要做的，就是把连续介质中的偏微分方程转化为可计算的代数方程组。这一节我们从热弹性本构出发，依次建立热传导方程、虚功原理、加权残量法，最终推导出耦合热-力问题的弱形式，并说明边界条件与初始条件如何设置。

热弹性本构关系

在热弹性问题中，总应变由两部分叠加而成：力学应变与热应变。以各向同性材料为例，小变形假设下总应变张量写为

ε_ij = ε_ij^e + ε_ij^th，其中热应变 ε_ij^th = α(T − T_ref)δ_ij

这里 α 是线膨胀系数，T_ref 是无应力参考温度，δ_ij 是克罗内克符号。需要特别强调：T_ref 的选择直接决定热应力的数值。工程上常取回流焊峰值温度或固化温度作为参考温度，因为在该温度下材料近似处于无应力状态。若参考温度设置错误，整个应力场会整体偏移一个常数项。

由胡克定律，应力与弹性应变的关系为

σ_ij = C_ijkl(ε_kl − α(T − T_ref)δ_kl)

对于各向同性材料，C_ijkl 由杨氏模量 E 和泊松比 ν 确定。封装中常见的塑封料、基板往往呈现正交各向异性，此时 C 矩阵需要 9 个独立弹性常数。本构方程是后续所有推导的物理基础。

热传导方程

瞬态温度场由能量守恒与傅里叶定律导出：

ρc_p ∂T/∂t = ∇·(k∇T) + Q

其中 ρ 是密度，c_p 是比热容，k 是热导率，Q 是内热源。封装分析中 Q 通常来自芯片功耗，可等效为体热源或面热源。稳态分析时左端项为零，方程退化为 ∇·(k∇T) + Q = 0。

热传导方程的边界条件有三类：第一类给定温度 T = T̄（如恒温台）；第二类给定热流密度 −k∂T/∂n = q̄（如芯片发热面）；第三类对流换热 −k∂T/∂n = h(T − T_∞)（如自然对流表面）。初始条件为 T(x,0) = T₀(x)。

虚功原理与加权残量法

力学平衡方程 ∇·σ + f = 0 在强形式下要求场变量二阶连续，这对复杂封装几何很难满足。虚功原理将其转化为积分形式：对任意虚位移 δu，有

∫_Ω σ : δε dΩ = ∫_Ω f·δu dΩ + ∫_Γ t̄·δu dΓ

这一步的物理含义是：内力虚功等于外力虚功。它把对 σ 的求导转移到了虚应变 δε 上，降低了对场函数连续性的要求。

加权残量法提供了更一般的框架。设控制方程记为 L(u) = 0，近似解 u^h 代入后产生残量 R = L(u^h)。选取权函数 w，令 ∫_Ω w·R dΩ = 0。当权函数取为形函数本身时即为伽辽金法，这也是有限元最常用的格式。伽辽金法结合分部积分，就能得到对称的弱形式。

耦合热-力弱形式推导

将热弹性本构代入虚功原理，得到力学弱形式：

∫_Ω δε : C : ε dΩ − ∫_Ω δε : C : α(T − T_ref)I dΩ = ∫_Ω f·δu dΩ + ∫_Γ t̄·δu dΓ

左端第二项就是热载荷项，它把温度场作为已知量引入力学方程。对热传导方程乘以权函数 δT 并分部积分，得到热弱形式：

∫_Ω ρc_p δT ∂T/∂t dΩ + ∫_Ω ∇δT·(k∇T) dΩ = ∫_Ω δT Q dΩ − ∫_Γ δT q̄ dΓ − ∫_Γ δT h(T − T_∞) dΓ

两个弱形式通过温度 T 耦合：温度场影响力学方程的右端热载荷，而力学变形对温度场的反馈在封装热机械分析中通常可忽略（即单向耦合）。因此实际求解时先算温度场，再把温度作为体载荷施加到力学模型上。这种单向耦合假设在大多数封装可靠性分析中是合理的。

边界条件与初始条件设置

边界条件的正确设置是仿真可信度的关键。力学方面，封装底部焊球通常施加对称约束或固定约束以消除刚体位移；对称面施加法向位移为零的条件。热学方面，需要区分芯片发热面（热流密度）、塑封料外表面（对流换热系数 h 通常取 5～25 W/(m²·K)）、以及底面焊球与 PCB 的传热路径。

初始条件方面，瞬态热分析必须给定初始温度场。典型回流焊曲线分析中，初始温度设为 25℃，然后按升温速率加载。参考温度 T_ref 建议取 175℃（无铅回流峰值附近）或材料固化温度，并在整个分析中保持一致。

综合来看，本节的核心逻辑链条是：热弹性本构提供应力-应变-温度关系 → 热传导方程描述温度场演化 → 虚功原理与加权残量法将强形式转化为弱形式 → 耦合弱形式通过温度项连接两个物理场 → 边界与初始条件封闭问题。理解这条链条，才能在有限元软件中正确设置材料参数、载荷步和边界条件，而不是机械地点击菜单。