賈劍平，徐鑫，李田雨

南昌大學(xué)先進(jìn)制造學(xué)院江西省機(jī)器人與焊接自動(dòng)化重點(diǎn)實(shí)驗(yàn)室 江西南昌 330031

1 序言

鋁合金是一種重要的輕量化材料，具有比強(qiáng)度高的特點(diǎn)，廣泛應(yīng)用于軌道交通、汽車和航空航天等領(lǐng)域[1，2]。由于鋁合金導(dǎo)熱系數(shù)和熱膨脹系數(shù)較大，若板材厚度較小，則在焊接高溫作用下容易產(chǎn)生嚴(yán)重的殘余應(yīng)力和焊接變形，不僅影響結(jié)構(gòu)美觀，也會(huì)降低其承載強(qiáng)度，無(wú)法滿足使用要求[3-5]。冷金屬過(guò)渡（CMT）技術(shù)通過(guò)換向送絲系統(tǒng)和數(shù)字式焊接控制系統(tǒng)，能夠降低焊接熱輸入量，為鋁合金薄板焊接提供了新的選擇[6-8]。

近年來(lái)，焊接數(shù)值模擬技術(shù)發(fā)展迅速，通過(guò)數(shù)值模擬能夠顯著減少試驗(yàn)量，降低生產(chǎn)成本。在CMT焊接數(shù)值模擬方面，AZAR等[9]首先提出了形狀參數(shù)隨時(shí)間正弦變化的新型雙橢球熱源模型，通過(guò)熱源形狀的周期性變化反映CMT焊接過(guò)程中電弧的瞬時(shí)變化；ZAPICO等[10]對(duì)AZAR提出的熱源模型進(jìn)行了簡(jiǎn)化，發(fā)現(xiàn)鋁合金CMT焊接過(guò)程中峰值溫度可達(dá)到1800℃，有利于去除氧化物；而整個(gè)焊接過(guò)程中平均溫度僅為約1300℃，有效控制了焊接熱輸入和熱影響區(qū)的擴(kuò)展；王遠(yuǎn)傳等[11]模擬了高強(qiáng)鋼的CMT焊接溫度場(chǎng)，分析了溫度對(duì)熱影響區(qū)寬度的影響；楊羅揚(yáng)等[12]通過(guò)模擬不銹鋼CMT增材制造過(guò)程，優(yōu)化了層間冷卻時(shí)間和增材成形路徑。張奇奇等[13]基于動(dòng)網(wǎng)格技術(shù)成功模擬了CMT焊過(guò)程中焊絲的送進(jìn)與回抽運(yùn)動(dòng)，最終模擬的電弧形態(tài)與實(shí)際拍攝的電弧形態(tài)一致。趙文勇等[14]對(duì)CMT電弧增材制造的熱-流場(chǎng)進(jìn)行研究，通過(guò)橢球熱源的周期性加載分析了熔池流動(dòng)行為、熔覆層形貌和熔池尺寸的影響因素。

本文采用ANSYS軟件建立三維仿真模型，通過(guò)焊接熱源的周期性加載，對(duì)5052鋁合金薄板CMT焊接溫度場(chǎng)進(jìn)行分析，根據(jù)數(shù)值模擬結(jié)果研究CMT焊接溫度場(chǎng)特點(diǎn)，并通過(guò)試驗(yàn)驗(yàn)證了溫度場(chǎng)模擬結(jié)果的準(zhǔn)確性。

2 有限元分析

2.1 材料屬性

焊接接頭形式為對(duì)接，焊接母材為5052鋁合金，尺寸為200m m×50m m×0.8m m，焊絲為ER5183鋁鎂焊絲，直徑為1.2mm。使用JMatPro軟件計(jì)算出5052鋁合金的高溫?zé)嵝阅?，如圖1所示。由于焊絲與母材成分相近，因此假設(shè)焊縫金屬熱物理性能與母材相同。同時(shí)，由于鋁合金電弧焊接中存在固態(tài)金屬熔化過(guò)程，所以在分析焊接溫度場(chǎng)時(shí)必須考慮熔化潛熱現(xiàn)象。根據(jù)相關(guān)文獻(xiàn)[15]，采用等效比熱容法處理鋁合金熔化潛熱。

圖1 JMatPro計(jì)算的材料熱性能參數(shù)

式中ce——等效比熱容[J/（g·℃）]；

c——材料本身的比熱容[J/（g·℃）]；

