肖翰林，姬書得，方喜風(fēng)，王 新，高雙勝 ，王玉曉

(1.沈陽航空航天大學(xué) 航空航天工程學(xué)部(院)，沈陽 110136；2.南車青島四方機車車輛股份有限公司 技術(shù)工程部，山東 青島 266111)

?

航空宇航工程

搭接接頭直線往返擺動焊溫度及應(yīng)力場的數(shù)值模擬

肖翰林1，姬書得1，方喜風(fēng)2，王 新2，高雙勝1，王玉曉1

(1.沈陽航空航天大學(xué) 航空航天工程學(xué)部(院)，沈陽 110136；2.南車青島四方機車車輛股份有限公司 技術(shù)工程部，山東 青島 266111)

在焊接工程中采用直線往返擺動焊接，利用有限元軟件SYSWELD進行了MIG焊接搭接接頭的溫度場與應(yīng)力場的數(shù)值模擬。結(jié)果表明，溫度場與應(yīng)力場的分布規(guī)律都與直線往返擺動的運絲方式有關(guān)，在擺動焊接過程中，焊縫的部分區(qū)域被周期性的二次加熱，一次加熱區(qū)的溫度峰值(875 ℃)要低于二次加熱區(qū)域的溫度峰值(967 ℃)，縱向殘余應(yīng)力在焊縫和熱影響區(qū)主要表現(xiàn)為拉應(yīng)力。縱向殘余拉應(yīng)力在焊縫區(qū)域呈周期性的波動，且后焊區(qū)域的值高于先焊區(qū)域。

直線往返擺動焊；搭接接頭；溫度場；殘余應(yīng)力場

為了使一次焊接獲得較寬的焊縫，減少焊道次數(shù)，在實際焊接中常常使用擺動焊接，這種焊接工藝可以大幅縮短焊接時間，提高生產(chǎn)效率，因此，此工藝在實際焊接中應(yīng)用廣泛，并且逐漸趨于成熟。對于熔化焊接過程的數(shù)值模擬，國內(nèi)外學(xué)者主要研究的工藝有兩種：非擺動(直線運絲)[1]和“Z”字型擺動[2-8]。姬書得等人[2]利用串狀熱源模型對手工擺動焊接過程的溫度及應(yīng)力場進行了模擬。蹤雪梅等人[4]根據(jù)電弧擺動的方向和速度特點，建立了階梯模型，并用它模擬擺動焊過程的溫度場。李晉梅等人[5]利用有限元軟件ANSYS，基于帶狀熱源模擬了“Z”字型擺動焊接過程的溫度場。但目前未見到關(guān)于模擬直線往返擺動工藝的相關(guān)報導(dǎo)。

本文利用焊接專業(yè)有限元軟件SYSWELD，以實際焊接中使用的直線往返擺動焊接方式為例，對MIG焊接搭接接頭過程的溫度場和殘余應(yīng)力場進行了數(shù)值模擬，并對模擬結(jié)果進行了詳細的分析，為實際焊接過程中運絲方式的選擇提供了理論依據(jù)。

1 MIG焊的有限元模型

1.1 模型材料及網(wǎng)格劃分

焊絲采用SAF-ER5356，板件材料為A7N01-T4鋁合金，接頭形式為搭接接頭，上下兩板尺寸相同，均為500mm×150mm×12mm，上板開有55°的坡口。同時，模擬用的兩種材料在不同溫度下的熱物理性能參數(shù)(如密度、熱導(dǎo)率和比熱容)和力學(xué)性能參數(shù)(如屈服強度、泊松比、彈性模量和線膨脹系數(shù))見文獻[9]。根據(jù)SYSWELD材料庫的結(jié)構(gòu)形式，把相應(yīng)的材料參數(shù)輸入到材料庫文件中，建立了焊絲SAF-ER5356和A7N01-T4鋁合金的材料庫文件，用SYSWELD軟件讀取此文件便可進行計算。

搭接接頭的三維實體模型使用CATIA建立，再用Hypermesh對其進行網(wǎng)格劃分，生成SYSWELD可識別的網(wǎng)格文件.ASC。在劃分網(wǎng)格過程中，網(wǎng)格劃分得越細，計算的精度越高，同時消耗的計算時間也越長。如果網(wǎng)格劃分得過于粗大，雖然可以縮短計算時間，但是會影響到計算結(jié)果的精度。由于焊縫及其附近區(qū)域的溫度梯度較大，而遠離焊縫區(qū)域的溫度梯度較小，為了平衡計算效率和計算精度之間的矛盾，在焊縫及其附近區(qū)域?qū)⒕W(wǎng)格細分，而遠離焊縫的區(qū)域網(wǎng)格劃分得比較粗大。劃分的網(wǎng)格如圖1所示。

