EAST超导托卡马克内真空室组件焊接数值模拟.docVIP

EAST 超导托卡马克内真空室组件焊接数值模拟 高大明、吴杰峰、肖唯物 摘要:焊接过程的数值模拟是焊接技术从“技艺”走向“工程科学”的重要标志之一。焊接的非平衡加热、冷却过程导致产生焊接应力与变形,严重影响焊接过程及结构的服役行为,是焊接结构研制过程中必须解决的关键技术问题之一。本文利用焊接传热理论和Goldak 双椭球热源模型,采用SYSWELD 分析软件对EAST 超导托卡马克内真空室组件焊接温度场及应力、应变分布情况进行了初步数值模拟,并与焊接工艺试验结果和实际情况进行了验证和分析。 关键词: 热源模型 焊接温度场 焊接应力与变形 数值模拟 0 序言 作为先进制造技术中的一个重要的组成部分,未来焊接技术的发展将从“技艺”走向“工程科学”,而焊接过程数值模拟技术的发展是其重要的标志[1]。目前在焊接领域的模拟对象主要包括温度、位移、应变、应力等。而焊接非平衡加热和冷却过程中产生的焊接应力和变形,不仅影响焊接结构的制造精度,而且对结构的服役性能产生重大影响,所以焊接应力和变形才是模拟的最终对象,也一直是焊接领域广泛关注和研究的重要课题之一。但是,焊接热过程的准确计算是焊接应力应变分析的基础,焊接热过程取决于外加热源的分布形式、被焊材料的热物理性能以及被焊件与周围的换热等因素。经典的Rykalin 公式以点、线、面分布热源来模拟各种焊接热源,它存在一定的局限性。焊接过程研究早期一直停留在二维水平上[2][3][4]。近年来,随着计算机技术和有限元技术的发展以及有关焊接数值模拟软件的开发,焊接过程的三维数值模拟的研究成为该领域的研究前沿。但由于焊接过程的复杂性,焊接过程的三维数值模拟仍停留在以典型接头作为研究对象,应用到实际焊接构件还存在一定的困难。鉴于此,我们采用了焊接专用分析软件SYSWELD 对EAST 超导托卡马克内真空室两个1/16扇形段焊接过程进行实际数值模拟,并将模拟结果与试验件和真实情况对比分析,提出类似焊接结构件数值模拟应关注的几个问题。 1 焊接过程物理模型的建立 1.1 控制方程 由于焊接热过程是典型的非线形瞬态热传导过程,该问题可作如下数学描述,不失一般性,我们考虑在三维情况下,对从物体中分割出来微元平行六面体作分析,应用傅立叶公式和能量守恒定律,建立的焊接热传导控制方程为: Q z z y T y x T x t T C p + ????+ ????+ ????=??λλλρ (1 式(1中,ρ 密度(kg/m 3; C P 比定压热容(J/kg.K ; T 温度(K ; t 时间(s ; λ 热导率(W/m.K ; x 、y 、z 坐标(m ; Q 内热源强度,取决于焊接热源的功率 图1 1/16段内真空室几何模型 1.2.EAST 内真空室有限元模型的建立和计算 1.2.1几何模型的建立 EAST 超导托卡马克内真空室是一个D 形截面,双层壁、双曲面、大直径、多弧段的全焊接结构件。它由16个22.5度扇形段组成一个直径6m ,高2.7m 的球环体,图1是两个1/16段对接的焊接装配模型,该模型是利用I-DEAS 建立的,这是由于内真空室形状较为复杂,I-DEAS 相对SYSWELD 软件在实体建模方面要方便和快捷些。在local 模块划分时,充分考虑了后续有限元模型中单 元的选择及边界条件的处理。在不影响整体计算结果的前提下,对不同弧段的局部模型结构作适当的分割和简化。如减少local 模块类型和便于与弧板的对接等。 1.2.2 焊接热源模型的建立 EAST 内真空室的焊接采用钨极氩弧焊(TIG ,而高斯热源模型是在二维计算时取得的面热源形式,它没有考虑电弧在熔深方向的加热作用。实际上,三维的热源模型计算才能真实地体现焊接热过程,Goldak 研究并提出了双椭球热源模型如图2所示 热源的前半部分是个1/4椭球,后半部分是另1/4椭球。前半部分椭球内热源分布为: 图2 双椭球形热源模型 表1 热源各参数取值 图4 1/16段内真空室 local 模块的换热边界条件 1 2 22.(.(3.3.31,,,(.........3.6c t v z b y a x f t z y x e e e c b a Q f q ??= τπ π (2 同样,后半部分椭球内热源分布为: 2 2 22.(.(3.3.32 ,,,(.........3.6c t v z b y a x r t z y x e e e c b a Q f q ??=τπ π (3 在式(2、(3中参数a ,b,c 1,c 2可取根据情况可取不同的值,它们是相互独立的,f f ,f r 分别为总的输入功率在熔池前、后两部分的分配系数,f f +f r =2,通常f f 取0.6,f r 取1.4,热源的其它参数值列于表1,它是根