趙延義，王澤武，范海貴，劉培啟

(大連理工大學(xué) 化工學(xué)院，遼寧大連 116024)

0 引言

第三代壓水堆AP1000發(fā)生破口失水事故時(shí)，低溫安注冷卻水會(huì)通過(guò)直接安注(Direct Vessel Injection，DVI)接管注入堆芯，從而使得DVI接管區(qū)域存在顯著的溫度梯度，并產(chǎn)生過(guò)大的熱應(yīng)力。與此同時(shí)，系統(tǒng)內(nèi)仍保持著較高的介質(zhì)壓力，即發(fā)生承壓熱沖擊(Pressurized Thermal Shock，PTS)現(xiàn)象，因此，DVI接管區(qū)域工作環(huán)境苛刻，且結(jié)構(gòu)復(fù)雜，屬于高應(yīng)力區(qū)[1]。當(dāng)核電廠運(yùn)行至壽期末時(shí)，也會(huì)由于反應(yīng)堆壓力容器(Reactor Pressure Vessel，RPV)長(zhǎng)期受到強(qiáng)烈的中子輻照，使得材料的斷裂韌性下降，如果此時(shí)發(fā)生PTS事件，有可能導(dǎo)致DVI接管與RPV筒體連接區(qū)形成裂紋缺陷，進(jìn)而引發(fā)安全事故[2]。由此可見(jiàn)，DVI接管裂紋強(qiáng)度的準(zhǔn)確計(jì)算及評(píng)估對(duì)于確保核電安全運(yùn)行具有重要意義，是當(dāng)前需要迫切解決的問(wèn)題之一。

杜青等[3]基于Schwartz-Neuman交替法研究了壓力容器內(nèi)接管表面裂紋的應(yīng)力強(qiáng)度因子分布，但研究限于線彈性斷裂力學(xué)理論分析。然而，對(duì)于韌性或彈塑性材料制成的承壓結(jié)構(gòu)，在斷裂發(fā)生之前，裂紋周邊會(huì)發(fā)生較大的塑性變形，應(yīng)力強(qiáng)度因子不適用于確定裂紋周圍的真實(shí)應(yīng)力狀態(tài)[4-5]，因此進(jìn)行J積分彈塑性斷裂力學(xué)分析是十分必要的。胡序春等[6]基于J積分理論，研究了不同工況、不同裂紋角度對(duì)裂紋尖端處J積分的影響，表明J積分計(jì)算對(duì)于判斷裂紋是否失穩(wěn)擴(kuò)展具有重要價(jià)值。

針對(duì)RPV斷裂行為，張麗屏等[1]進(jìn)行了溫度和壓力作用下RPV接管嘴內(nèi)隅角應(yīng)力強(qiáng)度因子計(jì)算研究，但未進(jìn)行彈塑性J積分計(jì)算；楊曉華[7]進(jìn)行了失水事故下RPV彈塑性分析，并計(jì)算了筒體內(nèi)表面軸向半橢圓形裂紋最深點(diǎn)的J積分值，但未考慮DVI接管對(duì)J積分值的影響。白鑫等[8]進(jìn)行了雙裂紋在核管道中應(yīng)力強(qiáng)度因子的相互影響研究；姚安林等[9]基于三維虛擬裂紋閉合技術(shù)研究了多裂紋輸氣管道附屬裂紋對(duì)于主裂紋前緣各點(diǎn)的干涉影響，但未在含有接管的結(jié)構(gòu)不連續(xù)區(qū)域進(jìn)行相應(yīng)的雙裂紋相互干涉效應(yīng)研究。

基于以前的研究，本文針對(duì)DVI接管與RPV筒體連接區(qū)的特殊結(jié)構(gòu)和高應(yīng)力區(qū)，以及RCC-M半橢圓形表面裂紋特征參數(shù)，研究熱力耦合作用下的高應(yīng)力區(qū)表面裂紋彈塑性斷裂強(qiáng)度J積分?jǐn)?shù)值計(jì)算方法，進(jìn)而探究溫差、裂紋參數(shù)和雙裂紋相互干涉效應(yīng)對(duì)J積分的影響規(guī)律，為RPV缺陷完整性評(píng)價(jià)提供理論參考。

