萬震宇， 周 霞， 張 昭,2

(1． 大連理工大學(xué) 運(yùn)載工程與力學(xué)學(xué)部 工程力學(xué)系工業(yè)裝備結(jié)構(gòu)分析國家重點實驗室，遼寧 大連116024;2. 大連理工大學(xué) 國際計算力學(xué)研究中心,遼寧 大連116024)

?

鋁合金攪拌摩擦焊后人工時效力學(xué)性能數(shù)值模擬

萬震宇1， 周 霞1， 張 昭1,2

(1． 大連理工大學(xué) 運(yùn)載工程與力學(xué)學(xué)部 工程力學(xué)系工業(yè)裝備結(jié)構(gòu)分析國家重點實驗室，遼寧 大連116024;2. 大連理工大學(xué) 國際計算力學(xué)研究中心,遼寧 大連116024)

采用KWN模型構(gòu)建攪拌摩擦焊接中Al-Mg-Si系鋁合金沉淀相演化計算模型，通過將屈服強(qiáng)度分為晶粒大小、固溶相和析出相三部分貢獻(xiàn)，可以計算平板攪拌摩擦焊后的屈服強(qiáng)度和硬度。進(jìn)一步研究不同焊后人工時效條件下，焊接平板力學(xué)性能變化的機(jī)理。結(jié)果表明：更長的焊后保溫時間有利于攪拌區(qū)力學(xué)性能的回復(fù)；較高的保溫溫度有利于攪拌區(qū)力學(xué)性能的快速回復(fù)，但是當(dāng)溫度高于200 ℃時，長時間保溫會使母材軟化，不利于力學(xué)性能回復(fù)；通過焊后人工時效不能明顯改善熱影響區(qū)的力學(xué)性能。

攪拌摩擦焊接；析出相；焊后熱處理；人工時效

作為一種固態(tài)連接技術(shù)，相比較于傳統(tǒng)熔焊中接頭狀態(tài)的改變而導(dǎo)致的扭曲變形以及較大的殘余應(yīng)力，攪拌摩擦焊接具有接頭力學(xué)性能好、污染低等優(yōu)點[1-2]，迅速在航空航天、船舶、高速列車、汽車等領(lǐng)域取得工業(yè)應(yīng)用[3-4]。在攪拌摩擦焊接中，焊接接頭被劃分為幾個不同的焊接區(qū)域：焊核區(qū)、熱影響區(qū)、熱機(jī)影響區(qū)和母材區(qū)[5-6]，不同焊接區(qū)域具有不同的力學(xué)性能[7-8]，主要是由于焊接過程導(dǎo)致不同區(qū)域具有不同的微觀組織結(jié)構(gòu)。建立微觀組織結(jié)構(gòu)數(shù)學(xué)模型對于了解攪拌摩擦焊接焊縫質(zhì)量控制機(jī)理具有重要意義。

結(jié)合實驗測試，F(xiàn)ratini等[7]提出了一種預(yù)測6082鋁合金攪拌摩擦焊接晶粒尺寸的經(jīng)驗公式，將晶粒尺寸與應(yīng)變、應(yīng)變率和溫度歷史相對應(yīng)；然而，該模型中相關(guān)系數(shù)與所采用的應(yīng)變度量方式和具體數(shù)學(xué)模型有較大關(guān)聯(lián)，對于不同的數(shù)學(xué)模型以及不同的應(yīng)變度量方式，需要重新測定或者通過嘗試確定相關(guān)參數(shù)。對于攪拌摩擦焊接微觀組織結(jié)構(gòu)的數(shù)值模擬，更廣泛的方法是基于Zener-Hollomon參數(shù)的度量方法。Azizieh等[8]確定了純Mg攪拌摩擦焊接平板晶粒尺寸和Z參數(shù)之間的關(guān)系，用于幫助分析所獲得的實驗結(jié)果。Zhang等[9]計算了6061鋁合金攪拌摩擦焊接平板的Z參數(shù)變化，通過Z參數(shù)和晶粒尺寸的關(guān)系式計算了攪拌區(qū)晶粒尺寸，并研究了攪拌頭尺寸和轉(zhuǎn)速對Z參數(shù)和攪拌區(qū)晶粒尺寸的影響。Gerlich等[10]則利用Z參數(shù)反推計算了7075鋁合金攪拌摩擦點焊過程中的應(yīng)變率。Z參數(shù)的使用和度量對于快速計算攪拌區(qū)晶粒尺寸具有重要意義。然而，為了進(jìn)一步反映攪拌摩擦焊接過程中的晶粒變化過程，Zhang等[11]發(fā)展了一種基于Monte Carlo方法的攪拌摩擦焊接過程的晶粒生長模型，可以反映不同溫度歷史作用下攪拌摩擦焊接平板不同部位的晶粒生長行為?；谶@一Monte Carlo模型，Grujicic等[12]發(fā)現(xiàn)攪拌摩擦焊接不同焊接區(qū)域晶粒變化是晶粒生長和晶粒細(xì)化過程相互競爭的結(jié)果。

