阮敬平 孫俊華 劉鵬

摘要：通過(guò)GAMBIT軟件建立穿孔型等離子弧焊焊接316不銹鋼的三維有限元模型并劃分網(wǎng)格，利用FLUENT軟件進(jìn)行迭代計(jì)算。通過(guò)編寫UDF程序?qū)崿F(xiàn)移動(dòng)熱源和力源的加載以及相的定義，并利用VOF方程追蹤相的界面，求解出橫截面、縱截面在不同時(shí)刻的溫度場(chǎng)云圖和熔池內(nèi)液態(tài)金屬的流動(dòng)情況。結(jié)果表明：隨著焊接時(shí)間的變化，工件在電弧熱的作用下逐漸被熔透，熔池內(nèi)液態(tài)金屬在等離子弧壓力和等離子流力及其反作用力的作用下發(fā)生動(dòng)態(tài)演變。對(duì)6 mm的316不銹鋼進(jìn)行焊接實(shí)驗(yàn)，并對(duì)比焊縫橫截面的數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果，兩者基本吻合，驗(yàn)證了數(shù)值模擬結(jié)果的可靠性。

關(guān)鍵詞：316不銹鋼;溫度場(chǎng);液態(tài)金屬;數(shù)值模擬

中圖分類號(hào)：TG409 ? ? ?文獻(xiàn)標(biāo)志碼：A ? ? ? ? 文章編號(hào)：1001-2003（2021）05-0056-05

DOI：10.7512/j.issn.1001-2303.2021.05.10

0 ? ?前言

等離子弧焊接（PAW）具有能量密度高、焊接速度快、應(yīng)力應(yīng)變小、焊縫成形美觀等優(yōu)點(diǎn)，特別適合焊接各種難熔、易氧化及熱敏感性強(qiáng)的金屬材料（如鎢、鉬、銅、鎳、鈦等）。相比于TIG、MIG焊，PAW在焊接中厚板時(shí)更具優(yōu)勢(shì);相比于激光焊、電子束焊等，等離子弧焊接具有成本低，對(duì)焊件要求低的優(yōu)點(diǎn)，是焊接中厚不銹鋼板常用的方法之一。

316不銹鋼是一種碳含量較低的奧氏體不銹鋼，由于它在海水或其他介質(zhì)中的耐腐蝕性能比其他不銹鋼更加優(yōu)異，因此常用于船舶制造、航空航天等領(lǐng)域。但316不銹鋼具有較低的熱導(dǎo)率、較高的電阻率以及較大的線膨脹系數(shù)，焊接難度較大[1-3]。目前，研究者對(duì)焊接316不銹鋼時(shí)的電弧特性、增材制造以及316與其他金屬材料的焊接等做了大量研究。尹玉祥[4]采用TIG往復(fù)焊對(duì)316不銹鋼進(jìn)行增材制造，得出了增材層數(shù)對(duì)電弧形態(tài)及電弧溫度的影響規(guī)律。顧偉[5]對(duì)316不銹鋼的焊接接頭組織和性能進(jìn)行研究，得到了均勻一致、質(zhì)量穩(wěn)定的焊接接頭。張俠洲[6]采用Q235R/316L復(fù)合基板，得到了焊接性能良好、抗拉強(qiáng)度高的焊接接頭。但針對(duì)等離子弧焊接316不銹鋼板的溫度場(chǎng)、流場(chǎng)的研究相對(duì)較少。

近年來(lái)，隨著數(shù)值模擬軟件SYSWELD、FLUENT等的使用，越來(lái)越多的研究者[7-11]通過(guò)數(shù)值模擬的方法來(lái)探究焊接過(guò)程中的溫度場(chǎng)和流場(chǎng)。武傳松教授團(tuán)隊(duì)[7-11]對(duì)等離子弧焊接304不銹鋼的溫度場(chǎng)、流場(chǎng)進(jìn)行數(shù)值模擬，采用VOF追蹤小孔界面，得到的數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果吻合，揭示焊接過(guò)程的內(nèi)在機(jī)理，并取得了豐碩的成果。文中采用FLUENT軟件分析316不銹鋼的溫度場(chǎng)和流場(chǎng)，并通過(guò)工藝實(shí)驗(yàn)驗(yàn)證了其正確性，為后續(xù)的工藝實(shí)驗(yàn)提供理論支撐。

