姬麗森，凌澤民，何 建

(重慶大學(xué)材料科學(xué)與工程學(xué)院， 重慶 400044)

?

Inconel 625鎳基合金管道焊接殘余應(yīng)力的數(shù)值模擬

姬麗森，凌澤民，何 建

(重慶大學(xué)材料科學(xué)與工程學(xué)院， 重慶 400044)

摘要：以Inconel 625鎳基合金管道的焊接接頭為研究對象，建立了熱-力學(xué)耦合的三維有限元模型，采用ANSYS軟件對該合金管道環(huán)焊縫對稱焊的殘余應(yīng)力進(jìn)行了數(shù)值模擬，分析了管道外表面軸向和環(huán)向殘余應(yīng)力分布，并進(jìn)行了試驗(yàn)驗(yàn)證；此外，還分析了預(yù)熱溫度對管道殘余應(yīng)力的影響。結(jié)果表明：該合金管道焊后外表面軸向與環(huán)向殘余應(yīng)力的模擬結(jié)果與試驗(yàn)結(jié)果在數(shù)值和分布趨勢上均比較吻合，證明了模型的準(zhǔn)確性；在焊縫及近焊縫區(qū)的管道外表面形成了軸向壓應(yīng)力和環(huán)向拉應(yīng)力，隨著距焊縫中心距離的增加，軸向壓應(yīng)力逐漸變?yōu)槔瓚?yīng)力，而環(huán)向拉應(yīng)力逐漸轉(zhuǎn)變?yōu)閴簯?yīng)力，并最終趨向于0；隨著預(yù)熱溫度的升高，管道外表面軸向和環(huán)向殘余應(yīng)力均降低。

關(guān)鍵詞：Inconel 625鎳基合金；管道；預(yù)熱；殘余應(yīng)力；數(shù)值模擬

0引言

Inconel 625鎳基合金具有較高的強(qiáng)度和塑性，優(yōu)異的疲勞和熱疲勞強(qiáng)度、氧化抗力、耐腐蝕性能,優(yōu)良的焊接性，廣泛應(yīng)用于航空航天、船舶、核能和石油化工等領(lǐng)域[1-3]。該合金制造的管道在應(yīng)用時需要進(jìn)行焊接，而焊接是一個局部不均勻加熱和冷卻的過程，易在焊接接頭處產(chǎn)生殘余應(yīng)力，影響到管道的強(qiáng)度、剛度和承載能力，并誘發(fā)應(yīng)力腐蝕[3]。因此，焊接殘余應(yīng)力受到了人們的廣泛關(guān)注。

有限元模擬是伴隨著計算機(jī)技術(shù)和數(shù)值分析方法發(fā)展而發(fā)展起來的結(jié)構(gòu)分析方法，被廣泛地應(yīng)用于焊接熱傳導(dǎo)分析、應(yīng)力和變形分析、焊接結(jié)構(gòu)斷裂分析等方面。通過有限元軟件可對焊接殘余應(yīng)力進(jìn)行精確模擬,從而優(yōu)化焊接結(jié)構(gòu)設(shè)計以及工藝參數(shù)，相比試驗(yàn)方法，節(jié)約了大量的時間和人力物力。國內(nèi)外學(xué)者對鎳基合金焊接殘余應(yīng)力的有限元模擬已開展了一些研究，Yegaie等[4]采用有限元方法研究了Monel 400板塊散熱鎢極惰性氣體保護(hù)焊(GTAW)焊接溫度場與殘余應(yīng)力分布；Purmohamad等[5]利用ANSYS有限元軟件研究了GTAW焊接Incoloy 800H管道對接環(huán)焊縫殘余應(yīng)力的分布；肖新華等[6]利用數(shù)值模擬方法預(yù)測了鎢極惰性氣體保護(hù)焊方法打底焊接大口徑Inconel 718鎳鉻合金管道的焊接溫度場和殘余應(yīng)力的分布。但對Inconel 625鎳基合金管道的焊接殘余應(yīng)力的研究報道相對較少。為此，作者采用ANSYS軟件對Inconel 625鎳基合金管道焊接過程進(jìn)行了有限元建模，模擬了管道環(huán)焊縫外表面軸向和環(huán)向殘余應(yīng)力的分布，并進(jìn)行了試驗(yàn)驗(yàn)證。

1試驗(yàn)方法

對外徑150 mm、壁厚7 mm、長300 mm 的Inconel 625鎳基合金管道分兩層進(jìn)行環(huán)焊縫焊接。第一層采用鎢極氬弧焊焊接，以ERNCrMo-3焊絲作為填充材料；第二層采用手工電弧焊焊接，填充材料為ENCrMo-3焊絲。管道的環(huán)縫對接焊如圖1所示，焊接工藝參數(shù)如表1所示。采取對稱焊的方法，焊接方向用時鐘位置表示為6→3→12，6→9→12,預(yù)熱溫度為20 ℃。

