反响堆冷却系统主管道研究.docx

反响堆冷却系统主管道研究（核动力工程杂志）2014年第三期1分析准则疲惫和应力的验收准则主要采用ASME规范NB-360098版，其应力和疲惫的验收准则见表1。表1中各式形式如下：假定管道部件作为承压边界，其主要失效特征为环向应力引起的塑性垮塌，即由于缺陷存在而产生一个轴向裂纹扩展趋势。为了控制环向应力作为静载荷引起管道失效，应当首先校核环向应力以保证知足相当的管道壁厚要求；知足要求后再分别校核一次应力强度要求，一次加二次应力强度范围要求若不知足还应当采用简化的弹塑性不连续性分析，通过计算峰值应力强度范围进而评价疲惫强度，以及关于热棘轮效应的安宁性评价。2管道热瞬态分析用于疲惫分析的管道瞬态主要包含RCS的压力、管道内流质流动的速度以及瞬态的温度。这些都是基于时间递增可拟合得到的函数曲线。2.1主管道内的对流换热对于管内单向流动的情况，其对流换热系数能够用Seider-Tate关系式计算适用范围是0.7Pr120；Re10000；L/D604，如式8所述。在各瞬态列表中看到的温度只是冷却剂的实时温度，该温度将通过饱和单相流热交换的方式传递到主管道内壁。通过对主管道内外表对流换热经过的研究能够看到，当确定了管道尺寸和流质特性等方面因素后，管道内对流换热的系数将成为一个与流体流速v相关的函数，主要通过表征管道内流体流通的状态，即Re来影响对流换热系数的取值。当流速上升，管道内流体更多地体现为层流特征，同理当流速下降时Re过小则更多地表现出自然对流产生的内外表湍流，在研究对流换热的经过中需要区别对待各种瞬态，然而大多数RCS瞬态将适用于Seider-Tate关系式。2.2热瞬态分析时间积分与时间步长为了保证时间积分与时间步长t的合理性，一般在计算经过中使用2种无量纲的因子Biot数Bi和Fourier数Fo来确定时间步长。Bi是不考虑尺寸的热阻对流和传导比例因子。2.3沿管道壁厚方向温度场的求解在有限元求解的时候，为了求得沿壁厚的温度分布，将模型建成一个同心圆平面，在内圆面上施加流体的对流换热系数和流体温度。考虑到主管道的外部包着隔热材料，因而保守以为管道外部为绝热，且初始温度设定为设计运行温度。即把整个模型设置成为温度内低外高的散热经过。之前已经介绍过，有限元分析后处理的目的在于提取1T和2T供本文第2节所述的公式计算使用。通过计算，在每一个瞬态时刻的温度分布范围均能够用一个函数Ty来表示图1，温度分布范围能够分成3部分：在疲惫计算中需要提取的是1T令1T=V和2T。这在后处理中能够通过管道壁厚方向定义途径，提取到离散节点的y值和对应的温度，进而进行积分求解。1T的计算得到的值也适用于ASMENB-3653.7热棘轮效应的影响校核，许用线性部分的温差1T应当不超过本文式7中规定的数值。同时应当注意到，根据NB-3611中的要求，在NB-3200中允许采用1.5Sm来替代公式7中的Sy值来进行热棘轮的安宁性评价。考虑到NB3600的应力计算结果能够被NB3200所包络，而通常NB3600的公式法结果更为保守，本文在最小壁厚的迭代经过中采用的是NB3600的有关阐述。3计算流程主管道疲惫分析的核心内容为沿壁厚方向的温差分布，即热瞬态分析。根据ASMEBPVCode第三卷NB分卷管道设计部分的内容，求解最小壁厚的流程如图2所示。以此为根据，采用循环累加的方法得到最小壁厚，该值可为以后可能出现的主管道壁厚减薄的情况提供审评根据。在求解主管道最小壁厚的经过中采用了优化算法，即设定一个目的函数以及多个等式或不等式边界条件，让目的函数趋近边界以得到最优解。流程展现了整个求解的思路和算法实现经过中的逻辑经过。详细的计算经过能够通太多种计算机语言实现。本文采用FORTRAN95语言编程实现标准规范公式法应力分析计算、疲惫雨流法计数和瞬态叠加的经过。在温度瞬态的计算经过中，则采用ANSYS表参数的方式对瞬态经过进行输入输入参数包括各瞬态下的冷却剂温度、RCS的压力、流速，进行求解和后处理，将温度结果导入FORTRAN95编写的程序进行处理。固然在核电的设计经过中不主张采取最优解的方法，但是本算例中的边界条件完全根据成熟且保守的标准规范，这些计算方法、参数选取、材料等均保守可接受。此外边界条件不牵涉制造、安装和检验经过中的技术条件，仅专注考虑在设计强度指标中的要求。因而在材料和制造等环节质量能够保证的前提下这种求解最值的方法是可靠的。4结束语本文从标准规范的要求出发，对管道分析的应力、疲惫和热棘轮情况进行了总结，采用有限元热瞬态分析代替有限差分方程来求解标准规范计算中需要的沿壁厚线性温度分布和非线性温度分布情况，同时采用优化方法进行迭代循环以不断趋近于能够知足ASME规范要求的最小壁厚值。本项研究可为管道设计、制造和安装经过中出现的偏差或不符合项提供参考，为安全审评提供根据。