何雄君, 宋 淵 ,肖 祥, 王進(jìn)軍

(1.武漢理工大學(xué) 交通學(xué)院，湖北 武漢 430063；2.中交二航局六公司，湖北 武漢 430014)

0 引 言

目前世界上已經(jīng)建成了很多斜拉橋，斜拉橋可以有各種各樣的形式，但主體都為3個(gè)部分：塔柱、斜拉索、主梁[1]。錨拉板作為連接主梁與斜拉索的裝置，在斜拉橋運(yùn)營(yíng)階段起著至關(guān)重要的作用。筆者以西固大橋?yàn)楸尘?，探究西固大橋錨拉板焊接殘余應(yīng)力與焊接熱輸入的關(guān)系。錨拉板由鋼材焊接而成，其主要的焊接形式為T型構(gòu)件多道焊焊接，因此，T型構(gòu)件多道焊焊接是錨拉板焊接中需要著重考慮的關(guān)鍵技術(shù)之一。焊接是鋼材之間結(jié)合時(shí)金屬材料迅速升溫隨即又迅速降溫的過程，焊接后不可避免地會(huì)在焊縫及焊接熱影響區(qū)產(chǎn)生焊接殘余應(yīng)力。焊接殘余應(yīng)力的存在會(huì)使錨拉板產(chǎn)生一定的變形，對(duì)后期安裝的精度產(chǎn)生影響；焊接殘余應(yīng)力還會(huì)影響錨拉板的承載能力與疲勞強(qiáng)度[2]。而錨拉板T型構(gòu)件多道焊中，焊接熱輸入是一個(gè)對(duì)焊接殘余應(yīng)力影響較大的因素。

目前，在錨拉板T型構(gòu)件多道焊中關(guān)于焊接熱輸入沒有一個(gè)明確的界定，基本是工人自身的工作經(jīng)驗(yàn)確定。而焊接熱輸入對(duì)焊接殘余應(yīng)力會(huì)產(chǎn)生很大的影響，因此，若能掌握焊接熱輸入對(duì)錨拉板T型構(gòu)件多道焊焊接殘余應(yīng)力的影響，對(duì)于斜拉橋的建設(shè)顯得尤為重要。

利用ANSYS軟件的數(shù)值模擬方法測(cè)量焊接殘余應(yīng)力具有成本低、勞動(dòng)量小的優(yōu)點(diǎn)，已有很多學(xué)者運(yùn)用該方法對(duì)焊接殘余應(yīng)力進(jìn)行了仿真模擬，但大部分是針對(duì)單道焊縫焊接殘余應(yīng)力的數(shù)值模擬[3]。然而實(shí)際錨拉板T型構(gòu)件是分道焊接而成，由于分道焊接時(shí)，后面焊層會(huì)對(duì)前面已經(jīng)施焊焊層產(chǎn)生熱影響，該影響又是不確定的，故多道焊縫比單道焊縫焊接殘余應(yīng)力分布要復(fù)雜得多。運(yùn)用ANSYS軟件揭示多道焊焊接而成的錨拉板T型構(gòu)件焊接殘余應(yīng)力與焊接熱輸入之間的關(guān)系，指導(dǎo)斜拉橋錨拉板焊接工藝，將具有很大的工程意義。

1 有限元模型

1.1 模型的建立

筆者以西固大橋?yàn)楸尘?，斜拉橋錨拉板的體積較大，運(yùn)用有限元分析軟件ANSYS模擬整個(gè)錨拉板的焊接對(duì)計(jì)算機(jī)的硬件要求很高，也需要很長(zhǎng)的計(jì)算時(shí)間[4]。因?yàn)殄^拉板的主要焊接形式為T型焊接，為了能較好的反應(yīng)錨拉板的焊接殘余應(yīng)力也避免對(duì)計(jì)算機(jī)硬件過高的要求，可以將錨拉板焊接的仿真分析轉(zhuǎn)化為T型構(gòu)件焊接的仿真分析。

