艦船磁場數(shù)值計算方法發(fā)展綜述

王德強(qiáng)，余 強(qiáng)

(1.中國人民解放軍92786部隊，廣東 湛江 520464;2.海軍工程設(shè)計研究局，北京 100070)

0 引言

現(xiàn)代艦船大多由鋼鐵材料建造，在地磁場和機(jī)械力的綜合作用下，成為一個巨大的浮動磁體。艦船在地球磁場的磁化下，其周圍空間具有一定的磁場，這個磁場被稱為艦船磁場。磁性水雷、磁性魚雷等武器的引信感應(yīng)到艦船磁場會在艦船附近爆炸，對艦船生命力構(gòu)成嚴(yán)重威脅，因此，艦船磁隱身是艦船隱身性能的重要指標(biāo)之一［1］。現(xiàn)代戰(zhàn)爭中利用艦船磁場對其實施探測和攻擊的武器發(fā)展十分迅速，如果不對艦船實施有效的磁防護(hù)措施將嚴(yán)重影響其生命力和戰(zhàn)斗力。

為提高艦船生命力，必須對艦船實施磁性防護(hù)。艦船磁性防護(hù)的目的是盡一切可能降低艦船的磁性特征，使其在周圍空間一定范圍內(nèi)產(chǎn)生的磁場小于磁引信武器或磁探測儀器所設(shè)定的動作閥值。艦船磁性防護(hù)的前提是全面、準(zhǔn)確地了解和掌握艦船空間磁場分布，而衡量磁性防護(hù)的指標(biāo)也必須了解艦船的空間磁場分布。

目前，一般通過磁傳感器測量和磁場數(shù)值建模來獲得艦船磁場。隨著電磁場數(shù)值計算技術(shù)的發(fā)展，可用于艦船磁場數(shù)值建模的方法已趨多元化［2］。鑒于艦船磁場數(shù)值建模的復(fù)雜性，如何建立一個精度和穩(wěn)定性都滿足艦船磁性防護(hù)需求的艦船磁場數(shù)值模型成為磁隱身技術(shù)中的一個重要課題。本文在廣泛查閱國內(nèi)外文獻(xiàn)的基礎(chǔ)上，對艦船磁場數(shù)值建模方法進(jìn)行系統(tǒng)總結(jié)，并分析和評價目前用于艦船磁場數(shù)值建模中等效源法、廣義多級技術(shù)、磁場積分法和磁場微分法等幾種典型方法。

1 艦船磁場分類及其特點

按照艦船磁場形成機(jī)理，可將艦船磁場分為感應(yīng)磁場和固定磁場兩部分。感應(yīng)磁場隨地磁場變化，而固定磁場在一般情況下保持相對穩(wěn)定。

1.1 感應(yīng)磁場

鋼質(zhì)艦船在地球磁場作用下，會產(chǎn)生感應(yīng)磁性，進(jìn)而在艦船周圍產(chǎn)生感應(yīng)磁場，感應(yīng)磁場是一種隨外界磁場 (主要是地球磁場)變化的磁場，它與艦船實時的位置、航向、姿態(tài)等因素有關(guān)，因而是一種瞬時效應(yīng)磁場。艦船感應(yīng)磁性的大小與下列因素有關(guān):

1)艦船航行緯度區(qū)地磁場的大小;2)艦船的航向、縱傾和搖擺;3)造船所用鋼材的磁性能;4)艦船形狀、尺寸及鐵磁設(shè)備在船上的分布情況。同一類型的艦船如果建造材料相同，結(jié)構(gòu)相同，則其感應(yīng)磁性相同。

1.2 固定磁場

