1、西安航空职业学院毕业论文ITER 校正场线圈氦冷却管焊接数值模拟姓 名: 专 业: 航空电子 班 级: 完成日期: 指导教师: 摘要 国际热聚变试验堆(ITER)校正场线圈(CC)氦冷却管和铠装电缆导体(CICC)铠甲焊接时,需控制超导缆的温度不高于 250,以免影响其超导性能。本文利用焊接有限元分析软件 SYSWELD 模拟了三种方案下氦冷却管和 CICC 铠甲焊接的温度场变化,在理论上找到了合理的焊接工艺参数,在此基础上进行了相应的焊接试验。试验结果和模拟结果的趋势较为吻合,证明了有限元分析的正确性,为 ITER 多种类似结构的预焊接工艺制定提供了依据;同时,通过不同方案的对比得到了最优焊
2、接方案。关键词 焊接 数值模拟 ITER 氦冷却管 SYSWELD国际热核聚变试验堆(I nternational Thermonuclear Experimental Reactor, ITER)校正场线圈(Correction Coils,CC)用于补偿 ITER 内部由于制造误差、安装误差、变形误差及其他误差引起的磁场分布误差,由 6 个 Top 线圈、6 个 Side 线圈和 6 个 Bottom 线圈组成 1,Side 线圈的局部结构如图 1 所示。线圈内部采用铠装电缆导体(Cable-in-Conduit Conductor,CICC)结构,其中的 NbTi 超导体运行于 10kA
3、 电流、4.2K 温度下 2。氦冷却管是超导体和提供低温环境的液氦管路的连接部分。图 1 Side 线圈局部结构 图 2 氦冷却管和 CICC 结构Fig. 1 Local structure of Side CC Fig. 2 Structure of He-inlet and CICCCICC 由铠甲、花包层和 NbTi 超导体组成。氦冷却管和 CICC 的结构如图 2 所示,采用钨极氩弧焊将液氦冷却管和铠甲连接起来,此焊缝极为重要,而且焊接难度大,主要体现在:首先,采用全位置焊,空间狭小,并要求全焊透;其次,在焊接过程中需控制超导缆处的温度不能高于 250,以保护超导缆的超导性能;再次,
4、需严格控制焊接质量以保证焊缝强度,使超导线圈工作时该焊缝在电磁力作用下不被破坏。所以,对焊缝开展焊接工艺参数的研究是非常有必要的。有限元法作为一种近似解法是解决工程实际问题的重要途径,其用离散单元组成的集合体代替原结构,用近似函数表示单元内的真实场变量,从而给出离散模型的数值解,能很好地适应复杂的几何形状、材料特性和边界条件 3。本文采用焊接有限元分析软件 SYSWELD 对氦冷却管和 CICC 铠甲的焊接进行了数值模拟;根据材料属性和实际焊接工艺试验得到焊缝截面进行了热源校核,根据设计的冷却工装加载边界条件,利用瞬态法进行温度场计算;通过有限元获得的焊接工艺参数,进行了相同参数的焊接试验,验
5、证了有限元模拟结果和焊接工艺参数的可靠性。1、几何模型的建立氦冷却管由圆管段、跑道形段及两段之间扫略而成的过渡段组成,如图 3。实际焊接时首先将跑道形段焊接在 CICC 上,然后依次焊接过渡段和圆管段。因仅跑道形段和 CICC 铠甲焊接的温度场可能影响内部 NiTi 超导体性能,故专门设计了水冷循环装置,并对此两部分的焊接温度进行了测量。因而氦冷却管的几何模型仅采用跑道形段,高度为 11.2mm,并将其阶梯型坡口简化为单边 V 型坡口。CICC 铠甲尺寸为 19.219.2,内侧是 0.5mm 厚的花包层。由于焊接试验时需在花包层内表面布置热电偶,故采用无导体的空心 CICC 电缆,其几何模型
6、也仅为铠甲和花包层,壁厚为二者之和,长度为 120mm。由此建立的总体几何模型如图 4 所示。图 3 氦冷却管组成和 CICC 导体截面 图 4 CICC 和氦冷却管几何模型Fig. 3 He-inlet components and CICC section Fig. 4 Geometrical model of He-inlet and CICC2、有限元分析2.1 有限元分析误差及减小措施有限元计算误差来源有建模误差、离散误差和数值误差。建模误差可通过改进数学模型来减少,离散化误差可通过运用更多的单元来减小;计算机表示的数据和处理过程引入的数值误差很小,但某些物理因素和不良的离散化会使误
7、差变大 4。本文基于此原理进行相应的误差控制。2.2 焊接传热问题有限元法基于加权余量法分析热传导的过程是:根据热传导问题的微分方程和边界条件,建立其等效积分形式;利用等效积分的一般形式加权余量法的原理建立近似解法并求解。有限元法区别于传统的加权余量法和求解泛函驻值的变分法,该法不是在整个求解域上假设近似函数,而是在单元上分片假设近似函数。瞬态温度场的场变量 ( )在直角坐标中应满足的微分方程为 5:,xyzt(在 内) 0xyzckkkQt 边界条件为:(在 边界上)= 1(在 边界上)xyzknknq2(在 边界上)()xyzah3式中, 为材料密度; 为比热容; 为材料沿三个主方向的导热
8、系数; 为内热源c,xyzk Q密度; 为边界外法线的方向余弦; 是在 边界上的给定温度;,xyzn=(,)t1是边界 上的给定热流密度; 为 边界上的外界环境温度(自然对()qt2a3流)或边界层绝热壁温度(强制对流) 。然后建立上述微分方程和边界条件的等效积分形式,利用伽辽金法选择任意函数,建立瞬态温度场的有限元求解方程: CKP式中,C 为热容矩阵,K 为热传导矩阵,P 为温度载荷向量, 为节点温度列阵, 为节点温度对时间的导数列阵。已将时间域和空间域的偏微分问题离散为 N 个节点温度 的常微分方程的初值问题,求i解即得温度场分布。2.3 热源数学模型采用数值模拟方法的关键是建立焊接热源
9、模型,即确定合理的焊接热流分布函数。三维空间温度分布随瞬间热源形态不同而有所差异 6。早期发展的点热源和线热源等热源分析模型不能用于有限元分析 3,常用于焊接热场有限元分析的分布式热源模型有双椭球热源、高斯锥体热源、高斯平面热源及组合热源。J. Goldak 提出的双椭球功率密度分布热源模型主要适用于有填充材料的 MIG、TIG 焊等,该模型设定体热源的前半部分为 1/4 椭球,后半部分为另 1/4 椭球,即沿焊接方向前后半椭球形状不相同。热源中心前面的区域温度梯度大,而后半部分温度梯度分布较缓,如图 5 所示。本文选择双椭球热源,并利用熔池边界准则对热源进行校核。图 5 双椭球热源模型Fig
10、. 5 Double ellipsoid heat source model前 1/4 椭球和后 1/4 椭球的功率密度函数 7分别为:(1)22fffhfh633(,)exp(-0yzqxyz xabcbc( ) ),(2)22rrrhrh(,)(-( ) ),其中, 分别为双半椭球体的半轴; 为前后半球体内热输入份额,frab、 、 、 cfr、。 fr+=22.4 网格模型由于稀疏网格的计算结果偏差较大,通常选取焊缝长度上包含至少 4 个网格,但过密的单元在提高精度的同时将大量增加计算时间,且网格细分到一定程度后,计算精度的提高将微乎其微;此外,单元最大边长和最小边长之比也影响计算精度
11、8。基于上述原则,利用 Visual-Mesh 划分了包含三道焊缝的全实体单元网格模型,并对边长比过大的局部网格进行了重新划分,如图 6 所示。图 6 氦冷却管和 CICC 的有限元网格模型Fig. 6 FEA Mesh model of He-inlet and CICC采用三种方案对从内到外的各层焊道进行了分段并定义网格组,具体分段方案如表 1。由于打底层焊缝对花包层内侧温度影响最大,故其分段数最多。表 1 分段方案方案 温度 打底层 填充层 盖面层方案一 13 12 11方案二 17 15 13方案三 24 22 202.5 数值模拟采用瞬态法对以上三种方案进行了数值模拟。根据对氦冷却管
12、跑道形段内壁设计的通水铜质冷却工装的进出口温度测量,发现由于流速较高,进出水口温度相差不大,因此对这段内壁的二维单元施加了恒温,其他表面和空气进行热交换。数值模拟温度场云图和热传感器预置点的温度变化曲线如下图所示,结果显示了三种方案的最高温度。(a)(b)(c)图 7 方案一(a) 、方案二( b) 、方案三(c)的温度场数值模拟结果Fig. 7 Numerical results of the three schemes3、试验验证对三种方案的数值模拟进行了试验验证,其中试验装夹条件、对流换热边界条件及施焊分段均与数值模拟时的加载一致。采用钨极氩弧焊进行焊接,在花包层背面布置热电偶并记录温度
13、变化,实际试验装置及焊接样件如图 8 所示。图 8 氦冷却管焊接试验及样件Fig. 8 Welding test and sample piece of He-inlet and CICC根据产品对最高温度的要求,试验最高温度及其与模拟结果的对比见表 2。表 2 模拟结果和试验结果Tab. 2 Results of FEA and test方案 温度 模拟值() 实际值()方案一 294.778 285.6方案二 279.751 262.2方案三 245.65 247.1由表 2 可以看到,焊接试验结果和数值模拟计算的趋势是一致的,但由于手工焊存在多方面不稳定因素如分段长度不均和停留时间灵活等
14、原因,试验结果相对于模拟值的差值有一定波动。另外,实测值相对于模拟值偏低或相当,可能原因是有限元模拟时跑道形段内壁加载的温度不合适或实际焊接时水温有所变化。4、结论采用有限元软件 SYSWELD 模拟了三种方案的氦冷却管和 CICC 铠甲焊接温度场变化,并进行了试验验证。通过模拟和试验得出了以下结论:(1 )三种方案的数值模拟结果显示,随着每层分段数的增加,铠甲内侧最高温度呈降低趋势;(2 )试验结果和数值模拟结果的趋势一致,表明了有限元模型及求解计算的正确性,为ITER 多种类似结构的预焊接工艺制定提供了依据。(3 )根据试验结果,方案三满足超导缆对铠甲内侧最高温度的要求。参考文献1曾文斌,
15、吴维越. ITER 校正场线圈的结构分析J,核聚变与等离子体物理,2009,29(3): 2642692 A. Foussat, N. Dolgetta, C. Jong, et al. From design to development phase of the ITER Correction Coils, IEEE T APPL SUPERCONJ, 2011, 21(3): 19601963 3张彦华. 焊接力学与结构完整性原理M. 北京:北京航空航天大学出版社,2007:26774罗伯特D库克,戴维S马尔库斯,迈克尔E普利沙,等. 有限元分析的概念与应用M. 第 4版. 西安:西安交通大学出版社,2007: 45王勖成. 有限单元法M,北京:清华大学出版社,2003: 4414486上田幸雄,村川英一,麻宁绪(日). 焊接变形和残余应力的数值计算方法与程序M. 成都:四川大学出版社. 2008: 87 7武传松. 焊接热过程与熔池形态M,北京:机械工业出版社,2008:278孔祥谦. 有限单元法在传热学中的应用M. 第三版. 北京:科学出版社,1998: 39 致 谢感谢我的同学和朋友,在我写论文的过程中给予我很多素材,还在论文的撰写和排版过程中提供热情的帮助。由于我的学术水平有限,所写论文难免有不足之处,恳请各位老师和学友批评和指正!