T型構(gòu)件的尺寸為：水平板400 mm×400 mm×36 mm，豎向板400 mm×400 mm×40 mm，模型主要采用8節(jié)點(diǎn)6面體單元，有少量的6節(jié)點(diǎn)5面體單元。運(yùn)用ANSYS軟件中的網(wǎng)格劃分命令對(duì)T型焊接構(gòu)件進(jìn)行網(wǎng)格劃分，在不影響計(jì)算結(jié)果的前提下可以對(duì)網(wǎng)格劃分進(jìn)行優(yōu)化，同時(shí)可以減少計(jì)算時(shí)間。 T型構(gòu)件可以分為3個(gè)區(qū)域：焊縫區(qū)、靠近焊接區(qū)、遠(yuǎn)離焊縫區(qū)。由于熱源對(duì)3個(gè)區(qū)域的影響作用依次減弱，故3個(gè)區(qū)域的網(wǎng)格密度也相應(yīng)地減小。劃分后的模型共有6 189個(gè)節(jié)點(diǎn)，13 851個(gè)單元，如圖1。

圖1 典型T型構(gòu)件ANSYS模型Fig. 1 ANSYS model of typical T-type component

1.2 溫度場(chǎng)

焊接是金屬材料局部迅速升溫，隨即又迅速降溫的過程。在焊接過程中隨著熱源的移動(dòng)，焊接構(gòu)件的溫度會(huì)由于時(shí)間和空間的不同而發(fā)生劇烈的變化，構(gòu)件的某些部位還會(huì)有受熱屈曲的現(xiàn)象。焊接過程中金屬的材料特性會(huì)發(fā)生變化，因此焊接過程中溫度場(chǎng)的分析是典型的非線性瞬態(tài)熱傳導(dǎo)問題，有限元中三維瞬態(tài)熱傳導(dǎo)問題的控制方程為[5]

(1)

式中：λ為金屬的導(dǎo)熱系數(shù)；T為溫度；x,y,z分別為T型構(gòu)件中垂直于焊縫方向、沿焊縫的方向和豎向板的方向；Q為熱源大??；ρ為金屬材料的密度；c為比熱容；t為焊接至金屬冷卻的時(shí)間。

導(dǎo)熱過程中通常運(yùn)用以下3類邊界條件[6]。

第1類為溫度邊界條件，表達(dá)式為

Ts=Ts(x,y,z,t)

(2)

第2類為熱流邊界條件，表達(dá)式為

(3)

第3類為對(duì)流換熱邊界條件，表達(dá)式為

(4)

1.3 熱源的選擇

西固大橋錨拉板的焊接形式為手工電弧焊，故選用高斯分布的熱源來進(jìn)行模擬實(shí)際熱源，高斯熱源的分布函數(shù)為[7]

(5)

式中：R為電弧有效加熱半徑；r為焊件上任意點(diǎn)至電弧中心的距離；η為焊接熱效率；U為焊接時(shí)的電壓；I為焊接時(shí)的電流。

運(yùn)用ANSYS軟件模擬一條焊縫的5道焊，即完整的一條焊縫是分5道焊接而成的，采用生死單元法來實(shí)現(xiàn)該過程，模型建立時(shí)將T型構(gòu)件以及5道焊縫全部建立起來。然后賦予第1道焊縫焊接材料并施加熱源，其余道焊縫則不賦予材料，不起焊縫的作用；當(dāng)賦予第2道焊縫材料并施加熱源時(shí)則第3道及后面的焊縫也不起焊縫的作用。按照此方法依次完成對(duì)5道焊縫和熱源的施加，即完成了一條完整的焊縫，實(shí)現(xiàn)5道焊的焊接工藝。

2 有限元及實(shí)際測(cè)量結(jié)果分析

2.1 有限元分析

筆者以T型構(gòu)件焊接為例(每條焊縫為5道焊)，研究焊接熱輸入對(duì)多道焊焊接殘余應(yīng)力的影響。焊接中選用的材料為Q370qE橋梁鋼，在焊接中，母材的相關(guān)參數(shù)會(huì)由于溫度的變化發(fā)生一定的變化，運(yùn)用ANSYS模擬的母材參數(shù)變化如表1。

表1 母材參數(shù)Table 1 Base material parameter

多道焊的焊接順序如圖2，圖中①、②表示兩條焊縫，每條焊縫有5道，焊縫為從內(nèi)到外依次施焊。

圖2 焊接順序Fig. 2 Welding sequence

實(shí)際上，隨著焊接的進(jìn)行，金屬材料的應(yīng)力應(yīng)變有一個(gè)復(fù)雜的變化過程。筆者首先將焊絲中輸入的熱量轉(zhuǎn)換為金屬材料的溫度場(chǎng)，再將溫度場(chǎng)轉(zhuǎn)換為應(yīng)力場(chǎng)，最后利用在有限元模型中設(shè)置的節(jié)點(diǎn)提取該節(jié)點(diǎn)的焊接殘余應(yīng)力。