船體由許多鋼材鉚焊而成，船上安裝了大量的鋼鐵設(shè)備，造船鋼板和鐵磁設(shè)備本身都具有一定固定磁性。艦船建造組裝后，這些鋼板和設(shè)備的固定磁性迭加起來就形成了艦船總的固定磁性的一部分。這些鋼板和設(shè)備在艦船建造過程中一直處于地球磁場中，在施工中受到長期機(jī)械震動，也必然會產(chǎn)生一定程度的不可逆磁化，從而形成固定磁性的另一部分。另外，在鋼板的焊接、切割等局部性的加熱、冷卻時鋼材內(nèi)部將產(chǎn)生應(yīng)力，這些應(yīng)力也會改變局部的磁特性。這樣，在艦船下水時就已經(jīng)形成了自己所特有的固定磁性。

艦船的固定磁性主要在艦船建造時期形成，它與下列因素有關(guān):

1)造船所用鋼材的磁特性。鋼材的磁特性不同，形成的固定磁性也將不同。

2)造船地區(qū)地磁場的大小。地磁垂向分量大的地方，形成較大的垂向固定磁性，地磁水平分量大的地方，形成較大的水平固定磁性。

3)造船船臺的方向。它決定了地磁場對艦船的水平磁化方向。

4)船體及主要鐵磁設(shè)備的形狀。

5)造船工藝等。

艦船在建造過程中，由于上述各種情況不完全相同，因而不同艦船甚至是同類型的艦船其固定磁性也將不同。

艦船固定磁性的固定不變相對于外界的磁化條件不發(fā)生大的變化而言。當(dāng)艦船的磁化條件有較大變化時，其固定磁性將會發(fā)生變化。艦船建造期間所形成的固定磁性在首次航行時就會發(fā)生一定的變化。這是因為艦船在建造時一直受到一個方向的磁化，而首次航行時艦船受到各不同方向的磁化，加上波浪沖擊和機(jī)器開動時振動等因素影響，使艦船在新的條件下達(dá)到新的穩(wěn)定磁化狀態(tài)。艦船經(jīng)過戰(zhàn)斗、炮火射擊，受到強(qiáng)烈震動，大風(fēng)浪襲擊以及船體經(jīng)過較大規(guī)模的修理之后，固定磁性也要發(fā)生變化。另外艦船更換基地，較長時間航行在另一個地區(qū)，其地磁場與原地區(qū)有明顯的差異，也會使固定磁性發(fā)生變化。

因此，由固定磁性在艦船周圍產(chǎn)生的固定磁場是在鋼材的磁滯效應(yīng)作用下由鋼材的歷史磁化形成的磁場，是一種相對穩(wěn)定的積累效應(yīng)磁場。

2 艦船磁場數(shù)值計算的特點

由于艦船的磁場是受地球磁場這個基本穩(wěn)定的磁場磁化產(chǎn)生的，則艦船在地球磁場中的磁化可以看作在弱磁場作用下鐵磁物質(zhì)的靜態(tài) (準(zhǔn)靜態(tài))磁化過程，而且其感應(yīng)磁場部分可以看作地磁場作用下的可逆磁化。另外，所研究的區(qū)域主要是艦船的外部區(qū)域，所以，可以將艦船磁場的數(shù)值計算歸屬到開域靜磁場問題的研究范疇。

簡而言之，艦船磁場數(shù)值建模可等效為考慮剩余磁化強(qiáng)度時鐵磁物體在外磁場作用下的靜磁建模問題。因此，電磁場數(shù)值技術(shù)中的靜磁建模方法均可用來進(jìn)行艦船磁場數(shù)值建模，但由于艦船結(jié)構(gòu)、尺寸和設(shè)備的復(fù)雜性及多樣性，將電磁場數(shù)值計算方法直接應(yīng)用于艦船磁場數(shù)值建模存在一定困難，艦船磁場數(shù)值建模的不利因素為:

1)艦船磁場數(shù)值建模為開域問題;

2)復(fù)雜的物理結(jié)構(gòu)給幾何建模和剖分增加了難度;

3)鐵磁設(shè)備和艦船殼體的尺寸相差很大，甚至達(dá)到10－4;

4)材料的多樣性給鐵磁材料磁性參數(shù)的確定增加了困難;

5)考慮剩余磁化強(qiáng)度的艦船磁場數(shù)值建模技術(shù)存在模型穩(wěn)定性問題;

