8.2 有限元法（FEM）

8.2 有限元法（FEM）

FDTD 用规则网格切空间，FEM 的回答更精细：把结构剖成任意形状的小块（单元），每块里用多项式插值逼近场，再把"让残差最小"的变分条件变成一个大型线性方程组。它对复杂几何的贴身贴合能力，使它成为 HFSS、COMSOL 等工业仿真软件的主力内核——机身、连接器、手机外壳的仿真几乎都跑在 FEM 上。

思路：从变分原理到线性方程组

以静电/亥姆霍兹问题为例。边值问题（如第 2 章的泊松方程 ∇²φ = −ρ/ε₀ + 边界条件）与一个泛函极值问题等价：

J[φ] = ∫[(1/2)|∇φ|² − ρφ/ε₀] dV → 取极小

（这不是巧合：第 2 章唯一性定理保证了两者同解；物理上泛函正是场能，"场自动选择能量极小的分布"。）

FEM 的三步走：

剖分：区域切成四面体/六面体单元，顶点为节点。三角形/四面体能贴着任意曲线边界走——这正是 FEM 几何能力的来源（对比 FDTD 的矩形台阶）。

插值：每个单元内 φ 用节点值的线性（或高阶）组合表示，φ(x) = ΣNᵢ(x)·φᵢ，Nᵢ 为形状函数（局部多项式）。整个域的未知量缩为节点值向量 {φᵢ}——连续问题离散成有限维。

组装与求解：把泛函对每个节点值求偏导并置零，得到 K·φ = f——K 为稀疏刚度矩阵（每个节点只与同单元的邻居耦合），f 为源项。解出 φ，再后处理 E = −∇φ、电容、损耗等。

这套"剖分-插值-变分"的流程对任何椭圆型方程（静电、静磁、频域亥姆霍兹）都成立；时谐电磁问题里得到的是复数方程组（解含相位），特征值问题版本给出波导的模式与谐振腔的谐振频率（第 6 章的"本征值骨架"在此被数值求解）。

与 FDTD 的性格对比

经验法则：几何复杂、频点不多 → FEM；频带很宽、时域现象重要 → FDTD；开放辐射/散射 → 下一节 MoM。三者不是竞争而是分工，商业软件里常常共存（HFSS 同时有 FEM 与积分方程求解器）。

演练：一维 FEM 手算静电场

# 一维泊松方程 d2φ/dx2 = -1（均匀源），两端 φ=0：解析解 φ=x(1-x)/2
import numpy as np
n = 5                       # 5 个内部节点，步长 h=1/6
h = 1.0/(n+1)
K = (2*np.eye(n) - np.eye(n, k=1) - np.eye(n, k=-1))/h   # 刚度矩阵（二阶差分）
f = np.ones(n)*h                                           # 源项积分（集中载荷）
phi = np.linalg.solve(K, f)
exact = np.array([(i+1)*h*(1-(i+1)*h)/2 for i in range(n)])
print("FEM 节点值:", np.round(phi, 4))
print("解析解   :", np.round(exact, 4))
# 输出两列一致到 4 位小数；加密网格 n=20 误差继续下降
​

五行核心代码就是完整 FEM 流程（剖分-插值-组装-求解的最简版）：K 的三对角结构就是"每个节点只连邻居"的稀疏性；右端 f 的 h 因子是载荷积分。把源换成 f = np.sin(np.linspace(0, np.pi, n+2)[1:-1])*h，解析解换成特解对照——换源不换框架，与第 1 章"格林函数哲学"遥相呼应。

网格收敛与自适应加密

FEM 的误差理论给出置信工具：单元尺寸 h 减半，线性元的误差约降为 1/4（二阶收敛）。工程实践靠自适应网格：先粗算一遍，找出误差指示子大的区域（场剧变处：拐角、缝隙、金属尖端——正是第 2 章"尖端 σ 大"的数值版），局部加密再算，迭代直到目标量（如谐振频率、S 参数）收敛到指定小数位。看一份仿真报告，第一件事是问"目标量随网格加密收敛了吗"——没有收敛性研究的单点结果只是"电脑画的画"。

易错点

单元质量差毁掉精度：细长扁平的四面体使插值误差爆炸，剖分不是"越密越好"而是"形状规整+关键区加密"；商业软件的网格检查（雅可比比、纵横比）必须过。

边界条件选错类型：狄利克雷（给定值）与诺依曼（给定导数）混用位置错了，解可能仍收敛但完全错误（第 2 章唯一性定理要求两类条件恰当地给在边界上）。

频域单点的"窄带错觉"：FEM 单频求解的成本结构诱使用户只算几个频点，漏掉窄带谐振；扫频策略（自适应抽取频点插值）要覆盖 S 参数的共振细节，第 6 章谐振腔的高 Q 特性对采样密度有硬要求。

补充辨析：三道进阶是非题

第一道："FEM 只能算静场与低频"——过时：频域 FEM 配合 PML 与积分方程边界，已覆盖到毫米波与光频的工业问题；限制来自问题的电尺寸而非频率标签。第二道："自适应网格收敛了，结果就可信"——只对了一半：自适应迭代消除的是离散误差，几何简化（倒角忽略、焊缝理想化）、材料参数误差（ε_r 随温漂）、边界简化（腔体开孔漏波）仍要靠工程师的判断与实验对标。第三道："高阶单元一定比线性单元好"——在光滑解区域是；但在拐角、尖端等场奇点附近，高阶多项式的优势消失，局部加密的线性元反而稳健；"单元阶数与网格密度的组合"才是真决策变量。三道题训练的是"误差预算"思维：把总误差拆成离散、模型、参数、测量四栏分别管理——这份预算表意识，是从"会用软件"到"会背书结果"的分水岭，也是第 9 章读仿真报告时反复要用的看家本领。

数字感受一组：连接器仿真典型 50 万未知量（稀疏矩阵，分钟级）、手机整机天线 100–300 万（小时级）、飞机级 RCS 全波需千万级（必须 MLFMA 或混合方法，天级）——三个量级对应三类商业软件的授权价位。理解算力账本，就理解了仿真部门的排期与预算为什么长成那个样子。