1、電磁場(chǎng)的數(shù)值計(jì)算方法物理系0702班 學(xué) 生 杜星星 指導(dǎo)老師 任麗英摘 要：數(shù)值計(jì)算方法是一種研究并解決數(shù)學(xué)問題數(shù)值近似解的方法， 廣泛運(yùn)用于電氣、軍事、經(jīng)濟(jì)、生態(tài)、醫(yī)療、天文、地質(zhì)等眾多領(lǐng)域。 本文綜述了電磁場(chǎng)數(shù)值計(jì)算方法的發(fā)展歷史、分類，詳細(xì)介紹了三種典型的數(shù)值計(jì)算方法有限差分法、有限元法、矩量法, 對(duì)每種方法的解題思路、原理、步驟、特點(diǎn)、應(yīng)用進(jìn)行了詳細(xì)闡述, 并就不同方法的區(qū)別進(jìn)行了深入分析, 最后對(duì)電磁場(chǎng)數(shù)值計(jì)算方法的應(yīng)用前景作了初步探討。關(guān)鍵詞：電磁場(chǎng)；數(shù)值計(jì)算；有限差分法；有限元法；矩量法引言自從1864年Maxwell 建立了統(tǒng)一的電磁場(chǎng)理論，并得出著名的Maxwell 方程以

2、來，經(jīng)典的數(shù)學(xué)分析方法是一百多年來電磁學(xué)學(xué)科發(fā)展中一個(gè)極為重要的手段, 圍繞電磁分布邊值問題的求解國(guó)內(nèi)外專家學(xué)者做了大量的工作。在數(shù)值計(jì)算方法之前, 電磁分布的邊值問題的研究方法主要是解析法，但其推導(dǎo)過程相當(dāng)繁瑣和困難,缺乏通用性,可求解的問題非常有限。上個(gè)世紀(jì)六十年代以來,伴隨著電子計(jì)算機(jī)技術(shù)的飛速發(fā)展,多種電磁場(chǎng)數(shù)值計(jì)算方法不斷涌現(xiàn),并得到廣泛地應(yīng)用,相對(duì)于解析法而言,數(shù)值計(jì)算方法受邊界形狀的約束大為減少,可以解決各種類型的復(fù)雜問題。但各種數(shù)值計(jì)算方法都有一定的局限性,一個(gè)復(fù)雜的問題往往難以依靠一種單一方法解決，因此如何充分發(fā)揮各種方法的優(yōu)勢(shì),取長(zhǎng)補(bǔ)短,將多種方法結(jié)合起來解決實(shí)際問題,即混

3、合法的研究和應(yīng)用已日益受到人們的關(guān)注。本文綜述電磁場(chǎng)的數(shù)值計(jì)算方法，對(duì)三種常用的電磁場(chǎng)數(shù)值計(jì)算方法進(jìn)行分類和比較。1電磁場(chǎng)數(shù)值計(jì)算方法的發(fā)展歷史在上世紀(jì)四十年代，就有人試探用數(shù)值計(jì)算的方法來求解具有簡(jiǎn)單邊界的電磁場(chǎng)問題，如采用Ritz法1，以多項(xiàng)式在整個(gè)求解場(chǎng)域范圍內(nèi)整體逼近二階偏微分方程在求解域中的解。五十年代，采用差分方程近似二階偏微分方程，誕生了有限差分?jǐn)?shù)值計(jì)算方法，開始是人工計(jì)算，后來采用機(jī)械式的手搖計(jì)算機(jī)計(jì)算，使簡(jiǎn)單、直觀的有限差分法得到應(yīng)用和發(fā)展，該方法曾在歐、美風(fēng)行一時(shí)。1964年美國(guó)加州大學(xué)學(xué)者Winslow以矢量位為求解變量，用有限差分法在計(jì)算機(jī)上成功地解算了二維非線性磁場(chǎng)1

