康玉寬 李眾立

(西南科技大學(xué)制造科學(xué)與工程學(xué)院 四川綿陽 621000)

?

兩同軸通電圓柱形螺線管間電磁作用力分析和計(jì)算

康玉寬 李眾立

(西南科技大學(xué)制造科學(xué)與工程學(xué)院 四川綿陽 621000)

通過建立數(shù)學(xué)模型，運(yùn)用電磁場理論、數(shù)值計(jì)算方法、微積分知識(shí)推導(dǎo)出兩通電圓柱形螺線管間互感，進(jìn)一步運(yùn)用虛功原理推導(dǎo)出兩同軸通電圓柱形螺線管間電磁作用力的表達(dá)式，著重分析和研究了兩同軸通電等大的圓柱形螺線管間電磁作用力的計(jì)算方法和思想。利用Matlab數(shù)學(xué)軟件對(duì)實(shí)際的等大的兩個(gè)同軸圓柱形螺線管電磁感線分布情況進(jìn)行了模擬，進(jìn)一步驗(yàn)證計(jì)算結(jié)果的正確性，并模擬了部分電磁作用力和變量的關(guān)系圖像。

電磁作用力 互感 虛功原理 關(guān)系圖像 電磁舉升系統(tǒng)

隨著新能源概念的提出，電能的有效利用在各個(gè)領(lǐng)域都有了新的發(fā)展，作為電能應(yīng)用的載體之一的電磁能應(yīng)用尤為突出，自上海磁懸浮列車通車運(yùn)行成功以來，人們對(duì)電磁能的研究再次達(dá)到新的高峰。分析研究兩同軸通電螺線管間電磁作用力計(jì)算方法，在實(shí)際問題和技術(shù)應(yīng)用中有著重要的理論和工程價(jià)值。

在工程自卸車的舉升系統(tǒng)中，由于液壓舉升系統(tǒng)技術(shù)成熟、實(shí)用性較強(qiáng)，所以應(yīng)用較為普遍，但是隨著石油能源不斷枯竭，加之液壓舉升系統(tǒng)的不穩(wěn)定性，在未來的一段時(shí)間內(nèi)，電磁舉升系統(tǒng)將逐漸取代液壓舉升系統(tǒng)在工程自卸車舉升系統(tǒng)中的地位。因此，對(duì)兩同軸通電螺線管間電磁作用力的研究就尤為重要。雖然實(shí)際應(yīng)用中有一些經(jīng)驗(yàn)計(jì)算或者近似估計(jì)也能滿足一般工程要求，但是這些由經(jīng)驗(yàn)得來的估算結(jié)果在很大程度上不能有效控制計(jì)算的精確度，而且不能滿足特殊的工程要求。

1 模型建立

1.1 通電圓柱形螺線管的自感

我們可以將通電圓柱形螺線管看成由許多同軸通電圓環(huán)組成，通過分析兩平行通電圓環(huán)互感來推導(dǎo)通電圓柱形螺線管的自感。首先由紐曼公式推導(dǎo)兩平行通電圓環(huán)互感，再根據(jù)微積分原理取圓柱形螺線管的微分單元，求解兩微分單元間的互感。最后在整個(gè)圓柱形螺線管上求積分，就可以準(zhǔn)確求得圓柱形螺線管的自感。

(1)

(2)

圖1 圓柱形螺線管Fig.1 Cylindrical solenoid

1.2 兩通電圓柱形螺線管的互感

分別在兩通電圓柱形螺線管上微分單元dZ1和dZ2,同求單個(gè)通電圓柱形螺線管的自感的微積分方法，再利用數(shù)值計(jì)算方法化簡，進(jìn)而求解兩通電圓柱形螺線管的互感。設(shè)兩圓柱形螺線管尾端和前端間距為b，半徑分別為R1，R2,在微分單元dZ1和dZ2上分別均勻地繞有線圈m，n匝，磁導(dǎo)率為μ。基本尺寸如圖2所示。則兩通電圓柱形螺線管的互感為[3]：

(3)

圖2 圓柱形螺線管Fig.2 Cylindrical solenoid

1.3 兩同軸通電圓柱形螺線管間電磁作用力計(jì)算

