• 匿名盲审
  • 学术期刊非营利性
  • 全球免费开放获取全文
  • 最新科研成果提供绿色通道

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

核电厂非能动安全壳冷却系统储水箱的振动特性分析

李晓蒙 何铮

李晓蒙, 何铮. 核电厂非能动安全壳冷却系统储水箱的振动特性分析[J]. 南方能源建设, 2015, 2(4): 102-106. doi: 10.16516/j.gedi.issn2095-8676.2015.04.018
引用本文: 李晓蒙, 何铮. 核电厂非能动安全壳冷却系统储水箱的振动特性分析[J]. 南方能源建设, 2015, 2(4): 102-106. doi: 10.16516/j.gedi.issn2095-8676.2015.04.018
Xiaomeng LI, Zheng HE. Dynamic Characteristic Analysis of the Passive Containment Cooling System Water Tank of Nuclear Power Plant[J]. SOUTHERN ENERGY CONSTRUCTION, 2015, 2(4): 102-106. doi: 10.16516/j.gedi.issn2095-8676.2015.04.018
Citation: Xiaomeng LI, Zheng HE. Dynamic Characteristic Analysis of the Passive Containment Cooling System Water Tank of Nuclear Power Plant[J]. SOUTHERN ENERGY CONSTRUCTION, 2015, 2(4): 102-106. doi: 10.16516/j.gedi.issn2095-8676.2015.04.018

核电厂非能动安全壳冷却系统储水箱的振动特性分析

doi: 10.16516/j.gedi.issn2095-8676.2015.04.018
基金项目: 

国家大型先进压水堆核电站重大专项关键设备设计分析技术研究课题 2012ZX06004-012

详细信息
    作者简介:

    李晓蒙(1988),女,北京人,工程师,硕士,主要从事核反应堆结构力学计算等(e-mail)lxmcau@126.com。

  • 中图分类号: TL364