1 焊接溫度場(chǎng)和流場(chǎng)的數(shù)值模擬

1.1 工藝實(shí)驗(yàn)

為了更好地了解等離子弧焊接316不銹鋼的工藝過(guò)程，驗(yàn)證數(shù)值模擬模型的正確性，焊接工藝實(shí)驗(yàn)采用以下參數(shù)：焊接電流150 A，電弧電壓19 V，焊接速度110 mm/min，保護(hù)氣流量18 L/min，等離子氣流量3 L/min，噴嘴距工件的高度為5 mm。將工件的橫截面拋光后用王水（濃鹽酸∶濃硝酸=3∶1）腐蝕，得到了焊縫橫截面熔合線形狀并進(jìn)行拍照。

1.2 建立有限元模型

采用前處理軟件GAMBIT建立三維有限元模型，模型尺寸為70 mm×20 mm×10 mm，包括中間6 mm的工件以及工件上下分別為2 mm的空氣層，有限元模型如圖1所示。焊接初始位置在坐標(biāo)原點(diǎn)位置，如圖1a所示，x軸為焊接速度方向，z軸為工件厚度方向，ABCD為volocity-inlet，AEDJ和BCGF為pressure-outlet，AMNB為symmetry，其余壁面為wall。在焊接過(guò)程中，由于焊縫區(qū)溫度梯度大，遠(yuǎn)離焊縫位置的溫度梯度小，因此為了節(jié)省計(jì)算資源，同時(shí)保證計(jì)算精度，采用非均勻網(wǎng)格，xOy面網(wǎng)格如圖1b所示。同時(shí)為了簡(jiǎn)化計(jì)算，假設(shè)導(dǎo)熱系數(shù)、粘度是溫度的函數(shù)，其余材料熱物性參數(shù)為常數(shù)。

1.3 “熱-力 ”模型的建立

等離子弧焊接涉及到復(fù)雜的傳熱、輻射、熔化與凝固等過(guò)程，因此建立合適的“熱-力”耦合模型是數(shù)值模擬的關(guān)鍵。根據(jù)等離子弧焊接熔池呈現(xiàn)倒喇叭狀的特點(diǎn)，采用雙橢球體+錐體熱源模型。雙橢球熱源方程為：

錐體熱源方程為：

式中 UPAW為電弧電壓;IPAW為焊接電流;η為焊接熱效率;vx為焊接速度;A前、 A后、B、C、re、ri、h1為熱源分布參數(shù);n上、n下 分別為雙橢球體熱源和錐體熱源的能量分配系數(shù)，且n上+n下=1。文中采用追蹤熔池最底部位置，實(shí)現(xiàn)對(duì)錐體熱源高度的實(shí)時(shí)控制。具體原理如圖2所示，h1=zi-C（其中zi為熱源總高度）。等離子弧焊接開始時(shí)只有上半部分熱源，此時(shí)n上=1，n下=0。隨著焊接過(guò)程的進(jìn)行，逐漸出現(xiàn)錐體熱源，此時(shí)設(shè)定n下=Dh1/L，n上=1-n下（式中，D為調(diào)節(jié)系數(shù)， L為工件厚度）。

在等離子弧焊接過(guò)程中，等離子弧壓力和等離子流力是形成小孔的重要因素。在焊接過(guò)程中，電弧力和電弧熱的分布趨勢(shì)相同，因此等離子弧壓力和等離子流力均采用雙橢圓分布模型，方程如下：

式中 ζ為等離子弧壓力調(diào)節(jié)系數(shù);μ0為真空磁導(dǎo)率;a11、a22、b11為壓力分布系數(shù)。

