韓 治,唐 暉,張春明,文 靜

(環(huán)境保護(hù)部核與輻射安全中心,北京 100082)

核電廠中存在大量的核安全級(jí)薄壁儲(chǔ)液容器,如輔助給水系統(tǒng)ASG水箱和換料系統(tǒng)PTR水箱等,其需要在地震情況下保持結(jié)構(gòu)完整性及其功能。GB50267-97《核電廠抗震設(shè)計(jì)規(guī)范》要求 “對(duì)于薄壁儲(chǔ)液容器計(jì)算應(yīng)計(jì)入器壁柔度的影響,并對(duì)壓應(yīng)力進(jìn)行儲(chǔ)液容器壁的失穩(wěn)校核”[1]。

大型儲(chǔ)液容器在商業(yè)、軍用、民用、核電廠,尤其是在石油化工業(yè)中十分常見,多年來(lái)各領(lǐng)域的科研人員針對(duì)儲(chǔ)液容器的地震響應(yīng)問(wèn)題建立了行之有效的分析方法,詳細(xì)發(fā)展歷程請(qǐng)參考綜述性文獻(xiàn)[2]。1957年,Housner將剛性儲(chǔ)液容器化為質(zhì)量彈簧系統(tǒng)的簡(jiǎn)化模型[3],模型簡(jiǎn)單實(shí)用,為六、七十年代儲(chǔ)液容器抗震計(jì)算所廣泛采用。然而,容器壁的剛性假設(shè),導(dǎo)致求出的脈動(dòng)液壓過(guò)小,按Housner理論設(shè)計(jì)的儲(chǔ)液容器在地震中仍大量損壞,從而使科研人員進(jìn)一步開展柔性壁錨固式儲(chǔ)液容器的抗震研究。ASCE 4-98適用于核設(shè)施所有核安全相關(guān)結(jié)構(gòu),包括且不限于地面及地下結(jié)構(gòu)、埋置管道、地面立式容器以及帶隔震系統(tǒng)的結(jié)構(gòu)[4],核電廠大型儲(chǔ)液容器屬于該適用范圍內(nèi)。ASCE 4-98標(biāo)準(zhǔn)的3.5節(jié)中給出了立式儲(chǔ)液容器的集中質(zhì)量彈簧模型與抗震計(jì)算方法,考慮了柔性壁影響,以樓層反應(yīng)譜為依據(jù)確定地震動(dòng)輸入?；诩匈|(zhì)量-彈簧模型,多位研究人員在核電廠儲(chǔ)液容器的地震響應(yīng)分析方面取得進(jìn)展[5-8]。

近年來(lái),隨著計(jì)算機(jī)硬件性能的提升和流固耦合計(jì)算方法的發(fā)展,利用流固耦合數(shù)值模擬方法,如任意拉格朗日歐拉方法 (ALE)、位移流體單元法、勢(shì)流體單元法、SPH粒子法和移動(dòng)粒子法 (MPS)等,進(jìn)行儲(chǔ)液容器地震響應(yīng)分析研究工作越來(lái)越多[9-14]。其中,陸道綱等[15]采用ALE方法、衛(wèi)媛媛等[12,16]采用勢(shì)流體單元法[12]、MPS方法[16]研究了快中子反應(yīng)堆主容器的自由液面晃動(dòng)問(wèn)題。

本文利用ASCE4-98建議的集中質(zhì)量-彈簧模型方法 (以下簡(jiǎn)稱MASS-SPRING方法)和ANSYS流固耦合方法 (以下簡(jiǎn)稱FSI方法)進(jìn)行核電廠儲(chǔ)液容器地震響應(yīng)分析,全面比較了兩者之間的差異,并分析了出現(xiàn)差異的原因。

1 儲(chǔ)液容器地震響應(yīng)分析方法

1.1 集中質(zhì)量-彈簧模型方法

