莊初立， 張永山， 汪大洋

(廣州大學(xué) 土木工程學(xué)院， 廣州 510006)

超設(shè)計(jì)地震作用下核島結(jié)構(gòu)三維震動(dòng)響應(yīng)與控制研究

莊初立， 張永山， 汪大洋

(廣州大學(xué) 土木工程學(xué)院， 廣州 510006)

建立了AP1000核島結(jié)構(gòu)(屏蔽廠房和鋼制安全殼)的ANSYS有限元模型，將ANSYS模態(tài)分析結(jié)果與既有研究成果和Midas Gen計(jì)算結(jié)果對(duì)比，驗(yàn)證了核島結(jié)構(gòu)有限元模型的正確性，基于此模型開(kāi)展安全停堆和超設(shè)計(jì)地震下的核島結(jié)構(gòu)三維震動(dòng)響應(yīng)研究，引入基礎(chǔ)隔震控制技術(shù)，研究隔震核島結(jié)構(gòu)在超設(shè)計(jì)地震作用下的震動(dòng)控制效果。研究結(jié)果表明，AP1000核島結(jié)構(gòu)在超設(shè)計(jì)地震作用下，水平和豎向加速度響應(yīng)劇烈，頂點(diǎn)最大水平和豎向加速度分別為安全停堆地震下的3.31倍和3.33倍，屏蔽廠房結(jié)構(gòu)應(yīng)力超過(guò)材料容許應(yīng)力進(jìn)入非線性狀態(tài)。采用基礎(chǔ)隔震控制后核島結(jié)構(gòu)在超設(shè)計(jì)地震下的相對(duì)水平加速度、上部結(jié)構(gòu)相對(duì)水平位移、構(gòu)件內(nèi)力和基底剪力等地震反應(yīng)大幅降低，基底剪力減小率均達(dá)78%以上，隔震控制效果理想，基礎(chǔ)隔震技術(shù)能有效地解決超設(shè)計(jì)地震下核島結(jié)構(gòu)的地震安全問(wèn)題。

核島結(jié)構(gòu)；三維震動(dòng)；超設(shè)計(jì)地震；基礎(chǔ)隔震

核電建設(shè)是我國(guó)能源建設(shè)的一項(xiàng)重要政策，隨著《核電安全規(guī)劃(2011-2020年)》和《核電中長(zhǎng)期發(fā)展規(guī)劃(2011-2020年)》相繼實(shí)施，我國(guó)核電站建設(shè)迎來(lái)新一輪高峰期。然而，地震不確定性和隨機(jī)性給核電站帶來(lái)安全威脅，2007年日本新潟縣發(fā)生強(qiáng)地震，柏崎刈羽核電站遭受超過(guò)原設(shè)計(jì)基準(zhǔn)的地震動(dòng)，其中2號(hào)機(jī)組基底記錄到的水平加速度最大值為606 gal，超過(guò)設(shè)計(jì)值的3.6倍，對(duì)廠址上的非安全相關(guān)建筑物和構(gòu)筑物造成了較為嚴(yán)重的破壞[1-2]。我國(guó)大部分區(qū)域處于地震區(qū)，近年來(lái)強(qiáng)震頻發(fā)給核電站帶來(lái)安全隱患，假想在類似汶川和玉樹等原屬于低烈度區(qū)(0.1～0.15g)上建有核電站，汶川地震時(shí)震中區(qū)水平加速度記錄達(dá)0.4～0.8g，那么核電站將遭遇兩倍于安全停堆地震的作用[3]，因此如何確保核電站在超設(shè)計(jì)地震作用時(shí)的安全，成為亟待解決的問(wèn)題。