4、，此后有限差分法在工程電磁場(chǎng)計(jì)算領(lǐng)域大為發(fā)展。 1965年，Winslow首先將有限元法從力學(xué)界引入電氣工程中，1969年加拿大MeGill大學(xué)P. P. Silvester運(yùn)用有限元法成功地進(jìn)行了波導(dǎo)的計(jì)算2；七十年代初，P. P. Silvester和M. V. K. Chari合作將有限元法應(yīng)用于二維非線性磁場(chǎng)的計(jì)算，成功地計(jì)算了直流電機(jī)、同步電機(jī)的恒定磁場(chǎng)。此后有關(guān)有限元法探討的論文越來越多，有限元法運(yùn)用的范圍由靜態(tài)場(chǎng)到渦流場(chǎng)到輻射場(chǎng)，由線性場(chǎng)到非線性場(chǎng)，由各向同性媒質(zhì)到各向異性、要考慮磁滯損耗，由工程電磁場(chǎng)到生物電磁場(chǎng)等等。有人認(rèn)為有限元法是求解工程電磁場(chǎng)的最有效最成功的方法。 有限

5、元法和有限差分法都是求解邊值問題的方法，屬于微分方程法。對(duì)于開區(qū)域或要求求解連續(xù)分布場(chǎng)量的區(qū)域，這類方法就會(huì)受到自身的限制。1972年英國(guó)盧瑟福實(shí)驗(yàn)室的等人提出了積分方程法的思想，給出了二維、三維場(chǎng)問題的離散形式2，由于此種方法只需離散源區(qū)，不需考慮邊界條件，所以它較好地解決了無界開域場(chǎng)和要求連續(xù)計(jì)算場(chǎng)量的問題。該方法計(jì)算精度高，但計(jì)算量很大。該實(shí)驗(yàn)室Sinkin等人又在積分方程法基礎(chǔ)上提出了邊界積分方程法(又稱邊界元法)，用此解決線性場(chǎng)的計(jì)算，計(jì)算量大為減小。此后該室的學(xué)者們將積分方程與微分方程法結(jié)合起來，提出了求解三維靜磁場(chǎng)的雙標(biāo)量位法等。 在解決天線輻射場(chǎng)、散射場(chǎng)問題中，矩量法是一個(gè)很重

6、要的數(shù)值計(jì)算方法。1968年R. F. Harrington發(fā)表了專著“Field computation by Moment Method”，對(duì)散射場(chǎng)、天線輻射場(chǎng)、波導(dǎo)場(chǎng)等方面的問題起了很好的推進(jìn)作用。除以上所介紹的方法外，隨著電磁場(chǎng)數(shù)值分析的不斷發(fā)展，各種新方法不斷涌現(xiàn)，如計(jì)算電場(chǎng)的模擬電荷法，最小二乘配點(diǎn)法，求解磁場(chǎng)的模擬電流法，以及計(jì)算場(chǎng)的圖論模型法，快速Fourier變換法、有限體元法、無網(wǎng)格計(jì)算法等等。各種方法互相配合，出現(xiàn)了一些混合方法，如：矩量法模擬電荷法、模擬電荷法有限元法、有限元法邊界元法等，有效地解決了一些實(shí)際問題。近年來人工神經(jīng)網(wǎng)絡(luò)，小波理論3等也引入了電磁場(chǎng)的數(shù)值計(jì)算

7、中，瞬態(tài)電磁場(chǎng)計(jì)算如時(shí)域有限差分法的應(yīng)用有了長(zhǎng)足的發(fā)展。總之隨著現(xiàn)有的電磁場(chǎng)數(shù)值計(jì)算方法的不斷深入發(fā)展、提高和完善，新的方法不斷產(chǎn)生。在電磁場(chǎng)的數(shù)值解法不斷發(fā)展的同時(shí)，人們并沒有忘記長(zhǎng)期以來所運(yùn)用的解析方法。解析法計(jì)算結(jié)果精確，且可以用解析式表達(dá)計(jì)算結(jié)果，受這些特點(diǎn)吸引，解析法與數(shù)值計(jì)算方法相結(jié)合形成的半解析法應(yīng)運(yùn)而生，也成為了一種主流解算方法，并還在不斷發(fā)展。 電磁場(chǎng)數(shù)值計(jì)算方法發(fā)展走向成熟的一個(gè)重要標(biāo)志是：成熟的方法越來越多地應(yīng)用于工程實(shí)際問題中，商業(yè)化通用軟件包不斷出現(xiàn)4。一個(gè)商業(yè)化軟件包通常由下面幾部分組成： 以上三部分中前、后處理占用了軟件包語(yǔ)句的90%以上，編程的主要工作量在此，而

