CN107612446B - 一种内置式永磁同步电机模型预测转矩控制方法 - Google Patents

一种内置式永磁同步电机模型预测转矩控制方法 Download PDF

Info

Publication number
CN107612446B
CN107612446B CN201711075880.9A CN201711075880A CN107612446B CN 107612446 B CN107612446 B CN 107612446B CN 201711075880 A CN201711075880 A CN 201711075880A CN 107612446 B CN107612446 B CN 107612446B
Authority
CN
China
Prior art keywords
synchronous motor
torque
permanent magnet
motor
coordinate system
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.)
Active
Application number
CN201711075880.9A
Other languages
English (en)
Other versions
CN107612446A (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.)
Tianjin Polytechnic University
Original Assignee
Tianjin Polytechnic 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 Tianjin Polytechnic University filed Critical Tianjin Polytechnic University
Priority to CN201711075880.9A priority Critical patent/CN107612446B/zh
Publication of CN107612446A publication Critical patent/CN107612446A/zh
Application granted granted Critical
Publication of CN107612446B publication Critical patent/CN107612446B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Ac Motors In General (AREA)
  • External Artificial Organs (AREA)

Abstract

一种内置式永磁同步电机模型预测转矩控制方法:通过系统采样和电机模型,补偿数字系统延时问题,并计算电机参考转矩Te *,使其经最大转矩电流比计算得到参考电流d、q轴分量;利用电机模型和备选空间电压矢量预测相应的永磁体转矩和磁阻转矩分量以及定子磁链d、q轴分量;对电机轻重载判断并求出最优电压矢量,即当电机在重载情况下使用永磁体转矩和磁阻转矩评价函数,在轻载时使用磁链矢量评价函数;进而将所求电压矢量作用于电机。本发明通过对内置式永磁同步电机模型预测转矩控制方法在最大转矩电流比控制下转矩的产生进行分析,设计了不需要系数整定的评价函数,该方法能够使电机在重载时具有更好的转矩控制效果。

Description

