基于ANSYS8.0的永磁直线电机的有限元分析及计算

浙江工业大学毕业论文 基于ANSYS8.0的永磁直线电机的有限元分析及计算 学生姓名：沈沉 指导教师：余佩琼 浙江工业大学信息 工程 路基工程安全技术交底工程项目施工成本控制工程量增项单年度零星工程技术标正投影法基本原理 学院电气工程系 摘 要 永磁直线电机是一种具有很高定位精度的新型电机。不同与其他励磁的直线电机，它采用永磁体作为励磁源。研究其磁场分布及力特性具有重要意义。相对于传统的解析法，有限元数值分析可以缩短电机的 设计 领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计 周期及减少设计成本，可对直线电机的磁场及力得出精确的分析。ANSYS8.0是一种在工程中广泛使用的有限元分析软件，采用该软件中的电磁场分析功能对永磁直线电机的磁场进行有限元的分析和计算，并在此分析的基础上对永磁直线电机的力场做进一步的计算和分析，对永磁直线电机的设计具有重要的工程意义。 通过电磁场的有限元数值分析方法，利用通用有限元分析软件ANSYS8.0建立平板型单边永磁直线电机的有限元模型，分析其2维静态磁场，得到初步的分析结果，并在这个分析的基础上对永磁直线电机的力场进行了进一步的分析，计算直线电机的推力和法向力，结合永磁直线电机的静态磁场，研究了永磁直线电机推力及法向力和电流变化的相互关系，对今后永磁直线电机的设计和研究具有一定的参考意义。 关键词 永磁直线电机、有限元、ANSYS、电磁场、推力、法向力 Finite Element Analysis and Calcultation of a Permanent Magnet Linear Motor Based on ANSYS8.0 Student: Chen Shen Advisor: Peiqiong Yu Department of Electric Engineering College of Information Engineering Zhejiang University of Technology Abstract The permanent magnet linear motor is a kind of new electrical engineering that has the very high fixed position accuracy.The differents between the permanent linear motor and the type of non-permanent is that it adopts thepermanent be the source of dlux Opposite in traditional resolution method,Finite element analysis can shorten the design period of the electrical engineering and reduce to the design cost,it also can get the analysis of a precision tu the magnetic field and fotce of the linear motor the ansys8.0 is a finitr element analvtical software. Throught the method of the Finite element for the electromagnetic analysis,we use ANSYS8.0 creat a model for the Permanent magnet lineat motor with finite element method.We analysis its 2-D setaic magnetic and get the first result. Then we analysis the force field by finish the analysis of its magetic field We calaulate the thrust and normal force combining the analysis of the permanent magnet linear motor,study the relationship between thecurrent and thrust,normal force.The work for this paper can give some help and advice to the study and design of the permanent linear motor Keywords PMLSM, ANSYS, FEM,Electromagnetic field, thrust, normal force 致 谢 在本次毕业设计的过程中，首先要衷心感谢我的指导老师余佩琼，在这一个学期的的毕业设计期间给于我很多的指导和帮助，在我遇到困难的时候，给了我许多有用的建议和提示。本文的顺利完成，离不开导师的关怀和悉心辅导，再次表示衷心的感谢。 同时，也感谢学校对毕业生学习工作的关心和支持，使我们能够顺利的进行毕业设计，感谢信息学院给我提供了良好的实验室，感谢学校图书馆及其网络资源，使我能查阅一些非常重要的参考资料。感谢我的朋友，同学对我支持，帮助。 最后感谢所有关心和支持毕业设计的学校领导和工作人员。 目 录 中文摘要 i 英文摘要 ii 致谢 iii 目录 iv 表列 x 图列 xii 第一章 绪论 1 1.1 直线电机的发展和应用 1 1.2 直线电机的分类 3 1.3 直线电机的工作原理 4 1.4 有限元分析及相关软件的介绍 5 1.4.1有限元分析简介 7 1.4.2通用有限元软件ANSYS8.0 7 1.5 课 题 快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题 的主要意义和论文的主要工作 7 1.5.1课题的主要意义 7 1.5.2论文的主要工作 8 第二章 ANSYS 8.0及其在电磁场有限元计算中的应用 9 2.1 电磁场的有限单元法 9 2.1.1 基本方程 9 2.1.2 位函数及边界条件 9 2.1.3 有限元法的求解 11 2.2 电磁力计算方法 12 2.2.1 虚功法 12 2.2.2 MAXWELL张量法 12 2.3 ANSYS有限元建模 14 2.4 网格的划分和单元的选定 17 2.4.1网格的划分 17 2.4.2单元的选定 17 2.5 边界条件的加载 18 第3章 永磁直线电机的电磁场有限元仿真分析……………………………………………………………………………………. 20 3.1 永磁直线电机的磁场分析方法概述 20 3.2 永磁直线电机的磁场有限元仿真模型 20 3.2.1永磁直线电机的结构 20 3.2.2永磁直线电机磁场有限元模型的建立 21 3.3 永磁直线电机的磁场分析……………………………………………………………………………………. ……………………………… 23 3.3.1永磁直线电机边端效应 23 　　 3.3.2永磁直线电机齿槽效应 24 3.3.3永磁直线电机磁场变化 24 第4章 永磁直线电机的力场特性的有限元仿真分析……………………………………………………………………… 25 4.1 永磁直线电机的矢量控制 25 4.2 推力和法向力与电流的大小的分析 27 4.2.1推力和法向力与电流的大小关系 27 4.2.2 推力系数 28 4.3 推力与法向力与电机电流相位的分析 28. 4.3.1推力与法向力电流相位关系 28 4.4 推力与法向力的波动问题 28 4.4.1永磁电机空载时的输出力 .29 4.4.2永磁电机额定电流下的输出力 31 第五章　结论与展望 34 参考文献 35 附录 ANSYS8.0使用说明 36 表列 表2－1 材料特性 14 表2－2 可选择单元类型 17 表3－1 永磁直线电机参数 21 表4－1 推力和法向力随电流大小变化数据 28 表4－2 推力和法向力随电流相位变化数据 30 表4－3 空载时推力和法向力波动变化数据 31 表4－4 额定电流下推力和法向力波动变化 32 图 列 图1—1 直线电机按工作原理的分类 3 图1－2 由旋转电机演变为直线电机的过程 4 图1－3 直线电机的工作原理 4 图1－4 有限元模型图 4 图1－5 最大切应力等应力线图 5 图2－1 永磁直线电机的2-D实体 15 图2－2 永磁直线电机的3-D实体模型 15 图2－3 定义完材料属性后的实体模型 16 图2－4 永磁直线电机的网络划分 17 图2－5 加载力标志和边界条件 17 图2－6 PMLSM的电流密度 17 图2－7 PLANE53 18 图2－8 边界条件示意图 19 图3－1 平板永磁直线电机的结构 21 图3－2 PMLSM绕组接线圈 21 图3－3 永磁直线电机的有限元模型 21 图3－4 铸铁的B-H曲线 22 图3－5 初级铁芯的B-H曲线 22 图3－6 网格分布（局部） 22 图3－7 永磁直线电机磁力线图 22 图3－8 永磁直线电机磁通密度矢量图 22 图3－9 永磁直线电机磁场强度矢量图 22 图3－10 永磁电机左边端效应 23 图3－11 永磁电机右边端效应 23 图3－12 永磁体单独作用的磁力线图 23 图3－13 额定电流单独作用下的磁力线图 23 图3－14 合成磁力线图 23 图4－1 矢量变换 26 图4－2 推力随电流大小变化曲线 27 图4－3 法向力随电流大小变化曲线 27 图4－4 推力随电流相位变化曲线图 28 图4－5 法向力随电流相位变化曲线 28 图4－6 初级移动的过程示意图 28 图4－7 推力波动 29 图4－8 法向力波动 29 图4－9 额定电流下推力波动 31 图4－10 额定电流下法向力波动 32 第一章 绪论 1.1直线电机的发展与应用 直线电机是一种将电能直接转化成直线运动机械能而不需要任何中间转换机构的传动装置。在结构上, 直线电机相当于旋转电机在顶上沿径向剖开并将其圆周拉直, 工作原理同旋转电机相似。因为它不需要任何中间环节就可直接驱动被控制对象产生直线轨迹运动, 所以我们把直线电机的传动方式称为“零传动”。 由于采用了“零传动”, 从而较传统的“旋转电机+ 机械变换环节”的传动方式有明显的优势, 如结构简单、无接触、无磨损、噪声低、速度快、精度高等。近年, 随着工业加工质量和运动定位精度等要求的提高, 直线电机受到了广泛的关注。在国外,直线电机驱动技术已进入工业化阶段, 但国内尚处于起步阶段。 近几年, 直线电机研究开发很快, 其结构设计、控制理论和电机原理等方面都有了很大的进展, 并已应用到国民经济的多个领域, 很多机械的直线运动均由直线电机完成，它是一种很有发展前途的新型电机。 直线电机的历史最早可追溯到1840 年惠斯登提出和制作雏形但不成功的直线电机, 至今已有160多年。其发展大致可分为三个阶段: 探索实验阶段(1840—1955)、开发应用阶段(1956—1970) 和实用商品化阶段(1971—)。 1．探索实验阶段(1840—1955) 从1840年到1955年的116年期间，直线电机从设想到实验到部分实验性应用，经历了一个不断探索，屡遭失败的过程。由于当时的制造技术，工程材料以及控制技术的现状，在经过断断续续20多年的顽强努力后，最终却未能获得成功。 2．开发应用阶段(1956—1970) 自1955年以来，直线电机进入了全面的开发阶段。这个时期，它可以说是直线电机的“文艺复兴时期”。特别是该时期的控制技术和材料的惊人发展，更加助长了这种势头。这种增长一方面固然是控制技术和材料促进了这种势头，另一方面，前些年的理论探讨和实验研究工作作为多年的技术贮备起到了重要的作用。 3．实用商品化阶段(1971—) 从1971年始到目前的这个阶段，直线电机终于进入了独立的应用时代，在这个时代，各类直线电机的应用得到了迅速的推广，制成了许多具有实用价值的装置和产品。 目前，直线电机在世界各国的应用大致可分为五个方面，即物流系统、工业设备、信息与自动化系统、交通与民用、军事及其它。 1.2 直线电机的分类 根据不同的 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 ,直线电机有不同的分类型式。直线电机按其结构型式主要可分为扁平型、圆筒型、圆盘型和圆弧型四种。 扁平型直线电机是一种扁平的矩形结构的直线电机，它有单边型和双边型之分。每种型式下又分别有短初级长次级或长初级短次级。圆筒型直线电机为一种外形如旋转电机的圆柱形的直线电机。这种直线电机一般均为短初级长次级型式。在需要的场合，我们还将这种电机做成既有旋转运动又有直线运动的旋转直线电机，至于旋转直线的运动体既可以是初级，也可以是次级。圆盘型直线电机，即该电机的次级是一个圆盘。其初级可以是单边型也可以是双边型。圆弧型电机，它的运动型式是旋转运动，且与普通旋转电机非常接近，然而它与旋转电机相比也是具有如圆盘型直线电机那样的优点，圆弧型与圆盘型的主要区别，在于次级的型式和初级对次级的驱动点有所不同。 直线电机，特别是直线感应电动机，按其功能用途主要可分为力电机，功电机和能电机。 1．力电机 力电机是指单位输入功率所能产生的推力，或单位体积所能产生的推力，主要用于在静止物体上或低速的设备上施加一定的推力的直线电机。它以短时运行、低速运行为主，例如阀门的开闭、门窗的移动、机械手的操作、推车等等。这种电机效率效低，甚至为零(如对静止物体上施加推力时，效率为零)，因此，对这类电机不能用效率这个指标去衡量它，而是用椎力／功率的比来衡量，即在一定的电磁推力下，其输入的功率越小，则说明其性能越好。 2．功电机 功电机主要作为长期进续运行的直线电机，它的性能衡量的指标与旋转电机基在一样，即可用效率、功率因数等指标来衡量其电机性能的优劣。例如高速磁悬浮列车用直线电机，各种高速运行的输送线等等。 3．能电机 能电机是指运动构件在短时间内所能产生的极高能量的驱动电机，它上要是在短时间、短距离内提供巨大的直线运动能，例如导弹、鱼雷的发射，飞机的起飞及冲击、碰撞等试验机的驱动等等。这类直线电机的主要性能指标是能效率(能效率＝输出的动能／电源所提供的电能)。 直线电机按其工作原理可分为两个大的方面，即直线电动机和直线驱动器。如图1-1所示。直线电动机包括交流直线感应电动机、交流直线同步电动机、直线直流电动机和直线步进(脉冲)电动机、混合式直线电动机等。直线驱动器包括直线振荡电动机、直线电磁螺线管电动机、直线电磁泵、直线超声波电功机等。 图1-1 直线电机按工作原理的分类 1.3 直线电机的工作原理 直线电机是一种将电能直接转换成直线运动机械能、而不需要任何中间转换机构的传动装置。它可以看成是将一台旋转电机沿径向剖开，并展成平面而成，如图1-2所示。由定子演变而来的一侧称为初级，由转子演变而来的一侧称为次级。在实际应用时，将初级和次级制造成不同的长度，以保证在所需行程范围内初级与次级之间的耦合保持不变。直线电机可以是短初级长次级，也可以是长初级短次级。考虑到制造成本、运行费用，目前一般均采用短初级长次级。 将图1－3a所示的旋转电机在顶上沿径向展开，并将圆周拉直，变得到了图1－3b所示的直线电机。在这台直线电机的三相绕组中通入三相对称正弦电流后，也会产生气隙磁场。当不考虑绕组中通入三相对称正弦电流后，也会产生气隙磁场。当不考虑由于铁心两端开断而引起的纵向边端效应时，这个气隙磁场的分布情况与旋转电机的相似，即可看成沿展开的直线方向呈正弦形分布。当三相电流随时间变化时，气隙磁场将按A，B，C相序沿直线移动。这个原理与旋转电机的相似，两者的差异是：这个磁场是平移的，而不是旋转的，因此称为行波磁场。 图1－3a旋转电机的工作原理 图1－3b直线电机的基本工作原理 1－定子 2－转子 3－磁场方向 1－初级 2－次级 3－行波磁场 显然，行波磁场的移动速度与旋转磁场在定子内圆表面上的线速度是一样的，即为 （m/s），称为同步速度，且 (1.1) 1.4 有限元分析及相关软件的简介 1.4.1有限元分析简介 目前在工程技术领域内常用到的数值模拟方法有:有限单元法,边界元法,离散单元法和有限差分法,但就其实用性和应用的广泛性而言,主要还是有限单元法。有限单元法的基本思想是将问题的求解域划分为一系列单元，单元之间仅靠节点连接。单元内部点的待求量可由单元节点量通过选定的函数关系插值求得。由于单元形状简单，易于由平衡关系或能量关系建立节点量之间的方程式，然后将各个单元方程“组集”在一起而形成总体代数方程组，计入边界条件后即可对方程组求解。单元划分越细。计算结果就越精确。 有限单元法的基本思想早在40年代初期就由人提出，但真正用于工程中则是在电子计算机出现后，“有限单元法”这个名称是1960年美国的克拉夫（Clough.R.W）在一篇题为“平面应力分析的有限单元法”论文中首先使用的。40年来，有限单元法的应用已由弹性力学平面问题扩展到空间问题,板壳问题，由静力平衡问题扩展到稳定性问题，动力问题和波动问题，分析的对象从弹性材料扩展到塑性，粘弹性，粘塑性和复合材料等，从固体力学扩展到流体力学、传热学、电磁学等领域。 图1－4和图1－5是用有限元分析法对直齿圆柱齿轮的轮齿进行的变形和应力分析，其中图1－4为有限元模型，图1－5是最大切应力等应力线图。在图1－5中采用八节点四边形等参单元把齿轮划分成网络，这些网络称为单元。网格间相互联接的交点称为节点，网格与网格的交界线称为边界。显然，图中的节点数是有限的，单元数目也是有限的，所以称为“有限单元”。这就是“有限元”一词的由来。 图1－4 有限元模型图 图1－5 最大切应力等应力线图 有限元分析计算的思路和作法可归纳如下： 1．物体离散化 将某个工程结构离散为由各种单元组成的计算模型。离散后单元与单元之间利用单元的节点互相连接起来；单元节点的设置、性质、数目等应视问题的性质，描述变形形态的需要和计算精度而定（一般情况，单元划分越细则描述变形情况越精确，即越接近实际变形，但计算量越大）。所以有限元法中分析的结构已不是原来的物体或结构物，而是同样材料的由众多单元以一定方式连接成的离散物体。这样，用有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理，则所获得的结果就与实际情况想符合。 2．单元特性分析 1）选择位移模式 在有限元法中，选择节点位移作为基本未知量时称为位移法；选择节点力作为基本未知量时称为力法；取一部分节点力和一部分节点位移作为基本未知量时称为混合法，位移法易实现计算机自动化，所以在有限元法中位移法应用范围最广。 当采用位移法时，物体或结构物离散化之后，就可把单元中的一些物理量如位移，应变和应力等由节点位移来表示。这时可以对单元中位移的分布采用一些能逼近原函数的近似函数予以描述。通常，有限元法中我们就将位移表示为坐标变量的简单函数。这种函数称为位移式或位移函数，如 ，其中 是待定系数， 是与坐标有关的某种函数。 2）分析单元的力学性质 根据单元的材料性质、形状、尺寸、节点数目、位置及其含义等，找出单元节点力和节点位移的关系式，这是单元分析中的关键一步。此时需要应用弹性力学中的几何方程和物理方程来建立力和位移的方程式，从而导出单元刚度矩阵，这是有限元法的基本步骤之一。 3）计算等效节点力 物体离散化后，假定力是通过节点从一个单元传递到另一个单元。但是，对于实际的连续体，力是从单元的公共边界传递到另一个单元中去的。因而，这种作用再单元边界上的表面力，体积力或集中力都需要等效地移到节点上去。也就是用等效的节点力来替代所有作用在单元上的力。 3．单元组集 利用结构力的平衡条件和边界条件把各个单元按原来的结构重新连接起来，形成整体的有限元方程 (1.2) 式中，K是整体结构的刚度矩阵：q是节点位移列阵;f是载荷列阵。 4．求解未知节点位移 解有限元方程式(1-3)得出位移。这里可以根据方程组的具体特点来选择合适的计算方法。通过上述分析，可以看出，有限元法的基本思想是“一分一合”，分是为了进行单元分析，合则是为了对整体结构进行综合分析。 1.4.2 通用有限元软件ANSYS 数值模拟技术通过计算机程序在工程中得到广泛的应用。到80年代初期，国际上较大型的面向工程的有限元通用程序达到几百种，其中著名的有：ANSYS, NASTRAN, ASKA, ADINA, SAP等。他们多采用FORTRAN语言编写，规模达几万条甚至几十万条语句，其功能越来越完善，不仅包含多种条件下的有限元分析程序而且带有功能强大的前处理和后处理程序。由于有限元通用程序使用方便，计算精度高，其计算结果已成为各类工业产品设计和性能分析的可靠依据。以ANSYS为代表的工程数值模拟软件，既有限元分析软件，不断吸取计算方法和计算机技术的最新进展，将有限元分析，计算机图形学和优化技术相结合，已成为解决现代工程学问题必不可少的有力工具。 ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件。由世界上最大的有限元分析软件公司之一的美国ANSYS开发，它能与多数CAD软件接口，实现数据的共享和交换，如Pro/Engineer，NASTRAN，Alogor，I－DEAS，AutoCAD等，是现代产品设计中的高级CAD工具之一。软件主要包括三个部分：前处理模块，分析计算模块和后处理模块。前处理模块提供了一个强大的实体建模及网格划分工具，用户可以方便地构造有限元模型；分析计算模块包括结构分析（可进行线性分析、非线性分析和高度非线性分析）、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析，可模拟多种物理介质的相互作用，具有灵敏度分析及优化分析能力；后处理模块可将计算结果以彩色等值线显示、梯度显示、矢量显示、粒子流迹显示、立体切片显示、透明及半透明显示（可看到结构内部）等图形方式显示出来，也可将计算结果以图表、曲线形式显示或输出。软件提供了100种以上的单元类型，用来模拟工程中的各种结构和材料。 1.5 课题的主要意义和论文的主要工作 1.5.1 课题的意义 今年来随着工业的发展.新型的高速、高精密、大行程、大推力的进给系统的需求越来越多.而传统的以“旋转电机 + 机械变换环节”方式进行的进给系统已不能满足现在工业的要求,新型的永磁直线电机具有的“零传动”的良好性能和驱动技术的成熟,使其得到了越来越多的关注。为了满足工程领域对新型高性能直线电动机不断增加的需求，许多新的技术和设计方法开始在直线电动机中得到应用。通过计算机辅助，采用有限元分析方法对直线电机进行设计具有良好的发展前景。 用传统的解析法对永磁直线电机进行设计和分析，需要简化很多条件，忽略很多实际因素的影响，往往不够精确，有限元分析作为近年来发展迅速的数值分析方法，在电磁场方面，能够精确的计算电机的各种参数，为此本课题选择ANSYS软件对永磁直线电机进行有限元的分析，希望通过对永磁直线电机进行有限元建模，来近一步提高对永磁直线电机的设计的研究。 1.5.2 论文的主要工作 论文的主要工作如下: 1．建立永磁直线电机样机在ANSYS8.0上的有限元模型。 2．利用ANSYS8.0 对2D有限元模型的静态磁场分析，给出2D磁力线图等数据 3．利用ANSYS8.0对2D有限元模型的力场进行分析计算，得到直线电机的推力及法向力，推力系数，直线电机的波动力等数据 第二章 ANSYS 8.0及其在电磁场有限元计算中的应用 2.1 电磁场的有限单元法 2.1.1 基本方程 1.麦克斯韦方程组 著名的麦克斯韦方程组是研究一切宏观电磁场问题的基础，也是电磁场有限元分析的依据和出发点。其微分形式如下： (2.1) (2.2) (2.3) (2.4) 式中，H——磁场强度矢量；B——磁通密度矢量；E——电场强度矢量；D——电位移矢量；J——传导电流密度矢量； ——自由电荷体密度。 在电磁场问题中的各向同性媒质中本构方程为 (2.5) 式中，ε——介电常数；υ——磁导率；σ——电导率。 在线性均匀及各向同性的媒质中，ε、υ和σ为常数。以稳定磁场为例，在这里所谓线性，就是指媒质中各点磁通密度B的大小与磁场强度H的大小成正比，所谓均匀，就是指媒质的组成情况处处相同，各点的导磁性能也一样；所谓各向同性，就是指沿着空间不同方向，媒质的导磁性能相同，因此磁通密度矢量与磁场强度矢量在空间有着同一方向；同范围，这是只讨论线性，均匀，各向同性媒质。 2.1.2 位函数及边界条件 2.1.2.1 位函数的微分方程 麦克斯韦方程组是场矢量之间的关系表达式，如果直接用来求解电磁场问题，在数学上存在较大困难。因此在分析电磁场问题时，常常引入一定的位函数作为求解的辅助量。 （1）标量电位的偏微分方程 静电场是无旋场,电场强度矢量的旋度处处为零,而对于任一标量函数,其梯度的旋度恒为零,因此在静电场中可以引入标量电位作为待求量 (2.6) 式中， 为标量电位。 方程（2.6）中的负号表示电场强度矢量的方向总是指向电位减小率最大的方向，将该式代入静电场的基本方程式中，即可导出标量电位满足的偏微分方程，这是一个泊松方程 (2.7) （2）标量磁位的偏微分方程 由麦克斯韦方程式（2.1－2.4）可知，在稳定磁场的无电流区域，磁场强度矢量的旋度为零。这时可以引入标量磁位作为待求量 (2.8) 式中， 为标量磁位。 将式（2.9）代入稳定磁场的基本方程式中，注意到这时电流密度矢量J 0，即可导出标量磁位满足的偏微分方程，这是一个拉普拉斯方程 (2.9) （3）矢量磁位的偏微分方程 在稳定磁场的有电流区域，磁场强度矢量的旋度不为零，因此不能采用标量磁位进行求解。但考虑到磁通密度矢量的散度恒为零，而对于任一矢量函数，其旋度的散度也是恒为零。因此可以引入矢量磁位来描述场域中有电流存在时的稳定磁场问题。 (2.10) 式中，A为矢量磁位 在平面磁场中，电流密度矢量J与矢量磁位A沿着Z轴方向，分别只有一个分量Jc与Az。在平面xoy上，Jc与Az是坐标x、y的函数。将上式代入稳定磁场的基本方程式中，得矢量磁位满足的偏微分方程为 (2.11) 从以上简明的推导可以看出，标量电位和矢量磁位满足泊松方程，标量磁位满足拉普拉斯方程，这些方程统称为泛定方程。 2.1.2.2 位函数的边界条件 边界条件是求解电磁场问题的关键。由于电磁场问题的复杂性，在很多情况下表现为边界条件的多变性。不同的问题，有不同的边界条件；同一问题在不同的情况下，也有不同的边界条件。在平面电磁场问题中，设求解区域为 ，它的边界为 ，计算所采用的变量为，边界条件通常分两种不同情况给出。 第一类边界条件为：在边界 上满足已知物理量 （2.12） 对于平面稳定电磁场问题，第一类边界条件用位函数给出有以下三种形式： 标量电位 (2.13) 标量磁位 (2.14) 矢量磁位 (2.15) 第二类边界条件为：满足物理量 在边界 上的法向导数 (2.16) 对于平面稳定电磁场问题，第二类边界条件用位函数给出，有以下三种形式： 标量电位： (2.17) 标量磁位： (2.18) 矢量磁位： (2.19) 式中 ——电位移矢量的法向分量； ——磁通密度矢量的法向分量； ——磁场强度矢量的切向分量； ——磁阻率，即磁导率 的倒数。 2.1.3 有限元法的求解 现以典型的三角形单元为例说明电磁场问题求解基本思路 1.单元分析 将求解域 离散成 个三角形单元，单元的三个顶点为i,j,m,选取单元位移函数 (2.20) 如前所述可得到单元内位移函数的表达式 （ ） (2.21) 在求解式中，总的能量泛函为单个能量泛函之和，对其中的面积分式和线积分式分别进行离散化处理。将线性插值函数代入面积分式就可以得到经过离散化后单元e的能量函数表达式，将该式对单元中每一顶点的位函数u求一阶偏导数。 2.整体分析 整个求解域的能量函数由每个单元的能量函数叠加而成， (2.22) 令上式为零，并代入单元分析中的结果中，就可以得到当能量函数达到极值时，位函数必须满足的矩阵方程：Ku＝R 最后，引入加强边界条件， ，求解。 2.2 电磁力计算方法 有限元计算电磁力的方法主要有三种：安培力定律、麦克斯韦张力法和、虚功原理法。以上方法都可以计算一个物体的总体受力情况，但麦克斯韦张力法和虚功原理法不能计算物体的受力分布。安培力定律只能应用于非磁性导体，能计算物体的受力分布。 在ANSYS有限元分析中，经常使用的有2种电磁力的计算方法，即电磁场的虚功法和MAXWELL张量法。 2.2.1 电磁场的虚功法 由于有限元的出发点是划分区域并使每个单元的能量泛函达到最小值,所以,通过虚位移使磁能变化而求力的虚功原理法非常合适有限元法分析.磁能使全局量,虚功原理不易受到不合适划分引起的局部误差的影响.首先,对给出的电流源和各种材料,用有限元法求解并计算出磁共能 。磁共能 可由下式求得： (2.23) (2.24) 然后，我们把需要计算力的单元移动微小的距离 ，在同样的电流下，重新用有限元法求解并计算出新位置的磁共能。 最后，得到物体受到的电磁力为 (2.25) 2.2.2 电磁场的MAXWELL张量法 根据麦克斯为的理论，作用于任意区域上的体积力可归为这一区域表面所受的张力。 首先，铁磁材料v能用面电流 和体电流 的分布替代，那么有： = (2.26) = (2.27) 其中 表示面积电位矢量的切向分量，下标t表示矢量的切向分量，磁矢量被定义为 (2.28) 通过以上公式，我们在任意区域（从理论上）都可以用面和体电流等效为磁性材料。这种思想与等效电流替代永磁体是一致的。然后，利用安培力定律，得力的密度 为: (2.29) (2.30) (2.31) 从式2－8可以得到 的微分形式 = (2.32) 其中T 为： (2.33) 则物体所受的力为: (2.34) 利用矢量微分理论，可以把T的体积分转化为面积分得： (2.35) 我们可以沿着面积的法向和切向方向对T进行积分，定义面积单位矢量的切向分量和法向分量为： (2.36) (2.37) 磁通密度的切向分量和法向分量为： （2.38) (2.39) 最后可得到用麦克斯伟张力法求解的公式，力的切向分量增量 和法向分量增量 ： (2.40) (2.41) 2.3 ANSYS有限元建模 利用ANSYS8.0建立永磁直线电机样机的有限元模型，具体流程如下： 1.设置有限元分析的物理环境，其中包括各种材料的性质的选定，定义模型的何尺寸及载荷的参数等 首先定义分析的物理环境，在进入程序前选择ANSYS多物理场模块，在进入程序后，改变任务题目对话框，输入名称：2-D model GUI: Utility Menu>File>Change Title>输入题目2-Dmodel>按OK 然后，定义单元类型，本次建模选择ANSYS提供的PLANE53单元。 GUI：Main Menu>preprocessor>Element Type>Add/Edit/Delete>按Add>选择Magnetic Vector and Vect Quad 8nod53(PLANE53)> 按close 设置计算时的电磁单位制是“米－安培—秒”制 GUI：Main Menu>Preprocessor>Material Props>Electromag Units>按OK,设置单位制为“MKS System” 定义材料特性：根据需要，定义10个材料，具体如表2－1所示。 GUI：Main Menu>Preprpcessor>Material Props>Constant>Isotropic>选择材料号1，再按OK>在“Relative permeability(MURX)”域中输入1>按Apply>新建材料号2-10，重复上述步骤，定义10个材料特性 表2－1材料特性 1号材料 空气 磁导率 1 2号材料 硅钢 B-H曲线 emagM54.SI_MPL 3号材料 线圈1 磁导率 1 4号材料 线圈2 磁导率 1 5号材料 线圈3 磁导率 1 6号材料 线圈4 磁导率 1 7号材料 线圈5 磁导率 1 8号材料 线圈6 磁导率 1 9号材料 永磁体 磁导率 矫顽力 1 804000 Oe 10号材料 铸铁 B-H曲线 emagSilicon.SI_MPL 其中Relative permeability（MURX）代表磁导率，Coercive force（MGXX）代表矫顽力。B-H曲线通过 Main Menu>Preprpcessor>Material Props>Read form File从ANSYS材料库中导入。 2 建立有限元分析模型，包括建立永磁直线电机的几何模型，给各区域赋于材料属性，对各区域进行划分及细化。 现在在ANSYS8.0中创建永磁直线电机实体模型，先创建电机铁轭的2－D实体模型，长度为510 mm，高度10 mm的矩形，在创建3－D实体模型时需要它的宽度参数为50 mm。 GUI>Preprpcessor>Modeling>Create>Areas>Rectangle>By 2 corners》在弹出框中输入 WPx = 0.255，Wpy = 0 ，Width = 0.51，Height = -0.01>按Apply，由此建立模型所需的所有平面。 然后在ANSYS中取布尔运算，使所以建立的面连成一个实体。 GUI>Main Menu>Preprpcessor>Modeling>Operate>Booleans>overlap>Areas>按pick all 到这里,直线电动机的2-D实体模型就建好了,如图2－2所示，图2-3显示了该电机的3-D实体模型 图2-1 永磁直线电机的2-D实体模型 图2-2 永磁直线电机的3-D实体模型 现在我们对电机2-D实体模型进行划分网格,建立它的有限元模型。划分网格前，先赋予实体模型材料属性。在定义永磁体的时候，需要特别定义它的NS极的方向。 GUI：Main Menu>Preprocessor>-Meshing>Mesh Attributes>Picked Areas >选择需要定义的面>选择单元和材料号。 图2－2显示了定义完材料属性后的实体模型,其中各种材料用不同的颜色来区分。 图2-3 定义完材料属性后的实体模型 定义完以后，就可以进行网格的划分了，此时需要考虑网格的划分类型和细化程度。在这里，选择细化程度为1（fine），采用自由网格划分方式。 GUI:MainMenu>Preprocessor>-Meshing-SizeCntrls>-SmartSize-Basic> GUI:Main Menu>Preprocessor>-Meshing->Mesh>-Areas-free>Pick All. 图2-5显示了永磁直线电机的网格划分的有限元模型。 图2-4 永磁直线电机的网格划分 3 加载荷和边界条件 在对永磁直线电机建模分析中，我们在模型上下边选择第一类边界条件，在左右边选择第二类边界条件。 GUI:MainMenu>Solution>-Loads->Apply>-Magnetic->Boundary>-Vector Porten-Flux Par’l-on lines> 为了求解推力和法向力，需要将永磁直线电机初级定义为一个组件，然后加上力标志。 GUI：Utility Menu>Select>Entities> Utility Menu>Select>Comp/Assembly>Creat Componet 图2－5 加载力标志和边界条件 在对永磁直线电机进行分析中，有时候要在线圈中加上电流载荷，它是以电流密度的形式赋予线圈的。如图2-6所示。 GUI MainMenu>Solution>-Loads->Apply>-Magnetic->Excitation>-Curr Density-On Elements, 图2-6 永磁直线电机加载后的示意图 现在，我们已经建立了一个永磁直线电机的有限元模型。接下去的工作就是对模型进行求解。求解需要注意分析类型和求解器的选择。具体可以参考附录中ANSYS8.0电磁场分析的使用说明。 2.4 网格的划分和单元的选择 2.4.1网格的划分 进行有限元分析，首先必须进行网格的划分，是有限元处理中主要工作，也是整个有限元分析的关键工作，它将集合模型转化为节点和单元构成的有限元模型。网格划分的质量和优劣将会对计算结果产生相当大的影响，其好坏直接影响到计算的准确性和计算进度。划分不合理甚至会导致计算不收敛。 为了建立正确，合理的有限元模型，在划分网格是应考虑以下基本原则： 1．网格数量 网格数量多少影响计算结果的精度和计算规模的大小，但要考虑计算时间的影响。在决定网格数量时也应考虑分析数据的类型。 2．网格疏密 网格疏密时指在结构不同部位采用大小不同的网格，这也是为了适应计算数据的分布特点。在计算数据变化梯度较大的部位（如应力集中处），为了较好地反映数据变化规律，需要采用比较密集的网格。 3．单元阶次 许多单元都具有线性，二次和三次等形式，其中二次和三次形式的单元称为高阶单元。选用高阶单元可提高计算精度，因为高阶单元的曲线或曲面边界能够更好的逼近结构的曲线和曲面边界，且高阶单元的节点数较多，在网格数量相同的情况下由高阶单元组成的模型规模要大的多，因此需要考虑计算精度和时间。 其他的原则还包括网格质量，网格分界面和位移协调性。 因为建立的永磁直线电机有限元的实体模型几何结构比较清楚，我们将采用自由划分网格的方式对永磁直线电机进行划分。 需要特别指出的是，在划分网格的疏密程度上，考虑到气隙层比较薄，太疏的的网格不能准确划分，因此我们对网格划分的疏密程度做了一定的要求。 2.4.2单元的选择 在ANSYS8.0的 2-D 静态磁场分析中，可供选择的单元类型如下表所列： 表2－2 可选择单元类型 类型 维数 形状 自由度 2-D实体单元 PLANE13 2－D 四变形，4节点 三角形，3节点 磁矢势（AZ），位移，温度或时间积分电势 PLANE53 2－D 四变形，8节点 三角形，3节点 磁矢势（AZ），位移，温度或时间积分电势 远场单元 INFIN9 2－D 线性 2节点 磁矢势（AZ）， INFIN10 2－D 四变形，4节点 或8节点 磁矢势（AZ），位移，温度 通用电路单元 CIRCU124 无 通用电路单元，6节点。 磁矢势（AZ）位移，温度或时间积分电势 本次设计建立的是永磁直线电机的2-D有限元模型，属于二维域的分析，一个二维区域可以用三角形或是四变形单元进行离散，同常，矩形单元最适合离散 规则 编码规则下载淘宝规则下载天猫规则下载麻将竞赛规则pdf麻将竞赛规则pdf 区域，三角形单元可用于离散不规则区域，如图2－1所示，建立的永磁直线电机的实体模型具有规则的区域，所以这次建模选定采用ANSYS8.0提供的PLANE53单元，它是四边形8节点的单位，具有很高的计算准确性。 PLANE53 单元说明 PLANE53 用于 2 维 (平面和轴对称) 磁场问题的建模。本单元有 8 个节点，每个节点最多 4 个自由度：磁矢量势的 z 分量 (AZ)、时间积分电标量势  (VOLT)、电流 (CURR) 和电动势降 (EMF)。 PLANE53 是以磁矢量势理论为基础的，可以用于以下低频磁场分析：静磁、涡流 (AC 时间谐波和瞬态分析)、电动力磁场 (voltage forced magnetic fields)(静态, AC时间谐波和瞬态分析) 以及电磁-电路耦合场 (静态，AC时间谐波和瞬态分析)。本单元具有非线性磁能力，可用于 B-H 曲线或永久磁体退磁曲线的建模。 图2-7 PLANE53单元 2.5 边界条件的加载 本次设计采用的2-D电磁场的分析，假定永磁直线电机的磁场沿Z轴方向，因此需要加载磁力线边界条件AZ＝0，使磁力线强制平行于永磁直线电机的上下表面。 因为需要分析电机的力场特性，所以需要加载力标志。需要注意的是要计算力的部分周围要包围一层空气单元。 在ANSYS中电流是通过电流密度的方式加载的。 图2－8显示了永磁直线电机边界条件加载的情况。上下边约束第一别界条件，左右边用第2边界条件。 图2－8 边界条件示意图 第三章 永磁直线电机的电磁场有限元仿真分析 3.1永磁直线电机的磁场分析方法概述 与其他电励磁直线电机不同，永磁直线电机采用永磁体励磁产生励磁磁场。永磁体在电机中即是磁源，又是磁路的组成部分。因为永磁体制造工艺，电机其他材料和磁路结构复杂，使永磁直线电磁计算的复杂性增加，计算结果的准确性降低。 有时为了简化分析计算，目前在许多工程问题中仍常采用场化路的方法，将空间实际存在的不均匀分布的磁场转化为转化成等效的多段磁路来进行分析。由于直线电机与传统直线电机在结构上存在较大的差别，加之其饱和特性，以路的观点进行电机性能理论的分析就又很大的局限性。相反，利用电机电磁场理论和有限元法进行直线电机电磁场分析与计算，以场的观点，全面，系统地分析电机的性能，以便进行电机的设计，性能分析及仿真计算，具有明显的优势。 电机电磁场数值分析主要采用有限元法、边界元法、有限差分法，其中，最有效，用途最广泛的时有限元法，与其他方法相比，有限元法具有以下优点： 1．处理第二类边界条件和内媒质边界条件非常方便，对于由多种材料组成，内部具有较多媒质分界面的电机电磁场来说，非常实用。 2．几何划分灵活，适合解决电机这类几何形状复杂的问题。 3．可较好的处理非线性问题。 对于永磁直线电机，存在永磁材料和非线性铁磁材料，结构较为复杂，采用普通电机分析方法很难得到准确的结果。而采用有限元分析方法进行研究，如果能建立适当的有限元模型，就能够得到准确的结果，如推力，感应电势等。 3.2永磁直线电机的磁场有限元仿真模型 3.2.1永磁直线电机的结构参数 本研究采用的永磁直线电机为短初级型单边结构。样机的相关参数如表3-1所示，图3-1为永磁直线电机的结构图。图3-2为绕组连接示意。 图3-1 永磁直线电机的结构图 根据图3-2，将电机线圈分成以下几组，以方便电流的加载。 A相：1,3,14,16 材料号3 Phase Angel 0度 2,4,13,15 材料号4 Phase Angel 180度 B相：9,11,22,24 材料号5 Phase Angel -120度 10,12,21,23 材料号6 Phase Angel 60度 C相：6,8,17,19 材料号7 Phase Angel 120度 5,7,18,20 材料号8 Phase Angel 300度 表3－1 永磁直线电机参数 项 目 符号 数值 初级 相数 m 3相 级数 p 6极 每相匝枢 N 160 槽数 Ns 24 齿距 (s 14.67mm 槽口宽度 s 6.67mm 槽深 d 29.2mm 铁心高度 I 43.2mm 铁心长度 L 360mm 次级 极距 16mm