剛度等效在地坑過濾器滯留籃抗震分析中的應用

蘭天寶，董嬌嬌，彭星銘，盛鋒

重大工程裝備

剛度等效在地坑過濾器滯留籃抗震分析中的應用

蘭天寶，董嬌嬌，彭星銘，盛鋒

（中國核電工程有限公司，北京 100840）

對地坑過濾器滯留籃的過濾板進行剛度等效計算。應用能量法和各向異性法對地坑過濾器滯留籃的過濾板進行等效剛度計算?；谧冃文芟嗟鹊哪芰糠ǖ玫降念l率遠大于地震反應譜的截止頻率，該方法適用于滯留籃模型的工程要求?；诟飨虍愋苑ǖ玫降闹黝l率與實際模型主頻率相差14%，該方法同樣適用于滯留籃的工程要求。2種方法均可滿足工程要求，各向異性方法更接近實際模型。最終采用基于能量法的模型和各項異性剛度等效過濾板都取得了較好的效果。

地坑過濾器；滯留籃；過濾板；等效；抗震分析

在壓水堆核電事故的分析中，一些假想的極限事故（如主管道破裂、主蒸汽管道破裂等）如果發(fā)生在安全殼內(nèi)，會引發(fā)安注、安全殼噴淋的動作，最終從主系統(tǒng)噴射出來的水、安注水、噴淋水都會進入安全殼底部，被收集在一個稱為地坑的大水池中[1-3]。為保證應急堆芯冷卻系統(tǒng)的正常運行，需要在安全殼底層設(shè)置地坑過濾器[4]，其主要功能是依靠自身濾網(wǎng)對事故工況下的碎片進行攔截，來保障冷卻系統(tǒng)和噴淋系統(tǒng)的正常運行。壓水堆核電站安全殼地坑過濾器在電廠嚴重事故工況下為防止堆芯融化，保證安全殼的完整性起著關(guān)鍵作用，核電規(guī)范設(shè)計要求對地坑過濾器事故工況下的應力有相應的要求[5-6]。

滯留籃是地坑過濾器的關(guān)鍵設(shè)備之一，破口或事故噴發(fā)的流體夾帶安全殼內(nèi)的碎渣全部流入滯留籃內(nèi)，經(jīng)過滯留籃過濾板的過濾，大部分碎渣留在滯留籃內(nèi)[7]。地坑過濾器滯留籃由鋼結(jié)構(gòu)框架和過濾板構(gòu)成，過濾板為三明治結(jié)構(gòu)，上下2層為帶穿孔的平板，由中間筋板連接固定，過濾板結(jié)構(gòu)如圖1所示。過濾板和鋼結(jié)構(gòu)框架尺寸細節(jié)相差較大，按照過濾板細節(jié)進行建模耗時、費力，因此在評估滯留籃抗震性能時，需要將過濾板進行簡化，等效為平板。采用反應譜法對滯留籃進行抗震分析計算，由于滯留籃浸沒在水中，若過于保守地簡化，較大的附加水質(zhì)量將會使得設(shè)計變得冗余笨重[8-9]，因此在計算時怎樣將過濾板進行簡化使其能逼近實際模型的剛度，是一個難點問題。

圖1 過濾板結(jié)構(gòu)

1 基于變形能等效的能量法

目前關(guān)于剛度等效的方法主要分2種[10]：第1種為平衡方程等效法，即直接利用等效后的板與筋板載荷和位移曲線協(xié)調(diào)一致，適合加筋板的結(jié)構(gòu)等效，工程上有相關(guān)經(jīng)驗公式，適合結(jié)構(gòu)簡單的模型，較為成熟；第2種方法為基于能量法的變形能一致的剛度等效。第2種方法應用廣泛，適合三明治結(jié)構(gòu)類型的等效。首先從變形能等效原理出發(fā)，導出將過濾板等效為無肋板的等效剛度計算公式。

1.1 剛度等效的適用

對于使用剛度等效法應充分考慮其適用情形，剛度等效法是為了簡化結(jié)構(gòu)的一種等效方法，無法保證與原結(jié)構(gòu)完全一致，需要考慮假設(shè)前提對該方法的影響，如等效后的結(jié)構(gòu)變形與原結(jié)構(gòu)相差較大，需要謹慎考慮該方法。如將一個突然變徑的桿進行等效成粗細均勻的桿（如圖2所示），其撓曲線是不同的，導致其臨界力不同，尤其是當變徑桿截面尺寸變化較大的情況下，即底部桿的剛度比上部細桿大很多，等效剛度法帶來的誤差更大[11]。

