浙江大学学报(工学版), 2024, 58(10): 2182-2191 doi: 10.3785/j.issn.1008-973X.2024.10.022

电子与电气工程

新型盘式行星永磁电机的等效动态磁网络模型

梁志文,, 刘欣,

天津工业大学 天津市现代机电装备技术重点实验室,天津 300387

Equivalent dynamic magnetic network model of novel disc planetary permanent magnet motor

LIANG Zhiwen,, LIU Xin,

Tianjin Key Laboratory of Modern Electromechanical Equipment Technology, Tiangong University, Tianjin 300387, China

通讯作者: 刘欣,女,教授,博导. orcid.org/0000-0002-1269-8464. E-mail:liuxin@tiangong.edu.cn

收稿日期: 2023-06-29  

基金资助: 国家自然科学基金资助项目(51875408);天津市研究生科研创新项目(2021YJSB229).

Received: 2023-06-29  

Fund supported: 国家自然科学基金资助项目(51875408);天津市研究生科研创新项目(2021YJSB229).

作者简介 About authors

梁志文(1999—),男,硕士生,从事新型盘式行星永磁电机电磁特性研究.orcid.org/0009-0000-5356-8343.E-mail:2131050710@tiangong.edu.cn , E-mail:2131050710@tiangong.edu.cn

摘要

为了解决传统永磁电机低速直驱能力弱和输出单一的问题,提出新型盘式行星永磁电机. 所提电机集盘式电机与行星驱动于一体,具有空间磁场结构复杂、低速直驱和多自由度动力输出的特点. 分析所提电机的电磁性能,研究盘式转子与行星定子不均匀气隙的动态磁路模型,利用径向分环法将三维磁场模型分环为二维模型. 考虑气隙不均、定子开槽和漏磁的影响,引入不均匀气隙位置分布函数,建立所提电机的等效气隙磁阻模型和等效漏磁阻模型. 划分盘式转子和行星定子间的磁场区域,将动态磁路模型转化为时变的等效磁网络拓扑结构,计算得到电机的磁链和空载反电势. 通过有限元和样机实验验证动态磁网络模型的有效性. 相比三维有限元仿真,等效磁网络提高了计算效率,便于后续的结构参数分析和优化,为盘式行星永磁电机的磁场建模和分析提供了参考.

关键词: 不均匀气隙 ; 盘式行星永磁电机 ; 径向分环 ; 气隙分布函数 ; 等效磁网络

Abstract

A novel disc planetary permanent magnet motor was proposed to resolve the problem of weak low-speed direct drive capability and the single output of traditional permanent magnet motors. Combining a disc motor with a planetary drive, the proposed motor has the characteristics of a complex spatial magnetic field structure, low-speed direct drive, and multi-degree freedom power output. The electromagnetic performance of the proposed motor and the dynamic magnetic circuit model of the uneven air gap between the disc rotor and the planetary stator were analyzed, and the 3D magnetic field model of the motor was divided into 2D models by the radial sectors-of-rings method. Considering the influence of uneven air gap, stator slots, and leakage flux, the equivalent air gap reluctance and equivalent leakage flux reluctance models of the proposed motor were established by introducing the uneven air gap position distribution function. The magnetic field area between the disc rotor and the planetary stator was divided, and the dynamic magnetic circuit model was transformed into a time-varying equivalent magnetic network topology. The flux linkage and no-load back electromotive force of the motor were calculated. The validity of the dynamic magnetic network model was verified by finite element analysis and prototype experiment. Compared with 3D finite element simulation, the equivalent magnetic network is convenient for the analysis and optimization of structure parameters, and it improves the calculation efficiency and provides reference for a magnetic field modeling and analysis of a disc planetary permanent magnet motor.

Keywords: uneven air gap ; disc planetary permanent magnet motor ; radial sectors-of-rings ; air gap distribution function ; equivalent magnetic network

