常规磁铁设计

1、常规磁铁设计目录常规磁铁设计1一、加速器中的磁铁1二、二极磁铁与四极磁铁的物理原理3三、二极磁铁工程设计6四、四极磁铁工程设计12五、最终设计结果17六、附录 程序代码19一、加速器中的磁铁加速器技术发展至今，对于带电粒子运动的控制一直都是由磁铁产生的电磁场来完成。带电粒子的输运、聚焦、引入引出等，无不需要磁铁。可以说磁铁的设计是加速器设计中最重要的部分之一。1930年Earnest O. Lawrence和他的学生建造第一台回旋加速器时，直径只有10cm，所用的磁铁也很简单。随着加速器技术的迅速发展，对磁铁的要求也越来越高，如前苏联的“杜布纳”加速器，需要36000吨磁铁，耗电量达到148M

2、W。在Betatron电子加速器中，要控制使中心处的磁场变化是轨道处磁场变化的二倍，即著名的“二比一定律”，而且同时要考虑带电粒子运动的聚焦，这就需要对磁铁进行特殊的设计。如今，现代的大型高能加速器对磁铁的要求也更加的高，其中为代表的类型是世界上普遍建造的同步加速器。本文对同步加速器中所需的电磁铁进行了初步的设计与讨论。同步加速器的基础是强聚焦原理。首先，在介绍强聚焦原理之前，有必要先介绍一下弱聚焦原理。带电粒子在圆形加速器中运动时，为使束流能稳定存在，必须对粒子束进行聚焦。建立柱坐标系如图一，定义n为磁场降落指数下标c代表粒子闭轨处，下标z表示轴向方向的磁感应强度。带电粒子沿轴向和和横向方向

3、的运动方程为分别为对这两个方程求解，其结果表明，当n0时，沿r方向运动稳定。其中习惯称沿r方向为沿径向方向，沿z方向为轴向方向。综合起来，只有当0n1时，整个横向运动才满足稳定。为使磁场降落指数满足上述条件，就需要对磁铁进行特殊处理，如将极面设计为倾斜的状态等。随着高能物理的发展，对粒子能量的追求也越来越高，弱聚焦原理越来越不能满足高能加速器的需求，其中磁铁过重，耗电过大均是很难解决的问题。1958年库兰特提出强聚焦原理，利用将n0的磁铁进行交替排布以实现聚焦，不再拘泥于0n1，推动了加速器的飞速发展。其中所用磁铁包括四极磁铁和二极磁铁。两种四极磁铁产生均匀梯度的磁场分布，进行交替排列，对径向

4、和轴向交替聚焦；二极磁铁产生匀强磁场，使粒子弯转。其中同步加速器的前期物理设计的主要工作就是设计出磁铁的排列方式，以使带电粒子稳定的被加速。二、二极磁铁与四极磁铁的物理原理对磁多极矩进行分析。在二维恒定磁场中，磁标势满足拉普拉斯方程如下将展开成多项式其中分别为i次项，写出一次项和二次项如下根据拉普拉斯方程，应有有=-。一次项和二次项分别代表了二极矩与四极矩磁场，即分别为匀强磁场和均匀梯度磁场。对于二极矩在加速器中，主环一般是水平放置的，我们通常希望只有垂直方向的磁场，即仅有。此时磁标势就可以为即可看出，此时磁标势等位线是和水平方向平行的，即二极磁铁的可以利用两个接近的极面，在极面间的间隙中近似

5、产生匀强磁场。对于四极矩再次求导得到磁感应强度梯度磁标势表达式即可看出四极矩可以由极面为双曲线的磁铁产生，且产生的磁场磁感应强度梯度为常数。类似的，我们还可以得到更高阶的六极矩、八极矩等。因为二极铁与四极铁相对与其他磁铁要为重要，同时又因为二者都有较强的对称性，故本文仅对二极磁铁与四极磁铁的横截面进行二维的设计。三、二极磁铁工程设计二极磁铁通常有C型、H型、框型等，由于C型磁铁一侧是开口的，便于其他部件的安装与维修，故设计采用C型磁铁。将上文lattice设计所得参数列入表中。如表1，表2。其中，好场区宽度指满足均匀度的磁铁极面宽度；磁间隙指极面间隙的高度。表1 内环二极磁铁参数好场区宽度80

