熱依汗古麗·木沙

（新疆大學機械工程學院，新疆烏魯木齊830047）

油田運輸管道中的彎管法蘭焊接仿真及優(yōu)化分析

熱依汗古麗·木沙

（新疆大學機械工程學院，新疆烏魯木齊830047）

基于Sysweld軟件平臺，采用雙橢球體熱源分布模型，實現(xiàn)對油田運輸管道中的彎管法蘭焊接過程的數(shù)值模擬仿真。使用UG建立彎管法蘭的三維焊接模型，利用Hypermesh、Visual-Weld和Sysweld等進行網(wǎng)格劃分、分組和模擬仿真，獲取了溫度場及應(yīng)力應(yīng)變場的分布情況。通過分析焊接溫度場分布和焊件形變得出溫度場對焊件殘余應(yīng)力的影響，優(yōu)化焊接裝卡條件后，改善焊件殘余應(yīng)力情況。

彎管法蘭；數(shù)值模擬仿真；溫度場；應(yīng)力應(yīng)變場；焊件殘余應(yīng)力

在利用管道運輸液、氣體時常有改變運輸方向的情況，管道與管道的連接一般采用法蘭盤連接，還可以通過螺栓預(yù)緊力來確保管道在額定壓強下運輸液、氣體的穩(wěn)定性。絕大多數(shù)管道運輸?shù)氖且?、氣體，而且管道一般都埋在環(huán)境比較惡劣的地方，因此，對彎管法蘭焊接質(zhì)量有較為嚴格的要求，為確保焊接質(zhì)量，采用有限元或有限差分的方法對焊接應(yīng)力進行預(yù)測，這樣可以縮短實驗周期和節(jié)約生產(chǎn)成本[1-2]。也就是說，采用計算機模擬仿真的方法來模擬焊接應(yīng)力，能夠帶來較大的經(jīng)濟效益。

本文對彎管法蘭焊接過程進行了模擬仿真，利用UG、Hypermesh、Visual-Weld、Sysweld等進行建模、劃分網(wǎng)格、對焊接組件及裝夾條件等進行分組，并設(shè)置載入材料屬性函數(shù)、熱源參數(shù)及網(wǎng)格等影響因素，進行求解觀察結(jié)果[3-5]。同時結(jié)合模擬仿真分析結(jié)果，改變焊接過程中加載約束，優(yōu)化彎管法蘭裝卡狀態(tài)，以獲得更高的焊接質(zhì)量。

1 彎管法蘭焊接的數(shù)值仿真

1.1 有限元模型的建立、建模、網(wǎng)格劃分、分組

焊接模型選用國標90°彎頭鋼管，材料Q235-A；用牌號為J422型號為E4303的焊條；選用與彎頭相對應(yīng)的型號為DN50平面法蘭盤，單面焊接，雙面螺母連接，對應(yīng)的工作壓力PN=1.8 MPa.選取焊腳尺寸為7 mm，彎管法蘭焊接中彎管壁厚為8 mm，焊腳厚度為7 mm.

此次焊接過程用雙焊槍對稱焊接，其焊接路徑如圖1所示。為了更好地保證結(jié)合點的焊接工藝，采用過八度焊接。

圖1 焊接路徑

UG建模后，導出的IGES格式模型裝入Hypermesh進行網(wǎng)格劃分：此次研究對象為簡化模型留下1/4體；并劃分密網(wǎng)格和疏網(wǎng)格區(qū)域，對模型進一步切割分塊，為保證網(wǎng)格的連續(xù)性，需在切割完后的各體之間采用過度連接的方法[6-7]。

在試驗中采用Visual Weld軟件對焊接件進行分組處理，在數(shù)控焊接機床上用兩只焊槍進行焊接，用Visual Weld軟件完成整個焊接過程中的各點、軌跡線、參考線及機械邊界定義等條件[8]。

1.2 焊接向?qū)?/p>

