1. 引言
质子是带1个正电荷的粒子,是原子核的组成部分。用于医学的质子来源于氢(H2),氢电离后成为质子(H+),经同步或回旋加速器加速到近光速后应用于治疗各种疾病。质子束的最大特征是它进入人体内形成的尖锐的Bragg峰,在形成峰之前的平坦段称为坪(Plateau),峰后则是一个突然减弱陡直的尾。由于Bragg峰太尖,所以可将它扩展后形成与肿瘤大小吻合的扩展Bragg峰(Spread Out Bragg Peak, SOBP) [1]。对于小肿瘤,可以根据它在体内的深度选择质子能量,使Bragg峰落在肿瘤位置上;对于较大的肿瘤,可以通过调整质子的能量使Bragg峰展宽到与肿瘤厚度相当。这样质子几乎不杀伤肿瘤后面和侧面的正常组织,对肿瘤前面的正常组织损伤也较小。作为癌症最有效的治疗手段之一,质子治疗以其独特的对人体正常组织的保护性和对病灶肿瘤细胞杀伤的精确性成为放射疗法的主要研究趋势,质子加速器是质子治疗系统中的核心装置。质子治疗加速器的能量要求在70~250 MeV之间,相应于穿透软组织或水3~40 cm的深度,束流强度则要求在nA量级。由于质子射程取决于质子的能量,因此质子治疗过程中要求质子的能量能够连续调节。
目前世界上的质子治疗仪加速器分为三种:直线加速器、同步加速器和回旋加速器,由于回旋加速器其占地小、部件少、运行便捷、束流稳定的优点广为应用,质子治疗回旋加速器又分为紧凑型回旋加速器和分离扇型回旋加速器。一台完整的回旋加速器主要由磁铁系统、高频系统、离子源、束流引出系统、束流诊断系统、靶系统、真空系统、冷却系统构成 [2],其中磁铁系统是回旋加速器中最主要的部件之一,他的作用是提供带电粒子在所控制的轨道中做圆周运动所需要的磁场强度,由磁铁、线圈、电源等设备组成。
本文主要基于磁路原理对230 Mev医用紧凑型回旋加速器的主磁铁系统进行设计。
2. 基本电磁参数的选择
磁铁形状选型:回旋加速器一般选用图示的方框型,粒子的运动平面是水平的,结构简单,真空室暴露出来的部分较少,粒子源及粒子引出设备安装调试方便。
Figure 1. Schematic diagram of structure of box magnet
图1. 方框型磁铁结构示意图
根据其竖直水平对称特性,取四分之一的模型结构为研究对象,图中各个参数及意义:R1为磁极半径,H1为气隙高度,R2为线圈外径距中心轴的距离,R3为磁轭内半径,R4为磁轭外半径,H2为线圈到中平面的距离,H3为磁极上部到中平面的距离,H4为上盖板到中平面的距离,R5和H5分别为除去导角的盖板半径和磁轭高度。
磁铁材料选择:对于回旋加速器的材料要求饱和磁感应强度高、矫顽力小、磁时效性小、磁导率高 [3],这里选用10号钢(C = 0.07%)。
气隙高度及磁峰区初始角的选择:参考IBA公司235 Mev质子回旋加速器的参数 [4],气隙高度H1选为0.08 m,初始峰区角α = 54˚。
磁场中平面平均磁感应强度B、磁密峰值BH及磁密谷值BV的选择:在磁刚度BR一定的情况下,
,而R越大,体积V越大,因此,要使磁铁小型化,应选择尽可能大的B。而回旋加速器的磁场在粒子加速过程中提供一个相当大的洛仑兹力,这个力有可能将粒子(一般为H−)的电子剥离 [5]。为了在较高能量时将束流损失限制在合理范围内,B应低于1.3 T,且如果工作在近饱和区,要耗费很大的励磁功率,这里我们初步选择B为1.0 T [6]。又根据B、BH、BV之间的关系:
(1)
可以初步选取峰区磁密
,谷区磁密
。
电流密度J的选择:综合考虑能耗成本、冷却方式等各种因素,电流密度J初步取为1.5 A/mm2 [7]。
磁极利用系数ξ取为0.86、线圈窗口的高和宽的比值λ取为1.5、导体有效截面与线圈窗口总截面积的比值η取为0.67 [7]。
3. 主磁铁结构尺寸及其他电磁参数的计算
1) 磁钢度BR计算:由相对论原理、质能方程和质子在加速器中的运动方程式 [7]:
(2)
(3)
(4)
可以计算出磁钢度:
(5)
式中W为质子的最终动能,即230 Mev,q为质子所带电荷(H−),即1.602117 × 10−19库仑,m0为质子的静止质量,值为1.672648 × 10−27。
2) 电压参考频率f计算:加速器主磁铁在一个圆周内采用4块磁极,即质子加速区也有四个,则加速电压得参考频率为质子回旋频率的4倍,即:
(6)
3) 根据磁刚度可以计算出磁极半径R1:
(7)
4) 中心磁感应强度B0的计算:根据等时性要求和相对论下质子的基本运动方程可以得到:
(8)
5) 根据中心磁感应强度B0可以求出离散半径r(i)上的平均磁密B(i):
(9)
6) 将(9)式带入(1)式可以求出每一半径r(i)处得磁极角度
:
(10)
7) 其他结构参数的计算:
H2由经验公式求得:
(11)
假定磁轭内半径与线圈外半径的距离之差为R23,R23的值由经验公式可得:
(12)
H3的值可参照图1由R3表示出来:
(13)
R4和H4均由经验公式得到 [7]:
(14)
(15)
根据全电流定律:
(16)
根据图1分析,该磁路分四段,即l1,l2,l3和气隙段,l1段的磁感应强度为峰区磁密BH,长度为
;l2段的磁感应强度由经验公式取为
,长度为
;l3段的磁感应强度按经验公式取为0.9 BH,长度为
;气隙段的磁感应强度为峰区磁密BH,长度为H1/2;根据10号钢的磁化曲线可以求出各段磁路的磁场强度。
励磁安匝数的表达式为 [8]:
(17)
由式(17)和(16)可以解出:
由R3的值进而求得其他值。
4. 根据设计参数建立磁铁模型
设计参数计算结果(表1):
Table 1. Structural parameters of the main magnet model
表1. 主磁铁模型结构参数
在AnsysMaxwell中建立主磁铁的2D有限元模型(1/2)如图2所示 [9]。图中,区域A1为磁极峰区(若为谷区模型时则为空气),A2盖板,A3为磁轭,A4为线圈,A5为空气,A6为气隙。
Figure 2. 2D model of main magnet (peak)
图2. 主磁铁2D有限元模型(峰区)
对模型添加Balloon边界条件,对线圈区域加载激励条件,初始电流密度取1.5 × 106 A/m2。主磁铁峰区和谷区模型的磁力线分布如图3(a)和图3(b)所示:
(a) (b)
Figure 3. (a) Magnetic field line distribution in peak region; (b) Magnetic field line distribution in valley region
图3. (a) 峰区磁力线分布;(b) 谷区磁力线分布
峰区和谷区中平面沿径向的磁场分布如图4(a)和图4(b)所示:
(a)(b)
Figure 4. (a) Mid-planar magnetic field distribution in peak region; (b) Mid-planar magnetic field distribution in valley region
图4. (a) 峰区中平面磁场分布;(b) 谷区中平面磁场分布
由图2和图3可以看出,主磁铁峰区中平面上的平均磁场为1.72 T左右,谷区中平面的平均磁场为0.2T左右,与设计值相差不大。
5. 结论
1) 本文通过磁路原理设计了230 Mev质子治疗仪回旋加速器主磁铁的主要结构和电磁参数,为后续对磁场的优化设计和电动力学计算奠定了基础。
2) 应用有限元分析软件ANSYSMaxwell对设计的主磁铁建立了2D的有限元模型。
3) 通过ANSYSMaxwell软件对设计的主磁铁进行了电磁场计算,计算结果表明:在磁铁峰区中平面磁通密度能达到1.72 T左右,在谷区中平面上磁通密度能达到0.2 T左右,均满足设计要求。