張廣剛， 趙強(qiáng)， 劉佳濤， 白瑞祥， 董紅剛

(1.大連理工大學(xué)，遼寧 大連 116024；2.中國(guó)航發(fā)沈陽(yáng)黎明航空發(fā)動(dòng)機(jī)有限責(zé)任公司，沈陽(yáng) 110043)

0 前言

慣性摩擦焊通過飛輪儲(chǔ)能，利用材料接觸面相對(duì)運(yùn)動(dòng)產(chǎn)生的摩擦熱使接頭發(fā)生塑性變形，在頂鍛力的作用下使工件達(dá)到原子間結(jié)合的精密焊接，是一種高效、環(huán)保、低熱輸入、易獲得優(yōu)質(zhì)接頭的固相連接方法[1-3]。FGH96高溫合金是基于損傷容限思想設(shè)計(jì)的第二代粉末高溫合金，在750 ℃下仍能保持優(yōu)異的服役性能，廣泛應(yīng)用于高性能發(fā)動(dòng)機(jī)渦輪盤和其它熱端部件的制造[4-6]。由于FGH96高溫合金的化學(xué)成分比較復(fù)雜，焊接區(qū)域的組織和性能對(duì)焊接溫度比較敏感，而在焊接過程中對(duì)溫度梯度和應(yīng)力進(jìn)行測(cè)量較為困難，因此采用數(shù)值模擬方法研究不同工藝參數(shù)對(duì)FGH96高溫合金慣性摩擦焊溫度場(chǎng)和應(yīng)力分布的影響，對(duì)于控制焊接過程和提升焊接質(zhì)量具有重要的工程價(jià)值。

Zhang等人[7]建立了二維熱-力耦合有限元模型，研究了GH4169高溫合金慣性摩擦焊接過程中的溫度場(chǎng)分布，計(jì)算得到的溫度曲線與熱電偶測(cè)得的數(shù)據(jù)吻合較好。朱海等人[8]建立了二維GH4169高溫合金軸對(duì)稱模型，分析了焊接條件對(duì)慣性摩擦焊溫度場(chǎng)和縮短量的影響，指出增大初始轉(zhuǎn)速、頂鍛力和轉(zhuǎn)動(dòng)慣量均能增加軸向縮短量，其中初始轉(zhuǎn)速對(duì)軸向縮短量的影響最為顯著。Wang等人[9]建立了鎳基高溫合金二維軸對(duì)稱模型，發(fā)現(xiàn)了軸向縮短量與初始轉(zhuǎn)速的平方成正比，與頂鍛力成對(duì)數(shù)的關(guān)系，同時(shí)指出軸向壓力對(duì)高溫寬度影響更加明顯。李文亞等人[10]基于ABAQUS/Explicit建立了三維GH4169高溫合金管剛塑性體有限元模型，計(jì)算了不同時(shí)刻飛邊的形貌和溫度場(chǎng)分布。卜文德等人[11]建立了Inconel 718合金三維熱-力耦合有限元模型，分析了慣性摩擦焊過程中溫度場(chǎng)變化特點(diǎn)，模擬所得溫度曲線與試驗(yàn)實(shí)測(cè)數(shù)據(jù)基本吻合。張利國(guó)等人[12]進(jìn)行了FGH96高溫合金慣性摩擦焊過程中飛邊形成規(guī)律的研究，指出增加初始轉(zhuǎn)速和頂鍛力均可增大飛邊尺寸和卷曲角度。鑒于計(jì)算機(jī)運(yùn)算能力的限制，以往對(duì)慣性摩擦焊的模擬多采用二維軸對(duì)稱或三維剛塑性體模型，這種簡(jiǎn)化模型在一定程度上可以揭示焊接過程中的熱-力耦合行為，但由于工件在焊接過程中會(huì)發(fā)生較大的塑性變形，這種簡(jiǎn)化模型相較于雙塑性體模型并不能很好的反映摩擦界面的熱力相互作用，隨著計(jì)算機(jī)運(yùn)算能力的提高，采用雙塑性體模型更符合實(shí)際的慣性摩擦焊過程。

文中基于ABAQUS軟件建立了FGH96高溫合金管慣性摩擦焊的三維雙塑性體模型，研究了慣性摩擦焊接過程中溫度場(chǎng)與軸向應(yīng)力場(chǎng)的演變規(guī)律及初始轉(zhuǎn)速、頂鍛力和轉(zhuǎn)動(dòng)慣量對(duì)慣性摩擦焊接過程的影響規(guī)律。

1 有限元模型的建立

