电磁场的数值计算方法

1、 电磁场的数值计算方法摘要：数值计算方法是一种研究并解决数学问题数值近似解的方法，广泛运用于电气、军事、经济、生态、医疗、天文、地质等众多领域。本文综述了电磁场数值计算方法的发展历史、分类，详细介绍了三种典型的数值计算方法有限差分法、有限元法、矩量法,对每种方法的解题思路、原理、步骤、特点、应用进行了详细阐述,并就不同方法的区别进行了深入分析,最后对电磁场数值计算方法的应用前景作了初步探讨。关键词：电磁场；数值计算；有限差分法；有限元法；矩量法引言自从1864年Maxwell建立了统一的电磁场理论，并得出著名的Maxwell方程以来，经典的数学分析方法是一百多年来电磁学学科发展中一个极为重要的

2、手段,围绕电磁分布边值问题的求解国内外专家学者做了大量的工作。在数值计算方法之前,电磁分布的边值问题的研究方法主要是解析法，但其推导过程相当繁琐和困难,缺乏通用性,可求解的问题非常有限。上个世纪六十年代以来,伴随着电子计算机技术的飞速发展,多种电磁场数值计算方法不断涌现,并得到广泛地应用,相对于解析法而言,数值计算方法受边界形状的约束大为减少,可以解决各种类型的复杂问题。但各种数值计算方法都有一定的局限性,一个复杂的问题往往难以依靠一种单一方法解决，因此如何充分发挥各种方法的优势,取长补短,将多种方法结合起来解决实际问题,即混合法的研究和应用已日益受到人们的关注。本文综述电磁场的数值计算方法，

3、对三种常用的电磁场数值计算方法进行分类和比较。1电磁场数值计算方法的发展历史在上世纪四十年代，就有人试探用数值计算的方法来求解具有简单边界的电磁场问题，如采用Ritz法1，以多项式在整个求解场域范围内整体逼近二阶偏微分方程在求解域中的解。五十年代，采用差分方程近似二阶偏微分方程，诞生了有限差分数值计算方法，开始是人工计算，后来采用机械式的手摇计算机计算，使简单、直观的有限差分法得到应用和发展，该方法曾在欧、美风行一时。1964年美国加州大学学者Winslow以矢量位为求解变量，用有限差分法在计算机上成功地解算了二维非线性磁场1，此后有限差分法在工程电磁场计算领域大为发展。1965年，Winsl

4、ow首先将有限元法从力学界引入电气工程中，1969年加拿大MeGill大学P.P.Silvester运用有限元法成功地进行了波导的计算；七十年代初,P.P.Silvester和M.V.K.Chari合作将有限元法应用于二维非线性磁场的计算，成功地计算了直流电机、同步电机的恒定磁场。此后有关有限元法探讨的论文越来越多，有限元法运用的范围由静态场到涡流场到辐射场，由线性场到非线性场，由各向同性媒质到各向异性、要考虑磁滞损耗，由工程电磁场到生物电磁场等等。有人认为有限元法是求解工程电磁场的最有效最成功的方法。有限元法和有限差分法都是求解边值问题的方法，属于微分方程法。对于开区域或要求求解连续分布场量

5、的区域，这类方法就会受到自身的限制。1972年英国卢瑟福实验室的C.W.Trowbridge等人提出了积分方程法的思想，给出了二维、三维场问题的离散形式2，由于此种方法只需离散源区，不需考虑边界条件，所以它较好地解决了无界开域场和要求连续计算场量的问题。该方法计算精度高，但计算量很大。该实验室Sinkin等人又在积分方程法基础上提出了边界积分方程法（又称边界元法），用此解决线性场的计算，计算量大为减小。此后该室的学者们将积分方程与微分方程法结合起来，提出了求解三维静磁场的双标量位在解决天线辐射场、散射场问题中，矩量法是一个很重要的数值计算方法。1968年R.F.Harrington发表了专著“

6、FieldcomputationbyMomentMethod”，对散射场、天线辐射场、波导场等方面的问题起了很好的推进作用。除以上所介绍的方法外，随着电磁场数值分析的不断发展，各种新方法不断涌现，如计算电场的模拟电荷法，最小二乘配点法，求解磁场的模拟电流法，以及计算场的图论模型法，快速Fourier变换法、有限体元法、无网格计算法等等。各种方法互相配合，出现了一些混合方法，如：矩量法模拟电荷法、模拟电荷法有限元法、有限元法边界元法等，有效地解决了一些实际问题。近年来人工神经网络，小波理论3等也引入了电磁场的数值计算中，瞬态电磁场计算如时域有限差分法的应用有了长足的发展。总之随着现有的电磁场数值

