激光焊接熔池動態(tài)行為數(shù)值模擬的進(jìn)展與分析

彭 進(jìn)，龍偉民，張永振

1.中機(jī)智能裝備創(chuàng)新研究院（寧波）有限公司，浙江 寧波 315700

2.河南科技大學(xué) 高端軸承摩擦學(xué)技術(shù)與應(yīng)用國家地方聯(lián)合工程實(shí)驗(yàn)室，河南 洛陽 471023

3.華北水利水電大學(xué)，河南 鄭州 450045

0 前言

大功率激光制造是一種以“光能源”和“光工具”作為材料加工手段的先進(jìn)技術(shù)，能夠很好地解決不同材料的加工、成形和精煉等問題，被稱為21世紀(jì)的綠色制造技術(shù)［1］。焊接是制造業(yè)的基礎(chǔ)工藝，在經(jīng)濟(jì)建設(shè)中發(fā)揮著重要作用［2］。激光焊接作為特種焊接，具有能量密度高、焊接變形小、焊接精度高、焊接接頭質(zhì)量與電子束焊接相差不大但不產(chǎn)生X射線等突出優(yōu)點(diǎn)［2-4］，廣泛應(yīng)用于航空航天、汽車制造、電力電子、船舶等領(lǐng)域［5-7］。

激光焊接常見的焊接缺陷有氣孔、咬邊、焊接變形等，而熔池的動態(tài)行為對焊接缺陷的形成有直接影響。對于熔池內(nèi)部的動態(tài)行為，國內(nèi)外主要采用實(shí)時監(jiān)測和數(shù)值模擬的方法進(jìn)行研究和分析。日本大阪大學(xué)在實(shí)時監(jiān)測熔池內(nèi)部動態(tài)行為研究方面具有代表性。隨著計(jì)算機(jī)技術(shù)的發(fā)展，激光焊接熔池行為的數(shù)值模擬技術(shù)得到了長足發(fā)展。通過數(shù)值模擬，可提前預(yù)測并采取相應(yīng)措施減少焊接缺陷的形成，這將有利于提高激光焊接工件的焊縫質(zhì)量。

文中對國內(nèi)外激光焊接熔池行為的數(shù)值模擬技術(shù)，如激光深熔焊、雙束激光焊接、真空激光焊接、磁場作用下的激光焊接、液態(tài)金屬作用下的激光焊接熔池動態(tài)行為數(shù)值模擬的研究現(xiàn)狀進(jìn)行了分析和概括。

1 激光深熔焊熔池動態(tài)行為數(shù)值模擬

激光深熔焊接以在熱源下方熔池內(nèi)部存在匙孔為主要特征，也被稱為匙孔模式激光焊接，整個焊接過程涉及的物理現(xiàn)象包括：激光材料交互、熱傳導(dǎo)、熔化、蒸發(fā)、羽輝形成、流體流動和自由表面形成等。在激光深熔焊熔池行為的數(shù)值模擬方面，Han等［8］通過建立的激光焊接模型分析了焊接驅(qū)動力的影響，發(fā)現(xiàn)反沖壓力在焊接過程中有明顯作用，當(dāng)不加載反沖壓力時，熔池中無法生成匙孔。Zhang等［9］建立了10 mm厚度鋼板激光全熔透焊接模型，發(fā)現(xiàn)在全熔透焊接中，熔池下表面可能比上表面更長、更不穩(wěn)定。全熔透過程中，激光的一部分能量通過匙孔直接從底部射出而不與液態(tài)金屬接觸，這與非熔透焊的激光能量作用機(jī)制有很大區(qū)別，熔池瞬態(tài)流場是反沖壓力與表面張力共同作用的結(jié)果，熔池波動周期約8 ms。

