AFAM 7050鋁合金的數值模擬

劉清賢，楊濱，東青，孫獎

1北京科技大學，北京，100083；2先進技術研究院，西雅圖，98033

引言

攪拌摩擦焊（FSW）是英國焊接技術研究所（TWI）于1991年發明的一種固相連接方法。由于它具有力學性能優異、效率高、污染小等諸多優點，適合鋁鎂等有色金屬及其合金的焊接，自發明初就得到了人們廣泛關注和深入研究。

攪拌摩擦焊增材制造（AFAM）是基于現有的攪拌摩擦焊技術，將其機理延伸和拓展到增材制造領域的新技術。在增材制造中，殘余應力嚴重影響材料的疲勞壽命和穩定性，同時溫度場是影響材料應力場的重要因素。由于溫度難以準確實時監測，無法精準評估其對應力場的影響，因此有必要對溫度場和殘余應力進行數值模擬。張昭[1]等人提出了雙軸肩攪拌頭產熱功率的計算方法，實現了對攪拌摩擦焊接T型接頭殘余狀態的預測。Khandkar[2]等人模擬了6061和2024兩種鋁合金的攪拌摩擦焊接后的殘余應力，實驗結果趨勢和Sutton[3]等人的大致吻合。本文基于完全熱力耦合模型，建立了熱輸入功率和轉速之間的關系，實現了溫度場的模擬，并在溫度場的基礎上進行了殘余應力的模擬。

1 數值模擬建模

1.1 有限元模型

采用八節點實體導熱單元solid70，基板尺寸為150×40×8mm，增材尺寸為150×10×8mm，攪拌頭軸肩直徑20mm，進行不均勻網格劃分，共劃分網格29899個，單元25600個。

1.2 材料的熱物性參數

選取實驗材料為7050鋁合金，基板尺寸150×40×8mm，增材寬度為10mm.攪拌頭為鎢鋼，實驗材料的具體熱物性參數如見表1：

泊松比為0.33，

表1 7050鋁合金不同溫度下熱物性參數

攪拌頭轉速750r/min，行進速度60 mm/min。

1.3 熱源模型

攪拌摩擦焊的主要熱源是摩擦生熱，由于AFAM中沒有攪拌針，故軸肩成為主要熱源，因此采用面熱源模型。熱流密度可以表示為：

其中R為軸肩半徑，0＜r＜R，Q為熱輸入功率。其中，扭矩和轉速之間關系如下式[4]：

其中w為攪拌頭轉速。

1.4 邊界條件

假設攪拌頭和基板結合面不發生熱量交換，其他和空氣接觸的面的換熱系數18 W＊（m2＊℃），墊板和空氣的換熱系數1800 W＊（m2＊℃），環境溫度10℃。

1.5 殘余應力模擬

進行殘余應力模擬前需要載入溫度場文件，然后施加載荷進行求解。由于焊接屬于大變形，在設定分析選項時，需要打開大變形選項。用Newton-Rahpson迭代法進行求解。

圖1 時間為6，80，156s時溫度場云圖

2 結果與討論

2.1 溫度場模擬結果

在攪拌頭轉速為750 r/min，前進速度為60 mm/min條件下，攪拌頭預熱時間6s，焊接時間150s，冷卻時間154s。焊接結束時溫度場如下圖所示。由圖1可知，溫度的分布呈現橢球形，最高溫度范圍在軸肩內，表明了軸肩是增材制造的主要熱源。

整個焊接過程可以分為：

（1）初始升溫過程，此時峰值溫度比較預熱結束約上升100℃；

（2）增材制造達到穩定后，峰值溫度隨之穩定；

（3）增材制造結束，溫度有所上升。結合實際增材制造過程，這是由于增材制造過程趨近于結束時，墊板散熱面積較初始下降，而產熱量不變，造成單位時間內產熱大于散熱，從而推動增材制造結束時溫度上升。同時，在接近結束時，墊板本身也存在熱量累積，也會導致溫度上升。

圖2 不同位置溫度變化曲線圖

圖3 不同位置殘余應力曲線圖

圖2為焊縫中心線方向上表面距離初始位置30、60、90、120mm處各點溫度隨時間變化情況。由圖2可知，隨攪拌頭向前移動，各處溫度依次達到峰值，然后迅速下降。

2.2 殘余應力模擬結果

圖3是增材樣件上表面垂直于焊縫中心線處縱向、橫向以及厚度方向殘余應力的數值模擬變化曲線。由圖3可知，厚度方向的殘余應力幾乎為0，可以忽略，主要殘余應力為縱向殘余應力。造成這一規律的主要原因是在增材制造過程中，鋁合金同時受到溫度場和攪拌頭軸肩的應力作用，當軸肩攪拌到此處時，產生熱量的擴散，由于溫度的下降容易發生材料的收縮，同時周圍的材料約束作用進而產生了縱向殘余應力。在增材制造過程中沒有發生相變，這種不穩定狀態隨著冷卻保留到室溫狀態。

從樣件上表面到下表面殘余應力逐漸遞減，從拉應力轉變為壓應力。此規律和溫度場的分布基本相同，驗證了溫度場產生的熱應力是增材制造的殘余應力的主要影響因素。由于單位長度的焊縫的熱輸入量和焊接速度二者呈負相關，即焊接速度減小，單位長度的熱輸入量增大[5]。利用此特性，可以有效減小增材制造過程中的局部區域受熱不均勻，從而減小增材制造的殘余應力。

3 結論

1.在轉速750 r/min，行進速度60mm/min時，7050鋁合金AFAM的峰值溫度在500℃范圍內，在固相線以下；

2.殘余應力主要為縱向殘余應力，在軸肩以內為拉應力，以外為壓應力；厚度方向數值很小，幾乎不變；

3.熱應力是殘余應力的主要影響因素，殘余應力的分布與熱應力有關。

參考文獻

[1]張昭，劉會杰.攪拌頭形狀對攪拌摩擦焊材料變形和溫度場的影響[J].焊接學報，2011，32（3）：5-8.

[2]Khandkar M Z H，Khan J A，Reynolds A P，et al. Predicting residual thermal stresses in friction stir welded metals[J]. Journal of Materials Processing Tech，2006，174（1-3）：195-203.

[3]Sutton M A，Reynolds A P，Wang D Q，et al. A Study of Residual Stresses and Microstructure in 2024-T3 Aluminum Friction Stir Butt Welds[J]. Journal of Engineering Materials &Technology，2002，124（4）：215-221.

[4]張正偉. 攪拌摩擦焊接構件殘余狀態和疲勞壽命研究[D]. 大連理工大學，2014.

[5]王大勇，馮吉才，王攀峰. 攪拌摩擦焊接熱輸入數值模型[J]. 焊接學報，2005，26（3）：25-28.