7、计算方法的不断深入发展、提高和完善，新的方法不断产生。在电磁场的数值解法不断发展的同时，人们并没有忘记长期以来所运用的解析方法。解析法计算结果精确，且可以用解析式表达计算结果，受这些特点吸引，解析法与数值计算方法相结合形成的半解析法应运而生，也成为了一种主流解算方法，并还在不断发展。电磁场数值计算方法发展走向成熟的一个重要标志是：成熟的方法越来越多地应用于工程实际问题中，商业化通用软件包不断出现4。一个商业化软件包通常由下面几部分组成：数据定义：几何尺寸、材料性能参数、边界条件前处理模拟化：空调剖分、网格自动产生、节点形成网格图形显示离散方程组系数矩阵形成数据处理求解代数方程组、非线性叠代按要

8、求输出计算结果场图显示（含线性媒质和非线性媒质区）后处理受力和损耗计算与图形显示局部场域分布的精细计算与显示以上三部分中前、后处理占用了软件包语句的90%以上，编程的主要工作量在此，而数据处理，也就是我们目前正在学习的数值计算方法仅占软件语句的10%以内，但它却是占用计算机内存量和消耗CPU时间的主要部分。电磁场数值计算方法的分类求解电磁问题的最终要求就是获得满足实际条件的Maxwell方程的解，借助于计算数学中的数值算法能够得到大多数电磁问题的近似解。数值算法的基本思想5就是把连续变量函数离散化，把微分方程化为差分方程；把积分方程化为有限和的形式，从而建立起收敛的代数方程组，然后利用计算机技

9、术进行求解。数值计算方法从求解方程的形式看，主要分为积分方程法和微分方程法两大类。积分方程法主要有矩量法和边界元法，微分方程法主要有有限差分法和有限元法。对两种方程法的比较，如表一所示。表一积分方程法和微分方程法的比较积分方程法微分方程法共性对场问题的处理是一致的，即需离散化场域，结果是数值解不同占八、离散域仅在场源区，无需对整个场域离散整个场域计算对象Q口场里先求位函数，再求场量求解域可在场域内某一局部区域求解，也可在全场域内求解全场域内求解计算程度较咼较低应用不适用边界区域复杂的场域边界形状复杂的场域较易处理联系两种方法的结合形式，可处理较复杂的电磁场问题几种重要的数值计算方法有限差分法在

10、电磁场数值计算方法中，有限差分法是应用最早的一种方法。有限差分法以其概念清晰，方法简单、直观，有大致固定的处理和计算模式，具有一定的通用性等特点，在电磁场数值分析领域内得到了广泛的应用。有限差分法的基本原理有限差分法的基本思想是把连续的定解区域用有限个离散点构成的网格来代替，这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上

11、的近似解。差分与差商设函数f(x)的自变量X有一小增量Ax二h，则f(x)的增量为纣(x)=f(X+h)-f(x)(3.1)Af(x)为函数f(x)的一阶差分。当增量h足够小，差分Af与微分df之间的差才足够小。一阶差分纣是自变量x的函数。按式(1)，计算Af(x)的差分A2f(x)称二阶差分，且A2f(x)二Af(x+h)-Af(x)(3.2)函数f(x)的一阶导数f(x)为fQ)=f=lim业dxAxAxtO应用差分，f(x)可表示为f,()Af(x)_f(x+h)f(x)f(x)_h故f(x)可表示为差分Af(x)除以有限小差分Ax的商，称为差商。同理，函数f(x)的二阶导数f(x)可表

12、示为Ax乎)x+Axdxxd2fdx23.4)11f(x+h)f(x)f(x)f(xh)hIhhf(x+h)2f(x)+f(xh)h2差分方程的构造现以二维静态电、磁场泊松方程的第一类边值问题为例，来具体阐明有限差分法的应用。设具有平行平面场特征的电磁场场域D,如图1所示,为一由闭合边界L所界定的平面域，其定解条件可表述为3.5)V2u6,y)_乜+乜_F6,y)G,y)eDQx2dy2uy)3.6)对于所给定的偏微分方程定解问题，应用有限差分法，首先需从网格剖分着手决定离散点的分布方式。原则上，可以采用任意的网格剖分方式，但这将直接影响所得差分方程的具体内容，进而影响解题的经济性与计算精度。

