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

1、燕山大學(xué)課程設(shè)計(jì)論文MATLAB模擬環(huán)形磁鐵磁場(chǎng)分布（By r oy ）摘要：和地球內(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 fielddistribution of ring magnetsAbstract： Similar to the earth's int

2、erior distribution of magnetic 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

3、parallel to, the more near the center more like a straight line extending to the polesKeyword: MATLAB . Ring magnet Line of magnetic induction、 Molecular electric currents 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ì)算和模

4、擬功能，利用安培分子電流假說(shuō)，建立起 數(shù)學(xué)模型，并利用MATLAB數(shù)學(xué)軟件進(jìn)行數(shù)學(xué)計(jì)算和形象的描述， 可以得到理論情況下直觀的電磁場(chǎng)的大小和方向的位置分布示意 圖。二 MATLAB pde 工具箱MATLAB是矩陣實(shí)驗(yàn)室(Matrix Laboratory)的簡(jiǎn)稱，和Mathe matica. 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è)基于MATL

5、AB的用來(lái)研究與求解偏微分方程 的工具集合。它為二維空間上的偏微分方程的研究與求解提供了一 個(gè)強(qiáng)大而靈活的工作環(huán)境。工具箱對(duì)偏微分方程采用有限元方法(F EM)進(jìn)行離散，提供如下功能：(1) 定義一個(gè)偏微分方程(PDE)問(wèn)題，包括確定二維研究區(qū)域、邊界 條件和方程系數(shù)；(2) 采用數(shù)值方法求解PDE問(wèn)題，包括在研究區(qū)域上生成剖分網(wǎng)格， 并在其上離散方程和產(chǎn)生近似解；(3) 圖形化顯示求解結(jié)果。PDE工具箱的基本方程是橢圓型偏微分方程：m三環(huán)形磁鐵磁場(chǎng)的計(jì)算環(huán)形磁鐵的形狀和結(jié)構(gòu)決定了其內(nèi)部磁場(chǎng)分布，由于環(huán)形磁鐵 的立體性使得它的磁場(chǎng)分布與環(huán)形線圈不同，但是它們二者又有很 多相似之處，我們可以仿照

6、環(huán)形線圈的磁場(chǎng)計(jì)算方法。利用安培分 子電流假說(shuō)和環(huán)路積分定理建立數(shù)學(xué)模型。磁感應(yīng)強(qiáng)度是描述磁場(chǎng) 性質(zhì)的重要物理量，磁感應(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)3燕山大學(xué)課程設(shè)計(jì)論文度的空間分布圖。4燕山大學(xué)課程設(shè)計(jì)論文#燕山大學(xué)課程設(shè)計(jì)論文圖1 環(huán)形磁鐵1）基本原理安培環(huán)路積分定理：L式中：真空磁導(dǎo)率 /o=4-xlO-7T.m/.磁介質(zhì)的相對(duì)磁導(dǎo)率 pr磁

7、介質(zhì)的絕對(duì)磁導(dǎo)率（簡(jiǎn)稱磁導(dǎo)率）分子電流的假說(shuō)是指在原子、分子等物質(zhì)微粒內(nèi)部，存在著一 種環(huán)形電流分子電流，分子電流使每個(gè)物質(zhì)微粒都成為微小的 磁體，它的兩側(cè)相當(dāng)于兩個(gè)磁極。仿效直流線圈的磁場(chǎng)計(jì)算公式我 們可以得到分子電流形成的磁場(chǎng)，一般磁鐵可以看成是很多個(gè)方向 一致的的分子磁場(chǎng)的疊加。直流導(dǎo)線圈產(chǎn)生磁場(chǎng)的基本規(guī)律：吐1，軸線二%IS上式為環(huán)形線圈軸線上的磁場(chǎng)分布情況，對(duì)于軸線以外的空間區(qū) 域，應(yīng)該用比奧一薩瓦爾定律:Idl xAtv r2通過(guò)積分來(lái)確定。2）模型建立物理模型由于環(huán)形磁鐵的磁場(chǎng)分布與磁化形狀有關(guān)，我們針對(duì)這樣的磁 鐵進(jìn)行研究，即磁鐵平放時(shí)上下為兩磁極。假設(shè)磁鐵被極化的很好，近全部

8、分子電流沿同一個(gè)方向（逆時(shí) 針），并且假設(shè)假設(shè)分子環(huán)形電流緊密挨著，分子電流之間電流相 互抵消，這樣可以簡(jiǎn)化為如下模型：環(huán)形磁鐵電流密度分層分布， 最外層為逆時(shí)針設(shè)為im，最內(nèi)層為in，之間由于抵消分子電流很小 設(shè)為120如圖3。圖3簡(jiǎn)化后電流分布計(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)為是空氣中。R4R7R1R5R6R2圖4環(huán)形磁鐵縱截面圖R4與R7電流相等，為h, R5與R6電流密度相等為I. ,R1與R2電流密度相等為L(zhǎng)靜磁場(chǎng)的控制方程用雙旋度方

