磁气矢量势
磁气矢量势的理论基础
矢量势$\mathbf{A}$
老师，磁气矢量势是什么？
磁通密度$\mathbf{B}$的「势」。$\mathbf{B}$的散度为零（$\nabla \cdot \mathbf{B} = 0$），因此必能表示为某个矢量场$\mathbf{A}$的旋度：
这对应于静电场的电势$\phi$吗？
对应关系：
| 静电场 | 静磁场 |
|---|---|
| $\mathbf{E} = -\nabla\phi$ | $\mathbf{B} = \nabla \times \mathbf{A}$ |
| $\phi$：标量 | $\mathbf{A}$：矢量（3个分量） |
| 泊松方程：$\nabla^2\phi = -\rho/\varepsilon$ | $\nabla^2\mathbf{A} = -\mu\mathbf{J}$（库仑规范） |
规范条件
$\mathbf{A}$无法唯一确定（加上任意标量场$\psi$的梯度后$\mathbf{B}$不变）。为了唯一性，需要引入规范条件：
- 库仑规范：$\nabla \cdot \mathbf{A} = 0$ — 静磁场标准
- 洛伦茨规范：$\nabla \cdot \mathbf{A} + \mu\varepsilon \partial\phi/\partial t = 0$ — 动态问题
2D中的$A_z$
在2D问题中（$z$方向均匀），$\mathbf{A} = A_z(x,y) \hat{z}$变成标量问题。磁力线是$A_z$的等高线。这是2D磁场有限元的最大优点。
总结
磁气矢量势——「无物理实体之量」为何支撑有限元？
磁气矢量势A（B=旋度A）被称为电磁场的「影子量」，无法直接测量但在计算中不可或缺。有限元磁场分析以A而非B为未知数，是因为与其处理旋度B=μ₀J的数值困难，不如自动满足散度B=0（磁通连续性）更重要。在A定式化中，散度B=0自动保证为旋度（旋度A）=0。1970年代Oliver C. Zienkiewicz等人将A定式化应用于有限元，如今已成为所有主要磁场有限元求解器的基础。
磁气矢量势的数值计算方法
A法的有限元定式化
弱形式：
$\mathbf{N}_i$：边单元形状函数（Nédélec单元）。
2D定式化
在2D中使用节点单元离散化$A_z$：
与静电场泊松方程形式完全相同。$\varepsilon \to \nu$、$\phi \to A_z$、$\rho \to J_z$。
总结
规范固定——AV定式化中「非唯一性」的数值技术解法
磁气矢量势A在规范变换（A→A+梯度φ）下不变，因此具有「规范自由度」无法唯一确定。有限元求解A需要固定规范，主要方法有「库仑规范（散度A=0）」和「树-余树（Tree-Cotree）规范」。库仑规范通过罚函数法或拉格朗日乘数法实现，被许多商业求解器采用。树-余树规范改善矩阵条件数，提高大规模问题的收敛性。规范固定的选择直接影响计算稳定性和效率，是重要的数值技术。
磁气矢量势的实际应用
实务
几乎所有磁场有限元分析都基于A法（矢量势法）。JMAG、Maxwell、COMSOL、FEMM的内部求解器均采用A法。
检查清单
「A定式化与磁通密度精度差异」——边单元vs节点单元的选择
磁场有限元中离散化A时，选择节点单元（标量自由度）或边单元（Nedelec单元）直接影响精度。节点单元不能保证切线分量连续性，在不同透磁率的界面处B的法向分量产生误差。边单元自然满足切线分量连续性，能抑制界面的数值跳跃。现代商业求解器几乎均采用边单元，但从许可成本和内存考虑，定制实现有时选用节点单元。在永久磁石和高透磁率材料界面处，边单元是精度重要性下的标准选择。
磁气矢量势的软件比较
工具
所有磁场有限元求解器都实现了A法。用户几乎无需特别考虑。FEMM可以直接可视化$A_z$，具有教学价值。
矢量势定式化工具——ANSYS Maxwell vs COMSOL AC/DC
采用磁气矢量势A（或AV）定式化的主要工具包括ANSYS Maxwell（A-φ定式化）、COMSOL AC/DC Module（A定式化，边单元标准配备）、Elmer FEM（开源，树-余树规范已实现）。Maxwell在大规模3D问题的AMG求解器（代数多重网格）优化方面进展显著。COMSOL通过LiveLink与MATLAB连接，便于自定义规范和材料模型的实现。OpenFOAM的electromagneticsics模块也采用A定式化，在流体-电磁耦合问题研究中应用广泛。
磁气矢量势的先进研究
先进
AV定式化的扩展——时间谐波分析与复矢量势
交流磁场分析中A被视为复量的「时间谐波（频域）定式化」广泛应用。以A=Re(Â exp(jωt))形式同时求解实部和虚部，定常周期解可一次完成。但含非线性BH曲线的铁芯无法用复A严格处理，需用「等效损耗电阻」近似或「时间步进（FDTD式）」。处理多频成分（逆变器高次谐波）时需多频率叠加计算或非线性时间步进，计算成本与频率数成正比增长。
磁气矢量势的故障排除
故障
「规范未固定导致有限元不收敛」——典型初学者陷阱
磁气矢量势A定式化不固定规范时，解无法唯一确定，反复求解器收敛缓慢或不收敛。此问题表现为「奇异矩阵（非正则系统）」，求解器报错为「不定矩阵」或「检出零特征值」。对策：①检查求解器规范设置（多数具「自动规范」设置）；②在边界条件中至少一处固定A的法向或切向分量。商业求解器自动处理规范固定，但定制有限元代码需明确实现规范，是容易被忽视的重要点。
价值
更详细
需报告