8、數(shù)據(jù)處理，也就是我們目前正在學(xué)習(xí)的數(shù)值計(jì)算方法僅占軟件語(yǔ)句的10%以內(nèi)，但它卻是占用計(jì)算機(jī)內(nèi)存量和消耗CPU時(shí)間的主要部分。2 電磁場(chǎng)數(shù)值計(jì)算方法的分類求解電磁問題的最終要求就是獲得滿足實(shí)際條件的Maxwell方程的解，借助于計(jì)算數(shù)學(xué)中的數(shù)值算法能夠得到大多數(shù)電磁問題的近似解。數(shù)值算法的基本思想5就是把連續(xù)變量函數(shù)離散化，把微分方程化為差分方程；把積分方程化為有限和的形式，從而建立起收斂的代數(shù)方程組，然后利用計(jì)算機(jī)技術(shù)進(jìn)行求解。數(shù)值計(jì)算方法從求解方程的形式看，主要分為積分方程法和微分方程法兩大類。積分方程法主要有矩量法和邊界元法，微分方程法主要有有限差分法和有限元法。對(duì)兩種方程法的比較，如表一

9、所示。表一 積分方程法和微分方程法的比較積分方程法微分方程法共性對(duì)場(chǎng)問題的處理是一致的，即需離散化場(chǎng)域，結(jié)果是數(shù)值解不 同 點(diǎn)離散域僅在場(chǎng)源區(qū)，無需對(duì)整個(gè)場(chǎng)域離散整個(gè)場(chǎng)域計(jì)算對(duì)象場(chǎng)量先求位函數(shù)，再求場(chǎng)量求解域可在場(chǎng)域內(nèi)某一局部區(qū)域求解，也可在全場(chǎng)域內(nèi)求解全場(chǎng)域內(nèi)求解計(jì)算程度較高較低應(yīng)用不適用邊界區(qū)域復(fù)雜的場(chǎng)域邊界形狀復(fù)雜的場(chǎng)域較易處理聯(lián)系兩種方法的結(jié)合形式，可處理較復(fù)雜的電磁場(chǎng)問題3 幾種重要的數(shù)值計(jì)算方法3.1 有限差分法在電磁場(chǎng)數(shù)值計(jì)算方法中，有限差分法是應(yīng)用最早的一種方法。有限差分法以其概念清晰，方法簡(jiǎn)單、直觀，有大致固定的處理和計(jì)算模式，具有一定的通用性等特點(diǎn)，在電磁場(chǎng)數(shù)值分析領(lǐng)域內(nèi)得

10、到了廣泛的應(yīng)用。3.1.1 有限差分法的基本原理有限差分法的基本思想是把連續(xù)的定解區(qū)域用有限個(gè)離散點(diǎn)構(gòu)成的網(wǎng)格來代替，這些離散點(diǎn)稱作網(wǎng)格的節(jié)點(diǎn)；把連續(xù)定解區(qū)域上的連續(xù)變量的函數(shù)用在網(wǎng)格上定義的離散變量函數(shù)來近似；把原方程和定解條件中的微商用差商來近似，積分用積分和來近似，于是原微分方程和定解條件就近似地代之以代數(shù)方程組，即有限差分方程組，解此方程組就可以得到原問題在離散點(diǎn)上的近似解。然后再利用插值方法便可以從離散解得到定解問題在整個(gè)區(qū)域上的近似解。 3.1.2 差分與差商設(shè)函數(shù)的自變量有一小增量，則的增量為 （3.1）為函數(shù)的一階差分。當(dāng)增量足夠小，差分與微分之間的差才足夠小。一階差分是自變量