1 J積分?jǐn)?shù)值計(jì)算模型開(kāi)發(fā)

1.1 AP1000壓力殼模型

AP1000的RPV結(jié)構(gòu)如圖1(a)所示，并采用如圖1(b)所示“四進(jìn)兩出”布置，入口、出口以及DVI接管兩兩對(duì)稱分布[10-11]。正常運(yùn)行工況下，冷卻劑從入口進(jìn)入，帶走堆芯產(chǎn)生的熱量，再?gòu)某隹诹鞒?。?dāng)壓水堆一回路冷卻劑系統(tǒng)發(fā)生破口失水事故時(shí)，反應(yīng)堆停堆，應(yīng)急冷卻水則從DVI接管注入，對(duì)RPV堆芯進(jìn)行緊急冷卻。

(a)主視圖

(b)俯視圖圖1 AP1000的RPV結(jié)構(gòu)示意Fig.1 Structural diagram of RPV for AP1000

1.2 有限元建模及表面裂紋模型

考慮到RPV結(jié)構(gòu)和載荷的對(duì)稱性，為了降低計(jì)算規(guī)模、提高計(jì)算效率，對(duì)圖1所示核壓力殼結(jié)構(gòu)進(jìn)行合理簡(jiǎn)化，忽略上、下封頭及其他零部件，并選取RPV結(jié)構(gòu)的1/4進(jìn)行三維有限元建模，如圖2所示。有限元模型包括RPV筒體、DVI接管、支座、冷卻劑入口接管以及出口接管。

根據(jù)RCC-M ZG1100規(guī)定，應(yīng)在壓力容器殼體、接管、封頭和焊縫處考慮典型的半橢圓表面裂紋，本文著重研究DVI接管與RPV筒體連接區(qū)的半橢圓表面裂紋。對(duì)核壓力殼整體結(jié)構(gòu)采用Solid 187單元進(jìn)行四面體網(wǎng)格劃分(如圖2所示)，并在DVI接管與RPV交界區(qū)域建立與接管軸向平行的半橢圓三維裂紋。為了保證裂紋計(jì)算精度，先采用“影響球”將DVI接管附近的網(wǎng)格進(jìn)行細(xì)化，如圖3(b)所示；然后在DVI接管圓角的中心區(qū)域建立局部坐標(biāo)系，使得X軸與過(guò)渡弧面法線平行，Z軸與裂紋前緣方向平行，Y軸為裂紋前緣擴(kuò)展法向，如圖3(c)所示；最后采用Solid 186六面體單元對(duì)裂紋局部區(qū)域進(jìn)行網(wǎng)格劃分，并設(shè)置半橢圓形裂紋的長(zhǎng)、短半軸長(zhǎng)度、積分區(qū)域半徑、裂紋前緣單元數(shù)、積分區(qū)域環(huán)向和徑向網(wǎng)格劃分?jǐn)?shù)，如圖3(d)所示。

圖2 RPV筒體與DVI接管三維有限元模型Fig.2 3D finite element model of the RPV shell andDVI nozzle

為了研究雙裂紋相互干涉作用，在主裂紋左側(cè)距離為l處按照?qǐng)D3所示相同的方法建立附屬裂紋，如圖4所示。同時(shí)，為了消除網(wǎng)格大小對(duì)計(jì)算結(jié)果的影響，進(jìn)行了多次細(xì)化網(wǎng)格試算，確定了整體結(jié)構(gòu)單元尺寸為46 mm，當(dāng)長(zhǎng)半軸長(zhǎng)度c=10 mm、短半軸深度a=5 mm,且僅存在主裂紋時(shí)，得到單元536 158個(gè)，節(jié)點(diǎn)786 348個(gè)。

圖3 表面裂紋模型示意Fig.3 Schematic diagram of the surface crack model

圖4 雙裂紋結(jié)構(gòu)示意Fig.4 Structural diagram of the double cracks

1.3 材料模型及邊界條件

1.3.1 彈塑性材料模型

