厚金屬板與任意方位圓柱線圈間的電磁力計算

羅垚

(武漢大學 電氣工程學院，湖北 武漢 430072)

?

厚金屬板與任意方位圓柱線圈間的電磁力計算

羅垚

(武漢大學 電氣工程學院，湖北 武漢 430072)

摘要:針對一個交流電激勵的圓柱線圈傾斜放置于金屬板上方時，它們之間的電磁作用力無法以傳統方法進行解析計算的問題，引入高階矢量勢求解該問題，并采用矢量分析和Struve函數化簡，可得到適用于任意方位圓柱線圈的二重積分表達式。所得結果表明任意方位圓柱線圈至多受到兩個方向力的作用。在線圈軸垂直或平行于導電板的特殊情形，表達式可進一步簡化為單重積分，而線圈將只受垂直力的作用。隨后以數值算例對給出公式的有效性進行驗證。數值結果表明，對于線圈軸的某個傾角，垂直力分量周期均值可以為零，而水平力分量周期均值將在某個傾角處達到最大值。當線圈通以直流電時，僅當金屬板的磁導率與空氣磁導率一致時線圈才不受電磁力的作用。

關鍵詞:圓柱線圈；電磁力；高階矢量勢；Bessel函數；Struve函數

0引言

當一個置于平板導體附近的矩形截面圓柱線圈通以交變電流時，由于線圈的交變電磁場將在導體中感應出渦流，線圈本身將受到這一渦流的電磁力作用。這樣的電磁結構在電工技術中有許多應用，如快速動作斷路器(fast acting circuit breaker,FACB)[1]。對這種結構的電磁作用力的計算，常見于線圈軸線與導電體表面垂直的情形[1-3]。若線圈軸線傾斜于導電體表面，則此時的電磁力分析傳統上只能借助于純數值方法，如有限元法(FEM)，邊界元法(BEM)。然而，考慮到純數值方法無法給出待求電磁力與電磁系統各結構參數之間的函數關系，故在這一問題中解析或半解析工作仍然具有不可取代的優(yōu)勢。在已有工作中，為了便于計算通常將磁通密度B表為磁矢勢A的旋度，但這一方法卻無法用以處理線圈軸與導電面不垂直的一般情形。為了對一般方位的線圈電磁力進行解析處理，我們將引入高階矢量勢并借助Fourier積分求得一般情形的電磁力表達式。所得的電磁力表達式包括3個函數，它們均含二重積分，在線圈軸垂直于導電面的特殊情形可化為單重積分表達式。

1導電板附近任意載流環(huán)的電磁力作用

設有電導率κ，磁導率μ=μrμ0(μ0=4π×10-7H/m)的無窮大金屬板，其厚度較之趨膚深度而言可視為無窮大，即可將該金屬板視為一半無窮導電空間。在其上方的非導電空間磁導率為μ0，有一任意方位和形狀的載流環(huán)c置于其中，而環(huán)c上載有電流i=Icosωt。由于i隨時間變化，故其將在下方的導電板中感應出渦流，而環(huán)c本身則將受到這一渦流在上半空間激發(fā)出的磁場的電磁力作用。傳統的勢論僅能處理載流環(huán)為水平的情形而無法對任意方位的載流環(huán)進行解析計算。然而，若考慮Coulomb規(guī)范

divA=0,

(1)

則可對A可引入高階矢勢P[4]

A=rotP,

(2)

(3)

其中：P1，P2為兩個僅與坐標有關的標量函數，ex為x方向的單位矢量。由此可以Fourier積分求得導電板中的感應電流在上半空間中激發(fā)的磁矢勢

(4)

其中

(5)

(6)

式(4)由Biot-Savart定律得出，故對環(huán)c的方位與形狀均無限制。另由式(4)，式(6)的形狀可以看出，無論載流環(huán)的形狀與方位如何，由導電板中感應電流在上半空間激發(fā)的磁矢勢均不含垂直分量。

得到磁矢勢A后即由B=rotA得到感應磁通密度。由Ampere定律得到隨時間變化的微分力為

dF(t)=I(t)drc×B(t)。

(7)

其中B(t)為導電板中感應電流在上半空間激發(fā)的磁感應強度。對于時諧激勵源，我們關心的通常是一個周期內的平均力

(8)

其中*表示復共軛量。從而載流環(huán)在一個周期內所受的平均力將為

(9)

2任意方位矩形截面圓柱線圈感應場的確定

設上半非導電空間中有一載流圓環(huán)圓心為A，半徑為R，圓環(huán)軸線為l，所載電流為i=Icosωt。若取過點A且垂直于導電平面的直線為z軸，將l在z=0上的投影取為y軸，并設l與z=0交于點B，則可設OA=H，AB=a，并記AB與z軸的夾角為γ，如圖1所示。在這樣的坐標選擇下，我們得到源點矢量為

