用matlab模擬環(huán)形磁鐵的磁場(chǎng)分布

1、matlab模擬環(huán)形磁鐵磁場(chǎng)分布摘要：和地球內(nèi)部的磁感線分布類似，環(huán)形磁鐵圓環(huán)中心的磁感線是垂直于環(huán)形平面的直線，其余的按距離環(huán)由近及遠(yuǎn)由環(huán)繞環(huán)的磁感線漸漸伸展成和中心平行的直線，越靠近中心的越像直線向兩極伸展。為了能夠形象的刻畫(huà)，我們使用matlab 強(qiáng)大的計(jì)算能力做了描述。關(guān)鍵字：matlab、環(huán)形磁鐵、磁感應(yīng)線、分子電流、安培環(huán)路定理matlab simulation of the magnetic field distribution of ring magnetsabstract: similar to the earths interior distribution of magn

2、etic induction lines, magnetic induction lines in the center of the ring magnet is perpendicular to the ring plane ,by the remaining distance from the near to the distant ring around the ring by the magnetic sense of line and centers gradually extended into a straight line parallel to, the more near

3、 the center more like a straight line extending to the poles.keyword: matlab、ring magnet line of magnetic induction、 molecular electric current、ampere ring circuit theorem一 引言作為一種人工磁化而制成的磁鐵，環(huán)形磁鐵有其本身特別的優(yōu)點(diǎn)和用處，研究它的磁場(chǎng)分布對(duì)了解環(huán)形磁鐵的性質(zhì)有著重要的意義。而為了對(duì)環(huán)形磁鐵磁場(chǎng)的分布情況進(jìn)行形象的描述，需要借助于計(jì)算機(jī)強(qiáng)大的計(jì)算和模擬功能，利用安培分子電流假說(shuō)，建立起數(shù)學(xué)模型，并利用mat

4、lab數(shù)學(xué)軟件進(jìn)行數(shù)學(xué)計(jì)算和形象的描述，可以得到理論情況下直觀的電磁場(chǎng)的大小和方向的位置分布示意圖。二 matlab pde工具箱matlab是矩陣實(shí)驗(yàn)室（matrix laboratory）的簡(jiǎn)稱，和mathematica、maple并稱為三大數(shù)學(xué)軟件。它在數(shù)學(xué)類科技應(yīng)用軟件中在數(shù)值計(jì)算方面首屈一指。matlab可以進(jìn)行矩陣運(yùn)算、繪制函數(shù)和數(shù)據(jù)、實(shí)現(xiàn)算法、創(chuàng)建用戶界面、連接其他編程語(yǔ)言的程序等，主要應(yīng)用于工程計(jì)算、控制設(shè)計(jì)、信號(hào)處理與通訊、圖像處理、信號(hào)檢測(cè)、金融建模設(shè)計(jì)與分析等領(lǐng)域。pde工具箱是一個(gè)基于matlab的用來(lái)研究與求解偏微分方程的工具集合。它為二維空間上的偏微分方程的研究與求

5、解提供了一個(gè)強(qiáng)大而靈活的工作環(huán)境。工具箱對(duì)偏微分方程采用有限元方法(fem)進(jìn)行離散，提供如下功能：(1)定義一個(gè)偏微分方程(pde)問(wèn)題，包括確定二維研究區(qū)域、邊界條件和方程系數(shù)；(2)采用數(shù)值方法求解pde問(wèn)題，包括在研究區(qū)域上生成剖分網(wǎng)格，并在其上離散方程和產(chǎn)生近似解；(3)圖形化顯示求解結(jié)果。pde工具箱的基本方程是橢圓型偏微分方程：三 環(huán)形磁鐵磁場(chǎng)的計(jì)算環(huán)形磁鐵的形狀和結(jié)構(gòu)決定了其內(nèi)部磁場(chǎng)分布，由于環(huán)形磁鐵的立體性使得它的磁場(chǎng)分布與環(huán)形線圈不同，但是它們二者又有很多相似之處，我們可以仿照環(huán)形線圈的磁場(chǎng)計(jì)算方法。利用安培分子電流假說(shuō)和環(huán)路積分定理建立數(shù)學(xué)模型。磁感應(yīng)強(qiáng)度是描述磁場(chǎng)性質(zhì)的

