3. 直流电机数学模型：电压平衡方程、转矩方程、机械特性方程推导
做电机控制这些年，我有个很深的体会：数学模型就是电机的“身份证”。你连它的脾气秉性都不了解，怎么去控制它？今天咱们就来聊聊直流电机最核心的三个方程——电压平衡、转矩方程、机械特性。搞懂了这些，后面的仿真和控制才有根。
核心观点：直流电机的数学模型是后续所有控制算法（PID、电流环、速度环）的理论基石。三个方程分别描述了电机的电气特性、力学特性和稳态性能。
3.1 电压平衡方程——电机的“电路灵魂”
先看最简单的等效电路。直流电机电枢回路，说白了就是一个电阻 R 和一个电感 L 串联，再串一个反电动势 E。反电动势是啥？就是电机转起来之后，电枢线圈切割磁感线产生的电压，方向跟外加电压相反。
根据基尔霍夫电压定律，我们写出：
U = E + I·R + L·(dI/dt)其中：
- U — 电枢两端外加电压 (V)
- E — 电枢反电动势 (V)
- I — 电枢电流 (A)
- R — 电枢回路总电阻 (Ω)
- L — 电枢回路总电感 (H)
反电动势 E 跟转速和磁通成正比：
E = Ce·Φ·n这里 Ce 是电动势常数，Φ 是每极磁通，n 是转速。
我的经验：在实际项目中，电阻 R 不是一成不变的。电机发热后，铜线电阻会上升，我见过有人仿真时忽略这个，结果电流环参数整定出来实际跑偏了20%。建议做高精度仿真时，把 R 写成温度的函数。
稳态时，电流不变，dI/dt = 0，方程简化为：
U = E + I·R这个简化形式在工程估算中非常常用。你想想看，堵转时 n=0，E=0，电流就是 U/R，这就是启动电流的来源。我遇到过有人选驱动器时没算这个，结果一启动就过流保护。
3.2 转矩方程——力从哪来？
电磁转矩是电机能拖动负载的根本原因。它由电枢电流和磁场相互作用产生：
Tem = CT·Φ·I其中：
- Tem — 电磁转矩 (N·m)
- CT — 转矩常数
- Φ — 每极磁通 (Wb)
- I — 电枢电流 (A)
这里有个重要的关系：CT 和 Ce 不是独立的。在 SI 单位制下，它们满足：
CT = 9.55 × Ce为什么会这样？因为电磁功率 Pem 既等于 E·I，又等于 Tem·ω。把 ω = 2πn/60 代进去一算，就出来了。
| 参数 | 符号 | 单位 | 典型值范围 |
|---|---|---|---|
| 电动势常数 | Ce | V/(rpm·Wb) | 0.01 ~ 0.5 |
| 转矩常数 | CT | N·m/(A·Wb) | 0.1 ~ 5.0 |
| 电枢电阻 | R | Ω | 0.1 ~ 10 |
| 电枢电感 | L | mH | 0.5 ~ 50 |
注意：转矩方程中的 I 是电枢电流，不是励磁电流。对于永磁直流电机，Φ 是常数，转矩就只跟电流成正比——这就是为什么电流环那么重要。我曾经在调试一个伺服系统时，发现转矩波动，最后查出来是电流采样有毛刺，滤波没做好。
3.3 机械特性方程——转速与转矩的“博弈”
机械特性描述的是稳态下转速 n 与电磁转矩 Tem 的关系。把前面两个方程联立起来：
从电压平衡方程：U = E + I·R，代入 E = Ce·Φ·n
从转矩方程：Tem = CT·Φ·I → I = Tem / (CT·Φ)
消去 I，得到：
n = U / (Ce·Φ) - [R / (Ce·CT·Φ²)] · Tem这就是直流电机的机械特性方程。它是一条直线：
- 斜率：-R / (Ce·CT·Φ²)，决定了特性曲线的“硬度”
- 空载转速：n0 = U / (Ce·Φ)，即 Tem=0 时的转速
- 堵转转矩：Tst = (CT·Φ·U) / R，即 n=0 时的转矩
说白了，这条线越平（斜率绝对值越小），电机的“硬特性”越好，负载变化时转速波动小。我做过一个张力控制项目，要求转速波动小于0.5%，最后选了低电阻、大磁通的电机，机械特性非常硬。
工程启示：机械特性方程告诉我们三个调节电机性能的手段：
- 调压 U — 改变空载转速，曲线平移
- 调磁 Φ — 改变斜率，弱磁升速
- 调电阻 R — 改变斜率，但效率会下降
3.4 知识体系总览
下面这张图把三个方程的关系串起来了。你看，电压方程和转矩方程通过电流 I 和反电动势 E 耦合，机械特性是它们的稳态结果。
3.5 三个方程的内在联系
我习惯把这三个方程看作一个整体：
- 电压方程管的是“电”的平衡——电压怎么分配
- 转矩方程管的是“力”的产生——电流怎么出力
- 机械特性管的是“运动”的结果——转速和转矩怎么配合
在动态过程中，这三个方程通过电流和转速互相耦合。你想想看，负载一变，转速就变，反电动势跟着变，电流跟着变，转矩又跟着变——这就是一个闭环。做仿真时，这三个方程就是核心的微分方程组。
避坑指南：我曾经在搭建仿真模型时，直接把机械特性方程当动态方程用，结果动态响应完全不对。记住：机械特性只适用于稳态，动态过程必须用电压方程和转矩方程的微分形式。