文中所采用的FGH96高溫合金管幾何模型及網(wǎng)格劃分如圖1所示。工件單管高度為15 mm，外徑60 mm，壁厚5 mm。由于摩擦界面附近存在較大的溫度梯度和塑性變形，所以對(duì)摩擦界面兩側(cè)5 mm區(qū)域進(jìn)行網(wǎng)格細(xì)化，單元尺寸設(shè)定為0.7 mm，遠(yuǎn)離摩擦界面的非重點(diǎn)區(qū)域采用梯度化網(wǎng)格過渡，網(wǎng)格類型為使用減縮積分及沙漏控制的C3D8RT八節(jié)點(diǎn)熱力耦合六面體網(wǎng)格，整個(gè)焊接模型共71 820個(gè)單元。為避免計(jì)算過程中網(wǎng)格畸變導(dǎo)致不收斂，在界面附近采用網(wǎng)格自適應(yīng)技術(shù)對(duì)網(wǎng)格進(jìn)行重新劃分。為引入慣性摩擦焊中飛輪的動(dòng)能，將旋轉(zhuǎn)端工件上部2 mm位置耦合至參考點(diǎn)(RP)，賦予RP點(diǎn)轉(zhuǎn)動(dòng)慣量和初始轉(zhuǎn)速以完成慣性摩擦焊接過程中的能量輸入。由于模擬的準(zhǔn)確性受材料本構(gòu)方程影響較大，所以選擇合適的本構(gòu)方程至關(guān)重要。Geng等人[13]分別采用Field-Backofen、Johnson-Cook和Arrhenius 3種材料本構(gòu)模型進(jìn)行了GH4169高溫合金線性摩擦焊飛邊形貌數(shù)值模擬的對(duì)照試驗(yàn)，結(jié)果表明采用雙曲正弦修正的Arrhenius方程所得到的計(jì)算結(jié)果與試驗(yàn)結(jié)果吻合較好，因此文中采用應(yīng)變補(bǔ)償型Arrhenius本構(gòu)方程。材料本構(gòu)方程及所涉及的相關(guān)熱物性參數(shù)，如密度、彈性模量、線膨脹系數(shù)和熱導(dǎo)率等均借鑒于文獻(xiàn)[14-15]。

圖1 工件的幾何模型與網(wǎng)格劃分

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

2.1慣性摩擦焊接過程

為描述慣性摩擦焊過程中溫度場(chǎng)和軸向應(yīng)力場(chǎng)的演變規(guī)律，根據(jù)經(jīng)驗(yàn)及參考文獻(xiàn)[16]，將焊接工藝參數(shù)初步設(shè)定為：初始轉(zhuǎn)速n=600 r/min、頂鍛力P=400 MPa和轉(zhuǎn)動(dòng)慣量I=2.4 kg·m2。數(shù)值計(jì)算得到的飛邊形貌與實(shí)際焊件形貌對(duì)比如圖2所示，實(shí)際FGH96高溫合金慣性摩擦焊接接頭的飛邊尺寸約2.0 mm，計(jì)算得到的飛邊尺寸為2.1 mm，誤差僅為5%，模擬結(jié)果與試驗(yàn)焊接件吻合較好，表明該建模方法符合焊接實(shí)際情況。

圖2 飛邊形貌對(duì)比

焊接過程中溫度場(chǎng)的分布如圖3所示，溫度場(chǎng)關(guān)于摩擦界面呈對(duì)稱分布。根據(jù)接合界面溫度場(chǎng)的變化情況，整個(gè)焊接過程可劃分為3個(gè)階段：① 0 ～0.7 s為初始摩擦階段；② 0.7～2.1 s為穩(wěn)定摩擦階段；③ 2.1～4.0 s為保壓冷卻階段。初始摩擦階段，界面中心溫度在0.1 s內(nèi)迅速升至1 226 ℃，這是由于FGH96高溫合金在常溫下抵抗變形的能力較強(qiáng)，工件在轉(zhuǎn)動(dòng)過程中，接合界面剪切粒子通過干摩擦產(chǎn)生大量熱，界面中心區(qū)域急速升溫，使得界面溫度徑向分布不均勻，呈中心溫度高、邊緣溫度低，如圖3a所示。隨著摩擦產(chǎn)熱不斷積累，高溫區(qū)域逐漸擴(kuò)展至整個(gè)接頭。0.7 s時(shí)，接頭高溫粘塑性金屬開始被擠出，形成飛邊，如圖3b所示，此時(shí)進(jìn)入穩(wěn)定摩擦階段，界面中心維持在1 335 ℃高溫平臺(tái)，持續(xù)高溫使得接合界面大量粘塑性金屬不斷被擠出從而使飛邊進(jìn)一步擴(kuò)展，如圖3c所示。隨后飛輪轉(zhuǎn)速下降，摩擦產(chǎn)熱速率低于散熱速率，界面溫度開始下降，至2.1 s時(shí)摩擦結(jié)束，此時(shí)界面最高溫度為1 127 ℃，如圖3d所示。工件停止轉(zhuǎn)動(dòng)后即進(jìn)入保壓冷卻階段，該階段已無(wú)熱輸入，熱量一部分輻射到環(huán)境中，一部分通過熱傳導(dǎo)向工件冷端傳遞，接頭溫度持續(xù)下降，由圖3e和圖3f可觀察到，軸向溫度梯度變小， 界面溫度分布逐漸均勻化。

