双丝共熔池GMAW焊接熔池流场和温度场的计算机模拟
- 格式:pdf
- 大小:1.37 MB
- 文档页数:6
《自动化技术与应用》2008年第27卷第7期28 | T echniques of Automation & Applications双丝共熔池GMAW 焊接熔池流场和温度场的计算机模拟*解生冕1,赵朋成2,黄石生2,文元美2(1.广州航海高等专科学校 计算机与信息工程系,广东 广州 510725;2.华南理工大学 机械工程学院,广东 广州 510641)摘 要:利用之前所建立的准稳态数值模型以及焊接热输入模型,对双丝共熔池GMAW焊接过程熔池成型、流场和温度场进行计算机模拟,并对模拟结果给予实验验证。
关键字:焊接熔池;流场;温度场;计算机模拟中图分类号:TP391.9 文献标识码:B 文章编号:1003-7241(2008)07-0028-06Computer Simulation of the Fluid Flow Field and the T empera-ture Field of a Twin-wire Co-Welding PoolXIE Sheng-mian 1, ZHAO Peng-cheng 2, HUANG Shi-sheng 2, WEN Yuan-mei 2(puter and Information Engineering Dep., Guangzhou Maritime College, Guangzhou 510725, China;2.Mechatronic Engineering College, South China Univ. of Tech., Guangzhou 510641, China)Abstract: The forming process of the welding pool, the fluid flow and the temperature field in the welding pool are simulated byusing the quasi-steady state numerical model and the welding heat model. The simulation the results are also verifieed validated by experiments.Keywords: welding pool; fluid flow field; temperature field; computer simulation*基金项目:国家自然科学基金(编号:50375054)收稿日期:2008-05-061 引言根据之前所建立的准稳态数值模型以及焊接热输入模型,对低碳钢的焊接未熔透熔池的三维形态、工件温度场和熔池内流体的流场进行数值模拟,对模拟结果需要给予实验验证。
2 计算过程和程序编制2.1 三维准稳态数值计算的主程序熔池表面变形的引入使得准稳态问题中流场和温度场的计算比稳态问题复杂。
首先要基于温度场的计算结果确定熔池形状,之后才能计算熔池的表面变形,而温度场和流场通过贴体坐标系与表面变形相耦合,同时流场又通过工件的热物性参数与温度场相耦合。
表面变形、流场和温度场之间的相互耦合关系如图1所示。
这是一个非线性、强耦合的数值求解问题。
根据SIMPLE方法的思想,采用表面变形、流场和热场分别求解的分离变量法。
每一个循环中,三个因变量的值依次改进一次,直到满足收敛条件。
模拟计算过程可以分以下几步:1.计算工件上的热场。
2.由热场的计算结果确定熔池的三维形状,根据熔池表面方程计算熔池的表面变形。
3.进行坐标系变换,计算贴体坐标系下熔池的流场,得到收敛的流场。
计算贴体坐标系下工件的热场,得到收敛的热场。
5.重复2 ̄4步,不断改进计算结果,直至相邻两次循环的表面变形、热场和流场满足一定的收敛条件。
图1 表面变形、热场和流场的耦合关系工业控制与应用Industry Control and Applications《自动化技术与应用》2008年第27卷第7期Techniques of Automation & Applications | 29模拟计算的主程序采用模块式设计,流程图如图2。
当流场和热场满足收敛条件时,迭代停止。
采用的收敛判据如下:005.0/1≤−+NNNTTT,005.0/1≤−+NNNUUU005.0/1≤−+NNNVVV,005.0/1≤−+NNNWWW设定收敛判据时,考虑了计算速度与计算精度之间的矛盾,在满足精度的前提下,充分考虑到了计算的经济性。
工件材料的比热、导热系数和液态金属的粘度等物性参数是温度的函数。
在本模型中,它们与温度的关系被做成子程序,在计算熔池的热场和流场程序中随时调用。
2.2 表面变形子程序采用Gauss-Seidel迭代法逐点计算熔池上下表面的变形量(图3)。
首先根据热场的计算结果确定熔池的边界,计算电弧压图2 计算流场和热场的主程序流程图工业控制与应用Industry Control and Applications图3 熔池表面变形计算子程序力。
用表面变形计算方程计算熔池自由表面的变形,判断是否收敛,重复计算表面变形,直到满足收敛条件。
2.3 熔池热场和流场的计算流场的计算采用SIMPLEC方法。
热场的计算过程只需求解热能守恒方程及其边界条件。
计算过程中,流场和热场的计算子程序采用较为严格的收敛条件:)(11KTTNN≤−+,001.0/1≤−+NNNUUU,001.0/1≤−+NNNVVV,001.0/1≤−+NNNWWW2.4 程序的适用性本程序是为运动电弧作用下的三维准稳态GMAW焊接过程的数值模拟设计的,适用性强。
设定焊丝熔化速度为零,可以模拟固定和运动电弧作用下GTAW焊准稳态过程。
3 焊接参数的设定3.1 试件材料、尺寸和物性参数本文模拟计算的焊接材料为Q235低碳钢(A3低碳钢)。
由于焊接工件沿焊接电弧运动方向的中心线对称,因此取其中的一半进行计算。
不均匀网格系统被固定到工件上,试件被划分了384×52×10个网格单元。
沿着y方向,在熔池的中心线附近,由于热输入集中,因此网格密度最大;在热影响区网格密度逐步减小;在试件的边缘,温度相对于中心较低,因此网格密度最小。
采用不均匀网格系统的主要目的是为了减少网格数目,从而减少计算机的内存使用量,加快运算速度。
对该材料的物性参数也进行了相应的考虑。
3.2 熔滴过渡参数的确定在GMAW焊接中,焊接电流是决定熔滴过渡方式的最基本参数,调整它可以得到需要的熔滴过渡形式,电流小时发生滴状过渡,电流相对高时发生射流过渡。
有研究表明[1],对于低碳钢焊丝,采用80%Ar-20%CO2保护时,随着焊接电流的增加,熔滴从滴状过渡到射滴过渡,然后到射流过渡。
熔滴过渡形式影响熔滴参数,建立不同条件下焊接工艺参数与熔滴参数的函数关系是对GMAW焊接熔池行为进行数值分析的重要问题,目前用于描述熔滴的过渡现象的理论有静力平衡理论和收缩不稳定理论[1,2,3]等。
但是由于影响熔滴过渡的因素较复杂,由这些理论建立的熔滴参数计算模型存在较多难于确定的因子,误差较大并且应用不方便。
本文根据目前的一些文献测定结果确定熔滴参数与焊接工艺参数的关系。
本文选用直径为1.2mm的低碳钢焊丝,保护气体为80%Ar-20%CO2,文献[4, 5]测定了不同焊接电流下低碳钢焊丝的熔滴参数,本文对测量数据进行统计分析得出熔滴半径rd和焊接电流I的关系式为:熔滴到达熔池表面的速度vd与焊接电流I可用下面的统计关系式表示,0626.1106377.1103829.9225d+×−×=−−IIv熔滴热焓量可根据由H Idr=+0551917151801..计算,文献[6]给出了低碳钢GMAW焊接熔池的平均热焓量Hv=1.《自动化技术与应用》2008年第27卷第7期30 | T echniques of Automation & Applications45~1.61(KJ/g),可得熔滴带入熔池的热焓量hw=Hdr–Hv。
焊丝的熔化速度Sm与I、V、We、d的关系由下式确定[7],()S I IV d I W dm e =×−×−×+×⋅−−−−8999710242981025198101525410101182324..././其中,We是焊丝干伸长,d是焊丝直径。
在进行数值分析时必须事先给定焊接基本参数I、V、We、d,可由上述公式计算出rd、Sm、vd、Hdr,其它计算所需的熔滴参数由下式求出,m r d d 3=43πρ,dm 4m S d 2ρπν=。
3.3 焊接参数的测量和确定双丝共熔池GMAW焊接参数包括脉冲电流、电弧电压、焊接速度、焊丝直径以及送丝速度等。
其中电弧电压和脉冲电流信号采集系统测量,焊接速度和送丝速度提前设定。
由于焊接电流是脉冲电流,在数值模拟过程中使用了加权平均法,将脉冲电流转换成了恒定电流。
焊接电流随时间的变化关系如图4所示。
则加权平均电流为:p b I t t t I t t t I 212211+++=。
另外,由于本文研究的是准稳态熔池,经过一个焊接周期以后前丝焊接电流对熔池的表面变形影响不大。
因此在计算表面变形时使用的是后丝电流,而计算熔池的流场和热场则采用前丝和后丝加权平均电流之和,即前丝和后丝的合并焊接电流。
4 模拟结果和讨论利用所建立的模型和算法,对双丝共熔池GMAW焊接过程中到达准稳态时熔池的三维形态、工件的温度场和熔池内液体金属的流场进行了数值模拟。
4.1 熔池三维形态的计算熔池三维形态包括熔池三维尺寸和熔池自由表面的变形。
取前丝在工件表面的投影点为原点,建立运动坐标系,即坐标系随着前丝同步移动。
计算结果都转换为移动坐标系下数据。
为了考察双丝GMAW焊接熔池的三维尺寸和熔池自由表面的变形,对低碳钢进行了堆焊试验,焊接参数分别为:(1)合并电流193A(前丝95A,后丝98A,电流频率均为100Hz);电弧电压前丝23.1V,后丝19.6V;焊接速度为0.5m/min。
(2)合并电流260A(前丝130A,后丝130A,电流频率均为100Hz);电弧电压前丝22.4V,后丝22.4V;焊接速度为0.5m/min。
不同焊接电流得到的熔池几何形状如图5和图6所示。
(a) 顶视图(b) 侧视图(c) 前视图(x =0mm)图5 焊接电流为193A 时熔池三维形状的计算结果(a) 顶视图(b) 侧视图(c) 正视图(x =0mm)图6 焊接电流为260A 时熔池三维形状的计算结果图4 脉冲电流的加权平均处理工业控制与应用Industry Control and Applications《自动化技术与应用》2008年第27卷第7期Techniques of Automation & Applications | 31由于双丝GMAW模型中使用的焊接电流是合并电流,所以图5和图6均为滴状过渡时熔池几何形状的计算结果。