CN106130429A - 无轴承永磁同步电机预测控制器及构造方法 - Google Patents

Abstract

本发明公开一种电力传动控制设备领域中的无轴承永磁同步电机预测控制器及构造方法，Clark变换器的输出端分别连接转矩绕组磁链与转矩初始观测模块以及悬浮力绕组磁链初始观测模块的输入端，悬浮力绕组磁链初始观测模块的输出端连接悬浮力绕组磁链预测模块的输入端，悬浮力绕组磁链预测模块的输出端连接悬浮力观测模块的输入端；转矩绕组磁链与转矩初始观测模块的输出端分别连接转矩绕组磁链与转矩预测模块和转矩与电压变换模块的输入端，转矩绕组磁链与转矩预测模块的输出端分别连接悬浮力观测模块、电压变换模块和悬浮力与电压变换模块的输入端，通过当前采样周期的磁链和转矩预测出下一周期磁链、转矩和悬浮力值，避免延时带来的问题。

Description

无轴承永磁同步电机预测控制器及构造方法

技术领域

本发明属于电力传动控制设备的技术领域，涉及一种多变量、非线性、强耦合、的无轴承永磁同步电机，无轴承永磁同步电机在密封泵、高速精密机械加工、航空航天、飞轮储能、生命科学、真空技术等领域均具有广泛的应用，适用于多种无轴承电机的高性能控制，通过对电机转矩和悬浮力进行解耦控制，能够实现电机转子稳定悬浮和运行。

背景技术

要实现无轴承永磁同步电机的精确控制，首先要对转矩与悬浮力进行解耦。采用矢量控制策略能够实现转矩与悬浮力的精确解耦，但同时也给控制系统带来大量的复杂的矢量变换计算，不利于控制系统性能的提高。直接转矩控制(DTC)有着算法简单、鲁棒性高等优点，能够弥补矢量控制策略存在的不足，它摒弃了矢量控制中系统过多依赖于电机参数以及矢量坐标变换的复杂性使得实际的控制效果难以达到理论分析的结果的控制思想，采用定子磁链定向和瞬时空间矢量理论，直接对电机定子磁链及转矩进行控制，是一种高性能变频调速控制方法，具有对电机模型参数的依赖程度小的优点。

传统的直接转矩和直接悬浮力控制策略是：定子磁链与电磁转矩是根据当前采样点的电压与电流值计算得到，并根据DTC策略得到相应的目标电压矢量，理想情况下该电压矢量在采样点k时刻施加，在k+1时刻使转矩与磁链达到给定值，但由于采样及数字计算的延迟，期望电压矢量实际上是在一个采样周期后，即k+l时刻施加，在k+2时刻才使转矩与磁链达到给定值，这使得转矩与磁链的响应始终滞后一个开关周期，降低了电机系统的运行性能。因此，传统的直接转矩和直接悬浮力控制策略中，由于采样和计算的延时，使目标电压矢量只能在下一个采样周期才能作用于系统，这样就使磁链、转矩和悬浮力的响应始终滞后一个周期，造成磁链、转矩和悬浮力的波动变大。

发明内容

本发明的目的是为了解决无轴承永磁同步电机采用传统的直接转矩和直接悬浮力控制存在的问题，提出一种新的无轴承永磁同步电机预测控制器及该预测控制器的构造方法，预测控制器能预测无轴承永磁同步电机下一时刻的磁链、转矩和悬浮力，减小由于采样和数字计算的延迟造成的控制误差，削弱控制延时造成转矩、磁链和悬浮力的波动。

本发明无轴承永磁同步电机预测控制器采用的技术方案是：包含电流传感器和Clack变换器，电流传感器的输出端连接Clark变换器的输入端，Clark变换器的输出端分别连接转矩绕组磁链与转矩初始观测模块以及悬浮力绕组磁链初始观测模块的输入端，悬浮力绕组磁链初始观测模块的输出端连接悬浮力绕组磁链预测模块的输入端，悬浮力绕组磁链预测模块的输出端连接悬浮力观测模块的输入端；所述转矩绕组磁链与转矩初始观测模块的输出端分别连接转矩绕组磁链与转矩预测模块和转矩与电压变换模块的输入端，转矩绕组磁链与转矩预测模块的输出端分别连接悬浮力观测模块、电压变换模块和悬浮力与电压变换模块的输入端，转矩与电压变换模块的输出端依次串接第一SVPWM和第一电压型逆变器；悬浮力与电压变换模块的输出端依次串接第二SVPWM和第二电压型逆变器。

