基于CFD數(shù)學(xué)模型的電磁攪拌技術(shù)對鋁合金激光焊填絲元素分的影響及計算機分析

利 莉，胡治宇

（景德鎮(zhèn)學(xué)院，江西景德鎮(zhèn)333000）

基于CFD數(shù)學(xué)模型的電磁攪拌技術(shù)對鋁合金激光焊填絲元素分的影響及計算機分析

利 莉，胡治宇

（景德鎮(zhèn)學(xué)院，江西景德鎮(zhèn)333000）

鋁合金激光焊接中，采用外加電磁場，操縱感應(yīng)體積力引起液相金屬流動，從而模擬含硅填料的元素分布。為了對熔池內(nèi)部液體動力學(xué)過程有更深刻的認(rèn)識，通過一個CFD模型借助計算機模擬熔體流動。本研究的主要內(nèi)容是，不同頻率外加磁場作用下，將填絲元素分布模擬結(jié)果與相同參數(shù)下的實驗結(jié)果進行對比。結(jié)果表明，兩種情況下，填料合金元素分布呈空間周期性。根據(jù)CFD模型可知，周期性感應(yīng)電磁體積力引起了熔融流體的調(diào)制，從而導(dǎo)致填料元素分布的變化。

鋁合金；激光焊接；電磁場；CFD模型

0 前言

磁場常見于各種應(yīng)用，尤其在金屬工業(yè)中應(yīng)用廣泛。使用磁場影響激光焊接過程一直都是焊接領(lǐng)域的研究重點之一，眾多研究結(jié)果表明，通過施加一個穩(wěn)定的磁場，可以顯著影響焊縫橫截面的形狀。含Si填充焊絲鋁合金激光焊接中，低頻交變磁場對熔珠稀釋有一定作用，整個焊縫均勻分布的Si含量大于2%時，可以顯著降低鋁合金（如6XXX-合金）熱裂紋的敏感性[1]。另有結(jié)果表明，磁場作用下進行激光焊接時主要是流體動力學(xué)過程[2]。對磁場參數(shù)與流體流動間的相互作用認(rèn)識越深刻，越能增加焊縫填充材料元素均勻分布的能力。

為了解決理論問題，建立了一個CFD模型，可以通過納維斯-托克斯方程處理由激光輻射感應(yīng)的流體流動現(xiàn)象和液態(tài)熔池與磁場間的相互作用這兩方面問題[3]。以往的模擬工作主要集中在兩個方向：一是對焊接現(xiàn)象進行詳細(xì)的模擬（如激光束吸收過程或污染的影響）；二是在可接受范圍內(nèi)簡化了所需的計算時間。

本研究的目的是為了更深刻地了解在含有合金元素的填充焊絲鋁合金全熔透激光焊接中，磁力攪拌對焊絲材料濃度分布的影響。研究重點是磁場和激光誘導(dǎo)液態(tài)熔池之間的復(fù)雜、瞬態(tài)的流體動力作用理論，其中熔池會導(dǎo)致填料凝固后特定的濃度分布。本研究建立了一個有關(guān)焊接過程的CFD模型。在適當(dāng)?shù)挠嬎銜r間范圍內(nèi)，該模型可用于模擬計算0.3 s的焊接順序。該模型可用于研究磁場不同頻率和磁通密度對元素分布的影響。為了進行比較，進行了全熔透激光Al99.5鋁焊接實驗，采用硅濃度高達(dá)18%的超共晶填絲。

1 實驗設(shè)備和實驗過程

焊接實驗采用輸出功率為5.5 kW的TruDisk 8002盤形激光器，光學(xué)系統(tǒng)的焦距為250 mm，準(zhǔn)直長度250 mm；纖維直徑200 μm，標(biāo)稱焦斑值為200μm。全熔透堆焊選擇厚3 mm的Al99.5薄板，焦點在工件上方2 mm處。填料為AlSi18焊絲（平均Si含量為18%），直徑1.2 mm。焊接速度v0=8 m/min，送絲速度vd=6 m/min。實驗為全熔焊，焊頭提供了不同磁通量密度B0和不同頻率f的磁場，如圖1a所示。磁場直流磁通密度分布與激光束同軸，隨著距焊頭軸向距離的增加，磁通密度幅值稍有減?。ㄒ妶D1b）。坐標(biāo)系的原點位于板材和激光束的交點。如圖1c所示，磁場脈沖t形狀為三角形。焊接后，對工件的縱截面取樣，采用EXD方法定量分析Si元素的分布。

圖1 焊接試驗

2 數(shù)值模型的建立

建立的帶有填充焊絲的全熔透堆焊CFD模型，可以用來計算在與激光束同軸的交變磁場（見圖1b）影響下的熔體流動。為了計算帶有填充焊絲下的熔池和焊縫形狀，采用如圖2所示的幾何模型。為了形成3 mm厚的鋁板，將模型設(shè)計成為一個3 mm厚的四邊形。在四邊形上連接一個直徑為0.6 mm的圓柱作為填充焊絲。兩個結(jié)構(gòu)通過旋轉(zhuǎn)對稱幾何形狀，即小孔形狀進行連通。這個簡化的小孔模型，通過在小孔表面附加蒸發(fā)溫度，可以模擬激光束的吸收。

