張成國,武小娟,孟凡玲

(沈陽理工大學(xué) 材料科學(xué)與工程學(xué)院,沈陽 110159)

1060鋁攪拌摩擦焊穩(wěn)態(tài)流場的數(shù)值模擬

張成國,武小娟,孟凡玲

(沈陽理工大學(xué) 材料科學(xué)與工程學(xué)院,沈陽 110159)

通過計算流體力學(xué)軟件Fluent構(gòu)建了1060鋁攪拌摩擦焊穩(wěn)態(tài)流場有限元模型,對塑型金屬材料的三維流場進(jìn)行了研究。數(shù)值模擬結(jié)果表明:在工件厚度方向上,塑形金屬的流動受軸肩影響明顯,速度隨到軸肩距離的增加而逐漸減小;縱截面上塑形金屬流動方向多處發(fā)生改變;橫截面塑性金屬的流動受速度合成不同的影響,并不對稱。前進(jìn)邊速度和變形范圍比返回邊的大。

攪拌摩擦焊;流場;數(shù)值模擬

攪拌摩擦焊接是由英國焊接研究所針對常規(guī)熔焊焊接性差的輕質(zhì)有色金屬開發(fā)的固相連接技術(shù)[1],隨著該技術(shù)的不斷發(fā)展,在工業(yè)各領(lǐng)域得到了廣泛應(yīng)用[2]。

攪拌摩擦焊過程中,塑形金屬的流動對于焊縫成形和焊接接頭質(zhì)量有重要影響[3],一直頗受國內(nèi)外學(xué)者的注意和研究。對塑型金屬流場的模擬,是一種高效率低成本的方法,能補足試驗方法難以準(zhǔn)確得到攪拌頭附近塑型金屬的流動情況和規(guī)律的短板。Colligan K等人[4]對異種金屬的攪拌摩擦焊研究表明,和同種材料焊接時比,材料都出現(xiàn)了強(qiáng)烈的機(jī)械混合和塑性變形,不對稱性特征明顯。Colegrove P A等人[5]使用Fluent模擬了攪拌摩擦焊塑性金屬材料的流動情況,比較了不同攪拌頭的流場流動情況。本文通過有限元軟件Fluent對攪拌摩擦焊穩(wěn)定階段的攪拌頭周圍流場情況進(jìn)行了模擬,并通過不同的截面展示和分析塑形金屬流動的情況。對優(yōu)化攪拌頭的參數(shù)、探尋攪拌摩擦焊焊接機(jī)理和確定缺陷位置具有一定的意義。

1 攪拌摩擦焊的三維有限元模型的建立

攪拌摩擦焊是個相當(dāng)復(fù)雜的過程,通常分解為攪拌頭旋轉(zhuǎn)擠進(jìn)工件階段、穩(wěn)定焊接階段和攪拌頭運動停止并退出板材的三個過程。為模擬順利進(jìn)行和縮減計算時間,本模型僅對焊接穩(wěn)定階段進(jìn)行模擬和分析,此階段工件與攪拌頭的相對運動保持不變,焊接橫截面的性能和形狀類似。將塑化金屬在穩(wěn)定焊接階段的流動看作為一種穩(wěn)態(tài)的流動。

1.1 模型的假設(shè)條件

攪拌摩擦焊過程中焊縫塑性金屬物理性能和流動性能與液體相近,具有流體運動的特性[6],為實現(xiàn)模擬的順利進(jìn)行和縮短模擬運算時間,提出以下假設(shè)條件:(1)攪拌頭與工件上表面垂直,忽略傾角對焊接過程的影響;(2)與工件直接接觸的攪拌頭表面符合固體無滑移前提,也就是攪拌頭周圍塑性材料隨著攪拌頭旋轉(zhuǎn)時,其遷移速度與攪拌頭旋轉(zhuǎn)線速度相等;(3)材料密度恒定,流體不可壓縮,將塑性金屬在焊接過程中的移動設(shè)定為穩(wěn)定狀態(tài)下流動。