进一步地，电流传感器检测到悬浮力绕组上的电流i_2a、i_2b、i_2c和转矩绕组上的电流i_1a、i_1b、i_1c，与第一电压型逆变器的直流电压U_DC1和第二电压型逆变器的直流电压U_DC2一起输入Clark变换器，Clark变换器输出α-β坐标系下的悬浮力绕组上的电流i_2α、i_2β、电压u_2α、u_2β和转矩绕组上的电流i_1α、i_1β、电压u_1α、u_1β；电流i_2α、i_2β和电压u_2α、u_2β输入悬浮力绕组磁链初始观测模块，悬浮力绕组磁链初始观测模块初始观测下的悬浮力绕组磁链幅值ψ_s2(k)及相位λ，悬浮力绕组磁链预测模块输出下一时刻的悬浮力绕组磁链幅值ψ_s2(k+1)及相位λ；电流i_1α、i_1β和电压u_1α、u_1β输入到转矩绕组磁链与转矩初始观测模块，转矩绕组磁链与转矩初始观测模块输出初始观测下的初始转矩T_e(k)、磁链ψ_s1(k)和相位角θ，转矩绕组磁链与转矩预测模块输出预测的下一时刻的转矩绕组磁链幅值ψ_s1(k+1)、转矩绕组合成气隙磁链幅值ψ_m1(k+1)、相位μ和预测的下一时刻转矩T_e(k+1)，转矩绕组磁链幅值ψ_s1(k+1)输入转矩与电压变换模块，转矩绕组合成气隙磁链幅值ψ_m1(k+1)及相位μ输入悬浮力观测模块和悬浮力与电压变换模块，悬浮力观测模块输出悬浮力的值F_α(k+1)、F_β(k+1)。

本发明无轴承永磁同步电机预测控制器的构造方法采用的技术方案是：首先构造转矩绕组磁链与转矩初始观测模块，再依次构造转矩绕组磁链与转矩预测模块、转矩与电压变换模块、构造悬浮力与电压变换模块、悬浮力观测器、悬浮力绕组初始磁链观测模块和悬浮力绕组磁链预测模块，最后由所构造的7个模块与第一、第二电压型逆变器、第一、第二SVPWM以及两个PID调节器、两个PI调节器一起共同构成无轴承永磁同步电机预测控制器。

本发明采用上述技术方案后的优点在于：

1、本发明采用转矩和悬浮力直接控制来实现无轴承永磁同步电机的高性能控制，将无轴承永磁同步电机复杂非线性强耦合系统的转矩和转子径向悬浮力控制问题转化为简单的直接转矩和直接悬浮力控制系统，很容易实现电机转矩和悬浮力之间的直接控制。具有控制结构简单，转速响应快，优良的动静态性能，对电机参数表现出强的鲁棒性等优点，摆脱了传统矢量控制过多依赖于电机参数、复杂的矢量坐标变换、电流跟踪型逆变器开关频率高的缺点。也摆脱了处于开环状态的悬浮力矢量控制对悬浮力控制精度和动态响应性能的限制。

2、本发明无轴承永磁同步电机预测控制器采用预测直接转矩和预测直接悬浮力控制策略，能通过当前采样周期的磁链和转矩预测出下一周期磁链、转矩和悬浮力值，利用初始观测的转矩、磁链、悬浮力的值，得出相应的目标电压矢量，在下一周期就能施加合适的电压矢量。避免延时带来的转矩、悬浮力脉动和磁链纹波等问题，有效弥补了时间延时对系统性能的影响。同时也保留了传统SVM调制直接转矩和直接悬浮力控制优良的动态特性，提高了系统性能和鲁棒性。

附图说明

图1是本发明无轴承永磁同步电机预测控制器的结构框图；

图2是图1中的转矩绕组磁链与转矩初始观测模块1的构造原理图；

图3是图1中的转矩绕组磁链与转矩预测模块2的构造原理图；

图4是图1中的转矩与电压变换模块3的构造原理图；

图5是图1中的悬浮力与电压变换模块4的构造原理图；

图6是图1中的悬浮力观测模块5的构造原理图；

图7是图1中的悬浮力绕组初始磁链观测模块6的构造原理图；

图8是图1中的悬浮力绕组磁链预测模块7的构造原理图。

图中：1.转矩绕组磁链与转矩初始观测模块；2.转矩绕组磁链与转矩预测模块；3.转矩与电压变换模块；4.悬浮力与电压变换模块；5.悬浮力观测模块；6.悬浮力绕组磁链初始观测模块；7.悬浮力绕组磁链预测模块；8.永磁同步电机预测控制器；9.Clack变换器；10、11.SVPWM；12、13.PI调节器；14、15.PID调节器；16.光电编码器；17、18.电涡流位移传感器；19.电压型逆变器；20.电压型逆变器；21.电流传感器；22.预测直接悬浮力控制器；23.预测直接转矩控制器。

具体实施方式

参见图1，本发明无轴承永磁同步电机预测控制器8包含电流传感器21、Clack变换器9、预测直接悬浮力控制器22和预测直接转矩控制器23。