Rai［10-13］建立了激光焊接熱流耦合模型，設(shè)定熔池表面為平面，匙孔處于準(zhǔn)穩(wěn)態(tài)，匙孔壁面的溫度為母材的沸點(diǎn)溫度。指出由于表面張力與溫度有直接關(guān)系，熔池表面的溫度梯度越大則表面張力梯度越大，熔池表面的流動方向?yàn)橛沙卓组_口處向熔池邊緣流動，進(jìn)而增大了熔池寬度和熔池長度。Ye等［14］建立了激光全熔透焊接傳熱和流動的三維模型，發(fā)現(xiàn)除了焊接速度外，Marangoni對流對工件溫度分布和熔池熔體流動也起著至關(guān)重要的作用，即使是薄板的全熔透焊接也不容忽視。Ribic等［15］研究了表面活性元素在激光深熔焊中的作用，發(fā)現(xiàn)表面活性元素會影響表面張力梯度的大小和方向，進(jìn)而引起Marangoni對流的變化，當(dāng)表面活性元素含量增大到一定程度，則會改變?nèi)鄢乇砻娴牧鲃臃较颍ㄈ鄢匾簯B(tài)金屬由熔池邊緣流向匙孔），進(jìn)而縮小熔池的尺寸。

Chang等［16］建立了鈦合金全熔透激光焊接模型，發(fā)現(xiàn)對于較低功率、較低速度的焊接過程，匙孔后面的流體流動是湍流和不穩(wěn)定的，形成渦流；對于更高功率、更高速度的焊接過程，流體流動湍流特性減弱，并且不會形成渦流。Sohail等［17］研究了激光焊接20 mm厚鋼板時熱輸入對熔池流場的影響，發(fā)現(xiàn)焊接速度會影響熔池尺寸，激光功率越高熔池渦流越明顯，這是因?yàn)樵谳^高的熱輸入作用下，匙孔吸收的激光能量增多，熔池的體積增多，熔池底部的流體流動趨勢增強(qiáng)。

Wang等［18］通過建立的熱源模型發(fā)現(xiàn)渦流存在于熔池頂部和底部表面，渦流的存在有利于熱傳遞和熔池的形成，匙孔的形貌與激光功率、焊接速度有關(guān)。當(dāng)焊接速度增加，匙孔前壁傾斜角度增大。Geiger等［19］建立的激光焊接模型中考慮了菲涅爾吸收、蒸氣壓力、表面張力、熔化和蒸發(fā)焓以及由于材料蒸發(fā)造成的能量損失，發(fā)現(xiàn)激光焊接鍍鋅板過程中，液態(tài)金屬沿匙孔前壁向下流動，匙孔的直徑尺寸出現(xiàn)周期性的變化。熔池中流體最高速度（約為0.5 m/s）出現(xiàn)在匙孔前壁的區(qū)域，在匙孔壁后部，熔體流動速度顯著降低。Zhao等［20］建立的模型考慮了等離子體、液態(tài)金屬、固態(tài)金屬的三種相，發(fā)現(xiàn)匙孔壁前方的液態(tài)金屬通過匙孔底部流向匙孔后方，熔池頂部的液態(tài)金屬從高溫區(qū)域流向低溫區(qū)域，隨著焊接時間的進(jìn)行匙孔是不斷波動的，當(dāng)匙孔深度突然降低，則在熔池底部形成焊接氣泡，這也是氣孔缺陷形成的主要原因，隨著激光持續(xù)照射進(jìn)匙孔，匙孔深度變大。

Li［21］和 Lu［22］等在模型中引入了凝固模型和氣泡模型，研究了T型接頭匙孔型氣孔的形成過程，發(fā)現(xiàn)匙孔的不穩(wěn)定性和熔池流動的復(fù)雜性導(dǎo)致了匙孔的倒塌以及氣泡的形成。熔池中不同部位的氣泡的移動軌跡不同，氣泡的移動主要受周圍流體及浮力的影響，一些氣泡移動過程中與其他氣泡合并。T型接頭的間隙由0 mm變?yōu)?.6 mm，焊縫氣孔率增大，當(dāng)T型接頭的間隙變?yōu)? mm時，焊縫氣孔率較低，這主要是由于在間隙附近匙孔保持張開的狀態(tài)。

