4. 直线电机数学模型

直线电机的伺服控制性能,最终取决于控制器对电机内部电磁关系的理解精度。与旋转电机类似,直线电机可以在三相静止的abc坐标系中建立原始方程,也可以通过Park变换与Clarke变换映射到随动子同步运动的dq坐标系中。两种坐标系描述的是同一台电机,但abc坐标系更贴近物理绕组,dq坐标系更适合控制器设计。本节围绕电压方程、磁链方程、推力方程和动态模型四条主线,梳理直线电机数学模型的完整脉络。

一、abc坐标系:从物理绕组出发的原始方程

直线电机的初级铁心开槽并嵌放三相绕组,通入三相对称电流后,在气隙中形成沿行程方向平移的行波磁场。在abc坐标系下,三相绕组的电压方程可写为:

u_a = R_s·i_a + dψ_a/dt
u_b = R_s·i_b + dψ_b/dt
u_c = R_s·i_c + dψ_c/dt

其中 R_s 为每相绕组电阻,ψ_a、ψ_b、ψ_c 为三相磁链。磁链由自感、互感以及永磁体磁链共同构成。对于永磁同步直线电机(PMSLM),三相磁链可表示为:

ψ_a = L_a·i_a + M_ab·i_b + M_ac·i_c + ψ_f·cos(θ_e)
ψ_b = M_ba·i_a + L_b·i_b + M_bc·i_c + ψ_f·cos(θ_e - 120°)
ψ_c = M_ca·i_a + M_cb·i_b + L_c·i_c + ψ_f·cos(θ_e + 120°)

式中 θ_e 为电角度,ψ_f 为永磁体等效磁链幅值。由于直线电机铁心开断、端部效应显著,三相自感和互感并非严格常数,而是随动子位置周期性变化,这也是abc坐标系模型复杂、难以直接用于控制器设计的根本原因。

abc坐标系的价值在于物理直观:它直接对应实际绕组电流与端电压,适合用于有限元仿真验证、故障分析和绕组不对称研究。

二、dq坐标系:解耦后的控制模型

通过Clarke变换将三相量投影到两相静止αβ坐标系,再经Park变换旋转到与动子磁场同步的dq坐标系,可消除时变电感带来的耦合。对于表贴式永磁同步直线电机,dq轴电感近似相等,即 L_d ≈ L_q = L_s。此时电压方程简化为:

u_d = R_s·i_d + L_d·di_d/dt - ω_e·L_q·i_q
u_q = R_s·i_q + L_q·di_q/dt + ω_e·(L_d·i_d + ψ_f)

磁链方程则变为:

ψ_d = L_d·i_d + ψ_f
ψ_q = L_q·i_q

其中 ω_e 为电角速度,与动子线速度 v 的关系为 ω_e = π·v / τ,τ 为极距。可以看到,dq坐标系下磁链不再显含位置角,电感成为常数,电压方程中仅剩旋转反电势项 ω_e·ψ_f 和交叉耦合项 ω_e·L_q·i_q,这为矢量控制提供了理想基础。

直线电机 abc → dq 坐标系变换与模型结构 abc 坐标系 u_a = R_s·i_a + dψ_a/dt u_b = R_s·i_b + dψ_b/dt u_c = R_s·i_c + dψ_c/dt 电感随位置变化 端部效应显著 Clarke 3/2 变换 αβ 坐标系 u_α = R_s·i_α + dψ_α/dt u_β = R_s·i_β + dψ_β/dt 两相静止 仍含位置信息 Park 旋转变换 dq 坐标系(同步旋转) u_d = R_s·i_d + L_d·di_d/dt − ω_e·L_q·i_q u_q = R_s·i_q + L_q·di_q/dt + ω_e·(L_d·i_d + ψ_f) ψ_d = L_d·i_d + ψ_f ψ_q = L_q·i_q F_e = (3π/2τ)·[ψ_f·i_q + (L_d−L_q)·i_d·i_q] 动态机械模型 F_e − F_L = M·dv/dt + B·v v = dx/dt ω_e = π·v / τ M:动子质量 B:粘滞摩擦系数