圖1　Inconel 625合金管道環(huán)焊縫焊接示意Fig.1　Schematic for the girth welding of Inconel 625 alloy tube

表1　焊接工藝參數(shù)Tab.1　Welding process parameters

采用X-350AL型X射線應(yīng)力測定儀測試Inconel 625鎳基合金管道外表面軸向和環(huán)向的主應(yīng)力(即殘余應(yīng)力)，采用鉻靶，Kα射線，選(211)面作為衍射晶面，入射角ψ0分別為0°和35°，衍射角2θ為151°～162°。殘余應(yīng)力測試點(diǎn)位于管道外表面軸線上(穿過圖1中9點(diǎn)位置)，距焊縫中心線的距離分別為0,7，17，27 mm，每個點(diǎn)處測6次，去掉最大值和最小值取平均值作為該點(diǎn)處的殘余應(yīng)力。測試前先對管道外表面測試點(diǎn)處進(jìn)行打磨，去除表面的氧化皮和油污等，再用NaCl飽和溶液進(jìn)行電解拋光。

2焊接過程的有限元模型

2.1　三維有限元模型

采用有限元法對管道環(huán)焊縫焊接過程進(jìn)行分析。為準(zhǔn)確模擬焊接的溫度場和應(yīng)力場，采用三維實(shí)體模型進(jìn)行計算，并對其進(jìn)行網(wǎng)格劃分。從焊縫中心到遠(yuǎn)離焊縫邊緣采取由緊密到稀疏的網(wǎng)格過渡方式劃分網(wǎng)格，這樣既可以節(jié)省計算時間，又可以提高計算精度。單元類型采用8節(jié)點(diǎn)的六面體單元，模型的網(wǎng)格劃分如圖2所示，模型節(jié)點(diǎn)數(shù)為29 100，單元數(shù)為21 900。

圖2　管道環(huán)焊縫焊接有限元模型Fig.2　Finite element model for the girth welding of the tube

2.2　溫度場計算模型

在計算中采用基于生死單元的焊接熱源加載模型，通過施加生熱率載荷模擬熱源的移動以及焊縫的形成。將全部焊接熱Q均勻分布在焊縫上[7]，假定在計算前所有焊縫單元不是激活的，即開始計算前將焊縫中所有單元“殺死”，相當(dāng)于焊前管道的組合裝配狀態(tài)，在計算過程中按照焊接順序?qū)⒈弧皻⑺馈钡膯卧凹せ睢?，以此來模擬焊縫金屬的填充過程；同時，給激活的單元施加生熱率(HGEN)，熱載荷的作用時間等于實(shí)際焊接時間。每一步計算完成后，刪除該步的生熱率，重新進(jìn)行下一步的加載運(yùn)算[8]。生熱率的計算公式為：

(1)

(2)

式中：HGEN為每個載荷步施加的生熱率；Q為熱源熱量；U為焊接電壓；I為焊接電流；η為電弧熱效率，取0.75；Aweld為焊縫的橫截面積；v為焊接速度；dt為每個載荷步的時間步長。

隨著焊接過程的逐步進(jìn)行，熱源也隨之在被焊工件上方不斷地向前移動，熱源的移動及數(shù)學(xué)模型通過ANSYS的APDL語言定義循環(huán)來實(shí)現(xiàn)。

計算過程中假定母材與焊縫金屬的熱物理性能參數(shù)相同，不同溫度下導(dǎo)熱系數(shù)λ[1]、比熱容c[1]和熱交換系數(shù)αc[9]如圖3所示。

圖3　Inconel 625合金的熱物理性能Fig.3　Thermophysical characteristic of Inconel 625 alloy

2.3　力學(xué)模型

采用ANSYS有限元軟件進(jìn)行應(yīng)力計算時，對溫度場和應(yīng)力場采用間接耦合方式，在溫度場計算結(jié)束后將模型的傳熱單元改為結(jié)構(gòu)單元，將瞬態(tài)溫度場(熱載荷)按單元節(jié)點(diǎn)一一對應(yīng)施加于應(yīng)力場計算過程中，從而實(shí)現(xiàn)應(yīng)力場與溫度場的耦合[10]。計算時設(shè)定材料模式為熱彈塑性，符合Von Mises屈服準(zhǔn)則，塑性區(qū)遵從流變法則，且各向同性硬化。Inconel 625合金在各個溫度下的彈性模量E、剪切模量G、線膨脹系數(shù)ɑ、屈服強(qiáng)度σs、泊松比ν等力學(xué)性能[1]如圖4所示。

圖4　Inconel 625合金的力學(xué)性能Fig.4　Mechanical properties of Inconel 625 alloy

3模擬結(jié)果與分析

3.1　溫度場