其中，预测直接悬浮力控制器22由电压型逆变器19、SVPWM11(即空间矢量脉宽调制模块11)、悬浮力与电压变换模块4、悬浮力绕组磁链初始观测模块6、悬浮力绕组磁链预测模块7、悬浮力观测模块5以及PID调节器14和PID调节器15组成。

预测直接转矩控制器23由转矩绕组磁链与转矩初始观测模块1、转矩绕组磁链与转矩预测模块2、电压型逆变器20、SVPWM10、转矩与电压变换模块3以及PI调节器12和PI调节器13组成。

电流传感器21检测无轴承永磁同步电机的悬浮力绕组上的电流i_2a、i_2b、i_2c和转矩绕组上的电流i_1a、i_1b、i_1c，电流传感器21的输出端连接Clark变换器9的输入端，Clark变换器9的输出端分别连接转矩绕组磁链与转矩初始观测模块1以及悬浮力绕组磁链初始观测模块6的输入端。同时将电压型逆变器20的直流电压U_DC1和电压型逆变器19的直流电压U_DC2也输入Clark变换器9。

Clark变换器9对输入信号进行处理，悬浮力绕组上的电流i_2a、i_2b、i_2c生成α-β坐标系下的电流i_2α、i_2β、电压u_2α、u_2β，转矩绕组上的电流i_1a、i_1b、i_1c生成α-β坐标系下的电流i_1α、i_1β、电压u_1α、u_1β。Clark变换器9输出电流i_2α、i_2β、电压u_2α、u_2β、电流i_1α、i_1β、电压u_1α、u_1β信号。

其中的电流i_2α、i_2β、电压u_2α、u_2β信号输入到悬浮力绕组磁链初始观测模块6，得到初始观测下的悬浮力绕组磁链幅值ψ_s2(k)及相位λ。悬浮力绕组磁链初始观测模块6的输出端连接悬浮力绕组磁链预测模块7的输入端，悬浮力绕组磁链预测模块7的输出端连接悬浮力观测模块5。浮力绕组磁链幅值ψ_s2及相位λ经过悬浮力绕组磁链预测模块7，得到观测出的下一时刻的悬浮力绕组磁链幅值ψ_s2(k+1)及相位λ。将悬浮力绕组磁链幅值ψ_s2(k+1)及相位λ输入到悬浮力观测模块5，将转矩绕组磁链与转矩预测模块2输出的预测出的转矩绕组合成气隙磁链幅值ψ_m1(k+1)及相位μ也输入悬浮力观测模块5，悬浮力观测模块5根据输入的悬浮力绕组磁链幅值ψ_s2(k+1)、相位λ以及转矩绕组合成气隙磁链幅值ψ_m1(k+1)、相位μ计算出预测实时的悬浮力的值F_α(k+1)和F_β(k+1)并输出。

其中的电流i_1α、i_1β、、电压u_1α、u_1β信号输入到转矩绕组磁链与转矩初始观测模块1，转矩绕组磁链与转矩初始观测模块1的输出端分别连接转矩绕组磁链与转矩预测模块2和转矩与电压变换模块3的输入端，转矩绕组磁链与转矩预测模块2的输出端分别连接悬浮力观测模块5、电压变换模块3和悬浮力与电压变换模块4的输入端，电压变换模块3的输出端连接依次串接SVPWM10和电压型逆变器20，电压型逆变器20控制无轴承永磁同步电机的。转矩绕组磁链与转矩初始观测模块1对输入的电流i_1α、i_1β、、电压u_1α、u_1β信号进行处理，得到初始观测下的初始转矩T_e(k)、磁链ψ_s1(k)和相位角θ，将初始转矩T_e(k)、磁链ψ_s1(k)与电压u_1α、u_1β输入到转矩绕组磁链与转矩预测模块2，生成预测的下一时刻的转矩绕组磁链幅值ψ_s1(k+1)、转矩绕组合成气隙磁链幅值ψ_m1(k+1)、相位μ和预测的下一时刻转矩T_e(k+1)。其中，转矩绕组合成气隙磁链幅值ψ_m1(k+1)和相位μ分别输入悬浮力观测模块5和悬浮力与电压变换模块4。

采用光电编码器16检测无轴承永磁同步电机转子的实际实际转速ω，将实际转速ω与转速指令值ω^*的比较得到的转速差值输入到PI调节器12中，调制后生成转矩指令值然后与转矩绕组磁链与转矩预测模块2输出的下一时刻转矩T_e(k+1)比较得到差值，将该转矩差值输入到PI调节器13中进行调制得到转矩绕组磁链相位角增量Δδ；将磁链相位角增量Δδ、磁链指令值下一时刻的转矩绕组磁链幅值ψ_s1(k+1)和相位角θ输入到转矩与电压变换模块3中，电压变换模块3生成两相静止坐标下的控制转矩的电压指令值和将电压指令值和变量输入到空间脉宽调制模块SVPWM10中得到逆变器开关信号S_1a、S_1b、S_1c，再输入给电压型逆变器20，电压型逆变器20控制无轴承永磁同步电机的转矩绕组。其中U_DC1为电压型逆变器20的直流电压，也被称为母线电压。