圖2 桿的等效示意

將過濾板結(jié)構(gòu)等效為單層板，無法展現(xiàn)模型的局部細節(jié)，但對于整體模型來說，并不會影響結(jié)構(gòu)的整體振型和鋼結(jié)構(gòu)應力、基礎(chǔ)載荷等。當需要對過濾板局部應力（如加強筋的應力）進行分析時，可建立局部細節(jié)模型分析。因此，采用剛度等效法適用于滯留籃整體模型的抗震分析。

1.2 能量法計算剛度等效

對于能量法計算等效剛度的研究已經(jīng)十分成熟[12-15]，相關(guān)研究給出了理論和公式，見式（1）、（2）。理論推理過程的基本假設(shè)是：具有相同平面尺寸的加肋板和無肋板在相同荷載的作用下具有相同的彈性變形能，根據(jù)虛功原理，外力所做的功恒等于加肋板變形后板內(nèi)積蓄的變形能a。同理，無肋板在外力所做的功恒等于板變形后板內(nèi)積蓄的變形能b，過程如圖3所示。

圖3 肋板和無肋板的等效

根據(jù)假設(shè)：ab，假設(shè)兩板具有相同的撓曲函數(shù)a(,)=b(,)=(,)，即肋板與無肋板的剛度相同，變形相同。對于a板，變形能由板和肋共同貢獻，見式（3）。

對于等效的b板，其彈性變形能為：

根據(jù)式（5）進一步計算板的頻率：

根據(jù)以上公式，求得等效后板的頻率為118 Hz。經(jīng)分析得，若單方向加肋板，在2個方向的抗彎剛度是不同的，僅通過修改等效板的厚度并不能表征這種特性。本文使用的基于變形能相等的方法平均了2個方向的抗彎剛度，并且假設(shè)了撓曲函數(shù)，而實際的模型變形形態(tài)可能比試函數(shù)更復雜，忽略了剪切變形能，因此會與實際模型存在較大誤差，低于實際模型的頻率[16]。但是根據(jù)等效板厚得出的頻率遠大于地震反應譜的截止頻率，該方法對于滯留籃模型可以滿足工程要求。

2 各向異性剛度矩陣法

彈性體如果在各個方向的彈性都相同，稱為各向同性體[17]；如果在各個方向的彈性并不完全相同，則稱為各向異性體。對于各向異性的薄板研究表明，板的2個方向的尺寸、兩向剛度及扭轉(zhuǎn)剛度對板的頻率均有直接影響[18]。對于加肋板，雖然可以用能量法計算它們的自然頻率和臨界載荷，但計算非常復雜[19]，將加肋條平板中肋條的剛度歸入平板的剛度，從而把具有肋條的各向同性板變換為不具有肋條的各向異性板，則計算大為簡化[20]。

為了進一步逼近實際模型，考慮板的橫向剪切的板殼廣義應力和應變的關(guān)系為：

因此各向異性板的薄膜和彎曲剛度矩陣表示為：

式中：1212為截面力；1212為截面彎矩。

根據(jù)參考文獻[21]，加肋板的抗彎剛度在沿著肋筋的方向可表示為：

式中：Ir為從肋板截面取出的某列（如圖4所示）沿板中心線的慣性矩。另一方向的慣性矩，考慮為三明治層合板，將中間層彈性模量設(shè)為接近0。根據(jù)文獻[21]，各向異性板的扭轉(zhuǎn)剛度H為：

式中：0為無肋時板的抗彎剛度；為剪切模量；為肋肋的抗扭慣性矩。

3 驗證與分析

為了對2種方法進行剛度等效后的計算結(jié)果進行驗證，建立滯留籃過濾板的實際模型[22]。由于實際模型模態(tài)存在局部細節(jié)的肋板振型，但參與質(zhì)量較小，因此選取參與質(zhì)量較大的模態(tài)振型分析，將實際模型得到的計算主頻率與以上2種方法得到的主頻率進行對比，見表1。

表1 主頻率對比

Tab.1 Comparison of main frequencies

實際模型主頻率模態(tài)振型如圖5所示，基于各向異性法對過濾板進行簡化后計算得到的主頻率模態(tài)振型如圖6所示。可以看出，實際模型和各向異性法等效后的模型計算結(jié)果整體趨勢相似，應力最大、最小值出現(xiàn)的位置相似。各向異性法計算的主頻率為263 Hz，相較實際模型計算的主頻率誤差為14%，計算的主頻率均已遠超過地震譜的截止頻率，認為該計算可以符合工程要求，且較基于變形能一致的能量法等效更接近于實際模型計算結(jié)果，該方法更適合在工程計算中用于對滯留籃過濾板進行剛度等效。

