3. 电磁场有限元基础与工具链

电磁执行器的动态特性，本质上由“磁场—电流—运动—机械负载”四者耦合决定。要在一台电磁阀、比例电磁铁或旋转执行器上复现吸合时间、电流拐点、涡流滞后与撞针速度，仅靠磁路法或静态有限元往往不够，必须进入瞬态磁场求解。本节围绕 Maxwell、Flux、MagNet 这类主流低频电磁场工具，说明瞬态求解的关键设置、网格与材料非线性处理，以及运动网格和带宽策略，使仿真结果具备工程可信度。

一、瞬态磁场求解的基本方程与工具链定位

瞬态磁场求解的核心是涡流方程（A-φ 或 A-V 形式）：

∇×(1/μ ∇×A) + σ(∂A/∂t + ∇φ) = J_s
∇·(σ(∂A/∂t + ∇φ)) = 0

其中 A 为磁矢位，φ 为电标量位，σ 为电导率，J_s 为源电流密度。与静态求解相比，瞬态求解多了 σ∂A/∂t 项，它同时带来涡流、集肤效应和磁通滞后。工具链上，Ansys Maxwell 适合二维/三维瞬态与运动带（Band）建模，并与 Simplorer/Twin Builder 做场路耦合；Altair Flux 在二维瞬态与旋转机械运动上成熟，参数化和脚本能力强；Infolytica MagNet（现属 Mentor/Siemens）以三维瞬态和运动网格见长。三者思路一致：几何—材料—激励—网格—运动—求解—后处理。

二、求解设置：时间步、激励与耦合

瞬态求解的时间步长直接决定能否捕捉电流拐点和吸合瞬间。经验做法是：在运动发生前用较大步长（如 0.5 ms），在衔铁即将吸合、气隙快速减小的区间加密到 0.05～0.1 ms，并开启自适应时间步。以 24 V 直流电磁阀为例，线圈电阻 12 Ω、电感约 30 mH，时间常数约 2.5 ms，若步长取 1 ms 会明显低估吸合时间。激励可用外电路（电压源+电阻+二极管续流）或直接给定电流波形；前者更接近真实驱动，后者便于快速扫描。三维模型还需设置求解域边界（通常为空气包，边界取磁绝缘或远场），并指定对称面以减少计算量。

三、网格划分与材料非线性

网格策略应“按物理梯度分配”：气隙、衔铁尖端、轭铁拐角处磁通密度梯度大，需局部加密；线圈和远场空气可粗化。二维模型常用三角形网格，气隙内至少 3～5 层单元；三维模型在运动带附近用扫掠或棱柱层。材料非线性方面，硅钢、纯铁、坡莫合金的 B-H 曲线必须输入真实数据，并注意：① 曲线需单调、无交叉；② 高场段可外推但不宜过度；③ 叠片铁心需用叠压系数修正等效 B-H 与电导率，以体现涡流损耗。若忽略非线性，饱和区磁通会被高估，吸力偏大、吸合时间偏短。

四、运动网格与带宽策略

电磁执行器的衔铁会运动，网格必须随动。主流做法是“运动带（Band）+ 气隙层重划”：在衔铁外包裹一层可变形空气带，带内网格随位移拉伸或压缩，带外网格保持静止。带宽应大于行程并留 20%～30% 余量，否则单元过度畸变会导致求解发散。二维模型可用滑移网格（Sliding Mesh），三维模型常用重划网格（Remeshing）或任意拉格朗日—欧拉（ALE）方法。运动速度较高时，还需考虑网格更新频率与时间步的匹配：若一步位移超过气隙单元的 1/3，应减小步长或增加气隙层数。此外，运动带内不宜放置非线性材料或细密结构，以免重划时破坏材料边界。

五、一个可复现的二维瞬态设置示例

以某比例电磁铁二维轴对称模型为例：线圈 500 匝，电阻 8 Ω，驱动电压 24 V，初始气隙 1.2 mm，衔铁质量 0.05 kg，复位弹簧刚度 2 N/mm。设置如下：求解类型 Transient，终止时间 20 ms，初始步长 0.2 ms，最小步长 0.02 ms，最大步长 0.5 ms；运动带包住衔铁，带宽 2 mm；材料为纯铁 B-H 曲线 + 叠压系数 0.96；激励用外电路电压源串联 8 Ω 电阻。运行后观察电流曲线：约 3.5 ms 处出现拐点，对应衔铁起动；约 7 ms 电流达到稳态 3 A，位移到达 1.2 mm。若把步长改为固定 1 ms，吸合时间会被算成约 5 ms，误差明显。这说明时间步与运动网格必须协同设置。

综上，瞬态磁场求解的可信度取决于四件事：时间步是否跟得上运动、网格是否跟得上梯度、材料是否跟得上饱和、运动带是否跟得上位移。把这四点与真实驱动电路、机械负载一起纳入模型，才能为后续电磁—机械动态联合仿真提供可靠的场量输入。