DEGMAW焊接温度场与热变形的有限元分析.pdf

《DEGMAW焊接温度场与热变形的有限元分析.pdf》由会员分享，可在线阅读，更多相关《DEGMAW焊接温度场与热变形的有限元分析.pdf（70页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、分类号：T G 4 4密级：论文题目：l单位代码：学号：东岁，蚤硕士学位论文y-19 8 2 6 1 主1 0 4 2 22 0 0 3 1 18 2 0D E G M A W 焊接温度场与热变形的有限元分析作者姓名专业张明贤材料加工工程指导教师姓名专业技术职务武传松教授2 0 0 6 年5 月1 2日山东大学硕士学位论文D E G M A W 焊接温度场与热变形的有限元分析摘要D E G M A W(d o u b l e e l e c t r o d eg a sm e t a la r cw e l d i n g)焊接工艺是一种新型的高效电弧焊工艺，它以M A G 焊为主体，以T I

2、 G 焊枪作旁路，利用旁路电弧的分流作用，分走了通过焊丝的部分电流，在焊接总电流不变的前提下减小了作用于母材的电流，在保证了熔敷率的同时，减小了对母材的热输入，可以实现高熔敷率和高速焊接。本文针对该新型焊接工艺的特点，建立D B G M A W 焊接热过程的数值分析模型，为实现其工艺参数的优化奠定基础。在熟练掌握S Y S W E L D 软件的基础上，对该软件自身功能进行了改进和二次开发。提出了余高的简单处理办法；利用单元死活技术，通过不同的网格划分来考虑实际的填充金属焊接过程中焊接熔池由于电弧压力、熔滴冲击力等力的作用引起熔池表面发生的变形，使用平面高斯热源与均匀体积热源相结合的复合热源作

3、用模式，建立了适用于D E-G M A W 焊接工艺的有限元模型，并对该工艺下的温度场进行了数值模拟。结果表明：利用本文所建立的热源作用模式和余高处理办法计算出的D E G M A W 焊缝断面形状尺寸与实验结果吻合良好。计算了通过焊丝的总电流相同时D E-G M A w M I G 焊接工艺温度场的分布情况以及焊缝断面熔合线形状。结果表明：D E-G M A W 焊接时(旁路电流8 0 A，下同)焊缝的正面熔宽、熔深与实验结果吻合良好；而M I G 焊接工艺时工件已经熔透，上表面的热影响区宽度也比D E G M A W 时的结果要大得多。利用所获得的温度场数据计算的应力、应变及变形的演变规律

4、与实际情况吻合较好，但D E G M A W 焊时相应单元的平均应力、等效应变以及相应节点的最大变形量均要比M I G 焊时的结果小。本文的研究在为D E G M A W 焊接工艺参数优化提供理论数据的同时进一步印证了该工艺的优越性。关键词：D E-G M A W，热源作用模式，焊接温度场，有限元计算美国国家科学基金会国际合作项目(D M!-0 3 5 5 3 2 4)A B S T R A C 下F i n i t eE l e m e n tA n a l y s i so f t h eW e l d i n gT e m p e r a t u r eP r o f i l e sa

5、n dT h e r m a lD i s t o r t i o ni nD E G【A WA B S T R A C TD o u b l e E l e c t r o d eG a sM e t a lA r cW e l d i n g(D E-G M A W)i san o v e la r cw e l d i n gp r o c e s sw i t hh i g he f f i c i e n c y I nD E G M A W,t h eM A Ga f ci st h em a i na r cw h i l et h eT I Ga r ci sab y p a

6、s sa r c B e c a u s et h eT I Ga r ct a k e sa w a ys o m ec u r r e n t,t h ec u r r e n tf l o w e do nt h eb a s em e t a li sd e c r e a s e db a s e do nt h ep r e m i s et h a tt h et o t a lc u r r e n ti sc o n s t a n t T h u s，t h ed e p o s i t i o nm t ec a nb eg u a r a n t e e da n dt

7、 h eh e a ti n p u tt ot h eb a s em e t a li sr e d u c e d T h eh i g l ld e p o s i t i o nr a t ea n dh i g hs p e e dw e l d i n gc o u l db er e a l i z e di nD E G M A W I nt h i ss t u d y,t h en u m e r i c a la n a l y s i sm o d e lo ft h eD E-G M A Wp r o c e s si sd e v e l o p e dw i t