1.2 熱源模型及焊接工藝參數(shù)

SYSWELD中提供了多種熱源模型，其中雙橢球熱源模型的熱流密度分布可以反映出沿深度方向焊接束流對焊接溫度的影響，因此本文采用雙橢球熱源模型來描述MIG焊接過程的熱輸入[10-15]，其數(shù)學(xué)表達式如下所示。

前半部分橢球內(nèi)的熱源分布為：

圖1 模擬用三維網(wǎng)格劃分

x≥0

x<0

ff與fr為前后部分的能量分配參數(shù)，其值分別為1.6和0.4。af,ar,b與c為高斯參數(shù)，af是前半部橢球的長度，ar是后半部橢球的長度，b為熔寬的1/2，c為熔深?；谠囼灲Y(jié)果，熱源模型中的熔寬、熔深與熱效率分別是10 mm、4 mm與0.7。如圖1(b)所示的焊縫是由兩層兩道組成的；第一層為非擺動(直線運絲)焊接；第二層為直線往返擺動焊接，其焊接過程中各層的電流與電壓見表1。本文主要分析第二道的擺動焊接，擺動方式如圖2。

表1 搭接接頭MIG焊的焊接工藝參數(shù)

1.3 其它條件

在模擬過程中，板件散熱的邊界條件主要考慮其自由表面與空氣的對流散熱和輻射散熱，設(shè)置暴露在空氣中的總散熱系數(shù)為25 W/m2℃。

圖2 直線往返擺動焊接示意圖

第一道與第二道的引弧與收弧處到下板兩端的距離均是10 mm；當(dāng)?shù)谝坏篮附油瓿赡M結(jié)果中的收弧處的溫度降到70℃時，才開始進行第二道的焊接模擬；第二道的擺動焊接工藝是采用的進10 mm退4 mm。

2 結(jié)果分析與討論

2.1 溫度場模擬結(jié)果

焊件上的溫度分布反映了復(fù)雜的焊接熱過程，它是影響焊接質(zhì)量的主要因素。溫度在焊件上的不均勻分布，引起熱應(yīng)變的不均勻，最終直接導(dǎo)致焊后焊件上出現(xiàn)殘余應(yīng)力。圖3給出的是直線往返擺動工藝下焊接搭接接頭第二層焊縫的溫度場云圖。通過分析可以看出，距離焊縫越遠，材料所經(jīng)歷的溫度峰值越低，這是由于熱源作用在焊縫區(qū)，遠離焊縫區(qū)材料的溫度是由熔池區(qū)的熱傳導(dǎo)引起的。由于接頭的第二層焊縫是采用直線往返擺動工藝焊接，焊接速度小，坡口處的部分點被重復(fù)加熱，這一現(xiàn)象被稱為二次加熱。擺動在第二層焊縫上是周期性的，全段都相繼按照擺動長度出現(xiàn)二次加熱，導(dǎo)致溫度升高。一次加熱區(qū)的溫度峰值為875 ℃(圖3(a)與圖3(b))，二次加熱區(qū)的溫度峰值會更高，達到967 ℃(圖3(c)與3(d))。綜上，擺動工藝會使焊接的溫度峰值在一定范圍內(nèi)上下波動，導(dǎo)致溫度場的周期性變化。

圖4是圖1(a)與圖1(b)中所示的A、B、C三點的溫度循環(huán)曲線，三點均位于下板件，距離焊趾2 mm；點A和點C距離下板件兩端的距離均是20 mm，點B位于線段AC的中點處。焊接過程中熱源依次經(jīng)過A、B、C三點。由于搭接接頭的焊縫是由兩層焊道組成的，第一道焊接結(jié)束，等到收弧處的溫度冷卻到70 ℃時，才開始第二道的焊接，因此，試板上點A、B與C都經(jīng)歷了兩個加熱和冷卻過程(見圖4)。

當(dāng)焊接開始時整個板件的溫度相對較低，在起弧處，熱源作用于板件時大量的熱量會向周圍擴散，導(dǎo)致熔池處的熱量變小，溫度峰值較低，甚至?xí)?dǎo)致起弧處不能完全熔化。隨著焊接過程的進行，板件的溫度逐漸升高，熱源作用在溫度較高的板件上時，向周圍擴散的熱量較少，因此導(dǎo)致了焊接過程熔池區(qū)的溫度峰值逐漸升高。收弧處，由于焊接熱輸入效率未改變，但散熱面積減少，導(dǎo)致更多的熱量集中于熔池區(qū)，造成收弧處的溫度峰值急劇上升。在焊接過程中，熔池區(qū)的熱量將向四周傳遞，使鄰近區(qū)域溫度升高。因此，高的熔池區(qū)溫度會使鄰近區(qū)域的材料也經(jīng)歷較高的溫度。上述就是點A、B、C三點的溫度峰值依次逐漸升高的原因(見圖4)。由于采用了直線往返擺動焊接工藝，在焊接第二道時部分區(qū)域(如圖4中的C點的溫度循環(huán)曲線)將在短時間內(nèi)出現(xiàn)兩個溫度峰值，在“Z”字形擺動的溫度場模擬結(jié)果中也出現(xiàn)了類似的情況[3]。

