基于FLUX軟件的艦船磁場計算

姜智鵬 莊飚 關(guān)濤 母海方

（92537部隊，北京 100161）

1 引言

艦載消磁系統(tǒng)是有效降低艦船磁場的重要措施之一，而要設(shè)計理想的消磁系統(tǒng)，必須首先對艦船磁場進行預(yù)測。目前常用的幾種艦船磁場預(yù)測方法分別是母型艦換算、經(jīng)驗公式計算、磁性船模試驗以及數(shù)值仿真計算。其中，母型艦法的使用條件比較苛刻，實際應(yīng)用時一般需要進行大量的修正，在修正過程中需要和其它方法結(jié)合使用，如經(jīng)驗公式、物理模擬試驗、仿真計算等；經(jīng)驗公式法主要依賴于計算人員的經(jīng)驗，受計算人員的人為干擾因素較大；磁性船模試驗法是目前國內(nèi)使用最廣泛的方法，該方法主要是對船模制造提出了很高的要求，此外其建造周期長、通用性差、試驗成本高。各種數(shù)值計算方法與計算機軟件技術(shù)相結(jié)合產(chǎn)生了數(shù)值仿真計算法，其具有快速、簡便、經(jīng)濟的優(yōu)點。隨著研究的不斷深入，數(shù)值仿真計算法在艦船磁場計算的應(yīng)用中越來越廣泛，特別是各種商用電磁場計算軟件的成熟，為艦船磁場的計算提供了基本的仿真平臺。本文以電磁場有限元軟件FLUX為基本的仿真平臺進行艦船磁場計算研究。

2 艦船磁場計算的求解方程

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

麥克斯韋電磁方程是電磁場分析的基礎(chǔ)，而靜態(tài)磁場方程由麥克斯韋電磁方程簡化而來，其基本方程為式(1)、(2)和(3)[1]。

磁介質(zhì)的高斯定律（磁通守恒定律）

麥克斯韋－安培方程（安培環(huán)路定律）

磁特性方程

靜態(tài)磁場問題根據(jù)研究對象可分為 2D模型和 3D模型，根據(jù)求解方程的狀態(tài)變量又可以分為矢量模型和標(biāo)量模型，對于 2D靜磁問題一般使用矢量磁位，而 3D靜磁問題則一般使用標(biāo)量磁位，因為標(biāo)量磁位在計算精度、所需內(nèi)存和解的唯一性上比矢量磁位有優(yōu)勢，這對于 3D問題很關(guān)鍵。艦船磁場的計算屬于 3D靜磁問題，因此使用標(biāo)量磁位，相比于矢量磁位只需要1/3的未知數(shù)，可大大減少計算時間及對計算機內(nèi)存的需要。

在靜態(tài)磁場問題的求解中引入標(biāo)量磁位φ，則

有多種選擇且每種選擇各對應(yīng)一種標(biāo)量磁位，每種標(biāo)量磁位又對應(yīng)不同的磁場求解方程。

有限元法靜磁場標(biāo)量模型的通用求解方程如下：

如下三種類型分別對應(yīng)不同的三種標(biāo)量磁位公式：沒有電流源、電流源不需要剖分、電流源需要剖分[1]。

要求解的方程為：

狀態(tài)變量為全標(biāo)量磁位φtot。

要求解的方程為：

要求解的方程為：

部分標(biāo)量磁位的使用，一方面提高了艦船感應(yīng)磁場的計算精度，另一方面使得無限域的模擬成為可能[2]。

3 艦船磁場計算

FLUX軟件被國外海軍廣泛應(yīng)用于艦船磁場的計算及消磁繞組的優(yōu)化設(shè)計，其基于有限元法[2-5]。針對艦船復(fù)雜的形體、船體大尺度與船殼薄壁之間的不一致、開域磁場的計算，該軟件都能較好的處理。

以某船為研究對象，使用FLUX軟件進行了艦船磁場的計算研究，現(xiàn)將基于FLUX軟件的艦船磁場計算的一般步驟和注意問題總結(jié)如下：