一种内置式永磁同步电机模型预测转矩控制方法
技术领域
本发明涉及一种电机模型预测转矩控制方法。特别是涉及一种内置式永磁同步电机模型预测转矩控制方法。
背景技术
随着电力电子技术的快速发展,人们对电能的使用更加广泛,电机作为能量转换的桥梁受到了前所未有的关注。其中,内置式永磁同步电机以其可靠性好,功率密度高等优点在工业制造、轨道交通、航空航天等领域广泛应用。传统的内置式永磁同步电机控制主要以获得电机最大的转矩-电流比为目标,控制策略主要包括空间矢量控制和直接转矩控制,其中,直接转矩控制具有结构简单,动态性能良好和鲁棒性强等优点。但由于直接转矩控制策略采用滞环控制器,因此电机的输出转矩波动较大。
模型预测控制策略具有动态响应快、目标函数配置灵活、易于处理约束优化问题等优点,特别是近年来,被广泛应用于电气驱动领域。模型预测转矩控制使用精确的内部模型取代转矩和磁链滞环控制器,在当前控制周期内对所有开关作用所产生的电机输出状态进行预测,通过评价函数使逆变器输出控制最优的电压矢量,从而使电机输出的转矩效果最佳。因此,控制器设计过程中,评价函数的设计对控制性能优劣的影响至关重要。
传统模型预测转矩控制的评价函数由转矩和磁链两个分量组成,但由于两个分量的单位不同,需要设计相应的系数来实现不同单位分量的同时控制。而该系数的选取一般缺乏通用的理论原则,需要通过大量的仿真和实验调试所取得的数据来确定,工作量较大。基于模型预测转矩控制的思想,将评价函数改为磁链矢量的方法改善了这个问题。但由于电机调速控制的直接因素是转矩,而该方法却并没能将转矩设置为直接控制量,所以导致电机的动态响应和转矩的控制效果变差,仍需要进一步改善。
传统模型预测转矩控制的评价函数由转矩和磁链两个分量组成,由于两个分量单位不同,一般需要设计相应的系数实现不同单位分量的同时控制,而该系数的选取较为繁琐。
发明内容
本发明所要解决的技术问题是,提供一种能够使电机在重载时具有更好的转矩控制效果的内置式永磁同步电机模型预测转矩控制方法。
本发明所采用的技术方案是:一种内置式永磁同步电机模型预测转矩控制方法,是用于内置式永磁同步电机控制系统的控制方法,包括如下步骤:
1)采集内置式永磁同步电机第k个控制周期开始时刻的电压电流信号以及转速信号,并进行坐标变换;
2)通过延时补偿模型,预测得到内置式永磁同步电机第(k+1)个控制周期开始时刻的旋转坐标系d-q下的电流分量;
3)将转速信号中的同步电机转子机械角速度ω(k)与给定转速ω*之间的差值,经PI控制器得到同步电机第k个控制周期的电磁转矩给定值Te *;根据电磁转矩给定值Te *,采用最大转矩电流比计算方法,得到旋转坐标系d-q下电流分量的给定值id *和iq *,并据电流分量的给定值id *和iq *求出同步电机定子磁链ψs的给定值和同步电机输出转矩的给定值;
4)将步骤2)预测得到的电流分量,以及内置式永磁同步电机控制系统中的逆变器中的8个电压矢量V0、V1、V2、……、V7,带入到电机预测模型中,预测出不同电压矢量Vn作用下的同步电机磁链和转矩,n=0,1,2……7,在预测过程中,由于电压矢量V0和电压矢量V7的作用效果相同,不考虑电压矢量V0的相关计算;
5)将预测出不同电压矢量Vn作用下的电机磁链和转矩,带入到评价函数中,从评价函数的计算结果中找到最小值,所述最小值所对应的逆变器电压矢量,即为最优电压矢量Vopt,并通过逆变器输出所述最优电压矢量。
步骤1)包括:
通过主控板上微处理器内部的模数转换接口采集第k个控制周期开始时刻的逆变器直流母线电压udc(k)和三相定子电流ia(k)、ib(k)、ic(k);
利用编码器得到第k个控制周期开始时刻的同步电机转子机械角速度ω(k)和电角速度ωe(k)以及转子位置电角度θ(k);
将采集到的三相定子电流ia(k)、ib(k)、ic(k),通过三相/两相静止坐标变换得到两相静止坐标系α-β下电流分量iα(k)和iβ(k),变换公式如下:
Figure BDA0001457022540000021
再通过两相旋转坐标变换得到旋转坐标系d-q下的电流分量id(k)和iq(k),变换公式如下:
Figure BDA0001457022540000022
式中ia(k)、ib(k)、ic(k)分别为同步电机三相定子电流,iα(k)、iβ(k)分别为同步电机两相静止坐标系α-β下的电流分量,id(k)、iq(k)分别为同步电机旋转坐标系d-q下的电流分量,θ(k)为同步电机转子位置电角度。
步骤2)包括:利用坐标变换后的电压电流信号以及转速信号包括电流分量id(k)、iq(k)和电角速度ωe(k),以及第k个控制周期施加的电压矢量V(k)在旋转坐标系d-q下的电压分量Vd(k)、Vq(k),通过延时补偿模型,预测得到第(k+1)个控制周期开始时刻的旋转坐标系d-q下电流分量
Figure BDA0001457022540000023
所述的延时补偿模型公式如下:
Figure BDA0001457022540000024
式中Ld,Lq分别为同步电机的直轴和交轴电感,Rs为定子电阻,Ts为控制周期的时长,ψf为转子永磁体磁链,ωe(k)为同步电机电角速度,Vd(k)、Vq(k)分别为第k个控制周期施加的电压矢量V(k)在旋转坐标系d-q下的电压分量,id(k)、iq(k)分别为同步电机旋转坐标系d-q下的电流分量,
Figure BDA0001457022540000031
分别为预测得到的第(k+1)个控制周期开始时刻的旋转坐标系d-q下电流分量,△Vd(k)、△Vq(k)分别由同步电机旋转坐标系d-q下的电流分量id(k)、iq(k)同上个周期预测得到的电流分量
Figure BDA0001457022540000032
做差经过PI控制器得到。
步骤3)所述的同步电机定子磁链ψs的给定值是采用下式计算:
Figure BDA0001457022540000033
式中,Ld,Lq分别为同步电机的直轴和交轴电感,ψd *和ψq *分别为定子磁链在旋转坐标系d-q下d轴的给定分量和q轴的给定分量,ψf为同步电机转子永磁体磁链。
步骤3)所述的同步电机输出转矩的给定值是采用下式计算:
Figure BDA0001457022540000034
式中,TM *为同步电机永磁体输出转矩的给定分量,TR *为同步电机磁阻输出转矩的给定分量,Ld,Lq分别为同步电机的直轴和交轴电感,ψf为同步电机转子永磁体磁链,p为同步电机的极对数。
步骤4)包括:
将预测得到的内置式永磁同步电机第(k+1)个控制周期开始时刻的旋转坐标系d-q下的电流分量
Figure BDA0001457022540000035
Figure BDA0001457022540000036
以及所述逆变器中的8个电压矢量V0、V1、V2、……、V7,带入到电机预测模型,预测出不同电压矢量Vn作用下,第(k+2)控制周期开始时刻同步电机定子磁链ψs在旋转坐标系d-q下的分量
Figure BDA0001457022540000037
Figure BDA0001457022540000038
及永磁体转矩
Figure BDA0001457022540000039
和磁阻转矩
Figure BDA00014570225400000310
电机预测模型公式如下,在预测过程中,由于电压矢量V0和电压矢量V7的作用效果相同,不考虑电压矢量V0的相关计算:
Figure BDA00014570225400000311
Figure BDA00014570225400000312
Figure BDA00014570225400000313
式中Ld,Lq分别为同步电机的直轴和交轴电感,Rs为定子电阻,Ts为控制周期的时长,p为同步电机的极对数,
Figure BDA0001457022540000041
分别为预测得到第(k+1)个控制周期开始时刻的旋转坐标系d-q下电流分量,n为和电压矢量Vn相关变量的的序号、Vdn、Vqn分别为逆变器电压矢量Vn在旋转坐标系d-q下的电压分量,
Figure BDA0001457022540000042
分别为电压矢量Vn作用后第(k+2)控制周期开始时刻同步电机定子磁链ψs在旋转坐标系d-q下的分量,
Figure BDA0001457022540000043
Figure BDA0001457022540000044
分别为电压矢量Vn作用后第(k+2)控制周期开始时刻同步电机永磁体输出转矩和同步电机磁阻输出转矩,ψf为同步电机转子永磁体磁链,由于电机的机械时间常数远大于电气时间常数,认定两个相邻控制周期中转速是不变的即式中ωe(k+1)=ωe(k)。
步骤5)包括:
将步骤4)求得的永磁体转矩
Figure BDA0001457022540000045
和磁阻转矩
Figure BDA0001457022540000046
与同步电机永磁体输出转矩的给定分量TM *和同步电机磁阻输出转矩的给定分量TR *的误差建立同步电机在重载情况衡量电压矢量Vn对控制效果的评价函数g(n)公式如下:
Figure BDA0001457022540000047
根据评价函数g(n)的计算结果,找到评价函数g(n)的最小值g(n)min,则与所述最小值g(n)min对应的电压矢量Vn为最优电压矢量Vopt
在同步电机轻载时需要使用磁链矢量评价函数以避免权重系数整定和电流波动,这时,评价函数g(n)是采用下式得到:
Figure BDA0001457022540000048
式中,ψd *和ψq *分别为定子磁链在旋转坐标系d-q下d轴的给定分量和q轴的给定分量,
Figure BDA0001457022540000049
Figure BDA00014570225400000410
为第(k+2)控制周期开始时刻同步电机定子磁链ψs在旋转坐标系d-q下的分量;
为使同步电机在重载情况下与轻载情况下的评价函数g(n)正常切换,设定常数TX为切换参考转矩,当同步电机的电磁转矩给定值Te *的绝对值大于TX时,认定为重载,当同步电机的电磁转矩给定值Te *的绝对值小于等于TX时认定为轻载。
为避免当同步电机的电磁转矩给定值Te *的绝对值接近于TX且上下波动时,导致同步电机在重载情况下与轻载情况下的评价函数g(n)切换频繁,故加入滞环比较器减少切换频率:设定环宽为σ,当同步电机的电磁转矩给定值Te *的绝对值大于设定常数TX与设定环宽σ之和,即同步电机的电磁转矩给定值Te *的绝对值处于滞环之上时,滞环输出δ(k)=1,为重载时的评价函数g(n);当同步电机的电磁转矩给定值Te *的绝对值小于设定常数TX与设定环宽σ之差时,即同步电机的电磁转矩给定值Te *的绝对值处于滞环之下时,滞环输出δ(k)=-1,为轻载时的评价函数g(n);当TX-σ≤|Te *|≤TX+σ,即同步电机的电磁转矩给定值Te *的绝对值处于滞环之内时,滞环输出等于上一个控制周期的输出,即δ(k)=δ(k-1),使用上一个控制周期使用的评价函数g(n);δ的初值δ(0)=-1;
从而,所述评价函数表示为:
Figure BDA0001457022540000051
式中
Figure BDA0001457022540000052
分别为电压矢量Vn作用后第(k+2)控制周期开始时刻同步电机定子磁链ψs在旋转坐标系d-q下的分量,
Figure BDA0001457022540000053
分别为电压矢量Vn作用后第(k+2)控制周期开始时刻同步电机永磁体输出转矩和同步电机磁阻输出转矩,ψd *和ψq *分别为定子磁链在旋转坐标系d-q下d轴的给定分量和q轴的给定分量,TM *为同步电机永磁体输出转矩的给定分量,TR *为同步电机磁阻输出转矩的给定分量,argming(n)为求数列{g(1),g(2),……,g(n)}中最小值所在列数的运算,即从评价函数的计算结果中找出最小值所对应的电压矢量。
本发明的一种内置式永磁同步电机模型预测转矩控制方法,通过对内置式永磁同步电机在最大转矩电流比控制下转矩的产生进行分析,设计了可直接控制转矩的评价函数,即将传统转矩和磁链分量的评价函数改为了永磁体转矩和磁阻转矩的评价函数,且该评价函数不需要系数的整定。同时结合磁链矢量控制的方法改善了所设计转矩评价函数在轻载下的缺点。该方法相比传统方法省去了繁琐的系数整定工作,相比磁链矢量控制的方法能够使电机在重载时具有更好的转矩控制效果。
附图说明
图1是本发明所用的内置式永磁同步电机控制系统;
图中
1:示波器 2:主控制板
3:DSP芯片 4:逆变器
5:电流传感器 6:电压传感器
7:同步电机 8:编码器
图2是发明一种内置式永磁同步电机模型预测转矩控制方法的流程图;
图3是发明一种内置式永磁同步电机模型预测转矩控制方法的框图;
图4是发明一种内置式永磁同步电机模型预测转矩控制方法中延时补偿模块的框图;
图5是发明所用空间电压矢量图。
具体实施方式
下面结合实施例和附图对本发明的一种内置式永磁同步电机模型预测转矩控制方法做出详细说明。
本发明的一种内置式永磁同步电机模型预测转矩控制方法,是用于图1所示的内置式永磁同步电机控制系统的控制方法,如图2、图3、图4所示,包括如下步骤:
1)采集内置式永磁同步电机第k个控制周期开始时刻的电压电流信号以及转速信号,并进行坐标变换;包括:
通过主控板上微处理器(DSP)内部的模数转换(A/D)接口采集第k个控制周期开始时刻的逆变器直流母线电压udc(k)和三相定子电流ia(k)、ib(k)、ic(k);
利用编码器得到第k个控制周期开始时刻的同步电机转子机械角速度ω(k)和电角速度ωe(k)以及转子位置电角度θ(k);
将采集到的三相定子电流ia(k)、ib(k)、ic(k),通过三相/两相静止坐标变换得到两相静止坐标系
Figure BDA0001457022540000061
下电流分量iα(k)和iβ(k),变换公式如下:
Figure BDA0001457022540000062
再通过两相旋转坐标变换得到旋转坐标系d-q下的电流分量id(k)和iq(k),变换公式如下:
Figure BDA0001457022540000063
式中ia(k)、ib(k)、ic(k)分别为同步电机三相定子电流,iα(k)、iβ(k)分别为同步电机两相静止坐标系α-β下的电流分量,id(k)、iq(k)分别为同步电机旋转坐标系d-q下的电流分量,θ(k)为同步电机转子位置电角度。
2)为解决数字控制系统控制滞后的问题,通过延时补偿模型,预测得到内置式永磁同步电机第(k+1)个控制周期开始时刻的旋转坐标系d-q下的电流分量;包括:利用坐标变换后的电压电流信号以及转速信号包括电流分量id(k)、iq(k)和电角速度ωe(k),以及第k个控制周期施加的电压矢量V(k)在旋转坐标系d-q下的电压分量Vd(k)、Vq(k),通过延时补偿模型,预测得到第(k+1)个控制周期开始时刻的旋转坐标系d-q下电流分量
Figure BDA0001457022540000064
所述的延时补偿模型公式如下:
Figure BDA0001457022540000065
式中Ld,Lq分别为同步电机的直轴和交轴电感,Rs为定子电阻,Ts为控制周期的时长,ψf为转子永磁体磁链,ωe(k)为同步电机电角速度,Vd(k)、Vq(k)分别为第k个控制周期施加的电压矢量V(k)在旋转坐标系d-q下的电压分量,id(k)、iq(k)分别为同步电机旋转坐标系d-q下的电流分量,
Figure BDA0001457022540000066
分别为预测得到的第(k+1)个控制周期开始时刻的旋转坐标系d-q下电流分量,△Vd(k)、△Vq(k)分别由同步电机旋转坐标系d-q下的电流分量id(k)、iq(k)同上个周期预测得到的电流分量
Figure BDA0001457022540000067
做差经过PI控制器得到。
3)将转速信号中的同步电机转子机械角速度ω(k)与给定转速ω*之间的差值,经PI控制器得到同步电机第k个控制周期的电磁转矩给定值Te *;根据电磁转矩给定值Te *,采用最大转矩电流比计算方法,得到旋转坐标系d-q下电流分量的给定值id *和iq *,并据电流分量的给定值id *和iq *求出同步电机定子磁链ψs的给定值和同步电机输出转矩的给定值;
所述的同步电机定子磁链ψs的给定值是采用下式计算:
Figure BDA0001457022540000071
式中,Ld,Lq分别为同步电机的直轴和交轴电感,ψd *和ψq *分别为定子磁链在旋转坐标系d-q下d轴的给定分量和q轴的给定分量,ψf为同步电机转子永磁体磁链。
所述的同步电机输出转矩的给定值是采用下式计算:
Figure BDA0001457022540000072
式中,TM *为同步电机永磁体输出转矩的给定分量,TR *为同步电机磁阻输出转矩的给定分量,Ld,Lq分别为同步电机的直轴和交轴电感,ψf为同步电机转子永磁体磁链,p为同步电机的极对数。
4)将步骤2)预测得到的电流分量,以及内置式永磁同步电机控制系统中的逆变器中的8个电压矢量V0、V1、V2、……、V7,带入到电机预测模型中,预测出不同电压矢量Vn作用下的同步电机磁链和转矩,n=0,1,2……7,在预测过程中,由于电压矢量V0和电压矢量V7的作用效果相同,不考虑电压矢量V0的相关计算;包括:
将预测得到内置式永磁同步电机第(k+1)个控制周期开始时刻的旋转坐标系d-q下的电流分量
Figure BDA0001457022540000073
Figure BDA0001457022540000074
以及所述逆变器中的8个电压矢量V0、V1、V2、……、V7,如图5和表1所示,带入到电机预测模型,预测出不同电压矢量Vn作用下,第(k+2)控制周期开始时刻同步电机定子磁链ψs在旋转坐标系d-q下的分量
Figure BDA0001457022540000075
Figure BDA0001457022540000076
及永磁体转矩
Figure BDA0001457022540000077
和磁阻转矩
Figure BDA0001457022540000078
电机预测模型公式如下,在预测过程中,由于电压矢量V0和电压矢量V7的作用效果相同,不考虑电压矢量V0的相关计算:
Figure BDA0001457022540000079
Figure BDA00014570225400000710
Figure BDA00014570225400000711
式中Ld,Lq分别为同步电机的直轴和交轴电感,Rs为定子电阻,Ts为控制周期的时长,p为同步电机的极对数,
Figure BDA00014570225400000712
分别为预测得到第(k+1)个控制周期开始时刻的旋转坐标系d-q下电流分量,n为和电压矢量Vn相关变量的的序号、Vdn、Vqn分别为逆变器电压矢量Vn在旋转坐标系d-q下的电压分量,
Figure BDA00014570225400000713
分别为电压矢量Vn作用后第(k+2)控制周期开始时刻同步电机定子磁链ψs在旋转坐标系d-q下的分量,
Figure BDA0001457022540000081
Figure BDA0001457022540000082
分别为电压矢量Vn作用后第(k+2)控制周期开始时刻同步电机永磁体输出转矩和同步电机磁阻输出转矩,ψf为同步电机转子永磁体磁链,由于电机的机械时间常数远大于电气时间常数,认定两个相邻控制周期中转速是不变的即式中ωe(k+1)=ωe(k)。
表1空间电压矢量表
Figure BDA0001457022540000083
表中S1、S3、S5分别表示两电平逆变器所对应的三相桥中上桥臂的开关管,S2、S4、S6分别表示两电平逆变器所对应的三相桥中下桥臂的开关管,为1时表示开通,为0时表示关断。
5)将预测出不同电压矢量Vn作用下的电机磁链和转矩,带入到评价函数中,从评价函数的计算结果中找到最小值,所述最小值所对应的逆变器电压矢量,即为最优电压矢量Vopt,并通过逆变器输出所述最优电压矢量;包括:
内置式永磁同步电机在最大转矩电流比控制目标下,转矩由永磁体转矩和磁阻转矩构成。为使转矩波动最小,则将步骤4)求得的永磁体转矩
Figure BDA0001457022540000084
和磁阻转矩
Figure BDA0001457022540000085
与同步电机永磁体输出转矩的给定分量TM *和同步电机磁阻输出转矩的给定分量TR *的误差建立同步电机在重载情况衡量电压矢量Vn对控制效果的评价函数g(n)公式如下:
Figure BDA0001457022540000086
根据评价函数g(n)的计算结果,找到评价函数g(n)的最小值g(n)min,则与所述最小值g(n)min对应的电压矢量Vn为最优电压矢量Vopt
由于当电机运行在轻载时,转矩评价函数容易造成较大的电流波动,所以在同步电机轻载时需要使用磁链矢量评价函数以避免权重系数整定和电流波动,这时,评价函数g(n)是采用下式得到:
Figure BDA0001457022540000087
式中,ψd *和ψq *分别为定子磁链在旋转坐标系d-q下d轴的给定分量和q轴的给定分量,
Figure BDA0001457022540000088
Figure BDA0001457022540000089
为第(k+2)控制周期开始时刻同步电机定子磁链ψs在旋转坐标系d-q下的分量;
为使同步电机在重载情况下与轻载情况下的评价函数g(n)正常切换,设定常数TX为切换参考转矩,当同步电机的电磁转矩给定值Te *的绝对值大于TX时,认定为重载,当同步电机的电磁转矩给定值Te *的绝对值小于等于TX时认定为轻载;
为避免当同步电机的电磁转矩给定值Te *的绝对值接近于TX且上下波动时,导致同步电机在重载情况下与轻载情况下的评价函数g(n)切换频繁,故加入滞环比较器减少切换频率:设定环宽为σ,当同步电机的电磁转矩给定值Te *的绝对值大于设定常数TX与设定环宽σ之和,即同步电机的电磁转矩给定值Te *的绝对值处于滞环之上时,滞环输出δ(k)=1,为重载时的评价函数g(n);当同步电机的电磁转矩给定值Te *的绝对值小于设定常数TX与设定环宽σ之差时,即同步电机的电磁转矩给定值Te *的绝对值处于滞环之下时,滞环输出δ(k)=-1,为轻载时的评价函数g(n);当TX-σ≤|Te *|≤TX+σ,即同步电机的电磁转矩给定值Te *的绝对值处于滞环之内时,滞环输出等于上一个控制周期的输出,即δ(k)=δ(k-1),使用上一个控制周期使用的评价函数g(n);δ的初值δ(0)=-1;
从而,所述评价函数表示为:
Figure BDA0001457022540000091
式中
Figure BDA0001457022540000092
分别为电压矢量Vn作用后第(k+2)控制周期开始时刻同步电机定子磁链ψs在旋转坐标系d-q下的分量,
Figure BDA0001457022540000093
分别为电压矢量Vn作用后第(k+2)控制周期开始时刻同步电机永磁体输出转矩和同步电机磁阻输出转矩,ψd *和ψq *分别为定子磁链在旋转坐标系d-q下d轴的给定分量和q轴的给定分量,TM *为同步电机永磁体输出转矩的给定分量,TR *为同步电机磁阻输出转矩的给定分量,argming(n)为求数列{g(1),g(2),……,g(n)}中最小值所在列数的运算,即从评价函数的计算结果中找出最小值所对应的电压矢量。

