基于計算機模擬技術的薄板GTAW焊接接頭溫度分布

楊 莎

（遵義師范學院，貴州遵義563002）

基于計算機模擬技術的薄板GTAW焊接接頭溫度分布

楊 莎

（遵義師范學院，貴州遵義563002）

使用類似K型熱電偶測量不同薄板焊縫接頭在GTAW焊接過程中熱影響區溫度的分布。對比通過有限元分析軟件建立的三維有限元模擬測量板材的溫度分布情況。由于熱電偶與熔合線接近且冷卻速度相似，證明實驗和仿真結果吻合，尤其表現在熱影響區的微觀組織上。另外，還可在焊接過程缺少填充材料時，觀察峰值溫度和不同材料的冷卻速度。

計算機模擬；GATW焊接；溫度分布

1 簡介

焊接是可靠、高效和實用的金屬填充過程，在各行業被廣泛使用，如航空航天、汽車、運輸和石油等。盡管這種填充過程具有許多優點，但會受到一些條件的影響。例如：熱周期等參數顯著影響焊縫的殘余應力、變形、微觀組織、熱影響區硬度等。導致的焊接缺陷直接影響焊縫性能。由于焊接是局部加熱過程，因此控制熱循環至關重要。

在此主要研究計算薄板焊接瞬態溫度分布。使用K-type熱電偶測量兩組相同材質和一組不同材質的GATW單道焊縫的焊接溫度分布。并運用有限元分析軟件，開發3D有限元模型，預測溫度循環，然后進行實驗模擬驗證。

2 試驗程序

試驗采用S304型不銹鋼和St37碳鋼材料，使用相同的對接焊縫接頭，焊接方法為鎢極惰性氣體保護焊（GTAW）。本試驗會將結果與有限元法獲得的結果進行比較。熱電偶位于工件水平表面中部孔上，從而測量焊縫兩邊不同位置的溫度。試件和熱電偶的位置示意如圖1所示。試件規格200mm×200mm× 3mm。在記錄測量溫度時，收集到的信號被轉移到數據記錄器和電腦上。數據記錄器設置10次/s讀數記錄熱電偶的信號。同時，使用Labview顯示溫度曲線。

圖1 熱電偶在工件表面的位置

當熱電偶與薄板連接時，熱電偶通過兩根導線接通并記錄數據，此時必須保證沒有其他導線與熱電偶導通，否則會導致測量值不準確。

工件焊接屬于沒有填充材料的單道焊，焊接保護氣為氬氣。在焊機上讀取電壓和電流。此外，記錄焊接時間，并根據時間計算出焊接速度。試驗電壓（U）、電流（I）、焊接速度（v）如表1所示。GTAW焊接的熱效率為50％，每毫米的熱輸入量為Q

表1 實驗參數和計算數據

下面將詳細討論瞬態溫度以及溫度隨距離焊縫熔合線遠近的變化情況。S304和St37化學成分如表2所示。

表2 S304不銹鋼和St37碳鋼的化學成分％

3 數值計算

在使用有限元分析軟件的基礎上，通過熱有限元計算三次單道焊接過程中溫度的分布，使用傳熱分析獲取溫度記錄。因材料的熱屬性不同，瞬態溫度場與工件的材料有關。

在GTAW焊接過程有限元模擬中，三種接頭同樣使用S304不銹鋼和St37碳鋼，模型尺寸與試驗工件相同，材料特性如表3所示，模型網格如圖2所示。在焊縫區及其附近劃分細網格。前期已有專家研究過網格尺寸大小的影響，其結果表明峰值溫度隨網格大小而變化。本次模擬中選擇劃分細網格，峰值溫度與實際值相差小于2％。DFLUX用戶子程序使用FORTRAN語言和模型來計算熱通量。建立一個近似于焊接熔池形狀和尺寸的雙橢球熱模型，模型中單個橢球內熱源分布方程為

后部分橢球內熱源分布方程為

式中x、y、z為焊接熔合線在雙橢球熱模型中的坐標；ff和fr分別為前后橢球熱量分布函數，ff+fr=2.0，本次實驗中ff=1.4，fr=0.6，可以看出在焊接過程中前半部分橢球內溫度斜率比后半部分大；Q為焊接熱源功率，其計算方法見式（1）。a1、a2、b、c為橢球形狀參數。通過實驗研究計算上述參數，從而調整焊接參數以得到理想的焊縫熔合區。

圖2 焊接模擬三維有限元模型

表3 St37碳素結構鋼和S304型不銹鋼的熱物理性質

由于在本次實驗中，GTAW焊接沒有使用焊縫填充金屬，其中一些焊接參數不再使用，如起弧、收弧等。當焊接電弧應用于工件，整個焊縫存在并且經過加熱。因此，上述技術的應用在仿真中引起了誤差，在這項工作的整個模擬時間模型中，移動熱量源靠近焊縫元件。

焊接中瞬態傳熱分析控制方程為

非線性各向同性傅里葉熱通量方程為

式中k是與溫度有關的熱導率。

