汽包应力分布的实时有限元分析.pdfVIP

V0119s 第19卷增刊 。豫之．怎口黧m uppl．199911 POWERENGINEERING 汽包应力分布的实时有限元分析 李斌 张聪师 段新会 071 (华北电力大学仿控所 003河北保定) 【摘要】本文利用当量管接头模型，采用一维有限元方法编制了汽包应力分布的实时肓限无监测程 序，并在电站仿真机中进行了应用。针对电站仿真机与实际参比机组任系统及显示参数一教的特点， 利用仿真机的运行参数对汽包狂热态启动过程中应力变化随时间的变化和应力分布进行r分析。 关键词汽包应力育限元 O前言 众所周知，汽包是锅炉最主要的厚肇承压部件。在机组启动、停机平¨负荷变化过程中， 汽包经常承受由于剧烈的温度变化引起的热应力及压力变化造成的机械应力的作用，这必然 加速汽包的疲劳寿命损耗。近年来汽包的低周疲劳寿命分析越来越受到人们的重视，运行人 员为了保证机组的安全运行迫切的需要一种监测汽包内应力变化的手段。为此，本文开发r 计算锅炉汽包应力的有限元实时监测程序，井在电站仿真机做了监测汽包应力实时分布的实 验。 1 计算模型的建立 汽包的最易损坏的部位是汽包F降管接头处，即通常所指的汽包纵向截面内转角和横向 截面外转角处，如果这两处不出现问题，整个汽包的寿命就可以得到保证。汽包下降管接头 是由两个直径不同的圆柱体组成的是三维结构．严格的讲，应该采用二维有限元方法进行分 析，但这种方法计算L作量很大，难以实现实时监测汽包应力的目的，因此要实现汽包应力 的实时监测必须对这一问题进行简化处理。 本文采瑚了当量管接头模型，将汽包r降管接头转化为球壳与圆+t体的组合结构，利用 二维轴对称有限元方法进行分析。当董管接头模l_是指在计算汽包温度分布时按汽包实际的 结构尺寸计算，当计算应力分布时将汽包壁的内径放大为原尺寸的3倍作为球壳的内径，其 他尺寸(包括汽包晕的厚度、下降管内径及r降管壁的厚度等)不变。这样计算山的』、i力分布 与采埘二维有限元方法的计算结果近似，可以满足{‘样计算的要求。这种模犁铂。保证计算精 度的前提r，减小了计算的I：作量，使程序可以满足实时监洲汽包麻力分布的要求， 计雉模刑的结构如图l，圈2所示。j刳中％为，e包内径．‘为当量管接头的内释。汽包 的儿何』t寸址表l所示。 35】 表1汽包的几何尺寸 单位：(m) l 汽包简体外直径 J 汽包简体内直径 l 下降管外直径 l 下降管内直径 l J 2149 J 1743 f 720 l 420 l 例l温度场计算有限元}呵格剖分圈 圈2．当量管接头有限元网格剖分图 2温度场计算 计算温度场时截取图l所示的汽包的横截面进行分析，同时考虑了汽包周向温度和径向 温度的分布不均。由于汽包可以沿横向自由变形，所以轴向温度的不均匀引起的热应力可以 略去不计。其导热方程为 n， zc窘+争一卢等=。 式中，=71(x，J，，r)为某一时刻r截面内的温度分布：五为金属的导热系数；P为金属 密度：c为金属的比热。 由f结构的对称性，边界sl为绝热边界条件，即g力=O；边界S2为第三类边界条件： 由于下降管接头取足够长，边界S3视为绝热边界条件处理，边界s4按取汽包保温层内壁温 度的第四类边界条件处理。初始条件给出初始时刻的温度分布。这样方程(1)的定解条件为 飘=。 ‘