目前，三向地震動(dòng)下核島結(jié)構(gòu)的動(dòng)力響應(yīng)問(wèn)題鮮有研究，周國(guó)良等[4]分析安全停堆地震下核島廠房動(dòng)力響應(yīng)時(shí)發(fā)現(xiàn)反應(yīng)譜分析結(jié)果比三向地震時(shí)程分析結(jié)果偏小，建議采用時(shí)程法進(jìn)行校核驗(yàn)算。賀秋梅等[5]對(duì)某高溫氣冷堆核電廠進(jìn)行了安全停堆三向地震時(shí)程分析，得出應(yīng)重點(diǎn)關(guān)注豎向地震動(dòng)對(duì)結(jié)構(gòu)地震反應(yīng)的影響。目前研究關(guān)注于安全停堆地震下核島結(jié)構(gòu)的三向地震反應(yīng)，而對(duì)于超設(shè)計(jì)地震下核島結(jié)構(gòu)的三向地震響應(yīng)研究并沒(méi)有涉及。超設(shè)計(jì)強(qiáng)震作用下核電站地震安全問(wèn)題的解決方案研究較為薄弱，而結(jié)構(gòu)隔震技術(shù)能有效減輕地震作用，可應(yīng)用于核電結(jié)構(gòu)提高其抗震安全性[6-7]，保護(hù)核電站在遭遇強(qiáng)震時(shí)結(jié)構(gòu)構(gòu)件、管道、設(shè)備等的安全性，避免發(fā)生核事故，本文以AP1000核島結(jié)構(gòu)為研究對(duì)象，研究核島結(jié)構(gòu)在安全停堆和超設(shè)計(jì)地震作用下的三向地震響應(yīng)規(guī)律，引入基礎(chǔ)隔震技術(shù)，解決突發(fā)強(qiáng)震時(shí)核電站的地震安全問(wèn)題。

1 AP1000核島結(jié)構(gòu)有限元建模與驗(yàn)證

1.1 模型的建立

AP1000是美國(guó)西屋電氣公司采用非能力安全系統(tǒng)的概念來(lái)改進(jìn)核電站的安全，并且簡(jiǎn)化設(shè)計(jì)來(lái)提高核電站的安全、建造、運(yùn)行和維護(hù)，屬于第三代核電技術(shù)。AP1000的反應(yīng)堆安全殼結(jié)構(gòu)(如圖1所示)由鋼制安全殼和屏蔽廠房組成，屏蔽廠房為鋼筋混凝土結(jié)構(gòu)，采用C45級(jí)混凝土，鋼制安全殼為薄殼結(jié)構(gòu)采用SA738B級(jí)高強(qiáng)鋼材。目前AP1000在美國(guó)尚無(wú)運(yùn)行業(yè)績(jī)，而我國(guó)在建的三門核電站和海陽(yáng)核電站均采用AP1000核電技術(shù)。

圖1 AP1000核島結(jié)構(gòu)Fig.1 AP1000 nuclear island structure

采用ANSYS軟件建立AP1000核島結(jié)構(gòu)精細(xì)化有限元模型，假定AP1000核島結(jié)構(gòu)位于硬巖石廠址，因此土壤與結(jié)構(gòu)的相互作用可忽略不計(jì)[8]，核島結(jié)構(gòu)采用固結(jié)的方式與基礎(chǔ)相連，底部邊界條件采用固定約束的方式。采用shell181單元，該單元適用于模擬薄殼至中等厚度殼結(jié)構(gòu)，有限元模型共劃分16 112個(gè)單元，13 868個(gè)節(jié)點(diǎn)，在屏蔽廠房進(jìn)氣孔處劃分500 mm，安全殼、屏蔽廠房和筏板基礎(chǔ)的劃分尺寸分別為1 600 mm、1 400 mm和800 mm，有限元模型如圖2所示。安全殼鋼材采用SA738型B級(jí)鋼，本構(gòu)采用理想彈塑性模型，屈服應(yīng)力為415 MPa，鋼材的彈性模量E=200 000 MPa，泊松比v=0.3。屏蔽廠房采用C45混凝土，泊松比v=0.2，混凝土本構(gòu)采用Rüsch建議的模型[9]，如圖3所示。

當(dāng)ε≤ε0時(shí)(上升段)

(1)

當(dāng)ε0≤ε≤εu時(shí)(下降段)

σ=fc

(2)

式中：fc為混凝土軸心抗壓強(qiáng)度，取fc=29.6 N/mm2；ε0為相應(yīng)于峰值應(yīng)力時(shí)的應(yīng)變，取ε0=0.002；εu為極限壓應(yīng)變，取εu=0.003 5。

圖2 AP1000核島結(jié)構(gòu)有限元模型Fig.2 Finite element model of nuclear island structure

圖3 混凝土應(yīng)力-應(yīng)變曲線Fig.3 Stress-strain curve of normal concrete

1.2 模型的驗(yàn)證

