高周波FEM
理论与物理
高频电磁场的FEM
老师，用FEM求解高频电磁场时的控制方程是什么？
亥姆霍兹方程（频域）：
$k_0 = \omega/c = 2\pi f/c$: 自由空间波数。使用棱边单元（Nédélec单元）离散化，求解复对称矩阵。
HFSS用的就是这种方法吧。
是的。Ansys HFSS是频域FEM的代表。四面体网格+二阶棱边单元可以高精度分析复杂形状。其自适应网格功能会自动迭代直至S参数收敛。
总结
- 亥姆霍兹方程 — 频域的波动方程
- 棱边单元 — 矢量场的自然离散化
- 复对称矩阵 — 使用直接法或迭代法求解
高频FEM的“数值色散”问题——波长越短误差越累积
高频FEM存在一个本质性的问题，即“数值色散”。这是在FEM近似中传播的电磁波的相速度偏离真实光速的现象。网格越粗、频率越高，这个误差就越大（相位误差随传播距离成比例累积）。对策有三：① 细化网格（单元边长控制在波长的1/6~1/10以下），② 使用高阶单元（二阶/三阶多项式），③ 采用数值色散较少的DG-FEM（不连续伽辽金有限元法）等方法。到了毫米波频段（30~300 GHz），这个问题会直接影响设计精度。
各项的物理意义
- 电场项 $\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial t$：法拉第电磁感应定律。随时间变化的磁通密度产生电动势。【日常例子】自行车发电机（发电机）通过旋转磁铁使附近的线圈产生电压——这是磁场随时间变化会感应出电场这一定律的直接应用。IH电磁炉也基于相同原理，高频磁场的变化在锅底感应出涡流，通过焦耳热加热。
- 磁场项 $\nabla \times \mathbf{H} = \mathbf{J} + \partial \mathbf{D}/\partial t$：安培-麦克斯韦定律。电流和位移电流产生磁场。【日常例子】电线通电时周围会产生磁场——这就是安培定律。电磁铁根据此原理工作，通过给线圈通电产生强磁场。智能手机的扬声器也是应用了电流→磁场→振膜受力这一原理。在高频（GHz频段天线等）下，位移电流 $\partial D/\partial t$ 不可忽略，它描述了电磁波的辐射。
- 高斯定律 $\nabla \cdot \mathbf{D} = \rho_v$：表明电荷是电通量的发散源。【日常例子】用垫板摩擦头发会产生静电使头发竖起——带电的垫板（电荷）放射状地发出电力线，对轻的头发施加力。设计电容器时，电极间的电场分布就是根据此定律计算的。ESD（静电放电）对策也基于高斯定律的电场分析。
- 磁通连续性 $\nabla \cdot \mathbf{B} = 0$：表示不存在磁单极子。【日常例子】将条形磁铁切成两半，也无法得到只有N极或只有S极的磁铁——N极和S极总是成对存在。这意味着磁力线描绘的是“没有起点和终点的闭合回路”。在数值分析中，为了满足此条件，采用矢量势 $\mathbf{B} = \nabla \times \mathbf{A}$ 的公式化方法，自动保证磁通连续性。
假设条件与适用范围
- 线性材料假设：磁导率、介电常数不依赖于磁场、电场强度（饱和区域需要非线性B-H曲线）
- 准静态近似（低频）：可忽略位移电流项（$\omega \varepsilon \ll \sigma$）。涡流分析中常用
- 2D假设（截面分析）：电流方向均匀且可忽略端部效应时有效
- 各向同性假设：对于各向异性材料（如硅钢板的轧制方向等）需要定义方向特性
- 不适用的情形：等离子体（电离气体）、超导体、非线性光学材料需要额外的本构关系
数值解法与实现
FEM的实现
请告诉我高频FEM实现上的注意事项。
频率扫描怎么做呢？
HFSS的快速频率扫描，通过少数几个直接解，用有理函数近似来插值整个频带的S参数。无需在每个频率点重新求解，因此非常高效。
总结
- $\lambda/5$〜$\lambda/10$网格 — 二阶单元的情况
- DDM — 区域分割并行实现大规模化
- 快速频率扫描 — 有理函数插值实现高速
棱边单元（Nedelec单元）——为何高频FEM中必不可少
在高频电磁场FEM中使用通常的节点单元（拉格朗日单元）会出现“伪解（虚假模式）”这种物理上不存在的虚假本征模式。为了防止这种情况而开发的就是棱边单元（也称为Nedelec单元、Whitney单元），它将自由度分配给单元的边而不是节点，自动满足∇·E=0（电场散度为零）条件。以HFSS为代表的主要高频FEM工具几乎都基于棱边单元。如果被问到“为什么HFSS使用Nedelec单元？”，最简单的答案就是“为了消除伪解”。
棱边单元（Nedelec单元）
专为电磁场分析设计的单元。自动保证切向分量的连续性，排除虚假模式。是3D高频分析的标准。
节点单元
用于标量势公式化。在静磁场的标量势法或静电场分析中有效。
FEM vs BEM（边界元法）
FEM: 对应非线性材料、非均匀介质。BEM: 自然处理无限域（开域问题）。混合FEM-BEM也有效。
非线性收敛（磁饱和）
使用牛顿-拉夫森法处理B-H曲线的非线性。残差标准: $||R||/||R_0|| < 10^{-4}$ 是通用的。
频域分析
通过时间谐波假设归结为稳态问题。需要进行复数运算，但宽带特性需通过时域分析获取。
时域的时间步长
需要最高频率成分的1/20以下的时间步长。隐式时间积分中也可使用更大的步长，但需注意精度。
频域与时域的区分使用
频域分析类似于“将收音机调到特定频率”——可以高效计算单一频率下的响应。时域分析类似于“同时录制所有频道”——可以再现包含所有频率成分的瞬态现象，但计算成本高。
实践指南
实务中的应用
波导滤波器、RF-MEMS、封装内的电磁场分析、连接器的S参数提取是典型应用。
实务检查清单
HFSS的自适应网格——“Delta S”的真面目
使用ANSYS HFSS时，在分析设置中有一个称为“最大Delta S”的收敛判定参数。这是指连续两次网格细化循环之间，S参数（反射系数、传输系数）的变化量低于此值时，即判定为“收敛”的阈值。默认值0.02（2%）对于多数情况足够，但如果需要精密评估滤波器的带外衰减特性或天线的宽带特性，则建议设置为0.005以下。不过，设置得越严格，通过的循环次数就越多，计算时间也会延长，因此根据要求的精度选择合适的值是实务的关键。
分析流程的比喻
电机的电磁场分析感觉上类似于“给吉他调音”。调整琴弦的粗细（线圈匝数）和琴桥的位置（磁铁配置），以引出最美的音色（高效的扭矩特性）。改变一个参数，整体的平衡就会改变——所以参数化研究很重要。
初学者容易陷入的误区
“空气区域？为什么要用网格划分空气？”——这是几乎所有初次接触电磁场分析的人都会产生的疑问。答案是“因为磁力线也会扩散到铁心之外”。如果将分析区域设置得紧贴铁心，无处可去的磁通会“撞上”边界壁面并反射，产生实际中不可能出现的磁通集中。想象一下房间太小，球在墙壁上不断弹跳的状态。
边界条件的思考方式
远场边界条件虽然不起眼但超级重要。需要在数值上表达“从这里开始是无限广阔的空间”。如果设置错误，磁通就会像撞上“看不见的墙壁”一样被反射回来。
软件比较
工具
| 工具 | 特点 |
|---|---|
| Ansys HFSS | 高频FEM的行业标准。自适应网格+DDM |
| COMSOL RF | 多物理场耦合。易于热-RF耦合 |
| CST Studio Suite | 也搭载FEM求解器（频域） |
| Altair FEKO | FEM+MoM混合。大规模问题 |
高频
相关主题
なった
詳しく
報告