CAE中的规范条件是什么？

规范条件是在电磁场解析等中，对标量势和矢量势表示所产生的数学自由度进行约束的附加条件。在不改变物理场（电场或磁场）的情况下，使数值解唯一确定，防止求解器发散或产生非物理振荡。

库仑规范的优点是什么？

库仑规范（∇・A=0）强制矢量势是无散度的。这使得方程的系数矩阵具有对角优势，大幅提高了迭代求解器（如ICCG法）的收敛性，使数值计算更加稳定。

洛伦兹规范在哪些解析中使用？

洛伦兹规范特别适用于包含时间变化的瞬态解析和电磁波辐射问题。应用此条件后，势的波动方程形式变为对称，相对论共变性也得以保持，使解析更加高效。

如果不设置规范条件，CAE解析会出现什么问题？

不设置规范条件会导致势不能唯一确定，数值解可能发散，产生非物理数值噪声（振荡），进而导致计算不稳定，无法获得准确结果。

规范条件 — CAE术语解释

分类：术语集 | 2026-01-15

CAE visualization for gauge condition - technical simulation diagram

规范条件

🙋

老师，我在电磁场解析教科书中看到"规范条件"这个词，但为什么要特意施加一个条件呢？

规范条件的理论基础

规范条件的物理意义

🙋

我经常看到"规范条件"这个术语。到底在"规范"什么呢？

🎓

好问题。这是电磁场解析和流体解析中的基础概念。例如，在电磁场的麦克斯韦方程中，标量势φ和矢量势A并不能唯一确定，存在数学上的"自由度"。这个多余的自由度需要用"规范条件"来约束。物理场（E或B）保持不变的情况下，通过固定势的取值表示形式来实现的。

🙋

如果有自由度，具体会产生什么不利影响？数值计算会变得不稳定吗？

🎓

完全正确。数值上解不能唯一确定，求解器会发散，或产生非物理振荡（数值噪声）。例如，库仑规范

$$\nabla \cdot \mathbf{A} = 0$$

强制矢量势A是无散度的（旋度场）。这样做，泊松方程的系数矩阵就具有对角优势，迭代求解器如ICCG法的收敛性大幅改善。在实务中，频域解析经常使用。

🙋

还有其他的规范条件吗？我听说过洛伦兹规范。

🎓

是的。另一个代表是洛伦兹规范。其表达式为

$$\nabla \cdot \mathbf{A} + \mu \epsilon \frac{\partial \phi}{\partial t} = 0$$

。这在包含时间变化的问题中特别有用，特别是瞬态解析和电磁波辐射问题。应用此条件后，势对应的波动方程形式变为对称，相对论协变性也得到保持。Ansys Maxwell的瞬态解析默认设置就是这个洛伦兹规范。

规范条件的数值计算手法

有限元法的导入方法

🙋

规范条件是如何融入有限元法的矩阵方程中的？只是简单的约束条件吗？

🎓

集成方式有两种主要方法："罚函数法"和"拉格朗日乘数法"。罚函数法是在能量函数中添加一个"惩罚"项，以处理不满足规范条件的情况。例如，库仑规范会添加

$$\frac{1}{2\mu_0} \int_\Omega (\nabla \cdot \mathbf{A})^2 d\Omega$$

这样的项。关键是罚函数系数的设置，系数过大会导致数值条件恶化。

🙋

拉格朗日乘数法呢？那好像更严格一些。

🎓

严格满足约束，但代价是未知数增加。对约束条件

$$\nabla \cdot \mathbf{A} = 0$$

引入拉格朗日乘数λ，求解扩展系统。最终的矩阵具有如下分块结构：

$$ \begin{bmatrix} \mathbf{K} & \mathbf{G}^T \\ \mathbf{G} & \mathbf{0} \end{bmatrix} \begin{bmatrix} \mathbf{A} \\ \boldsymbol{\lambda} \end{bmatrix} = \begin{bmatrix} \mathbf{b} \\ \mathbf{0} \end{bmatrix} $$

其中，G是散度算子的离散化。这类系统称为"鞍点问题"，需要特殊的求解器。

🙋

在实际软件中，用户能调整这些系数和方法吗？

🎓

高级软件可以。例如，COMSOL Multiphysics的AC/DC模块中，在"磁矢量势"物理场接口设置中，可以从"罚函数"或"拉格朗日乘数"中选择规范固定方法。罚函数系数的默认值通常与材料透磁率μ的倒数成比例（1/μ），但收敛困难时需要在1e-6到1e6范围内手动调整。

规范条件的实务应用

仿真设置指南

🙋

实际做解析时，在什么情况下需要特别注意规范条件？总是需要吗？

🎓

