圆筒型无铁心永磁直线电机空载磁场解析计算

永磁无刷直流电机空载气隙磁场和绕组反电势的解析计算摘要：该文利用许-克变换构造了考虑齿槽效应的气隙相对比磁导函数，该气隙相对比磁导函数反映了齿槽效应对气隙磁场分布的影响，且这种影响的程度随气隙中的径向位置而变化。

在忽略铁心饱和的情况下，结合偏微分方程的解析算法，提出了一种考虑齿槽效应的永磁无刷直流电机空载气隙磁场分布和相绕组反电动势的解析计算方法。

计算结果与二维有限元计算结果对比，其计算波形和大小吻合很好。

证明此方法是正确的、可靠的，为永磁无刷直流电机优化设计和性能分析提供了基本分析手段。

关键词：永磁无刷电机；气隙磁场；反电势；解析计算1 引言对于永久磁钢表面安装的永磁电机，由于定子铁心开有若干槽，使气隙磁导并非均匀值，从而导致电机气隙磁场并非理想的梯形波，其中含有幅值较大的齿谐波，当电机旋转时会引起相绕组交链磁链的波动，使相绕组反电势出现波动，进而导致绕组相电流的脉动，引起电磁转矩的波动，最终引起电机的振动和噪声。

可见，要准确计算永磁电机的电磁转矩波动，首先应准确计算电机气隙内的磁场分布，从而可以准确计算出电机相绕组的反电动势变化波形和电机的电磁转矩波动，以确定有效的改进措施和控制策略。

而准确计算永磁电机气隙内磁场分布的关键是如何考虑齿槽结构对气隙磁场分布的影响。

在气隙磁场的求解方法中，有限元数值计算方法可以准确计算出气隙磁场的分布波形，具有通用性强、适用于各种媒质的特点。

但其前处理过程复杂、计算时间较长，对使用者有较高的技术要求，在电机优化设计中不便采用。

解析方法可以较准确地计算气隙磁场分布波形，同时可以观察到气隙磁场分布与结构尺寸之间的关系，具有很大的工程实用价值。

文献[1]、[2]利用解析方法对气隙磁场进行计算，求解出气隙磁场的分布波形，但文中忽略了齿槽的影响。

文献[3]讨论了永磁电机中定子斜槽（或转子斜极）、永磁体磁化方式、气隙长度、转子半径和永磁体极弧系数对气隙磁场分布的影响，给出计算气隙磁场分布的经验公式，由此计算出相绕组反电动势变化波形，可是文中忽略了齿槽的影响，公式的通用性也较差。