8、 l lc o n s i d e r i n gi t sc h a r a c t e r i s t i c s A n dt h em o d e lc o u l db et h eb a s i sf o ro p t i m i z i n gt h ep a r a m e t e ro f t h en e ww e l d i n gp r o c e s s T h ef u n c t i o n so f t h es o t h v a r eS Y S W E L Da r ee x t e n d e d As i m p l i f i e dm e t h

9、o df o rp r o c e s s i n gt h ew e l dr e i n f o r c e m e n ti sp u tf o r w a r d T h e“d e a d-a c t i v e de l e m e n t s t e c h n i q u ea n dd i f f e r e n tm e s hp l o t t i n ga r ee m p l o y e dt oc o n s i d e rt h ed e f o r m a t i o no f t h ew e l dp o o ls u r f a c eu n d e rt

10、 h ea c t i o no ft h ea r cp r e s s u r e m o l t e nd r o p l e ti m p a c ta n d o n A2-DO a u s s i a nh e a ts o u r c ea n dav o l u m e t r i ch e a ts o u r c ea r ec o m b i n e dt o g e t h e rt os i m u l a t et h eh e a ti n p u ta sac o m p o s i t eh e a ts o u r c e Af m i t ee l e m

11、 e n ta n a l y s i s(F E A)m o d e lf o rD E-G M A Wi sd e v e l o p e dw h i c hi su s e dt 0m a k en u m e r i c a ls i m u l a t i o n so ft h et e m p e r a t u r ed i s t r i b u t i o nd u r i n gD E G M A W T h er e s u l t ss h o wt h a tt h ep r e d i c t e dw e l dd i m e n s i o n sa tc

12、r o s ss e c t i o n sm a t c hw i t l lt h ee x p e r i m e n to n e s W h e nt h et o t a lc u r r e n tf l o w i n gt h ew i r e(Di sc o n s t a n t,f o rM I Gw e l d i n g，t h eb y p a s sc u r r e n tf l o w i n gt h eT I Ot o r c h(I b D)i sZ e r o，w h i l ef o rD E-G M A W,Ic o n s i s t so f

13、t w op a r t s(t h eb y p a s sc u r r e n tI b pa n dt h eb a s em e t a lc u r r e n tI b m)，a n dI bi sn o n Z e r o F o rt h es t u d yc a s e1o fD E-G M A W,I-3 3 0 A，I b p=8 0 Aa n dI b m=2 5 0 A F o rt h es t u d yc a s e2o f M I C I=3 3 0 A I b p=O Aa n dI m=3 3 0 A T h en u m e r i c a la n

14、 a l y s i so f t e m p e r a t u r eI V 山东大学硕士学位论文f i e l d so nt h ew o r k p i e c e si sc o n d u c t e da n dt h ew e l dd i m e n s i o n sa tc r o s ss e c t i o n sa r eo b t a i n e d F o rc a s e1，t h ep r e d i c t e dw e l dw i d t ha n dp e n e t r a t i o na r ei ng o o da g r e e m e

15、n tw i t ht h ee x p e r i m e n t a lm e a s u r e m e n t s F o rc a s e2，t h ew o r k p i e c ei sf u l l yp e n e t r a t e db e c a u s eo fh i g h e rh e a ti n p u t,a n dt h ew i d t ho fb o t ht h ew e l da n dt h eH A Zi sm u c hg r e a t e rt h a nt h a ti nD E-G M A W B a s e do nt h ec

16、a l c u l a t e dt e m p e r a t u r ep r o f i l e s，t h ep r e d i c t e dt h e r m a ls t r e s sa n ds t r a i na n dt h ev a r i a l i O R So ft h e r m a ld i s t o r t i o no ft h ew e l d e dw o r k p i e c e sm a t c hw e l l、i t I lt h ep r a c t i c a ls i t u a t i o n s I nD E G M A Wp

17、r o c e s s t h em e a ns t r e s s，e q u i v a l e n ts t r a i no ft h ec o r r e s p o n d i n ge l e m e n ta n dt h em a x i m u mo ft h e r m a ld e f o r m a t i o no ft h ec o r r e s p o n d i n gn o d ea r el o w e rt h a nt h o s ei nM I GI tp r o v i d e sw i t hb a s i cd a t af o rt h

18、eo p t i m i z a t i o no fD E G M A Wa n df u r t h e rv e r i f i e st h ea d v a n t a g e so f D E-G M A Wp r o c e s s K e yw o r d s：D o u b l e-E l e c t r o d eG a sM e t a lA r cW e l d i n g(D E-G M A W)，m o d e lo f h e a ts o u r c e，t e m p e r a t u r ef i e l d s，f i n i t ee l e m e

19、n ta n a l y s i sf l E A)-V-主要符号表主要符号表余两形状参数r温度余高横截面面积L熔点比热容瓦环境温度板厚主弧电压焊丝直径L么旁路电弧电压杨氏模量v o焊接速度余高高度v，送丝速度电流干伸长度旁路电弧电流x，J，z移动坐标系主电弧电流x，f，：固定坐标系导热系数搿线性膨胀系数应变硬化指数口。初始温度的线性膨胀系数蒸发潜热常数瑾。综合散热系数蒸发率r h主电弧功率系数泊松比啦旁路电弧功率系数均匀体热源功率0加权系数表面热源功率O o余高形状参数均匀体热源五表面热源占总功率的份额对流和辐射散失热流密度p密度蒸发散失热流密度仉屈服强度电弧总有效热功率毋纵向应力平面高斯热