6)剖分單元較多的艦船幾何網(wǎng)格模型對計算時間提出了更高的要求。

上述不利因素限制了不同電磁場數(shù)值計算方法在艦船磁場數(shù)值建模技術(shù)中的應(yīng)用，如有限元法在處理開域問題時需對區(qū)域邊界進(jìn)行特殊處理、磁場積分法形成的系數(shù)矩陣為滿陣求解時需花費(fèi)大量的時間、等效源法存在建模穩(wěn)定性問題等［3］。

然而，艦船磁場數(shù)值建模也存在有利因素。在一定工程精度要求下，艦船磁場數(shù)值建模大多數(shù)情況不需考慮電流線圈的作用，激勵源僅為地磁場，與強(qiáng)電流產(chǎn)生的磁場相比，地磁場為弱磁場，弱磁作用下的磁化點處于鐵磁材料的線性磁化區(qū)，因此鐵磁材料磁性參數(shù)的線性化使得艦船磁場數(shù)值建模問題可大大簡化。

由于艦船磁場數(shù)值建模所具有的獨(dú)特特點，通用商業(yè)軟件如Ansys、Ansoft等直接應(yīng)用于艦船磁場建模相對困難，特別是考慮剩余磁化強(qiáng)度的磁場建模，采用磁場逆問題的思路來研究需對軟件進(jìn)行二次開發(fā)。因此，國內(nèi)外學(xué)者都從艦船磁場數(shù)值建模的特點出發(fā)，基于電磁場理論來對艦船磁場數(shù)值建模技術(shù)展開廣泛的研究。

3 艦船磁場數(shù)值計算方法

由國內(nèi)外艦船磁場數(shù)值建模技術(shù)資料來看，應(yīng)用較為廣泛的建模方法有等效源法、廣義多級技術(shù)、磁場積分法和磁場微分法。

3.1 等效源法

所謂等效源法［3］，即是在求解的場域之外，用一組虛設(shè)的模擬源來代替原問題中比較復(fù)雜的磁化物體，使模擬源在求解區(qū)域產(chǎn)生的場來等效實際磁化物體產(chǎn)生的場，如圖1所示。

圖1 等效源法求解艦艇磁場數(shù)值建模示意圖Fig.1 Equivalent source method for calculating warships magnetic field numerical modeling

在磁場等效源法中，常用的模擬體為磁荷、磁偶極子、均勻磁化橢球、磁化電流等。這些模擬體的分布參數(shù)一般需根據(jù)對問題的定性分析及計算經(jīng)驗人為確定，模擬體的磁性參數(shù)則由已知邊界上的情況或測量值來確定。當(dāng)模擬體的分布參數(shù)和磁性參數(shù)確定后，求解區(qū)域內(nèi)任意點的磁場值就可解析出來。

由于等效源法具有實施簡單靈活、編程易于實現(xiàn)、計算時間短且具有一定工程實用精度等特點，磁場等效源法已廣泛應(yīng)用于艦船磁場數(shù)值建模中，如艦船磁異常在不同深度的換算。然而等效源法在應(yīng)用中也表現(xiàn)出其不足之處，如模擬體的分布參數(shù)直接決定著等效源法計算精度及模型穩(wěn)定性，這給艦船磁場數(shù)值建模的通用化、實用化帶來了困難。目前一般通過以下措施增強(qiáng)等效源法的通用性:對模擬體分布參數(shù)優(yōu)化可減少人為因素的干擾，有效提高磁場模型的穩(wěn)定性［4］;采用正則化技術(shù)可降低測量噪聲對建模結(jié)果的影響［5］。

3.2 廣義多級技術(shù)

20世紀(jì)80年代中期，一種數(shù)值法和解析法的結(jié)合——半解析法應(yīng)運(yùn)而生，即廣義多級技術(shù)，也稱為級數(shù)法、多級子法、多級理論、圓 (球)形等效源法、新型等效源法［6］。廣義多級技術(shù)對計算機(jī)資源要求較少，既不需要離散場域，也不需要離散邊界，而且對于某一類問題可用統(tǒng)一格式獲得高階可導(dǎo)的連續(xù)解。其基本原理是:將偏微分方程中的未知函數(shù)，用一系列滿足該偏微分方程的本征函數(shù)特解展開，確定該通解中待定系數(shù)的方法采用點匹配法，并根據(jù)經(jīng)驗總結(jié)出一些使用規(guī)則。