圖3 擺動工藝下焊接時的溫度場云圖

圖4 點A、B與C的溫度循環(huán)曲線

2.2 應(yīng)力場模擬結(jié)果

圖5是直線往返擺動工藝下搭接接頭的縱向殘余應(yīng)力場云圖，圖6是焊縫表面上各點的縱向殘余應(yīng)力分布規(guī)律。通過分析可以得出，縱向殘余拉應(yīng)力主要集中在焊縫區(qū)和熱影響區(qū)，這與“Z”字形擺動的表面縱向殘余應(yīng)力模擬結(jié)果相似[6]；隨著到焊縫距離的增加，殘余拉應(yīng)力逐漸轉(zhuǎn)變?yōu)闅堄鄩簯?yīng)力，最后應(yīng)力值逐漸趨近于0。同時，由于溫度場的不均勻造成了焊縫的前后兩端為殘余壓應(yīng)力，中間為殘余拉應(yīng)力，殘余拉應(yīng)力和壓應(yīng)力在整個焊縫上趨近于平衡(圖5)。由于焊縫區(qū)域的材料經(jīng)歷了較大的溫度變化梯度，先被加熱到熔化狀態(tài)后再冷卻，在室溫時焊縫存在殘余壓縮塑性應(yīng)變，因此，焊縫區(qū)域有高的縱向殘余拉應(yīng)力。熱影響區(qū)的部分溫度峰值高于金屬材料的液相線，使原子的擴散能力變強，晶界遷移力大于阻力，因此晶粒獲得了很大的能量發(fā)生聚集和長大，這就導(dǎo)致了熱影響區(qū)發(fā)生再結(jié)晶現(xiàn)象生成大晶粒，長大的晶粒消除了焊縫區(qū)的部分拉應(yīng)力，導(dǎo)致焊縫區(qū)的殘余拉應(yīng)力低于熱影響區(qū)。

在擺動工藝下，焊縫區(qū)域所經(jīng)歷的溫度峰值在一定范圍內(nèi)呈周期性變化。一次加熱區(qū)的溫度峰值低于二次加熱區(qū)，溫度峰值的不同會導(dǎo)致焊后的殘余應(yīng)力值不同，因此，焊縫不同區(qū)域的殘余應(yīng)力峰值呈周期性波動，如圖6所示。同時，隨著焊接過程的進行，熔池區(qū)的溫度峰值逐漸升高，導(dǎo)致焊縫處的縱向殘余拉應(yīng)力值亦呈現(xiàn)逐漸上升的趨勢，即后焊區(qū)域的殘余拉應(yīng)力高于先焊區(qū)域(圖6)。而“Z”字形擺動的沿焊縫方向的縱向殘余應(yīng)力值在焊縫中部基本一致[2]，不存在波動和逐漸上升的現(xiàn)象。焊工根據(jù)焊接接頭形式、裝配間隙和焊條直徑等因素合理的選擇不同的擺動方式，直線往返擺動與“Z”字形擺動都有各自的優(yōu)點和適用情況，因此這兩種擺動方式在實際焊接中都被廣泛應(yīng)用。

圖5 搭接接頭的縱向殘余應(yīng)力場云圖

圖6 沿焊縫方向的縱向殘余應(yīng)力分布

3 結(jié)論

采用有限元軟件SYSWELD對直線往返擺動焊接過程進行了數(shù)值模擬，得到以下結(jié)論：

(1)由于采用直線往返擺動工藝，焊縫的部分區(qū)域經(jīng)歷了周期性的二次加熱，一次加熱區(qū)的溫度峰值為875℃，二次加熱區(qū)的溫度峰值為967℃；對于高熱傳導(dǎo)率的鋁合金來說，焊接溫度峰值隨著焊接過程的進行而逐漸升高。

(2)焊件的縱向殘余拉應(yīng)力主要集中在焊縫區(qū)和熱影響區(qū)；焊縫區(qū)的縱向殘余拉應(yīng)力分布呈現(xiàn)周期性波動，這是由于溫度場呈周期性變化；焊縫后焊區(qū)域的縱向殘余拉應(yīng)力值高于先焊區(qū)域。