由于能量法假定了變形撓曲函數(shù)，忽略了層間剪切能量，考慮綜合等效，因此相對過濾板實際模型的誤差較大。剛度等效將實際幾何模型包含的多層板簡單等效為一層板，會丟失一些模態(tài)振型，如上下兩層板的相離運動等一些振型，如圖7所示。這是由于模型簡化導致的，實際模型的質(zhì)量分布也和等效后的板存在差異。因此，當需要重點關(guān)注過濾板的應力狀態(tài)時[23]，則需要進行局部的模型應力分析，建立肋板符合實際的模型。工程上既要考慮進度更關(guān)注精度，在進行模型等效時，需要評估等效后帶來的誤差。

圖5 實際模型主頻率模態(tài)振型

圖6 各向異性法主頻率模態(tài)振型

圖7 過濾板的模態(tài)振型

滯留籃模型的整體尺寸約為7 000 mm×3 000 mm× 3 000 mm，由于板的剛度較大，對于采用等效板來建模獲取梁的彎矩、軸力、剪力等和基礎(chǔ)載荷是可行的。建立如圖8所示的滯留籃有限模型，應用譜分法[24-25]對滯留籃進行抗震分析。在整體模型中，用各向異性法對過濾板進行剛度等效簡化，過濾板的應力狀態(tài)采用局部模型進行分析。

圖8 滯留籃有限元模型

4 結(jié)論

本文闡述了剛度等效法的原理和適用條件，對比了能量法和各向異性2種方法在計算地坑過濾器滯留籃過濾板的結(jié)果區(qū)別，并根據(jù)實際模型對2種方法進行了驗證，結(jié)果表明：

1）基于變形能相等的能量法和各向異性法對滯留籃過濾板的剛度等效計算均可滿足工程要求，各向異性方法更接近實際模型。

2）由于幾何上的簡化和質(zhì)量的分布會導致等效后的模型與實際模型存在偏差，使用剛度等效法應評估這種影響，最終本文采用各向異性法完成了滯留籃的抗震分析。

[1] 李春, 張慶華, 常猛, 等. 先進壓水堆核電廠安全殼地坑濾網(wǎng)設(shè)計[J]. 核安全, 2012, 11(1): 61-64.

LI Chun, ZHANG Qing-hua, CHANG Meng, et al. The Containment Sump Screen Design in Advanced PWR NPPS[J]. Nuclear Safety, 2012, 11(1): 61-64.

[2] 路廣遙, 任俊生, 向文元, 等. 大型壓水堆核電站地坑過濾器發(fā)展及應用[J]. 電力與能源, 2011, 32(5): 389-391.

LU Guang-yao, REN Jun-sheng, XIANG Wen-yuan, et al. Development and Application of Containment Sump Strainer in PWR[J]. Power & Energy, 2011, 32(5): 389-391.

[3] 張衛(wèi), 龔釗, 朱京梅, 等. 秦山核電二期擴建工程安全殼地坑過濾器設(shè)計改進及遺留問題分析[J]. 核動力工程, 2013, 34(6): 128-131.

ZHANG Wei, GONG Zhao, ZHU Jing-mei, et al. Design Improvement and Open Item Analysis of Containment Sump Strainer in QINSHAN Phase II Extension Project[J]. Nuclear Power Engineering, 2013, 34(6): 128-131.

[4] 劉夏杰, 王德忠, 張臣剛, 等. 核電廠地坑過濾器附加質(zhì)量試驗研究[J]. 核動力工程, 2012, 33(1): 14-18.

LIU Xia-jie, WANG De-zhong, ZHANG Chen-gang, et al. Experimental Study on Added Mass of Sump Strainer in Nuclear Power Plants[J]. Nuclear Power Engineering, 2012, 33(1): 14-18.

[5] 陳俊. 基于變態(tài)模型的核電站新型安全殼地坑過濾器的研發(fā)[D]. 上海: 上海交通大學, 2015.

CHEN Jun. Design of a New Type of Rector Sump Strainer Based on Distorted Model[D]. Shanghai: Shanghai Jiao Tong University, 2015.

[6] 李石磊, 謝洪虎, 何英勇, 等. 安全殼地坑過濾器壓損數(shù)值模擬與試驗研究[J]. 核科學與工程, 2019, 39(6): 867-871.

