一种新的无轴承永磁同步电机径向悬浮力研究方法

内容需要下载文档才能查看

交流

河北

一种新的无轴承永磁

同步电机径向悬浮力研究方法

国网湖北省电力公司神农架供电公司

１引言

随着现代农业科技、电力电子技术、微电子技术、新型电机控制理论和稀土永磁材料的快速发展，使得永磁同步电动机在离心泵、机械加工、干燥机、高精数控机床等领域得以迅速的推广应用。永磁同步电机与传统的电励磁同步电机相比，效率高、损耗少、节电效果明显，这些优点在稀土永磁同步电机上表现的尤为显著，因而它成为了近年来广受研究并在各个领域中广泛应用的一种电动机［１］。无轴承永磁同步电机充分利用了磁轴承与电机产生电磁力原理的相似性，把磁轴承中产生径向力的绕组安装在电机定子上，通过解耦控制，独立控制了电机转矩和径向悬浮力。无轴承电机被广泛应用于需要免维修、长寿命运行，无菌、无污染以及有毒有害液体或气体的传输等场合究无轴承永磁同步电机是非常有价值的。

基于迦辽金加权余量法或变分原理的有限元法，最早用于力学计算中［４］，其具有以下优点：系数矩阵对称、正定且具有稀疏性，几何剖分灵活，有利于分析几何形状复杂的问题；适宜于处理非线性问题。本文在基于电磁场理论基础上推导出无轴承永磁同步电机径向悬浮力数学模型。基于有限元的优点，在Ａｎ－ｓｏｆｔ／Ｍａｘｗｅｌｌ有限元分析软件基础上，提出了一种简单的气隙磁密基波法对无轴承永磁同步电机径向悬浮力数学模型进行研究验证，仿真结果验证了径向悬浮力数学模型的正确性。

２径向悬浮力数学模型

忽略洛仑磁力和由于转子偏心而引起的径向作用力。不计电机磁饱和、涡流损耗与齿槽效应。假定铁芯和气隙的磁导率分并且μＦｅ远大于μ０，铁芯与气隙交界面上的法别为μＦｅ与μ０，

Ｈｍ，由于Ｈｍ≈０，作用向磁感应强度和切向磁场强度分别为Ｂｍ，在转子表面面积为元ｄＡ上的麦克斯韦力可表示为：

２２

ｄＦ１＝μＦｅ－μ（Ｂｍ２＋μＦｅμ０Ｈｍ）ｄＡ≈＝０００

［２－３］

汪泉汪明涛

将上式沿圆周积分，得径向悬浮力沿ｘ方向上的分量为：Ｆｘ＝πｌｒＢ１Ｂ２ｃｏｓ（λ－μ）

０

其中，Ｐ１＝１，Ｐ２＝２，ω１＝ω２。

同样的可得到径向悬浮力沿ｙ方向上的分量为：Ｆｙ＝πｌｒＢ１Ｂ２ｓｉｎ（λ－μ）０

正弦函数，令它们的幅值为：Ｆｒ＝πｌｒＢ１Ｂ２

０的值。

３电机的有限元模型与磁场分析

Ｐ２＝２，转矩绕组的额定电流为１０Ａ，径向悬浮力参数：Ｐ１＝１，绕组的额定电流为１０Ａ，定子槽数１２，定子铁芯外径１５５ｍｍ，定子铁芯内径９８ｍｍ，转子外径９３ｍｍ，转子内径５０ｍｍ，铁芯长度１３５ｍｍ，转矩绕组和径向悬浮力绕组每槽匝数为４０，转速３０００ｒｐｍ。

定子槽分成内层和外层两部分，分别嵌入径向悬浮力绕组和转矩绕组。图１为电机磁力线分布图。左边为只加转矩绕组电流时磁力线呈两极对称分布；中图表示的是只加径向悬浮力绕组电流时磁力线呈四极对称分布；右图为两绕组同时通过三相对称电流时合成磁力线分布，可以看出磁力线密集的位置位于ｙ轴的负方向上，而ｙ轴正方向上的磁力线则很稀疏，根据电磁场麦克斯韦力产生原理，可知电机沿ｙ轴负方向产生了径向悬浮力。

（６）（５）（４）

Ｆｘ正好为一余弦函数，而Ｆｙ为一分析Ｆｘ和Ｆｙ的公式可知，

分析（６）可得，若能获得气隙磁密幅值Ｂ１和Ｂ２便能得到Ｆｒ

研。因此，

（１）

ｌ为铁芯长度，ｒ为转子外半径，θ为式中，μ０＝４×１０－７Ｈ／ｍ，空间位置角。气隙中的合成磁密为：

Ｂ（θ，ｔ）＝Ｂ１ｃｏｓ（ω１ｔ－ρ１θ－μ）＋Ｂ２ｃｏｓ（ω２ｔ－ρ２θ－λ）

