扁平鋼箱梁U肋焊接殘余應力數(shù)值模擬分析

摘要：

已建成的扁平鋼箱梁橋在運營過程中發(fā)現(xiàn)一些與結(jié)構(gòu)、構(gòu)造和制造相關的病害，并且隨著服役時間的增長越來越突出。其原因是扁平鋼箱梁結(jié)構(gòu)采用薄板組焊而成，因此，扁平鋼箱梁在未受力之前，便已存在著由于焊接所產(chǎn)生的應力場和變形場， 導致扁平鋼箱梁普遍出現(xiàn)裂紋等病害。運用熱結(jié)構(gòu)耦合法對扁平鋼箱梁U肋與橋面板交接處焊接的溫度場和應力場進行數(shù)值模擬計算分析，并得到焊接殘余應力的分布規(guī)律，為進一步分析結(jié)構(gòu)病害與加固設計提供依據(jù)。

關鍵詞：

扁平鋼箱梁；焊接；殘余應力；裂紋；數(shù)值模擬

中圖分類號：

TU391

文獻標志碼：A

文章編號：16744764（2015）01010408

Numerical simulation for the residual stress of Urib welding onto the flat steel box girder

Wang Ruolin1，Li Gensen2

（1.School of Civil Engineering， Wuhan University， Wuhan 430072， P.R.China；

2.Yunnan Electric Power Design Institute， Kunming 650224， P.R.China）

Abstract：

The flat steel box girder is commonly used for longspan cablestayed bridges. the structural complexity， construction details and manufactures cause damage to the bridge and the damages seriously grow overtime. The main reason for these damages is that welding technology of flat steel box girder with thin plates and rib elements produces residual stress and strain fields during the welding process， and further more welding cracks commonly appear. In this study， thermalstructural couple methods was adopted to numerically simulate the welding high temperature fields and the stress fields at the joint of U rid and decks of a flat steel box. The welding residual stress distribution was studied as well for further analyzing on the structural damage and strengthening design strategies.

Key words：

flat steel box girder； welding； residual stress； crack； numerical simulation

扁平鋼箱梁因其具有良好的氣動性能而得到了廣泛的運用。中國近年來所建成的大跨度懸索橋和斜拉橋的主梁多采用該種結(jié)構(gòu)形式。但是，扁平鋼箱梁構(gòu)造復雜，加勁肋（U肋）、縱橫隔板布置密集，加之大跨度橋梁體量巨大[1]，并且在扁平鋼箱梁內(nèi)部，腹板與頂（底）板、加勁肋與頂（底）板都采用焊接連接，導致了箱梁內(nèi)部焊縫很多，而焊縫及其熱影響區(qū)由于焊接的不均勻熱過程，導致構(gòu)件中存在較大的焊接殘余應力、應變，幾何不連續(xù)性、力學性能不連續(xù)性和各種焊接缺陷等。當車輛從橋梁上經(jīng)過時，因橋面鋪裝缺損或者伸縮縫的不平整，造成了車輛振動而對橋梁作用交變荷載，在焊縫等薄弱環(huán)節(jié)，極易產(chǎn)生疲勞損傷，形成疲勞裂紋等，導致梁體剛度下降，加劇橋面板上鋪裝層的損傷。

王若林，等：扁平鋼箱梁U肋焊接殘余應力數(shù)值模擬分析

學者們對鋼結(jié)構(gòu)焊接殘余應力進行了大量的研究。在計算機技術推廣以前，因焊接工藝工程復雜多變，采用試驗手段作為基本研究方法，其模式為“理論試驗生產(chǎn)”[2]。到了近代，計算機技術得到推廣，也使數(shù)值模擬遍及了焊接的各個領域[2～20]，焊接研究的模式由“理論試驗生產(chǎn)”轉(zhuǎn)變?yōu)榱恕袄碚摂?shù)值模擬生產(chǎn)”。2000年，Masahito等針對以往計算焊接殘余應力的固有應變分析和熱彈塑性分析只能計算焊件表面應力的不足，提出了一種新的數(shù)值計算方法，對圓管對接焊厚度方向的焊接殘余應力進行了數(shù)值計算，并與中子衍射測得結(jié)果進行對比，實測結(jié)果與計算結(jié)果一致[5]。同年，Dong等針對焊接殘余應力對壓力管道疲勞強度的影響做了研究，并進行了有限元分析[6]； 近幾年對焊接殘余應力的研究除了傳統(tǒng)的試驗研究以外，運用有限元對其進行數(shù)值計算更是研究熱點并取得到長足的進步[710]。

