辅助槽对内置式永磁同步电机齿槽转矩的影响

陶彩霞，付名禄，高锋阳，万应利，刘汗青

(兰州交通大学 自动化与电气工程学院，兰州 730070)

永磁同步电机(permanent magnet synchronous motor, PMSM)因其微型化、功率因数高等优点，广泛应用于电动汽车、自动化、工业制造等行业。而电机齿槽转矩会影响电机性能，造成转矩波动，产生振动和噪声等[1]。因此，抑制齿槽转矩方法的研究具有重要的意义。

目前已有大量研究对永磁电机的齿槽转矩抑制方式开展了多层面的分析。Lin等[2]、Chen等[3]和鲍晓华等[4]采用规格不同的永磁体进行重新组合抑制齿槽转矩，但是需要对永磁体进行重新设计及重组，工艺复杂。郭有权等[5]、Ge等[6]在不同电机仿真中采用磁极偏移法削弱齿槽转矩，但是会引入新的谐波。李天元等[7]、Kang等[8]、杨玉波等[9]和周建等[10]用有限元法研究在内置式永磁同步电机转子上开规格不同的辅助槽对齿槽转矩产生的影响。唐旭等[11-12]归纳出异步起动永磁同步电机更加详细的定子参数分析公式，获得可以抑制齿槽转矩的参数取值原理,但未用有限元法做仿真分析。Ishikawa等[13]和Ren等[14]运用麦克斯韦张量法分析表贴式永磁同步电机的齿槽转矩。杨金歌等[15]和Bianchi等[16]分析开辅助槽的电机模型，研究了辅助槽求解子域，但是未全面分析辅助槽的各个参数。王轶楠等[17]、夏加宽等[18]和张登明等[19]用有限元法研究定子齿开辅助槽对表贴式永磁电机造成的影响，但辅助槽参数数据较少，取值有瑕疵，且分析开槽面积时未设置统一量。何庆领等[20]研究辅助槽对齿槽转矩的影响，但未分析开槽对电机其他性能如反电动势和气隙磁密的影响。现有研究报道大部分分析辅助槽对表贴式永磁电机齿槽转矩的影响，对内置式永磁电机齿槽转矩的研究较少。笔者采用定子齿开槽的方法，采用能量法解析有限元法分析内置式永磁电机。

笔者在文献[15-20]的基础上，以内置式永磁同步电机为样机，在开辅助槽的齿槽转矩表达式基础上，详细地解析了开槽数量、槽深和槽宽对齿槽转矩的影响，并分析了开辅助槽前后空载反电动势、气隙磁密谐波分析、波形畸变率等电机性能的变化；采用槽深槽宽分别为变量的有限元新方法，得到了齿槽转矩抑制效果最好的辅助槽参数，并首次提出以开槽面积为影响因素，在比较开槽形状时将其设置为统一量进行分析，分析结果更具有可比性。研究表明，合理的设计辅助槽参数能有效抑制永磁电机的齿槽转矩。

1 齿槽转矩的数学模型

齿槽转矩是电机不通电时的磁场能量W对定子与转子之间相对位置角α的负导数[21]，即

(1)

电机内存储的磁场能量近似为电机气隙和永磁体中的磁场能量：

(2)

式中：B(θ,α)为气隙磁密的分布函数；θ表示沿电机旋转方向变化的角度；μ0为空气和永磁体的磁导率；V为包括永磁体和气隙的积分区域；Wairgap+PM为气隙和永磁体中的磁场能量。

(3)

式中：Br(θ)为永磁体剩磁；hm(θ)为永磁体充磁方向长度沿圆周方向的分布；δ(θ,α)为有效气隙长度。

将式(3)代入(2)中，得

(4)

(5)

(6)

将式(4)～(6)代入式(1)计算，得到齿槽转矩表达式：

(7)

2 电机有限元分析

2.1 有限元模型

