基于新型权重因子校正和延时补偿的永磁同步电动机模型预测转矩控制*

时间：2024-07-28

王海良刘世民陈永清尹玉才

(1.北京科技大学能源与环境工程学院北京 100083；2.新兴际华集团技术中心北京 100070；3.中国五矿集团(唐山曹妃店)矿山控股有限公司唐山 063200；4.泰豪军工集团研究院北京 100166；5.新兴重工集团有限公司北京 100070)

1 引言

永磁同步电动机(Permanent magnet synchronous motor，PMSM)具有功率因数大、效率高和功率密度大等优点，已经在电梯、风电系统和电动汽车等工业领域得到了广泛的应用[1-5]；同时，永磁同步电动机控制技术也在不断改进与完善[6]。目前，广泛采用的高性能永磁同步电动机控制方法有矢量控制和直接转矩控制(Direct torque control，DTC)[7]。与矢量控制相比，直接转矩控制通过选择合适的空间电压矢量，控制定子磁链运动，通过迅速改变负载角，实现对转矩的直接控制。由于该控制方法具有对电动机参数依赖少、鲁棒性强、结构简单和动态响应快的优点，已经成为了学者们的研究热点[8]。然而，传统的DTC 方法采用迟滞比较器和开关表，使得开关频率和转矩脉动很大[9]。为解决以上问题，很多基于最优化控制理论的DTC 方案已经出现，如模型预测直接转矩控制(Model predictive direct torque control，MPC-DTC)和无差拍转矩控制等[10]。

由于MPC 具有动态响应快和易于实现的优点，已经被广泛应用于电动机和功率变换器控制领域中。MPC-DTC 以转矩和磁链为控制目标，通过代价函数选择使转矩和磁链跟踪误差最小的电压矢量，是一种最优化控制方法。目前，研究人员已经在MPC-DTC 领域取得了一些成果，然而在应用过程中仍存在许多问题。首先，MPC-DTC 代价函数中包含转矩和磁链两个物理量，由于它们的量纲和幅值范围等不同，需要设计权重因子来调整两参数在选择最优电压矢量时的重要性，目前权重因子校正主要依赖试凑和试验法，但是其过程比较繁琐。文献[11-12]提出了一种仅以转矩为控制目标的MPC-DTC，该方法将磁链预测过程与转矩计算相结合，省去了代价函数中的磁链项，可以避免权重因子的使用，但失去了跟踪磁链的功能。文献[13]研究了采用仿真和试验法校核权重因子的标准流程，通过仿真或试验参数化分析最佳参数，但由于参数化步长的限制，最佳权重因子依旧无法准确获得。除权重因子设计外，计算延时补偿是永磁同步电动机模型预测控制系统的另一挑战[8]。文献[14]针对三相逆变器模型预测电流控制提出了一种两步预测补偿策略，被广泛应用到电动机控制领域[15-17]，但迄今为止还没有更多的补偿策略出现，方法的单一性限制了模型预测控制理论体系的进一步发展。

为了校正MPC-DTC 代价函数中的权重因子，本文首先建立了永磁同步电动机离散预测模型，通过分析权重因子的本质，利用转矩和磁链对控制电压的响应速度不同，校核MPC 控制器中的权重因子；此外，为降低转矩和电流波动并扩展模型预测控制理论，提出了一种新型的计算延时补偿策略，通过在一个控制周期内进行两次电流采样实现延时时间预估和延时补偿，提高系统的控制性能。

2 PMSM 预测模型

在忽略铁心磁饱和且不计涡流和磁损耗的前提下，表贴式PMSM 在αβ坐标系下的数学模型为

式中，us和is分别为定子电压和电流矢量；Ls和Rs分别为定子绕组电感和电阻；p为电动机极对数；ωm为电动机转速；ψs和ψf分别为气隙磁链和永磁体磁链矢量，其中永磁体磁链满足

