核電反應(yīng)堆在地震和失水事故下的結(jié)構(gòu)動力響應(yīng)分析

伍時建，尚爾濤，劉 攀，金 挺，聶照宇

(中廣核工程有限公司，廣東 深圳 518124)

核電反應(yīng)堆在地震和失水事故下的結(jié)構(gòu)動力響應(yīng)分析

伍時建，尚爾濤，劉 攀，金 挺，聶照宇

(中廣核工程有限公司，廣東 深圳 518124)

以核電廠反應(yīng)堆在地震和失水事故下的結(jié)構(gòu)動力響應(yīng)分析為工程背景，對反應(yīng)堆結(jié)構(gòu)的動力響應(yīng)分析開展了研究。分析研究形成了-套可用于新核電站反應(yīng)堆結(jié)構(gòu)設(shè)計的完整設(shè)計分析系統(tǒng)，介紹了該方法所使用的分析模型，明確了動力學(xué)分析模型中考慮的重要因素。克服了結(jié)構(gòu)復(fù)雜性、載荷多樣性，結(jié)構(gòu)非線性、材料非線性、接觸非線性和流固耦合等多重非線性因素，更加精確的模擬了反應(yīng)堆結(jié)構(gòu)，提高了計算精度。該方法對我國自主開展反應(yīng)堆事故工況下的設(shè)計和安全分析具有良好的推廣和應(yīng)用前景。

核電廠反應(yīng)堆；地震和失水事故；動力學(xué)響應(yīng)分析

作為結(jié)構(gòu)設(shè)計的基礎(chǔ)，核電廠反應(yīng)堆在地震和失水事故下的動態(tài)反應(yīng)分析研究對于保證核安全至關(guān)重要。美國ASME規(guī)范[1]、法國RCC-M規(guī)范[2]和美國核管會標(biāo)準(zhǔn)審查大綱(SRP)[3]對反應(yīng)堆抗震及防失水事故(LOCA)動力效應(yīng)有關(guān)要求，在安全停堆地震(SSE)和LOCA聯(lián)合作用下，應(yīng)能保證反應(yīng)堆安全停堆，并維持在安全停堆狀態(tài)。

由于地震和LOCA是瞬態(tài)動力激勵，而核電站反應(yīng)堆的結(jié)構(gòu)形式又十分復(fù)雜，且廣泛存在諸如間隙、摩擦、接觸等非線性因素。使得對動力激勵的確定、結(jié)構(gòu)的模型化、動力分析的方法以及相關(guān)的計算軟件系統(tǒng)開發(fā)都非常困難[4]。

美國、法國和日本在核電廠主設(shè)備設(shè)計中涉及事故工況條件下的結(jié)構(gòu)動力反應(yīng)分析時，采用了整體結(jié)構(gòu)分析與部分試驗相結(jié)合、通用程序與專用程序相結(jié)合的方法。20世紀(jì)70年代，這些國家就開始了這方面的研究工作，特別是對堆內(nèi)構(gòu)件、燃料組件、驅(qū)動線、SG傳熱管等結(jié)構(gòu)在地震加失水下的結(jié)構(gòu)動力反應(yīng)作了大量的模型、甚至是實物的地震試驗和相應(yīng)的計算分析工作，積累了大量的試驗和研究數(shù)據(jù)[5]。

自20世紀(jì)80年代開始，國內(nèi)核電廠的設(shè)計分析中開展了主設(shè)備在地震和失水事故下的反應(yīng)分析和試驗研究工作，特別是對堆內(nèi)構(gòu)件、燃料組件、控制棒驅(qū)動線與蒸汽發(fā)生器四大主設(shè)備進(jìn)行了詳細(xì)專題研究[6]。

本文在國內(nèi)外反應(yīng)堆結(jié)構(gòu)動力響應(yīng)分析方法的基礎(chǔ)上，以反應(yīng)堆在地震和LOCA下反應(yīng)堆結(jié)構(gòu)動力響應(yīng)分析為工程背景，克服了結(jié)構(gòu)復(fù)雜性、載荷多樣性，結(jié)構(gòu)非線性、材料非線性、接觸非線性和流固耦合等多重非線性因素，更加精確地模擬反應(yīng)堆結(jié)構(gòu)，提高了計算精度。

