CEFR堆芯燃料組件抗震試驗方法研究

高付海，莫亞飛，王魯波

(中國原子能科學(xué)研究院 反應(yīng)堆工程技術(shù)研究部，北京 102413)

自2011年中國實驗快堆(CEFR)首次并網(wǎng)發(fā)電以來，為充分發(fā)揮其實驗平臺作用，CEFR堆芯燃料的設(shè)計和制造國產(chǎn)化逐步被提上日程。國產(chǎn)的CEFR燃料組件在材料選擇、結(jié)構(gòu)型式、制造裝配和焊接工藝上均與CEFR堆芯目前使用的組件存在差異，如結(jié)構(gòu)材料采用國產(chǎn)316Ti不銹鋼、燃料棒下端塞與小柵板配合改為蒜頭型設(shè)計方案等，國產(chǎn)燃料組件在地震下能否保證結(jié)構(gòu)完整性和氣密性，可通過抗震試驗的方式予以論證。針對快堆堆芯地震動力學(xué)響應(yīng)，國外特別是法國和日本進行了大量的試驗和數(shù)值計算研究工作，建立地震模擬振動臺，開展堆芯抗震試驗，開發(fā)專用計算程序(如法國CEA的CASTEM和日本PNC的FINAS)，研究目的在于準(zhǔn)確預(yù)測和正確評價地震作用下堆芯的結(jié)構(gòu)動力學(xué)響應(yīng)[1-6]。俄羅斯為了驗證CEFR組件的動態(tài)可靠性，開展了組件模擬件的耐振試驗考驗研究[7]。國內(nèi)中國原子能科學(xué)研究院[8-12]和華北電力大學(xué)[13-14]組織開展了持續(xù)的堆芯抗震研究工作，計算了CEFR堆芯在地震載荷下的動力學(xué)響應(yīng)，掌握了快堆堆芯抗震的數(shù)值計算方法。在試驗研究方面中國原子能科學(xué)研究院委托清華大學(xué)測試了CEFR組件模擬件在水中的固有頻率和阻尼比等動力學(xué)特性，積累了良好的工作基礎(chǔ)。

本文基于俄羅斯的組件耐振試驗方案，結(jié)合國內(nèi)機械設(shè)備抗震試驗法規(guī)標(biāo)準(zhǔn)要求，針對組件的結(jié)構(gòu)、堆內(nèi)安裝和變形特點，提出一套新的組件抗震試驗方法，新方法可保證對組件的抗震性能考驗是全面的，且更加符合組件在堆芯處的實際變形狀態(tài)。

1 俄羅斯CEFR組件耐振試驗臺架用于組件抗震試驗的適用性分析

組件耐振試驗的目的主要是考察主泵長時機械振動載荷對組件結(jié)構(gòu)性能的影響，載荷雖不同于地震，但其試驗方案對制定抗震試驗方案仍具有借鑒意義。試驗使用的臺架如圖1所示[1]，臺架模擬組件在堆芯中的安裝條件，試驗基本思想是取堆芯中由1根組件和其周圍介質(zhì)構(gòu)成的1個單元為試驗對象，以單根組件試驗來代替堆芯多排組件試驗。組件凸臺部位的定位通過螺釘來實現(xiàn)，保留燃料組件上、中兩個凸臺處間隙。臺架柱形容器內(nèi)充滿水，模擬液鈉的影響。

圖1所示的試驗臺架基本能反映組件在堆芯內(nèi)的實際安裝條件(圖2)，在振動載荷下組件管腳隨臺架一起運動，定位凸臺部位處受到臺架容器上設(shè)置的螺釘?shù)募s束，以模擬組件受到的相鄰組件凸臺約束作用。若該臺架方案直接用于組件抗震試驗，其與地震作用下處于堆芯中組件的響應(yīng)存在差異：模擬相鄰組件凸臺作用的螺釘是完全剛性的，會使得試驗組件受到的沖擊力高，頂端位移響應(yīng)卻較小；實際情況是相鄰組件并不是完全剛性的，是可變形的柔性體。因此可從定性角度判斷俄羅斯試驗臺架設(shè)計對組件約束得過硬，使得組件受到更為苛刻的高頻沖擊，以造成較大的機械振動疲勞損傷，實現(xiàn)耐振試驗的目的，但這與抗震試驗?zāi)康牟煌耆韧?/p>

圖1 試驗臺架與試驗件Fig.1 Test set-up and subassembly

圖2 堆芯組件的安裝狀態(tài)[8] Fig.2 Installation state of subassembly in core[8]

2 組件抗震試驗新方法

2.1 新方法的提出

