一种基于永磁同步电机转矩预测控制方法与流程

本发明涉及一种永磁同步电机技术领域，尤其是一种基于永磁同步电机转矩预测控制方法。

背景技术：

现有电机主要分为两类：直流电机和交流电机，根据结构和工作原理的不同，交流电机可以分为交流异步电机和同步电机，随着永磁材料性能的提高和价格的下降，永磁同步电机应用越来越广泛，并且由于永磁同步电机具有以高功率因素、高效率、低能耗、结构简单、体积小、较硬的机械特性、较快的动态响应速度、控制方法简单、运行可靠等优点，所以永磁同步电机使用范围广泛；

现有永磁同步电机的控制方法中矢量控制和直接转矩控制较为成熟，但是直接转矩控制是基于对磁链和转矩的滞环控制，根据滞环比较结果采用查表的方法选择基本电压矢量，这会使转矩和磁链产生较大脉动；

矢量控制转矩响应较慢，含有复杂的坐标变换，控制结构复杂等问题，随着数字信号处理器等控制芯片计算能力的大幅提高，现有的转矩预测控制显然不能满足需求。

技术实现要素：

针对现有技术的不足，本发明提供一种基于永磁同步电机转矩预测控制方法。

本发明的技术方案为：一种基于永磁同步电机转矩预测控制方法，其特征在于，包括以下步骤：

s1)、通过电流传感器获取k时刻永磁同步电机的三相定子电流isa(k)、isb(k)、isc(k)，并且通过电压传感器获取永磁同步电机直流侧电压udc，并且通过编码器获取永磁同步电机的转速ω；

s2)、通过将转速ω与参考转速ω^*做差后经pi调节器输出得到参考转矩指令te^*；

s3)、通过参考转矩指令te^*计算得到参考定子磁链ψs^*，其计算公式如下：

其中，p为极对数，ψf为转子磁链，|ψs^*|为参考定子磁链幅值，ls为定子电感；

s4)、通过计算参考定子磁链相位角θ，利用极坐标系得到参考定子磁链ψs^*在αβ轴下的分量ψsα^*、ψsβ^*，其中，参考定子磁链相位角θ计算公式如下：

θ＝arctan(lste^*/(1.5pψf²))+ωt，

其中，ω为转速，arctan(lste^*/(1.5pψf²))为转矩角，参考定子磁链ψs^*在αβ轴下的分量ψsα^*、ψsβ^*的计算公式如下：

ψsα^*＝|ψs^*|cosθ，

ψsβ^*＝|ψs^*|sinθ；

s5)、根据逆变器的驱动信号sa、sb、sc和直流侧电压udc计算得到k时刻永磁同步电机三相电压usa(k)、usb(k)、usc(k)，计算公式如下：

s6)、通过将永磁同步电机k时刻的三相电压usa(k)、usb(k)、usc(k)和永磁同步电机的三相定子电流isa(k)、isb(k)、isc(k)经过坐标变换后可以得到usα(k)、usβ(k)和isα(k)、isβ(k)，其计算公式如下：

s7)、通过usα(k)、usβ(k)和isα(k)、isβ(k)计算可以得到定子磁链ψs(k)在αβ轴下的分量ψsα(k)、ψsβ(k)，并根据定子磁链ψs(k)的分量ψsα(k)、ψsβ(k)计算得到转矩te(k)，计算公式如下：

其中，rs为定子电阻；

s8)、通过反pack坐标变换可以得到转子磁链ψf在静止坐标系下的分量ψfα和ψfβ，计算公式如下：

s9)、为使每个周期输出三个不同的电压矢量，通过三桥臂的8个开关状态，设定六个非零电压矢量(u1u2u3u4u5u6)和零电压矢量u0或u7，并将六个非零电压矢量(u1u2u3u4u5u6)和零电压矢量u0或u7组成12组电压矢量集，其中每组电压矢量集由两个相邻的非零电压矢量和零电压矢量u0或u7组成，具体分别为：(u0u1u2)、(u0u2u3)、(u0u3u4)、(u0u4u5)、(u0u5u6)、(u0u6u1)、(u7u1u2)、(u7u2u3)、(u7u3u4)、(u7u4u5)、(u7u5u6)、(u7u6u1)，并从这12组电压矢量集中任意选取一组电压矢量集作用于k+1时刻，从而使永磁同步电机快速达到额定转速，为了减少开关损耗，设置每组电压矢量左右对称；