采用电涡流位移传感器17和电涡流位移传感器18分别检测无轴承永磁同步电机转子的x轴与y轴的实时径向位移值x、y，将实时径向位移值x、y分别与转子位移指令值x^*、y^*分别对应地进行比较得到位移差值，将两组位移差值各输入到对应PID调节器中，即x轴位移差值输入到PID调节器14中，y轴位移差值输入到PID调节器15中，PID调节器14调制后生成x轴悬浮力指令值PID调节器15调制后生成y轴悬浮力指令值两个浮力指令值 分别与悬浮力观测模块5输出的实时悬浮力F_α(k+1)和F_β(k+1)进行比较得到对应的力的差值ΔF_α和ΔF_β。将力的差值信号ΔF_α和ΔF_β以及预测出的转矩绕组合成气隙磁链幅值ψ_m1(k+1)和相位μ输入到悬浮力与电压变换模块4，最终生成改变气隙平衡磁场的电压指令值和悬浮力与电压变换模块4的输出端依次串接空间脉宽调制模块SVPWM11和电压型逆变器19。然后将电压指令值和变量输入到SVPWM11中得到逆变器开关信号S_2a、S_2b、S_2c，再输入给电压型逆变器19，此电压型逆变器19控制无轴承永磁同步电机的悬浮力绕组。

无轴承永磁同步电机预测控制器8的具体构造方案分以下8步：

1、如图2所示，首先构造转矩绕组磁链与转矩初始观测模块1。初始观测下的转矩T_e(k)、磁链ψ_s1(k)和相位角θ的产生是由以下几个变量得到的：α-β坐标系下转矩绕组的电流i_1α、i_1β，电压u_1α、u_1β、定子绕组电阻R_s、转矩绕组极对数n_p，根据以下表达式构造转矩绕组磁链与转矩初始观测模块1：