[1]Li X Y,Wang L J,Yang L J,et al.Modeling of temperature field and pool formation during linear laser welding of DP1000 steel[J].Journal of Materials Processing Technology， 2014,214(9):1844-1851.

[2]姬書得，方洪淵，劉雪松，等.基于串狀熱源的手工擺動焊應(yīng)力場的數(shù)值模擬[J].焊接學(xué)報，2005,26(5):46-48.

[3]汪瓊.厚壁鋁合金擺動電弧窄間隙MIG焊工藝技術(shù)研究[D].哈爾濱：哈爾濱工業(yè)大學(xué)， 2013.

[4]蹤雪梅，吳斌，張立平，等.基于階梯模型的擺動焊接溫度場數(shù)值模擬[J].焊接學(xué)報，2014,35(11):9-12.

[5]李晉梅，雷毅.基于帶狀熱源的擺動焊接溫度場數(shù)值模擬[J].中國石油大學(xué)學(xué)報，2008,32(6):134-137.

[6]胡軍峰,方洪淵,楊建國，等.電弧擺動對焊接應(yīng)力場的影響[J].焊接學(xué)報,2007,28(1):65-68.

[7]張華軍,張廣軍,蔡春波，等.擺動焊接動態(tài)過程溫度場數(shù)值模擬[J].焊接學(xué)報，2008,29(2):69-76.

[8]Kim B H,Kang N H,Oh W T,et al.Effects of weaving laser on weld microstructure and crack for Al 6k21-T4 alloy[J].Journal of Marine Science and Technology,2011,27(1):93-96.

[9]朱瑞棟，董文超，陸善平，等.基于SYSWELD的A7N01鋁合金緩沖梁結(jié)構(gòu)焊接過程數(shù)值模擬[J].焊接,2014,(3):42-47.

[10]Piekarska W,Kubiak M.Three-dimensional model for numerical analysis of thermal phenomena in laser-arc hybrid welding process[J].International Journal of Heat and Mass Transfer,2011,54:4966-4974.

[11]李午申，白世武，嚴春研，等.9Ni鋼焊接溫度場的數(shù)值模擬[J].電焊機，2009,39(1):1-4.

[12]楊天冰，郭瑞杰.X70管道環(huán)焊縫接頭殘余應(yīng)力數(shù)值模擬[J].電焊機，2008,38(11):9-14.

[13]Manurung Y H P,Lidam R N,Rahim M R,et al.Welding distortionanalysis of multipass joint combination with different sequences using 3D FEM and experiment[J].International Journal of Pressure Vessels and Piping,2013,111-112:89-98.

[14]康惠，凌澤民，齊喜岑.基于SYSWELD對低合金鋼焊接接頭的數(shù)值模擬及實驗分析[J].熱加工工藝.2010,39(23)：200-202.

[15]張書權(quán)，王仲玨，代禮，等.基于SYSWELD的T型接頭溫度場的數(shù)值模擬[J].熱加工工藝,2011,40(7):133-135.

(責(zé)任編輯：吳萍 英文審校：劉敬鈺)

Numerical simulation of temperature and stress field of lap joint during the LRSW process

XIAO Han-lin1,JI Shu-de1,FANG Xi-feng2,WANG Xin2,GAO Shuang-sheng1,WANG Yu-xiao1

(1.Faculty of Aerospace Engineering,Shenyang Aerospace University,Shenyang 110136,China;2.Technical Engineering Department,CSR Qingdao Sifang Locomotive and Rolling Stock Co.,Ltd.,Qingdao 266111,China)

Linear reciprocating swing welding (LRSW),which has been widely used in the welding process,is simulated by the temperature field and stress field of MIG welded lap joint with the finite element software SYSWELD.The results of numerical simulation show that the distribution of temperature field and stress field is related to the wire-feeding mode of LRSW.The twice heating zones periodically appear along the weld in the process of LRSW and the peak temperature of once heating zone (875 ℃) is lower than that of twice heating zone (967 ℃).The tensile longitudinal residual stress mainly turns up in the weld and heating affected zone,while the value along the weld fluctuates periodically.The residual stress of rear welding region is higher than that of prior welding region.

linear reciprocating swing welding;Lap joint;temperature field;residual stress field

2015-01-06

國家自然科學(xué)基金(項目編號：51204111)；航空科學(xué)基金(項目編號：2013ZE54021，2014ZE54021)

肖翰林(1990-)，男，吉林白山人，碩士研究生,主要研究方向：焊接過程有限元仿真;E-mail: faithxhl@126.com；姬書得，男，河北滄州人，副教授，主要研究方向：焊接工藝及機理;E-mail:superjsd@163.com。

2095-1248(2015)04-0001-05

TG453

A

10.3969/j.issn.2095-1248.2015.04.001