VOF（流體體積函數(shù)）是常用的追蹤穩(wěn)態(tài)或瞬態(tài)氣液界面的方法之一，通過(guò)計(jì)算每個(gè)單元的體積分?jǐn)?shù)來(lái)追蹤兩個(gè)或者多個(gè)互不相融的流體界面。VOF的計(jì)算必須基于壓力求解器，在每個(gè)六面體網(wǎng)格中，氣液兩相的體積分?jǐn)?shù)和都為1。VOF控制方程為：

式中 vx、vy、vz分別為x、y、z方向的速度。

1.3 邊界條件

在焊接起始時(shí)刻，即t=0時(shí)：

在工件上表面：

在工件下表面：

在工件對(duì)稱面：

式中 T∞為環(huán)境溫度;n1、n2 分別為上、下表面的法向量;q為等離子弧傳遞到工件的熱流密度;qd為熱對(duì)流損失的熱流密度;qf為熱輻射損失的熱流密度;qz為蒸發(fā)損失的熱流密度。

2 結(jié)果分析

2.1 焊接溫度場(chǎng)和流場(chǎng)

利用建立的“ 熱-力 ”模型對(duì)6 mm的316不銹鋼進(jìn)行等離子弧焊接的溫度場(chǎng)和流場(chǎng)數(shù)值模擬，得到不同時(shí)刻縱截面以及橫截面的溫度場(chǎng)和流場(chǎng)變化圖，如圖3所示。在焊接的初始階段，等離子弧溫度高達(dá)10 000～20 000 K，工件表面迅速被加熱熔化形成熔池。隨著時(shí)間的進(jìn)行，等離子弧的熱量在熔池及其周圍發(fā)生熱傳導(dǎo)和熱輻射，使熔池周圍的固態(tài)金屬溫度升高，熔池內(nèi)的液態(tài)金屬在等離子弧壓力和等離子流力的作用下被擠壓到熔池四周，新的固態(tài)金屬又被暴露出來(lái)，繼續(xù)被等離子弧加熱，直至焊透工件。隨著焊接的進(jìn)行，焊槍向前移動(dòng)，電弧后方的溫度場(chǎng)產(chǎn)生后拖，逐步呈現(xiàn)出前窄后寬的雙橢球體形狀。溫度場(chǎng)在t2時(shí)刻熔寬大幅度變窄，出現(xiàn)錐體形狀，最終在t4時(shí)刻熔深最大，最終溫度場(chǎng)出現(xiàn)“ 倒喇叭狀”的形狀，隨后，焊接處于穩(wěn)定階段，溫度場(chǎng)僅出現(xiàn)小范圍波動(dòng)。

從焊縫縱截面流場(chǎng)可以看出，在t1階段，由于等離子弧壓力和等離子流力的作用，熔池中的液態(tài)金屬受到向下的作用力，液態(tài)金屬被擠壓開始下凹，電弧正下方的液態(tài)金屬由熔池表面流向熔池底部，但此時(shí)熔池下方的金屬?zèng)]有熔化，仍然以固體的形式存在，因此在固體金屬的反作用力下，熔池底部的液態(tài)金屬被擠壓向四周流動(dòng)，一部分液態(tài)金屬在上表面凝固，形成余高。隨著焊接過(guò)程的進(jìn)行，在t2～t3時(shí)刻，熔池內(nèi)的液態(tài)金屬繼續(xù)在等離子弧壓力和等離子流力的作用下被排開，從熔池內(nèi)壁流向四周。當(dāng)熔深繼續(xù)增大，出現(xiàn)熔寬大幅度減小的錐體熱源的位置時(shí)，熔池底部金屬向下流動(dòng)，熔池內(nèi)壁金屬繼續(xù)被擠壓自下向上流動(dòng)，熔池外壁金屬向下流動(dòng)，并在反作用力的作用下，沿熔池外壁流出。在t4～t6時(shí)刻，小孔貫穿整個(gè)工件，此時(shí)等離子弧壓力和等離子流力的反作用力大幅度降低，因此熔池內(nèi)中下部液態(tài)金屬向下流動(dòng)，熔池偏上部液態(tài)金屬向上流動(dòng)。