c0——潛熱比熱容[J/（g·℃）]；

L——材料的熔化潛熱（J/mol）；

TL——材料液相線溫度（℃）；

TS——材料固相線溫度（℃）。

2.2 模型建立

采用ANSYS 15.0有限元分析軟件建立幾何模型。為了減少單元與節(jié)點(diǎn)數(shù)量，提升計(jì)算效率，只建立單側(cè)的模型進(jìn)行計(jì)算，幾何模型如圖2所示。由圖2可知，幾何模型劃分為焊縫區(qū)、過(guò)渡區(qū)和遠(yuǎn)離焊縫區(qū)，模型整體尺寸為200mm×50mm×0.8mm。

圖2 幾何模型

為提升計(jì)算準(zhǔn)確度與計(jì)算效率，分別對(duì)焊縫區(qū)、過(guò)渡區(qū)和遠(yuǎn)離焊縫區(qū)進(jìn)行網(wǎng)格劃分。選用SOLID 70單元進(jìn)行溫度場(chǎng)模擬分析，模型中焊縫區(qū)單元尺寸設(shè)置為0.5mm×0.2mm×0.26mm，過(guò)渡區(qū)單元尺寸為0.5mm×0.8mm×0.26mm，遠(yuǎn)離焊縫區(qū)單元尺寸為0.5mm×5mm×0.26mm。最終生成的網(wǎng)格模型包含28800個(gè)單元、40100個(gè)節(jié)點(diǎn)，薄板網(wǎng)格劃分模型如圖3所示。

圖3 薄板網(wǎng)格劃分模型

2.3 初始條件與邊界條件

初始溫度設(shè)置為環(huán)境溫度30℃，在焊接過(guò)程中，鋁板表面存在熱輻射和熱對(duì)流[16]。根據(jù)式（3）[17]對(duì)工件表面換熱系數(shù)進(jìn)行處理，將處理后的等效換熱系數(shù)he通過(guò)熱流密度的形式施加到模型與空氣接觸的4個(gè)自由表面。由于焊接時(shí)采用金屬底板作為支撐，而金屬導(dǎo)熱系數(shù)較大，散熱較快，因此底面的換熱系數(shù)設(shè)定為100W/（m2·℃）。模型中心對(duì)稱面處設(shè)置絕熱邊界條件，即設(shè)置熱流密度為零。

式中he——等效換熱系數(shù)[W/（m2·℃）]；

T——表面溫度（℃）。

2.4 焊接熱源

由于焊接電弧的熱流分布具有一定對(duì)稱性，且本研究中焊接速度不大、母材板厚較小，所以本文使用高斯平面熱源（見(jiàn)圖4）作為熱源模型進(jìn)行模擬計(jì)算[18-20]。高斯熱源表達(dá)式為

圖4 高斯平面熱源模型

式中q(r)——距離熱源中心r處的熱流密度[J/（s·m2）]；

r——距熱源中心距離（mm）；

Qm——熱源加熱中心位置的最大熱流密度[J/（s·m2）]；

R——熱源的有效加熱半徑（mm），本文取有效加熱半徑R=2.5mm。

CMT技術(shù)的焊絲回抽頻率可達(dá)70Hz，因此溫度場(chǎng)計(jì)算的時(shí)間步長(zhǎng)設(shè)置為一個(gè)定值t，取t=0.015s。為了真實(shí)反映出CMT焊接過(guò)程熱輸入的周期性變化，通過(guò)ANSYS參數(shù)化設(shè)計(jì)語(yǔ)言（APDL）編寫(xiě)命令流，實(shí)現(xiàn)熱源函數(shù)的周期性加載，如圖5所示。圖5中將每個(gè)時(shí)間步t細(xì)分為一個(gè)加熱子步tH和一個(gè)冷卻子步tC，在加熱子步tH內(nèi)，對(duì)焊縫區(qū)域加載高斯熱源；在冷卻子步tC內(nèi)，設(shè)置焊縫區(qū)域熱流密度為零，以此反映CMT技術(shù)短路過(guò)渡時(shí)熱輸入的降低，每個(gè)冷卻子步占一個(gè)時(shí)間步長(zhǎng)的30%。為分析CMT焊接溫度場(chǎng)特點(diǎn)，同時(shí)對(duì)傳統(tǒng)MIG焊進(jìn)行模擬計(jì)算，計(jì)算中采用的三維模型、網(wǎng)格劃分、熱源參數(shù)和熱源效率均與CMT焊保持一致，但熱源加載方式為持續(xù)加載，不存在冷卻子步。