本研究中的分析样机为1台4极36槽的内置式异步起动永磁同步电机，电机相关参数如表1所示。使用Ansoft Maxwell软件建立电机模型，输入电机参数，赋材料属性，设置参数和自适应添加边界，再将其变换为Maxwell 2D模型，进行二维场仿真分析，得出有关参数。电机建模流程图如图1所示，电机有限元模型见图2。

表1 电机相关参数

图1 电机建模流程图Fig. 1 Motor modeling flow chart

图2 电机有限元模型Fig. 2 Motor finite element model

2.2 结果分析

图3 齿槽转矩波形图Fig. 3 Cogging torque waveform

3 辅助槽参数选择

3.1 辅助槽槽数对齿槽转矩的影响

样机为内置式，对其谐波次数进行分析。极数与槽数不同的配合形式可以使齿槽转矩的谐波次数不同，用fpn表示齿槽转矩的谐波次数：

(8)

式中NL为定子槽数z和极数2p的最小公倍数。

齿槽转矩由各次谐波组成，并且谐波的幅值同次数成反比，如果齿槽转矩的fpn上升，则幅值下降。因此，采用开辅助槽的措施改变NL(z,2p)以达到抑制齿槽转矩的作用。当开k个辅助槽时，则槽数z改变为(k+1)z，mNL(z,2p)增加，如果NL[(k+1)z,2p]/NL(z,2p)≠1，则齿槽转矩的fpn增高，齿槽转矩的幅值降低。本研究中的样机极数2p=4，槽数z=36，根据式(8)分析，fpn= 9m，开1个辅助槽时，fpn=18m；开2个辅助槽时，fpn=27m，辅助槽的数量对齿槽转矩的次数产生影响。当fpn增加，磁势和磁导谐波幅值降低，从而抑制齿槽转矩。然而过多的辅助槽数量导致气隙磁密降低，并且受电机结构和工艺复杂的制约，将矩形辅助槽的个数确定为2个。综上分析辅助槽的数量可以对齿槽转矩造成较大影响。

图4 辅助槽开槽的数量和位置Fig. 4 Locations and numbers of slots

图5 槽个数对齿槽转矩的影响Fig. 5 Effects of the number of slots on cogging torque

3.2 辅助槽槽深对齿槽转矩的影响

从改变Gn的原理分析槽深对齿槽转矩的影响，

(9)

式中：b为定子齿深；h为辅助槽深。

由式(9)可得，槽深变化将改变Gn，引起齿槽转矩发生变化；选取的h需符合要求，h过小抑制效果差，h会改变定子齿原有的磁路分布。

图6 不同槽深的齿槽转矩 Fig. 6 Cogging torque with different groove depth

图7 齿槽转矩随槽深的变化趋势Fig. 7 Cogging torque changing with groove depth

3.3 辅助槽槽宽对齿槽转矩的影响

从改变Gn的原理分析槽宽对齿槽转矩的影响：

(10)

式中：a为定子齿宽；l为辅助槽宽。

图8 不同槽宽的齿槽转矩Fig. 8 Cogging torque with different groove width

图9 齿槽转矩随槽宽的变化趋势Fig. 9 Cogging torque varies with groove width

3.4 辅助槽槽深槽宽同为变量对齿槽转矩的影响

表2 槽深槽宽同为变量的齿槽转矩

续表2

开辅助槽削弱齿槽转矩的原理实质上是齿槽转矩的基波周期发生变化，开槽后的齿槽转矩抵消了部分开槽前的齿槽转矩，因此总齿槽转矩减小。齿槽转矩的周期数并非遵循线性关系，而是服从傅里叶分解的c1sin(ωt+φ)关系，因此齿槽转矩的最优解并不完全满足槽深槽宽一个为定量，另外一个为变量的最优解。

3.5 辅助槽开槽形状对齿槽转矩的影响

图10 辅助槽形状Fig. 10 Different shapes of the auxiliary groove

