幾何拘束對鎳基合金裂紋擴展速率影響的有限元分析

,,

(1. 西安科技大學(xué) 理學(xué)院，西安 710054； 2. 西安科技大學(xué) 機械工程學(xué)院，西安 710054)

鎳基合金及其焊縫金屬被廣泛用于核電一回路壓力容器及異種金屬焊接接頭中，這些焊接接頭的應(yīng)力腐蝕開裂(SCC)是核電構(gòu)件的一種重要破壞形式[1]。為實現(xiàn)準確的焊接構(gòu)件完整性評價和壽命預(yù)測，精確預(yù)測在役構(gòu)件的SCC裂紋擴展速率(CGR)顯得尤為重要[2]。受到異種金屬焊接接頭材料力學(xué)不均勻性和殘余應(yīng)力的影響，用慢應(yīng)變速率試驗(SSRT)等方法精確測定焊接構(gòu)件的應(yīng)力腐蝕CGR非常困難。工程實踐表明，即使是同種材料，裂尖拘束效應(yīng)對構(gòu)件的CGR也有很大影響。通常將裂紋尺寸、試樣或結(jié)構(gòu)幾何與加載方式等引起的拘束效應(yīng)稱為幾何拘束，幾何拘束可進一步分為面外拘束和面內(nèi)拘束[3]。目前的研究主要針對裂尖拘束效應(yīng)對材料與結(jié)構(gòu)斷裂行為的影響[4]、拘束參數(shù)的定量化表征[5-7]、復(fù)合拘束效應(yīng)對高溫蠕變裂紋CGR的影響規(guī)律和機理[8-10]，關(guān)于幾何拘束效應(yīng)對SCC裂紋擴展速率影響的研究還較少，究其原因，是由于應(yīng)力腐蝕裂紋擴展的影響因素眾多，要從試驗中單獨剝離出幾何拘束對CGR的影響存在一定難度。

鑒于CGR試驗的費時費力和電子計算機計算能力的飛速發(fā)展，數(shù)值模擬技術(shù)被越來越多地應(yīng)用于裂紋擴展速率的分析，本工作以鎳基合金600為研究對象，建立了不同厚度和裂紋長度緊湊拉伸(CT)試樣的有限元模型，對不同拘束狀態(tài)的裂紋尖端力學(xué)場和表征SCC裂紋擴展速率的主要參量進行了研究，進一步分析了試樣厚度和裂紋長度這兩種幾何拘束條件對鎳基合金SCC裂紋擴展速率的影響，以期對長期服役的低拘束核電焊接構(gòu)件的壽命預(yù)測和完整性評價提供借鑒。

1 模型的建立

1.1 幾何和有限元模型

選用帶預(yù)制裂紋的CT試樣。按照美國ASTM E399-90標準，1T-CT試樣的幾何尺寸如圖1所示，寬度W=50 mm，通過改變厚度B(B=0.05W、0.1W、0.2W、0.3W、0.4W和0.5W)來改變面外拘束，通過改變裂紋長度a(a=0.5W、0.6W和0.7W)來改變面內(nèi)拘束。

圖1 1T-CT試樣的幾何尺寸(W=50 mm,a=0.5 W)Fig. 1 Geometric size of 1T-CT specimen(W=50 mm,a=0.5 W)

試樣全局模型和裂尖局部區(qū)域的有限元網(wǎng)格如圖2所示，均采用8節(jié)點六面體實體單元(C3D8)。為獲得較為詳細準確的裂尖應(yīng)力應(yīng)變，在裂尖處對網(wǎng)格進行了細化，以B=0.5W，a=0.5W的1T-CT試樣為例，其全局模型單元總數(shù)為177 408，裂尖局部區(qū)域單元總數(shù)為83 540。

在對各斷裂力學(xué)參量進行對比分析時，取恒定載荷作為各有限元模型的加載條件。通常對試樣施加恒定載荷[11]使裂尖的應(yīng)力強度因子KI=30 MPa·m1/2。為研究載荷水平對SCC裂紋擴展速率的影響，還取恒載荷KI=10 MPa·m1/2和KI=20 MPa·m1/2。