2.2 焊接殘余應(yīng)力的測(cè)量

對(duì)焊接殘余應(yīng)力的測(cè)量已存在多種測(cè)量方法[8]，總體可以分為機(jī)械測(cè)量方法和物理測(cè)量方法。機(jī)械測(cè)量方法主要有截條法、剝層法、盲孔法等。機(jī)械測(cè)量方法的原理是對(duì)結(jié)構(gòu)進(jìn)行破壞性試驗(yàn)，結(jié)構(gòu)的應(yīng)變會(huì)有相應(yīng)的變化，通過檢測(cè)這些變化即可求得殘余應(yīng)力的大小。機(jī)械測(cè)量方法的優(yōu)點(diǎn)是測(cè)量結(jié)果較為精確，缺點(diǎn)是會(huì)對(duì)構(gòu)件產(chǎn)生一定的損傷。物理測(cè)量方法[9]主要有X射線法、超聲波法、磁測(cè)法等，物理測(cè)量方法的原理是結(jié)構(gòu)受力會(huì)影響結(jié)構(gòu)的某些物理性能。而這些物理性能的改變可以由一定的參數(shù)變化反應(yīng)出來，通過測(cè)量這些參數(shù)變化經(jīng)過換算就可得到結(jié)構(gòu)的受力情況[10]。物理測(cè)量方法的優(yōu)點(diǎn)是不會(huì)對(duì)構(gòu)件有損害，但缺點(diǎn)是測(cè)量結(jié)果與構(gòu)件實(shí)際的受力狀況會(huì)有一定的偏差。

為保證試驗(yàn)的質(zhì)量，筆者運(yùn)用對(duì)被測(cè)構(gòu)件損傷較小且測(cè)量精度較高的盲孔法[11]來測(cè)量構(gòu)件的焊接殘余應(yīng)力。盲孔法的測(cè)量原理是在構(gòu)件的表面粘貼應(yīng)變花，然后在相應(yīng)的位置鉆取一定孔徑的圓孔，此時(shí)構(gòu)件內(nèi)部的殘余應(yīng)力會(huì)因?yàn)殂@孔而發(fā)生變化進(jìn)而有相應(yīng)應(yīng)變的變化。應(yīng)變花會(huì)識(shí)別應(yīng)變的變化，有了應(yīng)變，利用彈塑性力學(xué)等相關(guān)理論即可算得該殘余應(yīng)力的大小。實(shí)際中的焊接構(gòu)件及測(cè)量中測(cè)點(diǎn)布置如圖3。

圖3 焊接實(shí)物及測(cè)點(diǎn)布置Fig. 3 Welding object and measuring points layout

2.3 實(shí)測(cè)結(jié)果與軟件計(jì)算結(jié)果對(duì)比

實(shí)際焊接中運(yùn)用的相關(guān)工藝參數(shù)如表2。

表2 焊接參數(shù)Table 2 Welding parameters

采用實(shí)際中的焊接工藝，運(yùn)用盲孔法測(cè)得的沿焊縫方向的橫、縱殘余應(yīng)力與運(yùn)用ANSYS模擬實(shí)際焊接中調(diào)取的沿焊縫方向的橫、縱向應(yīng)力的比較如圖4。

圖4 焊接殘余應(yīng)力實(shí)測(cè)與數(shù)值模擬對(duì)比Fig. 4 Comparison of welding residual stress measurement and numerical simulation

從圖4中可以得到，用軟件計(jì)算的結(jié)果與盲孔法實(shí)測(cè)的結(jié)果雖然有一定的偏差，但總體比較吻合，焊接殘余應(yīng)力的趨勢(shì)一致，從而說明了運(yùn)用ANSYS軟件模擬該焊接工藝進(jìn)行殘余應(yīng)力檢測(cè)的可行性。圖4(a)顯示，沿著焊縫方向的橫向殘余應(yīng)力主要為壓應(yīng)力，拉應(yīng)力的范圍比較小，橫向殘余應(yīng)力基本以焊縫中間成對(duì)稱布置，拉應(yīng)力的最大值出現(xiàn)在距焊縫起點(diǎn)約50 mm處，拉應(yīng)力的大小約為4 MPa；圖4(b)顯示，沿焊縫方向的縱向殘余應(yīng)力在起弧、收弧的地方表現(xiàn)為壓應(yīng)力，其余地方表現(xiàn)為拉應(yīng)力，拉應(yīng)力的最大值約在焊縫中間出現(xiàn)，大約為135 MPa，縱向殘余應(yīng)力也基本以焊縫中間呈對(duì)稱布置。圖4的結(jié)果還可以得到，焊接材料的縱向殘余應(yīng)力要大于橫向殘余應(yīng)力。