20、源半径f加载、卸载判据时间口4qd办E矗，k七髟厶BC吼弦只，原创性声明本人郑重声明：所呈交的学位论文，是本人在导师的指导下，独立进行研究所取得的成果。除文中已经注明引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写过的科研成果。对本文的研究作出重要贡献的个人和集体，均己在文中以明确方式标明。本声明的法律责任由本人承担。论文作者签名：豕!翌国览日6 上7 乙关于学位论文使用授权的声明本人完全了解山东大学有关保留、使用学位论文的规定，同意学校保留或向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅：本人授权山东大学可以将本学位论文的全部或部分内容编入有关数据库进行检索，可

21、以采用影印、缩印或其他复制手段保存论文和汇编本学位论文。(保密论文在解密后应遵守此规定)，论文作者签名：磷，师签名：日期：乃鸯J 乙山东大学硕十学伊论文1 1 选题意义第一章绪论现代制造业对焊接生产加工提出了高质量和高效率的要求。由于电弧焊具有低成本和高适应性的特点，因此，对传统的电弧焊进行高效化改进，具有广泛的应用f；景。提高焊接效率的途径一般来说有两种：一是提高焊接速度；二是提高熔敷率。但当焊接速度提高到一定程度之后，焊缝成形会变差，出现咬边和驼峰等缺陷-。另一方面，对于传统的电弧焊接方法，电弧总是建立在电极和母材之间，焊接电流通过电弧全部作用于母材”1。如果要提高焊接速度或焊缝金属的熔敷

22、率，就势必要增大焊接电流；而随着焊接电流的增加，电弧对母材的热输入也相应增加，从而使焊接接头的热变形也随之增大，影响焊接质量。如果研发一种既能提高焊接速度和熔敷率、又能降低对母材的热输入、同时焊缝成形良好的新型电弧焊工艺，无疑具有广阔的应用I j 景和重要的实际意义。美国肯塔基大学张裕明教授研发出了这样一种工艺方法，即D E G M A W(d o u b l e e l e c t r o d eg a sm e t a la r cw e l d i n g)焊接工艺“。其工作原理如图1 所示。D E G M A W 工艺将一个T I G 焊枪与一个M A G 焊枪相组合，T I G 焊枪

23、构成旁路，M A G 焊枪与工件构成主路。流经焊丝的焊接电流，分为两部分，一是旁路电流l b。，二是施加到母材的电流碥。作用于焊丝上的焊接电流数值较高，有利于提高焊丝的熔化速度，从而提高熔敷率。T I G 焊枪构成的旁路，分流了一部分通过焊丝的焊接电流，在保证了焊缝熔敷率的同时，减小了作用于母材的热输入。由于D E G M A W 本质上属于电弧焊的改型，所以，它是低成本的高效焊接方法。美国肯塔基大学已经搭建起这种新工艺的焊接实验装置，并丌展了大量的焊接工艺实验。结果表明了D E G M A W 这种新工艺的可行性和可靠性。第一蕈绪论整lD E-G M A W 簿接互艺黎意塑F i g IS