焊接向?qū)椋翰牧希ǖ秃辖鸶邚姸冉Y(jié)鋼）函數(shù)（w_s355j2g3.mat）；熱源函數(shù)（hsf.fct），采用雙橢球型熱源模型；散熱函數(shù)（CONVECTIVE AND RADIATIVELOSSES mm），采用自然冷卻，邊界溫度30℃；導入網(wǎng)格（T_shapeDATA100.TIT）；用焊接起始點，焊接線，參考線等焊接路徑進行焊接操作；Ux、Uy、Uz三個方向的約束為邊界。

2 焊接溫度場分布

焊接時刻溫度云圖如圖2（a）、（b）所示，焊接速度為6 cm/min，總焊接時間16 s，空中自然冷卻300 s.從圖中可以看出最高溫度為1 500℃，熔融區(qū)在熱源中心靠近外部的區(qū)域，溫度約為950℃左右，靠內(nèi)區(qū)域為熱影響區(qū)域；熱源的最高溫度基本保持在1 600℃左右，從圖中還可以看出，隨冷卻過程的推移溫度會急劇降低，在300 s時，溫度基本降低到30℃左右。

圖2 （a）49 s時焊接時刻溫度云圖

圖2 （b）300 s時焊接時刻溫度云圖

在14 s左右時的所有單元溫度曲線如圖3所示，焊接中熱源的最高溫度為1 600℃左右，此時溫度曲線有兩個峰值（雙焊槍同時焊接）。焊接件中所有單元體在總焊接時間300 s時的溫度變化曲線如圖4所示，從圖中可以看出，在0~16 s左右時處于焊接狀態(tài)，從16~20 s左右時溫度突降到220℃左右，焊件因溫差較大其溫度迅速下降；20~300 s時為緩慢冷卻過程，當焊件冷卻至室溫整個焊接過程完成。

圖3 14s時所有單元溫度曲線

圖4 焊接單元體整個焊接過程的溫度曲線

從圖5、圖6中可以看到焊接過程中溫度變化曲線及各點的位置和溫度變化情況，還可以看出焊縫熱源邊上的點的最高溫度在110℃左右，距離焊縫較近的點的最高溫度約為270℃左右，在16 s時溫度曲線出現(xiàn)波峰；距離焊縫區(qū)域較遠處的溫度變化相對較小其差距不大。

圖5 選取點和對應(yīng)的溫度變化曲線情況圖

圖6 選取點和對應(yīng)的溫度變化曲線情況圖

3 應(yīng)力應(yīng)變場分析

焊接中經(jīng)過Mechanical Analysis求解可以得出應(yīng)力分布云圖。圖7為300 s時期應(yīng)力場分布情況，從圖中可以看出在熔融區(qū)與固態(tài)區(qū)域的過渡處出現(xiàn)應(yīng)力集中現(xiàn)象，選用低合金高強度結(jié)構(gòu)鋼，其抗拉強度為470~680 MPa，屈服強度大于等于350 MPa，結(jié)果在許用應(yīng)力范圍內(nèi)[9]。還可從圖中可知大多數(shù)點的應(yīng)力分布在0 MPa左右，最大應(yīng)力不超過390 MPa.本次對分析結(jié)果中Displacements進行了查看，結(jié)果如圖8所示，該圖反映在16 s、300 s兩時刻焊件的應(yīng)變情況，從圖中可以看出，冷卻后的最大變形位置、變形量的變化過程和在彎管處出現(xiàn)整體偏移的情況，偏移量雖然不大但會影響后續(xù)安裝過程，對彎管法蘭的密封性也有所影響。

圖7 300s時應(yīng)力場分布情況圖

圖8 單元變形情況圖

4 彎管法蘭焊接的優(yōu)化處理

對彎管法蘭焊接過程采取焊件重新定位、裝卡、法蘭盤底端施加約束，并在彎管端面施加裝卡條件來阻止剛性體運動等優(yōu)化處理。

優(yōu)化裝卡條件時，由于在彎管端面添加載荷，在應(yīng)力減小時會受軸方向影響，應(yīng)力分布曲線與之前曲線相差不大；經(jīng)過添加裝卡條件后，在彎管端面處應(yīng)變有所改善，降低到0.3 mm左右。