從焊縫橫截面流場(chǎng)可以看出，不同時(shí)刻的液態(tài)金屬流動(dòng)情況與縱截面相似，但橫截面左右兩側(cè)流體流動(dòng)比較均勻，而縱截面在電弧后方的流體流速較快。這是因?yàn)樵诤附舆^(guò)程中，等離子弧向后偏移，作用在電弧后方使焊接后方的等離子弧壓力的反作用力較大，促使其流速變大。

2.2 實(shí)驗(yàn)驗(yàn)證

將實(shí)驗(yàn)所得橫截面的焊縫熔合線與達(dá)到準(zhǔn)穩(wěn)態(tài)后的模擬結(jié)果進(jìn)行對(duì)比，如圖4所示，兩者基本吻合，驗(yàn)證了數(shù)值模擬結(jié)果的可靠性，這對(duì)優(yōu)化焊接工藝和節(jié)約成本具有重大意義。

3 結(jié)論

（1）劃分了既能保證精度又能節(jié)約計(jì)算資源的非均勻網(wǎng)格，建立了三維瞬態(tài)“熱-力”耦合模型。

（2）對(duì)等離子弧焊接316不銹鋼的溫度場(chǎng)和流場(chǎng)進(jìn)行數(shù)值模擬，結(jié)果表明：等離子弧焊接316不銹鋼的溫度場(chǎng)呈現(xiàn)“倒喇叭狀”，流體流動(dòng)主要是取決于電弧壓力和等離子流力及其反作用力。

（3）進(jìn)行了等離子弧焊接316不銹鋼的工藝實(shí)驗(yàn)，經(jīng)腐蝕后觀察焊縫橫截面的熔合線，并將其與模擬結(jié)果進(jìn)行對(duì)比，兩者基本吻合，驗(yàn)證了模擬結(jié)果的正確性，這為優(yōu)化等離子弧焊接316不銹鋼的焊接工藝提供了理論支撐。

參考文獻(xiàn)：

黃彥博，劉智慧. Q235B+TP316L不銹鋼復(fù)合板焊接工藝研究[J].金屬加工（熱加工），2020（9）：19-21.

姚祥宏，周琦，王克鴻. 基于焊道尺寸控制的電弧增材成形高氮奧氏體不銹鋼與316L不銹鋼交織結(jié)構(gòu)[J].材料工程，2020，48（1）：54-60.

錢強(qiáng).材料及材料的焊接行為[M]. 哈爾濱：WTI中德焊接國(guó)際，2019.

尹玉祥，肖笑，邱然鋒. 316L不銹鋼TIG焊增材電弧物理特性分析[J].電焊機(jī)，2019，49（12）：52-56.

顧偉，皇甫樂森，趙勇. 316L不銹鋼厚壁管全位置TIG焊接頭組織及性能[J].焊接技術(shù)，2018，47（11）：14-16.

張俠洲，王鳳會(huì)，陳延清. Q345R/316L復(fù)合板鎳基焊接接頭性能[J].電焊機(jī)，2017，47（8）：81-83，113.

霍玉雙，武傳松，陳茂愛. 等離子弧焊接小孔形狀和穿孔過(guò)程的數(shù)值分析[J]. 金屬學(xué)報(bào)，2011，47（6）：706-712.

Wu C S，Wang H G，Zhang Y M. A new heat source model for keyhole plasma arc welding in FEM analysis of the temperature profile[J]. Welding Journal，2006，85（12）：284-289.

張濤，武傳松. 穿孔等離子弧焊接熱場(chǎng)和流場(chǎng)的數(shù)值模擬[J]. 焊接學(xué)報(bào)，2011，32（7）：87-90.

張濤，武傳松，陳茂愛.穿孔等離子弧焊接熔池流動(dòng)和傳熱過(guò)程的數(shù)值模擬[J]. 金屬學(xué)報(bào)，2012，48（9）：1025-1032.

Li Tianqing，Wu Chuansong. An analytic formula describing the plasma arc pressure distribution[J].China Welding，2014，23（2）：7-11.