楊子酉,房玉超,何景山

(哈爾濱工業(yè)大學(xué) 先進(jìn)焊接與連接國家重點(diǎn)實(shí)驗(yàn)室，哈爾濱 150001)

電子束焊接因?yàn)榫哂心芰棵芏雀?、焊接速度快、焊后變形小、焊縫質(zhì)量高等優(yōu)點(diǎn)，被廣泛應(yīng)用于航空、航天、核能、交通運(yùn)輸及國防等高新技術(shù)制造領(lǐng)域[1-2]。如何得到高質(zhì)量的電子束焊縫，成為人們關(guān)注的焦點(diǎn)。隨著計(jì)算機(jī)能力的提升，數(shù)值分析逐漸成為人們研究熔池內(nèi)部溫度與流動(dòng)情況的重要手段。

選擇合適的焊接熱源模型是進(jìn)行高效準(zhǔn)確的數(shù)值模擬的關(guān)鍵。由于高能束焊接存在匙孔效應(yīng)，人們多采用雙橢球熱源模型[3]、高斯圓柱熱源模型[4]、高斯旋轉(zhuǎn)體熱源模型[5]和雙橢球與圓柱組合熱源模型[6]等。雖然這些體熱源和組合熱源可以在一定程度上反映高能束在空間中的分布規(guī)律，但不能真正體現(xiàn)高能束與匙孔壁面相互作用的動(dòng)態(tài)過程。近年來，高能束焊接熔池行為的研究取得了較大進(jìn)展，主要包括：熔池溫度場和流場耦合行為的研究[7-8]，匙孔形成過程的研究[9-10]，熔池驅(qū)動(dòng)力的研究[11-13]，焊接缺陷的研究[14-15]等。但很少有研究者對(duì)焊接進(jìn)入準(zhǔn)穩(wěn)態(tài)后熔池與匙孔的波動(dòng)與震蕩情況進(jìn)行分析。

本研究基于Fluent軟件，建立了電子束焊接三維瞬態(tài)模型，對(duì)非穿透焊接熔池進(jìn)行了數(shù)值模擬。分析了電子束焊接進(jìn)入準(zhǔn)穩(wěn)態(tài)后旋渦對(duì)熔池波動(dòng)的影響，并采用動(dòng)態(tài)面熱源首次將焊接熱源與匙孔形狀進(jìn)行耦合分析，總結(jié)出電子束焊接熔池的波動(dòng)規(guī)律。本研究可以為深入了解電子束焊接熔池行為提供前提和基礎(chǔ)。

1 有限元模型的建立

1.1 基本假設(shè)與控制方程

電子束焊接熔池流動(dòng)非常復(fù)雜，計(jì)算時(shí)必須進(jìn)行簡化，故提出以下假設(shè)：(1)熔池中液態(tài)金屬為不可壓縮牛頓流體，其流動(dòng)為層流；(2)忽略金屬蒸汽和匙孔壁面之間的摩擦剪切力；(3)金屬材料為各向同性，除比熱容、熱導(dǎo)率為溫度的函數(shù)之外，其余材料參數(shù)皆為常數(shù)；(4)試件初始溫度為300K。

在數(shù)值模擬中，焊接熔池要滿足計(jì)算流體力學(xué)三大守恒方程，即連續(xù)性方程、動(dòng)量守恒方程和能量守恒方程?？紤]到電子束焊接的物理過程，在動(dòng)量守恒方程中以表面張力、熱浮力、重力和金屬蒸汽反沖壓力作為源項(xiàng)，在能量方程中以電子束熱源、熔化焓和汽化焓作為源項(xiàng)。

1.2 邊界條件與數(shù)學(xué)模型

圖1為電子束焊接傳熱方式及受力過程的示意圖。對(duì)于熔池而言，金屬受高速電子束沖擊而受熱熔化，熔池表面金屬迅速蒸發(fā)，以蒸發(fā)和輻射的形式進(jìn)行散熱；液態(tài)金屬快速流動(dòng)，在熔池中以對(duì)流的方式進(jìn)行傳熱。對(duì)于未熔化的金屬而言，主要以輻射方式進(jìn)行散熱，因?yàn)殡娮邮附釉谡婵障逻M(jìn)行，所以不考慮工件表面對(duì)流換熱。