Claims (7)

1.一种内置式永磁同步电机模型预测转矩控制方法,是用于内置式永磁同步电机控制系统的控制方法,其特征在于,包括如下步骤:
1)采集内置式永磁同步电机第k个控制周期开始时刻的电压电流信号以及转速信号,并进行坐标变换;
2)通过延时补偿模型,预测得到内置式永磁同步电机第(k+1)个控制周期开始时刻的旋转坐标系d-q下的电流分量;
3)将转速信号中的同步电机转子机械角速度ω(k)与给定转速ω*之间的差值,经PI控制器得到同步电机第k个控制周期的电磁转矩给定值Te *;根据电磁转矩给定值Te *,采用最大转矩电流比计算方法,得到旋转坐标系d-q下电流分量的给定值id *和iq *,并据电流分量的给定值id *和iq *求出同步电机定子磁链ψs的给定值和同步电机输出转矩的给定值;
4)将步骤2)预测得到的电流分量,以及内置式永磁同步电机控制系统中的逆变器中的8个电压矢量V0、V1、V2、……、V7,带入到电机预测模型中,预测出不同电压矢量Vn作用下的同步电机磁链和转矩,n=0,1,2……7,在预测过程中,由于电压矢量V0和电压矢量V7的作用效果相同,不考虑电压矢量V0的相关计算;
5)将预测出不同电压矢量Vn作用下的电机磁链和转矩,带入到评价函数中,从评价函数的计算结果中找到最小值,所述最小值所对应的逆变器电压矢量,即为最优电压矢量Vopt,并通过逆变器输出所述最优电压矢量;包括:
将步骤4)求得的永磁体转矩
Figure FDA0002302516030000011
和磁阻转矩
Figure FDA0002302516030000012
与同步电机永磁体输出转矩的给定分量TM *和同步电机磁阻输出转矩的给定分量TR *的误差建立同步电机在重载情况衡量电压矢量Vn对控制效果的评价函数g(n)公式如下:
Figure FDA0002302516030000013
根据评价函数g(n)的计算结果,找到评价函数g(n)的最小值g(n)min,则与所述最小值g(n)min对应的电压矢量Vn为最优电压矢量Vopt
在同步电机轻载时需要使用磁链矢量评价函数以避免权重系数整定和电流波动,这时,评价函数g(n)是采用下式得到:
Figure FDA0002302516030000014
式中,ψd *和ψq *分别为定子磁链在旋转坐标系d-q下d轴的给定分量和q轴的给定分量,
Figure FDA0002302516030000015
Figure FDA0002302516030000016
为第(k+2)控制周期开始时刻同步电机定子磁链ψs在旋转坐标系d-q下的分量;
为使同步电机在重载情况下与轻载情况下的评价函数g(n)正常切换,设定常数TX为切换参考转矩,当同步电机的电磁转矩给定值Te *的绝对值大于TX时,认定为重载,当同步电机的电磁转矩给定值Te *的绝对值小于等于TX时认定为轻载;
为避免当同步电机的电磁转矩给定值Te *的绝对值接近于TX且上下波动时,导致同步电机在重载情况下与轻载情况下的评价函数g(n)切换频繁,故加入滞环比较器减少切换频率:设定环宽为σ,当同步电机的电磁转矩给定值Te *的绝对值大于设定常数TX与设定环宽σ之和,即同步电机的电磁转矩给定值Te *的绝对值处于滞环之上时,滞环输出δ(k)=1,为重载时的评价函数g(n);当同步电机的电磁转矩给定值Te *的绝对值小于设定常数TX与设定环宽σ之差时,即同步电机的电磁转矩给定值Te *的绝对值处于滞环之下时,滞环输出δ(k)=-1,为轻载时的评价函数g(n);当TX-σ≤|Te *|≤TX+σ,即同步电机的电磁转矩给定值Te *的绝对值处于滞环之内时,滞环输出等于上一个控制周期的输出,即δ(k)=δ(k-1),使用上一个控制周期使用的评价函数g(n);δ的初值δ(0)=-1;
从而,所述评价函数表示为:
Figure FDA0002302516030000021
式中
Figure FDA0002302516030000022
分别为电压矢量Vn作用后第(k+2)控制周期开始时刻同步电机定子磁链ψs在旋转坐标系d-q下的分量,
Figure FDA0002302516030000023
分别为电压矢量Vn作用后第(k+2)控制周期开始时刻同步电机永磁体输出转矩和同步电机磁阻输出转矩,ψd *和ψq *分别为定子磁链在旋转坐标系d-q下d轴的给定分量和q轴的给定分量,TM *为同步电机永磁体输出转矩的给定分量,TR *为同步电机磁阻输出转矩的给定分量,argming(n)为求数列{g(1),g(2),……,g(n)}中最小值所在列数的运算,即从评价函数的计算结果中找出最小值所对应的电压矢量。
2.根据权利要求1所述的一种内置式永磁同步电机模型预测转矩控制方法,其特征在于,步骤1)包括:
通过主控板上微处理器内部的模数转换接口采集第k个控制周期开始时刻的逆变器直流母线电压udc(k)和三相定子电流ia(k)、ib(k)、ic(k);
利用编码器得到第k个控制周期开始时刻的同步电机转子机械角速度ω(k)和电角速度ωe(k)以及转子位置电角度θ(k);
将采集到的三相定子电流ia(k)、ib(k)、ic(k),通过三相/两相静止坐标变换得到两相静止坐标系α-β下电流分量iα(k)和iβ(k),变换公式如下:
Figure FDA0002302516030000024
再通过两相旋转坐标变换得到旋转坐标系d-q下的电流分量id(k)和iq(k),变换公式如下:
Figure FDA0002302516030000025
式中ia(k)、ib(k)、ic(k)分别为同步电机三相定子电流,iα(k)、iβ(k)分别为同步电机两相静止坐标系α-β下的电流分量,id(k)、iq(k)分别为同步电机旋转坐标系d-q下的电流分量,θ(k)为同步电机转子位置电角度。
3.根据权利要求1所述的一种内置式永磁同步电机模型预测转矩控制方法,其特征在于,步骤2)包括:利用坐标变换后的电压电流信号以及转速信号包括电流分量id(k)、iq(k)和电角速度ωe(k),以及第k个控制周期施加的电压矢量V(k)在旋转坐标系d-q下的电压分量Vd(k)、Vq(k),通过延时补偿模型,预测得到第(k+1)个控制周期开始时刻的旋转坐标系d-q下电流分量
Figure FDA0002302516030000031
4.根据权利要求3所述的一种内置式永磁同步电机模型预测转矩控制方法,其特征在于,所述的延时补偿模型公式如下:
Figure FDA0002302516030000032
式中Ld,Lq分别为同步电机的直轴和交轴电感,Rs为定子电阻,Ts为控制周期的时长,ψf为转子永磁体磁链,ωe(k)为同步电机电角速度,Vd(k)、Vq(k)分别为第k个控制周期施加的电压矢量V(k)在旋转坐标系d-q下的电压分量,id(k)、iq(k)分别为同步电机旋转坐标系d-q下的电流分量,
Figure FDA0002302516030000033
分别为预测得到的第(k+1)个控制周期开始时刻的旋转坐标系d-q下电流分量,ΔVd(k)、ΔVq(k)分别由同步电机旋转坐标系d-q下的电流分量id(k)、iq(k)同上个周期预测得到的电流分量
Figure FDA0002302516030000034
做差经过PI控制器得到。
5.根据权利要求1所述的一种内置式永磁同步电机模型预测转矩控制方法,其特征在于,步骤3)所述的同步电机定子磁链ψs的给定值是采用下式计算:
Figure FDA0002302516030000035
式中,Ld,Lq分别为同步电机的直轴和交轴电感,ψd *和ψq *分别为定子磁链在旋转坐标系d-q下d轴的给定分量和q轴的给定分量,ψf为同步电机转子永磁体磁链。
6.根据权利要求1所述的一种内置式永磁同步电机模型预测转矩控制方法,其特征在于,步骤3)所述的同步电机输出转矩的给定值是采用下式计算:
Figure FDA0002302516030000036
式中,TM *为同步电机永磁体输出转矩的给定分量,TR *为同步电机磁阻输出转矩的给定分量,Ld,Lq分别为同步电机的直轴和交轴电感,ψf为同步电机转子永磁体磁链,p为同步电机的极对数。
7.根据权利要求1所述的一种内置式永磁同步电机模型预测转矩控制方法,其特征在于,步骤4)包括:
将预测得到的内置式永磁同步电机第(k+1)个控制周期开始时刻的旋转坐标系d-q下的电流分量
Figure FDA0002302516030000041
Figure FDA0002302516030000042
以及所述逆变器中的8个电压矢量V0、V1、V2、……、V7,带入到电机预测模型,预测出不同电压矢量Vn作用下,第(k+2)控制周期开始时刻同步电机定子磁链ψs在旋转坐标系d-q下的分量
Figure FDA0002302516030000043
Figure FDA0002302516030000044
及永磁体转矩
Figure FDA0002302516030000045
和磁阻转矩
Figure FDA0002302516030000046
电机预测模型公式如下,在预测过程中,由于电压矢量V0和电压矢量V7的作用效果相同,不考虑电压矢量V0的相关计算:
Figure FDA0002302516030000047
Figure FDA0002302516030000048
Figure FDA0002302516030000049
式中Ld,Lq分别为同步电机的直轴和交轴电感,Rs为定子电阻,Ts为控制周期的时长,p为同步电机的极对数,
Figure FDA00023025160300000410
分别为预测得到第(k+1)个控制周期开始时刻的旋转坐标系d-q下电流分量,n为和电压矢量Vn相关变量的的序号、Vdn、Vqn分别为逆变器电压矢量Vn在旋转坐标系d-q下的电压分量,
Figure FDA00023025160300000411
分别为电压矢量Vn作用后第(k+2)控制周期开始时刻同步电机定子磁链ψs在旋转坐标系d-q下的分量,
Figure FDA00023025160300000412
Figure FDA00023025160300000413
分别为电压矢量Vn作用后第(k+2)控制周期开始时刻同步电机永磁体输出转矩和同步电机磁阻输出转矩,ψf为同步电机转子永磁体磁链,由于电机的机械时间常数远大于电气时间常数,认定两个相邻控制周期中转速是不变的即式中ωe(k+1)=ωe(k)。
CN201711075880.9A 2017-11-03 2017-11-03 一种内置式永磁同步电机模型预测转矩控制方法 Active CN107612446B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711075880.9A CN107612446B (zh) 2017-11-03 2017-11-03 一种内置式永磁同步电机模型预测转矩控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711075880.9A CN107612446B (zh) 2017-11-03 2017-11-03 一种内置式永磁同步电机模型预测转矩控制方法