1.2 網(wǎng)格的劃分和邊界條件

用三維設(shè)計軟件SolidWorks構(gòu)建工件和攪拌頭的幾何模型,通過布爾減操作獲得流場模型。工件簡化成80mm×50mm×7mm同塊長方體,攪拌頭由軸肩和圓柱攪拌針構(gòu)成。攪拌針長度為6.3mm,直徑為8mm,軸肩直徑為20mm。攪拌頭位于工件中心,以降低邊界條件的設(shè)定與實際情況的差異帶來的偏差[7]。坐標(biāo)原點位于攪拌針端面圓心位置。三維坐標(biāo)中,縱截面上焊接前進(jìn)方向沿著X坐標(biāo)軸負(fù)方向,橫截面上焊縫前進(jìn)邊與Y坐標(biāo)軸負(fù)方向重合,攪拌針端面為Z=0mm截面。為避免構(gòu)建模型的復(fù)雜化和不運用動態(tài)網(wǎng)格,模型中軸肩簡化為一個平面,攪拌針表面無螺紋[8]。數(shù)據(jù)導(dǎo)入應(yīng)用前處理軟件GAMBIT,流體網(wǎng)格劃分如圖1所示。

圖1 網(wǎng)格的劃分

在實際攪拌摩擦焊焊接的過程中,工件被牢牢固定,攪拌頭邊旋轉(zhuǎn)邊前進(jìn),而模擬過程中依照常規(guī)的方式設(shè)置為攪拌頭旋轉(zhuǎn),工件做平移運動。圖2所示的為流場的邊界條件。模型中共設(shè)有三種邊界面,流體從速度進(jìn)口面流入計算區(qū)域,經(jīng)過旋轉(zhuǎn)的攪拌頭面,經(jīng)壓力出口面流出,塑性金屬在此密封的通道中流動[9]。圖中左側(cè)面邊界條件類型設(shè)置為速度進(jìn)口面,方向和焊接速度相反,大小相等。圖中右側(cè)面邊界條件類型設(shè)置為壓力出口面,此為相對壓力,大小為0,因為Fluent中壓力是相對于操作壓力而言。攪拌針壁面設(shè)置為動墻壁面,默認(rèn)無滑移,大小和方向同攪拌頭旋轉(zhuǎn)速度相同,為700rpm。軸肩壁面設(shè)置相同。左右和上下墻壁面設(shè)置為靜止壁面。

圖2 邊界條件

1.3 材料參數(shù)

本文研究內(nèi)容是攪拌摩擦焊接穩(wěn)定階段流場的模擬,密度設(shè)定恒定為2700kg/m3,粘度是重要參數(shù),測量難度較高,引用文獻(xiàn)[10]數(shù)值0.4kPa·s,其他參數(shù)[11]見表1。

表1 鋁1060不同溫度下的材料參數(shù)

2 模擬結(jié)果與分析

圖3顯示的是整體流場圖,能得出焊縫區(qū)域流場受到了攪拌頭旋轉(zhuǎn)帶動作用,攪拌頭附近塑性金屬隨著攪拌頭旋轉(zhuǎn),但攪拌頭影響范圍有限。同時各處流動速度是不相同的:水平方向上,在軸肩影響范圍內(nèi),材料的流動速度隨旋轉(zhuǎn)距離的增加而增大。軸肩輪廓外緣達(dá)到了最大值。

在軸肩影響范圍之外,材料的流動速度急劇減小。厚度方向上,能夠看出離工件軸肩越遠(yuǎn)的流場流速越小,離軸肩距離近的位置,流場范圍大,流速高。隨著離軸肩距離的增加,流場范圍逐漸減小,流速逐漸降低。

圖3 整體流場圖

模擬軟件的優(yōu)勢是可以在任意截面處觀察塑性金屬流動情況,下面選取幾處典型截面給予更加清晰的顯示。厚度方向上截取位于Z=6 mm、Z=4mm和Z=0mm水平截面。橫向和縱向的截面分別位于X=0mm和Y=0mm處。流場模擬結(jié)果如圖4所示。