(a) 1T-CT試樣(b) 裂尖區(qū)域圖2 1T-CT試樣的有限元網(wǎng)格Fig. 2 Mesh of the 1T-CT specimen (a) and the detail around crack tip (b)

1.2 材料模型

鎳基合金600屬于冪硬化材料，本構(gòu)關(guān)系通常用Ramberg-Osgood關(guān)系來表征

(1)

式中：ε為真實應(yīng)變，包含彈性和塑性應(yīng)變；σ為真實應(yīng)力；σ0和ε0分別為屈服應(yīng)力和屈服應(yīng)變；α和n分別為材料的應(yīng)變硬化系數(shù)和硬化指數(shù)。在340 ℃高溫高壓水環(huán)境中，鎳基合金600相關(guān)力學(xué)參數(shù)如下[2]:楊氏模量(E)為189 500 MPa,泊松比(v)為0.286,屈服應(yīng)力(σ0)為436 MPa,應(yīng)變硬化系數(shù)(α)為3.075,硬化指數(shù)(n)為6.495。

1.3 SCC裂紋擴展速率定量分析

根據(jù)滑移溶解理論建立的FORD-ANDRESRN模型[12-13]，裂紋擴展速率為

(2)

又有

(3)

式中:Mmol為金屬原子的摩爾質(zhì)量，ρ為金屬的密度，F(xiàn)為法拉第常數(shù)，z為溶解過程中的電荷變化，i0為裸露金屬表面的氧化電流密度，t0為環(huán)境中溶解電流衰減開始的時間常數(shù)。f為裂尖表面氧化膜的斷裂應(yīng)變。

對于式(2)，距離擴展裂尖前緣特征距離r0處的裂尖應(yīng)變εct可以用塑性應(yīng)變εp來代替，有

εct=εp|r=r0(4)

式中:r為裂紋尖端前緣的距離，當裂紋擴展時，r的變化會導(dǎo)致裂尖的塑性應(yīng)變增大，則有

(5)

式中:dεp/da是擴展裂紋前緣，裂尖特征距離r0處的塑性應(yīng)變率。

將式(5)代入式(2)，得到材料在高溫水環(huán)境中的SCC裂紋擴展速率為

(6)

目前，可以用有限元方法對特征距離r0處的塑性應(yīng)變率dεp/da進行求解[2]，進而得到含缺陷構(gòu)件的SCC裂紋擴展速率并對其進行壽命預(yù)測，而特征距離r0的選定和裂尖塑性區(qū)有很大關(guān)系，一般在幾十微米以內(nèi)[14]。

考慮到裂紋延展中的應(yīng)變重新分配，SHOJI等[15-16]根據(jù)裂尖應(yīng)變梯度理論提出

(7)

在應(yīng)力強度因子K不變的情況下，式(7)的第一項可忽略不計。即

(8)

從式(8)可以得到，特征距離r0處的塑性應(yīng)變率dεp/dr成為影響裂紋SCC行為的唯一力學(xué)因素，對于某一材料裂紋尖端固定的電化學(xué)環(huán)境，裂尖的塑性應(yīng)變率與裂紋的擴展速率實質(zhì)上是一致的。

2 結(jié)果與討論

2.1 試樣厚度對SCC裂尖塑性應(yīng)變率的影響