ASCE4-98基于Housner理論給出了立式儲(chǔ)液容器的抗震計(jì)算模型,即集中質(zhì)量-彈簧模型,見圖1。本文不對(duì)Housner計(jì)算模型做詳細(xì)介紹,僅列出與本文有關(guān)的一些假設(shè)[4]:流體無(wú)粘、無(wú)旋、不可壓;容器具有剛性的豎直側(cè)壁和水平底板 (由剛性容器導(dǎo)出的流體模型連接到彈性容器中,才能實(shí)現(xiàn)流體與彈性容器的耦合振動(dòng))。本文涉及的質(zhì)量-彈簧模型及相關(guān)參數(shù)的計(jì)算公式見下文,主要取自文獻(xiàn)[17]。

等效脈沖質(zhì)量公式如下

其中,MF是液體總質(zhì)量,R是容器半徑,h是液位高度。

圖1 柔性壁儲(chǔ)液容器水平地震響應(yīng)分析的流體質(zhì)量分布Fig.1 Distribution of fluid mass for horizontal seismic response analysis of storage tank with flexible walls

圖1中,MS表示儲(chǔ)液容器質(zhì)量;MI表示脈沖質(zhì)量;MC表示對(duì)流質(zhì)量;lI表示脈沖質(zhì)量施加的區(qū)域;lC表示對(duì)流質(zhì)量施加的區(qū)域。

脈沖質(zhì)量的質(zhì)心高度為

對(duì)流質(zhì)量公式為

對(duì)流質(zhì)量的質(zhì)心高度HC為

流體晃動(dòng)對(duì)應(yīng)的角頻率滿足:

1.2 ANSYS流固耦合有限元方法

對(duì)于地震中的儲(chǔ)液容器,工程上主要關(guān)注儲(chǔ)液容器的形變和應(yīng)力分布,通?？梢约僭O(shè)流體為無(wú)粘、理想可壓縮和小擾動(dòng)的,固體為線彈性的。采用流固耦合方法分析動(dòng)力學(xué)問(wèn)題時(shí),對(duì)固體結(jié)構(gòu)采用位移為基本變量,依據(jù)定義的流體和固體之間耦合作用傳遞方式的不同,流體可以采用位移、壓力以及各種勢(shì)函數(shù)作為基本變量,相應(yīng)的流固耦合有限元格式可以是位移-壓力格式、位移-位移格式和位移-勢(shì)函數(shù)格式。ANSYS的流體單元FLUID80對(duì)應(yīng)的是位移-位移格式,相應(yīng)的有限元方程為

其中,Ms是固體的質(zhì)量矩陣;Mf是流體的質(zhì)量矩陣;Ks是固體的剛度矩陣;Kf是流體的剛度矩陣;Us是固體節(jié)點(diǎn)位移向量;Uf是流體節(jié)點(diǎn)位移向量。

2 算例:核電廠ASG水箱地震響應(yīng)分析

2.1 ASG水箱結(jié)構(gòu)

應(yīng)急補(bǔ)水箱 (ASG)是核電廠余熱排出系統(tǒng)中的重要設(shè)備,在事故情況下確保向蒸汽發(fā)生器供應(yīng)足夠的水,以排出堆芯的衰變熱,正確地評(píng)價(jià)ASG水箱在地震中的響應(yīng)對(duì)核安全至關(guān)重要。

ASG水箱為立式儲(chǔ)液容器[5],由封頭、筒體、封底、筒體加強(qiáng)環(huán)等組成,結(jié)構(gòu)如圖2所示。筒體由五層不同厚度的筒體段焊接而成,厚度分別為16mm、12mm、10mm、8mm、8mm,高度分別為3000mm、3000mm、3000mm、3000mm和2600mm,筒體外直徑為9700mm,加強(qiáng)環(huán)外半徑為5350mm,內(nèi)半徑4850mm,厚度為10mm,液位標(biāo)高為14520mm。容器主要材料為20HR-B,密度7850kg·m-3,泊松比0.3,材料力學(xué)性能見表1。

圖2 ASG水箱結(jié)構(gòu)示意圖Fig.2 Schematic diagram of ASG tank

