4. 网格划分与自适应求解
在 HFSS 中,网格不是“画完模型后顺手点一下”的附属步骤,而是决定有限元方程组规模、求解时间与 S 参数精度的核心变量。HFSS 采用四面体网格对三维场域离散,每个四面体内部用一阶或二阶基函数逼近电场,网格越密,未知量越多,场分布刻画越细,但矩阵规模与求解时间也随之上升。理解网格划分,本质上是理解“用多少未知量换取多少精度”的工程折中。
一、HFSS 的网格控制类型
HFSS 提供多种网格操作,它们作用的对象和尺度不同,需要区分使用。
基于频率的网格(Frequency Based Mesh)是全局性控制。HFSS 在自适应迭代前会先按求解频率生成初始网格,其疏密由“每波长网格数”决定。默认情况下,HFSS 在介质中按约每波长 10~20 个网格单元布置,在金属表面附近自动加密。求解频率越高,波长越短,初始网格越密。因此,如果扫描范围是 8~12 GHz,而求解频率设为 10 GHz,那么 12 GHz 处的网格相对偏疏,高频端 S 参数误差通常更大。
长度网格(Length Based Mesh)用于对特定对象或区域强制指定最大单元边长。例如微带线导体宽度仅 0.7 mm,若全局网格边长达到 2 mm,导体边缘的电流分布就无法分辨。此时可对导体对象设置 Maximum Length 为 0.1~0.2 mm,强制在导体附近生成细网格。
表面网格(Surface Mesh)作用于对象表面,控制表面三角形的大小与曲率逼近。对波导圆角、同轴接头、弧形谐振腔等曲面结构,表面网格过粗会导致几何逼近误差,进而使端口阻抗和 S 参数出现系统性偏差。HFSS 的 Surface Deviation 与 Normal Deviation 参数即用于控制曲面离散精度。
自适应迭代(Adaptive Iteration)是 HFSS 自动加密网格的机制。每次迭代求解后,HFSS 计算各单元内的电场解误差,在误差最大的区域插入新网格点,重新求解,直到相邻两次迭代的 S 参数变化小于设定收敛判据(Delta S),或达到最大迭代次数。收敛判据通常设为 0.02,即 S 参数变化小于 0.02 时停止。
二、微带谐振器:长度网格与导体边缘加密
以一个 10 GHz 微带谐振器为例。介质基板 Rogers RO4003C,相对介电常数 3.55,厚度 0.508 mm,导体宽度约 1.1 mm,谐振长度约 5.2 mm。若仅依赖基于频率的全局网格,10 GHz 下介质中波长约 15.9 mm,按每波长 15 个单元计算,平均边长约 1.06 mm,与导体宽度相当。此时导体边缘的电流集中区域无法被准确捕捉,谐振频率可能偏移 1%~3%,回波损耗深度也会偏浅。
实际操作中,应对微带导体对象添加 Length Based Mesh,Maximum Length 设为 0.15 mm,约为导体宽度的 1/7。同时,对导体与介质交界面可启用表面网格,确保边缘电流路径被至少 3~4 个单元覆盖。自适应迭代通常需要 6~10 次才能将 Delta S 压到 0.02 以下。若迭代 4 次就“收敛”,往往说明初始网格过粗,收敛到的是局部粗糙解,需要检查是否因网格增长率过大而跳过细节。
三、波导器件:表面网格与曲面逼近
波导器件的几何特征与微带不同。以 X 波段矩形波导(WR-90,22.86 mm × 10.16 mm)为例,内部为规则矩形,但法兰连接处、圆角过渡、螺钉孔等曲面结构需要表面网格控制。若 Surface Deviation 设为 0.5 mm,而圆角半径仅 1 mm,曲面被离散为少数几个平面,端口阻抗计算会出现明显误差。
建议将 Surface Deviation 设为最小曲率半径的 1/10 以下,例如圆角半径 1 mm 时设为 0.1 mm。Normal Deviation 可设为 15°~30°,控制相邻表面法向变化。对于波导内部,基于频率的网格通常足够,因为场分布相对平滑;但波导窄边附近的电流集中区域仍建议用长度网格限制最大边长不超过窄边长度的 1/5,即约 2 mm。
四、收敛控制与结果判读
自适应迭代的收敛判据 Delta S 是判断网格是否足够的关键指标。Delta S 定义为相邻两次迭代中所有 S 参数幅度的最大变化量。设为 0.02 意味着 S 参数变化小于 0.02 时认为收敛。对于回波损耗要求较高的滤波器或谐振器,建议收紧到 0.005~0.01;对于一般传输结构,0.02 通常足够。
最大迭代次数建议设为 15~20,避免因振荡不收敛而无限循环。网格增长率(Refinement Per Pass)默认 30%,即每次迭代新增约 30% 的网格单元。若增长率过大,可能跳过误差局部极值;若过小,迭代次数增加但收敛更稳。对于微带谐振器,建议增长率设为 20%~25%。
五、网格密度对 S 参数精度的量化影响
以下为示例性对比,用于说明趋势而非实测数据。对同一微带谐振器,分别设置三种网格方案:
| 方案 | 导体最大边长 | Delta S | 迭代次数 | 谐振频率偏移 | 回波损耗深度 |
|---|---|---|---|---|---|
| A(粗) | 0.5 mm | 0.02 | 4 | +2.1% | -18 dB |
| B(中) | 0.15 mm | 0.01 | 8 | +0.4% | -32 dB |
| C(细) | 0.08 mm | 0.005 | 12 | +0.1% | -38 dB |
方案 A 虽然 Delta S 达到 0.02,但谐振频率偏移超过 2%,回波损耗深度仅 -18 dB,说明收敛判据本身不能替代网格合理性。方案 B 将导体最大边长降至 0.15 mm,频率偏移缩小到 0.4%,回波损耗深度改善到 -32 dB。方案 C 进一步加密,精度提升但收益递减,求解时间可能增加 3~5 倍。工程上通常选择方案 B 作为平衡点。
对波导器件,表面网格的影响更为突出。若 Surface Deviation 从 0.5 mm 降至 0.1 mm,端口驻波比计算误差可从 3% 降至 0.5% 以内,而网格量仅增加约 40%,因为曲面区域占比有限。这说明针对关键区域加密,比全局均匀加密更高效。
六、实操建议与检查清单
- 求解频率应设在扫描范围的中心或最高频,若关注高频端精度,可将求解频率设为扫描上限。
- 微带导体、缝隙、耦合区域必须添加长度网格,最大边长不超过特征尺寸的 1/5~1/10。
- 曲面结构检查 Surface Deviation 是否小于最小曲率半径的 1/10。
- 自适应迭代至少观察 6 次以上,Delta S 应单调下降或稳定收敛。
- 收敛后检查网格统计:单元总数、最大边长、最小边长,确认无异常尖角或薄片单元。
- 对关键 S 参数(如回波损耗深度、谐振频率)做网格无关性验证:加密一次后结果变化应小于 0.5%。
网格划分与自适应求解不是一次设置就结束的步骤,而是需要结合器件类型、频段、精度要求和计算资源反复调整的过程。微带谐振器重在导体边缘的长度网格,波导器件重在曲面的表面网格,而自适应迭代的收敛判据和增长率则决定了两者能否在合理时间内得到可信的 S 参数。掌握这些控制手段,才能在 HFSS 中把“算得出来”变成“算得准”。