6、 mm磁间隙g50 mm好场区磁感应强度B1.4 T好场区均匀度0.0001 T偏转角度7.5度长度1.0917 m曲率半径8.3403 m表2 外环二极磁铁参数好场区宽度80 mm磁间隙g50 mm好场区磁感应强度B1.224 T好场区均匀度0.0001 偏转角度7.5度长度1.2488 m曲率半径9.5403 m其中内环外环二极磁铁均为平行边磁铁。首先计算磁铁的安匝数。整个磁路由磁铁和一段气隙组成，由电磁学磁路定理可得安匝数F其中f为修正系数，f=1.051.10；B为好场区磁感应强度，单位为高斯；g为气隙高度，单位为cm。取f=1.1，可得安匝数F=61600安培由于二极磁铁的横截面是二

7、维对称的，所以仅研究上半部分即可，即取安匝数的一半。采用外截面为6X6mm2，内方孔4X4mm2，内通冷却水的线材，其横截面为20 mm2，在冷却水正常冷却情况下，允许通过的最大电流密度为100A/mm2，即每根导线允许通过的最大电流为200A。一般要求实际电流要小于最大电流。取设计电流为170A，则由上文所求安匝数可求得需要的线匝数N为为使绕制的线圈接线和冷却方便，故将线圈设计为绕制偶数排与偶数层，即10排18层共180匝。考虑导线之间的绝缘物件，导线绕制要间距1mm，线包外也要留2mm厚的绝缘层，再考虑留下0.5mm的安装空间。最终可得到线圈横截面的尺寸大小，尺寸设计为130X75mm2.

8、我们知道，在两个无限大平板磁极间是匀强磁场，但工程是不可能做到这一点的。但为了保证一定的磁感应强度均匀度，极面应在好场区宽度的基础上向外延长一定的宽度。根据国内外积累的建造加速器的经验，当均匀度达到万分之一时，即可满足加速器的要求。利用半经验公式，我们可估计所需极面要向外延伸的长度其中即为磁场均匀度，g为磁间隙高度。带入数据可估计极面应向外延长70mm，则极面宽度应为140+80=220mm。根据线圈的大小和极面的宽度，我们可以设计出磁铁的轮廓如下磁铁的材质采用硅钢片，其B-H数据见附录。由于硅钢的物理性质，当磁感应强度很强时，将会达到饱和状态。如下图为避免局部磁场感应强度饱和，在轮廓设计中，

9、在磁铁的内部转角处采用了圆角。而在磁铁的外围由于本身磁场较弱，设计圆角只会增加加工成本，故未采用圆角。接下来将磁铁的形状参数、安匝数、硅钢片B-H线数据写入possion程序的输入文件，进行计算。其中possion程序的数值计算方法为有限差分法，网格划分为三角形网格。将输入文件运行后，查看磁铁中是否有磁感应强度超过饱和值，然后对磁感应强度进行分析，确定好场区范围是否达到80mm，若未达到，则对极面进行垫补处理，最后调整安匝数，使磁感应强度为1.4T.本次设计中，第一次计算中有部分区域磁场饱和，进行调整圆角的半径大小；之后好场区仅有60mm左右，经多次调整，最终确定安匝数为30002，单根电线电

10、流为167A；增加垫补如图，单个垫补宽10mm，高0.5mm。最后计算结果如图所示其中红色的线为磁感线。分析计算数据，调出磁间隙的磁感应强度数据如图如图可见，在51cm-59cm的80mm宽度内，磁感应强度上下波动不超过万分之一，满足了好场区宽度的要求。对于外环设计与内环类似，鉴于减少工厂设计建造费用，故直接采用内环的磁铁尺寸，降低电流达到外环磁感应强度参数即可。计算结果如图可见基本达到设计要求。除安匝数为25310外，其他参数与内环相同。四、四极磁铁工程设计据lattice设计结果，列出内环散焦、聚焦磁铁，外环散焦、聚焦磁铁参数如下表表三 四极磁铁目标参数内环聚焦磁铁聚焦常数KF0.8725