表1 材料的力學(xué)性能參數(shù)Table 1 Mechanical properties of the material

2.2 地震動(dòng)輸入條件

本文的目的是研究 (滿水狀態(tài))儲(chǔ)液容器地震響應(yīng)分析方法,SSE地震的結(jié)構(gòu)阻尼比取4%,地震作用下流體阻尼比取0.5%,為便于計(jì)算結(jié)果之間的比較,將地震動(dòng)輸入條件簡(jiǎn)化為如表2所示。

表2 SSE地震水平方向反應(yīng)譜Table 2 seismic response spectrum of horizontal direction

2.3 模態(tài)分析

2.3.1 水箱模態(tài)分析 (MASS-SPRING方法)基于MASS-SPRING方法建立水箱有限元模型,如圖3所示,水箱采用殼單元SHELL281,質(zhì)量點(diǎn)采用質(zhì)點(diǎn)單元MASS21,彈簧采用彈簧單元COMBIN14。對(duì)水箱進(jìn)行模態(tài)分析,前兩階為對(duì)流質(zhì)量晃動(dòng)頻率,均為0.3048 Hz;水箱整體一階梁式振型的頻率為6.865 Hz。

圖3 集中質(zhì)量-彈簧有限元模型Fig.3 Finite element model based on Mass-Spring Method

圖4 振型模態(tài)參與系數(shù)比分布 (MASS-SPRING方法)Fig.4 Modal participation coefficient ratio based on MASS-SPRING Method

模態(tài)分析得到的振型中就只有前兩階為代表流體晃動(dòng)的振型,其他的都是水箱的振型。圖4給出了MASS-SPRING方法得到的振型模態(tài)參與系數(shù)比分布情況,模態(tài)參與系數(shù)比(M.C.RATIO)高于0.00001的振型有50個(gè),其中代表液體晃動(dòng)的振型2個(gè),整體殼式振型20個(gè),整體梁式振型28個(gè)。

2.3.2 水箱模態(tài)分析 (FSI方法)

基于ANSYS FSI方法建立水箱有限元模型,如圖5所示,水箱采用殼單元SHELL181,流體采用流體單元FLUID80。流體晃動(dòng)一階振型2個(gè),頻率均為0.3103 Hz,水箱整體一階梁式振型2個(gè),頻率均為9.195 Hz。

模態(tài)參與系數(shù)比高于0.00001的流體振型413個(gè),梁式振型7個(gè),殼式振型5個(gè),從圖6所示的振型模態(tài)參與系數(shù)比 (M.C.RATIO)分布圖,流體晃動(dòng)振型和容器整體梁式振型占據(jù)主導(dǎo)地位。

圖5 ANSYSFSI有限元模型Fig.5 Finite element model based on FSI Method

圖6 FSI方法得到的水箱振型模態(tài)系數(shù)比分布Fig.6 Modal participation coefficient ratio based on FSI Method

2.3.3 模態(tài)分析結(jié)果的對(duì)比

MASS-SPRING方法得到的流體一階振型的頻率為0.3048Hz,ANSYS FSI方法得到的流體一階振型的頻率為0.3103Hz,由公式(5)得到的流體一階振型的頻率為0.3068Hz,基本是一致的。

但是,比較圖4和圖6可以發(fā)現(xiàn),MASSSPRING方法得到的振型中代表流體晃動(dòng)的振型只有兩階,殼式振型占據(jù)了重要的地位;ANSYS FSI方法模態(tài)分析得到的振型中,流體振型和容器梁式振型占據(jù)主導(dǎo)地位。MASSSPRING方法得到的一階梁式振型的頻率為6.865Hz,ANSYS FSI方法得到的一階梁式振型的頻率為9.195Hz。兩種方法得到的一階梁式振型區(qū)別很大,這主要是因?yàn)閮煞N方法對(duì)于地震中流體固體之間的相互作用采用的建模方法差異很大,研究發(fā)現(xiàn)MASS-SPRING方法將流體對(duì)固體的脈沖作用簡(jiǎn)化為脈沖質(zhì)量,施加在水箱的下半部,脈沖質(zhì)量和水箱在地震中完全同步,也就是等價(jià)于水箱下半部的水體在地震中和水箱完全同步,然而真實(shí)地震中水體達(dá)不到和水箱完全的同步運(yùn)動(dòng),因此MASSSPRING方法的脈沖質(zhì)量的假設(shè)方法,引入了明顯的保守性。