中國對焊接殘余應力的研究，尤其是運用數(shù)值方法對其進行的研究起步較晚。最初是由西安交通大學的樓志文等人在上世紀70年代研究焊接溫度場和熱彈塑性應力場的分析中，引入了數(shù)值分析方法。上世紀80年代后，上海交通大學焊接教研室對焊接的熱傳導做了大量的數(shù)值分析研究，提出了變步長外推法方程用以求解非線性熱傳導。近年來，很多高校和研究機構(gòu)也對焊接展開了數(shù)值模擬的探討[2，5，710]。

筆者針對一實際橋梁扁平鋼箱梁U肋焊縫附近反復出現(xiàn)裂縫，對該焊接過程進行數(shù)值模擬，進而闡明焊接殘余應力的產(chǎn)生原因及其對鋼箱梁等鋼結(jié)構(gòu)所產(chǎn)生的影響。

1焊接過程分析建模

11有限元基本方程

焊接過程是一個局部被急速加熱，后又急速冷卻的過程，呈現(xiàn)的是一個典型的非線性瞬態(tài)熱傳導問題。該熱傳導過程的控制方程為式（1）和（2）。

xλxxTx+yλyyTy+zλzzTz+q=ρcdTdt（1）

dTdt=Vt+VxTx+VyTy+VzTz（2）

式中：c為材料比熱容；ρ為比材料密度；Vx、Vy、Vz為媒介傳導速率。

在對焊接溫度場進行計算時，通常采用以下幾種邊界條件。

第一類邊界條件，邊界上的溫度值已知

λxxTxnx+λyyTyny+λzzTznz=Tsx，y，z，t（3）

第二類邊界條件，邊界上熱流密度分布已知

λxxTxnx+λyyTyny+λzzTznz=qsx，y，z，t（4）

第三類邊界條件，邊界與周圍物體的熱交換已知

λxxTxnx+λyyTyny+λzzTznz=β（Tα-Ts）（5）

式中：qs為單位面積上的外部熱源輸入；β為表面換熱系數(shù)；Tα為邊界已知溫度；Ts為周圍介質(zhì)溫度。

12焊接過程應力場

焊接過程中，隨著溫度的變化，材料的彈性模量、屈服極限強度、線膨脹系數(shù)等都隨著溫度而不斷變化，呈現(xiàn)出明顯的非線性行為，其應力場的分析是以溫度為體荷載的復雜的非線性彈塑性分析。包括如下相關性分析。

1）路徑相關性。

2）率相關性。在對焊接的計算中，按率無關性考慮。

3）塑性分析相關準則。本文分析中，將材料視為均勻的，因此采用的屈服準則為Von Mises準則，對分析中所采用的強化準則，選取雙線性隨動強化準則。

13焊接過程模型簡化

焊接是一個涉及到電弧物理、傳熱、冶金和力學的復雜過程[2]。其應力應變場存在著材料非線性和幾何非線性等非線性問題，為簡化焊接的計算模型，將焊接熱應力場看作材料非線性瞬態(tài)問題，采用彈塑性力學模型。

在焊接計算中只考慮單向耦合，即焊接溫度場對焊接應力場的影響。

分析中主要考慮對流和熱傳導，而熱量的輻射則通過加大對流和熱傳導系數(shù)來近似代替。對焊接熱源的模擬，采用在焊縫相應單元節(jié)點處施加生熱率來進行模擬。

2扁平鋼箱梁U肋焊接過程有限元模擬分析

在對焊接過程的模擬計算中，采用ANSYS熱結(jié)構(gòu)耦合計算功能。

21熱結(jié)構(gòu)耦合

在分析中采用間接耦合法進行焊接殘余應力的模擬。即先進行熱分析，再將熱分析的結(jié)果作為體荷載加到模型中進行應力場分析，從而得到焊接后產(chǎn)生的殘余應力。在熱分析和結(jié)構(gòu)分析中采用相同的網(wǎng)格劃分方式，使單元具有相同的節(jié)點號碼，用ANSYS相關命令將熱單元轉(zhuǎn)換為結(jié)構(gòu)單元，并將熱分析結(jié)構(gòu)文件中計算所得的節(jié)點溫度直接施加到相應的結(jié)構(gòu)分析中的節(jié)點上。

22單元的生死

在焊接過程中，焊接前焊縫金屬是不存在的，隨著焊縫的延伸，焊縫處被金屬所填充，對此過程可采用ANSYS單元的生死來進行焊接過程的模擬：焊接前，將焊縫單元全部殺死，冷凝過程按施焊順序?qū)ⅰ氨粴⑺馈钡膯卧鸩郊せ?