① 艦船實體幾何模型建模。利用第三方軟件如SolidWorks等對艦船進行三維實體建模，模型建好后保存為 FLUX軟件可以識別的格式，如STEP格式。需要注意的是在艦船實體模型的建模過程中，對艦船的形體需要進行適當(dāng)?shù)暮喕雎砸恍┘?xì)微的地方，如很多倒角、圓角以及一些較小的部件等，因為過小的部件和細(xì)節(jié)在磁場分析中會導(dǎo)致剖分單元的困難，造成節(jié)點和單元數(shù)過多，甚至于無法進行剖分。

② 艦船幾何模型導(dǎo)入與糾錯。在SolidWorks中繪制完成艦船實體幾何模型后，就可以將生成的STEP文件導(dǎo)入FLUX軟件。在模型導(dǎo)入過程中需要注意的問題是導(dǎo)入坐標(biāo)系的選擇，此外，還會碰到模型無法導(dǎo)入的問題。由于生成的STEP文件含有的艦船幾何模型信息很多，受軟件所支持內(nèi)存量的限制，可能會造成導(dǎo)入操作無法成功。針對該問題，一勞永逸的辦法是增加軟件所能支持的內(nèi)存量，但在現(xiàn)有條件不能滿足的情況下，可考慮利用FLUX軟件的能夠支持幾何模型多次導(dǎo)入的功能，將原全船幾何模型分割成多個部分，分別保存為STEP格式文件。由于每個文件包含的幾何信息可以大為減少，能夠解決艦船實體幾何模型無法導(dǎo)入的問題。模型導(dǎo)入后需要對導(dǎo)入模型進行幾何檢查，查看有無幾何缺陷。如果有，需要對幾何模型進行修正和簡化。導(dǎo)入后的艦船實體幾何模型由于形體復(fù)雜，幾何缺陷會比較多，包括致命缺陷和非致命缺陷。雖然FLUX軟件具有自動修復(fù)功能，但自動修復(fù)功能無法滿足大多數(shù)幾何缺陷的修復(fù)要求，為此需要手動修復(fù)。

③ 模型的網(wǎng)格剖分。艦船形體特殊，結(jié)構(gòu)復(fù)雜，因此艦船三維實體模型網(wǎng)格剖分的基本考慮是盡量使用只有四面體單元的自動剖分器，這樣沒有六面體、棱柱等其他形式的單元混合，剖分會相對簡單一些。FLUX軟件是通過剖分單元尺寸和數(shù)量對網(wǎng)格剖分進行控制和調(diào)整的。其中，通過剖分點的設(shè)置對單元尺寸進行控制，即點調(diào)節(jié)；通過剖分線的設(shè)置對單元數(shù)量進行控制，即線調(diào)節(jié)。對于模型的絕大部分，使用點調(diào)節(jié)就可以，而對于一些細(xì)節(jié)的部分則建議進行線調(diào)節(jié)，這樣更利于剖分的成功[6]。歸納起來，成功實現(xiàn)滿足計算精度要求的剖分的基本原則是保證邊界盒有2～3層剖分單元（內(nèi)、外邊界盒之間），關(guān)鍵區(qū)域的剖分單元尺寸與磁性船模磁場實測點之間的最小距離相接近，同時盡量保證剖分的面元為等邊三角形。

④ 物理屬性的設(shè)置。對于屬于靜磁問題的艦船磁場計算問題，其在FLUX中使用靜態(tài)磁場模型、標(biāo)量磁位求解方程。本文中研究對象磁性船模由三種不同磁特性的鋼材制造，因此在其艦船磁場模型中有四種不同磁導(dǎo)率的材料（加上空氣），磁化曲線為線性，以相對磁導(dǎo)率作為唯一磁特性參數(shù)對材料進行定義。在模型的材料設(shè)置中，分別使用了volume regions和surface regions，并關(guān)聯(lián)到相關(guān)的體和面。材料屬性設(shè)置完成后，進行外加磁場和線圈的設(shè)置。外加磁場為空間分布均勻的地磁場，為了計算結(jié)果分析的方便，地磁場分量是分開設(shè)置的，即一次求解過程，只設(shè)一個地磁場分量（垂直分量或水平分量）。