對于Al-Mg-Si系鋁合金的攪拌摩擦焊接，僅僅晶粒尺寸不足以全面反映微觀組織結(jié)構(gòu)的變化過程，也不能完整解釋力學(xué)性能變化的機(jī)理，以6005鋁合金為例，攪拌摩擦焊接平板的強(qiáng)度不僅僅取決于晶粒，也取決于第二相的析出和基體中溶質(zhì)固溶，其中，析出相對力學(xué)性能的影響更大[13]。在攪拌摩擦焊接中，以Mg2Si為典型代表的第二相顆粒會經(jīng)歷一個完整的轉(zhuǎn)變過程：即過飽和固溶體→原子團(tuán)聚→G.P.區(qū)→ β″(Mg5Si3) → β′(Mg5Si3) → β(Mg2Si)，最終生成的Mg2Si作為析出相彌散在基體上，以細(xì)小彌散狀分布的析出相可以很大程度上強(qiáng)化原有的鋁合金，同時，焊后熱處理工藝可以明顯改變攪拌摩擦焊接平板的力學(xué)性能[14-16]。有必要發(fā)展一種結(jié)合第二相粒子固溶、析出機(jī)理的攪拌摩擦焊接力學(xué)性能計算模型，可以加深對攪拌摩擦焊接機(jī)理的了解和優(yōu)化控制相關(guān)的工藝參數(shù)以提高焊接質(zhì)量，以此為基礎(chǔ)，進(jìn)一步討論焊后人工時效工藝的設(shè)計與優(yōu)化，對攪拌摩擦焊接的進(jìn)一步應(yīng)用具有重要意義。

1 數(shù)學(xué)模型

對攪拌摩擦焊接過程的數(shù)值模擬采用網(wǎng)格重剖分模型[17]，有限元網(wǎng)格如圖1所示，焊接速率為100 mm/min，攪拌頭轉(zhuǎn)速為1000 r/min，攪拌頭移動過程持續(xù)135 s。6005鋁合金熱導(dǎo)率為180.2 W/(m·K)，比熱容為884.7 J/(kg·K)，與攪拌頭之間熱傳導(dǎo)系數(shù)為20000 W/(mm2·K)，與空氣熱交換系數(shù)為200 W/(mm2·K)。材料力學(xué)性能是溫度的函數(shù)，如圖2所示。

圖1 攪拌摩擦焊接有限元模型Fig.1 Finite element model of FSW

圖2 不同溫度下的力學(xué)性能曲線Fig.2 Mechanical properties used in simulation under different conditions (a)400 ℃;(b)500 ℃

接觸模型采用剪切摩擦模型[5]，

f=m·k

(1)

式中：m為摩擦剪切因子，取為0.8；k為剪切強(qiáng)度。

溫度和力學(xué)響應(yīng)在同一增量步內(nèi)分別進(jìn)行差分求解，采用如下有限元格式，

