麦克斯韦应力张量 — CAE术语解说

问得好。能量密度是一个标量，但麦克斯韦应力张量是一个描述"力流"的二阶张量。具体来说，它表示电磁场对物质施加的单位面积力（应力）。例如，真空中的张量为

这直接关系到网格质量和计算成本。例如，当计算旋转电机转子上的磁气扭矩时，可以使用仅依赖于转子内磁通密度分布的虚拟功法进行体积积分。但如果使用麦克斯韦应力张量，我们只需要"气隙"这一个面（实际上是一个壳层网格）上的数据来计算扭矩。这样就不需要对整个体积进行高精度网格划分，特别是可以避免边缘和角部奇异点的影响，这在实际工程中具有巨大意义。

这是从电磁场的动量守恒定律必然导出的形式。减法项

最关键的是在"具有均匀且足够细致网格的区域"设置积分面。例如，在Ansys Maxwell中计算永磁电机的扭矩时，应定义一个通过转子外周和定子内周中点的圆筒面作为积分面。这个面必须在气隙中磁通密度平稳变化的区域内。如果将面设置在铁心表面附近，透磁率不连续点处的数值误差会直接传入结果。

根据经验法则，为了捕捉磁通密度的变化，每个局部波长（或极对数的一个周期）的空间周期应至少配置10～20个网格单元。例如，对于4极电机，机械角180度是一个周期，所以积分面的圆周方向至少需要36～72个分割。此外，通过逐步细化网格（h-细化）来计算力或扭矩值是否收敛，这个"网格收敛性分析"在实际工作中是必须进行的步骤。

这是数值误差的主要来源之一。许多FEM求解器（如JMAG）会对积分面上的"高斯积分点"处的磁感应强度进行插值，并在该点计算张量进行积分。为了避免边界处的不连续问题，软件内部会采用一致地仅使用面"单侧"（通常是受力物体外侧）的场值的算法。用户需要注意的是确认积分面"法线方向"在软件内部的定义。如果方向相反，力的符号会反转。

这是一个普遍而关键的问题。一般的最佳实践方法如下：1) 首先，将虚拟功法（能量微分法）作为"基准值"。这是一种体积积分，基于整体能量变化，相对不易受到局部场数值误差的影响。2) 通过改变积分面位置和网格多次计算麦克斯韦应力法的结果，评估其变化范围。3) 如果两种方法的差异在目标精度范围内（例如1%以内），则认为该模型和网格是可信的。Ansys文档也推荐使用这两种方法进行相互验证。

有几个方法。第一是"力的平衡"。对于不含电流源或永磁体的空心区域周围的封闭曲面进行麦克斯韦应力张量面积分，理论上全力和全扭矩应该为零。计算这个值并检查它在数值上是否接近零。第二是与已知解析解的比较。例如，平行板电容器极板间的引力或平行电流线间的力都可以用简单公式表示。第三是检查输出的值是否物理合理。例如，对于直径约50mm的小型无刷直流电机，如果麦克斯韦应力法计算出的扭矩达到几十牛米，显然是单位系统设置错误。

需要区分时间平均力和瞬时力。例如，计算交流驱动感应电机的平均扭矩时，在时域过渡分析中计算数个周期后，对1个周期的瞬时扭矩进行时间平均。此时麦克斯韦应力张量也从瞬时电磁场计算，因此积分面在整个时间内需要固定（变形物体情况除外）。此外，当高次谐波分量较多时，如果时间步长不足，会产生混叠。时间步长的经验法则是，最高频率的至少一个周期应分为20步以上。

完全正确。主要有三点差异。1) **场的插值方法**：Ansys Maxwell从单元边定义的磁矢量势直接在积分点计算B。JMAG多使用在节点处插值的B。COMSOL允许用户选择"单元内评估点"还是"节点平均化"。2) **积分面的处理**：Maxwell自动生成"围绕对象的封闭曲面"，而JMAG需要用户明确定义"力计算面"。3) **张量的形式**：COMSOL在物质内（

不一定。强磁性体表面的力由于透磁率极高、磁场变化急剧，无论哪种方法数值误差都会很大。实际工作中，不过分关注铁心表面的力，而是关注通过气隙传递的整体力和扭矩。这种情况下任何软件的麦克斯韦应力法（配置恰当的话）都很有效。如果要评估铁心内部的磁致伸缩或磁力引起的变形和应力，通常使用COMSOL的多物理场耦合功能，或Ansys的Mechanical联合分析，将电磁解析的体积力密度传递给结构分析。

基本算法已经实现，但与商用软件的区别在于"验证支撑"和"用户支持"。FEMM专门用于二维静磁场分析，可以定义圆弧积分路径计算扭矩，但不支持三维和瞬态现象。Elmer是开源的，灵活性高，但积分面设置和结果后处理都需要用户自己用脚本编写，初学者难度大。商用软件（年许可证数百万日元）经过汽车和电机制造商的严格验证，特别是在转子扭矩脉动等计算中，其精度和再现性已为客户价值买单。

不仅仅是这个原因。第一个怀疑对象应该是"积分面离材料边界太近"。在透磁率不连续变化的界面（铁心表面），磁感应强度的法线分量是连续的，但切线分量是不连续的，数值解在那里也会急剧变化。如果积分面靠近这样的区域，仅移动一个网格单元的距离就会使场的值发生很大变化，导致积分结果不稳定。解决办法是把积分面至少设置在距气隙宽度20～30%的地方，即在磁场平稳变化的区域中央附近。

大多数是数值误差，特别是"网格通过噪声"。当转子旋转时，积分面保持固定，但通过它的转子侧网格单元的形状和配置在每个时间步都在变化。此时，网格界面处的场插值会产生不连续，导致力计算出现波动。缓解这个问题的方法有：1) 旋转部分和固定部分的网格尽可能细致且均匀，2) 在Ansys Maxwell中正确应用"Band"区域之类的滑动网格法，3) 对结果进行一定时间窗口（例如电气角一周期）的移动平均来平滑噪声。

首先检查"求解器残差"。如果静磁场分析还未完全收敛，场分布本身就含有误差。对于Ansys，检查"Residual"值是否低于设定阈值（默认1e-4等）。其次检查"模型对称性"。即使几何上完全对称，如果网格不对称，力平衡也会被破坏。特别是使用自动网格生成时容易出现。最后检查"积分面是否完全闭合"。在三维中，如果面有微小间隙或重叠，高斯散度定理就不成立，会产生误差。在完成这些检查后仍无解决，那就是该FEM离散化本身固有的数值误差水平，需要与虚拟功法的结果比较来判断是否在允许范围内。