从电枢电压到转子转角：手把手拆解直流电机数学模型，附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 基础模型搭建步骤

1. 创建新模型：启动Simulink，选择"Blank Model"模板

2. 添加关键模块：

   • 阶跃输入（Step）
   • 传递函数（Transfer Fcn）
   • 示波器（Scope）
   • 工作区输出（To Workspace）

3. 参数设置示例：

% 电机参数（以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 完整仿真模型结构


关键模块连接关系：

1. Step → Transfer Fcn → Scope
2. Transfer Fcn参数设置为[K],[tau 1]
3. 仿真时间设为2秒，步长0.001秒

4. 动态响应分析与工程启示

4.1 典型参数影响对比

观察结论：

• 转动惯量$J_m$增大 → 系统响应变慢，振荡加剧
• 摩擦系数$B_m$增大 → 阻尼效果增强，抑制振荡

4.2 实际调试技巧

1. 参数辨识方法：

   • 通过空载阶跃响应曲线拟合$\tau$和$K$
   • 测量稳态转速与电压比确定$K$
   • 测量转速下降曲线确定$B_m$

2. 控制器设计基础：

% PID控制器初步设计 Kp = 0.8*(2*Jm*wn)/K; Ki = 0.5*(wn^2)*Jm/K; Kd = 0.1*(2*Jm*wn)/K;

3. 常见问题排查：
   • 响应异常快：检查$J_m$是否设置过小
   • 持续振荡：验证$B_m$值是否合理
   • 稳态误差：确认是否忽略$L_a$导致模型失真

在完成基础仿真后，尝试修改Step信号的幅值从1V逐步增加到5V，观察系统在不同工作点的线性度表现。实际项目中，我们曾遇到某医疗设备驱动电机在低速段出现明显非线性，最终发现是电刷接触电阻变化导致$R_a$参数不恒定，这提醒我们简化模型有其适用边界。