尹文鋒,黃本生
(1.西南石油大學(xué) 工程訓(xùn)練中心,四川 成都 610500;2.西南石油大學(xué) 材料科學(xué)與工程學(xué)院,四川 成都 610500)
316L不銹鋼是一種超低碳奧氏體不銹鋼,大量應(yīng)用于防腐、防銹設(shè)備。但不銹鋼的價(jià)格偏高,往往采用不銹鋼和碳鋼的異鋼種焊接,以提高經(jīng)濟(jì)效益[1-2]。Q235為珠光體碳素結(jié)構(gòu)鋼,綜合性能較好,價(jià)格便宜,應(yīng)用廣泛。316L不銹鋼與Q235碳鋼的異種鋼焊接在設(shè)計(jì)制造中經(jīng)常被采用。但異種鋼的化學(xué)成分、物理性能存在較大差異,給焊接工藝制定帶來了困難[3]。因此,分析316L/Q235異種鋼焊接過程指導(dǎo)焊接工藝設(shè)計(jì)具有重要意義。隨著計(jì)算機(jī)技術(shù)的發(fā)展,焊接過程數(shù)值模擬為焊接過程分析、焊接工藝優(yōu)化提供了新途徑[4]。焊接溫度場對焊縫及熱影響區(qū)組織成分有很大影響,并且焊接溫度場是進(jìn)行焊接殘余應(yīng)力和焊接變形預(yù)測、評定焊接工藝合理性的重要依據(jù)[5]。
在此運(yùn)用Ansys軟件對316L/Q235異種鋼平板對接焊的動(dòng)態(tài)溫度場進(jìn)行數(shù)值模擬。分析了316L/Q235異種鋼平板對接焊的動(dòng)態(tài)焊接過程,探討了焊接溫度場的分布特點(diǎn)及焊接熱循環(huán)規(guī)律。
母材為316L不銹鋼和Q235碳鋼,焊接溫度場數(shù)值模擬需要的材料參數(shù)有:導(dǎo)熱系數(shù)、密度、比熱容、對流換熱系數(shù),它們與溫度成非線性關(guān)系,需要在Ansys軟件中建立在某些已知溫度下的材料參數(shù)表,通過插值法和外推法確定未知溫度下的參數(shù)值[6]。
模擬采用焊條電弧焊焊接溫度場。試件尺寸為100 mm×50 mm×8 mm,平板對接焊,開 V形坡口,坡口角度70°,焊縫采用兩道焊,第1道焊后冷卻時(shí)間90 s,第2道焊后冷卻時(shí)間240 s。焊接電流150~180 A,電弧電壓 25~30 V,焊接速度 3 mm/s。
由于焊接熱循環(huán)過程快速且不均勻,焊縫附近的溫度梯度很大,需采用較細(xì)密的網(wǎng)格,遠(yuǎn)離焊縫的地方,溫度梯度較小,可采用較稀疏的網(wǎng)格,既可以保證計(jì)算精度,又可以提高運(yùn)算速度。熱分析選擇三維熱實(shí)體單元solid70,它可用于三維穩(wěn)態(tài)和瞬態(tài)熱分析[7]??紤]到空氣對流對焊接溫度場的影響,在焊件表面創(chuàng)建熱表面效應(yīng)單元surf152,并將熱對流邊界條件則施加在表面效應(yīng)單元上。焊件劃分網(wǎng)格后如圖1所示。
圖1 焊件有限元網(wǎng)格模型Fig.1 Finite element mesh model of weldment
實(shí)踐證明,對于焊條電弧焊,采用高斯熱源模型可以得到比較滿意的模擬結(jié)果。選用單元內(nèi)部生熱的熱源模式[8]。
焊接母材初始溫度為25℃,在焊件模型外表面通過表面效應(yīng)單元施加對流邊界條件。為模擬熱源的移動(dòng),可將施加熱載荷的過程按時(shí)間分為n個(gè)載荷步,再將焊縫沿焊接方向等分為n段,在每個(gè)時(shí)間段施加對應(yīng)的熱流密度,當(dāng)下一段加載時(shí),刪除前面所加的高斯熱流密度,以前面所得的溫度場作為下一次加載的初始條件。加載過程應(yīng)用APDL編程實(shí)現(xiàn)[9]。焊后冷卻模擬將焊后溫度作為初始溫度,刪除焊接過程施加的高斯熱流密度,設(shè)定載荷步結(jié)束時(shí)間為冷卻時(shí)間進(jìn)行求解。
圖2為第1道焊剛開始時(shí)焊件的溫度分布云圖(t=2 s)。當(dāng)焊接熱源作用在焊縫單元上時(shí),熱源中心的溫度從室溫25℃升高到1664℃,達(dá)到焊材熔點(diǎn),形成熔池。圖3為第1道焊結(jié)束時(shí)溫度分布云圖(t=31 s)。熱源中心溫度達(dá)到1976℃比開始焊接時(shí)高,這是因?yàn)榍懊娴暮附舆^程存在預(yù)熱,與焊接熱源溫度的疊加。Q235的導(dǎo)熱系數(shù)大于316L,更利于熱量的傳導(dǎo),Q235焊板溫度上升速度快于316L,Q235焊板溫度總體要高于316L。熱源中心前方的等溫線比較密集,溫度梯度較大;熱源后方的等溫線分布比較疏散,溫度梯度稍小。
圖2 t=2 s時(shí)焊件溫度分布Fig.2 Temperature distribution of weldment at t=2 s
圖3 t=31 s時(shí)焊件的溫度Fig.3 Temperature distribution of weldment at t=31 s
圖4為第1道焊結(jié)束后,焊件開始冷卻時(shí),焊件的溫度分布云圖(t=36 s)。焊件高溫區(qū)域溫度梯度較大,溫度大幅降低;遠(yuǎn)離焊縫區(qū)域溫度梯度較小,溫度緩慢上升。Q235溫度變化較316L更快,這是因?yàn)楹缚p溫度明顯高于焊件及環(huán)境空氣溫度,焊接間的熱傳導(dǎo)較焊件與空氣間的熱對流占優(yōu)勢,因此,熱傳導(dǎo)系數(shù)較大的Q235焊板溫度更容易升高。圖5為第1道焊冷卻約70 s后焊件的溫度分布云圖(t=117 s)。焊件最高溫出現(xiàn)在第1道焊焊縫的末端及中間區(qū)域,其原因是焊縫末端在第1道焊接完成后溫度最高,且其冷卻時(shí)間比焊縫其他區(qū)域短。焊縫中間區(qū)域由于受到長時(shí)間熱源作用使該區(qū)域周圍溫度升高,溫度梯度較小,降溫過程也相對緩慢。經(jīng)過一段時(shí)間的冷卻,由于焊件間的熱傳導(dǎo),以及與空氣存在對流散熱,焊接各部分的溫度趨于均勻。
圖4 t=36 s時(shí)焊件溫度分布Fig.4 Temperature distribution of weldment at t=36 s
圖5 t=117 s時(shí)焊件的溫度分布Fig.5 Temperature distribution of weldment at t=117 s
為保證兩道焊焊接溫度相差不大,第2道焊采用較低的焊接電流。圖6為第2道焊開始時(shí)焊件的溫度分布云圖(t=123 s)。焊板整體溫度高于279℃,說明第一道焊對焊件產(chǎn)生了預(yù)熱效果。焊接熱源中心溫度為1602℃高于焊材熔點(diǎn),能夠形成熔池。圖7為第2道焊接過程中的焊件溫度分布云圖(t=138 s)。由于前面焊接過程的加熱,焊接熱源中心溫度升高2073℃。焊件最低溫度265℃,較第2道焊開始時(shí)的最低溫度279℃有所下降,這說明焊接過程中,焊件與空氣存在對流散熱,焊件的整體溫度場分布是焊接熱源加熱以及焊接表面散熱共同作用下的動(dòng)態(tài)平衡。
圖6 t=123 s時(shí)焊件溫度分布Fig.6 Temperature distribution of weldment at t=123 s
圖7 t=138 s時(shí)焊件的溫度Fig.7 Temperature distribution of weldment at t=138 s
圖8為第2道焊后自然冷卻約210 s后的焊件溫度場分布云圖(t=378 s)。Q235焊板溫度整體低于316L,這是因?yàn)镼235與空氣的對流換熱系數(shù)比316L與空氣的對流換熱系數(shù)大,當(dāng)沒有焊接熱源加熱,焊件與空氣的對流散熱就占主導(dǎo)地位,因此Q235焊板溫度就比316L降低得更快。
圖8 t=378 s時(shí)焊件溫度分布Fig.8 Temperature distribution of weldment at t=378 s
為分析焊縫及其附近熱影響區(qū)的焊接熱循環(huán),在焊件上選擇三條路徑,作出路徑上各節(jié)點(diǎn)的熱循環(huán)曲線。路徑布置如圖9所示。路徑1在316L焊板表面中央垂直于焊縫,各節(jié)點(diǎn)間距為10 mm。路徑2在焊板Q235表面平行于焊縫,距焊縫中央6 mm,各節(jié)點(diǎn)間距為15 mm。路徑3在兩塊焊板表面垂直于焊縫,各節(jié)點(diǎn)間距為20 mm。
圖9 路徑布置示意Fig.9 Path layout diagram
路徑1上節(jié)點(diǎn)在焊接過程中的熱循環(huán)曲線如圖10所示。焊接采用兩道焊,因此各點(diǎn)溫度均存在兩個(gè)峰值。當(dāng)焊接熱源距離較遠(yuǎn)時(shí),各點(diǎn)溫度均為室溫;當(dāng)熱源經(jīng)過特征點(diǎn)所在路徑時(shí),各點(diǎn)的溫度迅速上升,并且由于焊接加熱的局部性,越靠近焊縫中央其溫度上升越快,溫度值也更高。當(dāng)熱源遠(yuǎn)離該路徑后,各特征點(diǎn)溫度又迅速下降,并且溫度越高的點(diǎn)降溫越快。這是因?yàn)楹讣鞑糠珠g存在熱傳導(dǎo),并且溫度梯度很大,因此熱傳導(dǎo)及散熱迅速。冷卻階段,各點(diǎn)溫度變化趨勢基本相同,在開始冷卻時(shí)溫度下降較快,而后逐漸趨于穩(wěn)定,下降減慢。
圖10 路徑1上各點(diǎn)熱循環(huán)曲線Fig.10 Thermal cycle curves of path 1
圖11 路徑2上各點(diǎn)熱循環(huán)曲線Fig.11 Thermal cycle curves of path 2
路徑3上各點(diǎn)在第1道焊接過程中的熱循環(huán)曲線如圖12所示。a1點(diǎn)位于焊縫中心,位于Q235焊板上的b3、d3點(diǎn)距焊縫中心的距離分別與位于316L焊板上的點(diǎn)c3、e3相同。由于距離焊縫中心較遠(yuǎn),熱量主要來自于焊縫處局部熱源的熱傳導(dǎo),并且與空氣存在較大面積的對流散熱,b3、d3、c3、e3各點(diǎn)的溫度峰值遠(yuǎn)低于焊縫中心最高溫度,并且各點(diǎn)的溫度變化不是很劇烈。位于Q235焊板上的b3、d3節(jié)點(diǎn)溫度在相同時(shí)刻均高于316L焊板上的節(jié)點(diǎn)c3、e3;并且在焊接加熱階段,Q235焊板上節(jié)點(diǎn)溫度要先于316L開始上升,這主要是因?yàn)镼235熱導(dǎo)率及與空氣的對流散熱系數(shù)均比316L大,利于熱量的傳遞,溫度變化更迅速。
路徑2上點(diǎn)在第1道焊接過程中的熱循環(huán)曲線如圖11所示。其溫度變化總趨勢與路徑1上各點(diǎn)相似。路徑2上各點(diǎn)的溫度峰值均在1250℃以下低于母材熔點(diǎn),這是因?yàn)槁窂?上各點(diǎn)均在焊縫之外,是通過熱傳導(dǎo)將熔池?zé)崃總鬟f給各點(diǎn)的,并且由于先焊熔池的預(yù)熱作用,后焊區(qū)域節(jié)點(diǎn)的峰值溫度略高于先焊區(qū)域節(jié)點(diǎn)的峰值溫度。由于路徑2上各點(diǎn)到焊縫中心距離相等,因此所經(jīng)歷的熱循環(huán)過程變化規(guī)律基本相同,區(qū)別在于熔池加熱過程中各點(diǎn)溫度達(dá)到峰值的時(shí)間不同。
圖12 路徑3上各點(diǎn)熱循環(huán)曲線Fig.12 Thermal cycle curves of path 3
(1)運(yùn)用Ansys軟件模擬了316L/Q235異種鋼平板對接焊三維動(dòng)態(tài)焊接過程,得到了不同時(shí)刻焊接溫度場分布及焊件上不同節(jié)點(diǎn)的熱循環(huán)曲線。
(2)焊接溫度場高溫中心隨焊接熱源勻速向前移動(dòng),熱源前方的溫度梯度較大;后方溫度梯度比較緩和。Q235鋼較316L鋼熱傳導(dǎo)系數(shù)及對流散熱系數(shù)大,更利于熱量傳導(dǎo)和流動(dòng),Q235鋼焊板在焊接和冷卻過程中的溫度變化也更迅速。
(3)節(jié)點(diǎn)距焊縫中央熱源越近,熱源加熱作用越強(qiáng)烈,溫度升高越迅速且峰值越高;與焊縫中心等距離節(jié)點(diǎn)的熱循環(huán)曲線特征相似,區(qū)別是溫度達(dá)到峰值的時(shí)間先后有所不同;由于焊接熱源預(yù)熱作用,不同道次間焊接溫度隨著道次的增加而上升。
[1]李光福,李冠軍,方可偉,等.異材焊接件A508/52M/316L在高溫水環(huán)境中的應(yīng)力腐蝕破裂[J].金屬學(xué)報(bào),2011,47(07):797-803.
[2]Satyanarayanaa V V,Reddyb G M,Mohandas T.Dissimilar metal friction welding of austenitic-ferritic stainlesssteels[J].Journal of Materials Processing Technology,2005,160(2):128-137.
[3]劉中青,劉 凱.異種金屬焊接技術(shù)指南[M].北京:機(jī)械工業(yè)出版社,2000:30-80.
[4]宮大猛,雷 毅.數(shù)值模擬在焊接中的應(yīng)用分析[J].電焊機(jī),2012,42(6):58-62.
[5]劉 川,梁延德,鄧瑩瑩.超聲波塑料焊接溫度場有限元數(shù)值模擬[J].機(jī)械工程材料,2005,28(11):14-16.
[6]陳 浩,高傳玉,曹 敏,等.薄板鎂合金激光焊接溫度場分布的模擬研究[J].熱加工工藝,2012,41(3):154-159.
[7]遲露鑫,麻永林,邢淑清,等.壓力容器用鋼板焊接溫度場模擬與驗(yàn)證[J].焊接學(xué)報(bào),2012,33(5):61-65.
[8]肖 馮,米彩盈.T型角接頭焊接熱源模型研究[J].電焊機(jī),2010,40(6):41-45.
[9]董學(xué)偉,黎向鋒,左敦穩(wěn),等.7022鋁合金攪拌摩擦焊接全過程溫度場的數(shù)值模擬[J].機(jī)械工程材料,2012,36(10):92-96.