考虑磁桥非线性的内置式永磁...机空载电磁性能通用解析模型_李世奇.pdf

2023 年3 月 电 工 技 术 学 报 Vol.38 No.6 第 38 卷第 6 期 TRANSACTIONS OF CHINA ELECTROTECHNICAL SOCIETY Mar.2023 DOI:10.19595/ki.1000-6753.tces.220486 考虑磁桥非线性的内置式永磁同步电机 空载电磁性能通用解析模型 李世奇 佟文明 贾建国 唐任远（沈阳工业大学国家稀土永磁电机工程技术研究中心 沈阳 110870）摘要 该文提出一种将子域法及等效磁路法相结合的新型解析模型，可在计及磁桥饱和的情况下，计算具有任意极槽配合的 V 型/U 型/一字型内置式永磁电机的空载电磁性能。建模时首先根据转子结构将永磁体进行等效，以便于各子域的建立；为了考虑转子磁桥处的饱和，采用等效磁路法来计算磁桥处的磁导率，以便于相应子域边界条件方程的建立。电机无槽时的空载气隙磁通密度分布可以通过求解各子域交界面处边界条件方程组得到，然后通过保角映射的方法可考虑定子开槽的影响。该文所提模型可用于计算内置电机的空载气隙磁通密度分布、反电动势及齿槽转矩。将模型计算结果与实验测试及有限元计算结果进行对比，验证了所提出模型的有效性。与有限元法相比，该文所提模型具有建模速度快、耗时短，同时可达到近似精度的优点，这为相关电机的初始设计及优化带来了方便。关键词：内置式永磁电机 子域法 等效磁路法 磁桥饱和 中图分类号：TM351 0 引言 相比于表贴式永磁电机，内置式永磁电机（简称内置式电机）由于具有高功率密度、宽速度范围和高磁阻转矩等特点，被广泛用于工业场合与智能家居1-3。电机精确的磁场计算是预测其电磁性能和判断设计合理性的前提4。目前，计算永磁电机磁场的方法主要包括有限元法与解析法。有限元法可精确地预测出电机的磁场分布，但由于内置式电机转子结构复杂多变，加上计算时需要高精度的网格，这导致有限元法的计算时间较长，不利于电机的初始设计5-6。解析法是快速准确计算电机磁场分布的有效工具。已有国内外学者提出了几种预测永磁电机磁场分布的解析方法，例如，基于拉普拉斯-泊松方程的子域法7-13、基于磁通管理论的等效磁路法14-20，以及基于绕组函数理论的解析方法21-22等。子域法是一种通过不同子域之间的边界条件求解每个子域的场控制方程的方法。1993 年，Zhu Zhiqiang 教授在文献7-8中首次提出了这种方法来计算表贴式永磁电机的磁场分布；基于该方法，文献9中计算了轴向磁通永磁电机的磁场分布；考虑到内置式电机转子结构变化引起的边界条件变化，参考文献10-13中分别计算了表面嵌入式永磁电机、V 型内置式电机、U 型内置式电机以及切向式电机的空载气隙磁场。然而，上述文献应用子域法计算电机磁场时忽略了铁心饱和的影响，这将导致计算结果与实际情况间存在误差，尤其是对于磁桥严重饱和的内置式电机。为了考虑铁心饱和的影响，一些学者采用等效磁路法来计算永磁电机的气隙磁场分布。文献14-15中采用三维等效磁路法分别计算了直线电机与轴向磁通电机的气隙磁场分布。文献16-20通过等效磁路对不同转子结构的内置式电机的气隙磁场进行了计算。然而，等效磁路法的计算精度取决于划分节点的数量，节点过少会导致计算精度降低，节点过多会使计算更加复杂，增大计算时间。此外，由于定转子相对位置的变化，气隙中各节点的连接关系在不同时刻需要重新设置，这将导致计算进一 辽宁省“兴辽英才计划”（XLYC2007107）和辽宁省百千万人才工程资助项目。收稿日期 2022-04-02 改稿日期 2022-05-12 1422 电 工 技 术 学 报 2023 年 3 月 步复杂化。文献21-22中将绕组函数与等效磁路法相结合，可在考虑铁心饱和的情况下快速计算内置式电机的电枢磁场分布，节省了计算时间，但该方法只能得到电机的径向气隙密度，无法计算切向气隙磁通密度。为了在内置式电机磁场分布的计算精度与计算时间上达到一个平衡，参考文献23-25提出了一种改进的模型，可在考虑铁心饱和的情况下预测切向式永磁电机的气隙磁场分布。但是这种解析模型仅适用于切向式电机，由于转子结构的不同会使边界条件方程发生变化，故该模型不适用于求解其他类型的内置式电机气隙磁通密度。文献26-27分别使用混合模型获得了多层及 V形内置式电机的气隙磁场分布。文献28中提出一种将内置式电机等效为表贴式电机的解析模型，可对切向式、一字型、V 型内置式电机的空载磁场进行计算。但是上述模型在计算过程中，无论磁桥尺寸为何，均认为磁桥磁通密度为一恒定值，这会给计算结果带来误差。考虑到磁桥磁通密度随其尺寸的变化，文献29将 V 型内置式电机分为两种结构，并将两种结构的磁场叠加得到了 V 型电机的气隙磁场，但合成磁场的谐波含量与有限元结果不同，使计算的空载性能产生误差。本文在文献29中模型基础上进行改进，提出了一种结合子域法及等效磁路法的新型通用解析模型。该解析模型可在考虑磁桥非线性的情况下，直接计算具有任意极槽配合的 V 型/U 型/一字型内置式电机空载气隙磁场，计算结果无需叠加，在保证模型通用性情况下进一步提高计算精度。为便于子域的建立，文中将不同转子结构的电机进行相应的等效，这种等效为分析其他具有不规则永磁体排布结构的内置式电机提供了思路；模型计算时通过局部等效磁路来计算磁桥处的磁导率，以便于相应子域边界条件方程的建立，进而考虑磁桥非线性的影响。将所提模型计算结果与实验测试及有限元结果进行了对比，验证了所提出模型的有效性。1 建模假设及电机不同拓扑等效过程 本文在建模时做出如下假设：（1）忽略定子铁心及转子除磁桥外铁心的饱和，仅考虑磁桥处的饱和，并认为其饱和程度均匀。（2）忽略端部效应。（3）电枢绕组开路，电机磁场仅由永磁体激励。不同内置式电机的拓扑及其等效结构示意图如图 1 所示，图中将等效后的永磁体及空气区域分为径向结构与切向结构，并给出了转子的主要尺寸参数。模型建立在 r 坐标系中，规定逆时针方向为电机旋转的正方向。（a）V 型内置式电机及其等效结构 （b）U 型内置式电机及其等效结构 （c）一字型内置式电机及其等效结构 图 1 不同永磁体结构内置式电机拓扑及其 等效结构示意图 Fig.1 Topology of interior permanent magnet motors with different magnet structure and its equivalent schematic diagrams 本文将永磁体与转子外径间的磁桥称为磁桥1，将相同极性永磁体间的磁桥称为磁桥 2。各拓扑结构在等效时保证永磁体的厚度及跨度角不变。根据总磁通不变原理，对于 V 型电机可建立 mx1fb22xlf1m1l1222=2=2llRwwlRRwhRw=+=（1）进而得到，V 型拓扑中 Rf表达式为 b2ml12fml12=12wlRwwRhRw+（2）式中，hm与 lm分别为永磁体的厚度与宽度；Rl为磁桥 1 内半径；w1为永磁体与磁桥 1 间空气宽度；w2为永磁体与磁桥 2 间空气宽度；为磁极跨度角；Rf为等效后切向结构外半径；lx为等效后径向结构 第 38 卷第 6 期 李世奇等 考虑磁桥非线性的内置式永磁同步电机空载电磁性能通用解析模型 1423 中永磁体的宽度；与1分别为等效后径向及切向结构跨度角；wb2为磁桥 2 的宽度。根据总磁通不变原理，对于 U 型电机有 m1m2x1fxlf1122(2)=2lllRlRRw+=+=（3）进而得到，U 型拓扑中 Rf表达式为 m2m11lf2=12llwRR+（4）式中，lm1与 lm2分别为 U 型电机等效前永磁体径向及切向的宽度。根据总磁通不变原理，对于一字型电机有 m1f2b212(2)2=2lRww=+（5）进而得到，一字型拓扑中 Rf表达式为 b2m2f2=2wlwR+（6）可以发现，图 1 中各电机等效后的拓扑均具有磁桥结构、径向结构（径向永磁体或空气区域）与切向结构（切向永磁体或空气区域），本文针对以上三种等效拓扑结构，提出一种通用解析模型，对电机的空载电磁性能进行计算，模型建立时为便于各子域及其相应边界条件的建立，做出如下简化：电机径向结构均为永磁体，其宽度为 RlRf；切向分布结构均为永磁体，其跨度角为1+2，忽略磁桥 2的存在。所提出的通用解析计算模型如图 2 所示（后文称之为简化结构），图中，Rr为转子铁心外半径，Rs为定子铁心内半径，Rm为等效后切向结构内半径，p 为电机极对数，、与分别代表模型中的不同子域，Ax与 Hx 分别为不同子域处的矢量磁 图 2 所提出通用解析计算模型示意图（简化结构）Fig.2 Schematic diagram of the proposed general model 位与切向磁场强度（x=1,2,3,4），各子域间边界条件将在后文给出，为磁桥饱和区域的跨度角23，有 b1b1r2=hwR+（7）式中，hb1与 wb1分别为磁桥 1 的长度与宽度。为了验证所提等效方法的有效性，本文通过有限元对电机原始结构、等效结构及等效后简化结构（即所提出通用模型）的空载气隙磁通密度进行计算。由于 V型电机结构等效后与 U型电机结构相似，同时一字型电机结构又可看成 V 型电机结构永磁夹角为 180 时的特殊情况，故本文以具有代表性的 V型电机结构为例，通过有限元对所提出的等效方法进行了验证，所计算出的不同情况下电机的空载气隙磁通密度波形对比如图 3 所示。可以看出，原始电机与等效后电机的空载气隙磁通密度波形吻合较好，进而证明了所提出等效方法的有效性。图 3 不同情况下电机空载气隙磁通密度波形对比 Fig.3 Comparison of no-load air-gap flux density under different conditions 通过图 3 结果还可以看出，相比于原始电机气隙磁通密度波形，等效后简化结构的气隙磁通密度形状几乎未发生变化，区别仅表现在幅值的不同。气隙磁通密度形状几乎不变的原因是磁桥 2 的存在几乎不影响永磁体磁通进入气隙的路径，使得模型简化前后的气隙磁通具有相似的分布规律，故气隙磁通密度的形状变化较小；气隙磁通密度幅值发生变化的原因是由于模型简化后，永磁体的总量增多，使得总磁通增大，忽略磁桥 2 又会使漏磁减少，故使得计算出的气隙磁通密度幅值增大。U 型结构电机与一字型结构电机的气隙磁通密度波形同样具有上述分布规律。基于上述分析可知，通过对图 2 所示简化结构建立子域模型，求解各子域矢量磁位可得到空载气隙磁通密度分布，然后通过等效磁路法对计算结果 1424 电 工 技 术 学 报 2023 年 3 月 幅值进行修正，即可得到各原始电机的空载气隙磁通密度分布，具体实现过程见下文。2 解析模型的建立 2.1 解析建模 根据材料的不同特性，本文将图 2 中建立的模型分为四个子域：气隙（）、磁桥（）、径向永磁体（）与切向永磁体（）。考虑到模型对称性，本文假设磁极中心处为=0，并对 0,/2p内的模型进行解析，其中 p 为电机极对数。电机各子域矢量磁位满足 220r22211MMMrrrrrrr+=+AAA（8）式中，A 为矢量磁位；0为真空磁导率；Mr与 M分别为磁化强度的径向及切向分量。对于气隙子域与磁桥子域，磁化强度为零，故这两个子域满足 22222110rrrr+=AAA （9）22222110rrrr+=AAA （10）对于径向永磁体子域，以图 2 中 S 极为例，其满足 M=Brem/0，M/r=Mr/=0，其中 Brem为永磁体剩磁；对于切向永磁体子域，其满足 M=0，M/r=Mr/=0。故对于径向与切向永磁体子域有 22rem22211Brrrrr+=AAA （11）22222110rrrr+=AAA （12）考虑到模型的周期性和对称性，各子域的矢量磁位通解12,23-25可表示为 11nsr1(,)=sin()npnpnnrrrCDnpRR=|+|A（13）02022 2 122rl=1(,)ln2 cos2mmmmmrCDrrrmCDRR=+|+|A（14）03032 2 133lf1(,)ln2 cos2gggggrCDrrrgCDRR=+|+|A remBr （15）()()11(21)(21)2244fm0rem122211(,)42(21)(1)sin2(21)2kkkkkkrrrCDRRBkrk+=|=+|+|+A（16）式中，n、m、g、k 分别为气隙子域、磁桥子域、径向永磁体子域及切向永磁体子域的谐波次数；C1n、D1n、C02、D02、C2m、D2m、C03、D03、C3g、D3g、C4k及 D4k为各子域待定系数。不同子域间交界面处，矢量磁位 A 及切向磁场强度 H 是连续的。故在气隙子域与磁桥子域交界面处，满足如下边界条件 11rr11rr(,)(,)22(,)(,)22RRHRHR+=+=AA（17）在磁桥子域与径向分布永磁体子域交界面处，满足如下边界条件 11ll11ll(,)(,)22(,)(,)22RRHRHR=+=+AA（18）在径向分布永磁体子域与切向分布永磁体子域交界面处，满足如下边界条件 11ff11ff(,)(,)22(,)(,)22RRHRHR=+=+AA（19）另外，由于与磁导率无穷大的铁心相邻的子域切向磁场强度为零，可得到以下边界条件 s1m(,)002(,)002HRpHR=+（20）上述边界条件在模型中的位置如图 2 所示，由各边界条件得到的包含待定系数的方程将会在附录中给出，其中需要注意的是，方程中磁桥子域的磁导率未知，需要通过局部等效磁路法求得。待求解的方程中磁桥磁导率应与各原始电机磁桥 1 处的磁导率相同，为了计算该处磁导率，建立了如图 4a 所示的电机通用等效磁路，图中，Gm、第 38 卷第 6 期 李世奇等 考虑磁桥非线性的内置式永磁同步电机空载电磁性能通用解析模型 1425 Gas、Gb1、Gb2、Gg分别为永磁体等效内磁导、转子空气区域磁导、磁桥 1 磁导、磁桥 2 磁导以及气隙磁导；r为永磁体等效磁通源；m、as、b1、b2、以及g为各支路磁通，以上磁导及各支路磁通均可根据文献23,30中的公式及迭代法计算得到，进而可计算出方程中所需的磁桥磁导率。（a）原始电机 （b）简化结构 图 4 原始电机及简化结构通用等效磁路 Fig.4 Equivalent magnetic circuit of different model 在求解方程组得到各子域待定系数后，可根据气隙子域矢量磁位计算出电机径向及切向气隙磁通密度表达式为 r11sr1cos()npnpnnnnprrBCDnprRR=|=+|（21）1111rrss1sin()npnpnnnD npC nprrBnpRRRR=|=|（22）为了考虑模型简化对空载气隙磁通密度的影响，需要对式（21）、式（22）计算结果的幅值进行修正，修正后的径向及切向气隙磁通密度为 rmodrBKB=（23）modBKB=（24）其中 gmodg=K （25）式中，Kmod为修正系数；g为图 4a 所示原始电机磁路中的气隙支路磁通；g为图 4b 所示电机简化结构等效磁路中的气隙支路磁通，该磁路中各参数 计算方法与图 4a 中参数计算方法相同。2.2 定子开槽影响 本文采用保角映射的方法来考虑定子开槽对气隙磁通密度的影响。该方法的主要思想是通过几种保角变换将实际开槽的平面转换为无槽平面，然后可通过计算出的复相对气隙磁导表示开槽的影响31。考虑开槽后的径向及切向气隙磁通密度分别为 srrabBBB=+（26）sarbBBB=（27）式中，a与b分别为复相对气隙磁导的实部与虚部。2.3 空载电磁性能计算 根据磁链分布可以得到电机的空载反电动势为 ddEt=（28）其中 srdNNBS=（29）式中，为绕组磁链；N 为每相串联匝数；为相绕组磁通；S 为相绕组区域的面积。电机齿槽转矩可由麦克斯韦张量法计算得到 2acogsrs0dL rB B=（30）式中，La为铁心叠压长度。结合以上分析，图 5 给出了本文所提通用解析模型计算流程。图 5 内置电机空载性能解析计算流程 Fig.5 Flow chart of no-load performance analysis 3 解析计算结果及其验证 为了验证所提出解析模型的有效性，本文通过解析法与有限元法分别对 V 型结构、U 型结构与一字型结构电机的空载电磁性能进行了计算，并通过搭建实验平台对一台 6 极/36 槽 V 型结构电机与一台 4 极/18 槽一字型结构电机空载性能进行了测试。所计算电机的主要参数见表 1，所测试样机图片如 1426 电 工 技 术 学 报 2023 年 3 月 图 6、图 7 所示。表 1 电机主要参数 Tab.1 Main parameters of the motors 数 值 参 数 V 型电机U 型电机 一字型电机额定功率/kW 11 20 15 额定转速/(r/min)8 200 4 500 20 000 额定电压/V 380 189 380 额定电流/A 17.6 87.5 27.7 定子内径/mm 80 127 70 转子外径/mm 77.6 125.6 67.4 极对数 3 4 2 槽数 36 48 18 铁心长度/mm 80 75 110 定子槽口宽/mm 2 1.8 2 磁桥 1 宽 wb1/mm 1.5 2 1.5 磁桥 1 长 hb1/mm 3.6 4 4.5 磁桥 2 宽 wb2/mm 1.5 2 磁桥 2 长 hb2/mm 3.8 4 空气宽度 w1/mm 2 5 2.5 空气宽度 w2/mm 1 1 永磁体宽 lm/mm 14.4 8(lm1),27(lm2)15 永磁体厚 hm/mm 4 5 4 永磁体材料 N40UH SmCo24 N42UH （a）定子 （b）转子 （c）整机 图 6 V 型结构样机图片 Fig.6 Prototype pictures of V-shape motor （a）定子 （b）转子铁心 （c）整机 图 7 一字型结构样机图片 Fig.7 Prototype pictures of straight-shape motor 3.1 空载气隙磁通密度对比 通过有限元计算出的 V 型结构、U 型结构与一字型结构电机空载磁通密度云图如图 8 所示。解析模型计算出的各结构电机在气隙中心处的空载气隙磁通密度分布与有限元结果对比如图 9 所示，可以看出，解析计算结果与有限元结果吻合良好，说明了解析模型预测空载气隙磁通密度的有效性。（a）V 型结构电机 （b）U 型结构电机 （c）一字型结构电机 图 8 电机空载磁通密度云图 Fig.8 No-load flux density cloud map of motors （a）V 型结构电机 （b）U 型结构电机 第 38 卷第 6 期 李世奇等 考虑磁桥非线性的内置式永磁同步电机空载电磁性能通用解析模型 1427 （c）一字型结构电机 图 9 空载气隙磁通密度波形对比 Fig.9 Comparison of no-load air-gap flux density 3.2 空载反电动势对比 通过解析模型及有限元对 V 型、U 型与一字型电机在转速为 1 500 r/min、4 500 r/min 与 15 000 r/min情况下的空载线反电动势进行了计算，并将计算结果与实验结果进行了对比，如图 10 所示，其中还包括相应的谐波分析。可以看出，对于 V 型电机，有限元、解析及实验得到的反电动势基波幅值分别为85.5 V、86.2 V 和 84.8 V，相应的波形畸变率分别为 2.57%、1.95%及 3.02%；对于 U 型电机，有限元与解析得到的反电动势基波幅值分别为 213.7 V 和 （a）V 型结构电机 （b）U 型结构电机 （c）一字型结构电机 图 10 空载线反电动势波形及其谐波分析 Fig.10 Back electromotive force waveforms and harmonic analysis 215.9 V，相应的波形畸变率分别为 4.98%与 4.7%；对于一字型电机，有限元、解析及实验得到的反电动势基波幅值分别为 318.7 V、325.3 V 和 320.4 V，相应的波形畸变率分别为 1.09%、1.62%及 1.55%。可以看出，所提模型可有效计算电机空载反电动势。3.3 齿槽转矩对比 通过解析模型计算出了 V 型、U 型与一字型电机在一个齿距下的齿槽转矩波形，并将其与有限元计算结果进行对比，如图 11 所示。可以看出，对于V 型电机，有限元与解析得到的齿槽转矩幅值分别为 368 mNm 与 403 mNm，计算误差为 9.5%；对于 U 型电机，有限元与解析得到的齿槽转矩幅值分别为 609 mNm 与 584 mNm，计算误差为 4.1%；对于一字型电机，有限元与解析得到的齿槽转矩幅值分别为 57 mNm 与 51 mNm，计算误差为 10.5%。可以看出，解析计算结果与有限元结果吻合良好，计算误差均在工程允许范围内。图 11 齿槽转矩结果对比 Fig.11 Comparison of cogging torque 3.4 模型适用条件与计算时间 在列写附录中的方程时认为永磁体产生的磁通全 1428 电 工 技 术 学 报 2023 年 3 月 部沿方向进入磁桥 1 的边界，故当1=(+)2时，磁桥子域满足/=0A。随着磁桥饱和程度的增加，磁通不再全部沿方向进入磁桥，按照磁桥 子域满足/=0A计算会产生误差。文献29中总结了不同磁桥饱和程度下/=0A的适用条件，故 本文所提解析模型的适用条件与文献29中的 相同。在利用有限元计算电机性能时，需要求解 Nmesh Nmesh的非线性矩阵，其中 Nmesh是有限元中的节点数，这意味着有限元的计算时间主要取决于网格剖分数量。对于本文提出的解析模型，求解等效磁路以考虑磁桥饱和时仅需计算一个 NmNm的非线性矩阵，其中 Nm为图 4 中等效磁路的节点数，其值远小于 Nmesh；除此之外的计算时间可以忽略不计，因为附录中的方程均是线性的。故相比于有限元法，本文所提模型在计算时间上具有显著优势。4 结论 本文提出一种将子域法及等效磁路法相结合的通用解析模型，可在计及磁桥饱和的情况下，计算具有任意极槽配合的 V 型/U 型/一字型内置式永磁电机空载电磁性能。为便于模型建立，提出了一种转子结构等效方法，这种等效为分析其他具有不规则永磁体排布结构的内置式电机提供了思路。通过解析模型计算出的空载气隙磁通密度、线反电动势及齿槽转矩结果与有限元和实验测试结果吻合较好。解析模型计算出的线反电动势幅值及其畸变率与有限元和实验测量结果间的误差均在 3%以内；解析计算出的齿槽转矩幅值与有限元结果间的误差在 10%以内，满足工程需求。与有限元法相比，本文所提解析模型可在保证计算精度的前提下，实现 V 型/U 型/一字型内置永磁电机空载性能的快速计算，这为相关类型电机的初始设计及优化带来了方便。附 录 根据式（17），可得到如下方程 r0202r111s11lnnpnnnRCDRCDR=|+=+|（A1）2 lr22112rs=12=mnpmmnnnRRCDCDRR|+|（A2）r11s2 r021222rbl=122(1)1npnnmnmmmRn CDRmRDCDR|=|+|（A3）式中，rb为磁桥的相对磁导率，且 1121211=sin()d1=coscos22npnpnpnp +|（A4）当 np2m/时 11212211222=sin()cosd2cos()coscos22=12 1()mnpmnpnpmnpnp+|+|（A5）当 np=2m/时 11212211=sin()cosd222coscos22=4npnpnpnpnp+|+（A6）根据式（18），可得到如下方程 reml0303l0202l2 l223r=1ln=ln1mmmmBRCDRCDRRCDR+|+|（A7）2 2 ll33224fr=12=gmggmmmRRCDCDRR|+|（A8）0203rblrml11=DDRR （A9）2 l22rbr2 l033334rmf=12=22 mmmggggmRCDRgRDCDR|+|（A10）式中，rm为永磁体的相对磁导率，且 第 38 卷第 6 期 李世奇等 考虑磁桥非线性的内置式永磁同步电机空载电磁性能通用解析模型 1429 1121322=cosdsin2mmm+=|（A11）当mg 时 11211422 2=coscosd22mg+|222222sincos()=1mmgmgg （A12）当m=g 时 11221422=cosd=22g+|（A13）根据式（19），可得到如下方程 1remf0303f(21)2rem1f44f5222m10ln=4(2)1(1)(21)(2)kkkkkBRCDRBRCDRRk+=+|+|+|（A14）12 f33l(21)2rem1f44f6222m10=4(2)2(1)(21)(2)gggkkkkkRCDRBRCDRRk+=+|+|+|（A15）1(21)2f4411mfrem122212 f03533611l=1(21)(21)22(1)4(2)=(21)(2)442 22kkkkggggkkRCDRRBkgRDCDR+|+|+|（A16）其中 112115112(21)2(21)=sind=cos2(21)2(2)kkk+|+（A17）当 2(1+2)g(2k+1)时 1121612(21)2=sincosd22kg+|+111212(21)cos()cos(21)2(2)22121kgkgk+=+（A18）当 2(1+2)g=(2k+1)时 11216212 2=sincosd2()sin2ggg+|+=（A19）根据式（20），可得到如下方程 s11r=npnnRCDR （A20）1(21)2m144mfrem122212=(1)(21)4(2)(21)(2)kkkkRCDRRkBk+（A21）将式（A1）式（A3）、式（A7）式（A10）、（A14）式（A16）、式（A20）、式（A21）联立方程组，通过数学软件取有限阶次，建立矩阵对系数进行求解。参考文献 1 赵文祥,刘桓,陶涛,等.基于虚拟信号和高频脉振信号注入的无位置传感器内置式永磁同步电机MTPA控制J.电工技术学报,2021,36(24):5092-5100.Zhao Wengxiang,Liu Huan,Tao Tao,et al.MTPA control of sensorless IPMSM based on virtual signal and high-frequency pulsating signal injectionJ.Transactions of China Electrotechnical Society,2021,36(24):5092-5100.2 赵方伟,王秀和,赵文良,等.内置式永磁同步电机动态偏心故障下的轴电压解析分析和削弱J.电工技术学报,2022,37(4):837-848.Zhao Fangwei,Wang Xiuhe,Zhao Wenliang,et al.Analysis and reduction of shaft voltage in interior permanent magnet synchronous motors under dynamic eccentricity faultJ.Transactions of China Electrotechnical Society,2022,37(4):837-848.3 顾理成,陈前,赵文祥,等.五相永磁容错电机的相间短路容错控制J.电工技术学报,2022,37(8):1972-1981.Gu Licheng,Chen Qian,Zhao Wenxiang,et al.Inter-phase short-circuit fault-tolerant control for five-phase permanent magnet fault-tolerant motorsJ.Transactions of China Electrotechnical Society,2022,37(8):1972-1981.1430 电 工 技 术 学 报 2023年3月 4 石玉君,程子活,蹇林旎.两种典型的场调制型永磁电机的对比分析J.电工技术学报,2021,36(1):120-130.Shi Yujun,Cheng Zihuo,Jian Linni.Comparative analysis of two typical field modulated permanent-magnet machinesJ.Transactions of China Electro-technical Society,2021,36(1):120-130.5 王明杰,贾宛英,张志艳,等.永磁直线同步电机空载反电动势和推力的解析计算J.电工技术学报,2021,36(5):954-963.Wang Mingjie,Jia Wanying,Zhang Zhiyan,et al.Analytical calculation of no-load eletromotive force and thrust in permanent magnet linear synchronous motorsJ.Transactions of China Electrotechnical Society,2021,36(5):954-963.6 Liu Feng,Wang Xiuhe,Xing Zezhi,et al.Reduction of cogging torque and electromagnetic vibration based on different combination of pole arc coefficient for interior permanent magnet synchronous machineJ.CES Transactions on Electrical Machines and Systems,2021,5(4):291-300.7 Zhu Ziqiang,Howe D.Instantaneous magnetic field distribution in brushless permanent magnet DC motors,part I:open-circuit fieldJ.IEEE Transa-ctions on Magnetics,1993,29(1):124-135.8 Zhu Ziqiang,Howe D.Instantaneous magnetic field distribution in brushless permanent magnet DC motors,part:armature-reaction fieldJ.IEEE Transactions on Magnetics,1993,29(1):136-142.9 马霁旻,王杜,曲荣海,等.基于有取向硅钢的轴向磁通开关磁阻电机准三维解析分析与设计J.电工技术学报,2018,33(17):4069-4077.Ma Jimin,Wang Du,Qu Ronghai,et al.Quasi-three-dimensional analysis and design of an axial flux switched reluctance motor based on grain oriented silicon steelJ.Transactions of China Electro-technical Society,2018,33(17):4069-4077.10 Zhang Zhen,Xia Changliang,Wang Huimin,et al.Analytical field calculation and analysis of surface inset permanent magnet machines with high saliency ratioJ.IEEE Transactions on Magnetics,2016,52(12):1-12.11 Xu Lei,Zhang Chao,Zhu Xiaoyong,et al.Indirect analytical modeling and analysis of V-shaped interior PM synchronous machineJ.IEEE Access,2019,7:173786-173795.12 Hajdinjak M,Miljavec D.Analytical calculation of the magnetic field distribution in slotless brushless machines with U-shaped interior permanent mag-netsJ.IEEE Transactions on Industrial Electronics,2020,67(8):6721-6731.13 Pourahmadi-Nakhli M,Rahideh A,Mardaneh M.Analytical 2-D model of slotted brushless machines with cubic spoke-type permanent magnetsJ.IEEE Transactions on Energy Conversion,2018,33(1):373-382.14 赵玫,于帅,张华强.聚磁式横向磁通永磁直线电机的变磁导等效磁网络J.电机与控制学报,2020,24(4):12-22.Zhao Mei,Yu Shuai,Zhang Huaqiang.Variable permeability equivalent magnetic circuit network of flux-concentrated transverse flux permanent magnet linear machineJ.Electric Machines and Control,2020,24(4):12-22.15 He Mingjie,Li Weiye,Peng Jun,et al.Multi-layer quasi three-dimensional equivalent model of axial-flux permanent magnet synchronous machineJ.CES Transactions on Electrical Machines and Systems,2021,5(1):3-12.16 陈威,吴桂初,方攸同.基于绕组分布函数理论和动态磁网络的两种内置式永磁牵引电机解析建模方法J.电工技术学报,2020,35(增刊2):377-386.Chen Wei,Wu Guichu,Fang Youtong.Two analytical models based on winding function theory and dynamic reluctance mesh for interior permanent magnet traction machinesJ.Transactions of China Electrotechnical Society,2020,35(S2):377-386.17 Ding Ling,Liu Guohai,Chen Qian,et al.A novel mesh-based equivalent magnetic network for perfor-mance analysis and optimal design of