考慮地基不均勻性核電廠房結(jié)構(gòu)擬三維抗震安全分析

尹訓(xùn)強，羅 勇，王桂萱

(大連大學(xué) 土木工程技術(shù)研究與開發(fā)中心， 遼寧 大連 116622)

由于我國是多地震國家，地震活動頻繁、強度大、范圍廣，其中地震活動區(qū)占國土面積達60%以上，因此在對核電站設(shè)計時，需要考慮的首要問題是提高其抗震能力。在考慮局部場地特征及地震作用下，為了保證核電廠房結(jié)構(gòu)的及核電站設(shè)備完整性和運行的可行性，是核電抗震安全的重中之重。然而，在目前的核電站設(shè)計中考慮到設(shè)計的成本和核電站廠房結(jié)構(gòu)的安全性，往往都將核電站的設(shè)計方案和配套設(shè)施都進行統(tǒng)一化處理，即在核電廠房結(jié)構(gòu)設(shè)計過程中均采用標(biāo)準(zhǔn)化設(shè)計。而我國國土面積覆蓋面廣，各個地區(qū)地質(zhì)條件及復(fù)雜程度參差不齊，使得這種標(biāo)準(zhǔn)化設(shè)計面臨諸多問題。因此，出于對核電站廠房結(jié)構(gòu)的安全性考慮，在對核電站廠址的地基條件作判定時，對地基進行適應(yīng)性分析是必不可少的重要環(huán)節(jié)。而由于非均質(zhì)地基問題具有更好的代表性，進而在對核電廠房結(jié)構(gòu)進行地震響應(yīng)分析時，考慮地基不均勻性顯得非常有必要。

目前雖然對于土體非線性結(jié)構(gòu)最常用的方法是進行二維分析，但這種方法往往會低估結(jié)構(gòu)的抗震性能，另一方面，二維分析不能考慮土體結(jié)構(gòu)的平面外效應(yīng)。但是，對于真實現(xiàn)象來說，三維分析被認為是最為重視的。然而，創(chuàng)建復(fù)雜地基和動態(tài)的三維分析模型需要相當(dāng)長的時間且不夠經(jīng)濟[1-2]。為了解決上述問題，在對非均質(zhì)地基條件下的核電廠進行抗震安全分析時，對非均質(zhì)地基地條件下的核電廠房結(jié)構(gòu)進行擬三維分析及研究非均質(zhì)地基地條件下的核電廠房結(jié)構(gòu)擬三維地震響應(yīng)影響規(guī)律具有十分重要的意義。

因此，本文將以某建立在非均質(zhì)地基條件上內(nèi)陸核電廠CPR1000為研究對象，研究非均質(zhì)地基條件下核電廠房結(jié)構(gòu)擬三維震響應(yīng)影響規(guī)律，并與真三維模型的計算結(jié)果進行了對比分析，進一步對所建立擬三維計算模型的有效性及可靠性進行考察驗證。

1 擬三維動力計算基本理論

1.1 土-結(jié)構(gòu)相互作用運動方程

土-結(jié)構(gòu)相互作用運動方程：

(1)

式中：M為結(jié)構(gòu)的質(zhì)量矩陣；C為阻尼矩陣；K為剛度矩陣；u為位移向量；P為外力向量。

1.2 考慮平面外效應(yīng)的求解方程

考慮平面外的能量傳遞的擬三維求解運動方程：

(2)

式中：模型節(jié)點相對于剛性基礎(chǔ)的位移由u表示；切片平面應(yīng)變質(zhì)量矩陣由M表示，m是m與剛性基礎(chǔ)加速度a方向有關(guān)的向量；其中需要說明的是在由復(fù)模量構(gòu)成的剛度矩陣K中需要考慮材料的阻尼比，切片的平面厚度為單位厚度。

V是由切片平面內(nèi)兩側(cè)面的能量傳遞邊界產(chǎn)生的力，其表達式表示為：

(3)

F作用在切片的端部力，而不是作用在自由場的垂直面，可以表示為：

F=Guf

(4)

式中：G為簡單剛度矩陣，是自由場中與頻率無關(guān)的復(fù)模量構(gòu)成的。

T為能量傳遞有關(guān)的力，其表達式表示為：

T=(R+L)(u-uf)

(5)

其中R和L是由Lysmer等[3-4]引進的與頻率有關(guān)的邊界剛度矩陣。

2 非均質(zhì)地基條件下核電廠房擬三維地震分析模型

2.1 擬三維簡化模型的建立

為了便于理解，將提供簡單的示例進行說明。如圖1所示，在具有三維傾斜斜面的情況下，設(shè)定分析范圍當(dāng)如圖2所示存在三維傾斜斜面時，設(shè)定分析范圍，假定如圖3所示那樣分別獨立的二維橫截面。這里，將各區(qū)域稱為區(qū)域1、區(qū)域2、區(qū)域3。有效寬度取到各個區(qū)域的代表面的距離的一半?？紤]到本文研究對象核反應(yīng)堆上部結(jié)構(gòu)的復(fù)雜性，此方法主要針對下部結(jié)構(gòu)地基的簡化。

