洛伦兹力耦合分析 — 电磁力与结构变形的耦合方法

反而相反，结构工程师更应该了解这个课题。比如MRI超导磁石——起磁时线圈中流过几百安培的电流，在强磁场中受到洛伦兹力，线圈会产生数吨的膨胀力。这个力能否被支承结构承受，必须通过结构分析来评估。

更贴近的例子是变压器短路事故。短路电流上升到通常的10～25倍时，线圈间会瞬间产生数百kN的电磁力。这会导致线圈座屈或绝缘破坏。电力公司高度重视"变压器短路耐力试验"正是出于这个原因。此外还有直线电动机、电磁成形、轨道炮——只要有电流和磁场存在的地方，就必然存在洛伦兹力。

完全正确，那是出发点。包含电场和磁场的一般形式是：

在CAE耦合分析中，这需要扩展为连续体，用体积力密度的形式表达。涉及电流密度 $\mathbf{J}$（A/m²）与磁通密度 $\mathbf{B}$（T）的外积：

完全正确。与重力性质相同，作为体积力进入FEM的右端。结构运动方程式如下：

其中 $\boldsymbol{\sigma}$ 是应力张量，$\mathbf{f}_{\text{em}} = \mathbf{J} \times \mathbf{B}$ 为电磁体积力，$\rho$ 是密度，$\ddot{\mathbf{u}}$ 是加速度。静力分析时右侧为零。简言之，"通过电磁场分析求得 $\mathbf{J}$ 和 $\mathbf{B}$，计算其外积，将该力作为荷载施加到结构分析"——这就是洛伦兹力耦合的核心。

提得好。J × B 只能在有电流的区域定义力。但永久磁铁或铁氧体核这类磁性体没有电流流动，却受到磁场力的作用。此时使用Maxwell应力张量：

在物体表面对该应力张量进行面积分，得到作用在物体上的电磁力合力：

理论上，围绕物体的任意闭合曲面应该给出相同结果，但在数值上会因网格粗度而变化。实务中，"空气间隙积分"（将积分面离开物体表面一定距离）比较稳定。COMSOL默认采用这种方法。

眼光敏锐。磁歪（magnetostriction）是与洛伦兹力不同的电磁-结构耦合机制。磁性体被磁化后，伴随磁区旋转，晶格会发生微小变形。变压器铁心在50Hz/60Hz交流磁场中反复收缩，因此产生唱歌声。

磁歪表现为结构方程中的歪斜张量 $\varepsilon^{\text{ms}}$：

与热膨胀公式 $\sigma = C(\varepsilon - \alpha \Delta T)$ 形式相同，对吧？ 磁歪歪斜通常为 $10^{-6}$～$10^{-5}$ 数量级，电工钢大约 $\lambda_s \approx 5 \times 10^{-6}$。看似很小，但铁心整体累加后，成为变压器噪声和振动的主要来源。

完全正确。实际问题中这些往往叠加出现。比如电动机中，定子绕组受到洛伦兹力，同时铁心产生磁歪振动。做噪声分析时若不同时考虑两者，试验结果会无法匹配。

大多数情况下弱耦合足够。MRI磁石的静态评估、变压器短路耐力分析，结构变形都很小（毫米以下），对电磁场影响可忽略。

但强耦合必要的情况也存在。典型例子是电磁成形。向铝板瞬间施加涡电流，通过与磁场的相互作用实现高速变形的加工工艺。变形改变线圈与被加工材的间距，磁场分布和电流分布都会明显改变。单向法会高估力，从而给出错误的结果。直线电动机的空气间隙变动问题也是同理。

在低频电磁场分析中，从Maxwell方程忽略位移电流的准静态近似出发。用磁矢量势 $\mathbf{A}$ 和标量势 $\phi$ 的A-φ法进行定式化，是标准做法：

其中 $\mu$ 为透磁率，$\sigma$ 为导电率，$\mathbf{J}_s$ 为外部施加电流密度。由此可求出磁通密度 $\mathbf{B} = \nabla \times \mathbf{A}$，导体内诱导电流密度为 $\mathbf{J} = -\sigma(\partial \mathbf{A}/\partial t + \nabla \phi)$。这些的外积就是洛伦兹力。

完全一样的想法。但有重要区别。电磁FEM常用边元素（Nedelec元素）。标准节点元素无法保证磁通密度的切向连续性。与结构的Lagrange元素完全不同，连成时需特别注意元素类型不匹配。

主要有三种方法。各有所长，需根据用途选择：

概念完全相同。将物体虚拟微小位移 $\delta s$ 时，磁场能量变化求出力：