兩通電圓柱形螺線管間電磁作用力滿足牛頓第三定律，該電磁作用力僅有沿軸方向上的分量。

由電磁場理論可知，該系統(tǒng)的場能為[4-6]：

W=I1I2M

(4)

I1，I2分別是兩通電圓柱形螺線管的電流，M是兩通電圓柱形螺線管的互感。載流體系在外磁場中的受力計(jì)算沿軸方向上的分量應(yīng)用虛功原理則有

(5)

(6)

至此，討論兩同軸通電等大的圓柱形螺線管間電磁作用力，即R1=R2=R。

其中：

{h1∈(0,b),h2∈(b,c),h3∈(b,a+b),

h4∈(0,b),i=1,2,3,4}

(7)

完全橢圓積分[7-8]的導(dǎo)數(shù)分別為：

(8)

(9)

則：

(10)

(11)

(12)

(8)式和(9)式代入(12)式化簡得：

(13)

將(7)代入(13)化簡得[9-10]：

(14)

分析表達(dá)式(14)可以知道兩同軸通電圓柱形螺線管間電磁作用力是由4部分組成，而且這4部分的解析表達(dá)式格式相同，每一部分對(duì)應(yīng)一個(gè)h,而且h的取值范圍在表達(dá)式(7)中列出。在這里，當(dāng)F為正值時(shí)，表示斥力；當(dāng)F為負(fù)值時(shí)，表示吸引力。推導(dǎo)得到的結(jié)果是一個(gè)級(jí)數(shù)解，具體應(yīng)用中只能取級(jí)數(shù)的若干項(xiàng)進(jìn)行近似計(jì)算，因此，從這個(gè)角度來說，表達(dá)式得出的是近似解。

2 計(jì)算結(jié)果的數(shù)值分析和圖形模擬

2.1 計(jì)算結(jié)果的正確性驗(yàn)證

應(yīng)用Matlab軟件模擬實(shí)際的等大的兩個(gè)同軸圓柱形螺線管電磁感線分布情況，進(jìn)一步驗(yàn)證計(jì)算結(jié)果的正確性，如圖3-圖6所示。

圖3 螺線管中間截面的磁感線分布圖Fig.3 Distribution diagram of magnetic induction lines on solenoid’s middle section

圖4 螺線管中間截面的磁感線趨勢圖Fig.4 Tendency diagram of magnetic induction lines on solenoid’s middle section

圖5 螺線管上下兩個(gè)面和中間面的截面圖Fig.5 The sectional views of solenoid’s top and bottom sides and middle part

圖6 螺線管中間截面的磁感線立體圖Fig.6 Space diagram of magnetic induction lines of solenoid’s middle section

利用Matlab軟件計(jì)算兩個(gè)螺線管的磁感線分布，結(jié)合表達(dá)式(14)進(jìn)行如下分析。圖3為不考慮兩個(gè)螺線管電磁相互作用力后螺線管中間截面的磁感線分布圖，圖4為考慮兩個(gè)螺線管電磁相互作用力后螺線管中間截面的磁感線趨勢圖，圖中方框代表螺線管的位置。通過圖3、圖4的對(duì)比可以知道，考慮兩個(gè)螺線管電磁相互作用力后，圖4中磁感線的彎曲程度大于圖3中磁感線的彎曲程度。說明磁感線的排斥力加大，兩個(gè)螺線管之間中心區(qū)域磁感應(yīng)強(qiáng)度變小，圖5給出兩螺線管上下兩個(gè)面和中間面的截面圖，這個(gè)圖從立體的角度說明了兩個(gè)螺線管磁感線互斥，且中間磁感應(yīng)強(qiáng)度變小。圖6是圖4的立體圖，當(dāng)兩個(gè)管的磁感線相互接近時(shí)由于排斥力的作用會(huì)使磁感線的方向改變，再次驗(yàn)證了兩個(gè)相對(duì)螺線管中心磁感應(yīng)強(qiáng)度變小。綜上分析計(jì)算結(jié)果符合實(shí)際情況。

2.2 表達(dá)式的近似解和變量h的關(guān)系驗(yàn)證

