程虎 沈精虎 戴昆 劉蓉

摘要：? 針對地鐵枕梁焊接中存在的殘余應(yīng)力對列車性能和質(zhì)量產(chǎn)生的影響問題，本文基于Simufact. welding軟件，對B型鋁合金地鐵車體底架中的枕梁進(jìn)行焊接數(shù)值模擬分析，對焊接后的殘余應(yīng)力與變形分布進(jìn)行分析，并對不同的焊接接頭形式創(chuàng)建不同的熱源模型，對實際工況下的多層焊道進(jìn)行簡化處理，以單層焊道代替多層焊道，得到焊接的熱循環(huán)曲線。同時，以簡化處理得到的熱源模型代替常規(guī)的雙橢球熱源模型，對整個枕梁進(jìn)行數(shù)值模擬分析。仿真結(jié)果表明，由于焊縫附近位置節(jié)點(diǎn)被約束，當(dāng)焊縫焊接完成，進(jìn)入冷卻過程時，焊縫收縮。因遠(yuǎn)離焊縫的枕梁型材中部和弧形板處沒有約束，因此變形較大。該研究比較準(zhǔn)確地模擬了實際生產(chǎn)中焊接應(yīng)力的變化情況，為地鐵枕梁的焊接提供了技術(shù)支持。

關(guān)鍵詞：? Simufact.welding; 枕梁; 熱源模型; 熱循環(huán)曲線; 模擬分析

中圖分類號： TG404文獻(xiàn)標(biāo)識碼： A

作者簡介：? 程虎（1989-），男，碩士研究生，主要研究方向為焊接仿真分析。

通信作者：? 沈精虎（1965-），男，碩士生導(dǎo)師，主要研究方向為機(jī)械CAD/CAM技術(shù)。 Email： postmaster@laohu.net

地鐵屬于城市軌道交通的一部分，但與城市道路基本沒有交叉，其具有速度快、運(yùn)能大和時間準(zhǔn)[1]的優(yōu)點(diǎn)。近年來，隨著我國城市軌道交通技術(shù)的不斷發(fā)展，列車速度不斷提升，其安全性、可靠性、舒適性等問題得到了廣泛關(guān)注[23]。為提高列車的行進(jìn)速度，增大載客能力，提高車輛行駛的安全性，降低制造和運(yùn)營成本，地鐵輕量化是一種行之有效的措施[4]。目前，我國地鐵普遍采用鋁合金車體，鋁的密度約為碳鋼的1/3，與不銹鋼車體相比，鋁合金車體的自重一般可減少30%～40%[56]。國內(nèi)的鋁合金系主要有1000系純鋁系、2000系鋁銅系、3000系鋁錳系、4000系鋁硅系、5000系鋁鎂系、6000系鋁硅鎂系、7000系鋁鋅鎂系和以其他元素為主的8000系[7]。其中，6000系鋁硅鎂合金是國內(nèi)鋁合金牌號中的鍛鋁合金，具有良好的抗腐蝕性能和焊接性能，可熱處理強(qiáng)化，無應(yīng)力腐蝕傾向，具有極佳的壓力加工成型性[8]。由于有限元仿真技術(shù)已廣泛應(yīng)用于焊接問題的研究中，為降低動車組車體加工中的焊接殘余應(yīng)力，馬思群等人[9]對動車組KK端多道焊焊接殘余應(yīng)力進(jìn)行仿真研究;而為了提高鋼板焊接質(zhì)量，王士軍等人[10]采用Ansys軟件對鋼板填充焊接過程進(jìn)行了有限元模擬仿真分析;在不確定響應(yīng)下，Zheng Li等人[11]采用焊接模擬方法對結(jié)構(gòu)進(jìn)行了熱力耦合分析。鋁合金車體主要由底架、側(cè)墻、端墻和車頂?shù)却蟛考M成，其中車體枕梁位于底架中，是鋁合金車體重要的支撐和連接部件[12]，此部位與牽引梁等部件相連，受力較大且復(fù)雜[13]。枕梁中，各組件的連接方式為焊接，焊接方式均采用熔化極惰性氣體保護(hù)焊（metal inertgas welding，MIG）[14]，其焊縫較多，殘余應(yīng)力較大，影響車體的裝備精度和外形，降低車體的承載能力和使用壽命，易出現(xiàn)安全隱患，需要進(jìn)行校核，以確保列車行車安全[15]?；诖?，本文主要對B型鋁合金地鐵枕梁進(jìn)行焊接數(shù)值模擬分析，仿真較為準(zhǔn)確地模擬了實際生產(chǎn)中焊接應(yīng)力的變化情況。該研究對提高地鐵枕梁的焊接質(zhì)量具有重要意義。

