4. 直线电机数学模型
直线电机的推力控制与补偿，必须建立在一个能够同时描述电磁瞬态与机械运动的数学模型之上。与旋转电机不同，直线电机没有齿轮、丝杠等中间传动环节，动子直接驱动负载，因此其模型中的电磁推力、法向力、端部效应和参数非线性会直接反映到位置与速度输出上。本节围绕 dq 轴建模、电磁推力方程、机械运动方程、参数辨识与仿真验证五个环节，构建可用于控制器设计的完整模型链路。
一、dq轴建模：从三相到两轴
直线电机的初级绕组通常为三相分布绕组，若直接在三相静止坐标系下建立电压方程，电感矩阵随动子位置周期变化，耦合项复杂。通过 Clarke 变换与 Park 变换，可将三相变量投影到与动子磁场同步的 dq 坐标系中。对于永磁同步直线电机（PMSLM），忽略铁耗与磁路饱和时，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)其中 R_s 为初级电阻，L_d、L_q 为 dq 轴电感，ψ_f 为永磁体磁链，ω_e = π·v/τ 为电角速度，v 为动子线速度，τ 为极距。由于直线电机铁芯开断，d 轴与 q 轴磁路不对称程度通常高于旋转电机，L_d ≠ L_q 的凸极效应在推力波动分析中不可忽略。dq 轴建模的核心价值在于：将随位置变化的电感矩阵转化为直流量，使电流环可以按直流电机方式独立设计。
二、电磁推力方程与推力分量
在 dq 坐标系下，PMSLM 的电磁推力可写为：
F_e = (3π / 2τ) · [ψ_f · i_q + (L_d − L_q) · i_d · i_q]该式包含两个分量：永磁推力分量 ψ_f·i_q 与磁阻推力分量 (L_d−L_q)·i_d·i_q。对于表贴式永磁直线电机，L_d ≈ L_q，磁阻分量近似为零，采用 i_d = 0 控制即可使推力与 i_q 成正比，控制简单。对于内嵌式或凸极结构，适当施加负 i_d 可产生磁阻推力，提高推力密度，但会引入额外的铜耗与参数敏感性。工程中还需注意：直线电机的端部效应会使 ψ_f 沿行程轻微波动，齿槽力则表现为与位置相关的周期性扰动，二者在推力方程中通常作为附加扰动项处理。
三、机械运动方程与扰动项
动子沿导轨的直线运动由牛顿第二定律描述：
M · dv/dt = F_e − F_L − B_v · v − F_fric − F_detent(x) − F_end(x)其中 M 为动子与负载总质量，B_v 为粘滞摩擦系数，F_fric 为库仑摩擦，F_detent(x) 为齿槽力，F_end(x) 为端部效应力。位置与速度关系为 dx/dt = v。该方程揭示了推力控制的核心难点：控制器输出的 F_e 必须同时克服负载、惯性、摩擦和周期性扰动，才能获得平稳的速度与位置响应。若忽略 F_detent 与 F_end，在低速段会出现明显的推力波动与速度纹波。
四、参数辨识方法与仿真验证
模型精度取决于参数准确性。常用辨识方法如下表所示：
| 参数 | 辨识方法 | 激励条件 | 注意事项 |
|---|---|---|---|
| R_s | 直流注入法 | d 轴通直流，动子堵转 | 需等待温度稳定，避免温升导致阻值漂移 |
| L_d / L_q | 高频注入或电压阶跃响应 | d/q 轴分别注入小信号 | 注意磁饱和与交叉耦合影响 |
| ψ_f | 反电动势法 | 外力拖动动子匀速运动 | 测量空载线电压，取基波幅值 |
| M | 加减速实验 | 施加已知推力，记录加速度 | 需扣除摩擦与齿槽力影响 |
| B_v / F_fric | 多速度点稳态实验 | 不同速度下匀速运行 | 用最小二乘拟合推力-速度曲线 |
仿真验证建议：在 MATLAB/Simulink 或 PLECS 中搭建 dq 轴模型，将辨识参数回填，先验证空载反电动势波形，再验证阶跃推力下的速度响应，最后加入齿槽力与端部效应查表模型，对比实测电流与速度波形。若仿真与实测在低速段偏差较大，优先检查 F_detent 的相位与幅值是否准确。
常见误区：一是直接用旋转电机的参数代替直线电机，忽略端部效应导致的 ψ_f 波动；二是辨识 L_d、L_q 时未考虑动子位置对电感的影响，导致电流环增益失配；三是机械方程中遗漏 F_detent，使仿真速度曲线过于理想，无法复现实测纹波。建议在模型中加入位置相关的扰动查表，并在控制器中预留前馈补偿接口。
五、从模型到控制的衔接
完整的直线电机模型应形成“电压方程 → 电流环 → 推力方程 → 机械方程 → 位置/速度反馈”的闭环链路。dq 轴建模解决了电磁瞬态的直流量化问题，推力方程给出了 i_q 到 F_e 的映射，机械方程则决定了系统对外部扰动与负载变化的响应。参数辨识与仿真验证是保证模型可信度的必要步骤。只有模型中的 R_s、L_d、L_q、ψ_f、M、B_v 以及周期性扰动项均经过实验校准，后续的推力补偿策略才有可靠的被控对象基础。