為保證分析模型的正確性，將AP1000核島結(jié)構(gòu)ANSYS模型的模態(tài)分析結(jié)果與現(xiàn)有核電站核島結(jié)構(gòu)研究文獻(xiàn)進(jìn)行對(duì)比，如表1所示，從表1中可以看出，核電站安全殼結(jié)構(gòu)作為核反應(yīng)堆的最后一道屏障，安全性能十分重要，混凝土壁厚基本為1.0 m，安全殼結(jié)構(gòu)均具有較大的剛度。本文的分析模型第一階自振頻率為4.04 Hz，與表1中既有文獻(xiàn)的頻率數(shù)據(jù)較為接近，AP1000核島結(jié)構(gòu)高度比一般的核島結(jié)構(gòu)要高一些，所以其頻率與既有研究成果相比有所降低，核島結(jié)構(gòu)第一振型均為梁式平動(dòng)。此外，還采用第二個(gè)有限元軟件Midas Gen進(jìn)行校核驗(yàn)證，Midas Gen計(jì)算所得第一階頻率為3.85 Hz，與ANSYS結(jié)果誤差僅為4.70%，圖4給出了在安全停堆地震作用下兩個(gè)軟件的時(shí)程計(jì)算結(jié)果，從圖中可知其位移和加速度結(jié)果相當(dāng)吻合。因此，本文所建立的AP1000核島結(jié)構(gòu)有限元模型較為合理。

表1 核島結(jié)構(gòu)模型計(jì)算結(jié)果與既有結(jié)果對(duì)比

(a)安全殼頂點(diǎn)X向位移對(duì)比

(b)安全殼頂點(diǎn)X向加速度對(duì)比圖4 ANSYS與Midas Gen計(jì)算結(jié)果對(duì)比Fig.4 Comparison of calculation results between ANSYS and Midas

2 AP1000核電站地震動(dòng)輸入

AP1000核電站抗震I類SSCs必須能夠承受安全停堆地震，并保證在地震發(fā)生時(shí)執(zhí)行其預(yù)定安全功能的能力，根據(jù)美國(guó)核管委管理導(dǎo)則RG1.29的分類，AP1000抗震I類構(gòu)筑物包括屏蔽廠房、鋼安全殼及安全殼內(nèi)部結(jié)構(gòu)和輔助廠房，故而本文分析的AP1000核島結(jié)構(gòu)必須要能夠承受安全停堆地震(0.3g)。采用AP1000設(shè)計(jì)反應(yīng)譜為目標(biāo)譜擬合人工時(shí)程，AP1000核電站設(shè)計(jì)反應(yīng)譜，是基于美國(guó)核電管理導(dǎo)則RG1.60譜并研究高頻區(qū)域的放大影響而經(jīng)評(píng)價(jià)后確定的。

AP1000核電廠用于抗震I類結(jié)構(gòu)的地震動(dòng)為三個(gè)相互獨(dú)立正交、統(tǒng)計(jì)獨(dú)立的的人工合成加速度時(shí)程，擬合的設(shè)計(jì)時(shí)程其總持時(shí)不小于20 s，包括一個(gè)持續(xù)時(shí)間超過(guò)6 s，相對(duì)穩(wěn)定的、強(qiáng)烈的地面運(yùn)動(dòng)。所采用的設(shè)計(jì)時(shí)程是根據(jù)Taft實(shí)際地震記錄的綜合時(shí)程曲線進(jìn)行修改而生成的，采用Abrahamson等[15-16]提出的小波算法對(duì)實(shí)際強(qiáng)震記錄進(jìn)行修正，其方程函數(shù)如式(3)

(3)

Δt處的漸變余弦小波函數(shù)為

(4)

ψj(f)使用三線性模型來(lái)確定

(5)

式中,各參數(shù)含義詳見(jiàn)文獻(xiàn)[16]。

以AP1000標(biāo)準(zhǔn)反應(yīng)譜為目標(biāo)譜對(duì)Taft加速度時(shí)程進(jìn)行修正生成三條人工時(shí)程，計(jì)算中阻尼比采用0.05，總持時(shí)為30 s，步長(zhǎng)為0.01 s。此外還選取了Kery County和Qian An波作為時(shí)程分析輸入，人工時(shí)程和天然波與5%阻尼比AP1000設(shè)計(jì)反應(yīng)譜吻合較好，地震波的頻譜曲線如圖5所示。

(a)水平地震波譜曲線

(b)豎向地震波譜曲線圖5 地震波譜曲線Fig.5 The spectrum curves of earthquake waves

