曹昱澎 惠 虎 賀寅彪 李 輝

(1.華東理工大學(xué)機(jī)械與動(dòng)力工程學(xué)院承壓系統(tǒng)安全科學(xué)教育部重點(diǎn)實(shí)驗(yàn)室；2.上海核工程研究設(shè)計(jì)院工程設(shè)備所)

大量的壓力容器及管道等承壓設(shè)備都是用鐵素體鋼制成的。鐵素體鋼具有韌脆轉(zhuǎn)變現(xiàn)象，若因事故及過載等原因使設(shè)備運(yùn)行在韌脆轉(zhuǎn)變溫度區(qū)時(shí)，便有可能誘發(fā)嚴(yán)重的脆斷事故。鐵素體鋼在韌脆轉(zhuǎn)變區(qū)的斷裂過程非常復(fù)雜，失效機(jī)理與所處的溫度區(qū)域有關(guān)。在轉(zhuǎn)變區(qū)的中-上區(qū)域，一般延性撕裂與解理斷裂相互競爭；在轉(zhuǎn)變區(qū)的中-下區(qū)域至下平臺(tái)溫度區(qū)，一般發(fā)生純解理斷裂，在解理失穩(wěn)斷裂發(fā)生前沒有明顯的延性撕裂，宏觀表現(xiàn)為純脆斷的失效形式，在斷裂失效之前觀察不到結(jié)構(gòu)的明顯變形，這是非常危險(xiǎn)的。

夏比V型缺口沖擊試驗(yàn)由于試驗(yàn)簡便、費(fèi)用低廉而被廣泛用于鐵素體的韌脆轉(zhuǎn)變研究，至今已有上百年的歷史。但是，夏比沖擊功AKV不能直接關(guān)聯(lián)斷裂韌性KJc，所以無法建立缺陷尺寸與斷裂外載荷之間的關(guān)系。當(dāng)用于結(jié)構(gòu)的強(qiáng)度設(shè)計(jì)或完整性評定時(shí)，基于AKV對材料韌性的要求往往是根據(jù)經(jīng)驗(yàn)來確定的，不夠科學(xué)。因此，建立夏比沖擊功AKV與斷裂韌性的關(guān)系一直是學(xué)術(shù)界關(guān)注的研究熱點(diǎn)[1]。

在韌脆轉(zhuǎn)變區(qū)，夏比沖擊試樣在受沖擊的過程中，由于缺口尖端的應(yīng)力集中程度不及裂紋，通常容易在缺口尖端先發(fā)生一定程度的延性撕裂，延性裂紋的起裂和擴(kuò)展將改變原有缺口尖端的應(yīng)力應(yīng)變分布，進(jìn)而會(huì)對裂尖附近薄弱組織的解理起裂產(chǎn)生影響。作為研究夏比沖擊功與斷裂韌性參量的關(guān)聯(lián)的第一步，需要基于有限元數(shù)值模擬試樣在解理起裂前的延性損傷斷裂過程。筆者以石化設(shè)備常用的16MnR鋼(新國標(biāo)GB713-2008稱為Q345R鋼[2])為例，借助大型有限元分析軟件ABAQUS，考慮沖擊過程中高應(yīng)變率對基體材料的強(qiáng)化效應(yīng)，并耦合延性損傷GTN模型，對韌脆轉(zhuǎn)變區(qū)內(nèi)在下轉(zhuǎn)變溫度區(qū)(T=-86℃)和中-上轉(zhuǎn)變區(qū)(T=-50℃)的夏比沖擊試驗(yàn)進(jìn)行三維彈塑性數(shù)值模擬，重點(diǎn)分析延性裂紋的起裂擴(kuò)展。

1 有限元分析

1.1夏比V型缺口試樣有限元模型

由于試樣具有對稱性，模型取原試樣的1/4建立，模型尺寸符合文獻(xiàn)[3]中的規(guī)定(圖1)。缺口尖端最小單元尺寸為0.029mm×0.015mm×0.091mm。因?yàn)榭拷嚇颖砻婢惺刃。瑧?yīng)力梯度大，所以表面網(wǎng)格密；試樣中部拘束度大，應(yīng)力梯度小，因此試樣厚度中心處網(wǎng)格稀疏。沿試樣厚度方向從表面至試樣中部共劃分30層網(wǎng)格。模型共有36 690個(gè)C3D8R和C3D6R單元，40 961個(gè)節(jié)點(diǎn)，沖錘和支座設(shè)為剛體，沖錘的初始速度按照試驗(yàn)設(shè)定為5.24m/s。

