4、网格基础：网格类型（六面体、四面体、薄片）、网格剖分原理、网格质量指标

CST 的时域求解器（T 求解器）与频域求解器（F 求解器）均依赖网格离散 Maxwell 方程组。网格不仅是几何模型的数值近似，更直接决定了矩阵条件数、局部截断误差与时间步长稳定性。在“不收敛”的排障流程中，网格问题往往排在端口与边界之前——因为错误的网格会让端口激励无法正确耦合，也会让边界条件在非物理位置生效。本章聚焦三类核心网格：六面体、四面体、薄片网格，剖析其剖分原理与质量判定方法，帮助你在仿真报错或结果发散时快速定位网格根因。

4.1 三类网格的核心原理与适用场景

CST 中网格类型的选择并非随意为之，而是由求解器、电尺寸、几何复杂度共同决定。六面体网格是时域求解器的默认选择，其核心优势在于正交结构化拓扑，使得 FDTD（时域有限差分）更新方程可以高效迭代。六面体网格采用 PBA（Perfect Boundary Approximation）技术，允许网格边界部分填充材料，从而在不细化网格的情况下逼近曲面边界。但 PBA 的近似在尖锐边缘或薄介质层处会产生等效介电常数误差，这是导致谐振频率偏移的常见原因。

四面体网格主要用于频域求解器（F 求解器），其非结构化特性可以精确贴合任意复杂曲面，尤其适合圆角、倒角、螺旋结构。四面体网格的剖分基于 Delaunay 三角化算法，通过迭代优化点位置使最小角最大化。然而，四面体网格的缺点是低频时矩阵规模庞大，且对薄片结构（如 PCB 走线、屏蔽罩）会产生大量高长宽比单元，严重恶化矩阵条件数。

薄片网格（Thin Sheet）是 CST 针对电小厚度结构的专用处理方式。当结构厚度远小于工作波长（通常小于 λ/20）时，若仍用六面体或四面体剖分，会产生大量细长单元，导致时间步长急剧减小（时域）或矩阵病态（频域）。薄片网格将厚度方向压缩为单层等效表面，通过表面阻抗或传输条件建模，大幅减少网格数量。但薄片网格牺牲了厚度方向的场分布细节，对于需要精确计算趋肤效应或介质损耗的场景不适用。

4.2 网格剖分原理与关键参数控制

六面体网格的剖分过程分为三个步骤：首先根据模型外包围盒确定全局网格密度，然后通过“线”和“面”的局部加密设置细化关键区域，最后执行 PBA 修正。全局网格密度由“Cells per wavelength”参数控制，默认值为 10，即每个波长至少 10 个网格单元。对于高 Q 值结构或需要精确计算表面电流的模型，建议提升至 15-20。但盲目增加全局密度会线性增加内存消耗，并让时域时间步长按比例缩小，导致仿真时间立方级增长。

局部加密是解决“全局过密、局部不足”矛盾的核心手段。CST 提供三种局部加密方式：固定点数加密（在指定边或面上设置网格数）、曲率自适应加密（根据几何曲率自动加密）、场自适应加密（根据初步仿真结果加密高场强区域）。在排障实践中，最常用的是固定点数加密——例如在微带线馈电点附近设置 5-10 个网格，在缝隙天线槽口设置 3-5 个网格。需要特别注意的是，局部加密区域与周围网格的过渡必须平滑，否则会产生虚假反射。

四面体网格的剖分控制参数包括“Max. edge length”（最大边长）、“Min. edge length”（最小边长）和“Curvature refinement”（曲率细化）。最大边长通常设置为 λ/10 至 λ/15，最小边长则受限于几何最小特征尺寸。曲率细化系数控制曲面上的网格密度，默认值为 1.0，对于圆波导或球面结构建议提升至 2.0。四面体网格的剖分质量与几何模型的“清洁度”高度相关——导入的 STEP 文件中若存在微小碎面、重叠面或缝隙，会导致剖分器生成大量畸形单元。

// 六面体网格关键参数示例（CST 宏命令）
With Mesh
    .SetGlobalMeshType "Hexahedral"
    .SetCellsPerWavelength 15
    .SetLowerEdgeLimit 0.001
    .SetUpperEdgeLimit 0.5
    .SetLocalEdgeFactor "feed_line", 5
    .SetLocalEdgeFactor "slot_gap", 3
End With

4.3 网格质量指标与不收敛的关联

