焊接冷卻時間對錨拉管殘余應(yīng)力影響研究*

陳亞莉 何 佳 肖 祥 董 杰

(武漢理工大學交通學院1) 武漢 430063) (三峽大學土木學院2) 宜昌 443000)

焊接冷卻時間對錨拉管殘余應(yīng)力影響研究*

陳亞莉1)何 佳2)肖 祥1)董 杰1)

(武漢理工大學交通學院1)武漢 430063) (三峽大學土木學院2)宜昌 443000)

斜拉橋錨拉管焊接殘余應(yīng)力將會影響結(jié)構(gòu)的受力性能，為了研究焊接中冷卻時間對它的影響，采用數(shù)值模擬和試驗測量相結(jié)合的方法.運用磁測法和盲孔法對典型對接和T接焊縫模型的殘余應(yīng)力進行測量，同時采用ANSYS有限元軟件建立焊縫的三維模型，試驗測量值驗證了數(shù)值模擬方法的可行性.基于該數(shù)值模擬方法，研究在不同冷卻時間的作用下，斜拉橋錨拉管焊接殘余應(yīng)力的變化規(guī)律.結(jié)果表明，隨著錨拉管焊接冷卻時間的增長，殘余應(yīng)力的分布規(guī)律保持不變，但拉應(yīng)力峰值增大，壓應(yīng)力峰值減小.為了減小殘余應(yīng)力的不利影響，應(yīng)當縮短相鄰焊層間的冷卻時間或者采取預熱和保溫的措施.

道路工程；錨拉管；數(shù)值模擬；焊接殘余應(yīng)力；盲孔法；冷卻時間

0 引 言

在焊接的過程中會產(chǎn)生殘余應(yīng)力，對結(jié)構(gòu)造成不利的影響.目前，對于斜拉橋錨固區(qū)殘余應(yīng)力的試驗測量已有學者做過研究[1-3]，但是試驗測量數(shù)據(jù)有限，不能準確的顯示殘余應(yīng)力的變化情況，所以對焊縫進行數(shù)值模擬是非常有必要的.對焊接殘余應(yīng)力分析也只從焊接層數(shù)和道數(shù)、板件厚度、焊接順序以及焊接的速率、電流和電壓等方面進行參數(shù)分析[4-9].針對多層多道焊接時，每道焊縫焊接完成后的冷卻時間對最終殘余應(yīng)力的影響分析較少.在多層多道焊縫施焊的過程中，兩道焊縫之間的冷卻時間如果過短，道間溫度過高，容易引起熱影響區(qū)晶粒粗大，使焊縫強度及低溫沖擊韌性下降[10-11].如果冷卻時間過長，道間溫度過低，不重新預熱就施焊容易在層間產(chǎn)生冷裂紋[12].施焊工人往往是根據(jù)經(jīng)驗來控制冷卻時間，并沒有相對規(guī)范的時間控制.因此，有必要對焊接的冷卻時間進行研究.

本文以蘭州西固黃河大橋為研究背景，斜拉橋的錨固區(qū)起著連接主梁和斜拉索的作用，是結(jié)構(gòu)的主要受力部分.此處的焊縫密集，如果殘余應(yīng)力得不到有效的控制，對結(jié)構(gòu)受力會有較大影響.因此，對斜拉橋錨固區(qū)焊接殘余應(yīng)力進行預測和研究具有重要的意義.文中采用數(shù)值模擬和試驗測量相結(jié)合的方法，首先驗證了數(shù)值模擬的可行性，然后基于該數(shù)值模擬方法研究冷卻時間對斜拉橋錨拉管殘余應(yīng)力的影響，以便于采取有效的措施控制錨拉管的殘余應(yīng)力.

1 焊接模型試驗

試驗采用磁測法和盲孔法對比測量焊接試件的殘余應(yīng)力.盲孔法是半破壞性測試方法，若直接對錨拉管進行測量，錨拉管將無法投入到橋梁的使用中,因此，采用與該錨拉管相同材料、加工工藝以及同一批工人現(xiàn)場制作的典型對接和T接焊縫模型試件進行殘余應(yīng)力測量.

1.1 模型制作

試驗材料為Q370qE，采用手工電弧焊，焊接電壓為30 V，焊接電流為270 A，焊接速度為0.004 m/s.對接焊縫模型M1的橫截面尺寸、焊道布置見圖1a)，焊道焊接順序為1,2,3→4,5；T接焊縫模型M2的橫截面尺寸、焊道布置見圖1b)，焊道焊接順序為1,2→3,4,5.