AP1000壓力殼材料為SA508Gr.3低合金鋼，其隨溫度變化的材料性能參數(shù)如表1[12]所示。并根據(jù)文獻(xiàn)[12]，選用雙線性各向同性硬化彈塑性材料模型進(jìn)行彈塑性有限元數(shù)值計(jì)算。

表1 SA508Gr.3材料性能參數(shù)Tab.1 Material parameters of the SA508Gr.3

1.3.2 載荷與約束

在DVI接管注水過(guò)程中，當(dāng)RPV系統(tǒng)壓力降低至約4.9MPa時(shí)，蓄壓安注箱投入工作，此時(shí)注入的冷卻水溫度較低，而RPV內(nèi)部溫度較高，內(nèi)外壁最大溫差可達(dá)81 ℃[13]，是承壓熱沖擊最為劇烈的工況之一。為計(jì)算此最苛刻工況下的RPV筒體與DVI接管連接區(qū)裂紋斷裂特性，先在RPV筒體和DVI接管內(nèi)、外壁壁面各施加恒定溫度，內(nèi)、外壁溫度為116.85,197.85 ℃(溫差81 ℃)，從而得到整體結(jié)構(gòu)最苛刻時(shí)的溫度場(chǎng)；然后基于ANSYS Workbench平臺(tái)施加4.9 MPa流體壓力，采用間接法實(shí)現(xiàn)整體結(jié)構(gòu)熱力耦合數(shù)值計(jì)算。

由于對(duì)封頭、接管進(jìn)行了簡(jiǎn)化處理，因此需要在RPV頂部、底部、DVI接管入口、冷卻劑接管入口及出口按照式(1)分別施加等效軸向拉應(yīng)力pT,pB,pDVI,pc和ph，載荷施加方式如圖5所示。

圖5 載荷與約束示意Fig.5 Schematic diagram of the loads and constraints

(1)

式中,peq為等效軸向拉應(yīng)力pT,pB,pDVI,pc和ph，MPa;p為介質(zhì)壓力，MPa;Di為筒體和接管的內(nèi)徑，mm；Do為筒體和接管的外徑，mm。

同時(shí),在圖5中施加邊界約束條件，包括在支座底部端面施加軸向約束、在對(duì)稱面上施加對(duì)稱約束。

1.4 J積分計(jì)算模型

對(duì)于一個(gè)二維問(wèn)題，J積分的區(qū)域積分表示為:

(2)

q方向?yàn)榱鸭y尖端前方局部坐標(biāo)系的x軸方向，在沿積分路徑Γ的節(jié)點(diǎn)處取q向量為0，在積分路徑Γ內(nèi)部，除與Γ直接相連的中間節(jié)點(diǎn)外的所有節(jié)點(diǎn)取q向量為單位向量，這些帶有單位向量的節(jié)點(diǎn)稱為虛擬裂紋擴(kuò)展節(jié)點(diǎn)。

對(duì)于本文三維模型問(wèn)題，J積分的區(qū)域積分表示變成了體積積分，它同樣是在一組元素上計(jì)算的。虛擬裂紋擴(kuò)展節(jié)點(diǎn)，也被稱為裂紋尖端節(jié)點(diǎn)組件，是J積分計(jì)算中最重要的輸入數(shù)據(jù)元素之一，對(duì)于三維裂紋問(wèn)題，裂紋尖端的節(jié)點(diǎn)組件由沿裂紋前緣的節(jié)點(diǎn)組成。在裂紋有限元數(shù)值計(jì)算中，J積分計(jì)算主要基于SHIH等[14]的區(qū)域積分法，區(qū)域積分公式既適用于二維問(wèn)題的面積積分，也適用于三維問(wèn)題的體積積分，而且體積積分比輪廓積分和曲面積分具有更高的精度和計(jì)算效率。

2 計(jì)算結(jié)果與討論

2.1 模型驗(yàn)證

