祝小雨 劉勝道 周國華

(海軍工程大學電氣工程學院 武漢 430033)

?

一種磁性物體周邊磁場的估算方法*

祝小雨 劉勝道 周國華

(海軍工程大學電氣工程學院 武漢 430033)

在國防和航天工業(yè)中,有時需要對某設備的周邊磁場進行測量。而在實際操作中,單靠測量顯然無法較全面地反映出其周邊磁場的特點。為了推算出一般磁性物體的周邊磁場并提取其特征,從邊界元理論出發(fā),根據物體周邊合適點的磁場測量建立反演模型,求出正演所需要的面磁荷參數(shù)從而建立其正演模型,進而推算出其周邊磁場。用一般鐵磁性物體進行了驗證性實驗,其推算值與測量值吻合得比較好,誤差小,可應用于實際工程中。

磁場測量; 反演; 正演; 面磁荷

Class Number TM154.3

1 引言

在國防和航天工業(yè)中經常需要對某磁性設備或者大型船舶的周邊磁場進行測量,而在實際操作中,往往無法通過測量得到周邊所有的磁場數(shù)據,因此我們希望能建立磁性物體的數(shù)學模型并進行磁場推算。近年來,磁場建模方法得到了廣泛研究且取得了較多成果[1]。用于船舶磁場推算的數(shù)學模型主要有兩類:一類是依據嚴格的理論推導出來的,如大平面法、邊界元法等,其計算精度高,但需要整個場域邊界上的磁場數(shù)據,對測量要求較高[2];另一類是磁體模擬法,它基于在物體若干測量點處產生磁場的相似性,將整個物體磁場等效為若干分布于船舶所占空間內的模擬體的磁場[3],其只需少量的磁場測量數(shù)據,但缺乏理論上的嚴密性。面磁荷法是基于邊界元理論推導出的一種磁場建模方法[4],該方法結合了上述兩類數(shù)學模型的優(yōu)點,在經過嚴密的理論推導之后,截取部分場域邊界,使得該方法既能滿足工程上的計算精度要求,又能用少量的磁場測量數(shù)據建立模型。實際磁性設備的測量與計算表明,該方法所建立的模型具有計算精度高、穩(wěn)定性好及易于實現(xiàn)等優(yōu)點。

2 鐵磁性物質面磁荷建模方法

一般采用如圖1所示方式測量船舶的磁場,即點測法。所謂點測法,即在合適的場地將船舶固定在一定航向,然后將傳感器固定在水下某一深度呈一定的點狀排布,逐點測出其磁場數(shù)值[5]。

船舶磁場一般隨深度的增加而衰減,不同類型的船舶對應不同的標準深度。

當已知磁源下方某平面S上的等效面磁荷分布時,其下方任意點P處的磁感應強度可表示為[6～10]

(1)

其中rP為場點P矢徑,rQ為源點矢徑且位于面S上,σ(Q)為等效面磁荷密度。對于上式而言,σ(Q)未知,那么先根據測量出來的磁場數(shù)據對其進行反演,然后再對船舶磁場進行推算。式(1)中的S面理論上要求無限大,但是隨著距離增大時,磁場逐漸衰減到零。故可選擇適當?shù)腟T截面離散成若干面單元,當面單元足夠小的時候,其面磁荷密度則可近似看作常數(shù)。

圖1 船舶磁場測量示意圖

2.1 等效面磁荷磁場推算模型的推導

圖2 一般鐵磁物體磁場推算中等效面磁荷法計算原理示意圖

圖2所示為一般鐵磁物體磁場推算中等效面磁荷法計算原理示意圖,由邊界元理論,ABCD構成的區(qū)域閉合面S內任意一點P的標量磁位可表示為

(2)式中,|rPM|=|rP-rM|,rP為場點矢徑,rM為位于邊界S上的源點矢徑,n為外法線方向上的單位矢量。

將面AD、BC、CD延伸到無窮遠處,磁場隨著距離衰減在無窮遠處為0,故只需要在AB面上對式(2)進行積分?？紤]到在AB面上有n=-z,則有:

(3)

若場點選在點P關于AB面對稱的點Q處,由于點Q位于閉合區(qū)域外,則有:

(4)

式中|rQM|=|rQ-rM|。由于|rPM|=|rQM|且?(1/|rQM|)/?z=-?(1/|rPM|)/?z,那么式(4)可轉化為

(5)

式(3)加式(5)并對場點坐標求梯度的負數(shù)可得P點磁感應強度為

(6)

式中σM=-?φM(rM)/?z,可視為AB面上的等效面磁荷密度。也就是說,當已知AB面上的面磁荷密度分布時,由式(6)即可求解AB下空間任意點的磁感應強度。

2.2 等效面磁荷逆問題模型的建立及求解

我們知道,磁場隨著場點與磁源的距離增大而逐漸衰減,截取無窮大平面AB中有效區(qū)域ST,離散成若干足夠小的小面單元,此時,其面磁荷密度可視為常數(shù),由式(6)

(7)

其中i為離散單元數(shù),STi表示第i個單元面,σi表示第i個面單元等效面磁荷密度。