轴向充磁圆筒型永磁直线电机磁场解析黄克峰;李槐树;周羽【摘要】针对轴向充磁圆筒型永磁直线电机的磁场分布和气隙磁场的计算与径向充磁问题，提出了一种基于基本气隙磁场计算轴向充磁圆筒型永磁直线电机磁场的解析计算法．利用许-克变换计算出了该种结构的无槽和开槽后的基本气隙磁场强度，再计算出气隙最小处的气隙磁场强度，最终得出整个电机的气隙磁场．以一台10极9槽轴向充磁永磁直线样机为例来验证所提方法的有效性和准确性．结果表明：解析法和有限元法计算气隙磁场强度的变化规律完全一致，峰值分别为1．3 T和1．28 T，两者结果吻合较好．所提方法是一种快捷、有效、计算精度高的方法．%Aiming at the problem of the air⁃gap magnetic field distribution and calculation of the axially magnetized tubular permanent magnet linear motor(TPMLM), which is very different from the radially magnetized tubular per⁃manent magnet linear motor, a new analytic method of calculating air⁃gap magnetic field of the axially magnetized TPMLM was proposed for this paper, which was based on the basic air⁃gap magnetic field. This tool had calculated the basic air⁃gap magnetic field intensity of the slotless and slotted structure by using the Schwarz⁃Christoffel meth⁃od. And the air⁃gap magnetic field intensity of the motor was derived from the calculation of the magnetic field in⁃tensity at the minimum air⁃gap. A10⁃pole 9⁃slot axially magnetized TPMLM was studied to verify the validity and accuracy of the proposed method. The results showed that the peaks of the air⁃gap magnetic field calculated by the analytical method and finite element method were 1.3 T and 1.28 T. The variation principle of theair⁃gap magnetic field obtained from the two methods above had the same trends. The two results were in good agreement. The pro⁃posed method for calculating motor magnetic field is proven to be quick, effective, and accurate.【期刊名称】《哈尔滨工程大学学报》【年(卷),期】2013(000)007【总页数】6页(P883-888)【关键词】气隙磁场;轴向充磁圆筒型永磁直线电机;解析法;许-克变换【作者】黄克峰;李槐树;周羽【作者单位】海军工程大学电气工程系，湖北武汉430033;海军工程大学电气工程系，湖北武汉430033;海军工程大学电气工程系，湖北武汉430033【正文语种】中文【中图分类】TM352.9轴向充磁圆筒型永磁直线电机分析和设计的前提是对电机内的磁场分布进行准确计算.目前旋转永磁电机内的磁场计算方法有解析公式法［1-8］、有限元数值计算法［9-11］.利用有限元法能够考虑电机在实际运行中存在的磁路饱和、齿槽效应和绕组涡流等因素的影响.解析公式法具有物理概念清晰、计算量小、快捷等优点，便于直观考察结构尺寸、材质等对电机磁场的影响，这样能够快速有效地对电机的结构进行优化设计.对于永磁旋转电机在结构上的传统性利用解析法能够准确进行计算，但是轴向充磁圆筒型永磁直线电机是一种新的电机类型，并在结构上具有特殊性，利用解析法不能对该类电机进行有效计算，本文将采用基本气隙磁场法并利用许—克变换得出的隐函数公式进行准确计算.1 无槽轴向充磁磁势的解析模型1.1 基本气隙磁场的定义基本气隙磁场定义［12］如下:在圆筒型永磁直线电机磁场分布中最小气隙δmin处的磁感应强度密度定义为1T时的气隙磁场分布，即式中:Fδmin为最小气隙处的磁势，Bδmin为最小气隙处的磁场强度.1.2 无槽轴向充磁磁势的解析模型为了建立圆筒型永磁直线电机拓扑结构的磁场分布，假设:1)电机轴是无限长，因此磁场分布是轴对称并在x轴周期分布.2)电枢是无槽，铁心磁导是无穷大.无槽轴向充磁圆筒永磁直线同步电机(tubular including interior PM linearmotor，TIPMLM)结构如图1所示.该电机的特点是采用轴向充磁的永磁体，磁体在x轴方向交替更换极性，与高磁导铁心结合，形成若干个磁极，在圆柱气隙空间产生磁场，从而产生轴向电磁推力.永磁材料采用钕铁硼稀土材料，其去磁曲线如图2.图1中，τp为极距，h为转子侧铁心的长度，hm为永磁体充磁长度的一半.假定永磁体均匀磁化，其工作点为P( Hm，Bm )，回复线与退磁曲线重合，回复磁导率为μ0μr.其中，Hm为永磁体工作点的磁场强度，Bm为永磁体工作点的磁通密度.图1 无槽轴向充磁圆筒永磁直线同步电机结构示意Fig.1 The structural diagrammatic sketch of slotted TIPMLM由图2(a)，得到式中:Br为永磁材料的剩磁感应强度，μr为永磁材料相对磁导率，μ0为真空磁导率.轴向充磁圆筒型永磁直线电机中的定子内表面磁势分布可以有2种:1)如果将整个电机的磁势分布认为是由极间线性变化到极面，到极面下保持不变可以得到梯形波磁势如图3(a);2)如果将极间部分认为是非线性变化，从而等效成极间漏磁可以得到的磁势分布为方波磁势如图3(b).由图3可得出电机磁势分布的数学模型:式中:F1、F2分别为梯形波和方波的磁势，Hm为永磁的磁化强度，hm为永磁体充磁长度的一半，τp为极距;Fm为方波磁势的峰值，x为定子内侧的位置.图2 永磁体的退磁曲线Fig.2 The demagnetizing curve with PM图3 无槽轴向充磁定子侧磁势分布Fig.3 The scatter with magnetic potential diagram of stator in TIPM LM2 无槽轴向充磁气隙磁场的解析计算2.1 磁势为梯形波分布气隙磁场计算由于电机不开槽，气隙均匀可得将式(5)代入式(2)得式中:B1m为梯形波磁势的磁场强度，δ为不开槽下电机的均匀气隙.将式(6)代入式(3)，可得出磁势为梯形分布的气隙磁场磁场强度分布:2.2 磁势为方波的气隙磁场计算在1.2节假设的基础上增加2个假设:1)忽略导磁材料的磁压降影响;2)永磁体漏磁部分以漏磁系数修正.由图2(a)可得:B=Bmr，则式中:φr为永磁体的虚拟内禀磁通，φr=Br Am;φm为永磁体的输出磁通，φm=Bm Am;φmσ为永磁体的内部漏磁通，永磁体输出磁通与外磁路磁通相等，即式中:φ为磁路中的主磁通，φσ为外磁路漏磁通.由于外磁路结构复杂，较难准确计算φσ，为简便，本文用漏磁系数修正主磁通.即φm=σ0φ，其中σ0为漏磁系数.因此，只要计算磁路中主磁通φ即可，φ与定子表面磁场分布有关.下面分段计算定子表面磁场强度:极间部分的磁场强度可用极间漏磁进行等效，运用许—克变换进行计算;极面部分的磁场强度，气隙均匀可以用磁势进行计算.1)极间部分(图1中的AB段)磁场强度计算.计算AB段磁场强度时有2种方法，即用许—克变换得到的隐函数公式和近似计算的tan公式.隐函数公式:式中:αp为极弧系数，为定子侧任意位置的基本气隙磁场强度.tan 公式［12］:2)极面部分(图1中的BC段)磁场强度计算BC段的磁场强度计算:只要计算出电机在最小气隙处的磁场强度就可以求出整个电机的磁场分布.tan公式法、梯形波磁势法和隐函数法这3种方法计算结果如图4.隐函数法和有限元法的对比如图5，由图5可以看出:1)隐函数法和有限元法得到的结果规律是基本一致;2)2种方法得到的结果在峰值上也一致，都是1.212 T，只是在永磁体和动子铁心粘结处存在微小的差别.而从图4可看出隐函数法和其他2种方法存在较大差别，这也说明了利用隐函数法计算气隙磁场是有效的.图4 3种解析法计算结果对比Fig.4 Comparison of the results of three analytical methods图5 隐函数法与有限元法结果的对比Fig.5 Comparison of the implicit function with the FEM3 开槽轴向充磁气隙磁场的解析计算3.1 开槽后相对气隙磁场强度电机开槽后气隙磁场较为复杂，利用许—克变换可以很好地对这类问题进行计算(图6).把齿部气隙磁场强度值看作单位，利用许—克变换计算开槽后槽部相对齿部的相对气隙磁场强度.在计算时先假设:1)定、动子只有一面开槽，另一面光滑;2)铁心的磁导率为无穷大，其表面为等磁位面;3)槽深和槽节距都是无穷大.图6 开槽后的电机计算示意Fig.6 The calculated diagrammatic sketch of slotted TIPM LM按照上述假设，矩形槽在z平面和w平面的情况如图7 所示［13］.图7 单面开槽时气隙磁场的变换Fig.7 The transform course with single slotted sided air-gap field通过许—克变换可得到:式中，-1＜u＜1，b为定子侧槽口宽，a= 为开槽后的磁场强度与不开槽的磁场强度的相对值(如图8).图8 定子侧开槽后与无槽时的相对值Fig.8 The relative value of slotted and slot3.2 开槽后气隙磁场强度根据2.2节计算的不开槽气隙磁场强度，可以得出开槽后的气隙磁场强度:1)极间部分(图1中AB段):根据定子所处的位置，利用公式:分别求出和u，再利用u求出 .2)极面部分(图1中BC段):根据定子所处的位置，利用公式:求出u，再利用u求出 .计算得到开槽后的基本气隙磁场强度如图9，从图9中可以得到:1)在气隙最小处的基本气隙磁场强度仍为1;2)在极面与槽口相对时，此处的基本气隙磁场强度不再为1，而由于开槽的影响，大小随着所对槽位置的变化而变化.图9 定子侧开槽后的基本气隙磁场强度Fig.9 The basic air-gap field with slotted TIPM LM4 10极9槽电机解析计算气隙磁场的分布可以借助有限法(充分考虑磁饱和和电机的运行状况)进行验证.只要解析法计算得到的气隙磁场强度和有限元法计算的结果吻合就可以验证该方法的正确性.本文以10极9槽样机为算例进行验证分析.轴向充磁圆筒型永磁直线电机其主要结构尺寸参见表1.表1 轴向充磁圆筒型永磁直线电机主要结构参数Table 1 The key size of TLIPM linear motor参数值参数值不锈钢轴/mm 15 电机极数 10转子外径/mm 30 电机槽数 9气隙大小/mm 1.5 定子材料硅钢片定子外径/mm 48 充磁方向轴向极距 /mm 18 永磁体材料钕铁硼不开槽的有限元法和隐函数法的对比如图5;开槽后的有限元法和隐函数法结果如图10、11.图10轴向充磁圆筒型永磁直线电机开槽磁场强度云图Fig.10 The magnetic nephogram of slotted TPM LIM图11 许—克变换解析法和有限元法计算结果对比Fig.11 The results match the Schwarz-Christ of fel with the FEM从图11可得:1)2种方法计算结果的规律是一致的;2)隐函数法计算的气隙磁场强度的峰值为1.3 T，有限元法结果为1.28 T，两者的大小在误差范围内;3)图11中有2处位置的结果存在较大偏差，这个偏差引起的原因是在该位置槽基本上与极间部分相对，这样在电机运行中该处就存在聚磁效应，而运用隐函数法解析分析时该因素考虑不到.总体来说隐函数法和有限元法很吻合，存在的误差在工程允许范围内，也验证了本文中的隐函数法解析分析的正确性.5 结论本文研究结果表明所提方法具有如下特点:1)与有限元法相比，虽存在工程允许范围内的误差，但更为快捷方便、计算量小，非常有利于电机优化设计;2)解析法能够直观的认识电机结构参数影响电机磁场强度的规律;3)计算精度能满足工程设计需要.但仍然存在一些不足:该研究电机气隙磁场的解析方法只能得到x关于气隙磁场强度的函数，不能得到关于x的函数，这样将不利于研究电机结构参数对气隙磁场强度影响规律.因此很有必要对关于x的函数进行进一步研究.参考文献:【相关文献】［1］NICOLA B，SILVERIO B，DARIO D C.Tubular linear permanent magnet motors:an overall comparison ［J］.IEEE Transactions on Industry Applications， 2003， 39(2):466-475.［2］WANG JB，DAVID H.Tubularmodular permanentmagnet machines equipped with quasi Halbach magnetized magnets part I:magnetic field distribution，EMF，and thrust force［J］.IEEE Transactions on Magnetics，2005，41(9):2470-2478.［3］WANG JB，DAVID H.Tubular modular permanent magnet machines equipped with quasi Halbach magnetized magnets part II:magnetic field distribution，EMF，and Thrust Force［J］.IEEE Transactions on Magnetics，2005，41(9):2479-2489.［4］WANG J B，DAVID H，LIN Z Y.Design optimization of short stroke single phase tubular permanent magnet motor for refrigeration applications［J］.IEEE Transactions on Industry Electronic，2010，57(1):327-334.［5］WANG J B，WANG W Y，KAIS A.A linear permanent magnet motor for active vehicle suspension ［J］.IEEE Transactions on Vehicular Technology，2011，60(1):55-63. ［6］卢晓慧，梁加红.表面式永磁电机气隙磁场分析［J］.电机与控制学报，2011，15(7):14-20.LU Xiaohui，LIANG Jiahong.Air gapmagnetic field analysis of surface-mounted permanent magnet motors［J］.Electric Machines and Control，2011，15(7):14-20.［7］陈阳生，林友仰.永磁电机气隙磁密的分析计算［J］.中国电机工程学报，1994，14(5):17-26.CHEN Yangsheng，LIN Youyang.Calculation of air-gap flux density in permanent magnet machines［J］.Proceedings of the CSEE，1994，14(5):17-26.［8］王兴华，励庆孚，王曙鸿.永磁无刷直流电机磁阻转矩的解析计算方法［J］.中国电机工程学报，2002，22(10):104-108.WANG Xinghua，LI Qingfu，WANG Shuhong.Analytical calculation of cogging torque in brush less DC motor［J］.Proceeding of the CSEE，2002，22(10):104-108.［9］路义萍，丰帆，孙明琦.同步电机定子与气隙流场数值计算与分析［J］.电机与控制学报，2011，15(8):47-51.LU Yiping，FENG Fan，SUN Mingqi.Numerical calculation and analysis of fluid flow field of stator and air-gap of a synchronous machine［J］.Electric Machines and Control，2011，15(8):47-51.［10］王咏菲.圆筒型永磁直线同步电机的电磁场分析与研究［D］.天津:天津大学，2008:18-26.WANG Yongfei.Magnetic field analysis on tubular permanent magnet linear synchronous motors［D］.Tianjin:Journal of Tianjin University，2008:18-26.［11］陈卫宝，范承志，叶云岳.圆筒永磁直线电机齿槽力综合优化［J］.哈尔滨工程大学学报，2011，32(7):944-947.CHEN Weibao，FAN Chengzhi，YE prehensive optimization of the cogging force of a tubular permanent magnet linear motor ［J］.Journal of Harbin Engineering University，2011，32(7):944-947.［12］黄克峰，李槐树，金朝红.低速大推力圆筒永磁直线电机磁场分析［J］.微特电机，2012，40(2):4-6.HUANG Kefeng，LI Huaishu，JIN Chaohong.Magnetic field analysis of low speed and high thrust tubular permanent magnet linear motor［J］.Small＆ Special Electrical Machines，2012，40(2):4-6.［13］胡之光.电机电磁场的分析与计算［M］.北京:机械工业出版社，1984:37-42.HU Zhiguang.Calculation of flux density in machines［M］.Beijing:China Machine Press，1984:37-42.。