X80是國(guó)內(nèi)外輸氣管道的主導(dǎo)鋼材，文獻(xiàn)[15]基于ABAQUS對(duì)X80壓力管道表面裂紋進(jìn)行了J積分計(jì)算分析。為了驗(yàn)證本文開(kāi)發(fā)的數(shù)值模型的正確性，采用上述裂紋設(shè)置方法，基于ANSYS Workbench平臺(tái)在X80管道外表面建立與文獻(xiàn)[15]相同的半橢圓形表面裂紋，裂紋深度a=6 mm，半長(zhǎng)c分別為12，18 mm，即裂紋深長(zhǎng)比a/c分別為1/2,1/3，管道外徑為1 219 mm，壁厚18.4 mm，內(nèi)壓p=12 MPa。為了簡(jiǎn)便準(zhǔn)確地探究三維裂紋J積分的影響參數(shù)，以下分析時(shí)將裂紋前緣長(zhǎng)度s做歸一化處理。需要注意的是，文獻(xiàn)[15]中的計(jì)算結(jié)果是裂紋前緣J積分隨裂紋角度的變化，為了便于對(duì)比分析，需要將其轉(zhuǎn)化為裂紋前緣J積分隨歸一化裂紋前緣長(zhǎng)度s′的變化，具體轉(zhuǎn)換公式如下：

(3)

(4)

式中,s為θ弧度對(duì)應(yīng)的半橢圓形裂紋前緣弧長(zhǎng)，mm;t為積分變量;s′為歸一化裂紋前緣長(zhǎng)度;s周長(zhǎng)為半橢圓形裂紋的周長(zhǎng)，mm。

不同深長(zhǎng)比a/c下,本文與文獻(xiàn)[15]J積分分布如圖6所示?？梢钥闯?，本文得到的J積分分布與文獻(xiàn)[15]通過(guò)ABAQUS得到的計(jì)算結(jié)果吻合較好，裂紋最深點(diǎn)A處的最大J積分相對(duì)誤差僅為3.24%，在誤差接受范圍之內(nèi)，從而驗(yàn)證了本文建立的裂紋J積分計(jì)算模型的準(zhǔn)確性。

圖6 不同深長(zhǎng)比a/c下J積分分布曲線Fig.6 J integral distribution curves under the differentratios of a/c

2.2 裂紋評(píng)定

根據(jù)RCC-M《壓水堆核島機(jī)械設(shè)備設(shè)計(jì)和建造規(guī)則》，并考慮中子輻照對(duì)材料韌性的影響，得到SA508Gr.3材料對(duì)應(yīng)于出現(xiàn)失穩(wěn)擴(kuò)展和裂紋中止時(shí)的臨界應(yīng)力強(qiáng)度因子KⅠC和KⅠa分別為：

(5)

(6)

式中,T為材料溫度，℃；RTNDT為非延性轉(zhuǎn)變溫度，℃。

應(yīng)力強(qiáng)度因子計(jì)算式為：

(7)

式中,KⅠ為應(yīng)力強(qiáng)度因子，MPa·m1/2;E為彈性模量，GPa;υ為泊松比。

根據(jù)RCC-M ZG3210，對(duì)于接管圓角區(qū)域，取基準(zhǔn)缺陷為平面型半橢圓表面裂紋，深度為20 mm，長(zhǎng)度為40 mm，經(jīng)計(jì)算可得裂紋前緣最大J積分為7.45 mJ/mm2。由式(7)可得，裂紋前緣最大應(yīng)力強(qiáng)度因子KⅠ=39.13 MPa·m1/2，考慮中子輻照影響，以及在承壓熱沖擊分析時(shí)，RPV筒體材料的RTNDT為17.0 ℃[16]。由圖7可以看出，KⅠ小于材料斷裂韌性KⅠC和KⅠa值，表明在此工況及裂紋尺寸下不會(huì)發(fā)生裂紋擴(kuò)展。

圖7 承壓熱沖擊最為劇烈時(shí)的裂紋評(píng)定圖Fig.7 Evaluation diagram of the crack at the mostpressurized thermal shock

2.3 參數(shù)影響性分析

2.3.1 溫差載荷影響