圖2 填充焊絲工藝和近似小孔形狀的幾何模型

小孔周圍流體、熔池內(nèi)的馬蘭哥尼對流和浮力隨相變的影響（如潛熱或凝固和熔化潛熱的驅(qū)動力等）是通過增加動量和能量方程中的力與能量源來實現(xiàn)的。

根據(jù)式（1），通過增加動量方程中的體積力來研究浮力的影響

式中 cp為恒壓下的比熱容；T為控制體的溫度；T0為室溫；ΔH為熔化潛熱。β為液相比例，見式（3）

式中 TS和TL分別為固相和液相溫度。溫度范圍為34 K時，控制體的液相比例從0擴大到1。為了模擬動態(tài)凝固過程，通過增加動量方程中如式（4）所示的體積力來降低焊接速度為v0時，液相轉(zhuǎn)變?yōu)楣滔鄷r的流動。

式中 ▽T為局部溫度梯度；γT為表面張力溫度系數(shù)。這導(dǎo)致了液體沿著溫度梯度加速。負(fù)溫度系數(shù)情況下（此模型使用的鋁合金），液體向溫度降低的方向加速。磁場和速度場為的導(dǎo)電流體間的相互作用，引起了流動熔體和基體材料間的電流密度。電流密度和其他內(nèi)部或外部電場的表達(dá)式為

式中 t和t0分別為時間和時間偏量。電流和外部磁場（洛倫茲）力間的相互作用也可以添加到動量方程，如式（8）所示

用于這些計算流體特性的材料是純鋁（液體和固體），如表1所示?；w材料和焊絲的假設(shè)也是為了簡化計算。由焊絲中元素Si的分布可作為一個無質(zhì)量標(biāo)量φ，通過式（9）方程計算了其時間和空間分布

表1 模擬中的材料性能

計算中使用的動力學(xué)和熱邊界如表2所示。通過有限體積法求解控制方程，采用商業(yè)cfx5.1高分辨率CFD求解器平流方案，以及一階向后歐拉方程計算瞬時流量。因此幾何模型被離散成近百萬個四邊形控制體。無磁場的穩(wěn)態(tài)求解作為瞬態(tài)計算的初始條件。假定時間步長為1e-03 s。從準(zhǔn)靜態(tài)解開始，模擬焊接過程的總持續(xù)時間為0.3 s。

3 結(jié)果

3.1 頻率的影響

在焊接速度8 m/min，送絲速度vd=6 m/min的條件下進行計算。為了研究磁場頻率對元素分布的影響，選擇參考磁通量B0=160 mT。計算選擇的頻率為10 Hz、15 Hz和20 Hz。0.3 s焊接后幾何模型縱向?qū)ΨQ面濃度分布的計算機模擬結(jié)果如圖3所示。圖3中黑色線條表示熔池周圍基體材料和填絲的熔融等溫線。熔池內(nèi)浮力、小孔周圍的流體、馬蘭哥尼效應(yīng)和磁感應(yīng)體積力都導(dǎo)致熔池內(nèi)的強對流，對焊絲材料進行稀釋。在熔池表面附近，尤其是上表面，熔池被拉長。

表2 動態(tài)和熱邊界條件

圖3 熔池縱截面無質(zhì)量標(biāo)量的分布和熔融等溫線

3.2 磁通密度的影響

第二組計算用來研究參考磁通量B0對元素分布的影響。頻率固定10Hz，參考磁通量分別為160mT、260mT和360 mT。幾何模型縱向?qū)ΨQ平面計算濃度分布如圖5所示。

圖4 幾何模型上表面濃度分布的計算機模擬結(jié)果

圖5 幾何模型縱向?qū)ΨQ平面計算濃度分布

濃度分布的周期形狀隨磁通量的增加而改變。值得注意的是，磁通密度增加，焊縫上半部分濃度分布更加均勻。

圖6為頻率固定為10 Hz時，隨著磁通密度增加上表層濃度分布。由圖6可以看出，當(dāng)磁通密度從160 mT增加到260 mT時，對濃度分布影響不大，但可以發(fā)現(xiàn)當(dāng)磁通密度為360 mT時，高濃度分布圖略有縮小。

圖6 頻率固定為10 Hz時，磁通密度增加上表層濃度分布

3.3 實驗結(jié)果

為了與計算結(jié)果進行比較，對焊縫縱截面硅元素分布進行了EDX掃描分析，結(jié)果如圖7所示。白色區(qū)域為硅元素高濃度分布區(qū)，黑色區(qū)域主要為基體材料。

圖7 不同磁場下EDX測試的縱截面硅元素分布

4 討論

研究表明，洛倫茲力對熔池內(nèi)速度場有顯著影響，如果流體流動與磁場有一定夾角，洛倫茲力主要作用于原始熔體流動[5]。熔池內(nèi)的流動方向是從小孔前面流向熔池后部，幾乎垂直于磁場。