PDF (4241KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

梁志文, 刘欣. 新型盘式行星永磁电机的等效动态磁网络模型. 浙江大学学报(工学版)[J], 2024, 58(10): 2182-2191 doi:10.3785/j.issn.1008-973X.2024.10.022

LIANG Zhiwen, LIU Xin. Equivalent dynamic magnetic network model of novel disc planetary permanent magnet motor. Journal of Zhejiang University(Engineering Science)[J], 2024, 58(10): 2182-2191 doi:10.3785/j.issn.1008-973X.2024.10.022

永磁同步电机(permanent magnet synchronous motor,PMSM)具有高转矩密度、高功率因数和高过载能力等优点,被广泛应用于工业制造、航空航天以及轨道交通等领域[1-3],应用领域的不断扩展对电机系统高品质的运行性能提出更高的需求[4]. 在低转速高负载领域,传统的永磁同步电机要与变速装置匹配才能实现低速高负载的动力输出,传统永磁电机的单一输出增加了对多电机系统总成的依赖,不利于电机动力的高效利用. 直驱系统省去中间的齿轮传动系统,消除了齿轮带来的振动噪声、机械磨损和动力损耗,因此现代工业动力系统更倾向于直接驱动[5-7]. 为了适应新技术的发展,满足工业多样化的生产需求,国内外许多学者在不断探索新型拓扑下的驱动方式. Jungmayr等[8-9]提出新型的轴向串联磁齿轮复合电机,在结构上实现了磁齿轮与外转子永磁电机的一体化,减少了电机系统的体积和重量; Zhu等[10]基于行星磁齿传动提出混动汽车用磁性行星齿轮永磁无刷电机,该电机采用非接触的行星磁齿轮传动结构,具有功率密度高,调速范围广,体积紧凑等特点. 磁性齿轮与永磁电机的结合拓宽了传统永磁电机的驱动模式,但未解决传统永磁电机输出口单一、直驱能力弱以及功率和能量分配效率低的问题. Liu等[11]提出结合盘式电机与行星转子的多端口空间电机,能够实现多自由度动力输出;由于特殊的半包角式定子结构,该电机存在明显的端部效应. 分布式的驱动模式增加了电枢励磁绕组,导致散热性较差,并使控制难度增加.

本研究提出新型盘式行星永磁电机(disc planetary permanent magnet motor),该电机在驱动方式和拓扑结构上有别于以往的磁齿轮复合电机,且相比于传统的盘式行星永磁电机,具有轴向尺寸短、散热性好、结构简单和易于控制等优势. 永磁行星驱动具有直驱、低噪音及过载保护等优点,能够在有限的空间里实现多自由度的动力输出,为多源动力领域带来更多的驱动模式和功能. 本研究所提电机具有电磁驱动与磁齿轮传动2种动力传输模式,实现了机电磁的有机结合. 作为新型盘式行星永磁电机的重要指标表征的电磁性能是本文研究的重点. 该电机空间结构比较特殊,磁场呈现三维特性,传统的二维电机的分析方法不再适用,计算精确的三维有限元会消耗大量的计算资源和时间,给新型电机的初期设计带来不便[12]. 等效磁网络法通过构建区域磁通网络对磁场进行整体的等效计算,该方法结合磁路法计算速度快和有限元法求解精度高的优点,在电机的初期设计阶段具有很大的优势[13-14]. 曹东辉等[15]提出基于网格剖分的等效磁网络建模法,在考虑漏磁与局部饱和的情况下,精确地评估了转子磁障结构对降低永磁体涡流损耗的作用,提高了计算精度. 宫晓等[16]考虑漏磁形式,对轴向盘式永磁电机进行设计及优化,建立了这种电机的等效磁路网络模型,模型的计算结果与有限元以及实验结果吻合较好. 张淦等[17]建立针对磁通切换型永磁电机的非线性磁网络模型,通过引入深度饱和系数计算电机的磁链和反电势,该方法在电机的初期设计中可节省大量的计算时间. Ding等[18]提出新的基于等效磁网络建模方法,通过建立模块化的磁网格,对永磁电机的性能进行优化设计,该方法无需重复编程和建模,在电磁性能分析上具有很大优势.

本文针对新型盘式行星永磁电机复杂的空间磁场,研究定转子间磁场的时变性等效动态磁网络模型. 1)利用径向分环法,将该新型电机的三维磁路模型等效分环为二维模型,分析二维等效模型的磁路变化. 2)针对该电机气隙不均的特点,构建不均匀气隙位置分布函数,对各段气隙磁阻和漏磁磁阻进行等效计算,进一步建立该电机在单个电周期下的等效磁网络模型. 3)通过求解该等效磁网络模型得到该电机的磁链和空载反电势. 4)通过有限元仿真和样机实验验证该电机等效动态磁网络模型的有效性.

1. 所提电机的结构及工作原理

图1所示,新型盘式行星永磁电机主要由3个基本驱动部件构成:盘式转子、行星转子和行星定子. 行星定子由硅钢片叠压而成,绕组槽内放置通电线圈. 盘式转子和行星转子表面嵌有N-S极相间的永磁体,永磁体的布置由定转子间的磁力啮合关系决定,盘式转子上的永磁体径向布置,行星转子的永磁体内嵌于转子表面,并与转轴平行,保证在动力传递过程中磁力能够有效传递. 行星定子绕组通入三相交流电后,表面产生旋转的磁场,驱动上下2个盘式转子旋转;同时盘式转子上的永磁体与行星转子上的永磁体发生磁力耦合,带动行星转子转动. 新型拓扑结构下的盘式行星永磁电机有2种驱动方式: 1)通电的行星定子与盘式转子间的电磁驱动, 2)盘式转子与行星转子间的磁力驱动. 行星转子受到上下2个盘式转子的同步驱动,与盘式转子之间存在传动关系,两者满足一定传动比,在电机结构层面实现盘式电机与行星驱动系统的一体化设计. 在实际应用中,可根据不同的工况选择输出端口,提高电机的工作效率.

图 1

图 1   盘式行星永磁电机结构

Fig.1   Structure of disc planetary permanent magnet motor


2. 动态磁路模型分析

2.1. 分环等效模型

盘式转子上永磁体的间距沿着盘式转子的径向不断变化,使得盘式转子和行星定子间的气隙磁场在空间上呈现不均匀的三维分布. 由图1可知,2个盘式转子关于行星定子对称,且盘式转子与行星定子的磁路互不影响,进行任意一侧盘式转子与行星定子之间的空间磁场分析即可得到分析结果. 采用径向分环的思想,将盘式转子与行星定子间的不均匀三维磁路模型等效简化为若干个二维磁路模型的叠加[19],分环模型如图2所示. 图中,Rr, i为盘式转子第i个分环模型的分环半径.

图 2

图 2   单个盘式转子与行星定子的分环模型

Fig.2   Sectors-of-rings model of single disc rotor and planetary stator


2.2. 分环模型动态磁路分析