為了研究?jī)?nèi)、外壁溫差Δt所產(chǎn)生的熱應(yīng)力對(duì)J積分值的影響，在熱分析中保持RPV外壁溫度為197.85 ℃不變，對(duì)內(nèi)壁分別施加116.85,136.85,156.85,176.85 ℃，即內(nèi)、外壁溫差Δt分別為81,61,41,21 ℃。同時(shí)，結(jié)構(gòu)應(yīng)力分析時(shí)保持介質(zhì)壓力為4.9 MPa不變。橢圓形表面裂紋的短半軸長(zhǎng)度為5 mm，長(zhǎng)半軸長(zhǎng)度為20 mm恒定。圖8為溫差61 ℃時(shí)裂紋前緣J積分分布云圖，其中，“1”表示半橢圓形裂紋前緣起點(diǎn)，“2”表示半橢圓形裂紋前緣終點(diǎn)。

圖8 Δt=61 ℃時(shí)裂紋前緣J積分值分布云圖Fig.8 J integral distribution nephogram of thefront edge of crack at Δt=61 ℃

圖9示出內(nèi)、外壁不同溫差下裂紋前緣J積分值分布圖。可以看出，在相同裂紋尺寸和介質(zhì)壓力情況下，裂紋前緣J積分隨著溫差的增大而增大。一般來(lái)說(shuō)，裂紋更易在最深點(diǎn)發(fā)生啟裂，即J積分值最大處。當(dāng)內(nèi)外壁溫差為81 ℃時(shí)，裂紋前緣最深點(diǎn)的J積分可達(dá)2.97 mJ/mm2，而內(nèi)外壁溫差為21 ℃時(shí)，裂紋前緣最深點(diǎn)的J積分僅為0.43 mJ/mm2，即當(dāng)溫差上升2.86倍時(shí)，裂紋前緣最深點(diǎn)的J積分增大5.91倍。

圖9 不同溫差下裂紋前緣J積分分布曲線Fig.9 J integral curves of the front edge of crack underdifferent temperature differences

圖10示出裂紋前緣最深點(diǎn)J積分隨溫差變化曲線?？梢钥闯?，裂紋前緣最深點(diǎn)J積分隨溫差增大而增大，近似呈二次曲線關(guān)系?；谒脭?shù)據(jù)進(jìn)行函數(shù)關(guān)系擬合可得：

圖10 裂紋前緣最深點(diǎn)J積分隨溫差變化曲線Fig.10 Variation curve of J integral of the deepest pointon the front edge of crack with temperature difference

J′=3.077×10-4Δt2+1.093×10-2Δt

+6.574×10-2

(8)

式中,J′為裂紋前緣最深點(diǎn)J積分，mJ/mm2。

可見(jiàn)，當(dāng)外壁處于高溫而內(nèi)壁處于低溫時(shí)，RPV處于“外加熱”狀態(tài)，溫差應(yīng)力在內(nèi)壁面產(chǎn)生的拉應(yīng)力與內(nèi)壓引起的機(jī)械拉應(yīng)力疊加使得內(nèi)壁面處綜合應(yīng)力發(fā)生惡化，在經(jīng)歷長(zhǎng)時(shí)間中子輻照、材料斷裂韌性降低后，極易使得結(jié)構(gòu)不連續(xù)處的DVI接管高應(yīng)力區(qū)的軸向裂紋發(fā)生擴(kuò)展，因此在核電站壽期末應(yīng)注意控制安注水流量，盡量避免因過(guò)大的溫差引起裂紋快速擴(kuò)展，從而發(fā)生斷裂事故。

2.3.2 裂紋尺寸影響

裂紋深長(zhǎng)比a/c表示半橢圓形裂紋的短半軸與長(zhǎng)半軸的比值，為了研究裂紋尺寸對(duì)J積分的影響，保持裂紋短半軸為5 mm不變，長(zhǎng)半軸分別取20，15，10，6.25，5 mm，即裂紋深長(zhǎng)比為1/4，1/3，1/2，4/5，1。RPV內(nèi)外壁溫差為81 ℃，介質(zhì)壓力為4.9 MPa。

