静电场边界元法
静电场边界元法的理论基础
静电场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的核心中延续生命，这真是一个惊人的故事。
数值解法与实现
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$（距离的倒数）这种相对简单的形式，但对于各向异性电介质或多层介质，则会立刻变得复杂。在半导体封装的多层基板中，每个电介质层都需要不同的格林函数，为了精确进行其积分所需的数值技术也是研究的前沿领域。
静电场边界元法静电场边界元法实践指南
实务
输电线路的电场分析、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的原因之一。