嚴(yán)浩東,馬思群，于佳田，蒼 松

(1.大連交通大學(xué) 機(jī)車車輛工程學(xué)院，遼寧 大連 116028； 2.遼寧軌道交通職業(yè)學(xué)院，遼寧 沈陽(yáng) 110023)

0 引 言

塔式起重機(jī)作為特種設(shè)備之一，主要應(yīng)用在交通運(yùn)輸、傳統(tǒng)機(jī)械、核工業(yè)等領(lǐng)域，是一種典型的焊接結(jié)構(gòu)，其安全性一直備受行業(yè)的高度重視。近年來(lái)，隨著計(jì)算機(jī)技術(shù)的不斷發(fā)展，數(shù)值模擬技術(shù)在焊接過(guò)程中的力學(xué)分析、能量分析、組織性能分析方面均得到越來(lái)越多的應(yīng)用[1]。武漢理工大學(xué)的張凱[2]在熱彈塑性有限元法的基礎(chǔ)下，針對(duì)船體結(jié)構(gòu)典型T型接頭角焊縫進(jìn)行了焊接變形預(yù)測(cè)，研究了構(gòu)件平面尺寸、焊接速度、焊角尺寸、板材的厚度以及焊接順序的變化等因素對(duì)焊接變形的影響。上海交通大學(xué)的朱琦峰[3]將熱彈塑性有限元法進(jìn)行簡(jiǎn)化，分析了不同焊接順序下箱體模型的溫度場(chǎng)、應(yīng)力場(chǎng)以及焊接變形數(shù)值，對(duì)比不同方案下三者的變化規(guī)律。焊接產(chǎn)生的殘余應(yīng)力會(huì)影響整座塔吊的安全應(yīng)用，筆者基于熱彈塑性理論，采用雙橢球熱源模型，建立了二維網(wǎng)格模型，使用Hypermesh軟件和SYSWELD軟件對(duì)適合該雙臂塔式起重機(jī)焊接數(shù)值模擬的網(wǎng)格尺寸、熱影響區(qū)寬度進(jìn)行了確定，還對(duì)焊接電流對(duì)焊接件溫度場(chǎng)、熔池形狀與焊后殘余應(yīng)力的影響規(guī)律做了研究。

1 研究對(duì)象

塔尖是主要承力機(jī)構(gòu)，用來(lái)承受臂架和平衡臂拉繩傳來(lái)的上部荷載。文中以某型雙臂塔式起重機(jī)的塔尖為研究對(duì)象。如圖1所示，該塔尖由A字承力結(jié)構(gòu)、鋼架和外部吊架等構(gòu)件組焊而成。為了研究方便，首先對(duì)塔尖兩端起承力作用的A字結(jié)構(gòu)開展焊接仿真的相關(guān)研究。該A字結(jié)構(gòu)所用材料為Q345，由3 mm的鋼板、6 mm的中空?qǐng)A管和7 mm的中空外架通過(guò)焊縫1～6焊接而成。焊縫1～6全為角焊縫且為單道焊接，焊縫5和焊縫6為環(huán)焊縫，其結(jié)構(gòu)如圖2所示。

圖1 塔尖整體外貌圖 圖2 A字結(jié)構(gòu)焊縫示意圖

2 熱源模型選取

熱彈塑性法是焊接數(shù)值模擬的一種基礎(chǔ)方法。熱彈塑性有限元法焊接分為以下兩步：首先分析焊接傳熱，隨后進(jìn)行應(yīng)力變形分析[4]。在使用熱彈塑性法進(jìn)行數(shù)值模擬分析時(shí)，最為關(guān)鍵和基礎(chǔ)的是熱源模型的選取。常用的熱源模型有平面高斯熱源模型、雙橢球分布熱源模型、3D高斯圓錐形熱源模型等。相對(duì)于其它熱源模式，雙橢球熱源模型擁有更為復(fù)雜的熱流密度函數(shù)和較多的焊接定義參數(shù)，能夠較好地反應(yīng)熱源分布情況，一般適用于手工電弧焊、熔化極氣體保護(hù)焊、熔化極活性氣體保護(hù)焊等焊接方法。

在實(shí)際焊接過(guò)程中，金屬材料在電流作用下熔化形成熔池，熔池內(nèi)的能量因電弧移動(dòng)呈不對(duì)稱分布，存在前后差異。以焊接方向?yàn)檎较颍娀『蟛康哪芰看笥陔娀∏安康哪芰浚p橢球熱源模式對(duì)這一現(xiàn)象進(jìn)行了很好地模擬，其由前后兩個(gè)不同的1/4橢球體組成，如圖3所示。

