电机设计方法.docVIP

第2章 电磁场有限元分析简介 电磁场的边值问题实际上是求解给定边界条件下的麦克斯韦（Maxwell）方程组及由方程组深化出的其他偏微分方程问题。从求解问题的技术手段上来说，它可以分为解析求解和数值求解两大类。对于简单模型，有时可以得到方程的解析解。若模型复杂度增加，则往往很难获得模型的解析解。随着计算工具，特别是高速大容量电子计算机的发展，电磁场数值分析已深入到工业生产各个领域，解决问题的面越来越广，分析的问题也日趋复杂。电磁场数值分析是一门综合性的学科，涉及电磁场理论、数值分析、计算方法、计算机基础知识及高级语言等多个方面，但在计算上存在着共性。有限元法是一种常用的数值方法，并有相应的电磁软件问世，其中ANSOFT公司的Maxwell 3D/2D就是非常优秀的电磁分析软件。 本章将对电磁场的基本理论、电磁场有限元的求解及ANSOFT公司的Maxwell 3D/2D作简单的介绍。至于完整的电磁理论描述，读者可以参考许多教科书。如果读者已熟悉电磁理论，完全可以略过本章，直接从第2章开始学习如何使用Maxwell电磁软件。 1.1电磁场基本理论 1.1.1麦克斯韦方程组 在19世纪中叶，麦克斯韦在总结前人工作的基础上，提出了适用于所有宏观电磁现象的数学模型，称之为麦克期韦方程组。它是电磁场理论的基础，也是工程电磁场数值分析的出发点。 麦克斯韦方程组包括微分和积分两种形式，在此仅给出它们的微分形式，通过它们可以导出能用有限元处理电磁问题的微分方程。 麦克斯韦方程组为 法拉第电磁感应定律 麦克斯韦-安培定律 高斯电通定律 高斯磁通定律 电荷守恒定律 式中，E为电场强度，V/m；D为电通量密度，C/m；H为磁场强度，A/m；B为磁通量密度，T；J为电流密度，A/㎡；P为电荷密度C/m3。 上面5个方程中包含两个旋度方程式（1.1）、式（1.2）和3个散度方程式（1.3）、式（1.4）和式（1.5）。  1.1.2麦克斯韦方程组各方程之间的关系 上面提到的麦克斯韦方程组的5个方程中，只有3个方程是独立的，另外两个相关方程可以从独立方程中导出。其中两个旋度方程肯定是独立方程，另外一个独立方程可以在散度方程式（1.3）和式（1.5）中任选一个，方程式（1.4）只能作为相关方程。读者可以参考表1.1。 表1.1 麦克斯韦方程组中的独立方程与相关方程 独 立 方 程 相 关 方 程 1 2 3 (1.1) (1.2) (1.3) (1.1) (1.2) (1.5) 1 2 (1.4) (1.5) (1.4) (1.3) 1、方程式（1.1）与式（1.4）的关系 对方程式（1.1）两边取散度，有 根据矢量恒等式，可知式（1.6）左端恒等于零 设在场域内B关于时间和场点二阶混合偏导数连续，则式（1.7）可以化为 即 C是与时间无关的常数。同理，也是与时间无关的常数，只要在初始时刻t=0取C=0，则在t＞0以后的任意时刻恒有 由此，可以看出方程式（1.1）与式（1.4）是相关的，由方程式（1.1）可以推导出式（1.4）。 2、方程式（1.2）、式（1.3）与式（1.5）之间的关系 对方程式（1.2）两边取散度，有 显然，如果仅仅利用方程式（1.2）不能同时导出方程式（1.3）和式（1.5）。这时，要私将方程式（1.3）设为独立方程，联合方程式（1.11）推导出方程式（1.5）；要么将方程式（1.5）设定为独立方程，联合方程式（1.11）推导出方程式（1.3）。 1）方程式（1.11）与式（1.3）联合推导式（1.5） 将方程式（1.3）代入式（1.11）有 2）方程式（1.11）与式（1.5）联合推导式（1.3） 将方程式（1.5）代入式（1.11）有 即 这里，C为与时间无关的常数，那么只要在初始时刻t=0取C=0，则在t＞0以后的任意时刻恒有 1.1.3 本构关系 场量E、D、B、H之间的关系，由媒质的特性决定，对于线性介质，本构关系为 D= B= J= 式中，ε为介质的介电常数，F/m；μ为介质的磁导率，H/m；σ为介质的电导率，S/m。 还需要说明的是，对于各向同性介质，ε、μ和σ是标量；对于各向导性介质，它们是张量。 如果希望得到电磁场问题的惟一解，除了上述方程外，还需要配备定解条件；对于瞬变场，需要配备边界条件和初始条件；对于静态场、稳态场、时谐场，只需配备边界条件。 1.1.4 二阶电磁场微分方程 在实际有限元计算中，通常并不针对麦克斯韦方程组中的一阶方程，常常先将方程化为二阶方程，然后针对二阶方程进行有限元数值求解。实际上，比较方便的做法是根据场的基本性质，引入辅助的计算量，如标量电势?、矢量磁位A