姚振寧，劉大明，朱興樂

(海軍工程大學 電氣與信息工程學院，湖北 武漢，430033)

在艦船磁性防護領域，通過測量艦船周圍空間磁場的磁感應強度，檢測艦船磁性狀態(tài)以及評估艦船磁性防護水平。對于鋼鐵建造成的艦船，若周圍沒有其他鐵磁物體，則其空間磁場只由艦船的固定磁場和被地磁場磁化產(chǎn)生的感應磁場組成。若艦船周圍存在如拖船、鐵磁設備等鐵磁物體且距離較近，則艦船與這些鐵磁物體相互磁化，產(chǎn)生附加的感應磁場，因而艦船周圍空間磁場的磁感應強度測量值并不等于艦船磁場與鐵磁物體磁場的磁感應強度的疊加。由于附加感應磁場的存在，空間磁場的磁感應強度測量值不能如實反映艦船磁性，如何由空間磁場的磁感應強度測量值計算艦船磁場的磁感應強度成為艦船磁防護工作中可能碰到的技術難題，解決這個問題必須求出附加感應磁場的磁感應強度。目前，許多文獻通過有限元法或積分方程法對磁性物體的感應磁場的磁感應強度計算進行了大量研究[1-7]，這2 種方法需要對物體或場域進行剖分，剖分一般較復雜，且剖分的疏密程度直接影響計算精度。另外，這些研究大多是計算地磁場磁化產(chǎn)生的感應磁場的磁感應強度，由物體間相互磁化產(chǎn)生的附加感應磁場的磁感應強度計算研究卻很少。為此，本文作者考慮附加感應磁場的影響，建立磁模型，再把物體的磁化等效為其磁模型中各磁偶極子的磁化，給出一種求解附加感應磁場的磁感應強度的計算方法。

1 基于磁體模擬法的磁場數(shù)學模型

艦船與鐵磁物體形狀不規(guī)則、磁化不均勻，難以用解析法求解其磁場的磁感應強度。一般采用磁體模擬法建立磁模型，即用若干個具有特定磁矩的磁性模擬體(磁偶極子、旋轉(zhuǎn)橢球體)所產(chǎn)生的磁場模擬不規(guī)則物體的磁場，然后由磁性模擬體計算場域空間的磁場的磁感應強度[8]。這里采用磁偶極子陣列進行建模。

假設磁偶極子數(shù)為n，艦船在第j 個測量點產(chǎn)生磁場三分量為hxj，hyj和hzj，則有方程組：

式中：mxi，myi和mzi分別為第i 個磁偶極子沿x，y 和z 軸方向的磁矩分量；fxji，fyji，fzji，gxji，gyji，gzji，exji，eyji和ezji為與測量點坐標、磁偶極子位置有關的常數(shù)，由它們組成了方程組的系數(shù)矩陣。在實際艦船磁場的磁感應強度測量中，在艦船下方某個平面測量磁場的磁感應強度，根據(jù)式(1)組成矛盾方程組，求其最小二乘解可得各個磁偶極子的磁矩，進而可求場域空間磁場的磁感應強度，即建立了磁模型。

磁偶極子原則上可以布設在所模擬物體區(qū)域的任何位置，為了簡便，一般布設在水線面上。為了提高模型精度，一般將磁偶極子個數(shù)取得較多，但太多的磁偶極子使系數(shù)矩陣列相關性增強，其條件數(shù)變得非常大，以上方程組將呈現(xiàn)病態(tài)，因而模型不再穩(wěn)定。系數(shù)矩陣條件數(shù)與磁偶極子位置有關，為了建立穩(wěn)定的磁模型，可以利用智能優(yōu)化算法來確定磁偶極子的分布位置[9-11]。

本文以系數(shù)矩陣條件數(shù)作為目標函數(shù)，利用改進的蟻群算法-蟻群系統(tǒng)[12-14]對磁偶極子位置進行優(yōu)化，如圖1 所示，考慮實際艦船與吃水線相交情況，把艦船水線面形狀近似看作橢圓，假設其磁偶極子個數(shù)為M, 將水線面沿長度方向也分成M 等分，則第i個磁偶極子在陰影范圍內(nèi)優(yōu)化取值(i=1, 2, …, M)。對于艦船周圍鐵磁物體，其磁偶極子位置的優(yōu)化按照此方法進行。