1焊接過程數(shù)值模擬分析

在計算鋁合金車體枕梁結(jié)構(gòu)應(yīng)力應(yīng)變時，采用SIMUFACT公司旗下的Simufact. welding軟件進(jìn)行焊接工藝仿真計算，并對焊接過程進(jìn)行數(shù)值模擬分析。

1.1幾何模型和網(wǎng)格劃分

枕梁結(jié)構(gòu)幾何模型寬1 138 mm，長2 580 mm。共有焊縫17條，最短處焊縫238 mm，最長焊縫2 580 mm，焊接類型為MIG焊，枕梁結(jié)構(gòu)幾何模型如圖1所示。

枕梁結(jié)構(gòu)的網(wǎng)格劃分采用六面體結(jié)構(gòu)單元，在Hypermesh軟件中進(jìn)行網(wǎng)格劃分，網(wǎng)格單元103 760個，節(jié)點(diǎn)數(shù)139 801。焊縫處的網(wǎng)格尺寸約為5 mm，熱影響區(qū)的網(wǎng)格尺寸為10～15 mm，其母材網(wǎng)格尺寸為15～20 mm，網(wǎng)格尺寸呈漸變增加趨勢[16]。枕梁結(jié)構(gòu)劃分后的網(wǎng)格模型如圖2所示。

1.2材料數(shù)據(jù)庫的建立

在鋁合金車體中，枕梁材料為6005A-T6，焊槍填充材料為AlMg4.5MnZr，材料化學(xué)成分如表1所示。

為準(zhǔn)確模擬實際工況下的焊接過程，需要在Simufact. welding中建立6005A-T6和AlMg4.5MnZr的數(shù)據(jù)庫，使用JMatPro軟件可以較為準(zhǔn)確地查找出材料在不同溫度下的熱力學(xué)性能，材料比熱隨溫度變化曲線如圖3所示，材料熱導(dǎo)率隨溫度變化曲線如圖4所示，材料密度隨溫度變化曲線如圖5所示。材料的這三種熱力學(xué)性能是影響焊接分析結(jié)果最主要的三個參數(shù)，在Simufact. welding中需要對這三項參數(shù)重新進(jìn)行設(shè)定。

1.3熱源模型

熱源模型是作用在焊接件上，在時間域和空間域上熱輸入分布特點(diǎn)的一種數(shù)學(xué)表達(dá)式[17]。目前，用于焊接數(shù)值分析中所有的焊接熱源模型大都不隨時間的變化而變化，也就是靜態(tài)熱源[18]。Simufact. welding支持常規(guī)型和激光焊兩個預(yù)定義的熱源，常規(guī)型用于氣焊、手工電弧焊、MSG、WSG和UP焊接;激光焊用于激光和電子束焊接（正太分布區(qū)域熱源和恒定圓柱熱源組合）。由于枕梁處的焊接是MIG焊，屬于常規(guī)焊接的一種，采用雙橢球熱源模型，雙橢球熱源模型如圖6所示。

1.4焊接順序

按照實際加工順序設(shè)置焊接順序，枕梁焊接順序如圖7所示。圖中帶S的焊縫為相同焊縫的對稱位置，按照如下順序進(jìn)行焊接：

1.5焊接接頭

在實際工藝中，鋁合金地鐵枕梁的焊接形式有兩種：一種是兩個10 mm板對接，另一種是兩個25 mm板對接。對于薄板（10 mm）焊接，實際焊接工藝采用3層焊道，對于厚板（25 mm）焊接，采用10～14層焊道。對不同形式的焊接接頭形式創(chuàng)建不同的熱源模型。

兩個10 mm板對接接頭，焊接效率取0.8～0.9，焊接接頭形式如圖8所示，圖8中，t1=10 mm，t2=10 mm，b≤1 mm，c=2～3 mm，α=70°。焊接順序如圖9所示，其焊接工藝規(guī)程如表2所示。

兩個25 mm板對接接頭，焊接效率取0.8～0.9，焊接接頭形式如圖10所示，圖10中，t1=25 mm，t2=25 mm，b≤1 mm，c=2～3 mm，α=70°。焊接順序如圖11所示，其焊接工藝規(guī)程如表3所示。