⑤ 求解及后處理。FLUX軟件可進行參數(shù)化求解，在艦船磁場的求解中，可將各材料的相對磁導(dǎo)率設(shè)為變量，對每個變量賦多個值，如對三種材料的相對磁導(dǎo)率各賦三個值，則通過一次求解過程就可以得出 27種不同磁導(dǎo)率組合的計算結(jié)果。求解完成后，可直接在軟件中對某一物理量的值進行三維顯示，還可定義平面矩陣，將該平面矩陣上特定物理量的值輸出到文本。

4 計算結(jié)果的驗證和分析

對某地地磁環(huán)境下A（淺深度）、B（較深深度）兩個不同深度的艦船磁場的Zix、Ziy、Ziz分量進行了仿真計算，并與磁性船模磁場的實測結(jié)果進行比較。以Ziz分量為例，其A、B兩個深度的磁場比較圖見圖1和圖2。

從圖1、2中可以看出，計算曲線與測量曲線的趨勢基本一致。經(jīng)統(tǒng)計,所有點計算值和實測值的相對誤差均小于20%，其中絕大部分點的誤差在5%以內(nèi)，小部分點的誤差在5%～15%之間，只有個別點的誤差超過15%，在15%～20%之間。產(chǎn)生誤差的原因有很多，經(jīng)過分析，認(rèn)為主要原因有以下幾種：

① 幾何模型的簡化所產(chǎn)生的誤差。由于艦船形體復(fù)雜，在幾何建模時，對其流線型的船體以及軸等附件進行了一定程度的簡化，這勢必引入計算誤差。

② 艙室設(shè)置與實際分布情況不盡相同所產(chǎn)生的誤差。由于對磁性船模艙室分布情況沒有完全掌握，艙壁是根據(jù)對相關(guān)資料的解讀和判斷而設(shè)置的，不能避免與實際分布情況有出入。

圖 1 垂向感應(yīng)磁場的垂向分量 Ziz（A深度）

圖 2 垂向感應(yīng)磁場的垂向分量 Ziz（B深度）

③ 有限元法數(shù)值計算的固有計算誤差。

④ 磁導(dǎo)率測試的誤差。磁導(dǎo)率是根據(jù)原材料樣品的磁特性測試確定的，受樣品加工方式的影響較大，樣品切割方式的不同會導(dǎo)致磁導(dǎo)率變化的不同。

⑤ 測磁系統(tǒng)固有的測量誤差。

⑥ 測量點定位引起的測量誤差。若某一測量點定位不夠精確，則不可避免地引起一定的測量誤差，將對應(yīng)點的計算值與它相比較，就會出現(xiàn)不一致的情況。

上述原因引起的誤差可通過采取措施進行控制。

5 結(jié)束語

本文以專用的電磁場有限元分析軟件 FLUX為仿真計算的平臺，對艦船磁場的有限元仿真計算問題進行了研究和探索。通過對艦船實體幾何模型的建模、模型的有限元剖分等進行的深入研究，實現(xiàn)了艦船磁場的有限元計算。仿真計算結(jié)果經(jīng)磁性船模實測磁場值驗證，誤差小于20%，可應(yīng)用于工程實踐。

[1] Phsical Applications: Magnetic, Electic, Thermal. [M].FLUX 9.30 User’s guide. 2005, 3:15-40.

[2] Brunotte X, Meunier G, Bongiraud J P. Ship Magnetizations Modelling by the Finite Element Method [J]. IEEE Trans. Mag, 1993, 29(2):1970-1975.

[3] Dr. Santana-Diaz E, Tims R. A Complete Underwater Electric and Magnetic Signature Scenario Using Computational Modeling [R].

[4] Pierre Dallet J. Ship Magnetic Modelization by Finite Element[J]. MARELEC, 2001, 35:5.

[5] Dr. Birsan M. User Experience: Defence Research for Development Canada[J]. FLUX MAGAZINE, 2006,50:11.

[6] General Tools Geometry and Mesh [M]. FLUX 9.30 User’s guide. 2005, 1:131-224.