11、的函數(shù)。按式（1）,計(jì)算的差分稱二階差分，且 （3.2）函數(shù)的一階導(dǎo)數(shù)為 應(yīng)用差分，可表示為 （3.3）故可表示為差分除以有限小差分的商，稱為差商。 同理，函數(shù)的二階導(dǎo)數(shù)可表示為 （3.4）3.1.3 差分方程的構(gòu)造現(xiàn)以二維靜態(tài)電、磁場(chǎng)泊松方程的第一類邊值問題為例，來具體闡明有限差分法的應(yīng)用。設(shè)具有平行平面場(chǎng)特征的電磁場(chǎng)場(chǎng)域,如圖1所示,為一由閉合邊界所界定的平面域，其定解條件可表述為 （3.5） （3.6） 對(duì)于所給定的偏微分方程定解問題，應(yīng)用有限差分法，首先需從網(wǎng)格剖分著手決定離散點(diǎn)的分布方式。原則上，可以采用任意的網(wǎng)格剖分方式，但這將直接影響所得差分方程的具體內(nèi)容，進(jìn)而影響解題的經(jīng)濟(jì)性與

12、計(jì)算精度。為簡(jiǎn)化問題，通常采用完全有規(guī)律的分布方式，這樣在每個(gè)離散點(diǎn)上就能得出相同形式的差分方程，有效地提高解題速度，因而經(jīng)常采用正方形的網(wǎng)格的剖分方式。現(xiàn)即以這種正方形網(wǎng)格剖分場(chǎng)域，也就是說，用分別與兩坐標(biāo)軸平行的兩簇等距網(wǎng)格線來生成正方形網(wǎng)格，即xL03MDyhh124圖1 正方形網(wǎng)格劃分 為步長(zhǎng)，網(wǎng)格線的交點(diǎn)稱為節(jié)點(diǎn)，這樣域就離散化為由網(wǎng)格節(jié)點(diǎn)標(biāo)成的離散點(diǎn)得集合。對(duì)場(chǎng)域中節(jié)點(diǎn)是一典型節(jié)點(diǎn)，它與周圍的1，2，3和4點(diǎn)構(gòu)成一個(gè)對(duì)稱星型。設(shè)這些離散點(diǎn)上待求函數(shù)的近似值記為 ，則式（6）可近似離散化為（3.7）即 （3.8） 若式（6）=0,則節(jié)點(diǎn)上函數(shù)的值等于其四周相鄰點(diǎn)函數(shù)值的平均。因?yàn)椴罘?/p>

13、方程（7），（8）只出現(xiàn)待求函數(shù)在點(diǎn)及其四個(gè)臨近點(diǎn)的值，故稱之為五點(diǎn)差分格式6，根據(jù)差分方程組解出各離散點(diǎn)處的待求函數(shù)值。3.2 有限元法傳統(tǒng)的變分法在20世紀(jì)二三十年代為其新型時(shí)期，理論上發(fā)展很快，各種變分問題的最后求解都可歸結(jié)為解尤拉方程的邊值問題，然而只有在一些特殊情況下尤拉方程才能求出精確解，在大多數(shù)情況下，尤拉方程的精確解無法求出。四五十年代，隨著計(jì)算機(jī)的出現(xiàn)，使其在實(shí)際應(yīng)用中逐漸為比較靈活、通用的有限差分法所替代。但是，有限差分法在理論上沒有以變分原理為基礎(chǔ)，因而其收斂性和數(shù)值穩(wěn)定性往往得不到保證。隨后發(fā)展形成的有限元法正是變分法與有限差分法相結(jié)合的成果，它取長(zhǎng)補(bǔ)短地在理論上以變分