對于大型焊接仿真，熱源模型一般不會選擇雙橢球熱源，而是選擇簡化的熱循環(huán)曲線作為熱源模型[21]。如果焊縫焊道過多，熱循環(huán)曲線的獲取會變得困難，即使得到了熱循環(huán)曲線，大型結(jié)構(gòu)的仿真計算也變得耗時巨大，實際工程中的經(jīng)濟(jì)效益會變低。因此，可進(jìn)一步簡化多層焊道的焊縫，用單層焊道代替多層焊道，從而得出簡化的熱循環(huán)曲線。將兩個10 mm板對接接頭分別做單層焊、雙層焊和多層焊處理，得到焊接后的殘余應(yīng)力。同理，將兩個25 mm板對接接頭分別做單層焊、雙層焊和多層焊處理。焊接殘余應(yīng)力隨距焊縫中心位置變化曲線如圖12所示。

由圖12可以看出，采用不同層數(shù)的焊接，殘余應(yīng)力的大小與分布趨勢相似，因此可以把多層焊焊接用單層焊焊接簡化處理，得到單層焊接的熱循環(huán)曲線。焊接溫度隨時間變化曲線如圖13所示。由圖13可以看出，該熱循環(huán)曲線分別提取了2個10 mm板對接接頭和2個20 mm板對接接頭焊縫沿著焊接方向中間截面的節(jié)點(diǎn)從開始的初始溫度20 ℃到焊縫熔化并升溫到最高溫度隨后冷卻到200 ℃以下的溫度變化。

2結(jié)果與分析

采用簡化后的熱循環(huán)曲線作為熱源進(jìn)行焊接仿真。在鋁合金車體枕梁的整體計算中，枕梁母材材料為6005A-T6，焊縫填充材料為AlMg4.5MnZr，熔化溫度范圍為580～620 ℃。按照圖7梁焊接順序進(jìn)行焊接，并加載熱循環(huán)曲線，當(dāng)焊縫溫度升到最高時，溫度場最高為677.96 ℃。枕梁焊接完成后，位移最大為2.68 mm，發(fā)生在枕梁型材中間位置，枕梁整體位移變形如圖14所示。

由于焊縫附近位置節(jié)點(diǎn)被約束，當(dāng)焊縫焊接完成進(jìn)入冷卻過程時，焊縫收縮，而遠(yuǎn)離焊縫的枕梁型材中部和弧形板處沒有約束，因此變形較大。采用熱循環(huán)曲線作為熱源模型，加載完畢后，提取枕梁有限元模型中應(yīng)力最大焊縫和應(yīng)力最小焊縫縱向焊接殘余應(yīng)力，殘余應(yīng)力隨距焊縫中心位置變化曲線如圖15所示。由圖15可以看出，焊縫處的殘余應(yīng)力的軌跡是相似的，且符合一般焊接過程的分布規(guī)律[17]。

3結(jié)束語

本文主要通過Simufact. welding軟件，對B型鋁合金地鐵枕梁焊接數(shù)值模擬仿真分析。分析枕梁結(jié)構(gòu)在焊接過程的溫度場、應(yīng)力和變形情況。在實際工藝下，通過對比多層焊道與簡化后單層焊道的殘余應(yīng)力可知，單層焊、三層焊和多層焊殘余應(yīng)力趨勢相似，數(shù)值差別不大，多層焊道簡化為單層焊道是可行的，可以大大減少工作時間;在焊接應(yīng)力場的分析中，焊接裝夾位置處會出現(xiàn)較大的等效應(yīng)力，產(chǎn)生應(yīng)力集中，除裝夾位置外，應(yīng)力較大的區(qū)域出現(xiàn)在焊縫熱影響區(qū)。平行于焊縫方向的應(yīng)力（縱向應(yīng)力）會大于垂直焊縫方向的應(yīng)力（橫向應(yīng)力），焊接順序?qū)堄鄳?yīng)力的大小產(chǎn)生一定的影響，兩個距離較近的焊縫，先焊接焊縫周圍的殘余應(yīng)力小于后焊接焊縫周圍的殘余應(yīng)力，在實際焊接過程中，可以把重要的部位或者容易集中的部位先進(jìn)行焊接，以減小殘余應(yīng)力。該研究為地鐵枕梁的焊接提供了理論依據(jù)。

參考文獻(xiàn)：

[1]孫寧.我國城市軌道交通行業(yè)發(fā)展現(xiàn)狀與對策[J]. 城市軌道交通， 2008， （4）： 51-55.

[2]王文靜， 惠曉龍， 馬紀(jì)軍. 高速列車設(shè)備艙支架疲勞裂紋機(jī)理研究[J]. 機(jī)械工程學(xué)報， 2015， 51（6）： 142-147.