圖1 磁偶極子位置的優(yōu)化Fig.1 Magnetic dipoles’ optimization

2 周圍有鐵磁物體存在時的艦船磁場的磁感應強度計算

2.1 基本思想

當艦船周圍有鐵磁物體存在時，它們之間相互磁化產(chǎn)生各自的附加感應磁場。在實際的艦船磁場的磁感應強度測量中，如果只能測得其下方某一測量面上總磁場的磁感應強度h 與鐵磁物體磁場的磁感應強度h2，假設艦船磁場磁感應強度為h1，則有

式中：h1′為艦船的附加感應磁場的磁感應強度； h2′ 為鐵磁物體的附加感應磁場的磁感應強度。可見，由于艦船與鐵磁物體相互磁化產(chǎn)生了附加感應磁場，總磁場的磁感應強度與鐵磁物體磁場的磁感應強度之差并不等于真實的艦船磁場的磁感應強度。

由式(2)可知，要得到真實的艦船磁場的磁感應強度必須計算出附加感應磁場的磁感應強度。一個物體的附加感應磁場是由另一物體對其磁化產(chǎn)生的，這種磁化是非均勻磁場對非規(guī)則物體的非均勻磁化，無法直接用解析式進行磁場的磁感應強度計算，所以必須等效為對規(guī)則物體的均勻磁化。結(jié)合磁模型，不妨認為對物體的磁化等效為對其磁模型中各磁偶極子的磁化，這樣就可以直接用解析法計算磁場的磁感應強度。

2.2 附加感應磁場的磁感應強度的計算

設艦船磁模型中磁偶極子的三分量磁矩為mxj，myj，mzj(用M 個磁偶極子模擬，j=1, 2, …, M)，設鐵磁物體磁模型中磁偶極子的三分量磁矩為mxj，myj，mzj(用N-M 個磁偶極子模擬，j=M+1, M+2, …, N)，則鐵磁物體對艦船的第i 個磁偶極子產(chǎn)生的附加磁化場三分量為(i=1, 2, …, M)

艦船對鐵磁物體的第i 個磁偶極子產(chǎn)生的附加磁化場三分量為(i=M+1, M+2, …, N)

艦船、鐵磁物體的第i 個磁偶極子的附加感應磁矩三分量為(i=1, 2, …, N)

式中：λi為第i 個磁偶極子的附加感應磁矩與其附加磁化場之間所設的比例常數(shù)。式(4)中艦船磁偶極子的磁矩由h-h2根據(jù)式(1)近似求出，式(3)中鐵磁物體磁偶極子的磁矩由h2根據(jù)式(1)求出。

比例常數(shù)λi可借助于均勻地磁場對艦船、鐵磁物體的磁化而求出，即

式中：HD為地磁場的磁感應強度水平分量；Z1yk和Z2yk分別為艦船、鐵磁物體在第k 個測量點由HD橫向磁化產(chǎn)生的感應磁場的磁感應強度垂直分量，一般通過測量計算求出，作為已知數(shù)據(jù)。將Z1yk和Z2yk代入式(6)即可組成關于λi的矛盾方程組，求其最小二乘解可得λi。

根據(jù)式(3)～(6)可得附加感應磁矩 m′xi，m ′yi和 m′zi，進而可得附加感應磁場的磁感應強度為

3 計算實例