11、15磁感应强度梯度10.1879 T/m梯度均匀度0.0001磁铁长度0.6 m内环散焦磁铁聚焦常数KD-0.846984磁感应强度梯度9.8897 T/m梯度均匀度0.0001磁铁长度0.6 m外环聚焦磁铁聚焦常数0.840794磁感应强度梯度9.8174 T/m梯度均匀度0.0001磁铁长度0.6 m外环散焦磁铁聚焦常数-0.817091磁感应强度梯度-9.5407 T/m梯度均匀度0.0001磁铁长度0.6 m其中，磁感应强度梯度为聚焦常数与磁刚度的乘积，根据前文，磁刚度数值为当四极磁铁极头的横截面顶部曲线严格满足双曲线时，磁铁梯度也严格为常数，但与二极磁铁的设计类似，工程上是无法实现严

12、格的双曲线的，故在保证一定精度下，采用双曲线和直线结合的几何结构。首先对安匝数进行计算，由磁路定理求得单个磁极的安匝数其中1.1为计算过程中的修正系数。代入数据可得安匝数=11000，同二极磁铁的设计过程，利用此公式计算出的安匝数还需要进一步的微调才能满足设计要求。设计线圈结构为6层10排共60匝，单根导线电流约为180A。依然采用前述导线。则线圈横截面尺寸即为60mmX36mm。接下来结合线圈尺寸设计磁铁几何结构。先由孔径大小50mm确定四极铁中心内接圆的大小，为保证一定精度，取内接圆半径为100mm。双曲线断开的位置即B点坐标由经验公式x=（1.31.5）a来确定，双曲线断开后设计为一段切

13、线段BC加一段平行线段CD。其中平行线段CD可减少双曲线两端的顶端效应，类似于二极磁铁的垫补。切线段长设计为5mm，平行线段长度由经验取6mm。同样，在磁铁内侧设计圆角以避免磁场饱和。初步设计好磁铁几何尺寸后，写入输入文件进行计算，由于四极磁铁的对称性，所以仅取八分之一部分进行计算即可。根据程序计算结果微调几何结构，经过多次调整切线与平行线的长度及安匝数，最终确定四极磁铁的几何尺寸如图软件的计算结果如图调出四种四极磁铁磁感应强度梯度图线如图以上各图中左端当接近于x=0时，出现数据发散的情况，这是因为程序在对磁感应强度求梯度时，当分母接近于0时误差被放大所致。实际情况下，越靠近图像的左端，梯度是

14、越均匀的。由以上各图，可见在025mm内，磁感应强度梯度的均匀度达到了万分之一标准。为减少设计成本，四个四极磁铁采用同一几何结构，调整线圈电流改变安匝数即可分别达到各个磁铁对磁感应强度梯度大小的要求。五、最终设计结果二极磁铁内环弯转磁铁安匝数30002 A导线电流166.67 A内环弯转磁铁安匝数25310 A导线电流140.61 A四极磁铁内环聚焦磁铁安匝数10930 A导线电流182.17 A内环散焦磁铁安匝数10608 A导线电流176.80 A外环聚焦磁铁安匝数10535 A导线电流175.58 A外环散焦磁铁安匝数10235 A导线电流170.58 A参考文献1、李泉凤.2002.

15、电磁场数值计算与电磁铁设计 M.北京：清华大学出版社，149-184.2、裴元吉 王相綦.2012. 加速器物理学.合肥：中国科学技术大学课程讲义，63-94.六、附录 程序代码二极铁® kprob=0, ! Poisson or Pandira problemmode=0, ! Some materials have variable permeabilityxreg1=4.5,kreg1=2,xreg2=26.5, kreg2=58,xreg3=29,kreg3=68,xreg4=42, kreg4=124,xreg5=44.5, kreg5=134,xreg6=66.5, kreg

16、6=190,kmax=200,yreg1=3, lreg1=20,yreg2=11, lreg2=54,yreg3=14, lreg3=74,yreg4=35.5, lreg4=138,lmax=140 &po x=0,y=0 &po x=87,y=0 &po x=87,y=43&po x=0,y=43 &po x=0,y=0 ® mat=1,cur=30002 &po x=29,y=3.5 &po x=29,y=11 &po x=42,y=11 &po x=42,y=3.5 &po x=29,y=3.5® mat=1,cur=-30002 &po x=66,y=3.5&po x=