式中，ψf为永磁体磁链幅值；θ为转子位置角度。

PMSM 的电磁转矩方程为

式中，Te为输出电磁转矩；p为电动机极对数。

在k时刻，利用一阶欧拉法对式(1)和式(2)进行离散化，得到电动机磁链和电流的未来状态预测表达式

式中，T为开关周期；us(k)为待选电压矢量(控制量)；is(k)和ψs(k)为k时刻的状态量；is(k+1)和ψs(k+1)分别为k+1 时刻的状态量。根据式(5)可知，当预测下一时刻磁链时，需要已知当前时刻磁链值ψs(k)，但电动机中一般不会安装磁链传感器，因此，需要构建观测器实时检测电动机磁链值。将式(1)在k-1 时刻利用欧拉法离散化，得到

与式(6)不同，us(k-1)为k-1 时刻施加的最佳电压矢量，根据式(8)可知，ψs(k-1)又可根据前一时刻状态量计算得到。因此，当磁链初始值ψs(0)给定后，以后每一时刻的磁链都可观测。在t=0 时，假设系统电流为0，气隙磁链只包含永磁体磁链，根据式(3)可知，电动机磁链初始值为

式中，θ0为转子初始位置。

3 PMSM 模型预测转矩控制

3.1 MPC-DTC 实现方法

3.1.1 系统原理

PMSM 模型预测转矩控制原理框图如图1 所示。电动机转速经过转速控制器(Automatic speed regulator，ASR)调节后产生参考转矩Te*，而参考磁链通过最大转矩电流比(Maximum torque per ampere，MTPA)算法计算而得，其中

图1 模型预测转矩控制原理框图

在一个采样周期内，根据当前电动机转速、位置、电流和上一周期状态信息，将7 个待选电压矢量(u1，u2，…，u7)分别代入预测模型得到预测磁链幅值ψs(k+1)和转矩Te(k+1)，并计算出对应的代价函数值，并通过比较选择使代价函数最小的电压矢量及开关状态发送该周期的控制信号，其中，以磁链和转矩为控制目标的代价函数为

式中，λ为权重因子。

3.1.2 问题描述

首先，与传统的控制器(PI)不同，MPC-DTC 控制器是基于最优控制理论的非线性控制器，目前，仍缺乏简单、高效的参数设计方法以保证系统的稳定性，而传统的试凑法具有过程复杂和可靠性低的缺点，因此，研究模型预测转矩控制器权重因子的本质和校核方法具有十分重要的工程价值。

其次，以转矩跟踪为例，模型预测转矩控制过程如图2 所示，Tr为电动机真实输出转矩，To为任一周期内优化选择的最佳电压矢量，Tp表示与7 个电压矢量对应的下一时刻转矩估算值。在tk－1～tk周期内，由图2 可知，在不考虑计算延时情况下，即采样与最佳电压矢量施加都发生在tk－1时刻，u4会使代价函数值最小，且当在该周期内发送对应的开关状态控制信号后，电动机转矩会具有最佳跟踪特性。然而，电压矢量与开关状态选择需要循环执行7 次预测算法，加上排序，庞大的计算会导致信息采样与电压施加之间存在延时td。若在tk－1+td时刻施加u4，一个采样周期后，转矩控制无法达到预期的效果。相比之下，若u2或u3作为该周期内的控制电压，转矩跟踪特性会更加优良。因此，计算延时会影响最佳控制电压选择过程，如果不采取有效的补偿措施，计算延时效应将持续存在于整个控制过程(例如tk和tk+1时刻)，降低系统控制性能。

图2 模型预测转矩控制过程及延时问题

3.2 权重因子设计方法

3.2.1 权重因子本质