對(duì)核島結(jié)構(gòu)進(jìn)行地震時(shí)程分析時(shí)設(shè)置了0.2g、0.3g、0.4g、0.6g、0.8g、1.0g共六組加速度峰值工況，進(jìn)行三向地震時(shí)程分析，結(jié)構(gòu)阻尼采用常用的瑞利阻尼模型，其中阻尼比取0.05[17]。考慮到實(shí)際所記錄到的強(qiáng)震記錄其豎向加速度峰值有接近或大于水平加速度峰值的情況[18]，故本文分析時(shí)豎向峰值地面加速度保守地假定與水平方向取值相同。

3 AP1000核島結(jié)構(gòu)三維地震響應(yīng)研究

3.1 加速度響應(yīng)分析

圖6給出了在安全停堆地震和超設(shè)計(jì)地震作用下沿屏蔽廠房高度分布的最大水平和豎向加速度峰值。由圖6知，三組地震波作用下的屏蔽廠房水平加速度沿結(jié)構(gòu)外壁高度變化規(guī)律一致，最大水平加速度發(fā)生在屏蔽廠房頂點(diǎn)處，隨著地震輸入強(qiáng)度的增加，屏蔽廠房水平加速度顯著增大，超設(shè)計(jì)地震作用下屏蔽廠房頂點(diǎn)加速度過(guò)大，其中1.0g工況時(shí)最大水平加速度達(dá)52.54 m/s2(人工波)，為安全停堆地震下的3.31倍。超設(shè)計(jì)地震作用下核島結(jié)構(gòu)遭遇的動(dòng)力響應(yīng)為設(shè)計(jì)地震動(dòng)下的數(shù)倍，威脅結(jié)構(gòu)構(gòu)件及內(nèi)部設(shè)備管道的使用狀態(tài)，應(yīng)當(dāng)引起足夠的重視。

三向地震作用下屏蔽廠房沿外壁高度分布的豎向加速度在結(jié)構(gòu)高度約60 m處突變，原因在于屏蔽廠房約60 m處由圓錐形混凝土壁支承著上部達(dá)2 800 m3容積的冷卻水箱，結(jié)構(gòu)豎向剛度有所減小且上部冷卻水箱質(zhì)量較大，造成豎向加速度突變，放大了豎向加速度響應(yīng)，Kern County、QianAn和人工波作用下，頂點(diǎn)處的豎向加速度峰值放大為屏蔽廠房60 m處的倍數(shù)分別為2.94、2.19和1.70。超設(shè)計(jì)地震1.0g工況下屏蔽廠房頂點(diǎn)最大豎向加速度達(dá)41.13 m/s2(QianAn波)，為安全停堆地震下的3.33倍，超設(shè)計(jì)地震作用下的強(qiáng)烈加速度響應(yīng)對(duì)核島結(jié)構(gòu)構(gòu)件、管道和設(shè)備構(gòu)成極大的安全威脅。可見(jiàn)，AP1000核島結(jié)構(gòu)設(shè)計(jì)應(yīng)該著重分析豎向地震動(dòng)所帶來(lái)的不利影響。

(a)Kern County波 (b)QianAn波 (c)人工波 圖6 屏蔽廠房最大加速度響應(yīng)Fig.6 Maximum acceleration responses of the shield structure

3.2 位移響應(yīng)分析

圖7給出了安全停堆地震和超設(shè)計(jì)地震下屏蔽廠房沿外壁高度分布的位移響應(yīng)。由圖7知，AP1000屏蔽廠房外壁水平位移峰值隨結(jié)構(gòu)高度逐漸增大，無(wú)明顯突變現(xiàn)象，安全停堆地震下最大水平位移為20.30 mm(人工波)，超設(shè)計(jì)地震1.0g工況下最大水平位移為67.14 mm(人工波)，為安全停堆地震下的3.31倍，但由于屏蔽廠房結(jié)構(gòu)整體剛度大，在安全停堆和超設(shè)計(jì)地震下其水平位移峰值較小。

屏蔽廠房豎向位移同樣在圓錐形混凝土壁處發(fā)生突變，Kern County、QianAn和人工波作用下頂點(diǎn)處的豎向位移峰值放大為屏蔽廠房60 m處的倍數(shù)分別為1.47、2.29和1.40。在超設(shè)計(jì)地震1.0g工況下屏蔽廠房的最大豎向位移為11.11 mm(人工波)，雖然在沒(méi)有考慮冷卻水箱儲(chǔ)水的情況下，在三向地震波作用下其豎向位移峰值不大，但考慮到上部冷卻水箱結(jié)構(gòu)的響應(yīng)放大及水箱容積之大，其豎向地震不容忽視。