在盘式行星永磁电机运行过程中,盘式转子和行星定子之间的磁通路径随着盘式转子的转动不断变化,磁路模型也随之改变. 根据磁路的周期性,将盘式转子与行星定子间的磁路划分为4个区域,分别分析每个区域里的磁路. 在盘式转子与行星定子的磁路分环模型中,假设盘式转子永磁体等间隔. 定、转子分环模型参数如图3所示. 图中,τ为盘式转子齿间距,wm为盘式转子齿宽,c为行星定子槽口宽度,Lstp为行星定子极靴宽度,hpm为永磁体厚度,Rs为行星定子半径,$ \gamma $为分环模型下盘式转子的转动路径. 各分环模型中盘式转子永磁体的齿间距τi满足:

图 3

图 3   定、转子的分环模型参数

Fig.3   Sectors-of-rings model parameters of stator and rotor


$ {\tau _i} = {\text{π }}{R_{{\text{r}},i}}/p - {w_{\text{m}}}. $

式中:Rr,i为盘式转子第i个分环的平均转动半径(i=1,2,3,···),p为盘式转子磁齿的极对数.

为了分析与行星定子A相绕组匝链的磁通,以A相绕组磁链为零的位置为起始位置,建立分环模型下定转子的等效磁路模型,起始位置处的等效磁路关于行星定子中心线对称,磁路模型如图4所示. 图中,Φ1Φ2分别为流向两侧的主磁通,RδRδi均为主磁通等效气隙磁阻,RstiRsyi分别为行星定子齿和轭部等效磁阻,Rryi为盘式转子轭部等效磁阻(i=1,2···),Rpm为永磁体内部磁阻,Fm为永磁体磁动势. 当盘式转子转动时,盘式转子上永磁体与行星定子极靴的相对位置发生变化,磁路的气隙长度改变,盘式转子上的磁通不断与绕组匝链. 在时刻t,分环模型中盘式转子的转动路径为

图 4

图 4   起始位置的等效磁路模型

Fig.4   Equivalent magnetic circuit model of initial position


$ \gamma = {R_{{\rm{r }},i}}\omega t. $

式中:$ \omega $为盘式转子的角速度. 为了简化磁路模型,仅分析与A相绕组匝链的磁路,对变化的磁路按盘式转子转过的区域分别进行分析,具体有以下4种区域的磁路模型.

1) 区域1:$ 0 < \gamma < {w_{\text{m}}}/4 $ .在区域1中,起始位置中的磁路1(图4)开始与A相绕组匝链,流进A相绕组中的磁通开始增加,主磁路如图5所示. 在该区域下,随着盘式转子的转动,盘式转子磁极1与行星定子齿2的相对面积开始增大,使得起始位置中的磁路1有一部分开始流向齿2,在经过齿2后穿过气隙流入到磁极2,这部分磁路穿过气隙与A相电枢绕组匝链,形成主磁通Φ3.

图 5

图 5   区域1中的磁通路径及等效磁路模型

Fig.5   Flux path and equivalent magnetic circuit model in region 1


2)区域2:$ {w_{\text{m}}}/4 \leqslant \gamma < {\tau _1} - {w_{\text{m}}}/2 $ .区域2中除了磁路3外,还有一部分为磁极1本身形成的闭合磁路Φσ,磁路模型如图6所示. 这是由于磁极1与定子齿2的相对面积增加,与定子齿1的相对面积减少,气隙长度变长,在磁路3的基础上,一部分磁通经过行星定子齿2后,开始分流到定子齿1,并且回到磁极1,形成闭合回路,这部分磁路Φσ也与A相绕组匝链.

图 6

图 6   区域2中的磁通路径及等效路模型

Fig.6   Flux path and equivalent magnetic circuit model in region 2


3)区域3:$ {w_{\text{m}}}/2 \leqslant \gamma < {\tau _1}/2 $ .在区域3中,随着行星定子齿1与磁极1的远离,之间的气隙长度增加. 齿1与磁极2的相对间距在缩短,啮合面积增加,气隙长度减小,磁通Φσ不再存在,区域3中的磁路模型如图7所示. 该区域产生流经齿1穿过气隙流向磁极2的磁通Φ4. 在此区域内,当磁极1的左端面与齿2极靴左端面重合时,齿2极靴与磁极1的相对面积最大,由于行星定子铁心的磁导率远大于空气,流进齿2的磁通量达到最大,流向A相绕组的磁通达到最大值.

图 7

图 7   区域3中的磁通路径及等效磁路模型

Fig.7   Flux path and equivalent magnetic circuit model in region 3


4)区域4:$ {\tau _1}/2 \leqslant \gamma < ({\tau _1}+w)/2 $. 在区域4中,磁极2与行星定子齿3间的气隙长度增加,相对面积减小,与行星定子齿1的气隙长度缩短,相对面积增大,使得磁通 Φ3(图7)流经定子齿3的气隙磁阻变大,根据磁阻最小原理,定子齿3不再有磁通流过. 这时磁通在穿过气隙后进入定子齿2,流经定子轭部流入相近的定子齿1,之后穿过气隙流入磁极2,形成一条闭合磁通Φ5Φ5磁路模型如图8所示.

图 8

图 8   区域4中的磁通路径及等效磁路模型

Fig.8   Flux path and equivalent magnetic circuit model in region 4