圖3 焊接過程溫度場(chǎng)分布云圖

圖4為中心軸線上不同位置點(diǎn)的溫度變化曲線，取點(diǎn)位置分別距離摩擦界面0.7 mm，2.1 mm，3.5 mm和4.9 mm。由a曲線可以看出，節(jié)點(diǎn)升溫速率較快，溫度在0.5 s內(nèi)可升至1 000 ℃。此時(shí)界面溫度可達(dá)1 318℃，摩擦界面上逐漸形成一層高溫塑性金屬，界面材料流動(dòng)性增強(qiáng)，摩擦扭矩減小，產(chǎn)熱功率降低，從而導(dǎo)致升溫速率減緩。在1.0～1.8 s之間，摩擦產(chǎn)生的熱量與通過熱傳遞及飛邊擠出帶走的熱量達(dá)到動(dòng)態(tài)平衡，a點(diǎn)溫度維持在1 200 ℃。隨后由于轉(zhuǎn)速下降，摩擦產(chǎn)熱不足，接頭溫度開始下降，至2.1 s時(shí)飛輪能量耗盡，進(jìn)入保壓冷卻階段。隨著離摩擦界面距離的增大，節(jié)點(diǎn)升溫速率明顯變緩，且峰值溫度逐漸降低。這是由于a點(diǎn)距離摩擦界面近，產(chǎn)熱機(jī)制主要為摩擦生熱，所以其升溫速率較快，峰值溫度更高，達(dá)到峰值溫度所用時(shí)間也最短[17]。對(duì)于b點(diǎn)、c點(diǎn)和d點(diǎn)而言，升溫速率主要受材料熱導(dǎo)率的影響。由于傳熱的遲滯效應(yīng)和熱量耗散，軸向各節(jié)點(diǎn)的峰值溫度和升溫速率均隨距離增大而顯著降低，且距離摩擦界面越遠(yuǎn)的位置達(dá)到最高溫度所用的時(shí)間越長(zhǎng)，溫度梯度越小。在保壓冷卻過程中，4個(gè)節(jié)點(diǎn)的溫差逐漸減小，接頭溫度場(chǎng)趨于一致。4.9 mm位置的d點(diǎn)最高溫度為628 ℃，低于FGH96高溫合金的服役溫度750 ℃，由此可見焊接過程對(duì)FGH96高溫合金母材組織和性能影響較小。

圖4 中心軸線上特征點(diǎn)位置溫度曲線

焊接過程中工件軸向應(yīng)力分布隨時(shí)間的演變過程如圖5所示。由于旋轉(zhuǎn)端約束條件不同于滑移端，旋轉(zhuǎn)端的應(yīng)力值略大于相同位置的滑移端，單側(cè)工件應(yīng)力分布不均勻，旋轉(zhuǎn)端與滑移端的軸向應(yīng)力關(guān)于摩擦界面呈近似對(duì)稱分布。遠(yuǎn)離界面位置軸向應(yīng)力為正值，表示其處于拉應(yīng)力狀態(tài)，工件周邊軸向應(yīng)力為負(fù)值，表示其處于壓應(yīng)力狀態(tài)。摩擦界面中心位置附近軸向壓應(yīng)力值最大，應(yīng)力值從中心到邊緣逐漸降低。初始摩擦階段，由于頂鍛力的作用，摩擦界面主要受壓應(yīng)力。0.7 s時(shí)高溫區(qū)域擴(kuò)展至整個(gè)摩擦界面，接頭屈服強(qiáng)度降低，生成的粘塑性金屬流動(dòng)性增強(qiáng)，在頂鍛力作用下被擠出形成飛邊。隨著飛邊量增加，軸向壓力沿飛邊擠出方向發(fā)生分解，界面粘塑性金屬受力發(fā)生改變，界面邊緣應(yīng)力得到釋放，其應(yīng)力狀態(tài)也隨之轉(zhuǎn)變?yōu)槔瓚?yīng)力，隨著焊接過程進(jìn)行，拉應(yīng)力值逐漸從186 MPa降至77 MPa。隨著轉(zhuǎn)速下降，界面升溫速率低于散熱速率，接頭溫度開始下降，界面拉應(yīng)力逐漸升高。保壓冷卻階段，溫度場(chǎng)逐漸均勻化，接頭軸向壓應(yīng)力值減小，飛邊根部則由于受到擠壓，形成較高的壓應(yīng)力集中區(qū)。

