Comsol Multiphysics 在脉冲强磁场中的应用

Radius(cm)
套上不锈钢之后的应力状态
放入液氮中冷却之后的应力状态
21
B. 放电过程中和放电结束时的应力状态
4
Stress(GPa)
3 2 1 0 -1
r hoop z von Mises
A
1.0 0.8 0.6 0.4 0.2 0.0 -0.2 -0.4 1 2 3 4 5 6 7 8 Radius(cm)
⑥
放电过程结束之后，因为磁体导体材料进入了塑性变形，而加固材料还 在弹性阶段，有一个回弹的过程，因此还有残余应力的存在。 v 脉冲磁体多次放电，从而疲劳寿命问题也是我们需要考虑的另一个方面。
18
弯曲应变方程
εH0 =
r−ρ
ρ
ρ
tension O compression
r
19
A. 磁体绕制过程
250 225 200 175 150 125 100 75 0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07
R a dius (m )
17
三、脉冲强磁体应力分析
v 一个脉冲磁体一次实验过程中的应力状态有： ① 在磁体绕制过程中，磁体导体以及加固材料要承受因弯曲变形而导致的 应力和一定的预应力； ② ③ ④ ⑤ 通常还会在外加固层套上一层加热的不锈钢，然后冷却到室温； 最后的组装过程中，还会给轴向方向紧上螺栓，从而施加了一定的压力； 磁体做好之后，准备开始实验前，要放入液氮中冷却到77K，这样就又 有一个冷却过程； 冷却充分之后，开始进入实验阶段，由于磁体处于放电过程中，就有电 磁力的作用，加上磁体导体材料发热及传导，从而在磁体导体和加固层 中就会导致热应力的产生；
T im e (m s )
Time (ms)
整个磁体模型中的电流密度动态演示
中平面处电流密度分布
5.00E+009 4.00E+009 3.00E+009
1 ms 5 ms 20 ms
Current density (A/m )
2
2.00E+009 1.00E+009 0.00E+000 -1.00E+009 -2.00E+009 -3.00E+009 0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07
300
T e m pe ra ture
250
200
150
100
50
13
不同近似情况下的磁场波形
90 80 70 60
E ddy c urre nt+ Ma g ne tore s is ta nc e Ma g ne tore s is ta nc e only N o e ddy c urre nt a nd no m a g ne tore s is ta nc e
σ
∂Aφ ∂t
∂T ρC p + ∇ ⋅ ( −k∇T ) = Q ∂t 1 1 σ ( B, T ) = = 1.1 ρ ( B, T ) ⎛ ⎡ ρ ( 273) ⎤ ⎞ −3 ρ (T ) ⎜1 + 10 ⎢ B ⎥ ⎟ ⎜ ⎣ ρ (T ) ⎦ ⎟ ⎝ ⎠
+ ∇ × ( µ0 −1µr −1∇ × Aφ ) = J φe
current (A)
-1.5x10 0
4 .8 m s 1 0 m s 2 0 m s 3 5 m s
5
30
35
5
me (ms)
Time (ms)
T e m pe ra ture (T )
r hoop
Stress(GPa)
Stress(GPa)
z von Mises
z von Mises
1
2
3
4
5
6
1
2
Radius(cm)
3 4 Radius(cm)
5
6
7
磁体弯曲之后的应力状态
磁体绕制之后的应力状态
20
0.4 0.3 0.2 0.1 0.0 -0.1 -0.2 -0.3 -0.4 0.3 0.2 0.1 0.0 -0.1 -0.2 -0.3
COMSOL Multiphysics 在脉 冲强磁场中的应用
宋运兴 华中科技大学武汉强磁场中心
1
提纲 一、脉冲强磁场概述 二、脉冲强磁体放电过程分析 三、脉冲强磁体应力分析 四、脉冲强磁体快速冷却分析
2
一、脉冲强磁场概述
80
强磁场产生 强大的电动力
Magnetic Field (T)
PMT50S PMT60S PMT70S PMT80S
10
磁通渗透过程
1.0 0.8
normalized Field
0 ms 0.35 ms 1 ms 4.8 ms
0.6 0.4 0.2 0.0 0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07
Radius (m)
11
中心场磁场波形和dB/dt波形
40 30 20 10 0 -10 0 5 10 15 20 25 30 80 70
60
40
20
0 0 20
Time (ms)
40
60
制约场强提高的主要因素：大电流和强磁场的相互作用产生强大的电动力； 焦耳热和众多电磁效应(磁致电阻效应和涡流效应)在磁体内引起的温升。 如何发挥现有材料性能设计出尽可能高参数的磁场是个极具挑战的课题!
3
短脉冲磁体
单线圈 ＋ 单电容器电源
电触头
电路设计图
Time (ms)
400 350 300
Current in ss (A)
Time (ms)
Temperature (K)
Tss0 Tss2 Tss5 Tss10 Tss20
250 200 150 100 50 0.00
0.02
0.04
0.06
0.08
0.10
16
不锈钢使用在内加固层
90 80 70 60
0.4 0.3 0.2 0.1 0.0 -0.1 -0.2 -0.3 -0.4 0.3 0.2 0.1 0.0 -0.1 -0.2 Stress(GPa) Stress(GPa)
r hoop
0.4 0.3 0.2 0.1 0.0 -0.1 -0.2 -0.3 -0.4 0.3 0.2 0.1 0.0 -0.1 -0.2
电阻和电感随时间的变化
100 1 .2 90 1 .1 80
dB/dt (T/ms)
R es is ta nc e (m Ω)
Field (T)
50 40 30 20 10 0 35
70
Induc ta nc e
60
1 .0
0 .9 50
40
0 .8
30 0 5 10 15 20 25 30 35 0 .7
8
二、脉冲强磁体放电过程分析
忽略不对称因素，认为脉冲磁体是轴对称结构
9
脉冲磁场放电过程涉及到电场、磁场以及温度场， 是一个复杂的多物理场耦合过程，而COMSOL Multiphysics具备耦合场分析上的优势，这是我选择 它的主要原因。
电路方程 电磁方程 热传导方程 电导率函数
diL ( Rex + Rcoil )iL + ( Ls _ in + Ls _ out + LM + Lex ) × − uC = 0 dt
60
Pick-up coil voltage (V)
12 10 8 6 4 2 0 -2 0
40
10 kv
30 20 10 0 20
pick-up coil voltage
5
10
15
Time (ms)
Magnetic field (T)
50
7
PMDS（ Pulsed Magnet Design Software ）
diSS ( j ) diL n − uC = 0 ( Re x + R ) iL + L + ∑ M1, j dt j =2 dt
15
不锈钢使用在外加固层
Current in conductor (A)
4x10 3x10 2x10 1x10 0
4 4 4 4
80
60
Field (T)
B0 Bss2 Bss5 Bss10 Bss20
环氧树脂 不锈钢筒
磁场波形
Cu导线 Glass纤维
80 PM T50S PM T60S PM T70S PM T80S
Zylon纤维 CuAg导线 Zylon纤维
60
( T)
连接头
40 20 0 0
20
40
60
(m s)
4
长脉冲磁体
双线圈 ＋ 脉冲发电机电源
10 kV 母线
6.9 kV 母线
内线圈
Stress(GPa)
r hoop z von Mises
B
9
10
(A) 磁场峰值时候的应力状态 (B) 放电结束后的残余应力状态
22
A 脉冲结束后的等效塑性应变分布 (A) 考虑弯曲应变 (B) 不考虑弯曲应变
0.05
0.06
0.07
Radius (m)
400
350
tra ns port c urre nt,m a g ne tore s is ta nc e a nd e ddy c urre nt tra ns port c urre nt a nd m a g ne tore s is ta nc e tra ns port c urre nt only