2.4 地震響應(yīng)分析

2.4.1 水箱地震響應(yīng)分析 (MASS-SPRING方法)

由于儲(chǔ)液容器在地震中主要的破壞形式是容器側(cè)壁底部區(qū)域附近的屈曲變形,所以實(shí)際工程問(wèn)題關(guān)注水箱底部區(qū)域側(cè)壁的主應(yīng)力S1和垂直軸向應(yīng)力Sy。

圖7 水箱側(cè)壁面第一主應(yīng)力S1Fig.7 First principal stress of side wall of tank

考慮靜水壓力、自重和水平X方向地震載荷,質(zhì)量-彈簧模型方法得到如圖7、圖8所示應(yīng)力結(jié)果,最大主應(yīng)力S1出現(xiàn)在底部,為56.8MPa,最大的垂直軸向應(yīng)力Sy同樣出現(xiàn)在底部,為56.8 MPa。比較主應(yīng)力S1和軸向應(yīng)力Sy的分布,可以看出在水平X方向地震載荷激勵(lì)下,水箱側(cè)壁面上Sy占據(jù)主導(dǎo),符合地震中屈曲失穩(wěn)破壞的容器的應(yīng)力分布特點(diǎn)。

圖8 水箱側(cè)壁面的垂直軸向應(yīng)力SyFig.8 Vertical axial stress of side wall of tank

2.4.2 水箱地震響應(yīng)分析 (FSI方法)

基于ANSYS FSI方法的水箱地震響應(yīng)分析,首先進(jìn)行水箱靜力學(xué)分析,再進(jìn)行預(yù)應(yīng)力模態(tài)分析和譜分析。從圖9、圖10應(yīng)力分布圖可以看出,最大主應(yīng)力S1出現(xiàn)在水箱底部,達(dá)到22.3 MPa,最大軸向應(yīng)力Sy出現(xiàn)在水箱底部,達(dá)到22.1 MPa,對(duì)比S1和Sy可以看出,水箱底部的應(yīng)力Sy占據(jù)主導(dǎo)地位。圖11、圖12給出重力作用下水箱的應(yīng)力結(jié)果,最大周向應(yīng)力為26.2MPa,最大軸向應(yīng)力為7.88MPa,靜力學(xué)分析得到的應(yīng)力分布是預(yù)應(yīng)力模態(tài)分析的初始輸入數(shù)據(jù)。水箱在地震激勵(lì)下真實(shí)的應(yīng)力分布需要由譜分析應(yīng)力結(jié)果迭加靜力學(xué)分析應(yīng)力結(jié)果來(lái)得到。

圖9 水箱側(cè)壁面的第一主應(yīng)力S1Fig.9 First principal stress of side wall of tank

圖10 水箱側(cè)壁面的軸向應(yīng)力SyFig.10 Vertical axial stress of side wall of tank

圖11 重力載荷作用下水箱的周向應(yīng)力Fig.11 Circumferential stress of tank under gravity loading

圖12 重力載荷作用下水箱的軸向應(yīng)力Fig.12 Axial Stress of Water Tank under Gravity Load

2.4.3 地震響應(yīng)分析結(jié)果的對(duì)比

