开关磁阻电机磁场有限元分析

开关磁阻电机磁场有限元分析班级：电控（研）13-1姓名：万冬学号：4713206881摘要本文以开关磁阻电机（SR电机）磁场的有限元分析及稳态特性研究为主题开展研究工作，以期促进开关磁阻电机的推广应用及开关磁阻电机磁场理论的发展。针对开关磁阻电机定、转子极存在显著的边缘效应和高度的局部饱和特点，本文在系统地总结电机电磁场以及工程电磁场数值计算理论的基础上，基于大型有限元分析软件ANSYS，采用全场域二维电磁场有限元分析，对开关磁阻电机的磁场分布、静态特性等进行了大量的计算。通过仿真的结果，一方面得出了开关磁阻电机在几个典型转子位置下的磁场分布图，并总结了电机内磁场的分布规律；另一方面根据计算结果绘制出了开关磁阻电机的磁化曲线族，即i关系曲线，电感曲线和静态转矩曲线等参数，这些工作为分析开关磁阻电机的工作原理及开关磁阻电机的进一步开发和应用，建立开关磁阻电机合理的非线性模型提供了理论基础和可靠依据。21前言开关磁阻电机（SwitchedReluctanceMotor，SRM）是继变频调速系统、无刷直流电动机调速系统之后发展起来的最新一代无级调速系统，是集现代微电子技术、数字技术、电力电子技术、红外光电技术及现代电磁理论、设计和制作技术为一体的光、机、电一体化高新技术。它具有调速系统兼具直流、交流两类调速系统的优点。英、美等经济发达国家对开关磁阻电动机调速系统的研究起步较早，并已取得显著效果，产品功率等级从数w直到数百kw，广泛应用于家用电器、航空、航天、电子、机械及电动车辆等领域。开关磁阻电机是一种新型调速电机，调速系统兼具直流、交流两类调速系统的优点，是继变频调速系统、无刷直流电动机调速系统的最新一代无极调速系统。它的结构简单坚固，调速范围宽，调速性能优异，且在整个调速范围内都具有较高效率，系统可靠性高。主要有开关磁阻电机、功率变换器、控制器与位置检测器四部分组成。控制器内包含功率变换器和控制电路，而转子位置检测器则安装在电机的一端。利用电机电磁场理论和有限元法进行SR电机磁场分析与计算，在SR电机的研究中占据十分重要的地位，它是整个电机设计和运行性能分析的基础。由于SR电机结构与传统的交直流电机具有很大的差别，加之其过饱和特性，以路的观点进行电机性能的理论分析便显出很大的局限性；相反，以场的观点，全面、系统地分析电机性能，以便进行电机设计、性能分析及仿真计算，显示出极大的优越性。电机磁化曲线在SR电机电磁场分析与计算中占据十分重要的位置。通过这一特性的计算与分析，可以清楚地了解电机能量转换的方式与大小以及电机内部的饱和情况；同时，它是计算电感、磁通、转矩和功率的基础，亦是电机优化与仿真的依据。本文所借助的主要数学方法和分析手段是有限元法。目前，电磁场有限元分析大体上出现了以下几个方面的发展趋势：一是对原有方法不断地完善和改进。如有限元周期边界的新处理方法、网格快速可靠全自动自适应生成、后验误差估计与自适应新方法以及有限元多项式方法等。二是与其它数值计算方法相互耦合。除了传统的FEM-BEM法、有限元—模拟电荷法、有限元—积分方程法外，近年又出现有限元—级数耦合法及H棱边有限元—边界元耦合算法等。耦合法能实现不同方法的优势互补，解决多子域、多连通域的复杂问题。三是新方法的开发应用以及新技术的不断融入。如棱边有限元法、叠层有限元法、有限元的外推插值法、无限元法等。值得注意的是，由于神经网络和小波分析在电磁场中应用日益增多，已出现了“小波—Galerkin有限元法等新方法。32开关磁阻电机磁场的有限元分析由于SR电机的定、转子均为凸极结构，使它在运行过程中存在磁通周期性变化和严重的局部饱和现象，以路的观点进行电机性能的理论分析便显出很大的局限性，并且绕组的磁链和电流之间有着高度的非线性关系，同时磁链随转子位置的变化而变化，所以SR电机必须计算一组不同转子位置下的磁化曲线族),(i。由此可见，在开关磁阻电机的研究中，电磁场的计算具有重要的意义，它是整个SR电机甚至是整个SR电机驱动系统分析计算的基础，也就是说磁场计算的准确程度将直接影响到电机的转矩特性、磁链特性、电感特性的准确程度，将影响到对SR电机稳态运行性能分析的准确性以及控制策略的可靠性。因此对SR电机磁场的计算分析是很有必要的。有限元法与其它的电磁场数值计算方法相比，有着无可比拟的优势，它使得复杂结构、复杂边界情况的定解问题得到解答，并能胜任非线性问题以及时变场的计算和分析，因此，它是目前工程技术领域中实用性最强、应用最为广泛的数值计算方法。它的基本思想是将问题的求解域划分为一系列单元，单元之间仅靠节点连接，单元内部点的待求物理量可由单元节点物理量通过选定的函数关系插值求得。由于单元形状简单，易于由平衡关系和能量关系建立节点量之间的方程式，然后将各个单元方程“组装”在一起而形成总体代数方程组，加入边界条件后即可对方程组求解，从而得到求解域内节点上的物理量值。本文的主要内容之一即是采用二维有限元法对开关磁阻电机的磁场进行分析。本章主要介绍电磁场分析的有限元方法以及开关磁阻电机二维有限元模型的建立与求解。2.1开关磁阻电机磁场的有限元分析开关磁阻电机运行时，在定子和转子铁心中存在着整体饱和及严重的局部饱和，计算这种非线性磁场，可以采用有限元法通过计算磁场中的矢量磁位A来完成。本文所分析的开关磁阻电机样机的结构参数见表2-1。表2-1样机的结构参数定子极数8铁心长度180mm定子极弧21°转子极数6每相绕组匝数136转子极弧24°定子外径210mm气隙0.4mm定子轭高13.5mm转子外径115mm转子极高16mm转子轭高17.5mm2.1.1开关磁阻电机二维有限元模型的建立1.基本假设由于磁场随转子位置而时刻变化，采用部分场域分析的方法，边界条件较难确定，因此本文采用全场域分析的方法。为了建立SR电机内部磁场的微分方程，确定求解区域和有限元求解的边界条件，作如下假设：4（1）忽略电机端部磁场效应，磁场沿轴向均匀分布，矢量磁位A和电流密度J只有轴向分量Az和Jz，故磁感应强度只有Bx和By分量（以下将A简写为Az，J简写为Jz）；（2）磁滞效应的影响忽略不计，认为铁心冲片材料的磁化曲线是单值的；（3）磁场仅被限制于电机的内部，定子的外部边界及转子的内部边界认为是零矢量磁位线；（4）不计交变磁场在导电材料中如定子绕组及机座中的涡流反应，因此开关磁阻电机的磁场可作为非线性似稳磁场来处理。2.数学模型及边界条件根据以上假设，针对开关磁阻电机的全场域二维平面磁场，用矢量磁位表述的边值问题为01121zzzzAJyAyxAx（2-1）式中：为材料的磁导率；21，分别为定子外圆周和转子内圆周。上式等价于以下变分问题0min212122zzzzzzAdxdyJAyAxAAW（2-2）式中：为SR电机的求解场域。二维有限元法就是从式（2-2）出发，将求解区域剖分离散，在离散单元内构造矢量磁位的插值函数，然后利用插值方法将式（2-2）的条件变分问题，离散化为多元函数的极值问题，即化为一组关于各个节点矢量磁位的代数方程组，求解之得到矢量磁位的数值解。2.1.2开关磁阻电机二维有限元模型的求解采用三角形剖分单元，在单元e内采用线性插值，得到矢量磁位的线性插值函数3,13,131ANANAiii（2-3）式中：下标1，2，3表示三角形剖分单元e的顶点，它们之间是逆时针排列的。基函数N为3,2,121,iycxbayxNiiii，将前面的线性插值函数对x，y分别求偏导数，则得53322113322112121AcAcAceyAAbAbAbexA（2-4）已知整个求解区域的三角形单元的总和为E，那么能量泛函可以表示为求解区域内所有单元泛函之和。EezzzzzEeezdxdyJAyAxAAWAW122121（2-5）单元e的能量泛函可表示为32123322112332211381AAAJAcAcAcAbAbAbAWeZeze（2-6）求上式对三节点的一阶偏导eeeeeeeeeeeeeeeeeeeeeeeePAKPPPAAAKKKKKKKKKAWAWAW321321333231232221131211321（2-7）其中2,31s,41,，，lccbbKslslesl3,2,13elJPl，对e个单元进行总体合成，则有PAKPPPAAAKKKKKKKKKAWAWAWeeeeeeeNNeNeNeNeeN321321212222111211321（2-8）式中：N为节点总数。根据泛函极值条件可知0PAK（2-9）即6PAK（2-10）在SR电机中，铁心导磁材料的磁化曲线呈非线性，即式中μ不是常数，为磁感应强度的非线性函数，由（2-6）式可知，上式所示的代数方程组为一非线性方程组，系数[K]是未知量矢量磁位{A}的非线性函数。非线性方程组的求解一般归结为多次迭代求解同阶的线性代数方程组。常用的解法有欠松弛迭代法、最速下降法、牛顿-拉夫逊迭代法、共轭梯度法、波前法等。本文采用波前法求解非线性方程组，波前法仍然是高斯消元的一种形式，只是该方法中组装单元刚度矩阵与消元是同时进行的。在实际消元过程中，它不是等到所有单元全部组装进整体刚度矩阵后再消元，而是将与此节点相关的所有单元（即连接此节点的单元）组装进整体刚度矩阵就开始消元了，把已完成消元的结果写入文件，释放占用的内存，再存放其它数据，再进行消元，最后从文件中读出消元结果求节点位移。这种交替进行的方法，占用内存小，可以求解大规模问题，在求解电机内部电磁场时，它具有求解时间短、收敛性好的特点。2.1.3开关磁阻电机的解后处理本文对开关磁阻电机的有限元分析仅仅求出磁势的分布是远远不够的，还要得到其它物理参数或工程参数，例如磁通密度、磁场储能、磁链、电磁力等。通常把这个从势函数到具体物理量的计算过程称为“有限元的解后处理”。下面就本文所涉及的物理量介绍解后处理的方法。1.等势线图为了直观地和定性地表达电机场分布，需要用图形来表示场分布的性质和状况，等势线图便是最为通用的一种。将求解区域中磁势值相同的点用线联结起来，便形成了等势线图。一般情况下，可以令相邻两条等势线的势差相同。这样等势线的疏密就表示了场磁通量密度的强弱，也就是说，等势线密集的地方其磁通量密度也强。2.磁通量密度磁场的强弱程度通常可用磁通量密度B来表示。有限元方法直接得到的数值解为矢量磁势A在每个节点上的值，这些势函数在某个单元上的分布函数由该单元的几个节点上势函数值的插值表示。而磁通量密度则由势函数的方向导数求得，即xAyyAxABee（2-11）3.磁场储能由于开关磁阻电机定、转子由硅钢片叠压而成，工作时都会进入不同程度的饱和状态，对于这种非线性材料，磁场储能密度（即单位体积的储能）由下式定义BmHdBw0（2-12）这说明磁场储能密度由材料的B/H特性和工作磁通密度来决定。整个电机的总磁场储7能Wm写成以上储能对体积的积分，对二维问题来说体积积分被简化为面积积分再乘以铁心长度l，即dHdBldwlWEeBEemmee10122（2-13）4.磁通量穿过每一个曲面S的磁通量φm（