6、重要物理量，磁感應(yīng)線能形象地表示磁場(chǎng)在空間的分布。若能通過(guò)多媒體手段將磁感應(yīng)線形象地展現(xiàn)出來(lái)，就能將抽象轉(zhuǎn)為直觀，更有助我們研究環(huán)形磁鐵的性質(zhì)和應(yīng)用。為了解決這一難題，我們運(yùn)用 matlab 強(qiáng)大的數(shù)學(xué)運(yùn)算和繪圖功能，通過(guò)對(duì)環(huán)形磁鐵的磁感應(yīng)強(qiáng)度矢量分布圖進(jìn)行計(jì)算機(jī)仿真，形象地展現(xiàn)出了磁感應(yīng)強(qiáng)度的空間分布圖。圖1 環(huán)形磁鐵1） 基本原理安培環(huán)路積分定理：式中：真空磁導(dǎo)率 磁介質(zhì)的相對(duì)磁導(dǎo)率 磁介質(zhì)的絕對(duì)磁導(dǎo)率（簡(jiǎn)稱磁導(dǎo)率） 分子電流的假說(shuō)是指在原子、分子等物質(zhì)微粒內(nèi)部，存在著一種環(huán)形電流分子電流，分子電流使每個(gè)物質(zhì)微粒都成為微小的磁體，它的兩側(cè)相當(dāng)于兩個(gè)磁極。仿效直流線圈的磁場(chǎng)計(jì)算公式我們可以得

7、到分子電流形成的磁場(chǎng)，一般磁鐵可以看成是很多個(gè)方向一致的的分子磁場(chǎng)的疊加。直流導(dǎo)線圈產(chǎn)生磁場(chǎng)的基本規(guī)律：上式為環(huán)形線圈軸線上的磁場(chǎng)分布情況，對(duì)于軸線以外的空間區(qū)域，應(yīng)該用比奧薩瓦爾定律：通過(guò)積分來(lái)確定。2）模型建立物理模型由于環(huán)形磁鐵的磁場(chǎng)分布與磁化形狀有關(guān)，我們針對(duì)這樣的磁鐵進(jìn)行研究，即磁鐵平放時(shí)上下為兩磁極。假設(shè)磁鐵被極化的很好，近全部分子電流沿同一個(gè)方向（逆時(shí)針），并且假設(shè)假設(shè)分子環(huán)形電流緊密挨著，分子電流之間電流相互抵消，這樣可以簡(jiǎn)化為如下模型：環(huán)形磁鐵電流密度分層分布，最外層為逆時(shí)針設(shè)為im，最內(nèi)層為in，之間由于抵消分子電流很小設(shè)為iz。如圖3。圖2 環(huán)形磁鐵俯視圖圖3簡(jiǎn)化后電流分

8、布計(jì)算模型 沿線圈軸向取出一個(gè)截面來(lái)分析，如圖4所示。在這個(gè)截面中，包含個(gè)區(qū)域，r1+r4+r5和r6+r2+r7為環(huán)形磁鐵區(qū)域邊緣處電流密度較大，中間較小。其他區(qū)域?yàn)樗诳臻g?，F(xiàn)在認(rèn)為是空氣中。圖4 環(huán)形磁鐵縱截面圖r4與r7電流相等，為im，r5與r6電流密度相等為in ,r1與r2電流密度相等為iz靜磁場(chǎng)的控制方程用雙旋度方程來(lái)描述：1a=j在平面問(wèn)題中我們假定電流方向平行于z軸，于是僅有磁矢位存在z分量：a=(0，0，a)，j=(0，0，j)雙旋度方程可以簡(jiǎn)化為標(biāo)量橢圓型偏微分方程：1a=j式中j=j（x y）。這樣在一個(gè)截面內(nèi)的磁場(chǎng)線分布，經(jīng)過(guò)旋轉(zhuǎn)360后，形成的三維分布就是正真的環(huán)

