电机设计方法

《电机设计方法》由会员分享，可在线阅读，更多相关《电机设计方法（18页珍藏版）》请在金锄头文库上搜索。

1、第2章电磁场有限元分析简介电磁场的边值问题实际上是求解给定边界条件下的麦克斯韦(Maxwell )方程组及由方程组深化出的其他偏微分方程问题。从求解问题的技术手段上来说，它可以分为解析求解和数值求解两大类。对于简单模型，有时可以得到方程的解析解。若模型复杂度增加，则往往 很难获得模型的解析解。 随着计算工具，特别是高速大容量电子计算机的发展，电磁场数值分析已深入到工业生产各个领域，解决问题的面越来越广， 分析的问题也日趋复杂。电磁场数值分析是一门综合性的学科，涉及电磁场理论、 数值分析、计算方法、计算机基础知识及高级语言等多个方面，但在计算上存在着共性。有限元法是一种常用的数值方法，并有相应的

2、电磁软件问世，其中 ANSOFT公司的Maxwell 3D/2D就是非常优秀的电磁分析软件。本章将对电磁场的基本理论、电磁场有限元的求解及ANSOFT公司的Maxwell 3D/2D作简单的介绍。至于完整的电磁理论描述，读者可以参考许多教科书。如果读者已熟悉电磁 理论，完全可以略过本章，直接从第2章开始学习如何使用 Maxwell电磁软件。1.1电磁场基本理论1.1.1麦克斯韦方程组在19世纪中叶，麦克斯韦在总结前人工作的基础上，提出了适用于所有宏观电磁现象 的数学模型，称之为麦克期韦方程组。 它是电磁场理论的基础， 也是工程电磁场数值分析的 出发点。麦克斯韦方程组包括微分和积分两种形式，在此

3、仅给出它们的微分形式，通过它们可以导出能用有限元处理电磁问题的微分方程。麦克斯韦方程组为法拉第电磁感应定律3B?xB=at麦克斯韦-安培定律dD高斯电通定律高斯磁通定律F - B = L电荷守恒定律式中，E为电场强度，V/m ; D为电通量密度，C/m ; H为磁场强度，A/m ; B为磁通量 密度，T; J为电流密度，A/m2； P为电荷密度 C/m3。上面 5 个方程中包含两个旋度方程式（ 1.1 ）、式（ 1.2 ）和 3 个散度方程式（ 1.3 ）、式 （1.4 ）和式（ 1.5 ）。1.1.2麦克斯韦方程组各方程之间的关系上面提到的麦克斯韦方程组的5个方程中，只有3个方程是独立的，另

4、外两个相关方程可以从独立方程中导出。其中两个旋度方程肯定是独立方程，另外一个独立方程可以在散度方程式（1.3 ）和式（1.5 ）中任选一个，方程式（1.4 ）只能作为相关方程。读者可以参考 表 1.1。表1.1麦克斯韦方程组中的独立方程与相关方程独立方程相关方程1231 2(1.1)(1.2)(1.3)(1.4)(1.5)(1.1)(1.2)(1.5)(1.4)(1.3)1、方程式（1.1）与式（1.4）的关系对方程式（1.1）两边取散度，有0BV Vx E） = -V at根据矢量恒等式，可知式（1.6）左端恒等于零SB设在场域内B关于时间和场点二阶混合偏导数连续，则式（1.7）可以化为d-

5、V-B = OSt即V B = CC是与时间无关的常数。同理，普朋也是与时间无关的常数，只要在初始时刻t=0取C=0,则在t 0以后的任意时刻恒有V-B = O由此，可以看出方程式（1.1）与式（1.4）是相关的，由方程式（1.1 ）可以推导出式（1.4）。2、方程式（1.2）、式（1.3）与式（1.5）之间的关系对方程式（1.2）两边取散度，有V * （V x H） = V-J+V-D=fl显然，如果仅仅利用方程式（1.2）不能同时导出方程式（1.3）和式（1.5） o这时，要 私将方程式（1.3）设为独立方程，联合方程式（1.11）推导出方程式（1.5）;要么将方程式（1.5 ）设定为独立