LI Shi-lei, XIE Hong-hu, HE Ying-yong, et al. Simulation and Experimental Study on Head Loss of Containment Sump Strainer[J]. Nuclear Science and Engineering, 2019, 39(6): 867-871.

[7] 侯建飛, 王慶禮, 司恒遠. 核電廠事故后安全殼內(nèi)置換料水箱碎片傳輸性能分析[J]. 核動力工程, 2021, 42(5): 76-80.

HOU Jian-fei, WANG Qing-li, SI Heng-yuan. Analysis of Post-Accident Debris Transport Performance in In-Containment Refueling Water Storage Tank in Nuclear Power Plants[J]. Nuclear Power Engineering, 2021, 42(5): 76-80.

[8] 于肖宇. 核電站地坑過濾器附加質(zhì)量的研究[D]. 上海: 上海交通大學, 2010.

YU Xiao-yu. Study of Added Mess the Containment Sump Strainer in Nuclear Plant[D]. Shanghai: Shanghai Jiao Tong University, 2010.

[9] 謝洪虎, 李石磊, 張峰, 等. 核電廠安全殼內(nèi)置換料水箱過濾系統(tǒng)過濾性能及阻力特性研究[J]. 核動力工程, 2019, 40(6): 130-134.

XIE Hong-hu, LI Shi-lei, ZHANG Feng, et al. Assessment of Filtering and Pressure Drop Performances of IRWST Sump Strainers for Nuclear Power Plants[J]. Nuclear Power Engineering, 2019, 40(6): 130-134.

[10] SERTI? J, KOZAK D, SAMARD?I? I. Calculation of Reaction Forces in the Boiler Supports Using the Method of Equivalent Stiffness of Membrane Wall[J]. The Scientific World Journal, 2014, 2014: 392048.

[11] 黃義仿, 賈汝民. 關(guān)于等效剛度法適用條件的討論[J]. 力學與實踐, 1996, 18(2): 59-60.

HUANG Yi-fang, JIA Ru-min. Discussion on the Applicable Conditions of Equivalent Stiffness Method[J]. Mechanics and Practice, 1996, 18(2): 59-60.

[12] 張煜, 包太, 秦永留. 加肋板在橫向載荷和板中面力共同作用下的撓度分析[J]. 貴州大學學報(自然科學版), 2009, 26(1): 111-114.

ZHANG Yu, BAO Tai, QIN Yong-liu. Analysis of Deflection of Stiffened Plates with Rids under Transverse Load Together with the Neutral Plane Forces[J]. Journal of Guizhou University (Natural Science Edition), 2009, 26(1): 111-114.

[13] 丁婷, 丁圣果, 姜宇. 基于能量原理的肋型板等效剛度[J]. 西安建筑科技大學學報(自然科學版), 2013, 45(6): 858-863.

DING Ting, DING Sheng-guo, JIANG Yu. Equivalent Stiffness of Rib Board Based on the Energy Principle[J]. Journal of Xi’an University of Architecture & Technology (Natural Science Edition), 2013, 45(6): 858-863.

[14] 李綺文, 丁婷, 丁圣果. 計算肋形板自振頻率的能量法[J]. 山西建筑, 2012, 38(26): 33-34.

LI Qi-wen, DING Ting, DING Sheng-guo. Energy Method of Calculating Natural Vibration Frequency of Rib Plate[J]. Shanxi Architecture, 2012, 38(26): 33-34.

[15] 孟維迎, 孫鵬文, 張?zhí)m挺, 等. 基于能量法的復合材料層合板性能等效算法研究與應用[J]. 太陽能學報, 2014, 35(8): 1505-1510.

MENG Wei-ying, SUN Peng-wen, ZHANG Lan-ting, et al. Research and Application of Property Equivalence Algorithm for Composite Lamination Plate Based on Energy Method[J]. Acta Energiae Solaris Sinica, 2014, 35(8): 1505-1510.

[16] 吳曉. 用改進能量法計算梁屈曲荷載和固有頻率[J]. 工程與試驗, 2019, 59(2): 1-3.

WU Xiao. Calculation of the Buckling Loads and Natural Frequency of Beams by Improved Energy Method[J]. Engineering & Test, 2019, 59(2): 1-3.

[17] 徐芝綸. 彈性力學(下)[M]. 北京: 高等教育出版社, 1980.

XU Zhi-lun. Elastic Mechanics (II)[M]. Beijing: Higher Education Press, 1980.

[18] 丁圣果. 分析結(jié)構(gòu)力學[M]. 貴陽: 貴州科技出版社, 2010.