大型薄壁儲(chǔ)液容器在地震中的主要破壞形式是屈曲失穩(wěn),軸向應(yīng)力Sy是屈曲與否的關(guān)鍵。從圖8、圖10、圖12中選取底部附近水箱側(cè)壁面的垂直軸向應(yīng)力Sy,形成圖13所示的軸向應(yīng)力Sy對(duì)比結(jié)果?？梢钥闯鰞煞N方法得到的軸向應(yīng)力Sy趨勢(shì)相同,底部軸向應(yīng)力Sy最大,隨垂直坐標(biāo)Y的升高,水箱側(cè)壁面軸向應(yīng)力Sy不斷下降,在水箱壁厚變化的結(jié)構(gòu)不連續(xù)處Sy出現(xiàn)波動(dòng)現(xiàn)象。

圖13 水箱下部壁面垂直方向應(yīng)力Sy計(jì)算結(jié)果對(duì)比圖Fig.13 Comparison of the vertical stress in bottom of tank

但是,兩種方法得到的軸向應(yīng)力存在明顯的差異,MASS-SPRING方法得到的應(yīng)力Sy比ANSYS FSI方法得到的應(yīng)力Sy高80%左右。

究其原因,由表3所示的振型有效質(zhì)量、振型有效質(zhì)量與總有效質(zhì)量之比這兩個(gè)參數(shù)可以計(jì)算出,MASS-SPRING方法建立的有限元模型的總有效質(zhì)量為1319937kg,而ANSYS FSI方法建立的有限元模型的總有效質(zhì)量為1172651kg(約等于水箱的真實(shí)總質(zhì)量)。這說(shuō)明,MASS-SPRING方法在建模過(guò)程中確定的對(duì)流質(zhì)量和脈沖質(zhì)量偏保守,計(jì)算中發(fā)現(xiàn),脈沖質(zhì)量對(duì)模態(tài)分析得到的梁式振型的頻率和地震響應(yīng)分析得到的應(yīng)力結(jié)果起著決定性的影響,計(jì)算結(jié)果對(duì)對(duì)流質(zhì)量不敏感。MASS-SPRING方法將脈沖質(zhì)量施加在水箱的下半部,脈沖質(zhì)量和水箱在地震中完全同步,使得地震響應(yīng)分析得到的軸向應(yīng)力Sy結(jié)果偏保守。

表3 重要模態(tài)有效質(zhì)量和模態(tài)參與系數(shù)比Table 3 Effective mass and modal participation coefficient ratio of significant mode

3 結(jié) 論

本文分別利用MASS-SPRING方法和ANSYS FSI方法對(duì)核電廠ASG水箱進(jìn)行地震響應(yīng)分析,兩種方法得到的地震響應(yīng)分析結(jié)果差異明顯。

兩種方法得到的模態(tài)分析和地震響應(yīng)分析結(jié)果出現(xiàn)差異的原因如下:

(1)MASS-SPRING方法把流體對(duì)水箱的作用等效為集中質(zhì)量-彈簧系統(tǒng),施加在水箱上部,脈沖質(zhì)量均布在水箱下半部,從而將流固耦合問(wèn)題轉(zhuǎn)化為傳統(tǒng)的結(jié)構(gòu)動(dòng)力學(xué)問(wèn)題。模態(tài)分析結(jié)果表明,MASS-SPRING方法建立的有限元模型的總有效質(zhì)量為1319937 kg,而ANSYS FSI方法建立的有限元模型的總有效質(zhì)量為1172651 kg(約等于水箱的真實(shí)總質(zhì)量),MASS-SPRING方法建模過(guò)程中的保守假設(shè),使得有限元模型的有效質(zhì)量明顯超出真實(shí)總質(zhì)量。

(2)MASS-SPRING方法將流體對(duì)固體的脈沖作用簡(jiǎn)化為脈沖質(zhì)量,施加在水箱的下半部,脈沖質(zhì)量和水箱在地震中完全同步,也就是等價(jià)于水箱下半部的水體在地震中和水箱完全同步,然而真實(shí)地震中水體達(dá)不到和水箱完全的同步運(yùn)動(dòng),因此MASSSPRING方法的脈沖質(zhì)量的假設(shè)方法,引入了明顯的保守性,相應(yīng)的,地震響應(yīng)分析得到的應(yīng)力結(jié)果偏保守。