圖11示出裂紋不同深長(zhǎng)比a/c下的裂紋前緣J積分分布曲線。可以看出，盡管裂紋前緣J積分仍然以裂紋最深點(diǎn)為中心，兩邊呈近似對(duì)稱分布，但當(dāng)裂紋a/c接近4/5時(shí)，裂紋前緣大部分區(qū)域近似呈現(xiàn)為直線；而當(dāng)a/c處于1/4～4/5區(qū)間時(shí)，J積分分布呈現(xiàn)中間高、兩邊低的特點(diǎn)；當(dāng)a/c處于4/5～1區(qū)間時(shí)，J積分則呈現(xiàn)出中間低、兩邊高的規(guī)律?？傮w看來(lái)，裂紋最深點(diǎn)處的J積分隨著a/c的增大而減小，而裂紋前緣起點(diǎn)和終點(diǎn)處的J積分隨著a/c的增大而增大；當(dāng)a/c相對(duì)較小時(shí)，裂紋更易在裂紋前緣最深點(diǎn)發(fā)生擴(kuò)展，且a/c越小，結(jié)構(gòu)發(fā)生破壞的概率越高；而當(dāng)a/c相對(duì)較大時(shí)，裂紋更易在裂紋前緣起點(diǎn)和終點(diǎn)處發(fā)生擴(kuò)展。

圖11 不同a/c時(shí)裂紋前緣J積分分布曲線Fig.11 J integral curves of the front edge of crackunder different ratios of a/c

2.3.3 雙裂紋相互干涉效應(yīng)

為了研究DVI接管圓角區(qū)域雙裂紋相互作用對(duì)裂紋前緣J積分的影響，保持溫差81 ℃和介質(zhì)壓力4.9 MPa載荷不變，在原有a/c=1/2的主裂紋左側(cè)l=18.40 mm處建立一個(gè)如圖4所示的附屬裂紋。附屬橢圓形表面裂紋的長(zhǎng)度保持10 mm，深度分別取10，8，5 mm，即分別取附屬裂紋深長(zhǎng)比a/c為1,4/5,1/2。

圖12示出附屬裂紋a/c=1時(shí)，單裂紋及雙裂紋相互作用下的主裂紋J積分分布。可以看出，當(dāng)存在雙裂紋相互作用時(shí)，主裂紋的裂紋前緣J積分要略高于僅存在單裂紋時(shí)的J積分。

圖12 雙裂紋相互作用時(shí)主裂紋J積分分布曲線Fig.12 J integral curves of the main crack underinteraction of double cracks

為了便于表征附屬裂紋對(duì)主裂紋的干涉效應(yīng)，參考相關(guān)文獻(xiàn)[9]，定義裂紋相互干涉作用因子η:

η=Jm/Js

(9)

式中,Jm為DVI接管圓角區(qū)域含有雙裂紋時(shí)的J積分，mJ/mm2;Js為DVI接管圓角區(qū)域含有相同尺寸單裂紋時(shí)的J積分，mJ/mm2。

圖13示出裂紋相互干涉作用因子η與附屬裂紋深度對(duì)應(yīng)關(guān)系曲線，可以看出,裂紋相互干涉作用因子η都大于1。因此，當(dāng)附屬裂紋位于主裂紋左側(cè)18.40 mm處時(shí)，附屬裂紋對(duì)主裂紋前緣的J積分起到了增強(qiáng)效應(yīng)。裂紋相互干涉作用因子η隨著附屬裂紋深度的增大而有所增加。當(dāng)附屬裂紋深度確定時(shí)，裂紋相互干涉作用因子η隨歸一化裂紋前緣長(zhǎng)度s′近似呈指數(shù)形式增加，這是因?yàn)榱鸭y前緣終點(diǎn)B與裂紋前緣起點(diǎn)C相比更加靠近附屬裂紋，因此當(dāng)雙裂紋以圖4所示形式布置時(shí)，靠近附屬裂紋的裂紋前緣終點(diǎn)B處的J積分變化更應(yīng)該引起重視。

圖13 裂紋相互干涉作用因子η與附屬裂紋深度對(duì)應(yīng)曲線Fig.13 Corresponding curves between the crack interactionfactor η and the depth of accessory crack

