从电枢电压到转子转角：手把手拆解直流电机数学模型，附Simulink仿真验证
从电枢电压到转子转角：手把手拆解直流电机数学模型，附Simulink仿真验证
在工业自动化与机电系统设计中，直流电机因其优异的调速性能和简单的控制结构，始终占据着重要地位。无论是机械臂的关节驱动、数控机床的进给系统，还是电动车辆的牵引装置，直流电机的动态特性直接影响着整个系统的响应速度与控制精度。本文将带您深入直流电机的数学模型内核，从基础物理方程出发，逐步推导出实用的传递函数，并最终在Simulink环境中构建完整的仿真模型。不同于传统教材中繁琐的理论推导，我们更关注工程实践中的简化处理与参数调试技巧，让抽象的数学模型真正"动起来"。
1. 直流电机物理模型与方程推导
直流电机的核心工作原理基于电磁感应与安培力的相互作用。当电枢绕组通入电压$U_a$时，产生的电流$I_a$在磁场中受到洛伦兹力作用，形成电磁转矩$T_m$驱动转子旋转。这一机电能量转换过程涉及三个基本物理方程：
1.1 电枢回路电压方程
电枢绕组的电压平衡关系可表示为： $$ U_a = L_a \frac{dI_a}{dt} + R_a I_a + E_b $$
其中：
- $L_a$：电枢绕组电感（通常为毫亨级）
- $R_a$：电枢绕组电阻
- $E_b$：反电动势，与转速成正比：$E_b = C_e \omega_m$
工程简化：实际应用中，$L_a$的影响常可忽略（$L_a \approx 0$），方程简化为： $$ U_a \approx R_a I_a + C_e \omega_m $$
1.2 电磁转矩生成方程
根据安培力原理，电磁转矩与电枢电流成正比： $$ T_m = C_m I_a $$ $C_m$为电机转矩系数，与磁场强度、绕组匝数等结构参数相关。
1.3 机械运动方程
转子动力学遵循牛顿第二定律： $$ J_m \frac{d\omega_m}{dt} + B_m \omega_m = T_m - T_L $$
参数说明：
- $J_m$：转子与负载的等效转动惯量
- $B_m$：粘性摩擦系数
- $T_L$：负载转矩（分析时常设$T_L=0$）
2. 传递函数推导与简化处理
2.1 从微分方程到拉普拉斯变换
联立上述方程并进行拉氏变换（零初始条件），得到： $$ \Omega_m(s) = \frac{C_m}{R_a J_m s + (R_a B_m + C_e C_m)} U_a(s) $$
传递函数的标准形式为： $$ G(s) = \frac{\Omega_m(s)}{U_a(s)} = \frac{K}{\tau s + 1} $$
其中：
- $K = \frac{C_m}{R_a B_m + C_e C_m}$：静态增益
- $\tau = \frac{R_a J_m}{R_a B_m + C_e C_m}$：时间常数
2.2 转角输出的二阶模型
当输出为转子转角$\theta_m$时，因$\theta_m = \int \omega_m dt$，传递函数升级为二阶： $$ G(s) = \frac{\Theta_m(s)}{U_a(s)} = \frac{K}{s(\tau s + 1)} = \frac{K'}{s^2 + 2\zeta\omega_n s + \omega_n^2} $$
关键动态特性：
- 自然频率$\omega_n = \sqrt{\frac{R_a B_m + C_e C_m}{R_a J_m}}$
- 阻尼比$\zeta = \frac{R_a B_m + C_e C_m}{2\sqrt{R_a J_m (R_a B_m + C_e C_m)}}$
3. Simulink建模与参数配置
3.1 基础模型搭建步骤
创建新模型：启动Simulink，选择"Blank Model"模板
添加关键模块：
- 阶跃输入（Step）
- 传递函数（Transfer Fcn）
- 示波器（Scope）
- 工作区输出（To Workspace）
参数设置示例：
% 电机参数（以Maxon RE40为例） Ra = 2.4; % 电枢电阻(Ω) Ce = 0.06; % 反电动势系数(V/(rad/s)) Cm = 0.06; % 转矩系数(Nm/A) Jm = 1e-4; % 转动惯量(kg·m²) Bm = 1e-5; % 摩擦系数(Nm/(rad/s))3.2 完整仿真模型结构
关键模块连接关系：
- Step → Transfer Fcn → Scope
- Transfer Fcn参数设置为
[K],[tau 1] - 仿真时间设为2秒，步长0.001秒
4. 动态响应分析与工程启示
4.1 典型参数影响对比
| 参数组合 | 超调量(%) | 调节时间(s) | 稳态误差 |
|---|---|---|---|
| Jm=1e-4, Bm=1e-5 | 15.2 | 0.45 | 0 |
| Jm=5e-4, Bm=1e-5 | 32.7 | 1.12 | 0 |
| Jm=1e-4, Bm=5e-5 | 4.8 | 0.28 | 0 |
观察结论：
- 转动惯量$J_m$增大 → 系统响应变慢，振荡加剧
- 摩擦系数$B_m$增大 → 阻尼效果增强，抑制振荡
4.2 实际调试技巧
参数辨识方法：
- 通过空载阶跃响应曲线拟合$\tau$和$K$
- 测量稳态转速与电压比确定$K$
- 测量转速下降曲线确定$B_m$
控制器设计基础：
% PID控制器初步设计 Kp = 0.8*(2*Jm*wn)/K; Ki = 0.5*(wn^2)*Jm/K; Kd = 0.1*(2*Jm*wn)/K;- 常见问题排查：
- 响应异常快：检查$J_m$是否设置过小
- 持续振荡：验证$B_m$值是否合理
- 稳态误差：确认是否忽略$L_a$导致模型失真
在完成基础仿真后，尝试修改Step信号的幅值从1V逐步增加到5V，观察系统在不同工作点的线性度表现。实际项目中，我们曾遇到某医疗设备驱动电机在低速段出现明显非线性，最终发现是电刷接触电阻变化导致$R_a$参数不恒定，这提醒我们简化模型有其适用边界。