控制棒下落时间计算模型.PDFVIP

第37 卷 第11 期 核 技 术 Vol.37, No.11 2014 年11 月 NUCLEAR TECHNIQUES November 2014 控制棒下落时间计算模型 刘 新 陈先龙 张 修 于晓雷 （深圳中广核工程设计有限公司 深圳 518172 ） 摘要 反应堆驱动线是反应堆控制的“生命线”，执行重要的安全功能。控制棒下落时间不仅是驱动线设计考 核的重要指标，也是核电站安全分析的关键参数。基于压水堆驱动线的结构特性，对控制棒组件在下落过程进 行受力分析，提出落棒时间计算的数学模型，采用有限差分法求解运动方程，并且在此基础上编制计算程序。 计算结果最终表明：程序计算值与岭澳核电站二期机组热态落棒试验值符合较好，基于此计算方法编制的落棒 时间计算程序可以用于压水堆停堆时的安全评价。 关键词 控制棒组件，驱动线，落棒时间 中图分类号 TL364 DOI: 10.11889/j.0253-3219.2014.hjs.37.110604 控制棒组件在反应堆中的落棒时间对于反应堆 的安全运行至关重要。反应堆驱动线由运动组件和 1 理论模型 固定组件两部分组成。运动组件包括控制棒棒束组 控制棒运动组件不仅受到重力还受到阿基米德 件(Rod Cluster Control Assembly, RCCA)和驱动杆， 浮力、水力压力、水力摩擦力等，其中有些力随时 固定组件是运动组件的下降通道，主要包括：控制 间进行变化，为了计算落棒时间，我们需要对计算 棒驱动机构(Control Rod Drive Mechanism, CRDM)、 模型进行简化假设，并且要把运动过程按照时间步 热套管、控制棒导向筒组件和燃料组件导向管。反 长进行离散，假设在单个时间步长过程中力是恒定 应堆驱动线各部件结构复杂，落棒时间计算需要考 的，然后对时间进行积分计算。 虑反应堆内冷却剂的流速、温度和压力的影响，也 要考虑外界因素（如地震）的影响。 1.1 简化假设 控制棒在下落过程中，由于通道流体流速的变 (1) 在控制棒导向筒内，作用在星型架和控制 化，导致受到的水力学力（水力摩擦力和水力压力） 棒的水力学力被忽略； 会不断的变化，而且变化值与下落位置是非线性关 系，从而增加了计算落棒时间的难度[1] 。窦一康等[2] (2) 导向管内的水温为常数且等于堆芯中冷却 剂的平均温度； 提出了地震工况下用有限差分法对振动方程作离 (3) 在计算的任何时刻，流体应为单向的。 散，对时间历程则用Wilson-θ法进行逐步分析，对 由碰撞引起的非线性采用修正的Newton 法作迭代 1.2 控制棒受力分析 处理的计算方法。但是，常微分方程的 Wilson-θ 控制棒下落过程受到的力包括恒力和变力（如 求解方法比较复杂，因此，我们采用有限差分方法， 图1 所示），其中，恒力有重力F 、浮力F 、额定 将下落时间进行离散，同时利用圆管中流体压差计