劉向臻，劉龍杰

（200093 上海市 上海理工大學(xué) 機(jī)械工程學(xué)院）

0 引言

焊接是一種非穩(wěn)定的局部加熱過程，焊接過程中的高溫加熱使焊縫區(qū)域的金屬迅速熔化，鄰近焊縫的金屬受熱產(chǎn)生不均勻膨脹。高溫作用下，焊接過程中的母材金屬屈服強(qiáng)度下降，焊縫附近母材發(fā)生塑性變形。焊接完成后，熔池冷卻，焊縫開始收縮，已經(jīng)發(fā)生塑性變形的部位無法恢復(fù)到原來的狀態(tài)。冷卻到室溫后，焊縫冷卻收縮造成焊縫附近母材縱向壓縮，焊縫縱向拉伸，這種無法恢復(fù)的變形造成焊接構(gòu)件內(nèi)部存在大量無法釋放的應(yīng)力，稱為結(jié)構(gòu)的殘余應(yīng)力。焊接完成后的殘余應(yīng)力是影響工件性能的主要原因之一，要研究和預(yù)測焊接結(jié)構(gòu)在焊接完成后的焊接變型以及焊接完成后的殘余應(yīng)力，其前提必然是對焊接過程中的瞬態(tài)溫度場有充分準(zhǔn)確的描述和計算。在進(jìn)行焊接過程的仿真分析時，熱源模型的選擇是進(jìn)行瞬態(tài)溫度場分析的基礎(chǔ)，不同熱源模型作為輸入會直接影響輸出的焊接模擬結(jié)果。本文以T 型焊接樣件為研究對象，以包括半球狀熱源、旋轉(zhuǎn)高斯體熱源、雙橢球熱源在內(nèi)3 種不同形式的體熱源為研究重點(diǎn)，分析比對了不同體熱源焊接模擬結(jié)果的準(zhǔn)確性，能夠進(jìn)一步指導(dǎo)焊接殘余應(yīng)力的分析。

目前，在進(jìn)行焊接仿真分析過程中，有多種熱源模型可供選擇，其中有Pavelic V[1]分析提出的高斯表面熱源模型，Goldak J[2]通過驗證提出的雙橢球熱源模型等。這些在焊接仿真過程中構(gòu)建的熱源加載形式，一方面能將熱能以熱流密度分布的形式限定其分布函數(shù)的區(qū)域，以模擬熔池和焊縫的形態(tài)，另一方面能以熱傳導(dǎo)的方式將熱量傳遞到函數(shù)區(qū)域以外。并且當(dāng)材料的熱物理參數(shù)發(fā)生變化后，還可以此為依據(jù)分析材料的熱物理性質(zhì)和隨溫度變化的相變潛熱。

國內(nèi)外學(xué)者以上述不同形式的熱源模型為依據(jù)，也對焊接數(shù)值模擬進(jìn)行了進(jìn)一步分析。孟慶國[3]等選取雙橢球熱源模型作為熱量輸入，對低碳鋼多層多道焊過程進(jìn)行了模擬，對溫度場進(jìn)行了分析后證明結(jié)果符合預(yù)期；華一品[4]在建立了T型接頭焊接有限元分析數(shù)值計算模型的基礎(chǔ)上，以T 型接頭埋弧焊為研究對象，計算了瞬時溫度場和殘余應(yīng)力場；Flint T F[5]等將雙橢球分布推廣到雙橢球錐形熱功率密度模型中，準(zhǔn)確地描述了更廣泛的幾何形狀和焊接過程的瞬態(tài)溫度場；李紅濤[6]等以Q345 中厚板為研究對象，將雙橢球和均勻體兩種熱源分布方式作為輸入，研究焊接其對溫度場的影響。研究結(jié)果表明：當(dāng)焊接初始條件相同時，二者瞬態(tài)溫度場分析結(jié)果相似，后者的結(jié)果更為理想。但均勻體熱源加載時焊縫起始位置溫度較低。因而提出在焊接仿真分析中焊接的起始位置采用雙橢球熱源，其它位置采用均勻體熱源的加熱方式。

本文選取規(guī)格分別為100.0 mm×40.0 mm×3.2 mm 與100.0 mm×60.0 mm×3.2 mm 的Q345板進(jìn)行焊接，在ANSYS 軟件中建立對應(yīng)的簡化模型，分別以半球狀、旋轉(zhuǎn)高斯體、雙橢球3 種熱源形式作為輸入，分析輸出結(jié)果以研究不同形式的體熱源對焊接過程中的瞬時溫度場影響。

1 熱源模型及參數(shù)

1.1 半球狀熱源模型

在電弧引起的沖擊較小、熔池的深度較淺時，平面分布的熱源模型還能較為準(zhǔn)確地反應(yīng)焊接狀態(tài)，但此時沒有考慮到焊接時熔深的作用，為了將焊接時熔深方面的影響考慮在內(nèi)，發(fā)展出更加有效的均勻體熱源模型[7]。由于均勻體熱源模型過于簡化，有學(xué)者提出更為符合實際情況的半球狀熱源分布函數(shù)，其表達(dá)式為：