圖1 夏比V型缺口試樣的有限元模型

1.2材料參數(shù)

16MnR鋼在-86、-50℃時(shí)的真應(yīng)力-真塑性應(yīng)變曲線如圖2所示。彈性模量E=200GPa，泊松比ν=0.3。

圖2 16MnR的真應(yīng)力-真塑性應(yīng)變曲線

試樣在受沖擊過程中快速變形，各點(diǎn)的應(yīng)變率不同，缺口尖端的應(yīng)變率甚至高達(dá)1 000/s數(shù)量級[4]。16MnR屬于應(yīng)變率敏感材料，屈服強(qiáng)度和抗拉強(qiáng)度隨應(yīng)變率的升高而明顯升高。本研究測試了16MnR在340、1 860/s應(yīng)變率下的沖擊拉伸性能。在數(shù)值模擬時(shí)采用Cowper-Symonds過應(yīng)力模型來考慮材料的應(yīng)變率強(qiáng)化效應(yīng)，以便插值取得各應(yīng)變率下的材料本構(gòu)關(guān)系[5]：

(1)

式中p——需要擬合的材料參數(shù)；

σ0——準(zhǔn)靜態(tài)下的真應(yīng)力。

(2)

(3)

式中fc——臨界孔洞體積分?jǐn)?shù)；

fF——失效時(shí)的孔洞體積分?jǐn)?shù)；

f*——等效孔洞體積分?jǐn)?shù)；

q1，2，3——材料的本構(gòu)參數(shù)；

σeq——宏觀Mises等效應(yīng)力；

σm——宏觀靜水應(yīng)力；

σy——材料屈服應(yīng)力。

損傷演化方程由新孔洞形核和孔洞長大兩部分函數(shù)組成：

(4)

df=f0

式中fN——適合孔洞形核的顆粒體積分?jǐn)?shù)；

f0——初始孔洞體積分?jǐn)?shù)；

SN——分布的標(biāo)準(zhǔn)差；

εN——平均孔洞形核應(yīng)變；

εp——臨界塑性應(yīng)變；

GTN模型共有9個(gè)參量需要標(biāo)定，分為兩類，第一類屬于本構(gòu)參量，包括q1、q2、q3，通常q1=1.50，q2=1.00，q3=(q1)2=2.25；第二類屬于材料參數(shù)，形核參數(shù)εN和SN通常取0.3、0.1。初始空洞體積分?jǐn)?shù)f0可以通過材料的化學(xué)成分用Franklin公式估算[9]。fN、fc、fF需要經(jīng)過模擬與試驗(yàn)結(jié)果之間反復(fù)比較試算得到。經(jīng)反復(fù)調(diào)試標(biāo)定得到16MnR的GTN參量(表1)。

表1 16MnR鋼的GTN模型標(biāo)定參量

2 結(jié)果與討論

2.1模擬的載荷-撓度曲線與延性撕裂

有限元模擬的Mises云圖如圖3所示，圖中缺口尖端的圓弧區(qū)域?yàn)镚TN模型預(yù)測的延性裂紋擴(kuò)展，裂紋從試樣中截面(厚度方向)的缺口尖端起裂并擴(kuò)展，剩余韌帶將承擔(dān)沖錘繼續(xù)下壓的載荷。夏比沖擊試樣受沖擊的過程可分為彈性變形、整體屈服、延性裂紋起裂和延性裂紋穩(wěn)定擴(kuò)展、突然的解理斷裂和最后剪切唇的形成。

圖3 夏比沖擊模擬的Mises云圖