14、原理為基礎(chǔ)，在具體方法構(gòu)造上又利用了有限差分法網(wǎng)格離散化處理的思想7。3.2.1 有限元法的基本原理dxdsA(x1, y1)B(x2, y2)xy圖2 最速降線問題O有限元法是以變分原理為基礎(chǔ)，將要求解的微分方程型數(shù)學(xué)模型邊值問題，首先轉(zhuǎn)化為相應(yīng)的變分問題，即泛函求極值問題；然后，利用剖分插值將變分問題離散化為普通多元函數(shù)的極值問題，最終歸結(jié)為一組多元的代數(shù)方程組，求解該方程組，從而獲得邊值問題的數(shù)值解。3.2.2 泛函、變分問題簡(jiǎn)介 在微積分學(xué)形成初期，以數(shù)學(xué)物理問題為背景，與多元函數(shù)的極值問題相對(duì)應(yīng)，已在幾何、力學(xué)上提出了若干個(gè)求解泛函極值的問題。如圖2中的質(zhì)點(diǎn)最速降線問題所述，質(zhì)點(diǎn)A從

15、定點(diǎn)自由下滑到定點(diǎn)B，試求使滑行時(shí)間最短的質(zhì)點(diǎn)下滑軌道。圖示滑行弧段所需時(shí)間為 滑行總時(shí)間為 （3.9）（3.9）式不僅取決于積分端點(diǎn)和，而且取決于的選取。取決于，所以是函數(shù)的函數(shù)，稱之為的泛函，記作。于是所述之最速降線問題，在數(shù)學(xué)上就歸結(jié)為研究泛函的極值問題，即 （3.10）泛函的極值問題就稱為變分問題。對(duì)一般問題而言，可導(dǎo)出下列對(duì)應(yīng)于一個(gè)自變量、單個(gè)函數(shù)及其導(dǎo)數(shù)的已知函數(shù) （3.11）式中為、和的已知函數(shù)。泛函的自變量不是一般的自變量，而是一個(gè)或幾個(gè)函數(shù)所屬的函數(shù)族。在端點(diǎn)和上分別等于給定值的無數(shù)個(gè)函數(shù)中，僅有一個(gè)能使定積分達(dá)到極小值，此函數(shù)稱為極值函數(shù)。因此，變分問題就在于尋求使泛函達(dá)到

16、極值的該極值函數(shù)，即分析研究泛函的極值問題。3.2.3 泛函的變分與尤拉方程泛函變分問題的經(jīng)典解法有兩種，一種稱之為直接解法，另一類是間接解法。直接解法是直接把泛函的極值問題近似地轉(zhuǎn)化為一般多元函數(shù)的極值問題，用有限維子空間中的函數(shù)去逼近無窮維空間中的極值函數(shù)，從而近似求得泛函的極值。間接解法是將變分問題轉(zhuǎn)化為尤拉方程（微分方程）的定解問題，即邊值問題來求解。以式（3.11）這種最簡(jiǎn)形式來推導(dǎo)尤拉方程。設(shè)函數(shù)稍有變化，記作，稱之為的變分，它反映了整個(gè)函數(shù)的變化量。這樣泛函的值也應(yīng)隨之變動(dòng)，相應(yīng)于變分的泛函增量為 （3.12）將（3.12）式由多元函數(shù)的泰勒公式展開 （3.13）式中作為泛函增量

