BLDC电机控制：PI与ADRC方案对比与实现

1. 直流无刷电机控制概述

直流无刷电机（BLDC）凭借其高效率、长寿命和低维护成本等优势，已成为现代工业运动控制领域的核心执行器件。与传统有刷直流电机相比，其通过电子换相取代机械换向，消除了电刷火花和磨损问题。但在控制层面，这种结构变化也带来了新的挑战——需要精确控制三相逆变器的开关时序来实现转子位置跟踪。

在实际工程应用中，我们通常采用转速-电流双闭环控制架构。这种分层设计源于一个基本控制思想：电流环作为内环负责电磁转矩的快速调节，转速环作为外环处理机械动态。这种结构能有效解耦电机的电磁和机械时间常数，电流环的带宽通常设计为转速环的5-10倍。

关键设计原则：内环（电流）的响应速度必须显著快于外环（转速），否则会出现控制耦合导致系统振荡。这个原则同样适用于其他多环控制系统。

2. 传统PI控制方案实现

2.1 电流环PI设计要点

电流环的核心任务是实现q轴电流（转矩分量）的精确跟踪。在Simulink中，我们首先需要建立电机的电压方程模型：

matlab复制function dI = motor_current(V, I, L, R, Ke, omega)
    % V: 输入电压
    % I: 当前电流
    % L: 绕组电感
    % R: 绕组电阻
    % Ke: 反电动势常数
    % omega: 转子角速度
    dI = (V - R*I - Ke*omega)/L;
end

PI参数整定采用工程常用的"带宽法"：

1. 确定期望闭环带宽ω_bw（通常取开关频率的1/10）
2. 比例系数 Kp = L·ω_bw
3. 积分时间常数 Ti = L/R

实际调试时需注意：

• 电感L会随电流饱和而变化，需留有余量
• 采样延迟会限制实际可达带宽
• 过高的Kp会导致PWM开关噪声被放大

2.2 转速环PI的特殊考量

转速环处理的是机械动态，其PI参数整定更为复杂。推荐采用"模最优"整定方法：

1. 将电流环简化为惯性环节：1/(1 + 2Ts)
2. 电机机械方程：J·dω/dt = Kt·Iq - B·ω - Tl
3. 典型参数选择：
   • Kp = (2·J)/(3·Kt·TΣ)
   • Ti = 4·TΣ

其中TΣ为系统总滞后时间（包含电流环等效滞后）。

实测经验：转速PI的输出限幅应设置为电机允许的最大电流值，同时需要加入抗积分饱和逻辑。我曾在一个伺服项目中，因忽略这点导致电机过流烧毁MOSFET。

3. ADRC控制方案实现

3.1 ADRC核心思想解析

自抗扰控制器(ADRC)的精妙之处在于其将系统内部动态和外部扰动统一视为"总扰动"，通过扩张状态观测器(ESO)实时估计并补偿。对于我们的BLDC系统：

1. 二阶ESO设计：

   matlab复制function [z1, z2, z3] = eso(y, u, h, beta1, beta2, beta3)
       persistent z1_prev z2_prev z3_prev
       e = y - z1_prev;
       z1 = z1_prev + h*(z2_prev + beta1*e);
       z2 = z2_prev + h*(z3_prev + beta2*e + b0*u);
       z3 = z3_prev + h*beta3*e;
       % 更新状态
       z1_prev = z1; z2_prev = z2; z3_prev = z3;
   end
   

   其中z3就是对总扰动的估计。

2. 非线性反馈组合：
   采用最速控制综合函数fhan()代替线性PID，其数学表达为：

   matlab复制function u = fhan(x1, x2, r, h)
       d = r*h^2;
       a0 = h*x2;
       y = x1 + a0;
       a1 = sqrt(d*(d+8*abs(y)));
       a2 = a0 + sign(y)*(a1-d)/2;
       sy = (sign(y+d) - sign(y-d))/2;
       a = (a0+y-a2)*sy + a2;
       sa = (sign(a+d) - sign(a-d))/2;
       u = -r*(a/d - sign(a))*sa - r*sign(a);
   end
   

