胡美娟，何浩華，齊麗華，李為衛(wèi)，陳宏遠(yuǎn)

(1.中國(guó)石油集團(tuán)石油管工程技術(shù)研究院，石油管材及裝備材料服役行為與結(jié)構(gòu)安全國(guó)家重點(diǎn)實(shí)驗(yàn)室 陜西 西安 710077；2.長(zhǎng)慶油田物資供應(yīng)處(物資管理部) 陜西 西安 710016)

0 引 言

由于焊接時(shí)快速的加熱、冷卻及溫度分布的不均勻性，焊接接頭表現(xiàn)為一個(gè)連續(xù)變化并具有非常陡峭組織梯度的區(qū)域。根據(jù)焊接方法和工藝的不同，焊接熱影響區(qū)的大小從零點(diǎn)幾毫米到幾十毫米。在這個(gè)寬度范圍內(nèi)包括有幾個(gè)組織和性能不同的特定溫度區(qū)。利用實(shí)際的焊接接頭來(lái)研究測(cè)定特定區(qū)域的力學(xué)性能和焊接特性是十分困難的。焊接熱模擬利用熱-力學(xué)模擬試驗(yàn)機(jī)在試樣上重現(xiàn)焊接熱影響區(qū)在焊接時(shí)的熱、應(yīng)力及應(yīng)變循環(huán)，使得試樣在較大的尺寸范圍內(nèi)獲得焊接熱影響區(qū)某一特定溫度區(qū)的均勻組織，進(jìn)而實(shí)現(xiàn)對(duì)其組織及性能的研究[1-2]。

焊接熱模擬試驗(yàn)的關(guān)鍵就是制定符合實(shí)際的焊接熱-力過(guò)程模擬方案。焊接熱模擬方案的制定就是模擬區(qū)域在焊接時(shí)所經(jīng)歷的溫度-時(shí)間，應(yīng)力-時(shí)間和應(yīng)變-時(shí)間三種類型曲線的確定。其中，溫度隨時(shí)間的變化即焊接熱循環(huán)曲線的制定不僅是焊接領(lǐng)域物理模擬的基礎(chǔ)和主要內(nèi)容，也是預(yù)測(cè)和評(píng)定應(yīng)力應(yīng)變條件下焊接接頭性能的基礎(chǔ)。美國(guó)DSI(Dynamic Systems Inc.)公司生產(chǎn)的Gleeble 3系列熱-力學(xué)模擬試驗(yàn)機(jī)中配備的焊接熱模擬軟件以實(shí)測(cè)的和經(jīng)典的數(shù)學(xué)解析表達(dá)式為基礎(chǔ)，包含了F(s,d)表格；Hannerz；Rykalin-2D；Rykalin-3D；Rosenthal和指數(shù)冷卻(Exponential cooling)6種不同的數(shù)學(xué)模型[3-4]。本文就針對(duì)這些模型對(duì)焊接熱循環(huán)的計(jì)算方法和應(yīng)用等問(wèn)題進(jìn)行深入的探討。

1 計(jì)算焊接熱循環(huán)的數(shù)學(xué)模型

1.1 F(s,d)表格模型

基于F(s,d)表格模型的熱循環(huán)曲線是通過(guò)對(duì)前人實(shí)際測(cè)得的不同材料以時(shí)間和距離為函數(shù)的數(shù)據(jù)表進(jìn)行插值所得。其公式如式(1)所示：

(1)

式中：T0為預(yù)熱溫度或?qū)娱g溫度， ℃；s為時(shí)間，s；d為距焊縫中心線的距離，cm；InputFactor為根據(jù)預(yù)熱溫度所得輸入系數(shù)。

軟件目前提供的焊接材料有25.4 mm厚鋼板、25.4～12.7mm厚不銹鋼、12.7 mm厚鋼板、12.7 mm厚鈦板、6.35 mm 厚鈦 a55板和6.35 mm 厚鋯2板。

在焊接熱循環(huán)計(jì)算時(shí)，F(xiàn)(s,d)表格模型需要輸入峰值溫度Tmax、預(yù)熱溫度T0和焊接線能量E，并選擇焊接材料。如果輸入的數(shù)值不在F(s,d)表格的范圍之內(nèi)，系統(tǒng)會(huì)自動(dòng)提示直到輸入的值達(dá)到要求為止。焊接熱模擬軟件根據(jù)輸入的參數(shù)，自動(dòng)計(jì)算出該點(diǎn)距焊縫中心線的距離和經(jīng)歷的熱循環(huán)曲線。