(2)

Ku=P

(3)

式中：K為剛度陣；u為位移列陣；P為載荷列陣；CT為熱容矩陣；KT為熱傳導(dǎo)矩陣；PT為熱載荷列陣。

在攪拌摩擦焊接過程及焊后熱處理中，析出相的形核率表征為，

(4)

(5)

粒子長大速率為[18]：

(6)

式中：Cp為析出相溶質(zhì)濃度，Ci為析出相粒子與基體界面溶質(zhì)濃度，D為擴(kuò)散系數(shù)，r為析出相粒子半徑。

界面溶質(zhì)濃度Ci根據(jù)Gibbs-Thomson方程可表述為[19]：

(7)

式中：γ為界面能，取為0.26 J/m2，Vm為析出相單位摩爾體積，取為7.62×10-5m3/mol。

粒子處于即不增長也不溶解臨界狀態(tài)的計算公式為[18]：

(8)

為了描述完整過程的粒子尺寸分布的時間演化過程，將粒子半徑連續(xù)分布離散為一個大量的子集集合，并且子集的范圍為Δr，每一個子集由它所對應(yīng)的半徑ri和沉淀相數(shù)量密度Ni所決定，結(jié)合文獻(xiàn)[18-20]對KWN模型[21]的修正，采用Fortran77編制求解程序，求解流程圖如圖3所示，

圖3 求解流程圖Fig.3 Flow chart of numerical solution

Al-Mg-Si鋁合金材料的強(qiáng)度分為純鋁基體晶體的貢獻(xiàn)、固溶強(qiáng)化和沉淀強(qiáng)化：

σy=σ0+σss+σp

(9)

式中：σ0為純鋁基體晶體的貢獻(xiàn)，取為10 MPa；固溶強(qiáng)化σss由Mg和Si的濃度確定。

(10)

式中：kMg，kSi分別為溶質(zhì)元素Mg和Si的固溶強(qiáng)化系數(shù)；CMg，CSi分別為基體中所含Mg和Si的濃度。析出相強(qiáng)化σp由式(11)確定：

(11)

Fsh=kGbr

(12)

Fbp=2βGb2

(13)

6系鋁合金的維氏硬度HV與屈服強(qiáng)度σy存在線性關(guān)系[19-20]，

HV=Aσy+B

(14)

式中：A，B為常數(shù)。

采用Deform-3D軟件模擬攪拌摩擦焊接過程，提取有限元模型中的溫度歷史作為輸入數(shù)據(jù)，結(jié)合基于Fortran77自行編制的微觀結(jié)構(gòu)演化和力學(xué)性能預(yù)測程序，計算第二相的演化和基于第二相演化的力學(xué)性能。

2 結(jié)果與分析

圖4和圖5所示為攪拌摩擦焊接的溫度場(水平焊接進(jìn)行至90 s)和塑性變形(水平焊接進(jìn)行至135 s)。由此，可以得到較為完整的攪拌摩擦焊接過程的升溫和降溫曲線，以便于進(jìn)行后續(xù)的析出相計算。

圖4 攪拌摩擦焊接溫度場Fig.4 Temperature field in friction stir welding

圖5 攪拌摩擦焊接塑性變形Fig.5 Plastic deformation in friction stir welding