17、66,y=11&po x=79,y=11&po x=79,y=3.5&po x=66,y=3.5® mat=2,mtid=1,mshape=0 &po x=5, y=0&po x=27,y=0&po x=27, y=7&po nt=2,x0=33,y0=7,x=0,y=6&po x=38,y=13&po nt=2,x0=38,y0=7,x=6,y=0 &po x=44,y=4.5 &po nt=2,x0=46,y0=4.5,x=0,y=-2 &po x=47,y=2.5 &po x=47,y=2.45 &po x=48,y=2.45 &po x=48,y=2.5 &po x=62,y=2

18、.5 &po x=62,y=2.45 &po x=63,y=2.45 &po x=63,y=2.5 &po x=64, y=2.5 &po nt=2,x0=64,y0=4.5,x=2,y=0 &po x=66, y=35&po x=5, y=35&po x=5, y=0&mt mtid=1bgam=0.00000 0.0017513135900 0.001747079950 0.0017417421000 0.0017354981050 0.0017283091100 0.001720141150 0.0017109631200 0.0017007531250 0.0016894941300

19、0.0016771741350 0.0016637862800 0.0010806942850 0.0010680512900 0.0010561422950 0.0010449123000 0.0010343093050 0.0010242893100 0.0010148093150 0.0010058283200 0.0009973123250 0.0009892263300 0.0009815393350 0.0009742224000 0.0009049524500 0.0008567985000 0.0008184935500 0.0007880856000 0.0007642026

20、500 0.0007458637000 0.0007323767500 0.0007232618000 0.0007182098500 0.0007170549000 0.0007197589500 0.00072641110000 0.00073723110500 0.00075259410578 0.000756258011319 0.000795102211940 0.000837520912451 0.000883470312912 0.000929368013313 0.000976467113654 0.001025325513935 0.001076426314216 0.001

21、125492414447 0.001176747514618 0.001231360314789 0.001284686515020 0.001331557915131 0.001387925115252 0.001442377015432 0.001491201915594 0.001538935115705 0.001591849716180 0.001854255516840 0.002375296917150 0.002915451917360 0.003456619417620 0.003972983717830 0.004486316718200 0.005494505518950

22、 0.007917656419500 0.010256410320200 0.014858841020650 0.019379845020950 0.023866348421600 0.037037037021900 0.045662100523000 0.086956521723386 0.100281000023850 0.118163000024408 0.138742000025079 0.162246000025885 0.188858000026854 0.218695000028019 0.2517840000 &四极铁® kprob=0, ; Poisson or Pan

23、dira problemmode=-1 ; Materials have fixed permeabilitydx=.05,dy=.05, ; Mesh intervals yminf=0,ymaxf=15 ; Fixed Y for field interpolationxminf=0,xmaxf=15 ; X range for field interpolation; The next 6 terms refer to the harmonic analysis:ktype=4, ; Quadrupole symmetrynterm=10, ; Number of coefficient

24、snptc=10, ; Number of arc points for interpolationrint=7.0, ; Radius of the arc for interpolationangle=45, ; Angular extent of arc (default start = 0)rnorm=1.0 & ; Aperture radius for normalization&po x=0.0,y=0.0 &po x=34,y=0.0 &po x=34,y=34 &po x=0.0,y=0.0 ® mat=3,mtid=3 &po x=3.5355,y=3.5355 &p

25、o nt=3,r=5,x=7,y=1.7857 &po x=8.049,y=1.627 &po x=10.355,y=1.627 &po x=19.137,y=5.146&po nt=2,x0=20.253,y0=2.362,x=3,y=0 &po x=23.253,y=0.0 &po x=32.253,y=0 &po x=32.253,y=12.362 &po x=22.375,y=22.375&po x=3.5355,y=3.5355 ® mat=1,cur=10930 &po x=15.819,y=3.6 &po x=21.819,y=3.6&po x=21.819,y=0&po

26、x=15.819,y=0&po x=15.819,y=3.6 ® ibound=0 &po x=34,y=34 &po x=3.5355,y=3.5355 &po x=0.0,y=0.0 &mt mtid=3bgam=0.00000E+00 0.0017513135 ! Start of B,Gamma data0.11420E+04 0.00175131350.29530E+04 0.00101595040.51140E+04 0.00078216660.84760E+04 0.00070786440.96670E+04 0.00072411300.10578E+05 0.00075625800.11319E+05 0.00079510220.11940E+05 0.00083752090.12451E+05 0.00