劉軼暄，李慧琴，麻永林，邢淑清

(內(nèi)蒙古科技大學(xué) 材料與冶金學(xué)院，內(nèi)蒙古 包頭 014010)

Monel－400合金環(huán)焊溫度場(chǎng)和殘余應(yīng)力場(chǎng)的數(shù)值模擬

劉軼暄，李慧琴，麻永林，邢淑清

(內(nèi)蒙古科技大學(xué) 材料與冶金學(xué)院，內(nèi)蒙古 包頭 014010)

采用ANSYS大型通用有限元分析軟件，利用APDL語言模擬計(jì)算了Monel－400合金環(huán)焊溫度場(chǎng)。繪制出焊縫熔合線及其附近的熱循環(huán)曲線，將熱分析得到的節(jié)點(diǎn)溫度作為體載荷進(jìn)行應(yīng)力場(chǎng)的數(shù)值模擬計(jì)算。結(jié)果表明：外表面節(jié)點(diǎn)x方向的應(yīng)力剛開始時(shí)為拉應(yīng)力，12 s后轉(zhuǎn)變?yōu)閴簯?yīng)力，而中間節(jié)點(diǎn)和內(nèi)表面節(jié)點(diǎn)x方向的應(yīng)力始終為壓應(yīng)力，并大于外表面的應(yīng)力值；外表面節(jié)點(diǎn)在y方向的應(yīng)力為175 MPa，該應(yīng)力值是屈服強(qiáng)度的76%，中間節(jié)點(diǎn)和內(nèi)表面節(jié)點(diǎn)的最高應(yīng)力為拉應(yīng)力，z方向的應(yīng)力分布規(guī)律與x方向的應(yīng)力分布基本相同。

Monel－400合金；溫度場(chǎng)；應(yīng)力場(chǎng)，數(shù)值模擬

0 前言

Monel－400是Ni－Cu固溶體的單相奧氏體合金，在中高溫的強(qiáng)酸、強(qiáng)堿中具有良好的耐腐蝕性，使其在一些化工和核工業(yè)設(shè)備制造中占有重要的位置[1]。Monel－400合金焊接性較好，其焊接結(jié)構(gòu)在國外一些國家早已廣泛采用。我國近年來也在開發(fā)Monel系列合金及其在特殊化工設(shè)備中的應(yīng)用。由于Mone－400合金的熱物理性能與普通鋼材有較大的差異，在焊接現(xiàn)場(chǎng)中容易出現(xiàn)一些焊接缺陷[2]，其中焊接殘余應(yīng)力的大小直接影響焊接質(zhì)量的好壞。

為了能夠較好地預(yù)測(cè)焊接殘余應(yīng)力，在此利用ANSYS有限元軟件對(duì)Monel－400合金圓筒環(huán)焊縫進(jìn)行了數(shù)值模擬計(jì)算。先是計(jì)算溫度場(chǎng)，然后將溫度場(chǎng)作為應(yīng)力場(chǎng)計(jì)算的初始條件進(jìn)行焊接殘余應(yīng)力的計(jì)算，這為研究實(shí)際的焊接殘余應(yīng)力提供了參考依據(jù)。

1 有限元模型的建立

圓筒為Monel－400合金，熱物理參數(shù)如表1所示。實(shí)體模型尺寸：外半徑50 mm，內(nèi)半徑40 mm，壁厚10 mm，筒高100 mm。在坡口底部進(jìn)行打底焊，焊接道次為三道，焊接速度 6.89 mm/s。

表1 Monel－400的熱物理性能

為了節(jié)省計(jì)算時(shí)間，取圓筒1/4，并且將實(shí)體模型分成三個(gè)區(qū)：焊縫區(qū)、過渡區(qū)和母材區(qū)，其中母材區(qū)距離焊縫較遠(yuǎn)，并對(duì)母材區(qū)線性化處理，即遠(yuǎn)離焊縫的母材區(qū)的熱物理參數(shù)為常數(shù)，不隨溫度變化而變化，只有焊縫區(qū)和過度區(qū)在焊接過程中的溫度較高，熱物理參數(shù)變化較大[3]。

為了計(jì)算精確，劃分網(wǎng)格時(shí)焊縫區(qū)和過渡區(qū)細(xì)密(網(wǎng)格尺寸約1.76 mm)，母材區(qū)相對(duì)稀疏，采用八節(jié)點(diǎn)六面體的單元類型，使用映射方法。整個(gè)模型的單元數(shù)為21852個(gè)，節(jié)點(diǎn)數(shù)為25839。有限元網(wǎng)格模型如圖1所示。