如圖4a所示，-86℃下夏比示波沖擊試驗(yàn)記錄的載荷-撓度曲線表現(xiàn)出明顯的純脆斷特征，在曲線的彈性段試樣就發(fā)生突然失穩(wěn)的解理斷裂，試樣斷面上未發(fā)現(xiàn)明顯的延性裂紋擴(kuò)展，而圖4b所示的是典型的夏比沖擊韌脆轉(zhuǎn)變區(qū)中-上區(qū)域的載荷-撓度曲線，試樣經(jīng)歷了彈性變形、明顯塑性變形、延性裂紋起裂擴(kuò)展后，載荷-撓度曲線突降，表明發(fā)生了脆斷。-86、-50℃下模擬的夏比沖擊試樣的載荷-撓度曲線與對應(yīng)溫度下解理起裂前的試驗(yàn)曲線吻合得較好。對于-86℃下的有限元計(jì)算，在對應(yīng)試驗(yàn)試樣脆斷的載荷步下，模擬預(yù)測沒有發(fā)生延性起裂，與試驗(yàn)結(jié)果一致。

a.-86℃

b.-50℃

-50℃下延性裂紋擴(kuò)展與試樣撓度的關(guān)系如圖5所示，3個(gè)試樣斷裂時(shí)的撓度分別為3.5、4.0、5.0mm，斷面上延性裂紋最大擴(kuò)展長度分別為0.60、0.77、0.88mm。在對應(yīng)的撓度下，GTN模型預(yù)測的擴(kuò)展長度分別為0.75、1.03、1.55mm。考慮到計(jì)算機(jī)處理能力有限，試樣韌帶方向上的有限元網(wǎng)格不可能劃分得非常細(xì)，而且模型是以逐個(gè)刪除單元的方式漸進(jìn)地模擬延性裂紋的擴(kuò)展，因此，可以認(rèn)為對延性裂紋擴(kuò)展的預(yù)測結(jié)果是可接受的。

圖5 延性裂紋長度與撓度的關(guān)系

以上計(jì)算結(jié)果表明三維彈塑性有限元分析結(jié)合局部損傷的GTN模型可以很好地預(yù)測延性裂紋的起裂與擴(kuò)展。由于在韌脆轉(zhuǎn)變溫度區(qū)解理斷裂的發(fā)生屬于一種概率事件，所以有限元不能直接模擬解理斷裂。

2.2夏比沖擊試樣中的應(yīng)變率

有限元計(jì)算發(fā)現(xiàn)，夏比沖擊試樣在受到?jīng)_錘高速?zèng)_擊的過程中，在延性裂紋即將起裂前，試樣中截面缺口尖端附近的應(yīng)變率可達(dá)1 900/s，試樣表面處應(yīng)變率略低，在缺口尖端附近的應(yīng)變率約1 100/s，如圖6所示。應(yīng)變率達(dá)到上百的高應(yīng)變率區(qū)集中在缺口尖端附近1mm的范圍內(nèi)，遠(yuǎn)離缺口尖端應(yīng)變率下降得很快。缺口尖端是材料發(fā)生延性損傷或解理起裂的斷裂過程區(qū)，可見模擬夏比沖擊時(shí)考慮應(yīng)變率對材料的強(qiáng)化效應(yīng)是非常有必要的。

圖6 張開應(yīng)變率沿韌帶的分布

2.3延性裂紋起裂和擴(kuò)展對裂尖場的影響

夏比V型缺口試樣中間平面上，沿韌帶的最大主應(yīng)力分布(距離初始缺口位置)如圖7所示。從圖7可以看出，一旦延性裂紋從缺口尖端起裂對裂尖的最大主應(yīng)力分布影響很大，最大主應(yīng)力峰值由約1 600MPa升高到1 900MPa。比較延性裂紋擴(kuò)展長度分別為0.30、1.05、1.50mm時(shí)的裂紋前端的最大主應(yīng)力分布，發(fā)現(xiàn)在延性裂紋擴(kuò)展過程中沿韌帶上的最大主應(yīng)力峰值略有降低，它的位置逐漸略微遠(yuǎn)離裂紋尖端，高應(yīng)力區(qū)略微增大。因?yàn)榻饫頂嗔寻l(fā)生的概率大小受應(yīng)力大小和高應(yīng)力控制區(qū)域大小的共同作用的影響，所以從力學(xué)場的角度分析，在筆者所計(jì)算的范圍內(nèi)，夏比沖擊試樣的延性裂紋起裂將明顯促使解理斷裂的發(fā)生，一旦延性裂紋起裂后若未隨之發(fā)生解理斷裂，隨著延性裂紋擴(kuò)展，裂尖發(fā)生解理斷裂的概率受延性撕裂的影響不大，即延性裂紋的擴(kuò)展不會(huì)顯著地促進(jìn)或抑制解理起裂。

