向祖權(quán)　劉　彬　杜洪文　劉春雨

(武漢理工大學(xué)高性能船舶技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室　交通學(xué)院　武漢　430063)

?

翼板厚度差異對T形接頭焊接變形的影響研究

向祖權(quán)劉彬杜洪文劉春雨

(武漢理工大學(xué)高性能船舶技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室交通學(xué)院武漢430063)

摘要：T形接頭焊接結(jié)構(gòu)形式廣泛應(yīng)用于船舶建造和修造的各個工序中，針對翼板厚度差異引起的T形接頭焊接溫度場、橫向收縮量、縱向收縮量以及角變形的研究規(guī)律問題，利用SYSWELD有限元分析軟件，對其焊接過程進(jìn)行數(shù)值模擬仿真，研究得到隨著翼板厚度的增加，焊接過程中翼板上下表面的峰值溫度差值逐漸增大，焊接完成后翼板的橫向收縮量、縱向收縮量以及角變形逐漸減小，這與實(shí)際生產(chǎn)規(guī)律相符，為造船及修船過程中由于翼板厚度差異引起的T形接頭焊接變形的預(yù)測，以及精度控制提供了前提條件.

關(guān)鍵詞：T形接頭；翼板厚度；橫向收縮量；縱向收縮量；角變形

0引言

在船舶建造過程中，從各種零部件的裝焊到分段的裝配，以及總段的合攏過程都需要大量的T形接頭焊接結(jié)構(gòu)件，相對于對接焊焊接工序來說，T形接頭焊接工序過程中的溫度場以及變形場要復(fù)雜很多，對于船廠在焊接過程中的精度控制要求也更高.目前針對國內(nèi)眾多中小型船廠在進(jìn)行T形接頭焊接過程中的精度控制問題，比如在不同翼板厚度下由于焊接過程中T形焊接構(gòu)件溫度的急劇變化造成的板材收縮變形[1]以及角變形降低了船舶建造過程中的焊接質(zhì)量，為了能夠控制這種情況，研究和探索由于翼板厚度的差異引起的T形接頭焊接溫度場、橫向收縮量、縱向收縮量以及角變形的規(guī)律，對于提高造船精度有著十分重要的意義.

早期的對于T形接頭焊接變形的研究，國內(nèi)外的許多專業(yè)人員做了大量研究工作[2-6]，對于T形接頭焊接變形的模擬研究工作出現(xiàn)了解析法、固有應(yīng)變法、熱-彈-塑性有限元法、基于人工神經(jīng)網(wǎng)絡(luò)的焊接變形分析法等方法.基于SYSWELD有限元分析軟件，對不同翼板厚度下T形接頭的焊接變形進(jìn)行了數(shù)值模擬仿真研究，總結(jié)了由于翼板厚度的差異引起的T形接頭焊接溫度場、橫向收縮變形、縱向收縮變形以及角變形的影響規(guī)律，從而為有效地預(yù)測和控制T形接頭的焊接變形以及提高造船精度提供了依據(jù).

1焊接物理模型的建立

1.1材料特性與幾何模型

母材所選用的材料為Q345低合金高強(qiáng)鋼，材料的主要化學(xué)成分見表1[7-8].T形構(gòu)件的尺寸：翼板為200 mm×100 mm×8 mm，腹板為200 mm×100 mm×8 mm.關(guān)于T形接頭幾何模型的創(chuàng)建，以及網(wǎng)格劃分過程是在Visual-mesh中實(shí)現(xiàn)的.

表1　Q345低合金高強(qiáng)鋼材料主要化學(xué)成分質(zhì)量分?jǐn)?shù)/%

主要研究T形接頭翼板厚度t在8，10，12，14 mm下焊接溫度場,以及變形場的不同情況，T形構(gòu)建模型見圖1.

圖1　T形構(gòu)件模型

1.2有限單元模型與裝卡條件

有限單元模型網(wǎng)格的劃分見圖2，采用SYSWELD中的8節(jié)點(diǎn)六面體單元進(jìn)行網(wǎng)格劃分，而且為了兼顧有限單元模型計算精度和計算速度之間的平衡在靠近焊縫區(qū)域的網(wǎng)格劃分比較細(xì)密，遠(yuǎn)離焊縫區(qū)域的網(wǎng)格劃分比較稀疏，劃分后的有限單元模型共有18 201個節(jié)點(diǎn)，24 250個單元(這對應(yīng)于翼板厚度為8 mm時，其他翼板厚度下有限單元模型的節(jié)點(diǎn)數(shù)與單元數(shù)目做適當(dāng)修改).

在實(shí)際焊接過程中，T形接頭焊接構(gòu)件是沒有外部夾持力即沒有位移和載荷約束，為了防止有限單元模型產(chǎn)生剛體位移和轉(zhuǎn)動，采用3支點(diǎn)自由約束條件[9]，即在節(jié)點(diǎn)A處施加x,y,z方向的剛性約束，在節(jié)點(diǎn)B施加y,z方向上的剛性約束，在節(jié)點(diǎn)C施加z方向上的剛性約束，裝卡條件見圖2.