圖5為預(yù)熱溫度為20 ℃、按6→9→12方向焊接時，9點(diǎn)位置處焊縫中心上一個采樣點(diǎn)的溫度隨時間變化的曲線。由圖5可知，該采樣點(diǎn)共經(jīng)過兩次升降溫過程，焊接第一層時，受到加熱作用其峰值溫度達(dá)到1 642 ℃，超過了Inconel 625合金的熔點(diǎn)，可見該采樣點(diǎn)在第一道焊縫內(nèi)；第二層焊在管外壁進(jìn)行，熱源未直接作用于該采樣點(diǎn)，因此峰值溫度在700 ℃左右。這反映了多層焊熱循環(huán)的特點(diǎn)。

圖5　20 ℃預(yù)熱溫度下焊縫中心線上采樣點(diǎn)處的焊接熱循環(huán)曲線Fig.5　Welding thermal cycle curve of the sampling point at theweld center line with preheating temperature of 20 ℃

3.2　應(yīng)力場

因管道上6點(diǎn)和12點(diǎn)處經(jīng)歷了兩次加熱過程，殘余應(yīng)力分布異常[11]，而且對稱焊接方式使6→3→12和6→9→12兩個半周管道殘余應(yīng)力分布趨勢相同，因此選取9點(diǎn)附近區(qū)域進(jìn)行了軸向和環(huán)向殘余應(yīng)力的測試與模擬。圖6為不同預(yù)熱溫度下管道焊后殘余應(yīng)力的分布，預(yù)熱溫度分別為20,150,300 ℃，試驗(yàn)值是在室溫下測得的殘余應(yīng)力。

圖6　管道外表面軸向和環(huán)向殘余應(yīng)力模擬與試驗(yàn)結(jié)果Fig.6　Axial (a) and hoop (b) residual stress distributions by simulation and experiment on outer surface of tube

由圖6可見，模擬管道外表面軸向和環(huán)向殘余應(yīng)力的分布趨勢與試驗(yàn)值較一致；軸向殘余應(yīng)力在焊縫及近縫區(qū)表現(xiàn)為壓應(yīng)力，模擬和試驗(yàn)得到的軸向最大壓應(yīng)力分別為492 MPa和480 MPa，隨著距焊縫距離增加，壓應(yīng)力逐漸減小，并轉(zhuǎn)變?yōu)閿?shù)值很小的拉應(yīng)力，最終應(yīng)力趨于0；環(huán)向殘余應(yīng)力在焊縫及近縫區(qū)為拉應(yīng)力，模擬和試驗(yàn)得到的環(huán)向最大拉應(yīng)力分別為280 MPa和290 MPa，隨著距焊縫距離增加，拉應(yīng)力快速轉(zhuǎn)變?yōu)閴簯?yīng)力，模擬和試驗(yàn)得到的環(huán)向最大壓應(yīng)力分別為316 MPa和297 MPa，隨后壓應(yīng)力逐漸減小至0。這是由焊接時管道受到的局部非均勻加熱和冷卻所導(dǎo)致的。隨著焊接過程的進(jìn)行，焊接熱源沿環(huán)焊縫焊接線移動，焊縫金屬熔化后依次凝固，造成管道內(nèi)部約束的變化；焊后環(huán)焊縫沿徑向進(jìn)行環(huán)向收縮，在焊縫與近縫區(qū)形成向內(nèi)凹的收縮變形，從而在管道軸線方向引起彎曲應(yīng)力。因此，管道焊接后在焊縫及近焊縫區(qū)的外表面形成軸向壓應(yīng)力和環(huán)向拉應(yīng)力，隨著遠(yuǎn)離焊縫，軸向壓應(yīng)力逐漸變?yōu)槔瓚?yīng)力，而環(huán)向拉應(yīng)力逐漸轉(zhuǎn)變?yōu)閴簯?yīng)力，并且隨著距焊縫中心距離的增加軸向和環(huán)向殘余應(yīng)力趨于0。

管道焊后外表面軸向與環(huán)向殘余應(yīng)力的模擬結(jié)果與試驗(yàn)結(jié)果在數(shù)值和分布趨勢上均比較吻合，說明建立的有限元模型可以比較精確地模擬管道對接焊后外表面的軸向與環(huán)向殘余應(yīng)力。由圖6還可發(fā)現(xiàn)，在不同預(yù)熱溫度下，管道軸向與環(huán)向殘余應(yīng)力的分布趨勢基本一致；隨預(yù)熱溫度升高，軸向殘余應(yīng)力和環(huán)向殘余應(yīng)力的最大拉、壓應(yīng)力均逐漸降低，預(yù)熱溫度為20,150,300 ℃的軸向最大壓應(yīng)力分別為492，450，337 MPa;環(huán)向最大拉、壓應(yīng)力分別由20 ℃的280，316 MPa 依次降低到150 ℃下的254，290 MPa和300 ℃下的223，170 MPa，最大壓應(yīng)力降低幅度較最大拉應(yīng)力的更明顯。焊前對管道進(jìn)行整體預(yù)熱，使焊接過程中焊縫區(qū)與遠(yuǎn)離焊縫區(qū)的溫差減小，可有效降低熱應(yīng)力。