圖6所示為攪拌摩擦焊接焊后和采用人工時效處理后的橫截面力學(xué)性能分布，左側(cè)縱坐標(biāo)為硬度，右側(cè)縱坐標(biāo)為屈服強(qiáng)度。從圖6可以看出，焊后攪拌區(qū)的屈服強(qiáng)度及硬度明顯降低，主要是由于在焊接過程中，攪拌區(qū)溫度較高，大量析出相回溶，導(dǎo)致攪拌區(qū)析出相體積分?jǐn)?shù)偏低，大約為0.05%，而在焊前，母材的析出相體積分?jǐn)?shù)為0.8%，析出相的體積分?jǐn)?shù)減少是攪拌區(qū)力學(xué)性能降低的主要原因，而在焊后進(jìn)行人工時效，185 ℃保溫0.5 h和185 ℃保溫5.5 h結(jié)果差異巨大，較長時間的保溫會使攪拌區(qū)內(nèi)的析出相體積分?jǐn)?shù)增加，從而通過新增析出相的彌散增強(qiáng)使力學(xué)性能得到部分恢復(fù)，甚至可以達(dá)到母材硬度的97%以上。0.5 h保溫也可以部分回復(fù)攪拌區(qū)的力學(xué)性能，但是效果較差，主要是由于析出相析出需要時間，所以短時間保溫析出相體積分?jǐn)?shù)較低。析出相體積分?jǐn)?shù)和析出相平均尺寸如圖7所示，從圖7可以直觀看出析出相體積分?jǐn)?shù)隨焊接熱處理工藝不同所產(chǎn)生的變化。然而，在熱影響區(qū)，由于析出相平均尺寸較大，由于Orowan效應(yīng)，大尺寸的析出相顆粒在保溫過程中繼續(xù)增大，從而使焊后熱處理對熱影響區(qū)力學(xué)性能的回復(fù)不利，同樣繼續(xù)增加保溫時間對熱影響區(qū)的影響也較弱，這說明焊后人工時效可以有效回復(fù)攪拌區(qū)的力學(xué)性能，但是不能明顯改善熱影響區(qū)的力學(xué)性能。

圖6 焊后以及人工時效后截面力學(xué)性能分布Fig.6 Mechanical property distributions after weld and after artificial ageing

圖7 焊后及人工時效截面的析出相體積分?jǐn)?shù)和平均半徑 (a)體積分?jǐn)?shù)；(b)平均半徑Fig.7 Volume fractions and mean radius of precipitated phase after weld and after artificial ageing (a) volume fraction；(b) mean radius

為了說明計算結(jié)果的有效性，選取同樣材料的焊后熱處理實驗數(shù)據(jù)[18]進(jìn)行對比，如圖8所示。采用185 ℃保溫5.5 h可以使焊接核心區(qū)域的硬度達(dá)到母材的98%，與當(dāng)前材料攪拌摩擦焊接采用同樣的焊后人工時效工藝硬度回復(fù)程度相當(dāng)，回復(fù)程度的差異不超過1%，從而能夠驗證本工作模型所反映的焊后人工時效影響的正確性和有效性。如圖8(b)所示為計算攪拌摩擦焊接以及焊后熱處理演化模型得到的焊接橫截面硬度分布；圖8(a)為所參考的相同材料進(jìn)行熔焊熱循環(huán)過程實驗測試得到的硬度分布，上坐標(biāo)為熔焊過程產(chǎn)生的不同熱循環(huán)達(dá)到的最高峰值溫度沿焊縫橫截面所處的不同距離。需要說明的是，實驗參考值所展示的熱循環(huán)過程為焊縫中心處溫度以50 ℃/s從室溫升至峰值溫度并且保持60 s后再以30 ℃/s從峰值溫度冷卻至室溫的過程?？梢钥闯?，攪拌摩擦焊接數(shù)值計算中的熱循環(huán)過程與熔焊在頂峰溫度的保留時間大不相同。由此結(jié)果可知，相較于攪拌摩擦焊接，在熔焊過程中熱影響區(qū)經(jīng)歷的軟化時間更長，達(dá)60 s，宏觀上所測得到的硬度值更小。

圖8 實驗[18]和模擬的焊后熱處理對力學(xué)性能影響 (a)實驗參考；(b)數(shù)值結(jié)果Fig.8 Effect of PWHT on mechanical property in experiment [18] and numerical model (a)experimental tests;(b)numerical results