(10)

運用積分

(11)

(12)

其中I1(α)為第一類一階變形Bessel函數，我們得到矢量激勵普函數為

圖1　任意位置圓環(huán)的坐標系選取Fig.1　　　　Selection of coordinate system for arbitrarily　　　 located circular loop

(13)

其中

u=nsinγ-jpcosγ。

(14)

從而

(-exn+eym)dmdn。

(15)

矩形截面圓柱線圈的感應場可由對單個載流圓環(huán)的感應場進行積分得到。設有一導電半無窮空間上方的圓柱線圈，其匝數為N，內外半徑分別為R1，R2，高為2h，坐標系取法如圖1，這時線圈軸與z軸夾角為γ，線圈中心點A與坐標原點間的距離為H，AB=d。將式(15)對軸向與徑向積分，并運用積分

(16)

其中Ln(x)為n階變形Struve函數，則線圈產生的感應磁通密度為

(e2hpcosγ-1)e-j(mx+ny)-pzg1(u)dmdn。

(17)

(e2hpcosγ-1)e-j(mx+ny)-pzg1(u)dmdn。

(18)

(e2hpcosγ-1)e-j(mx+ny)-pzg1(u)dmdn。

(19)

3矩形截面圓柱線圈的受力計算

首先有

(20)

運用式(11)，式(12)得

(21)

(22)

(23)

(24)

(25)

其中

v=psinγ-jncosγ。

(26)

從而

(27)

(28)

(29)

(30)

(31)

將式(21)～式(25)中的I置換為NI/2h(R2-R1)，我們得到矩形截面圓柱線圈所受的力為

ez(-wy+jwxcosγ)]。

(32)

也就是說，導電板上方任意方位的圓柱線圈將只受到y與z方向的力，這是因為圓柱線圈關于x軸對稱。對于軸線平行于導電平面的圓柱線圈，有γ=π/2，式(32)中各函數均將趨于0/0的不定式情形，顯然此時僅有z向的力不為零，經過極限步驟

(33)

得到其值為

(34)

式(34)表明水平圓柱線圈的受力與線圈本身的長度無關。這是可以理解的，因為對無窮大金屬板而言水平軸線圈的每匝都處于同等位置，每匝線圈與金屬板的電磁相互作用特性是一樣的，因此電磁力與匝數有關而與線圈長度無關。

4檢驗：線圈軸垂直于導電平面的特殊情形

現以線圈軸垂直于導電平面的特殊情形來檢驗式(32)的正確性。在這種情形有γ=0。對這種軸對稱的場，應只有力的垂直分量而無水平分量。引入

(35)

(36)

從而Fy=0。另一方面有

(37)

此式與文獻[3]中給出的表達式一致。

5數值計算與比較

作為算例，取線圈參數H=0.1 m，h=0.03 m，R1=0.03 m，R2=0.08 m，N=500，激勵電流參數為I=5 A，ω=2 000π rad/s，金屬板參數κ=5.8×107S/m，μr=1(銅)，以式(32)計算Fy與Fz隨夾角γ的變化值并繪制曲線圖，并將此模型以有限元軟件CST Studio仿真的結果進行驗證，在仿真中分別取γ=0°，15°，30°，45°，60°，75°，90°。其結果如圖2所示，圖中的離散點即為CST的仿真計算結果。

圖2　圓柱線圈受力隨夾角γ的變化Fig.2　　　　Forces acting on the cylindrical coils　　　 dependent on the axial angle γ

由圖2可以看到，在我們所取的各項參數下，γ=0，π/2時Fy=0，且Fy有一最大值出現在γ=70°附近，隨著γ的變化，Fy將保持方向不變。而Fz在γ=60°附近有一零點，過零之后Fz將反向。

保持其他各參數不變，取κ=1×107S/m(普通鐵磁材料電導率)，我們觀察式(34)所描述的水平軸圓柱線圈的受力隨激勵電流角頻率ω變化的規(guī)律。分別繪制μr=1，μr=100，μr=1 000時的線圈受力曲線，并以將此模型以有限元軟件CST Studio仿真的結果進行驗證，在仿真中分別取ω=1 000，3 000，5 000，7 000，9 000 rad/s，其結果如圖3所示，圖中離散點為CST仿真結果。在此處的計算中假定μr為常數，即不考慮磁材料飽和以及磁滯的影響。

圖3　平行軸圓柱線圈受力隨電流角頻率ω的變化Fig.3　　　　Forces acting on the cylindrical coils 　　　dependent on the axial angle ω