24、c h e m a t i c s k e t c ho f t h e D E-G M A W p r o c e s s为了蜜现D E G M A W 工艺豹芏激应爱，必须系统解决这静新型焊接方法豹关键工艺理论问题和工蕊参数优化问题。如果单凭焊接工艺试验来优化工艺参数，不仅簧花费大量的人力物力，其结果也仅适斌于特定的试骏条件。因此，对D E G M A W 焊接过程遴行计算祝模拟，再辅瑷多豢豹工艺试验验证，就能够宛量分析相应的焊接温度场与焊缝形状尺寸，将为这种新工艺的参数优化提供基础数据。露丽疑蠢卡分重要的理论意义和忑程实用影 馕。为达成这露标。美国嚣塔基大学翩遗研究中心焊接研究室与出东大

25、学材辩连接按术研究掰一起合作，程撮美国N S F 资助下，深入系统地丌展这方面的研究工作。本文就是这项国际合作礤究瀑题的部分内容。1 2 高效焊接工艺研究现状一黢认为，褰效滓缓技术胃班势必褥类：一楚镑对浑叛黪荔熔敷率浮臻方法，二是针对薄板的高速度焊接方法，两炎方法是相互联系的。综合备种为提高焊接生产率而提出的焊接工艺，大致可分为以下几个方颟蛳；(I)增加送丝速度；(2)采爝多丝、多程霹接；(3)在簿谈避羧孛添掇金满耪；(4)减，j、接头截瑟哭寸；(5)采用复合热源技术，等。山东人宁砀十宁何论文l、增加送丝速度以增加送丝速度柬提高焊接生产率的典型代表是T I M E(T r a n s f e

26、r r e dI o n i z e dM o l t e nE n e r g y)焊接工艺，该工艺由加拿大人J o h nC h u r c h 于二十世纪八十年代研究成功。它是一种高性能M A G 焊工艺。主要从高速和高熔敷率两个方面来实施。其中高速度用于薄板焊接，高熔敷率用于厚板焊接。与传统M A G 相比，T I M E 工艺具有送丝速度快(从1 6 m m i n 提高到5 0 m m i n)，焊丝伸出长度长(从1 0 1 5 m m 提高到2 0 3 5 m m)、采用四元保护气体(6 5 A r,2 6 5 H e，8 C 0 2，O 5 0 2)等特点”1。将焊丝熔敷率提高

27、了2 3 倍。文献【1 1】指出：焊丝伸出长度是影响熔敷率的主要因素；采用大焊丝伸出长度匹配无氮混合气体可以大幅度提高焊丝送丝速度；实现高熔敷率焊接的核心问题是：在连续大电流下获得稳定的熔滴过渡形式。因此，各国焊接研究人员开始研究以少氦或无氦混合气体做保护，实现高熔敷率焊接的高效M A G 焊接工艺。T I M E 焊接工艺在保护气体中加入了氦气，氦气具有较高的电离能，使电弧受到压缩，电弧电压提高，弧压可达到4 8 V。同时c 0 2 使弧压进一步提高，提高了射流过渡临界电流值，使得在大电流下得到稳定的熔滴过渡过程保证焊缝成形良好。但是由于氦气在焊接过程的使用，使焊接成本高于传统M A G 焊

28、工艺“”，因此阻碍了此工艺的普及推广。2、多丝多弧焊工艺在多丝多弧焊接新工艺方面，同本、瑞士、德国等国公司在多根焊丝配以单个或多个电源方面开展了大量的焊接研究工作，在提高焊接生产速度和金属熔敷率方面取得了一些实用化的成果。例如日本的藤村浩史开发的多丝焊接系统，可用于角焊缝的高速焊接，焊速可以达到1 S m m i n“”。多丝多弧焊能增加焊接效率的原理是很容易理解的，它适用于多种弧焊方法，其中已熔化极焊接方法的应用居多“”1。对于熔化极气体保护焊，多丝多弧焊分为两种类型，一种是每根焊丝各自具有自己的焊接电源、送丝机构及调节机构：另一种是所有焊丝共用一个焊接电源，从一个导电嘴向外输送“。多丝多弧

29、焊工艺在实施过程中要求电弧之间距离近一些，这样有利于提高热源的能量密度以及热源的利用率，但是电弧之日J 距离减小后，由于电弧产生的电磁场之间差生互相干扰，使焊接过程变得不稳定，所以要解决这一问题，需要有专用的双丝焊第一章绪论接电源。同时，多丝多弧焊接提高焊接熔敷率同时也增加了对母材的热输入，因而其在生产率的提高是以牺牲焊接质量为前提的“”。3、在焊接过程中添加金属粉在焊接过程中添加金属粉的工艺，主要是用于埋弧焊，也有应用在熔化极气体保护焊的报道汹洲。在不增加电弧能量的前提下，添加会属粉，可提高熔敷效率3 0 一5 0，在实际中应用还有将多丝埋弧焊与添加会属粉结合起来工艺嘲1。由于金属粉颗粒非常