圖9為不同保溫溫度對力學(xué)性能的影響，左側(cè)縱坐標(biāo)為硬度，右側(cè)縱坐標(biāo)為屈服強(qiáng)度。采用140 ℃保溫5 h的熱處理效果明顯低于180 ℃保溫5 h的情況，隨著保溫溫度的不斷增加，攪拌區(qū)的力學(xué)性能回復(fù)明顯增加，這說明較高的保溫溫度會有利于攪拌區(qū)力學(xué)性能的回復(fù)；但是當(dāng)溫度更高時，由于析出相粒子的快速長大，會使焊后熱處理后的攪拌摩擦焊接平板力學(xué)性能變差，當(dāng)采用220 ℃保溫5 h時，攪拌區(qū)的硬度為72 HV，遠(yuǎn)低于180 ℃保溫5 h的情況，母材區(qū)的硬度也出現(xiàn)了明顯的降低，這顯然不利于攪拌摩擦焊結(jié)構(gòu)件的服役安全。

圖9 不同保溫溫度對攪拌摩擦焊接平板力學(xué)性能的影響Fig.9 Effect of holding temperature on mechanical properties of friction stir weld

3 結(jié) 論

(1)通過KWN模型建立了攪拌摩擦焊接過程中的析出相演化和力學(xué)性能預(yù)測模型，進(jìn)一步研究了焊后熱處理對攪拌摩擦焊接平板力學(xué)性能的影響機(jī)理，與實驗結(jié)果對比驗證了模型的有效性。

(2)更長的焊后保溫時間有利于攪拌區(qū)力學(xué)性能的回復(fù)，在人工時效180 ℃和5.5 h保溫可以使攪拌區(qū)力學(xué)性能回復(fù)到母材的90%。

(3)較高的保溫溫度有利于攪拌區(qū)力學(xué)性能的快速回復(fù)，但是當(dāng)溫度高于200 ℃時，長時間保溫會使母材軟化，不利于力學(xué)性能回復(fù)。

(4)通過焊后熱處理不能明顯改善熱影響區(qū)的力學(xué)性能。

[1] MISHRA R S, MA Z Y. Friction stir welding and processing[J]. Materials Science and Engineering：R,2005,50:1-78.

[2] 薛鵬，張星星，吳利輝，等. 攪拌摩擦焊接與加工研究進(jìn)展[J]. 金屬學(xué)報，2016,52(10):1222-1238.

(XUE P，ZHANG X X，WU L H，etal．Research progress on friction stir welding and processing[J]．Acta Metallurgica Sinica,2016，52(10)：1222-1238．)

[3] NANDAN R，DEBROY T，BHADESHIA H K D H．Recent advances in friction-stir welding-process, weldment structure and properties[J]. Progress in Materials Science,2008,53:980-1023.

[4] HE X C, GU F S, BALL A. A review of numerical analysis of friction stir welding[J]. Progress in Materials Science,2014,65:1-66.

[5] 張昭，張洪武. 攪拌摩擦焊的數(shù)值模擬[M]. 北京：科學(xué)出版社，2016.

[6] ZHANG Z H, LI W Y, LI J L,etal. Effective predictions of ultimate tensile strength, peak temperature and grain size of friction stir welded AA2024 alloy joints[J]. The International Journal of Advanced Manufacturing Technology,2014,73:1213-1218.

[7] FRATINI L, BUFFA G. CDRX modelling in friction stir welding of aluminium alloys[J]. International Journal of Machine Tools & Manufacture,2005,45:1188-1194.

[8] AZIZIEH M, BAHADORI R, ABBASI M,etal. Effect of friction stir processing on the microstructure of pure magnesium castings[J]. International Journal of Cast Metal Research,2015,28(6):345-351.

[9] ZHANG Z, WU Q. Numerical studies of tool diameter on strain rates, temperature rises and grain sizes in friction stir welding[J]. Journal of Mechanical Science and Technology,2015,29(10):4121-4128.

[10] GERLICH A, AVRAMOVIC-CINGARA G, NORTH T H. Stir zone microstructure and strain rate during Al 7075-T6 friction stir spot welding[J]. Metallurgical and Materials Transactions A,2006,37:2773-2786.