圖3 雙橢球熱源模型

a、b、cb、cf均為熱源參數(shù)；a影響熔池熔寬；b影響熔深；cb和cf的確定則需要一定的經(jīng)驗(yàn)和一定的時(shí)間精力。ff、fr為前后橢球的熱量輸入分布系數(shù)，利用雙橢球式可以得出前后橢球的熱流分布為[5]：

x≥0

(1)

x<0

(2)

前半部分的熱輸入為：

(3)

對(duì)于后半部分同理為:

(4)

由于：

(5)

所以：

(6)

3 網(wǎng)格尺寸及熱影響區(qū)對(duì)仿真的影響

3.1 網(wǎng)格尺寸對(duì)仿真結(jié)果的影響

圖4所示為本次焊接數(shù)值模擬的技術(shù)路線，由圖可知焊接仿真的結(jié)果與幾何模型的還原度、熱源的類型、材料性能等多重因素相關(guān)。劃分的網(wǎng)格在整個(gè)過(guò)程中起著決定性的作用。網(wǎng)格尺寸的大小會(huì)影響仿真的精度，也會(huì)影響計(jì)算的迭代次數(shù)以及仿真的所用時(shí)間。根據(jù)已有研究發(fā)現(xiàn)，網(wǎng)格尺寸越小，仿真精度越高，同時(shí)仿真時(shí)間增加。

圖4 焊接數(shù)值模擬技術(shù)路線

為了使仿真結(jié)果更加準(zhǔn)確，需要研究網(wǎng)格大小對(duì)焊接數(shù)值模擬的影響。首先，使用Pro/E軟件對(duì)塔尖上的A字結(jié)構(gòu)進(jìn)行幾何建模。由于工藝設(shè)計(jì)的需要，此結(jié)構(gòu)中存在著很多的倒角，這些倒角的刪除與否對(duì)焊接仿真的結(jié)果不會(huì)有影響，但是會(huì)增加網(wǎng)格劃分的難度，也會(huì)在一定程度上影響網(wǎng)格劃分質(zhì)量。因此，在建模軟件中直接建立出了幾何清理后的中面模型。

其次，在HyperMesh中將建立完的中面模型劃分成兩種不同尺寸的網(wǎng)格，一種為5 mm的網(wǎng)格，如圖5所示，另一種為10 mm的網(wǎng)格，采用四節(jié)點(diǎn)四邊形單元為主，三節(jié)點(diǎn)三角形單元為輔的劃分方式。焊縫單元可以獨(dú)立于母材而被劃分成實(shí)體單元或者殼單元兩種形式，殼固耦合模型和純殼單元模型的仿真結(jié)果基本吻合。因此，在Visual-Weld和Visual-Mesh中對(duì)焊縫以及熱影響區(qū)進(jìn)行二維單元的劃分。經(jīng)過(guò)網(wǎng)格劃分，5 mm模型的單元數(shù)量為25 806,節(jié)點(diǎn)數(shù)為25 416；10 mm模型的單元數(shù)量為6 651，節(jié)點(diǎn)數(shù)為6 444。

圖5 5 mm網(wǎng)格離散圖

為了研究方便，對(duì)模型進(jìn)行了以下合理的簡(jiǎn)化：①熔池的液態(tài)金屬設(shè)定為不可壓縮性粘性流體，流動(dòng)方式為層間流動(dòng)[6]；②不對(duì)焊件進(jìn)行預(yù)熱處理，設(shè)置焊接初始溫度為20 ℃；③除了與溫度有關(guān)的參數(shù)以外，材料的其他參數(shù)都不隨溫度的變化而變化；④該塔吊采用的材料為Q345低碳鋼，熔化溫度為1 400～1 500 ℃，SYSWELD現(xiàn)有的材料庫(kù)中S355J2G3的性能與所需材料相似，直接采用該材料定義焊接母材以及焊縫單元。

最后，采用控制變量法進(jìn)行網(wǎng)格尺寸對(duì)焊接數(shù)值模擬的影響研究。在進(jìn)行焊接導(dǎo)向設(shè)置的時(shí)候，需要保證熱影響區(qū)的范圍、模型的夾持條件、冷卻環(huán)境、焊接電流與焊接速度等焊接參數(shù)需保持一致(焊接電壓為25 V，焊接電流為250 A，焊接速度為10 mm/s，焊接熱效率為0.75)，最終得到圖6中的兩種網(wǎng)格尺寸下的焊接仿真準(zhǔn)穩(wěn)態(tài)溫度場(chǎng)圖。由于焊縫周圍的金屬在1 500 ℃時(shí)已完全融化。因此，為了觀察方便，文中的溫度場(chǎng)上限全部調(diào)整為1 500 ℃。