為了檢驗上述艦船磁場計算方法的可行性，利用磁場的磁感應強度測量值進行計算驗證。實驗磁場模型如圖2 所示。鐵質(zhì)的艦船模型長4.8 m、最大寬度0.54 m，其水線面近似為橢圓形，鐵磁物體為1 個長2.1 m、半徑0.2 m 的鐵質(zhì)空心圓筒，其水線面為長方形，坐標原點取船模水線面中心。在下方測量平面放置高精度的磁通門傳感器陣列(對應船模左舷外、左舷、龍骨、右舷、右舷外正下方5 條測量線，每條測量線31 個測量點)，傳感器縱向間距0.24 m、橫向間距0.27 m。首先，測量船模與圓筒在一起的總磁場的磁感應強度以及船模不存在時的圓筒磁場的磁感應強度，然后，根據(jù)上述算法計算船模磁場的磁感應強度；最后，測量船模磁場的磁感應強度，將其測量值與計算值進行比較。

圖2 實驗磁場模型示意圖Fig.2 Sketch map of experimental magnetic model

取磁偶極子所在水線面寬度為實際的1/4。用30個磁偶極子模擬船模磁場，用15 個磁偶極子模擬圓筒磁場，則磁偶極子通過蟻群算法優(yōu)化后的分布位置如圖3 所示。

圖3 磁偶極子分布位置Fig.3 Distribution of magnetic dipoles

在實際的艦船磁性檢測站，通常只分析研究艦船磁場的磁感應強度的垂直分量。圖4 所示為船模磁場的磁感應強度垂直分量在左舷、龍骨、右舷下方的測量值與2 種計算值的對比，其中直接計算值由總磁場的磁感應強度直接減去圓筒磁場的磁感應強度求出，且磁感應強度已經(jīng)過歸一化處理。由圖4 可見：算法計算值與測量值吻合較好，相對均方差為2.0%，最大相對誤差為5.4%，說明該計算方法能夠有效地計算有其他鐵磁物體干擾時的艦船磁場的磁感應強度；而不考慮附加感應磁場的磁感應強度的直接計算值在(1,3)區(qū)段誤差較大，相對均方差為5.2%，最大相對誤差為13.4%，這是因為船模與圓筒接觸區(qū)域的相互磁化比較嚴重，在其正下方(1,3)區(qū)段產(chǎn)生的附加感應磁場的磁感應強度較大。

圖4 船模磁場測量值與計算值對比Fig.4 Comparisons between calculated magnetic field and measured one

從上述計算過程來看，該計算方法也存在一定誤差：磁偶極子的附加磁化場是由下方測量面磁場的磁感應強度換算得到，且磁場的磁感應強度測量數(shù)據(jù)有限，這屬于遠場到近場的換算問題[15-16]，磁偶極子的附加磁化場存在誤差；另外，該算法把物體的磁化近似等效為其磁模型中各磁偶極子的磁化，沒有考慮磁偶極子之間相互干擾，導致磁偶極子的附加感應磁矩與其附加磁化場之間的比例常數(shù)λi存在誤差。但是，由計算結(jié)果可知：以上2 種誤差并不大，總的來說本文所提算法是有效可行的。

4 結(jié)論

(1) 提出了一種艦船與鐵磁物體在一起時的艦船磁場計算方法。

(2) 該算法首先分別建立艦船與鐵磁物體的磁模型，并利用蟻群算法優(yōu)化磁偶極子的分布位置，然后把物體的磁化等效為其磁模型中各磁偶極子的磁化，求出附加感應磁場的磁感應強度，進而得出真實的艦船磁場的磁感應強度。測試結(jié)果表明該計算方法可行。

[1] Bartusek K, Fiala P, Mikulka J. Numerical modeling of magnetic field deformation as related to susceptibility measured with an MR system[J]. Radioengineering, 2008, 17(4): 113-118.

[2] Nguyen T S, Guichon J M, Chadebec O, et al. Ships magnetic anomaly computation with integral equation and fast multipole method[J]. IEEE Transactions on Magnetics, 2011, 47(5):1414-1417.

[3] 郭成豹, 張曉鋒, 肖昌漢, 等. 地磁場中空心鐵磁圓柱體縱向磁化狀態(tài)數(shù)值分析[J]. 哈爾濱工程大學學報, 2006, 27(6):886-890.GUO Chengbao, ZHANG Xiaofeng, XIAO Changhan, et al. The numerical simulation of the longitudinal magnetic state of a hollow ferromagnetic cylinder in the geomagnetic field[J].Journal of Harbin Engineering University, 2006, 27(6): 886-890.

