毕业论文范文——ITER校正场线圈氦冷却管焊接数值模拟

西安航空职业学院毕业论文ITER校正场线圈氦冷却管焊接数值模拟姓 名： 专 业： 航空电子 班 级： 完成日期： 指导教师: 摘要 国际热聚变试验堆（ITER）校正场线圈（CC）氦冷却管和铠装电缆导体（CICC）铠甲焊接时，需控制超导缆的温度不高于250，以免影响其超导性能。本文利用焊接有限元分析软件SYSWELD模拟了三种方案下氦冷却管和CICC铠甲焊接的温度场变化，在理论上找到了合理的焊接工艺参数，在此基础上进行了相应的焊接试验。试验结果和模拟结果的趋势较为吻合，证明了有限元分析的正确性，为ITER多种类似结构的预焊接工艺制定提供了依据；同时，通过不同方案的对比得到了最优焊接方案。关键词 焊接 数值模拟 ITER 氦冷却管 SYSWELD国际热核聚变试验堆（International Thermonuclear Experimental Reactor，ITER）校正场线圈（Correction Coils，CC）用于补偿ITER内部由于制造误差、安装误差、变形误差及其他误差引起的磁场分布误差，由6个Top线圈、6个Side线圈和6个Bottom线圈组成1，Side线圈的局部结构如图1所示。线圈内部采用铠装电缆导体（Cable-in-Conduit Conductor，CICC）结构，其中的NbTi超导体运行于10kA电流、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，以保护超导缆的超导性能；再次，需严格控制焊接质量以保证焊缝强度，使超导线圈工作时该焊缝在电磁力作用下不被破坏。所以，对焊缝开展焊接工艺参数的研究是非常有必要的。有限元法作为一种近似解法是解决工程实际问题的重要途径，其用离散单元组成的集合体代替原结构，用近似函数表示单元内的真实场变量，从而给出离散模型的数值解，能很好地适应复杂的几何形状、材料特性和边界条件3。本文采用焊接有限元分析软件SYSWELD对氦冷却管和CICC铠甲的焊接进行了数值模拟；根据材料属性和实际焊接工艺试验得到焊缝截面进行了热源校核，根据设计的冷却工装加载边界条件，利用瞬态法进行温度场计算；通过有限元获得的焊接工艺参数，进行了相同参数的焊接试验，验证了有限元模拟结果和焊接工艺参数的可靠性。1、几何模型的建立氦冷却管由圆管段、跑道形段及两段之间扫略而成的过渡段组成，如图3。实际焊接时首先将跑道形段焊接在CICC上，然后依次焊接过渡段和圆管段。因仅跑道形段和CICC铠甲焊接的温度场可能影响内部NiTi超导体性能，故专门设计了水冷循环装置，并对此两部分的焊接温度进行了测量。因而氦冷却管的几何模型仅采用跑道形段，高度为11.2mm，并将其阶梯型坡口简化为单边V型坡口。CICC铠甲尺寸为19.219.2，内侧是0.5mm厚的花包层。由于焊接试验时需在花包层内表面布置热电偶，故采用无导体的空心CICC电缆，其几何模型也仅为铠甲和花包层，壁厚为二者之和，长度为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 有限元分析误差及减小措施 有限元计算误差来源有建模误差、离散误差和数值误差。建模误差可通过改进数学模型来减少，离散化误差可通过运用更多的单元来减小；计算机表示的数据和处理过程引入的数值误差很小，但某些物理因素和不良的离散化会使误差变大4。本文基于此原理进行相应的误差控制。2.2 焊接传热问题有限元法基于加权余量法分析热传导的过程是：根据热传导问题的微分方程和边界条件，建立其等效积分形式；利用等效积分的一般形式加权余量法的原理建立近似解法并求解。有限元法区别于传统的加权余量法和求解泛函驻值的变分法，该法不是在整个求解域上假设近似函数，而是在单元上分片假设近似函数。瞬态温度场的场变量（）在直角坐标中应满足的微分方程为5： （在内） 边界条件为： （在边界上） （在边界上） （在边界上）式中，为材料密度；为比热容；为材料沿三个主方向的导热系数；为内热源密度；为边界外法线的方向余弦；是在边界上的给定温度；是边界上的给定热流密度；为边界上的外界环境温度（自然对流）或边界层绝热壁温度（强制对流）。然后建立上述微分方程和边界条件的等效积分形式，利用伽辽金法选择任意函数，建立瞬态温度场的有限元求解方程：式中，C为热容矩阵，K为热传导矩阵，P为温度载荷向量，为节点温度列阵，为节点温度对时间的导数列阵。已将时间域和空间域的偏微分问题离散为N个节点温度的常微分方程的初值问题，求解即得温度场分布。2.3 热源数学模型采用数值模拟方法的关键是建立焊接热源模型，即确定合理的焊接热流分布函数。三维空间温度分布随瞬间热源形态不同而有所差异6。早期发展的点热源和线热源等热源分析模型不能用于有限元分析3，常用于焊接热场有限元分析的分布式热源模型有双椭球热源、高斯锥体热源、高斯平面热源及组合热源。J. Goldak提出的双椭球功率密度分布热源模型主要适用于有填充材料的MIG、TIG焊等，该模型设定体热源的前半部分为1/4椭球，后半部分为另1/4椭球，即沿焊接方向前后半椭球形状不相同。热源中心前面的区域温度梯度大，而后半部分温度梯度分布较缓，如图5所示。本文选择双椭球热源，并利用熔池边界准则对热源进行校核。 图5 双椭球热源模型 Fig. 5 Double ellipsoid heat source model前1/4椭球和后1/4椭球的功率密度函数7分别为： （1） （2）其中，分别为双半椭球体的半轴；为前后半球体内热输入份额，。 2.4 网格模型由于稀疏网格的计算结果偏差较大，通常选取焊缝长度上包含至少4个网格，但过密的单元在提高精度的同时将大量增加计算时间，且网格细分到一定程度后，计算精度的提高将微乎其微；此外，单元最大边长和最小边长之比也影响计算精度8。基于上述原则，利用Visual-Mesh划分了包含三道焊缝的全实体单元网格模型，并对边长比过大的局部网格进行了重新划分，如图6所示。 图6 氦冷却管和CICC的有限元网格模型Fig. 6 FEA Mesh model of He-inlet and CICC采用三种方案对从内到外的各层焊道进行了分段并定义网格组，具体分段方案如表1。由于打底层焊缝对花包层内侧温度影响最大，故其分段数最多。 表1 分段方案方案温度 打底层填充层盖面层方案一13 1211方案二171513方案三2422202.5 数值模拟采用瞬态法对以上三种方案进行了数值模拟。根据对氦冷却管跑道形段内壁设计的通水铜质冷却工装的进出口温度测量，发现由于流速较高，进出水口温度相差不大，因此对这段内壁的二维单元施加了恒温，其他表面和空气进行热交换。数值模拟温度场云图和热传感器预置点的温度变化曲线如下图所示，结果显示了三种方案的最高温度。 （a） （b） （c）图7 方案一（a）、方案二（b）、方案三（c）的温度场数值模拟结果Fig. 7 Numerical results of the three schemes3、试验验证对三种方案的数值模拟进行了试验验证，其中试验装夹条件、对流换热边界条件及施焊分段均与数值模拟时的加载一致。采用钨极氩弧焊进行焊接，在花包层背面布置热电偶并记录温度变化，实际试验装置及焊接样件如图8所示。 图8 氦冷却管焊接试验及样件Fig. 8 Welding test and sample piece of He-inlet and CICC根据产品对最高温度的要求，试验最高温度及其与模拟结果的对比见表2。表2 模拟结果和试验结果Tab. 2 Results of FEA and test方案温度模拟值（）实际值（）方案一294.778285.6方案二279.751262.2方案三245.65247.1由表2可以看到，焊接试验结果和数值模拟计算的趋势是一致的，但由于手工焊存在多方面不稳定因素如分段长度不均和停留时间灵活等原因，试验结果相对于模拟值的差值有一定波动。另外，实测值相对于模拟值偏低或相当，可能原因是有限元模拟时跑道形段内壁加载的温度不合适或实际焊接时水温有所变化。4、结论 采用有限元软件SYSWELD模拟了三种方案的氦冷却管和CICC铠甲焊接温度场变化，并进行了试验验证。通过模拟和试验得出了以下结论：（1）三种方案的数值模拟结果显示，随着每层分段数的增加，铠甲内侧最高温度呈降低趋势；（2）试验结果和数值模拟结果的趋势一致，表明了有限元模型及求解计算的正确性，为ITER多种类似结构的预焊接工艺制定提供了依据。（3）根据试验结果，方案三满足超导缆对铠甲内侧最高温度的要求。参考文献1曾文斌，吴维越. 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上田幸雄，村川