3. HFSS求解器与网格剖分原理
在高速互连的三维电磁仿真中，HFSS 的核心并不是“把模型画出来然后点求解”，而是把连续的麦克斯韦方程组离散成有限维代数方程组，再用迭代方式逼近真实场分布。理解这一过程，才能判断结果是否可信、网格是否足够、求解设置是否合理。本节围绕有限元法（FEM）、自适应网格、收敛判据与求解设置四条主线展开，重点说明它们之间的因果关系与工程判断依据。
一、有限元法：从连续场到离散矩阵
HFSS 采用基于四面体单元的有限元法。其基本思想是：将求解区域剖分为大量小四面体，在每个单元内用多项式基函数近似电场 E，再通过伽辽金加权残量法把矢量波动方程转化为稀疏线性方程组：
[K]{E} = {b}
K：刚度矩阵（含材料、边界、波数信息）
E：待求电场自由度向量
b：激励与边界条件右端项HFSS 默认使用一阶或二阶四面体棱边元（edge element），未知量定义在棱边上而非节点上，这样能自然满足切向连续、抑制伪解。单元尺寸越小、阶数越高，未知量越多，矩阵越大，但场分辨率越高。对高速互连而言，关键区域是导体边缘、过孔反焊盘、参考平面开槽、连接器引脚附近——这些位置场梯度大，需要更密网格；而均匀传输线段可以适当放宽。
二、自适应网格：误差驱动的迭代加密
HFSS 不要求用户一次性给出完美网格，而是采用自适应迭代：先给一个较粗的初始网格求解，再根据场解误差估计，在误差大的单元处加密，重新求解，直到满足收敛判据。其流程可概括为“求解—估计—加密—再求解”的闭环。
自适应网格的关键在于误差估计器。HFSS 常用基于 S 参数变化量 ΔS 的判据：每次迭代后比较当前与上一次的 S 参数，若所有频点、所有端口的 |ΔS| 最大值小于设定阈值，则认为收敛。也可结合场量变化或能量误差。加密时并非全局均匀加密，而是按误差密度选择性细分，因此效率远高于一次性全局细网格。
三、收敛判据：判断“算够了没有”
收敛判据是自适应迭代的停止条件，直接决定结果可信度与计算成本。HFSS 中常用设置包括：
- Maximum Delta S：相邻两次迭代 S 参数最大变化量，默认 0.02，高速互连建议 0.005～0.01。
- Maximum Number of Passes：最大迭代次数，默认 6，复杂结构可设 10～15。
- Minimum Number of Passes：最少迭代次数，防止偶然收敛，通常设 2。
- Delta S per Pass 曲线：在求解信息中查看收敛趋势，若持续下降后趋于平缓，说明收敛正常。
工程判断上，不能只看“是否达到阈值”。若 ΔS 在阈值附近震荡、或某频点始终不收敛，往往说明网格在关键区域仍不足、端口设置不当、或存在尖锐几何。此时应检查端口尺寸、材料边界、是否遗漏参考平面，而不是盲目增加迭代次数。
四、求解设置：从频域到工程结果
HFSS 默认采用频域求解，适合高速互连的宽带 S 参数提取。关键设置包括：
| 设置项 | 典型值/建议 | 影响 |
|---|---|---|
| 求解频率 | 中心频率或最高关注频率 | 决定初始网格密度与自适应加密基准 |
| 扫频方式 | 插值扫频 / 快速扫频 | 插值适合窄带，快速适合宽带；点数影响精度 |
| 频率范围 | 如 0.1～20 GHz | 覆盖信号带宽与谐波关注范围 |
| 端口激励 | 波端口 / 集总端口 | 波端口适合传输线，集总端口适合芯片引脚 |
| 最大迭代次数 | 6～15 | 不足则未收敛，过多则耗时 |
| Delta S 阈值 | 0.005～0.02 | 越小越准，计算量越大 |
综上，HFSS 的求解过程是“有限元离散—自适应加密—收敛判据—求解设置”的闭环。理解每一步的物理与数值含义，才能在高速互连仿真中既保证精度，又控制计算成本。实际工程中，建议把收敛判据、端口设置与网格统计作为每次仿真的必查项，形成可复现的求解配置习惯。