圖1 有限元網(wǎng)格模型

2 理論基礎(chǔ)

焊接是一個(gè)局部快速加熱到高溫，并隨后快速冷卻的過程。因此焊接溫度場(chǎng)分析屬于典型的非線性瞬態(tài)熱傳導(dǎo)問題。非線性瞬態(tài)熱傳導(dǎo)問題的控制方程為

式中 c為材料比熱容；ρ 為材料密度；λ為導(dǎo)熱系數(shù)；T為溫度場(chǎng)分布函數(shù)；為內(nèi)熱源強(qiáng)度；t為傳熱時(shí)間。這些參數(shù)中λ，ρ，c都隨溫度變化[4]。

對(duì)Monel－400合金管道結(jié)構(gòu)，選用的彈塑性應(yīng)力應(yīng)變本構(gòu)關(guān)系為：材料的屈服服從Von－Mises屈服準(zhǔn)則，塑性區(qū)符合流變法則，并假設(shè)各向同性硬化，應(yīng)力應(yīng)變?cè)谖⑿〉臅r(shí)間增量?jī)?nèi)呈線性變化[5]。

焊接應(yīng)力應(yīng)變計(jì)算以熱傳導(dǎo)分析為基礎(chǔ)，同時(shí)考慮焊接區(qū)域組織轉(zhuǎn)變時(shí)應(yīng)力應(yīng)變場(chǎng)帶來的影響。材料屬于彈性或塑性狀態(tài)的應(yīng)力應(yīng)變關(guān)系為

式中 [D]為彈性或彈塑性矩陣；｛C｝為與溫度有關(guān)的向量。

3 加載計(jì)算

在模擬計(jì)算過程中，采用完全的牛頓—拉夫森法進(jìn)行瞬態(tài)熱分析。設(shè)置母材區(qū)的初始溫度為25℃，整個(gè)模型的外表面施加對(duì)流邊界條件，由于本研究焊后圓筒處于自然冷卻狀態(tài)，對(duì)流換熱系數(shù)取為9 W/(m2·℃)[4]。

熱源的移動(dòng)依靠分布循環(huán)加載來實(shí)現(xiàn)[5]。本研究熱源的移動(dòng)是通過ANSYS的APDL語言的循環(huán)語句*DO－*ENDDO來實(shí)現(xiàn)。利用生死單元來模擬焊縫金屬的填[6]。

4 計(jì)算結(jié)果和分析

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

在整個(gè)模擬計(jì)算過程中，焊接一道次時(shí)間為12 s，然后冷卻10 s。圖2為第三道次焊接完成時(shí)的整個(gè)焊件的溫度分布云圖，此時(shí)熱源溫度最高為789.256℃，最低溫度 46.946℃。圖3為焊接完成冷卻10 s后的溫度分布云圖，此時(shí)焊件最高溫度為336.241℃，等溫線范圍擴(kuò)大。圖4為冷卻大約9 min后整個(gè)焊件的溫度分布云圖，溫度分布呈橢圓形，焊接接頭的最高溫度為137.226℃，最低溫度為134.233℃，且分布在筒的兩端。

在1/4圓筒的中間截面選取處的一個(gè)面，在該截面上取五個(gè)位置，其中位置3在焊縫熔合線上，位置2和位置4位于位置3兩側(cè)1.7 mm處，位置1和位置5位于位置3兩側(cè)3mm處。其熱循環(huán)曲線如圖5所示。

圖2 56 s時(shí)的溫度場(chǎng)

圖3 66 s時(shí)的溫度場(chǎng)

圖4 566 s時(shí)的溫度場(chǎng)

由圖5可知，剛開始時(shí)，五個(gè)位置的溫度相同都為環(huán)境溫度25℃。6s時(shí)，位置1、位置2、位置3的溫度急速上升，達(dá)到最高為焊接溫度1 600℃，然后開始急速下降，但下降速度較上升速度慢，而位置4的最高溫度達(dá)到490℃，位置5的溫度僅僅達(dá)到150℃。直到28 s時(shí)，五個(gè)位置的溫度趨于相同，為120℃。第二道次時(shí)只有位置3的溫度最高達(dá)到焊接溫度1 600℃，第三道次時(shí)最高溫度為590℃，然后趨于一致。