按照文獻[15]，試驗程序是根據(jù)設(shè)備安裝位置處的運行基準(zhǔn)地震(SL1)和安全停堆地震(SL2)樓層響應(yīng)譜，采用單向、雙向或三向人工時程對設(shè)備先做5次SL1地震模擬試驗，再做1次SL2地震模擬試驗?？茖W(xué)合理的組件抗震試驗方法須處理好兩個基本問題：1) 在滿足試驗要求的前提下，從降低試驗難度角度考慮，可否以單向地震試驗替代三向試驗；2) 能否設(shè)計一合理的試驗臺架，準(zhǔn)確有效地模擬地震作用下堆芯組件的安裝條件和碰撞接觸相互作用狀態(tài)。

關(guān)于能否采用單向激振進行抗震試驗，《核設(shè)備抗震鑒定試驗指南》中規(guī)定：當(dāng)三向中的每兩向相互耦合很小，或無耦合，單向試驗是適用的；若是單向試驗，應(yīng)分別對x、y、z軸進行。組件的結(jié)構(gòu)基本對稱，存在3個正交軸，各正交方向在地震下的響應(yīng)基本上無耦合，符合單向試驗的條件。組件的縱向剛度偏大，響應(yīng)可忽略不計。組件響應(yīng)以水平兩個方向響應(yīng)為主，組件之間的碰撞主要體現(xiàn)在水平方向，鑒于此，組件的抗震試驗可忽略縱向試驗，只進行水平方向的試驗。

關(guān)于試驗臺架設(shè)計，嚴(yán)格講為了準(zhǔn)確模擬堆芯地震響應(yīng)組件的抗震試驗臺架應(yīng)能容納全堆芯的所有組件，但這對試驗費用和振動試驗臺提出了不必要的過高要求。俄羅斯采用的以單組件振動試驗來替代全堆芯試驗的臺架設(shè)計思路是可借鑒的，它大幅降低了全堆芯組件同時振動試驗的難度，且基本能反映出振動時組件的狀態(tài)，但在對可能出現(xiàn)碰撞作用的相鄰組件的處理上采用螺釘近似，該處理過于剛性，不能模擬出相鄰組件的柔性特征，試驗結(jié)果可能會對組件受到的碰撞沖擊力估計偏高，而對組件產(chǎn)生的變形和應(yīng)力水平估計偏低。為敘述方便，將該種類型臺架稱為剛性臺架。為了模擬出相鄰組件的柔性，也可將相鄰組件簡化成一柔性元件，其剛度約等于相鄰組件的抗彎剛度，組件間隙同樣保留，將該種臺架稱為柔性臺架。一個科學(xué)合理的抗震試驗方法必須兼顧考慮以上兩種臺架方案對組件的不同影響，本文提出的組件抗震試驗新方法即同時考慮以上兩種臺架試驗。

2.2 抗震試驗新方法的計算分析

使用專業(yè)計算軟件CASTEM，分別建立安裝于剛性臺架上的單組件(方案1)、安裝于柔性臺架上的單組件(方案2)及接近真實狀態(tài)的堆芯單排組件(方案3)抗震計算模型，分析計算地震時3種方案下組件的位移和沖擊力響應(yīng)，通過計算結(jié)果的對比驗證單組件抗震試驗方法的可行性。

1) 抗震計算模型

在地震情況下燃料組件主要表現(xiàn)出梁的特性，因此計算模型中采用變截面梁模擬單根組件，采用虛擬的液體附加質(zhì)量考慮液鈉與組件的流固耦合作用。組件下部管腳與管座之間存在間隙，使用等效彈簧元模擬，以等效的線性邊界條件代替非線性的間隙條件，等效彈簧的剛度均為8.8×105N/m。組件凸臺與其他組件凸臺之間、組件凸臺與外緣圍板之間在地震時可能發(fā)生碰撞，均采用間隙-彈簧單元模擬。

圖3 地震作用下的堆芯加速度時程輸入Fig.3 Core acceleration history input for seismic action

地震時堆芯柵板處加速度時程如圖3所示，時間步長為0.01 s，共3 200步，即總時間為32 s。

(1) 方案1，單組件在剛性臺架上地震下的位移和沖擊力響應(yīng)

組件與臺架之間的接觸剛性碰撞通過間隙彈簧單元模擬，如圖4所示，間隙值與實際一致，碰撞彈簧認(rèn)為是剛性的，所以要求其剛度取值較大，計算設(shè)定為3.0×109N/m。

圖4 單組件處于剛性臺架上的計算模型Fig.4 Computation model for single subassembly at rigid set-up

(2) 方案2，單組件在柔性臺架上SL1地震下的位移和沖擊力響應(yīng)

圖5 單組件處于柔性臺架上的計算模型Fig.5 Computation model for single subassembly at flexible set-up

