6.3 有限元法与变分原理

6.3 有限元法与变分原理

矩量法统治开放域，但遇到填充复杂介质的谐振腔、含各向异性材料的波导过渡、人体组织这样的非均匀场景，积分方程开始吃力。本节的有限元法（FEM）换一条哲学路线：不去逐点满足微分方程，而是宣告"真实场分布使能量泛函取极小"——把物理问题变成能量最小化问题，再用四面体把任意几何切成碎块，每块内用简单函数逼近场。学完本节，你将理解变分原理为何天然兼容复杂材料、自适应加密为何总能把算力花在刀刃上，以及 FEM 与 MoM 如何在混合框架里交接班。下一节收尾本章：高频渐近方法与工程软件把值班室的工具箱补齐。

一、变分原理：把方程换成能量

电磁场的微分方程组（1.1 节）与一个能量积分的极值条件等价：对谐变场，泛函 F(E) = ½∫(1/μ|∇×E|² − ω²ε|E|²)dV 在满足边界条件的所有可能场中，真实解使 F 取驻定值。这层"最小能量"外衣带来三个礼物：其一，复杂材料只是改变被积函数——各向异性、非均匀、有耗介质统统不必特殊处理；其二，边界条件分两类，强加（金属面场值固定）与自然（对称面、端口自动满足），开放边界的处理反而干净；其三，离散化的误差有能量范数的全局控制，收敛性有理论保底。

剖分离散后（四面体是三维标配，任意几何都能贴合），每个单元内的场用低阶多项式形函数展开，泛函极值条件化为一个大型稀疏矩阵方程 K·e = b。稀疏是 FEM 的财务优势：矩阵非零元只出现在几何相邻单元之间，内存与求解量随未知量近似线性增长，千万级未知量在普通工作站上可解——这与 MoM 的稠密账单形成鲜明对照。

二、自适应加密与谐振问题

FEM 最聪明的机制是自适应网格加密：先粗网格跑一遍，估计每单元误差，只在高梯度区（缝隙、尖角、介质突变处）局部加密，循环直到误差达标。算力精准流向难啃的区域，而不是均匀浪费在平静地带。谐振与传输问题的求解各有讲究：本征模问题（找腔体谐振频率与场型）归结为广义特征值问题；端口驱动的散射参数问题用向量拟合或扫频策略，射频器件的 S 参数曲线往往扫一次要数小时——自适应频点采样（在谐振附近自动加密采样）因此成为商用软件的标配。

# 一维亥姆霍兹方程的最简有限元：线性元解腔内驻波
import numpy as np

k = 2*np.pi/1.0            # 波长一米的腔
L, N = 1.0, 40             # 腔长一米 剖四十段
x = np.linspace(0, L, N+1)
h = x[1]-x[0]

K = np.zeros((N-1, N-1))
for i in range(N-1):
    K[i, i] = 2/h - (k**2)*h*(2/3)
    if i < N-2:
        K[i, i+1] = -1/h - (k**2)*h/6
        K[i+1, i] = K[i, i+1]

b = np.zeros(N-1); b[0] = 1.0        # 左端激励 单位场
e = np.linalg.solve(K, b)
x_in = x[1:-1]
exact = np.sin(k*x_in)/np.sin(k*L)   # 驻波解析解
err = np.max(np.abs(e-exact))
print(f"最大节点误差 {err:.2e}（线性元量级，随网格平方收敛）")
# 把 N 翻倍重跑：误差应降到约四分之一——二阶收敛的直接验证
​

把 N 加倍重跑、观察误差缩到约四分之一，就是数值方法课里"二阶收敛"最省事的现场验证——值班室自己的仪器自检。

三、棱边元与吸收边界：两处容易翻车的现场