30、细小容易被电弧吹向四周，因此它的应用主要局限在埋弧焊，在熔化极气体保护焊的应用非常少。4、双面电弧焊接工艺根据所采用焊接电源的形式，可分为双电源和单电源两类汹1。与单面焊接相比，双面电弧焊接技术有如下特点删：能量集中、热影响区窄，焊缝质量好；可一次焊透中等厚度板、减小焊接工序、提高焊接生产率：电弧能量利用率高、节约能源。但是双面电弧焊接工艺的可达性差，要求工件焊接区域两侧要有焊枪的活动空间，这使得该工艺的广泛应用受到局限。5、复合热源的应用所谓复合热源是指将两种以上的热源叠加起来，以求得更强的能量密度，如在等离子束中加激光、在电弧中加激光等嘲州。国内外在这方面的研究较多。日本学者T o m o

31、 h i k oS H I D A 等人研究了利用C 0 2 激光接铝合金时复合 f i g 或M I G电弧的工艺，使得焊接过程中熔池吸收的激光能量增加。天津大学胡绳荪等人曾研制了一套激光电弧复合热源焊接装置，采用该装置在激光功率为8 0 0 W时，配以9 0 A 的焊接电弧，可与5 k W 激光焊机的焊接能力相当3。复合热源的焊接工艺尤其是激光电弧复合工艺主要是针对铝激光焊的。虽然在一定程度上能使激光能量利用率提高，但是同时会增大对工件的热输入、热影响区变宽，焊接接头质量未必能比只用激光焊时高。从这个意义上说，能量利用率的提高，也是以焊接接头质量下降为代价的“”。为了解决增大熔敷率的同时增

32、加了对母材的热输入这一矛盾，山东大学邹增大教授提出了双芯电弧焊接工艺啪1。双芯单弧焊中，工件不接电源，双芯单弧焊条的两焊芯分别接电源的两极。两焊芯与工件短路或用碳棒短路引燃后，形成单山东大学硕十学付论文一电弧。电弧可以在空日J 的任何位置引燃和燃烧。利用电弧和熔滴的热量来加热熔化母材，达到牢固的冶金结合。这种工艺主要用在材料的堆焊表面改性方面。1 3 弧焊热过程数值模拟的研究进展焊接热过程是影响焊接质量和生产率的主要因素之一。焊接热过程的准确计算和测定，是焊接冶金分析、应力变形分析以及焊接过程进行控制的I j 提。焊接热过程是极其复杂的，其复杂性主要表现为乜9 1：(1)焊接热源能量密度分布梯

33、度大，加热极不均匀；(2)焊接加热速度极快，在很短的时间内热源把大量的热能传递给焊件：(3)焊接热源相对于工件移动，传热过程不稳定；(4)焊接熔池中会属强烈运动，内部进行一系列的物理化学反应。由于焊接过程涉及的变量数目繁多，单凭积累工艺试验数据来深入了解和控制焊接过程，既不切实际又成本昂贵和费时费力。随着计算机技术的发展，通过一组描述焊接基本物理过程的数学方程柬模拟焊接过程，采用数值模拟方法求解以获得焊接过程的定量认识，即焊接过程的计算机模拟，成了一种强有力的手段。采用科学的模拟技术和少量的试验验证，以替代过去一切都要通过大量重复试验的方法，不仅节省了大量的人力和物力，而且还可以通过数值模拟解

34、决一些目前无法在实验室里进行直接研究的复杂问题。在制造业，计算机模拟与仿真可以增加材料利用率2 5，节约生产成本3 0，产品设计至实际投产的时间缩短4 0。因此计算机模拟方法为焊接科学技术的发展创造了有力的条件啪。计算机模拟方法为焊接科学技术的发展创造了有力的条件”。从上世纪5 0 年代丌始，人们开始用数值法解决传热学中的温度分布问题。随着计算机的应用和发展，数值法求解热传导微分方程向两个方向发展，即差分法和有限元法。差分法的长处是对于具有规则的几何特征和均匀的材料特征问题，它的程序设计和计算过程简单，收敛性也比较好。差分法的缺点是往往局限于规则的差分网格，且只看到节点的作用。而忽略了把节点连