17、的線性主部為 （3.14）稱為泛函的一次變分（簡(jiǎn)稱變分）。而、分別是函數(shù)變分及其導(dǎo)數(shù)的二次、三次齊次式等的積分，依次稱為二次變分，三次變分令變分問題的解為，且設(shè)極值解稍有變動(dòng)，令 （3.15）式中為任意給定的微量實(shí)參數(shù)，值就確定了函數(shù)族中的某一曲線，進(jìn)而確定泛函之值；而是定義于區(qū)間且滿足齊次邊界條件的可微函數(shù)。于是泛函=就成為變量的函數(shù)，且當(dāng)時(shí)獲極值函數(shù)的解。在時(shí)取得極值的必要條件是 （3.16）利用分部積分，并根據(jù)變分與微分順序可互換原理，(3.16）式可寫為 （3.17）由于（3.17）對(duì)任意均成立，故有 （3.18）方程（3.18）就稱為泛函（3.11）的極值問題的尤拉方程。綜上所述，有

18、限元法的基本特點(diǎn)是：（1）離散化過程保持了明顯的物理意義。這是因?yàn)椋兎衷砻枋隽酥湮锢憩F(xiàn)象的物理學(xué)中的最小作用原理（如力學(xué)中的最小勢(shì)能原理、靜電學(xué)中的湯姆遜定理等）。（2）優(yōu)異的解題能力。 與其他數(shù)值計(jì)算方法相比較，有限元法在適應(yīng)場(chǎng)域邊界幾何形狀，以及媒質(zhì)物理性質(zhì)變異情況復(fù)雜的問題求解上，有突出的優(yōu)點(diǎn)。（3）從數(shù)學(xué)理論意義上講，有限元法作為應(yīng)用數(shù)學(xué)的一個(gè)重要分支，很少有其他方法應(yīng)用的這樣廣泛。它使微分方程的解法和理論面目一新，推動(dòng)了泛函分析與計(jì)算方法的發(fā)展。3.3 矩量法矩量法，是近年來在天線、微波技術(shù)和電磁波散射等方面廣泛應(yīng)用的一種方法。從這些實(shí)際問題涉及開域、激勵(lì)場(chǎng)源分布形態(tài)較為復(fù)雜等

19、特征出發(fā)，矩量法是將待求的積分方程問題轉(zhuǎn)化為一個(gè)矩陣方程的問題7，借助于計(jì)算機(jī)，求得其數(shù)值解，從而在所得激勵(lì)源分布的數(shù)值解基礎(chǔ)上，即可算出輻射場(chǎng)的分布及其波阻抗等特性參數(shù)。3.3.1矩量法的基本原理先選定基函數(shù)對(duì)未知函數(shù)進(jìn)行近似展開，代入算子方程，再選取適當(dāng)?shù)臋?quán)函數(shù)，使在加權(quán)平均的意義下方程的余量等于0，由此將連續(xù)的算子方程轉(zhuǎn)換為代數(shù)方程。原則上，矩量法可用于求解微分方程和積分方程，但用于微分方程時(shí)所得到的代數(shù)方程組的系數(shù)矩陣往往是病態(tài)的，故在電磁場(chǎng)問題中主要用于求解積分方程。3.3.2 加權(quán)余量法設(shè)給定邊值問題的場(chǎng)方程統(tǒng)一表述為如下的算子方程，即 （3.19）已知邊界條件為 （3.20） （

20、3.21）其中: 是線性算子,是待求函數(shù), 是已知的源。若為精確解，則方程(3.19)和邊界條件(3.20),(3.21)應(yīng)該完全滿足。但大多數(shù)情況下，不能得到的精確解，只能通過數(shù)值方法進(jìn)行估計(jì)。構(gòu)造一個(gè)由有限個(gè)線性無關(guān)函數(shù)(=1,2,)所組成的基函數(shù)集合，借以展開待求函數(shù)的近似解為 （3.22）將代入式(3.19)中必然存在誤差，即 (3.23) 取一個(gè)歸屬于試探函數(shù)的權(quán)函數(shù)集合，令 (=1,2, ,) （3.24）式(3.24)由個(gè)方程構(gòu)成的方程組，它等價(jià)于人為地強(qiáng)制近似解，使其因不能精確地滿足場(chǎng)方程而導(dǎo)致的誤差在平均的含義上等于零。按式(3.24)展開，所構(gòu)成的各種求解積分或微分方程近似