從國內(nèi)外文獻(xiàn)可以看出，廣義多級技術(shù)已用來解決三維靜態(tài)電磁場問題［7－9］，且國外已將其應(yīng)用于艦船磁場數(shù)值建模技術(shù)領(lǐng)域，如基于磁多級子成像技術(shù)來重構(gòu)艦船磁性狀態(tài)等［10－12］。

最近，美國科學(xué)家Alexander V.Kildishev等人研究了一系列潛艇磁場多極成像方面的課題，主要是將磁源的標(biāo)量磁勢在扁平旋轉(zhuǎn)橢球體坐標(biāo)系下展成諧波形式，其諧波系數(shù)即描述了磁源的多極特征，然后用外部磁場的測量數(shù)據(jù)來求得諧波系數(shù)，進(jìn)而可以得到用級數(shù)表示的磁場多極展開式。這些研究在潛艇外部磁場的分析方面取得一定進(jìn)展，但是建立在假設(shè)艇體為扁旋轉(zhuǎn)橢球體近似的基礎(chǔ)之上，對于水面艦船的分析還沒有相關(guān)資料。

從總結(jié)出的規(guī)律看，由于缺乏理論指導(dǎo)，使求解結(jié)果受制于使用者的經(jīng)驗，因此，廣義多級技術(shù)在艦船磁場數(shù)值建模中還沒有得到廣泛的應(yīng)用。

3.3 磁場積分法

磁場積分法主要包括積分方程法和邊界積分法兩類。

1)積分方程法

就積分方程法而言，它從宏觀的角度來描述場，場區(qū)中每點場值僅取決于所有場源對它的影響，即一次場源 (激勵源)和二次場源 (磁化源)，再通過畢奧－薩伐定律將場點和源點聯(lián)系在一起，建立以場點場量為未知量的代數(shù)方程組，求解上述該方程組即可獲得整個求解區(qū)域的解，如圖2所示。

圖2 積分方程法建模示意圖Fig.2 Integral-equation methods for modeling

選擇不同的場點場量，如磁化強(qiáng)度、磁場強(qiáng)度、磁標(biāo)勢 (簡化勢和全標(biāo)勢)、磁矢勢等，可得到不同的代數(shù)方程。根據(jù)特定的問題，選擇不同場點場量作為求解變量可有效簡化問題的求解。由于積分方程法將磁化介質(zhì)視為二次場源，所以其具有離散僅限于磁性介質(zhì)區(qū)、便于處理開域問題且不需考慮邊界條件的優(yōu)點［13］。

國內(nèi)外資料表明，積分方程法已在艦船磁場數(shù)值建模技術(shù)中得到了廣泛的應(yīng)用，如法國LEG公司以閉環(huán)消磁技術(shù)為對象，將積分方程法成功應(yīng)用于求解艦船薄殼體的磁場數(shù)值建模［14］;我國也對積分方程法在艦船磁場數(shù)值建模中的應(yīng)用展開了較為廣泛的研究，且取得了較多成果［15－16］。