綜上所述,基于兩種方法進(jìn)行核電廠儲(chǔ)液容器的地震響應(yīng)分析,ANSYS FSI方法得到的應(yīng)力結(jié)果明顯小于MASS-SPRING方法得到的應(yīng)力結(jié)果,MASS-SPRING方法得到的應(yīng)力結(jié)果是偏保守的。對(duì)于ANSYS FSI方法得到的應(yīng)力結(jié)果將進(jìn)一步進(jìn)行試驗(yàn)驗(yàn)證。

[1]國(guó)家地震局.GB50267-97核電廠抗震設(shè)計(jì)規(guī)范[S].北京:中華人民共和國(guó)建設(shè)部,1997.

[2]溫德超,鄭兆昌,孫煥純.儲(chǔ)液罐抗震研究的進(jìn)展[J].力學(xué)進(jìn)展,1995,25(1):60-76.

[3]HOUSNER G W.Dynamic pressures on accelerated fluid containers[J].Bulletin of the Seismological Society of America.1957,47(1):15-35.

[4]ASCE.ASCE 4-98 Seismic Analysis of Safe-related Nuclear Structures and Commentary[S].Reston:ASCE,1999.

[5]文靜,路燕,徐宇,等.核電廠輔助給水系統(tǒng)水箱抗震分析研究[J].原子能科學(xué)技術(shù),2015,49(5):897-902.

[6]李楠,韓治.采用組合質(zhì)量-彈簧模型的快堆臥式貯鈉罐晃動(dòng)分析[J].原子能科學(xué)技術(shù),2015,49(9):1642-1647.

[7]WEN Jing,LU Daogang,GAO Xiaoan,et al.Mass-spring model used to simulate the sloshing of fluid in the container under the earthquake[C].18th International Conference on Structural Mechanics in Reactor Technology(SMiRT-18).Beijing:Atomic Energy Press,2005:3636-3644.

[8]劉嘉一,白欣然,李曉軒.質(zhì)量-彈簧模型在儲(chǔ)液容器抗震分析中的應(yīng)用[J].原子能科學(xué)技術(shù),2013,47(6):947-951

[9]李彥民,徐剛,任文敏,等.儲(chǔ)液容器流固耦合動(dòng)力響應(yīng)分析計(jì)算[J].工程力學(xué),2002(4):29-32.

[10]徐剛,任文敏,張維,等.儲(chǔ)液容器的三維流固耦合動(dòng)力特性分析[J].力學(xué)學(xué)報(bào),2004(3):328-335.

[11]陸道綱,張立殷.核電廠中流固耦合現(xiàn)象數(shù)值模擬研究綜述[J].原子能科學(xué)技術(shù),2009,43(S1):76-83.

[12]衛(wèi)媛媛,陸道綱,梁亞平.與流體晃動(dòng)相耦合的快堆主容器地震響應(yīng)的數(shù)值模擬[J].原子能科學(xué)技術(shù),2008,42(S2):472-476.

[13]王建軍,李其漢,朱梓根,等.自由液面大晃動(dòng)的流固耦合數(shù)值分析方法研究進(jìn)展[J].力學(xué)季刊,2001(4):447-454.

[14]GOUDARZI M A.Seismic Behavior of a Single Deck Floating Roof Due to Second Sloshing Mode[J].Journal of Pressure Vessel Technology,2012,135(1):244-252.

[15]陸道綱.CEFR主容器內(nèi)正弦三波激勵(lì)下液面晃動(dòng)響應(yīng)[J].核科學(xué)與工程,2003,23(4):306-310.

[16]衛(wèi)媛媛,陸道綱.基于移動(dòng)粒子法的快堆自由表面流體對(duì)容器頂蓋沖擊現(xiàn)象的數(shù)值模擬.原子能科學(xué)技術(shù),2009,43(10):910-914.

[17]居榮初,曾心傳.彈性結(jié)構(gòu)與液體的耦聯(lián)振動(dòng)理論[M].北京:地震出版社,1983:115-123.