電子束焊接過程中會(huì)產(chǎn)生匙孔效應(yīng)，在匙孔壁面法線方向上，液態(tài)金屬主要受反沖壓力、液態(tài)金屬靜壓力以及彎曲液面產(chǎn)生的附加壓力作用；在匙孔壁面切線方向上，液態(tài)金屬主要受表面張力所產(chǎn)生的切向力以及金屬蒸汽與匙孔壁面摩擦產(chǎn)生的剪切力作用。在數(shù)值計(jì)算中，利用VOF算法追蹤焊接熔池的自由表面。

圖1 電子束焊接示意圖Fig.1 Sketch map of electron beam welding

在匙孔壁面上，電子束焊接熱源采用高斯面熱源，金屬蒸汽反沖壓力采用Semak提出的反沖壓力模型[16]，它們分別為：

(1)

(2)

式中：qr是焊接熱源能量密度；Q為電子束功率；R0是電子束聚焦半徑；(x2+y2)為熔池中任意點(diǎn)距離電子束中心的徑向距離；Pr為反沖壓力；p0為標(biāo)準(zhǔn)大氣壓；ΔHv為汽化潛熱；Tv為材料汽化溫度；R為理想氣體常數(shù)；T為熱力學(xué)溫度。

1.3 網(wǎng)格劃分與材料屬性

焊接熔池是熱力耦合的結(jié)果，而金屬的熱物性參數(shù)對(duì)焊接溫度場和熔池驅(qū)動(dòng)力有著重要的影響[17]，因此材料熱物性參數(shù)的選取十分重要。本研究采用10mm厚的2219鋁合金進(jìn)行電子束焊接，材料參數(shù)取自文獻(xiàn)[18]。電子束加速電壓60kV，電子束流15mA，焊接速率10mm/s，焊接方向沿x軸自左向右。模型尺寸為30mm×15mm×16mm，其中在金屬板上、下各設(shè)置3mm的自由液面變化區(qū)域。利用ICEM對(duì)模型進(jìn)行網(wǎng)格劃分，在近焊縫區(qū)采用加密網(wǎng)格，在遠(yuǎn)離焊縫區(qū)采用疏松網(wǎng)格，最小網(wǎng)格尺寸為0.125mm，網(wǎng)格總數(shù)為771,120，節(jié)點(diǎn)總數(shù)為801,477，計(jì)算區(qū)域外圍采用outflow邊界條件。計(jì)算域初始溫度300K，環(huán)境壓力為0Pa。

2 電子束焊接熔池的周期性波動(dòng)

2.1 焊接熔池溫度場和流場耦合過程的分析

圖2是電子束焊接過程中熔池縱截面和橫截面溫度場與流場隨時(shí)間演變圖，其中S和R分別代表表面張力和反沖壓力所引起的流動(dòng)。焊接初始階段，母材受電子束轟擊而迅速熔化形成熔池，液態(tài)金屬溫度繼續(xù)升高到沸點(diǎn)而蒸發(fā)為金屬蒸汽，其產(chǎn)生的反沖壓力使熔池表面下凹，熔池產(chǎn)生R1和R2方向的流動(dòng)，如圖2(a)所示。從圖2(b)，(c)中可以看出，隨著匙孔深度的增加，匙孔底部的液態(tài)金屬被金屬蒸汽反沖壓力向下壓，形成R1方向的流動(dòng)；隨著匙孔深度和壁面面積的增加，匙孔溫度上低下高，此時(shí)，匙孔壁面附近的液態(tài)金屬沿著匙孔壁面向S1方向流動(dòng)；在位于熔池表面的負(fù)表面張力梯度作用下，匙孔中的液態(tài)金屬向后排，受到固液界面的阻礙而回流，形成S2方向的流動(dòng)。在匙孔深度增加的過程中，匙孔底部的金屬蒸汽反沖壓力使液態(tài)金屬向R1方向流動(dòng)；R1方向上的液態(tài)金屬受到固液界面的阻礙而改變方向，然后被匙孔壁面的表面張力牽引形成S1方向的流動(dòng)； S2方向上和S1方向上的液態(tài)金屬在熔池中部匯合，共同流向匙孔壁面，如圖2(d)，(e)所示。