從仿真結(jié)果可以看出，網(wǎng)格的尺寸會(huì)嚴(yán)重影響焊接仿真結(jié)果。觀察圖6(b)中的10 mm網(wǎng)格的準(zhǔn)穩(wěn)態(tài)溫度場(chǎng)發(fā)現(xiàn)，采用10 mm網(wǎng)格劃分的模型，最低溫度為-12.584 6 ℃，低于設(shè)置的室溫20 ℃；最高溫度為1 880.45 ℃，與所測(cè)的實(shí)際最高溫度相差較大。而且，觀察該網(wǎng)格尺寸下的整個(gè)焊接仿真過(guò)程發(fā)現(xiàn)，溫度負(fù)峰值始終會(huì)出現(xiàn)在熱源的前面，跟隨著熱源的移動(dòng)而動(dòng)，并且熔池形狀與采用的熱源類型所對(duì)應(yīng)形成的熔池不同。因此，10 mm的網(wǎng)格無(wú)法準(zhǔn)確地對(duì)該塔吊進(jìn)行焊接仿真。

圖6(a)所示為5 mm的溫度場(chǎng)圖，準(zhǔn)穩(wěn)態(tài)峰值溫度為1 538.97 ℃，結(jié)果與實(shí)測(cè)結(jié)果相符，誤差在標(biāo)準(zhǔn)范圍之內(nèi)。由此看來(lái)，選擇合理的單元尺寸尤為重要，在以后的焊接仿真研究工作中，需要根據(jù)模型的大小、仿真精度要求來(lái)選擇合適正確的網(wǎng)格尺寸。而且對(duì)此模型，10 mm單元和5 mm單元的計(jì)算所需時(shí)長(zhǎng)相差不大。因此，對(duì)于一定規(guī)格的模型，在綜合考慮數(shù)值模擬的同時(shí)，可以優(yōu)先考慮使用小網(wǎng)格，小網(wǎng)格單元質(zhì)量更高，仿真結(jié)果更為準(zhǔn)確。

圖6 兩種網(wǎng)格尺寸的焊接溫度場(chǎng)圖

通過(guò)觀察溫度場(chǎng)，還可以發(fā)現(xiàn)當(dāng)達(dá)到焊接準(zhǔn)穩(wěn)態(tài)時(shí)，電弧前方的溫度梯度要大于后方的溫度梯度。這是因?yàn)殡娀⊙刂附臃较蛞苿?dòng)時(shí)，對(duì)電弧后方的焊接件會(huì)有類似預(yù)熱的過(guò)程，所以電弧前方的焊接件相對(duì)于后方的加熱時(shí)間短，熱輸入量小，因此電弧前方的溫度梯度較大，其前后的溫度分布不對(duì)稱，并且峰值溫度隨著電弧的前移而前移。

3.2 熱影響區(qū)的寬度對(duì)仿真結(jié)果的影響

焊接熱影響區(qū)的性能對(duì)焊接件的使用壽命有著很大的影響。對(duì)于工業(yè)生產(chǎn)常用的鋼材，根據(jù)每個(gè)區(qū)域的峰值溫度，可以將熱影響區(qū)劃分為粗晶熱影響區(qū)、細(xì)晶熱影響區(qū)、臨界熱影響區(qū)和亞臨界熱影響區(qū)。熱影響區(qū)的寬度與焊接方法、焊接工藝、母材厚度等多種因素有關(guān)，該塔吊采用的電弧焊的熱影響區(qū)寬度一般為6.0～8.5 mm[7]。

為了研究熱影響區(qū)寬度對(duì)焊件溫度變化的影響。在焊接電流250 A，焊接電壓25 V等其他焊接參數(shù)統(tǒng)一的情況下，建立了3 mm、6 mm、8.5 mm、10 mm、15 mm五種寬度的熱影響區(qū)。觀察各熱影響區(qū)寬度下的焊接仿真過(guò)程，發(fā)現(xiàn)當(dāng)焊接時(shí)間為40 s時(shí)，焊接早已達(dá)到準(zhǔn)穩(wěn)態(tài)。因此，以該時(shí)刻的電弧中心到達(dá)的焊縫單元上的節(jié)點(diǎn)為觀測(cè)點(diǎn)，如圖7所示。