13、为简化问题，通常采用完全有规律的分布方式，这样在每个离散点上就能得出相同形式的差分方程，有效地提高解题速度，因而经常采用正方形的网格的剖分方式。现即以这种正方形网格剖分场域D,也就是说，用分别与x,y两坐标轴平行的两簇等距网格线来生成正方形网格，即x二x二ihi(i二0,l,2)y二y.二jhj(j二0,l,2)h为步长，网格线的交点O(,y)称为节ij点，这样D域就离散化为由网格节点标成的离散点得集合。对场域D中节点OC,y)是一典型ij节点，它与周围的1，2,3和4点构成一个对称星型。设这些离散点上待求函数的近似值记为x图1正方形网格划分u=u(i,j)uu(i+1,/)uu(i,j+1)

14、u=u(i1,j)uu(i,j1)01234则式(6)可近似离散化为-1L(i+1,j)一2u(i,j)+u(i1,j)+丄L(i,j+1)一2u(i,j)+u(i,j1)=F(3.7)h2h2u(i+1,j)+u(i,j+1)+u(i1,j)+u(i,j1)4u(i,j)h2F(3.8)若式(6)F=0,则节点O上函数u的值等于其四周相邻点函数值的平均。因为差分方程(7),(8)只出现待求函数u在点OC,y)及其四个临近点的值，故ij称之为五点差分格式6，根据差分方程组解出各离散点处的待求函数值。有限元法传统的变分法在20世纪二三十年代为其新型时期，理论上发展很快，各种变分问题的最后求解都可

15、归结为解尤拉方程的边值问题，然而只有在一些特殊情况下尤拉方程才能求出精确解，在大多数情况下，尤拉方程的精确解无法求出。四五十年代，随着计算机的出现，使其在实际应用中逐渐为比较灵活、通用的有限差分法所替代。但是，有限差分法在理论上没有以变分原理为基础，因而其收敛性和数值稳定性往往得不到保证。随后发展形成的有限元法正是变分法与有限差分法相结合的成果，它取长补短地在理论上以变分原理为基础，在具体方法构造上又利用了有限差分法网格离散化处理的思想7。3.2.1有限元法的基本原理有限元法是以变分原理为基础，将要求解的微分方程型数学模型边值问题，首先转化为相应的变分问题，即泛函求极值问题；然后，利用剖分插值

16、将变分问题离散化为普通多元函数的极值问题，最终归结为一组多元的代数方程组，求解该方程组，从而获得边值问题的数值解。泛函、变分问题简介在微积分学形成初期，以数学物理问题为背景，与多元函数的极值问题相对应，已在几何、力学上提出了若干个求解泛函极值的问题如图2中的质点最速降线问题所述，质点A图2最速降线问题从定点(x,y)自由下滑到定点B(x,y)，试求1122使滑行时间最短的质点下滑轨道y二y(x)。图示滑行弧段ds所需时间为Id_ds_secadx_寸1+y2dxdt_vJ2gyv2gy滑行总时间为Jy(x)_Ty(x)_JTdt_Jx2工2dx(3.9)0 x1x12gy(3.9)式J_Jy(

17、x)不仅取决于积分端点x活口x2，而且取决于y_y(x)的选取。J取决于y(x)，所以J是函数y(x)的函数，称之为y(x)的泛函，记作Jy(x)。于是所述之最速降线问题，在数学上就归结为研究泛函Jy(x)的极值问题，即Jy(x)_J七芈匸dx_min、x1(3-10)、y(xi)_0y(x2)_y2泛函的极值问题就称为变分问题。对一般问题而言，可导出下列对应于一个自变量x、单个函数y(x)及其导数y(x)的已知函数Jy_Jx2F(x,y,y)dx(3.11)x式中F为x、y和yr的已知函数。泛函Jy的自变量不是一般的自变量，而是一个或几个函数所属的函数族y(x)。在端点x和x上分别等于给定值

18、的无数12个函数y(x)中，仅有一个亍(x)能使定积分Jy达到极小值，此函亍(x)数称为极值函数。因此，变分问题就在于寻求使泛函达到极值的该极值函数y(x)，即分析研究泛函的极值问题。泛函的变分与尤拉方程泛函变分问题的经典解法有两种，一种称之为直接解法，另一类是间接解法。直接解法是直接把泛函的极值问题近似地转化为一般多元函数的极值问题，用有限维子空间中的函数去逼近无穷维空间中的极值函数，从而近似求得泛函的极值。间接解法是将变分问题转化为尤拉方程(微分方程)的定解问题，即边值问题来求解。以式(3.11)这种最简形式来推导尤拉方程。设函数y(x)稍有变化，记作y+5y，8y称之为y(x)的变分，它