圖4 焊縫不同厚度處的水平截面流場圖

取位于焊縫上表面較近處水平截面,觀察流場情況,如圖4a所示。此位置距離軸肩0.3mm,受軸肩影響較大。能夠看出塑性金屬流動基本呈現(xiàn)和軸肩相同形狀大小的形貌,但區(qū)域比軸肩略大些。最大流動速度出現(xiàn)在軸肩影響范圍邊緣位置。前進(jìn)邊和返回邊兩側(cè)的流線密度不相同。塑性金屬流動在前方焊縫中心位置出現(xiàn)分流,軸肩帶動部分材料從前進(jìn)側(cè)順著攪拌頭外緣到達(dá)軸肩后方區(qū)域。另外部分塑性金屬繼續(xù)隨攪拌針旋轉(zhuǎn),并繞過焊縫中心線,繼續(xù)遷移一定距離落在返回側(cè)。

取焊縫厚度中間處水平截面,觀察流場情況,如圖4b所示。此位置距軸肩2.3mm,受軸肩影響已經(jīng)較弱。相比于圖4a塑性金屬流動劇烈區(qū)域明顯減小,僅略大于攪拌針?biāo)浇孛?。能夠得出其主要受攪拌針作用。流體最高速度值也有所減小,能夠得出軸肩對流體流動的帶動效果比攪拌針要好。兩者相同的是前進(jìn)邊和返回邊兩側(cè)的流線密度不相同。流動速度上,前進(jìn)邊大于返回邊。這是由于返回邊流動速度方向和攪拌針旋轉(zhuǎn)方向重合,兩種作用互相疊加作用的結(jié)果[12]。

取攪拌針端面處水平截面,觀察流場情況,如圖4c所示?？梢钥闯龃私孛嫠苄谓饘俚牧鲃虞S肩已經(jīng)沒有影響,僅受攪拌針端面影響,但作用也明顯減弱。塑形金屬流動區(qū)域小,流動速度低,材料雖有遷移,但范圍較小,流動不充分,這也是容易出現(xiàn)焊接缺陷的位置。

圖5為焊縫縱向截面流場情況圖,能夠發(fā)現(xiàn),最大切向速度出現(xiàn)在軸肩外緣位置。塑形金屬向前流動中遇攪拌頭的阻擋和旋轉(zhuǎn)帶動,軸肩下方附近塑形金屬呈漩渦狀,中間大部分材料繞著攪拌針旋轉(zhuǎn),還有小部分材料從攪拌針端面下方流過到達(dá)攪拌針后方。攪拌針后方附近的塑性金屬出現(xiàn)分流,貼緊攪拌針的跟隨攪拌針旋轉(zhuǎn),靠近攪拌針的向上遷移趨勢大,并且越接近軸肩位置,此部分材料范圍越大。離攪拌針距離略大的,向后遷移趨勢大。

圖5 焊縫縱截面流場情況圖

圖6為焊縫橫截面流場情況圖,由圖6可看出,在軸肩帶動影響下,攪拌頭附近塑形金屬流動方向明顯發(fā)生改變。部分塑性金屬在軸肩下方流動遷移形成漩渦,部分塑形金屬隨著軸肩旋轉(zhuǎn)遷移,流動到后方填補移動造成的瞬時空腔。圖中左側(cè)是流場中的前進(jìn)邊,塑性金屬的流動速度方向與軸肩的旋轉(zhuǎn)速度方向一致。右側(cè)是流場中的返回邊,塑性金屬的流動速度方向與軸肩的旋轉(zhuǎn)速度方向正好相反。這兩種速度效果相互疊加或抵消,造就了前進(jìn)邊塑性金屬的最高速度和變形區(qū)域比返回邊的大。

圖6 焊縫橫截面流場情況圖

3 結(jié)論