圖7 最大主應(yīng)力在裂紋前緣的分布

2.4能量分離

常規(guī)的夏比沖擊試驗(yàn)給出的沖擊吸收功AKV包含了材料彈性應(yīng)變能、塑性變形能、斷裂能和形成剪切唇所需的能量，這些能量很難通過試驗(yàn)將它們彼此分離出來。借助示波沖擊試驗(yàn)也只能從載荷-位移曲線上粗略地將彈性能和試樣整體屈服后的吸收功區(qū)分開。借助ABAQUS有限元分析，以-50℃下的夏比沖擊試驗(yàn)為例，可以將試樣在沖擊過程中的各部分能量分離出來，如圖8所示。隨著試樣變形撓度的增加，彈性能保持約2J不變，所占比例最小，塑性變形能占總功的比例最大，余下的為斷裂能等其他能量。由此可見，總吸收功AKV大部分都消耗在了塑性變形上，而其他能量占AKV的比重都很小。

圖8 試樣的總吸收功、塑性能和彈性能隨撓度的變化

3 結(jié)論

3.1考慮應(yīng)變率對材料的強(qiáng)化效應(yīng)和GTN延性損傷模型，借助ABAQUS有限元軟件可以較準(zhǔn)確地模擬韌脆轉(zhuǎn)變溫度區(qū)的夏比沖擊試驗(yàn)，預(yù)測延性裂紋的起裂和擴(kuò)展。對比夏比沖擊試樣斷面發(fā)現(xiàn)，GTN模型預(yù)測的延性裂紋擴(kuò)展長度是可接受的。延性裂紋從試樣中截面缺口尖端處起裂，裂紋起裂使原缺口尖端的主應(yīng)力明顯增大，對解理起裂有促進(jìn)作用，在延性裂紋擴(kuò)展過程中，裂尖的應(yīng)力分布變化不大。

3.2夏比沖擊試樣的缺口根部附近的應(yīng)變率可達(dá)1 000/s，高應(yīng)變率區(qū)集中在缺口尖端附近很小的區(qū)域內(nèi)。

3.3模擬分離得到夏比沖擊功AKV包含的各部分能量，發(fā)現(xiàn)在沖擊過程中塑性功消耗占絕大部分。

[1] 秦江陽，王印培，柳曾典.JIC和AKV之間的關(guān)系研究[J].材料工程，2001，18(2)：15～19.

[2] GB 713-2008，鍋爐和壓力容器用鋼板[S].北京：中國標(biāo)準(zhǔn)出版社，2008.

[3] GB/T 229-2007，金屬夏比缺口沖擊試驗(yàn)方法[S].北京：中國標(biāo)準(zhǔn)出版社，2007.

[4] Wilshaw T R.The Deformation and Fracture of Mild Steel Charpy Specimens[J].Journal of Iron and Steel Institute，1966，204(9)：936～942.

[5] Bodner S，Symonds P.Plastic Deformations in Impact and Impulsive Loading of Beams[C].Plasticity Proceedings of the Second Symposium on Naval Structural Mechanics.Rhode Island：Proceedings of the Second Symposium on Naval Structural Mechanics，1960：488～500.

[6] Gurson A L.Continuum Theory of Ductile Rupture by Void Nucleation and Growth：Part 1 Yield Criteria and Flow Rules for a Porous Ductile Media[J].Journal of Engineering Materials and Technology，1977，99(1)：2～15.

[7] Tvergaard V.On Localization in Ductile Materials Containing Spherical Voids[J].International Journal of Fracture，1982，18(4)：237～252.

[8] Tvergaard V，Needleman A.Analysis of the Cup-cone Fracture in a Round Tensile Bar[J].Acta Metallurgica，1984，32(1)：157～169.

[9] Franklin A G.Comparison Between a Quantitative Microscope and Chemical Methods for Assessment of Non-metallic Inclusions[J].Journal of Iron Steel and Institute，1969，207(2)：181～186.