JMAG和Ansys Maxwell等工具会自动计算并输出该虚功法力作为节点力。精度与稳定性平衡很好，在转矩计算等场景已成为事实标准。

这是耦合分析的关键。网格映射有三种方法：

电磁侧计算的合力与结构侧转移后的合力必须一致。这个性质崩坏的话，结构模型就会莫名其妙地产生或消失力量。变压器短路力分析中，转移的力保守性偏离5%以上，应力评价就无法定量使用。所以转移后必须确认合力一致的习惯很重要。

这是过渡耦合的难点。电磁场以表皮效应时常数（毫秒～微秒）变化，结构固有振动周期是毫秒～秒。时间尺度常相差2～3个数量级。

有效对策是亚循环法。电磁场用细时间步 $\Delta t_{\text{em}}$ 求解，结构用粗步 $\Delta t_{\text{str}} = n \cdot \Delta t_{\text{em}}$ 求解。力的转移在结构每步处进行，该区间的电磁力时间平均值转移给结构。变压器短路解析典型是电磁场 $\Delta t = 0.1 \text{ms}$，结构 $\Delta t = 1 \text{ms}$ 这样的组合。

讲解弱耦合（单向）情况的标准流程。大多实务情况都能用这个应对：

常见做法。比如用JMAG做电磁分析，力数据输出给Nastran或Abaqus。JMAG有Nastran格式力数据输出功能。在Ansys阵营则用Ansys Maxwell和Mechanical由Workbench连接最方便。异软件间连携时要注意力文件格式（CSV、UNV、.bdf等）和网格映射方式。

电磁网格最关键的是表皮深度的分辨率。表皮深度 $\delta$ 为：

导体表面深度 $\delta$ 以内集中了63%的电流，所以这个区域必须有最少3层、最好5层要素。比如铜在50Hz的表皮深度约9.4mm，但短路事故初期会有数kHz的高次谐波，$\delta < 1$ mm。

这是结构工程师容易踩坑的地方。电磁场存在于空气中，所以周围空间也要网格化。导体和磁性体周边的空气要离开足够远，不然磁通受到边界扭曲，力计算精度下降。对策是用无限元素（Infinite Element）或PML（完全匹配层）来处理。

列举典型例子：

ITER中，等离子体突然崩溃（失稳）时，毛毯结构在数毫秒内承受数MN电磁力。结构的过渡响应评估缺一不可。Abaqus电磁-结构耦合模块或COMSOL都被用上了。

主要工具比较如下。要点是"电磁场求解质量"和"结构耦合通路成熟度"两方面：

完全正确。COMSOL的最大优势是在一个GUI、一套网格、一个求解器内处理多物理耦合。网格映射的问题完全消除，强耦合也能无缝进行。缺点是大规模问题求解速度不如Ansys或JMAG。

JMAG可输出.bdf文件（Nastran格式）或.csv，供Nastran或Abaqus读取。通过JMAG-RT与MBD仿真（Adams等）连接，可做马达电磁激励力→车身振动→噪声的全链路耦合。日本汽车厂商马达开发中这套手法是主流。

引人关注的几个前沿方向：

有用FEM数据做监督学习的方法，也有PINNs靠支配方程残差无数据学习的。但目前精度还达不到FEM，多用于"设计初期筛选"或"实时监测代理模型"。

正是最前沿课题。电力变压器数字孪生构想中，运行中的负荷电流作输入，电磁力→线圈应力实时推定，累积疲劳监视。完整FEM计算来不及，事先学习电磁力→结构响应映射，用模态展开法或ROM（降阶模型）高速计算。MRI对超导失效检知与结构健全性评估的数字孪生也在研究中，需秒级反应。

耦合特有的问题可归纳为三大类。逐个看下去。

电磁力空间振荡（噪声）是最常遇到的问题。原因几种：

对策：先把电磁网格细化到表皮深度1/5以下。还有噪声的话，用空间滤波（高斯光滑）。COMSOL的"Volume Average"运算子，Maxwell的力光滑选项。但过度光滑会歪曲物理，转移后合力要确认未变。

松弛有效，但先找原因。双向耦合发散主要两种：

Aitken松弛更新式：

$\mathbf{r}_k$ 为第 $k$ 迭代的残差向量。COMSOL的Segregated Solver内装此法。Ansys Workbench需手动设松弛系数，初期保守 $\omega = 0.3$ 开始逐步提高较安全。

常见的汇总如下：

典型案例。耦合分析要"动手前纸上估算一遍力的数量级"。粗估 $F = J \times B \times V$（$V$ 导体体积），对比FEM结果。1个数量级以上偏差就有问题。