圖1 焊縫截面尺寸及焊道布置

1.2 試驗方法

磁測法是基于鐵磁性材料的磁致伸縮效應(yīng)，通過測量磁性變化來測定鐵磁材料中的應(yīng)力大小.首先在對接和T接焊縫試件上打磨拋光，在測量區(qū)域畫出5行5列共25個測點，相鄰測點間距為70 cm.然后采用JH-30殘余應(yīng)力檢測儀對測點進行測量.再采用盲孔法測量進行對比分析.盲孔法是利用JHZK型打孔設(shè)備在平面測點上鉆一個小孔，在該點貼上三向應(yīng)變計R1，R2和R3，通過應(yīng)變計測量出測點在鉆孔過程中釋放的應(yīng)變大小，再根據(jù)相關(guān)理論公式計算出殘余應(yīng)力大小.盲孔法共5個測點，測點縱向位于試件中心，與焊縫方向垂直.相關(guān)理論公式為

(1)

(2)

(3)

式中：ε1，ε2，ε3分別為三個應(yīng)變計R1，R2，R3的測量的應(yīng)變；A，B為應(yīng)變釋放系數(shù)，由試驗測定；θ為殘余應(yīng)力σ1方向與R1軸向的夾角.

2 有限元數(shù)值模擬

為了研究錨拉管殘余應(yīng)力的分布規(guī)律以及不同冷卻時間對錨拉管殘余應(yīng)力的影響，首先建立與實際對接和T接焊縫試件具有相同尺寸和焊道布置的幾何模型，模擬焊接的過程.對比數(shù)值模擬與試驗測量值，發(fā)現(xiàn)兩者結(jié)果較吻合.基于此種模擬方法，建立錨拉管的有限元模型，僅改變冷卻時間參數(shù)的大小，得到最終殘余應(yīng)力的分布和變化規(guī)律.

2.1 模型建立

采用ANSYS軟件建立焊接模型，利用sweep網(wǎng)格劃分單元.焊縫區(qū)源面運用映射方法劃分成加密的網(wǎng)格，遠離焊縫區(qū)源面劃分成較稀疏的網(wǎng)格，過渡區(qū)的單元采用自由網(wǎng)格劃分方式，使兩邊的單元均勻銜接.最終能夠在縮短運算周期的基礎(chǔ)上得到一個準確的計算結(jié)果，模型網(wǎng)格劃分結(jié)果見圖2.

圖2 模型網(wǎng)格劃分結(jié)果圖

該計算模型運用間接耦合單元，即首先利用溫度場單元對焊接模型進行熱分析，再將該單元轉(zhuǎn)化為對應(yīng)結(jié)構(gòu)場單元，進行應(yīng)力應(yīng)變的分析計算.根據(jù)該原則，選用SOLID70單元作為溫度場單元，SOLID95單元作為結(jié)構(gòu)分析單元.

2.2 材料特性

Q370qE結(jié)構(gòu)鋼的熱物理性資料比較少，參考文獻[13]中焊接材料的性能，用于本文的數(shù)值模擬.其導熱系數(shù)、密度、比熱容、泊松比、線性膨脹系數(shù)、屈服應(yīng)力和彈性模量參數(shù)隨溫度變化情況見表1.

2.3 熱源模擬

采用生死單元法模擬焊接填充過程中熱源的輸入，焊接熱量均勻的分布在焊縫上.在建立好焊接幾何模型之后，首先將焊縫單元“殺死”，然后根據(jù)焊接的順序，依次“激活”焊縫單元，給激活的單元施加生熱率.生熱率的計算公式為

表1 Q370qE鋼熱物理性能

(4)

式中：Al為焊縫的橫截面積；dt為每個荷載步的時間步長；k為焊接熱電弧效率，本文取值為0.7；U為焊接電壓；I為焊接電流；v為焊接速度.

2.4 邊界條件

進行溫度場分析時，施加熱力學邊界條件，其中包括恒溫、熱對流和輻射邊界條件.恒溫邊界條件模擬焊接模型周圍空氣環(huán)境溫度，本文取值為25 ℃；熱對流和輻射邊界條件是模擬焊接模型與外界氣體在接觸面上的熱量交換，為了便于計算，考慮總的換熱系數(shù).進行應(yīng)力應(yīng)變場分析時，施加位移邊界條件，目的是防止模型在模擬焊接計算過程中發(fā)生剛體位移，施加的約束不能影響焊件在計算過程中的自由變形.