由圖3可以看到，當μr=100，激勵電流角頻率ω約為4 000 rad/s時水平線圈受力平均值可以為零。μr=1，ω=0即線圈通以直流電時線圈受力為零，但在μr≠1的情形，ω=0時線圈受力并不為零。這實際上反映了鐵磁材料對通電線圈的磁力作用。

6結論

借助高階矢量勢，我們能夠解析地處理厚金屬板上方的任意放置的矩形截面圓柱線圈的受力問題。除了線圈軸垂直或平行于導電平面的特殊情形，該線圈將受到兩個方向的電磁作用力，其中水平力沿著線圈軸在z=0平面的投影方向，另一個力沿著垂直于導電平面的方向。當線圈軸垂直或平行于導電平面時，線圈將僅受到垂直于導電平面方向的作用力。隨著線圈軸與z軸夾角γ的變化，垂直方向力分量的時間平均值將在某個特定γ處變?yōu)榱恪τ讦蘲>1的金屬板，線圈所受平均電磁力可能對某個特定電流頻率變?yōu)榱恪ｋ娏鹘穷l率ω=0即激勵電流為直流時，僅當μr=1時線圈所受電磁力才會為零。

參 考 文 獻:

[1]RAJOTTE R J,DROUET M G.Experimental analysis of a fast acting circuit breaker mechanism—electrical aspects[J].IEEE Transactions on Power Apparatus and Systems,1975,94(1): 89-96.

[2]LIPINSKI W.über die abschirmung station?ter magnetfelder von zylinderspuler[J].Archiv für Elektrotechnik,1976,58: 9-14.

[3]LIPINSKI W,SIKORA R.Verluste und kraftangriff von wirbelstr?men im leitenden halbraum bei zylindersymmetrischen magnetischen wechselfeld[J].Archiv für Elektrotechnik,1975,57:173-175.

[4]HANNAKAM L.Wirbelstr?me im leitenden halbraum bei beliebiger form der erregenden leiterschleife[J].Archiv für Elektrotechnik,1972,54: 253-257.

[5]HANNAKAM L.Transienter skineffekt in einer platte endlicher dicke bei beliebiger form der erregenden leiterschleife[J].Archiv für Elektrotechnik,1978,60: 361-371.

[6]THEODOULIDIS T,SKARLATOS A.Eddy current interaction of an arbitrarily positioned probe coil with a conductive cylinder[J].IEEE Transactions on Magnetics,2012,48(8): 2392-2394.

[7]THEODOULIDIS T,DITCHBURN R J.Mutual impedance of cylindrical coils at an arbitrary position and orientation above a planar conductor[J].IEEE Transactions on Magnetics,2007,43(8): 3368-3370.

[8]THEODOULIDIS T.Analytical model for tilted coils in eddy-current nondestructive inspection[J].IEEE Transactions on Magnetics,2005,41(9): 2447-2454.

[9]MICOLAU G,PICHENOT G,PREMEL D,et al.Dyad-based model of the electric field in a conductive cylinder at eddy-current frequencies[J].IEEE Transactions on Magnetics,2004,40(2): 400-409.

[10]MIORELLI R.Eddy current modeling of narrow cracks in planar-layered metal structures[J].IEEE Transactions on Magnetics,2012,48(10): 2551-2559.

(編輯：劉琳琳)

Calculations for electromagnetic forces between thick metal plates and cylindrical coils of arbitrary location

LUO Yao

(School of Electrical Engineering,Wuhan University,Wuhan 430072,China)

Abstract:For an inclined circular coil with alternating current excitation and placed above a thick metal plate,analytical solution of electromagnetic forces between them can not be found by traditional methods.Using the high-order vector potential,the problem was solved analytically,and vector analysis and Struve functions were employed to simplify the final results.The expressions for the electromagnetic forces acting on the cylindrical coil of arbitrary position were obtained.The results indicate that there are at most two electromagnetic forces acting on the coils.When the axes of coils are perpendicular or parallel to the conducting surface,the results can be simplified to 1-D integrals and only perpendicular component exists.The numerical calculations were implemented to verify the obtained results.The results show that the time mean value of the vertical force component can vanish for certain inclined angles of coil axis,and the horizontal force component can attain maximum value for some other angles.If direct currents flow through the coil,the electromagnetic forces will vanish only when the plate has the same permeability with the air.

Keywords:cylindrical coils; electromagnetic force；high-order vector potential; Bessel functions; Struve functions

收稿日期:2013-05-04

基金項目：國家自然科學基金青年基金(51507115)

作者簡介:羅垚(1983—)，男，博士，研究方向為電磁場計算及數學物理方程。

通訊作者:羅垚

DOI:10.15938/j.emc.2016.05.004

中圖分類號:TM 12

文獻標志碼:A

文章編號:1007-449X(2016)05-0023-05