19、反映了整个函数的变化量。这样泛函Jy的值也应随之变动，相应于变分5y的泛函增量为AJ二Jy+5y-Jy二I冬f(x,y+6y,y+y)F(x,y,y)dx(3.12)x1将(3.12)式由多元函数的泰勒公式展开FOF1B2FB2FB2FAJJF5y+F+2磐(8y)2+2等+笋y)2+血ByBy2By2ByByBy2=J+52J+53J+(3.13)式中作为泛函增量AJ的线性主部为BFBF5J=Ix2一5y+5yrdx(3.14)x1ByBy5J称为泛函Jy的一次变分(简称变分)。而52j、53j分别是函数变分5y及其导数5y，的二次、三次齐次式等的积分，依次称为二次变分，三次变分令变分问题的

20、解为y=y(x)，且设极值解y=y(x)稍有变动y+5y，令5y二眄(x)(3.15)式中s为任意给定的微量实参数，s值就确定了y二y(x,e)函数族中的某一曲线，进而确定泛函Jy(x,s)之值；而耳(x)是定义于区间x,x且满足12)=0齐次边界条件的可微函数。于是泛函J(y+sn)=Jy(x,s)=(s)12就成为变量s的函数，且当s二0时获极值函数的解。(s)在s二0时取得极值的必要条件是dF（）L=（0）=Jx2：5y+：5ydx二0內Qy利用分部积分，并根据变分与微分顺序可互换原理，（3.16）式可写为3.16）=0 x13.17）3Fd3F1X2石dx（而）5曲=0由于（3.17）

21、对任意5y均成立，故有3.18）方程（3.18）就称为泛函（3.11）的极值问题的尤拉方程。综上所述，有限元法的基本特点是：（1）离散化过程保持了明显的物理意义。这是因为，变分原理描述了支配物理现象的物理学中的最小作用原理（如力学中的最小势能原理、静电学中的汤姆逊定理等）。（2）优异的解题能力。与其他数值计算方法相比较，有限元法在适应场域边界几何形状，以及媒质物理性质变异情况复杂的问题求解上，有突出的优点。（3）从数学理论意义上讲，有限元法作为应用数学的一个重要分支，很少有其他方法应用的这样广泛。它使微分方程的解法和理论面目一新，推动了泛函分析与计算方法的发展。矩量法矩量法，是近年来在天线、微

22、波技术和电磁波散射等方面广泛应用的一种方法。从这些实际问题涉及开域、激励场源分布形态较为复杂等特征出发，矩量法是将待求的积分方程问题转化为一个矩阵方程的问题7，借助于计算机，求得其数值解，从而在所得激励源分布的数值解基础上，即可算出辐射场的分布及其波阻抗等特性参数。3.3.1矩量法的基本原理先选定基函数对未知函数进行近似展开，代入算子方程，再选取适当的权函数，使在加权平均的意义下方程的余量等于0，由此将连续的算子方程转换为代数方程。原则上，矩量法可用于求解微分方程和积分方程，但用于微分方程时所得到的代数方程组的系数矩阵往往是病态的，故在电磁场问题中主要用于求解积分方程。加权余量法设给定边值问题

23、的场方程统一表述为如下的算子方程，即L(f)二g(3.19)已知边界条件为u=u(r)(3.20)ssb鋼=q(r)(3.21)dns2sb其中：L是线性算子，f是待求函数，g是已知的源。若u为精确解，贝V方程(3.19)和边界条件(3.20),(3.21)应该完全满足。但大多数情况下，不能得到u的精确解，只能通过数值方法进行估计。构造一个由有限个线性无关函数N(i=1,2,，n)所组成的基函数集合in，借以展开待求函数u的近似解为=YNu=Inb(3.22)iii=1将代入式(3.19)中必然存在误差，即R()=L-g(3.23)取一个归属于试探函数的权函数集合b，令iW(Lg)dV=0(j

24、=1,2,，n)(3.24)jv式(3.24)由n个方程构成的方程组，它等价于人为地强制近似解，使其因不能精确地满足场方程而导致的误差在平均的含义上等于零。按式(3.24)展开，所构成的各种求解积分或微分方程近似解的方法可被统称为加权余量法8。因为按给定权函数W展开式的式(3.24)，即意味着余量R=Lu-g对W取矩的一组平衡式，故式(3.24)的构造亦就被称为矩量法。基于加权余量式(3.24)，进行移项处理，便得(j=1,2,n)3.25)iWLudV=iWgdVvjvj将式(3.22)代入式(3.25)的左端，有vjNu11丿i=1=uiWL(N)dVijii=1v3.26)为了书写方便，