式中：q（x，y，z）——點(diǎn)（x，y，z）處熱流密度；Q——熱輸入率；c——半球體的半徑。

1.2 旋轉(zhuǎn)高斯體熱源模型

在分析過程中，平面高斯熱源模型是一種典型的平面熱源分布函數(shù)。此模型是在平面內(nèi)將熱源能量按照高斯函數(shù)的分布，在給定半徑的圓內(nèi)施加熱量，其具體公式為：

式中：q（r）——與圓心中相距r 處的熱流密度；σq——高斯熱源分布參數(shù)；Q——熱輸入率；r——與圓心距離。

但是在計算過程中，平面高斯熱源認(rèn)為焊接時的熱源分布密度具有對稱性的特點(diǎn)。當(dāng)焊接過程處在低速階段時，可以看作是近似合理的；但當(dāng)焊接速度較大時，熱源的分布并不對稱，平面高斯熱源不能準(zhǔn)確表達(dá)熱源的分布過程。為準(zhǔn)確地模擬出熔池形狀，吳甦[8]提出了旋轉(zhuǎn)高斯體熱源模型，該熱源模型為Gauss 曲線繞對稱軸旋轉(zhuǎn)生成的曲面圍繞而成的曲面體，如式（3）、式（4）所示，其形式如圖1 所示。

圖1 旋轉(zhuǎn)體高斯熱源模型Fig.1 Gaussian heat source model of revolution

式中：H——熱源高度；Q——熱輸入率；R0——熱源半徑。

1.3 雙橢球熱源模型

John Goldak 提出的雙橢球熱源分布函數(shù)綜合考慮了焊接時前端溫度梯度較大而后端溫度梯度較小的特點(diǎn)，沿焊接方向，將焊接熱輸入分為前后兩個橢球來處理。前后兩部分橢球的熱流密度分布函數(shù)分別如式（5）、式（6）所示，其形式如圖2 所示。

圖2 雙橢球熱源模型Fig.2 Double ellipsoid heat source model

沿x 軸方向雙橢球正方向熱源分布函數(shù)為

沿x 軸方向雙橢球負(fù)方向熱源分布函數(shù)為

式中：qf（x，y，z）——前半球熱流密度；qr（x，y，z）——后半球熱流密度；a1，a2，b，c——前后橢球的形狀參數(shù)；Q——焊接熱輸入；ff，fr——前后部分所占總輸入量的比例。

2 有限元模型的建立及分析

2.1 有限元模型

試驗采用的是惰性氣體保護(hù)焊方法，最后完成的焊接試樣如圖3 所示。相對應(yīng)，在ANSYS軟件中建立該T 型樣件的簡化模型，最后建立的簡化模型如圖4 所示。該模型有196 539 個節(jié)點(diǎn)，40 400 個單元。

圖3 T 型焊接試樣Fig.3 T-type welding specimen

圖4 ANSYS 構(gòu)建的簡化模型Fig.4 Simplified model built in ANSYS

2.2 材料參數(shù)

材料在不同溫度下物理性能存在差異，并且材料的熱物理參數(shù)對模擬結(jié)果會產(chǎn)生較大影響。本文采用的材料為Q345 結(jié)構(gòu)鋼，由文獻(xiàn)[9]知其性能隨溫度變化規(guī)律如表1 所示，缺失溫度下的值則通過線性插值法獲得。

表1 Q345 結(jié)構(gòu)鋼的物理性能Tab.1 Physical properties of Q345 structural steel

（續(xù)表）

2.3 熱源參數(shù)

本文在進(jìn)行實際工件焊接時采用惰性氣體保護(hù)焊，工作參數(shù)：焊接電流150 A，電壓20 V，焊接速度5 mm/s，焊接效率為0.83。將焊接時的實際參數(shù)作為焊接熱源參數(shù)的選擇標(biāo)準(zhǔn)，得到的不同熱源參數(shù)如表2 所示。

表2 不同體熱源參數(shù)選擇Tab.2 Parameter selection of different body heat sources

3 3 種熱源分析結(jié)果對比

3.1 焊接溫度場分析結(jié)果對比

將3 種不同的熱源模型作為熱載荷輸入，用ANSYS 進(jìn)行分析，結(jié)果如圖5—圖7 所示。圖5為不同時刻半球狀熱源模型作用下的溫度場結(jié)果。圖6 與圖7 分別為對應(yīng)時刻下旋轉(zhuǎn)高斯體熱源模型與雙橢球熱源模型作用下的溫度場結(jié)果。由溫度場分析結(jié)果可知，焊接熱源中心沿著焊縫不斷移動，距焊接熱源中心處溫度最高，熱量梯度呈橢球狀分布。隨著焊接的結(jié)束，由于對流換熱，溫度逐漸降低，焊接完成100 s 后，最高溫度小于50 ℃；焊接完成2 000 s 后，恢復(fù)至室溫。

圖5 半球狀熱源模型溫度場結(jié)果Fig.5 Temperature field results of hemispherical heat source model