1.2 Hannerz模型

根據(jù)瑞典斯德哥爾摩皇家技術(shù)學(xué)院N.E.Hannerz教授的研究成果，焊接熱循環(huán)的計(jì)算公式為式(2)～式(4)：

(2)

A=Δt×(500-T0)2(800-T0)2

(3)

B=1 300-2T0

(4)

式中：T0為預(yù)熱溫度， ℃；Tmax為峰值溫度， ℃；Δt為從800～500 ℃的冷卻時(shí)間，s；t為時(shí)間，s。

在焊接熱循環(huán)計(jì)算時(shí)，Hannerz模型需要輸入峰值溫度Tmax，預(yù)熱溫度T0以及從800～500 ℃的冷卻時(shí)間Δt，其中峰值溫度和預(yù)熱溫度的值必須在100～2 000 ℃間變化。根據(jù)輸入的參數(shù)，焊接熱模擬軟件自動(dòng)計(jì)算出所經(jīng)歷的熱循環(huán)曲線。

1.3 Rykalin-2D和Rykalin-3D模型

基于雷卡林對(duì)焊接熱過(guò)程的研究，當(dāng)焊件為二維導(dǎo)熱時(shí)，采用Rykalin-2D模型計(jì)算焊接熱循環(huán)的公式為式(5)～式(9)：

(5)

(6)

b=4π×k×c×p

(7)

(8)

(9)

當(dāng)焊件為三維導(dǎo)熱時(shí)，采用Rykalin-3D模型計(jì)算焊接熱循環(huán)的公式為式(10)～式(14)：

(10)

(11)

(12)

(13)

(14)

式中：Q為焊接線能量，J/cm；c為比熱，J/(g·℃)；ρ為密度，g/cm3；k為熱導(dǎo)率，W/(cm·℃)；d為板厚，cm；T1，T2用于定義冷卻時(shí)間， ℃；t0為預(yù)熱結(jié)束時(shí)間，s；Δt為從溫度T2冷卻到溫度T1的時(shí)間，s。

在焊接熱循環(huán)的計(jì)算時(shí)，Rykalin-2D模型需要輸入材料的物性參數(shù)比熱c、密度ρ、熱導(dǎo)率k，焊接時(shí)的加熱速率ωh、峰值溫度Tmax、預(yù)熱溫度T0、焊接線能量E、峰值溫度停留時(shí)間，定義冷卻時(shí)間的溫度T1，T2、冷卻時(shí)間Δt和和焊件厚度。Rykalin-3D模型則不需要輸入板厚。其中，峰值溫度停留時(shí)間必須在0.01～100 s間變化，焊接線能量E和冷卻時(shí)間Δt兩參數(shù)只需輸入其中一個(gè)，焊接熱模擬軟件會(huì)自動(dòng)計(jì)算出另一個(gè)的值和所經(jīng)歷的熱循環(huán)曲線。

1.4 Rosenthal模型

基于美國(guó)學(xué)者Dr.Daniel Rosenthal 對(duì)厚板焊件三維導(dǎo)熱的研究，采用Rosenthal模型計(jì)算焊接熱循環(huán)的公式為式(15)～式(17)：

(15)

(16)

(17)

式中：Q為焊接電弧的功率，J/s；η為焊接熱效率；v為焊接速度，cm/s；c為比熱，J/(g·℃)；ρ為密度，g/cm3；k為熱導(dǎo)率，W/(cm·℃)；r為距電弧中心的距離，cm；d為某點(diǎn)距焊縫中心線的垂直距離，cm。

在焊接熱循環(huán)的計(jì)算時(shí)，Rosenthal模型需要輸入材料的物性參數(shù)比熱c、密度ρ、熱導(dǎo)率k和焊接時(shí)的預(yù)熱溫度T0、焊接線能量E、焊接速度v，距焊縫中心線的垂直距離d以及焊接熱效率η。根據(jù)輸入的參數(shù)，焊接熱模擬軟件自動(dòng)計(jì)算出所經(jīng)歷的熱循環(huán)曲線。

