1、电磁学课程论文课程名称电磁学论文题目对通电螺线管磁场分布的若干问题的研究姓名学院系别专业年级2012级学号任课教师2013年12月1日1对通电螺线管磁场分布的若干问题的研究（中国科学技术大学近代力学系，安徽合肥230027）摘要本文对通电螺线管外部磁场是否为零进行了讨论，又对非密绕螺线管的磁场分布进行了计算。关键词通电螺线管磁场分布ARESEARCHINTOSOMEQUESTIONSOFTHEMAGNETICFIELDDISTRIBUTIONOFASOLENOIDWITHELECTRICCURRENTDEPARTMENTOFMODERNMECHANICS,UNIVERSITYOFSCIENCE
2、ANDTECHNOLOGYOFCHINA,ANHUIHEFEI230027ABSTRACTTHEARTICLEGIVESADISCUSSIONONWHETHERTHEMAGNETICFIELDOUTSIDEASOLENOIDISZEROANDICALCULATEDTHEDISTRIBUTIONOFANONETIGHTLYWOUNDSOLENOIDKEYWORDSASOLENOIDWITHELECTRICCURRENTTHEDISTRIBUTIONOFMAGNETICFIELD1引言螺线管的磁场是电磁学和电工技术必须面对的问题，一直有人关注。同时螺线管是用的最多的一种基本线圈形式，广泛用于军事、
3、经济、生态医疗、天文地质等众多领域中。所以对于螺线管空间磁场分布求解是有重要意义的。但一般电磁学教材对通电螺线管的磁场分布的讨论不够深入，对通电螺线管外部磁场是否为零的讨论比较模糊，而且一般仅限于无穷长密绕螺线管磁场的分布。对非密绕螺线管磁场分布的讨论很少。基于此，本文对通电螺线管的磁场分布进一步做了讨论。2无穷长密绕通电螺线管外磁场是否为零对于螺线管外的磁场，不少电磁学教材武2断得认为管外磁场为零。下面是一种论证无穷长通电管外磁场为零的方法。其思路是如图一所示的矩形回路ABCD,根据AMPERE环路定理，有0ABCDABBCCDDABDLBDLBDLBDLBDLNIAB已知轴线上（AB段）的
4、0BNI，故0ABBDLNIAB又因0,0BCDABDLBDL故要求CDBDLBCD管外管外结论是0B管外但是，上述论证是不严密的，因为他预先假定了CD路径上的B与DL同方向，事实上，由0BDL管外并不能得出B管外0，最多只能说明管外磁场的轴向分量为零。也有人提出了如下的证明方法。下面，让我们看看是如何证明的。1管外的磁场既不是Z的函数，也不是的函数，故管外磁场不可能是均匀磁场，只能是径向坐标R的函数。2管外磁场B管外的磁场的径向分量BR0，为了证明此结论，做如图二的同轴圆柱面。根据磁场的GAUSSS定理，通过上述圆柱面的磁通量为零，即DSDSDSDS0BBBB左底面右底面侧面图一由对称性，通
5、过左右底面的磁通量应互相抵消，故RDS0BDSBBDS侧面侧面式中，BR是管外磁场B的径向分量，故有BR0图2ABCD5ABCDZ33管外磁场的Z分量应为零，即BZ0,为了论证这一结论，做如图2所示的矩形回路ABCD，其中AB边和CD边与轴线的距离可任意。由对称性，因无净电流穿过次闭合回路，根据AMPER环路定理，有120ZZABABCDBCCDABBDLBDLBDLBDLBDL式中BZ1和BZ2分别是AB边和CD边处磁场的轴向分量。因在BC边和DA边的管内部分的磁场B与DL垂直，故对积分无贡献；而在BC边和DA边管外部分处，管外磁场的径向分量BR为零故对积分也无贡献。于是上式右端第二项和第四
6、项均为零。所以120ZZABCDBDLBDL又因AB边和CD边离轴线的距离是任意的，故管外磁场的Z分量只能是BZ1BZ2BZ0那么这种证发是否就无懈可击呢仔细分析，不难发现这种证法将螺线管的螺距忽略了，忽略了螺距，就不能叫螺线管。在张之翔编的电磁学千题解的第408页中的5218题给出了螺距不能忽略的螺线管模型，结果是，电流既有轴向的分量，也有垂直轴向的分量，垂直轴向的分量在管外激发的磁场为零，然而沿着轴向的分量在管外激发的磁场并不为零图3下面我们再来看一下螺距的大小与轴向电流的关系。首先我们看一下产生螺纹的一种方法。如图3所示，取一条倾角为的直角三角形条，绕在半径为R的圆筒上，纸条的斜边就构成
7、了螺旋线。当纸条的底边绕圆筒一周时，螺旋线推进的距离即为螺距。显然有TAN2DR单位长度的线圈数为1COT2NDR（1）如图4所示，对于密绕线圈，螺距即为导线的直径，在密绕条件下，相当于螺线管表面有一均匀的表面电流层。定义通过单4位长度线段（该单位长度位于电流面且与电流方向垂直）的电流强度为电流面密度J。（见课本P118最下面）由图4知J的大小为1COSJNI（2）图4式中N为轴向单位长度包含的线圈数，I为导线中的电流强度。如图4所示，J可以分为J轴向和J环向两个分量，前者沿螺线管的轴向，后者沿螺线管横截圆的切向，把（1）式带入（2）式中，得2SINIJR并有JSIN2JCOSCOT2IJRI
8、JR轴向环向轴向总电流为J2IRI轴向轴向通过单位长度线段的环向电流为JCOTN2IIIR环向环向由此可见，轴向电流I轴向与螺距无关，螺距再小I轴向也不会为零，且总等于I，所以忽略螺距得出的结果是不能使人信服的。螺线管外的磁场相当于轴向一根电流强度为I的导线产生的磁场，方向为环绕轴向电流的方向。下面让我们带入实际的数字来比较螺线管外部磁场和内部磁场大小设电流为1A，螺线管半径为01M，单位长度的线圈数为500。则80402102R6283210IBBTBNIT环管外内可以发现B内是远大于B管外的。因此在一定的条件下B管外可以忽略不计，但是螺线管外的磁场却不等于零。3有限长非密绕通电螺线管产生的
9、磁场在轴线上的分布对于无穷长密绕通电螺线管大家都有了比较清晰的了解，下面我们来讨论非密绕通5电螺线管磁场的空间分布。如图5所示是非密绕螺旋导线，螺旋半径为A，轴向单位长度的圈数为N,螺距为，显然有1N图5设通过螺旋导线的电流强度为I，为了求轴上观察点O的磁感应强度，用如下方法设置直角坐标系以观察点O为坐标原点，落线管的中心轴为Z轴，用过O点并与Z轴垂直的平面与螺旋线的交点为Q，连接O和Q两点的直线为X轴，XY平面与Z轴垂直。在螺旋导线的P处取线元DL,其位矢为R，R在XY平面内的投影在图5中用虚线表示，该虚线与X轴的夹角为，于是，P点的坐标即R（X,Y,Z为COSSIN22XAYAZN线元DL
10、的三个分量即DL（DLX,DLY,DLZ为SINCOS22DLXDXADDLYDYADDLZDZN根据BIOTSAVRT定律，O点的磁场为20314PIPDLRBR（3）式中P1和P2分别代表螺线管的左端和右端，I是螺线管中的电流强度，R是电流元IDL到O点所引的矢量，它是电流元位矢R的关系为RR即COSSIN2XYZRXARYARZN因COSSIN22SINCOS2XZYYZDLRRDLRDLDADANNADNCOSSIN22COSSIN2YXZZXDLRRDLRDLDAADNNADNXYZQIDLP1IP2622222SINCOSZYXXYDLRRDLRDLADADAD1122222221
11、2222122RXYZANNN代入（3）式，得出磁场B的三个分量为2121212130X223/220223/220223/2320223/2SINCOS2422SINCOS2COSSIN2242YZABNDNNAINADNABINADNAIDBNANA式中，1和2为螺旋线左右端的方位角，由上式知，若已知1和2则可以求出O点的磁场。下面我们以一个实际的模型进行计算，设螺旋半径A005M,单位长度的线圈匝数为50，真空中的磁导率为04107亨利/米，螺线管中的电流强度为1A。设1222，由于求出上面三式的解析解比较困难，故我们采取数值积分的方法。利用MATLAB编写龙贝格求积分的程序，并进行计算
12、，最后得出结果为2223/22SINCOS500D0XB2223/22COSSIN50002689D）2Y6000026891327110BINA2223/221500030D325020003023335104ZIBNAT下面让我们看一下当螺线管直径不变时随着角的增大BY和BZ是如何变化的。为此，我们令1和2的大小相等，都等于。用MATLAB画出BY和BZ与的关系如下。05010015020025015105005115X106BY与的关系曲线YB图一与的关系曲线7ZB图二与的关系曲线Y,ZBB图三与角大小的关系曲线4无穷长密绕通电螺线管在螺线管表面产生的磁场大小。由于面电流所在处的任一点的
13、磁感强度B，等于从该面两边趋于该点时磁感强度的极限值（B和B）的平均。即1B2BB（这是一个普遍规律，这个规律与面电荷所在处电场强度的规律相同。【参见张之翔，电磁学教学札记（高等教育出版社，1987）】又由于螺线管外部磁场的大小与内部磁场的大小相比可以忽略。因此螺线管磁场表面的磁场近似等于内部磁场大小的1/2。即01N2BI表面5结论由于任何无穷长通电螺线管的螺距都不能忽略，因此，通电螺线管便会存在沿着轴线的电流分布，从而导致通电螺线管外部的磁场不为零，本文还给出了非密绕通电螺线管内轴线上某一点的磁场强度的三个分量的积分表示，并通过给相关参数赋一组特殊值求得BX,BY,BZ,的数值解。最后绘图
14、表示了BY,BZ,与的关系曲线图。附龙贝格积分的MATLAB源代码FUNCTIONR,QUAD,ERR,HROMBERF,A,B,N,DELTAF是被积函数A,B分别是积分的上下限N1是T数表的列数DELTA是允许误差R是T数表QUAD是所求积分值M1HBAERR1J0RZEROS4,4R1,1HFEVALF,AFEVALF,B/2WHILEERRDELTAHH/2S0FORP1M050100150200250101234567X105BY与的关系曲线BZ与的关系曲线05010015020025001234567X105BZ与的关系曲线8XAH2P1SSFEVALF,XENDRJ1,1RJ,1
15、/2HSM2MFORK1JRJ1,K1RJ1,KRJ1,KRJ,K/4K1ENDERRABSRJ,JRJ1,K1ENDQUADRJ1,J1求磁场BY,BZ随角变化的MATLAB源代码YZEROS1,101ZZEROS1,101T1FORI12201TYTVPAINTCOSXXSINX/5PI2X23/2,I,IYTYT4PI107PI2500005EYTZTVPAINT1/5PI2X23/2,I,IZTZT1072PI5030052ZTTT1END参考文献【1】胡友秋程福臻叶邦角电磁学与电动力学上册（科学出版社）【2】张之翔电磁学千题解（科学出版社）【3】张之翔电磁学教学札记（高等教育出版社）【4】陈秉乾舒幼生胡望雨电磁学专题研究（高等教育出版社）【5】赵凯华,陈熙谋电磁学上册（高等教育出版社）【6】哈里德D,瑞斯尼克R物理学第二卷、第二册（科学出版社）【7】张德丰MATLAB数值计算方法（机械工业出版社）9