圖7 觀測(cè)點(diǎn)選取位置圖

在同一熱影響區(qū)寬度下，提取各個(gè)被選節(jié)點(diǎn)的溫度變化曲線，并求其平均值，進(jìn)行數(shù)據(jù)整理之后，得到圖8中的不同熱影響區(qū)寬度下的溫度變化曲線。

圖8 不同熱影響區(qū)寬度的溫度變化柱狀圖

從圖8可以發(fā)現(xiàn)在焊接數(shù)值模擬仿真時(shí)，熱影響區(qū)的選取寬度會(huì)對(duì)焊接溫度場(chǎng)產(chǎn)生一定的影響。熱影響區(qū)選的寬度不夠大，即達(dá)不到熱影響區(qū)最低理論值時(shí)，會(huì)使仿真溫度峰值過(guò)大，無(wú)法正確指導(dǎo)實(shí)際焊接。當(dāng)選取的熱影響區(qū)寬度大于或等于8.5 mm時(shí)，溫度場(chǎng)的峰值溫度不再發(fā)生變化，即當(dāng)熱影響區(qū)范圍大于焊件匹配的熱影響區(qū)寬度的理論值時(shí)，就不會(huì)對(duì)數(shù)值模擬的溫度場(chǎng)產(chǎn)生影響，同一位置的溫度變化曲線相同。并且，將熱影響區(qū)寬度為6 mm的溫度變化曲線與8.5 mm以上的對(duì)比，發(fā)現(xiàn)各時(shí)刻的溫度值相差在10 ℃以內(nèi)。這是因?yàn)闊嵊绊憛^(qū)寬度與母材厚度有關(guān)，工件越薄，熱影響區(qū)寬度越接近理論值下限。因此，在確保結(jié)果準(zhǔn)確性的基礎(chǔ)上，同時(shí)考慮建模的方便，對(duì)該A字結(jié)構(gòu)進(jìn)行焊接仿真時(shí)，焊接熱影響區(qū)統(tǒng)一選取為10 mm的寬度。在對(duì)轉(zhuǎn)臺(tái)、起重臂等其他雙臂塔式起重機(jī)部位進(jìn)行焊接仿真時(shí)，熱影響區(qū)的寬度也可以統(tǒng)一選取為10 mm，能夠在結(jié)果可靠的基礎(chǔ)上大幅度地降低離散模型的難度。

4 焊接電流對(duì)焊接的影響

焊接過(guò)程會(huì)受焊接電流、焊接電壓、焊接速度等多個(gè)參數(shù)的綜合影響。由于篇幅限制，文中先對(duì)焊接電流這一焊接參數(shù)開展研究，并在特定的焊接電壓、焊接速度的前提下，確定適合該塔尖焊接電流大小。

4.1 焊接電流對(duì)溫度場(chǎng)的影響

改變焊接電流會(huì)使焊接熱輸入量發(fā)生改變。如表1所列，在焊接電壓(25 V)、焊接速度(10 mm/s)、焊接熱效率(0.75)等其它焊接參數(shù)保持不變的前提下，確定了三組焊接方案，開展焊接電流對(duì)焊件焊縫周圍溫度變化影響規(guī)律的研究。同樣以圖7中的節(jié)點(diǎn)為觀察點(diǎn)，計(jì)算出每一焊接電流下的平均溫度變化值，得到圖9中的不同焊接電流下的溫度變化曲線。

表1 不同焊接電流對(duì)應(yīng)的焊接方案

圖9 不同焊接電流下的溫度變化曲線圖

從圖中可以發(fā)現(xiàn)，隨著焊接電流的不斷增大，電弧前后的溫度梯度變大，溫度場(chǎng)的最大溫度不斷升高。當(dāng)電流達(dá)到240 A時(shí)，焊接的穩(wěn)態(tài)溫度達(dá)到1 465.63 ℃，已達(dá)到焊件所用材料的融化溫度范圍。但為了保證焊縫周圍的金屬能夠完全融化，需繼續(xù)提高焊接電流。當(dāng)焊接電流達(dá)到250 A時(shí)，準(zhǔn)穩(wěn)態(tài)峰值溫度達(dá)到1 538.97 ℃，已超過(guò)碳鋼融化溫度的上限，確保焊件能夠完全融化。

4.2 焊接電流對(duì)熔池的影響