三、推力方程:电磁推力与参数关系

在dq坐标系下,直线电机的电磁推力表达式为:

F_e = (3π / 2τ) · [ψ_f·i_q + (L_d − L_q)·i_d·i_q]

对于表贴式永磁同步直线电机,L_d ≈ L_q,推力简化为 F_e = (3π/2τ)·ψ_f·i_q,即推力与q轴电流成正比。这意味着只要控制 i_q,就能线性控制推力;而 i_d 通常被控制为零以实现最大推力电流比。对于内置式或凸极效应明显的直线电机,磁阻推力分量 (L_d − L_q)·i_d·i_q 不可忽略,此时需要采用最大转矩电流比策略分配 i_d 与 i_q。

模型要素 abc坐标系 dq坐标系
电压方程 三相独立,含时变电感 两轴解耦,含旋转反电势
磁链方程 自感+互感+永磁磁链 ψ_d = L_d·i_d + ψ_f,ψ_q = L_q·i_q
推力表达 需先求三相磁链再合成 F_e = (3π/2τ)[ψ_f·i_q + (L_d−L_q)·i_d·i_q]
控制适用性 适合仿真与故障分析 适合矢量控制与电流环设计
参数特性 电感随位置周期变化 电感近似常数,便于线性化

四、动态模型与工程案例

将电磁推力方程与机械运动方程联立,即得到直线电机的完整动态模型:

F_e = (3π/2τ)·[ψ_f·i_q + (L_d − L_q)·i_d·i_q]
F_e − F_L = M·dv/dt + B·v
v = dx/dt
ω_e = π·v / τ

该模型表明:电流环输出 i_q 产生推力,推力克服负载 F_L、惯性力 M·dv/dt 和粘滞摩擦力 B·v 后决定动子速度,速度再通过极距换算为电角速度,反馈回电压方程中的反电势项。这是一个典型的机电耦合双闭环结构。

案例:某平板型永磁同步直线电机,τ = 32 mm,ψ_f = 0.12 Wb,M = 8 kg,B = 15 N·s/m。当 i_q = 5 A、i_d = 0 时,电磁推力 F_e = (3π/2×0.032)×0.12×5 ≈ 88.4 N。若负载为 40 N,则加速度 a = (88.4 − 40 − 15v)/8,初始加速度约 6.05 m/s²,随着速度上升,粘滞摩擦项增大,加速度逐渐减小,最终趋于稳态速度。

五、常见误区与复盘

  • 误区一:直接把旋转电机的dq模型套用到直线电机,忽略极距τ与线速度v的换算关系,导致电角速度计算错误。
  • 误区二:认为所有直线电机都有 L_d = L_q。只有表贴式结构近似成立,内置式或初级开槽较深时凸极效应明显。
  • 误区三:忽略端部效应和法向力。直线电机铁心开断导致三相电感不对称,实际推力存在波动,高精度场合需在模型中补偿。
  • 误区四:把 i_d = 0 控制当作万能策略。当需要弱磁扩速或利用磁阻推力时,必须重新分配 i_d 与 i_q。
工程复盘建议:在控制器参数整定前,先用有限元或实验测取 ψ_f、L_d、L_q、R_s 随位置的变化曲线,再决定采用线性dq模型还是含位置扰动的扩展模型。模型精度直接决定电流环带宽和推力波动抑制效果。

综上,abc坐标系提供了直线电机最原始的物理描述,dq坐标系则通过坐标变换实现了磁链与推力的解耦表达。电压方程、磁链方程、推力方程和机械动态方程共同构成了直线电机伺服控制的理论基石。理解两套坐标系之间的映射关系与各自适用边界,是后续矢量控制、电流环设计和推力波动抑制的前提。