保守考慮，組件凸臺高度處可能發(fā)生的最大單側(cè)位移取堆芯單排組件的最大總間隙為9.7 mm。抗震計算模型通過兩組間隙-彈簧單元模擬組件之間的相互作用，如圖5所示。一組模擬組件的最大單側(cè)位移限位，間隙為9.7 mm，彈簧碰撞剛度為3.0×109N/m；另一組模擬相鄰組件的柔性碰撞，間隙值跟實際一致，彈簧剛度設(shè)定為相鄰組件的抗彎剛度，有限元軟件計算出上部凸臺彈簧剛度為5.4×103N/m，中部凸臺彈簧剛度為1.1×104N/m。

(3) 方案3，單排組件在地震下的位移和沖擊力響應(yīng)

組件與組件之間的接觸剛性碰撞通過間隙彈簧單元模擬，如圖6所示，間隙值與實際一致，碰撞彈簧剛度取墊塊局部結(jié)構(gòu)剛度，均設(shè)定為3.0×107N/m。

圖6 單排組件計算模型Fig.6 Computation model for single array of subassemblies

2) 抗震計算結(jié)果對比分析

圖7 單根組件處于剛性臺架上的 沖擊力響應(yīng)與單排組件沖擊力響應(yīng)對比Fig.7 Impact force response result comparison of single subassembly at rigid set-up and single array of subassemblies

對于單組件剛性臺架計算模型，組件在定位凸臺處由于剛性碰撞作用，與單排組件計算模型獲得的計算結(jié)果相比，組件在剛性臺架上凸臺受到的碰撞沖擊力更大(圖7)，碰撞力基本可包絡(luò)單排組件的結(jié)果；但組件在剛性臺架上產(chǎn)生的頭部變形由于受到限制，變形遠小于單排組件的計算結(jié)果(圖8)；由此可判斷剛性臺架使得組件在地震時受到更大的局部高頻沖擊載荷作用。對于單組件柔性臺架計算模型，與單排組件計算模型獲得的計算結(jié)果相比：組件在柔性臺架上凸臺受到的碰撞沖擊力遠小于單排組件的結(jié)果(圖9)，但組件頭部產(chǎn)生的變形基本包絡(luò)單排組件的結(jié)果(圖10)，由此可判斷柔性臺架必然使得組件受到更多的低頻振動，而低頻振動更易與柵板聯(lián)箱處作為地震輸入作用到組件的反應(yīng)譜的峰值頻段重疊、引發(fā)組件結(jié)構(gòu)一次應(yīng)力破壞。單排組件的計算結(jié)果顯示：堆芯中組件的實際響應(yīng)既有高頻沖擊，又有低頻振動。以上結(jié)果表明如果對單組件在剛性臺架和柔性臺架依次分別進行抗震試驗，其試驗結(jié)果既偏保守，又有代表性，可保證在相同地震輸入下單組件的變形(應(yīng)力)和沖擊振動包絡(luò)堆芯組件群中的組件響應(yīng)。

圖8 單根組件處于剛性臺架上的 位移響應(yīng)與單排組件位移響應(yīng)對比Fig.8 Displacement response result comparison of single subassembly at rigid set-up and single array of subassemblies

圖9 單根組件處于柔性臺架上的 沖擊力響應(yīng)與單排組件沖擊力響應(yīng)對比Fig.9 Impact force response result comparison of single subassembly at flexible set-up and single array of subassemblies

圖10 單根組件處于柔性臺架上的 位移響應(yīng)與單排組件位移響應(yīng)對比Fig.10 Displacement response result comparison of single subassembly at flexible set-up and single array of subassemblies

3 結(jié)論

本文以組件耐振試驗方案為基礎(chǔ)，基于有限元計算結(jié)果，結(jié)合國內(nèi)抗震試驗法規(guī)要求，保守考慮提出了一種新的組件抗震試驗方法，得到如下結(jié)論。

1) 俄羅斯組件耐振試驗雖不同于抗震試驗，但其試驗臺架設(shè)計方案可借鑒用于組件抗震試驗。俄羅斯的組件耐振試驗臺架設(shè)計對組件約束過硬，使得組件受到的偏苛刻的高頻沖擊，但對于頻率成分以低頻為主的地震激勵，使得組件的位移響應(yīng)偏小，低頻應(yīng)力偏低，從應(yīng)力強度角度判斷基于該臺架方案的抗震試驗結(jié)果不保守。

2) 按照《核設(shè)備抗震鑒定試驗指南》中的試驗要求，分別依次開展單根組件位于本文提出的剛性臺架和柔性臺架上的開展抗震試驗，其結(jié)果是偏保守的，可保證在相同地震輸入下單組件位移響應(yīng)、沖擊響應(yīng)基本可包絡(luò)位于堆芯組件群中的各組件的響應(yīng)。因此本文提出的基于單根組件的抗震試驗新方法要求組件模擬件依次在剛性臺架和柔性臺架上完成抗震試驗，試驗后檢查組件結(jié)構(gòu)完整性和氣密性。