s10)、根据与计算得到下一时刻k+1的定子磁链分量的预测值ψsα(k+1)和ψsβ(k+1)，其中，由于定子电阻rs的阻值很小，可以忽略，计算公式如下：

其中，kαi＝uαi、kβi＝uβi、kαj＝uαj、kβj＝uβj为不同非零电压矢量(ui、uj)在αβ轴下作用的磁链斜率，t1和t2为相对应非零电压矢量(ui、uj)作用的时间长，其中，t1和t2的计算式为：

其中，δψsα(k)＝ψsα(k+1)-ψsα(k)＝ψsα^*(k)-ψsα(k)；δψsβ(k)＝ψsβ(k+1)-ψsβ(k)＝ψsβ^*(k)-ψsβ(k)；另外零电压矢量u0或u7作用的时间为t0＝ts-t1-t2；

s11)、由于转矩含有电流和磁链乘积，从而下一时刻k+1的电流预测值的计算公式如下：

由于可将电流预测值中含电压部分替换，得到：

则可得下一时刻k+1的转矩的预测值te(k+1)为：

其中，由于isα0、isβ0均为k时刻的值，故无需预测下一时刻k+1的电流值；

s12)、通过在价值函数g中加入开关次数f，从而区别零电压矢量u0与u7，并且依次将步骤s10)和s11)中在12组不同电压矢量集作用下的定子磁链分量的预测值ψsα(k+1)和ψsβ(k+1)与转矩的预测值te(k+1)输入到价值函数g中，分别计算得到12个不同的价值函数g值，计算公式如下：

g＝|te(k+1)-te^*(k)|+λ|ψs(k+1)-ψs^*(k)|+f，

其中，λ为加权系数，f为开关次数；

s13)、选取步骤s12)中最小的价值函数gmin值所对应的电压矢量集，以及相对应的电压矢量作用的时间，得到三个不同脉冲信号作用于逆变器，从而控制三个上下管的开通，实现逆变器的控制。

本发明的有益效果为：设计合理、实用性强、操作简单，并且无需计算下一时刻的电流预测值，从而避免了由于坐标变换导致的计算量大、计算过程复杂的问题，从而减少了计算量，大大简化了计算工作，另外，还避免了由于滞环控制导致的转矩和磁链产生脉动较大、控制转矩响应较慢等问题，

附图说明

图1为本发明的流程示意图；

图2为本发明的框架结构图；

图3为本发明中8个电压矢量和驱动信号之间的关系示意图；

图4为本发明中12组不同的电压矢量及其相应电压矢量作用时间的脉冲示意图；

具体实施方式

下面结合附图对本发明的具体实施方式作进一步说明：

如图1和图2所示，一种基于永磁同步电机转矩预测控制方法，其特征在于，包括以下步骤：

s1)、通过电流传感器获取k时刻永磁同步电机的三相定子电流isa(k)、isb(k)、isc(k)，并且通过电压传感器获取永磁同步电机直流侧电压udc，并且通过编码器获取永磁同步电机的转速ω；

s2)、通过将转速ω与参考转速ω^*做差后经pi调节器输出得到参考转矩指令te^*；

s3)、通过参考转矩指令te^*计算得到参考定子磁链ψs^*，其计算公式如下：

其中，p为极对数，ψf为转子磁链，|ψs^*|为参考定子磁链幅值，ls为定子电感；

s4)、通过计算参考定子磁链相位角θ，利用极坐标系得到参考定子磁链ψs^*在αβ轴下的分量ψsα^*、ψsβ^*，其中，参考定子磁链相位角θ计算公式如下：

θ=arctan(lste^*/(1.5pψf²))+ωt，

其中，ω为转速，arctan(lste^*/(1.5pψf²))为转矩角，参考定子磁链ψs^*在αβ轴下的分量ψsα^*、ψsβ^*的计算公式如下：

ψsα^*＝|ψs^*|cosθ，

ψsβ^*＝|ψs^*|sinθ；