如圖2(f)所示，隨著焊接過程的進(jìn)行，匙孔深度達(dá)到準(zhǔn)穩(wěn)態(tài)，但熔池體積仍在持續(xù)增大、熔池拖尾長度也在不斷增加。在匙孔沿著焊接方向移動(dòng)的過程中，不停地產(chǎn)生被匙孔排開的液態(tài)金屬，由于匙孔的產(chǎn)生是一個(gè)由淺入深的過程，所以將被排開的液態(tài)金屬可以分為兩個(gè)部分：對(duì)于匙孔上半部分，液態(tài)金屬向四周排開，產(chǎn)生S2方向的流動(dòng)，流動(dòng)到熔池邊界的液態(tài)金屬回流，并在匙孔壁面表面張力的牽引下形成S3方向的流動(dòng)；對(duì)于匙孔下半部分，液態(tài)金屬被金屬蒸汽向下壓，最終產(chǎn)生流經(jīng)固液界面、最終匯入匙孔表面的S1方向的流動(dòng)。S1和S3方向上的液態(tài)金屬在熔池中間部位匯合，流向匙孔表面，緊接著再產(chǎn)生上述流動(dòng)，周而復(fù)始。

0.8s之后，電子束熱源繼續(xù)移動(dòng)，熔池流動(dòng)行為與0.8s基本保持一致，改變的僅僅是熔池體積和熔池拖尾長度。圖3是電子束焊接過程中熔池長度和寬度隨時(shí)間變化曲線。從圖3可以看出：焊接0.5s以后，熔池寬度趨于穩(wěn)定；焊接1s以后，熔池拖尾長度趨于穩(wěn)定，可以認(rèn)為焊接過程從此刻開始進(jìn)入準(zhǔn)穩(wěn)態(tài)。

圖3 熔池尺寸隨時(shí)間變化的模擬曲線Fig.3 Simulation curves of molten pool size vs time

2.2 焊接熔池中渦流的變化規(guī)律及其產(chǎn)生原因

熔池進(jìn)入準(zhǔn)穩(wěn)態(tài)后，盡管熔池的尺寸和形狀基本保持不變，但其仍然在小幅度震蕩，說明熔池內(nèi)部仍然處于一個(gè)周期性波動(dòng)的狀態(tài)。為了更好地研究熔池波動(dòng)情況，提取焊接進(jìn)入準(zhǔn)穩(wěn)態(tài)后，不同時(shí)刻熔池中的流線，如圖4所示。從圖4可以看出，匙孔前方熔池中的液態(tài)金屬流向匙孔后方熔池，為匙孔后方熔池提供了液態(tài)金屬。不難發(fā)現(xiàn)，匙孔后方熔池中存在渦流，這些渦流的大小和位置并不是固定不變的。

焊接進(jìn)入準(zhǔn)穩(wěn)態(tài)后，在熔池上方的逆時(shí)針渦流并排存在，如圖4(a)～(e)所示，這些渦流在熔池表面引起了波浪，波浪前進(jìn)的方向?yàn)檠厝鄢乇砻嫦蚝?。從圖4(a)，(e)中不難看出，熔池下方存在一個(gè)較大的順時(shí)針渦流。這個(gè)渦流在液態(tài)金屬的牽引下，逐漸移動(dòng)到匙孔壁面，如圖4(b)，(f)所示。渦流有繼續(xù)向匙孔移動(dòng)的趨勢，但受到匙孔壁面的阻礙而逐漸變??；渦流在撞擊匙孔壁面的過程中，由于順時(shí)針旋轉(zhuǎn)，導(dǎo)致渦流中的液態(tài)金屬流向匙孔壁面，渦流使匙孔壁面某一位置產(chǎn)生了一個(gè)大的凸起，阻礙了電子束與匙孔底部的直接作用，此時(shí)的匙孔處于坍塌的臨界狀態(tài)，如圖4(c)，(g)所示。匙孔表面凸起的部位與電子束直接發(fā)生作用，受到向下的金屬蒸汽反沖壓力，并在表面張力的驅(qū)使下坍塌，如圖4(d)，(h)所示。通過分析可知，熔池下方經(jīng)歷了兩個(gè)循環(huán)：(a)～(d)和(e)～(h)。在兩個(gè)循環(huán)中，熔池下方的渦流經(jīng)歷了4個(gè)過程：渦流的形成、渦流的移動(dòng)、渦流與匙孔相互作用、渦流的破壞。