Dynamic Characteristic Analysis of the Passive Containment Cooling System Water Tank of Nuclear Power Plant

  • 摘要: PCS储水箱是AP系列核电站非能动安全壳冷却系统的重要组成部分。水体和水箱的耦合作用会影响结构的整体性能,并可能造成水箱破坏,因此有必要在考虑流固耦合前提下对储水箱进行振动分析。首先给出了研究水体晃动的理论模型,得到了储水量分别为60%、70%、80%时水箱的1阶固有频率。之后运用有限元软件ANSYS对空水箱模型及含三种不同储水量的水箱模型进行了模态分析,研究了储水量对结构振动特性的影响,并将数值仿真结果与理论计算结果进行了对比。最后,对空箱及含60%储水量的水箱模型进行了地震反应谱分析,结果显示水体由于固有频率低,在地震载荷下易被激起,在设计中必须加以考虑。
  • 图  1  Housner理论计算模型示意图

    Fig.  1  Housner's Theoretical Calculation Model

    图  2  等式(5)的前四阶根与Ri/Ro的关系图

    Fig.  2  The Diagram of the Relationship Between the Previous Four Roots and Ri/Ro of Equation (5)

    图  3  等效深度示意图

    Fig.  3  The Equivalent Depth Schematic

    图  4  空箱模型和60%储水箱模型示意图

    Fig.  4  Models of Empty Tank and 60% Capacity Tank

    图  5  约束后的水箱模型图

    Fig.  5  PCS Tank Model with All Restrictions

    图  6  空箱和60%储水箱模态图

    Fig.  6  The First and Fifth Mode Figures of Water

    图  7  60%、70%、80%储水箱的模态曲线图

    Fig.  7  Frequency Diagram of 60%, 70% and 80% Capacity Tank

    图  8  空箱模型和60%储水模型水体位移图

    Fig.  8  Water Displacement of Empty Tank and 60% Capacity Tank

    表  1  理论计算结果

    Tab.  1.   Results of Theoretical Calculations

    贮水量 60% 70% 80%
    ha/ m 4.8 5.6 6.4
    f / Hz 0.125 0.133 0.140
    下载: 导出CSV

    表  2  固有频率表

    Tab.  2.   Table of the Natural Frequency

    阶数 空箱模型频率 60%储水箱模型频率
    1 27.687 0.120
    2 27.687 0.120
    3 31.905 0.160
    4 31.905 0.208
    5 36.475 0.209
    6 36.475 0.277
    7 39.948 0.277
    8 39.948 0.324
    9 41.197 0.325
    10 41.197 0.356
    11 44.371 0.356
    12 44.371 0.361
    13 46.611 0.361
    14 46.611 0.374
    15 46.800 0.374
    16 46.800 0.391
    17 53.015 0.392
    18 53.015 0.394
    19 54.586 0.399
    20 54.586 0.399
    下载: 导出CSV

    表  3  三种模型固有频率表

    Tab.  3.   The Natural Frequency of Three Models

    阶数 60%储水箱模型频率/ Hz 70%储水箱模型频率/ Hz 80%储水箱模型频率/ Hz
    1 0.120 0.129 0.137
    2 0.120 0.129 0.137
    3 0.160 0.177 0.186
    4 0.208 0.220 0.228
    5 0.209 0.220 0.228
    6 0.277 0.284 0.287
    7 0.277 0.284 0.287
    8 0.324 0.327 0.327
    9 0.325 0.328 0.329
    10 0.356 0.357 0.359
    下载: 导出CSV

    表  4  仿真结果与理论计算对比表

    Tab.  4.   Comparison Between Simulation Results and the Theoretical Results

    储水量 60% 70% 80%
    fr 0.125 0.133 0.140
    fs 0.120 0.129 0.137
    偏差σ 4% 3% 2%
    下载: 导出CSV

    表  5  楼层反应谱

    Tab.  5.   Floor Response Spectrum

    频率 / Hz 0.100 0.135 0.138 0.174 0.522 1.035 100
    谱值 / g 0.023 0.043 0.045 0.087 0.299 0.573 1.024
    下载: 导出CSV
  • [1] 林诚格,郁祖盛. 非能动安全先进压水堆核电技术(中册) [M]. 北京:原子能出版社,2010, 509-510.
    [2] 程旭东. 大型储罐的抗震可靠性分析及其震害预测 [D]. 北京:中国石油大学,2010:4-5.
    [3] 高晓安,陆道纲,朱玉巧,等.考虑流固耦合的储液容器地震反应计算 [J]. 工程抗震,2001, 9(3):41-44.
    [4] 徐刚,任文敏,张维,等. 储液容器的三维流固耦合动力特性分析 [J]. 力学学报,2004,36(3):328-335.
    [5] 赵冯. 两种圆形储液容器的动力学分析 [D]. 西安:西安电子科技大学,2013.
    [6] Housner G W. Dynamic Pressure on Accelerated Fluid Containers [J]. Bull.Seism.Soc.Am, 1957, 47(1):15-35.
    [7] 鲁丽,杨翊仁,金建明. 反应堆吊篮在空气和静水中的振动特性分析[J]. 西南交通大学学报,2004, 39(1):82-85.
    [8] ESMAILZADEH M, LAKIS A A, THOMAS M. Three-dimensional Modeling of Curved Structures Containing and/or Submerged in Fluid [J]. Finite Elements in Analysis and Design, 2008(44): 334-345.
  • 加载中
图(8) / 表 (5)
计量
  • 文章访问数:  148
  • HTML全文浏览量:  43
  • PDF下载量:  11
  • 被引次数: 0
出版历程
  • 收稿日期:  2015-11-14
  • 刊出日期:  2015-12-25

核电厂非能动安全壳冷却系统储水箱的振动特性分析

doi: 10.16516/j.gedi.issn2095-8676.2015.04.018
    基金项目:

    国家大型先进压水堆核电站重大专项关键设备设计分析技术研究课题 2012ZX06004-012

    作者简介: 作者简介:李晓蒙(1988),女,北京人,工程师,硕士,主要从事核反应堆结构力学计算等(e-mail)lxmcau@126.com。

  • 中图分类号: TL364

摘要: PCS储水箱是AP系列核电站非能动安全壳冷却系统的重要组成部分。水体和水箱的耦合作用会影响结构的整体性能,并可能造成水箱破坏,因此有必要在考虑流固耦合前提下对储水箱进行振动分析。首先给出了研究水体晃动的理论模型,得到了储水量分别为60%、70%、80%时水箱的1阶固有频率。之后运用有限元软件ANSYS对空水箱模型及含三种不同储水量的水箱模型进行了模态分析,研究了储水量对结构振动特性的影响,并将数值仿真结果与理论计算结果进行了对比。最后,对空箱及含60%储水量的水箱模型进行了地震反应谱分析,结果显示水体由于固有频率低,在地震载荷下易被激起,在设计中必须加以考虑。

English Abstract

