反应堆压力容器及堆内构件整体大规模三维有限元地震分析.docxVIP

第32卷第2期核动力工程Vol.32.No.22011年4月NuclearPowerEngineeringApr.2011

文章编号：0258-0926(2011)02-0044-05

反应堆压力容器及堆内构件整体大规模三维有限元地震分析

姜乃斌1，臧峰刚1，张利民2，张传勇2

（1.中国核动力研究设计院核反应堆系统设计技术重点实验室，成都，6100412.日本泰科诺斯达株式会社北京代表处，北京，100098）

摘要：运用适用于大规模有限元模型求解的TechnoStarVENUS程序，对反应堆压力容器、堆内构件和燃料组件进行了整体大规模的三维有限元地震分析。分析中分别考虑了阻尼效应和非线性接触，在有限元模型自由度数超过2×107的情况下，分别给出了这2种情况在运行基准地震（OBE）激励下的时程响应。实现了现有计算机硬件条件下的整体大规模三维有限元地震分析，为今后核动力工程中的复杂设备或系统的整体动力分析探索了一条新的途径。

关键词：反应堆结构力学；地震分析；有限元；大规模中图分类号：O327 文献标识码：A

1引 言反应堆地震分析是核能领域结构力学分析工

作中重要内容。由于反应堆压力容器及堆内构件结构复杂、非线性因素众多，三维实体有限元分析十分困难，国内外通常采用二维简化方法对反应堆及堆内构件进行地震分析。虽然使用二维简化方法计算速度快，对计算机软硬件要求低，但由于模型的建立需要大量的工程经验和试验参数，不利于推广。本文试图探索一种反应堆压力容器及堆内构件的整体三维有限元地震分析方法，为今后核动力工程中的复杂设备或系统的整

体动力分析提供技术储备。

起主要支承作用的堆内构件和燃料组件。起主要支承作用的堆内构件包括：上部支承板构件、上支承柱、吊篮、上堆芯板、下堆芯板、堆芯支承柱和堆芯支承板。对于堆内测量系统、控制棒驱动机构、控制棒组件以及反应堆内冷却剂以及其他堆内构件，只通过折算密度的方式考虑质量的影响，而忽略其刚度的影响。按照燃料组件的外围尺寸，将燃料组件简化为213.5mm×213.5mm×4057.0mm的均匀方柱。

首先建立分析部件的三维实体模型，再进行有限元网格划分。模型离散采用4节点4面体单元，整个模型单元数超过4.38×106，节点

数大于6.98×106，总自由度数高达2×107。各

2有限元分析 个部件的单元和节点数量如表1所示，网格图见2.1分析软件 图1。

为实现反应堆压力容器及堆内构件整体大规模的三维有限元地震分析，采用适用于大规模有限元模型求解的新一代分析软件——TechnoStarVENUS（TSV）程序。该程序的求解器采用CoarseGridBasedConjugateGradient（CGCG）算法，可

以大幅度增强和提高分析求解的能力和速度。

模型中定义上支承板裙筒外表面与吊篮筒体内表面之间的接触以及吊篮筒体出口与压力容器出口端面之间的接触为非线性接触关系。定义堆芯支承板与压力容器筒体连接节点的水平自由度耦合，以模拟径向支承键的作用；耦合上堆芯板

与吊篮筒体连接节点的水平节点自由度，但竖直

2.2模型简化与网格划分 方向自由；模型中其他各部件接触位置处均定义建立的三维实体模型包括反应堆压力容器、 为共节点的固定连接。

收稿日期：2010-01-11；修回日期：2010-08-02

姜乃斌等：反应堆压力容器及堆内构件整体大规模三维有限元地震分析 45

表1各部件的单元和节点数量 表2反应堆各部件的材料属性

Table1ElementNumberandNodeNumberofEachComponent

Table2MaterialPropertiesofEachReactorComponent

部件 弹性模量/MPa 泊松比密度/kg·m3

部件

上部支承板构件吊篮

上支承柱

上堆芯板

单元数 节点数121263 198923340061 545263470125 79