4結(jié)論

(1) 建立了Inconel 625鎳基合金管道環(huán)焊縫對稱焊的三維實(shí)體模型，模擬得到焊后管道外表面軸向與環(huán)向殘余應(yīng)力分布，模擬結(jié)果與試驗(yàn)結(jié)果在數(shù)值和分布趨勢上均比較吻合，模型較準(zhǔn)確。

(2) 焊縫及近焊縫區(qū)的管道外表面軸向殘余應(yīng)力為壓應(yīng)力，環(huán)向殘余應(yīng)力為拉應(yīng)力；隨著距焊縫中心距離的增加軸向壓應(yīng)力逐漸轉(zhuǎn)為拉應(yīng)力,環(huán)向拉應(yīng)力轉(zhuǎn)為壓應(yīng)力，并且最終趨于0。

(3) 隨著預(yù)熱溫度的升高，Inconel 625鎳基合金管道外表面的軸向和環(huán)向殘余應(yīng)力降低。

參考文獻(xiàn):

[1]Special Metals. Inconel alloy 625[Z/OL]. (2003-08-13)[2014-09-15]. http://www.specialmetals.com/documents/Inconel%20alloy%20625.pdf.

[2]孟慶森,張柯柯.金屬焊接性基礎(chǔ)[M].北京:化學(xué)工業(yè)出版社, 2010.

[3]陸文江,中尾嘉邦,筱崎賢二.鎳基合金焊接熱影響區(qū)的液化裂紋敏感性[J].焊接學(xué)報, 1993,14(3):186-194.

[4]YEGAIE Y S, KERMANPUR A, SHAMANIAN M. Numerical simulation and experimental investigation of temperature and residual stresses in GTAW with a heat sink process of Monel 400 plates[J].Journal of Materials Processing Technology, 2010, 210:1690-1701.

[5]PURMOHAMAD H, KERMANPUR A, SHAMANIAN M.Numerical simulation and experimental investigation of residual stresses in the circumferential butt GTAW of incoloy 800H pipes[J].Journal of Materials Engineering and Performance, 2010, 19(1):13-21.

[6]肖新華,李輝燕.鎳鉻基合金管結(jié)構(gòu)GTAW打底焊接變形和殘余應(yīng)力的有限元預(yù)測[J]. 熱加工工藝, 2013,42(13):171-173.

[7]王長利.焊接溫度場和應(yīng)力場的數(shù)值模擬[D].沈陽:沈陽工業(yè)大學(xué), 2005.

[8]陳家權(quán),沈煒良,尹志新,等.基于單元生死的焊接溫度場模擬計算[J].熱加工工藝, 2005,34(7):64-65.

[9]RYBICKI E F. Computation of residual stresses due to multipass welds in piping systems[J]. Journal of Pressure Vessel Technology, 1979,101:149-154.

[10]李燦,傅定發(fā),王冠,等.焊接順序?qū)?061-T6鋁合金矩形焊縫焊接殘余應(yīng)力及變形量的影響[J].機(jī)械工程材料, 2012,36(7):88-92.

Numerical Simulation of Welding Residual Stress

in Inconel 625 Nickel-based Alloy Tube

JI Li-sen, LING Ze-min, HE Jian

(School of Material Science and Engineering, Chongqing University, Chongqing 400044, China)

Abstract:Based on the butt joint of Inconel 625 nickel-based alloy tube, the three-dimensional thermal and mechanical finite element model was established. The ANSYS software was used to simulate the residual stress of the nickel-based alloy tube after girth welded by the symmetry welding method, and the simulated distribution of axial and hoop residual stress on outer surface of the tube were analyzed and verified by experiment. In addition, the effect of the preheating temperatures on the residual stresses was investigated. The results show that the axial and hoop residual stress distribution on outer surface of the tube obtained by simulation coincided well with the experimental results, indicating the accuracy of the simulation model. The axial residual compressive stresses and hoop residual tensile stresses formed on outer surface of the weld and near weld zone. With the increase of the distance from the weld center line, the axial compressive stress gradually turned to tensile stress and hoop tensile stress to compressive stress, and then tended to zero. The axial and hoop residual stresses both decreased with the increase of the preheating temperature.

Key words:Inconel 625 nickel-based alloy; tube; preheating; residual stress; numerical simulation

中圖分類號：TG404

文獻(xiàn)標(biāo)志碼：A

文章編號：1000-3738(2015)12-0067-04