3 數(shù)值模擬與試驗結(jié)果對比

待對接和T接焊縫試件冷卻至室溫25 ℃之后，選取對接焊縫上表面路徑1路徑上和T接焊縫路徑2路徑上各測點為研究對象，采用磁測法和盲孔法共同測量其殘余應(yīng)力大小，同時提取數(shù)值模擬方法計算出的殘余應(yīng)力大小，圖3～4分別是對接和T接焊縫試驗測量值和數(shù)值模擬計算結(jié)果的比較.其中橫向為垂直于焊縫方向，縱向為平行于焊縫方向.

圖3 對接焊縫殘余應(yīng)力分布

圖4 T接焊縫殘余應(yīng)力分布

由圖4可知，橫向殘余應(yīng)力呈現(xiàn)出W形，該方向應(yīng)力是由焊縫及附近塑形區(qū)域橫向和縱向收縮共同引起的，在焊縫中心處達到拉應(yīng)力峰值52.98 MPa.在過渡區(qū)應(yīng)力值逐漸減小，達到最大壓應(yīng)力值63.95 MPa，該區(qū)域應(yīng)力值不連續(xù)主要是過渡區(qū)單元劃分不均勻?qū)е碌?隨著離焊縫的距離增大，應(yīng)力值逐漸增大，最后趨近于零.縱向殘余應(yīng)力呈現(xiàn)出V形，在焊縫及附近的塑形區(qū)表現(xiàn)為壓應(yīng)力，之后迅速增大再緩慢減小，拉應(yīng)力峰值為154.50 MPa，壓應(yīng)力峰值為275.19 MPa.由圖5可知，T接焊縫在橫向和縱向均呈現(xiàn)出M形，焊縫中心受壓，在焊料和母材連接焊趾區(qū)域處出現(xiàn)最大殘余拉應(yīng)力.這是由于焊料在受高溫膨脹后，隨著焊件的冷卻，焊縫區(qū)收縮變形受到約束，在焊趾處達到最大拉應(yīng)力數(shù)值.應(yīng)力值隨著離焊縫距離的增加而減小，然后緩慢增加，最后趨于零.T接焊縫橫向和縱向的拉應(yīng)力峰值分別為98.20，156.85 MPa，壓應(yīng)力峰值分別為118.93，149.27 MPa.以上的結(jié)論符合一般殘余應(yīng)力的變化規(guī)律.綜合上述對接和T接焊縫的應(yīng)力結(jié)果來看，縱向殘余應(yīng)力的數(shù)值和變化幅度更大，對結(jié)構(gòu)更加不利.

由圖4～5可知，磁測法和盲孔法試驗測量數(shù)據(jù)的變化趨勢和數(shù)值大小均與計算模擬值較吻合.說明該計算模擬方法可以反映模型的真實受力狀態(tài)，模擬值還可以補充試驗中未測點殘余應(yīng)力的分布，所以利用該模擬方法研究冷卻時間對錨拉管殘余應(yīng)力的影響是可行的.

4 冷卻時間對殘余應(yīng)力的影響

斜拉橋錨拉管焊縫橫截面尺寸及焊道布置見圖5，焊縫采用多層多道焊接，焊接順序為1→2→…→10.通過溫度場的模擬結(jié)果發(fā)現(xiàn)，每道焊縫焊接完成之后冷卻16 min時，焊件形成均勻的溫度場，而且與外界環(huán)境溫度接近，此時構(gòu)件內(nèi)部的應(yīng)力可以認為是最終的殘余應(yīng)力，不再隨著冷卻時間的增長而變化[13].因此，選取冷卻時間分別為0，4，8，12和16 min時對應(yīng)的殘余應(yīng)力值進行比較，研究不同冷卻時間對最終焊接殘余應(yīng)力的影響.基于上節(jié)中縱向殘余應(yīng)力對結(jié)構(gòu)受力更為不利的結(jié)論，錨拉管主要研究冷卻時間對縱向殘余應(yīng)力的影響.分別選取路徑3和路徑4路徑上各測點為研究對象，圖5～6分別為測點縱向應(yīng)力值隨位置和冷卻時間的變化情況.

圖5 錨拉管焊縫橫截面尺寸

圖6 縱向殘余應(yīng)力

