1、控制与应用技术EMCA2023，50(7)收稿日期:2023-03-14;收到修改稿日期:2023-04-03作者简介:郭培遥(1999)，男，硕士研究生，研究方向为永磁同步电机及其控制技术。梁建伟(1979)，男，硕士，副教授，研究方向为永磁同步电机及其控制技术。(通信作者)辐条式磁齿轮电机齿槽转矩抑制郭培遥，梁建伟，王昕华，尹利英(江西理工大学 电气工程与自动化学院，江西 赣州341000)摘要:磁齿轮电机具有转矩密度高、转子结构坚固的特点。针对外转子辐条式磁齿轮电机齿槽转矩高等问题，以一台三相 18 槽 28 极外转子辐条式磁齿轮电机作为研究对象，利用磁场调制理论对气隙进行解析计算，然后
2、基于麦克斯韦张量法分析外转子辐条式磁齿轮电机齿槽转矩产生的机理，并通过调节槽口角度和转子削角对磁齿轮结构进行优化。有限元分析结果表明，优化后电机的反电动势三次谐波被削弱了 53%、齿槽转矩幅值被削弱了 934%、输出转矩脉动下降了 737%。有限元分析结果验证了调节槽口角度和转子削角的结构优化方法的有效性。关键词:磁齿轮电机;磁场调制;麦克斯韦张量法;齿槽转矩中图分类号:TM351文献标志码:A文章编号:1673-6540(2023)07-0066-08doi:1012177/emca2023074Cogging Torque Suppression of Spoke Magnetic Gea
3、r MotorGUO Peiyao，LIANG Jianwei，WANG Xinhua，YIN Liying(School of Electrical Engineering and Automation，Jiangxi University of Science and Technology，Ganzhou 341000，China)Abstract:Magnetic gear motor has the characteristics of high torque density and strong rotor structure To solvethe problem of high
4、cogging torque of external rotor spoke magnetic gear motor，a three-phase 18-slot 28-pole externalrotor spoke magnetic gear motor is taken as the research object The air gap is calculated analytically using themagnetic field modulation theory，and then the mechanism of cogging torque of the external r
5、otor spoke magnetic gearmotor is analyzed based on Maxwell tensor method The magnetic year structure is optimized by adjusting the grooveangle and rotor cutting angle The finite element analysis results show that the third harmonic of the oplimized motorsback electromotive force is reduced by 53%，th
6、e amplitude of cogging torque is reduced by 934%，and the outputtorque ripple is reduced by 737%The finite element analysis results verify the effectiveness of the structuraloptimization method by adjusting the groove angle and rotor cutting angleKey words:magnetic gear motor;magnetic field modulatio
7、n;Maxwell tensor method;cogging torque0引言随着电气自动化和稀土行业的发展1，高功率密度和宽调速范围的永磁同步电机在电动汽车、航空航天、风力发电等领域都有着优异的表现2。其中，研制的基于磁场调制的外转子辐条式磁齿轮电机拥有更高的转矩密度，在低速、大扭矩等场合有明显的优势，因此优化电机输出转矩脉动、削弱齿槽转矩，对于提高电机性能具有十分重要的意义。针对磁场调制同步磁齿轮电机，一方面其基于磁场调制理论延伸出众多的平面拓扑结构3，其定子、转子、绕组和永磁体均有丰富的结构。磁齿轮电机最早由香港大学的 Chau 等4 结合磁齿662023，50(7)控制与应用技术E
8、MCA轮提出，其将一台高速永磁同步电机置于磁齿轮内部，将电动机和传动装置相结合。Wang 等5 将磁齿轮内部高速转子更改为相同极数的带绕组定子，相较于上述电机转矩提升了 20%5。QU等6 证实了磁齿轮电机和表贴式游标电机磁场耦合原理相同的特性，并对游标电机反电动势和输出转矩的产生机理进行阐述7。程明等89 基于单定子齿轮电机提出了一种切向充磁的外转子聚磁电机，其定子绕组采用低极数方案替代高转速游标电机内转子，相较于传统永磁同步电机其转矩提升了 80%。另一方面由于磁场调制电机聚磁效应导致齿槽转矩幅值大10，王培欣等11 采用余弦削极的转矩脉动抑制方法，解析了谐波产生的转矩脉动在输出转矩中的占
9、比，该方法有效地削弱了齿槽转矩。此外还有定子齿开槽12、磁极偏移13、偏心距磁极14、永磁体分段15、转子表面开槽等齿槽转矩削弱方法。文献 16通过对气隙磁场进行解析，采用定子齿肩削角的方案，其对齿槽转矩有明显的削弱效果17。本文基于磁场调制理论18，对一台 28 极 18槽的外转子辐条式磁齿轮电机气隙磁场进行解析。首先基于麦克斯韦张量法推导出电机产生齿槽转矩的气隙谐波极对数，然后采用优化定子开槽角和转子削角参数的方法优化气隙磁密，通过有限元二维仿真分析对比了优化前后电机的气隙空载磁场、转矩脉动和反电动势，对方案的可行性进行了验证。1电机磁场调制分析基于外转子辐条式磁齿轮电机磁场调制原理，分析
10、了其工作转矩和转矩脉动的形成机理和产生转矩脉动的主要谐波次数，推导出影响产生转矩的因素，为削弱转矩脉动和齿槽转矩的方法提供理论依据。11气隙磁场调制机理本文以一台 28 极 18 槽的外转子辐条式磁齿轮电机为例，定子和转子均采用硅钢片叠压，绕组为 8 极 18 槽的分数槽绕制，永磁体采用钕铁硼材料切向充磁来提高气隙的磁场强度，外转子辐条式磁齿轮电机主要参数如表 1 所示，电机结构如图 1 所示。表 1辐条式磁齿轮电机基本参数参数名称参数值转子内径/mm189定子外径/mm188转子外径/mm220永磁体宽度/mm35定子槽角/()5电机轴长/mm300定子槽数18永磁体极对数14绕组极对数4气
11、隙长度/mm05图 128 极 18 槽外转子辐条式磁齿轮电机结构由麦克斯韦张量法可以得到转矩产生的机理，外转子辐条式磁齿轮电机的谐波转矩为同极数且不同转速谐波相对耦合所产生，其表达式为T=l0(o i)SrBnBtandS(2)式中:l 为电机轴向长度;o为外转子内径;i为定子外径;r 为力臂;Bn、Btan分别为气隙曲面S 单元上径向和切向磁通密度矢量;0为真空磁导率。图 2永磁体作用下等效气隙磁动势波形当定子不开槽时，其永磁体在气隙中产生的磁感应强度分布可以等效为如图 2 所示的方波。电机气隙和永磁体内的磁矢位分别满足拉普拉斯方程和泊松方程，将气隙磁感应强度进行离散傅76控制与应用技术E
12、MCA2023，50(7)里叶变换，得到气隙磁感应强度随转子转动时，不同半径、不同 下气隙位置观测点处的径向分量Bn(r，)和切向分量 Btan(r，)的表达式为19 Bn(r，)=i=135B(r)n(i)cos ip(t)+ip0(3)Btan(r，)=i=135B(r)tan(i)sin ip(t)+ip0(4)式中:为转子旋转的转速;0为转子初始位置角;p 为转子极对数;B(r)j、B(r)i为气隙磁动势谐波傅里叶系数，且只与计算点的气隙半径 r 有关，其表达式为B(r)n(i)=B004Kniicosi2()ip(ip)2 1ro()ip1io()ip+1+ir()ip+1(5)B(
13、r)tan(i)=B004Ktaniisini2()ip(ip)2 1ro()iP1io()ip+1+ir()ip+1(6)式中:为 12 m;Kni、Ktani为关于谐波次数的系数;B0为永磁体剩余磁化强度。当定子开槽，引入定子齿作为调磁块后，由于铁心磁导率(r2000)远大于空气磁导率，磁力线大多会流入调磁铁心中，气隙中的相对磁导率分布如图 3 所示。图 3简化气隙磁导波形从定性角度分析，引入调磁块的过程可以由无开槽的气隙磁密乘以一个相对磁导来表示20 B=B(r，)(7)式中:B为开槽后的气隙磁场分布函数;B 为开槽的气隙磁场分布函数;(r，)为气隙磁导函数。通过保角变换，将开槽后的相对
14、气隙磁导率转换为一个复数函数:=r+j(8)式中:r、j分别为气隙磁导径向、切向分量。槽间距的复磁导率分布的一个周期波形可以用傅里叶级数的形式表示，得到整个气隙的复磁导率函数，即:r=0+n=1rncos(mN)(9)=n=0nsin(mN)(10)式中:0为径向磁导0 阶分量;rn、n为谐波幅值。将式(5)、(6)、(9)、(10)代入式(7)中可得到引入调磁块后电机气隙磁密表达式，其中永磁体气隙磁密调制后径向和切向分量分别为Bn=Bn(r，)r+Btan(r，)=0i=1，3，5Bn(i)cos ip(t)+ip0+i=0q=1n=1Bn(iqn)cos(nN+qip)qiptnN+qip
15、()+qi0(11)Btan=Bn(r，)+Btanr=0i=1，3，5Btan(i)sin ip(t)+ip0+i=0q=1n=1Btan(iqn)sin(nN+qip)qiptnN+qip()+qi0(12)式中:m、n 为气隙磁导的谐波次数，i、j 为气隙磁密谐波次数。12齿槽转矩产生机理由式(11)和式(12)可知，经过定子铁心的调制作用，气隙磁场含有大量高频谐波，根据麦克斯韦张量法即式(2)可知，产生转矩脉动的必要条件为两个磁动势波形的周期数相同且转速不同，即应满足如下式所示条件:n1N+qi1p=n2N+qi2p(13)引入 Q 为 N、p 的最小公倍数移项可得如下表达式aQn1+
16、n2=N(14)862023，50(7)控制与应用技术EMCAaQi1+i2=p(15)式中:a 为常数。各次谐波对应的转速为iqn=qipnN+qip(16)将案例样机 28 极 18 槽参数代入上式，由于永磁体磁场谐波幅值和气隙磁导谐波幅值随着谐波次数的增加而减小，产生转矩脉动的气隙磁密谐波主要由 14 对、238 对、266 对极磁场谐波构成，气隙磁导谐波主要由 18 次、252 次谐波和平均分量组成。根据式(13)和式(16)分别对磁场调制后 238对和 266 对磁场谐波产生因素和转速进行分析，见表 2 和表 3。表 2238 对气隙谐波组成分量原磁场极对数原磁场轴线位置磁导谐波极对数调制后磁场极对数调制磁场轴线位置谐波转速142522381701723802380表 3266 对气隙谐波组成分量原磁场极对数原磁场轴线位置磁导谐波极对数调制后磁场极对数调制磁场轴线位置谐波转速142522661901926602660由表 2 和表 3 可知，14 对磁场谐波分别和252 次磁导谐波耦合出 238 对和 266 对气隙谐波，转子 14 对基波旋转一周期同调制后的 238对、26