圖1 三維傾斜斜面

圖2 擬三維解析模型圖

圖3擬三維分析中二維平面示意圖

2.2 土-結(jié)構(gòu)相互作用計算模型

目前，在核電廠房結(jié)構(gòu)動力響應(yīng)分析及抗震設(shè)計中，考慮土-結(jié)構(gòu)動力相互作用的影響是必不可少的一個環(huán)節(jié)。雖然目前無限地基動力相互作用計算模型種類較多，但在眾多土-結(jié)構(gòu)動力相互作用模型中應(yīng)用比較成熟的主要有四種：人工邊界法[4-6]、邊界元法[7]、比例邊界有限元法[8]和阻尼抽取法[9-10]。相較于其他子結(jié)構(gòu)模型，人工邊界法具有其自身獨特的優(yōu)勢：一方面，在人工邊界法中可以截取的地基有限區(qū)域邊界處設(shè)置不同形式的能量邊界條件，進而能夠有效準(zhǔn)確地?zé)o限地基遠場的波動效應(yīng)；另一方面，在地基處理中可方便地考慮土體非線性和土-結(jié)構(gòu)相互作用的影響。其中黏彈性邊界的主要思想是在人工邊界節(jié)點處設(shè)置一系列的彈簧-阻尼器單元(如圖4所示)，基于此黏彈性邊界單元不僅能有效地模擬無限地基輻射阻尼效應(yīng)，同時還能有效地模擬人工邊界外半無限空間介質(zhì)的彈性恢復(fù)性能，此外還具有良好的頻率穩(wěn)定性和較高的精度以及較好的魯棒性[11]。因此黏彈性邊界在相關(guān)的科研和工程中得到更多應(yīng)用。

圖4黏彈性外邊界數(shù)值模型示意圖

3 擬三維黏彈性人工邊界的二次開發(fā)

3.1 地震波的輸入方法

在對結(jié)構(gòu)-地基動力相互作用問題進行分析時，為了在黏彈性人工邊界處實現(xiàn)地震波的輸入，借助于劉晶波[12]介紹的一種地震波的直接輸入方法，即在邊界節(jié)點處采用力學(xué)中脫離體的概念，其主要思想是將輸入地震動轉(zhuǎn)化為作用于人工邊界節(jié)點上的等效荷載的方法來實現(xiàn)地震波動輸入。人工邊界節(jié)點處的等效荷載可表達為：

(6)

3.2 擬三維分析法中用戶自定義單元開發(fā)流程及嵌入

在ANSYS軟件中的二次開發(fā)工具主要有4個：APDL、UPFs、UIDL、TclTk，其中UPFs是ANSYS二次開發(fā)的核心工具，UPFs是基于FORTRAN程序的二次開發(fā)，其主要功能是實現(xiàn)對ANSYS功能的擴充，即對ANSYS單元庫的擴充。同時ANSYS軟件也為用戶提供了諸多子程序，其中，UserElem.f子程序不僅可以幫助用戶傳遞自定義單元所需的的必要數(shù)據(jù)，同時還可以在動力計算過程中自動更新標(biāo)準(zhǔn)ANSYS中的數(shù)據(jù)庫和相關(guān)文件。需要說明的是UserElem.f子程序是基于ANSYS單元水平上的接口。基于UserElem.f子程序，用戶可以幾乎創(chuàng)建各種單元類型而無需直接訪問ANSYS數(shù)據(jù)庫和文件。因此,利用UPFs以及UserElem的特性,并結(jié)合二維及三維黏彈性人工邊界有關(guān)理論,用戶自定義擬三維黏彈性人工邊界的開發(fā)具體流程詳見參考文獻[13]。

為了實現(xiàn)將用戶自定義擬三維黏彈性人工邊界單元嵌入到ANSYS中方便于工程應(yīng)用，需要將接口子程UserElem與標(biāo)準(zhǔn)程序進行連接，其具體連接過程詳見參考文獻[13]，進而可以實現(xiàn)在ANSYS中對用戶自定義單元的嵌入。再根據(jù)2.1節(jié)所介紹的方法完成對下部結(jié)構(gòu)地基的簡化，簡化后的下部結(jié)構(gòu)成為二維片狀結(jié)構(gòu)，但上部結(jié)構(gòu)仍為三維結(jié)構(gòu)，由于用戶自定義單元無法對下部二維片狀結(jié)構(gòu)進行網(wǎng)格剖分，故采用MESH200單元對簡化后的下部二維平面進行有限元網(wǎng)格剖分，再利用ANSYS內(nèi)置命令將MESH200修改替換為用戶自定義的擬三維黏彈性邊界單元即可。

4 算例驗證

4.1 計算模型和地基計算參數(shù)

