直流电机数学模型与状态空间分析——精炼版

直流电机数学模型与状态空间分析

? 核心提要：核心提要
直流电机模型是所有电机控制算法（包含 PMSM FOC）的物理基石。本文从基础的机电耦合方程出发，推导其连续状态空间表达式，对比不同的离散化方案，并落脚于工程实际中的参数辨识与控制策略设计。
前置/关联节点：[[现代控制理论基础]]、[[坐标变换与FOC原理]]、[[滑模观测器 SMO]]

1. 核心物理模型与基础方程

直流电机的动态特性由机械域、电气域以及连接两者的机电耦合关系共同决定。

1.1 机械运动方程

描述转子在电磁转矩与外部阻力作用下的动力学行为：

\[T_e - T_L - B \omega_m = J \frac{d\omega_m}{dt} \]

• \(T_e\): 电磁转矩 (N·m)
• \(T_L\): 外部负载转矩 (N·m)
• \(B\): 黏性摩擦系数 (N·m·s)
• \(J\): 系统总转动惯量 (kg·m²)，包含电机本体 \(J_{motor}\) 与折算负载 \(J_{load}\)
• \(\omega_m\): 机械角速度 (rad/s)

1.2 电气电压方程 (电枢平衡方程)

将电枢绕组等效为 RL 串联与反电动势的回路：

\[u_a = R_a i_a + L_a \frac{di_a}{dt} + e_b \]

• \(u_a\): 端电压 (V)
• \(i_a\): 电枢电流 (A)
• \(R_a\): 电枢电阻 (\(\Omega\))
• \(L_a\): 电枢电感 (H)
• \(e_b\): 反电动势 (V)

1.3 机电耦合方程与能量守恒

电机本质是能量转换装置，其核心桥梁由以下两式建立：

1. 转矩生成：$$T_e = K_t i_a$$
2. 反电势生成：$$e_b = K_e \omega_m$$

? 工程洞察：工程洞察：\(K_t\) 与 \(K_e\) 的等价性
在 SI 国际单位制下，基于理想状态下的能量守恒原则（电磁功率 = 机械功率，即 \(e_b i_a = T_e \omega_m\)），转矩常数与反电动势常数在数值上必然相等：\(K_t = K_e\)。
微观本质上，\(K_e = \frac{p N}{2\pi a} \Phi\)，由极对数、绕组结构及单极磁通决定。

---

2. 系统动态响应与稳态特性

2.1 阶跃负载扰动响应 (开环特性)

当系统处于稳态 (\(u_a\) 恒定) 且突加负载 \(T_L\) 时，系统的演变逻辑如下：

1. 合力矩突变：\(T_e < T_L\)，加速度 \(\frac{d\omega_m}{dt} < 0\)，转速 \(\omega_m\) 跌落。
2. 反电势下降：由于 \(e_b = K_e \omega_m\)，\(e_b\) 随之降低。
3. 电流被动攀升：由 \(\frac{di_a}{dt} = \frac{u_a - e_b - R_a i_a}{L_a}\) 可知，压差增大导致电流 \(i_a\) 上升。
4. 新稳态平衡：\(T_e\) 随电流增加直至重新等于 \(T_L\)。最终表现为：转速下降，电流增大（直流电机的软机械特性）。

? 笔记：闭环控制的必要性
纯物理反馈无法在维持恒压 \(u_a\) 的同时抵抗负载扰动并保持转速。必须引入 PI 速度环，通过动态调节 \(u_a\) (PWM 占空比) 主动注入电流以抵抗扰动，打破物理死锁。

2.2 双时间常数解耦

• 电气时间常数：\(\tau_e = \frac{L_a}{R_a}\) (通常为 ms 级别)
• 机械时间常数：\(\tau_m = \frac{J R_a}{K_t K_e}\) (通常为 10~100 ms 级别)
  工程结论：电气动态响应速度远高于机械响应速度。这构成了串级控制（外环速度环，内环电流环）的理论基础。

---

3. 状态空间表达 (面向算法实现)

为了进行现代控制策略（如 LQR、状态观测器设计）分析，需将连续微分方程转化为标准状态空间形式：\(\dot{x} = Ax + Bu + Ed\)。

选取电流与转速（储能元件对应的物理量）作为状态变量：\(x = [i_a, \omega_m]^T\)。

\[\begin{bmatrix} \dot{i}_a \\ \dot{\omega}_m \end{bmatrix} = \begin{bmatrix} -\frac{R_a}{L_a} & -\frac{K_e}{L_a} \\ \frac{K_t}{J} & -\frac{B}{J} \end{bmatrix} \begin{bmatrix} i_a \\ \omega_m \end{bmatrix} + \begin{bmatrix} \frac{1}{L_a} \\ 0 \end{bmatrix} u_a + \begin{bmatrix} 0 \\ -\frac{1}{J} \end{bmatrix} T_L \]

• 矩阵 A (系统矩阵)：反映电机固有的衰减与振荡特性（包含全部出厂物理参数）。
• 矩阵 B (输入矩阵)：表明控制量 \(u_a\) 只能直接干预电流，必须通过矩阵 A 的耦合作用才能间接控制转速。

---

4. 离散化方案与 MCU 部署

在 MCU (如 DSP/ARM) 中实现数字控制或状态观测器时，必须对连续模型 \(A, B\) 进行离散化 (\(A_d, B_d\))，采样周期为 \(T_s\)。

4.1 前向欧拉法 (Forward Euler)

直接使用差分近似导数 \(\dot{x} \approx \frac{x_{k+1} - x_k}{T_s}\)。

• \(A_d = I + A T_s\)
• \(B_d = B T_s\)

?? 警告：适用场景与局限
仅适用于高开关频率 (\(T_s\) 极小，如 FOC 的 10kHz~20kHz 中断)。本质上是矩阵指数的一阶泰勒截断，若采样率不足，离散极点将跑出单位圆导致观测器发散爆炸。

4.2 零阶保持器法 (Zero-Order Hold, ZOH)

基于 MCU 的 PWM 保持特性，通过对连续解析解进行积分，保留所有泰勒高阶项：

• \(A_d = e^{A T_s}\)
• \(B_d = (\int_{0}^{T_s} e^{A \eta} d\eta) B\)

? 建议：算法落地建议
在实际编写无感 FOC (Sensorless FOC) 或卡尔曼滤波代码时，通常在 MATLAB 中通过 c2d 函数离散化算好 \(A_d, B_d\) 的数值，以宏定义常数的形式固化在 C 代码中，将复杂的微积分运算降维为 MCU 极易处理的线性乘加运算。

---

5. 总结：通往 FOC 的终极映射

深刻理解直流电机模型，是掌握交流伺服核心算法的必经之路。在永磁同步电机 (PMSM) 的磁场定向控制 (FOC) 中，利用 Park 变换将三相静止坐标系投影到随转子同步旋转的 \(d-q\) 坐标系后，其 \(q\) 轴电压方程为：

\[u_q = R i_q + L_q \frac{di_q}{dt} + \omega_e \psi_f \]

对比直流电机的 \(u_a = R_a i_a + L_a \frac{di_a}{dt} + K_e \omega_m\)：

• \(i_q\) 等效于电枢电流 \(i_a\)
• \(\omega_e \psi_f\) 等效于反电动势 \(e_b\)

FOC 的本质，正是利用数学上的坐标变换，在算法层面上将一台复杂的非线性交流电机，完美解耦重塑成一台纯粹、线性、极易控制的直流电机。