25、令工uiWL(N)dV三:W,L(N)：和iWgdV=fW,gi=1ivjijivjj代入(3.26)式，则可写成3.27)工uJWL(N)dV=:W,g；：(j二1,2,n)ivjiji=1这样，即展开成含n个未知数u的n个方程。若用矩阵形式表示，则有ilKJ=g(3.28)综上所述，矩量法的特点是：矩量法将连续方程离散化为代数方程组,既适用于求解微分方程,又适用于求解积分方程。它的求解过程简单,求解步骤统一应用起来比较方便，然而需要一定的数学技巧,如离散化的程度、基函数与权函数的选取,矩阵求解过程等。另外必须指出的是,矩量法可以达到所需要的精确度、解析部分简单,可计算量很大,即使用高速大容

26、量计算机,计算任务也很繁琐。电磁场数值计算的应用前景电磁场数值计算方法近六十年来发展如此之快，除由于从事这方面的科研人员的努力之外，主要是其研究成果迅速被电机、电器、变压器、加速器、微波器件、计算机磁头等领域采用9，对改善产品性能、降低生产成本，起着越来越大的作用。众所周知，任何电磁器件，包括国民经济中应用极其广泛的电机与变压器，其能量转换都是通过电磁场来实现的。但传统的电磁产品的设计方式由于客观条件与手段的限制，把场的实际分布参数当作集中参数处理，不可避免地带来相当大的误差。在迫不得已的情况下，只能用模拟、实验等方法处理，其耗费大、周期长，可借鉴的经验不多，而实验的模拟也不都是有条件的，此外

27、，工农业和日常生活中所用电磁装置越来越多，生产和销售竞争剧烈，因此有效的设计方法显然受到重视。在采用电磁场数值计算以前的任何方法，即使是十分精巧的代数解析方法，也只能适用于特别简单的几何结构以及一些特殊的假定模型，有时模型甚至简化到不能容忍的程度，但除此之外别无它法。因为实际的电磁场问题，特别是电机电磁场，其边界情况十分复杂，加上铁的饱和以及导体中的涡流效应，解析方法是无能为力的。而数值方法可以模拟复杂的形状，可以适应非线性问题，可以进行涡流的分布计算。电磁场的数值计算与流体力学分析、温度分析、机械应力分析、电磁力分析和生产计划等等有机联系起来，由计算机完成全部设计，构成所谓CAE1(系统(C

28、omputerAidedEngineering)。CAE系统基本可包括设计与制造的全部过程，产品的设计、制造方案、准备零件、草图、计算、成本核算、生产、机床数控和试验、模拟和产品的自动测试都可能列入到CAE系统中。目前CAE系统在西方国家发展很快，估计每年增长40%，尤其在机械行业中，形成了所谓CADMAT即计算机辅助设计、加工与试验。从长远观点看，电磁场数值计算要发挥更实际的作用，必须与其他多种学科相结合以形成CAE系统，CAE的发展必然给工业结构上带来巨大的变化。结束语数值计算是一门计算的艺术，电磁场的数值计算横跨了多个学科，是数学理论、电磁理论和计算机应用能力的完美组合。通过本次论文设计

29、，我对差分、变分、泛函、加权余量法等数学知识有了深刻的认识，掌握了有限差分法、有限元法、矩量法的基本原理，并能进行简单的运算。但若用数值计算方法解决复杂的电磁场问题，则在很大程度上依赖于数学知识及计算机编程能力，这就需要进一步系统学习相关的理论知识。参考文献文舸一.计算电磁学的进展与展望J.电子学报，1995,23(10)：62-69.刘圣民.电磁场的数值方法M.武汉：华中理工大学出版社，1991：23-45.邓东皋，彭立中.小波分析J.数学进展，1991,20(3)：294-310.洪伟.计算电磁学研究进展J.东南大学学报(自然科学版)，2002,32(3)：335-339.倪光正，杨仕友，

30、钱秀英等.工程电磁场数值计算M.北京：机械工业出版社，2003：123-161.盛剑霓.工程电磁场数值分析M.西安：西安交通大学出版社，1991：69-75.方静，汪文秉.有限元法和矩量法结合分析背腔天线的辐射特性J.微波学报，2000,16(2)：139-143连汉雄.电磁场理论的数学方法M.北京：北京理工大学出版社，1990:10-15.张平文，刘法启，张宇.小波函数值的计算J.计算数学，1995(3)173-185.楼仁海，符果行，袁敬闳.电磁理论M.成都：电子科技大学出版社，1996：73-102.ElectromagneticnumericalmethodAbstract：Numericalcalculation,whichiswidelyusedinmanyfields,suchaselectric,militaryaffairs,economy,ecology,medicaltreatment,astronomy,geologyandsoonisamethodofgettingthenumericalapproximatesolutioninthestudyingand