35、接起来单元的贡献。有限元法可以用任意形状的网格划分区域，还可以根掘具体问题有疏有密地布置节点。另外，有限元法是用统一观点对区域内节点及边界节点列出计算格式，使各节点在精度上比较协调。有限元法具有很大的灵活性和适应性，特别适用于具有复杂形状和边界条件的问题。在热传导问题中有限元法得到广泛应用的另一个重第一章绪论要原因是：在实际应用中，温度场的计算往往服务于热应力场的计算，例如计算焊接热应力应变的动态过程及残余应力和变形，第一步就是必须进行焊接传热分析。在这种情况下，采用有限元法便于把两者统一起来。真J 下意义上的焊接热过程的数值分析开始于2 0 世纪7 0 年代，1 9 7 5 年加拿大的Z P

36、 a l e y 和H D H i b b e r t”1 用有限差分法编制了可以分析非矩形截面以及常见的单层、双层U、V 型坡口的焊接热传导的计算机程序，考虑了材料物理性能与温度的关系，并将熔化区内的单元作为加热的热源来处理，但忽略了向周围环境的散热损失，并假设工件为无限长。美国M I T 的K M a s u b c h i 等人用有限元研究了水下焊接传热问题，美国的GW K m t z y 于1 9 7 6 年的博士论文中也用有限元建立了二维焊接温度场的计算模型，并考虑了相变潜热问题“。该模型考虑了导热系数和比热随温度的变化，边界条件中也考虑了试件与周围环境介质J 日J 的辐射和对流换热

37、，但没有说明焊接热源的处理方法，忽略了在电弧运动方向上的传热。加拿大学者J G o l d a k 等应用有限元对焊接温度场进行了比较详细的分析，提出一种新的双椭球形焊接热模型m”1。该模型考虑了熔池内液体的流动和电磁力作用下，内部的磁流体动力学情况，比较清楚的给出了熔化焊时能量密度分布的结果。该模型不仅适用于焊条电弧焊接，而且还可以用于埋弧焊。但模型中忽略了辐射传热，使得该模型仍然存在一些不足。S K o u m。”建立了模拟厚板激光或等离子焊接堆焊的准稳态三维温度场的差分计算模型。在国内，1 9 8 1 年，西安交通大学唐慕尧等首先用有限元法计算了薄板准稳态焊接温度场m 3，之后，上海交通

38、大学在焊接热传导数值分析方面做了许多工作，提出了求解非线性热传导方程的变步长外推法，考虑了材料热物理性能随温度的变化以及表明的散热情况，建立了焊接温度场的有限元计算模型和相应的计算机程序。并在脉冲T I G 焊接温度场等问题进行了成功的实例分析m 蜊。哈尔滨工业大学曹振宁建立了三维T I G M I G 焊接熔池流场和热场的准稳态数值模型H 0 1。把电弧作用下的热场、流场和磁场作用进行了有效的耦合，克服了以往T I G 焊接熔池流场和热场不能处理熔透的局限性。郑炜建立了脉冲T I G 焊接熔池流场与热场动态过程的瞬念数值模型“。建i l l 彖夫学碜+学帝论文立了一套遁合于该模裂非稳念、非线

39、性、多区域、强祸合特点的数值模拟方法。张逢缨建立了基予凌态薅接的黪稳惑T K 3 辫攘热袭导摸型泔。闫最溘建立了基于P H O E N I C S 的三维瞬态T I G 焊接熔池数值模拟模型“”。赵朋成建立众熔透焊接熔池形态瞬时行为的数值模拟模型“”，作者采用了双椭网平面高斯热源模墅，劳羹讨途了熔涟臻整的翔摇。2 0 馓纪9 0 年代焊接热传导分析可以说已达到较成熟的地舻。随着计算机容量的扩大和速度的增长，三维瞬态焊接温度场的分析已不存在实质性的困难呻蝴4迄今为止，相对予T I G(钨极氟弧焊)焊接，M I G M A G 焊接过程数值模拟方面的研究工作开展得较少，主要怒由于M I G M A