核電廠房上部結(jié)構(gòu)采用集中質(zhì)量模型，模型相關(guān)參數(shù)詳見文獻[14]，地基動態(tài)參數(shù)如表1所示。CPR1000 反應(yīng)堆核島廠房(RX)主要由三大部分組成，這三大部分主要包括：筏板基礎(chǔ)、安全殼結(jié)構(gòu)、內(nèi)部結(jié)構(gòu)。反應(yīng)堆核島廠房結(jié)構(gòu)整體采用集中質(zhì)量-梁單元模擬，結(jié)構(gòu)的質(zhì)量和轉(zhuǎn)動慣量均集中在各節(jié)點上，而兩相鄰節(jié)點間的幾何慣性矩和剪切面積用連接節(jié)點的梁單元來進行模擬，計算模型見圖5，圖中節(jié)點編號為1～12，梁單元編碼為①～，結(jié)構(gòu)在平面上即X和Y向是對稱的。核電廠房結(jié)構(gòu)核反應(yīng)堆模型如圖6所示。

圖5 ANSYS三維有限元分析模型

圖6核島廠房結(jié)構(gòu)核反應(yīng)堆模型

選取計算模型如圖7、圖8所示(不同地層剖面方向)模型地基長為120 m，寬為120 m，高為40 m。為了考慮地震動力響應(yīng)中的地基-結(jié)構(gòu)相互作用的影響，分析模型同時包括地基與核反應(yīng)堆廠房結(jié)構(gòu)模型，利用大型通用有限元軟件ANSYS建立了X-Y平面和Y-Z平面模型，平面單元采用Mesh200單元來進行離散，進而實現(xiàn)對真實三維ANSYS有限元模型(見圖5)的模擬，計算中選擇X-Z平面和Y-Z平面進行動力響應(yīng)分析，為了考慮地震動力響應(yīng)中地震方向隨機性的影響，在計算中將地震波分成兩次來進行輸入，兩次地震動的輸入方向均為水平方向和豎直向。地基計算參數(shù)見表1。

圖7 擬三維X-Z平面

圖8 擬三維Y-Z平面

4.2 地震動的輸入

地震動加速度時程曲線，如圖9所示，其中，水平向X、Y向地面運動加速度峰值均為0.932 m/s2，豎直向Z向地面運動加速度峰值為0.540 m/s2，總持時為25 s，時間步長為0.01 s。三維有限元模型及擬三維有限元模型均在地基部分邊界處施加有黏彈性人工邊界，地震動輸入方式按照3.1節(jié)方法輸入，從基礎(chǔ)底部垂直輸入。

圖9輸入地震動時程曲線

4.3 計算結(jié)果對比分析

考慮到集中質(zhì)量點1、8、12位置的特殊性，即在整個計算模型的頂部，相對于下部結(jié)構(gòu)節(jié)點的地震響應(yīng)而言，這三個觀測點的地震響應(yīng)較大，因此，提取節(jié)點1、8、12三個觀測點三個方向阻尼比為5%的加速度反應(yīng)譜，如圖10、圖11所示，利用ANSYS軟件對X-Y平面及Y-Z平面采用擬三維動力分析法計算加速度反應(yīng)譜與真三維計算結(jié)果進行對比分析。表2為反應(yīng)堆廠房結(jié)構(gòu)典型位置位移峰值。

從圖10、圖11可以看出，隨著各觀測點高程的增加，節(jié)點1、12、8的加速度及位移響應(yīng)均呈現(xiàn)出逐漸增大的趨勢，同時各節(jié)點加速度反應(yīng)譜曲線及位移變化趨勢也基本一致。將真三維計算結(jié)果與擬三維計算結(jié)果對比來看，加速度反應(yīng)譜曲線變化趨勢基本相同，雖然后者的計算結(jié)果整體量值有所減小，減小幅度均在15%～20%范圍內(nèi)，最大減小幅度約為15%，從表2中也可以看出，在擬三維下計算的典型位置的位移峰值相對于真三維計算結(jié)果而言，整體量值有所減小，減小幅度均在10%范圍內(nèi)，基本能滿足工程實際運用要求，綜上所述，可以驗證在ANSYS中建立的不均勻性巖性地基條件下的擬三維動力計算模型具有一定的可靠性與準(zhǔn)確性。

圖10 節(jié)點1加速度反應(yīng)譜對比

5 結(jié) 論

(1) 通過有效寬度選出典型二維獨立截面，可以實現(xiàn)對真實三維模型的簡化，進而實現(xiàn)對非均質(zhì)巖性地基條件下核電廠房結(jié)構(gòu)擬三維的地震響應(yīng)分析。

圖11核島上部結(jié)構(gòu)樓層加速度反應(yīng)譜

(2) 通過三個觀測點的加速度反應(yīng)譜，可以看出，擬三維計算結(jié)果與真三維計算結(jié)果變化趨勢基本一致，僅幅值有一定差距，但差距并不大，驗證了本文所建立的擬三維模型的準(zhǔn)確性與可靠性。

(3) 在對大型復(fù)雜實際工程作地震響應(yīng)分析時，通過擬三維法對核電廠房結(jié)構(gòu)作地震響應(yīng)分析，可以獲得更加合理的經(jīng)濟效益，同時相對于真三維的建模時間及計算效率而言，擬三維動力分析法都有著自身明顯的優(yōu)勢。