不一定总是需要意识到。许多商业软件在后台已经自动应用了合适的规范。但在两种主要情况下，用户需要明确改变设置。第一，求解器不收敛，或结果出现非物理数值噪声（例如，磁通密度在区域内随机微小振荡）。第二，涉及开放区域（无限远）的辐射问题，或没有导体的静磁场解析。

🙋

为什么开放区域问题需要特别注意？

🎓

在开放区域，无法施加一个自然的边界条件（狄利克雷条件），即势值在无限远处为零。此时，规范条件在确定无限远处势的行为中起着重要作用。例如，仅有外部电流源的静磁场问题中，如果不固定规范，A会在无限远处发散，导致数值计算失败。因此需要用"树网格"进行规范固定，或设置特定点的"点规范"。

🙋

我听说过"树规范"。那是做什么的？

🎓

树规范（或树接地）是通过断开网格的边来消除环路，从而唯一定义矢量势A的技术。在整个网格中构造一个"树"结构，在作为"主干"的边上强制设置A的切向分量为0。这样保持了磁通，但消除了势的自由度。JMAG和Flux等软件在静磁场解析设置中提供了"树规范创建"选项。对于复杂的3D模型特别有效，但需要注意树的选择方式可能对结果产生微小影响。

规范条件的软件对比

各软件的实现与默认值

🙋

主流CAE软件在规范条件的处理上有差异吗？

🎓

差异很大。Ansys Maxwell（3D）的静磁场求解器默认使用"迭代求解器"，其内部自动进行规范处理。用户可选择"自动""点规范"或"树规范"等。而在瞬态和频率响应解析中，如前所述，洛伦兹规范被内置为默认。

🙋

Abaqus的电磁场解析怎么处理？

🎓

Abaqus/Standard的电磁模块主要处理涡流解析和静磁场解析。采用A-φ法（电磁矢量势和电气标量势），规范固定通过"罚函数法"隐式进行。手册中没有详细说明，但存在"稳定化参数"调整选项来改善收敛性。该值由单元大小和材料常数自动计算，但发散时需在1e-12到1e-3范围内手动修改。

🙋

我听说开源的Elmer和GetDP允许更细致的控制。

🎓

完全同意。Elmer FEM的"MagnetoDynamics"求解器中，可在求解器设置文件（SIF）中明确写入"Gauge Fixing = Logical True"，并进一步选择"Gauge Fixing Method"为"coulomb"或"lorentz"。GetDP要求用户自己编写弱形式，因此需要将规范条件作为附加方程显式定义。例如库仑规范时，需为所有测试函数q添加

$$\int_\Omega \mathbf{A} \cdot \nabla q \, d\Omega = 0$$

这样的式子。自由度高，但需要深厚的理论理解。

规范条件的故障排除

常见错误及对策

🙋

与规范相关的具体错误信息有哪些？

🎓

取决于软件。Ansys Maxwell中，当规范不足导致系数矩阵奇异时，会出现"The matrix solver failed due to a singular matrix. Please check the model and boundary conditions."这样的消息。COMSOL中，瞬态解析开始时可能出现"Failed to find consistent initial values."的错误，这可能表示初始条件与规范的不一致。

🙋

被告知"奇异矩阵"时，首先应该检查什么？

🎓

第一步检查所有边界条件是否都已设置。特别是磁矢量势A是否在某些区域没有施加任何狄利克雷条件（固定值）。如果所有边界都只有自然边界条件（诺伊曼条件），那么规范根本没有被固定，矩阵必然奇异。对策是在物理上合理的点（例如，预期磁通为零的远场或对称面）添加一个"A = 0"的边界条件。如果难以这样做，启用软件的"自动规范固定"选项。

🙋

计算收敛但磁通密度分布随网格细节波动的情况呢？

🎓

这是"数值规范污染"的典型症状。规范固定可能过弱（罚函数系数太小），或网格太粗，无法表达高频分量。首先尝试细化网格。如果仍未改善，在Ansys Maxwell中，提高"Solver Settings"的"Advanced"标签中"Gauge Penalty"系数，从默认的1.0改为10.0或100.0。在COMSOL中，将罚函数法系数改为"手动"，设置为材料透磁率倒数的100到1000倍（例如空气中为1e6到1e7），然后重新计算。

🙋

瞬态解析中，每个时间步求解器收敛性变差，这也是规范的问题吗？

🎓

有可能。特别是在使用洛伦兹规范的瞬态解析中，时间离散化方法（隐式还是显式）与规范条件的相容性会成为问题。对策是减小时间步长。或者，某些软件提供"规范重置"选项，该功能在数个步长后强制修正因误差积累而偏离规范条件的势。这能提高长时瞬态解析的稳定性。