1、 装 备 环 境 工 程 第 20 卷 第 2 期 96 EQUIPMENT ENVIRONMENTAL ENGINEERING 2023 年 2 月 收稿日期：20210413；修订日期：20220214 Received：2021-04-13；Revised：2022-02-14 作者简介：兰天宝（1990），男，硕士。Biography：LAN Tian-bao(1990-),Male,Master.通讯作者：盛锋（1980），男。Corresponding author：SHENG Feng(1980-),Male.引文格式：兰天宝,董娇娇,彭星铭,等.刚度等效在地坑过滤器滞留篮抗震分
2、析中的应用J.装备环境工程,2023,20(2):096-101.LAN Tian-bao,DONG Jiao-jiao,PENG Xing-ming,et al.Application of Stiffness Equivalent In Seismic Analysis of Containment Sump Filter Retention BasketJ.Equipment Environmental Engineering,2023,20(2):096-101.刚度等效在地坑过滤器滞留篮抗震 分析中的应用 兰天宝，董娇娇，彭星铭，盛锋（中国核电工程有限公司，北京 100840）摘要：
3、目的目的 对地坑过滤器滞留篮的过滤板进行刚度等效计算。方法方法 应用能量法和各向异性法对地坑过滤器滞留篮的过滤板进行等效刚度计算。结果结果 基于变形能相等的能量法得到的频率远大于地震反应谱的截止频率，该方法适用于滞留篮模型的工程要求。基于各向异性法得到的主频率与实际模型主频率相差 14%，该方法同样适用于滞留篮的工程要求。结论结论 2 种方法均可满足工程要求，各向异性方法更接近实际模型。最终采用基于能量法的模型和各项异性刚度等效过滤板都取得了较好的效果。关键词：地坑过滤器；滞留篮；过滤板；等效；抗震分析 中图分类号：TL362 文献标识码：A 文章编号：1672-9242(2023)02-00
4、96-06 DOI：10.7643/issn.1672-9242.2023.02.013 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)ABSTRACT:The work aims to calculate the
5、 stiffness equivalent of the filter plate of the containment sump filter retention bas-ket.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 equa
6、l 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 metho
7、ds 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.KEY WORDS:containment sump filter,retention basket,filter plate,equivalent,seismic an
8、alysis 在压水堆核电事故的分析中，一些假想的极限事故（如主管道破裂、主蒸汽管道破裂等）如果发生在重大工程装备 第 20 卷 第 2 期 兰天宝，等：刚度等效在地坑过滤器滞留篮抗震分析中的应用 97 安全壳内，会引发安注、安全壳喷淋的动作，最终从主系统喷射出来的水、安注水、喷淋水都会进入安全壳底部，被收集在一个称为地坑的大水池中1-3。为保证应急堆芯冷却系统的正常运行，需要在安全壳底层设置地坑过滤器4，其主要功能是依靠自身滤网对事故工况下的碎片进行拦截，来保障冷却系统和喷淋系统的正常运行。压水堆核电站安全壳地坑过滤器在电厂严重事故工况下为防止堆芯融化，保证安全壳的完整性起着关键作用，核电规
9、范设计要求对地坑过滤器事故工况下的应力有相应的要求5-6。滞留篮是地坑过滤器的关键设备之一，破口或事故喷发的流体夹带安全壳内的碎渣全部流入滞留篮内，经过滞留篮过滤板的过滤，大部分碎渣留在滞留篮内7。地坑过滤器滞留篮由钢结构框架和过滤板构成，过滤板为三明治结构，上下 2 层为带穿孔的平板，由中间筋板连接固定，过滤板结构如图 1 所示。过滤板和钢结构框架尺寸细节相差较大，按照过滤板细节进行建模耗时、费力，因此在评估滞留篮抗震性能时，需要将过滤板进行简化，等效为平板。采用反应谱法对滞留篮进行抗震分析计算，由于滞留篮浸没在水中，若过于保守地简化，较大的附加水质量将会使得设计变得冗余笨重8-9，因此在计
10、算时怎样将过滤板进行简化使其能逼近实际模型的刚度，是一个难点问题。图 1 过滤板结构 Fig.1 Structure diagram of filter plate 1 基于变形能等效的能量法 目前关于刚度等效的方法主要分 2 种10：第 1种为平衡方程等效法，即直接利用等效后的板与筋板载荷和位移曲线协调一致，适合加筋板的结构等效，工程上有相关经验公式，适合结构简单的模型，较为成熟；第 2 种方法为基于能量法的变形能一致的刚度等效。第 2 种方法应用广泛，适合三明治结构类型的等效。首先从变形能等效原理出发，导出将过滤板等效为无肋板的等效刚度计算公式。1.1 刚度等效的适用 对于使用刚度等效法应
11、充分考虑其适用情形，刚度等效法是为了简化结构的一种等效方法，无法保证与原结构完全一致，需要考虑假设前提对该方法的影响，如等效后的结构变形与原结构相差较大，需要谨慎考虑该方法。如将一个突然变径的杆进行等效成粗细均匀的杆（如图 2 所示），其挠曲线是不同的，导致其临界力不同，尤其是当变径杆截面尺寸变化较大的情况下，即底部杆的刚度比上部细杆大很多，等效刚度法带来的误差更大11。图 2 杆的等效示意 Fig.2 Equivalent diagram of rod 将过滤板结构等效为单层板，无法展现模型的局部细节，但对于整体模型来说，并不会影响结构的整体振型和钢结构应力、基础载荷等。当需要对过滤板局部应
12、力（如加强筋的应力）进行分析时，可建立局部细节模型分析。因此，采用刚度等效法适用于滞留篮整体模型的抗震分析。1.2 能量法计算刚度等效 对于能量法计算等效刚度的研究已经十分成熟12-15，相关研究给出了理论和公式，见式（1）、（2）。理论推理过程的基本假设是：具有相同平面尺寸 ab的加肋板和无肋板在相同荷载 q 的作用下具有相同的弹性变形能，根据虚功原理，外力所做的功恒等于加肋板变形后板内积蓄的变形能 Ua。同理，无肋板在外力 q 所做的功恒等于板变形后板内积蓄的变形能 Ub，过程如图 3 所示。图 3 肋板和无肋板的等效 Fig.3 Equivalent diagram of ribbed
13、slab and unribbed slab 98 装 备 环 境 工 程 2023 年 2 月 ()aad dqWx,yx y=U (1)()bbd dqWx,yx y=U(2)根据假设：Ua=Ub，假设两板具有相同的挠曲函数 Wa(x,y)=Wb(x,y)=W(x,y)，即肋板与无肋板的刚度相同，变形相同。对于 a 板，变形能由板和肋共同贡献，见式（3）。222aa2200222220000ddd d2dd11dd22bammbbiiyyiiDWWUx yxyWWEIyGIyx yy=+|+|(3)式中：Da为 a 板的抗弯刚度，()3a aa212 1E tD=-；EIy为一根肋的抗弯刚
14、度，3a12yE bhEI=（矩形截面，b为宽度，h 为高度）；GIy为一根肋的扭转刚度。对于等效的 b 板，其弹性变形能为：2bb00ddd ddd22ba22DWWU=x y2xy+|(4)式中：Db为 b 板的抗弯刚度，()3b bb212 1E tD=-。最终根据变形能相等，便可求出等效后板的厚度 tb。过滤板的尺寸 a=840 mm，b=720 mm，肋板高 h=60 mm，肋宽 w=4 mm，10 个肋等间距分布在 x 方向。取挠曲函数()sincosxyW x,yAab=，最终得到等效厚度为 13 mm。根据式（5）进一步计算板的频率：22211Df=mab+|(5)式中：a 和
15、 b 为板的长宽尺寸；D 为板的抗弯刚度；m为板的均布质量。根据以上公式，求得等效后板的频率为 118 Hz。经分析得，若单方向加肋板，在 2 个方向的抗弯刚度是不同的，仅通过修改等效板的厚度并不能表征这种特性。本文使用的基于变形能相等的方法平均了 2 个方向的抗弯刚度，并且假设了挠曲函数，而实际的模型变形形态可能比试函数更复杂，忽略了剪切变形能，因此会与实际模型存在较大误差，低于实际模型的频率16。但是根据等效板厚得出的频率远大于地震反应谱的截止频率，该方法对于滞留篮模型可以满足工程要求。2 各向异性刚度矩阵法 弹性体如果在各个方向的弹性都相同，称为各向同性体17；如果在各个方向的弹性并不完
16、全相同，则称为各向异性体。对于各向异性的薄板研究表明，板的 2 个方向的尺寸、两向刚度及扭转刚度对板的频率均有直接影响18。对于加肋板，虽然可以用能量法计算它们的自然频率和临界载荷，但计算非常复杂19，将加肋条平板中肋条的刚度归入平板的刚度，从而把具有肋条的各向同性板变换为不具有肋条的各向异性板，则计算大为简化20。为了进一步逼近实际模型，考虑板的横向剪切的板壳广义应力和应变的关系为：=D(6)mmbbssDDD=DDD|=|(7)式中：Dm、Db、Ds分别为薄膜、弯曲和剪切刚度，直接影响薄膜和弯曲应力。对于各向同性材料，mD=2101011002Et|-|-|，2bm12tDD=，s2(1)EtDk=+。各向异性的材料，可直接在各向同性的基础上修改另一方向的属性。因此各向异性板的薄膜和弯曲刚度矩阵表示为：11121122223312121112112222331212000000000000sym0AANANANDDMDMDM|=|(8)式中：N1、N2、N12为截面力；M1、M2、M12为截面弯矩。根据参考文献21，加肋板的抗弯刚度在沿着肋筋的方向可表示为：r11EIDa=(9)式