用matlab模拟环形磁铁的磁场分布综述

1、燕山大学课程设计论文MATLAB模拟环形磁铁磁场分布摘要：和地球内部的磁感线分布类似，环形磁铁圆环中心的磁感线是垂 直于环形平面的直线，其余的按距离环由近及远由环绕环的磁感线渐渐伸 展成和中心平行的直线，越靠近中心的越像直线向两极伸展。为了能够形 象的刻画，我们使用matlab 强大的计算能力做了描述。关键字：MATLAB、环形磁铁、磁感应线、分子电流、安培环路定理MATLAB simulation of the magnetic fielddistribution of ring magnetsAbstract: Similar to the earths interior distribu

2、tion of magnetic inductionlin es, magn etic in ducti on lines in the cen ter of the ring magnet is Perpendicular to the ring plane ,By the remaining distanee from the near to the distant ring around the ring by the magneti c sense of line and centers gradually extended into a straight line parallel

3、to, the more near the center more like a straight li ne exte nding to the poles.Keyword : MATLAB 、 Ring magnet Line of magn eticin duct ion、Molecular electric current、 Ampere ring circuit theorem一引言作为一种人工磁化而制成的磁铁，环形磁铁有其本身特别的优 点和用处，研究它的磁场分布对了解环形磁铁的性质有着重要的意 义。而为了对环形磁铁磁场的分布情况进行形象的描述，需要借助 于计算机强大的计算和模拟功能

4、，利用安培分子电流假说，建立起 数学模型，并利用MATLAB数学软件进行数学计算和形象的描述， 可以得到理论情况下直观的电磁场的大小和方向的位置分布示意 图。MATLAB pde 工具箱MATLAB是矩阵 实验室(Matrix Laboratory )的 简称，和 Mathe matica、Map le并称为三大数学软件。它在数学类科技应用软件中 在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数 和数据、实现算法、创建用户界面、连接其他编程语言的程序等， 主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信 号检测、金融建模设计与分析等领域。PDE工具箱是一个基于MATLA

5、B的用来研究与求解偏微分方程 的工具集合。它为二维空间上的偏微分方程的研究与求解提供了一 个强大而灵活的工作环境。工具箱对偏微分方程采用有限元方法(F EM)进行离散，提供如下功能：定义一个偏微分方程(PDE)问题，包括确定二维研究区域、边界 条件和方程系数； 采用数值方法求解PDE问题，包括在研究区域上生成剖分网格， 并在其上离散方程和产生近似解；图形化显示求解结果。PDE工具箱的基本方程是椭圆型偏微分 方程：三环形磁铁磁场的计算环形磁铁的形状和结构决定了其内部磁场分布，由于环形磁铁 的立体性使得它的磁场分布与环形线圈不同，但是它们二者又有很 多相似之处，我们可以仿照环形线圈的磁场计算方法。

6、利用安培分 子电流假说和环路积分定理建立数学模型。磁感应强度是描述磁场 性质的重要物理量，磁感应线能形象地表示磁场在空间的分布。若 能通过多媒体手段将磁感应线形象地展现出来，就能将抽象转为直 观，更有助我们研究环形磁铁的性质和应用。为了解决这一难题， 我们运用MATLAB强大的数学运算和绘图功能，通过对环形磁铁的 磁感应强度矢量分布图进行计算机仿真，形象地展现出了磁感应强 度的空间分布图。图1 环形磁铁1）基本原理安培环路积分定理：B d =忖I内L式中：真空磁导率 =4二10=Tm/A磁介质的相对磁导率磁介质的绝对磁导率（简称磁导率）卩=亠0分子电流的假说是指在原子、分子等物质微粒内部，存在

7、着一 种环形电流分子电流，分子电流使每个物质微粒都成为微小的 磁体，它的两侧相当于两个磁极。仿效直流线圈的磁场计算公式我 们可以得到分子电流形成的磁场，一般磁铁可以看成是很多个方向 一致的的分子磁场的叠加。直流导线圈产生磁场的基本规律：B中心=皿2R2n R2x212上式为环形线圈轴线上的磁场分布情况，对于轴线以外的空间区 域，应该用比奥萨瓦尔定律：d B4 n r2通过积分来确定。2）模型建立物理模型由于环形磁铁的磁场分布与磁化形状有关，我们针对这样的磁 铁进行研究，即磁铁平放时上下为两磁极。假设磁铁被极化的很好，近全部分子电流沿同一个方向（逆时 针），并且假设假设分子环形电流紧密挨着，分子