[4] 周國華, 肖昌漢, 閆輝, 等. 一種弱磁作用下鐵磁物體感應磁場的計算方法[J]. 哈爾濱工程大學學報, 2009, 30(1): 91-95.ZHOU Guohua, XIAO Changhan, YAN Hui, et al. A method to calculate the induced magnetic field of ferromagnetic objects in a weak magnetic field[J]. Journal of Harbin Engineering University, 2009, 30(1): 91-95.

[5] Le Dorze F, Bongiraud J P. Modeling of degaussing coils effects in ships by the method of reduced scalar potential jump[J]. IEEE Transactions on Magnetics, 1998, 34(5): 2477-2480.

[6] Pechenkov A N. A mathematical algorithm for solving the inverse problem of magnetostatic flaw detection[J]. Russian Journal of Nondestructive Testing, 2005, 41(11): 714-718.

[7] Pechenkov A N. Numerical simulation of an inverse problem of magnetostatic flaw detection[J]. Russian Journal of Nondestructive Testing, 2005, 41(11): 719-723.

[8] 周耀忠, 張國友. 艦船磁場分析計算[M]. 北京: 國防工業(yè)出版社, 2004: 135-136.ZHOU Yaozhong, ZHANG Guoyou. Analysis and computing of marine vessels’ magnetic field[M]. Beijing: National Defense Industry Press, 2004: 135-136.

[9] 劉勝道, 劉大明, 肖昌漢, 等. 基于遺傳算法的磁性目標磁模型[J]. 武漢理工大學學報, 2008, 32(6): 1017-1020.LIU Shengdao, LIU Daming, XIAO Changhan, et al. Magnetic model for magnetic objects based on genetic algorithm[J].Journal of Wuhan University of Technology, 2008, 32(6):1017-1020.

[10] 徐杰, 劉大明, 周國華, 等. 一種基于遺傳優(yōu)化算法的潛艇高空磁場換算方法[J]. 艦船科學技術, 2009, 31(1): 156-159.XU Jie, LIU Daming, ZHOU Guohua, et al. Application of genetic algorithm in continuation of magnetic field of submarines in high areas[J]. Ship Science and Technology, 2009,31(1): 156-159.

[11] 張朝陽, 肖昌漢, 徐杰. 基于微粒群優(yōu)化算法的艦船磁模型分析[J]. 華中科技大學學報, 2010, 38(11): 124-128.ZHANG Zhaoyang, XIAO Changhan, XU Jie. Analysis of magnetic model for ships using particle swarm optimization method[J]. Journal of Huazhong University of Science and Technology, 2010, 38(11): 124-128.

[12] Dorigo M, Gambardella L M. Ant colony system: A cooperative learning approach to the traveling salesman problem[J]. IEEE Trans on Evolutionary Computation, 1997, 1(1): 53-56.

[13] Dorigo M, Caro G D, Gambardella L M. Ant algorithms for discrete optimization[J]. Artificial Life, 1999, 5(2): 137-172.

[14] 付宇. 蟻群優(yōu)化算法的改進及應用[D]. 上海: 上海海事大學物流工程學院, 2006: 20-22.FU Yu. The improvement and application of ant colony optimization[D]. Shanghai: Shanghai Maritime University.College of Logistics Engineering, 2006: 20-22.

[15] 徐世浙. 地球物理中的邊界單元法[M]. 北京: 科學出版社,1995: 40-75.XU Shizhe. BEM in geophysics[M]. Beijing: Science Press,1995: 40-75.

[16] 閆輝, 肖昌漢, 殷克全, 等. 遞推算法在船舶磁場遠場向近場換算中的應用[J]. 兵工學報, 2010, 31(9): 1200-1203.YAN Hui, XIAO Changhan, YIN Kequan, et al. The application of recursive algorithm on ship’s magnetic field extrapolation[J].Acta Armamentarii, 2010, 31(9): 1200-1203.