不同厚度1T-CT試樣裂紋尖端的塑性應(yīng)變率沿裂紋前端的分布如圖3所示。計算時，用裂尖的法向塑性應(yīng)變ε22來代替式(5)中的裂尖塑性應(yīng)變εp。隨著裂尖恒載荷KI的增大，裂尖塑性應(yīng)變率dεp/dr逐漸增大，即載荷越大，裂紋擴展速率越大，這與LU等[16]研究的冷加工316L不銹鋼的應(yīng)力腐蝕裂紋擴展速率結(jié)果是一致的。低載荷水平條件下，試樣厚度引起的裂尖拘束效應(yīng)對裂紋擴展速率的影響較明顯。特征距離范圍內(nèi)，r0<100 μm時，不同厚度試樣的法向塑性應(yīng)變率幾乎相等。當KI≥20 MPa·m1/2時，在高載荷水平條件下，不同厚度試樣裂尖的裂紋擴展速率趨于一致，拘束效應(yīng)對裂紋擴展速率幾乎沒有什么影響。從某種意義上來說，裂尖塑性應(yīng)變率屬于裂紋擴展驅(qū)動力的一種，是促進裂紋擴展的；但是隨著拘束和載荷的減少，淺裂紋結(jié)構(gòu)的斷裂韌性甚至能增大好幾倍，即低載荷薄試樣的裂紋擴展阻力比厚試樣的大。因此，對于特征距離r0處差別不大的塑性應(yīng)變率，厚試樣的裂紋擴展速率大于薄試樣的。這與談建平等[9-10]關(guān)于高溫轉(zhuǎn)子鋼蠕變裂紋擴展速率和304不銹鋼疲勞裂紋擴展速率的試驗結(jié)果吻合。當試樣的厚寬比B/W≥0.3時，裂尖塑性應(yīng)變率比較接近，說明厚試樣的裂尖拘束效應(yīng)對裂紋擴展速率的影響不明顯。

(a) KI=10 MPa·m1/2(b) KI=20 MPa·m1/2(c) KI=30 MPa·m1/2圖3 不同KI條件下試樣厚度對裂尖塑性應(yīng)變率的影響Fig. 3 Effect of specimen thickness on the normal plastic strain rates at crack tips with different KI

2.2 裂紋長度對SCC裂尖應(yīng)變率的影響

由于低載荷水平下試樣厚度所引起的拘束效應(yīng)比較明顯，研究了KI=10 MPa·m1/2時，裂紋長度對1T-CT試樣裂尖塑性應(yīng)變率的影響。由圖4可見:特征距離范圍內(nèi)，r0<100 μm時，不同厚度深裂紋試樣的塑性應(yīng)變率比較接近,裂紋長度a越大，裂尖應(yīng)變率越小，長裂紋的擴展速率較小。試樣較薄時，裂紋長度對SCC裂尖塑性應(yīng)變率幾乎沒有影響；試樣較厚時，裂紋長度所引起的拘束效應(yīng)才顯現(xiàn)出來[17]。對比圖4和圖3，對于試樣厚度和裂紋長度兩種不同的拘束條件，試樣厚度的拘束水平對裂紋擴展速率的影響大于裂紋長度對其的影響。

(a) B=0.05 W(b) B=0.3 W(c) B=0.5 W圖4 試樣厚度不同時,裂紋長度對裂尖法向塑性應(yīng)變的影響Fig. 4 Effect of crack length on the normal plastic strain rates at crack tips with different thickness of sample

3 結(jié)論

(1) 不同厚度試樣裂尖特征距離范圍內(nèi)的法向塑性應(yīng)變率非常接近，厚試樣的裂紋擴展速率高于薄試樣的。低載荷水平條件下，試樣厚度引起的裂尖拘束效應(yīng)對裂紋擴展速率的影響較明顯。高載荷水平條件下，不同厚度試樣的裂尖拘束度趨于一致。B/W≥0.3時，試樣厚度引起的裂尖拘束效應(yīng)對CGR幾乎沒有影響。

(2) 與薄試樣相比，試樣較厚時，裂紋長度所引起的拘束效應(yīng)較明顯。對于試樣厚度和裂紋長度兩種不同的拘束條件，試樣厚度引起的拘束效應(yīng)對裂紋擴展速率的影響大于裂紋長度對其的影響。

(3) 核電焊接構(gòu)件的SCC裂尖長期處于低拘束低載荷狀態(tài)，要對其裂紋擴展速率進行準確的定量預(yù)測，裂尖拘束效應(yīng)對裂紋擴展速率的影響不容忽視。