网格质量直接决定矩阵条件数与数值色散误差。对于六面体网格，最重要的质量指标是“长宽比”（Aspect Ratio），即单元最长边与最短边的比值。理想值为 1，工程上允许小于 5。当长宽比超过 10 时，时域求解器的 Courant 稳定性条件会强制时间步长降至最小边对应的值，导致仿真时间爆炸。更严重的是，高长宽比单元在 PBA 修正时会产生等效介电常数误差，引发谐振频率偏移。

四面体网格的核心质量指标是“偏斜度”（Skewness），定义为单元内切球半径与外接球半径之比。偏斜度越接近 1，单元越接近正四面体。CST 中偏斜度低于 0.2 的单元被视为“退化单元”，会导致频域矩阵奇异。另一个重要指标是“体积比”（Volume Ratio），即最大单元与最小单元的体积比。在包含微带线（宽度 0.1mm）与接地板（尺寸 100mm）的模型中，体积比可能达到 10⁶，此时频域求解器的迭代收敛速度会急剧下降。

薄片网格的质量评估相对特殊，主要关注“表面网格分辨率”与“厚度方向采样率”。薄片网格的默认厚度采样为 1 层，对于需要考虑趋肤效应的导体，应设置为 2-3 层。但每增加一层，网格数量呈线性增长，且厚度方向的网格尺寸必须远小于趋肤深度，否则会引入非物理的损耗。

4.4 案例拆解：微带滤波器不收敛的网格诊断

背景与约束：某 5.8GHz 微带带通滤波器，基板为 Rogers RT5880（εr=2.2，厚度 0.508mm），微带线宽 1.5mm，滤波器尺寸 30mm×20mm。使用 CST 时域求解器，初始网格采用默认设置（10 cells/λ），仿真结果 S21 在 5.5-6.0GHz 范围内出现剧烈振荡，且能量守恒曲线在 0.5ns 后偏离 1。

初步现象：S21 曲线在通带内出现 ±3dB 的波纹，且网格剖分报告显示最大长宽比为 18.7，位于微带线与接地板之间的介质区域。进一步检查发现，基板厚度方向仅划分了 2 层网格，而微带线宽度方向仅 4 个网格。

分析判断：长宽比 18.7 意味着时间步长被介质层最小网格限制，导致仿真时间步长仅为理论值的 1/5。更关键的是，微带线宽度方向的 4 个网格无法准确表征横向电流分布，导致特征阻抗计算偏差约 8%，进而引发谐振频率偏移和带内波纹。

操作步骤：① 将全局网格密度提升至 15 cells/λ；② 在微带线宽度方向设置局部加密为 8 个网格；③ 在基板厚度方向设置 4 层网格；④ 对滤波器耦合缝隙区域设置曲率自适应加密（系数 1.5）；⑤ 重新剖分后，最大长宽比降至 6.2，仿真时间步长提升 2.3 倍。

结果对比：优化后 S21 波纹降至 ±0.5dB 以内，能量守恒曲线在 0.3ns 后稳定在 0.999 以上。仿真时间从 45 分钟降至 28 分钟（得益于时间步长增大），内存占用增加 35%。

复盘经验：微带类结构的网格瓶颈往往不在全局密度，而在“薄介质层 + 窄导体”的局部区域。排障时应优先检查长宽比报告，而非盲目增加全局网格。此外，PBA 修正对微带线边缘的等效介电常数影响显著，建议在关键传输线区域使用“Edge refinement”而非单纯增加网格数。

4.5 常见误区与总结复盘

误区一：网格越密越精确。实际上，当网格密度超过一定阈值后，精度提升趋于饱和，而数值色散误差可能因累积效应反而增大。对于时域求解器，网格密度超过 20 cells/λ 后，继续加密只会增加计算负担，不会显著改善 S 参数精度。

误区二：忽略几何清理直接剖分。导入的 CAD 模型常包含微小倒角、短边、重叠面，这些特征会迫使剖分器生成大量小尺寸单元，导致长宽比恶化。建议在建模阶段使用“Heal”功能修复几何缺陷，或将不影响电磁性能的微小特征（如小于 λ/100 的圆角）删除。

误区三：薄片网格可以替代所有薄结构。薄片网格假设厚度方向场均匀，对于介质谐振器、铁氧体器件等厚度方向存在驻波的结构，薄片网格会完全错误。此时应使用六面体或四面体网格，并确保厚度方向至少 3 层网格。

误区四：只关注网格数量不关注网格分布。两个模型网格数量相同，但一个均匀分布、一个在关键区域集中加密，仿真精度可能相差数倍。网格剖分的核心原则是“在电场变化剧烈的地方加密，在均匀区域稀疏”。