通過Sysweld對管道中彎管法蘭焊接件焊接過程進行模擬仿真，可以有效的反應(yīng)焊接過程的溫度場的變化，并由此獲得了熱影響區(qū)的分布[10-12]。分析其溫度場和應(yīng)力應(yīng)變場情況，發(fā)現(xiàn)在使用雙焊槍180°對角焊的方式下，冷卻過程中殘余應(yīng)力最大為389.87 MPa，應(yīng)變最大為0.3 mm，經(jīng)過優(yōu)化裝卡條件后，分析發(fā)現(xiàn)冷卻過程中殘余應(yīng)力減小約20 MPa，應(yīng)變降到最低，此時保證了彎管處的尺寸和角度要求，為后續(xù)安裝做好了準備，而且法蘭盤裝卡孔的位置的形變量也在要求的范圍內(nèi)，對法蘭連接過程也沒有影響，初步達到優(yōu)化效果。

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

[2]高軍義，周建平，吳金強，等.T形焊接接頭溫度場及力場的三維數(shù)值模擬[J].機械設(shè)計與制造，2012（11）：13-15.

[3]王運炎，朱莉.機械工程材料[M].北京：機械工業(yè)出版社，2012.

[4]Chen，Yu-Hua；Wang，Yong.Numerical simulation of ther mal cycle of in-service welding onto active pipeline based on SYSWELD[J].Hanjie Xuebao/Transactions of the China Welding Institution，2007，28（1）：85-88.

[5]李金閣.Q345/2Cr13鋼異種材料焊接過程有限元分析[D].重慶：重慶大學，2012.

[6]Bilenko，G.A.Study of the structure and mechanical proper ties of a weld in a heat exchanger by computer modeling[J]. Metallurgist，2013，57：573-578.

[7]黎超文，王勇，李立英，等.T形接頭的焊接溫度場三維動態(tài)有限元模擬[J].焊接學報，2011，32（8）：33-36.

[8]SYSWELD.Engineering guide of training and Toolbox[M]. France：ESI，Group，2007.

[9]沈松泉，黃振仁，顧竟成.壓力管道安全技術(shù)[M].南京：東南大學出版社，2001.

[10]Syahroni，N；Hidayat，M.I.P.Numerical simulation of weld ing sequence effect on temperature distribution，residual stresses and distortions of T-joint fillet welds[J].Advanced Materials Research，2011：254-259.

[11]歐先印.基于SYSWELD TIG焊數(shù)值模擬[D].南昌：南昌航空大學，2012.

[12]董克權(quán)，劉超英，陳英俊.雙橢球熱源模型加載算法研究[J].機械設(shè)計與制造，2008（11）：60-62.

The Welding Simulation Analysis and Optimization to Pipe Flange in Oil Transportation

Reyihanguli·Musha
（School of Mechanical Engineering，XinJiang University，Urumqi 830047，China）

Based on the sysweld software platform，the double ellipsoid heat source distribution model is adopted to realize the numerical simulation of the welding process of the bend flange in the oil transportation pipeline.The three-dimensional welding model of elbow flange is established by using UG.Meshing，grouping and simulation are carried out by using Hypermesh，Visual-Weld and Sysweld.The distribution of temperature field and stress strain field are obtained.By analyzing the distribution of welding temperature field and the deformation of weldment，the influence of temperature field on residual stress of weldment was obtained，and the residual stress of weldment was optimized after the welding condition was optimized.

pipe flange；numerical simulation；temperature field；stress field；welding residual stress

TP319

A

1672-545X（2017）06-0108-03

2017-03-20

基于數(shù)控加工技術(shù)的螺旋錐齒輪齒面電化學光整加工工藝基礎(chǔ)研究（項目批準號：51665055）作者簡介：熱依汗古麗·木沙（1969-），女，新疆烏魯木齊人，學士，主要研究方向為機械設(shè)計。