6、方程，联合方程式（1.11）推导出方程式（1.3 ）o1）方程式（1.11）与式（1.3）联合推导式（1.5）将方程式（1.3）代入式（1.11）有2）方程式（1.11）与式（1.5）联合推导式（1.3）将方程式（1.5）代入式（1.11）有訶 D-p) = QV-D-p = C这里，C为与时间无关的常数， 那么只要在初始时刻 t=0取C=0,则在t0以后的任意时刻恒有1.1.3本构关系场量E、D、B、H之间的关系，由媒质的特性决定，对于线性介质，本构关系为D=.lB= i：J=j 1.式中，&为介质的介电常数，F/m ;卩为介质的磁导率，H/m ;为介质的电导率，S/m。还需要说明的是，对于

7、各向同性介质，&、卩和b是标量；对于各向导性介质，它们是张量。如果希望得到电磁场问题的惟一解，除了上述方程外，还需要配备定解条件；对于瞬变 场，需要配备边界条件和初始条件；对于静态场、稳态场、时谐场，只需配备边界条件。1.1.4二阶电磁场微分方程在实际有限元计算中，通常并不针对麦克斯韦方程组中的一阶方程，常常先将方程化为二阶方程，然后针对二阶方程进行有限元数值求解。实际上，比较方便的做法是根据场的基本性质，引入辅助的计算量，如标量电势？、矢量磁位A等。Maxwell常用的求解方程有二维、三维静电场求解器所满足的泊松方程V (V$) = -p二维稳恒电场求解器所满足的拉普拉斯方程V (crV(|

8、) = 0二维交变电场求解器所满足的复数拉普拉斯方程V- (o+jwe)V|) = 0二维静磁场求解器所满足的非齐次标量波动方程二维涡流场求解器所满足的波动方程组V7 X x A) = (V(|) jtj)A)(u+ jci)It 二I dfl = J (Vcj) ja)A)(cj + k口3二维轴向磁场涡流求解器所满足的齐次波动方程VxH +jwpH = O三维静磁场和涡流求解器所满足的齐次波动方程组（去 *xH）+ j 硼1. 1. 5电磁场求解的边界条件电磁场问题求解中，有各种各样的边界条件，结合Maxwell 3D / 2D，归结起来可概 括为 6 类。1、自然边界条件自然边界条件是软

9、件系统的默认边界条件， 不需要用户指定， 是不同媒质交界面场量 的切向和法向边界条件。2诺伊曼边界条件电磁场教科书中常常称诺伊曼边界条件为第二类边界条件， 它规定了边界处势的法向 导数分布。 Maxwell 所提到的是齐次诺伊曼边界，即法向导数为零。它是 Maxwell 系统默认 边界条件，不需要用户指定。3狄利克莱边界条件 电磁场教科书中常常称狄利克莱边界条件为第一类边界条件， 有限元计算领域， 常常 称其为约束边界条件，或本质边界条件。它规定了边界处势的分布，势是边界位置的函数， 也可以是常数和零。4对称边界条件 对称边界条件包括奇对称和偶对称两大类。奇对称边界可以模拟一个设备的对称面，

10、在对称面的两侧电荷、 电位、 电流等满足大小相等， 符号相反。 偶对称边界可以模拟一个设 备的对称面，在对称面的两侧电荷、电位、电流等满足大小相等，符号相同。采用对称边界 条件可以减小模型的尺寸，有效地节省计算资源。5匹配边界条件 匹配边界条件是模拟周期性结构的对称面，使主边界和从边界场量具有相同的幅度 （ 对 于时谐量还有相位 ），相同或相反的方向。6气球边界条件气球边界条件是 Maxwell 2D 求解器常见的边界条件，常常指定在求解区的外边界处 , 用于模拟绝缘系统等。除此之外，有一些求解器中还有各自“特色”的边界条件如交变电场中的电阻边界、 涡流场中的阻抗边界，主要用来模拟很薄的介质层

11、。12 电磁场求解的有限元方法所谓的有限元法， 就是将整个区域分割成许多很小的子区域， 这些子区域通常称为 “单 元”或“有限元”，将求解边界问题的原理应用于这些子区域中，求解每个小区域，通过选取 恰当的尝试函数， 使得对每一个单元的计算变得非常简单， 经过对每个单元进行重复而简单 的计算， 再将其结果总和起来， 便可以得到用整体矩阵表达的整个区域的解， 这一整体矩阵又常常是稀疏矩阵，可以更进一步简化和加快求解过程。由于计算机非常适合重复性的计算F面就以一个简单的和处理过程，因此整体矩阵的形成过程很容易使用计算机处理来实现。例子说明有限元法的基本原理。1. 2. 1 一维有限元法1.【例11】

12、问题的描述考虑一个两极板电容器的静电场分布问题。极板间充有密度为-的自由电荷，即自由电荷密度恰好等于介质的介电常数。前面已经介绍静电场所满足方程式(1.19)。考虑到极板间只有一种介质，可以导出本例中静电场满足的方程为= 1假定极板都接在 0.5V电源端，极板的间距为 2,如图1.1所示。由于电容器的激励和几何形状都关于Y轴对称，只要求解整个区域的一半即可，而另外一半可由对称关系得出。从边界条件上看，这种对称结构导致电力线垂直穿过Y轴，使电势在该对称轴上沿 x方向的变化率为零，即对称面可以用齐次诺伊曼条件表示， 为简化起 见，这里没有考虑实际的物理单位。受到狄里克莱和齐次诺伊曼边界条件的约束，

13、即。：；二:=0x=0假定电容器极板的尺寸远远大于极板间的距离t那么电容器的电势分布问题简化为一维边值问题。2.有限元求解用有限元求解问题的第一步就是划分单元。一般说来.单元数越多，则近似解的精度就越高，当然计算量也就越大，越费时。所以单元数应该足够多，以保证精度。对于例【1.11所考虑的问题，整个区域(0, 1)被分割成4个单元，记为单元 el、单元e2、单元e3和单元 e4。这些单元大小可以相同，也可以不同。在实际问题中，根据场分布的疏密程度，有限元 可以具有不同的尺寸，这样便于处理复杂的几何结构和激励源。分割后的区域由4个单元和5个点组成，这些点称为节点”，对应于这5个点的电 势记为？1

14、、?2、?3、？4、?5，每一个单元都由相邻的两个节点所限定，如图1 . 2所示。对于一维空间来说， 一个单元只是一个线段。 对二维空间来说， 有限元可以有各种形 状，如三角形、矩形等。作为一种数值计算方法，有限元并非用来寻求问题的解析解。实际上许多工程问题目前都无法找到解析解。有限元的作用就在于求解分布场的势函数在每个节点上的近似值，而势函数在单元其他位置的值，可以用插值方法获得。 如果采用线性插值方法表示势函数，则称为一阶有限元。如果采用高阶插值表示势函数，则称为高阶有限元。本节只介绍一阶有限丿元。在图1.3中，将任意单元记为e”，对应于这一有限元有两个节点：却和险1，这两个节点上的电势分

15、别记为 :和群他显然它们为待定的未知量，称为自由度。对于一阶有限元，由于采用线性插值，如果将单元上的电势分布用图形表示，实际上就成为一条连接两个节点电势值的线段。这一分布函数记为？e，一旦节点上的电势被求出，在单元上的其他各点的电势值即可由线性关系得到。显然，求解电势分布的关键是找到节点上的电势值。采用加权余数法（伽辽金法）或变分法（里兹法）可以得到以下代数方程组式中，匚为尝试函数；.为各待定系数。可以令待定系数 Ci为各节点上的电势值 ?i,这样一旦解出各待定系数，也就获得了节点上的电势值，也就是说，求解待定系数和求解节点电势成为一个统一的计算过程t这是有限元法的巧妙处之一。另外，还需要设法使方程的近似解满足边界处的狄里克莱条件。实际上，采用有限元法，满足这一边界条件并不困难， 只要令狄里克莱边界上的节点电势为给定的值即可，同