（２）

其中，Ｂ１，Ｂ２分别为永磁体、转矩绕组共同产生的气隙磁密幅值和径向悬浮力绕组单独产生的气隙磁密基波幅值；μ、λ为ω２为转矩绕组和径向悬浮力绕组的电角频率。初始相位角；ω１、

ｄＦ（θ，ｔ）＝可将转子表面单位面积上的径向悬浮力写成：［Ｂｃｏｓ（ω１ｔ－ρ１θ－μ）＋Ｂ２ｃｏｓ（ω２ｔ－ρ２θ－λ）］２ｄθ１

０

（３）

图1

磁力线分布

４径向悬浮力的有限元值和数学模型值

同时给转矩绕组和径向悬浮力绕组通以１０Ａ额定三相对称交流电时，在Ａｎｓｏｆｔ／Ｍａｘｗｅｌｌ中Ｆｘ和Ｆｙ的有限元波形曲线分别如图２所示。由图可知Ｆｘ为一余弦曲线，而Ｆｙ为一正弦曲线，与

46

2014年第8

期

内容需要下载文档才能查看

径向悬浮力数学模型分析结果一致。从图中读取其幅值记作Ｆｔ＝２８０Ｎ。

将表１可得曲线图５。

河北

交流

图2有限元波形

由前面分析得知，径向悬浮力的数学模型值的求取可转换成对Ｂ１和Ｂ２的求取，结合公式便能得到径向悬浮力的数学模型幅值。

当转矩绕组电流为额定值１０Ａ，径向悬浮力绕组不通电时，气隙磁密波形和进行ＦＦＴ分解后的波形分别如图３所示。

图5

径向悬浮力的比较

由图５中的曲线对比可知，随着径向悬浮力绕组电流的增两者数值基本相等，误差极小。且大，Ｆｔ和Ｆｒ基本上呈线性增加，Ｆｔ比Ｆｒ略小，造成这种现象的原因是Ａｎｓｏｆｔ综合考虑了磁饱和，定子磁势谐波、涡流损耗和齿槽效应等因素。通过有限元分析证明了径向悬浮力数学模型是正确的，为无轴承永磁同步电机的优化设计提供了理论依据。

６结论

本文使用了有限元的研究方式，控制绕组电流，进而分析了

图3

转矩绕组气隙磁密波形

无轴承永磁同步电机的内部磁场，证明了径向悬浮力产生原理是正确的。此外，还运用简单的基波幅值法研究了径向悬浮力的数学模型，与有限元仿真数值与数学模型计算值相结合，从而使径向悬浮力数学模型的正确性和有效性得到了验证。为无轴承永磁同步电机的进一步优化提供了借鉴。

参考文献：

[1]吴钦木,韦书龙,李捍东等.永磁同步电机驱动系统效

率优化控制参数变化影响研究[J].电机与控制应用,2012,39(6):18-23.

[2]A.Chiba,J.Asama.Influenceofrotorskewininductiontypebearinglessmotor[J].IEEETransactionsonMagnetics,2012,48(11):4646-4649.

[3]D.C.Fei,H.Q.Zhu.Studyondecouplingcontrolofbearinglesspermanentmagnetsynchronous

motorsbased

on

in-versesystemtheory[J].EngineeringScience,2005,7(11)48-54.

[4]朱熀秋，张涛.无轴承永磁同步电机有限元分析[J].中国电机工程学报,2006,26(3):136-140.

根据图３，电机气隙磁密的基波幅值Ｂ１＝０．１８Ｔ。

当径向悬浮力绕组电流为额定值１０Ａ，转矩绕组电流不通电时，气隙磁密波形和进行ＦＦＴ分解后的波形分别如图４所示。

图4各次谐波波形

根据图４，电机气隙磁密的基波幅值为Ｂ２＝０．２Ｔ。Ｂ２与其它已知量代入公式（６）中，将Ｂ１、得Ｆｒ＝２８２．５Ｎ。５径向悬浮力的的比较分析

当转矩绕组和径向悬浮力绕组分别通以１０Ａ额定三相对称交流电时，径向悬浮力的有限元模型值Ｆｔ和数学模型值Ｆｒ之间的误差很小，它们基本相等。即Ｆｘ、Ｆｙ的有限元波形与数学模型是一致的。接下来分析一般情况：

当转矩绕组电流保持额定值１０Ａ不变，Ｆｔ和Ｆｒ随径向悬浮力绕组电流的变化情况如下表所示。

表1径向悬浮力随径向悬浮力绕组电流变化的取值径

向

悬浮力绕组电流/A

0246810

Ft/N057115168225280

Fr/N057.5116169.5227282.5

2014年第8

期

47