上述4个区域为盘式行星永磁电机前半个电周期内磁路的变化情况,根据磁路的对称性可知,后半个电周期的磁路模型与前半个电周期的4个区域的磁路模型对称,在分析这4个磁路模型后,由磁路对称原则[20]可以得到盘式行星永磁电机在单个电周期的等效磁路模型. 为了验证动态磁路划分的合理性及各磁路流经的区域,通过有限元软件对该电机二维模型进行动态仿真,不同临界位置的磁密分布如图9所示. 可以看出,随着盘式转子的转动,流经行星定子齿的磁通先是关于行星定子中心对称分布,而后在转动过程中不断地向一侧流入,并不断与绕组匝链. 磁密云图的分布情况验证了上述磁路划分的合理性,为后续的动态磁网络模型的搭建确立了基础.

图 9

图 9   不同临界位置的磁密云图

Fig.9   Magnetic dense cloud at different critical positions


3. 等效动态磁网络模型建立

在盘式转子转动过程中,定转子间的气隙在不断变化,气隙长度不是常数,而是关于盘式转子转动位置的函数. 气隙长度分布函数是建立等效气隙磁阻模型和构建等效动态磁网络的关键.

3.1. 不均匀气隙位置分布函数

为了得到随盘式转子位置变化的气隙长度,建立不均匀气隙位置分布函数δ(X). 如图10所示,以盘式转子和行星定子一侧的气隙分布为分析对象,行星定子沿x轴展开为等效的直线定子模型,在定转子磁通回路中,实际的气隙长度δ(虚线部分)可以等效为磁力线流出时与行星定子上的交点在展开直线定子过程中以O点为圆心半径为R(X)的圆扫过的圆弧δ(X)(实线部分). 由不均匀气隙长度等效模型得到盘式行星永磁电机的等效气隙长度为

图 10

图 10   不均匀气隙长度的等效模型

Fig.10   Equivalent model of uneven air gap length


$ \delta = \delta (X)+{\delta _0}. $

式中:$ {\delta _0} $为最小气隙长度,$ \delta $为等效气隙长度.不均匀气隙位置函数表达式为

$ \delta (X) = \frac{{{\text{π}} \theta }}{{{{180}^\circ}}}\sqrt {{X^2}+{{(\sqrt {{R_{\mathrm{s}}}^2 - {X^2}} +{R_{\mathrm{s}}})^2}}} . $

在建立定转子不均匀气隙长度等效模型时,未考虑行星定子开槽情况. 由于定子槽口区气隙与铁心磁导率相差大,导致槽口区的气隙磁密产生间隙突变,磁路会发生改变[21]. 为了得到实际气隙长度分布情况,以分环模型中半径为Rr,1处的定转子为对象,分析行星定子槽口处的气隙路径,建立槽口处的气隙磁通路径模型. 槽口处的气隙路径模型如图11所示.

图 11

图 11   行星定子齿槽气隙磁路模型

Fig.11   Air gap magnetic circuit model of planetary stator slot


$ \left.\begin{split} &\delta =\left\{ \begin{array}{*{20}{l}}\delta (X)+{\delta }_{0},& \text{非槽口区域};\\\delta (X)+{\delta }_{0}+\dfrac{{\text{π}} (c/2-\Delta X)}{2},& \text{槽口1区域};\\\delta (X)+{\delta }_{0}+\dfrac{{\text{π}} (m-n)}{2{\mathrm{cos}}\;\beta },& \text{槽口2区域}.\end{array}\right.\\&\qquad \delta (X) = \theta {\text{π}} R(m)/{180^\circ} ,\; n = {R_{\text{s}}}\sin \alpha \text{;} \\&\qquad \alpha = \frac{{{\text{π}} ({L_{{\mathrm{stp}}}}+1.5c)}}{{({L_{{\mathrm{stp}}}}+c)p}} \text{,} \\&\qquad \beta = \arctan \left(\frac{{\sqrt {{R^2}(m) - {m^2}} +{\delta _0}}}{m}\right) .\end{split}\right\}$

由式(5)得到定转子在半径Rr,1处的气隙长度,对比等效无槽行星定子的气隙长度,结果如图12所示. 可以看出,行星定子在槽口处气隙长度有明显的突变,此突变会影响流经槽口的气隙磁阻大小. 须用槽口气隙长度分布函数求解槽口处的气隙磁阻.

图 12

图 12   考虑槽口时的气隙长度对比

Fig.12   Air gap length comparison considering notch


3.2. 等效气隙磁阻模型

由动态磁路的分析可知,气隙磁通模型随着盘式行星永磁电机的转动而不断变化. 为了求得所提电机变化的气隙磁阻,建立盘式转子与行星定子之间的气隙磁通模型,示意图如图13所示. 根据气隙磁通的流向,气隙磁通模型可分为齿对磁通和槽间磁通,槽间磁通由于行星定子槽口的影响,磁路产生畸变,使得槽间磁通形状与齿对磁通有所不同.

图 13

图 13   气隙磁通模型(局部)

Fig.13   Air gap flux model (local)


图14所示,通过分析磁通模型可以得到4种典型的气隙磁阻模型,其中齿对磁通形状等效为扇形磁通和不规则磁通. 在进行齿对磁通求解时,扇形磁通两侧的气隙长度根据盘式转子与行星定子的相对位置并由不均匀气隙位置分布函数求得;将不规则气隙模型(阴影1、2)等效成扇形磁通模型(阴影2、3),便可按照扇形磁通求解气隙磁阻Rg2. 槽间磁通1为矩形加1/4圆环,槽间磁通2为矩形加1/4半圆,槽间磁通的气隙长度由有槽口气隙长度分布函数得到. 由上述气隙磁通模型可以得到4种典型气隙磁通的磁导表达式,分别为

