基于谐波分析的永磁电机齿槽转矩抑制

2019年第47卷第1期
D
设计分析esign and analysis
杜晓彬等 基于谐波分析的永磁电机齿槽转矩抑制
37
　收稿日期:2017-07-26
基于谐波分析的永磁电机齿槽转矩抑制
杜晓彬,黄开胜,黄　信
(广东工业大学,广州510006)
摘　要:齿槽转矩是永磁电机特有的问题之一,也是设计制造过程中必须考虑的关键问题㊂为抑制电机的齿槽转矩,采用解析法与有限元法相结合的方式,研究了一种开辅助槽抑制电机齿槽转矩的方法㊂运用傅里叶级数对电机的齿槽转矩进行谐波分析,得到齿槽转矩的基波以及各次谐波的大小;将开辅助槽电机的齿槽转矩简化看作辅助槽与磁钢产生的齿槽转矩与原有电机的齿槽转矩的叠加,利用解析法确定辅助槽的开槽位置,使得特定位置的辅助槽与磁钢产生的齿槽转矩的谐波分量与电机原有的齿槽转矩的特定的谐波分量互相抵消,从而抑制齿槽转矩的基波以及低次数谐波分量,使得齿槽转矩幅值下降,并给出了具体的辅助槽位置公式㊂结合有限元方法,采用与普通槽相等的槽口宽度,并对辅助槽的深度进行仿真分析,得到最优的辅助槽深度㊂仿真结果表明,在特定的位置上开特定个数的辅助槽,能对齿槽转矩的基波与低次数谐波起到明显的削弱作用,降低齿槽转矩幅值㊂
关键词:永磁电机;齿槽转矩;辅助槽;有限元法
中图分类号:TM351 文献标志码:A 文章编号:1004-7018(2019)01-0037-04
Reduction of Cogging Torque of Permanent Magnet Machines
Based on Harmonic Analysis
DU Xiao⁃bin ,HUANG Kai⁃sheng ,HUANG Xin
(Guangdong University of Technology,Guangzhou 510006,China)
Abstract :Cogging torque is one of the special problems of permanent magnet machines,which is also a key problem to
be considered in the process of design and manufacture.In order to reduce the cogging torque of motor,a method of open⁃ing the auxiliary slots was presented to reduce the cogging torque of the motor by means of the combination of analytical method and finite element method.Fourier series was used to analyze the cogging torque of the motor,the magnitude of fun⁃damental harmonic and each harmonic were obtained.The cogging torque of the motor was simplified as the superposition of the torque between the cogging torque for auxiliary slots and cogging torque for original motor.The slot position of the auxil⁃
iary slots was determined by the analytical method.The harmonic component of the cogging torque of auxiliary slots in the special position was counteracted with the specific harmonic component of the original cogging torque of the motor.Thus,
the fundamental and low harmonic components of the cogging torque were suppressed and the amplitude of cogging torque decreased.The location formula of auxiliary slots was also bined with the finite element method,the depth of the auxiliary slots was simulated and analyzed so that optimum depth of the auxiliary slots was obtained,which the auxiliary slot width was equal to the ordinary slots.The simulation results show that,the auxiliary slots with a certain number in the particular position can reduce the fundamental and low harmonic of the cogging torque and the amplitude of cogging torque.
Key words :permanent magnet machine,cogging torque,auxiliary slots,finite element method
0　引　言
随着高性能永磁材料的出现,以及现代电机控制技术的发展,永磁电机在国民经济各部门的运用越来越广泛㊂永磁电机采用磁钢产生机电能量转换所需要的磁场,具有结构简单㊁体积小㊁运行可靠等特点
[1-2]
㊂然而,由于电枢开槽,当定转子产生相对
运动时,磁钢与定子齿部的相互作用会引起电机内磁场储能的变化,从而产生齿槽转矩,导致输出转矩
波动,引起电机振动和噪声,并影响电机的控制精度,因此如何有效地削弱齿槽转矩一直以来都是专
家学者研究的热点之一㊂
削弱电机齿槽转矩的方法有多种,例如采用斜极或者斜槽的方法㊁采用不等槽口宽度的方法㊁采用磁极偏移的方法㊁对偏心距或者极弧系数进行优化的方
法等[3-7],其中,开辅助槽削弱电机的齿槽转矩是一种有效便捷的方法㊂传统的开辅助槽的方法一般采用均匀开槽,即辅助槽在定子齿冠上均匀分布,每个电枢齿冠开1~3个辅助槽;或者采用有限元仿真的方式对辅助槽的位置㊁深度㊁宽度进行单变量参数优化扫描分析,得出最优的辅助槽参数,进行开槽㊂传统的方法难以确定开槽个数和最优的开槽位置,且当采用有限元法对位置进行确定时,花费时间较长㊂
D设计分析
esign and analysis 2019年第47卷第1期
杜
晓
彬
等 基于谐波分析的永磁电机齿槽转矩抑制 38　
本文采用解析法与有限元分析相结合的方法确定辅助槽参数㊂通过对齿槽转矩进行傅里叶分析,得到齿槽转矩的各次谐波,利用解析法确定辅助槽的开槽位置和个数,使得特定位置和个数的辅助槽产生的齿槽转矩谐波分量能有效抵消原有电机齿槽转矩的基波和低次数谐波㊂辅助槽槽口宽度采用与普通槽槽口宽度一致,并利用有限元方法分析了最优的辅助槽深度,从而达到抑制齿槽转矩的目的㊂以一台12槽10极的永磁同步电机为例进行仿真分析,验证了本文方法的有效性和正确性㊂
1　电机的齿槽转矩分析及辅助槽开槽方法
根据电机学原理,当定转子产生相对运动时,由于定子齿与磁钢产生相互作用,磁钢极弧中部与定子齿之间磁导基本不变,而磁钢两侧与定子齿之间磁导变化较大,导致电机磁场储能变化,从而产生了齿槽转矩㊂齿槽转矩被定义:
T cog=-∂W∂α(1)式中:W为电机磁共能;α为定转子相对位置角㊂当电机不斜槽或者斜极时,可将齿槽转矩展开为傅里叶表达式[8-9]:
T cog=∑¥n=1T n sin(N s nα)(2)式中:T n为齿槽转矩的傅里叶系数;N s为电机旋转一周齿槽转矩的周期数,其值为电枢槽数与极数的最小公倍数,即N s=LCM(2p,Q)㊂
当电机定子齿上开辅助槽时,相当于改变了电机的极槽配合,同样也影响了N s的大小㊂由式(2)可知,齿槽转矩的周期取决于N s,随着齿槽转矩的周期数变大,则频率越高,振幅越小,选择合适的辅助槽数,能有效抑制齿槽转矩的幅值㊂
为了简化分析,在分析时认为辅助槽的影响与普通定子槽影响一样,且本文只讨论矩形槽㊂由于电机总的齿槽转矩可以认为是每个磁钢所对应的齿槽转矩的叠加[10],故本文将开辅助槽的电机分为2个次级结构进行分析㊂如图1㊁图2所示,以一台12槽10极电机在齿冠上开一个辅助槽为例,第一部分为电机定子槽部分的次级结构,第二部分为电机辅助槽部分的次级结构,两部分则可以合成一个完整的开槽电机的定子铁心㊂当电机开辅助槽时,电机总齿槽转矩为第一部分次级结构对应的齿槽转矩与第二部分次级结构对应的齿槽转矩的叠加㊂则第一㊁二部分次级结构对应的齿槽转矩:
T slot=∑¥n=1T n sin(N s nα)(3)
T aux=∑¥n=1T n sin[N s n(α+θ)](4)式中:θ为沿着电机旋转方向普通槽组与辅助槽组相差的位置角㊂当电机开辅助槽时,电机总的齿槽转矩:
T cog=T slot+T aux (5)
图1　
第一部分次级机构
图2　第二部分次级机构
为了使得辅助槽对应的齿槽转矩的一次谐波能对电机原有的齿槽转矩的一次谐波产生抵消作用,令θ取值:
θ1=(12+k)2πN
s
(6)式中:k为非负整数,且由于θ1取值不能超过一个齿距,故k≤N s/Q㊂则当n为1时,将式(6)代入式(5),可得到T cog一次谐波分量为0,电机原有齿槽转矩的一次谐波被完全抵消,只剩下二次及以上的谐波㊂为了抑制n次谐波,可选在不同的位置开一组辅助槽,令θ取值:
θn=(12+k)2πnN
s
(7)式中:k为小于等于nN s/Q的非负整数㊂则将式(7)代入式(5),可得到T cog的n次谐波分量为0,电机原有齿槽转矩的n次谐波被完全抵消㊂
由于不均匀的槽口宽度以及不均匀的槽分布会引进新的谐波[11]㊂对于式(6)㊁式(7),k取值有多个,即抑制n次谐波辅助槽的位置可以有多种选择,应该按照如下原则选取:同一个齿冠上开的辅助槽,应使得辅助槽开口依齿冠中心线呈对称;对于电机上所有的辅助槽与普通槽,应使得辅助槽槽口与普通槽开口在铁心表面尽量均匀分布㊂而对于辅助槽槽口宽,应选取与普通槽槽口宽一致㊂为了简化分析,分析时将辅助槽与普通槽的影响等效,而实际上,辅助槽的影响并不等同于普通槽的影响,故可以通过有限元方法对辅助槽深度进行参数优化扫描分析,选取能抑制原有齿槽转矩的基波以及n次谐波的最优辅助槽深度㊂
2019年第47卷第1期
D
设计分析esign and analysis
杜晓彬等 基于谐波分析的永磁电机齿槽转矩抑制
39
　2　有限元仿真分析
为了进一步分析和验证本文方法的正确性与有效性,本文以一台12槽10极的电机为例进行有限元仿真分析,通过在合适的位置上开辅助槽,抑制电机原有齿槽转矩的基波以及二次谐波,从而抑制齿槽转矩幅值㊂
2.1　永磁同步电机齿槽转矩分析
本文采用的12槽10极电机基本参数如表1所示,采用ANSYS Maxwell 2D 软件建立电机模型,并对电机原有的齿槽转矩进行仿真分析,由图3可知,齿槽转矩的幅值为52.7844mN㊃m㊂对齿槽转矩进行傅里叶分析,得到各次谐波幅值分布如图4所示,可以看出,较大的谐波分量主要集中在基波和二次谐波,第三至六次分量很小,其中,基波分量幅值为49.1862mN㊃m㊂
表1　电机模型基本参数
参数
数值参数
数值
额定功率P e /
W
170定子外径Φso /mm 88.8转子外径Φro /mm 107定子内径Φsi /mm 33转子内径Φri /mm 90
铁心长度l Fe /mm 45槽口宽b /mm 2.5槽口深h /mm 2槽肩深h′/mm 3极弧系数αp
0.8磁钢厚度δ/mm
5.6图3　优化前电机齿槽转矩波形
图4　优化前电机齿槽转矩各次谐波分量
2.2　辅助槽位置分析
采用本文的抑制方法进行分析,对于12槽10极的分数槽电机,槽数与极数的最小公倍数N s 为
60,即定子旋转一周产生的齿槽转矩周期数为60,且每个定子齿对应的机械角度为2π/Q =30°㊂对于基波,采用式(6)进行计算,且注意θ取值不能超过一个齿宽,则θ取值为3°,6°,15°,21°,27°,即每一组的辅助槽开槽位置应与普通槽偏移以上计算度数㊂以一个齿冠为例子,如图5所示,可以看出,开槽的位置是关于齿槽中心线呈现对称分布的,且15°对应的辅助槽位置刚好位于齿冠中心线上㊂图5　抑制基波的辅助槽开槽位置
采用辅助槽槽口宽与定子槽口宽度一致,槽深为1mm,由于在3°,27°时辅助槽与定子槽重合,仿
真分析在6°,15°,21°时3个不同位置开辅助槽时电机总齿槽转矩,即每一组的辅助槽开槽位置与普通槽偏移以上度数,如图6所示㊂可以看出,3个波形基本一致,说明开槽效果基本相同,只有细微不同㊂其中,15°位置的总齿槽转矩最低,为15.8702mN㊃m㊂
对波形进行傅里叶分析,得到各次谐波幅值分布如
图7所示㊂
可以看出,电机原齿槽转矩基波被明显抑制,然而,二次谐波明显升高,这是因为辅助槽对应齿槽转矩的二次谐波对电机原齿槽转矩二次谐波
起到叠加作用,使其增大㊂对比不同开槽位置下的谐波分析结果,发现15°位置的总齿槽转矩前三次谐波幅值皆比其他2个位置开槽时的总齿槽转矩前三次谐波幅值低,其中,15°位置基波分量辅助为
1.6292mN㊃m,比未开槽前有明显的下降,说明比均匀开槽以及依齿冠中心线开槽效果要好㊂故对于齿槽转矩的基波,当开一组槽时,可选取15°位置进行开槽;当开两组槽时,可选取6°和21°位置同时进行开槽,但是这种开槽方式使得普通槽与辅助槽在铁心表面分布不均匀,效果较差;当开三组槽时,可在6°,15°,21°位置同时进行开槽㊂
图6　不同位置开槽的电机总齿槽转矩波形
图7　不同位置开槽的总齿槽转矩谐波分布
2.3　齿槽转矩基波及二次谐波抑制分析
根据上述的分析结果,为了简单方便,选取开一组槽,且在15°位置开槽㊂为了进一步确定最优的
辅助槽深度,采用有限元软件参数优化分析的方法
D
设计分析esign and analysis 2019年第47卷第1期
杜晓彬等 基于谐波分析的永磁电机齿槽转矩抑制
40
　对辅助槽深度进行扫描,槽深取值范围为0.1~1.5
mm,得出最优的辅助槽深度为1.1mm,如图8所示㊂可以看出,电机总齿槽转矩为15.1257mN㊃m㊂对波形进行傅里叶分析,得到各次谐波幅值分布如图9所示,基波分量幅值为0.0319mN㊃m,电机原齿槽转矩基波分量基本被抵消,二次谐波分量为
14.9133mN㊃m㊂
图8　
抑制基波后电机总齿槽转矩波形
图9　抑制基波后电机总齿槽谐波分布
为了进一步抑制电机原齿槽转矩二次谐波,采用式(7)进行计算,n 取2,且注意θ取值不能超过一个齿宽,则θ取值为1.5°,4.5°,7.5°,10.5°,13.5°,16.5°,19.5°,22.5°,25.5°,28.5°,按照尽量均匀开槽和依齿冠中心线开槽的原则,选取7.5°和22.5°位置开槽,并使新的辅助槽与15°位置辅助槽
槽深一致,采用有限元分析参数优化方法对槽深重新进行扫描分析,取值范围为0.1~1.5mm,得出最优的辅助槽深度为1.2mm,电机总齿槽转矩取最低幅值为1.5198mN㊃m,比电机原齿槽转矩幅值下降了97.12%㊂图10和图11分别为优化前后电机齿槽转矩波形对比和傅里叶分析㊂从图10
可以看
图10　
优化前后电机齿槽转矩波形对比
图11　优化前后电机总齿槽转矩波形谐波分布对比
出,齿槽转矩幅值明显下降㊂对波形进行傅里叶分析,得到各次谐波幅值分布,如图11所示,基波分量以及二次谐波分量幅值分别为0.0113mN ㊃m,
1.0208mN㊃m,相较于电机原齿槽转矩,分别下降了99.97%,87.56%,前两次谐波得到明显抑制㊂3　结　语
本文对齿槽转矩进行了谐波分析,采用解析法与有限元法相结合的方式,研究了一种抑制齿槽转矩的方法:对电机原齿槽转矩进行傅里叶分析;利用辅助槽对应的齿槽转矩谐波分量抵消电机原齿槽转矩的基波和低次数谐波,并给出了抵消基波和n 次
谐波的开槽位置的具体公式和选取原则;辅助槽槽口宽采用普通槽槽口宽度,并利用有限元法分析了最优的辅助槽深度㊂以一台12槽10极电机为例进行仿真分析,结果表明,采用本文的方法能有效抑制齿槽转矩的基波和低次数谐波,使得齿槽转矩幅值明显下降㊂参考文献
[1]　谭建成.永磁无刷直流电机技术[M].北京:机械工业出版社,2011.
[2]　王秀和.永磁电机[M].2版.北京:中国电力出版社,2011.
[3]　石松宁,王大志,ShiSongning,等.消弱永磁驱动器齿槽转矩的斜极优化方法[J].电工技术学报,2015,30(22):27-32.[4]　杨玉波,王秀和,陈谢杰,等.基于不等槽口宽配合的永磁电动
机齿槽转矩削弱方法[J].电工技术学报,2005,20(3):40-44.[5]　杨玉波,王秀和,张鑫,等.磁极偏移削弱永磁电机齿槽转矩方法[J].电工技术学报,2006,21(10):22-25.
[6]　李延升,窦满峰,樊鑫.表贴式永磁电动机气隙磁场及齿槽转
矩解析计算[J].微特电机,2012,40(12):10-15.
[7]　王明星,王爱元,李轶华.一种优化齿槽转矩抑制永磁同步电
机振动和噪声的方法[J].电机与控制应用,2017,44(2):110-114.
[8]　BIANCHI N,BOLOGNANI S,CAPPELLO A D F.Reduction of
cogging force in PM linear motors by pole-shifting[J].IEE Pro⁃ceedings-Electric Power Applications,2005,152(3):703-709.
[9]　Zhu Z Q,Howe D.Influence of design parameters on cogging
torque in permanent magnet machines[J].IEEE Transactions on Energy Conversion,2000,15(4):407-412.
[10]　HUANG S.A method for reducing cogging torque by slot-opening
shift in permanent magnet motors[J].Transactions of China Elec⁃trotechnical Society,2013,28(3):99-106.
[11]　刘婷,欧阳红林,黄守道.基于重复单元削弱永磁风力发电机
齿槽转矩[J].电工技术学报,2011,26(12):43-48.
作者简介:杜晓彬(1991 ),男,硕士研究生,研究方向为新型特种电机设计及其控制㊂。