📚 课程目录第4章电磁场理论基础

4. 电磁场理论基础

电磁制动器的多物理场仿真，本质上是在求解一组耦合的偏微分方程：线圈电流产生磁场，运动导体切割磁感线感应出涡流，涡流磁场又反过来阻碍原磁场变化，最终以洛伦兹力的形式表现为制动力。理解这条链条上的每一个理论环节，是正确设置COMSOL物理场接口、边界条件和材料参数的前提。本章从麦克斯韦方程组出发，逐层推导到可直接用于仿真的涡流方程、集肤深度公式与电磁力表达式。

一、麦克斯韦方程组与准静态近似

麦克斯韦方程组是电磁场理论的总纲，其微分形式为：

∇·D = ρv          （高斯定律）
∇·B = 0           （磁通连续性）
∇×E = -∂B/∂t      （法拉第定律）
∇×H = J + ∂D/∂t   （安培-麦克斯韦定律）

在电磁制动器的工作频率（通常为直流或几十至几百赫兹）下，位移电流密度 ∂D/∂t 远小于传导电流密度 J，可以忽略，这就是磁准静态近似（MQS）。此时安培定律简化为 ∇×H = J，电场与磁场通过法拉第定律耦合。COMSOL的“磁场（mf）”接口正是基于这一近似建立的，而“磁场和电场（mef）”接口则保留了完整形式，适用于高频或强电容耦合场景。

磁通连续性方程 ∇·B = 0 是引入磁矢势的直接依据：任何无散场都可以表示为一个矢量场的旋度，因此令 B = ∇×A，该方程自动满足。

二、磁场本构关系

本构关系描述介质对磁场的响应，是连接场量与材料属性的桥梁：

B = μ0 μr H + Br
J = σ (E + v × B)

其中 μ0 = 4π×10⁻⁷ H/m 为真空磁导率，μr 为相对磁导率，Br 为剩余磁通密度（永磁体用），σ 为电导率，v 为导体运动速度。电磁制动器中通常涉及三类材料：空气/非磁性区（μr ≈ 1）、铁磁材料（如低碳钢，μr 可达 1000～5000，但存在饱和）、永磁体（如钕铁硼，μr ≈ 1.05，Br ≈ 1.2 T）。

常见误区：把 μr 当作常数输入。铁磁材料在制动器气隙附近往往已进入非线性区，必须使用 B-H 曲线（磁化曲线）并勾选“磁饱和”选项，否则计算出的制动力会明显偏大。COMSOL中可通过“磁场”接口的“磁通密度模”后处理检查是否超过饱和值（低碳钢约 1.6～2.0 T）。

三、磁矢势与磁标势

直接求解 B 和 H 需要处理多个分量和约束条件，引入势函数可以降维并自动满足约束。两种势函数各有适用场景：

势函数	定义	控制方程	适用场景
磁矢势 A	B = ∇×A	∇×(μ⁻¹∇×A) + σ∂A/∂t = Js	含涡流、含电流源、二维/三维通用
磁标势 Vm	H = -∇Vm	∇·(μ∇Vm) = 0	无电流区域、永磁体与铁芯的静态磁场

磁矢势 A 是涡流仿真的核心变量，因为感应电场 E = -∂A/∂t 直接由它给出，涡流密度 J = σE 也随之确定。在二维轴对称或平面模型中，A 只有一个分量 Az，计算量大幅降低。磁标势 Vm 则常用于永磁体-铁芯的静态预分析，但要求区域内无传导电流，否则 ∇×H = J ≠ 0 会破坏标势的存在条件。

实操建议：在COMSOL中建立电磁制动器模型时，可先用“磁场，无电流（mfnc）”接口配合磁标势快速估算气隙磁密分布，确认磁路设计合理后，再切换到“磁场（mf）”接口加入涡流与运动效应，这样能显著缩短调试时间。

四、涡流效应与集肤深度

当导体（如制动盘）在磁场中运动或处于交变磁场中时，磁通变化会在导体内感应出闭合电流，即涡流。涡流产生的磁场方向总是阻碍原磁通的变化（楞次定律），宏观上表现为制动力。涡流方程可由法拉第定律与安培定律联立得到：

∇×(μ⁻¹∇×A) + σ(∂A/∂t + v·∇A) = Js

其中 v·∇A 项代表导体运动引起的对流效应，在COMSOL中通过“平移运动”或“旋转机械”功能引入。涡流并非均匀分布，而是集中在导体表层，其衰减特征由集肤深度 δ 描述：

δ = √(2 / (ω μ σ)) = √(1 / (π f μ σ))

以铜制动盘为例，σ = 5.8×10⁷ S/m，μ = μ0，当 f = 50 Hz 时，δ ≈ 9.3 mm；若频率升至 500 Hz，δ 降至约 2.9 mm。这意味着高频下涡流几乎只存在于薄表层，网格划分必须在表层加密，否则会严重低估涡流损耗与制动力。

五、电磁力计算与仿真验证

制动力的体积力密度由洛伦兹力给出 f = J × B，总制动力为体积分 F = ∫(J × B)dV。在COMSOL中，可通过“磁场”接口的“力计算”功能直接得到导体域上的总力，也可用“体积分”算子对 J×B 后处理。对于旋转制动盘，还需将力乘以半径得到制动力矩。

一个典型的验证流程如下：

建立二维轴对称模型，设置线圈电流密度 Js 与制动盘电导率 σ。
使用“旋转机械”功能设定制动盘转速 ω，或使用“平移运动”设定线速度 v。
在制动盘表层按 δ/3 的尺寸加密网格，至少布置 3～5 层单元。
求解后提取 J×B 体积分，与解析公式 F ≈ σ v B² A δ 对比，误差应控制在 10% 以内。
扫描转速或频率，观察制动力随速度的变化曲线，确认是否出现饱和趋势。

关键结论：集肤深度决定了网格密度与求解精度，涡流方程决定了物理场耦合方式，而洛伦兹力体积分则是制动器性能的最终输出。三者环环相扣，任何一环设置不当都会导致仿真结果偏离工程实际。

六、常见误区复盘

忽略磁饱和：铁芯 μr 取常数导致气隙磁密虚高，制动力偏大。
网格未按集肤深度加密：涡流被“抹平”，制动力偏小且不收敛。
误用磁标势：在含电流区域使用 Vm 会得到错误结果，应改用磁矢势 A。
遗漏运动项：静止模型无法体现速度对涡流的影响，必须加入 v×B 或旋转机械耦合。
单位不一致：σ 用 S/m、μ 用 H/m、ω 用 rad/s，混用 Hz 与 rad/s 会导致 δ 计算错误。

掌握以上理论后，下一阶段的建模工作就有了明确的物理依据：知道为什么选“磁场”接口、为什么设置运动、为什么在表层加密网格，以及如何判断结果是否可信。