Amara［23］、Zhang［24］等在模型中引入惰性氣體噴吹在熔池上方，模型考慮了熔池金屬蒸汽與液態(tài)金屬的相互作用以及熔池液態(tài)金屬與空氣的相互作用，通過在熔池上方吹入惰性氣體可提高熔池金屬的流動行為進(jìn)而得到質(zhì)量較好的焊接接頭。噴射的氣流不僅推動熔融金屬向熔池后方流動，同時在熔池表面形成適當(dāng)?shù)膲毫τ兄诰S持匙孔上部出口周圍的凹面。在噴射氣流的情況下推動匙孔開口上方周圍隆起的液態(tài)金屬向熔池后方流動，有利于擴(kuò)大匙孔開口尺寸，維持匙孔的穩(wěn)定性。

通過對以上文獻(xiàn)的分析可知，對于激光深熔焊熔池動態(tài)行為的數(shù)值模擬，目前主要集中在對焊接驅(qū)動力（表面張力、反沖壓力等）、焊接參數(shù)（焊接速度、激光功率等）、外部條件（表面活性元素、惰性氣體噴吹在熔池上方）對熔池流動、匙孔穩(wěn)定性的影響，以及匙孔失穩(wěn)與焊縫氣孔率的關(guān)系。研究相對全面，有利于后期激光深熔焊接實(shí)驗(yàn)工作的開展。

2 雙束激光焊接熔池動態(tài)行為數(shù)值模擬

在雙束激光焊接熔池行為的數(shù)值模擬方面，Pang等［25］建立了雙光束激光焊接模型，發(fā)現(xiàn)當(dāng)以串行排布模式進(jìn)行雙束激光焊接時，沿焊接方向后一束激光形成的匙孔深度大于前一束激光形成的匙孔深度，這主要是由于前束激光的預(yù)熱效應(yīng)造成的。在相同的焊接熱輸入條件下，與單束激光焊接相比，雙束串行激光焊接熔池會產(chǎn)生較小的湍流流動。陳偉東等［26］建立了并行雙束激光焊接模型，并對其匙孔行為進(jìn)行研究，發(fā)現(xiàn)提高焊接速度，匙孔深度減小，同時匙孔形貌的波動行為減弱；隨著雙光束激光斑點(diǎn)距離的增大，匙孔深度減小，匙孔開口尺寸增大，當(dāng)雙束激光光斑距離過大時，則在熔池中出現(xiàn)兩個獨(dú)立的匙孔，這對焊接過程也會產(chǎn)生很大程度的影響。

Yang等［27-30］建立了T型接頭雙側(cè)激光焊接模型，發(fā)現(xiàn)T型接頭雙側(cè)激光焊接過程中存在能量耦合的情況，熔池最大長度位于雙束激光的交匯部位，匙孔底部出現(xiàn)渦流流動行為，而且在匙孔上方存在對流。T型接頭雙側(cè)激光束是否貫通以及是否同步焊接對熔池流動行為影響較大，這將影響焊接氣泡的逃逸行為。當(dāng)保持T型接頭雙側(cè)匙孔貫通、同步焊接時有利于焊接氣泡逃逸出熔池，提高焊接速度同樣有利于焊接氣泡逃逸出熔池。

通過對以上雙束激光焊接熔池動態(tài)行為數(shù)值模擬的分析，發(fā)現(xiàn)與單束激光焊接相比，雙束激光焊接有利于提高匙孔、熔池的穩(wěn)定性。筆者認(rèn)為，考慮到雙束激光焊接的優(yōu)勢，后期可以將雙束激光焊接模型引入單束激光焊接過程容易產(chǎn)生焊接缺陷的“特殊工件”中（如異型接頭等），以期揭示其提高焊接穩(wěn)定性的機(jī)理。

3 真空激光焊接熔池動態(tài)行為數(shù)值模擬

在真空激光焊接熔池動態(tài)行為數(shù)值模擬方面，Pang及Fabbro等［31-33］建立了真空環(huán)境激光焊接熱源模型，研究發(fā)現(xiàn)相同工藝參數(shù)下真空環(huán)境激光焊接熔深明顯大于大氣環(huán)境下的激光焊接熔深，這主要是真空環(huán)境下材料的沸點(diǎn)溫度降低導(dǎo)致的。隨著焊接速度的增加，匙孔前壁平均溫度逐漸升高，而且與大氣環(huán)境相比，在真空環(huán)境下匙孔前壁溫度增加的幅度明顯較大。

