其中所用磁铁包括四极磁铁和二极磁铁。
两种四极磁铁产生均匀梯度的磁场分布,进行交替排列,对径向和轴向交替聚焦;二极磁铁产生匀强磁场,使粒子弯转。
其中同步加速器的前期物理设计的主要工作就是设计出磁铁的排列方式,以使带电粒子稳定的被加速。
二、二极磁铁与四极磁铁的物理原理
对磁多极矩进行分析。
在二维恒定磁场中,磁标势
满足拉普拉斯方程如下
将φ展开成多项式
其中
分别为i次项,写出一次项和二次项如下
根据拉普拉斯方程,应有有
=-
。
一次项和二次项分别代表了二极矩与四极矩磁场,即分别为匀强磁场和均匀梯度磁场。
对于二极矩
在加速器中,主环一般是水平放置的,我们通常希望只有垂直方向的磁场,即仅有
。
此时磁标势就可以为
即可看出,此时磁标势等位线是和水平方向平行的,即二极磁铁的可以利用两个接近的极面,在极面间的间隙中近似产生匀强磁场。
对于四极矩
再次求导得到磁感应强度梯度
磁标势表达式
即可看出四极矩可以由极面为双曲线的磁铁产生,且产生的磁场磁感应强度梯度为常数。
类似的,我们还可以得到更高阶的六极矩、八极矩等。
因为二极铁与四极铁相对与其他磁铁要为重要,同时又因为二者都有较强的对称性,故本文仅对二极磁铁与四极磁铁的横截面进行二维的设计。
三、二极磁铁工程设计
二极磁铁通常有C型、H型、框型等,由于C型磁铁一侧是开口的,便于其他部件的安装与维修,故设计采用C型磁铁。
将上文lattice设计所得参数列入表中。
如表1,表2。
其中,好场区宽度指满足均匀度的磁铁极面宽度;磁间隙指极面间隙的高度。
表1内环二极磁铁参数
好场区宽度
80mm
磁间隙g
50mm
好场区磁感应强度B
1.4T
好场区均匀度
≤0.0001T
偏转角度
7.5度
长度
1.0917m
曲率半径
8.3403m
表2外环二极磁铁参数
好场区宽度
80mm
磁间隙g
50mm
好场区磁感应强度B
1.224T
好场区均匀度
≤0.0001
偏转角度
7.5度
长度
1.2488m
曲率半径
9.5403m
其中内环外环二极磁铁均为平行边磁铁。
首先计算磁铁的安匝数。
整个磁路由磁铁和一段气隙组成,由电磁学磁路定理可得安匝数F
其中f为修正系数,f=1.05~1.10;B为好场区磁感应强度,单位为高斯;g为气隙高度,单位为cm。
取f=1.1,可得安匝数F=61600安培
由于二极磁铁的横截面是二维对称的,所以仅研究上半部分即可,即取安匝数的一半。
采用外截面为6X6mm^2,内方孔4X4mm^2,内通冷却水的线材,其横截面为20mm^2,在冷却水正常冷却情况下,允许通过的最大电流密度为100A/mm^2,即每根导线允许通过的最大电流为200A。
一般要求实际电流要小于最大电流。
取设计电流为170A,则由上文所求安匝数可求得需要的线匝数N为
为使绕制的线圈接线和冷却方便,故将线圈设计为绕制偶数排与偶数层,即10排18层共180匝。
考虑导线之间的绝缘物件,导线绕制要间距1mm,线包外也要留2mm厚的绝缘层,再考虑留下0.5mm的安装空间。
最终可得到线圈横截面的尺寸大小,尺寸设计为130X75mm^2.
我们知道,在两个无限大平板磁极间是匀强磁场,但工程是不可能做到这一点的。
但为了保证一定的磁感应强度均匀度,极面应在好场区宽度的基础上向外延长一定的宽度。
根据国内外积累的建造加速器的经验,当均匀度达到万分之一时,即可满足加速器的要求。
利用半经验公式,我们可估计所需极面要向外延伸的长度
其中
即为磁场均匀度,g为磁间隙高度。
带入数据可估计极面应向外延长70mm,则极面宽度应为140+80=220mm。
根据线圈的大小和极面的宽度,我们可以设计出磁铁的轮廓如下
磁铁的材质采用硅钢片,其B-H数据见附录。
由于硅钢的物理性质,当磁感应强度很强时,将会达到饱和状态。
如下图
为避免局部磁场感应强度饱和,在轮廓设计中,在磁铁的内部转角处采用了圆角。
而在磁铁的外围由于本身磁场较弱,设计圆角只会增加加工成本,故未采用圆角。
接下来将磁铁的形状参数、安匝数、硅钢片B-H线数据写入possion程序的输入文件,进行计算。
其中possion程序的数值计算方法为有限差分法,网格划分为三角形网格。
将输入文件运行后,查看磁铁中是否有磁感应强度超过饱和值,然后对磁感应强度进行分析,确定好场区范围是否达到80mm,若未达到,则对极面进行垫补处理,最后调整安匝数,使磁感应强度为1.4T.本次设计中,第一次计算中有部分区域磁场饱和,进行调整圆角的半径大小;之后好场区仅有60mm左右,经多次调整,最终确定安匝数为30002,单根电线电流为167A;增加垫补如图,单个垫补宽10mm,高0.5mm。
最后计算结果如图所示
其中红色的线为磁感线。
分析计算数据,调出磁间隙的磁感应强度数据如图
如图可见,在51cm-59cm的80mm宽度内,磁感应强度上下波动不超过万分之一,满足了好场区宽度的要求。
对于外环设计与内环类似,鉴于减少工厂设计建造费用,故直接采用内环的磁铁尺寸,降低电流达到外环磁感应强度参数即可。
计算结果如图
可见基本达到设计要求。
除安匝数为25310外,其他参数与内环相同。
四、四极磁铁工程设计
据lattice设计结果,列出内环散焦、聚焦磁铁,外环散焦、聚焦磁铁参数如下表
表三四极磁铁目标参数
内环聚焦磁铁
聚焦常数KF
0.872515
磁感应强度梯度
10.1879T/m
梯度均匀度
≤0.0001
磁铁长度
0.6m
内环散焦磁铁
聚焦常数KD
-0.846984
磁感应强度梯度
9.8897T/m
梯度均匀度
≤0.0001
磁铁长度
0.6m
外环聚焦磁铁
聚焦常数
0.840794
磁感应强度梯度
9.8174T/m
梯度均匀度
≤0.0001
磁铁长度
0.6m
外环散焦磁铁
聚焦常数
-0.817091
磁感应强度梯度
-9.5407T/m
梯度均匀度
≤0.0001
磁铁长度
0.6m
其中,磁感应强度梯度为聚焦常数与磁刚度的乘积,根据前文,磁刚度数值为
当四极磁铁极头的横截面顶部曲线严格满足双曲线时,磁铁梯度也严格为常数,但与二极磁铁的设计类似,工程上是无法实现严格的双曲线的,故在保证一定精度下,采用双曲线和直线结合的几何结构。
首先对安匝数进行计算,由磁路定理求得单个磁极的安匝数
其中1.1为计算过程中的修正系数。
代入数据可得安匝数
=11000,同二极磁铁的设计过程,利用此公式计算出的安匝数还需要进一步的微调才能满足设计要求。
设计线圈结构为6层10排共60匝,单根导线电流约为180A。
依然采用前述导线。
则线圈横截面尺寸即为60mmX36mm。
接下来结合线圈尺寸设计磁铁几何结构。
先由孔径大小50mm确定四极铁中心内接圆的大小,为保证一定精度,取内接圆半径为100mm。
双曲线断开的位置即B点坐标由经验公式x=(1.3~1.5)a来确定,双曲线断开后设计为一段切线段BC加一段平行线段CD。
其中平行线段CD可减少双曲线两端的顶端效应,类似于二极磁铁的垫补。
切线段长设计为5mm,平行线段长度由经验取6mm。
同样,在磁铁内侧设计圆角以避免磁场饱和。
初步设计好磁铁几何尺寸后,写入输入文件进行计算,由于四极磁铁的对称性,所以仅取八分之一部分进行计算即可。
根据程序计算结果微调几何结构,经过多次调整切线与平行线的长度及安匝数,最终确定四极磁铁的几何尺寸如图
软件的计算结果如图
调出四种四极磁铁磁感应强度梯度图线如图
以上各图中左端当接近于x=0时,出现数据发散的情况,这是因为程序在对磁感应强度求梯度时,当分母接近于0时误差被放大所致。
实际情况下,越靠近图像的左端,梯度是越均匀的。
由以上各图,可见在0~25mm内,磁感应强度梯度的均匀度达到了万分之一标准。
为减少设计成本,四个四极磁铁采用同一几何结构,调整线圈电流改变安匝数即可分别达到各个磁铁对磁感应强度梯度大小的要求。
五、最终设计结果
二极磁铁
内环弯转磁铁
安匝数
30002A
导线电流
166.67A
内环弯转磁铁
安匝数
25310A
导线电流
140.61A
四极磁铁
内环聚焦磁铁
安匝数
10930A
导线电流
182.17A
内环散焦磁铁
安匝数
10608A
导线电流
176.80A
外环聚焦磁铁
安匝数
10535A
导线电流
175.58A
外环散焦磁铁
安匝数
10235A
导线电流
170.58A
参考文献
1、李泉凤.2002.电磁场数值计算与电磁铁设计[M].北京:
清华大学出版社,149-184.
2、裴元吉王相綦.2012.加速器物理学.合肥:
中国科学技术大学课程讲义,63-94.
六、附录程序代码
二极铁
®kprob=0,!
PoissonorPandiraproblem
mode=0,!
Somematerialshavevariablepermeability
xreg1=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,kreg6=190,
kmax=200,
yreg1=3,lreg1=20,
yreg2=11,lreg2=54,
yreg3=14,lreg3=74,
yreg4=35.5,lreg4=138,
lmax=140&
&pox=0,y=0&
&pox=87,y=0&
&pox=87,y=43&
&pox=0,y=43&
&pox=0,y=0&
®mat=1,cur=30002&
&pox=29,y=3.5&
&pox=29,y=11&
&pox=42,y=11&
&pox=42,y=3.5&
&pox=29,y=3.5&
®mat=1,cur=-30002&
&pox=66,y=3.5&
&pox=66,y=11&
&pox=79,y=11&
&pox=79,y=3.5&
&pox=66,y=3.5&
®mat=2,mtid=1,mshape=0&
&pox=5,y=0&
&pox=27,y=0&
&pox=27,y=7&
&pont=2,x0=33,y0=7,x=0,y=6&
&pox=38,y=13&
&pont=2,x0=38,y0=7,x=6,y=0&
&pox=44,y=4.5&
&pont=2,x0=46,y0=4.5,x=0,y=-2&
&pox=47,y=2.5&
&pox=47,y=2.45&
&pox=48,y=2.45&
&pox=48,y=2.5&
&pox=62,y=2.5&
&pox=62,y=2.45&
&pox=63,y=2.45&
&pox=63,y=2.5&
&pox=64,y=2.5&
&pont=2,x0=64,y0=4.5,x=2,y=0&
&pox=66,y=35&
&pox=5,y=35&
&pox=5,y=0&
&mtmtid=1
bgam=0.000000.0017513135
9000.001747079
9500.001741742
10000.001735498
10500.001728309
11000.00172014
11500.001710963
12000.001700753
12500.001689494
13000.001677174
13500.001663786
28000.001080694
28500.001068051
29000.001056142
29500.001044912
30000.001034309
30500.001024289
31000.001014809
31500.001005828
32000.000997312
32500.000989226
33000.000981539
33500.000974222
40000.000904952
45000.000856798
50000.000818493
55000.000788085
60000.000764202
65000.000745863
70000.000732376
75000.000723261
80000.000718209
85000.000717054
90000.000719758
95000.000726411
100000.000737231
105000.000752594
105780.0007562580
113190.0007951022
119400.0008375209
124510.0008834703
129120.0009293680
133130.0009764671
136540.0010253255
139350.0010764263
142160.0011254924
144470.0011767475
146180.0012313603
147890.0012846865
150200.0013315579
151310.0013879251
152520.0014423770
154320.0014912019
155940.0015389351
157050.0015918497
161800.0018542555
168400.0023752969
171500.0029154519
173600.0034566194
176200.0039729837
178300.0044863167
182000.0054945055
189500.0079176564
.010*******
202000.0148588410
206500.0193798450
.023*******
216000.0370370370
219000.0456621005
230000.0869565217
233860.1002810000
238500.1181630000
244080.1387420000
250790.1622460000
258850.1888580000
268540.2186950000
280190.2517840000&
四极铁
®kprob=0,;PoissonorPandiraproblem
mode=-1;Materialshavefixedpermeability
dx=.05,dy=.05,;Meshintervals
yminf=0,ymaxf=15;FixedYforfieldinterpolation
xminf=0,xmaxf=15;Xrangeforfieldinterpolation
;;Thenext6termsrefertotheharmonicanalysis:
ktype=4,;Quadrupolesymmetry
nterm=10,;Numberofcoefficients
nptc=10,;Numberofarcpointsforinterpolation
rint=7.0,;Radiusofthearcforinterpolation
angle=45,;Angularextentofarc(defaultstart=0)
rnorm=1.0&;Apertureradiusfornormalization
&pox=0.0,y=0.0&
&pox=34,y=0.0&
&pox=34,y=34&
&pox=0.0,y=0.0&
®mat=3,mtid=3&
&pox=3.5355,y=3.5355&
&pont=3,r=5,x=7,y=1.7857&
&pox=8.049,y=1.627&
&pox=10.355,y=1.627&
&pox=19.137,y=5.146&
&pont=2,x0=20.253,y0=2.362,x=3,y=0&
&pox=23.253,y=0.0&
&pox=32.253,y=0&
&pox=32.253,y=12.362&
&pox=22.375,y=22.375&
&pox=3.5355,y=3.5355&
®mat=1,cur=10930&
&pox=15.819,y=3.6&
&pox=21.819,y=3.6&
&pox=21.819,y=0&
&pox=15.819,y=0&
&pox=15.819,y=3.6&
®ibound=0&
&pox=34,y=34&
&pox=3.5355,y=3.5355&
&pox=0.0,y=0.0&
&mtmtid=3
bgam=0.00000E+000.0017513135!
StartofB,Gammadata
0.11420E+040.0017513135
0.29530E+040.0010159504
0.51140E+040.0007821666
0.84760E+040.0007078644
0.96670E+040.0007241130
0.10578E+050.0007562580
0.11319E+050.0007951022
0.11940E+050.0008375209
0.12451E+050.0008834703
0.12912E+050.0009293680
0.13313E+050.0009764671
0.13654E+050.0010253255
0.13935E+050.0010764263
0.14216E+050.0011254924
0.14447E+050.0011767475
0.14618E+050.0012313603
0.14789E+050.0012846865
0.15020E+050.0013315579
0.15131E+050.0013879251
0.15252E+050.0014423770
0.15432E+050.0014912019
0.15594E+050.0015389351
0.15705E+050.0015918497
0.16180E+050.0018542555
0.16840E+050.0023752969
0.17150E+050.0029154519
0.17360E+050.0034566194
0.17620E+050.0039729837
0.17830E+050.0044863167
0.18200E+050.0054945055
0.18950E+050.0079176564
0.19500E+.010*******
0.20200E+050.0148588410
0.20650E+050.0193798450
0.20950E+.023*******
0.21600E+050.0370370370
0.21900E+050.0456621005
0.23000E+050.0869565217
0.23386E+050.1002810000
0.23850E+050.1181630000
0.24408E+050.1387420000
0.25079E+050.1622460000
0.25885E+050.1888580000
0.26854E+050.2186950000
0.28019E+050.2517840000&