(a)Kern County波 (b)QianAn波 (c)人工波 圖7 屏蔽廠房最大位移響應(yīng)Fig.7 Maximum displacement responses of the shield structure

3.3 應(yīng)力響應(yīng)分析

表2給出了不同地震動(dòng)峰值工況下AP1000屏蔽廠房結(jié)構(gòu)的最大壓應(yīng)力。從表2中可以看出，在安全停堆地震(0.3g)作用下屏蔽廠房最大壓應(yīng)力為11.8 MPa，在容許應(yīng)力范圍之內(nèi)，超設(shè)計(jì)地震作用下有達(dá)到混凝土應(yīng)力-應(yīng)變關(guān)系的峰值29.6 MPa，此時(shí)結(jié)構(gòu)進(jìn)行非線性狀態(tài)。日本柏崎刈羽核電站在遭遇超設(shè)計(jì)地震作用時(shí)造成了非安全相關(guān)建、構(gòu)筑物的嚴(yán)重破壞，因此，核島結(jié)構(gòu)抗震設(shè)計(jì)中需要考慮超設(shè)計(jì)基準(zhǔn)的保守性。

表2 各工況下屏蔽廠房最大壓應(yīng)力

4 AP1000核島結(jié)構(gòu)隔震控制研究

4.1 隔震設(shè)計(jì)

隔震結(jié)構(gòu)通過(guò)在基礎(chǔ)結(jié)構(gòu)與上部結(jié)構(gòu)間設(shè)置隔震支座，使上部結(jié)構(gòu)與地震動(dòng)水平成分隔離，其中隔震支座由夾層薄鋼板和薄橡膠片相互交錯(cuò)疊置制作而成，橡膠片受壓時(shí)的壓縮變形較大，橡膠片向水平方向膨脹，因此能降低承載能力。針對(duì)AP1000核島結(jié)構(gòu)在超設(shè)計(jì)地震作用下其加速度反應(yīng)強(qiáng)烈，危及到結(jié)構(gòu)構(gòu)件、管道和設(shè)備安全以及屏蔽廠房進(jìn)入非線性問(wèn)題，引入基礎(chǔ)隔震技術(shù)，進(jìn)行結(jié)構(gòu)震動(dòng)控制研究。隔震支座建模時(shí)采用COMBIN14單元模擬豎向性能，水平方向性能采用COMBIN40單元模擬，通過(guò)在結(jié)構(gòu)底部與基礎(chǔ)間加入隔震單元，固結(jié)隔震單元兩端節(jié)點(diǎn)，構(gòu)成隔震體系。

核島結(jié)構(gòu)隔震設(shè)計(jì)采用高阻尼橡膠支座，高阻尼橡膠隔震支座是一種具有較小水平剛度、良好濾波和能量耗散能力的裝置，采用雙線性模型模擬高阻尼隔震支座的滯回曲線(見(jiàn)圖8)。在核島結(jié)構(gòu)底部布置了217個(gè)高阻尼橡膠隔震支座，該支座豎向極限承載力為7 700 kN，等效水平剛度(100%)為1.25 kN/mm，支座高度374 mm，其中橡膠總厚度為204 mm，最大變形為280 mm。采用基礎(chǔ)隔震技術(shù)后結(jié)構(gòu)周期延長(zhǎng)為3.37 s。圖9較明顯地展示了非隔震與隔震核島結(jié)構(gòu)的變形差異，隔震后上部結(jié)構(gòu)在水平方向上做平動(dòng)運(yùn)動(dòng)。

圖8 高阻尼橡膠支座雙線性模型Fig.8 Bilinear hysteretic model of high damping rubber bearings

(a)非隔震 (b)隔震圖9 隔震與非隔震核島結(jié)構(gòu)第一階振型Fig.9 The first mode of isolated and non-isolated structure

4.2 加速度響應(yīng)對(duì)比

以人工波為例，對(duì)比非隔震與隔震核島結(jié)構(gòu)的水平加速度響應(yīng)(見(jiàn)圖10)，隔震后屏蔽廠房沿高度分布的水平加速度變化不明顯，水平加速度大幅度減小，1.0g工況時(shí)頂點(diǎn)水平加速度減幅高達(dá)77.54%，隔震后有效地控制了核島結(jié)構(gòu)的動(dòng)力響應(yīng)，上部結(jié)構(gòu)水平加速度沿高度基本不變，相比于隔震前核島結(jié)構(gòu)由搖擺運(yùn)動(dòng)轉(zhuǎn)變?yōu)槠絼?dòng)。