積分方程法從宏觀的角度來描述場，場區(qū)中每點的磁場值取決于所有場源對它的影響，場點和源點的聯(lián)系是通過畢奧－薩伐定律實現(xiàn)的。離散只需在非線性鐵磁區(qū)域中進(jìn)行即可，這使得數(shù)據(jù)輸入和網(wǎng)格剖分大為簡化，所以它能夠較好地解決開域問題及連續(xù)場的計算問題。但是，積分方程法雖很易于處理開域問題，但也有不足之處，用積分方程法求解非線性問題時，考慮到船體比較大，離散的單元可能比較多，而且確定物質(zhì)磁化狀態(tài)的離散方程的系數(shù)矩陣是非對稱滿陣，加上每一剖分單元重心上的場參數(shù)是用向量來描述的，在三維場中每一單元重心上形成3個未知數(shù)，因此需要相當(dāng)大的內(nèi)存來存儲系數(shù)矩陣的元素。此外，單元間耦合系數(shù)的求取也要消耗大量機(jī)時，而且鐵區(qū)中的介質(zhì)參數(shù)要反復(fù)迭代求解，要耗占大量的CPU時間和內(nèi)存，當(dāng)艦船復(fù)雜鐵磁材料離散單元數(shù)越多時，其求解時間和所需內(nèi)存將成指數(shù)增加，所以一般只能兼顧計算精度和計算時間來取適當(dāng)?shù)碾x散單元數(shù)，這也一定程度上制約了積分方程法的應(yīng)用。在艦船磁場數(shù)值計算中若單獨(dú)應(yīng)用積分方程法最好采用一些方法來對其進(jìn)行改進(jìn)，比如利用小波矩量法和快速多極技術(shù)等方法。

2)邊界積分法

邊界積分法也稱邊界元法。在場域內(nèi)媒質(zhì)為各向同性、線性和均勻的條件下，可得直接邊界積分法和間接邊界積分法 (包括單層源法和雙層源法)兩類;在分域媒質(zhì)分界面上可得分界面上的邊界積分方程［3］。直接邊界元法從等效原理出發(fā)，將閉合分界面內(nèi)的源在場點的磁效應(yīng)用閉合分界面上的等效源代替;間接邊界元法是從疊加原理出發(fā)，在邊界上設(shè)置一組單層源或雙層源，用該組源疊加后在場點產(chǎn)生的磁效應(yīng)來等效替代原始場。間接邊界元法與直接邊界元法相類似，前者可從后者推導(dǎo)而得，唯一的區(qū)別在于前者同一點上只有1種等效源，而后者有2種等效源。

分域媒質(zhì)上的邊界積分法與積分方程法有許多類似之處。二者的區(qū)別在于邊界積分法采用了分部積分，從而把描述磁場的第二類Fredholm積分方程在一定條件下轉(zhuǎn)化為邊界積分方程。除與積分方程法同樣的優(yōu)點之外，其離散只需要在邊界進(jìn)行，因此可降低問題的維數(shù)。

國內(nèi)外文獻(xiàn)表明，用邊界積分法進(jìn)行艦船磁場數(shù)值建模已展現(xiàn)出其方便性，因此也得到了廣泛的應(yīng)用，如文獻(xiàn) ［17－19］采用邊界積分法對艦船殼體感應(yīng)磁場進(jìn)行了數(shù)值建模，計算結(jié)果與實驗結(jié)果表現(xiàn)出很好的一致性。

雖然邊界積分法在艦船磁場數(shù)值建模中具有離散簡單、可降低維數(shù)等優(yōu)點，但其形成的系數(shù)矩陣也為不對稱滿陣，當(dāng)離散單元較多時，需花費(fèi)較多內(nèi)存和求解時間。

由磁場積分法的上述過程可以看出，磁場積分法是一種適用的艦船磁場數(shù)值建模方法。從目前最新研究進(jìn)展可以發(fā)現(xiàn)，快速多級子技術(shù)［20］、小波稀疏化矩陣技術(shù)［6］等可用來解決磁場積分法中存在的不足之處，因此，磁場積分法將在艦船磁場數(shù)值建模中具有廣闊的應(yīng)用前景。

3.4 磁場微分法

在磁場微分法中最典型的數(shù)值計算方法為有限元法。有限元法是以變分原理為基礎(chǔ)的一種數(shù)值計算方法，它應(yīng)用變分原理把所要求解的邊值問題轉(zhuǎn)化為相應(yīng)的變分問題，對求解區(qū)域進(jìn)行離散、插值，離散化后的變分問題為普通多元函數(shù)的極值問題，進(jìn)而得到一組多元的代數(shù)方程組，求解代數(shù)方程組即可得到邊值問題的數(shù)值解。一般需經(jīng)過如下步驟［3］:

1)給出與待求邊值問題相應(yīng)的泛函及其變分問題;

2)將連續(xù)域離散成剖分單元之和，選擇合適的插值函數(shù);將未知的連續(xù)函數(shù)離散成有限項函數(shù)之和，即將無限個自由度的問題離散成有限個自由度問題;

3)求解泛函極值，離散出矩陣方程;

4)用直接法或迭代法等方法求解有限元方程。

由上述有限元法求解邊值問題步驟可以看出，找出與邊值問題相對應(yīng)的變分問題是有限元法的關(guān)鍵。利用有限元法求解邊值問題過程中，形成的系數(shù)矩陣是稀疏陣，隨著求解大型稀疏方程組數(shù)值技術(shù)的發(fā)展，有限元法在眾多領(lǐng)域得到了較多的應(yīng)用。文獻(xiàn) ［21－24］表明，由于艦船磁場數(shù)值建模屬于開域靜磁場建模問題，直接應(yīng)用有限元法進(jìn)行求解需對區(qū)域邊界進(jìn)行特殊處理。

目前國外一些商業(yè)有限元軟件也可以用來計算艦船感應(yīng)磁場，比如OPERA有限元軟件與TOSCA分析軟件被多個國家的海軍用于艦船感應(yīng)磁場建模;另外，由CEDRAT開發(fā)的FLUX3D有限元軟件也已經(jīng)被用于研究艦船磁特征，該軟件可以采用薄殼單元對薄鋼板外殼的船體結(jié)構(gòu)進(jìn)行感應(yīng)磁場建模，并且可以計算消磁線圈的繞組效率。但是這些商用軟件的價格十分昂貴，功能也較為有限，而且難以實現(xiàn)利用測量值對計算結(jié)果進(jìn)行修正的功能，這在一定程度上限制了這些軟件的推廣。

有限元法是用有限個自由度來近似描述一個連續(xù)體。在開域問題中，研究區(qū)域為整個空間，則必須把邊值為0的邊界取到相當(dāng)遠(yuǎn)處，這會使網(wǎng)格和節(jié)點數(shù)變得非常多 (據(jù)稱，該軟件對艦船的剖分單元數(shù)達(dá)到100萬以上)，這使計算機(jī)內(nèi)存和CPU時間消耗非常大，如此巨大的數(shù)值求解系統(tǒng)可能難以用于艦船磁場測量的數(shù)據(jù)處理工作。所以在要研究的開域問題中單獨(dú)應(yīng)用有限元法需要進(jìn)行一些等效邊界條件的處理。通過保角變換將開域問題映射為閉域問題可方便利用有限元法來進(jìn)行艦船磁場數(shù)值建模研究。此外，逆矢徑變換、有限元法與其他方法混合建模等可彌補(bǔ)有限元法的不足。

4 艦船固定磁場數(shù)值計算的理論難點

艦船固定磁場的求取一般只能用反演的方法，可以歸屬到磁場逆問題的計算范疇，磁場逆問題在理論上有以下兩大難點:

一是解的非唯一性，即有無數(shù)艦船區(qū)域內(nèi)磁源分布的解對應(yīng)于給定的磁場分布，艦船區(qū)域內(nèi)的磁源分布不能由船外的磁場唯一確定，從而這個問題在本質(zhì)上是不適定的。在磁場逆問題中通常很難獲得唯一解，必須加入很多限制條件，才有可能得到磁源參數(shù)的唯一解。