Tang等［34］建立了真空環(huán)境鋁合金激光焊接的三維數(shù)學(xué)模型，發(fā)現(xiàn)在較低的環(huán)境壓力下，激光能量在匙孔壁上分布更均勻，產(chǎn)生氣泡的可能性較小。環(huán)境壓力也會影響熔池的形狀和大小以及氣泡移動的距離，進(jìn)一步影響氣泡能否逃逸出熔池。

Li等［35］建立了真空環(huán)境激光焊接數(shù)學(xué)模型，發(fā)現(xiàn)在真空環(huán)境下，匙孔的尺寸變大而且匙孔壁面的凸起現(xiàn)象減弱，熔池內(nèi)部漩渦流動現(xiàn)象減弱，匙孔壁的液態(tài)金屬流動速度變快，這是由于真空環(huán)境下母材金屬蒸發(fā)速率提升導(dǎo)致的。

通過對以上文獻(xiàn)的分析可知，真空環(huán)境激光焊接可提高匙孔的穩(wěn)定性、改變?nèi)鄢氐牧鲃?。但模型中沒有考慮匙孔內(nèi)部金屬蒸汽的可壓縮動力學(xué)，以及其與匙孔內(nèi)壁間相互作用。在后期的研究中，可編制相應(yīng)的程序?qū)霐?shù)值模擬軟件中，研究匙孔內(nèi)部金屬蒸汽與匙孔內(nèi)壁間相互作用機(jī)理。

4 磁場作用下的激光焊接熔池動態(tài)行為數(shù)值模擬

在磁場作用下的激光焊接熔池動態(tài)行為數(shù)值模擬方面，Zhou等［36］建立了磁場作用下的脈沖激光焊接模型，研究電磁場對瞬態(tài)熔池金屬流動、匙孔動態(tài)行為、氣孔形成的影響。研究發(fā)現(xiàn)，氣孔的形成與液態(tài)金屬的凝固速度、回填速度有關(guān)。如果匙孔上方的液態(tài)金屬回填速度不夠快，不能及時到達(dá)匙孔底部，則容易產(chǎn)生氣孔缺陷。在激光焊接過程中，通過施加電磁場可在匙孔倒塌時提高液態(tài)金屬的回填速度，進(jìn)而降低焊縫氣孔缺陷。

Chen等［37-40］研究了磁場作用下的鋁合金激光全熔透的熱電電流和熱電磁效應(yīng)，發(fā)現(xiàn)熱電電流密度在熔池達(dá)到穩(wěn)定前處于持續(xù)波動狀態(tài)。在磁場作用下，隨著磁通密度從0 mT增加到238 mT，鋁合金焊縫的成形質(zhì)量得到提高，尤其是焊縫表面寬度明顯減小。當(dāng)外加的磁場由200 mT增至500 mT時，將減弱熱毛細(xì)對流效應(yīng)，有利于改變焊縫形貌。

Bachmann等［41-45］建立了電磁場作用下的激光焊接三維熱流耦合模型，模型中考慮了熱毛細(xì)對流、自然對流、相變潛熱的情況，發(fā)現(xiàn)在鋁合金激光焊接過程中，沿焊接方向在橫向磁場作用下，熔池底部的液態(tài)金屬流動速度明顯降低，并且由于熔池流體的作用使熱量沿焊縫橫向傳遞，這樣焊縫呈現(xiàn)V形的趨勢增強(qiáng)。

綜上可知，磁場能夠影響熔池流動行為、抑制缺陷，進(jìn)而提升焊縫性能。后期可通過編制程序建立宏微觀耦合的多尺度仿真模型，模擬磁場對全尺寸熔池凝固過程行為的影響，進(jìn)一步揭示磁場影響下的熔池流動、組織演變的相關(guān)性。

5 液態(tài)金屬作用下的激光焊接熔池動態(tài)行為數(shù)值模擬

在液態(tài)金屬作用下的激光焊接熔池動態(tài)行為數(shù)值模擬方面，Rao等［46］指出，激光-GMAW復(fù)合焊接時，當(dāng)焊絲熔化的液態(tài)熔滴滴入熔池時，熔滴所攜帶的熱能和動能進(jìn)入熔池中，會造成熔滴進(jìn)入熔池的接觸面區(qū)域溫度升高，熔池獲得了額外的熱量，其凝固速度將降低。Gao等［47］建立了激光-MIG復(fù)合焊接模型，考慮了反沖壓力、電磁力等焊接作用力的影響，發(fā)現(xiàn)激光焊接反沖壓力以及焊絲熔化后的熔滴滴入熔池均會造成焊縫熔深的增加。