21、解的方法可被統(tǒng)稱為加權(quán)余量法8。因?yàn)榘唇o定權(quán)函數(shù)展開式的式(3.24)，即意味著余量對(duì)取矩的一組平衡式，故式(3.24)的構(gòu)造亦就被稱為矩量法。基于加權(quán)余量式(3.24)，進(jìn)行移項(xiàng)處理，便得 (=1,2, ,) （3.25） 將式(3.22)代入式(3.25)的左端，有 （3.26）為了書寫方便，令 和代入(3.26)式，則可寫成 （3.27）這樣，即展開成含個(gè)未知數(shù)的個(gè)方程。若用矩陣形式表示，則有 （3.28）綜上所述，矩量法的特點(diǎn)是：矩量法將連續(xù)方程離散化為代數(shù)方程組, 既適用于求解微分方程, 又適用于求解積分方程。它的求解過程簡(jiǎn)單, 求解步驟統(tǒng)一, 應(yīng)用起來比較方便，然而需要一定的數(shù)學(xué)技

22、巧, 如離散化的程度、基函數(shù)與權(quán)函數(shù)的選取, 矩陣求解過程等。另外必須指出的是, 矩量法可以達(dá)到所需要的精確度、解析部分簡(jiǎn)單, 可計(jì)算量很大, 即使用高速大容量計(jì)算機(jī), 計(jì)算任務(wù)也很繁瑣。4 電磁場(chǎng)數(shù)值計(jì)算的應(yīng)用前景電磁場(chǎng)數(shù)值計(jì)算方法近六十年來發(fā)展如此之快，除由于從事這方面的科研人員的努力之外，主要是其研究成果迅速被電機(jī)、電器、變壓器、加速器、微波器件、計(jì)算機(jī)磁頭等領(lǐng)域采用9，對(duì)改善產(chǎn)品性能、降低生產(chǎn)成本，起著越來越大的作用。眾所周知，任何電磁器件，包括國(guó)民經(jīng)濟(jì)中應(yīng)用極其廣泛的電機(jī)與變壓器，其能量轉(zhuǎn)換都是通過電磁場(chǎng)來實(shí)現(xiàn)的。但傳統(tǒng)的電磁產(chǎn)品的設(shè)計(jì)方式由于客觀條件與手段的限制，把場(chǎng)的實(shí)際分布參數(shù)

23、當(dāng)作集中參數(shù)處理，不可避免地帶來相當(dāng)大的誤差。在迫不得已的情況下，只能用模擬、實(shí)驗(yàn)等方法處理，其耗費(fèi)大、周期長(zhǎng)，可借鑒的經(jīng)驗(yàn)不多，而實(shí)驗(yàn)的模擬也不都是有條件的，此外，工農(nóng)業(yè)和日常生活中所用電磁裝置越來越多，生產(chǎn)和銷售競(jìng)爭(zhēng)劇烈，因此有效的設(shè)計(jì)方法顯然受到重視。在采用電磁場(chǎng)數(shù)值計(jì)算以前的任何方法，即使是十分精巧的代數(shù)解析方法，也只能適用于特別簡(jiǎn)單的幾何結(jié)構(gòu)以及一些特殊的假定模型，有時(shí)模型甚至簡(jiǎn)化到不能容忍的程度，但除此之外別無它法。因?yàn)閷?shí)際的電磁場(chǎng)問題，特別是電機(jī)電磁場(chǎng)，其邊界情況十分復(fù)雜，加上鐵的飽和以及導(dǎo)體中的渦流效應(yīng)，解析方法是無能為力的。而數(shù)值方法可以模擬復(fù)雜的形狀，可以適應(yīng)非線性問題，可