圖2　有限單元模型及裝卡條件

1.3焊接熱源模型

整個焊接過程采用半自動焊，由于在焊接過程中熱流密度分布的不均勻性導(dǎo)致后方的加熱區(qū)域比前方的加熱區(qū)域大，使得整個加熱區(qū)域?yàn)殡p半橢球體并且前后2個半橢球體是不相同的，因此焊接熱源模型選用,見圖3的雙橢球體積分布熱源.

圖3　雙橢球體積熱源

圖4　校核完成的熱源模型橫截面圖

前后半橢球的熱流密度分布由式(1)、式(2)確定.

(1)

(2)

式中：af為橢球的前半軸長度；ar為后半軸長度；b為熔寬；c為熔深；q為熱流密度；Q為熱輸入功率，Q=ηUI.采用的焊接方式為半自動焊，填充的角焊縫為等腰直角三角形，雙面單層焊接，保護(hù)氣體為CO2，氣體流量大小為14 L/min，電流180 A，電壓26 V，采取的焊接速度為5 mm/s，焊接的熱效率η取0.8，經(jīng)過SYSWELD校核得到的T形接頭焊接熱源模型見圖4.

在進(jìn)行T形接頭焊接時，先焊接第一道焊縫，見圖2，待焊后溫度冷卻至100 ℃左右之后再進(jìn)行第二道焊縫的焊接.在實(shí)際裝焊的過程中，焊縫處的填充焊絲須滿足與待焊母材相容且具有相似的熱物理性能的條件，因此在SYSWELD模擬計算過程中假定焊縫處的填充焊絲與母材為同一種材料，并且利用熱-彈-塑性非線性的傳熱方程來描述電弧熱在被焊試件內(nèi)部的傳導(dǎo)方程[10]，即

(3)

式中：ρ為材料的密度；c為材料的比熱容；T=T(x,y,z,t)為節(jié)點(diǎn)(x,y,z)在t時刻的溫度；λx，λy，λz為材料沿x，y，z方向上的熱傳導(dǎo)系數(shù)；q為材料內(nèi)部的熱流密度.在進(jìn)行有限元模型計算時，利用牛頓法則和波爾茨曼定律分別考慮T形焊接工件與外部環(huán)境對流和輻射的邊界條件[11].

2數(shù)值模擬結(jié)果及分析

2.1焊接溫度場結(jié)果及討論

2.1.1翼緣板上表面與下表面溫度

圖5為在翼板厚度在8，10，12，14 mm下焊接過程中垂直于焊縫中央截面的翼緣板上表面與下表面節(jié)點(diǎn)(如圖1的A，B節(jié)點(diǎn))的溫度變化曲線，可以看到A，B節(jié)點(diǎn)的溫度變化曲線都呈現(xiàn)一種升溫速率遠(yuǎn)遠(yuǎn)大于冷卻速率的規(guī)律，主要原因是在冷卻過程中節(jié)點(diǎn)處會受到后面熔池的再熱作用減緩了其冷卻速率，還可以看到A，B節(jié)點(diǎn)在350 s左右時溫度會再次上升，主要是節(jié)點(diǎn)開始受到第二道焊縫焊接過程的影響.

表2為在節(jié)點(diǎn)A，B處的峰值溫度，以及其峰值溫度差值.根據(jù)圖5及表2知，隨著T形構(gòu)件翼板厚度的增加，焊接過程中翼緣板上下表面節(jié)點(diǎn)的溫度差值逐漸增大，這是因?yàn)樵谀M焊接過程中所選用的熱源是相同的，即熱輸入相同的前提下由于翼板厚度的增加導(dǎo)致在傳熱過程中翼緣板上下表面的溫度差加大.

圖5　節(jié)點(diǎn)A，B處的溫度變化曲線

表2　節(jié)點(diǎn)A，B處的峰值溫度，以及溫度差值

2.2焊接變形結(jié)果及討論

2.2.1焊接整體變形

焊接完成后翼板厚度在8,10,12,14 mm下的T形焊接構(gòu)件的整體變形云圖見圖6，為了方便觀察將變形結(jié)果放大5倍.T形焊接構(gòu)件的整體變形是由其橫向收縮變形、縱向收縮變形,以及角變形等綜合影響的結(jié)果.

2.2.2焊接收縮變形

對于T形焊接構(gòu)件翼板的橫向收縮變形以及縱向收縮變形，取載荷步在3 600 s時翼板厚度在8,10,12,14 mm下的焊接變形數(shù)值，見圖7.

圖6　不同翼板厚度下的T形焊接構(gòu)件的整體變形云圖

圖7　收縮量變化曲線

由圖7可知，隨著T形焊接構(gòu)件翼板厚度的增加，焊后翼板的橫向收縮變形以及縱向收縮變形逐漸減小.

2.2.3焊接角變形