Wu等［48］通過建立的激光-MIG復(fù)合焊接模型，發(fā)現(xiàn)復(fù)合焊接過程中匙孔的形成分為三個階段：初始狀態(tài)的穩(wěn)定階段、匙孔快速擴(kuò)展階段以及熔池匙孔振蕩/動態(tài)穩(wěn)定階段。Cho、Na等［49-50］建立了激光-GMA復(fù)合焊接三維模型，引入了電磁力、表面張力以及電弧力，通過光線追蹤技術(shù)實(shí)現(xiàn)了激光能量在匙孔壁面的多次反射。研究發(fā)現(xiàn)，在焊絲熔化后以液態(tài)熔滴方式進(jìn)入熔池過程中，匙孔出現(xiàn)了張開和閉合情況，在反沖壓力的作用下匙孔底部的金屬向下流動，當(dāng)流動到熔池底部的固/液界面時，則向熔池后方流動，在浮力作用下產(chǎn)生向上的流動趨勢，最終在匙孔后方產(chǎn)生一個流動漩渦，此時熔池內(nèi)部流動行為較為復(fù)雜。Cho、Na等［51］建立了激光-GMA復(fù)合焊接合金元素分布模型，發(fā)現(xiàn)焊縫上部的Cr/Ni含量較多，且比較均勻，這是由于在匙孔后方存在順時針的流動漩渦，流動漩渦的作用有利于合金元素的均勻分布。

Peng等［52］建立了熔滴填充作用下的激光焊接熱流耦合模型，研究發(fā)現(xiàn)，當(dāng)熔滴填充進(jìn)入熔池時，由于熔滴距離匙孔較近容易造成對匙孔的沖擊，影響匙孔的穩(wěn)定性。熔滴填充也會對熔池流動速度產(chǎn)生影響，熔滴進(jìn)入熔池過程中匙孔底部的流體流動速度變大。

通過對以上文獻(xiàn)的分析可知，液態(tài)金屬進(jìn)入熔池會對匙孔、熔池造成一定程度的沖擊，影響焊接的穩(wěn)定性。但液態(tài)金屬進(jìn)入熔池并與其進(jìn)行傳熱傳質(zhì)，同樣會影響合金元素在焊縫的分布。后期可編制相應(yīng)程序，與合金元素分布和組織演變相結(jié)合，進(jìn)一步揭示對焊縫質(zhì)量影響的機(jī)理。

6 結(jié)論

（1）激光深熔焊接過程中，反沖壓力、由表面張力梯度引起的Marangoni對流對熔池的動態(tài)行為影響較大。焊接過程中，匙孔處于實(shí)時波動狀態(tài)，當(dāng)匙孔倒塌后，則產(chǎn)生焊接氣泡。在熔池上方噴吹惰性氣體可提高匙孔以及熔池流場的穩(wěn)定性。

（2）相比單束激光焊接，雙束串行激光焊接熔池產(chǎn)生較小的湍流流動，提高焊接速度，匙孔深度減小，同時匙孔形貌的波動行為減弱。T型接頭雙側(cè)匙孔貫通、同步焊接時有利于焊接氣泡逃逸出熔池。

（3）真空激光焊接可提高熔池流場的穩(wěn)定性，熔池內(nèi)部漩渦流動現(xiàn)象減弱，產(chǎn)生氣泡的可能性較小。

（4）通過施加電磁場可在匙孔倒塌時提高液態(tài)金屬的回填速度，進(jìn)而降低焊縫氣孔缺陷。磁場將減弱熱毛細(xì)對流效應(yīng)，有利于改善焊縫成形質(zhì)量。

（5）焊絲熔化后的熔滴滴入熔池會增加焊縫熔深、降低熔池凝固速度。同時也會影響熔池的流場，造成對匙孔的沖擊，匙孔易出現(xiàn)閉合情況。