24、以進(jìn)行渦流的分布計(jì)算。電磁場(chǎng)的數(shù)值計(jì)算與流體力學(xué)分析、溫度分析、機(jī)械應(yīng)力分析、電磁力分析和生產(chǎn)計(jì)劃等等有機(jī)聯(lián)系起來，由計(jì)算機(jī)完成全部設(shè)計(jì)，構(gòu)成所謂CAE10系統(tǒng)(Computer Aided Engineering)。CAE系統(tǒng)基本可包括設(shè)計(jì)與制造的全部過程，產(chǎn)品的設(shè)計(jì)、制造方案、準(zhǔn)備零件、草圖、計(jì)算、成本核算、生產(chǎn)、機(jī)床數(shù)控和試驗(yàn)、模擬和產(chǎn)品的自動(dòng)測(cè)試都可能列入到CAE系統(tǒng)中。目前CAE系統(tǒng)在西方國(guó)家發(fā)展很快，估計(jì)每年增長(zhǎng)40%，尤其在機(jī)械行業(yè)中，形成了所謂CADMAT即計(jì)算機(jī)輔助設(shè)計(jì)、加工與試驗(yàn)。從長(zhǎng)遠(yuǎn)觀點(diǎn)看，電磁場(chǎng)數(shù)值計(jì)算要發(fā)揮更實(shí)際的作用，必須與其他多種學(xué)科相結(jié)合以形成CAE系統(tǒng)，CA

25、E的發(fā)展必然給工業(yè)結(jié)構(gòu)上帶來巨大的變化。結(jié)束語(yǔ)數(shù)值計(jì)算是一門計(jì)算的藝術(shù)，電磁場(chǎng)的數(shù)值計(jì)算橫跨了多個(gè)學(xué)科，是數(shù)學(xué)理論、電磁理論和計(jì)算機(jī)應(yīng)用能力的完美組合。通過本次論文設(shè)計(jì)，我對(duì)差分、變分、泛函、加權(quán)余量法等數(shù)學(xué)知識(shí)有了深刻的認(rèn)識(shí)，掌握了有限差分法、有限元法、矩量法的基本原理，并能進(jìn)行簡(jiǎn)單的運(yùn)算。但若用數(shù)值計(jì)算方法解決復(fù)雜的電磁場(chǎng)問題，則在很大程度上依賴于數(shù)學(xué)知識(shí)及計(jì)算機(jī)編程能力，這就需要進(jìn)一步系統(tǒng)學(xué)習(xí)相關(guān)的理論知識(shí)。參考文獻(xiàn)1 文舸一計(jì)算電磁學(xué)的進(jìn)展與展望J電子學(xué)報(bào)，1995，23(10)：62-692 劉圣民電磁場(chǎng)的數(shù)值方法M武漢：華中理工大學(xué)出版社，1991：23-453 鄧東皋，彭立中小波

26、分析J數(shù)學(xué)進(jìn)展，1991， 20(3)： 294-3104 洪偉計(jì)算電磁學(xué)研究進(jìn)展J東南大學(xué)學(xué)報(bào)（自然科學(xué)版），2002，32（3）：335-3395 倪光正，楊仕友，錢秀英等工程電磁場(chǎng)數(shù)值計(jì)算M北京：機(jī)械工業(yè)出版社，2003： 123-1616 盛劍霓工程電磁場(chǎng)數(shù)值分析M西安：西安交通大學(xué)出版社，1991： 69-757 方靜，汪文秉有限元法和矩量法結(jié)合分析背腔天線的輻射特性J微波學(xué)報(bào)，2000，16(2)：139-1438 連漢雄電磁場(chǎng)理論的數(shù)學(xué)方法M北京：北京理工大學(xué)出版社， 1990：10-15 9 張平文，劉法啟，張宇小波函數(shù)值的計(jì)算J計(jì)算數(shù)學(xué)，1995(3)：173-18510 樓

27、仁海，符果行，袁敬閎電磁理論M成都：電子科技大學(xué)出版社，1996：73-102 Electromagnetic numerical method Department of Physics 0702 Student Du XingXing Tutor Ren LiYing Abstract：Numerical calculation, which is widely used in many fields, such as electric, military affairs, economy, ecology, medical treatment, astronomy, geology and so on is a method of getting the numerical approximate solution in the studying and solving of mathematical problems. This paper reviews the development history and classification of the electromagnetic numerical calculation, and explicitly introduces three typical numerical calculati