圖5 CMT熱輸入示意

3 溫度場(chǎng)結(jié)果分析

3.1 CMT焊接溫度場(chǎng)特點(diǎn)

模擬采用的焊接參數(shù)見(jiàn)表1。CMT焊接過(guò)程中三個(gè)連續(xù)時(shí)刻的溫度場(chǎng)云圖如圖6所示。圖6a中，電弧正處于燃燒狀態(tài)，在電弧熱輸入作用下，薄板溫度升高，熔池中心最高溫度達(dá)到1043℃；圖6b中，熔滴與熔池接觸形成短路，電弧熄滅，同時(shí)焊機(jī)自動(dòng)將焊接電流降低至極小值，此時(shí)焊接熱輸入幾乎為零，熔池中心最高溫度降低至896℃，并且薄板整體溫度均呈現(xiàn)降低趨勢(shì)；圖6c中，焊絲回抽，熔滴順利向熔池過(guò)渡，焊接電弧重新被引燃，繼續(xù)向薄板輸送熱量。此時(shí)，薄板的整體溫度上升，溫度場(chǎng)迅速恢復(fù)到降溫之前的狀態(tài)。由此可看出，在CMT周期性變化的熱輸入下，薄板溫度場(chǎng)也呈現(xiàn)有規(guī)律的變化，整個(gè)焊接過(guò)程由若干個(gè)這樣的周期組成。

圖6 CMT焊接溫度場(chǎng)云圖

表1 模擬采用的焊接參數(shù)

與CMT焊相比，傳統(tǒng)MIG焊工藝中不存在焊絲回抽運(yùn)動(dòng)和熱輸入降低的過(guò)程，因此整個(gè)焊接過(guò)程中的熱輸入相對(duì)穩(wěn)定。為了與CMT數(shù)值模擬結(jié)果進(jìn)行比較，選取MIG焊過(guò)程中三個(gè)連續(xù)時(shí)刻的溫度場(chǎng)分布云圖進(jìn)行分析，如圖7所示。由圖7可看出，MIG焊過(guò)程中溫度場(chǎng)分布較為穩(wěn)定，熔池中心的最高溫度保持在1480℃左右。

圖7 MIG焊接溫度場(chǎng)云圖

如圖8所示，在焊縫中段沿焊接方向等距選取11個(gè)節(jié)點(diǎn)（n1～n11），通過(guò)*GET命令提取各節(jié)點(diǎn)在整個(gè)焊接過(guò)程中的峰值溫度進(jìn)行分析。圖9所示為不同送絲速度下CMT焊和MIG焊過(guò)程中各節(jié)點(diǎn)峰值溫度及其平均值。保持焊接速度不變，當(dāng)送絲速度為2.3m/min時(shí)，CMT焊各節(jié)點(diǎn)的平均峰值溫度為881.4℃，MIG焊各節(jié)點(diǎn)平均峰值溫度為1266℃；當(dāng)送絲速度為2.6m/min時(shí)，CMT焊各節(jié)點(diǎn)的平均峰值溫度為1039.5℃，MIG焊各節(jié)點(diǎn)平均峰值溫度為1480.9℃；當(dāng)送絲速度為2.9m/min時(shí)，CMT焊各節(jié)點(diǎn)的平均峰值溫度為1206.5℃，MIG焊各節(jié)點(diǎn)平均峰值溫度為1717.4℃；當(dāng)送絲速度為3.2m/min時(shí)，CMT焊各節(jié)點(diǎn)的平均峰值溫度為1376.1℃，MIG焊各節(jié)點(diǎn)平均峰值溫度為1966.4℃。由此可看出，在焊接速度保持恒定時(shí)，隨著送絲速度增大，CMT焊和MIG焊模擬結(jié)果中各節(jié)點(diǎn)的峰值溫度均有明顯上升。在各組參數(shù)下，CMT焊模擬結(jié)果中各節(jié)點(diǎn)的峰值溫度與平均溫度均遠(yuǎn)低于MIG焊模擬溫度，體現(xiàn)出CMT技術(shù)熱輸入較低的特點(diǎn)。

圖8 節(jié)點(diǎn)選取示意

圖9 不同送絲速度下節(jié)點(diǎn)的峰值溫度與平均值