3 焊接熱輸入對(duì)焊接殘余應(yīng)力影響

第2節(jié)中實(shí)測(cè)結(jié)果與軟件分析結(jié)果比較吻合，說明用ANSYS軟件模擬該焊接工藝進(jìn)行殘余應(yīng)力檢測(cè)的可行性。筆者將實(shí)際中的焊接熱輸入做上下調(diào)動(dòng)，即形成了其余參數(shù)不變，焊接熱輸入分別為1 226、1 716 、2 860 J/mm的3組試驗(yàn)，3組試驗(yàn)編號(hào)為Ⅰ、Ⅱ、Ⅲ，焊縫熱輸入分組如表3。

表3 焊縫熱輸入分組Table 3 Group of welding heat input

運(yùn)用ANSYS軟件對(duì)上述3組試驗(yàn)進(jìn)行仿真分析，下面提取3組試驗(yàn)沿焊縫方向的橫、縱向殘余應(yīng)力的結(jié)果進(jìn)行對(duì)比分析來說明焊接熱輸入對(duì)焊接殘余應(yīng)力的影響。運(yùn)用ANSYS模擬的3組不同熱輸入的應(yīng)力云圖如圖5，沿焊縫方向的橫、縱向殘余應(yīng)力的對(duì)比如圖6。

圖5 不同焊接熱輸入的殘余應(yīng)力云圖Fig. 5 Residual stress nephogram with different welding heat input

圖6 不同焊接熱輸入下數(shù)值模擬結(jié)果對(duì)比Fig. 6 Comparison of numerical simulation results with different welding heat input

由圖6(a)可知，當(dāng)焊接熱輸入增加時(shí)，沿焊縫方向的橫向殘余應(yīng)力中的拉應(yīng)力區(qū)域會(huì)變小，橫向殘余應(yīng)力有向壓應(yīng)力轉(zhuǎn)化的趨勢(shì)。在焊縫的起、收弧處的橫向殘余應(yīng)力變化很小，原因可能是在焊縫起、收弧處焊接材料有自由面，熱量可以與外界交換，故其受焊接熱輸入的影響不明顯。當(dāng)焊接熱輸入增加時(shí)，橫向殘余應(yīng)力有減小的趨勢(shì)，總體來說，焊接熱輸入對(duì)沿焊縫方向的橫向殘余應(yīng)力影響較小。由圖6(b)可知，縱向焊接殘余應(yīng)力峰值基本處于焊縫的中間，隨著焊接熱輸入的增加，縱向焊接殘余應(yīng)力的峰值在增大。當(dāng)焊接熱輸入由1 226 J/mm增加到1 716 J/mm時(shí)，縱向焊接殘余應(yīng)力最大值增加了約22 MPa；焊接熱輸入由1 716 J/mm增加到2 860 J/mm時(shí)，縱向焊接殘余應(yīng)力最大值增加了約8 MPa，故縱向焊接殘余隨著焊接熱輸入的增加會(huì)增大，但應(yīng)力增大的幅度會(huì)變小。

4 結(jié) 論

通過對(duì)T型構(gòu)件焊接情況的數(shù)值模擬與實(shí)際盲孔測(cè)試法的比較以及對(duì)不同焊接熱輸入下T型構(gòu)件多道焊的焊接情況進(jìn)行數(shù)值模擬，經(jīng)過比較、分析，可得到如下的結(jié)論：

1)運(yùn)用ANSYS軟件對(duì)焊接構(gòu)件進(jìn)行數(shù)值模擬得到的焊接殘余應(yīng)力值與測(cè)量精度較高的盲孔法測(cè)量殘余應(yīng)力值相符，說明了運(yùn)用此軟件模擬該焊接工藝的可行性，這可以節(jié)省很大的財(cái)力與物力。

2)焊接構(gòu)件起、收弧與外界接觸處的有效殘余應(yīng)力小于其余區(qū)域的有效殘余應(yīng)力，說明焊接的邊界對(duì)焊接的結(jié)果會(huì)有一定的影響。