在這種假設(shè)下，計算的濃度分布，尤其是周期性形狀的一個合理解釋是，熔池內(nèi)速度場的周期控制引起熔體流動產(chǎn)生偏移。從式（6）和式（8）可以理解洛倫茲力主要作用于水平方向，引起流體阻力。熔體流動在水平方向的偏移可以減小磁場與熔體間的夾角，從而減小流動阻力。調(diào)制的熔體流動主要決定了元素的分布，而本研究提出的模擬忽略了擴散等影響。

由于體積力與局部速度成比例，因此高速度能引起高體積力。從計算結(jié)果可知，由于馬蘭哥尼加速度的最高速度在上下表面。此外，還應(yīng)考慮磁場在徑向和軸向的不均勻性。最高磁通密度位于熔池上表面附近。因此，考慮到馬蘭哥尼對流強大的加速度，此區(qū)域會產(chǎn)生最高的感應(yīng)力。這個假設(shè)也與其他結(jié)果相符。

考慮的第二個效應(yīng)是后熔池邊界溶體流動的偏移也會引起熔體從熔池表面向內(nèi)部流動。這種流動模式也是通過感應(yīng)力調(diào)制，也能影響濃度分布。

對頻率影響的計算結(jié)果顯示，濃度分布呈周期性調(diào)制（尤其是頻率為10 Hz和20 Hz）。由磁場和熔體流動偏轉(zhuǎn)共同引起的脈沖力也可以很容易理解。有個重要問題需要注意，調(diào)制不總是嚴(yán)格的周期性，可以看到15 Hz沒有明顯的調(diào)制。原因可能是溶體流動的瞬態(tài)失穩(wěn)，并且沒有任何外加磁場。導(dǎo)致頻率為15 Hz時濃度分布結(jié)果的原因是，流體不穩(wěn)定性與磁場作用下流體調(diào)制共同作用的一個瞬態(tài)疊加，造成了干擾，阻礙了周期性。

實驗結(jié)果證實了濃度周期性分布，只是周期頻率與計算結(jié)果不同。原因是模型的簡化，尤其是小孔不穩(wěn)定性的影響。磁通密度對元素分布的均勻性有顯著影響。從計算結(jié)果可以看出，隨著磁通密度的增加，高濃度和低濃度模式越來越模糊，這表明濃度的均勻性增加。

5 結(jié)論

（1）磁場頻率和磁通密度對填料合金元素濃度分布有顯著影響。

（2）沒有特定調(diào)制時，磁場頻率主要決定局部濃度分布，導(dǎo)致周期性濃度分布。

（3）磁通密度決定了感應(yīng)力的絕對值，對濃度分布的均勻性有顯著影響，隨著磁通密度的增加，合金元素濃度均勻性增強。

[1]周世杰，李云濤，宋繼順，等.激光焊接熱源和焊速對溫度場影響的數(shù)值模擬[J].電焊機，2014，44（10）：85-89.

[2]李志剛，賈慧芳，李剛，等.水下焊接穩(wěn)弧裝置的磁場及其雙流道流場[J].電焊機，2014，44（3）：55-59.

[3]劉海青，劉秀波，孟祥軍，等.金屬基體激光熔覆陶瓷基復(fù)合涂層的裂紋成因及控制方法[J].材料導(dǎo)報，2013（11）：60-63.

[4]Tang Z，Gatzen M.Influence on the dilution by laser welding of aluminium with magnetic stirring.In：Laser Assisted Net ShapeEngineering(LANE2010)Physics Procedia 5/2，eds：M.Schmidt，F(xiàn).Vollertsen，M.Geiger.Meisenbach Bamberg，2010：125-137.

[5]陳淑仙，李明偉.CZ法砷化鎵單晶生長中熔體流動狀態(tài)轉(zhuǎn)換[J].無機材料學(xué)報，2007（1）：15-20.

Effect and computer analysis of electromagnetic stirring technique on the filler wire element distribution in laser welding of aluminium based on the CFD mathematical model

LI Li，HU Zhiyu
（Jingdezhen University，Jingdezhen 333000，China）

In order to simulate the element distribution resulting from the flow of the liquid metal due to induced volume forces during laser welding of aluminium with silicon containing filler wire，additional external electromagnetic fields are used.A computational fluid dynamics(CFD)model is developed to numerically analyse the flow of the melt for a more profound understanding of fluid dynamics process inside the pool.The mainly purpose of this work is the simulating element distribution result compared with the experimental result of filler wire material under a coaxial magnetic field with different frequencies.The element distribution of the filler material has a spatial periodicity in the two cases.It is shown that the periodic induced electromagnetic volume force results in a modulation of the melt flow and then the change of the distribution of the filler material according to the CFD model.

aluminum alloy；laser welding；electromagnetic field；CFD mathematical model

TG249.2

：A

1001-2303（2015）09-0118-05

10.7512/j.issn.1001-2303.2015.09.26

2015-03-31

利 莉（1983—），女，江西景德鎮(zhèn)人，講師，碩士，主要從事計算機應(yīng)用與計算機網(wǎng)絡(luò)的研究工作。