二是解的不定性或解不連續(xù)依賴數(shù)據(jù) (又稱解不穩(wěn)定、病態(tài)特性)，即使輸入極小的噪聲或擾動都會引起解的振蕩，這種病態(tài)趨勢隨解的參量數(shù)目增加而上升，即表征磁源磁化特征的待定變量越多，得到的矩陣方程的病態(tài)越嚴(yán)重，它使許多常規(guī)求解方法變得不穩(wěn)定從而失效，會導(dǎo)致磁源參數(shù)的解非常不穩(wěn)定，而且這種不穩(wěn)定與算法無關(guān)，僅提高數(shù)據(jù)準(zhǔn)確性是難以克服這些難點的，只有通過增加信息來修正不穩(wěn)定，即在求解磁場逆問題時對源區(qū)的結(jié)構(gòu)及解的空間參量有某些限制與約束，才能得到在一定范圍內(nèi)合理的唯一解。

在艦船固定磁場的逆問題求取中，面臨的一個重要問題是:測量得到的固定磁場垂向分量總和數(shù)據(jù)很難提供足夠的信息來確定一個唯一、穩(wěn)定的解，想要從理論上精確地求出其剩余磁化強(qiáng)度是非常困難的，只有通過研究固定磁場的變化特征及船體的磁化特點，根據(jù)經(jīng)驗或試驗結(jié)果加入一些限制條件 (包括測量信息和先驗信息)才可能分離出固定磁場分量。

4 結(jié)語

電磁場數(shù)值計算技術(shù)的發(fā)展，已為艦船磁場數(shù)值計算提供了技術(shù)基礎(chǔ)。鑒于艦船磁場數(shù)值計算其獨(dú)特的特點，選擇合適的靜磁建模方法來分析研究艦船磁場顯得十分必要。根據(jù)上述分析可知，各種不同的艦船磁場數(shù)值計算方法都有各自優(yōu)缺點，采用單一的數(shù)值計算方法難以徹底解決艦船磁場數(shù)值建模問題。艦船磁場數(shù)值計算應(yīng)用背景的多樣化，要求必須根據(jù)艦船磁隱身技術(shù)中特定的應(yīng)用背景合理選擇滿足工程精度要求的艦船磁場數(shù)值計算方法。

［1］葉平賢，龔沈光.艦船物理場［M］.北京:兵器工業(yè)出版社，1992.YE Ping-xian，GONG Shen-guang.Ship physics field［M］.Beijing:Ordnance Industry Press，1992.

［2］周耀忠，張國友.艦船磁場分析計算［M］.北京:國防工業(yè)出版社，2004.ZHOU Yao-zhong，ZHANG Guo-you.Ship magnetic field analysis and calculation［M］.Beijing:National Defence Industry Press，2004.

［3］盛劍霓.工程電磁場數(shù)值分析［M］.西安:西安交通大學(xué)出版社，1991.SHENG Jian-ni.Numerical analysis of engineering electromagnetic fields［M］.Xi'an:Xi'an Jiaotong University Press，1991.

［4］L Shen-dao，X Chang-han，Z Guo-hua.Magnetic model of magnetic objects based on genetic algorithm ［C］.ICEF，2008.

［5］KAMONDETDACHA R，KILDISHEV A V，NYENHUIS J A.Multipole characterization of a magnetic source using a truncated SVD［J］.IEEE Trans.Magn.2004，40(4):2176－2178.

［6］閻照文.電磁場新型等效源法［M］.北京:北京工業(yè)大學(xué)出版社，2003.YAN Zhao-wen.Electromagnetic equivalent source method［M］.Beijing:Beijing Industrial University Press，2003.

［7］鄭勤紅，盛劍霓.三維靜電場分析的多級理論［J］.華北電力大學(xué)學(xué)報，1996，10(4):14 －19.ZHENG Qin-hong，SHENG Jian-ni.Multi-level theory of three-dimensional electrostatic field analysis［J］.North China Electric Power University，1996，10(4):14 －19.

［8］鄭勤紅，曾華，解福瑤.用多極理論計算三維靜磁場問題［J］.中國電機(jī)工程學(xué)報，2005，10(25):348 －351.ZHENG Qin-hong，ZENG Hua，XIE Fu-yao.Multipole theory for three-dimensional static magnetic field［J］.Proceedings of the CSEE，2005，10(25):348 －351.

