一种利用形状函数估算进行一次性磁场垫补的方法技术

本发明专利技术提出了一种利用形状函数估算进行一次性磁场垫补的方法，包括如下步骤：实测各个半径处的平均磁场，该各个半径处的平均磁场与理论的等时性平均磁场之间的差为

【技术实现步骤摘要】
一种利用形状函数估算进行一次性磁场垫补的方法


[0001]本专利技术属于回旋加速器磁场垫补
，尤其涉及一种利用形状函数估算进行一次性磁场垫补的方法
。

技术介绍

[0002]磁场垫补是等时性回旋加速器中必不可少的环节
。
传统方法的磁场垫补通常需要垫补3次以上，每次切削镶条后由于角宽度变化了
(
所述角宽度为含有镶条的磁极的角宽度，有可能是磁极的一侧含有镶条，也有可能是磁极的两侧含有镶条
)
都需要重新计算形状函数，如此，完成整个垫补过程至少需要计算3次以上的形状函数计算
。
所述需要计算3次以上的形状函数的原因是，镶条的垫补决不允许“过垫补”而必须采用“逐渐逼近”的方法，这是由于切削垫补过程中，只能切削镶条而不能增加镶条厚度
。
如果只按照初始形状函数
dB
i
进行一步到位的垫补量计算，当垫补量较大时，有可能出现“过垫补”，即镶条切削过量，造成垫补失败
。
所以只能采取每次垫补量都比较小
、
重复形状函数的计算
、
重复将镶条到加工厂切削等步骤
、
以保证不会出现“过垫补”的现象
。
然而多次垫补就需要多次计算形状函数
。
形状函数的形成如图3所示，其横坐标为加速器半径，纵坐标为磁场，图3中的初始形状函数
dB
i
是由多个带有具有独立波峰的形状函数叠加而成的，据统计，用有限元软件计算图2其中一个独立波峰的形状曲线就需要2个多小时，如果将图2所有具有独立波峰的形状函数计算完毕就需要将近2天的时间，并且这只是一次计算形状函数的时间，如果是3次以上的垫补就需要几天的计算形状函数的时间，导致按照传统的垫补方法需要高的计算资源与较长的磁场垫补时间周期
。

技术实现思路

[0003]本专利技术为解决现有技术存在的问题，提出一种利用形状函数估算进行一次性磁场垫补的方法，目的在于解决按照传统磁场垫补方法需要较长的磁场垫补时间周期的问题
。
[0004]本专利技术为解决其技术问题提出以下技术方案：
[0005]一种利用形状函数估算进行一次性磁场垫补的方法，包括如下步骤：
[0006]步骤一
、
实测各个半径处的平均磁场，该各个半径处的平均磁场与理论的等时性平均磁场之间的差为
B
shim
；若在
M
个半径处进行了平均磁场的测量，则
B
shim
为
M
维的向量，
B
shim
＝
[B
shim1
,B
shim2
,B
shim3
,
…
B
shimM][0007]步骤二
、
未切削镶条时，经过有限元软件模拟计算得到初始形状函数
dBi
；
[0008]步骤三
、
将各个半径处的磁场差值
B
shim
分为
N
份，进行
N
次计算，每次一计算预计磁场变化量为
B
shim
/N
＝
[
B
shim1
/N,B
shim2
/N,B
shim3
/N,
…
B
shimM
/N
]；
[0009]步骤四
、
第
k
次计算：预计磁场变化量为
B
shim
/N
；
k
取值范围为0到
N
‑1；
[0010]步骤五
、
建立基于估算函数的形状函数
dB
k
；
[0011]步骤六
、
通过垫补算法得到第
k
次计算的垫补量
X
k
，
X
k
＝
[X
k2
,X
k4
,
…
X
k2n
],n
为切除三角形计算磁场变化的次数；为第
k
次计算时第
2j
个折点处的垫补量，
j
取值范围为1到
n。
[0012]步骤七
、
完成上述
N
次计算后，第
2j
个折点处总的垫补量为
N
次垫补量的和即：以上共计算了
R2，
R4，
…
，
R
2n
处，共
n
个折点处的垫补量，其他
n+1
个折点垫补量通过线性插值得到
。
[0013]其特征在于：所述基于估算函数的形状函数
i
＝
0。
[0014]进一步地，所述步骤四的第
k
次计算，当
k
大于0时，其他各次计算均使用根据
dB
i
所估算的形状函数
dB
k
。
[0015]进一步地，所述建立基于预估算法的形状函数
dB
k
，既是利用形状函数
dBi
获得预估镶条不同厚度时的形状函数其中，
dBk
为除了基础形状函数
dB
i
以外的第
k
个形状函数，具体过程如下：
[0016]1)
计算镶条初始角宽度为
A
i
(R)
时的相对角宽度变化量其中，
A
i
(R)
为镶条初始角宽度向量，
A
i
(R)
＝
[A
i
(R2),A
i
(R4),
…
,A
i
(R
2n
)],
[0017]n
为形状函数切割三角形的个数；为与半径为
R
处镶条当前边缘曲线相距
h
处的角宽度向量；
f
i
(R)
为相对角宽度变化量向量；
f
i
(R)
＝
[f
i
(R2),f
i
(R4),
…
,f
i
(R
2n
)]；
[0018]2)
第
k
次计算，镶条角宽度变为
A
k
(R)
时的相对角宽度变化量：
[0019]其中，
A
k
本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.
一种利用形状函数估算进行一次性磁场垫补的方法，包括如下步骤：步骤一
、
实测各个半径处的平均磁场，该各个半径处的平均磁场与理论的等时性平均磁场之间的差为
B
shim
；若在
M
个半径处进行了平均磁场的测量，则
B
shim
为
M
维的向量，
B
shim
＝
[B
shim1
,B
shim2
,B
shim3
,
…
B
shimM
]
步骤二
、
未切削镶条时，经过有限元软件模拟计算得到初始形状函数
dBi
；步骤三
、
将各个半径处的磁场差值
B
shim
分为
N
份，进行
N
次计算，每次一计算预计磁场变化量为
B
shim
/N
＝
[B
shim1
/N,B
shim2
/N,B
shim3
/N,
…
B
shimM
/N]
；步骤四
、
第
k
次计算：预计磁场变化量为
B
shim
/N
；
k
取值范围为0到
N
‑1；步骤五
、
建立基于估算函数的形状函数
dB
k
；步骤六
、
通过垫补算法得到第
k
次计算的垫补量
X
k
，
X
k
＝
[X
k2
,X
k4
,
…
X
k2n
],n
为切除三角形计算磁场变化的次数；为第
k
次计算时第
2j
个折点处的垫补量，
j
取值范围为1到
n
；步骤七
、
完成上述
N
次计算后，第
2j
个折点处总的垫补量为
N
次垫补量的和即：以上共计算了
R2，
R4，
…
，
R
2n
处，共
n
个折点处的垫补量，其他
n+1
个折点垫补量通过线性插值得到；其特征在于：所述基于估算函数的形状函数
2.
根据权利要求1所述一种利用形状函数估算进行一次性磁场垫补的方法，其特征在于：所述步骤四的第
k
次计算，当
k
大于0时，其他各次计算均使用根据
...

【专利技术属性】
技术研发人员：边天剑，安世忠，管锋平，魏素敏，邢建升，王飞，陈忻禹，陆锦荣，凌丽，
申请(专利权)人：中国原子能科学研究院，
类型：发明
国别省市：