3.2 Simulink实现技巧

在搭建ADRC的Simulink模型时，有几个关键细节需要注意：

1. ESO的离散化实现：

   • 采用Tustin双线性变换而非欧拉法，可提高数值稳定性
   • 添加输出限幅防止积分饱和
   • 初始状态设置为电机启动时的典型工况

2. 参数整定经验公式：

   • 观测器带宽ωo ≈ (3~5)ωc（ωc为控制器带宽）
   • b0取系统增益的标称值（如Kt/J）
   • 非线性因子r的取值与转速量程相关

3. 工程实用技巧：

   • 添加扰动观测输出滤波（一阶低通）
   • 对z3（扰动估计）进行死区处理
   • 实现平滑的ADRC-PI切换逻辑

4. 对比仿真与结果分析

4.1 测试工况设计

为全面评估两种控制策略，我们设计了三类测试场景：

1. 阶跃响应测试：

   • 空载启动至额定转速
   • 50%负载突加
   • 转速指令突变

2. 抗扰测试：

   • 模拟负载转矩波动（±20%阶跃）
   • 电源电压跌落（80%→60%）
   • 参数失配（转动惯量增加50%）

3. 鲁棒性测试：

   • 电机参数漂移（R+30%，Kt-15%）
   • 测量噪声注入（10dB SNR）
   • 采样延时增加（从100μs到1ms）

4.2 性能指标对比

我们定义以下量化评估指标：

4.3 典型波形对比

1. 突加负载工况：

   • PI控制：转速跌落85rpm，恢复过程出现2次振荡
   • ADRC：转速跌落仅35rpm，无振荡平滑恢复

   关键差异源于ADRC的扰动前馈补偿能力，其ESO在负载变化的5ms内就准确估计出了扰动转矩。

2. 参数失配测试：
   当故意将模型中的转动惯量设置偏差+50%时：

   • PI控制：超调增至25%，出现持续小幅振荡
   • ADRC：性能指标几乎不变，展现出卓越的鲁棒性

5. 工程实施经验分享

5.1 PI控制调试心得

1. 参数整定步骤：

   • 先调电流环：将Kp从零逐渐增大至出现轻微振荡，然后回退30%
   • 再调转速环：固定Ti=4TΣ，调整Kp使阶跃响应临界阻尼
   • 最后微调Ti：在±20%范围内优化调节时间

2. 常见问题处理：

   • 高频振荡：检查PWM死区时间与采样同步性
   • 稳态误差：确认积分器未饱和，适当增大Ti
   • 响应迟缓：检查速度反馈滤波是否过度

5.2 ADRC实施要点

1. 降低实现难度的技巧：

   • 先用线性ADRC（LADRC）验证基本功能
   • ESO初始阶段采用较低带宽，稳定后逐步提高
   • 添加手动扰动注入测试观测器性能

2. 参数整定流程：

   matlab复制% 简易整定脚本示例
   wc = 2*pi*50;   % 期望控制带宽(Hz)
   wo = 3*wc;      % 观测器带宽
   b0 = Kt_nom/J_nom;
   beta1 = 3*wo;
   beta2 = 3*wo^2;
   beta3 = wo^3;
   r = 100*wc;     % 非线性因子
   

3. 数字实现注意事项：

   • 采用Q15格式定点数提高DSP运算效率
   • ESO状态变量需要定期抗饱和复位
   • 添加安全逻辑防止异常输出

在实际的电动舵机项目中，我们将ADRC与传统PI进行对比测试：在相同规格电机下，ADRC使位置跟踪误差减小了62%，同时在-20℃低温启动时成功率达100%（PI方案有15%的启动失败率）。这充分证明了ADRC在复杂工况下的优势。