［9］鄭勤紅，盛劍霓，解福瑤，等.用二階矢量位計算三維矢量磁位靜磁場問題的多極理論［J］.中國電機(jī)工程學(xué)報，1998，18(4):279 －382.ZHENG Qin-hong，SHENG Jian-ni，XIE Fu-yao，et al.Multipole theory of using second order vector portential to solve three-dimensional static magnetic field［J］.Proceedings of the CSEE，1998，18(4):279 －382.

［10］KILDISHEV A V，NYENHUIS J A，MORGAN M A.Multipole analysis of an elongated magnetic source by a cylindrical sensor array.IEEE Trans.Magn，2002，38(5):2465－2467.

［11］KILDISHEV A V，KAMONDETDACHA R，NYENHUIS J A.Prediction of the magnetic field below an axissymetrical planar sensor array with a magnetic source located above the measurement plane［J］.DT09，IEEE，2002.

［12］KILDISHEV A V，NYENHUIS J A.External magnetic characterization of marine vehicles［M］.IEEE，2000:1145－1147.

［13］樊明武，顏威利.電磁場積分方程法［M］.北京:機(jī)械工業(yè)出版社，1988.FAN Ming-wu，YAN Wei-li.Electromagnetic field integral equation method［M］.Beijing:MechanicalIndustry Press，1988.

［14］CHADEBEC O，COULOMB J L，BONGIRAUD J P，et al.Recent improvements for solving inverse magnetostatic problem applied to thin shells［J］.IEEE Trans，Magn.，2002，38(2):1005 －1008.

［15］郭成豹，何明，周耀忠.積分方程法計算艦船感應(yīng)磁場［J］.海軍工程大學(xué)學(xué)報，2001(12):72 －75.GUO Cheng-bao，HE Ming，ZHOU Yao-zhong.Integral equation method for calculating ship induction magnetic field［J］.Naval University of Engineering，2001(12):72 －75.

［16］周國華，肖昌漢，閆輝，等.基于單元表面積分和微粒群優(yōu)化算法的鐵磁物體感應(yīng)磁場求解方法［C］.全國電工理論與新技術(shù)2007學(xué)術(shù)會議論文集，2007:322－326.ZHOU Guo-hua，XIAO Chang-han，YAN Hui，et al.A method to calculate the induced magnetic field of ferromagnetic objects in a weak magnetic field［C］.National Electrical Theory and New Technology 2007 Conference Papers，2007:322 －326.

［17］KAHLER G R，DELLA TORRE E.Minimizing the deformation of a static magnetic field by the presence of a ferromagnetic body［J］.IEEE Trans.Magn.1991，27(6):5025 －5027.

［18］RIOUX-DAMIDAU F，BANDELIER B，PENVEN P.Perturbation of a static magnetic field by a thin shell ferromagnetic device［J］.IEEE Trans.Magn.1993，29(6):2434 －2436.

［19］RIOUX-DAMIDAU F，BANDEFIER B，PENVEN P.A fast and precise determination of the static magnetic field in the presence of thin iron shells［J］.IEEE Trans.Magn.，1995，31(6):3491－3493.

［20］盛新慶.計算電磁學(xué)要論［M］.北京:科學(xué)出版社，2004.SHENG Xin-qing.Computing generality for electromagnetism［M］.Beijing:Science Press，2004.

［21］RODGER D，LEONARD P J，LAI H C.Surface elements for modelling 3D fields around thin iron sheets［J］.IEEE Trans.Magn.，1993，29(2):1483 －1486.

［22］BRUNOTTE X，MEUNIER G.Line element for efficient computation of the magnetic field created by thin iron plates［J］.IEEE Trans.Magn.，1990，26:2196 －2198.

［23］BRUNOTTE X，MEUNIER G，Jean-Paul Bongiraud.Ship magnetizations modelling by the finite element method［J］.IEEE Trans.Magn.，1993，29(2):1970 －1975.

［24］DORZE F L，BONGIRAUD J P.Modeling of degaussing coils effects in ships by the method of reduced scalar potential jump［J］.IEEE Trans.Magn.，1998，34(5):2477 －2480.