圖6 旋轉(zhuǎn)高斯體熱源模型溫度場結(jié)果Fig.6 Temperature field results of rotating Gaussian body heat source model

圖7 雙橢圓熱源模型溫度場結(jié)果Fig.7 Temperature field results of double elliptical heat source model

通過比較相同時間下3 種熱源模型的瞬時溫度場，3 種體熱源模型的溫度場分布趨勢基本一致，熔池的形狀也基本相似。雖然半球狀熱源模型與旋轉(zhuǎn)高斯體熱源模型分別假設(shè)焊接熱量是按照半圓體和高斯體分布，但是由于存在熱傳導(dǎo)以及熱源移動的作用，從表面看，焊接中心附近的溫度場分布也近似雙橢圓狀。從溫度場結(jié)果看，3 種體熱源的分布狀況基本相似，因此要確切比對3 種熱源模型的瞬時溫度場結(jié)果差異還需要提取溫度測試點(diǎn)進(jìn)行更加詳細(xì)的比對。

3.2 焊趾處溫度結(jié)果對比

在焊趾部位取等距離的3 個點(diǎn)，提取3 個點(diǎn)的溫度數(shù)據(jù)，提取部位如圖8 所示，半球狀熱源模型、旋轉(zhuǎn)高斯體熱源模型與雙橢球熱源模型分析結(jié)果分別如圖9—圖11 所示。由于在100 s 到2 000 s 之間溫度由50 ℃變化到20 ℃，溫度變化緩慢，因而僅提取前100 s 數(shù)據(jù)進(jìn)行分析。

圖8 有限元模型中溫度點(diǎn)取值Fig.8 Temperature point value in finite element model

圖9 半球狀熱源模型中三點(diǎn)溫度場變化圖Fig.9 Temperature field of three points in hemispherical heat source model

圖10 旋轉(zhuǎn)高斯體熱源模型中三點(diǎn)溫度場變化圖Fig.10 Temperature field of three points in rotating Gaussian heat source model

圖11 雙橢球熱源模型中三點(diǎn)溫度場變化圖Fig.11 Temperature field of three points in double ellipsoid heat source model

由以上溫度變化曲線可知，無論是哪一種焊接熱源模型，焊趾處三點(diǎn)都是在熱源中心移動到附近處時溫度達(dá)到最大值；峰值過后，由于對流換熱，溫度迅速降低。

由圖9—圖11 對比可知，在焊接過程中，采用旋轉(zhuǎn)高斯體熱源時，溫度最高，其最高溫度為A 點(diǎn)溫度227.93 ℃，B 點(diǎn)溫度627.77 ℃，C 點(diǎn)溫度550.70℃；其次是半球狀熱源，最高溫度為A 點(diǎn)溫度246.42 ℃，B 點(diǎn)溫度573.06 ℃，C 點(diǎn)溫度515.52 ℃；雙橢球熱源溫度最低，A 點(diǎn)最高溫度為221.03 ℃；B 點(diǎn)最高溫度為457.71 ℃；C 點(diǎn)最高溫度為399.99 ℃。

3.3 熔池形貌對比分析

為比較3 種體熱源在T 型焊接工件的惰性氣體保護(hù)焊有限元分析中的真實性，對焊接工件進(jìn)行熔深測量，熔池截面形貌如圖12 所示。

圖12 熔池截面形貌Fig.12 Cross-section morphology of molten pool

分別提取3 種熱源模型的熔池與實際焊接工件的熔池做對比分析，3 種不同熱源模型的熔池形貌如圖13 所示。

圖13 3 種不同熱源模型的熔池形貌Fig.13 Pool morphology of three different heat source models

由圖13 可知，由于半球體模型與旋轉(zhuǎn)體高斯熱源模型溫度較高，其熔池深度較高，而雙橢球熱源模型能有效反應(yīng)焊接過程中熔池形態(tài)。因而，采用雙橢球熱源模型進(jìn)行焊接模擬仿真具有一定的準(zhǔn)確性。

4 結(jié)論

本文以T 型焊接樣件為研究對象，在ANSYS中使用半球狀熱源、旋轉(zhuǎn)高斯體熱源、雙橢球熱源分別對T 型焊接工件有限元模型進(jìn)行分析，得到以下結(jié)論：

（1）在焊接條件一致的條件下，3 種體熱源模型所得到的瞬態(tài)溫度場趨勢一致，熔池形狀均為橢球狀，但3 種熱源模型所得到的溫度場存在差異，取焊趾處溫度進(jìn)行分析，旋轉(zhuǎn)高斯體熱源產(chǎn)生的溫度最高，其次是半球狀熱源，雙橢球熱源溫度最低。

（2）經(jīng)過熔池形貌的對比發(fā)現(xiàn)，相對于旋轉(zhuǎn)高斯體熱源和半球狀熱源，雙橢球熱源模型能有效模擬焊接過程中熔池形態(tài)。因而，在本次焊接中采用雙橢球熱源模型進(jìn)行焊接模擬仿真具有一定的準(zhǔn)確性。