9、形磁鐵的磁場(chǎng)立體分布。四matlab程序設(shè)計(jì)及實(shí)現(xiàn)通過(guò)使用matlab pde工具箱解的的二維平面圖5所示：圖5 環(huán)形磁鐵磁場(chǎng)分布的橫截面圖由圖中可以直觀的看出截面磁場(chǎng)的分布，圖中用磁感應(yīng)線描述磁感應(yīng)線分度的密度也表征了磁感應(yīng)強(qiáng)度的大小，而其方向可以由圖中箭頭描述。旋轉(zhuǎn)360后類似下圖（理想圖）圖3 三維理想效果圖與三維情況下的理想效果圖比較可知模擬結(jié)果能較好的符合真實(shí)情況。五結(jié)束語(yǔ) 本次模擬中模型的建立過(guò)程中用到了一些猜想，和理想化處理。因此所刻畫(huà)的環(huán)形磁鐵磁場(chǎng)分布并非真實(shí)情況，在此僅作一種理想的求解介紹。并且環(huán)形磁鐵的磁化方式有多種，磁場(chǎng)分布都會(huì)有很大的不同，因此仍需繼續(xù)努力，找到更好的辦

10、法，來(lái)細(xì)膩的刻畫(huà)。參考文獻(xiàn):1 朱衛(wèi)娟.matlab在物理教學(xué)中的應(yīng)用.2 陳懷琛.matlab及其在理工課程中的應(yīng)用指南.西安電子科技大學(xué)出版社.2003年3 陳熙謀、趙凱華.電磁學(xué).高等教育出版社.1985年6月，第2版4 matlab工具箱應(yīng)用/蘇金明, 張蓮花, 劉波等編附錄一：編程代碼：function pdemodelpde_fig,ax=pdeinit;pdetool(appl_cb,6);set(ax,dataaspectratio,8 5 1);set(ax,plotboxaspectratio,1 1 1);set(ax,xlim,-8 8);set(ax,ylim,-5

11、5);set(ax,xtickmode,auto);set(ax,ytickmode,auto); % geometry description:pderect(-2 -1 0.25 -0.25,r1);pderect(1 2 0.25 -0.25,r2);pderect(-8 8 5 -5,r3);pderect(-2 -1.8999999999999999 0.25 -0.25,r4);pderect(-1.1000000000000001 -1 0.25 -0.25,r5);pderect(1 1.1000000000000001 0.25 -0.25,r6);pderect(1.899

12、9999999999999 2 0.25 -0.25,r7);set(findobj(get(pde_fig,children),tag,pdeeval),string,r1+r2+r3+r4+r5+r6+r7) % boundary conditions:pdetool(changemode,0)pdesetbd(9,.neu,.1,.0,.0)pdesetbd(8,.neu,.1,.0,.0)pdesetbd(2,.neu,.1,.0,.0)pdesetbd(1,.neu,.1,.0,.0)% mesh generation:setappdata(pde_fig,hgrad,1.3);se

13、tappdata(pde_fig,refinemethod,regular);setappdata(pde_fig,jiggle,char(on,mean,);pdetool(initmesh)pdetool(refine)pdetool(refine)% pde coefficients:pdeseteq(1,.1./(1.0)!1./(1.0)!1./(1.0)!1./(1.0)!1./(1.0)!1./(1.0)!1./(1.0),.0.0!0.0!0.0!0.0!0.0!0.0!0.0,.0!-8!8!2!-0.75!-2!0.75,.1.0!1.0!1.0!1.0!1.0!1.0!1.0,.0:10,.0.0,.0.0,.0 100)setappdata(pde_fig,currparam,.1.0!1.0!1.0!1.0!1.0!1.0!1.0;.0!-8!8!2!-0.75!-2!0.75 ) % solve parameters:setappdata(pde_fig,solveparam,.str2mat(0,35280,10,pdeadworst,.0.5,longest,0,1e-4,fixed,inf)% plotflags and user data strings:setappd