CN105071717B - 利用电流谐波抑制表贴式永磁同步电机转矩波动的方法 - Google Patents
利用电流谐波抑制表贴式永磁同步电机转矩波动的方法 Download PDFInfo
- Publication number
- CN105071717B CN105071717B CN201510493654.7A CN201510493654A CN105071717B CN 105071717 B CN105071717 B CN 105071717B CN 201510493654 A CN201510493654 A CN 201510493654A CN 105071717 B CN105071717 B CN 105071717B
- Authority
- CN
- China
- Prior art keywords
- harmonic
- current
- electromotive force
- electromagnetic torque
- amplitude
- 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.)
- Expired - Fee Related
Links
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 69
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000004804 winding Methods 0.000 claims abstract description 34
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 8
- 238000001228 spectrum Methods 0.000 claims abstract description 8
- 230000009471 action Effects 0.000 claims description 14
- 230000002401 inhibitory effect Effects 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 5
- 238000002347 injection Methods 0.000 claims description 4
- 239000007924 injection Substances 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 239000010410 layer Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000002356 single layer Substances 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 241001270131 Agaricus moelleri Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002829 reductive effect Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Landscapes
- Control Of Motors That Do Not Use Commutators (AREA)
Abstract
本发明公开了一种利用电流谐波抑制表贴式永磁同步电机转矩波动的方法,包括:将永磁同步电机反电动势波形进行傅里叶分解,获得反电动势谐波频谱和各次谐波的幅值大小;将电流表示成傅里叶级数形式,基于三相反电动势和电流表达式求解电磁转矩模型;将电磁转矩根据谐波成分不同进行分类,求解各类电磁转矩波动最小时电流谐波满足的条件;根据反电动势谐波特点判断电磁转矩波动主要来源,确定电流谐波次数、幅值和相位;将含有特定谐波的电流通入永磁同步电机三相绕组中,抑制电机电磁转矩波动。本方法能够有效抑制表贴式永磁同步电机电磁转矩波动,并为永磁电机转矩波动抑制提供了一种新的思路,可进一步用于其他类型永磁电机的转矩波动抑制中。
Description
技术领域
本发明涉及表贴式永磁同步电机转矩波动抑制技术领域,尤其涉及一种利用电流谐波抑制表贴式永磁同步电机转矩波动的方法。
背景技术
表贴式永磁同步电机具有较高的效率、功率因数和转矩密度,因而被广泛应用到电梯、风电、混合动力汽车等众多工业领域中。对于表贴式永磁同步电机来说,理论上如果电机的反电动势和电流均是理想的正弦波,则电机可产生恒定的电磁转矩。然而实际上,一方面由于电机本身所采用绕组的分布形式、齿槽效应、磁极布置等因素,使得电机气隙磁场发生畸变,产生磁链谐波,并进一步导致电机反电动势产生谐波;另一方面,由于变频器所采用的电力电子器件的非线性特性,如开关管的管压降、死区时间等,导致电机产生电流谐波。反电动势谐波与电流谐波相互作用,导致永磁同步电机产生转矩波动,严重影响电机系统的控制精度和运行平稳性,因而成为电机设计和控制重点解决问题之一。
为了解决永磁同步电机转矩波动的问题,学者们进行了大量研究。1)针对电机本体结构进行优化设计,尽可能地改善气隙磁场分布,减少电机本身产生的反电动势谐波,如优化转子磁极结构、改变永磁体充磁方式、采用斜槽或斜极等方法。但该类方法不仅增加了电机加工的难度和成本,而且由于加工误差的存在往往难以达到理想效果,无法完全消除反电动势谐波。2)针对变流器非线性特性产生的电流谐波通过控制策略补偿或抑制,如无死区开关控制、时间补偿、电流反馈型电压补偿等方法。而这类方法需要增加新的硬件,算法较为复杂,且由于电流过零点检测不准,容易导致误补偿。
近年来,人们渐渐意识到通过单纯对反电动势谐波和电流谐波进行抑制难以完全消除永磁同步电机的转矩波动,且在抑制谐波的过程中往往需要牺牲系统的某些其他性能,因此有专家提出从谐波利用的角度对电机系统的性能进行提升,为永磁同步发电机的转矩波动抑制提供了一种新的思路。但永磁同步电机转矩波动的产生是反电动势和电流谐波共同作用的结果,且反电动势和电流的谐波成分较为复杂,若要利用电流谐波来达到抑制转矩波动的目的,就必须对两种谐波之间的作用机理进行深入分析。由于目前关于谐波作用机理方面的研究较少,导致电流谐波利用技术受到了较大的限制。
发明内容
本发明基于表贴式永磁同步电机三相绕组反电动势与电流谐波之间的作用机理,提供了一种利用电流谐波抑制表贴式永磁同步电机电磁转矩波动的方法。
本发明采用的技术方案如下:
利用电流谐波抑制表贴式永磁同步电机转矩波动的方法,包括以下步骤:
1)对表贴式永磁同步电机反电动势波形进行傅里叶分解,获得反电动势谐波的分布频谱和各次谐波的幅值大小;
2)将非正弦电流表示成傅里叶级数形式,基于三相反电动势和电流傅里叶级数表达式求解电磁转矩模型;
3)根据电磁转矩谐波成分不同进行分类,求解各类电磁转矩波动最小时电流谐波满足的条件;
4)根据反电动势谐波特点判断电磁转矩波动的主要来源,确定电流谐波次数、幅值和相位;
5)将含有所述步骤4)确定的含有谐波的三相绕组电流通入表贴式永磁同步电机三相绕组中,抑制电机电磁转矩波动。
前述的步骤1)中,对表贴式永磁同步电机反电动势波形进行傅里叶分解,得到:
其中,分别为表贴式永磁同步电机定子三相绕组对应的反电动势;E1为反电动势基波幅值;ω为电机电角速度;t为时间;v为反电动势谐波次数,且v≠1;Ev为v次反电动势谐波分量对应的幅值;αv为v次反电动势谐波分量的初始相位角;kev为v次谐波反电动势向量旋转方向的系数,其值为1 或-1,kev=1代表v次谐波反电动势向量旋转方向与基波旋转方向相同,kev=-1代表v次谐波反电动势向量旋转方向与基波旋转方向相反。
前述的步骤2)中,将非正弦电流表示成傅里叶级数形式如下:
其中,分别为三相绕组中的电流;I1为电流基波幅值;ψ1为电机内功率因数角;u为电流谐波次数,且u≠1;Iu为u次电流谐波分量对应的幅值;βu为u次电流谐波分量的初始相位角;kiu表示 u次谐波电流相量旋转方向的系数,其值为1或-1,kiu=1代表u次谐波电流相量旋转方向与基波旋转方向相同,kiu=-1代表u次谐波电流相量旋转方向与基波旋转方向相反。
前述的步骤2)中,电磁转矩模型为:
其中,Tem为电磁转矩,ωr为电机旋转机械角速度,ωr=2πf/p;f为绕组电压和电流频率;p为极对数。
前述的步骤3)中,根据电磁转矩谐波成分不同进行分类,如下:
将永磁同步电机非正弦的三相反电动势和电流代入电磁转矩表达式式(3)中,可得:
Tem=T1+T2+T3+T4 (4)
其中,
即将电磁转矩分为T1,T2,T3,T4四类,T1是由电机基波反电动势和基波电流作用产生的恒定的电磁转矩;T2是由谐波反电动势和基波电流作用产生的电磁转矩分量;T3是基波反电动势与谐波电流作用产生的电磁转矩分量;T4是由谐波反电动势和谐波电流产生转矩分量。
前述的求解各类电磁转矩波动最小时电流谐波满足的条件,如下:
由于T2和T3特征相似,且单独不可能为0,因此将T2+T3作为一组,T4单独作为一组单独进行分析:
1)若令T2+T3=0,即
可知,当任意v=u且kev=kiu时,若
且βu=±π-ψ1+αv
即当电流所含谐波分量的次数与反电动势谐波次数相同,电流谐波幅值比例与反电动势谐波幅值比例一致,且任意次电流谐波初始相位角满足βu=±π-ψ1+αv时,T2与T3两项产生的转矩波动可完全消除;
2)若令T4=0且kevv≠kiuu,即
当kevv=kiuu时,T4中的各分量均不为0,无法自行消除,因此仅能通过对应谐波产生的转矩分量进行相互抵消,即令任意两次反电动势谐波v1、v2分量和相应两次电流谐波u1、u2分量引起的转矩波动之和为 0,如下:
当任意u1=v1、v2=u2且kev1=kiu1、kev2=kiu2时,若满足条件
且
即当任意两次电流所含谐波分量的次数与相应两次反电动势谐波次数相同,该两次电流谐波幅值比例与反电动势谐波幅值比例一致,且该两次电流谐波初始相位角满足上述关系时,T4产生的转矩波动可被完全消除,此时电磁转矩的平均值将保持不变。
前述的步骤4)根据反电动势谐波特点判断电磁转矩波动的主要来源,确定电流谐波次数、幅值和相位,具体如下:
由于永磁同步电机三相绕组结构对称,在定子圆周上互相错开120°电角度分布,反电动势中不含有偶数次以及3的倍数次谐波,因此v=5、7、11、13……,即反电动势仅含有6k±1次谐波,其中k=1,2, 3……,且6k-1次谐波反电动势向量的旋转方向与基波相反,kev=-1,而6k+1次谐波反电动势向量的旋转方向与基波相同,kev=1;同样,电流谐波次数也应为u=5、7、11、13……,且当u=6k-1时kiu=-1,当 u=6k+1,时kiu=1;
由于各电磁转矩波动分量大小均与反电动势谐波幅值和电流谐波幅值的乘积成正比,而一般情况下,反电动势和电流基波的幅值要远大于谐波幅值,且谐波中5、7次谐波幅值最大,因此由基波与5、7次谐波产生的电磁转矩波动分量要大于其余电磁转矩波动分量,所以判断反电动势中5、7次谐波和基波产生电磁转矩波动为主要来源,
在此基础上,给定三相绕组电流为:
其中,
I5=E5I1/E1;β5=π-ψ1+α;I7=E7I1/E1;β7=π-ψ1+α7。
前述的步骤5),将含有谐波的三相绕组电流通入表贴式永磁同步电机三相绕组中,具体为,
将式(12)进行坐标变换,得到d轴电流分量Id与q轴电流分量Iq,分别为:
然后,将式(13)给定到电机控制器中,再通过电流闭环控制使得电机定子电流中注入适当5次和7 次谐波,实现谐波电流注入。
本发明的有益效果是:本发明建立了非正弦反电动势和电流作用下表贴式永磁同步电机电磁转矩数学模型,根据基波与谐波组合的不同将电磁转矩分为4种类型的分量,研究电磁转矩波动最小时反电动势谐波和电流谐波满足的条件,并结合表贴式永磁同步电机反电动势谐波分布的特点,重点对5、7次谐波产生的电磁转矩波动分量进行抑制,简化电流谐波方案。本发明为永磁同步电机研究人员分析非正弦反电动势情况下电机电磁转矩特性提供了一种研究手段,为表贴式永磁同步电机转矩波动抑制提供了一种新的方法,并为电流谐波利用技术奠定了一定理论基础。
附图说明
图1为本发明的利用电流谐波抑制表贴式永磁同步电机转矩波动方法的流程图;
图2为本发明的实施例中一台4极6槽双层绕组表贴式永磁同步电机有限元模型;
图3(a)为图2的4极6槽永磁同步电机的三相反电动势波形图;图3(b)为图2的4极6槽永磁同步电机经傅里叶分解后的反电动势谐波频谱图;
图4(a)为图2的4极6槽永磁同步电机正弦电流与注入5、7谐波的电流对比图;图4(b)为图2 的4极6槽永磁同步电机在图4(a)两种电流波形下产生的电磁转矩波形对比图。
图5为本发明的实施例中一台4极12槽单层绕组表贴式永磁同步电机有限元模型;
图6(a)为图5的4极12槽永磁同步电机的三相反电动势波形图;图6(b)为图5的4极12槽永磁同步电机经傅里叶分解后的反电动势谐波频谱图;
图7(a)为图5的4极12槽永磁同步电机正弦电流与注入5、7谐波的电流对比图;图7(b)为图5 的4极12槽永磁同步电机在图7(a)两种电流波形下产生的电磁转矩波形对比图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面结合附图和具体实施方式对本发明作进一步地详细描述。
如图1所示,本发明的利用电流谐波抑制表贴式永磁同步电机转矩波动的方法包括以下步骤:
(1)对表贴式永磁同步电机反电动势波形进行傅里叶分解,获得反电动势谐波的分布频谱和各次谐波的幅值大小;具体为:
式中,分别为表贴式永磁同步电机定子三相绕组对应的反电动势;E1为反电动势基波幅值;ω为电机电角速度,ω=2πf;t为时间;v为反电动势谐波次数,且v≠1;Ev为v次反电动势谐波分量对应的幅值;αv为v次反电动势谐波分量的初始相位角;kev为v次谐波反电动势向量旋转方向的系数,其值可为1或-1,kev=1代表v次谐波反电动势向量旋转方向与基波旋转方向相同,kev=-1代表v次谐波反电动势向量旋转方向与基波旋转方向相反。
(2)将非正弦电流表示成傅里叶级数形式,基于三相反电动势和电流傅里叶级数表达式求解电磁转矩模型;具体为:
参照反电动势谐波表达式,非正弦电流可表示为:
式中,分别为三相绕组中的电流;I1为电流基波幅值;ψ1为电机内功率因数角;u为电流谐波次数,且u≠1;Iu为u次电流谐波分量对应的幅值;βu为u次电流谐波分量的初始相位角;kiu表示 u次谐波电流相量旋转方向的系数,其值可为1或-1,kiu=1代表u次谐波电流相量旋转方向与基波旋转方向相同,kiu=-1代表u次谐波电流相量旋转方向与基波旋转方向相反。
基于三相绕组反电动势和电流的傅里叶级数表达式得到永磁同步电机电磁转矩Tem模型如下:
式中,ωr为电机旋转机械角速度,ωr=2πf/p;f为绕组电压和电流频率;p为极对数。
(3)根据电磁转矩谐波成分不同进行分类,求解各类电磁转矩波动最小时电流谐波满足的条件。具体为:
将永磁同步电机非正弦的三相反电动势和电流代入电磁转矩表达式中,可得:
Tem=T1+T2+T3+T4 (4)
其中,
式中,根据电磁转矩谐波组成不同将其分量分为T1,T2,T3,T4四类。T1是由电机基波反电动势和基波电流作用产生的恒定的电磁转矩;T2是由谐波反电动势和基波电流作用产生的电磁转矩分量;T3是基波反电动势与谐波电流作用产生的电磁转矩分量;T4是由谐波反电动势和谐波电流产生转矩分量。
进一步对电磁转矩波动中的T2、T3和T4分量进行分析,探索非正弦反电动势情况下使得电磁转矩波动分量最小时电流波形满足条件。由于T2和T3特征较为相似,且单独不可能为0,因此将(T2+T3)作为一组,T4单独作为一组单独进行分析。
1)若令T2+T3=0,即
可知,当任意v=u且kev=kiu时,若
且βu=±π-ψ1+αv
即当电流所含谐波分量的次数与反电动势谐波次数相同,电流谐波幅值比例与反电动势谐波幅值比例一致,且任意次电流谐波初始相位角满足βu=±π-ψ1+αv时,T2与T3两项产生的转矩波动可完全消除。
2)若令T4=0(kevv≠kiuu),即
当kevv=kiuu时,T4中的各分量均不为0,无法自行消除,因此仅能通过对应谐波产生的转矩分量进行相互抵消,即令任意两次反电动势谐波v1、v2分量和相应两次电流谐波u1、u2分量引起的转矩波动之和为 0。
当任意u1=v1、v2=u2且kev1=kiu1、kev2=kiu2时,若满足条件
且
即当任意两次电流所含谐波分量的次数与相应两次反电动势谐波次数相同,该两次电流谐波幅值比例与反电动势谐波幅值比例一致,且该两次电流谐波初始相位角满足上述关系时,T4产生的转矩波动可被完全消除,此时电磁转矩的平均值将保持不变。
(4)根据反电动势谐波特点判断电磁转矩波动的主要来源,确定电流谐波次数、幅值和相位。具体如下:
步骤(3)中的两种情况均希望电流谐波的次数与反电动势谐波次数相同,且电流谐波幅值大小与反电动势谐波幅值成比例。由于永磁同步电机三相绕组结构对称,在定子圆周上互相错开120°电角度分布,反电动势中不含有偶数次以及3的倍数次谐波,因此v=5、7、11、13……,即反电动势仅含有6k±1 次谐波,其中k=1,2,3……。且6k-1次谐波反电动势向量的旋转方向与基波相反,kev=-1;而6k+1次谐波反电动势向量的旋转方向与基波相同,kev=1。因此,可得电流谐波次数也应为u=5、7、11、 13……,且当u=6k-1时kiu=-1;当u=6k+1,时kiu=1。各次电流谐波幅值大小与反电动势谐波幅值成正比。
但是,由于一般情况下ψ1不能为±π/2,否则T1为0,电机将无法产生恒定的电磁转矩。对比步骤(3) 中的两种情况可以发现,电磁转矩波动分量最小时电流初始相位角并不一致,无法同时满足。另外,由于反电动势中含有大量的谐波,若要根据所有反电动势谐波分量去调制出对应的电流谐波分量,实现难度较大,若调制方法不当甚至还可能会产生相反的效果,恶化永磁同步电机的运行性能。
为解决这两个问题,进一步对电磁转矩各分量进行分析。可以发现各电磁转矩波动分量大小均与反电动势谐波幅值和电流谐波幅值的乘积成正比。而一般情况下,反电动势和电流的基波幅值要远大于谐波幅值,且谐波中5、7次谐波幅值较大,因此由基波与5、7次谐波产生的电磁转矩波动分量要大于其他电磁转矩波动分量。为简化谐波电流方案,本发明仅针对反电动势中5、7次谐波和基波产生电磁转矩波动进行抑制,给定三相绕组电流应为
其中,
I5=E5I1/E1;β5=π-ψ1+α;I7=E7I1/E1;β7=π-ψ1+α7
(5)将含有特定谐波的电流通入表贴式永磁同步电机三相绕组中,使电机电磁转矩波动得到有效抑制。具体为:
将步骤(4)中获得的含有谐波的三相绕组电流进行坐标变换得到d轴电流分量Id与q轴电流分量Iq,分别为:
将上式给定到电机控制器中,再通过电流闭环控制使得电机定子电流中注入适当5次和7次谐波,实现谐波电流注入方案,达到抑制电磁转矩波动的目的。
虽然,所描述的实施例仅是针对表贴式永磁同步电机,而本发明也可应用到含有其他类型永磁电机中。
实施例1
101:图2为一台4极6槽双层绕组表贴式永磁同步电机有限元模型。电机额定频率为50Hz,额定转速1500r/min,额定电流有效值为3A,即I1=4.24A。表贴式永磁同步电机通常采用id=0控制方式进行控制,即ψ1=0。
102:通过有限元仿真可以获得该电机的三相反电动势波形,如图3(a)所示。进一步通过傅里叶变换获得反电动势各次谐波的频谱分布以及各次谐波幅值,如图3(b)所示。
103:由图3(b)可以看到反电动势波形中5、7次为其主要谐波,而且他谐波的幅值均小于基波幅值的2%,因此可忽略其他谐波,重点对5、7谐波产生的电磁转矩波动进行抑制。在考虑5、7次谐波初始相位角的基础上,该电机的三相反电动势可表示为:
从图3(b)中,可以得出,反电动势基波幅值E1为161.5,5次谐波分量对应的幅值E5为21.0,7次谐波分量对应的幅值E7为8.2。
104:根据本发明步骤2和步骤3中的分析,为了对5、7次反电动势谐波与电流基波产生的电磁转矩波动进行抑制,需要满足电流所含谐波分量的次数与反电动势谐波次数相同,电流谐波幅值比例与反电动势谐波幅值比例一致,且电流谐波初始相位角满足β5=±π-ψ1+α5、β7=±π-ψ1+α7的条件。则可得电流表达式应为:
105:将电流从三相静止坐标系经坐标变换转换到d-q坐标系下,可得
可以看到,本发明所提的方法需要在永磁同步电机控制器电流闭环中注入6倍频的谐波。对比图4中理想正弦电流波形和注入5、7次谐波后的电流波形以及各自产生的电磁转矩可以发现,注入5、7次谐波后,电流变为尖顶波,但电磁转矩波动仅为正弦电流波形下的42%,有效抑制的永磁同步电机的电磁转矩波动。
实施例2
101:图5为一台4极12槽单层绕组表贴式永磁同步电机有限元模型。电机额定频率为50Hz,额定转速1500r/min,额定电流有效值为3A,即I1=4.24A。表贴式永磁同步电机通常采用id=0控制方式进行控制,即ψ1=0。
102:通过有限元仿真可以获得该电机的三相反电动势波形,如图6(a)所示。进一步通过傅里叶变换获得反电动势各次谐波的频谱分布以及各次谐波幅值,如图6(b)所示。
103:由图6(b)可以看到反电动势波形中5、7次为其主要谐波,而且他谐波的幅值均小于基波幅值的2%,因此可忽略其他谐波,重点对5、7谐波产生的电磁转矩波动进行抑制。在考虑5、7次谐波初始相位角的基础上,该电机的三相反电动势可表示为:
从图6(b)中,可以得到从图3(b)中,可以得出,反电动势基波幅值E1为184.4,5次谐波分量对应的幅值E5为24.7,7次谐波分量对应的幅值E7为8.1。
104:根据本发明步骤2和步骤3中的分析,为了对5、7次反电动势谐波与电流基波产生的转矩波动进行抑制,需要满足电流所含谐波分量的次数与反电动势谐波次数相同,电流谐波幅值比例与反电动势谐波幅值比例一致,且电流谐波初始相位角满足β5=±π-ψ1+α5、β7=±π-ψ1+α7的条件。则可得三相电流表达式应为:
105:将电流从三相静止坐标系经坐标变换转换到d-q坐标系下,可得
可以看到,本发明所提的方法需要在永磁同步电机控制器电流闭环中注入6倍频的谐波。对比图7中理想正弦电流波形和注入5、7次谐波后的电流波形以及各自产生的电磁转矩可以发现,注入5、7次谐波后,电流变为平顶波,但电磁转矩波动仅为正弦电流波形下的47%,有效抑制的永磁同步电机的电磁转矩波动。
以上以本发明的实施例为中心,详细介绍了本方法的具体应用过程,所描述的计算流程或某些特征的具体体现,应当理解为本说明书仅仅是针对给出实施例的电机特性来描述本发明,实际上对于不同结构的表贴式永磁同步电机进行反电动势谐波分析和谐波电流注入时某些细节上会有所变化,这些变化应该属于本发明范围内。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.利用电流谐波抑制表贴式永磁同步电机转矩波动的方法,其特征在于,包括以下步骤:
1)对表贴式永磁同步电机反电动势波形进行傅里叶分解,获得反电动势谐波的分布频谱和各次谐波的幅值大小;
2)将非正弦电流表示成傅里叶级数形式,基于三相反电动势和电流傅里叶级数表达式求解电磁转矩模型;
3)根据电磁转矩谐波成分不同进行分类,求解各类电磁转矩波动最小时电流谐波满足的条件;根据电磁转矩谐波成分不同进行分类,如下:
将永磁同步电机非正弦的三相反电动势和电流代入电磁转矩模型中,可得:
Tem=T1+T2+T3+T4 (4)
其中,
即将电磁转矩分为T1,T2,T3,T4四类,T1是由电机基波反电动势和基波电流作用产生的恒定的电磁转矩;T2是由谐波反电动势和基波电流作用产生的电磁转矩分量;T3是基波反电动势与谐波电流作用产生的电磁转矩分量;T4是由谐波反电动势和谐波电流产生转矩分量;Tem为电磁转矩,E1为反电动势基波幅值;ω为电机电角速度;ωr为电机旋转机械角速度,t为时间;v为反电动势谐波次数,且v≠1;Ev为v次反电动势谐波分量对应的幅值;αv为v次反电动势谐波分量的初始相位角;kev为v次谐波反电动势向量旋转方向的系数,其值为1或-1,kev=1代表v次谐波反电动势向量旋转方向与基波旋转方向相同,kev=-1代表v次谐波反电动势向量旋转方向与基波旋转方向相反;I1为电流基波幅值;ψ1为电机内功率因数角;u为电流谐波次数,且u≠1;Iu为u次电流谐波分量对应的幅值;βu为u次电流谐波分量的初始相位角;kiu表示u次谐波电流相量旋转方向的系数,其值为1或-1,kiu=1代表u次谐波电流相量旋转方向与基波旋转方向相同,kiu=-1代表u次谐波电流相量旋转方向与基波旋转方向相反;
所述求解各类电磁转矩波动最小时电流谐波满足的条件,如下:
由于T2和T3特征相似,且单独不可能为0,因此将T2+T3作为一组,T4单独作为一组单独进行分析:
1)若令T2+T3=0,即
可知,当任意v=u且kev=kiu时,若
且βu=±π-ψ1+αv
即当电流所含谐波分量的次数与反电动势谐波次数相同,电流谐波幅值比例与反电动势谐波幅值比例一致,且任意次电流谐波初始相位角满足βu=±π-ψ1+αv时,T2与T3两项产生的转矩波动可完全消除;
2)若令T4=0且kevv≠kiuu,即
当kevv=kiuu时,T4中的各分量均不为0,无法自行消除,因此仅能通过对应谐波产生的转矩分量进行相互抵消,即令任意两次反电动势谐波v1、v2分量和相应两次电流谐波u1、u2分量引起的转矩波动之和为0,如下:
当任意u1=v1、v2=u2且kev1=kiu1、kev2=kiu2时,若满足条件
且
即当任意两次电流所含谐波分量的次数与相应两次反电动势谐波次数相同,该两次电流谐波幅值比例与反电动势谐波幅值比例一致,且该两次电流谐波初始相位角满足上述关系时,T4产生的转矩波动可被完全消除,此时电磁转矩的平均值将保持不变;
4)根据反电动势谐波特点判断电磁转矩波动的主要来源,确定电流谐波次数、幅值和相位;具体如下:
由于永磁同步电机三相绕组结构对称,在定子圆周上互相错开120°电角度分布,反电动势中不含有偶数次以及3的倍数次谐波,因此v=5、7、11、13……,即反电动势仅含有6k±1次谐波,其中k=1,2,3……,且6k-1次谐波反电动势向量的旋转方向与基波相反,kev=-1,而6k+1次谐波反电动势向量的旋转方向与基波相同,kev=1;同样,电流谐波次数也应为u=5、7、11、13……,且当u=6k-1时kiu=-1,当u=6k+1,时kiu=1;
由于各电磁转矩波动分量大小均与反电动势谐波幅值和电流谐波幅值的乘积成正比,而一般情况下,反电动势和电流基波幅值要远大于谐波幅值,且谐波中5、7次谐波幅值最大,因此由基波与5、7次谐波产生的电磁转矩波动分量要大于其余电磁转矩波动分量,所以判断反电动势中5、7次谐波和基波产生电磁转矩波动为主要来源,
在此基础上,给定三相绕组电流为:
其中,
I5=E5I1/E1;β5=π-ψ1+α;I7=E7I1/E1;β7=π-ψ1+α7;
5)将含有所述步骤4)确定的含有谐波的三相绕组电流通入表贴式永磁同步电机三相绕组中,抑制电机电磁转矩波动,具体为,
将式(12)进行坐标变换,得到d轴电流分量Id与q轴电流分量Iq,分别为:
然后,将式(13)给定到电机控制器中,再通过电流闭环控制使得电机定子电流中注入适当5次和7次谐波,实现谐波电流注入。
2.根据权利要求1所述的利用电流谐波抑制表贴式永磁同步电机转矩波动的方法,其特征在于,所述步骤1)中,对表贴式永磁同步电机反电动势波形进行傅里叶分解,得到:
其中,分别为表贴式永磁同步电机定子三相绕组对应的反电动势;E1为反电动势基波幅值;ω为电机电角速度;t为时间;v为反电动势谐波次数,且v≠1;Ev为v次反电动势谐波分量对应的幅值;αv为v次反电动势谐波分量的初始相位角;kev为v次谐波反电动势向量旋转方向的系数,其值为1或-1,kev=1代表v次谐波反电动势向量旋转方向与基波旋转方向相同,kev=-1代表v次谐波反电动势向量旋转方向与基波旋转方向相反。
3.根据权利要求1所述的利用电流谐波抑制表贴式永磁同步电机转矩波动的方法,其特征在于,所述步骤2)中,将非正弦电流表示成傅里叶级数形式如下:
其中,分别为三相绕组中的电流;I1为电流基波幅值;ψ1为电机内功率因数角;u为电流谐波次数,且u≠1;Iu为u次电流谐波分量对应的幅值;βu为u次电流谐波分量的初始相位角;kiu表示u次谐波电流相量旋转方向的系数,其值为1或-1,kiu=1代表u次谐波电流相量旋转方向与基波旋转方向相同,kiu=-1代表u次谐波电流相量旋转方向与基波旋转方向相反。
4.根据权利要求2或3所述的利用电流谐波抑制表贴式永磁同步电机转矩波动的方法,其特征在于,所述步骤2)中,电磁转矩模型为:
其中,Tem为电磁转矩,ωr为电机旋转机械角速度,ωr=2πf/p;f为绕组电压和电流频率;p为极对数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510493654.7A CN105071717B (zh) | 2015-08-13 | 2015-08-13 | 利用电流谐波抑制表贴式永磁同步电机转矩波动的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510493654.7A CN105071717B (zh) | 2015-08-13 | 2015-08-13 | 利用电流谐波抑制表贴式永磁同步电机转矩波动的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105071717A CN105071717A (zh) | 2015-11-18 |
CN105071717B true CN105071717B (zh) | 2018-03-09 |
Family
ID=54501022
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510493654.7A Expired - Fee Related CN105071717B (zh) | 2015-08-13 | 2015-08-13 | 利用电流谐波抑制表贴式永磁同步电机转矩波动的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105071717B (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105406775A (zh) * | 2015-12-02 | 2016-03-16 | 江苏科技大学 | 一种永磁同步电机电磁转矩的辨识方法 |
CN105450099B (zh) * | 2015-12-28 | 2017-11-24 | 浙江兆丰机电股份有限公司 | 一种电动汽车的轮毂电机电磁减震装置 |
DE102016211406A1 (de) * | 2016-06-24 | 2017-12-28 | Robert Bosch Gmbh | Verfahren zur Regelung einer Synchronmaschine und Regelvorrichtung für eine Synchronmaschine |
CN107664746B (zh) * | 2016-07-29 | 2020-11-20 | 杭州三花研究院有限公司 | 步进电机堵转检测方法 |
CN106385206B (zh) * | 2016-11-30 | 2019-01-18 | 上海卫星工程研究所 | 高精度伺服控制系统纹波力矩辨识及抑制方法 |
CN107517030B (zh) * | 2017-06-30 | 2020-06-30 | 广东美的暖通设备有限公司 | 电机控制系统及其电机振动抑制方法和装置以及空调器 |
CN107425779B (zh) * | 2017-07-24 | 2020-07-03 | 广东美的暖通设备有限公司 | 电机的控制方法、装置、空调器及存储介质 |
CN107241048B (zh) * | 2017-07-24 | 2020-07-03 | 广东美的暖通设备有限公司 | 电机的控制方法、装置、空调器及存储介质 |
CN107425778B (zh) * | 2017-07-24 | 2020-07-03 | 广东美的暖通设备有限公司 | 电机的控制方法、装置、空调器及存储介质 |
CN107453675B (zh) * | 2017-07-24 | 2020-07-03 | 广东美的暖通设备有限公司 | 电机的控制方法、装置、空调器及存储介质 |
CN111256281B (zh) * | 2018-11-30 | 2021-10-22 | 广东美的制冷设备有限公司 | 运行控制方法及系统、压缩机和空调器 |
CN111525713B (zh) * | 2020-04-22 | 2021-12-28 | 东南大学 | 一种集中绕组外转子磁场调制电机的转矩脉动削弱方法 |
CN111740653A (zh) * | 2020-06-18 | 2020-10-02 | 浙江理工大学 | 一种表贴式永磁同步电机转矩波动系数计算方法及装置 |
CN114070157B (zh) * | 2020-08-04 | 2023-06-16 | 美的威灵电机技术(上海)有限公司 | 基于感应电动势的电机的控制方法、电机和存储介质 |
CN113612423B (zh) * | 2021-08-02 | 2024-03-01 | 上海数明半导体有限公司 | 基于反电动势的最大转矩电流比控制方法、装置 |
CN114583702B (zh) * | 2022-05-05 | 2022-08-05 | 湖北工业大学 | Lcl型有源电力滤波器的并联电流控制系统及控制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101515780A (zh) * | 2009-04-03 | 2009-08-26 | 东南大学 | 一种注入电流谐波补偿永磁电机定位力矩的控制方法 |
CN101860300A (zh) * | 2010-06-09 | 2010-10-13 | 东南大学 | 基于空间矢量调制的永磁电机转矩脉动抑制方法 |
CN102201770A (zh) * | 2011-05-30 | 2011-09-28 | 重庆大学 | 一种注入谐波电压抑制永磁同步电机谐波电流的控制方法 |
-
2015
- 2015-08-13 CN CN201510493654.7A patent/CN105071717B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101515780A (zh) * | 2009-04-03 | 2009-08-26 | 东南大学 | 一种注入电流谐波补偿永磁电机定位力矩的控制方法 |
CN101860300A (zh) * | 2010-06-09 | 2010-10-13 | 东南大学 | 基于空间矢量调制的永磁电机转矩脉动抑制方法 |
CN102201770A (zh) * | 2011-05-30 | 2011-09-28 | 重庆大学 | 一种注入谐波电压抑制永磁同步电机谐波电流的控制方法 |
Non-Patent Citations (5)
Title |
---|
Design of currents to reduce torque ripple in brushless permanent magnet motors;Hung JY等;《IEE Proceedings B - Electric Power Applications》;19930731;第140卷(第4期);全文 * |
Minimization of Torque Ripple in Brushless DC Motor Drives;HOANG LE-HUY等;《 IEEE Transactions on Industry Applications》;19860731;第IA-22卷(第4期);全文 * |
Torque ripple suppression of surfacemounted permanent magnet synchronous motor using harmonic injected currents;M. Rahimi等;《Power Electronic & Drive Systems & Technologies Conference (PEDSTC), 2010 1st》;20100524;全文 * |
一种抑制非理想反电势的BLDC电机转矩脉动的方法;肖伟等;《电工电能新技术》;20070123;第26卷(第1期);第51-54页 * |
用谐波注入抑制永磁同步电机转矩脉动;廖勇等;《中国电机工程学报》;20110725;第31卷(第21期);第119-127页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105071717A (zh) | 2015-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105071717B (zh) | 利用电流谐波抑制表贴式永磁同步电机转矩波动的方法 | |
Miyama et al. | Vibration reduction by applying carrier phase-shift PWM on dual three-phase winding permanent magnet synchronous motor | |
Niguchi et al. | Cogging torque analysis of magnetic gear | |
US9391546B2 (en) | Direct torque control method for inhibiting torque ripples | |
Parsa et al. | Interior permanent magnet motors with reduced torque pulsation | |
Zhou et al. | Influence of magnet shape on the cogging torque of a surface-mounted permanent magnet motor | |
Xia et al. | A novel direct torque and flux control method of matrix converter-fed PMSM drives | |
Chabchoub et al. | PMSM cogging torque reduction: Comparison between different shapes of magnet | |
Ramezani et al. | The modeling and position-sensorless estimation technique for a nine-phase interior permanent-magnet machine using high-frequency injections | |
Brahim et al. | Cogging torque minimization of surface-mounted permanent magnet synchronous machines using hybrid magnet shapes | |
Hua et al. | Improved model‐predictive‐flux‐control strategy for three‐phase four‐switch inverter‐fed flux‐reversal permanent magnet machine drives | |
Yan et al. | The torque ripple reduction in PMAREL machine using time-space harmonics analysis of air-gap flux density | |
Scuiller et al. | Design of a bi-harmonic 7-phase PM machine with tooth-concentrated winding | |
Xuelei et al. | Research on field-weakening control of multiphase permanent magnet synchronous motor | |
Fangyang et al. | Study on flux-weakening control for PMSM | |
Scuiller | Design of a 7-phase surface-mounted pm machine with tooth-concentrated winding | |
Wang et al. | Analysis of harmonic current in permanent magnet synchronous motor and its effect on motor torque | |
Dai et al. | The control strategy of open-winding permanent magnet synchronous motor drive system based on five-leg inverter | |
Yang et al. | Research on Winding MMF Harmonics of the PMSM with the Same Number of Poles and Slots | |
Zhang et al. | Comparative analysis of electromagnetic force inverter fed pmsm drive using field oriented control (foc) and direct torque control (dtc) | |
Junaid et al. | Improved model of the iron loss for the permanent magnet synchronous motors | |
Han et al. | Two-segment three-phase PMSM drive with carrier phase-shift PWM | |
Lu et al. | Compensation for harmonic flux and current of permanent magnet synchronous motor by harmonic voltage | |
Faiz et al. | Torque ripple and switching frequency reduction of interior permanent magnet brushless direct current motors using a novel control technique | |
Zhang et al. | Research on a sensorless SVM-DTC strategy for induction motors based on modified stator model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180309 |
|
CF01 | Termination of patent right due to non-payment of annual fee |