2. 封装互连的电磁场物理基础
封装互连的全波电磁场分析，本质上是在求解一个带边界条件的偏微分方程问题。要理解仿真工具在做什么、为什么某些设置会显著影响结果，必须回到麦克斯韦方程组及其派生形式。本节从物理方程出发，逐步过渡到全波求解的数学本质，为后续建模与结果判读建立基础。
一、麦克斯韦方程组：全波分析的出发点
在封装互连的尺度下（特征尺寸从微米到毫米，信号频率从几百 MHz 到数十 GHz），电磁场随时间变化，必须使用完整形式的麦克斯韦方程组。微分形式如下：
∇ × E = -∂B/∂t （法拉第定律）
∇ × H = J + ∂D/∂t （安培-麦克斯韦定律）
∇ · D = ρ （高斯定律）
∇ · B = 0 （磁通连续性）
其中 E 为电场强度（V/m），H 为磁场强度（A/m），D 为电位移矢量（C/m²），B 为磁感应强度（T），J 为电流密度（A/m²），ρ 为自由电荷体密度（C/m³）。前两个方程描述电场与磁场的相互激发，是波动传播的核心；后两个方程描述场的散度约束。
在封装互连中，我们关心的物理量通常是电压、电流、S 参数和阻抗，它们都可以由 E 和 H 通过积分或比值得到。因此，全波求解的目标就是先求出空间中的 E 和 H 分布，再提取这些工程参数。
二、本构关系：材料如何影响场
麦克斯韦方程组本身不完整，还需要本构关系将 D 与 E、B 与 H、J 与 E 联系起来。对于封装中常见的线性、各向同性介质：
D = εE = ε₀εᵣE
B = μH = μ₀μᵣH
J = σE
其中 ε₀ ≈ 8.854×10⁻¹² F/m，μ₀ = 4π×10⁻⁷ H/m，εᵣ 为相对介电常数，μᵣ 为相对磁导率，σ 为电导率（S/m）。封装互连中，介质层（如 FR4、BT 树脂、ABF 膜）的 εᵣ 通常在 3.0～4.5 之间，损耗角正切 tanδ 在 0.005～0.025 之间；铜导体 σ ≈ 5.8×10⁷ S/m。
本构关系决定了场在材料中的传播速度、衰减和波长。例如，在 εᵣ = 4 的介质中，电磁波速度约为真空中光速的一半，波长也相应缩短。这意味着同样物理长度的互连，在介质中对应的电长度更大，更容易出现谐振和相位延迟。
三、边界条件：不同材料交界处的约束
封装结构由导体、介质和空气等多种材料组成，场在不同材料交界面上必须满足边界条件。设 n 为从介质 1 指向介质 2 的单位法向量，切向分量为 t：
n × (E₂ - E₁) = 0 （电场切向连续）
n × (H₂ - H₁) = Jₛ （磁场切向跳变等于面电流）
n · (D₂ - D₁) = ρₛ （电位移法向跳变等于面电荷）
n · (B₂ - B₁) = 0 （磁感应法向连续）
对于理想导体（PEC）表面，E 的切向分量为零，H 的切向分量由面电流决定。对于实际铜导体，在 GHz 频段趋肤深度 δ = √(2/(ωμσ)) 约为 1～2 μm，远小于导体厚度，因此常近似为有耗导体边界或表面阻抗边界。边界条件的正确处理直接决定仿真结果的准确性，尤其是导体表面粗糙度和介质界面处的场集中效应。
四、波动方程：从方程组到可求解形式
将本构关系代入麦克斯韦方程组，消去 H 或 E，可得到只含一个未知矢量的二阶偏微分方程。在无源区域（J = 0，ρ = 0），电场满足齐次矢量波动方程：
∇²E - μ₀ε₀εᵣ ∂²E/∂t² = 0
对于时谐场（e^{jωt} 约定），令 k₀² = ω²μ₀ε₀，得到亥姆霍兹方程：
∇²E + k₀²εᵣE = 0
有源区域则右端出现激励项 -jωμ₀J。封装互连中，激励通常通过端口面上的等效电流或电压源引入。波动方程的意义在于：它把“场如何传播”转化为一个边值问题——给定几何、材料、边界和激励，求满足方程和边界条件的场分布。
五、全波求解的数学本质
全波求解的核心步骤可以概括为：将连续的偏微分方程离散为有限维矩阵方程，再求解该矩阵方程得到离散场系数。以有限元法（FEM）为例：
- 网格剖分：将封装互连的三维区域划分为四面体或六面体单元。网格尺寸通常取 λ_min/10～λ_min/5，其中 λ_min 为最高频率下介质中的最小波长。例如，在 εᵣ=4、f=10 GHz 时，λ ≈ 15 mm，网格尺寸约 1.5～3 mm；但在导体边缘和过孔附近，需局部加密到 0.1～0.3 mm。
- 基函数展开：在每个单元内用矢量基函数（如棱边元）展开 E，将连续场表示为有限个系数的线性组合。
- 加权残差：用伽辽金法对波动方程做加权积分，将微分方程转化为弱形式，降低对场连续性的要求。
- 组装矩阵方程：将所有单元的贡献组装成全局稀疏矩阵方程
Ax = b，其中 A 为刚度矩阵，x 为待求场系数，b 为激励向量。 - 求解与后处理：用直接法或迭代法求解 x，再通过端口积分提取 S 参数、特性阻抗和损耗。
FDTD 则采用时域差分，将空间和时间都离散为网格，通过蛙跳格式逐步推进场量；MoM 则基于积分方程，只离散导体表面，适合开放区域问题。三种方法的数学本质相同：把无限维的场问题转化为有限维的代数问题。
六、从方程到工程参数的映射
求解得到 E 和 H 后，工程参数通过以下方式提取：
| 工程参数 | 提取方式 | 关键条件 |
|---|---|---|
| S 参数 | 端口面上对 E×H 积分得到入射/反射/传输功率 | 端口需满足单模传输，参考面一致 |
| 特性阻抗 Z₀ | 由电压 V=∫E·dl 与电流 I=∮H·dl 的比值得到 | 需明确电压积分路径和电流回路 |
| 导体损耗 | P_loss = ∫σ|E|²dV 或表面阻抗积分 | 需正确设置导体电导率和粗糙度模型 |
| 介质损耗 | P_loss = ωε₀εᵣtanδ∫|E|²dV | 需输入频率相关的 tanδ |
理解这一映射关系，才能在仿真结果异常时判断是场求解出了问题，还是后处理设置不当。例如，S 参数在某个频点出现非物理的尖峰，可能是网格谐振或端口模式不纯；插入损耗随频率单调上升但斜率异常，可能是导体粗糙度模型或介质损耗参数设置错误。
本节建立的方程—本构—边界—波动方程—离散求解这条链条，是后续所有建模、仿真和结果分析的基础。只有清楚每一步的数学含义和物理约束，才能在全波分析中做出合理假设、设置正确参数并准确判读结果。