40、 G 焊缝余高的存在给数值模拟豢寒一定瓣难。徐文立诗算了毫强锻会会双丝爆接嬖专瓣温疫场“”，餐其诗算逶翟中没有考虑填充金属的影响。张建勋、王军等建赢了C 0 2 气体保护焊接焊缝形状的神缎网络预测模型，将焊缝横截面上、下轮廓线简化处趱为两条抛物线“，弱弱享拳绞溺终垂赛瑟耱经颈溅浮缝形获。涤霉等奁磅究裹逮簿揍戆溽缝皎迭形成机理时，建立了平檄表面熔敷念麟形状方程，用来处理C 0 2 高速焊接时焊缝形状的问题。E P a r d o 脚1 采用有限元法建立了M I G 焊接三维稳态热传导模型，计算了M I G 簿菝溶缝懿忍蒋形获襄簿缝余裹。模黧没有考虑熔滚熬表瑟交形，挺熔池表面作为平面处理。P T

41、e k r i w a l”等也采用有限元方法建立了M I G 焊接传热模型，通过在步距上加入温度为2 3 0 0 K 的熔滴，模拟熔滴带入熔池的热量，这些萃元懿形凌、漫凄怒先确定，竣蠢考虑浮缝豹佘裹穗溶逢表嚣交形。武传松和美国威斯康辛大学的K T s a o 首次研究T M I G M A G 焊接熔池中的流体流动和传热过程”州，提出了殿定电弧二维M I G M A G 焊接对流、传热模型。1 9 9 2 年，戴传裣潜建立了运番毫蔽律羯下三缍掰l G 嬲矗G 簿羧落逸串懿滚褡漉动和传热过程的准稳态模型。在文献 5 6-5 7】中，武传松和德国的L D r o n 建立了运动电弧M I G 焊

42、接三维潞度场的数值分析模型。孙俊受拇1 建立赘A 茬i G 嬲A G 簿羧络泡流场，激发场黪模型，考虑了电弧热流密度在变形熔池表面的分御模式，以及熔滴热焓慧在熔池内部的分布。第一鼋绪论!s!s,1 1 1 1 4 1,l!i l l l l,ii l l l,i i,l l l,i*目1 4 焊接应力场数值模拟的研究现状焊接威力和变形的存在，是导数焊接裂纹和接头强度与性能降低的耋嚣因素。对焊接残余应力的产生机理及焊接过程中瞬态热应力应变场的研究，一盥是焊接工终卷关洼鲍阕题。关于媾接疲宠夔磅究骞鼹耱途径，一耱是菱毳追溯工馋的温度历史，仅基于接汰区域的力学性能指标，通过获得的非协调应交来得别焊接残

43、余应力。这种方法太多需要在实体上或模型上直接测量，不仅耗费时目J 长，费受毽器爨。在秘学抉遽菱震、薮楗攀喜不錾阂墼懿今天，久簌瓣簿接接头黪鞭裂韧性、材料可焊性、焊接裂纹的力学行为和焊接热模拟方面的研究工作，已经远不能满足于只对焊后残余应力应变状态的冶金力学特征的认识，而希望能够给出海接接头煞应力应交逡鼷懿实际嚣受戬及每一释瓣瓣凌态霉秀。雯一魏途径怒跟踪焊接过程的温度变化过程，从焊接热应力的产生入手对整个焊接过程中成力应变的演变过程进行分析计算。它不仪可以得到焊接残余应力与变形的结果，而且还戆获霉游接菝妥应力疲交j 童程静实际瑟变臻及簿一夸瓣嚣豹麓鑫露为等窍关信息，更清楚地掌握焊接残余应力和变形

44、的产生机理。因此，对焊接过程中的瞬态热应力研究在焊接力学中占有重要的地位。”。磅究挥接应力窝澎变豹数篷方法窍燕粪塑瞧蠢羧元分辑、鬻有痘交法、秧弹塑性分析，考虑相变与热应力祸合效应等。有限冗方法是焊接络构热力学分析最重要的分析手段之一，但它又是极其复杂的分析。2 G 鏊纪7 0 年代疆柬，上霹幸罐帮挺瑙荚一疆凌了考虑耪褥力学往麓与懑发有关的二维和三维焊接热弹塑性有限单元法1，并发展成为-f 7 新的学科“计算焊接力学”。E。嚣R y b i c k i 蕊酬镣入提密了繇缝对接浮瞽残余应力与交彤豹有限元模黧，模型中考虑了材料力学性能参数的非线性，材料幽弹塑性状念发生线弹性卸载，每一道焊聪几何形状改