9、程來(lái)描述:Vx在平面問(wèn)題中我們假定電流方向平行于z軸，于是僅有逆矢位存在Z分量：A = (0, 0, A), J= (0, 0, J)雙旋度方程可以簡(jiǎn)化為標(biāo)量橢圓型偏微分方程:式中 J二J (x y)o這樣在一個(gè)截面內(nèi)的磁場(chǎng)線分布，經(jīng)過(guò)旋轉(zhuǎn)360。后，形成的三維 分布就是正真的環(huán)形磁鐵的磁場(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。后類似下圖（理想圖）與三維情況下的理想效

10、果圖比較可知模擬結(jié)果能較好的符合 真實(shí)情況。五.結(jié)束語(yǔ)本次模擬中模型的建立過(guò)程中用到了一些猜想，和理想化處 理。因此所刻畫(huà)的環(huán)形磁鐵磁場(chǎng)分布并非真實(shí)情況，在此僅作一種 理想的求解介紹。并且環(huán)形磁鐵的磁化方式有多種，磁場(chǎng)分布都會(huì) 有很大的不同，因此仍需繼續(xù)努力，找到更好的辦法，來(lái)細(xì)膩的刻 畫(huà)C參考文獻(xiàn)：1 朱衛(wèi)娟.MATLAB在物理教學(xué)中的應(yīng)用2 陳懷琛.« MATLAB及其在理工課程中的應(yīng)用指南.西安電子科技大學(xué) 出版社.2003年3 陳熙謀、趙凱華.電磁學(xué).高等教育出版社.1985年6月，第2版MATLAB I具箱應(yīng)用/蘇金明，張蓮花，劉波等編10燕山大學(xué)課程設(shè)計(jì)論文附錄一：編程代

11、碼：function pdemode1 pde_figz ax=pdeinit;pdetool(1appl_cb1,6);set(axz ataAspectRatio1, 8 5 1);set ( ax z 1PlotBoxAspectRat i o1z 1 1 1 );set(axXLimJ-8 8 );set(axYLimJ-5 5 );set (axz 1 XTickMode1z * auto 1 );set( axz 1YT ickMode1, 1 auto 1 );% Geometry description:pderect(-2 -1 0.25 -0.25,Rl1);pderect

12、(1 2 0.25 -0.25,1R21);pderect( - 8 8 5 -5 , 1R31);JR4;,*R5')；'R6')；'R7 »);1PDEEval1 ) / 1 Strpderect(-2 -1.8999999999999999pderect( -1.1000000000000001 -1 pderect(1 1.10 0 0 0 0 0 0 0 0 0 0 0 0 01 0 .pderect(1.8999999999999999 2 0.25 025, set (findobj (get(pde_fig，Children* ) ,

13、1 Tag 1 z ing1,1R1+R2+R3+R4+R5+R6+R71)% Boundary conditions:pdetool ( 1 changemode' 9 0)pdesetbd(9 1 neu 0 J1 0 1 )pdesetbd(8,* neu '1 / 1 0 1 z . . 0 1 )pdesetbd (2 1neu ' , 1, 'O' 0 1 )pdesetbd (11neu '1, 0mW w w 11IJI% Mesh generation:'regular ');on * f 1 mean 9 f

14、' 9 );setappdata(pde_f igz 1Hgrad1 , 1.3 );setappdata(pde_f ig, 're f ineme thod1 , setappdata(pde_f igz 'j iggle' ,char ( pdetool(1 initmesh1) pdetool(1 ref ine1 ) pdetool ( 1refine 1 )% PDE coef ficients :10燕山大學(xué)課程設(shè)計(jì)論文1./ (1.0) !1./ (1.0) ! 1 . / (1 . 0) ! 1 . / (1.0) !1 . / (1.0) !

15、1./(1.0) !1./ (1.0) ' z00!00!00!0.0!00!00!00J 0!-8!8!2!-0 75!-2!0 75 ',1.0!1.0!1.0!1.0!1.0!1.0!1.0, / * 0 : 101, 0 0 ' 0 0 ' 1 0 100 1 ) setappdata(pde_f igz 1currparam1 f1.0!1.0!1.0!1.0!1.0!1.0!1.0,;. 0!-8!8!2!-0.75!-2!0.751 ) % Solve parameters:setappdata(pde_figz 'solveparam 1.s t r2mat ( 1 0 1 , 1 3 52 8 0 1 z,pdeadworstl / .1 0.5 1 z 1 longest 1 z 101 , 1 IE-4 1 z 1 1 z 1 fixed1 z 1 Inf 1 )% Plotflag