在焊接過程中，熱輻射和焊縫表面的對流會引發熱量流失。在焊縫附近高溫區主要是熱輻射損失，遠離焊縫的低溫區主要是熱對流導致的熱量損失。因此使用聯合導熱邊界條件，從而避免輻射建模的難度。同樣，在模型中使用熱對流系數。對流表面傳熱系數如表4所示。

表4 對流表面傳熱系數

4 結果和討論

4.1 St37碳鋼對接焊縫

沿著焊接方向的等溫線和焊縫接頭的有限元仿真模型如圖3所示。工件表面上離焊縫熔合線距離不同的各個點的溫度分布如圖4所示。其中距離熔合線3mm處的溫度最高，達到了550°C。由圖4可知，隨著距離的增加，溫度出現非線性趨勢的降低。其原因是焊槍的局部加熱和材料的熱性能與溫度的變化是非線性關系。圖5是上述各點實驗結果和有限元模擬結果對比的情況。兩種方法結果吻合良好，區別在于實驗初始階段工件溫度上升更快。這可能是由于試驗會有熔合物飛濺以及仿真模型熱源與實際焊接不同造成的。焊縫熔線周圍各點的峰值溫度如圖6所示。

圖3 仿真模型中熔合線的等溫線

4.2 S304型不銹鋼對接焊縫

通過實驗測量焊縫接頭左邊薄板的各點溫度，分布如圖7所示。其最高溫度約470°C。對比圖4和圖7的結果，由于不同材料的熱屬性（熱導率）不同，St37比S304冷卻更快。實驗和仿真模擬過程中焊縫熔線周圍各點的峰值溫度如圖8所示。

圖4 E1實驗各點的溫度分布

圖5 實驗各點溫度和仿真模擬結果對比

圖6 E1實驗和仿真模擬各點的峰值溫度

圖7 E2實驗各點的溫度分布情況

圖8 E2實驗和仿真模擬各點的峰值溫度

圖9 E3實驗中St37各點的溫度分布情況

圖10 E3實驗中S304各點的溫度分布情況

4.3 不同金屬的對接焊縫

由St37碳鋼和S304型不銹鋼兩種不同材料的薄板組成的對接焊縫溫度分布如圖9、圖10所示。可以看出S304的溫度高于St37。實驗測量和有限元模擬結果如圖11、圖12所示。圖13為兩種材料的冷卻曲線，明顯S304比St37冷卻更快。在焊接時工件中沒有填充材料，同時考慮到焊接的對稱性，得出兩種材料薄板的峰值溫度和冷卻速率主要取決于它們的熱屬性，特別是熱導率。峰值溫度出現差異的另一個原因是由于St37碳鋼受到相變的影響，因為St37碳鋼固態相變時會吸收一定比例的熱源能量。

圖11 St37碳鋼實驗和仿真模擬時各點的溫度對比

圖12 S304碳鋼實驗和仿真模擬時各點的溫度對比

圖13 兩種材料在實驗和仿真模擬中的冷卻曲線

5 結論

使用S304型不銹鋼和St37型碳鋼薄板進行了三組GTAW焊接實驗，通過三維有限元進行計算機模擬，并將實驗測量的溫度分布和仿真模擬的結果進行對比。其結果表明該種溫度模型在焊接過程中具有良好的預測能力。

（1）各點的焊接峰值溫度分布和溫度降低速度與其到熔合線的距離呈非線性關系。通過預測冷卻速率與距離的關系，可用于熱影響區微觀結構的預測。

（2）由于熱傳導率不同，在焊接過程中，在焊縫熔合線周圍S304的峰值溫度高于St37。

[1]高侖.鎂合金成形技術的開發與應用[J].輕合金加工技術，2004（3）：5-12.

[2]邢靜忠，李軍.ANSYS的建模方法和網格劃分[J].中國水運（學術版），2006，6（9）：16-18.

[3]何良菊,李培杰.中國鎂工業現狀與鎂合金開發技術[J].鑄造技術，2003，24（3）：161-162.

[4]Kou S，Le Y.Improveweldqualitybylow frequencyarcoscillation[J].Welding Journal，1985，64（3）：51-55.

Tem perature distribution of sheet GTAW welding joint based on computer simulation technology

YANG Sha
（ZunyiNormalCollege，Zunyi 563002，China）

Use similar type K thermocouple tomeasure different sheetweld joint in heat affected zone temperature distribution in the process of GTAW welding.The temperature distribution of three dimensional finite element simulation ofmeasuring plate at the same time by comparing the finite element analysis software to establish.The thermocouple and the fusion line near and cooling rate are similar，prove that the experimental and simulation results agree，especially themicrostructure in the heat affected zone.In addition，also the lack of fillingmaterial in the welding process，the observation of the peak temperature and cooling rate of differentmaterials.

computer simulation；GATW welding；temperature distribution

TG409

A

1001-2303（2015）07-0101-06

10.7512/j.issn.1001-2303.2015.07.22

2014-11-17；

2015-07-07

楊莎（1982—），女，貴州銅仁人，講師，碩士，主要從事電子技術和網絡應用技術的研究。