由圖6a)可知，在路徑3上，焊縫及附近熱影響區(qū)殘余應(yīng)力表現(xiàn)為“山”形，呈現(xiàn)出“拉-壓”交替式變換，在焊縫中心和焊趾處出現(xiàn)拉應(yīng)力峰值，最大值為232.98 MPa，接近材料的屈服強度.同時從殘余應(yīng)力云圖可知，在焊趾處有應(yīng)力集中現(xiàn)象，這主要是因為焊縫區(qū)溫度急劇變化，材料發(fā)生較大的膨脹和收縮變形，在過程中受到周圍區(qū)域的約束，產(chǎn)生了塑形變形，從而引起了較大的殘余應(yīng)力.由圖6b)可知，路徑4上各點在焊縫中部形成了較穩(wěn)定的殘余拉應(yīng)力區(qū)域，在焊縫的邊緣由于邊界效應(yīng)的影響，拉應(yīng)力逐漸變?yōu)閴簯?yīng)力，而且由于焊接順序的先后不同，殘余應(yīng)力的變化程度也有所差異.

由圖6可知，在不同冷卻時間的作用下，錨拉管縱向殘余應(yīng)力的分布規(guī)律保持不變，但應(yīng)力峰值發(fā)生了改變.圖7為縱向最大殘余拉應(yīng)力和壓應(yīng)力隨時間的變化圖，從由圖7可知，拉應(yīng)力峰值隨著冷卻時間的增長而增長，壓應(yīng)力峰值隨著冷卻時間的增長而減小.拉應(yīng)力是引起結(jié)構(gòu)受力不利的主要因素，因此冷卻時間越長，對結(jié)構(gòu)受力越不利.冷卻時間越長，構(gòu)件溫度越低，新填焊料與已焊構(gòu)件的溫差越大，就會引起較大的變形和約束作用，使殘余拉應(yīng)力越大.相反，冷卻時間越短，新填焊料與已焊構(gòu)件形成一個相對溫差較小的溫度場，可以減小材料的收縮和變形，使殘余拉應(yīng)力越小.用Ri/j(i,j=0,4,8,12,16)表示冷卻時間從i變?yōu)閖時，對應(yīng)最大殘余應(yīng)力的變化率，見表2.由表2可見,拉應(yīng)力峰值的變化率R0/4>R4/8>R8/12>R12/16，壓應(yīng)力峰值的變化率在46%左右波動.在第一個4 min內(nèi)，拉應(yīng)力峰值的變化率最大，兩個路徑分別達到19%和81%.這是因為焊件在冷卻的過程中，溫度越高，熱量擴散越快，因此溫度下降的速度越快.所以在第一個4 min的冷卻過程中，已焊構(gòu)件溫度變化最大，與新填焊料的溫差變化率最大，因此殘余應(yīng)力的變化率也最大.綜上所述，冷卻時間的長短會影響結(jié)構(gòu)的溫度，從而引起了殘余應(yīng)力的變化.冷卻時間越長，殘余拉應(yīng)力峰值越大，對結(jié)構(gòu)越不利.可以采取減少冷卻時間，或者適當?shù)谋睾皖A熱措施，來減少新填焊料與已焊構(gòu)件的溫差，從而達到減小殘余拉應(yīng)力的目的.

圖7 應(yīng)力隨時間變化曲線

%

5 結(jié) 論

1) 通過模擬方法，得到對接和T接焊縫的殘余應(yīng)力分布規(guī)律，模擬得到的數(shù)據(jù)與試驗值基本吻合，試驗結(jié)果驗證了數(shù)值模擬方法的可行性.對比對接和T接焊縫縱橫向殘余應(yīng)力的變化規(guī)律，發(fā)現(xiàn)縱向殘余應(yīng)力的數(shù)值和變化幅度更大，對結(jié)構(gòu)受力更為不利.

2) 經(jīng)過模擬數(shù)據(jù)分析，獲得了斜拉橋錨拉管縱向殘余應(yīng)力的分布規(guī)律.垂直于焊縫方向測點的縱向殘余應(yīng)力呈現(xiàn)“拉—壓”交替式分布，在焊縫中心處達到拉應(yīng)力峰值，其值接近材料的屈服強度；位于焊縫上的測點的縱向殘余應(yīng)力在中部形成穩(wěn)定的拉應(yīng)力區(qū)域，兩邊由于邊界效應(yīng)影響，逐漸變?yōu)閴簯?yīng)力.

3) 分析得到冷卻時間對斜拉橋錨拉管殘余應(yīng)力的影響.隨著冷卻時間的增長，殘余應(yīng)力的分布規(guī)律保持不變，但拉應(yīng)力峰值逐漸增大，壓應(yīng)力峰值逐漸減小，對結(jié)構(gòu)越為不利.