1.5 指數(shù)冷卻模型

指數(shù)冷卻模型是指試樣在加熱到峰值溫度后，在800 ℃和500 ℃間以指數(shù)方式按指定的速率冷卻。指數(shù)冷卻模型焊接熱循環(huán)的計(jì)算公式為式(18)：

(18)

式中：Tmax為峰值溫度， ℃；Δt為從800 ℃冷卻到500 ℃的時(shí)間，s；t為時(shí)間，s。

在焊接熱循環(huán)的計(jì)算時(shí)，指數(shù)冷卻模型需要輸入峰值溫度Tmax，加熱速率ωh和從800～500 ℃間的冷卻時(shí)間Δt。根據(jù)輸入的參數(shù)，焊接熱模擬軟件自動(dòng)計(jì)算出所經(jīng)歷的熱循環(huán)曲線。

2 焊接熱循環(huán)計(jì)算分析實(shí)例

2.1 不同數(shù)學(xué)模型焊接熱循環(huán)的計(jì)算實(shí)例

本文以12.7 mm厚X80管線鋼板材埋弧焊接的熱模擬為例，分析不同焊接熱循環(huán)計(jì)算模型的特點(diǎn)和應(yīng)用范圍。X80管線鋼的物理性能參數(shù)取室溫值，密度、比熱和熱傳導(dǎo)系數(shù)分別為7.871 g/cm3、0.481 J/(g·℃)和0.653 W/(cm·℃)。焊接工藝參數(shù)為：焊接電壓32 V，焊接電流500 A，焊接速度為11 mm/s，預(yù)熱溫度為100 ℃。加熱速度取為130 ℃/s，峰值溫度取為1 300 ℃，高溫停留時(shí)間為1 s，焊接后從800～500 ℃的冷卻時(shí)間采用德國(guó)鋼鐵學(xué)會(huì)推薦的D.Vwer等人提出的工程應(yīng)用公式，如式(19)～式(21)，在三維熱傳導(dǎo)條件下[5-8]：

(19)

在二維熱傳導(dǎo)條件下：

(20)

臨界厚度的判別公式為：

(21)

式中：η′為不同焊接方法的相對(duì)熱效率，F(xiàn)3和F2為三維和二維熱傳導(dǎo)條件下的接頭系數(shù)，具體見(jiàn)文獻(xiàn)[5]。計(jì)算出材料的臨界板厚為14.73 mm，屬于二維熱傳導(dǎo)形式，采用式(20)計(jì)算從800～500 ℃的冷卻時(shí)間為22.72 s。

利用Gleeble 3系列熱-力學(xué)模擬試驗(yàn)機(jī)中配備的焊接熱模擬軟件，輸入本文所研究的焊接工藝條件，6種數(shù)學(xué)模型計(jì)算所得的熱循環(huán)曲線如圖1所示。由圖可知，在相同的焊接工藝條件下，采用不同的數(shù)學(xué)模型計(jì)算，所得到六條熱循環(huán)曲線存在較大的差異， 6種模型從800～500 ℃的冷卻時(shí)間依次為：16.13、22.81、18.92、4.08、3.38和22.38 s。因此，在進(jìn)行焊接熱模擬試驗(yàn)以前，應(yīng)根據(jù)實(shí)際研究的焊接工藝，選擇適合的數(shù)學(xué)模型計(jì)算焊接熱循環(huán)。六種模型中Rykalin-3D模型和Rosenthal模型計(jì)算所得的冷卻速度對(duì)比其他四種差異性較大，這主要是因?yàn)檫@兩種模型是針對(duì)厚板三維熱傳導(dǎo)形式。因此在選擇數(shù)學(xué)模型計(jì)算焊接熱循環(huán)前，應(yīng)首先明確所研究焊接工藝的實(shí)際傳熱方式。

圖1 6種數(shù)學(xué)模型計(jì)算所得的焊接熱循環(huán)曲線

2.2 數(shù)學(xué)模型的可控性和適用范圍