23單元類型

焊接殘余應力的有限元模擬，單元類型的選擇考慮了以下幾點：具有間接耦合功能；具有單元的生死功能；可以對焊縫處進行規(guī)則劃分；還應具有熱和結(jié)構(gòu)單元的轉(zhuǎn)換功能。綜合以上幾點，選定Solid70單元進行溫度場分析，并選用與之對應的Solid185進行應力場分析。單元如圖1所示。

24模型尺寸

頂板厚12 mm，U肋厚8 mm，上口320 mm，下口204 mm，高260 mm，U肋和焊縫如圖2所示。因其對稱性，只取其中心一半建立模型，取面板寬320 mm，長300 mm。對構(gòu)件模型縱向劃分120個單元，橫向在焊縫及近縫區(qū)采用較密的網(wǎng)格劃分，遠離焊縫的位置則采用較粗的網(wǎng)格，劃分網(wǎng)格后，模型包含122 570個單元，如圖3所示。

圖1焊接分析單元

Fig.1Welding analysis elements

圖2扁平鋼箱梁U肋及焊縫（單位：mm）

Fig.2Urib and welding seam of the flat box girder

圖3焊縫有限元模型

Fig.3Welding seam FE modeling

25材料參數(shù)

橋梁U肋及橋面板采用16Mnq（Q345q）鋼。高溫下的鋼材物理性能相關文獻很少，在本文中的計算，需要考慮的鋼材熱物理性能，即材料彈性模量，屈服應力，熱傳導系數(shù)，熱膨脹系數(shù)，比熱等參數(shù)與溫度的變化關系，均參考歐洲規(guī)范Eurocode3[20] Part12相關規(guī)定取用。

26空氣對流系數(shù)

經(jīng)查閱相關資料，常溫下空氣的對流系數(shù)一般為5～25 W/（m2·T），在計算中統(tǒng)一選取為10 W/（m2·T），即10-6 W/（mm2·T）。

2.7焊接熱源

焊接時室溫取為25℃，焊接熱效率取080。相關參數(shù)如表1。

表1焊接熱源參數(shù)

Table 1Heat resources parameters of welding

焊材電流/A電壓/V焊速/（m·h-1）CO2氣流量/（L·min-1）

E501T1280±2028±223±215～20

28焊接過程溫度場分析

281焊接加熱過程溫度場分析使用生熱率沿焊縫進行加熱，加熱過程共分為120個荷載步，每個荷載步又分為兩個子步，收斂精度為0001。運用單元的生死功能沿焊縫按熱源移動方向順序加熱各段中間焊縫單元的節(jié)點，在對下一步施加生熱率時刪除上一步所施加的生熱率。該過程用循環(huán)來實現(xiàn)。

焊接加熱過程溫度場的變化如圖4 所示。

圖4U施焊加熱過程溫度場變化

Fig.4Temperature variation during welding heating process

由圖4可看出，焊接過程在一開始時焊縫處溫度變化較大，急劇升溫，一段時間后穩(wěn)定在一定溫度，本文計算中焊縫穩(wěn)定溫度為2 317℃，并且溫度帶隨著焊接熱源的移動而移動，相對于熱源位置，其溫度分布并不改變。

282焊接冷卻過程溫度場分析焊縫冷卻時的溫度梯度較加熱時小很多，因此在本文計算中對焊縫冷卻過程設定時間為2 h，每個時間步為60 s，每個步長再分為10個子步，其計算結(jié)果如圖5 所示，前3圖分別為冷卻開始1、2、3 min的溫度場云圖，后3圖分別為冷卻1、15和2 h時的溫度分布云圖。由圖中可看出焊接結(jié)束后的幾分鐘之內(nèi)構(gòu)件溫度急劇下降，以后溫度變化越來越慢。冷卻2 h后，最高溫與最低溫相差不超過4℃，此時可認為構(gòu)件降溫已完成。

圖5焊縫冷卻過程溫度場變化

Fig.5Temperature variation during welding cooling process

焊接從加熱到冷卻過程中，距離焊縫不同距離的金屬其溫度變化不同，如圖6所示。

由圖4和圖6可看出，焊縫處溫度一旦達到穩(wěn)定后，其溫度在構(gòu)件橫向幾乎不再產(chǎn)生變化，只是沿著焊縫方向移動而已。并且隨著距焊縫距離的增加，溫度也急劇下降，溫度梯度很大。