其中T形焊接構(gòu)件角變形量以z方向上的節(jié)點(diǎn)位移變化來衡量，圖8表示按照圖1所示的a～b路徑各節(jié)點(diǎn)焊接完成冷卻至室溫后角變形的分布特征.對于T形焊接構(gòu)件的角變形應(yīng)該從兩方面進(jìn)行考慮:一方面是隨著翼板厚度的增加使得翼緣板上下表面的溫差加大(由圖5及表2得到)，從而會加大角變形的產(chǎn)生；另一方面，由于翼板厚度增加使得板材的剛度增大，會阻礙其角變形的產(chǎn)生.結(jié)果表明，隨著T形焊接構(gòu)件翼板厚度的增加，翼板的角變形逐漸減小，說明翼板厚度即剛度增加對角變形阻礙的影響作用大于溫度差值驅(qū)動角變形增加的影響作用.

圖8　a～b路徑上各節(jié)點(diǎn)的角變形變化曲線

3結(jié) 束 語

利用SYSWELD焊接模擬軟件對在船體結(jié)構(gòu)不同翼板厚度下T形焊接構(gòu)件焊接變形數(shù)值模擬研究的結(jié)果發(fā)現(xiàn)：在溫度場中，隨著T形焊接構(gòu)件翼板厚度的增加，翼板邊緣上下表面的溫度差值逐漸增大；變形場中，隨著T形焊接構(gòu)件翼板厚度的增加，翼板的橫向收縮量、縱向收縮量以及角變形逐漸減小.

根據(jù)數(shù)值模擬結(jié)果發(fā)現(xiàn)T形焊接構(gòu)件翼板的角變形是隨著翼板厚度的增加逐漸減小，說明翼板剛度增加對于角變形的抑制力在一定程度上大于翼緣板上下表面溫度差值對于角變形的驅(qū)動力，在整體上呈現(xiàn)出一種角變形隨著翼板厚度的增加逐漸減小的規(guī)律，與實(shí)際焊接規(guī)律相符.

參 考 文 獻(xiàn)

[1]向祖權(quán)，劉彬，郭純軒，等.板材焊縫間隙差異對溫度場及應(yīng)力場的影響[J].船舶工程,2015,37(10):79-83.

[2]DENG D, OGAWA K, KIYOSHIMA S, et al. Prediction of residual stresses in a dissimilar metal welded pipe with considering cladding, buttering and post weld heat treatment [J]. Computational Materials Science, 2009,47:398-408．

[3]DENG D, LIANG W, MURAKAWA H. Determination of welding deforma tion in fillet-welded joint by means of numerical simulation and comparison with experimental measurements [J]. Journal of Materials Processing Technology, 2007,183:219-225．

[4]遲鳴聲.鎂合金復(fù)合熱源焊接工藝與數(shù)值模擬研究[D].大連:大連理工大學(xué),2005.

[5]徐琳,嚴(yán)仁軍.T形焊接接頭殘余應(yīng)力與變形的三維數(shù)值模擬[J].江蘇船舶,2007(1):5-8.

[6]SPRARAGEN W, ETTINGER W G. Shringkage Distortion in Welding [J]. Welding Journal，2014,23(11):544-559.

[7]李振江.基于SYSWELD的焊接接頭溫度場和殘余應(yīng)力場研究[D].北京:北京交通大學(xué),2010.

[8]馬慶芳.實(shí)用熱物理性質(zhì)手冊[M].北京:中國農(nóng)業(yè)機(jī)械出版社,1986.

[9]黃道業(yè).高強(qiáng)度鋼T形結(jié)構(gòu)焊接變形固有應(yīng)變數(shù)值模擬[J]. Hot Working Technology, 2013,42(11):218-224.

[10]鄧德安，童彥剛，周中玉.薄壁低碳鋼管焊接變形的數(shù)值模擬[J].焊接學(xué)報,2011(2):81-85.

[11]周一俊，鄧德安，馮可，等.低碳鋼薄板單道堆焊焊接變形的數(shù)值模擬[J].焊接學(xué)報,2013(12):101-105.

Research on T-joint Welding Deformation in Different Flange Thickness

XIANG ZuquanLIU BinDU HongwenLIU Chunyu

(KeyLaboratoryofHighPerformanceShipTechnologyandSchoolofTransportation,WuhanUniversityofTechnology,Wuhan430063,China)

Abstract：T-joint welding structure is widely used in ship construction and repairing, while there is almost no studies on the effect of different flange thickness on temperature field, transverse shrinkage, longitudinal shrinkage and angular deformation. This paper applies the software SYSWELD to simulate the welding process of T-joint. With the increase of flange thickness, the peak temperature difference between the flange up and flange down is found to increase gradually in the process of welding, while the transverse shrinkage, longitudinal shrinkage and angular distortion are reduced gradually. This is consistent with the actual process and provides the premise for the forecast and control of T-joint deformation resulting from the difference of flange thickness in the process of ship construction and repairing.

Key words：T-joint; flange thickness; transverse shrinkage; longitudinal shrinkage; angular deformation

doi:10.3963/j.issn.2095-3844.2016.02.012

中圖法分類號：U671.3

收稿日期：2016-02-10

向祖權(quán)(1973- )：男，博士，副教授，主要研究領(lǐng)域?yàn)榇跋冗M(jìn)制造技術(shù)