静电场边界元法
理论与物理
静电场BEM
老师，静电场也能用BEM吗？
开放空间问题BEM最合适。FEM需要将周围空间网格化，而BEM只需处理导体/电介质的表面。
格林函数 $G = 1/(4\pi r)$（3D）。以表面电位和法向电场为未知数求解。
只用表面网格就能处理无限区域！
是的。在输电线路周围电场、雷击浪涌屏蔽效果、EMC问题的开放空间分析中能发挥强大作用。
FEM与BEM的比较
总结
- 仅用表面网格即可求解开放空间 — BEM的最大优点
- 与FEM互补 — 封闭区域用FEM，开放空间用BEM
- FEM-BEM耦合 — 内部FEM + 外部BEM的混合方法
边界积分方程的起源——19世纪数学家格林的孤独发现
作为静电场BEM理论基础之一的“格林定理”的发表者乔治·格林（1793〜1841），是诺丁汉一位面包师的儿子，自学成才的数学家。他于1828年自费出版的论文生前几乎被忽视，但在其去世后被开尔文勋爵重新发现，成为电磁学的基石。近200年过去了，其格林函数仍在BEM的核心中延续生命，这真是一个惊人的故事。
各项的物理意义
- 电场项 $\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假设（截面分析）：电流方向均匀，可忽略端部效应时有效
- 各向同性假设：各向异性材料（如硅钢板的轧制方向等）需要定义方向特性
- 不适用的案例：等离子体（电离气体）、超导体、非线性光学材料需要额外的本构关系
数值解法与实现
BEM的离散化
表面用三角形/四边形单元离散化。以各单元的电势$\phi$和法向电通量密度$D_n$为未知数构建联立方程:
影响矩阵$[H], [G]$是格林函数的面积分。是稠密矩阵，因此内存为$O(N^2)$。
FMM加速
大规模问题中使用FMM（快速多极子法）可加速至$O(N\log N)$。FastCap和COMSOL的BEM模块中已实现。
总结
- $[H]\{\phi\} = [G]\{D_n\}$ — BEM的基本方程
- 稠密矩阵 → 用FMM加速
BEM的“格林函数”并非易事
边界元法的核心是格林函数——用数学公式表示“某点放置的点电荷对周围电势的影响”。静电场情况下是 $1/r$（距离的倒数）这种相对简单的形式，但对于各向异性电介质或多层介质，则会立刻变得复杂。在半导体封装的多层基板中，每个电介质层都需要不同的格林函数，为了精确进行其积分所需的数值技术也是研究的前沿领域。
边元（Nedelec元）
专用于电磁场分析的单元。自动保证切向分量的连续性，排除伪模式。是3D高频分析的标准。
节点元
用于标量势公式化。在静磁场的标量势法或静电场分析中有效。
FEM 对比 BEM（边界元法）
FEM: 对应非线性材料、非均匀介质。BEM: 可自然处理无限区域（开域问题）。混合FEM-BEM也有效。
非线性收敛（磁饱和）
用牛顿-拉弗森法处理B-H曲线的非线性。残差标准: $||R||/||R_0|| < 10^{-4}$ 为通用标准。
频域分析
通过时间谐波假设归结为稳态问题。需要复数运算，但宽带特性需通过时域分析获取。
时域的时间步长
需要最高频率成分的1/20以下的时间步长。隐式时间积分中也可使用更大的步长，但需注意精度。
频域与时域的区分使用
频域分析类似于“调谐到收音机的特定频率”——可以高效计算单一频率下的响应。时域分析类似于“同时录制所有频道”——可以再现包含所有频率成分的瞬态现象，但计算成本高。
实践指南
实务
输电线路的电场分析、EMC屏蔽效果、雷击脉冲分析是主要的应用方向。
检查清单
- [ ] 表面网格的质量（三角形纵横比）是否足够
- [ ] 奇异积分的处理是否正确（自单元的积分）
- [ ] 观测点是否在表面上（奇异点处值会发散）
高压输电绝缘子设计中选用BEM的理由
输电铁塔的绝缘子（绝缘体）为了防止污损（灰尘或海盐附着）导致的“闪络”，精密的电位分布分析不可或缺。绝缘子形状是复杂的曲面，周围是无限广阔的大气——这正是BEM擅长的场景。实际上，在275kV或500kV的超高压绝缘子设计中，会使用BEM进行形状优化，以使表面的电场强度峰值控制在20kV/mm以下。
分析流程的比喻
电机的电磁场分析类似于“给吉他调音”。调整弦的粗细（线圈匝数）和琴桥的位置（磁铁配置），以引出最美妙的音色（高效的扭矩特性）。改变一个参数，整体的平衡就会改变——所以参数化研究很重要。
初学者容易陷入的陷阱
“空气区域？为什么要用网格切割空气？”——这是几乎所有初次接触电磁场分析的人都会产生的疑问。答案是“因为磁力线也会扩散到铁心之外”。如果将分析区域恰好限定在铁心边界，无处可去的磁通会“撞上”边界壁并反射，产生实际中不可能出现的磁通集中。想象一下房间太小，球在墙上不断弹跳的状态。
边界条件的思考方式
远方的边界条件虽然不起眼但超级重要。需要在数值上表达“从这里开始是无限广阔的空间”。如果设置错误，磁通就会像撞上“看不见的墙”一样被反射回来。
软件比较
工具
| 工具 | 特点 |
|---|---|
| COMSOL (BEM模块) | 支持FEM-BEM耦合 |
| FastCap（MIT） | BEM静电容。开源 |
| Ansys Q3D | 内部使用BEM |
| Integrated Engineering Software | BEM专用。Coulomb, Amperes |
“专注于BEM的工具”意外稀少的原因
纵观市场，专注于静电场BEM的商业工具意外地少，很多情况下是以“为FEM添加BEM功能”的形式存在。原因是实现的难度——完全的BEM会生成稠密矩阵，因此在大规模问题中内存和计算成本会爆炸性增长。直到2000年代FMM（Fast Multipole Method，快速多极子法）实用化后，大规模BEM才变得现实。这也是工具供应商偏向FEM的原因之一。
选型时最重要的三个问题
- “要解决什么问题”：所需的物理模型、单元类型是否支持静电场边界元法。例如，流体中是否有LES支持，结构中是否有接触、大变形的支持能力会造成差异。
- “由谁使用”：新手团队适合GUI丰富的工具，有经验者适合脚本驱动的灵活工具。类似于汽车的自动挡（GUI）和手动挡（脚本）的区别。
- “要扩展到什么程度”：基于未来分析规模扩大（HPC支持）、向其他部门扩展、与其他工具的联动等长远考虑进行选择，有助于长期降低成本。