圖5 同一個(gè)表面圓周上的熱循環(huán)曲線

4.2 應(yīng)力場(chǎng)結(jié)果分析

在圓筒的中間截面選取一個(gè)面，在熔合線的外側(cè)選三個(gè)節(jié)點(diǎn)，節(jié)點(diǎn)10032圓筒的外表面，節(jié)點(diǎn)12814位于內(nèi)部，節(jié)點(diǎn)10653位于圓筒內(nèi)部，節(jié)點(diǎn)分布如圖6所示。

圖6 節(jié)點(diǎn)分布

圖7 為三個(gè)節(jié)點(diǎn)x方向的應(yīng)力分布圖。0～12s隨著熱源逐漸接近節(jié)點(diǎn)，x方向的應(yīng)力逐漸增大，都表現(xiàn)為拉應(yīng)力，從12 s時(shí)拉應(yīng)力急速下降，三個(gè)節(jié)點(diǎn)由拉應(yīng)力轉(zhuǎn)變?yōu)閴簯?yīng)力，外表面壓應(yīng)力最小為32 MPa，內(nèi)表面壓應(yīng)力最大為70 MPa，第一道次焊接完成后開始冷卻，壓應(yīng)力逐漸變小，外表面節(jié)點(diǎn)壓應(yīng)力減少緩慢，內(nèi)表面節(jié)點(diǎn)壓應(yīng)力降低較快。28 s時(shí)外表面節(jié)點(diǎn)x方向應(yīng)力由壓應(yīng)力轉(zhuǎn)變?yōu)槔瓚?yīng)力，中間節(jié)點(diǎn)應(yīng)力表現(xiàn)不明顯，內(nèi)表面節(jié)點(diǎn)x方向應(yīng)力先增大后降低。34 s時(shí)，三節(jié)點(diǎn)應(yīng)力同時(shí)增加，表現(xiàn)為壓應(yīng)力，外表面壓應(yīng)力最小為65 MPa，內(nèi)表面節(jié)點(diǎn)最大為110 MPa，中間節(jié)點(diǎn)為100 MPa。冷卻過程中，三節(jié)點(diǎn)的x方向應(yīng)力都減小。50s時(shí)，三個(gè)節(jié)點(diǎn)應(yīng)力變化和第二道次的應(yīng)力變化相似。56s時(shí)焊接完成開始冷卻，三節(jié)點(diǎn)都表現(xiàn)為壓應(yīng)力。最終該截面處外表面節(jié)點(diǎn)應(yīng)力為拉應(yīng)力，中間節(jié)點(diǎn)和內(nèi)表面節(jié)點(diǎn)均為壓應(yīng)力，其中內(nèi)表面節(jié)點(diǎn)壓應(yīng)力最大為50 MPa。

圖7 x方向瞬時(shí)應(yīng)力分布

y方向瞬時(shí)應(yīng)力分布如圖8所示。由圖8可見，在焊接第一道次時(shí)，內(nèi)表面節(jié)點(diǎn)y方向瞬時(shí)應(yīng)力即徑向瞬時(shí)應(yīng)力表現(xiàn)為拉應(yīng)力，內(nèi)表面節(jié)點(diǎn)表現(xiàn)為壓應(yīng)力，中間節(jié)點(diǎn)變化范圍較小，開始時(shí)為拉應(yīng)力，然后變?yōu)閴簯?yīng)力，并逐漸變小，直到28 s時(shí)轉(zhuǎn)變?yōu)槔瓚?yīng)力，開始冷卻時(shí)，又表現(xiàn)為壓應(yīng)力，一直到焊接完成。外表面節(jié)點(diǎn)的y方向瞬時(shí)應(yīng)力表現(xiàn)為壓應(yīng)力，28 s時(shí)急速增加，增加到最大為175 MPa；到34 s時(shí)(第二道次焊接完成開始冷卻時(shí))壓應(yīng)力急速降低，轉(zhuǎn)變?yōu)槔瓚?yīng)力達(dá)到90 MPa后開始變小。焊接完成后開始冷卻，在該過程中應(yīng)力表現(xiàn)不大。內(nèi)表面節(jié)點(diǎn)表現(xiàn)為壓應(yīng)力，且變化不大，直到焊接完成開始冷卻時(shí)，由拉應(yīng)力轉(zhuǎn)變?yōu)閴簯?yīng)力。外表面和中間節(jié)點(diǎn)應(yīng)力很小。

圖8 y方向瞬時(shí)應(yīng)力分布