45、变对应力的影响等。K M a s u b u c h i 汹嘲农焊接残余成力和交形懿颈溺和控制方疆进行了许多研究工作，讨论了影响分析精度的若干因索：焊接现象的复杂性：岗温下重要物瑷参数的缺乏以及系列分极辩鲍误差积累簿。翔拿穴的J，G o l d a k 辩“等对扶络赢銎j 室温时的焊接热应力_ i 薤行了分析研究，提出了各个温度段的本构方程：温度低于0 5 熔点时，材料为弹塑性；温度在0 5 0 8山东大学硕十学忙论文倍熔点时，材料为弹一粘塑性；温度超过0 8 倍熔点时，材料为线粘塑性。瑞典的L K a r l s s o n o“删等对大板拼接的焊接变形和应力进行了分析研究，特别是分析了焊缝

46、前端B J 隙的变化和点固焊的影响。另外还采用同样方法对薄壁管子的焊接残余应力和变形进行了分析。计算中考虑了材料性能的温度依赖性和相变膨胀影响在2 0 世纪8 0 年代初，国内的西安交通大学和上海交通大学等就开始了关于焊接热弹塑性理论记载数值分析方面的研究工作。西安交通大学唐慕尧旧1 与沪东造船场合作，对大板拼接单面焊终端裂纹的产生机理和防止进行了实验和数值研究。通过对热弹塑性理论本构关系的研究，用虚功原理建立了有限元方程，并编制了相应的平面问题程序，程序中考虑了材料的物理、力学参数对温度的依赖性。他还就焊接过程的力学行为的数值研究方法进行了论述，认为用热弹塑性理论的有限元法来进行整个焊接过程

47、的力学行为模拟是完全有效的。关桥等呻1 学者也编制了用于进行平板轴对称焊接应力应变分析的有限差分计算机程序和有限元程序，对于薄板氩弧焊点状热源加热的应力应变过程进行了计算分析，并讨论了有限差分法、有限元法和微分方程数值解法之自J 的差别。近些年来上海交通大学与闩本大阪大学对三维焊接应力与变形问题进行了共同研究，提出了改善计算精度和收敛性的若干途径，发展了有关三维焊接分析程序“。天津大学单平忡1 用A B A Q U S 程序对薄壁球形结构焊接残余应力和变形进行了有限元分析。清华大学鹿安理等“7 1 针对实际结构应力与变形的数值模拟，研究了动态可逆的自适应网格技术、焊缝熔敷金属填充的处理、并行计

48、算、材料性能在高温时的处理、降阶积分等关键问题。哈尔滨工业大学魏艳红等采用单元生死技术，消除了焊接构件中熔池变形对熔池尾部应力应变场的影响；通过加大材料线膨胀系数的方法，考虑凝固收缩对熔池尾部应力应变场的影响，从而建立了一种计算凝固裂纹驱动力的有效方法”。1 5 数值模拟软件的现状早期对熔池流场和热场进行数值模拟的研究者一般采用个人开发的程序进行计算。其特点是小而灵活，针对性强，输入输出系统简洁，能够迅速解决特定第一章绪论的研究问题。T Z a c h a r i am 1 和J D o m e y”分别使用套名蔓J W E L D E Rc o d e 的程序，用于模拟三维运动电弧的T I

49、G 焊接过程和重力对熔池形状的影响。利用该软件T Z a c h a r i a 还对6 0 6 1 铝的自动和非自动焊接以及对具有复杂外型的异型6 0 6 1 铝件的焊接进行了模拟。郑炜“”自行开发-j F H P T I G(F l u i d f l o w a n d h e a t t r a n s f e r i np u l s e d c u r r e n tT I Gw e l dp 0 0 1)程序模拟了脉冲T I G 焊的流场和热场。绝大多数研究者在公开发表的论文中没有提及其程序的详细内容，因而个人开发的程序的主要结构和模式均无法得知。二十世纪八十年代以来，计算技术的

50、飞速发展促进了商业软件的出现。不同于个人丌发的软件，商业软件的特点是通用性强，内容庞大，应用领域宽，有方便的模块接口，完善的输入输出系统和计算数掘的后处理系统。目前适用于T I O焊接熔池流体流动和传热过程数值模拟的商业软件”2 1 有P H O E N I C S、F L U E N T、C F X 等，其中P H O E N I C S 软件应用广泛。C h o o 使用P H O E N I C S 软件对固定电弧T I G焊接未熔池的传热和流动问题进行了数值模拟，但假定熔池的表面是不变形的。闫风洁“”应用P H O E N I C S 软件进行了二次开发，建立了运动电弧作用下三维瞬态未