圖5 焊接過程軸向應(yīng)力云圖

2.2初始轉(zhuǎn)速對(duì)焊接過程的影響

為研究初始轉(zhuǎn)速對(duì)界面峰值溫度的影響規(guī)律，將轉(zhuǎn)動(dòng)慣量設(shè)為I=2.4 kg·m2，頂鍛力設(shè)為P=350 MPa，初始轉(zhuǎn)速分別設(shè)為550 r/min，600 r/min和650 r/min，對(duì)應(yīng)的工件停止轉(zhuǎn)動(dòng)時(shí)間分別為1.80 s，2.35 s和2.55 s。不同轉(zhuǎn)速下距離摩擦界面0.7 mm位置的溫度變化曲線如圖6所示。可以看出在不同初始轉(zhuǎn)速下，界面中心溫度變化趨勢(shì)類似，均為快速升溫后保持穩(wěn)定。在550 r/min時(shí)，高溫持續(xù)時(shí)間為1.2 s，曲線峰值溫度也隨之降至1 130 ℃。當(dāng)初始轉(zhuǎn)速為650 r/min時(shí)，整個(gè)焊接過程中超過1 000 ℃的高溫持續(xù)時(shí)間延長(zhǎng)至2.1 s，曲線峰值溫度可達(dá)1 218 ℃；即隨著初始轉(zhuǎn)速的提高，高溫持續(xù)時(shí)間和峰值溫度都有所增加。這是因?yàn)樵谙嗤霓D(zhuǎn)動(dòng)慣量下，初始轉(zhuǎn)速越大，飛輪儲(chǔ)能越多，高溫持續(xù)時(shí)間越長(zhǎng)。同時(shí)，在較高的初始轉(zhuǎn)速下，兩側(cè)工件相對(duì)運(yùn)動(dòng)速度更大，摩擦產(chǎn)熱效率有所增加，從而其峰值溫度更高。

圖6 不同轉(zhuǎn)速特征點(diǎn)溫度曲線

不同初始轉(zhuǎn)速下形成的焊接接頭飛邊形貌如圖7所示。當(dāng)初始轉(zhuǎn)速為550 r/min時(shí)，界面溫度為1 070 ℃，飛邊量較少；當(dāng)初始轉(zhuǎn)速為650 r/min時(shí)，界面溫度為1 142 ℃，飛邊量明顯增加。當(dāng)初始轉(zhuǎn)速較低時(shí)，熱輸入不足，材料流動(dòng)性差，界面材料達(dá)到粘塑性狀態(tài)的量較少。在實(shí)際焊接過程中，飛邊擠出量不足，會(huì)使得工件表面金屬氧化物難以清除，影響接頭質(zhì)量。當(dāng)初始轉(zhuǎn)速升高時(shí)，熱輸入增加，摩擦界面將生成足夠多的粘塑性金屬，飛邊擠出量增多，使得界面金屬氧化物完全清除。但當(dāng)轉(zhuǎn)速過大時(shí)會(huì)導(dǎo)致熱輸入過大，熱影響區(qū)擴(kuò)大，從而導(dǎo)致工件性能下降，并且較大的熱輸入還會(huì)使得工件過度縮短，不利于后續(xù)加工[18]。

圖7 不同轉(zhuǎn)速摩擦停止時(shí)刻的飛邊形貌

2.3頂鍛力對(duì)焊接過程的影響

為研究頂鍛力對(duì)接頭飛邊形貌的影響規(guī)律，將初始轉(zhuǎn)速和轉(zhuǎn)動(dòng)慣量分別設(shè)為n=600 r/min和I=2.4kg·m2，頂鍛力以50 MPa為增量從300～450 MPa變化，摩擦結(jié)束時(shí)間分別為2.24 s，2.13 s，2.05 s和2.00 s。不同轉(zhuǎn)速下距離摩擦界面1.4 mm位置的溫度曲線如圖8所示?？梢钥闯?，不同頂鍛力下各節(jié)點(diǎn)溫度曲線變化一致。當(dāng)壓力較低時(shí)，節(jié)點(diǎn)達(dá)到峰值溫度所用時(shí)間較長(zhǎng)，高溫持續(xù)時(shí)間較短；當(dāng)壓力增加時(shí)，界面摩擦效率明顯增大，更快達(dá)到峰值溫度，峰值溫度也隨之變高。

