3、网格剖分原理：有限元法（FEM）与时域有限差分法（FDTD）中网格的作用，网格类型（四面体、六面体、三角形）及其适用场景
天线仿真中，网格剖分是将连续几何模型离散为有限数量单元的过程，其质量直接决定求解精度、内存开销与计算时长。有限元法（FEM）和时域有限差分法（FDTD）作为两类主流全波算法，对网格的依赖方式截然不同：FEM基于非结构化网格适应复杂曲面，FDTD基于结构化网格依赖空间步长稳定性。理解网格在两种算法中的角色差异，以及四面体、六面体、三角形等单元类型的几何特性与适用边界，是进行扫参分析前必须夯实的理论基础。本章将从算法原理出发，剖析网格作用机制，并给出工程选型建议。
3.1 核心原理：FEM 与 FDTD 对网格的依赖机制
FEM 的网格作用：有限元法从变分原理出发，将 Maxwell 方程组转化为泛函极值问题，在离散单元内使用基函数（如 Whitney 矢量基函数）近似电磁场分布。网格单元是基函数的支撑域，单元形状与尺寸决定了插值精度。FEM 采用非结构化网格，四面体可以贴附任意弯曲表面，因此对天线辐射体、介质分层、馈电结构等复杂几何具有天然适应性。求解过程最终形成稀疏矩阵方程，网格数量直接决定矩阵维度，影响内存与求解时间。
FDTD 的网格作用：时域有限差分法直接在时空域对 Maxwell 旋度方程进行中心差分近似，采用 Yee 元胞（正交六面体）将空间离散为规则网格。电场与磁场分量在元胞棱边和面中心交替分布，时间步进推进。FDTD 的网格必须满足 Courant 稳定性条件（Δt ≤ 1/(c√(1/Δx²+1/Δy²+1/Δz²))），且网格尺寸通常取最高工作频率波长的 1/10~1/20。由于采用正交结构化网格，曲面结构只能通过阶梯近似描述，需引入共形网格技术（如 CFDTD）改善精度。
3.2 网格类型对比与适用场景
| 网格类型 | 维度 | 典型算法 | 优势 | 局限 | 适用场景 |
|---|---|---|---|---|---|
| 四面体 | 3D 体 | FEM | 任意几何自适应、局部加密灵活 | 数值色散较大、内存开销高 | 复杂曲面天线、共形天线、介质谐振天线 |
| 六面体 | 3D 体 | FDTD / FEM | 网格方向性精度高、计算效率高 | 曲面拟合困难、生成算法复杂 | 微带天线、阵列天线、波导结构、平面周期结构 |
| 三角形 | 2D 表面 | FEM 边界元 / 矩量法 | 表面建模灵活、与体网格衔接好 | 仅覆盖表面，需配合体网格 | 贴片天线表面、PCB 走线、孔径耦合面 |
| 棱柱（三棱柱） | 2.5D 拉伸 | FEM | 薄层结构效率高、过渡连接好 | 仅适用于拉伸几何 | 介质基板、薄涂层、多层 PCB |
四面体网格的工程选择：当仿真对象包含曲率半径小于波长的结构（如螺旋天线、锥形喇叭、介质透镜）时，四面体网格能以较少单元数量获得较高几何拟合度。但需注意，四面体的数值色散各向异性较强，扫参时若频率范围较宽，需在最高频点验证网格收敛性。
六面体网格的工程选择：FDTD 中六面体网格是唯一选择，其正交性保证了差分格式的简洁性。对于平面型天线（如贴片天线、Vivaldi 天线），六面体网格可以沿结构主方向对齐，减少阶梯误差。FEM 中也可使用六面体网格，但生成高质量全六面体网格对几何要求苛刻，通常仅用于规则波导或腔体结构。
三角形表面网格的衔接作用：在 FEM 中，三角形网格常作为四面体网格的表面边界条件载体，用于定义端口、辐射边界或阻抗边界。在多层结构中，三角形表面网格与棱柱网格配合，可有效处理薄介质层。
3.3 案例拆解：贴片天线扫参中的网格策略对比
背景与约束：某 2.4 GHz 矩形微带贴片天线，基板厚度 1.6 mm，介电常数 4.4，需扫参分析贴片长度 L 从 28 mm 至 32 mm 变化时谐振频率的偏移规律。仿真工具分别采用 FEM 求解器（四面体网格）与 FDTD 求解器（六面体网格）。
初始现象与目标：FEM 默认自适应网格在 L=30 mm 时谐振频率为 2.41 GHz，FDTD 默认网格（Δ=1 mm）结果为 2.38 GHz，两者偏差 1.3%。目标是通过网格收敛性验证，确定扫参结果的可靠性。
关键输入数据与操作步骤：
- FEM 网格收敛验证：设置四面体最大边长分别为 λ/8、λ/12、λ/16（λ 为 2.4 GHz 自由空间波长 125 mm），对比谐振频率变化。结果：λ/12 与 λ/16 的谐振频率差小于 0.2%，判定 λ/12 为收敛网格。
- FDTD 网格收敛验证：设置网格步长 Δ=2 mm、1 mm、0.5 mm，观察谐振频率变化。结果：Δ=1 mm 与 Δ=0.5 mm 的频率差为 0.4%，但 Δ=1 mm 时阶梯误差导致基板厚度方向仅 1~2 个网格，需采用共形网格技术修正。
- 扫参执行：FEM 使用 λ/12 网格，每个 L 值重新剖分；FDTD 使用 Δ=1 mm 固定网格，仅修改贴片边界位置。
结果对比与复盘：FEM 扫参曲线平滑，谐振频率随 L 增大线性下降（约 0.8 MHz/mm）；FDTD 曲线出现轻微台阶效应，尤其在 L=30 mm 附近偏差 0.5%。原因在于 FDTD 阶梯近似使贴片有效电长度随网格对齐状态跳变。最终以 FEM 结果为准，FDTD 用于验证趋势一致性。
3.4 常见误区与总结复盘
误区一：网格越密精度越高。实际上，当网格密度达到一定程度后，精度提升趋于饱和，而内存与时间开销急剧增加。更关键的是，过密网格可能引入数值噪声，尤其在 FDTD 中，网格过密导致时间步长减小，累积舍入误差上升。
误区二：忽略网格质量指标。四面体网格的雅可比行列式、长宽比、扭曲度直接影响 FEM 矩阵条件数。扫参时若几何变化导致网格质量恶化，求解可能不收敛或产生伪解。建议每次剖分后检查最小网格质量，低于 0.2 时需调整剖分参数。
误区三：FDTD 中直接使用非均匀网格处理曲面。非均匀网格虽然能局部加密，但会引入额外的数值反射。对于曲面结构，应优先使用共形网格技术（如 Yu-Mittra 共形方案），而非简单加密。
误区四：扫参时网格策略不一致。FEM 中若不同参数点采用不同自适应迭代次数，会导致网格密度不一致，使扫参曲线出现虚假波动。应固定自适应迭代次数或固定目标网格数量。
总结复盘：网格剖分是天线仿真中“一分精度、一分代价”的典型环节。FEM 与 FDTD 的网格策略差异源于算法本质：FEM 以几何拟合换取矩阵求解的灵活性，FDTD 以规则网格换取时间步进的简洁性。工程实践中，应根据天线结构复杂度、扫参范围、计算资源综合选型。对于多参数扫参，建议先进行网格收敛性验证，确定基准网格密度，再执行批量仿真，以确保扫参结果的可比性与可信度。