為了研究附屬裂紋與主裂紋之間的距離l對(duì)裂紋相互干涉作用因子η的影響，保持溫差81 ℃和介質(zhì)壓力4.9 MPa載荷不變，分別在主裂紋左側(cè)l為20.1，18.40,16.7 mm處建立一個(gè)附屬裂紋，主裂紋和附屬裂紋尺寸相同，取a=5 mm，c=10 mm。

圖14示出裂紋相互干涉作用因子η與附屬裂紋和主裂紋之間的距離對(duì)應(yīng)關(guān)系曲線。

圖14 不同雙裂紋距離下裂紋相互干涉作用因子η變化曲線Fig.14 Variation curves of crack interaction factor ηunder different double crack distances

從圖14可以看出，裂紋相互干涉作用因子η隨著附屬裂紋和主裂紋之間距離的減小而增加，即距離越小、附屬裂紋對(duì)主裂紋前緣的J積分增強(qiáng)效應(yīng)越顯著，且距離越小，裂紋相互干涉作用因子η隨歸一化裂紋前緣長(zhǎng)度的增加而增長(zhǎng)得越快，在靠近附屬裂紋的裂紋前緣終點(diǎn)B處達(dá)到最大值。

為了研究附屬裂紋角度φ對(duì)裂紋相互干涉作用因子η的影響，保持溫差81 ℃和介質(zhì)壓力4.9 MPa 載荷不變，在主裂紋左側(cè)l=16.7 mm處，分別建立角度φ為0°,30°,45°,60°,90°的附屬裂紋，主裂紋和附屬裂紋的尺寸相同，取a=5 mm,c=10 mm。

圖15示出裂紋相互干涉作用因子η與附屬裂紋角度φ的對(duì)應(yīng)關(guān)系曲線?？梢钥闯觯鸭y相互干涉作用因子η隨著附屬裂紋角度φ的減小而增大，即附屬裂紋角度越小，附屬裂紋對(duì)主裂紋前緣的J積分增強(qiáng)效應(yīng)越顯著，而當(dāng)φ=90°時(shí)，附屬裂紋對(duì)主裂紋前緣的J積分增強(qiáng)效應(yīng)幾乎可以忽略不計(jì)。

圖15 不同附屬裂紋角度φ下裂紋相互干涉作用因子η變化曲線Fig.15 Variation curves of the crack interaction factor ηunder the different angles φ of accessory crack

3 結(jié)論

針對(duì)目前RPV筒體和DVI接管連接區(qū)的裂紋擴(kuò)展特性研究不足的情況，本文基于ANSYS Workbench熱力耦合平臺(tái)和RCC-M規(guī)范，開(kāi)發(fā)了DVI接管與RPV筒體連接高應(yīng)力區(qū)含半橢圓形表面裂紋的彈塑性斷裂J積分?jǐn)?shù)值計(jì)算模型，主要結(jié)論如下。

(1)基于RCC-M《壓水堆核島機(jī)械設(shè)備設(shè)計(jì)和建造規(guī)則》，對(duì)最苛刻溫差工況下裂紋進(jìn)行了安全評(píng)定。

(2)通過(guò)多參數(shù)對(duì)比分析，表明溫差載荷、裂紋尺寸和多裂紋干涉對(duì)J積分值具有顯著影響。裂紋前緣最深點(diǎn)J積分隨溫差增大而增大，隨深長(zhǎng)比的增大而減??；裂紋相互干涉作用因子隨附屬裂紋深度的增加、附屬裂紋與主裂紋距離的減小以及附屬裂紋角度的減小而有所增加，且在靠近附屬裂紋的裂紋前緣終點(diǎn)達(dá)到最大值。

(3)需要重點(diǎn)關(guān)注RPV筒體和DVI接管連接區(qū)裂紋分布情況，進(jìn)行科學(xué)的評(píng)定。若承壓件發(fā)生熱沖擊事件時(shí)，應(yīng)控制好DVI接管區(qū)溫差，避免過(guò)高的熱應(yīng)力引起表面裂紋快速擴(kuò)展。