圖10 人工波作用下屏蔽廠房最大水平加速度Fig.10 Maximum horizontal acceleration of the shield structure under artificial wave

4.3 位移響應(yīng)對(duì)比

圖11所示為非隔震與隔震核島結(jié)構(gòu)在Kern County波作用下核島結(jié)構(gòu)水平方向沿高度分布的最大位移曲線。非隔震核島結(jié)構(gòu)在各工況下呈現(xiàn)明顯的剪切變形，而隔震核島結(jié)構(gòu)位移集中于隔震層，上部結(jié)構(gòu)基本上近似剛體平動(dòng)，結(jié)構(gòu)相對(duì)水平位移很小。在安全停堆地震作用下隔震支座變形在限值280 mm之內(nèi)，當(dāng)超設(shè)計(jì)地震達(dá)0.6g時(shí)支座變形超過(guò)限值，因此在進(jìn)行核島結(jié)構(gòu)隔震設(shè)計(jì)時(shí)，應(yīng)考慮采用大直徑隔震支座及采取相應(yīng)的限位措施，以保證核島結(jié)構(gòu)在超設(shè)計(jì)地震下的性能。

4.4 應(yīng)力響應(yīng)對(duì)比

設(shè)置高阻尼橡膠隔震支座后，人工波作用時(shí)屏蔽廠房混凝土最大壓應(yīng)力在各工況下均減小80%左右(見(jiàn)圖12)，且隔震后屏蔽廠房混凝土最大壓應(yīng)力均在容許應(yīng)力范圍內(nèi)，結(jié)構(gòu)處于彈性狀態(tài)，有效地保護(hù)了核島結(jié)構(gòu)構(gòu)件在超設(shè)計(jì)地震作用下不致破壞。

4.5 基底剪力響應(yīng)對(duì)比

表3所示為Kern County波作用下非隔震與隔震核島結(jié)構(gòu)最大基底剪力對(duì)比。在各工況下，隔震核島結(jié)構(gòu)基底剪力減小率都達(dá)78%以上。

5 結(jié) 論

建立了AP1000核島結(jié)構(gòu)的ANSYS有限元模型，與現(xiàn)有研究成果和Midas Gen計(jì)算結(jié)果對(duì)比，驗(yàn)證了所建立核島結(jié)構(gòu)有限元模型的正確性，基于此模型，分析安全停堆和超設(shè)計(jì)地震作用下核島結(jié)構(gòu)的三維地震響應(yīng)，對(duì)核島結(jié)構(gòu)進(jìn)行隔震控制，研究基礎(chǔ)隔震技術(shù)解決超設(shè)計(jì)地震下響應(yīng)過(guò)大問(wèn)題的有效性，主要結(jié)論如下：

(1)三向地震作用下，屏蔽廠房水平加速度、水平位移峰值沿結(jié)構(gòu)高度均勻變化且最大值發(fā)生在屏蔽廠房頂部，豎向加速度、豎向位移峰值均在冷卻水箱處發(fā)生突變快速增大，最大值同樣發(fā)生在屏蔽廠房頂部，應(yīng)對(duì)超設(shè)計(jì)地震下核島結(jié)構(gòu)的豎向響應(yīng)引起重視。

(2)超設(shè)計(jì)地震作用下，其加速度響應(yīng)為安全停堆地震下的數(shù)倍，頂點(diǎn)最大水平、豎向加速度為安全停堆地震下的3.31倍和3.33倍。

(3)超設(shè)計(jì)地震作用下對(duì)核島結(jié)構(gòu)構(gòu)件、管道和設(shè)備構(gòu)成極大的安全威脅，屏蔽廠房結(jié)構(gòu)應(yīng)力超過(guò)材料容許應(yīng)力，結(jié)構(gòu)進(jìn)入非線性狀態(tài)。