完整的設(shè)計分析流程，如圖1所示。

圖1 設(shè)計分析流程圖Fig.1 Design analysis flowchart

1 反應(yīng)堆結(jié)構(gòu)

反應(yīng)堆結(jié)構(gòu)由多個部件組成，包括壓力容器(RPV)、堆芯、堆內(nèi)構(gòu)件(RVI)、控制棒驅(qū)動機(jī)構(gòu)、堆頂結(jié)構(gòu)、RPV支承環(huán)等。RVI在RPV內(nèi)支承和固定堆芯組件(見圖2)，其作用是提供流道，支承和固定堆芯，對控制棒進(jìn)行保護(hù)和導(dǎo)向，屏蔽堆芯輻射等。下部RVI包括吊籃(CB)、堆芯圍筒組件(CS)、下支承板(LSP)和流量分配組件(HDD)等；上部RVI包括上部支承組件(上支承組件)、堆芯上板(UCP)、控制棒導(dǎo)向筒(CRGT)、支承柱和壓緊彈簧等。

圖2 RPV和RVIFig.2 RPV and RVI

反應(yīng)堆堆芯由177個尺寸相同、截面為正方形的燃料組件排列而成，位于UCP和LSP之間，外圍被CS包圍。

2 地震和LOCA分析模型

反應(yīng)堆系統(tǒng)結(jié)構(gòu)十分復(fù)雜，并且需要考慮結(jié)構(gòu)間非線性連接剛度、間隙、阻尼和流固耦合等多種非線性因素的影響，分析時必須作合適的簡化。通常將反應(yīng)堆系統(tǒng)簡化為水平和垂直兩個獨立的非線性分析模型，水平模型表示系統(tǒng)在通過容器中軸線的平面內(nèi)的平移和轉(zhuǎn)動特性，而垂直模型表示系統(tǒng)在垂直方向的動態(tài)特性。

將主回路系統(tǒng)簡化，主回路系統(tǒng)對反應(yīng)堆的作用通過一個剛度矩陣和質(zhì)量矩陣來模擬，施加在RPV進(jìn)出口接管位置。根據(jù)國內(nèi)外經(jīng)驗，堆頂結(jié)構(gòu)和CRDM對反應(yīng)堆動力響應(yīng)的影響很小，堆頂結(jié)構(gòu)和CRDM簡化為集中質(zhì)量施加在RPV上封頭頂部。

2.1 水平分析模型

使用ANSYS軟件進(jìn)行建模。對于水平分析，反應(yīng)堆系統(tǒng)簡化成由二維梁單元(或管單元)和集中質(zhì)量單元組成的模型，各部件之間通過非線性彈簧、水動力質(zhì)量矩陣、剛度矩陣或剛性梁連接，并引入阻尼、間隙，水平分析模型如圖3所示。

圖3 水平分析模型Fig.3 horizontal Analysis Model

2.1.1 模型簡化

壓力容器(RPV)和吊籃(CB)采用管單元模擬，相關(guān)部件簡化成集中質(zhì)量和轉(zhuǎn)動慣量施加在相應(yīng)位置的節(jié)點上。RPV和CB內(nèi)的冷卻劑以集中質(zhì)量的形式施加在相應(yīng)位置的節(jié)點上。

支承柱(COLUMN)、控制棒導(dǎo)向筒(CRGT)和燃料組件(FA)簡化成具有等效質(zhì)量和等效剛度的梁單元，保證前幾階頻率和模態(tài)與有限元模型相符，并考慮冷卻劑的附加質(zhì)量。

壓緊彈簧簡化為集中質(zhì)量施加于上支承法蘭對應(yīng)的節(jié)點。USP、下支承板(LSP)和堆芯上板(UCP)在水平方向視為剛體。

2.1.2 剛度計算

RPV與土建結(jié)構(gòu)之間的連接剛度包含三部分：土建剛度、RPV支承環(huán)剛度和RPV進(jìn)出口接管剛度。將RPV與土建結(jié)構(gòu)之間的連接簡化成剛度矩陣施加在反應(yīng)堆模型。在RPV進(jìn)出口接管位置引入主系統(tǒng)的剛度矩陣和質(zhì)量矩陣來考慮主系統(tǒng)對反應(yīng)堆的影響。

