第5期2021年5月
机械设计与制造
Machinery Design&Manufacture217
常态运行的ITER杜瓦结构分析
汪献伟谢飞2，金环3，何庆1
(1.江苏理工学院机械工程学院，江苏常州213001；.顺德职业技术学院电子信息工程学院，广东佛山528300；
3.中国科学院等离子体物理研究所，安徽合肥230031)
摘要：作为一个封装托卡马克的真空压力容器，杜瓦将为冷屏和磁体系统提供真空环境，同时传递来自于自身和线圈系统的载荷，在正常和非正常运行的过程中杜瓦将会遭受地震载荷、电磁力、测试压力、自身净重等载荷的冲击，其中净重和压力载荷是影响杜瓦结构安全性的关键因素，该研究首先基于无力矩理论计算了杜瓦在压力载荷、自身净重以及两者联合作用下的应力状态，获得了杜瓦的危险区域，在此基础上通过应力线性化，提取了危险区域的薄膜、弯曲和峰值应力，最后基于ASME的分析设计对各种不同的应
力进行了评定，并分析了局部应力集中的原因，这对于杜瓦的设计优化以及中国聚变工程实验堆杜瓦的设计制造都具有一定的指导作用。
关键词：托卡马克；顶盖；真空容器；杜瓦；应力线性化
中图分类号:TH16；TL351.1文献标识码:A文章编号：1001-3997(2021)05-0217-04
The Structural Analyses of ITER Cryostat During Normal Operation
WANG Xian-wei1,XIE Fei2,JIN Huan3,HE Qing1
(l.School of Mechanical Engineering,Jiangsu University of Technology,Jiangsu Changzhou213001,China;2.School of Electronics and Information Engineering,Shunde Polytechnic,Guangdong Foshan528300,China;3.Institute of Plasma Physics,Chinese Academy of Sciences,Anhui Hefei230031,China)
Abstract：As a vacuum vessel enveloping the Tokamak machine,cryostat will provide the vacuum environment for thermal shield and magnet system and also transfer the loads derived from itself and magnetic coils.During the normal and off-normal operation,cryostat may encounter the impact of seis mic load,electromagnetic force,t est pressure,atmospheric pressure and its dead weight.A nd dead wei
ght and pressure are the key factors that impact the structural safety of cryostat.Firstly,based on the no torque theory,we calculate the cryostat stress state under the dead weight,external pressure and their combination, then we obtain the dangerous region.On this basis the membrane stress,bending stress and peak stress in the dangerous area are extracted by stress linearization.In the end,stress assessment of cryostat is carried out based on the design-by-analysis method specified in ASME.Meanwhile,the reasons of local stress intensity are analysed.The study will provide somewhat guidance for the cryostat design optimization as well as for the design and fabrication of cryostat of China Fusion Engineering Test Reactor.
Key Words:Tokamak；Top Lid；Vacuum Vessel；Cryostat；Stress Linearization
1引言
杜瓦是国际热核实验堆的核心部件之一，它是一个内压为10-^Pa的抽真空密封容器。杜瓦包裹托卡马克装置主机，从而为冷屏和磁体系统提供真空运行环境。在杜瓦的周向开设有大量的窗口可以为磁体馈线、水冷管、检测引线等提供安装通道，此外杜瓦可以利用其支撑结构将托卡马克主机和杜瓦自身在正常运行以及特定的事故状态下产生的热载荷传递到地面，以减小热载荷产生的危害。整个杜瓦系统主要由四大部分组成：最上方的顶盖、最下方的平底支撑、以及中间区域的上环体、下环体，如图1所示。
由于杜瓦的结构尺寸庞大，杜瓦的安装按照自下而上的顺序进行垂直装配，装配完毕后的杜瓦高度约29m，直径约28m，壁厚约50mm|1-3]，杜瓦的主要参数，如表1所示。杜瓦的支撑系统主要包含三大部分：裙座支撑、基座环和径向托架，裙座作为杜瓦筒体的外延部分，通过地脚螺栓固定在杜瓦生物屏蔽层的地基上，该安装方式确保了裙座具有足够的强度来支撑杜瓦的重量，同时承受地震载荷的冲击。基座环由两组同心的圆柱形板壳组成，板壳
来稿日期:2020-07-24
基金项目：广东省自然科学基金(2017A030310640)
作者简介：汪献伟,(1984-)，男，洛阳人，博士研究生，讲师，主要研究方向：聚变堆磁体线圈的力学特性分析
218汪献伟等：常态运行的ITER杜瓦结构分析第5期
内部设置有径向的加强筋，每相邻两组加强筋之间的周向间隔为4。，基座环的下方与18个沿周向布置的支撑柱相连接，基座环的主要作用是传递各种类型的载荷，诸如重力载荷、地震载荷、电磁载荷以及由于装置热变形产生的支反力载荷等，托卡马克真空室、环向场线圈、极向场线圈在运行状态下产生的各种载荷都可以通过基座环传递给重力支撑系统。径向托架沿圆周方向焊接在杜瓦的外部，每相邻两个
托架的周向间隔为20。，托架可以紧密的锁合在生物屏蔽层的托架槽中，从而抵抗地震载荷引起的结构晃动。的安全裕度。杜瓦作为抽真空压力容器，在正常运行状态下必须要承受外部的大气压力载荷，此外，杜瓦净重作为一种惯性载荷，贯穿杜瓦运行的始终。由于杜瓦的结构尺寸庞大、窗口众多，重力和压力载荷可能导致杜瓦的结构大变形和局部应力集中，它们是杜瓦的结构设计中必须考虑的载荷类型。
2.2边界条件reactor pressure vessel
杜瓦的边界条件施加取决于载荷状态和支撑方式，在冷却液泄露进入杜瓦的反常运行状态下，支撑裙座可以允许杜瓦的径向运动|7］。据此杜瓦裙座径向自由度被保留而其余的自由度均被抑制。对于基座环，其底部和杜瓦的支撑柱刚性连接，而支撑柱通过地脚螺栓固定在地基上，因此基座环既没有平动自由度也没有转动自由度。为了建模分析的方便，在有限元分析中可以忽略支撑柱模型，但必须在支撑柱和基座环的交界面施加相应的固定约束。对于支撑托架，在杜瓦反常冷却条件下其可以径向运动，从而缓解杜瓦由于局部热膨胀差异引起的热应力，但是支撑托架紧密的锁扣在生物屏蔽层的配合插槽中，从而限制杜瓦的圆周方向转动。20。杜瓦的计算模型和约束形式，杜瓦在压力和重力载荷下的力学状态分析将基于该约束模型，如图2、图3所示。
图1杜瓦的结构
Fig.1The Structure of Cryostat
表1杜瓦的参数
Tab.1Parameters of Cryostat
参数项数值
杜瓦的外径28.54m
杜瓦的高度29.25m 壁厚(40~180)mm 筒体厚度50&60mm 设计基准压力10-4Pa 表面积~3400m2顶盖重量656ton
上下筒体重量1409ton 总重~3500ton
图220。杜瓦计算模型
Fig.2The Calculation Model of20。Cryostat
2基于无力矩理论的杜瓦力学计算
2.1杜瓦的载荷状态
杜瓦在运行过程中要遭受各种不同载荷事件的冲击，主要包括：惯性载荷、等离子体电流破裂事件、等离子体垂直位移事件、磁体系统快速放电事件、真空度降低事件以及正常运行状态下的大气压力载荷。其中等离子体电流破裂、等离子体垂直位移、磁体系统快速放电将会诱发电磁力，对杜瓦产生电磁冲击，电磁分析表明，电磁载荷在杜瓦内部产生的最大屈斯卡应力大约5
MPa远远小于由其它载荷事件引起的内应力，因此在杜瓦的结构强度计算中可以忽略电磁载荷的影响1地震载荷作为一种惯性载荷，将会导致杜瓦的稳定性破坏，文献16的分析表明杜瓦的最大屈曲值远小于地震加速度，也就是说杜瓦不可能发生屈曲失效，此外，由地震载荷引起的最大米塞斯应力远小于杜瓦的设计应力阈值，因此地震载荷也不会导致杜瓦的结构破坏。冷却氦泄漏、空气或者冷却水进入杜瓦都会改变真空环境，导致杜瓦的真空度降低，为了避免冷却液泄露引起的杜瓦反常失效，在杜瓦试车之前需要进行压力测试实验，测试压力为设计压力的1.15倍，测试压力主要用来模拟验证杜瓦在真空度降低的情况下是否具有足够
Fig.3The Constraints of Cryostat
2.3压力和惯性载荷下的杜瓦结构分析
杜瓦采用全焊接圆柱型增强单层腔体，其主要材料为不锈钢304/304L,304/304L的材料特性，如表2所示。考虑到杜瓦结构尺寸庞大，开展详细的全模型分析不仅对硬件配置要求极高，而且耗时耗力。鉴于
杜瓦的近似轴对称性，可创建循环对称模型来开展结构分析。由于在杜瓦不同区域窗口的布置形式不尽相同，为了保证分析结果的可靠性，本分析选择杜瓦的最危险区域，创建20。计算模型，并在模型的高低边界施加循环对称约束以模拟杜瓦整模在外载荷下的力学状态。
根据厚径比杜瓦属于薄壁压
No.5
May.2021
机械设计与制造219
力容器的范畴，在有限元分析中杜瓦的腔体和加强筋全部采用壳
单元进行模拟。杜瓦的圆柱型筒体在测试压力或大气压力下的应 力状态可以进一步分解为轴向应力、周向应力和径向应力，基于
无力矩理论杜瓦筒体的周向应力最大，轴向应力次之，而径向应 力最小。由于开孔和窗口的存在，导致杜瓦的应力分布发生改变,
其中窗口区域附近的应力较大，如图4所示。杜瓦筒体在重力和
压力联合作用下的最大屈斯卡应力约83MPa ，远低于材料的设计
应力强度，最大变形约4.6mm，不会与其它内部部件发生干涉。
表2材料特性
Tab.2 The  Material  Property
材料特性值( 304/304L )单位最小拉伸强度517/483
MPa 最小屈服强度207/172MPa 弹性模量195000MPa 设计应力强度
137
MPa 泊松比
0.3-
剪切模量
75000MPa 密度7.78x10-5
N/mm 3
a)屈斯卡应力
(b)总变形
5
rwxML  scumoH
RTFF-I, sue TIFC-L
PawwCrapt-krE  ESKET-L
DMX  J44.77M 5MM  -5HK  M4.OT4
WOini
H  9W7
(c )屈斯卡应力 (d)总变形
图4净重和压力载荷下的应力和变形
Fig.4 Stress  and  Deformation  Under  the  Dead  Weight  and  Pressure  Load
在净重和大气压力联合作用下的应力和变形，如图4(a)、图
4(b)所示。在净重，大气压力以及测试压力联合作用下的应力和
变形，如图4(c)、图4(d)所示。
杜瓦的顶盖采用碟形封头，由三部分组成：球形区域、筒体 区域和过渡区域。球型区域处于两向应力状态，且轴向应力与周
向应力相同，筒体区域的周向应力为轴向应力的2倍，而过渡区 域的轴向和周向应力分别，如式(1)、式(2 )所示罔。实际上考虑到
杜瓦顶盖在压力载荷和重力载荷下可能岀现的塑性破坏，在杜瓦 顶盖上布置了一系列的径向和环向加强筋以增强顶盖的稳定性，
加强筋的存在将明显改变顶盖的应力应变状态，分析结果表明最
大应力岀现过渡区壁面和径向加强筋的交界点，最大变形岀现在
椭圆区域的中心位置，总变形约44mm 。实际上顶盖的环向加强
筋1#,#,#之间存在间隔分布的径向加强筋，该加强筋将有助
于减小顶盖椭圆区域的总变形，如图 5 所示。
(a)有限元模型
(b)完整实体模型
图5顶盖结构
Fig.5 Structure  of  Top  Lid
r  D  _ i
滓亠 +2-r 1 (1)
渍2t  [ r 1+石n^
蓡
滓电L +
灯s 蓡
⑵
1 sin 渍
」[ r 1 sm 渍
式中:P —介质压力；s —碟形封头壁厚；；一过渡圆弧之半径；;—
碟形封头之平均半径；一所求应力点第二曲率半径与回转
轴的夹角。
3应力线性化及结构评估
3.1顶盖的应力线性化
顶盖由于曲率不连续、厚度不一致，在外部压力载荷的作用
下除了产生薄膜应力还存在较大的弯曲应力和峰值应力。为了详
细分析厚度方向上各种不同类型应力的分布状态，可以通过应力 线性化提取弯曲应力和峰值应力。对于任意壁厚为S ，中面曲率
半径为R 的筒体，假定其实际应力、平均应力、等效弯曲应力和
峰值应力分别为s,s ”,sb,s p 。根据静力等效、静弯矩等效以及弯曲
应力沿截面的合力为零的原则，可以得到沿壁厚方向各种不同类 ( 3)
⑷
( 5)
b 型应力的分布状态|9]:
s m
=
(
I 1 +kI 2)/
S
s b  =] _ S  + ^^ 蔀 f  I 2+-------
(12厶一
S  I 1
)
b  \
s  S  丿丨
12-k  S s  =s -s m  -s (h  h
h  2
滓dt ,2
= |滓tdt ,3
二|滓t  dt ,,—中面的曲率和通过 -h
丿-h
丿
-
h
中面的法向坐标。
由于应力主要沿厚度方向变化，分析采用plane  182单元来
模拟顶盖的横截面。 由于顶盖几何结构的轴对称性，可只创建横
截面的1/2模型，应力线性化路径选择曲率不连续的过渡区域和 厚度不连续的交界区域，如图6所示。
图6联合载荷下的顶盖应力线性化路径
Fig.6 Stress  Linearized  Paths  of  Top  Lid  Under  Combination
Load
No.5
May.2021 220机械设计与制造
应力线性化的结果，如图7、图8所示。在靠近内壁一侧，路径1上的峰值应力很小，而弯曲应力远大于
总体薄膜应力，是导致顶盖岀现局部大应力的关键因素，这与文献1的分析结果基本一致。在路径2上弯曲应力较小，而在靠近筒体外壁一侧，各径上的峰值应力急剧增升高，这是由该区域筒体厚度的阶跃型变化所导致的。
5.01315.03925.06538.09145.117
路径1长度方向
图7路径1的应力线性化
Fig.7Stress Linearization Along Path1
060.004120.008180.012240.016
30.00290.006150.01210.014270.018
路径2长度方向
图8路径2的应力线性化
Fig.8Stress Linearization Along Path2
3.2杜瓦的应力评估
杜瓦的应力评估基于ASME的分析设计方法|11]，在分析设计中应力被进一步划分为总体薄膜应力、局部薄膜应力、弯曲应力、二次应力、峰值应力等不同类型，并针对不同的应力类型采用不同的评定标准，如表3所示。杜瓦在净重和压力载荷下产生的应力主要由总体薄膜应力、局部薄膜应力、弯曲应力以及峰值应力组成。不锈钢的设计应力强度为172MPa，由图4可知杜瓦的总体薄膜应力和局部薄膜应力均小于各自的应力强度阈值，不会导致结构的总体破坏，但是由图7可知顶盖的过渡区域存在较大的弯曲应力是杜瓦顶盖设计需要重点关注的区域，也是杜瓦结构优化设计的重点位置。
表3各种应力的限制
Tab.3The Limitation of Different Kinds Stress
应力类型符号阈值
一次总体薄膜应力滓
m
s”
一次局部薄膜应力滓
L
1.5S m
一次薄膜加弯曲应力滓
L+滓b
1.5S m
一次加二次应力滓
l+滓b+Q3S m
总应力滓
l+滓b+Q+滓p S a 4结论
采用有限元分析软件创建20。杜瓦保守计算模型计算了杜瓦在净重、压力载荷联合作用下的杜瓦力学状
态，分析结果表明：无论是在净重、大气压力载荷作用下，还是在净重、测试压力、大气压力联合作用下，杜瓦的最大屈斯卡应力总是岀现在顶盖的过渡区域，杜瓦的最大变形总是岀现在顶盖球型区域的中心位置；在净重和压力载荷作用下，杜瓦的总体薄膜应力和局部薄膜应力均未超岀各自的设计应力阈值，不会导致结构的整体破坏；杜瓦顶盖的过渡区域由于曲率不连续在压力载荷的作用下存在较大的弯曲应力，杜瓦顶盖与筒体的交界区域由于厚度突变存在较大的峰值应力，为了保证杜瓦的安全运行，以上区域是杜瓦结构优化设计需要重点考虑的区域。
参考文献
[1]Doshi Bharat,Zhou Caipin,Ioki Kimihiro.ITER cryostat-an overview and
design progress J]Fusion Engineering and Design,2011,86(9—11):1924 -1927.
[2]Sekachev Igor,Meekins Michael,Sborchia Carlo.The cryostat and subsys­
tems development at ITER[C].//25th International Cryogenic Engineer­ing Conference and International Cryogenic Materials Conference.Ensch—ede：Elsevier Science fBV,2015:294-301.
[3]Doshi.B,Xie Han,Zhou Caipin.Design and manufacture of the ITER cr—
yostat[C]./2013IEEE25th Symposium on Fusion Engineering^San Fra—ncisco：IEEE，2013：13846204.
[4]Alekseev A.,Arslanova,D.,Belov,  A..Computational models for electro­
magnetic transients in ITER vacuum vessel,cryostat and thermal shield [J[.Fusion Engineering and Design,2013,88(9-10)：1904—1907.
[5]葛剑，赵韩，宋云涛.等离体破裂情形下ITER杜瓦的电磁研究J].低
温与超导,2012,40(10)：37-41.
(Ge Jian,Zhao Han,Song Yun—tao.Electromagnetic studies on ITER cr­yostat during plasma disruption[j].Cryo.&Supercond.,2012,40(10)：37-41.)
[6]Liang Chao,Ye Min—you,Yao Da—mao.The structure analysis of ITER
cryostat based on the finite element method[j].Fusion Engineering and Design,2013,88(1)：42-45
[7]Manish Kumar Pandey,Saroj Kumar Jha,Girish Kumar Gupta.Alternate
design of ITER cryostat skirt support system C].//10th Asia plasma and fusion association conference.Gandhinagar: Asia Plasma&Fus Assoc, 2015：UNSP012038
[8]虞新建，曾涛，罗云蓉.碟形封头与筒体连接处边缘应力的研究J].机
械,2008,35(1)：8-10.
(Yu Xin—jian,Zeng Tao,Luo Yun—rong.Research on stress between dish cover and cylinder joint J[.Mechanical,2008,35(1)：8-10.)
[9]龚曙光，谢桂兰.压力容器分析设计中的应力分类方法J].化工装备
技术,2000,21(3)：27-31.
(Gong Shu—guang,Xie Gui—lan.The stress classification method in ana­lysis design of pressure vessel[J[.Chemical Equipment Technology, 2000,21(3)：27-31.)
[10]范章，左树春.碟形封头在外压作用下的屈曲分析J].压力容器，
2014,31(1)：3439.
(Fan Zhang,Zuo Shu—chun.Buckling analysis of dished head under external pressure J[.Pressure Vessel,2014,31(1):34—39.)
[11]G.Vitupier,M.Meekins,C.Sborchia.ITER cryostat structural analysis
[C]〃2015IEEE26th Symposium on Fusion Engineering.Austin:IEEE,
2015：
1605261.

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。