(4) 采用隔震控制后，核島結(jié)構(gòu)加速度、位移、應(yīng)力和基底剪力等得到較好地控制，有效地解決了超設(shè)計(jì)地震下核島結(jié)構(gòu)響應(yīng)強(qiáng)烈的問(wèn)題。此外進(jìn)行核島結(jié)構(gòu)隔震設(shè)計(jì)時(shí)應(yīng)盡可能采用大直徑隔震支座及采取相應(yīng)的限位措施，以保證核島結(jié)構(gòu)在超設(shè)計(jì)地震作用下的性能。

[ 1 ] 潘蓉. 日本柏崎刈羽核電廠新潟地震震害經(jīng)驗(yàn)[J]. 核安全, 2008(1): 20-23.

PAN Rong. The lessons learned from the damage of kashiwazaki-kariwa nuclear power plant caused by chuetu-oki earthquake[J]. Nuclear Safety, 2008(1): 20-23.

[ 2 ] 潘蓉, 李小軍, 胡勐乾. 實(shí)際地震下核電站場(chǎng)地地震反應(yīng)非線性研究[J]. 工業(yè)建筑, 2014, 44(增刊1)： 608-620.

PAN Rong, LI Xiaojun, HU Mengqian. Study on nonlinear seismic response of nuclear power station field under real earthquakes[J]. Industrial Construction, 2014, 44(Sup 1): 608-620.

[ 3 ] 夏祖諷. 核電廠的抗震設(shè)計(jì)輸入及AP1000核島隔震課題簡(jiǎn)介[J]. 中國(guó)工程科學(xué), 2013, 15(4): 52-61.

XIA Zufeng. Seismic input of NPP & topic of seismic-isolated research for AP1000 nuclear island[J]. Engineering Sciences, 2013, 15(4): 52-61.

[ 4 ] 周國(guó)良, 唐暉, 魏超, 等. 新型壓水堆核電廠核島廠房結(jié)構(gòu)地震反應(yīng)分析[C]// 第22屆全國(guó)結(jié)構(gòu)工程學(xué)術(shù)會(huì)議論文集第III冊(cè). 烏魯木齊:新疆大學(xué), 2013:332-337.

[ 5 ] 賀秋梅, 李小軍, 張江偉, 等. 某高溫氣冷堆核電廠結(jié)構(gòu)地震反應(yīng)分析[J]. 震災(zāi)防御技術(shù), 2014, 9(3): 454-461.

HE Qiumei, LI Xiaojun, ZHANG Jiangwei, et al. Seismic response analysis of high temperature gas cooled reactor nuclear power plant[J]. Technology for Earthquake Disaster Prevention, 2014, 9(3): 454-461.

[ 6 ] 謝禮立, 翟長(zhǎng)海. 核電工程應(yīng)用隔震技術(shù)的可行性探討[J]. 地震工程與工程振動(dòng), 2012, 32(1): 1-10.

XIE Lili, ZHAI Changhai. A prospective study on applicability of base isolation in nuclear power plants[J]. Journal of Earthquake Engineering and Engineering Vibration, 2012, 32(1): 1-10.

[ 7 ] 周福霖, 譚平, 劉德穩(wěn). 核電站安全的地震威脅及對(duì)策分析[J]. 中國(guó)工程科學(xué), 2013, 15(4): 36-40.

ZHOU Fulin, TAN Ping, LIU Dewen. The earthquake threat and countermeasures to the safety of nuclear power plants[J]. Engineering Sciences, 2013, 15(4): 36-40.

[ 8 ] 林誠(chéng)格. 非能動(dòng)安全先進(jìn)壓水堆核電技術(shù)[M]. 北京: 原子能出版社, 2010.

[ 9 ] 江見(jiàn)鯨, 陸新征. 混凝土結(jié)構(gòu)有限元分析[M]. 北京: 清華大學(xué)出版社, 2013.

[10] 凌云, 王曉雯, 袁芳. CNP1000安全殼抗震分析[J]. 核技術(shù), 2010, 33(2): 138-142.

LING Yun, WANG Xiaowen, YUAN Fang. CNP1000 containment earthquake resistance analysis[J]. Nuclear Techniques, 2010, 33(2): 138-142.

[11] 段安, 錢稼茹. CNP1000核電廠安全殼模型結(jié)構(gòu)抗震安全分析[J]. 工程力學(xué), 2009, 26(4): 153-157.

DUAN An, QIAN Jiaru. Aseismic safety analysis of a containment vessel model for CNP1000 nuclear power plant[J]. Engineering Mechanics, 2009, 26(4):153-157.