圆环形永磁片产生的磁场的解析计算及其磁化强弱的检查徐建成1,徐建民2(11首都师范大学物理系,北京　100037;　21北京矿冶研究总院,北京　100044) 摘　要:介绍沿轴向充磁的圆环形永磁片产生的磁场的解析计算方法,得出与圆环形永磁片同轴的圆所包围的磁通量和场内各点磁感应强度的计算公式,并与现有公式做了比较。

这些公式为通过测量圆环形永磁片产生的磁场来检查它磁化强弱的工作提供了理论依据。

关键词:选矿工程;圆环形永磁体;磁场计算;磁通量计算;磁选中图分类号:TD92411;TD457;O44112　文献标识码:A 文章编号:1001-0211(2006)04-0076-05收稿日期:2006-05-15作者简介:徐建成(1949-),男,北京市人,副教授。

沿轴向充磁的圆环形永磁体在空间产生的磁场计算有实际意义,然而商用计算软件编制的基础是否准确合理,对计算误差的处理是否适于特定情况,对计算结果会有很大影响。

有些文献报道中,公式推导不够严格[1-2]。

圆环形永磁片是永磁辊式磁选机的重要部件,它的磁性直接影响磁选机的性能,从一批充过磁的圆环形永磁片里挑出磁性大的片制造磁辊能保证磁辊的性能,下面的计算有助于提高这种挑选的可靠性,计算中采用国际单位制。

1　圆环形永磁片产生磁场的计算图1是过圆环形永磁片对称轴的截面图,其上下两个平面平行,内外两个圆柱面同轴,充磁方向与其对称轴平行。

计算中用柱坐标系(ρ,Z ,θ),Z 轴与圆环形永磁片对称轴重合,坐标原点在圆环形永磁片中心。

计算中假设圆环形永磁片是均匀磁化的,还假设永磁材料的相对磁导率等于1。

图1　过圆环形永磁片对称轴的截面图Fig 11　Section through symmetric axis of annularpermanent magnetic sheet111　磁化的圆环形永磁片的物理模型按照安培分子电流观点,磁化的永磁体外部任一点的磁场,与因永磁体磁化而在永磁体的侧面造成的一圈束缚面电流所产生的磁场一样[3-4]。