4. 磁路等效模型建立方法

电磁执行器的磁路设计，本质上是在有限空间内安排铁芯、永磁体、气隙和线圈，使磁通尽可能高效地穿过工作气隙并产生所需电磁力。直接对三维几何做场求解虽然精确，但计算量大、参数扫描慢，难以在方案阶段快速迭代。磁路等效模型把连续分布的磁场离散成若干段“磁阻”，用类似电路的方式求解磁通分布，是执行器设计中性价比最高的分析手段。本节围绕“如何把实际结构转化为磁阻网络”展开，重点讲清气隙、铁芯、永磁体和漏磁路径的磁阻计算，以及节点法与回路法的求解思路。

一、从磁场到磁路：等效的前提与约定

磁路法的理论依据是安培环路定律与磁通连续性。当铁磁材料磁导率远高于周围介质、且磁通主要沿截面近似均匀的路径闭合时，可以把磁场问题近似为“磁通在磁阻网络中流动”的问题。此时定义磁动势 F（单位 A，安匝）与磁阻 R_m（单位 H^-1），磁路欧姆定律写作：

F = N·I = Φ · R_m
R_m = l / (μ · A)

其中 l 为磁路长度（m），A 为有效截面积（m²），μ = μ_r·μ₀ 为材料磁导率，μ₀=4π×10^-7 H/m。需要特别强调方向约定：本文统一取线圈电流按右手螺旋定则产生的磁通方向为正，永磁体以磁化方向从 S 极指向 N 极（内部）为正，所有磁阻压降沿磁通方向取正。参照系一旦确定，全文与图示保持一致，避免符号混乱。

等效建模的核心判断是：哪些路径必须保留，哪些可以合并。通常把执行器拆成“主磁路”和“漏磁路”两部分。主磁路承担有效工作磁通，漏磁路则包括气隙边缘扩散、铁芯侧面之间的旁路、永磁体端部自闭合等。漏磁并非可以忽略，在气隙较大或永磁体较薄时，漏磁系数可达 1.2～1.5，直接影响推力估算精度。

二、各段磁阻的计算方法

气隙磁阻是最关键也最容易被低估的一段。理想情况下 R_g=l_g/(μ₀A_g)，但实际气隙磁通会向外扩散，等效面积大于几何面积。工程上常用经验修正：对矩形极面，取 A_g=(a+l_g)(b+l_g)，即每边外扩约半个气隙长度；对圆柱极面，取 A_g=π(r+l_g/2)²。若气隙 l_g=0.5 mm、极面 20 mm×20 mm，则几何面积 400 mm²，修正后约 (20.5)²≈420 mm²，磁阻下降约 5%，在精密推力计算中不可忽略。

铁芯磁阻需考虑材料非线性。硅钢或纯铁在低磁密时 μ_r 可达数千，但接近饱和（B≈1.5～1.8 T）时急剧下降。因此铁芯磁阻不能简单用常数计算，而应结合 B-H 曲线迭代：先假设磁通，查得 μ_r，算出磁阻与压降，再回代修正。若铁芯工作磁密控制在 1.2 T 以下，可近似取 μ_r=2000～4000 做初算。铁芯分段时，每段用其平均长度和平均截面积，转角处可适当增加等效长度。

永磁体磁阻的处理方式与普通磁阻不同。永磁体既是磁通源又是磁阻，常用“磁动势源 + 内磁阻”串联模型：

F_m = H_c · l_m
R_m = l_m / (μ_r · μ_0 · A_m)

其中 H_c 为矫顽力，l_m、A_m 为磁化方向长度与截面积，μ_r 为回复磁导率（钕铁硼约 1.05）。注意永磁体的工作点由外磁路负载线决定，不能孤立计算。

漏磁路径通常按几何近似估算。例如两铁芯侧面间距 d、重叠长度 L，可把漏磁路径看作一系列并联微元，用 R_leak≈d/(μ₀·L·w) 估算，w 为等效宽度。更精细时可用分段积分或经验漏磁系数修正主磁通。

三、节点法与回路法求解磁通分布

把各段磁阻连接成网络后，求解方法与电路完全对应。节点法以磁通连续性为依据，对每个节点列 ΣΦ=0，未知量为节点磁位（磁动势），适合节点少、支路多的网络；回路法以安培环路定律为依据，对每个独立回路列 ΣF=ΣΦR_m，未知量为回路磁通，适合回路清晰、串联为主的磁路。实际执行器磁路往往既有串联主路又有并联漏磁，两种方法可混用：先按回路法求主磁通，再用节点法分配漏磁。

下面这张图给出一个典型电磁执行器的磁阻网络等效结构，包含永磁体支路、铁芯段、工作气隙和两条主要漏磁路径，并标注了节点与回路编号，便于对照后续求解。

以图中网络为例，节点法可对节点 A、B、C 列磁通连续性方程，节点 D 取为参考零磁位。回路法则取回路 I（主磁路）和回路 II（永磁体—漏磁旁路）列安培环路方程。两种方法得到的未知量不同，但结果应一致，可用于交叉校验。

四、一个可复算的示例

设某执行器永磁体为钕铁硼，l_m=8 mm，A_m=300 mm²，H_c=890 kA/m，μ_r=1.05；工作气隙 l_g=0.5 mm，极面 20 mm×20 mm；铁芯总长 120 mm，平均截面积 350 mm²，初取 μ_r=3000。先算各段磁阻：

R_pm = 0.008 / (1.05 × 4π×10⁻⁷ × 300×10⁻⁶) ≈ 2.02×10⁷ H⁻¹
R_g  = 0.0005 / (4π×10⁻⁷ × 420×10⁻⁶) ≈ 9.47×10⁵ H⁻¹
R_c  = 0.12 / (3000 × 4π×10⁻⁷ × 350×10⁻⁶) ≈ 9.09×10⁴ H⁻¹
F_m  = 890000 × 0.008 = 7120 A

忽略漏磁时主磁通 Φ = F_m/(R_pm+R_g+R_c) ≈ 7120/(2.02×10⁷+0.947×10⁶+0.091×10⁶) ≈ 3.36×10^-4 Wb，气隙磁密 B_g=Φ/A_g≈0.80 T。若计入漏磁系数 1.25，则有效气隙磁通降至约 2.69×10^-4 Wb，B_g≈0.64 T。这个差异说明漏磁路径必须显式建模，否则推力会被高估约 25%。

五、建模流程与常见误差来源

完整的等效建模可按以下顺序推进：先按几何划分主磁路区段，逐段计算磁阻并标注节点；再识别主要漏磁路径，用经验公式或分段积分估算漏磁磁阻；然后选择节点法或回路法列方程，求解磁通与磁位；最后用铁芯 B-H 曲线校核磁密，若饱和则修正磁阻并迭代。误差主要来自三处：气隙边缘扩散系数取值不当、铁芯非线性未迭代、漏磁路径遗漏或重复计算。工程上建议用有限元对等效模型做一次标定，之后即可用等效网络做快速参数扫描。

掌握磁阻网络的建立与求解后，就能在方案阶段快速判断永磁体尺寸、气隙大小和铁芯截面是否匹配，为后续有限元精细仿真提供合理的初值和边界条件。