(1)使用三維造型軟件Solidworks構(gòu)建了流場的三維幾何模型,使用前處理軟件Gambit通過非均勻四面體網(wǎng)格方法劃分了網(wǎng)格,使用計算模擬軟件Fluent模擬出了攪拌摩擦焊的流場情況。運用整體和不同截面對流場情況進(jìn)行詳細(xì)說明。

(2)在1060攪拌摩擦焊的流場中,在工件水平方向上,塑形金屬的流動速度隨到攪拌頭旋轉(zhuǎn)中心的距離的增加而減小;在工件厚度的方向上,塑形金屬的流動受軸肩影響明顯,速度隨到軸肩距離的增加而逐漸減小。

(3)縱截面上塑形金屬流動方向多處發(fā)生改變,有分流現(xiàn)象;橫截面塑性金屬的流動受速度合成不同的影響,并不對稱。前進(jìn)邊速度和變形范圍比返回邊的大。

[1] 王國慶,趙衍華.鋁合金的攪拌摩擦焊接[M].中國宇航出版社,2010:1-3.

[2] 關(guān)橋.攪拌摩擦焊—未來的連接技術(shù)[C]//北京：攪拌摩擦焊國際會議，2012:10-11.

[3] 王善林,柯黎明,邢麗.攪拌頭形狀對焊縫塑化金屬流動行為的影響[J].南昌航空大學(xué)學(xué)報:自然科學(xué)版,2005,19(1):62-66.

[4] Colligan K.Material flow behavior during friction stir welding of aluminum[J].Welding Journal,1999,172(2):220.

[5] Colegrove P A,Shercliff H R.3-Dimensional CFD modelling of flow round a threaded friction stir welding tool profile[J].Journal of Materials Processing Technology,2005,169(2):320-327.

[6] 欒國紅.攪拌摩擦焊流變特性研究[J].航空制造技術(shù),2003(11):22-25.

[7] 張文斌.2A12鋁合金攪拌摩擦焊溫度場及塑性材料流動的數(shù)值分析[D].濟(jì)南：山東大學(xué),2010:32-33.

[8] 王鑫.7A52鋁合金攪拌摩擦焊實驗研究及流動仿真[D].北京：清華大學(xué),2009:48-49.

[9] 姬書得,劉偉,張利國，等.基于湍流計算模型的TC4鈦合金攪拌摩擦焊過程的流場分析[J].熱加工工藝,2012,41(13):172-174.

[10] 李寶華.攪拌頭形狀與工藝參數(shù)對塑性金屬遷移行為的影響[D].南昌：南昌航空大學(xué),2008:73-74.

[11] 王艷輝.1060鋁攪拌摩擦焊溫度場數(shù)值模擬[D].鎮(zhèn)江：江蘇科技大學(xué),2012:39-40.

[12] 金延野.攪拌摩擦焊接過程塑性流動特征研究[D].沈陽：沈陽航空航天大學(xué),2013:33-34.

(責(zé)任編輯：王子君)

NumericalSimulationof1060AluminumStirringFrictionWeldingSteadyFlowField

ZHANG Chengguo,WU Xiaojuan,MENG Fanling

(Shenyang Ligong University,Shenyang 110159,China)

By computational fluid dynamics software Fluent,1060 aluminum FSW steady flow finite element model is established,the numerical simulation of three-dimensional flow field of metal plastic material is studied.The results show that,in the direction of the thickness of the workpiece,the flow of the plastic is influenced by the shoulder,and the velocity decreases with the increase of the distance from the shoulder;the flow direction of the plastic is changed in many directions in the longitudinal section;the flow of the plastic is affected by velocity synthesis and not symmetrical in the cross section.The A_side velocity and deformation range are larger than the R_ side.

FSW;flow field;numerical simulation

TB391.9

A

2016-09-23

張成國(1991—)，男，碩士研究生；通訊作者：武小娟(1973—)，女，副教授，博士，研究方向：材料成型連接技術(shù)。

1003-1251(2017)05-0049-04