利用C程序設(shè)計(jì)、Matlab數(shù)學(xué)軟件[11]計(jì)算并模擬變量作出圖像。兩同軸通電圓柱形螺線管間電磁作用力由解析表達(dá)式(14)可以看出，作用力F與該兩同軸通電圓柱形螺線管間距離h有關(guān)，在其他條件不變的情況下，作用力F隨著h的變化而變化。

圖7 兩圓柱形螺線管間電磁力和距離的仿真圖Fig.7 The simulation diagram of the electromagnetic force and distance between two cylindrical solenoids

由于兩同軸通電圓柱形螺線管間電磁作用力的4部分的解析表達(dá)式格式相同，可知兩同軸通電等大的圓柱形螺線管間電磁作用力的其他三部分在距離h不同取值情況下，電磁作用力F和距離h的關(guān)系利用Matlab數(shù)學(xué)軟件模擬的圖形和圖7所示相似,在此不作具體Matlab軟件模擬。

3 結(jié)論

文章有效利用前人在螺線管間互感以及通電平行環(huán)的電磁力計(jì)算方面已取得的研究成果,基于磁場的疊加原理,分析研究螺線管的空間磁場分布，理論推導(dǎo)出同軸通電圓柱形螺線管間電磁作用力解析表達(dá)式。利用Matlab軟件模擬實(shí)際的等大的兩個(gè)同軸圓柱形螺線管電磁感線分布情況，進(jìn)一步驗(yàn)證了計(jì)算結(jié)果符合實(shí)際情況。

[1] 於乾英,富知愚.圓柱型螺線管磁體的最佳形狀[J].上海大學(xué)學(xué)報(bào):自然科學(xué)版,2002,8(1)：76-79.

[2] 張永吉.兩個(gè)共軸螺線管間互感的計(jì)算[J].大學(xué)物理,2000,5(5)：15-20.

[3] 洪正平.共軸密繞串聯(lián)螺線管之間的互感系數(shù)[J].山東師范大學(xué)學(xué)報(bào):自然科學(xué)版,2005,20(2)：101-102.

[4] 劉成林.大學(xué)物理[M].南京：南京大學(xué)出版社,2014.150-240.

[5] 趙凱華, 陳熙謀.電磁學(xué)，第2版[M].北京:高等教育出版社, 2006.494-502.

[6] 侯昭武.通電螺線管產(chǎn)生的電磁場[J].廣西教育學(xué)院學(xué)報(bào),2008,(5)：154-156.

[7] 田碩,田睿.完全橢圓積分在物理學(xué)中的應(yīng)用[J].中國西部科技,2009,8(7)：90-92.

[8] ZHOU Li-ming，QIU Song-liang,WANG Fei. Estimates for Complete Elliptic Integral by Trigonometric Functions with Applications[N]. Journal of Zhejiang Sci-Tech University，2012，29(5):722-726.

[9] 朱思華, 楊萬民.兩共軸平行圓環(huán)電流之間的電磁作用力[J].大學(xué)物理,2005,24(10)：24-26.

[10] 康兆敏，方秀男.數(shù)值計(jì)算方法,第二版[M].清華大學(xué)出版社,2013.

[11] (美)穆爾.MATLAB實(shí)用教程,第二版[M].高會(huì)生，劉童娜，李聰聰譯. 電子工業(yè)出版社,2010.

Analysis and Calculation of Electromagnetic Force between Two Coaxial Electric Cylindrical Solenoids

KANG Yu-kuan, LI Zhong-li

(SchoolofManufacturingScienceandEngineering,SouthwestUniversityofScienceandTechnology,Mianyang621010,Sichuan,China)

By establishing a mathematical model and using the electromagnetic field theory, the numerical computation method and calculus technique, etc., the mutual inductance between the two solenoids was obtained. Furthermore, the formula of expressing the above-said electromagnetic force was achieved by adopting the virtual work principle, with the computing method and ideology mainly analyzed and researched. By using the Matlab mathematical software, the distribution of the two coaxial electric cylindrical solenoids’ electromagnetic lines was simulated and the computed result was further verified. The image of relation between the electromagnetic force and their corresponding variables was also simulated.

Electromagnetic force; Mutual inductance; Virtual work principle; Image of relationship; Electromagnetic lifting system

2014-08-11

康玉寬，男，碩士研究生。E-mail:1173020132@qq.com

O369

A

1671-8755(2015)01-0085-05