3)焊接熱輸入對(duì)焊縫方向橫向殘余應(yīng)力的影響較小，對(duì)縱向殘余應(yīng)力的影響較大；隨著焊接熱輸入的增加縱向殘余應(yīng)力會(huì)增大，橫向殘余應(yīng)力會(huì)減小，但應(yīng)力變化的幅度會(huì)減小。

參考文獻(xiàn)(References)：

[1] 陳明憲.斜拉橋的發(fā)展與展望[J].中外公路,2006,26(4):76-86.

CHEN Mingxian. The development and prospect of cable-stayed bridge [J].JournalofChina&ForeignHighway, 2006, 26(4): 76-86.

[2] 劉小渝,劉秀偉. 斜拉橋鋼錨拉板區(qū)域焊接應(yīng)力消除試驗(yàn)研究[J].公路交通科技,2008(11):82-86,91.

LIU Xiaoyu, LIU Xiuwei. Experimental study on releasing welding residual stress in steel tensile anchor plate area of cable-stayed bridges [J].JournalofHighwayandTransportationResearchandDevelopment, 2008(11): 82-86, 91.

[3] 趙東升,劉玉君,孫敏科,等.碳鋼與304不銹鋼焊接殘余應(yīng)力的計(jì)算[J].焊接學(xué)報(bào),2012,33(1):93-95,117-118.

ZHAO Dongsheng, LIU Yujun, SUN Minke, et al. Calculation welding of residual stress of carbon steel and 304 stainless steel [J].TransactionsoftheChinaWeldingInstitution, 2012, 33(1): 93-95, 117-118.

[4] 馬思群,袁永文,馮良波,等.焊接速度對(duì)鋁合金多道焊焊接殘余應(yīng)力影響研究[J].鐵道學(xué)報(bào),2014,36(1):16-21.

MA Siqun, YUAN Yongwen, FENG Liangbo, et al. Research on effect of welding speed on aluminum alloy multi-pass welding residual stress [J].JournaloftheChinaRailwaySociety, 2014, 36(1): 16-21.

[5] 汪建華.焊接數(shù)值模擬技術(shù)及其應(yīng)用[M]. 上海:上海交通大學(xué)出版社,2003: 16-20.

WANG Jianhua.WeldingNumericalSimulationTechnologyandItsApplication[M]. Shanghai: Profile of Shanghai Jiaotong University Press, 2003: 16-20.

[6] 趙秋,吳沖.U肋加勁板焊接殘余應(yīng)力數(shù)值模擬分析[J].工程力學(xué),2012,29(8):262-268.

ZHAO Qiu, WU Chong. Numerical analysis of welding residual stress of U-RIB stiffened plate[J].EngineeringMechanics, 2012, 29(8): 262-268.

[7] 肖馮,米彩盈.T型角接頭焊接熱源模型研究[J].電焊機(jī),2010,40(6): 41-45.

XIAO Feng, MI Caiying. Study on heat source model for T-joint fillet weld [J].ElectricWeldingMachine, 2010, 40(6): 41-45.

[8] 張根樹,張穎.常用焊接殘余應(yīng)力測(cè)量方法的探究[J].科技信息,2012,32: 637.

ZHANG Genshu, ZHANG Ying. Research on ordinary measure method of welding residual stress[J].Science&TechnologyInformation, 2012, 32: 637.

[9] 沈軍,林波,遲永剛,等.殘余應(yīng)力物理法測(cè)量技術(shù)研究狀況[J].材料導(dǎo)報(bào),2012,26(1):120-125.

SHEN Jun, LIN Bo, CHI Yonggang, et al. Research status of residual stress physical method measurement techniques[J].MaterialsReview, 2012, 26(1): 120-125.

[10] 王慶明,孫淵.殘余應(yīng)力測(cè)試技術(shù)的進(jìn)展與動(dòng)向[J].機(jī)電工程,2011,28(1):11-15,41.

WANG Qingming, SUN Yuan. Research development on the test methods of residual stress[J].Mechanical&ElectricalEngineeringMagazine, 2011, 28(1): 11-15, 41.

[11] 王超逸,呂丹,蘇洪英,等.測(cè)定鋼板殘余應(yīng)力的改進(jìn)盲孔法[J].理化檢驗(yàn)(物理分冊(cè)),2015,51(11):781-783

WANG Chaoyi, LV Dan, SU Hongying, et al. Improved blind-hole method of measuring residual stress of steel plates[J].PhysicalTestingandChemicalAnalysisPartA(PhysicalTesting), 2015, 51(11): 781-783.