RPV法蘭與USP法蘭之間和RPV法蘭與CB法蘭之間的非線性碰撞剛度和轉(zhuǎn)動剛度，USP法蘭與CB法蘭之間壓緊彈簧的垂直和轉(zhuǎn)動剛度，通過有限元計算得到(見圖4)。

圖4 RPV與USP法蘭，CB法蘭連接示意圖Fig.4 RPV and USP flanges，GB flange connection diagram

CB出口管嘴和RPV出口接管內(nèi)凸臺之間的剛度，LSP與RPV之間徑向支承鍵的剛度，CB與UCP之間的剛度和CB與CS之間的剛度也通過有限元計算得到。

RVI之間以及RVI與RPV之間的間隙均采用反應(yīng)堆穩(wěn)態(tài)運行間隙。

2.2 垂直分析模型

對于垂直分析，反應(yīng)堆系統(tǒng)簡化成由彈簧單元和集中質(zhì)量單元組成的模型，垂直分析模型如圖5所示。

圖5 垂直分析模型Fig.5 Vertical Analysis Model

將下部堆內(nèi)構(gòu)件、燃料組件、上部堆內(nèi)構(gòu)件分成中部和邊緣部分(按下支承板的二階振型進(jìn)行劃分)，這樣可以更好地描述上部和下部堆內(nèi)構(gòu)件的變形及其對堆芯的影響，并提供更為詳細(xì)的地震和LOCA載荷分布。

3 動力學(xué)分析模型中考慮的因素

3.1 阻尼

采用等效黏性阻尼來體現(xiàn)系統(tǒng)中的阻尼效應(yīng)，在系統(tǒng)運動微分方程中，假設(shè)Rayleigh阻尼矩陣如下：

[C]=α[M]+β[K]

式中：[M]為質(zhì)量矩陣；[K]為剛度矩陣；α、β為Rayleigh阻尼常數(shù)，可通過以下方程求解：

ξi=α/(2ωi)+βωi/2

式中：ξi為結(jié)構(gòu)的阻尼比，ωi為結(jié)構(gòu)頻率(圓頻率)，一般可選取結(jié)構(gòu)的兩個最重要的頻率值來求解α、β。反應(yīng)堆各構(gòu)件的重要頻率多在5～40Hz之間，因此選擇5Hz和40Hz這兩個頻率來計算Rayleigh阻尼常數(shù)。對于反應(yīng)堆地震和LOCA分析，阻尼比取4%，計算α、β。阻尼比與頻率的關(guān)系曲線如圖6所示。

圖6 Rayleigh阻尼Fig.6 Rayleigh damping

3.2 流固耦合

RPV中充滿冷卻劑流體，RVI和燃料組件完全浸沒于流體中，需要考慮流固耦合作用。

RPV與CB之間、上支承組件裙筒與CB筒體之間、燃料組件與CS之間的流固耦合作用采用Fluid38單元模擬，支承柱和CRGT的流固耦合作用以附加水質(zhì)量的形式模擬，反應(yīng)堆中其他區(qū)域的冷卻劑質(zhì)量以集中質(zhì)量的形式施加在RPV和CB節(jié)點上。

Fluid38單元形式如圖7所示。

圖7 流固耦合單元Fig.7 Fluid-structure coupling element

上支承組件裙筒與CB筒體之間的環(huán)形水隙長度較短，流固耦合的附加水質(zhì)量根據(jù)Fritz給出的公式進(jìn)行計算[7]：

燃料組件相互之間的流固耦合作用很小，只考慮燃料組件與圍板之間的流固耦合作用。

4 地震和LOCA作用下的結(jié)構(gòu)響應(yīng)分析

4.1 地震和LOCA輸入

地震分析中，采用土建結(jié)構(gòu)地震分析得到的RPV支承處的地震加速度時程作為地震載荷輸入[8]，如圖8所示。

圖8 RPV支承處地震加速度Fig.8 Seismic acceleration of RPV support