圖6距焊縫不同距離節(jié)點溫度變化（1/4、1/2截面）

Fig.6Note temperature variation of different distances from welding seam（respectively，1/4、1/2 crosssection as example）

圖7溫度荷載

Fig.7Temperature loading

3扁平鋼箱梁U肋焊接殘余應力有限元模擬分析

31焊接熱應力

在對構(gòu)件進行溫度場分析的基礎之上，將計算所得溫度作為體荷載加在構(gòu)件上，使用ETCHG命令將熱分析單元轉(zhuǎn)換為相應的結(jié)構(gòu)分析單元進行應力場的分析。因焊接時焊縫及熱影響區(qū)處于高溫狀態(tài)，導致了高度的材料非線性行為，導致計算量大、收斂困難等問題[910]。為此，文中計算采用加大荷載步數(shù)等方法來加強其收斂。計算采用每個荷載步分為40個子步計算。根據(jù)文獻[7] 所述，當溫度高于1 000 ℃時，其對焊接殘余應力的影響并不大，將溫度高于1 000 ℃的單元賦予其溫度為1 000 ℃，以加強收斂，溫度荷載如圖7所示。計算開始時，由于焊縫處于高溫狀態(tài)，產(chǎn)生膨脹對周圍區(qū)域造成壓縮而使焊縫熱影響區(qū)呈現(xiàn)較高的應力，如圖8 所示。由圖8可看出，焊接所產(chǎn)生的熱應力高應力區(qū)集中在熱影響區(qū)，而焊縫處則表現(xiàn)為應力趨近于0。依據(jù)上述簡化原則，焊接結(jié)束時計算所得的熱應力Mises應力分布如圖9所示。由圖中可以看出，此時的最大應力達到了384709 MPa，已超出了鋼材常溫下的屈服極限。

圖8焊縫區(qū)熱應力

Fig.8Heat stress of welding zone

圖9焊接結(jié)束（45 s）時Mises應力云圖

Fig.9Mises stress cloud at the end of welding（45 s）

規(guī)定沿焊縫長度方向為縱向，沿此方向的應力為縱向應力，以σz表示；沿橋面板寬度方向為橫向，沿此方向的應力為橫向應力，以σx表示。構(gòu)件厚度方向上的應力變化很小，圖10 所示為焊接時沿焊接方向的縱向應力和橫向應力，圖中橫坐標為距離（縱向），縱坐標為應力，單位均為MPa。由圖10，在焊縫方向上，熔池部分應力較小且趨于零，熔池前端承受壓應力，而熔池尾部則承受拉應力，隨著熱源的移動，焊縫處各點的應力首先表示為壓應力，而后成為拉應力，并且其最大應力均超過了鋼材在該溫度下的屈服極限，因此必然會產(chǎn)生塑性應變，為焊接殘余應力的產(chǎn)生提供了必要的條件。

圖10沿焊縫應力變化

Fig.10Stress variation along welding seam

32焊接殘余應力

從理論上分析，焊接殘余應力是由于焊接區(qū)以遠高于周圍區(qū)域的溫度迅速加熱，使焊縫區(qū)的材料膨脹，熱膨脹受到周圍較冷金屬的約束，因此產(chǎn)生了熱應力，然而受熱區(qū)溫度升高后屈服極限下降，熱應力便可超過該溫度下鋼材的屈服應力，從而形成了塑性熱壓縮。當冷卻時，溫度急劇下降，造成焊接區(qū)比周圍區(qū)域相對縮短、變窄、縮小，因此，該區(qū)域就呈現(xiàn)出拉應力，而其周圍區(qū)域則呈現(xiàn)壓應力。圖11 為冷卻10 min后的Mise殘余應力分布云圖，由圖中可看出，該時刻在焊縫處存在著較大的殘余應力，其最大值34683 MPa。該時刻縱向殘余應力分布云圖如圖12 所示。由圖中可看出，雖然沿板厚方向殘余應力值不相同，但縱向殘余應力均表現(xiàn)為中間拉應力，兩邊壓應力在焊縫區(qū)和近焊縫區(qū)殘余應力最大，并且距焊縫越近，應力值越大。

圖11冷卻10 min后Mises應力分布云圖

Fig.11Mises stress cloud after 10 min cooling

圖12冷卻后縱向殘余應力分布云圖

Fig.12Longitudinal stress cloud after cooling

冷卻后橫向殘余應力分布云圖如圖13所示.由圖中也可明顯看出冷卻后在焊縫區(qū)存在著較大的殘余拉應力，其值一般都在100 MPa以上，造成了該區(qū)域的應力集中，也為該處的破壞埋下了隱患。