4) 冷卻時間越短，殘余拉應(yīng)力峰值的變化率越大.為了減小殘余拉應(yīng)力的大小，應(yīng)減少相鄰焊層間的冷卻時間或者采用預熱和保溫的措施.

[1]劉小渝，王彥博.觀音巖大橋錨拉管區(qū)域焊接應(yīng)力試驗[J].重慶交通大學學報(自然科學版),2008,27(增刊1):875-879.

[2]王彥博.斜拉橋錨拉管區(qū)域焊接殘余應(yīng)力試驗研究分析[D].重慶：重慶交通大學碩，2008.

[3]XU G X, WU C S, MA X Z, et al. Numerical analysis of welding residual stress and distortion in laser+gmaw hybrid welding of aluminum alloy t-joint[J]. Acta Metallurgica Sinica,2013,26(3):352-360.

[4]郭彥林，陳航，袁星.厚鋼板對接焊接三維有限元數(shù)值模擬與分析[J].建筑科學與工程學報,2014,31(1):90-97.

[5]王艷飛，鞏建鳴，蔣文春.焊縫層數(shù)對特厚度管板焊接殘余應(yīng)力與變形影響的有限元分析[J].上海交通大學學報,2013,47(11):1675-1679.

[6]YANG X H, LI Y N, ZhAO W Z. Numerical simulation of the welding deformation for the side sill of the bogie frame based on local-global method[J].China Welding,2007,16(4):11-16.

[7]衛(wèi)軍，汪東明，劉曉春,等.鋼梁橋焊接殘余應(yīng)力的計算機模擬[J].橋梁建設(shè),2015,45(4):94-99.

[8]錢海盛，張宇，鄧紅川，等.焊接工藝參數(shù)對Q235鋼焊接殘余應(yīng)力的影響[J].熱加工工藝,2015,44(5):169-171.

[9]鄧德安，清島祥一.焊接順序?qū)癜搴附託堄鄳?yīng)力分布的影響[J].焊接學報,2011,32(12):55-60.

[10]王憲，杜義，張京海.道間溫度對10CrNi5MoV鋼氣體保護焊接頭力學性能的影響[J].材料開發(fā)與應(yīng)力,2014,29(2):37-41.

[11]HU M J, LIU J H. Effects of zonal heat treatment on residual stresses and mechanical properties of electron beam welded TC4 alloy plates[J]. Transactions of Nonferrous Metals Society of China,2009,19:324-329.

[12]中華人民共和國住房與城鄉(xiāng)建設(shè)部.鋼結(jié)構(gòu)焊接規(guī)范：GB 50661-2011[S].北京：中國建筑出版社，2011.

[13]馮艷輝.基于ANSYS的焊接梁殘余應(yīng)力分析[D].西安：西安建筑科技大學，2008.

Research on the Influence of Welding Cooling Time on the Residual Stress of Tensile Anchor Tube

CHEN Yali1)HE Jia1)XIAO Xiang1)DONG Jie1)

(SchoolofTransportation,WuhanUniversityofTechnology,Wuhan430063,China)1)(SchoolofCivilEngineering,ChinaThreeGorgesUniversity,Yichang443000,China)2)

The welding cooling time of tensile anchor tube on the Cable-Stayed Bridge has great influence on the residual stress of the structure. Therefore, research on the cooling time is considered of great importance. The influence of cooling time on the residual stress is studied by means of both the numerical simulation and experimental measurement. Firstly, the residual stress of typical butt and T welds are measured by magnetic and blind hole method. Meanwhile, a three-dimensional model of welding seam is established by ANSYS finite element software. The test results verify the feasibility of the numerical simulation method. In addition, based on the numerical simulation method, variation of residual stress of tensile anchor tube on the Cable-Stayed Bridge under different cooling time is researched. Results show that with the growth of welding cooling time of tensile anchor tube, the distribution law of residual stress remains unchanged, but the peak tensile stress increases and the peak compressive stress decreases. In order to reduce the adverse effects of residual stress, measures of shortening the cooling time or preheating and heat preservation should be taken.

road engineering; tensile anchor tube; numerical simulation; welding residual stress; blind hole method; cooling time

2017-05-18

*甘肅省交通運輸廳科技項目資助(201674)

U448.27

10.3963/j.issn.2095-3844.2017.04.032

陳亞莉(1993—)：女，碩士生，主要研究領(lǐng)域為橋梁結(jié)構(gòu)耐久性設(shè)計與安全評估