FEM 的坑集中在两处，都值得提前踩过。第一处叫伪解：如果直接用节点插值（把电场三个分量当普通标量场插值），在波导与腔体的本征模求解里会算出一批物理上不存在的"模"，它们是离散格式违反了场量散度约束的产物。解药是棱边元（Nédélec 单元）：把自由度安插在四面体的棱边上，插值对象改为沿棱边的切向场分量，切向连续天然保证，伪解从源头消失。所有正经的高频 FEM 软件默认用棱边元，理解这一点，就理解了为什么同一模型在不同软件间移植时自由度数目对不上。第二处是开放域截断：FEM 要求计算域有界，人工截断面上 outward 行波必须被"吸收"而不是反射回来。一阶吸收边界只对法向入射近似有效，斜入射反射明显；工程主流是完美匹配层——在截断面外裹一层人为设计的各向异性有耗介质，理论上对任意入射角、任意极化的外行波阻抗完全匹配，数值上层数取十层上下即可把反射压到实用水平（介质薄层贴太近会让场尚未衰减就触底反弹，这与 FDTD 的吸收边界纪律相通）。

谐振腔调频是 FEM 日常出勤里最常见的活，灵敏度可以先手算估计：

# 圆柱腔 TE011 模的谐振频率与端盖调谐灵敏度
import math

c = 299792458.0
x01 = 1.8412                      # 贝塞尔函数导数第一个零点，TE01 模专用
def te011(radius, depth):
    return c / (2*math.pi) * math.sqrt((x01/radius)**2 + (math.pi/depth)**2)

f0 = te011(0.05, 0.05)
f1 = te011(0.05, 0.051)           # 可动端盖外拉一毫米
print(f"标称频率 {f0/1e9:.3f} GHz，端盖外拉 1 毫米后 {f1/1e9:.3f} GHz")
print(f"调谐灵敏度约 {(f0-f1)/1e6:.0f} MHz/mm，量程约覆盖几十兆赫兹")
​

一毫米行程换来约五十兆赫兹的频率搬移，这类手算灵敏度是建模前的"预飞行检查"：仿真值若与它差出量级，先怀疑边界条件与模式判读，而不是急着加密网格。

四、与 MoM 的交接与混合框架

FEM 的软肋恰是 MoM 的强项：FEM 要求计算域有限（要人为截断并配吸收边界或完美匹配层），无截断的开放辐射问题做起来别扭；而 MoM 天生开放域、却不擅复杂介质。自然的结论是混合：目标表面用 MoM 的未知电流，复杂介质区用 FEM 的体网格，两者在交界面上通过等效原理耦合（1.2 节边界条件的数值化身影）。天线带罩分析（天线 MoM + 罩体介质 FEM）、埋地散射（地表 MoM + 地下 FEM）是经典组合。混合框架的交接纪律要严：交界面上场连续性条件的离散精度决定整体精度，两个域的网格密度要匹配，否则短板效应让贵的那半白算。

把三大数值方法的分工摊在一张桌上对照，选型时不必临场犹豫：

表格读法：没有全能冠军，只有任务匹配。宽带、瞬态找时域；腔体、复杂介质找有限元；开放域的高精度阻抗与散射找矩量法；三者互相交接的混合框架，正是现代大型仿真软件的卖点所在——调度值班室从不押宝单一工具，而是按订单派单。

本节要点回顾

• 变分哲学：真实场使能量泛函驻定，复杂材料只改被积函数，天然兼容各向异性与非均匀。
• 稀疏红利：四面体剖分带来稀疏矩阵，内存与求解近似线性增长，与 MoM 的稠密账单互补。
• 自适应加密：误差估计驱动局部加密，算力流向缝隙与尖角，收敛有能量范数保底。
• 本征与扫频：谐振找模用特征值问题，散射参数靠自适应频点采样，谐振附近自动加密。
• 混合耦合：MoM 管开放表面、FEM 管复杂介质，交界面连续性是整体精度的命门。

下一节补齐值班室最后一批工具：电大尺寸的高频渐近捷径、算法混编的交接纪律，以及把这一切装配成生产力的大型仿真软件。