3. 电磁仿真基础与封装应用

封装电磁仿真的核心任务，是在封装几何与材料分布已知的前提下，求解结构内部及周围的电场与磁场分布，并据此提取电阻、电感、电容、阻抗、S参数、辐射方向图等工程参数。要理解仿真软件在做什么，必须先回到麦克斯韦方程组这一物理源头，再理解高频下导体内部电流如何重新分布，最后把这两件事落到封装互连与天线两类典型场景上。

一、麦克斯韦方程组：仿真器真正在解的方程

麦克斯韦方程组由四个方程构成，分别描述电场与电荷、磁场与电流、时变电场与磁场之间的耦合关系。在封装电磁仿真常用的频域形式下，可写为：

∇ × E = -jωB          （法拉第定律：时变磁场产生电场）
∇ × H = J + jωD       （安培定律：传导电流与位移电流共同产生磁场）
∇ · D = ρ             （高斯定律：电荷是电场的源）
∇ · B = 0             （磁通连续性：不存在磁单极子）

其中 E 为电场强度，H 为磁场强度，B 为磁感应强度，D 为电位移矢量，J 为传导电流密度，ρ 为自由电荷体密度，ω 为角频率。对封装仿真而言，最关键的是前两个旋度方程：它们决定了电磁能量如何在互连、介质和空间中传播与耦合。

实际仿真不会直接求解上述微分形式，而是根据问题类型选择不同的数值方法。常见的有有限元法（FEM）、时域有限差分法（FDTD）、矩量法（MoM）和边界元法（BEM）。选择依据主要看三点：结构是封闭还是开放、是否需要宽频带扫描、以及电尺寸大小。例如封装内电源分配网络（PDN）的阻抗提取常用 FEM 或 MoM，而天线辐射问题更常用 FDTD 或 MoM。

二、涡流与趋肤效应：高频电流为何只走表面

当导体处于时变磁场中时，导体内部会感应出闭合的涡旋电流，即涡流。涡流的方向总是阻碍原磁通的变化，这既是电感在高频下损耗增加的来源，也是电磁屏蔽和涡流制动的工作基础。在封装中，涡流主要出现在焊球、铜柱、引线框架和金属散热盖等大块导体中。

趋肤效应则是涡流在导体自身电流上的体现。交变电流产生的磁场在导体内部感应出反向涡流，使电流向导体表面集中。电流密度从表面向内部按指数衰减，衰减到表面值 1/e 处的深度称为趋肤深度 δ：

δ = sqrt( 2 / (ω μ σ) ) = sqrt( 1 / (π f μ σ) )

以铜为例，σ ≈ 5.8×10⁷ S/m，μ ≈ μ₀ = 4π×10⁻⁷ H/m。代入可得：

这张表说明一个关键判断：当铜线厚度大于趋肤深度时，电流只在一层薄壳内流动，有效截面积减小，交流电阻显著高于直流电阻。仿真中若仍按均匀电流分布计算，会低估高频损耗。因此高频封装互连建模时，必须让网格在导体表面附近足够密，或采用表面阻抗边界条件来等效趋肤效应，否则结果会偏离实际。

三、封装互连场景：从直流电阻到高频阻抗

封装互连包括焊球、铜柱、再分布层（RDL）、基板走线和键合线。低频时，互连主要表现为电阻；频率升高后，电感、电容和趋肤效应逐渐主导。一个典型的封装互连电磁仿真流程如下：

1. 几何简化：保留影响电磁场分布的关键结构，如焊球阵列、参考平面、介质层，去除对结果影响很小的圆角、文字标记和测试焊盘。
2. 材料赋值：铜导体设 σ = 5.8×10⁷ S/m，焊料设 σ ≈ 7×10⁶ S/m，介质按频变模型输入 εr 和 tanδ。
3. 端口设置：在互连两端设置波端口或集总端口，端口阻抗通常设为 50 Ω，并明确参考导体。
4. 频率扫描：例如从 10 MHz 扫到 20 GHz，观察插入损耗 S21 和回波损耗 S11 随频率的变化。
5. 结果提取：由 S参数换算得到等效 R、L、C、G，或直接查看阻抗曲线 Z(f)。

以一个简化示例说明判断依据：某封装走线在 1 GHz 时仿真得到 S21 = -0.8 dB，S11 = -18 dB；在 10 GHz 时 S21 = -3.5 dB，S11 = -8 dB。这说明高频下反射和损耗同时恶化，可能原因是走线特征阻抗偏离 50 Ω，或参考平面不连续导致回流路径变长。此时应检查走线宽度、介质厚度和参考平面开槽位置，而不是盲目增加导体厚度。

四、天线场景：封装天线与电磁辐射

封装天线（AiP）将天线集成在封装基板或芯片上方，工作频率通常在 24 GHz、60 GHz 或 77 GHz 等毫米波频段。此时封装尺寸与波长可比，电磁场不再局限于导体内部，而是向空间辐射。仿真目标也从阻抗提取转变为方向图、增益、效率和波束宽度。

天线仿真与互连仿真的关键区别在于边界条件。互连仿真常用理想匹配层（PML）或辐射边界来吸收向外传播的波，避免反射污染结果；天线仿真则必须使用开放边界，并保证天线到边界距离大于四分之一波长。以 60 GHz 为例，自由空间波长约 5 mm，封装天线到边界距离应大于 1.25 mm，实际工程中常取 3～5 mm 以留出余量。

天线场景中还需要关注以下参数：

• 回波损耗 S11：通常要求在工作频段内低于 -10 dB，表示 90% 以上功率被辐射。
• 辐射效率：考虑导体损耗和介质损耗后的实际辐射功率占比，毫米波封装天线常要求高于 70%。
• 增益与方向图：用于判断波束是否覆盖目标区域，以及旁瓣是否过高。
• 封装材料影响：塑封料和基板介质的 εr 与 tanδ 会改变谐振频率和效率，仿真中必须使用频变实测参数。

一个常见的教学示例是 60 GHz 贴片天线阵列：4×4 阵列，单元间距约 2.5 mm（半波长），基板 εr = 3.0，厚度 0.2 mm。仿真中若发现谐振频率偏移到 58 GHz，首先检查介质 εr 是否设置偏高，其次检查馈线长度和地平面尺寸是否引入额外寄生电容。调整馈线长度或增加匹配枝节后，S11 应回到 -15 dB 以下。

五、把方程、效应与场景串起来

麦克斯韦方程组给出物理约束，涡流与趋肤效应决定导体内部电流如何分布，封装互连和天线则是这些规律在具体结构上的体现。仿真时，先判断工作频率与结构电尺寸，再决定用哪种数值方法和网格密度；然后检查材料参数是否频变、端口和边界是否合理；最后用 S参数、阻抗、方向图等指标判断结果是否可信。只有把物理判断和仿真设置对应起来，才能避免“算得出结果、却解释不了结果”的常见问题。