s5)、根据逆变器的驱动信号sa、sb、sc和直流侧电压udc计算得到k时刻永磁同步电机三相电压usa(k)、usb(k)、usc(k)，计算公式如下：

s6)、通过将永磁同步电机k时刻的三相电压usa(k)、usb(k)、usc(k)和永磁同步电机的三相定子电流isa(k)、isb(k)、isc(k)经过坐标变换后可以得到usα(k)、usβ(k)和isα(k)、isβ(k)，其计算公式如下：

s7)、通过usα(k)、usβ(k)和isα(k)、isβ(k)计算可以得到定子磁链ψs(k)在αβ轴下的分量ψsα(k)、ψsβ(k)，并根据定子磁链ψs(k)的分量ψsα(k)、ψsβ(k)计算得到转矩te(k)，计算公式如下：

其中，rs为定子电阻；

s8)、通过反pack坐标变换可以得到转子磁链ψf在静止坐标系下的分量ψfα和ψfβ，计算公式如下：

s9)、为使每个周期输出三个不同的电压矢量，通过三桥臂的8个开关状态，设定六个非零电压矢量(u1u2u3u4u5u6)和零电压矢量u0或u7，具体如图3所示，其中，u0(000)表示电压矢量u0的对应的驱动信号为(000)，并将六个非零电压矢量(u1u2u3u4u5u6)和零电压矢量u0或u7组成12组电压矢量集，其中每组电压矢量集由两个相邻的非零电压矢量和零电压矢量u0或u7组成，具体分别为：(u0u1u2)、(u0u2u3)、(u0u3u4)、(u0u4u5)、(u0u5u6)、(u0u6u1)、(u7u1u2)、(u7u2u3)、(u7u3u4)、(u7u4u5)、(u7u5u6)、(u7u6u1)，并从这12组电压矢量集中任意选取一组电压矢量集作用于k+1时刻，从而使永磁同步电机快速达到额定转速，为了减少开关损耗，设置每组电压矢量左右对称；

s10)、根据与计算得到下一时刻k+1的定子磁链分量的预测值ψsα(k+1)和ψsβ(k+1)，其中，由于定子电阻rs的阻值很小，可以忽略，计算公式如下：

其中，kαi＝uαi、kβi＝uβi、kαj＝uαj、kβj＝uβj为不同非零电压矢量(ui、uj)在αβ轴下作用的磁链斜率，t1和t2为相对应非零电压矢量(ui、uj)作用的时间长，其中，t1和t2的计算式为：

其中，δψsα(k)＝ψsα(k+1)-ψsα(k)＝ψsα^*(k)-ψsα(k)；δψsβ(k)＝ψsβ(k+1)-ψsβ(k)＝ψsβ^*(k)-ψsβ(k)；另外零电压矢量u0或u7作用的时间为t0＝ts-t1-t2；

s11)、由于转矩含有电流和磁链乘积，从而下一时刻k+1的电流预测值的计算公式如下：

由于可将电流预测值中含电压部分替换，得到：

则可得下一时刻k+1的转矩的预测值te(k+1)为：

其中，由于isα0、isβ0均为k时刻的值，故无需预测下一时刻k+1的电流值；

s12)、通过在价值函数g中加入开关次数f，从而区别零电压矢量u0与u7，并且依次将步骤s10)和s11)中在12组不同电压矢量集作用下的定子磁链分量的预测值ψsα(k+1)和ψsβ(k+1)与转矩的预测值te(k+1)输入到价值函数g中，分别计算得到12个不同的价值函数g值，计算公式如下：

g＝|te(k+1)-te^*(k)|+λ|ψs(k+1)-ψs^*(k)|+f，

其中，λ为加权系数，f为开关次数；

s13)、选取步骤s12)中最小的价值函数gmin值所对应的电压矢量集，以及相对应的电压矢量作用的时间，得到三个不同脉冲信号作用于逆变器，从而控制三个上下管的开通，实现逆变器的控制，其中，12组不同电压矢量集作用时间和相应的电压矢量集形成的脉冲信号如图4所示，为使驱动信号左右对称，调整作用时间的先后顺序。

上述实施例和说明书中描述的只是说明本发明的原理和最佳实施例，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。