[11] ZHANG Z, WU Q, GRUJICIC M,etal. Monte Carlo simulation of grain growth and welding zones in friction stir welding of AA6082-T6[J]. Journal of Materials Science,2016,51:1882-1895.

[12] GRUJICIC M, RAMASWAMI S, SNIPES J S,etal. Prediction of the grain-microstructure evolution within a friction stir welding (FSW) jointviathe use of the Monte Carlo simulation method[J]. Journal of Materials Engineering and Performance, 2015, 24:3471-3486.

[13] SIMAR A, BRéCHET Y, DE MEESTER B,etal. Sequential modeling of local precipitation, strength and strain hardening in friction stir welds of an aluminum alloy 6005A-T6[J]. Acta Materialia,2007,55:6133-6143.

[14] SIDHAR H, MARTINEZ N Y, MISHRA R S,etal. Friction stir welding of Al-Mg-Li 1424 alloy[J]. Materials & Design,2016,106:146-152.

[15] LI B, SHEN Y F, LUO L,etal. Effects of processing variables and heat treatments on Al/Ti-6Al-4V interface microstructure of bimetal clad-plate fabricatedviaa novel route employing friction stir lap welding[J]. Journal of Alloys and Compounds,2016,658:904-913.

[16] HUANG Y X, WAN L, LYU Z L,etal. Microstructure and microhardness of aluminium alloy friction stir welds with heat treatment[J]. Science and Technology of Welding and Joining,2016,21:638-644.

[17] ZHANG Z, WAN Z Y. Predictions of tool forces in friction stir welding of AZ91 magnesium alloy[J].Science and Technology of Welding and Joining,2012,17(6):495-500.

[18] GRONG ?. Metallurgical modelling of welding [M]．2nd ed.London: The Institute of Materials,1997.

[19] MYHR O R, GRONG ?, ANDERSEN S J. Modelling of the age hardening behaviour of Al-Mg-Si alloys[J]．Acta Materialia,2001,49(1):65-75.

[20] MYHR O R, GRONG ?, FJ?R H G,etal. Modelling of the microstructure and strength evolution in Al-Mg-Si alloys during multistage thermal processing[J]．Acta Materialia，2004，52：4997-5008.

[21] WAGNER R, KAMPMANN R. Materials science and technology—a comprehensive treatment[M]. Weinhem: Wiley-VCH,1991.

(責(zé)任編輯：徐永祥)

NumericalSimulationofMechanicalPropertyofPostFrictionStirWeldArtificialAgeingofAluminumAlloy

WAN Zhenyu1， ZHOU Xia1， ZHANG Zhao1,2

(1．State Key Laboratory of Structural Analysis for Industrial Equipment, Department of Engineering Mechanics, Dalian University of Technology, Dalian 116024, Liaoning China；2．International Research Center for Computational Mechanics, Dalian University of Technology, Dalian 116024, Liaoning China)

KWN model was used to establish the precipitation evolution model of friction stir welding of Al-Mg-Si alloy. The yield strength was divided into three parts: the contribution from grain size, the contribution from solid solution and the contribution from the precipitations. Based on this model, the yield strength and hardness of friction stir weld was predicted. The effect of post weld artificial ageing on mechanical properties of friction stir weld was further investigated. The results indicate that longer holding time can be beneficial to the recovery of mechanical properties in the stirring zone. Higher temperature can lead to quick recovery of mechanical properties in the stirring zone, but when the holding temperature is higher than 200 ℃, longer holding time can lead the base metal softened, which is harmful to the service of friction stir welds. The mechanical property in the heat affected zone cannot be improved by post weld artificial ageing.

friction stir welding; precipitation; post weld heat treatment; artificial ageing

2017-01-09；

2017-03-27

國家自然科學(xué)基金項目(11572074； 11232003)

張昭(1979—)，男，博士，教授，研究方向為制造工藝力學(xué)，(E-mail) zhangz@dlut.edu.cn。

10.11868/j.issn.1005-5053.2016.000003

TG402

： A

： 1005-5053(2017)04-0019-06