4. 网格剖分策略与收敛性

在毫米波频段，电路尺寸往往只有几个毫米，而金属边缘、缝隙、过孔等细节却处在几十微米量级。全波仿真器把连续场离散成有限未知量，网格就是这一离散化的载体。网格剖分策略决定了矩阵规模、求解时间，也直接决定结果是否可信。本节围绕电大与电小结构的不同处理方式，讨论自适应网格、曲边单元、网格密度控制与收敛判据，并用微带不连续性案例说明如何验证网格无关性。

电大与电小结构对网格的不同诉求

毫米波结构常同时包含电大与电小特征。以 60 GHz 微带线为例，自由空间波长约 5 mm，介质中有效波长约 3.3 mm；而导体厚度可能只有 2 μm，缝隙宽度可能只有 50 μm。若按最小特征统一加密，网格量会爆炸；若按最大尺寸统一粗化，边缘电流又无法解析。合理做法是分区对待：对电大区域（如均匀传输线段、大块介质）允许较粗网格，对电小区域（如不连续处、耦合缝隙、过孔壁）局部加密。

判断网格是否足够，不能只看单元尺寸，还要看每波长采样点数。工程上常要求介质中每波长至少 10～20 个网格，边缘处再额外加密。对于曲边结构，如圆形过孔、弧形耦合线，还需用曲边单元贴合几何边界，否则阶梯近似会引入虚假反射。

自适应网格与曲边单元

自适应网格的核心是“先粗后细、按误差加密”。典型流程是：先用较粗网格求解，估计各单元的场变化或残差，在误差大的区域细分单元，再求解，直到误差指标或 S 参数变化满足阈值。误差估计常用场梯度、电流密度跳变或后验残差。对毫米波微带不连续处，误差通常集中在导体边缘和介质界面，自适应加密会优先落在这些位置。

曲边单元用于描述弯曲边界。以圆形过孔为例，若用直边三角形逼近，过孔周长被阶梯化，等效半径和电感都会偏移。采用二阶曲边单元后，边界由二次曲线描述，几何误差从 O(h) 降到 O(h²)，在相同网格量下能显著改善过孔阻抗和回波损耗的精度。对毫米波频段，曲边单元对耦合器圆弧、滤波器圆角等结构尤其重要。

网格密度与收敛判据

网格密度不能只靠“每波长多少点”一刀切。对毫米波结构，建议同时控制三类参数：全局最大单元尺寸、介质中每波长采样点数、导体边缘局部加密层数。一个可操作的起点是：介质中每波长 15 个单元，导体边缘额外 2～3 层加密，曲边结构采用二阶曲边单元。之后通过自适应迭代或手动加密，观察 S 参数是否稳定。

收敛判据通常分三层。第一层是 S 参数稳定：相邻两次加密的 |ΔS11| 小于 0.05 dB、|ΔS21| 小于 0.1 dB。第二层是能量守恒：无耗结构应满足 |S11|²+|S21|²≈1，偏差小于 0.5%。第三层是场量收敛：导体边缘电流密度峰值变化小于 2%。三层判据同时满足，才可认为网格无关。

微带不连续性案例：网格无关性验证

下面以一个 60 GHz 微带阶梯不连续结构为例，说明验证过程。结构参数为：介质 εr=3.0，厚度 h=0.127 mm，导体厚度 t=2 μm，输入线宽 W1=0.28 mm，输出线宽 W2=0.18 mm，阶梯跳变位于结构中部。端口采用波端口，频率扫描 55～65 GHz。

验证时固定几何与材料，仅改变网格密度。第一组采用粗网格，全局最大单元 0.08 mm，介质中每波长约 8 个单元，网格量约 1.2 万。第二组加密到每波长 12 个单元，网格量约 3.5 万。第三组每波长 18 个单元，并在阶梯跳变处额外加密，网格量约 8 万。第四组每波长 24 个单元，网格量约 15 万。

表中数据为教学示例，用于说明趋势而非实测结果。可以看到，G1 到 G2 变化较大，说明粗网格尚未解析边缘电流；G3 到 G4 的 ΔS11 已小于 0.05 dB，能量偏差也降到 0.2%，可判定网格无关。若继续加密到 G5，S11 变化通常小于 0.02 dB，但网格量翻倍，求解时间显著增加，因此 G3 或 G4 是精度与效率的折中点。

实操建议与常见误区

实际项目中，建议先做一次快速粗网格扫描，确认端口、材料和几何无误，再启动自适应加密。自适应迭代次数可设为 3～5 次，每次加密比例控制在 20%～40%，避免单次加密过猛导致网格量突增。对毫米波电大结构，可先用均匀粗网格加局部加密，再对关键区域做曲边处理。

常见误区包括：只按全局尺寸加密而忽略边缘；用一阶直边单元逼近圆弧；只看 S11 不看能量守恒；以及把“网格量最大”等同于“结果最准”。网格无关性验证的本质是证明结果不再随网格变化，而不是追求无限加密。掌握这一判据，才能在毫米波全波仿真中既保证精度，又控制计算成本。