图 14

图 14   4种典型的气隙磁通模型

Fig.14   Four typical air gap flux models


$ \frac{1}{{{R_{{\text{g}}1}}}} = \frac{{{\mu _0}{l_{\text{a}}}\Delta {X_i}}}{{[\delta ({X_1})+\delta ({X_2})]/2}}, $

$ \frac{1}{{{R_{{\text{g}}2}}}} = \frac{{{\mu _0}{l_{\text{a}}}\Delta {X_i}}}{{\delta ({X_3})}}, $

$ \frac{1}{{{R_{{\text{g}}3}}}} = \frac{{2{\mu _0}{l_{\text{a}}}}}{{\text{π}} }\ln \left(1+\frac{{{\text{π}} {r_1}}}{{2\delta (X)}}\right), $

$ \frac{1}{{{R_{{\text{g}}4}}}} = \frac{{{\mu _0}{l_{\text{a}}}k}}{{\delta (X)+k{\text{π}} /4}}. $

式中:$ {\mu _0} $为真空磁导率,$ {l_{\mathrm{a}}} $为分环模型轴向长度.

3.3. 漏磁模型建立

为了简化磁网络模型,建立动态磁路模型时只分析与绕组匝链的主磁通Φi和有效漏磁通Φσ,没有考虑未与电枢绕组匝链的漏磁通Φσm. 对盘式行星永磁电机在运行过程中定转子间产生的漏磁进行分析计算[22],具体漏磁形式如图15所示. 图中,ΦσbmΦσim分别为永磁体端部漏磁通和定转子齿之间的极间漏磁通. 可以看出,盘式行星永磁电机的漏磁[23]主要2个部分组成: 永磁体的端部漏磁, 盘式转子与行星定子齿之间的极间漏磁. 在盘式转子与行星定子间的磁通中,由于不均匀气隙的存在,盘式转子上永磁体齿与行星定子极靴有相对错位,导致永磁体一部分磁通在流经行星定子极靴后未经过行星定子轭部又流回盘式转子永磁体齿,这一部分漏磁通被称为极间漏磁. 通过分析发现,极间漏磁仅存在于磁齿与行星定子齿1、2之间,并随着永磁体与定子极靴的重合面积增加而不断减少;端部漏磁是由永磁体的端部而产生的漏磁,是永磁电机不可避免的漏磁情况.

图 15

图 15   气隙漏磁模型

Fig.15   Air gap leakage flux model


为了得到漏磁通的气隙磁阻,建立盘式行星永磁电机定转子间漏磁模型如图16所示,其中RσigRσbg分别为极间漏磁和端部漏磁对应的气隙漏磁磁阻,Rσist为极间漏磁经过的定子齿部磁阻,Rrym为漏磁通经过的转子轭部磁阻. 漏磁区域中各部等效气隙磁阻及磁导的表达式分别为

图 16

图 16   漏磁的等效磁网络

Fig.16   Equivalent magnetic network of leakage flux


$ {R_{\sigma {\text{ig}}}} = \frac{{[\delta ({x_1})+\delta ({x_2})]/2}}{{{\mu _0}{l_{\mathrm{a}}}({x_1} - {x_2})}}, $