Publications (2)

Publication Number Publication Date
CN107612446A CN107612446A (zh) 2018-01-19
CN107612446B true CN107612446B (zh) 2020-06-02

Family

ID=61085393

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711075880.9A Active CN107612446B (zh) 2017-11-03 2017-11-03 一种内置式永磁同步电机模型预测转矩控制方法

Country Status (1)

Country Link
CN (1) CN107612446B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108448976B (zh) * 2018-03-21 2019-09-27 华中科技大学 一种永磁同步电机最大转矩电流比控制装置
FR3089368A1 (fr) * 2018-11-30 2020-06-05 IFP Energies Nouvelles Procédé de commande d’une machine tournante triphasée et système de commande associé
CN109617481A (zh) * 2018-12-25 2019-04-12 南京越博电驱动系统有限公司 一种永磁同步电机预测转矩控制方法
CN110445441B (zh) * 2019-06-29 2021-07-27 天津大学 一种永磁同步电机预测转矩控制方法
CN110880891A (zh) * 2019-11-25 2020-03-13 沈阳工程学院 一种基于模型预测的多绕组永磁风力发电机效率优化方法
CN112751513B (zh) * 2020-12-24 2022-10-11 珠海格力电器股份有限公司 一种电机的控制方法、装置、电机、存储介质及处理器
CN113809958B (zh) * 2021-08-05 2023-11-03 杭州电子科技大学 一种内置式永磁同步电机的最大再生功率制动控制方法
CN113640667B (zh) * 2021-08-20 2023-12-19 东风汽车集团股份有限公司 一种电机eol下线零点自动校准方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103516286A (zh) * 2013-08-09 2014-01-15 天津大学 一种可改善输入输出性能的矩阵变换器直接转矩控制方法
CN103944459A (zh) * 2014-04-09 2014-07-23 天津大学 一种适用于多电机系统的直接转矩控制方法
CN104953914A (zh) * 2015-07-10 2015-09-30 泉州装备制造研究所 一种异步电机无传感器预测转矩控制系统及方法
CN105548792A (zh) * 2015-12-28 2016-05-04 中南大学 基于预测控制的矩阵变换器开关开路故障诊断方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103516286A (zh) * 2013-08-09 2014-01-15 天津大学 一种可改善输入输出性能的矩阵变换器直接转矩控制方法
CN103944459A (zh) * 2014-04-09 2014-07-23 天津大学 一种适用于多电机系统的直接转矩控制方法
CN104953914A (zh) * 2015-07-10 2015-09-30 泉州装备制造研究所 一种异步电机无传感器预测转矩控制系统及方法
CN105548792A (zh) * 2015-12-28 2016-05-04 中南大学 基于预测控制的矩阵变换器开关开路故障诊断方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Analysis of One Step Delay in Direct Torque Control of Permanent Magnet Synchronous Motor and Its Remedies;Yongchang Zhang, Jianguo Zhu, Wei Xu;《International Conference on Electrical Machines and Systems. IEEE》;20101231;792-797 *
电动汽车用永磁同步电机改进预测电流控制;陈晨 等;《新型工业化》;20170920;第7卷(第9期);22-30 *