$\left\{\begin{array}{c}{\mathrm{\ψ}}_{s1\mathrm{\α}}=\∫(u_{1\mathrm{\α}}-R_s{i}_{1\mathrm{\α}})dt\\ {\mathrm{\ψ}}_{s1\mathrm{\β}}=\∫(u_{1\mathrm{\β}}-R_s{i}_{1\mathrm{\β}})dt\\ \left|{\mathrm{\ψ}}_{s1}\right(k\left)\right|=\sqrt{{\mathrm{\ψ}}_{s1\mathrm{\α}}^2+{\mathrm{\ψ}}_{s1\mathrm{\β}}^2}\\ \mathrm{\θ}=\arctan ({\mathrm{\ψ}}_{s1\mathrm{\β}}/{\mathrm{\ψ}}_{s1\mathrm{\α}})\end{array}\right..$

T_e(k)＝3n_p(ψ_s1ai_1β-ψ_s1βi_1α)/2

2、如图3所示，构造转矩绕组磁链与转矩预测模块2。为了获取预测的下一时刻的转矩T_e(k+1)、转矩绕组磁链幅值ψ_s1(k+1)、转矩绕组合成气隙磁链幅值ψ_m1(k+1)和相位μ。需要由转矩绕组磁链与转矩初始观测模块1提供的初始观测的转矩T_e(k)、磁链ψ_s1(k)以及α-β坐标系下转矩绕组的电压u_1α、u_1β，下一时刻的磁链由初始磁链经过一个采样周期T_s推算得到；下一时刻d-q轴的转矩绕组电流i_1d(k+1)、i_1q(k+1)由初始转矩绕组的电流i_1d(k)、i_1q(k)、电压u_1d、u_1q、初始转速ω(k)、绕组电感L_d、L_q、永磁磁链ψ_f经过一个采样周期T_s推算得到；这样就可以根据磁链和电流的关系，估算出下一时刻的转矩T_e(k+1)。其计算过程如下：

$\left\{\begin{array}{c}{\mathrm{\ψ}}_{s1\mathrm{\α}}(k+1)={\mathrm{\ψ}}_{s1\mathrm{\α}}\left(k\right)+\[u_{1\mathrm{\α}}-R_s{i}_{1\mathrm{\α}}\left(k\right)\]T_s\\ {\mathrm{\ψ}}_{s1\mathrm{\β}}(k+1)={\mathrm{\ψ}}_{s1\mathrm{\β}}\left(k\right)+\[u_{1\mathrm{\β}}-R_s{i}_{1\mathrm{\β}}\left(k\right)\]T_s\\ \left|{\mathrm{\ψ}}_{s1}(k+1)\right|=\sqrt{{\mathrm{\ψ}}_{s1\mathrm{\α}}{(k+1)}^2+{\mathrm{\ψ}}_{s1\mathrm{\β}}{(k+1)}^2}\end{array}\right.$

$\left\{\begin{array}{c}{i}_{1d}(k+1)=i_{1d}\left(k\right)+\frac{u_{1d}-R_s{i}_{id}\left(k\right)+\mathrm{\ω}\left(k\right)L_q{i}_{1q}\left(k\right)}{L_d}{T}_s\\ i_{1q}(k+1)=i_{1q}\left(k\right)+\frac{u_{1q}-R_s{i}_{1q}\left(k\right)-\mathrm{\ω}\left(k\right){\mathrm{\ψ}}_f-\mathrm{\ω}\left(k\right)L_d{i}_{1d}\left(k\right)}{L_q}{T}_s\end{array}\right.,$

$T_e(k+1)=\frac{3}{2}{n}_p\[{\mathrm{\ψ}}_{s1\mathrm{\α}}(k+1)i_{1\mathrm{\β}}(k+1)-{\mathrm{\ψ}}_{s1\mathrm{\β}}(k+1)i_{1\mathrm{\α}}(k+1)\]$

转矩绕组合成气隙磁链幅值ψ_m1(k+1)由预测的转矩绕组的磁链减去漏感L_1l产生的磁链，由下式可得：

$\left\{\begin{array}{c}{\mathrm{\ψ}}_{m1\mathrm{\α}}(k+1)={\mathrm{\ψ}}_{s1\mathrm{\α}}(k+1)-L_{1l}{i}_{1\mathrm{\α}}(k+1)\\ {\mathrm{\ψ}}_{m1\mathrm{\β}}(k+1)={\mathrm{\ψ}}_{s1\mathrm{\β}}(k+1)-L_{1l}{i}_{1\mathrm{\β}}(k+1)\\ \left|{\mathrm{\ψ}}_{m1}\right(k+1\left)\right|=\sqrt{{\mathrm{\ψ}}_{m1\mathrm{\α}}^2(k+1)+{\mathrm{\ψ}}_{m1\mathrm{\β}}^2(k+1)}\\ \mathrm{\μ}=\arctan ({\mathrm{\ψ}}_{m1\mathrm{\β}}\left(k+1\right)/{\mathrm{\ψ}}_{m1\mathrm{\α}}\left(k+1\right))\end{array}\right..$

3、如图4所示，构造转矩与电压变换模块3。转矩与电压变换模块3由转矩绕组预测磁链ψ_s1(k+1)与相位θ、给定参考磁链和磁链相位角增量Δδ，得到两相静止坐标下的控制转矩的电压指令值和按下式构造转矩与电压变换模块3：

$\left\{\begin{array}{c}{u}_{s\mathrm{\α}}^{*}=R_s{i}_{1\mathrm{\α}}(k+1)+\frac{\left|{\mathrm{\ψ}}_{s1}^{*}\right|\cos (\mathrm{\θ}+\mathrm{\Δ}\mathrm{\δ})-\left|{\mathrm{\ψ}}_{s1}(k+1)\right|\cos \mathrm{\θ}}{T_s}\\ u_{s\mathrm{\β}}^{*}=R_s{i}_{1\mathrm{\β}}(k+1)+\frac{\left|{\mathrm{\ψ}}_{s1}^{*}\right|\sin (\mathrm{\θ}+\mathrm{\Δ}\mathrm{\δ})-\left|{\mathrm{\ψ}}_{s1}(k+1)\right|\sin \mathrm{\θ}}{T_s}\end{array}\right.,$

其中，i_1α(k+1)、i_1β(k+1)是α-β轴的转矩绕组电流，R_s是定子绕组电阻R_s。

4、如图5所示，构造悬浮力与电压变换模块4。悬浮力与电压变换模块4由悬浮力指令值和和预测值F_α(k+1)和F_β(k+1)的比较差值ΔF_α和ΔF_β、转矩绕组合成气隙磁链幅值ψ_m1(k+1)和相位μ计算得到：

$\left\{\begin{array}{c}{u}_{2\mathrm{\α}}^{*}=R_s{i}_{2\mathrm{\α}}+K_M^{-1}{\mathrm{\ψ}}_{m1}^{-1}(k+1)(\mathrm{\Δ}{F}_{\mathrm{\α}}\cos \mathrm{\μ}+\mathrm{\Δ}{F}_{\mathrm{\β}}\sin \mathrm{\μ})\\ u_{2\mathrm{\β}}^{*}=R_s{i}_{2\mathrm{\β}}+K_M^{-1}{\mathrm{\ψ}}_{m1}^{-1}(k+1)(\mathrm{\Δ}{F}_{\mathrm{\α}}\sin \mathrm{\μ}-\mathrm{\Δ}{F}_{\mathrm{\β}}\cos \mathrm{\μ})\end{array}\right.,$

其中，i_2α、i_2β是α-β坐标系下的电流。

5、如图6所示，构造悬浮力观测器5。由于输入的磁链均为预测的下一时刻的悬浮力绕组磁链幅值ψ_s2(k+1)和相位λ、转矩绕组合成气隙磁链幅值ψ_m1(k+1)和相位μ计算得到，所以该悬浮力的值也为预测的下一时刻的值F_α(k+1)和F_β(k+1)。这与x轴与y轴位移的指令值与测量值的差经PID调制后得到的悬浮力指令值和进行比较得到差值ΔF_α和ΔF_β，采用下式构造出悬浮力观测器5：

$\left\{\begin{array}{c}{F}_{\mathrm{\α}}(k+1)=K_M{\mathrm{\ψ}}_{m1}(k+1){\mathrm{\ψ}}_{s2}(k+1)cos(\mathrm{\λ}-\mathrm{\μ})\\ F_{\mathrm{\β}}(k+1)=K_M{\mathrm{\ψ}}_{m1}(k+1){\mathrm{\ψ}}_{s2}(k+1)sin(\mathrm{\λ}-\mathrm{\μ})\end{array}\right.,$

其中K_M的表达式为：式中：p₁、p₂分别为转矩绕组磁场与悬浮力绕组磁场等效极对数；l为电机转子铁心有效长度；r为转子半径；n₁、n₂为转矩绕组和悬浮力绕组每相串联的有效匝数；μ₀为空气磁导率。

6、如图7所示，构造悬浮力绕组初始磁链观测模块6。电流传感器21检测到悬浮力绕组上的电流i_2a、i_2b、i_2c，加上电压型逆变器19的直流电压U_DC2，经过Clark变换器9生成α-β坐标系下的电流电压i_2α、i_2β、u_2α、u_2β，将这些电流电压信号输入到该模块中，经如下推导：

$\left\{\begin{array}{c}{\mathrm{\ψ}}_{s2\mathrm{\α}}\left(k\right)=\∫(u_{2\mathrm{\α}}-R_s{i}_{2\mathrm{\α}})dt\\ {\mathrm{\ψ}}_{s2\mathrm{\β}}\left(k\right)=\∫(u_{2\mathrm{\β}}-R_s{i}_{2\mathrm{\β}})dt\\ \left|{\mathrm{\ψ}}_{s2}\right(k\left)\right|=\sqrt{{\mathrm{\ψ}}_{s2\mathrm{\α}}^2\left(k\right)+{\mathrm{\ψ}}_{s2\mathrm{\β}}^2\left(k\right)}\\ \mathrm{\λ}=\arctan ({\mathrm{\ψ}}_{s2\mathrm{\β}}\left(k\right)/{\mathrm{\ψ}}_{s2\mathrm{\α}}\left(k\right))\end{array}\right.,$

得到初始观测下的悬浮力绕组磁链幅值ψ_s2及相位λ。

7、如图8所示，构造悬浮力绕组磁链预测模块7。悬浮力绕组初始磁链观测模块6得出的悬浮力绕组磁链幅值ψ_s2(k)在α-β坐标下的磁链ψ_s2α(k)和ψ_s2β(k)及相位λ作为输入量，加上定子绕组电阻R_s、采样周期T_s，可以预测出下一时刻的悬浮力绕组磁链幅值ψ_s2(k+1)：

$\left\{\begin{array}{c}{\mathrm{\ψ}}_{s2\mathrm{\α}}(k+1)={\mathrm{\ψ}}_{s2\mathrm{\α}}\left(k\right)+(u_{2\mathrm{\α}}-R_s{i}_{2\mathrm{\α}})T_s\\ {\mathrm{\ψ}}_{s2\mathrm{\β}}(k+1)={\mathrm{\ψ}}_{s2\mathrm{\β}}\left(k\right)+(u_{2\mathrm{\β}}-R_s{i}_{2\mathrm{\β}})T_s\\ \left|{\mathrm{\ψ}}_{s2}(k+1)\right|=\sqrt{{\mathrm{\ψ}}_{s2\mathrm{\α}}^2(k+1)+{\mathrm{\ψ}}_{s2\mathrm{\β}}^2(k+1)}\end{array}\right..$

8、最终由上述7个大模块，再与电压型逆变器20、电压型逆变器19、SVPWM10、SVPWM11以及两个PID调节器14、15、两个PI调节器12、13一起共同组成完整的无轴承永磁同步电机预测控制器8，如图1所示，即预测直接转矩和悬浮力直接控制器，可根据不同的控制要求采用不同的硬件或软件来实现。

本发明所述的无轴承永磁同步电机预测控制器8对电机控制时，对于转矩部分，采用Clark变换器9实时计算电机的转矩绕组相电流与电压型逆变器直流电压，然后转矩绕组磁链与转矩初始观测模块1检测出磁链和转矩的初始值，经由转矩绕组磁链和转矩预测模块2来获取直接转矩和直接悬浮力控制所需的无轴承永磁同步电机绕组磁链、转矩和悬浮力的相关信息；再将所预测的绕组实时磁链应用于悬浮力观测模块5、转矩与电压变换模块3和悬浮力与电压变换模块4。其次将转速指令值和实时转速进行比较得到转速差值，通过PI调制后生成转矩指令值。然后将转矩指令值与转矩实时预测值进行比较得到转矩差值，通过PI调制后生成转矩绕组磁链相位角增量。再根据磁链幅值指令值、实时绕组转矩分量预测磁链幅值和相位在转矩与电压变换模块3中生成电压静止坐标指令值，将得到两组电压静止坐标指令值通入空间矢量脉宽调制模块SVPWM10调制后得到电压型逆变器开关信号，驱动电压型逆变器实现无轴承永磁同步电机转矩的预测控制。对于悬浮力部分，采用同一个Clark变换器9实时计算悬浮力绕组定子相电流和电压型逆变器直流电压，然后悬浮力绕组磁链初始观测模块6检测出磁链的初始值，通过悬浮力绕组磁链预测模块7获取下一时刻的磁链信息。悬浮力观测模块5将预测悬浮力绕组磁链的幅值和相位，结合预测转矩绕组合成气隙磁链的幅值和相位实时计算出预测悬浮力静止坐标分量的值。其次将转子位置指令值与电涡流位移传感器17、18获得的实际位移值分别进行比较得到位移差值。再经PID调制后生成悬浮力指令值，然后将悬浮力指令值与实时预测到的下一时刻的悬浮力比较得到悬浮力差值通入悬浮力与电压变换模块4生成电压静止坐标指令值。将得到两组电压静止坐标指令值通入空间矢量脉宽调制模块SVPWM11调制后得到电压型逆变器19的开关信号，驱动逆变器实现无轴承永磁同步电机悬浮力的预测控制。该种预测控制的结构简单，响应速度快，确保电机转子悬浮的稳定性和运行的鲁棒性。通过当前采样周期的磁链和转矩值来预测出下一周期磁链和转矩值，得出相应的目标电压矢量，在下一周期就能施加合适的电压矢量，避免延时带来的纹波问题。既可使无轴承永磁同步电机具有优良的动、静态控制性能，抗电机参数变化及抗负载扰动能力强，又能有效地提高无轴承永磁同步电机的各项控制性能指标，如动态响应速度、稳态跟踪精度。

Claims (6)

1.一种无轴承永磁同步电机预测控制器，包含电流传感器(21)和Clack变换器(9)，其特征是：电流传感器(21)的输出端连接Clark变换器(9)的输入端，Clark变换器(9)的输出端分别连接转矩绕组磁链与转矩初始观测模块(1)以及悬浮力绕组磁链初始观测模块(6)的输入端，悬浮力绕组磁链初始观测模块(6)的输出端连接悬浮力绕组磁链预测模块(7)的输入端，悬浮力绕组磁链预测模块(7)的输出端连接悬浮力观测模块(5)的输入端；所述转矩绕组磁链与转矩初始观测模块(1)的输出端分别连接转矩绕组磁链与转矩预测模块(2)和转矩与电压变换模块(3)的输入端，转矩绕组磁链与转矩预测模块(2)的输出端分别连接悬浮力观测模块(5)、电压变换模块(3)和悬浮力与电压变换模块(4)的输入端，转矩与电压变换模块(3)的输出端依次串接第一SVPWM(10)和第一电压型逆变器(20)；悬浮力与电压变换模块(4)的输出端依次串接第二SVPWM(11)和第二电压型逆变器(19)。

2.根据权利要求1所述无轴承永磁同步电机预测控制器，其特征是：电流传感器(21)检测到悬浮力绕组上的电流i_2a、i_2b、i_2c和转矩绕组上的电流i_1a、i_1b、i_1c，与第一电压型逆变器(20)的直流电压U_DC1和第二电压型逆变器(19)的直流电压U_DC2一起输入Clark变换器(9)，Clark变换器(9)输出α-β坐标系下的悬浮力绕组上的电流i_2α、i_2β、电压u_2α、u_2β和转矩绕组上的电流i_1α、i_1β、电压u_1α、u_1β；电流i_2α、i_2β和电压u_2α、u_2β输入悬浮力绕组磁链初始观测模块(6)，悬浮力绕组磁链初始观测模块(6)初始观测下的悬浮力绕组磁链幅值ψ_s2(k)及相位λ，悬浮力绕组磁链预测模块(7)输出下一时刻的悬浮力绕组磁链幅值ψ_s2(k+1)及相位λ；电流i_1α、i_1β和电压u_1α、u_1β输入到转矩绕组磁链与转矩初始观测模块(1)，转矩绕组磁链与转矩初始观测模块(1)输出初始观测下的初始转矩T_e(k)、磁链ψ_s1(k)和相位角θ，转矩绕组磁链与转矩预测模块(2)输出预测的下一时刻的转矩绕组磁链幅值ψ_s1(k+1)、转矩绕组合成气隙磁链幅值ψ_m1(k+1)、相位μ和预测的下一时刻转矩T_e(k+1)，转矩绕组磁链幅值ψ_s1(k+1)输入转矩与电压变换模块(3)，转矩绕组合成气隙磁链幅值ψ_m1(k+1)及相位μ输入悬浮力观测模块(5)和悬浮力与电压变换模块(4)，悬浮力观测模块(5)输出悬浮力的值F_α(k+1)、F_β(k+1)。

3.根据权利要求2所述无轴承永磁同步电机预测控制器，其特征是：光电编码器检测电机转子的实际转速ω，转速ω与转速指令值ω^*比较得到转速差值输入到第一PI调节器(12)，生成转矩指令值T_e ^*与所述预测的下一时刻转矩T_e(k+1)比较得到转矩差值，该转矩差值输入第二PI调节器(13)得到的转矩绕组磁链相位角增量Δδ并输入转矩与电压变换模块(3)。

4.根据权利要求2所述无轴承永磁同步电机预测控制器，其特征是：电涡流位移传感器检测电机转子的x轴、y轴的实时径向位移值x、y，实时径向位移值x、y与转子位移指令值x^*、y^*分别对应地进行比较得到位移差值，将两组位移差值各输入到对应的PID调节器中生成x轴、y轴的悬浮力指令值悬浮力指令值分别与所述悬浮力的值F_α(k+1)、F_β(k+1)比较得到对应的力的差值ΔF_α、ΔF_β，力的差值ΔF_α、ΔF_β输入到悬浮力与电压变换模块(4)。

5.一种如权利要求1所述无轴承永磁同步电机预测控制器的构造方法，其特征是：首先构造转矩绕组磁链与转矩初始观测模块(1)，再依次构造转矩绕组磁链与转矩预测模块(2)、转矩与电压变换模块(3)、悬浮力与电压变换模块(4)、悬浮力观测器(5)、悬浮力绕组初始磁链观测模块(6)和悬浮力绕组磁链预测模块(7)，最后由所构造的7个模块与第一电压型逆变器(20)、第二电压型逆变器(19)、第一SVPWM(10)、第二SVPWM(11)以及两个PID调节器、两个PI调节器一起共同构成无轴承永磁同步电机预测控制器。

6.根据权利要求1所述无轴承永磁同步电机预测控制器的构造方法，其特征是：

根据公式构造转矩绕组磁链与转矩初始观测模块(1)；

T_e(k)＝3n_p(ψ_s1ai_1β-ψ_s1βi_1α)/2

$\left\{\begin{array}{c}{\mathrm{\ψ}}_{s1\mathrm{\α}}(k+1)={\mathrm{\ψ}}_{s1\mathrm{\α}}\left(k\right)+\[u_{1\mathrm{\α}}-R_s{i}_{1\mathrm{\α}}\left(k\right)\]T_s\\ {\mathrm{\ψ}}_{s1\mathrm{\β}}(k+1)={\mathrm{\ψ}}_{s1\mathrm{\β}}\left(k\right)+\[u_{1\mathrm{\β}}-R_s{i}_{1\mathrm{\β}}\left(k\right)\]T_s\\ \left|{\mathrm{\ψ}}_{s1}(k+1)\right|=\sqrt{{\mathrm{\ψ}}_{s1\mathrm{\α}}{(k+1)}^2+{\mathrm{\ψ}}_{s1\mathrm{\β}}{(k+1)}^2}\end{array}\right.$

根据公式和

$T_e(k+1)=\frac{3}{2}{n}_p\[{\mathrm{\ψ}}_{s1\mathrm{\α}}(k+1)i_{1\mathrm{\β}}(k+1)-{\mathrm{\ψ}}_{s1\mathrm{\β}}(k+1)i_{1\mathrm{\α}}(k+1)\]$

构造转矩绕组磁链与转矩预测模块(2)；

根据公式构造转矩与电压变换模块(3)，i_1α(k+1)、i_1β(k+1)是α-β轴的转矩绕组电流，R_s是定子绕组电阻R_s；

根据公式构造悬浮力与电压变换模块(4)；

根据公式构造悬浮力观测器(5)，p₁、p₂分别为转矩绕组磁场与悬浮力绕组磁场等效极对数；l为电机转子铁心有效长度；r为转子半径；n₁、n₂为转矩绕组和悬浮力绕组每相串联的有效匝数；μ₀为空气磁导率；

根据公式构造悬浮力绕组初始磁链观测模块(6)；

根据公式构造悬浮力绕组磁链预测模块(7)。