$ R{'_{\sigma {\text{ig}}}} = \frac{{[\delta ({k_1})+\delta ({k_2})]/2}}{{{\mu _0}{l_{\mathrm{a}}}({k_2} - {k_1})}}, $

$ \frac{1}{{{R_{\sigma {\text{bg}}}}}} = \int_{\Delta x}^{{k_\sigma }} {\frac{{{\mu _0}{l_{\text{a}}}{\mathrm{d}}z}}{{{\text{π}} z}}} . $

3.4. 动态等效磁网络建立

由动态磁路模型分析可知,每支与A相绕组匝链的磁路都要经过行星定子齿2. 将划分的动态磁路进行整合并简化成具有时变性的二维等效磁网络模型,得到一对极下的等效动态磁网络模型,如图17所示. 根据基尔霍夫定律列出等效动态磁网络的磁动势方程矩阵. 各磁通在半个电周期内的磁网络磁动势表达式为

图 17

图 17   一对极下的等效磁网络模型

Fig.17   Equivalent magnetic network model for one pair of poles


$ {{{\boldsymbol{R}}}^{ - 1}}{{\boldsymbol{F}} = }{\boldsymbol{\varPhi }}; $

$ {{\boldsymbol{F}}} = [2{F_{\text{m}}}\quad{\text{ 0 }}\quad{F_{\text{m}}}\quad{F_{\text{m}}}\quad{{ 0\quad 0\quad 0 \quad0}}]; $

$ {\boldsymbol{\varPhi }} = {[{\varPhi _{1,2}} \quad{\varPhi _3} \quad {\varPhi _4} \quad {\varPhi _\sigma } \quad {\varPhi _{\sigma {\text{im}}}} \quad {\varPhi _{\sigma {\text{bm}}}} \quad {\varPhi _\delta } \quad {\varPhi _{\text{m}}}]^{\mathrm{T}}}; $

$ \;{\begin{split} & \boldsymbol{R}= \\ &\qquad \left[ \begin{matrix} {{R}_{1,2}} & 0 & 0 & 0 & 0 & 0 & 0 & 0 \\ 0 & {{R}_{3}} & -{{R}_{4}} & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & {{R}_{\sigma \text{ }}} & -{{R}_{\text{st}}} & 0 & {{R}_{\delta }}+{{R}_{\text{st}}} & {{R}_{\text{rym}}}+{{R}_{\text{pm}}} \\ 0 & 0 & 0 & 0 & 0 & {{R}_{\sigma \text{bm}}} & 0 & {{R}_{\text{rym}}}+{{R}_{\text{pm}}} \\ 0 & 0 & 0 & 0 & {{R}_{\sigma \text{i}}} & -{{R}_{\sigma \text{bg}}} & {{R}_{\delta }} & 0 \\ 0 & 0 & {{R}_{4}} & -{{R}_{\sigma }} & 0 & 0 & 0 & {{R}_{\text{ry}}} \\ 0 & 1 & 1 & 1 & 1 & 0 & -1 & 0 \\ 0 & 0 & 0 & 0 & 0 & 1 & 1 & -1 \\\end{matrix} \right]; \\ \end{split}} $

式中:Φ为各支路中的磁通向量,F为支路中的磁动势列向量,R为等效磁阻系数矩阵.

4. 有限元分析及实验验证

4.1. 电机三维有限元分析

盘式行星永磁电机的磁场呈三维分布,分析盘式转子与行星定子的磁场时,可对盘式转子和行星定子进行三维有限元仿真[24]. 在有限元仿真时,通过设置盘式转子的转速,得到A相绕组的磁链和空载反电势在单个电周期内的波形. 行星定子A相绕组磁链最大值与最小值时定转子相对位置及磁场强度分布情况如图18所示.

图 18

图 18   磁链极值的行星定子磁场分布

Fig.18   Planetary stator magnetic field distribution for flux linkage extreme values


4.2. 磁网络解析计算及实验验证

为了提高电机初期设计效率,通过有限元分析,在误差允许范围内,将盘式转子和行星定子模型的分环数最终确定为6个,利用等效磁网络法分别求解磁链和空载反电势,进行结果叠加,再与有限元和样机实验进行对比. 三维有限元和等效磁网络模型计算得到该电机A相的磁链波形如图19所示. 由于空气磁导率远小于铁心,使得气隙磁阻两端的磁压降远大于铁心磁阻,在忽略铁心磁阻,不考虑磁饱和的情况下,动态磁网络求解结果与有限元仿真得到的磁链波形基本吻合. 由图可知,动态磁网络计算结果的峰值与三维有限元仿真的平均相对误差为6.97%. 误差的主要原因是未考虑电机空间漏磁,使得磁网络计算结果偏大. 2种方法计算的磁链波形在波峰处有一定位置偏离,原因是在对盘式转子与行星定子分环时,将分环后的盘式转子等效成直线运动下的转子,行星定子与等效直线转子之间的各磁路区域与分环前实际定转子各磁路区域存在位置上的偏离.

图 19

图 19   一对极在单个电周期下的磁链波形对比

Fig.19   Comparison of flux linkage waveform for one pair of poles in one electrical cycle


等效磁网络和三维有限元所用的时间Ts以及对应的剖分网格数Ng表1所示,由表可知,采用等效磁网络法求解磁场的精度相对较高,耗时少,效率高,适用于新型电机设计的初期阶段.

表 1   不同有限元分析方法的计算资源比较

Tab.1  Comparison of calculation resources for different finite element analysis methods

方法名称NgTs/s
等效磁网络法32318
三维有限元法16 1571661

新窗口打开| 下载CSV


行星定子四极12槽/盘式转子8极样机如图20所示,其设计参数与有限元仿真一致,主要参数如表2所示. 在实验平台中,通过原动机拖动样机,测得样机的空载反电势波形. 通过等效磁网络求得该电机的磁链后,对A相磁链求导,得到盘式行星永磁电机A相空载反电势的波形.

图 20

图 20   样机及实验平台

Fig.20   Prototype and experimental platform


表 2   电机设计的主要参数

Tab.2  Main parameters of motor design

参数数值参数数值
行星定子外径Dsout/mm60行星轮轴向长度la/mm40
盘式转子外径Drout/mm180盘式转子磁齿数z18
行星转子外径Dpout/mm60行星轮磁齿数z24
永磁体厚度hpm/mm5行星定子极数4
永磁体长度lm/mm40绕组匝数100
行星轮4最小气隙长度δ0/mm1

新窗口打开| 下载CSV


$ {E_{\text{A}}} = - \frac{{{\text{d}}{{\varPsi }_{\text{A}}}}}{{{\text{d}}t}}. $

式中:EAA相绕组空载反电势,ΨAA相绕组磁链. 等效磁网络法与样机实验及三维有限元仿真得到的A相绕组一对极在一个电周期内的空载反电势波形如图21所示. 可以看出,利用磁网络法解析计算的空载反电势波形与样机实验及三维有限元的变化基本一致,能够反映该电机的空载反电势变化趋势和幅值大小. 磁网络求解得到的空载反电势波形相对于样机实验和有限元有一定的误差,这是求导对磁链偏差有放大作用导致的. 此外,A相绕组磁链波峰的偏离使得对磁链求导后空载反电势的波形受到有一定影响.

图 21

图 21   一对极在单个电周期下的反电势波形对比

Fig.21   Comparison of back electromotive force waveform for one pair of poles in one electric cycle


5. 结 论

(1)利用分环的思想,将盘式转子与行星定子间的三维磁路模型简化为二维磁路模型. 该方法能够节省大量的分析时间,提高所提新型电机初期设计效率.

(2)针对所提新型电机气隙结构不均匀的特点,在考虑行星定子槽口影响下构建不均匀气隙位置分布函数. 计算气隙磁阻时引入该函数,得到随位置变化的动态气隙磁阻模型,进一步提高气隙磁阻计算的精度.

(3)考虑不均匀气隙长度和漏磁的影响,建立所提新型电机的等效气隙磁阻和漏磁磁阻. 通过划分动态磁路区域,构建电机动态等效磁网络模型. 等效磁网络能够简化磁场的计算模型,有利于提高磁场计算效率.

(4)通过对比三维有限元仿真与等效磁网络模型的磁链计算结果以及实验样机的空载反电势结果,验证了所建等效磁网络模型的有效性和正确性. 等效动态磁网络模型可以节省计算时间,提高所提电机的初期设计效率,为后续该电机的电磁参数优化奠定基础.

新型盘式行星永磁电机仍处于初期研究阶段,通过进一步优化该电机结构参数来提高其气隙磁密,改善反电势的正弦性是未来研究的方向和重要课题.

参考文献

ZHU X, XIANG Z, QUAN L, et al

Multimode optimization research on a multiport magnetic planetary gear permanent magnet machine for hybrid electric vehicles

[J]. IEEE Transactions on Industrial Electronics, 2018, 65 (11): 9035- 9046

DOI:10.1109/TIE.2018.2813966      [本文引用: 1]

马伟明, 王东, 程思为, 等

高性能电机系统的共性基础科学问题与技术发展前沿

[J]. 中国电机工程学报, 2016, 36 (8): 2025- 2035

[本文引用: 1]

MA Weiming, WANG Dong, CHENG Siwei, et al

Common basic scientific problems and development of leading-edge technology of high performance motor system

[J]. Proceedings of the CSEE, 2016, 36 (8): 2025- 2035

[本文引用: 1]

张东, 邹国棠, 江建中, 等

新型外转子磁齿轮复合电机的设计与研究

[J]. 中国电机工程学报, 2008, 28 (30): 67- 72

DOI:10.3321/j.issn:0258-8013.2008.30.011      [本文引用: 1]

ZHANG Dong, CHAU K T, JIANG Jianzhong, et al

Design and research of a novel magnetic-geared outer-rotor compact machine

[J]. Proceedings of the CSEE, 2008, 28 (30): 67- 72

DOI:10.3321/j.issn:0258-8013.2008.30.011      [本文引用: 1]

ATALLAH K, CALVERLEY S D, HOWE D

Design, analysis and realization of a high-performance magnetic gear

[J]. Electric Power Applications, 2004, 151 (2): 135- 143

DOI:10.1049/ip-epa:20040224     

鲍晓华, 刘佶炜, 孙跃, 等

低速大转矩永磁直驱电机研究综述与展望

[J]. 电工技术学报, 2019, 34 (6): 1148- 1160

[本文引用: 1]

BAO Xiaohua, LIU Jiwei, SUN Yue, et al

Review and prospect of low-speed high-torque permanent magnet machines

[J]. Transactions of China Electrotechnical Society, 2019, 34 (6): 1148- 1160

[本文引用: 1]

JUNGMAYR G, MARTH E, SEGON G. Magnetic-geared motor in side-by-side arrangement-concept and design [C]// IEEE International Electric Machines and Drives Conference . San diego: IEEE, 2019: 847– 853.

[本文引用: 1]

NIELSEN S S, HOLM R K, RASMUSSEN P O

Conveyor system with a highly integrated permanent magnet gear and motor

[J]. IEEE Transactions on Industry Applications, 2020, 56 (3): 2550- 2559

DOI:10.1109/TIA.2020.2977877      [本文引用: 1]

ZHU X, XIANG Z, QUAN L, et al

Multimode optimization research on a multiport magnetic planetary gear permanent magnet machine for hybrid electric vehicles

[J]. IEEE Transactions on Industrial Electronics, 2018, 65 (1): 9035- 9046

[本文引用: 1]

LIU X, ZHENG X, WANG X

Modeling of uneven air gap magnetic field for disc planetary permanent magnet machine with segmented gap Halbach array

[J]. IEEE Transactions on Energy Conversion, 2023, 38 (1): 693- 702

DOI:10.1109/TEC.2022.3206403      [本文引用: 1]

LIU X, ZHENG X, WANG X

Analysis of end-effect torque and multi-objective optimization for multi-port spatial permanent magnet synchronous motor

[J]. IEEE Transactions on Magnetics, 2022, 58 (9): 8206711

[本文引用: 1]

RAFAQ M S, JUNG J W

A comprehensive review of state-of-the-art parameter estimation techniques for permanent magnet synchronous motors in wide speed range

[J]. IEEE Transactions on Industrial Informatics, 2020, 16 (7): 4747- 4758

DOI:10.1109/TII.2019.2944413     

赵玫, 于帅, 张华强

聚磁式横向磁通永磁直线电机的变磁导等效磁网络

[J]. 电机与控制学报, 2020, 24 (4): 12- 22

[本文引用: 1]

ZHAO Mei, YU Shuai, ZHANG Huaqiang

Variable permeability equivalent magnetic circuit network of flux-concentrated transverse flux permanent magnet linear machine

[J]. Electric Machines and Control, 2020, 24 (4): 12- 22

[本文引用: 1]

韩力, 韩雪峰, 欧先朋, 等

无刷双馈电机运行特性的等效磁网络分析法

[J]. 电机与控制学报, 2019, 23 (1): 27- 34

[本文引用: 1]

HAN Li, HAN Xuefeng, OU Xianpeng, et al

Operating performances analysis of brushless doubly-fed machine using magnetic equivalent circuit

[J]. Electric Machines and Control, 2019, 23 (1): 27- 34

[本文引用: 1]

吴迪, 辜承林

新型横向磁通永磁电机三维磁路分析

[J]. 中国电机工程学报, 2010, 30 (33): 90- 95

[本文引用: 1]

WU Di, GU Chenglin

3D magnetic circuit analysis of a novel transverse flux permanent magnet machine

[J]. Proceedings of the CSEE, 2010, 30 (33): 90- 95

[本文引用: 1]

曹东辉, 赵文祥, 吉敬华, 等

基于网格剖分法的磁通切换永磁电机等效磁网络建模与分析

[J]. 中国科学: 技术科学, 2021, 51 (8): 970- 980

DOI:10.1360/SST-2020-0388      [本文引用: 1]

CHAO Donghui, ZHAO Wenxiang, JI Jinghua, et al

Modeling and analysis of a mesh-based equivalent magnetic circuit for a flux switching PM machine

[J]. Scientia Sinica: Technologica, 2021, 51 (8): 970- 980

DOI:10.1360/SST-2020-0388      [本文引用: 1]

宫晓, 徐衍亮

轴向磁场盘式永磁电机等效磁路网络及气隙漏磁的分析计算

[J]. 电机与控制学报, 2013, 17 (10): 59- 64

[本文引用: 1]

GONG Xiao, XU Yanliang

Lumped parameter magnetic circuit analysis of axial flux permanent magnet motor and its analytical calculation of air gap leakage

[J]. Electric Machines and Control, 2013, 17 (10): 59- 64

[本文引用: 1]

张淦, 花为, 程明, 等

磁通切换型永磁电机非线性磁网络分析

[J]. 电工技术学报, 2015, 30 (2): 34- 42

DOI:10.3969/j.issn.1000-6753.2015.02.005      [本文引用: 1]

ZHANG Gan, HUA Wei, CHENG Ming, et al

Analysis of nonlinear magnetic network models for flux-switching permanent magnet machines

[J]. Transactions of China Electrotechnical Society, 2015, 30 (2): 34- 42

DOI:10.3969/j.issn.1000-6753.2015.02.005      [本文引用: 1]

DING L, LIU G, CHEN Q, et al

A novel mesh-based equivalent magnetic network for performance analysis and optimal design of permanent magnet machines

[J]. IEEE Transactions on Energy Conversion, 2019, 34 (3): 1337- 1346

DOI:10.1109/TEC.2019.2900263      [本文引用: 1]

朱熙, 范瑜, 吕刚, 等

单边盘式感应电机的数学模型与转矩分析

[J]. 中国电机工程学报, 2010, 30 (24): 69- 74

[本文引用: 1]

ZHU Xi, FAN Yu, LÜ Gang, et al

Modeling and torque analysis of a disc induction motor

[J]. Proceedings of the CSEE, 2010, 30 (24): 69- 74

[本文引用: 1]

张昌锦, 陈志辉, 梅庆枭, 等

横向磁通永磁同步电机非线性动态等效磁网络模型分析

[J]. 中国电机工程学报, 2019, 39 (1): 307- 314

[本文引用: 1]

ZHANG Changjin, CHEN Zhihui, MEI Qingxiao, et al

Nonlinear dynamic equivalent magnetic network model analysis of transverse flux permanent magnet synchronous motor

[J]. Proceedings of the CSEE, 2019, 39 (1): 307- 314

[本文引用: 1]

孙明灿, 唐任远, 韩雪岩, 等

轴向磁通永磁电机空载气隙磁场建模分析

[J]. 中国电机工程学报, 2018, 38 (5): 1525- 1533

[本文引用: 1]

SUN Mingcan, TANG Renyuan, HAN Xueyan, et al

Analysis and modeling for open circuit air gap magnetic field prediction in axial flux permanent magnet machines

[J]. Proceedings of the CSEE, 2018, 38 (5): 1525- 1533

[本文引用: 1]

QU R, LIPO T A

Analysis and modeling of air-gap and zigzag leakage fluxes in a surface-mounted permanent-magnet machine

[J]. IEEE Transactions on Industry Applications, 2004, 40 (1): 121- 127

DOI:10.1109/TIA.2003.821790      [本文引用: 1]

ZHAO J, MA T, LIU X, et al

Performance analysis of a coreless axial-flux PMSM by an improved magnetic equivalent circuit model

[J]. IEEE Transactions on Energy Conversion, 2021, 36 (3): 2120- 2130

DOI:10.1109/TEC.2020.3040009      [本文引用: 1]

颜建虎, 冯奕

聚磁式横向磁通永磁盘式风力发电机设计与分析

[J]. 中国电机工程学报, 2017, 37 (9): 2694- 2700

[本文引用: 1]

YAN Jianhu, FENG Yi

Design and analysis of a flux-concentrated transverse flux permanent magnet disk wind generator

[J]. Proceedings of the CSEE, 2017, 37 (9): 2694- 2700

[本文引用: 1]

/