李晓蒙, 何铮. 核电厂非能动安全壳冷却系统储水箱的振动特性分析[J]. 南方能源建设, 2015, 2(4): 102-106. doi: 10.16516/j.gedi.issn2095-8676.2015.04.018
引用本文: 李晓蒙, 何铮. 核电厂非能动安全壳冷却系统储水箱的振动特性分析[J]. 南方能源建设, 2015, 2(4): 102-106. doi: 10.16516/j.gedi.issn2095-8676.2015.04.018
Xiaomeng LI, Zheng HE. Dynamic Characteristic Analysis of the Passive Containment Cooling System Water Tank of Nuclear Power Plant[J]. SOUTHERN ENERGY CONSTRUCTION, 2015, 2(4): 102-106. doi: 10.16516/j.gedi.issn2095-8676.2015.04.018
Citation: Xiaomeng LI, Zheng HE. Dynamic Characteristic Analysis of the Passive Containment Cooling System Water Tank of Nuclear Power Plant[J]. SOUTHERN ENERGY CONSTRUCTION, 2015, 2(4): 102-106. doi: 10.16516/j.gedi.issn2095-8676.2015.04.018
  • AP系列核电站是美国西屋公司开发的第三代非能动型压水反应堆核电站,具有相对简单性、安全性和可靠性等特点。其中非能动安全壳冷却系统(Passive Containment Cooling System, PCS)是其主要设计特点之一。非能动安全壳冷却系统由一台与安全壳屏蔽构筑物结构合为一体的冷却水储水箱(以下简称PCS储水箱)、管道以及相关的仪表、阀门等构成。PCS储水箱位于屏蔽构筑物顶部,除了安全壳除热功能外、还起到了代替乏燃料水池以及在地震安全停堆后为消防用储水箱补水等重要作用[1]

    由于PCS储水箱具有举足轻重的地位,在设计上按安全停堆地震(SSE)等级,抗震I类设计[1]。然而,在地震作用下,水箱中水体的振动与固体的振动不同,水箱内带有自由表面的水体不但随着水箱水平运动而且由于惯性和重力作用,其表面会出现升降而发生晃荡。同时,水体晃荡会对限制其运动的固体壁面产生晃荡力,从而反过来影响结构的运动。因此,PCS储水箱在地震作用下的动力响应非常复杂,它不仅受到结构的振动和水体的晃荡影响,还必须考虑结构与水体的相互作用。而工程上有很多实例都说明了:地震作用导致水体晃动而产生的晃荡力是导致储水结构破坏的主要原因[2-4]。因此,研究PCS水箱的振动特点及在地震下的动力响应,是十分必要的。

    文中首先给出了描述水箱流体-固体耦合振动问题的简化理论模型,之后运用大型有限元软件ANSYS对PCS储水箱的振动模态进行了计算,分析了不同储水量对水箱固有频率的影响,并与理论计算结果进行了对比。最后对PCS储水箱进行了抗震反应谱分析,得到了水体的晃动作用对水箱结构整体的抗震性能有很大影响,在设计中必须将其考虑在内的结论。

    • 流体动力学是分析储液容器内液体晃动问题的理论基础,一般的理论方法是对Navier-Stokes方程直接进行理论分析。然而理论分析求解往往比较复杂,且对于形状不规则的流动或是非线性问题不能合理简化,大多数都不能得到解析结果[5]。为了规避求解Navier-Stokes方程,Housner等人提出了简化的集中质量法[6]。这一方法将液动压力分为两部分,一部分是随容器作同步运动的液体惯性力产生的脉冲压力,另一部分是由容器内液体自由晃动而产生的对流压力。Housner将流体部分简化成不同的质量单元和弹簧单元,如图1所示。

      图  1  Housner理论计算模型示意图

      Figure 1.  Housner's Theoretical Calculation Model

      图1中脉冲质量M1,固接在箱体上;第一阶振动等效质量M2弹性连接箱体上;k1为第一阶等效弹簧刚度;h0为等效质量M1的作用位置;h1为等效质量M2的作用位置;f1为储罐中液体晃动的第一阶模态频率。对于圆柱形储液容器有:

      ((1))

      但是PCS储水箱并不是规则的圆柱形,因此需要进行近似处理。液体各阶模态固有频率的理论计算公式见式(2)~式(4):

      ((2))
      ((3))
      ((4))

      式中:n为流体振动模态的阶数;g为重力加速度;νn是下面等式(5)的第n个根。

      式(5)前四阶根与Ri/Ro的函数关系曲线如图2所示,可得v1=1.397 [5]

      图  2  等式(5)的前四阶根与Ri/Ro的关系图

      Figure 2.  The Diagram of the Relationship Between the Previous Four Roots and Ri/Ro of Equation (5)

      ((5))

      式中:J1(x)为第一类一阶贝塞尔函数,Y1(x)为第二类一阶贝塞尔函数。

      流体的等效流体深度ha指:将同体积的流体充入截面为矩形且有与水箱内径、外径相同的环形水箱后的流体深度,如图3所示。

      图  3  等效深度示意图

      Figure 3.  The Equivalent Depth Schematic

      设PCS储水箱外半径Ro为11.58 m,内半径Ri为5.3 m,水箱内侧壁高5.96 m,外侧壁高10.68 m,壁厚0.9 m,水箱体积约为2 672 m3。为考察贮水容量对液体频率的影响,分别计算贮水量为60%、70%及80%时液体一阶固有频率。得到的理论计算结果见表1

      表 1  理论计算结果

      Table 1.  Results of Theoretical Calculations

      贮水量 60% 70% 80%
      ha/ m 4.8 5.6 6.4
      f / Hz 0.125 0.133 0.140
    • 虽然Housner简化处理方法已被国内、外学者广泛采用,但是这种方法具有一定的局限性,并未能真正考虑流-固之间的动力耦合关系。而ANSYS有限元程序在数值模拟流-固耦合方面具有较强的实力[8]。本文将采用ANSYS程序对PCS储水箱的振动模态进行计算分析。

    • 分别建立空箱及储水量为60%的水箱有限元模型。为了考察水体晃动对结构整体振动特性的影响,在空箱模型中加入了60%的静水质量。

      箱内液体假设为无旋、不可压缩的理想流体,采用Fluid80单元来模拟;箱体材料为混凝土,采用shell63单元进行模拟。水箱的材料特性取值如下:弹性模量E = 24.9 GPa,泊松比μ = 0.17,密度ρ = 2 432 kg/m3

      一般取流体的弹性模量为结构弹性模量的1%,泊松比取相对很小的值[7]。基于这一原则,水体的材料特性在常温(20 ℃)下取值如下:弹性模量E = 2.1 GPa;泊松比μ = 0.3×10-4;密度ρ = 1 000 kg/m3,热膨胀系数α = 5.5×10-6K-1。两种模型的剖面图见图4

      图  4  空箱模型和60%储水箱模型示意图

      Figure 4.  Models of Empty Tank and 60% Capacity Tank

      对结构和液体共用节点或临近节点进行自由度耦合,具体规则为:水箱底部共用节点Z方向自由度耦合,水箱壁面临近节点法向自由度耦合,切向自由度自由。水体上表面为自由液面,下表面为固定壁条件。

      选取水箱节点的Ux自由度为主自由度,选取水体自由表面处的Uz自由度为主自由度。施加约束后的水箱模型如图5所示。

      图  5  约束后的水箱模型图

      Figure 5.  PCS Tank Model with All Restrictions

    • 采用ANSYS模态分析中的“缩减法”进行计算,分别得到两种模型的前20阶固有频率,如下表2所示。

      表 2  固有频率表

      Table 2.  Table of the Natural Frequency

      阶数 空箱模型频率 60%储水箱模型频率
      1 27.687 0.120
      2 27.687 0.120
      3 31.905 0.160
      4 31.905 0.208
      5 36.475 0.209
      6 36.475 0.277
      7 39.948 0.277
      8 39.948 0.324
      9 41.197 0.325
      10 41.197 0.356
      11 44.371 0.356
      12 44.371 0.361
      13 46.611 0.361
      14 46.611 0.374
      15 46.800 0.374
      16 46.800 0.391
      17 53.015 0.392
      18 53.015 0.394
      19 54.586 0.399
      20 54.586 0.399

      从频率表中可看出,由于空箱模型分别关于XOZ平面和YOZ平面对称,故空箱模型的每阶模态在XY方向上会出现重频现象。而有60%储水量的水箱模型由于加入了流体的关系,固有频率降低很多,同样,由于模型的对称性,也逐渐出现了重频现象。

      根据频率表,绘出两种模型的频率曲线图,如图6所示。可看出,空箱模型的固有频率的增加趋势呈线性,而60%储水箱模型的固有频率的增加趋势呈二次曲线型。

      图  6  空箱和60%储水箱模态图

      Figure 6.  The First and Fifth Mode Figures of Water

      上图为60%储水箱模型中水体的1阶、5阶模态振型图。在前20阶振型中,水箱结构中没有发生明显的变形,振型全部为流体的晃动。同时还发现,水体的z方向的晃动模态比较密集且频率较小,所以地震时储液的晃动比较容易被激发,这也是导致水箱顶部失效的主要原因。

    • 为了考虑水箱内储水量对水箱自振频率的影响,分别计算储水量为60%、70%以及80%时,水箱的前20阶自振频率。由于篇幅有限,这里仅列出三种储水量模型的前10阶频率,如表3所示。

      表 3  三种模型固有频率表

      Table 3.  The Natural Frequency of Three Models

      阶数 60%储水箱模型频率/ Hz 70%储水箱模型频率/ Hz 80%储水箱模型频率/ Hz
      1 0.120 0.129 0.137
      2 0.120 0.129 0.137
      3 0.160 0.177 0.186
      4 0.208 0.220 0.228
      5 0.209 0.220 0.228
      6 0.277 0.284 0.287
      7 0.277 0.284 0.287
      8 0.324 0.327 0.327
      9 0.325 0.328 0.329
      10 0.356 0.357 0.359

      将上述三种储水量的模态频率综合对比分析,见图7,可以看出,随着储水量的增多,各阶对应的模态频率一直变大。不同储水量下,前8阶的固有频率相差较大,之后的高阶频率差距逐渐减小并越来越接近。

      图  7  60%、70%、80%储水箱的模态曲线图

      Figure 7.  Frequency Diagram of 60%, 70% and 80% Capacity Tank

    • 将仿真计算的三种不同储水量的1阶固有频率与理论计算结果进行对比,记理论计算结果为fr,仿真计算结果为fs,两者对比如表4所示。其中fsfr的偏差计算公式为:

      表 4  仿真结果与理论计算对比表

      Table 4.  Comparison Between Simulation Results and the Theoretical Results

      储水量 60% 70% 80%
      fr 0.125 0.133 0.140
      fs 0.120 0.129 0.137
      偏差σ 4% 3% 2%
      ((6))

      通过表4发现,数值仿真与理论计算结果相差不多,偏差在可接受范围内。仿真计算的结果要大于理论计算结果,且随着储水量的增多,两者的偏差越来越小。造成偏差主要有两点:(1)PCS水箱的形状不规则,文章中采用的理论公式是在环形水箱的基础上进行的近似;(2)理论公式中没有考虑流-固耦合之间的作用,同时有许多简化,因此会与数值仿真结果存在一定偏差。

    • PCS储水箱位于屏蔽构筑物顶端,距地面高度约为90米,因此需要用对应的楼层反应谱对其进行抗震分析。阻尼比为0.5%,90米处的竖直方向的楼层反应谱,如表5所示:

      表 5  楼层反应谱

      Table 5.  Floor Response Spectrum

      频率 / Hz 0.100 0.135 0.138 0.174 0.522 1.035 100
      谱值 / g 0.023 0.043 0.045 0.087 0.299 0.573 1.024

      分别对PCS空水箱模型及60%储水量的水箱模型进行反应谱分析,结果如图8所示。

      图  8  空箱模型和60%储水模型水体位移图

      Figure 8.  Water Displacement of Empty Tank and 60% Capacity Tank

      空箱模型的最大位移发生在箱体顶部,最大位移为0.58 mm;而含60%储水量的水箱模型结果则有较大不同。首先反应谱分析下,水体的晃动较为剧烈,最大晃动高度达到2.82 m,而水箱发生的位移量则较小仅为0.53 mm,与空箱模型接近。这是因为水体的固有频率小,从而很容易被激起振动,而水箱结构的固有频率高,相对稳定。综上,在地震作用下,水体晃动非常明显,在水箱设计中,必须要将其影响考虑进去,并在实际工程中对水箱水位进行严格监控,以防止水箱失效。

    • 本文利用有限元软件ANSYS对核电厂屏蔽构筑物上方的PCS储水箱进行了模态分析及抗震反应谱分析,研究了不同储水量对水箱整体振动特性的影响。计算表明:水体由于固有频率较低,在地震载荷作用下很容易被激起,并产生较剧烈的晃动及晃动高度。因此在PCS水箱的设计中,必须将水体的晃动因素加以考虑。

  • 参考文献 (8)

    目录

      /

      返回文章
      返回