CN111030542A - 一种永磁同步电机预测电流控制方法和装置 - Google Patents

一种永磁同步电机预测电流控制方法和装置 Download PDF

Info

Publication number
CN111030542A
CN111030542A CN201911392796.9A CN201911392796A CN111030542A CN 111030542 A CN111030542 A CN 111030542A CN 201911392796 A CN201911392796 A CN 201911392796A CN 111030542 A CN111030542 A CN 111030542A
Authority
CN
China
Prior art keywords
vector
ref
current
permanent magnet
region
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201911392796.9A
Other languages
English (en)
Other versions
CN111030542B (zh
Inventor
刘平
刘超
黄守道
苗轶如
童林林
陈常乐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hunan Lanhai Electrical Engineering Co ltd
Original Assignee
Hunan University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hunan University filed Critical Hunan University
Priority to CN201911392796.9A priority Critical patent/CN111030542B/zh
Publication of CN111030542A publication Critical patent/CN111030542A/zh
Application granted granted Critical
Publication of CN111030542B publication Critical patent/CN111030542B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/22Current control, e.g. using a current control loop
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/14Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P2205/00Indexing scheme relating to controlling arrangements characterised by the control loops
    • H02P2205/01Current loop, i.e. comparison of the motor current with a current reference
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P2207/00Indexing scheme relating to controlling arrangements characterised by the type of motor
    • H02P2207/05Synchronous machines, e.g. with permanent magnets or DC excitation

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Ac Motors In General (AREA)

Abstract

本发明公开了一种永磁同步电机预测电流控制方法和装置,该预测电流控制方法中通过对当前时刻永磁同步电机的直流母线电压、三相电流、转子旋转电角速度与转子电角度进行采样,计算得到静止坐标系的参考电压矢量,随后提出了针对基本电压矢量扇区的多区域划分方法,在考虑了过调制情况下选取出了最佳矢量,最后给出了各个基本电压矢量的最佳作用时间的精确简化表达式;故本发明的方法能够改进已有过调制优化方法且减少三角函数的使用,大大减少了计算量,并同时进一步减小了永磁同步电机交直轴电流纹波,减小转矩脉动和定子电流总谐波畸变率,提高了系统的稳态性能。

Description

一种永磁同步电机预测电流控制方法和装置
技术领域
本发明涉及电机控制方法领域,更具体地说,涉及一种永磁同步电机预测电流控制方法和装置。
背景技术
永磁同步电机以其功率密度高、效率高、可靠性好等优点,在电动汽车等电力传动领域得到了广泛应用。电机控制策略的优劣直接决定了驱动系统性能的好坏。随着电气传动产业化的迫切需求,对电驱系统的高效率、高密度与高可靠性要求也日益严苛。当前交流电机控制方式应用较广泛的主要是磁场定向矢量控制(FOC)和直接转矩控制(DTC)。传统FOC虽然具有稳态性能好、转矩脉动小等优点,但由于其电流内环一般基于比例-积分(PI)控制器进行设计,存在积分饱和PI参数整定难,交直轴(dq轴)电流控制相互影响/解耦算法较复杂以及系统的约束不好处理等问题,导致其电流环动态响应能力受限,无法进一步提高,难以在宽速度范围内同时获得较优的动态和稳态性能;DTC虽然动态性能极佳,但其普遍问题是转矩脉动较大,仍然不能很好地满足工业上对交流电机的要求,并且低速性能较差。
模型预测控制(MPC)被广泛认为是继FOC和DTC之后的新一代高性能控制策略,易于集成系统非线性约束和实现对电力电子变换器和电机系统多目标控制,具有动态响应快、电流控制性能好等特点,且不需要对硬件做改动,控制灵活、适用性强,具有非常好的应用前景。与DTC相比,模型预测控制在矢量选择上更加准确有效,它通过对电机状态进行预测,来选取当前时刻的最优电压矢量,从而能够获得更好的稳态性能。
传统模型预测电流控制在一个采样周期内只有一个电压矢量作用,而单矢量预测控制需要遍历所有的基本电压矢量来预测系统下一时刻对应的状态,计算量较大;当矢量作用时间大于采样周期时,参考电压矢量处于过调制区,在过调制区的处理方法为令矢量作用时间等于采样周期,因此单矢量预测控制导致较大的电流脉动。
为了提高系统的稳态性能,无差拍和空间矢量调制的思想被引入。三矢量预测电流控制能有效减小转矩脉动和定子电流总谐波畸变率,但是在三矢量预测控制过调制优化方法中,矢量作用时间表达式包含三角函数,计算较为复杂。相比于二矢量预测控制,三矢量预测控制开关频率增大,变换器的损耗明显增加。总体来说,二矢量预测电流控制有比传统预测电流控制更好的稳态性能,比三矢量预测电流控制更低的开关频率。
在二矢量预测控制中,存在例如“CN201610321625.7-基于无差拍优化与双矢量模型预测的电机控制方法”的专利利用转矩和磁链同时无差拍求得参考电压矢量,利用转子位置信息得到静止坐标系下的参考电压矢量和位置角,根据位置角判断电压矢量所处的扇区,选择相应的备选电压矢量进行计算,得到使目标函数值最小的预测电压矢量;计算得到预测电压矢量和零矢量的作用时间。但是,在CN201610321625.7专利中,当矢量作用时间大于采样周期时,采用的处理方法为将作用时间线性缩小至采样周期以内,仍存在一定误差,且非零矢量作用时间计算较为复杂;在计算矢量作用时间之后需要判定是否过调制,存在计算量大的问题。
此外,在“CN201810787417.5-一种基于几何法的永磁同步电机两矢量预测控制系统及方法”的专利中,其第65-70段以及图2中的确定备选电压矢量和双电压矢量选择的计算中存在大量三角函数,提高了基于占空比的作用时间的计算的复杂度和计算量,其针对过调制的控制方法存在误差,且没对III+区域细分,求作用时间后必须判断t>Ts、t<0(以参考矢量和AC重合为例,当参考矢量长度≥2AC,在BC上的投影
Figure BDA0002345457960000031
此时作用时间大于等于Ts,若t>Ts则令t=Ts,若t<0则令t=0)。此外,CN201810787417.5专利中求作用时间是用投影长度d除以该边长,该专利的图3(b)中d为Vop在BC边上的投影长度,这种求解方法存在较大误差,比如当Vop和AO重合且p在BC上时,d=0,u1和u2某一个作用时间为0,误差为BC/2,而按照现有技术的误差计算方法,其误差等于0,故其求解方法存在较大误差。
综上所述可知,在已有的现有技术中,用于两电平逆变器驱动永磁同步电机的二矢量预测电流控制仍采用传统处理方法;已出现在三矢量预测控制中的过调制优化方法计算依然复杂,对于永磁同步电机预测控制,仍然缺少一种能够同时减小转矩脉动和定子电流总谐波畸变率、提高稳态性能、改进了已有过调制优化方法、减少三角函数的使用以减少电机控制算法计算量的PMSM预测电流控制方法或装置。
发明内容
(一)技术问题
基于上述的技术缺陷,本发明提供一种永磁同步电机预测电流控制方法和装置,其实际上属于一种考虑过调制优化的永磁同步电机二矢量预测电流控制方法,该方法能够进一步减小永磁同步电机交直轴电流纹波,减小转矩脉动和定子电流总谐波畸变率,提高稳态性能,改进简化了已有过调制优化方法,并同时减少三角函数的使用以减少电机控制方法的计算量。
(二)技术方案
本发明提供的一种永磁同步电机预测电流控制方法,其特征在于,该电流控制方法包括如下步骤:
步骤1:对当前时刻永磁同步电机的直流母线电压udc(k)、三相电流ia(k),ib(k),ic(k)、转子旋转电角速度ωe(k)与转子电角度θe(k)进行采样;
步骤2:计算dq轴的参考电压矢量udref和uqref
步骤3:将udref和uqref转换为静止坐标系;根据两相旋转到两相静止的反Park变换,计算出静止坐标系的参考电压矢量uαref,uβref,uαref和uβref分别为静止坐标系下参考电压矢量对应的定子电压实部和虚部分量,两者共同合成参考电压矢量uref
步骤4:计算静止坐标系的参考电压矢量uαref和uβref所在扇区si
步骤5:判断uref在扇区si对应区域并选择最佳矢量,根据已推导的最佳矢量作用时间的表达式直接计算作用时间;该步骤5中具体包括:考虑过调制区和线性调制区,将扇区分成多个区域,选择与参考电压矢量uref相邻的两个非零矢量和零矢量作为备选矢量,选择线性调制区内距离参考电压矢量最近的矢量作为最佳矢量,并根据对应区域内最佳矢量与备选矢量的关系计算备选矢量中各个基本电压矢量的最佳作用时间。
步骤6:应用基本电压矢量的最佳作用时间,以完成对永磁同步电机的控制。
进一步的,所述步骤2中具体包括以下步骤:
步骤2.1:根据采集的三相电流以及转子电角度计算两相旋转坐标系下的dq轴电流id(k)和iq(k):
Figure BDA0002345457960000051
步骤2.2:永磁同步电机电压方程为:
Figure BDA0002345457960000052
Figure BDA0002345457960000053
其中,Rs、ψf分别表示电机定子电阻、转子永磁体磁链;Ld,Ld分别表示旋转坐标系下d轴、q轴上的电感分量;
采用一阶欧拉离散方程:
Figure BDA0002345457960000061
其中Ts为采样周期,由公式(3)-(4)得k+1时刻电流预测值为:
Figure BDA0002345457960000062
Figure BDA0002345457960000063
步骤2.3:通过以下目标函数选取最优电压矢量:
Figure BDA0002345457960000064
其中,
Figure BDA0002345457960000065
分别为d轴和q轴参考电流,电流预测控制的目标是dq轴电流预测值跟随参考电流,理想情况下有
Figure BDA0002345457960000066
Figure BDA0002345457960000067
代入到公式(2)~(3),可计算同步旋转坐标系上d轴、q轴上的参考电压矢量udref,uqref
Figure BDA0002345457960000068
进一步的,所述步骤3中具体满足下式:
Figure BDA0002345457960000069
进一步的,所述步骤4包括以下步骤:
步骤4.1:计算参考电压矢量uref的空间位置角θref
Figure BDA00023454579600000610
公式(10)所求空间位置角θref的取值范围为[-π,π],将其转换到[0,2π]:如果uβref大于等于0,θref不变;如果uβref小于0,θref变为θref+2π;
步骤4.2:对于三相永磁同步电机,对应于两电平三相电压源逆变器,共有8个基本电压矢量,包括2个零矢量与6个非零矢量;非零矢量在空间上从0度开始,逆时针旋转,间隔为π/3,至2π共有6个基本电压矢量均匀分布,依次定义为ui,其中i=1,2,3,…,6,零矢量定义为u0和u7;将8个基本电压矢量所在平面划分为12个扇区,即扇区si,下标i为1到12之间的自然数,每个扇区si的圆心角为30°;在每个扇区非零矢量的一端,作两条与非零矢量夹角为30°的线,将每个扇区划分为3个区域,每个扇区包含区域1-3,离起点O距离最近的区域1为线性调制区,离起点O距离次近的区域2包含线性调制区和过调制区,离起点O距离最远的区域3为过调制区,根据公式10中的空间位置角θref确定参考电压矢量uref所在扇区si
Figure BDA0002345457960000071
上式中,函数floor()为向下取整函数。
进一步的,所述步骤5中具体包括:
以参考电压矢量uref在第2扇区为例,其它扇区以此类推,当参考电压矢量uref在第2扇区时,与其相邻的两个非零基本矢量为u1、u2,零矢量为u0、u7;基本矢量u1、u2、u0、u7合成的矢量起点为O,终点落在OA、AB、OB上,在第2扇区,AB左下角的区域1为线性调制区,AB右上角的区域3为过调制区,区域2包含线性调制区和过调制区;备选矢量最佳作用时间为t0~t2;其中,t1表示u1、u3、u5的作用时间,t2表示u2、u4、u6的作用时间,t0表示零矢量的作用时间,u0、u7的作用时间分别等于0.5t0,t0、t1、t2中至少有一个等于0;
(a)当参考电压矢量uref位于区域1时,uref需满足uαref≤|u|/2;此时线性调制区uOD到uref的距离最近,最佳矢量为uOD;在静止坐标系中,联立OA和CD的表达式,如公式(12)所示:
Figure BDA0002345457960000081
可得D点的横坐标为:
Figure BDA0002345457960000082
其中,|u|为备选电压矢量中非零电压矢量的幅值,也是非零电压u2的幅值;
OA是u2的幅值,长度为|u|,则
Figure BDA0002345457960000083
OD=2xD
则u2的作用时间为:
Figure BDA0002345457960000084
u1的作用时间为t1=0;零矢量的作用时间为t0=Ts-t1-t2
(b)当参考电压矢量uref位于区域2时,uref需满足uαref>0.5|u|且
Figure BDA0002345457960000085
无论uref在线性调制区还是过调制区,均有线性调制区uOF到uref的距离最近,矢量uOF的一端点为O,另一端点在AB上,则此时最佳矢量为uOF;在静止坐标系中,联立AB和EF的表达式,如公式(15)所示:
Figure BDA0002345457960000091
可得F点的坐标为:
Figure BDA0002345457960000092
由t1u1+t2u2=uOF,得
Figure BDA0002345457960000093
得u2的作用时间为:
Figure BDA0002345457960000094
u1的作用时间为t1=Ts-t2,u0的作用时间为t0=0;
(c)当参考电压矢量uref位于区域3时,uref需满足
Figure BDA0002345457960000095
由于此时线性调制区u2到uref的距离最近,则此时最佳矢量为u2,此时u2的作用时间为t2=Ts,u1的作用时间为t1=0,u0的作用时间为t0=Ts-t1-t2
进一步的,所述步骤6中具体包括:将基本矢量的作用时间输入到脉冲发生器,通过控制逆变器桥臂的通断施加给电机对应的定子电压,以完成对永磁同步电机的控制。
在另外一个方面,本发明还公开了一种永磁同步电机预测电流控制装置,包括:至少一个处理器;以及与所述处理器通信连接的至少一个存储器,其中:
所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令能够执行如上述任一项所述的电流控制方法。
在另外一个方面,本发明还公开了一种非暂态计算机可读存储介质,所述非暂态计算机可读存储介质存储计算机指令,所述计算机指令使所述计算机执行如上述任一项所述的电流控制方法。
(三)有益效果
相对于现有技术,本发明的预测电流控制方法中具体给出了优化了12扇区和区域1-3的划分方法,在考虑了过调制情况下选取出了最佳矢量,最后给出了各个矢量的最佳作用时间t0~t2的精确简化表达式,该方法改进了已有过调制优化方法且减少了三角函数的使用,大大减少了计算量,并同时进一步减小了永磁同步电机交直轴电流纹波,减小转矩脉动和定子电流总谐波畸变率,提高了系统的稳态性能。
附图说明
图1是本发明永磁同步电机预测电流控制系统的控制结构图。
图2是本发明的预测电流控制方法的流程图。
图3是本发明的预测电流控制方法的12扇区划分图。
图4是本发明的预测电流控制方法的12扇区细分图。
图5是参考矢量uref位于区域1时的位置图。
图6是参考矢量uref位于区域2时的位置图
图7是参考矢量uref位于区域3时的位置图。
图8是采用本发明的预测电流控制方法下电机转速为1500rpm突加15Nm负载时的三相定子电流、交直轴电流和转矩波形图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
图1为本发明的系统结构图,主要包括直流电源、三相两电平电压源型逆变器、控制器和永磁同步电机等。其中,永磁同步电机作为受控单元,三相电压源型逆变器作为功率模块以驱动受控单元,控制器从受控单元采集电机系统的机械量和电气量,再采用本发明的控制方法来选择最优的控制矢量通过调制单元输入到逆变器,实现电机的电流控制。
本发明方法的控制流程图如图2所示。其永磁同步电机预测电流控制方法包括如下步骤:
步骤1:对当前时刻永磁同步电机的直流母线电压udc(k)、三相电流ia(k),ib(k),ic(k)、转子旋转电角速度ωe(k)与转子电角度θe(k)进行采样;
步骤2:计算dq轴的参考电压矢量udref和uqref
上述步骤2中具体包括以下步骤:
步骤2.1:根据采集的三相电流以及转子电角度计算两相旋转坐标系下的dq轴电流id(k)和iq(k):
Figure BDA0002345457960000111
步骤2.2:永磁同步电机电压方程为:
Figure BDA0002345457960000112
Figure BDA0002345457960000121
其中,Rs、ψf分别表示电机定子电阻、转子永磁体磁链;Ld,Ld分别表示旋转坐标系下d轴、q轴上的电感分量;
采用一阶欧拉离散方程:
Figure BDA0002345457960000122
其中Ts为采样周期,由公式(3)-(4)得k+1时刻电流预测值为:
Figure BDA0002345457960000123
Figure BDA0002345457960000124
步骤2.3:通过以下目标函数选取最优电压矢量:
Figure BDA0002345457960000125
其中,
Figure BDA0002345457960000126
分别为d轴和q轴参考电流,电流预测控制的目标是dq轴电流预测值跟随参考电流,理想情况下有
Figure BDA0002345457960000127
Figure BDA0002345457960000128
代入到公式(2)(3),可计算同步旋转坐标系上d轴、q轴上的参考电压矢量udref,uqref
Figure BDA0002345457960000129
步骤3:将udref和uqref转换为静止坐标系;根据两相旋转到两相静止的反Park变换,计算出静止坐标系的参考电压矢量uαref,uβref,uαref和uβref分别具体为静止坐标系下参考电压矢量对应的定子电压实部和虚部分量,两者共同合成uref,其具体满足下式9:
Figure BDA0002345457960000131
步骤4:计算静止坐标系的参考电压矢量uαref和uβref所在扇区si
所述步骤4包括以下步骤:
步骤4.1:计算参考电压矢量uref的空间位置角θref
Figure BDA0002345457960000132
公式(10)所求空间位置角θref的取值范围为[-π,π],将其转换到[0,2π]:如果uβref大于等于0,θref不变;如果uβref小于0,θref变为θref+2π。
步骤4.2:对于三相永磁同步电机,对应于两电平三相电压源逆变器,共有8个基本电压矢量,包括2个零矢量与6个非零矢量;非零矢量在空间上从0度开始,逆时针旋转,间隔为π/3,至2π共有6个基本电压矢量均匀分布,依次定义为ui,其中i=1,2,3,…,6,零矢量定义为u0和u7
将8个基本电压矢量所在平面划分为12个扇区,即扇区si,下标i为1到12之间的自然数,每个扇区si的圆心角为30°,如图3所示;在每个扇区非零矢量的一端,作两条与非零矢量夹角为30°的线,将每个扇区划分为3个区域,每个扇区包含区域1-3,离起点O距离最近的区域1为线性调制区,离起点O距离次近的区域2包含线性调制区和过调制区,离起点O距离最远的区域3为过调制区,如图4和图6-7所示。可根据公式10中的空间位置角θref确定参考电压矢量uref所在扇区si
Figure BDA0002345457960000141
上式11中,函数floor()为向下取整函数。
步骤5:判断uαβref在扇区si对应区域并选择最佳矢量,并根据最佳矢量计算最佳作用时间t0~t2
步骤5中的具体包括:
考虑过调制区和线性调制区,将扇区分成多个区域,选择线性调制区内距离参考电压矢量最近的矢量作为最佳矢量,在考虑区域所属类型的情况下(即线性调制区和过调制区),并根据对应区域内最佳矢量与备选矢量的关系计算备选矢量中各个基本电压矢量的最佳作用时间t0~t2;其中,t1表示u1、u3、u5的作用时间,t2表示u2、u4、u6的作用时间,t0表示零矢量的作用时间,u0、u7的作用时间分别等于0.5t0,t0、t1、t2中至少有一个等于0。
以参考电压矢量uref在第2扇区为例,说明如何在划分特定区域的情况下选取最佳矢量和计算简化后的作用时间。
当参考电压矢量uref在第2扇区时,与其相邻的两个非零基本矢量为u1、u2,零矢量为u0、u7;基本矢量u1、u2、u0、u7合成的矢量起点为O,终点落在OA、AB、OB上,在第2扇区,AB左下角的区域1为线性调制区,AB右上角的区域3为过调制区,区域2包含线性调制区和过调制区,设t1表示u1、u3、u5的作用时间,t2表示u2、u4、u6的作用时间,t0表示零矢量的作用时间,图5-7具体标注了扇区2的区域1-3的划分,在每个扇区非零矢量的一端,作两条与非零矢量夹角为30°的线,分成3个区域。
(a)当参考电压矢量uref位于区域1时,如图5所示,uref需满足uαref≤0.5|u|;此时线性调制区uOD到uref的距离最近,最佳矢量为uOD。在静止坐标系中,联立OA和CD的表达式,如公式(12)所示:
Figure BDA0002345457960000151
可得D点的横坐标为:
Figure BDA0002345457960000152
其中,|u|为备选电压矢量中非零电压矢量的幅值,也是非零电压u2的幅值;
OA是u2的幅值,长度为|u|,则
Figure BDA0002345457960000153
OD=2xD
则u2的作用时间为:
Figure BDA0002345457960000154
u1的作用时间为t1=0,零矢量的作用时间为t0=Ts-t1-t2
(b)当参考电压矢量uref位于区域2时,如图6所示,uref需满足uαref>0.5|u|且
Figure BDA0002345457960000155
无论uref在线性调制区还是过调制区,均有线性调制区uOF到uref的距离最近,矢量uOF的一端点为O,另一端点在AB上,则此时最佳矢量为uOF,在静止坐标系中,联立AB和EF的表达式,如公式(15)所示:
Figure BDA0002345457960000161
可得F点的坐标为:
Figure BDA0002345457960000162
由t1u1+t2u2=uOF,得
Figure BDA0002345457960000163
得u2的作用时间为:
Figure BDA0002345457960000164
u1的作用时间为t1=Ts-t2,u0的作用时间为t0=0。
(c)当参考电压矢量uref位于区域3时,如图7所示,uref需满足
Figure BDA0002345457960000165
由于此时线性调制区u2到uref的距离最近,则此时最佳矢量为u2,此时u2的作用时间为t2=Ts,u1的作用时间为t1=0,u0的作用时间为t0=Ts-t1-t2
同理,可计算参考矢量uref位于其他扇区时基本矢量的作用时间,并求出在考虑了过调制优化情况下所有的uref对应的12组最佳作用时间t0~t2
值得一提的是,由于每个扇区的位置不一样,则最佳作用时间t0~t2的公式和对应uref需满足的条件也会有所不一样,但是其简化时间t0~t2的计算方式都与上述第2扇区的计算原理相同,此处限于篇幅,故不再赘述其它11个扇区的计算方式,因为这些对本领域技术人员来说都是不用再付出创造性劳动就能够得到的结果。其12个扇区的作用时间表参见表1。
表1 12个扇区的作用时间表
Figure BDA0002345457960000171
Figure BDA0002345457960000181
Figure BDA0002345457960000191
Figure BDA0002345457960000201
Figure BDA0002345457960000211
上述表1中,区域1和区域3内的零矢量作用时间为t0=Ts-t1-t2
步骤6:应用最佳作用时间t0~t2;将基本矢量的作用时间输入到脉冲发生器,通过控制逆变器桥臂的通断施加给电机对应的定子电压,以完成对永磁同步电机的控制。
为验证本发明控制方法的优势,与id和iq无差拍且无过调制优化的二矢量预测电流控制策略进行对比,测试工况为电机转速为1500r/min在t=0.5s时突加15Nm负载。本发明控制方法下的系统工作波形如图8所示,图8中的图a、b、c分别表示三相定子电流、交直轴电流和转矩的波形。采用标准差计算得到两种控制策略直交轴电流脉动、转矩脉动及输出电流总谐波畸变率(THD)如表2所示。
表2电机转速为1500r/min突加15Nm负载时两种策略的仿真结果对比
Figure BDA0002345457960000212
此外,值得一提的是,对于专利CN201810787417.5的方法,经过发明人的实际仿真,发现其谐波和转矩的波动都很大,且其方法存在较大误差,在仿真方面,该专利的效果比不上传统的单矢量预测控制,更不用说传统的二矢量预测控制和本发明的方法。
需要说明的是,上述的电流控制方法可以转换为程序指令,既可以使用包括处理器和存储器的电流控制装置来运行实现,也可以通过非暂态计算机可读存储介质中存储的计算机指令来实现。
最后,本发明的方法仅为较佳的实施方案,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种永磁同步电机预测电流控制方法,其特征在于,该电流控制方法包括如下步骤:
步骤1:对当前时刻永磁同步电机的直流母线电压udc(k)、三相电流ia(k),ib(k),ic(k)、转子旋转电角速度ωe(k)与转子电角度θe(k)进行采样;
步骤2:计算dq轴的参考电压矢量udref和uqref
步骤3:将udref和uqref转换为静止坐标系;根据两相旋转到两相静止的反Park变换,计算出静止坐标系的参考电压矢量uαref,uβref,uαref和uβref分别为静止坐标系下参考电压矢量对应的定子电压实部和虚部分量,两者共同合成参考电压矢量uref
步骤4:计算静止坐标系的参考电压矢量uαref和uβref所在扇区si
步骤5:判断uref在扇区si对应区域并选择最佳矢量,根据已推导的最佳矢量作用时间的表达式直接计算作用时间;该步骤5中具体包括:考虑过调制区和线性调制区,将扇区分成多个区域,选择与参考电压矢量uref相邻的两个非零矢量和零矢量作为备选矢量,选择线性调制区内距离参考电压矢量最近的矢量作为最佳矢量,并根据对应区域内最佳矢量与备选矢量的关系计算备选矢量中各个基本电压矢量的最佳作用时间。
步骤6:应用基本电压矢量的最佳作用时间,以完成对永磁同步电机的控制。
2.根据权利要求1所述的永磁同步电机预测电流控制方法,其特征在于,所述步骤2中具体包括以下步骤:
步骤2.1:根据采集的三相电流以及转子电角度计算两相旋转坐标系下的dq轴电流id(k)和iq(k):
Figure FDA0002345457950000021
步骤2.2:永磁同步电机电压方程为:
Figure FDA0002345457950000022
Figure FDA0002345457950000023
其中,Rs、ψf分别表示电机定子电阻、转子永磁体磁链;Ld,Lq分别表示旋转坐标系下d轴、q轴上的电感分量;
采用一阶欧拉离散方程:
Figure FDA0002345457950000024
其中Ts为采样周期,由公式(3)-(4)得k+1时刻电流预测值为:
Figure FDA0002345457950000025
Figure FDA0002345457950000026
步骤2.3:通过以下目标函数选取最优电压矢量:
Figure FDA0002345457950000027
其中,
Figure FDA0002345457950000028
分别为d轴和q轴参考电流,电流预测控制的目标是dq轴电流预测值跟随参考电流,理想情况下有
Figure FDA0002345457950000029
Figure FDA00023454579500000210
代入到公式(2)~(3),可计算同步旋转坐标系上d轴、q轴上的参考电压矢量udref,uqref
Figure FDA0002345457950000031
3.根据权利要求1所述的永磁同步电机预测电流控制方法,其特征在于,所述步骤3中具体满足下式:
Figure FDA0002345457950000032
4.根据权利要求1所述的永磁同步电机预测电流控制方法,其特征在于,所述步骤4包括以下步骤:
步骤4.1:计算参考电压矢量uref的空间位置角θref
Figure FDA0002345457950000033
公式(10)所求空间位置角θref的取值范围为[-π,π],将其转换到[0,2π]:如果uβref大于等于0,θref不变;如果uβref小于0,θref变为θref+2π;
步骤4.2:对于三相永磁同步电机,对应于两电平三相电压源逆变器,共有8个基本电压矢量,包括2个零矢量与6个非零矢量;非零矢量在空间上从0度开始,逆时针旋转,间隔为π/3,至2π共有6个基本电压矢量均匀分布,依次定义为ui,其中i=1,2,3,…,6,零矢量定义为u0和u7;将8个基本电压矢量所在平面划分为12个扇区,即扇区si,下标i为1到12之间的自然数,每个扇区si的圆心角为30°;在每个扇区非零矢量的一端,作两条与非零矢量夹角为30°的线,将每个扇区划分为3个区域,每个扇区包含区域1-3,离起点O距离最近的区域1为线性调制区,离起点O距离次近的区域2包含线性调制区和过调制区,离起点O距离最远的区域3为过调制区,根据公式10中的空间位置角θref确定参考电压矢量uref所在扇区si
Figure FDA0002345457950000041
上式中,函数floor()为向下取整函数。
5.根据权利要求4所述的永磁同步电机预测电流控制方法,其特征在于,所述步骤5中具体包括:
以参考电压矢量uref在第2扇区为例,其它扇区以此类推,当参考电压矢量uref在第2扇区时,与其相邻的两个非零基本矢量为u1、u2,零矢量为u0、u7;基本矢量u1、u2、u0、u7合成的矢量起点为O,终点落在OA、AB、OB上,在第2扇区,AB左下角的区域1为线性调制区,AB右上角的区域3为过调制区,区域2包含线性调制区和过调制区;备选矢量最佳作用时间为t0~t2;其中,t1表示u1、u3、u5的作用时间,t2表示u2、u4、u6的作用时间,t0表示零矢量的作用时间,u0、u7的作用时间分别等于0.5t0,t0、t1、t2中至少有一个等于0;
(a)当参考电压矢量uref位于区域1时,uref需满足uαref≤|u|/2;此时线性调制区uOD到uref的距离最近,最佳矢量为uOD;在静止坐标系中,联立OA和CD的表达式,如公式(12)所示:
Figure FDA0002345457950000042
可得D点的横坐标为:
Figure FDA0002345457950000051
其中,|u|为备选电压矢量中非零电压矢量的幅值,也是非零电压u2的幅值;
OA是u2的幅值,长度为|u|,则
Figure FDA0002345457950000052
OD=2xD
则u2的作用时间为:
Figure FDA0002345457950000053
u1的作用时间为t1=0;零矢量的作用时间为t0=Ts-t1-t2
(b)当参考电压矢量uref位于区域2时,uref需满足uαref>0.5|u|且
Figure FDA0002345457950000054
无论uref在线性调制区还是过调制区,均有线性调制区uOF到uref的距离最近,矢量uOF的一端点为O,另一端点在AB上,则此时最佳矢量为uOF;在静止坐标系中,联立AB和EF的表达式,如公式(15)所示:
Figure FDA0002345457950000055
可得F点的坐标为:
Figure FDA0002345457950000056
由t1u1+t2u2=uOF,得
Figure FDA0002345457950000057
得u2的作用时间为:
Figure FDA0002345457950000061
u1的作用时间为t1=Ts-t2,u0的作用时间为t0=0。
(c)当参考电压矢量uref位于区域3时,uref需满足
Figure FDA0002345457950000062
由于此时线性调制区u2到uref的距离最近,则此时最佳矢量为u2,此时u2的作用时间为t2=Ts,u1的作用时间为t1=0,u0的作用时间为t0=Ts-t1-t2
6.根据权利要求1所述的永磁同步电机预测电流控制方法,其特征在于,所述步骤6中具体包括:将基本矢量的作用时间输入到脉冲发生器,通过控制逆变器桥臂的通断施加给电机对应的定子电压,以完成对永磁同步电机的控制。
7.一种永磁同步电机预测电流控制装置,其特征在于,包括:
至少一个处理器;以及与所述处理器通信连接的至少一个存储器,其中:
所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令能够执行如权利要求1至6任一项所述的预测电流控制方法。
8.一种非暂态计算机可读存储介质,其特征在于,所述非暂态计算机可读存储介质存储计算机指令,所述计算机指令使所述计算机执行如权利要求1至6任一项所述的预测电流控制方法。
CN201911392796.9A 2019-12-30 2019-12-30 一种永磁同步电机预测电流控制方法和装置 Active CN111030542B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911392796.9A CN111030542B (zh) 2019-12-30 2019-12-30 一种永磁同步电机预测电流控制方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911392796.9A CN111030542B (zh) 2019-12-30 2019-12-30 一种永磁同步电机预测电流控制方法和装置

Publications (2)

Publication Number Publication Date
CN111030542A true CN111030542A (zh) 2020-04-17
CN111030542B CN111030542B (zh) 2021-03-19

Family

ID=70195886

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911392796.9A Active CN111030542B (zh) 2019-12-30 2019-12-30 一种永磁同步电机预测电流控制方法和装置

Country Status (1)

Country Link
CN (1) CN111030542B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111464095A (zh) * 2020-05-14 2020-07-28 长沙丹芬瑞电气技术有限公司 一种永磁同步电机控制方法及系统
CN111464100A (zh) * 2020-05-14 2020-07-28 长沙丹芬瑞电气技术有限公司 一种基于模型预测控制的永磁同步电机控制方法及系统
CN111654064A (zh) * 2020-06-23 2020-09-11 广东电网有限责任公司电力科学研究院 一种虚拟同步发电机的控制方法和相关装置
CN111682817A (zh) * 2020-06-29 2020-09-18 潍柴动力股份有限公司 内嵌式永磁同步电机的控制方法及装置
CN111830435A (zh) * 2020-07-27 2020-10-27 北京航空航天大学 一种六相永磁容错电机系统功率管开路故障诊断方法
CN112054736A (zh) * 2020-09-11 2020-12-08 南通大学 一种优化分区调制的永磁同步电机模型预测电流过调制控制方法
CN113179065A (zh) * 2021-05-12 2021-07-27 西南交通大学 一种永磁同步电机模型预测脉冲序列控制方法
CN118694233A (zh) * 2024-08-26 2024-09-24 江西理工大学 一种电力牵引系统及其控制方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107134966A (zh) * 2017-04-18 2017-09-05 天津大学 一种基于有限状态预测的开绕组永磁电机电流控制方法
CN108880374A (zh) * 2018-07-02 2018-11-23 华中科技大学 直线感应电机含电流约束双矢量模型预测控制方法
CN109039189A (zh) * 2018-07-17 2018-12-18 东南大学 基于几何法的永磁同步电机两矢量预测控制系统及方法
CN109586637A (zh) * 2018-12-04 2019-04-05 南京工程学院 一种基于改进双矢量mpc的pmslm电流控制方法
CN109802615A (zh) * 2019-01-10 2019-05-24 北方工业大学 双矢量绕组开路永磁同步电机电流预测控制方法和设备

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107134966A (zh) * 2017-04-18 2017-09-05 天津大学 一种基于有限状态预测的开绕组永磁电机电流控制方法
CN108880374A (zh) * 2018-07-02 2018-11-23 华中科技大学 直线感应电机含电流约束双矢量模型预测控制方法
CN109039189A (zh) * 2018-07-17 2018-12-18 东南大学 基于几何法的永磁同步电机两矢量预测控制系统及方法
CN109586637A (zh) * 2018-12-04 2019-04-05 南京工程学院 一种基于改进双矢量mpc的pmslm电流控制方法
CN109802615A (zh) * 2019-01-10 2019-05-24 北方工业大学 双矢量绕组开路永磁同步电机电流预测控制方法和设备

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
孙晓: "基于电流预测的永磁同步电机矢量控制算法研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 *
孙翀: "《开绕组永磁同步电机驱动系统优化模型预测控制技术研究》", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111464095A (zh) * 2020-05-14 2020-07-28 长沙丹芬瑞电气技术有限公司 一种永磁同步电机控制方法及系统
CN111464100A (zh) * 2020-05-14 2020-07-28 长沙丹芬瑞电气技术有限公司 一种基于模型预测控制的永磁同步电机控制方法及系统
CN111464095B (zh) * 2020-05-14 2022-03-18 长沙丹芬瑞电气技术有限公司 一种永磁同步电机控制方法及系统
CN111654064A (zh) * 2020-06-23 2020-09-11 广东电网有限责任公司电力科学研究院 一种虚拟同步发电机的控制方法和相关装置
CN111682817A (zh) * 2020-06-29 2020-09-18 潍柴动力股份有限公司 内嵌式永磁同步电机的控制方法及装置
CN111830435A (zh) * 2020-07-27 2020-10-27 北京航空航天大学 一种六相永磁容错电机系统功率管开路故障诊断方法
CN111830435B (zh) * 2020-07-27 2021-07-13 北京航空航天大学 一种六相永磁容错电机系统功率管开路故障诊断方法
CN112054736A (zh) * 2020-09-11 2020-12-08 南通大学 一种优化分区调制的永磁同步电机模型预测电流过调制控制方法
CN113179065A (zh) * 2021-05-12 2021-07-27 西南交通大学 一种永磁同步电机模型预测脉冲序列控制方法
CN113179065B (zh) * 2021-05-12 2022-08-19 西南交通大学 一种永磁同步电机模型预测脉冲序列控制方法
CN118694233A (zh) * 2024-08-26 2024-09-24 江西理工大学 一种电力牵引系统及其控制方法

Also Published As

Publication number Publication date
CN111030542B (zh) 2021-03-19

Similar Documents

Publication Publication Date Title
CN111030542B (zh) 一种永磁同步电机预测电流控制方法和装置
US11356045B2 (en) Double virtual voltage vectors predictive torque control method without weighting factor for five-phase permanent magnet synchronous motor
CN110995076A (zh) 一种永磁同步电机模型预测电流控制方法
CN106936356B (zh) 矢量筛选和占空比结合的电机模型预测控制系统及方法
CN106788027B (zh) 一种基于最优二矢量组合的模型预测控制方法及系统
WO2021031685A1 (zh) 一种大功率直驱永磁同步电机控制调制方法
CN109347387A (zh) 基于模型预测的电机控制方法及控制装置
CN112217437A (zh) 一种永磁同步电机三矢量模型预测电流控制电路和方法
CN111654225B (zh) 一种开绕组永磁游标直线电机的模型预测电流控制方法
Xu et al. Three-vector-based model predictive current control with zero-sequence current suppression for open-winding LPMVM drives
Zhang et al. Multi-stage series model predictive control for PMSM drives
Saeed et al. Generalized deadbeat solution for model predictive control of five-phase PMSM drives
CN115021636A (zh) 一种基于无差拍的五相永磁同步电机模型预测容错控制方法
CN115833671A (zh) 一种双三相电机高精度模型预测电流控制系统及控制方法
Chen et al. Band-based multi-step predictive torque control strategy for PMSM drives
Li et al. Model predictive current control algorithm based on joint modulation strategy for low-inductance PMSM
CN111900907B (zh) 一种基于切换点优化的永磁电机模型预测磁链控制方法
JP2010172133A (ja) Ipmモータ制御装置及び制御方法
Riccio et al. Modulated model-predictive integral control applied to a synchronous reluctance motor drive
CN113992095B (zh) 一种低复杂度的双三相永磁同步发电机pmsg模型预测直接功率控制方法
CN115276501B (zh) 一种永磁同步电机无差拍预测电流控制方法
Ma et al. Dichotomy solution based model predictive control for permanent magnet linear synchronous motors
CN112467776A (zh) 电流源型变流器系统、控制方法及空间矢量调制方法
Hadla Predictive load angle and stator flux control of SynRM drives for the full speed range
CN109831136B (zh) 一种三相交流电机的驱动方法、系统及装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230810

Address after: Room 705, Guishui Building, No. 469, Section 1, Wanjiali Middle Road, Gaoqiao Street, Yuhua Block, Changsha City, Hunan Province, 410199

Patentee after: Hunan Tianzheng Youxuan Enterprise Planning Co.,Ltd.

Address before: Yuelu District City, Hunan province 410082 Changsha Lushan Road No. 2

Patentee before: HUNAN University

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20231013

Address after: Room 401, 4th Floor, Building G1-1, Lushan Science and Technology Innovation Park, No. 966 Lushan South Road, Science Village Community, Yuelu Street, Yuelu District, Changsha City, Hunan Province, 410000

Patentee after: Hunan Lanhai Electrical Engineering Co.,Ltd.

Address before: Room 705, Guishui Building, No. 469, Section 1, Wanjiali Middle Road, Gaoqiao Street, Yuhua Block, Changsha City, Hunan Province, 410199

Patentee before: Hunan Tianzheng Youxuan Enterprise Planning Co.,Ltd.