二维时谐分析 Ansys工程电磁场有限元分析 华科电气
基于ANSYS软件的电机电磁场有限元分析解读
基于ANSYS软件的电机电磁场有限元分析发表时间：2007-9-11 作者: 黄劭刚夏永洪张景明来源: 万方数据关键字: APDL语言同步发电机电磁场有限元介绍了应用ANSYS自带的APDL编程语言进行软件开发，将该软件应用于同步发电机空载磁场分析中，在电机的电磁场计算中实现了电机的自动旋转、自动施加载荷的功能，使用、修改方便，并且计算速度快。
通过对电磁场计算结果的后处理，得出了同步发电机的旋转磁场波形和电压波形。
样机测试结果验证了分析结果的正确。
1 前言ANSYS软件是一个功能强大、灵活的，融结构、热、流体、电磁、声学于一体的大型通用有限元分析软件。
广泛用于核工业、石油化工、航空航天、机械制造、土木工程等一般工业及科学研究领域的设计分析。
在实际的电机电磁场分析中，电机的转子磁极形状、定子齿槽形状、气隙大小以及铁磁材料均已确定，但是当转子相对十定子齿槽的位置不同时一，其计算结果也不相同。
为了分析电机电磁场问题，若把定、转子相对位置固定不变进行求解，再对电磁场计算结果进行傅立叶级数分解来计算电机绕组的电势则误差太大。
为此，需要对定、转子不同位置时一分别进行计算，然后通过电磁场的计算结果求出电机何个定子齿部磁通随转角变化的关系，然后根据磁通的变化率求出电机基波绕组的电势。
ANSYS软件是目前应用最为广泛、使用最方便的通用有限元分析软件之一，应用ANSYS软件来分析电机电磁场是非常有效的。
但是当采用ANSYS软件的图形用户界面( GUI)操作方式时，每次定、转子之间的旋转、网格剖分、施加载荷进行求解、查看计算结果等都需要人工进行重复操作，使用起来非常繁琐，并且效率低。
为此，木文采用ANSYS软件的APDL语言编写的软件对同步发电机的空载磁场进行研究，实现了电机定、转子之间的自动旋转，自动网格剖分，自动施加载荷以及自动求解的功能。
整个电磁场分析过程无需人工进行干预，使用方便，便于修改，并且大大提高了计算速度。
Ansoft Maxwell 2D 工程电磁场有限元分析
1.3.1Maxwell 2D 的边界条件
静磁场有以下几种边界条件： 矢量磁位边界条件 对称边界条件 气球边界条件 主边界条件 从边界条件 • 1.3.1.1 Default Boundary Conditions 自然边界条件，也称纽曼边界条件，可以用来描 述两个相接触的物体，在接触面上，磁场强度H 的切向分量和磁感应强度B 的法向分 量保持连续。
Maxwell 2D 基础
1.3 Maxwell 2D 的边界条件和激励源
边界条件和激励源方式按照不同的求解器来设定。 按照计算模型所需的求解器不同，主要可以分为以下 6 大类： 求解器 可计算的执行参数 静磁场 矩阵（电容）、力、转矩 涡流场 矩阵（电感）、力、转矩、磁通量 瞬态磁场 矩阵（阻抗）、力、转矩、磁通量 静电场 导纳、电流 交变电场 电导、电流 直流传导电场 注： 瞬态磁场是指被求解问题随时间做一定有规则的运动， 以及所加载激励是时间、位置、或者速度的函数关系，
软件默认的参数变量为_t，在X、Y、Z 三个方向上都可以设置为_t的函数，而在 Start_t 和End_t 中设置参数_t 的起始和终止范围，通过Points 项可以设置由多少个点 组成该参数曲线，若设置为0 则表示由软件默认的点数组成，此时的曲线较为光滑 ，若该项设置过少则曲线将有多段直线组成。
Maxwell 2D 基础
1.2 Maxwell 2D 的材料管理
1.2.2 常用硅钢片50W600的添加
以硅钢片50W600材料为例，先要了解该材料的特性，找到相应的相对磁化曲线 表，它的磁化曲线是非线性的；电导率在2e+6 S/m 左右。 添加步骤： 1、材料命名：50W600 2、选定坐标系：Cartesian 3、设置相关参数 Relative Permeability --相对磁导率设为 非线性曲线，点击右侧“Bh Curve”，进入 磁化曲线表， 输入相应的数据，点击OK 。 Bulk Conductivity --电导率设为2000000。 后两项默认即可。
23.1 电磁场有限元分析简介_ANSYS 有限元分析从入门到精通_[共3页]
![23.1 电磁场有限元分析简介_ANSYS 有限元分析从入门到精通_[共3页]](https://img.taocdn.com/s3/m/ba975fdbba1aa8114531d91f.png)
第23章
ANSYS电磁场分析
23.1 电磁场有限元分析简介
1．磁场分析对象
利用ANSYS/Emag或ANSYS/Multiphysics模块中的电磁场分析功能，ANSYS可分析计算下列的设备中的电磁场，例如：
●电力发电机；
●磁带及磁盘驱动器；
●变压器；
●波导；
●螺线管传动器；
●谐振腔；
●电动机；
●连接器；
●磁成像系统；
●天线辐射；
●图像显示设备传感器；
●滤波器；
●回旋加速器；
●电解器；
●磁悬浮装置。
ANSYS磁场分析的有限元公式由磁场的Maxwell方程组导出，通过将标量势或矢量势
等引入Maxwell方程组中并考虑其电磁性质关系，开发出适合于有限元分析的方程组。
ANSYS软件的其他一些功能增强了软件的电磁分析能力和灵活性。
例如，用户可方便地选择MKS、CGS或其他一些单位制作为电磁场分析的单位制。
作为标准的Frontal求解器的替代者，。
ANSYS有限元分析二维静态磁场仿真
一周总结报告一、ANSYS学习1.学习情况目前正在边看书籍边操作ANSYS系统，已经了解了ANSYS的基本操作系统以及ANSYS 分析过程的三大步骤，大体上知道了它的整个工作流程。
目前正在深入仔细学习每一部分的详细步骤。
现在已经学习了ANSYS有限元分析典型步骤、实体建模、网格划分、创建有限元模型，正在学习加载和求解这一部分。
2.理论知识(1)网格划分与创建有限元模型①设置单元属性，包括：a.选择单元类型，如常用的有PLANE13，PLANE53，INFIN110；在Element Type中设置；b.设置单元实常数，如线圈横截面积、匝数、导体填充率等；c.设置材料属性，如泊松比、材料密等；d.设置单元坐标系统。
②通过网格划分工具设置网格划分属性包括：a.单元属性分配设置，作用是在网格划分之前为模型(包括实体和有限元模型)分配单元属性；b.智能划分水平控制；c.单元尺寸控制，单元尺寸的意思是单元边的长度。
③实体模型的划分ANSYS有两种方式对实体模型进行网格划分。
映射网格划分方法：最大特点就是必须使用形状规则的单元划分，对于面对象必须使用三角形单元或四边形单元，对于体对象只能使用六面体单元。
故划分对象必须形状规则。
不是任何形状的对象都能用映射网格划分。
(2)加载和求解有限元分析的主要目的在于得到系统在特定激励源和边界条件下的响应。
这些激励以及边界条件统称为载荷。
所以载荷包括边界条件和激励。
磁场分析中常见的载荷有磁势、磁通量边界条件等。
载荷分为六大类：自由度约束、集中力载荷、面载荷、体载荷、惯性载荷以及耦合场载荷。
关于载荷步、子步和平衡迭代，通过阅读理论知识自己的理解的总结是：一个实际加载过程需要多次施加不同的载荷才能满足要求，每一步就称为一个载荷步。
一个载荷步可以通过多个子步来逐渐施加。
平衡迭代用于考虑收敛的非线性分析。
3.仿真结果目前按照教程的步骤将ANSYS从建立模型到加载求解再到查看后处理器的整个分析过程大体操作了一遍，目的就是先通过简单模型熟练ANSYS的整体操作。
电磁应力 Ansys工程电磁场有限元分析 华科电气
1.电磁张Leabharlann 与电磁力：电磁张量定义：
有：
应用散度定理，得到：
式中， 。
上述公式反映了法拉第和Maxwell对电磁力的看法，他们认为，电磁力是靠连续的弹性媒质（即以太）来传递的。后来的研究表明，以太的假设不是必要的，但是上述公式给出了计算电磁力的正确结果。
对于包围在体积V内的静止电荷和恒定电流，有下列关系式：
2.虚功原理：
二维静态磁场分析：
有限元ansys电磁场分析详解PPT课件
• 选择OK
第4页/共33页
• 重复这些步骤，定义定子磁体材料3
• 为转子磁体平行磁化方向定义11号局部坐标系 • 水平方向反时针30度（总体坐标+X 轴） • 局部坐标系原点与总体坐标系一致 Utility>workplane>local coord. systems>create local CS>at specified location
• 选择 OK • 与前面一样重新设置衔铁的关联 • 对除有约束方程的节点外的所有外部节点重新施加平行条件 • 执行求解
第29页/共33页
• 显示磁通密度和磁力线迭加图 – 由于衔铁位置改变，磁力线随 着变化 – 定子内最大磁密BSUM增大 – 模型交界处磁场连续
BSUM (T)
第30页/共33页
谢谢您的观看！
第33页/共33页
2.5-33
第17页/共33页
• 模拟有许多磁极的电机，周期性边界 条件非常有用
• 右图显示的是一个10极永磁电机
• 模拟转子的运动。当转子转动时，电 流会变化。
• 定子槽内显示电流密度
• 本模型也允许转子和定子相互独立
• 观看动画，可执行动画文件:
mach2d.avi
定子
第18页/共33页
转子
约束方程—不相同网格
• 将定子一侧边界上的节点建立组件. • 选择定子模型边界上线段 • 选择STATOR组件 • 再选择边界上线段 • 选择所选线段上的全部节点 • 建立单节点组件CE_N
第22页/共33页
定子内半径 全部节点
• 选择衔铁组件ARMATURE • 选择节点组件 CE_N • 应用约束方程生成器
Preproc>coupling/ceqn>adjacent regions
ANSYS电磁场分析指南
ANSYS电磁场分析指南（共17章）ANSYS电磁场分析指南第一章磁场分析概述：ANSYS电磁场分析指南第二章2-D静态磁场分析：ANSYS电磁场分析指南第三章２-Ｄ谐波（ＡＣ）磁场分析：ANSYS电磁场分析指南第四章２-Ｄ瞬态磁场分析：ANSYS电磁场分析指南第五章３-Ｄ静态磁场分析（标量法）：ANSYS电磁场分析指南第六章３-Ｄ静态磁场分析（棱边元方法）：ANSYS电磁场分析指南第七章３-Ｄ谐波磁场分析（棱边单元法）：ANSYS电磁场分析指南第八章３-D瞬态磁场分析（棱边单元法）：ANSYS电磁场分析指南第九章3-D静态、谐波和瞬态分析（节点法）：ANSYS电磁场分析指南第十章高频电磁场分析：ANSYS电磁场分析指南第十一章磁宏：ANSYS电磁场分析指南第十二章远场单元：ANSYS电磁场分析指南第十三章电场分析：ANSYS电磁场分析指南第十四章静电场分析（h方法）：ANSYS电磁场分析指南第十五章静电场分析（P方法）：ANSYS电磁场分析指南第十六章电路分析：ANSYS电磁场分析指南第十七章其它分析选项和求解方法：第一章磁场分析概述1.1磁场分析对象利用ANSYS/Emag或ANSYS/Multiphysics模块中的电磁场分析功能，ANSYS可分析计算下列的设备中的电磁场，如：·电力发电机·磁带及磁盘驱动器·变压器·波导·螺线管传动器·谐振腔·电动机·连接器·磁成像系统·天线辐射·图像显示设备传感器·滤波器·回旋加速器在一般电磁场分析中关心的典型的物理量为：·磁通密度·能量损耗·磁场强度·磁漏·磁力及磁矩· S-参数·阻抗·品质因子Q·电感·回波损耗·涡流·本征频率存在电流、永磁体和外加场都会激励起需要分析的磁场。
Ansys电磁场分析简介
3、什么时候要用二维分析？
2-D平面分析：
忽略终端效应
模型位于X-Y平面
电流方向只沿X-Y面的 法线方向（Z方向）
磁场只有X-Y面内的分 量
2-D轴对称分析：
模型位于X-Y平面
电流方向只沿X-Y面的 法线方向（圆周Z方向）
磁场只有X-Y面内的分 量
第12页/共22页
4、什么时候要用三维分析？
ANSYS 软件是融结构、热、流体、电磁、声学于一体的大型通用 有限元软件，可广泛的用于核工业、铁道、石油化工、航空航天、机 械制造、能源、汽车交通、国防军工、电子、土木工程、生物医学、 水利、日用家电等一般工业及科学研究。该软件提供了不断改进的功 能清单，具体包括：结构高度非线性分析、电磁分析、计算流体力学 分析、设计优化、接触分析、自适应网格划分及利用ANSYS 参数设计 语言扩展宏命令功能。
6、ANSYS磁场分析的类型
2-D静态磁场分析，分析直流电或永磁体所产生的磁 场，用矢势法。
2-D谐波磁场分析，分析低频交流电流或交流电压所 产生的磁场，用矢势法。
2-D瞬态磁场分析，分析随时间任意变化的电流或外 场所产生的磁场，用矢势法。
第15页/共22页
ANSYS磁场分析的类型（续1）
ANSYS 介绍
ANSYS是一种广泛的商业套装工程分析软件。所谓工程分析软件， 主要是在机械结构系统受到外力负载所出现的反应，例如应力、位移、 温度等，根据该反应可知道机械结构系统受到外力负载后的状态，进 而判断是否符合设计要求。一般机械结构系统的几何结构相当复杂， 受的负载也相当多，理论分析往往无法进行。想要解答，必须先简化 结构，采用数值模拟方法分析。由于计算机行业的发展，相应的软件 也应运而生，ANSYS 软件在工程上应用相当广泛，在机械、电机、土 木、电子及航空等领域的使用，都能达到某种程度的可信度，颇获各 界好评。使用该软件，能够降低设计成本，缩短设计时间。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第三章二维时谐分析（2-D Harmonic (AC) Analysis）分析对象：正弦交变电流产生的效应∙Eddy currents涡流∙Skin effects集肤效应∙Power loss due to eddy currents涡流损耗∙Forces and torque力和力矩∙Impedance and inductance阻抗和自感∙Two contacting bodies with dissimilar meshes不同网孔（如转子/定子气隙）典型应用：∙Transformers变压器∙Induction machines电感器∙Eddy-current braking systems涡流刹车系统∙Most electromagnetic devices that work on AC交流电磁装置不能有永磁体，不考虑磁滞效应。
3.1 线性分析与非线性分析∙严格上讲，时谐分析只适用于线性分析。
∙中等饱和的非线性问题，如果不关心波形畸变，只关心时间平均的量，则可分析。
∙B-H曲线为等效值。
∙严格的非线性分析只能通过瞬态分析完成。
3.2 所用的单元实体单元：远场单元：电路单元：3.3 创建2D 时谐分析物理环境自由度选项∙ AZ ：无外加电压；用于短路导体∙ AZ-VOLT ：允许外加电压，可模拟多种状况t ϕ∂=--∇∂A Ej j V ωω∇=--E ANote ：d V t ϕ=⎰ (time-integrated potential)，（单位：V s ）。
同一断面上V 是常数。
使用时需要对所有相关节点进行耦合。
∙ AZ-CURR ：用于电压源驱动的线圈（线圈不计涡流）模型物理特性设置∙短路导体∙开路导体：耦合电位；无净电流∙载流块导体：耦合电位；有净电流鼠笼转子等∙载压线圈：耦合电流；外加电压∙载流线圈：外加电流，均匀电流密度建模与剖分∙ 集肤效应的考虑集肤深度：δ==每个集肤深度至少划分一层单元。
施加边界条件与载荷∙ 边界条件与载荷可以加在实体（线、面）上，也可以直接加在有限元模型（节点、单元）上。
∙ 周期性边界条件：Command(s): P ERBC2DGUI: Main Menu> Preprocessor> Loads> Define Loads> Apply> Magnetic> Boundary> Vector Poten> Periodic BCs奇对称（odd symmetry）选项表示半个周期对称条件；偶对称（even symmetry）表示全周期对称。
幅值、相位与频率相位（单位为°）：只在多相分析时才有效，单相时令为0。
频率（单位为Hz）：Command(s): H ARFRQGUI: Main Menu> Solution> Load Step Opts> Time/Frequenc> Freq and Substps求解∙定义时谐分析类型Main Menu> Solution> Analysis Type> New Analysis，选择Harmonic analysis.∙定义分析类型（默认Full）Main Menu> Solution> Analysis Options ∙选择方程求解器（2D默认波前求解器）Main Menu> Solution> Analysis Options∙非线性问题迭代选项（默认）Main Menu> Solution> Load Step Opts> Nonlinear> Convergence Crit∙设置频率：Main Menu> Preprocessor> Loads> Load Step Opts> Time/Frequenc> Freq and Substps∙启动求解：线性问题：Main Menu> Solution> Solve> Current LS非线性问题：Command(s): H MAGSOLVGUI: Main Menu> Solution> Solve> Electromagnet> Harmonic Analys> Opt&Solv∙跟踪收敛曲线∙结束求解Command(s): F INISHGUI: Main Menu> Finish查看结果获得的数据:∙Nodal magnetic flux density (BX, BY, BSUM)∙Nodal magnetic field intensity (HX, HY, HSUM)∙Nodal magnetic forces (FMAG: components X, Y, SUM) ∙Nodal Lorentz magnetic forces (FMAG; components X, Y, SUM)∙Nodal reaction current segments (CSGZ)∙Joule heat per unit volume (JHEAT)示例：载压线圈时谐分析命令流：!: 设置环境/batch,list/prep7/title, Voltage-fed thick stranded coil in free space/com, Calculate the current in the coil from an applied AC voltage/com,/com, Parameters/com, n = 500 coil turns/com, s = .02 (coil winding width and depth) meters /com, r = 3*s/2 meters (mean radius of coil)/com,/com,/nopr!: 选择单元et,1,53,,,1 ! airet,2,53,2,,1 ! voltage forced coilet,3,110,,,1 ! far-fieldemunit,mks!:设置材料属性mp,murx,1,1mp,murx,2,1mp,rsvx,2,3.00e-8 ! resistivity of coil!:定义参数s=.02n=500r=3*s/2!: 几何建模rectng,s,2*s,0,s/2pcirc,0,6*s,0,90pcirc,0,12*s,0,90!:布尔操作aovlap,all!: 设置属性、剖分网格asel,s,area,,1aatt,2,1,2asum ! 所选单元的面积求和*get,a,area,,area ! area of 1/2 coil cross-sectionasel,s,area,,4aatt,1,1,3asel,allcsys,1 !: 选择使用柱坐标系lsel,s,loc,x,9*s !: 选取半径为9s的线段lesize,all,,,1esize,,8mshape,0,2d ! mapped mesh with quadsmshkey,1amesh,1,4,3 ! mesh far-field and coilsmrtsize,2esize,s/4mshape,1,2d ! specifytriangle elementsmshkey,0 ! free meshamesh,5 ! mesh air region!: 定义线圈实常数r,1,2*a,500,,1,1 ! coil constants：序号1，面积2a，! 匝数500，电流1A，填充因子1!: 选择线圈单元，耦合电流esel,s,mat,,2 ! get coil elementsnsle,s !: 选择节点cp,1,curr,all ! couple curr dof in coilallsel,all !: 选择全部对象!: 设置无限远边界csys,1nsel,s,loc,x,12*ssf,all,inf ! set infinite surface flag csys,0 !: 换回直角坐标系!: 设置对称轴（y轴）上的边界条件nsel,s,loc,x,0d,all,az,0 ! flux-parallel condition allsel,allfinish!: 建模及剖分结束。
以下进入求解过程/soluantype,harm !: 定义分析类型为谐分析esel,s,mat,,2 !: 给线圈单元施加电压载荷bfe,all,vltg,,12 ! 12 volts load (real) esel,allharfrq,60 ! 60 Hz.solvefinish!: 选择求解结束，进入后处理器/post1set,1 !: 选择实部plf2d ! Plot real flux lines!: 获取线圈电流n1=node(s,0,0) ! get anode on the coil*get,ireal,node,n1,currset,1,1,,1 读取虚部结果plf2d ! Plot imaginary flux lines *get,imag,node,n1,curr*status ! Imaginary and real components of currentFinishU=12 VI=1.1921976-j1.62066033 AZ=R+jX=U/I=3.53429+j4.80448 ΩR=3.53429 Ω, X= j4.80448 ΩL=0.0127443 H 基于直流计算的电阻：500 3.534292lR S ρ==Ω 视在功率：*7.153199.723962UI P j P jQ ==+=+ （有功功率和无功功率） 基于ANSYS 单元计算的功率损耗：3.57659327.153186W ⨯=（等于P ） 基于ANSYS 单元计算的存贮能量：0.005624320.0112486J ⨯=如果套用212e E I L =的公式，I e 为有效值，得到 22110.012897W 24e E I L I L === 比基于单元计算的结果略大10%。
原因不清楚。
将线圈高度延长4倍，保持线圈匝数、厚度、填充因子等其他参数不变，问电感、电阻、损耗等将如何变化？U=12 VI=1.629333-j4.41289096 AZ=R+jX=U/I=0.883589+j2.39306 ΩR=0.883589 Ω, X= j2.39306 ΩL=0.00634778 H 基于直流计算的电阻：5000.883573lR S ρ==Ω 能量损耗：*9.776326.47732UI P j P jQ ==+=+ 基于ANSYS 单元计算的存贮能量：0.01423120.028462J ⨯=换算成平均功率（j ω⨯）：基于ANSYS 单元计算的工具损耗：4.88829.776W ⨯=作业：2D非线性时谐分析（教材实例2）一厚度5mm的无限大钢板位于磁场强度Hm=2644.1A/m的线圈中，频率50Hz。