图11 不同槽型齿槽转矩Fig. 11 Cogging torque with different groove shapes

4 开辅助槽后模型的分析

在定子齿开设矩形辅助槽，槽深为0.7 mm，槽宽为1.2 mm，面积为0.84 mm2，完成开槽前后齿槽转矩、气隙磁密、谐波分析、空载反电动势电机性能对比分析如图12～15所示。

图12 齿槽转矩对比Fig. 12 Cogging torque comparison

图13 气隙磁密(B)对比Fig. 13 Air gap magnetic density contrast

气隙磁密是电机重要的物理量，在定子齿开辅助槽会对电机的气隙磁密产生影响，因此需要研究开辅助槽后气隙磁密是否引入新的谐波。电磁力由切向力和径向力两部分组成，其中切向力致使转子运行,径向力则造成电机定转子形变震动，这会对电机的耐久度产生负面影响,并且会产生噪声，因此需要分析电机的径向磁密，研究气隙磁密中的谐波含量。因为电机磁路构造的特征，气隙磁密含有大量的谐波不能达到理想的正弦波分布，这会使电机在运行过程中产生杂声与抖动，导致较大的转矩波动，可靠性降低[22]。

为了研究电机气隙磁密，可以对电机在静磁场完成分析，当电机在额定转速1 500 r/min运行时，设置气隙中心位置R=(Rm+Rs)/2的气隙磁密，其中Rm为磁钢外径和Rs分别为定子内径。气隙磁密中径向磁密的计算公式[23]如下：

B=Bxcosθ+Bysinθ。

(11)

式中：B为径向磁密；Bx为磁密x轴分量；By为磁密y轴分量；θ为柱坐标中的相应角度。

气隙磁密波形畸变率用谐波畸变率RTHD表示[24]，RTHD越小则谐波含量越少，波形越近似正弦波。

(12)

式中：Bn为第n次气隙磁密谐波幅值，B1为气隙磁密基波幅值。开槽前后气隙磁密波形谐波分析结果如表3～4所示。

图14 各次谐波对比 Fig. 14 Harmonic comparison

图15 空载反电动势(E)对比Fig. 15 No-load back electromotive force comparison

表3 开槽前气隙磁密波形谐波分析结果

表4 开槽后气隙磁密波形谐波分析结果

在Matlab中求气隙磁密波形畸变率，开槽前的基波幅值为0.879 1 T，气隙磁密波形畸变率为21.23%；开槽后基波幅值为0.832 6 T，气隙磁密波形畸变率为17.78%。磁密幅值略微降低，不会对电机出力产生较大影响。谐波分析表明5、7次谐波削弱效果明显，因此开辅助槽在确保电机气隙磁密基波幅值不影响电机出力的情况下，有效地抑制了磁通密度中谐波含量，使空载气隙磁密波形更近似正弦波，抑制了齿槽转矩。

5 结 论

笔者采用有限元分析法，分析了在永磁同步电机的定子齿冠开辅助槽时，辅助槽数量、槽深、槽宽、开槽面积和开槽形状对齿槽转矩抑制的影响，并根据仿真结果分析了开槽前后反电动势、气隙磁密、各次谐波与畸变率的变化。

1)开槽数量会影响齿槽转矩，开2个辅助槽比开1个辅助槽抑制效果更明显。

2)槽深会影响齿槽转矩，齿槽转矩随着开槽深度增加先减小至最优值，然后增大；槽宽比槽深对齿槽转矩的影响更明显，齿槽转矩随着开槽宽度增加先减小至临界值，然后增大。

5)当槽宽与开槽面积相同时，矩形辅助槽比三角形和半圆形对齿槽转矩的抑制效果更好。

6)开辅助槽前后空载反电动势变化不大，气隙磁密基波幅值降低，5、7次谐波减小程度显著，气隙磁密波形畸变率下降。开辅助槽抑制齿槽转矩的同时使气隙磁密波形更接近于正弦波。