代价函数中必须采用权重因子调节磁链和转矩重要性的原因包括：①两变量的量纲和幅值范围不同。首先，磁链的量纲为Wb 而转矩的量纲为N·m，这导致两个参数属性具有本质不同；其次，电动机磁链幅值与其变化范围一般很小，特别是对于中小型电动机[9,18-19]，其磁链范围在0.01～1.00 Wb，而电动机输出转矩的变化范围达到0～20 N·m。因此，在MPC-DTC 算法中，需要首先消除磁链和转矩量纲和幅值范围之间的差异；② 磁链和转矩对控制变量的响应速度不同。该因素还可以表述为当同样的电压施加于电动机后，磁链和转矩在相同时间内的变化量不同，此时，变量的动态特性会影响其在代价函数中的重要性。此时，需要分析磁链和转矩响应速度之间的差异，并进行平衡，便可实现权重因子校正。

3.2.2 标准化

借助标幺化控制思想，为消除变量间的量纲与幅值差异，在把预测的未来磁链ψs(k+1)和转矩Te(k+1)代入代价函数选择最优控制电压之前可通过除以最大值的方式进行标准化，即

式中，ψsmax和Temax分别是电动机最大磁链和最大转矩，其中，Temax可设置为额定值，即Temax=Terated，而ψsmax为电动机输出额定转矩时的磁链值，通过试验测得；ψsn*和ψsn(k+1)为标准化的磁链，Ten*和Ten(k+1)为标准化的转矩，它们的量纲统一为p.u.，且幅值范围为[0，1]。通过标准化操作，代价函数中的两变量具有了相同的量纲和幅值范围，只需要平衡它们对电压变化的响应速度即可完成参数校正。

3.2.3 响应速度平衡

为简化分析，将PMSM 电动机模型转换到旋转轴坐标系下

式中，usr和isr分别为定子电压和电流矢量；usr_q和isr_q分别为q轴电压和电流；ψsr为旋转坐标系下气隙磁链。离散化模型为

在tk～tk+1周期内，假设施加电压发生扰动为Δusr(k)=(Δusr_d(k)，Δusr_q(k))，根据式(17)和式(18)，磁链和电流在相同时间(一个周期)内由该扰动造成的变化量为

式中，Δusr_d(k)和Δusr_q(k)数量级相当，假设两者相等，即

根据式(19)可知，PMSM 的电磁转矩在T内由Δusr(k)引起的变化量为

根据式(22)和式(23)定义转矩磁链响应速度比δ为

为平衡转矩磁链响应速度在代价函数中的影响，权重因子λ设置为

3.3 延时补偿策略

在tk－1～tk周期内(图2)，若在tk－1时刻获得tk－1+td时刻的磁链和转矩值，并利用它们进行未来状态预测，MPC-DTC 执行过程可等效为采样与最佳电压矢量施加都发生在tk－1+td时刻的无延时控制。由式(5)～(7)可知，为得到延时td后的磁链与转矩值，应预先估计该时刻的电流值。在此基础上，本文提出的延时补偿过程应包含三部分：①延时时间计算；② 补偿电流预测；③补偿磁链和转矩预测。以α轴电流iα为例，假设电流在一个采样周期内线性变化，补偿策略的具体实施过程如下。

3.3.1 延时时间计算

传统的MPC-DTC 算法，在每个周期执行的初始阶段，进行一次电动机状态信息采样，且采集的信息直接用于预测算法中，导致工程师无法评估延时时间的长短以及其对电动机控制性能的影响，为此，图3 设计了一种基于“单周期双采样”的延时时间估算方法。在无其他补偿算法执行时，在每一个采样周期内，电流采样操作会被执行两次：第一次仍然在每个周期的初始阶段，第二次则发生在发送开关状态信号之前。其中，在tk－1～tk周期内，第一次采样电流为iar_1(k－1)，第二次采样电流为iar_2(k－1)；在tk～tk+1周期内，第一次和第二次采样电流分别为iar_1(k)和iar_2(k)。根据假设，计算延时时间可由式(26)计算得到

图3 延时时间计算