[12] 張有佳, 李小軍, 賀秋梅. 基于分層殼單元的安全殼結(jié)構(gòu)有限元分析[J]. 建筑科學(xué), 2014, 30(5): 35-40.

ZHANG Youjia, LI Xiaojun, HE Qiumei. Finite element analysis of a containment vessel structures based on the multi-layer shell elements[J]. Building Science, 2014, 30(5): 35-40.

[13] 孫鋒, 潘蓉, 王威, 等. 某核電站安全殼隔震動(dòng)力響應(yīng)規(guī)律初探[J]. 工程抗震與加固改造, 2014, 36(3): 27-30.

SUN Feng, PAN Rong, WANG Wei, et al. Simulation analysis of seismic isolation dynamic response lawson concrete containmentin nuclear power plants[J]. Earthquake Resistant Engineering and Retrofitting, 2014, 36(3): 27-30.

[14] 劉立. 某核安全殼振動(dòng)臺(tái)試驗(yàn)數(shù)值模擬研究[D]. 哈爾濱: 哈爾濱工程大學(xué), 2009.

[15] ABRAHAMSON N A. Non-stationary spectral matching[J]. Seismological Research Letters, 1992, 63(1):30-39.

[16] HANCOCK J, WATSON-LAMPREY J, ABRAHAMSON N A, et al. An improved method of matching response spectra of recorded earthquake ground motion using wavelets[J]. Journal of Earthquake Engineering, 2006, 10:67-89.

[17] 王國(guó)波, 王亞西, 于艷麗, 等. 基礎(chǔ)提離對(duì)核電站結(jié)構(gòu)地震響應(yīng)的影響分析[J]. 振動(dòng)與沖擊, 2015, 34(7):228-233.

WANG Guobo, WANG Yaxi,YU Yanli,et al. Analysis on the influence of foundation uplift on structure seismic response[J]. Journal of Vibration and Shock, 2015, 34(7): 228-233.

[18] 周正華, 周雍年, 趙剛. 強(qiáng)震近場(chǎng)加速度峰值比和反應(yīng)譜統(tǒng)計(jì)分析[J]. 地震工程與工程振動(dòng), 2002, 22(3): 15-18.

ZHOU Zhenghua, ZHOU Yongnian, ZHAO Gang. Statistical analysis of peak ratios and response spectra of near-field accelerograms[J]. Earthquake Engineering and Engineering Vibration, 2002, 22(3): 15-18.

A study on three-dimensional seismic vibration response and vibration control of a nuclear island structure under beyond-design basis earthquake

ZHUANG Chuli, ZHANG Yongshan, WANG Dayang

(School of Civil Engineering, Guangzhou University, Guangzhou 510006, China)

A three-dimensional finite element model of AP1000 nuclear island structure (shield building and containment vessel) was built by using ANSYS program. In order to ensure the accuracy of the finite element model, the model was analyzed in two finite element analysis software sets and the results from the two analyses were compared. Moreover, the ANSYS results were compared with the results in literature. The finite element model in ANSYS was used to study the three-dimensional seismic vibration response of nuclear island structure under safe shutdown earthquake and beyond-design basis earthquake. The seismic isolation technology was introduced, and analysis was performed to investigate its effectiveness under beyond-design basis earthquake. The seismic responses of the non-isolated nuclear island structure increased remarkably under the beyond-design basis earthquake compared to the safe shutdown earthquake; the top horizontal and vertical acceleration increased by 3.31 and 3.33 times, respectively, and the seismic stresses in the shield building structure exceeded the allowable stresses of materials. Compared to the non-isolated nuclear island structure, the base seismic isolation system had largely reduced horizontal accelerations, superstructure displacements, seismic stresses and base shear force under the beyond-design basis earthquake. With the base isolation, the rate of reduction in the base shear force was more than 78%. Therefore the base seismic isolation system can effectively ensure the safety of a nuclear island structure under the beyond-design basis earthquake.

nuclear island structure; three-dimensional seismic vibration; beyond-design basis earthquake; base seismic isolation

國(guó)家自然科學(xué)基金項(xiàng)目(51378135；51408140)；教育部博士點(diǎn)基金項(xiàng)目(20134410120003)；廣州市科技計(jì)劃項(xiàng)目(201510010291)

2016-04-05 修改稿收到日期： 2016-06-30

莊初立 男，博士生, 1990年生

張永山 男，博士，教授，1964年生

TU352.1

A

10.13465/j.cnki.jvs.2017.16.037