圖4 熔池縱截面流線的演變過程(a)t=1.505s；(b)t=1.524s；(c)t=1.528s；(d)t=1.530s；(e)t=1.534s；(f)t=1.561s；(g)t=1.582s；(h)t=1.584sFig.4 Calculated streamlines from molten pool in longitudinal-section side views(a)t=1.505s；(b)t=1.524s；(c)t=1.528s；(d)t=1.530s；(e)t=1.534s；(f)t=1.561s；(g)t=1.582s；(h)t=1.584s

依照上述分析，可以根據(jù)熔池中液態(tài)金屬的流動(dòng)把熔池分為3個(gè)區(qū)域：匙孔前方(Zone Ⅰ)、匙孔后方上半部分(Zone Ⅱ)和匙孔后方下半部分(Zone Ⅲ)，具體分界線如圖5所示。在熔池的3個(gè)不同區(qū)域中，液態(tài)金屬呈現(xiàn)出不一樣的流動(dòng)情況，渦流的旋轉(zhuǎn)方向和移動(dòng)方向也不盡相同。通過文獻(xiàn)[19]可知，當(dāng)液體流動(dòng)方向和流速大小迅速發(fā)生改變時(shí)，就有可能產(chǎn)生渦流，下面根據(jù)此理論對(duì)這些不同區(qū)域渦流的產(chǎn)生原因進(jìn)行分析。

圖5 熔池不同分區(qū)示意圖Fig.5 Diagram of different zones of molten pool

區(qū)域Ⅰ是匙孔前方的熔池，由于電子束熱源剛剛與之接觸，因此該區(qū)域體積較小。電子束更多地與匙孔前壁面接觸，以產(chǎn)生新的匙孔。匙孔前壁面上的液態(tài)金屬受較大的反沖壓力，使液態(tài)金屬向下流動(dòng)，越過匙孔底端，最終流向熔池后方，如圖6(a)中的R方向。該區(qū)域的存在為匙孔后方熔池體積的穩(wěn)定起到了重要作用。

在區(qū)域Ⅱ中，靠近熔池表面會(huì)有多個(gè)逆時(shí)針渦流產(chǎn)生，其形成過程示意圖如圖6(b)所示。該區(qū)域靠近熔池表面的液態(tài)金屬在表面張力的作用下向V1方向流動(dòng)，因受到固液界面的阻礙而回流(V2)；一部分液態(tài)金屬在改變流動(dòng)方向的過程中與沿V1方向遷移的液態(tài)金屬相遇，另一部分流向S方向；沿V2流動(dòng)的液態(tài)金屬經(jīng)歷了固液界面的剪切摩擦后流速變小，小于V1方向的液態(tài)金屬流速，因此會(huì)形成逆時(shí)針渦流。S方向上的液態(tài)金屬一部分撞擊了匙孔壁面而改變方向，向匙孔開口處流動(dòng)，緊接著撞擊熔池上表面而改變方向(V3)，匯入V1。液態(tài)金屬流動(dòng)方向的連續(xù)改變，形成了一個(gè)逆時(shí)針渦流。兩個(gè)渦流在液態(tài)金屬的帶動(dòng)下不斷移動(dòng)。

在區(qū)域Ⅲ中，區(qū)域Ⅰ的液態(tài)金屬直接流入該區(qū)域。其中一部分液態(tài)金屬由區(qū)域Ⅰ進(jìn)入?yún)^(qū)域Ⅲ時(shí)，被匙孔壁面下方表面張力吸引，沿V1方向流動(dòng)；另一部分液態(tài)金屬沿固液界面流動(dòng)較短距離后被匙孔壁面上方表面張力吸引，形成V2方向的流動(dòng)；剩余大部分液態(tài)金屬沿固液界面向熔池后方流動(dòng)，與熔池上方向下流動(dòng)的金屬相遇，產(chǎn)生S方向的流動(dòng)。匙孔壁面從下到上溫度呈遞減趨勢、表面張力呈遞增趨勢，導(dǎo)致V1方向上的液態(tài)金屬流速較V2大。這兩個(gè)方向上的液態(tài)金屬撞擊匙孔壁面而迅速改變流動(dòng)方向并相遇。由于液態(tài)金屬流動(dòng)方向的連續(xù)改變，并且V1方向上流速快，所以產(chǎn)生了順時(shí)針的渦流。該區(qū)域的渦流形成過程如圖6(c)所示。