通過(guò)比較溫度分布云圖與焊縫中段各節(jié)點(diǎn)的峰值溫度可得出，使用CMT工藝進(jìn)行薄板焊接能夠顯著減小焊接熱輸入量，薄板整體溫度更低。根據(jù)節(jié)點(diǎn)平均溫度峰值的變化發(fā)現(xiàn)：隨著送絲速度增大，兩種工藝模擬的節(jié)點(diǎn)溫度峰值之間的差值逐漸增大，CMT工藝的低熱輸入優(yōu)勢(shì)更加明顯。

3.2 溫度場(chǎng)試驗(yàn)驗(yàn)證

如圖10所示，為了驗(yàn)證溫度場(chǎng)模擬結(jié)果的準(zhǔn)確性，采用K形熱電偶測(cè)量鋁合金薄板表面P1、P2和P3三點(diǎn)（距離焊縫中心分別為15mm、20mm和25mm）的實(shí)時(shí)溫度并記錄，繪制三點(diǎn)的熱循環(huán)曲線。在有限元模型中相應(yīng)位置選取節(jié)點(diǎn)，在后處理模塊中將節(jié)點(diǎn)溫度定義為時(shí)間的變量，提取各節(jié)點(diǎn)的模擬熱循環(huán)曲線，與實(shí)測(cè)熱循環(huán)曲線進(jìn)行比較。選用A2參數(shù)進(jìn)行測(cè)溫試驗(yàn)。

圖10 測(cè)溫點(diǎn)示意

P1、P2和P3三點(diǎn)的實(shí)測(cè)熱循環(huán)曲線和模擬熱循環(huán)曲線如圖11所示。在實(shí)際焊接過(guò)程中，隨著焊接電弧的移動(dòng)，熱量從焊縫中心部位向周圍金屬傳遞，測(cè)量點(diǎn)的溫度不斷上升并迅速達(dá)到溫度峰值；電弧離開(kāi)后，測(cè)量點(diǎn)的溫度逐漸下降，最終恢復(fù)到室溫值。試驗(yàn)中，P1點(diǎn)距離焊縫最近，峰值溫度最高；P3點(diǎn)距離焊縫較遠(yuǎn)，峰值溫度較低。由此可看出，薄板上表面距離焊縫中心不同位置處節(jié)點(diǎn)的模擬熱循環(huán)曲線與試驗(yàn)測(cè)量的熱循環(huán)曲線變化趨勢(shì)基本一致。分析發(fā)現(xiàn)，模擬溫度值略高于實(shí)際溫度，可能是因材料參數(shù)和散熱邊界條件設(shè)置而與實(shí)際存在一定誤差。

圖11 熱循環(huán)曲線對(duì)比

在垂直焊縫方向切取試樣，對(duì)金相試樣進(jìn)行打磨與拋光后，使用凱勒試劑對(duì)焊縫試樣進(jìn)行腐蝕，拍攝焊縫截面圖，并將其與模擬的熔池截面圖進(jìn)行比較，如圖12所示。由圖12可知，模擬結(jié)果與實(shí)際焊縫輪廓比較接近。

圖12 焊縫截面示意

通過(guò)熱電偶測(cè)溫試驗(yàn)與焊縫截面輪廓對(duì)比，可認(rèn)為溫度場(chǎng)模擬結(jié)果與實(shí)際焊接溫度場(chǎng)分布情況接近，采用該模型進(jìn)行鋁合金薄板CMT焊接有限元分析具有可行性，能夠用于后續(xù)的分析研究。

4 結(jié)束語(yǔ)

本文使用ANSYS軟件建立鋁合金焊接溫度場(chǎng)三維分析模型，考慮材料熱物理參數(shù)隨溫度的非線性變化，通過(guò)命令流實(shí)現(xiàn)了焊接熱源的周期性加載，得出以下結(jié)論。

1）溫度場(chǎng)模擬結(jié)果與試驗(yàn)結(jié)果較為吻合，表明該模型能夠反映CMT焊接溫度場(chǎng)特點(diǎn)，可用于CMT焊接數(shù)值模擬研究。

2）CMT焊接過(guò)程中存在多個(gè)冷卻階段，能夠顯著降低焊接熱輸入，對(duì)于鋁合金薄板焊接十分有利。

3）與傳統(tǒng)MIG焊工藝相比，CMT焊接薄板整體溫度較低，且送絲速度越大，CMT技術(shù)的冷卻作用越顯著。