Also Published As

Publication number Publication date
CN107612446A (zh) 2018-01-19

Similar Documents

Publication Publication Date Title
CN107612446B (zh) 一种内置式永磁同步电机模型预测转矩控制方法
CN107317532B (zh) 基于滑模的永磁同步电机预测电流控制方法和系统
CN106936356B (zh) 矢量筛选和占空比结合的电机模型预测控制系统及方法
CN106788075B (zh) 基于改进欧拉法的快速矢量筛选预测转矩控制方法
Lin et al. Improved model-free predictive current control for synchronous reluctance motor drives
CN107154763B (zh) 永磁同步电机无差拍直接转矩控制系统及控制方法
CN106788045B (zh) 一种永磁同步电机模型预测pi动态权重并行控制方法
CN108448982B (zh) 一种基于空间电压矢量预测的直接转矩控制方法
CN109617481A (zh) 一种永磁同步电机预测转矩控制方法
CN110995072B (zh) 一种电机转子位置估测方法
CN108512473B (zh) 三相四开关永磁同步电机调速系统的直接转矩控制方法
Abu-Rub et al. Speed observer system for advanced sensorless control of induction motor
CN113992098A (zh) 一种永磁同步电机稳态性能的提升方法、系统
CN112910359A (zh) 一种改进型永磁同步直线电机模型预测电流控制方法
CN110601631A (zh) 一种基于占空比调制的磁通切换型永磁直线电机直接推力控制方法
Tatemarsu et al. New approaches with sensorless drives
Schmidt et al. Sensorless control of a permanent magnet synchronous machine down to near zero speed applied to position motion control
CN116073720A (zh) 一种改善永磁同步电机预测控制系统性能的装置及方法
CN115514279A (zh) 一种电机脉宽调制优化驱动方法及系统
Yin et al. Sliding mode control for servo motors based on the differential evolution algorithm
Wang et al. Sensorless predictive control for an induction machine
CN110601629B (zh) 一种提高pmsm伺服系统电流环响应速度并且减小震荡的方法
Chang et al. Improved rotor position detection speed based on optimized BP neural network for permanent magnet synchronous motor
CN110808704B (zh) 一种高开关频率逆变器的低控制频率控制方法
Geethanjali et al. Direct torque control of induction motor drive using tms320f28335 digital signal controller

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