轴向磁通电机驱动装置用永磁电机磁场特性仿真研究发表时间:2025-12-03 16:07 抽象的 轮缘驱动装置(RDD)将电机和叶轮集成在一起,可实现涡轮机械系统的无轴、模块化和一体化运行,具有广阔的应用前景。为了减小RDD的轴向长度和径向厚度,本文设计了一种采用薄轭宽齿分数槽集中绕组定子和无铁芯哈尔巴赫永磁体阵列转子的电机。对其气隙磁场特性进行了理论分析和有限元仿真。结果表明,对于薄轭宽齿分数槽集中绕组永磁电机,磁极产生的谐波磁场主要应考虑磁极谐波磁动势与恒定气隙比磁导率相互作用产生的磁场分量,以及磁极基波磁动势与基波和二阶谐波气隙比磁导率相互作用产生的磁场分量。电流产生的谐波磁场主要考虑由电流产生的小极对数(NOPP)大幅谐波磁动势与恒定气隙比磁导率相互作用产生的磁场分量。与径向磁通密度相比,切向磁通密度具有相同的NOPP和频率分量,相位差为90°。二者之间的基波幅值差异较大,而谐波幅值差异较小。 关键词: 轮缘驱动装置;电机;分数槽集中绕组;哈尔巴赫永磁体阵列;磁场特性 1. 引言 轮缘驱动装置(RDD),也称为轮缘驱动螺旋桨(RDP)或轮缘驱动推进器(RDT),其中的“轮缘”通常指的是叶轮轮缘。叶轮是机械设备中广泛应用的能量转换部件,例如泵叶轮和螺旋桨叶轮[
图 1. 按应用分类,RDD 包括船舶推进螺旋桨和管道输送泵,分别如图 2a、2b 所示。
图 2. 按结构分类,RDD 包括带支撑轴的 RDD 和不带支撑轴的 RDD,分别如图 3a、3b 所示。
图 3. 以电机为原动机的无轴轮缘驱动推进器(RDD)最早由Kort于1940年在德国专利DE688114中提出[ RDD关键技术的研究主要包括导流槽和叶轮水动力性能的优化[ 本文提出了一种用于 RDD 的永磁电机设计,该电机结合了薄轭宽齿分数槽集中绕组定子和无铁芯 Halbach 永磁体阵列转子,并对其磁场特性进行了仿真研究。 2. RDD用永磁电机的设计 当RDD用作管道输送泵时,仅存在内部流场而无外部流场。然而,当其用作船舶推进螺旋桨时,则同时存在内部流场和外部流场。无论在何种应用场景下,RDD都应具有相对较小的径向厚度和轴向长度。 2.1. 薄轭宽齿分数槽集中绕组定子 普通电机的电枢绕组端为喇叭形,占据相对较大的轴向空间,如图4a所示。分数槽集中绕组电机的端绕组结构简单,电机的轴向长度相对较小[
图 4. 分数槽集中绕组末端较短,使其非常适合用于RDD永磁电机。在确定电机槽极组合方案时,目标是获得较大的基波绕组系数、较小的谐波磁动势分量、较小的谐波反电动势幅值以及较小的齿槽转矩。综合考虑上述因素以及叶轮外径(本文中叶轮外径为250 mm),最终选择了24槽20极的绕组方案。 三相绕组的布置如图5a所示。图中,线圈边 A 和 X 属于 A 相绕组,线圈边 B 和 Y 属于 B 相绕组,线圈边 C 和 Z 属于 C 相绕组。线圈边 X、Y 和 Z 中的电流大小相等,方向与线圈边 A、B 和 C 中的电流方向相反。由于 RDD 在水中工作,因此绕组用胶水密封,如图5b所示。
图 5. 图 6显示了 24 槽 20 极分数槽集中绕组的电枢谐波磁动势的单位值,其中未知节点类型:字体表示磁动势极对的数量,未知节点类型:字体未知节点类型:字体未知节点类型:字体未知节点类型:字体未知节点类型:字体未知节点类型:字体表示单位磁动势的值。
图 6. 2.2. 无芯哈尔巴赫永磁体阵列转子 普通永磁电机的磁极通常由极性相反、径向磁化的交替永磁体组成,如图7a所示。通过将切向磁化的永磁体(如图7b所示)嵌入到普通永磁电机的永磁体排列中,可以形成哈尔巴赫永磁体阵列[
图 7. RDD的永磁电机转子通常在水介质中运行。叶轮叶片和轮缘一般采用耐腐蚀不锈钢制成,为保证转子强度,不锈钢轮缘需要一定的厚度。不锈钢的磁导率远低于铁磁材料。同时,为了保证RDD的平稳运行,转子的厚度应尽可能小。因此,本文的转子未采用传统的铁磁材料,而是结合单侧磁聚集效应的哈尔巴赫永磁体阵列和不锈钢轮缘,形成了一种无铁芯哈尔巴赫永磁体阵列转子方案。每个永磁体均并联磁化。永磁体的排列方式如图8所示,箭头指示磁化方向。
图 8. 最终,RDD 的永磁电机的设计参数如下表 1所示。 表 1.
3. 基于一维磁路法的气隙磁场谐波分析 一维线性磁路法是电机设计和分析中最早采用的方法。它将电机气隙磁通密度视为气隙磁势与气隙比磁导率的乘积。虽然该方法对复杂的几何形状和材料非线性考虑不足,但其概念和逻辑清晰[ 3.1. 径向气隙磁通密度 在空载条件下,磁极磁动势可表示为[ 公式中,正值或负值表示电枢磁动势的旋转方向。正值表示电枢磁动势的旋转方向与转子相同,而负值表示电枢磁动势的旋转方向与转子相反。 考虑谐波电流时,气隙合成磁动势为 当忽略铁磁材料的磁阻时,气隙磁通密度可以表示为气隙磁动势与气隙比磁导率的乘积。需要强调的是,磁通密度本质上是一个矢量,而磁动势和比磁导率均为标量。磁动势和比磁导率在切向方向上没有解析表达式。磁动势乘以比磁导率的方法是一种一维磁场计算方法。由此得到的磁通密度仅为气隙磁通密度矢量的径向分量,而无法考虑气隙磁通密度矢量的切向分量,也无法考虑气隙径向方向上径向磁通密度幅值的变化。气隙磁通密度仅是圆周角位移和时间的函数,即 1 (19) 公式中,1 在空载条件下,将式 (14) 和 (16) 代入式 (19) 并进行三角函数乘法运算,即可得到气隙的径向磁通密度。 1 (20) 在空载条件下,忽略比磁导率的谐波分量,径向气隙磁通密度的来源、空间中的NOPP以及基频倍数(简称MOFEF)如表2所示。 表 2.
在负载条件下,仅考虑基波电流,将式 (14) 和 (17) 代入式 (19) 并进行三角函数的乘法和运算,即可得到气隙的径向磁通密度。 1 (21) 在考虑基波电流时,忽略了气隙比磁导率的谐波分量。与空载工况相比,径向气隙磁通密度的新来源 NOPP 和 MOFEF 如表 3所示。 表 3.
考虑谐波电流,将式 (14) 和 (18) 代入式 (19),并对三角函数进行乘法和运算,即可得到气隙的径向磁通密度。 1 (22) 在考虑谐波电流时,忽略了气隙比磁导率的谐波分量。与仅考虑基波电流相比,径向气隙磁通密度的新来源 NOPP 和 MOFEF 如表 4所示。 表 4.
3.2.切向气隙磁通密度 一维磁路法无法考虑切向气隙磁通密度。为了考虑切向气隙磁通密度,必须采用二维磁场法。许多学者对电机二维磁场的求解进行了研究。 对于磁极的磁场,参考文献[ ⎧⎩⎨2D先生( (23) 公式中,rm 在考虑槽效应时,参考文献[ ⎧⎩⎨2D先生( (24) 公式中,rm 对于电枢磁场,参考文献[ ⎧⎩⎨2Dar( (25) 公式中,rℎ和tℎ表示气隙径向和切向磁通密度各谐波分量的幅值。它们的值与磁场点的径向坐标分量、气隙谐波磁通密度的NOPP、谐波电流的MOFEF、定子铁芯内径、转子铁芯外径、绕组匝数、各谐波电流幅值、谐波节距系数和谐波分布系数相关。可以看出,在不考虑槽效应的情况下,对于表面贴装式永磁电机电枢气隙磁场,径向和切向磁通密度各谐波分量的幅值不同,但径向和切向磁通密度谐波分量的NOPP和MOFEF相同,且具有相同NOPP和MOFEF的径向和切向磁通密度谐波分量的相位角仍然相差90°。 考虑槽效应时,相应电枢气隙磁场的径向和切向分量的解析表达式为: ⎧⎩⎨2Dar( (26) 公式中,rℎ 在负载条件下,气隙的合成磁场是磁极磁场和电枢磁场的叠加。由以上分析可知,考虑电枢电流的各种谐波和槽效应时,径向和切向磁通密度各谐波分量的幅值不同,但径向和切向磁通密度各谐波分量的净正交极化率(NOPP)和磁极磁场强度因子(MOFEF)相同,且具有相同NOPP和MOFEF的径向和切向磁通密度谐波分量的相位角仍相差90°。基于此结论,可以对基于一维磁路法的气隙磁场解析模型进行补充,得到气隙磁通密度切向分量的表达式。 1 (27) 公式中的系数是与具有相同 NOPP 和 MOFEF 的径向气隙磁通密度谐波幅值成正比的切向气隙磁通密度谐波幅值的系数。 4. 基于二维有限元法的气隙磁场仿真分析 本文采用有限元法[ 4.1 空载 图9展示了空载工况下磁场线和磁场内磁通密度的分布情况。可以看出,尽管转子没有铁芯,但采用哈尔巴赫结构后,磁场主要集中在靠近转子气隙的一侧。转子内侧也存在一些磁线,但磁通密度值相对较小。哈尔巴赫永磁体阵列的单侧磁聚集效应十分明显。定子齿和轭的磁通密度约为1.6 T,因此定子铁芯工作在铁磁材料磁化曲线的“拐点”附近,充分利用了铁磁材料的磁导率而不发生过饱和。这有利于在保证电磁性能的同时提高电机的功率密度。空载工况下转子内侧以及定子齿和轭的磁通密度分布验证了本文电机设计的合理性。不考虑磁极极性时,由于槽数和极数的最大公约数为4,电机沿圆周方向的磁场分布重复4次。因此,可以看出电机的整体磁场在空间上呈4个周期分布。
图 9. 在空载条件下,径向气隙磁通密度的波形和谐波含量直接影响电机反电动势的波形和谐波含量。图10显示了无铁芯哈尔巴赫阵列永磁电机在特定时刻空载条件下的径向气隙磁通密度波形和频谱。
图 10. 可以看出,空载条件下哈尔巴赫阵列永磁电机径向气隙磁通密度的3次、5次和7次谐波含量相对较高,这是因为哈尔巴赫阵列每极仅有2个磁块,且各永磁块的磁化角度差为90°。如果持续增加每极磁块的数量并减小各永磁块的磁化角度差,气隙磁通密度将趋于正弦波形,但这会使制造工艺更加复杂。 虽然图10反映了电机气隙磁通密度的空间分布特征,但它无法反映气隙磁通密度随时间的变化。所得到的气隙磁通密度分布仅仅是某一时刻特定极槽相对位置处的气隙磁通密度分布。实际上,由于磁极与定子齿的相对位置随时间变化,气隙磁场分布也会随之变化。 为了全面分析电机实际工况下气隙磁通密度随时间和空间的变化,本文计算了一个电周期内不同转子位置的气隙磁场。所得气隙磁通密度径向分量和切向分量的时空分布如图11所示。
图 11. 从图11可以看出,在气隙圆周360°范围内,气隙磁通密度的径向分量和切向分量均有10个周期,反映了24槽20极电机的基本磁场。在0.006 s内,气隙磁通密度的径向分量和切向分量均经历了两个半周期的正负值,反映了电机的一个完整电周期。气隙磁通密度的最大径向分量接近1 T,而最大切向分量约为0.3 T,波形畸变较大。 气隙磁通密度径向和切向分量的时空分布图只能粗略确定气隙磁通密度径向分量的基频值,不利于分析气隙磁通密度径向和切向分量的谐波值,而气隙磁通密度的谐波分量是引起电机电磁振动的重要因素。对不同时刻磁通密度空间分布进行一维傅里叶变换得到的谐波幅值各不相同,且随时间变化。具有相同NOPP的空间谐波包含不同转速的分量。本文对气隙磁通密度径向和切向分量的时空分布进行二维傅里叶变换,得到的时空频谱如图12和图13所示。图中,将空间360°圆周内的磁通密度分布的一次波动视为空间中一对磁极的磁通密度。图 12b和图 13b分别是图 12a和图 13a的放大局部视图。
图 12.
图 13. 从图12a可以看出,空载条件下径向气隙磁通密度谐波主要与转子磁极(10对)有关,谐波磁通密度的NOPP为(2 通过比较图 12和图 13,可以发现空载时气隙磁通密度径向和切向分量的 NOPP 和频率分量完全相同,验证了前文对气隙径向和切向磁通密度关系的分析——即切向磁通密度谐波总是与径向磁通密度谐波同时出现,它们的相位差仅为 90°,但 NOPP 和 MOFEF 相同。此外,通过比较图 12a和图 13a,可以发现基波(NOPP 基于有限元法获得不同时刻的磁场分布后,即可通过相磁链随时间变化的差值商计算电机反电动势随时间的变化。为验证本文有限元法计算电机电磁场的有效性,对一台24槽20极轮缘驱动样机进行了实验。图14展示了实测反电动势与有限元法计算反电动势的对比。
图 14. 如图14所示,有限元法计算得到的反电动势与实测值总体吻合度较高,验证了本文所用有限元法的有效性。由于分数槽集中绕组电机的谐波绕组系数通常较低,因此其反电动势的谐波含量也相对较低。此外,实测反电动势频谱中奇次谐波的幅值低于有限元计算结果,且存在分数次谐波。这主要是由于电机实际制造过程中永磁体安装公差的累积,导致磁极偏移和转子磁极分布不完全对称,如图15所示。
图 15. 多极分数槽集中绕组表面贴装式永磁电机的电磁气隙相对较大。电机制造过程中,转子偏心和气隙不均匀的影响相对较小。然而,由于磁极块数量较多,容易累积安装公差,导致磁动势和反电动势的分数谐波产生。实际制造过程中应注意这一点。 4.2. 仅包含基本电流 在负载条件下,电枢电流会对气隙磁场产生影响。当仅考虑基波电流时,轮缘驱动装置中24槽20极分数槽集中绕组无铁芯哈尔巴赫阵列永磁电机的气隙磁通密度径向和切向分量的时空频谱如图16和图17所示。
图 16.
图 17. 从图 16a可以看出,考虑基波电流后,磁极产生的谐波磁通密度的变化与NOPP 4.3. 含有谐波电流 当考虑谐波电流时,轮缘驱动装置的 24 槽 20 极分数槽集中绕组无铁芯 Halbach 阵列永磁电机的气隙磁通密度径向分量和切向分量的时空频谱分别如图 18和图 19所示。
图 18.
图 19. 考虑谐波电流后,图 16a和图 18a的对比未观察到明显差异。然而,从图 16b和图 18b的对比可以看出,谐波电流并未引入新的谐波分量,而是增大了磁通密度的幅值。MOFEF 5. 结论 对于用于 RDD 的薄轭宽齿分数槽集中绕组永磁电机,可以得出以下结论: 对于空载条件下的谐波磁场,主要磁场谐波分量应考虑磁极谐波磁动势与恒定气隙比磁导率相互作用产生的分量,以及磁极基波磁动势与基波和二阶谐波气隙比磁导率相互作用产生的磁场分量。 对于电机基波和谐波电流产生的谐波磁场,主要应考虑由基波和谐波电流产生的小幅值谐波磁动势与恒定气隙比磁导率相互作用而产生的磁场分量。其他幅值较小的谐波磁场分量可以忽略不计。 切向磁通密度谐波总是与径向磁通密度谐波同时出现;它们具有相同的 NOPP 和 MOFEF,它们之间的相位差为 90°,它们的幅值差很小。 |