圖13冷卻后橫向殘余應力分布云圖

Fig.13Transversal residual stress distribution after cooling

4結(jié)論

1）運用焊接溫度場分析的基本理論，包括熱分析中經(jīng)典的傳熱方程和有限元基本方程，利用熱結(jié)構(gòu)間接耦合方法，建立了具有位移邊界條件、初始溫度條件、散熱條件的U肋加勁板焊接有限元模型，用ANSYS單元的生死模擬扁平鋼箱梁U肋焊接過程焊縫金屬的熔化和填充，可以得到焊接熱應力和殘余應力場。

2）溫度場的分析結(jié)果表明，在焊接過程中焊縫處形成了穩(wěn)定的溫度場，其變化過程也比較符合實際情況。隨著熱源的移動，焊縫處各點的應力首先表示為壓應力，而后成為拉應力，并且其最大應力均超過了鋼材在該溫度下的屈服極限，必然產(chǎn)生塑性變形。

3）焊接殘余應力的計算結(jié)果表明：殘余應力的高應力區(qū)集中在焊縫及其熱影響區(qū)。

4）對同一截面，無論在縱向和橫向，其焊接殘余應力沿板厚方向分布均有變化。

參考文獻：

[1]劉毅，李愛群，丁幼亮，等.潤揚懸索橋扁平鋼箱梁局部應力的測試與分析[J].東南大學學報，2007，37（2）：275279.

Liu Y，Li A Q，Ding Y L，et al. Computational and experimental studies on stress analysis of flat steelboxgirder of Runyang Suspension Bridge [J]. Journal of Southeast University：Natural Science Edition，2007，37（2）：275279

[2]趙銳.焊接殘余應力的數(shù)值模擬及控制消除研究[D].遼寧大連：大連理工大學，2006

[3]上田幸雄，村川英一，麻寧緒（日）.焊接變形和殘余應力的數(shù)值計算方法與程序[M].成都：四川大學出版社，2008.

[4]Mahin K W，Winters W，Holden T M. Prediction and measurement of residual elastic strain distributions in gas tungsten arc welds [J].Welding Journal（USA），1991，70（9）：245260.

[5]Masahito M，Makoto H，Toshio H.Numerical analysis of welding residual stress and its verification using neutron diffraction measurement [J].Engineering Materials and Technology，2000，122：98103.

[6]Dong P，Brust F W.Welding residual stresses and effects on fracture in pressure vessel and piping components a millennium review and beyond [J].Pressure Vessel Technology，2000，122：329388.

[7]張陜鋒.正交異性板扁平鋼箱梁若干問題研究[D].南京：東南大學，2006.

[8]呂建民，陳懷寧，林泉洪.一種消除小直徑管道焊接殘余應力技術的數(shù)值分析[J].焊接學報，2003（4）：8386

Lyu J M， Chen H N， Lin Q H. Numerical analysis of a method for relieving welding stresses of girthweld pipes with small diameters [J]. Transactions of the China Welding Institution，2003（4）：8386.

[9]周宏，羅宇，蔣志勇，等.基于ANSYS的船舶復雜結(jié)構(gòu)焊接變形預測研究[J].中國船舶，2010，51（2）：5866

Zhou H， Luo Y， Jiang Z Y， et al. Research of complicated ship structure welding distortion based on ANSYS [J]. Ship Building of China，2010，51（2）：5866

[10]陳建波，羅宇，龍哲.大型復雜結(jié)構(gòu)焊接變形熱彈塑性有限元分析[J].焊接學報，2008，29（4）：6972.

Chen J B， Luo Y， Long Z. Analysis on welding distortion of large complicated structure by thermal elasticplastic finite element method [J]. Transactions of the China Welding Institution，2008，29（4）：6972.

[11]Kim S H，Young R C. Effects of alloying elements on mechanical and fracture properties of base metals and simulated heat affected zones of SA508 steels [J].Metallurgical and Materials Transactions，2001，3（4）：903911.

[12]Dong P， Hong J K. The master sn curve approach to fatigue of piping and vessel welds[J]. Welding in the World，2008，48：2836

[13]Lee C H， Chang K H. Finite element computation of fatigue growth rates for model I cracks subjected to welding residual stresses [J].Engineering Fracture Mechanics，2011，78：25052520

[14]張國智，胡仁喜，陳繼剛，等.ANSYS100熱力學有限元分析實例知道教程[M].北京：機械工業(yè)出版社，2007.

（編輯呂建斌）