焊接熱輸入與焊接電流的大小有著密切的關(guān)系，其他參數(shù)不變，焊接電流一定時(shí)，焊接熱輸入即為固定值。把熔池看成一個(gè)由焊接熱輸入和熔池散熱引起的熱行為整體，則熔池的散熱量只與其本身的大小、溫度以及其周圍的溫度有關(guān)，且散熱的速度與熔池大小和溫度成正比，與周圍溫度成反比。在焊接起始階段，熔池較小且其內(nèi)部溫度較低，散熱量小于吸熱量，熔池隨著電弧的前移而長(zhǎng)大，散熱量也不斷增加，熔池內(nèi)部的熱量增加速度和熔池長(zhǎng)大速度變緩，直到熔池形狀達(dá)到一定的尺寸形狀，散熱量等于吸熱量，熔池形狀確定。熔池形狀變化如圖10所示。

圖10 不同焊接電流下的熔池形狀

根據(jù)圖10所示，穩(wěn)態(tài)下的熔池的寬度和長(zhǎng)度隨著電流的增加不斷增大。這是因?yàn)樵黾雍附与娏?，?huì)使焊接熱輸入增加，焊接件上的焊縫處的各點(diǎn)的熱輸入也相應(yīng)增加，導(dǎo)致該焊接件的溫度場(chǎng)分布范圍變大，熔池也逐漸增大。

4.3 焊接電流對(duì)殘余應(yīng)力的影響

焊接電流是影響焊接殘余應(yīng)力產(chǎn)生的重要因素，而焊后殘余應(yīng)力會(huì)引起焊接件發(fā)生疲勞斷裂、失穩(wěn)，是導(dǎo)致工件使用壽命減短的主要原因之一[8]。因此，在對(duì)塔尖該部位的焊接溫度場(chǎng)模擬的基礎(chǔ)上，還對(duì)該焊接件的應(yīng)力場(chǎng)進(jìn)行了有限元計(jì)算，得到如圖11所示的殘余應(yīng)力云圖。

圖11 不同焊接電流下的殘余應(yīng)力

從應(yīng)力云圖中可以看出焊接電流為250 A的平均應(yīng)力最大為315.03 MPa，260 A的平均應(yīng)力最大為317.21 MPa，大的焊縫殘余應(yīng)力發(fā)生在焊縫中心處，而且在其他焊接參數(shù)不變的情況下，焊后殘余應(yīng)力會(huì)隨著焊接電流的增大而增大。而當(dāng)焊接電流為250 A時(shí)，焊接準(zhǔn)穩(wěn)態(tài)峰值溫度就已達(dá)到Q345的熔點(diǎn)。

因此，在焊接速度為10 mm/s、焊接電壓25 V的實(shí)際焊接情況下，采用250 A的焊接電流較為合適。焊接完成后，可以通過(guò)整體熱處理、錘擊焊縫等方法降低焊接殘余應(yīng)力。

5 結(jié) 論

文中針對(duì)某型雙臂塔式起重機(jī)塔尖的焊接問(wèn)題，基于熱彈塑性理論，通過(guò)使用HyperMesh和SYSWELD軟件進(jìn)行了數(shù)值仿真，得出了適合本塔式起重機(jī)仿真分析的網(wǎng)格尺寸與熱影響區(qū)寬度，并在特定的焊接電壓和焊接速度的情況下，確定了最佳的焊接電流，為該雙臂塔式起重機(jī)的焊接仿真以及實(shí)際生產(chǎn)提供了指導(dǎo)。

(1) 對(duì)焊接數(shù)值模擬的工件進(jìn)行前處理時(shí)，要選取適合的單元尺寸和熱影響區(qū)寬度。5 mm是離散該塔尖較為理想的單元尺寸，10 mm的熱影響區(qū)寬度可以應(yīng)用到整個(gè)塔吊的焊接模擬。

(2) 在焊接電壓為25 V，焊接速度10 mm/s的情況下，250 A焊接電流的峰值溫度已達(dá)到材料熔點(diǎn)，殘余應(yīng)力也在材料允許值內(nèi)。因此，采取250 A的焊接電流對(duì)該塔尖進(jìn)行焊接較為合適。

(3) 當(dāng)其他焊接尺寸不變時(shí)，焊接熔池的尺寸與殘余應(yīng)力都會(huì)隨著焊接電流的增大而增大。

(4) 可使用SYSWELD對(duì)整體塔吊等其他大型焊接結(jié)構(gòu)進(jìn)行仿真模擬，能夠大大縮短產(chǎn)品的設(shè)計(jì)周期，從而降低生產(chǎn)成本。