LOCA分析中，考慮余熱排除出系統(tǒng)接管的雙端斷裂(RHR)和安全注入系統(tǒng)接管的雙端斷裂(SIL)兩種破口形式下的載荷。

水平方向的LOCA水力載荷如圖9所示。

圖9 水平方向LOCA水力載荷Fig.9 horizontal LOCA hydraulic load

考慮RVI壓緊彈簧預(yù)緊力和燃料組件壓緊彈簧的預(yù)緊力作用。

4.2 地震和LOCA分析結(jié)果

使用ANSYS軟件，采用直接積分方法，進(jìn)行瞬態(tài)動力學(xué)分析，最終得到反應(yīng)堆內(nèi)各主要部位的位移和作用力。

計算得到用于堆頂結(jié)構(gòu)和CRDM地震分析的RPV上封頭的位移時程如圖10所示。

圖10 RPV上封頭位移(三個方向)Fig.10 RPV head displacement(three directions)

5 結(jié)論

本文介紹了華龍一號反應(yīng)堆在地震和LOCA下的結(jié)構(gòu)動力響應(yīng)分析研究。在國內(nèi)外反應(yīng)堆結(jié)構(gòu)動力響應(yīng)分析方法的基礎(chǔ)上引進(jìn)消化吸收，更加精確模擬反應(yīng)堆結(jié)構(gòu)，提高了計算精度。已應(yīng)用在華龍一號反應(yīng)堆結(jié)構(gòu)設(shè)計和安全評審中，對我國自主開展百萬千瓦級先進(jìn)壓水堆核電廠反應(yīng)堆在地震和失水事故下的設(shè)計和安全分析具有良好的推廣和應(yīng)用前景。

[1] ASME Boiler and Pressure Vessel Codel[S].American Society of Mechanical Engineers，1997.

[2] Design and construction rules for mechanical components of PWR nuclear islands[S].AFCEN，2007.

[3] Standard Review Plan for the Review of Safety Analysis Reports for Nuclear Power Plants：LWR Edition[S]. NUREG-0800.

[4] 毛慶. 核反應(yīng)堆地震和失水事故下結(jié)構(gòu)動力分析研究[J].核動力工程，2002，23：93-98.

[5] 姚偉達(dá). 核電廠主設(shè)備在地震加失水事故下的結(jié)構(gòu)反應(yīng)分析研究[J].核動力工程，2002，23：135-138.

[6] 劉文進(jìn). 反應(yīng)堆結(jié)構(gòu)三維非線性抗震分析研究[J].核動力工程，2013，34：24-26.

[7] Fritz R. J. The Effect of Liquids on the Dynamic Motions of Immersed Solids[J]，ASME Journal of Engineering for Industry，1972.

[8] 姜乃斌. 反應(yīng)堆冷卻劑系統(tǒng)地震分析的地面時程輸入方法研究[C]//第14屆全國反應(yīng)堆結(jié)構(gòu)力學(xué)會議論文集，2006.

DynamicResponseAnalysisofNuclearPowerReactorunderSeismicandLOCA

WUShi-jian，SHANGEr-tao，LIUPan，JINTing，NIEZhao-yu

(China Guangdong Nuclear Power Engineering Co. Ltd.，Shenzhen of Guangdong Prov. 518124，China)

Based on the structural dynamic response analysis of a nuclear power plant reactor under earthquake and LOCA, the dynamic response analysis of reactor structure is carried out. A complete design and analysis system for reactor structure design of new nuclear power plant has been formed. The analytical model used in this method has been introduced, and the important factors considered in the dynamic analysis model have been defined. It overcomes many factors, such as structural complexity, load diversity, structural nonlinearity, material nonlinearity, contact nonlinearity and fluid-structure coupling. It simulates the reactor structure more accurately and improves the calculation accuracy. This method has a good prospect of popularization and application in the accident condition design and safety analysis of China's independent reactor.

Nuclear reactor；Seismic and LOCA conditions；Dynamic response analysis

2017-01-19

伍時建(1986—)，男，四川人，工程師，碩士，現(xiàn)主要從事核電站設(shè)計工作

O34

A

0258-0918(2017)06-0943-05