[3]李廣權(quán)， 劉明志， 王文靜， 等. 高速動車組齒輪箱疲勞裂紋機(jī)理研究[J]. 機(jī)械工程學(xué)報， 2017， 53（2）： 99-105.

[4]曹蘭勝. 國內(nèi)外地鐵文化建設(shè)之于合肥的經(jīng)驗借鑒與發(fā)展路徑[J]. 重慶交通大學(xué)學(xué)報： 社會科學(xué)版， 2020， 20（2）： 1-2.

[5]施仲衡. 建設(shè)符合我國國情的地下鐵路[J]. 都市快軌交通， 2013， 26（2）： 1-2.

[6]海邦軍. 鋁合金車體設(shè)計研究[J]. 鐵道車輛， 2003， 41（10）： 26-29.

[7]朱信科. 客車車體結(jié)構(gòu)強(qiáng)度和常用材料探討[J]. 鐵道車輛， 2009， 47（6）： 33-35.

[8]李曉敏. 中國鋁型材市場及未來發(fā)展趨勢[J]. 四川有色金屬， 2010（4）： 1-5.

[9]馬思群， 金輝， 謝素明， 等. 高速動車組KK端多道焊焊接殘余應(yīng)力仿真研究[J]. 鐵道學(xué)報， 2015， 37（7）： 27-31.

[10]王士軍， 孫增光， 周永鑫， 等. 基于ANSYS的焊接溫度場仿真分析[J]. 熱加工工藝， 2020， 49（9）： 120-121， 129.

[11]Zheng L， Launert B， Hartmut P. Application of the stochastic finite element method in welding simulation [J]. Welding in the World， 2018， 62（5）： 905-912.

[12]牛煒霖. 論軌道交通鋁合金車體設(shè)計[J]. 上海電氣技術(shù)， 2020， 13（1）： 34-37.

[13]王元良， 周友龍， 胡久富， 等. 鋁合金運(yùn)載工具輕量化及其焊接技術(shù)的發(fā)展[J]. 鋁加工， 2005， 35（9）： 14-18.

[14]陳祝年. 焊接工程師手冊[M]. 北京： 機(jī)械工業(yè)出版社， 2002.

[15]岳彩昂， 李東風(fēng)， 葛懷普， 等. 全焊接鋁合金B(yǎng)型地鐵車體焊接變形的研究[J]. 金屬加工， 2016（20）： 40-42.

[16]王亞男. 高速動車組鋁合金車體T形焊縫接頭的SYSWELD數(shù)值模擬研究[J]. 焊接技術(shù)， 2019， 48（1）： 32-36.

[17]方洪淵. 焊接結(jié)構(gòu)學(xué)[M]. 北京： 機(jī)械工業(yè)出版社， 2008.

[18]程久歡， 陳俐， 于有生. 焊接熱源模型的研究進(jìn)展[J]. 焊接技術(shù)， 2004， 33（1）： 13-15.

[19]莫春立， 錢百年， 國旭明， 等. 焊接熱源計算模式的研究進(jìn)展[J]. 焊接學(xué)報， 2001， 22（3）： 93-96.

[20]肖馮， 米彩盈. T型角接頭焊接熱源模型研究[J]. 電焊機(jī)， 2010， 40（6）： 40-45.

[21]劉陽， 于云海， 祁建偉， 等. 120A型地鐵轉(zhuǎn)向架側(cè)梁焊接變形控制工藝研究[J]. 焊接技術(shù)， 2019， 48（4）： 80-82.

Abstract：? Aiming at the influence of residual stress on train performance and quality in the bolster welding， based on Simufact. welding software， the numerical simulation of welding was carried out on the bolster in the metro vehicle aluminum car body to solve practical engineering problems and the residual stress and deformation distribution after welding were analyzed. Different heat source model is created according to different modes of welded joints and the simplification of actual working conditions of multilayer welding bead， multilayer welding bead is replaced by a single welding bead to get the welding thermal cycle curve. The simplified heat source model is used instead of the conventional double ellipsoid heat source model to simulate the whole bolster， and the welding temperature field and residual stress are obtained. The simulation results show that due to the constraint of the position nodes near the weld， the weld shrinks when the welding is completed and enters the cooling process. Because there is no constraint at the middle of the bolster profile and the arc plate far away from the weld， the deformation is large. The research accurately simulates the change of welding stress in actual production， and provides technical support for the welding of bolster.

Key words： simufact. welding; bolster; heat source model; thermal cycling curve; simulation analysis