圖8 不同頂鍛力定點(diǎn)位置溫度曲線

不同頂鍛力下接頭飛邊形貌如圖9所示。當(dāng)頂鍛力為300 MPa時(shí)，飛邊擠出較少，不利于接頭氧化物的自清理。隨著頂鍛力的增大，摩擦界面及其附近區(qū)域的金屬流動(dòng)性增強(qiáng)，飛邊擠出量逐漸增多。當(dāng)頂鍛力為450 MPa時(shí)，產(chǎn)生的飛邊形貌不規(guī)則，旋轉(zhuǎn)端的飛邊根部位置出現(xiàn)應(yīng)力集中，在該區(qū)域形成環(huán)形凹陷；材料的屈服強(qiáng)度隨著溫度升高而降低，在該頂鍛力作用下，工件根部的應(yīng)力集中區(qū)域會(huì)產(chǎn)生屈服變形。對(duì)于管狀焊件而言，飛邊有內(nèi)外兩部分，內(nèi)飛邊不易去除，導(dǎo)致服役時(shí)內(nèi)飛邊殘余拉應(yīng)力較高的區(qū)域容易開裂，因此適當(dāng)?shù)膲毫?duì)形成良好性能的接頭具有重要意義[19]。

圖9 不同頂鍛力條件下2.00 s時(shí)的飛邊形貌

2.4轉(zhuǎn)動(dòng)慣量對(duì)焊接過程的影響

為研究轉(zhuǎn)動(dòng)慣量對(duì)接頭飛邊形貌的影響規(guī)律，將初始轉(zhuǎn)速和頂鍛力設(shè)分別為n=600 r/min和P=400 MPa，轉(zhuǎn)動(dòng)慣量設(shè)置為1.8 kg·m2，2.4 kg·m2，2.8 kg·m2和3.2 kg·m2，計(jì)算得到的工件停止轉(zhuǎn)動(dòng)時(shí)間分別為1.48 s，2.06 s，2.29 s和2.50 s，接頭飛邊形貌如圖10所示。由于轉(zhuǎn)動(dòng)慣量較小，摩擦階段持續(xù)時(shí)間短，接頭熱輸入不足，此時(shí)界面最高溫度僅為895 ℃，導(dǎo)致飛邊較小，不利于接頭進(jìn)行自清潔；當(dāng)轉(zhuǎn)動(dòng)慣量增加，界面摩擦產(chǎn)熱時(shí)間變長(zhǎng)，熱輸入增加，相同時(shí)刻接頭界面的峰值溫度升高，飛邊形成量顯著增大。

圖10 不同轉(zhuǎn)動(dòng)慣量條件下2.00 s時(shí)的飛邊形貌

3 結(jié)論

(1)建立了三維雙塑性體慣性摩擦焊有限元模型，采用自適應(yīng)技術(shù)解決網(wǎng)格畸變問題。FGH96高溫合金管在慣性摩擦焊過程中，摩擦表面中心處加熱功率最大，峰值溫度可達(dá)1 335 ℃，軸向溫度梯度較大。在保壓冷卻過程中，接頭軸向溫度梯度減小，溫度場(chǎng)逐漸均勻化。

(2)工件兩側(cè)軸向應(yīng)力呈對(duì)稱分布，摩擦過程界面軸向壓應(yīng)力變化復(fù)雜，局部軸向壓應(yīng)力可達(dá)900 MPa，摩擦界面附近的壓應(yīng)力值從中心向邊緣位置減小，擠出的飛邊應(yīng)力狀態(tài)從壓應(yīng)力變?yōu)槔瓚?yīng)力。降溫過程中，飛邊末端出現(xiàn)拉應(yīng)力集中，飛邊根部出現(xiàn)壓應(yīng)力集中。

(3)增加初始轉(zhuǎn)速、頂鍛力和轉(zhuǎn)動(dòng)慣量均能增大飛邊量。提高初始轉(zhuǎn)速和轉(zhuǎn)動(dòng)慣量可以增加焊接熱輸入，延長(zhǎng)摩擦?xí)r間，其中轉(zhuǎn)速增加可以提高界面峰值溫度；增大頂鍛力能提高能量轉(zhuǎn)換效率，縮短摩擦?xí)r間，增加飛卷曲程度。