3.3.2 补偿电流预测

由图4 可知，以tk～tk+1周期为例，获得计算延时时间后，在利用模型预测未来状态之前，可根据式(27)算得补偿电流

则在tk+td的电流为

同理，β轴补偿电流也可得到。

图4 补偿电流预测策略

3.3.3 补偿磁链和转矩预测

考虑权重因子校正和延时补偿的模型预测转矩控制原理框图如图5 所示。

图5 考虑权重因子校正和延时补偿的MPC-DTC 原理框图

4 试验验证

首先，试验是为了验证所设计的MPC-DTC 控制器权重因子能够保证系统的动态以及稳态性能，其次，通过比较延时补偿前后的系统性能，验证本文提出的权重因子和延时补偿设计与分析方法是有效的。本文针对一台表贴式永磁同步电动机进行算法验证，其系统参数如表1 所示。试验控制平台以DSP TMS320F2812 为主控芯片，逆变器采用智能功率模块(IPM)PM75RLA12，开关电动机位置检测采用光电编码器，加载和记录转速由Magtrol 公司的10 kW 测功台完成，磁链和电磁转矩数据由数字处理器DSP 记录，并通过RS-485发送至上位机。

根据式(24)和式(25)可得试验中磁链权重因子为λ=171.58。

表1 永磁同步电动机控制系统参数

4.1 稳态性能

利用本文提出的MPC-DTC 算法控制电动机旋转于3 000 r/min，图6 和图7 分别为电动机空载和负载(4.5 N·m)时的稳态特性。首先，电动机能按照给定参考转速旋转，稳态误差为0，说明所设计的磁链权重因子满足稳态性能要求；其次，空载状态下，平均气隙磁链约为0.057 6 Wb，负载时，磁链值略大于空载值(0.057 9 Wb)；最后，虽然电动机转矩存在波动，但其平均输出值展现出了良好的跟踪性。

图6 带延时补偿空载稳态特性

图7 带延时补偿负载稳态特性

4.2 动态性能

利用本文提出的方法控制电动机从0 加速至500 r/min，在 0.3 s 时刻，设置参考转速为3 000 r/min，并保持至0.6 s，然后，控制电动机降速至1 000 r/min。图8 为电动机的加减速动态特性图。首先，在电动机加减速过程中，电动机磁链都会变大，当电动机转速再次达到稳定后，磁链恢复至0.056 Wb；其次，加减速过程中输出最大电磁转矩分别约为±20 N·m，这是由转速控制器限幅输出导致的，电动机稳定运行时，输出电磁转矩略大于0；当电动机从500 r/min 加速至3 000 r/min，上升时间约为80 ms，超调非常小(小于10 r/min)，且降速特性与加速特性非常相似，这表明所设计的权重因子及MPC-DTC 算法具有良好的动态性能。

图8 带延时补偿加减速动态特性

为验证电动机的抗负载干扰特性，控制电动机稳定运行至2 000 r/min，在0.25 s 时刻，突加4.5 N·m 负载转矩，图9 为突加负载时的电动机动态特性。突加负载时，电动机转速大范围下降(94 r/min)，但电动机转矩很快恢复至2 000 r/min，调节时间约为 30 ms。试验表明，所设计的MPC-DTC 参数具有良好的抗负载干扰特性。

图9 带延时补偿突加负载动态特性

4.3 延时补偿前后性能对比

当采用无延时补偿的MPC-DTC 算法控制电动机工作于3 000 r/min 时，电动机的空载稳态特性如图10 所示，可以看出，磁链和转矩脉动大小分别为0.021 Wb 和4.2 N·m，与图6 相比(磁链和转矩脉动分别为0.0178 Wb 和3.6 N·m)，脉动增长约为17%，所以，所提出的延时补偿策略能够降低控制磁链和转矩脉动，提高系统的性能。

图10 无延时补偿空载稳态特性