焊接熱循環(huán)由4個(gè)主要特征參數(shù)決定：加熱速度、峰值溫度、高溫停留時(shí)間、冷卻速度或某一溫度區(qū)間的冷卻時(shí)間。針對(duì)油氣領(lǐng)域常用的管線鋼，加熱速度影響奧氏體均質(zhì)化和碳化物的溶解過(guò)程，進(jìn)而影響到冷卻過(guò)程最終的相變產(chǎn)物和性能。峰值溫度決定了所模擬焊接熱影響區(qū)的特定區(qū)域。高溫停留時(shí)間則是指相變溫度以上的停留時(shí)間，包括加熱和冷卻階段，這主要是影響奧氏體晶粒長(zhǎng)大和合金元素的溶解。但是在Gleeble系列模擬軟件中，高溫停留時(shí)間指的是在峰值溫度的停留時(shí)間。冷卻速度是影響焊接后最終轉(zhuǎn)變組織和性能的重要因素，也是焊接熱循環(huán)模擬的重要變量。一般試驗(yàn)要求的冷卻速度為從800～500 ℃的冷卻速度，也就是冷卻時(shí)相變發(fā)生溫度區(qū)間的冷卻速度。如果考慮馬氏體轉(zhuǎn)變則應(yīng)考慮800～300 ℃的冷卻速度。預(yù)熱溫度也是焊接時(shí)的常見(jiàn)工藝參數(shù)，但是對(duì)于焊接物理模擬熱循環(huán)計(jì)算的影響較小。

6種模型中，F(xiàn)(s,d)模型根據(jù)實(shí)際的測(cè)量結(jié)果，采用分段模擬加熱，因此對(duì)各個(gè)階段焊接熱循環(huán)特征參數(shù)都是可以控制的。Hannerz模型和Rosenthal模型無(wú)法對(duì)加熱速度進(jìn)行設(shè)置。Rosenthal模型無(wú)法直接定義峰值溫度，而是通過(guò)調(diào)整距焊縫中心線的位置獲得不同的峰值溫度。針對(duì)高溫停留時(shí)間，僅僅Rykalin-2D和Rykalin-3D模型考慮了峰值溫度停留時(shí)間對(duì)熱影響區(qū)性能的影響。所有模型都可實(shí)現(xiàn)從800～500 ℃冷卻速度的控制，但是針對(duì)其他溫度段冷卻速度，則只有Rykalin-2D和Rykalin-3D進(jìn)行了考慮。表1總結(jié)了6種模型對(duì)焊接熱循環(huán)特征參數(shù)的輸入控制能力，在焊接熱模擬前可根據(jù)現(xiàn)有的工藝參數(shù)，選擇合適的模型計(jì)算焊接熱循環(huán)。

表1 6種數(shù)學(xué)模型對(duì)焊接熱循環(huán)特征參數(shù)的控制

總的來(lái)說(shuō)，可以實(shí)現(xiàn)所有焊接熱循環(huán)特征參數(shù)的設(shè)定，在模擬材料的物性參數(shù)已知，焊接工藝參數(shù)明確的情況下是計(jì)算焊接熱循環(huán)參數(shù)的最佳選擇，兩者的選擇只需確定材料焊接工藝的傳熱方式。Hannerz模型和Rosenthal模型則分別是在焊接熱循環(huán)控制參數(shù)要求不高的情況下，相對(duì)簡(jiǎn)單的二維和三維傳熱方式計(jì)算公式。指數(shù)冷卻模型是材料物性參數(shù)未知的情況下，計(jì)算焊接熱循環(huán)的選擇。

3 結(jié) 論

本文基于DSI公司Gleeble 3系列熱-力學(xué)模擬試驗(yàn)機(jī)中配備的焊接熱模擬軟件，深入探討了F(s,d)表格，Hannerz，Rykalin-2D，Rykalin-3D，Rosenthal和指數(shù)冷卻6種不同焊接熱循環(huán)計(jì)算模型。在選擇數(shù)學(xué)模型計(jì)算焊接熱循環(huán)前，應(yīng)首先明確所研究焊接工藝的實(shí)際傳熱方式。在焊接材料和焊接工藝明確的情況下，選用Rykalin-2D和Rykalin-3D模型計(jì)算溫度-時(shí)間曲線可以實(shí)現(xiàn)更多焊接熱循環(huán)特征參數(shù)的控制。