由圖2,從S內取合適的測量面測量其若干點磁感應強度值,若所選取的場點P為對應傳感器位置處,那么就能建立以單元面磁荷密度為未知量的等效面磁荷逆問題離散模型。

(8)

其中,i為測量點個數(shù),式(8)的矩陣形式為Ax=b。

3 實驗驗證

對某設備的六個面進行編號記為A、B、C、D、E、F。將設備置于某南北朝向軌道上。軌道上每5cm標記一個點,共81個點。將傳感器置于第41號點正下方。如圖3所示。

將設備置于較遠位置測出背景磁場,然后將設備放在軌道上由南至北81個點處,依次測出每點的磁場。

通過面磁荷法由上一層傳感器的數(shù)據推算下一層傳感器的磁場,然后將推算出的數(shù)據與下一層傳感器的實測數(shù)據進行比較。

推算數(shù)據與實測數(shù)據的比較如圖4～9所示。(圖中162個點分別為12、13號傳感器數(shù)據)。

圖3 設備周邊磁場測量示意圖

圖4 設備A面下方推算與實測數(shù)據比較圖

圖5 設備B面下方推算與實測數(shù)據比較圖

圖6 設備C面下方推算與實測數(shù)據比較圖

圖7 設備D面下方推算與實測數(shù)據比較圖

圖8 設備E面下方推算與實測數(shù)據比較圖

圖9 設備F面下方推算與實測數(shù)據比較圖

由圖可知,用面磁荷法得到的推算值與實際測量值吻合較好。

通過建立的模型對此設備周圍的磁場分布進行推算,選取距離C面300mm大小為570mm×570mm的面,研究其上均勻分布的31×31點陣列的磁場。

圖10 距C面300mm處面上Bx分布圖

圖11 距C面300mm處面上By分布圖

圖12 距C面300mm處面上Bz分布圖

圖13 距C面300mm處面上總場B分布圖

從圖10～13中可以看出,Bx的最大值為61nT,最小值為-31nT;By的最大值為61nT,最小值為-26nT;Bz的最大值為98nT,最小值為-48nT;總場B的最大值為101nT,最小值為25nT。

4 結語

針對實際操作中單靠測量無法較全面地反映磁性設備周邊磁場特點的問題,從邊界元理論出發(fā),在嚴密的理論推導之后,通過截取部分場域邊界,得到了一種磁性物體周邊磁場的估算方法。該方法既能滿足工程上的計算精度要求,又能用少量的磁場測量數(shù)據建立模型,從而可推算出其周邊磁場。實際磁性設備的測量與計算表明,該方法不但所需的磁場測量數(shù)據較少,而且能滿足工程上的計算精度要求,穩(wěn)定性高,易于實現(xiàn)。

[1] 周耀忠,張國友.艦船磁場分析計算[M].北京:國防工業(yè)出版社,2004:135-136.

[2] 周國華,肖昌漢,劉勝道,等.一種艦艇磁隱身中磁場推算的新方法[J].兵工學報,2009,30(7):951-957.

[3] 郭成豹,張曉鋒,肖昌漢,等.地磁場中空心鐵磁圓柱體縱向磁化狀態(tài)數(shù)值分析[J].哈爾濱工程大學學報,2006,27(6):886-890.

[4] 周國華,肖昌漢,閆輝,等.一種弱磁作用下鐵磁物體感應磁場的計算方法[J].哈爾濱工程大學學報,2009,30(1):91-95.

[5] 周國華,肖昌漢,劉勝道,等.船舶磁場的面磁荷建模方法及穩(wěn)定性分析[J].哈爾濱工程大學學報,2010,31(5):553-558.

[6] 陳杰,魯習文.基于磁荷面分布的艦船磁場預測方法[J].物理學報,2009,58(6):3839-3844.

[7] Norgren M, Sailing H 2000 IEEE Trans. Magn. 34 370.

[8] Tejedor M, Rubio H, Elbaile L, lglesias R 1995 IEEE Trans. Magn. 25 800.

[9] Frumkis L, Benzion K 1999 IEEE Trans. Magn. 30 4155.

[10] Wingo R A,Holmes J J, Lackey M H 1992 Nav. Eng. J.100 230.

Estimation Method of Magnetic Field Around an Magnetic Object

ZHU Xiaoyu LIU Shengdao ZHOU Guohua

(College of Electrical Engineering, Naval University of Engineering, Wuhan 430033)

In the national defence industy and aerospace industries, sometimes a device around the magnetic field needs to be surveyed and drawed. Actually it’s too hard to describe it only by measurements. In order to calculate and describe the magnetic field around the general magnetic objects, the suitable point around the object of the magnetic field is measured to derive the inverse model from the boundary theory, then the parameters are obtained to derive the direct model. After that the magnetic field around the object can be described. Confirmatory tests are completed, the projected values matches well with the measured values. The method can be applied in engineering.

measurements, inverse model, direct model, surface magnetic charge

2014年10月2日,

2014年11月27日

祝小雨,男,研究生,研究方向:磁場計算。

TM154.3

10.3969/j.issn1672-9730.2015.04.043