圖9 z方向瞬時(shí)應(yīng)力分布

z方向瞬時(shí)應(yīng)力分布如圖9所示。由圖9可見，z方向的瞬時(shí)應(yīng)力分布和x方向的應(yīng)力分布規(guī)律完全一致，但大小不同，z方向的瞬時(shí)應(yīng)力最大為內(nèi)表面節(jié)點(diǎn)，表現(xiàn)為壓應(yīng)力，大小為120MPa。中間節(jié)點(diǎn)也為壓應(yīng)力最大值為105 MPa；而外表面節(jié)點(diǎn)變化不大，剛開始表現(xiàn)出較小的拉應(yīng)力，第一道次焊接完成開始冷卻后，轉(zhuǎn)變?yōu)閴簯?yīng)力，最大值為60 MPa。在第三道次焊接完成開始冷卻時(shí)內(nèi)表面節(jié)點(diǎn)和中間節(jié)點(diǎn)壓應(yīng)力幾乎相等為55 MPa，并逐漸減小，外表面節(jié)點(diǎn)壓應(yīng)力為20 MPa并逐漸減小。

5 結(jié)論

(1)靠近焊縫處的節(jié)點(diǎn)溫度梯度最大，熱影響區(qū)次之，母材區(qū)的溫度梯度最小。

(2)利用單元生死技術(shù)可以有效的模擬焊縫的形成過程和焊接熱載荷的輸入。

(3)外表面節(jié)點(diǎn)x方向的應(yīng)力剛開始時(shí)為拉應(yīng)力，然后又轉(zhuǎn)變?yōu)閴簯?yīng)力，而中間節(jié)點(diǎn)和內(nèi)表面節(jié)點(diǎn)x方向的應(yīng)力始終為壓應(yīng)力，并大于外表面的應(yīng)力值。y方向的應(yīng)力在外表面變化最大，最大可以達(dá)到175 MPa。z方向的應(yīng)力分布規(guī)律和x方向的應(yīng)力分布基本相同。

[1]王貴生.蒙乃爾合金的焊接[J].焊管，2000，23(4)：32－34.

[2]武廣昭.蒙乃爾合金管件的焊接[J].熱加工工藝，2007，36(19)：89－91.

[3]董巷海，劉建華，杜漢斌，等.激光薄板拼焊過程溫度場(chǎng)和應(yīng)力應(yīng)變場(chǎng)的數(shù)值模擬[J].電焊機(jī)，2004，34(9)：21－24.

[3]張文鉞.焊接傳熱學(xué)[M].北京：機(jī)械工業(yè)出版社，1987：4－7.

[4]陳家權(quán)，沈煒良，尹志新，等.基于單元生死的焊接溫度場(chǎng)模擬計(jì)算[J].熱加工工藝，2005(7)：64－65.

[5]胡紅軍，楊明波，張丁非.ANSYS10.0材料工程有限元分析實(shí)例教程[M].北京：電子工業(yè)出版社，2008：147－149.

[6]龔曙光，謝桂蘭.ANSYS操作命令與參數(shù)化編[M].北京：機(jī)械工業(yè)出版社，2004：128－128.

Numerical simulation of monel400 alloy welding temperature field and residual stress field

LIU Yi-xuan，LI Hui-qin，MA Yong-lin，XING Shu-qing
(School of Material and Metallurgy，Inner Mongolia University of Science ＆ Technology，Baotou 014010，China)

The article analyses Monel－400 alloy cylinder circle weld temperature field and residual stress field by finite element，using ANSYS finite element analyse of large-scale general-purpse software，and draw the thermal cycling curves of the weld fusion line and near location.Then the temperature of node in thermal analysis as the result load on the body of the numerical simulation calculation of stress results show that：The stress of the surface direction node of x beginning to tensile stress，then transform for the compressive stress,but the nodes of the intermediate and the inside surface in x direction is always compressive stress，And the stress value more than appearance；The stress in the node of the surface is the largest in the y direction is 175 MPa，it is 76%on the yield strength，the highest tensile stress of intermediate nodes and the inside surface nodes，The distribution of the stress in z direction and x direction are the same.

Monel－400 alloy；temperature field；residual stress field；numerical simulation

TG457.19

A

1001－2303(2011)03－0075－04

2010－11－09

劉軼暄(1982—)，男，內(nèi)蒙古人烏蘭察布市人，碩士，主要從事焊接的數(shù)值模擬工作。