高斯定律（静电场）
高斯定律（静电场）的理论基础
概述
教授，高斯定律是麦克斯韦方程之一吧？在静电场分析中如何使用？
高斯定律是"通过闭合曲面的电通量总量等于其内部全部电荷"的定律。对应麦克斯韦方程的第一式。
微分形式为
其中 $\mathbf{D} = \varepsilon \mathbf{E}$ 是电通量密度，$\rho$ 是体积电荷密度。
这与FEM静电场分析的关系是什么？
将 $\mathbf{E} = -\nabla\phi$ 代入可得泊松方程。
FEM求解器通过离散化求解该方程。也就是说，高斯定律是FEM静电场分析的根本方程。
对称性利用的解析解
从高斯定律可以直接求得电场的情况也存在吧？
仅在具有高度对称性的情况下才能求解。
| 对称性 | 高斯面 | 电场 |
|---|---|---|
| 球对称（点电荷） | 同心球面 | $E = Q/(4\pi\varepsilon r^2)$ |
| 圆柱对称（线电荷） | 同轴圆柱面 | $E = \lambda/(2\pi\varepsilon r)$ |
| 平面对称（面电荷） | 直方体 | $E = \sigma/(2\varepsilon)$ |
这些解析解对FEM结果的有效性验证至关重要。无论是COMSOL还是Ansys Maxwell，首先应在对称问题上与理论值进行比较，然后再进行复杂问题的分析，这是实务的基本原则。
"闭合曲面的形状无关紧要"——高斯定律的巧妙自由度
高斯定律的妙处在于"无论闭合曲面是什么形状，其内部的电荷就能决定积分值"。球形、立方体、奇怪的形状都可以。利用这一点，对于对称性强的问题（球状电荷或无限长圆柱电荷）计算会大大简化。比如高压电缆的截面——中心有芯线，外边是屏蔽层——这种几乎完全圆柱对称的结构中，只需将高斯面设为同轴圆柱形，就能一下子求出电场分布。FEM当然能数值求解，但如果想"验证解析解与数值解是否一致"时，这种直观的解法在实际工作中就派上大用场了。
高斯定律（静电场）的数值计算手法
数值解法的详细
高斯定律导出的泊松方程用FEM如何离散化？
Galerkin法的弱形式化是基础。用试验函数 $w$ 相乘然后分部积分可得
离散化后得到 $[K]\{\phi\} = \{f\}$。形式与结构FEM的刚度方程完全相同。
介电常数对应杨氏模量？
正是如此。因此，有结构FEM经验的人很容易理解静电场FEM。
高斯定律的电荷计算
从FEM结果后处理中求得导体上的电荷如何操作？
利用高斯定律对导体表面的电通量密度进行积分。
COMSOL具备表面积分的后处理功能，可以直接计算。由此值推导出静电容量 $C = Q/V$。
容量矩阵的提取
多导体系统的容量矩阵怎样求得？
将导体 $j$ 设定为1V，其他所有导体设为0V，计算各导体上的感应电荷。对所有导体重复此操作就能得到容量矩阵 $C_{ij}$。Ansys Q3D完全自动化了这个操作，在PCB布线的寄生容量提取中被广泛使用。
高斯面的选取方式对FEM精度的影响
数值分析中利用高斯定律验证电荷量时，"在哪里设置高斯面"会大大影响结果精度。如果高斯面穿过网格粗糙的区域，面积分的误差就会增大，电荷收支会不平衡。实务中的做法是"在电极表面与分析空间之间插入网格充分细密的区域后再设置高斯面"。将形状接近球体或直方体可以简化积分，解析后用高斯定律进行验证能有效检查"FEM是否真的正确计算了电荷"。
高斯定律（静电场）的实务应用
实践指南
实务中进行静电场分析的流程请教教我。
以高压设备的绝缘设计为例来讲解。
Step 1：模型构建
- 电极形状的CAD导入
- 绝缘材料的介电常数设定（环氧树脂：$\varepsilon_r=3.5$、油：$\varepsilon_r=2.2$、SF6：$\varepsilon_r=1.0$）
- 空气/周边区域的追加
Step 2：边界条件
- 高压电极的Dirichlet条件（$\phi = V_0$）
- 接地电极设为 $\phi = 0$
- 外部边界设无穷远条件
Step 3：网格与求解
- 电极边缘部细网格（曲率半径的1/5以下）
- 采用二阶单元（电场精度考量）
电场强度的评估基准是什么？
与材料各自的绝缘耐压值进行比较。
| 材料 | 绝缘耐压值 [kV/mm] | 安全系数 |
|---|---|---|
| 空气（1atm） | 3.0 | 2～3 |
| 环氧树脂 | 20～30 | 3～5 |
| SF6（0.5MPa） | 8.9 | 2.0 |
三相接触点（三重点）的电场集中特别要注意吧。
完全同意。导体、绝缘体、气体相交的三重点是电场集中的温床。在COMSOL中通过对三重点附近进行网格细化处理是实务的标准做法。
高压送电线绝缘子设计中使用高斯定律的故事
在高压送电线中电压可达数百kV，因此电杆上附装的"绝缘子"的绝缘设计至关重要。绝缘子的形状（伞片数量和间隔）改变时，表面的电场分布也会改变，在有雨水和污垢时的绝缘性能会大幅波动。实务中利用基于高斯定律的静电场分析来预测绝缘子周围电场的最大值，进行设计。雷电浪涌发生时的瞬间电场集中在哪里，污垢条件下的闪络概率是多少——这类问题都用静电场CAE来求解。虽然不太起眼，但对送电线的稳定运行来说是不可或缺的技术。
高斯定律（静电场）的软件对比
商用工具对比
请比较静电场分析工具。
按功能别来整理主要工具。
| 功能 | COMSOL | Maxwell | FEMM | Q3D |
|---|---|---|---|---|
| 2D静电场 | ○ | ○ | ○ | - |
| 3D静电场 | ○ | ○ | - | ○ |
| 自动自适应网格 | ○ | ○ | - | ○ |
| 无穷远边界 | ○ | ○ | △ | ○ |
| 容量矩阵自动提取 | ○ | ○ | △ | ○ |
| 多物理耦合 | ○ | △ | - | - |
COMSOL与Maxwell的选用如何区分？
COMSOL可以自由追加自定义方程，适合非线性介电体和空间电荷注入等研究型问题。Maxwell的自动自适应网格功能强大，初学者也能少费工夫进行网格设计。
FEMM在实务中也能用吗？
限于2D问题范围内完全够用。同轴构造或导线周围的电场计算用免费软件就能得到高精度结果。有Lua脚本API可以进行自动化。但3D问题需要COMSOL或Maxwell。
高斯定律验证功能——工具选择的隐形指标
在评估商用静电场分析工具时，"是否具有基于高斯定律的电荷收支验证功能"是个微妙但重要的指标。将解析结果的电通量密度 $\mathbf{D}$ 对闭合曲面进行积分，验证是否与内含电荷量相符的工具信度更高。ABAQUS和COMSOL能在后处理中进行此验证，但有些工具可能没有列为标准功能或操作比较复杂。高压绝缘设计现场中，有的团队设定了"用高斯定律验证电荷收支在1%以内"这样的质量标准，足见该验证在实务中的重要程度。
高斯定律（静电场）的前沿研究
前沿话题
关于高斯定律相关的最新研究动向请教一下。
有几个重要课题。
非线性介电体的电场解析
在强铁电体（BaTiO3、PZT）中，$\mathbf{D}$ 与 $\mathbf{E}$ 的关系呈非线性且显示迟滞。高斯定律本身不变，但本构关系变成非线性需要迭代求解。在COMSOL中实装Preisach模型的研究正在推进。
概率论绝缘设计
听说从决定论安全系数向概率论设计转变。
用Weibull分布统计处理绝缘破坏。基于FEM电场分析结果，计算各单元的破坏概率。
在GIS（气体绝缘开关柜）和HVDC电缆设计中开始纳入IEC规格。
电气流体力学（EHD）耦合
绝缘油中的空间电荷由对流输送，影响电场分布。变压器油中的电场分析需考虑EHD耦合。用COMSOL的CFD模块与AC/DC模块进行耦合求解的方法备受瞩目。
为什么半导体仿真中高斯定律不可或缺
在半导体器件（晶体管、二极管）的仿真中，必须联立求解载流子（电子和空穴）的输运方程与泊松方程（出自高斯定律）。即使器件尺寸小到几十纳米，从高斯定律来的 $\nabla \cdot \mathbf{D} = \rho$ 也是支配方程的中心。在最先进的2nm工艺节点中，栅长仅有数十个原子那么小，"经典连续体电磁学是否仍然成立？"的研究在进行。但目前TCAD（工艺仿真工具）的主流仍基于高斯定律，在芯片设计中发挥实时作用。
高斯定律（静电场）的故障排除
故障排除
静电场分析常见的故障请教一下。
归纳一下代表性问题。
1. 高斯面积分电荷不匹配
导体表面的积分电荷与施加电荷不符。
对策：细化网格。至少要将导体面的单元尺寸设为特征尺寸的1/20以下。COMSOL中用表面积分的"es.Dn"指定导体全表面进行计算。
2. 浮动导体电位不定
既不加电压也不接地的导体置入后收敛失败。
浮动导体需"Floating Potential"边界条件。COMSOL和Ansys Maxwell都有同名设定。忘记这个设置会导致矩阵奇异。
3. 分析区域大小的影响
改变区域后答案也变了。
对策：用COMSOL的"无穷远单元域"或Maxwell的"Radiation Boundary"。FEMM则应采用对象20倍以上的区域。扩大区域至2倍后结果在1%以内不变就可确认收敛。
4. 介电体界面处电场不连续
不同介电常数的边界处电场跳变是否正常？
这是正常的。法向电通量密度 $D_n = \varepsilon E_n$ 连续，但电场 $E_n$ 不连续。满足 $\varepsilon_1 E_{n1} = \varepsilon_2 E_{n2}$。切向电场连续。为了正确反映这种物理现象，网格应与界面相对齐。
"电荷消失了！"——电荷收支不平衡的原因追查
进行静电场分析时用高斯定律进行电荷收支检查，结果"投入的电荷量与积分值相差很大"的情况时有发生。原因多数是网格粗糙，但陷阱有时也很隐蔽。比如穿过介电体边界的面单元会使 $\mathbf{D}$ 的法向成分产生不连续，收支就会错乱。使用周期边界条件时也容易把边界面重复计算。"用高斯定律检查收支不平时＝模型或网格有问题"这个信号可用来进行有效的故障排除，使整个分析过程更有效率。
帮助
更详细
错误