有限元分析电场 — 静电场与高频电磁场仿真的数值方法复盘

一、背景：电场仿真的特殊数值挑战

电场仿真是有限元分析中一个比较特殊的领域——它不涉及经典的”力-位移”关系，而是求解Maxwell方程组或其简化形式。电场仿真的挑战来自几个方面：电场在导体边缘的奇异性（尖角处电场趋于无穷）、介质界面的法向电场不连续性、以及高频下电磁波在空间中的传播与衍射需要高波数分辨率。

我实际处理过的电场仿真项目包括：高压绝缘子电场分布优化、PCB板间电磁兼容（EMC）分析、变压器油纸绝缘结构设计、以及射频滤波器的S参数提取。不同类型的问题需要不同的求解策略。

二、静电场分析

2.1 控制方程与边界条件

静电场分析求解拉普拉斯方程或泊松方程：

∇·(ε∇φ) = -ρ

其中φ是电势，ε是介电常数，ρ是自由电荷密度。边界条件有三类：

• Dirichlet（第一类）：指定导体表面的电势（如高压端100 kV，接地端0 V）
• Neumann（第二类）：指定表面的电场法向导数/表面电荷密度（如对称面的零法向电场）
• Floating Potential：导体为等势体但总电荷量固定（如悬浮金属颗粒）

**漂浮导体（Floating Conductor）**在有限元分析电场中是一个容易被忽视但很重要的边界条件。绝缘结构中的微小金属杂质（如变压器油中的金属颗粒）在电场中感应出净电荷，成为一个”漂浮导体”——它的电势不是人为给定的，而是通过总电荷为零的条件求解得到的。如果不处理为漂浮导体，简单设为0 V会严重低估局部电场强度（偏小50-80%）。

2.2 电场奇异性与网格细化

导体尖角处的电场理论上趋于无穷（数学奇异性）。有限元网格无法正确解析奇异场——网格越细，电场越高，不会收敛到一个有限值。这与力学中的应力奇异是同一个问题。

处理方案：

1. 圆角替代尖角：将实际的微小圆角（如0.1-0.5 mm）几何建模出来。这是一个”物理修正”——现实中不存在数学意义上的完美尖角
2. 不直接使用奇异点的场值：在奇异点以外的位置（如距离尖角>1mm）评估电场强度
3. 子模型技术：用全局模型得到电势分布，用局部子模型在尖角区域加精细圆角后重新计算

2.3 自适应网格与误差估计

ANSYS Maxwell和COMSOL都支持基于电场梯度的自适应网格细化（h-adaptivity）。关键设置是自适应循环的终止条件：

• 能量误差：全局能量误差<1%（COMSOL默认<0.1%）
• 局部场误差：每个单元的局部电场误差<5%
• 自适应步数：通常3-5次自适应循环即可收敛

经验：如果经过5次自适应后误差仍未收敛——通常是几何中存在未处理的尖角，即使无限细化也无法收敛。回到几何模型，找出并修正尖角。

三、交流/准静态电场

3.1 准静态假设的适用边界

当电场的频率满足一定条件时，可以使用准静态近似（忽略位移电流中的磁感应项）。对于电力设备（50/60 Hz），波长为6000/5000 km，远大于设备尺寸，准静态近似完全有效。

但对于高频设备（如500 MHz的RF器件，波长≈0.6 m），波长与器件尺寸可比，必须使用全波求解。

判断准则：如果器件最大尺寸<波长的1/10，准静态近似误差<3%，可以接受。

3.2 复介电常数与损耗正切

交流电场下，介质的介电常数表达为复数：

ε* = ε’ – jε”

其中ε’是介质储存电能的能力，ε”代表介质损耗（tanδ=ε”/ε’）。有限元分析电场中的介质损耗直接转化为热：

p_loss = ωε”|E|²/2

这是评估介质是否会发生热击穿的关键参数。对于环氧树脂（tanδ≈0.001-0.01），50 Hz下的介质损耗通常可以忽略。但对于油浸纸（tanδ≈0.003-0.02），在高压直流（HVDC）场合，由于谐波含量高，介质损耗可能造成显著的局部温升。

四、高频电磁场（全波分析）

4.1 网格密度与波长

高频电磁场的有限元分析（全波求解）中，网格密度的要求比静电场严格得多——每个波长至少需要10个二阶单元（或20个一阶单元）来正确解析波的相位变化。这被称为”10点每波长”规则。

对于10 GHz的电磁波（波长30 mm），如果仿真域为100×100×50 mm³，二阶四面体的特征单元尺寸为3mm，需要的单元数约为50万——计算量对于单台工作站完全可控。但如果频率升到60 GHz（波长5 mm），同样的物理空间中需要的单元数量增加216倍（体积与波长的立方成反比），单元数超过1亿——需要集群。

4.2 完美匹配层（PML）

与声学仿真类似，高频电磁场中也用PML吸收外向传播的电磁波。PML的性能取决于三个参数：厚度、层数和多项式阶数。COMSOL中推荐设置：

• PML厚度：1/2波长
• PML层数：8层
• 多项式阶数：3阶

在这些设置下，PML对法向入射波的反射系数<10⁻⁴。

4.3 端口设置：波端口 vs 集总端口

RF器件的激励通常通过端口（Port）设置。ANSYS HFSS和COMSOL提供了两种端口类型：

• 波端口（Wave Port）：假设端口处为无限长的均匀传输线结构（如波导、同轴线），通过2D特征模求解得到端口处的场分布。精度高但要求端口处的传输线足够长（>3倍截面尺寸）。
• 集总端口（Lumped Port）：在端口处施加一个集总电压/电流激励。不需要端口延伸，但精度略低于波端口（特别是对于非TEM模的传输线）。

对于微带线和共面波导，集总端口通常足够；对于矩形波导和同轴线，波端口精度更高。

五、多物理场耦合

5.1 电-热耦合

电流流过导体产生焦耳热（电阻损耗）和介质中的介质损耗都会导致温升。电-热耦合分析的难点在于材料属性的温度依赖性——电导率和介电常数都随温度变化。

耦合策略：电-热弱耦合（顺序耦合）适用于大多数场景——先跑电场分析得到损耗分布，再将损耗作为热源跑热分析。如果温度变化导致电导率变化超过30%，需要迭代耦合（多次电→热→更新电导率→重新电→…）。

5.2 电-力耦合

高压设备中的电场力（Maxwell应力张量）可能引起结构的机械振动。对于交流设备（50 Hz），电场力的频率是100 Hz（两倍基频）——如果这个频率接近结构的固有频率，会发生谐振。

对于GIS（气体绝缘开关设备）的母线，电动力引起的振动幅值通常为μm量级，但在谐振条件下可能放大10-50倍，产生可听噪声（嗡嗡声）。

六、复盘

有限元分析电场问题的核心挑战因频率而异：

• DC/工频：奇异性处理（避免在尖角处取场值）+ 浮动导体边界
• 准静态：介质损耗评估 + 复介电常数
• 高频：波长分辨率（10点/波长）+ PML + 端口设置

电场不是”算一下就出来”——每个频率段有每个频率段的坑。