高斯定律（静电场）
理论与物理
概述
老师，高斯定律是麦克斯韦方程组之一吧？在静电场分析中怎么使用呢？
高斯定律是“穿过闭合曲面的电通量总和等于其内部的总电荷”的定律。对应麦克斯韦方程组的第一式。
微分形式为
这里 $\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求解当然也能数值得到，但在“想确认是否与解析解一致”时，这种直观的解法在现场非常有用。
各项的物理意义
- 电场项 $\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公式化从高斯定律导出的泊松方程？
基于Galerkin法的弱形式化是基础。乘以权函数 $w$ 并进行分部积分得到
离散化后得到 $[K]\{\phi\} = \{f\}$。与结构FEM的刚度方程形式完全相同。
介电常数对应杨氏模量对吧。
没错。所以如果有结构FEM的经验，就很容易理解静电场FEM。
基于高斯定律的电荷计算
如何通过后处理从FEM结果求导体上的电荷？
使用高斯定律对导体表面的电位移矢量进行积分。
在COMSOL中，可以使用表面积分的后处理功能直接计算。根据此值可以算出静电容量 $C = Q/V$。
电容矩阵提取
多导体系统的电容矩阵怎么求？
将导体 $j$ 设为1V，其他所有导体设为0V，计算各导体上的感应电荷。对所有导体重复此操作即可得到电容矩阵 $C_{ij}$。Ansys Q3D中此操作完全自动化，广泛用于PCB布线的寄生电容提取。
高斯面的选择影响FEM精度
在数值分析中使用高斯定律验证电荷量时，“高斯面取在哪里”会极大地影响结果的精度。如果高斯面穿过网格粗糙的区域，面积分的误差会变大，导致电荷收支不平衡。在实际工作中，常采用“在电极表面和分析空间之间，设置网格足够精细的区域后再设定高斯面”的技巧。有时将形状接近球体或长方体会使积分更简单，在分析后的验证中熟练运用高斯定律，可以检查“FEM是否真的计算出了正确的电荷”。
边单元（Nedelec单元）
专用于电磁场分析的单元。自动保证切向分量的连续性，排除伪模式。是3D高频分析的标准。
节点单元
用于标量势公式化。在静磁场的标量势法或静电场分析中有效。
FEM vs BEM（边界元法）
FEM: 对应非线性材料·非均匀介质。BEM: 能自然处理无限域（开域问题）。混合FEM-BEM也有效。
非线性收敛（磁饱和）
用牛顿-拉夫森法处理B-H曲线的非线性。残差标准: $||R||/||R_0|| < 10^{-4}$ 是通用的。
频域分析
通过时间谐波假设归结为稳态问题。需要复数运算，但宽带特性需通过时域分析获取。
时域的时间步长
需要最高频率成分的1/20以下的时间步长。隐式时间积分中也可使用更大的步长，但需注意精度。
频域与时域的使用区分
频域分析类似于“调谐到收音机的特定频率”——可以高效计算单一频率下的响应。时域分析类似于“同时录制所有频道”——能再现包含所有频率成分的瞬态现象，但计算成本高。
实践指南
实践指南
请告诉我实际工作中进行静电场分析的流程。
以高压设备的绝缘设计为例来说明。
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以下）
- 使用二阶单元（为了提高电场精度）
电场强度的评价标准是？
三重点（导体·绝缘体·气体交汇点）的电场集中需要特别注意吧。
没错。导体、绝缘体、气体交汇的三重点是电场集中的温床。在COMSOL中，通过网格细化将单元集中在三重点附近是实际工作中的常规做法。
高斯定律用于高压输电线路绝缘子设计的故事
高压输电线路的电压高达数百kV，因此安装在电线杆上的“绝缘子”的绝缘设计非常重要
なった
詳しく
報告