圖6 熔池不同區(qū)域渦流形成過程 (a)Zone Ⅰ；(b)Zone Ⅱ；(c)Zone ⅢFig.6 Process of vortex formation in different region (a)Zone Ⅰ；(b)Zone Ⅱ；(c)Zone Ⅲ

2.3 電子束與匙孔壁面耦合行為

從圖4中可以看出，電子束焊接過程中，匙孔的形狀并不是固定不變的。區(qū)域Ⅲ中渦流的移動(dòng)使熔池變形，進(jìn)而影響匙孔的變化，提取圖4中熔池形狀和電子束在匙孔壁面上的作用位置，如圖7所示，圖中灰色部分代表熔池，綠色部分代表電子束與匙孔的作用位置。不難看出，電子束主要作用在匙孔壁面凸起的位置上。對(duì)于匙孔前壁面：匙孔壁面上凸起的部位接收電子束的能量，溫度的迅速升高使之產(chǎn)生蒸氣反沖壓力，由凸起變?yōu)榘枷?，如此反?fù)，匙孔前壁面的液態(tài)金屬快速波動(dòng)。對(duì)于匙孔后壁面：熔池的翻滾不斷地阻礙電子束與匙孔后壁面下側(cè)接觸，使其溫度低于金屬沸點(diǎn)而不能穩(wěn)定存在；在匙孔后方液態(tài)金屬表面張力和凹液面附加壓力的作用下，匙孔底部被液態(tài)金屬補(bǔ)充。由于焊接熱源以恒定的速度向前移動(dòng)，電子束不停地挖掘新的匙孔前壁面，而匙孔后壁面也不斷地被液態(tài)金屬補(bǔ)充，這就造成了匙孔一直是傾斜狀態(tài)，即匙孔底部具有一定的滯后性。

圖7 電子束與匙孔壁面的相互作用 (a)t=1.505s；(b)t=1.524s；(c)t=1.528s；(d)t=1.530sFig.7 Interaction between electron beam and keyhole (a)t=1.505s；(b)t=1.524s；(c)t=1.528s；(d)t=1.530s

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

由于在實(shí)際焊接時(shí)熔池的流動(dòng)行為很難檢測，因此，人們多采用數(shù)值模擬的手段來研究焊接熔池的動(dòng)態(tài)過程，并用對(duì)比焊縫橫截面形貌的方法來進(jìn)行驗(yàn)證數(shù)值模擬的有效性[20]。本研究模擬的焊縫熔化線形狀與實(shí)際焊縫熔化線形狀對(duì)比如圖8所示。從圖8可以看出，模擬的焊縫熔化線形狀與尺寸和實(shí)際焊縫熔化線接近，因此可以驗(yàn)證本研究的模擬結(jié)果是合理的。

圖8 焊縫橫截面試驗(yàn)與計(jì)算結(jié)果對(duì)比Fig.8 Comparison of simulated weld shape with actual weld joint morphology

3 結(jié)論

(1)綜合考慮了電子束與匙孔壁面的耦合作用，建立了電子束焊接三維瞬態(tài)模型，對(duì)10mm厚2219鋁合金的電子束焊接熔池周期性波動(dòng)過程進(jìn)行了數(shù)值模擬。

(2)焊接進(jìn)入準(zhǔn)穩(wěn)態(tài)后，熔池尺寸并不是恒定不變的，而是在小范圍內(nèi)震蕩，焊接熔池呈周期性波動(dòng)。

(3)根據(jù)液態(tài)金屬流動(dòng)規(guī)律可以將熔池分為3個(gè)區(qū)域，其中在區(qū)域Ⅱ和區(qū)域Ⅲ會(huì)產(chǎn)生渦流。區(qū)域Ⅰ中的液態(tài)金屬起到維持熔池體積穩(wěn)定的作用，區(qū)域Ⅱ中的逆時(shí)針渦流起到維持熔池表面積的作用，區(qū)域Ⅲ中的順時(shí)針渦流促進(jìn)熔池坍塌，使匙孔變得不穩(wěn)定。

(4)匙孔壁面上的電子束并不是均勻分布的，兩者的相互作用使匙孔底部具有滯后性。