電位分布解析
理论与物理
概述
老师，电势和电场有什么关系呢？
电势 $\phi$ 是一个标量场，电场是其负梯度。
因为可以将3个分量的矢量场归结为1个标量场，所以在FEM中以求解电势为标准方法。
未知数减少到三分之一，这影响很大啊。
电势满足的方程是泊松方程。
在无电荷区域则变为拉普拉斯方程 $\nabla^2 \phi = 0$。
边界条件
电势分析的边界条件有哪些种类？
有三种基本的边界条件。
边界条件 数学表达 物理意义 Dirichlet $\phi = V_0$ 电势固定（导体表面） Neumann $\partial\phi/\partial n = -\sigma_s/\varepsilon$ 法向电场指定 对称边界 $\partial\phi/\partial n = 0$ 与电力线平行的面
导体是等势体，所以 $\phi = \text{const}$。在COMSOL中，可以直观地设置为“Electric Potential”和“Ground”边界条件。
静电能
如何从电势求出静电能呢？
静电能通过对电场的平方进行体积积分求得。
$$ W_e = \frac{1}{2}\int_\Omega \varepsilon |\nabla\phi|^2\,d\Omega $$
这与电容器的能量 $W = \frac{1}{2}CV^2$ 一致。在COMSOL中可以使用“Volume Integration”直接计算。
Coffee Break 闲谈
“电势是标量，为什么比矢量的电场更方便？”
电场 $\mathbf{E}$ 是矢量量（具有x,y,z分量），而电势 $\phi$ 是标量量（只有一个数字）。在FEM求解时，求解标量可以将自由度减少到三分之一，从而大幅降低计算量。“直接求解电场”不如“先求解电势，然后用 $\mathbf{E} = -\nabla\phi$ 计算电场”来得高效。这就是电势分布分析成为CAE标准方法的原因。另一方面，在电势无法定义（积分值随路径变化）的时变电磁场中，会使用“矢量势”来代替电势。选择标量还是矢量决定了计算成本——这是作为CAE工程师需要重视的视角。
各项的物理意义
- 电场项 $\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假设（截面分析）：电流方向均匀且可忽略端部效应时有效
- 各向同性假设：对于各向异性材料（如硅钢板的轧制方向等）需要定义方向特性
- 不适用的案例：等离子体（电离气体）、超导体、非线性光学材料需要额外的本构关系
电势分析的边界条件有哪些种类？
有三种基本的边界条件。
| 边界条件 | 数学表达 | 物理意义 |
|---|---|---|
| Dirichlet | $\phi = V_0$ | 电势固定（导体表面） |
| Neumann | $\partial\phi/\partial n = -\sigma_s/\varepsilon$ | 法向电场指定 |
| 对称边界 | $\partial\phi/\partial n = 0$ | 与电力线平行的面 |
导体是等势体，所以 $\phi = \text{const}$。在COMSOL中，可以直观地设置为“Electric Potential”和“Ground”边界条件。
如何从电势求出静电能呢？
静电能通过对电场的平方进行体积积分求得。
这与电容器的能量 $W = \frac{1}{2}CV^2$ 一致。在COMSOL中可以使用“Volume Integration”直接计算。
“电势是标量，为什么比矢量的电场更方便？”
电场 $\mathbf{E}$ 是矢量量（具有x,y,z分量），而电势 $\phi$ 是标量量（只有一个数字）。在FEM求解时，求解标量可以将自由度减少到三分之一，从而大幅降低计算量。“直接求解电场”不如“先求解电势，然后用 $\mathbf{E} = -\nabla\phi$ 计算电场”来得高效。这就是电势分布分析成为CAE标准方法的原因。另一方面，在电势无法定义（积分值随路径变化）的时变电磁场中，会使用“矢量势”来代替电势。选择标量还是矢量决定了计算成本——这是作为CAE工程师需要重视的视角。
各项的物理意义
- 电场项 $\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公式化。
将泊松方程的弱形式离散化。用形状函数 $N_i$ 近似电势，则单元刚度矩阵为
因为是正定对称矩阵，所以CG法最优。预处理使用AMG有效。COMSOL中会自动选择MUMPS和AMG。
从电势到电场的精度提升
有什么提高电场精度的技巧吗？
电势连续，但 $\mathbf{E} = -\nabla\phi$ 在单元边界处可能不连续。
- SPR法: 基于超收敛点处的梯度值，以面片为单位重构电场
- 二次单元: 若电势为二次，则电场呈一次变化，精度提高
- 平滑处理: 在COMSOL中可作为后处理选项使用
使用二次单元是实务的最低要求呢。
没错。一次单元中电场在单元内为恒定值，无法正确捕捉电场的集中。Ansys Maxwell的自适应网格也基于二次单元运行。
输电线路绝缘子设计者“均等化”电位分布的巧思
高压输电线路的悬式绝缘子（多个串联使用的碗状物），实际上串联的每个绝缘子上的电压并不是均匀分配的。靠近电线侧的绝缘子电场更强，会发生“电压分配不均”。通过CAE的电位分布分析确认这一点，并通过优化屏蔽环或金具的形状、配置来均等化电位分布。过去曾使用“电解槽法”实验来确认——将缩尺模型放入装满水的槽中测量电极间的电位，这是一种巧妙利用拉普拉斯方程类比性的方法。现在虽然FEM一键搞定，但了解这些实验的创意巧思，对分析会更有帮助。
边单元（Nedelec单元）
专用于电磁场分析的单元。自动保证切向分量的连续性，排除伪模式。是3D高频分析的标准。
节点单元
用于标量势公式化。在静磁场的标量势法或静电场分析中有效。
FEM vs BEM（边界元法）
FEM: 对应非线性材料、非均匀介质。BEM: 可自然处理无限域（开域问题）。混合FEM-BEM也有效。
非线性收敛（磁饱和）
用牛顿-拉夫森法处理B-H曲线的非线性。残差标准: $||R||/||R_0|| < 10^{-4}$ 为一般标准。
频域分析
通过时间谐波假设归结为稳态问题。需要进行复数运算，但宽带特性需通过时域分析获取。
时域的时间步长
需要最高频率成分的1/20以下的时间步长。隐式时间积分中可使用更大的步长，但需注意精度。
频域与时域的使用区分
频域分析类似于“将收音机调到特定频率”——可以高效计算单一频率下的响应。时域分析类似于“同时录制所有频道”——可以再现包含所有频率成分的瞬态现象，但计算成本高。
实践指南
实践指南
请告诉我电位分布分析的具体应用实例。
高压套管示例
请展示具体步骤。
变压器套管的电位分析中，
1. 创建2D轴对称模型（旋转对称结构）
2. 芯线设为 $\phi = 275\text{kV}/\sqrt{3}$，法兰设为 $\phi = 0$
3. 对电容锥设置浮动电位
4. 设置各介电层的 $\varepsilon_r$
等位线密集的地方电场越强。锥端部和法兰附近是关键部位。
和等高线地图的读法一样呢。
正是如此。优化电容锥的片数和位置，使等位线间隔均匀。这是绝缘设计的核心。
现场“不会触电”多亏了等电位作业
进行带电作业（高压输电线带电状态下作业）的电力作业人员，通过“等电位作业”的方法来防止触电。作业者穿着导电服接触电线，使身体与电线电位相同，消除身体与电线之间的电位差，从而避免电流流过身体。这正是电位分布分析的实际应用，利用了“在等位面上就不会受到电场力”的原理。在CAE中仔细观察电位云图，会逐渐培养出“沿着这个等位面移动就是安全路径”的感觉。分析结果直接关系到实际安全保障，这正是静电场的魅力所在。
分析流程的比喻
电机的电磁场分析类似于“给吉他调音”。调整弦的粗细（线圈匝数）和琴桥位置（磁铁配置），以引出最美的音色（高效的扭矩特性）。改变一个参数，整体的平衡就会改变——因此参数化研究很重要。
初学者容易陷入的误区
“空气区域？为什么要用网格划分空气？”——这是几乎所有初次接触电磁场分析的人都会产生的疑问。答案是“磁力线也会扩散到铁心之外”。
なった
詳しく
報告