8、电流之间电流相 互抵消，这样可以简化为如下模型：环形磁铁电流密度分层分布， 最外层为逆时针设为i m，最内层为in，之间由于抵消分子电流很小 设为i z。如图3。图2环形磁铁俯视图图3简化后电流分布计算模型沿线圈轴向取出一个截面来分析，如图4所示。在这个截面中， 包含个区域，R1+R4+R5和R6+R2+R7为环形磁铁区域边缘处电流密 度较大，中间较小。其他区域为所在空间。现在认为是空气中。R4R7R1R5R6R2图4环形磁铁纵截面图R4与R7电流相等，为Im, R5与R6电流密度相等为I n ,R1与R2电流密度相等为I静磁场的控制方程用双旋度方程来描述：在平面问题中我们假定电流方向平行于z

9、轴，于是仅有磁矢位存在z分量：A=(0, Oi A), = (0, 0, I)双旋度方程可以简化为标量椭圆型偏微分方程:式中 J=J （ x y ）o这样在一个截面内的磁场线分布，经过旋转360。后，形成的三维 分布就是正真的环形磁铁的磁场立体分布。四MATLAB程序设计及实现通过使用matlab PDE 工具箱解的的二维平面图5所示：CdrlEdur. A VacEdn Hie-ld.日OQ-图5 环形磁铁磁场分布的横截面图由图中可以直观的看出截面磁场的分布，图中用磁感应线描述 磁感应线分度的密度也表征了磁感应强度的大小，而其方向可以由 图中箭头描述。旋转360。后类似下图（理想图）图3 三

10、维理想效果图与三维情况下的理想效果图比较可知模拟结果能较好的符合 真实情况。五结束语本次模拟中模型的建立过程中用到了一些猜想，和理想化处 理。因此所刻画的环形磁铁磁场分布并非真实情况，在此仅作一种 理想的求解介绍。并且环形磁铁的磁化方式有多种，磁场分布都会 有很大的不同，因此仍需继续努力，找到更好的办法，来细腻的刻 画。参考文献:1 朱卫娟.MATLAB 在物理教学中的应用.2 陈怀琛.MATLAB 及其在理工课程中的应用指南西安电子科技大学出版社.2003年3 陈熙谋、赵凯华.电磁学.高等教育出版社.1985年6月，第2版4 MATLAB 工具箱应用/苏金明，张莲花，刘波等编附录编程代码：f

11、un cti onpdemodelpde_fig,ax=pdeinit;pdetool(appl_cb,6);set(ax,DataAspectRatio,8 5 1);set(ax,卩 lotBoxAspectRatio,1 1 1);set(ax,XLim，-8 8);set(ax,YLim，-5 5);set(ax,XTickMode, auto);set(ax,YTickMode, auto);% Geometry descripti on:pderect(-2 -1 0.25 -0.25,R1);pderect(1 2 0.25 -0.25,R2)i;pderect(-8 8 5 -

12、5,R3);pderect(-2 -1.8999999999999999 0.25-0.25,R4);pderect(-1.10000000000000011 0.25 -0.25,R5);pderect(1 1.1000000000000001 0.25-0.25,R6);pderect(1.8999999999999999 2 0.25-0.25,R7);set(fi ndobj(get(pde_fig,Childre n),Tag,PDEEval),String , R1+R2+R3+R4+R5+R6+R7)% Bou ndary con diti ons:pdetool(cha nge

13、mode,0)pdesetbd(9,.neu1,0,.0)pdesetbd(8,.neu,.1,. .O,.O)pdesetbd(2,.neu,.1,. .O,.O)pdesetbd(1,.neu,.1,. .O,.O)Hgrad ,1.3);refi nemethodjiggle,char(),regular);on , mean， );% Mesh gen erati on: setappdata(pde_fig, setappdata(pde_fig, setappdata(pde_fig, pdetool(in itmeshpdetool(refi nepdetool(refi ne%

14、 PDE coefficie nts:pdeseteq(1,1./(1.0)!1./(1.0)!1./(1.0)!1./(1.0)!1./(1.0)!1./(1.0)!1.,Inf )1 1 1 0 0 0 1 0 1）；）；/(1.0) ,.0.0!0.0!0.0!0.0!0.0!0.0!0.00!-8!8!2!-0.75!-2!0.751.0!1.0!1.0!1.0!1.0!1.0!1.00:10 ,.0.0,.0.0,.0 100 )setappdata(pde_fig,1.0!1.0!1.0!1.0!1.0!1.0!1.00!-8!8!2!-0.75!-2!0.75% Solve parameters:setappdata(pde_fig,str2mat( 0, 352800.5, Iongest, 0% Plotflags and user data stri ngs:setappdata(pde_fig,0 1 0 1);seta