DING Sheng-guo. Mechanics of Analysis Structural[M]. Guiyang: Guizhou Science and Technology Publishing Press, 2010.

[19] 丁圣果, 劉廣寧, 丁綺文. 肋型樓蓋的一種有限元算法及其內(nèi)力分布[J]. 建筑結(jié)構(gòu), 1999(6): 46-50.

DING Sheng-guo, LIU Guang-ning, LI Qi-wen. A Finite Element Method of Ribbed Floor System and Its Internal Distribution[J]. Building Structure, 1999(6): 46-50.

[20] 楊章, 童根樹, 張磊. 對稱布置2根單側(cè)加勁肋的有效剛度[J]. 浙江大學學報(工學版), 2016, 50(8): 1446-1455.

YANG Zhang, TONG Gen-shu, ZHANG Lei. Effective Rigidity of Two One-Side Stiffeners Arranged Symmetrically[J]. Journal of Zhejiang University (Engineering Science), 2016, 50(8): 1446-1455.

[21] BLEVINS R D. Formulas for natural frequency and mode shape[M]. New York: Van Nostrand Reinhold Co, 1979.

[22] 謝貽權(quán), 何福保. 彈性和塑性力學中的有限單元法[M]. 北京: 機械工業(yè)出版社, 2000.

XIE Yi-quan, HE Fu-bao. Finite Element Method of Elasticity and Plastic Mechanics[M]. Beijing: China Machine Press, 2000.

[23] 陳學潮, 劉永仁, 金瑞椿. 加肋管板的變形和應力分析[J]. 上海力學, 1980, 1(3): 13-27.

CHEN Xue-chao, LIU Yong-ren, JIN Rui-chun. The Analysis of Deformations and Stresses of Ribbed Tube-Plate[J]. Chinese Quarterly of Mechanics, 1980, 1(3): 13-27.

[24] 錢浩, 鄧晶晶, 徐定耿, 等. 安全殼再循環(huán)地坑濾網(wǎng)應力分析[J]. 力學季刊, 2007, 28(1): 165-169.

QIAN Hao, DENG Jing-jing, XU Ding-geng, et al. Stress Analysis of Recirculating Pit Filter in Containment[J]. Chinese Quarterly of Mechanics, 2007, 28(1): 165-169.

[25] 朱賀. 循環(huán)水浸泡影響下的新型地坑過濾器抗震分析[J]. 核安全, 2014, 13(4): 90-94.

ZHU He. Anti-Seismic Analysis on a New Type of Sump Filter in Circulating Water Immersion Condition[J]. Nuclear Safety, 2014, 13(4): 90-94.

Application of Stiffness Equivalent In Seismic Analysis of Containment Sump Filter Retention Basket

LAN Tian-bao, DONG Jiao-jiao, PENG Xing-ming, SHENG Feng

(China Nuclear Power Engineering Co., Ltd., Beijing 100840, China)

The work aims to calculate the stiffness equivalent of the filter plate of the containment sump filter retention basket.Energy method and anisotropy method were used to calculate the equivalent stiffness of the filter plate of the containment sump filter retention basket. The frequency obtained by the energy method based on equal deformation energy was much larger than the cutoff frequency of seismic response spectrum, which was suitable for the engineering requirements of the stranded basket model. There was a 14% difference between the principal frequency obtained by anisotropy method and that of the actual model. Both methods can meet the engineering requirements, and the anisotropic method is closer to the actual model. Finally, the model based on the energy method and the anisotropic stiffness equivalent filter plates has good results.

containment sump filter, retention basket, filter plate, equivalent, seismic analysis

2021-04-13；

2022-02-14

LAN Tian-bao (1990-), Male, Master.

盛鋒（1980—），男。

Corresponding author：SHENG Feng (1980-), Male.

蘭天寶, 董嬌嬌, 彭星銘, 等. 剛度等效在地坑過濾器滯留籃抗震分析中的應用[J]. 裝備環(huán)境工程, 2023, 20(2): 096-101.

TL362

A

1672-9242(2023)02-0096-06

10.7643/ issn.1672-9242.2023.02.013

2021–04–13；

2022–02–14

蘭天寶（1990—），男，碩士。

LAN Tian-bao, DONG Jiao-jiao, PENG Xing-ming, et al. Application of Stiffness Equivalent In Seismic Analysis of Containment Sump Filter Retention Basket[J]. Equipment Environmental Engineering, 2023, 20(2): 096-101.

責任編輯：劉世忠