ABAQUS对高温管道裂纹尖端应力场的有限元分析模拟.pdfVIP

Abaqus 对高温管道裂纹尖端应力场的有限元分析模拟 王阳杰 （华东理工大学机械与动力工程学院 上海） 摘 要：蠕变是高温管道的主要破坏形式，为了对实际含缺陷的高温管道实际寿命的 预测，定量化研究蠕变的应力重分布时间及裂尖应力场是十分必要的。本文使用Abaqus 有限元软件对某超超临界再热段蒸汽管道进行了模拟。模拟结果与实验室试样结果趋势一 致。 关键词：蠕变；Abaqus；管道；应力重分布 1 引言： 蠕变是指在载荷和温度不变的条件下，材料的变形随时间的增加而缓慢增大的现 象，其过程经历裂纹的萌生、扩展，直至构件发生断裂[1]。随着火力发电技术的发展， 蒸汽管道也承受着越来越高的内压与温度，对其进行高温结构完整性评定便尤为重要， 搞清蠕变裂纹尖端的应力-应变场是这一切的基础。对建立准确的寿命预测有着重要的 实际意义。 本文使用 Abaqus 有限元软件对含缺陷的直管进行了应力应变场进行了详细的有 限元模拟，模型中使用标准椭圆模拟了管道的轴向表面裂纹，通过对不同时间下的蠕变 应力应变场的分析能够对含缺陷的高温管道应力应变场随时间的分布有了更清楚的了 解，以此更好地对高温结构进行完整性评定。 2 有限元模型 2.1 材料 在有限元模拟使用弹性-幂率蠕变材料模型[2]，蠕变应变率由下式给出： n n ε =A σ=ε ( σ/ σ) c 0 0 o （1） 材料选用2.25Cr1Mo钢，其在565 C 下的弹性和蠕变参数为[3]：弹性模量 E=154GPa，泊松比υ 0.3 ，蠕变常数 A=1.733E-26MPa-nh-1，应力指数n=9.732，标称应力取  -1 σ=330MPa ，相应的蠕变应变率 =0.0531h 。 0 ε0 2.2 直管模型几何尺寸 高温构件选用文献中的管道尺寸进行数值模拟[4] 。管道外半径 R=418mm，管壁 厚t=38mm，裂纹尺寸a/t=0.2、a/c=0.4 。在内壁和裂纹面上同时施加内压P=10Mpa。 为了避免边界效应的影响，管道长L=30t=1140mm[4] 。几何形状如图1 所示。裂纹简 化处理为椭圆形。 图1 裂纹参数表征 为了减少有限元计算的规模，根据管道构件的对称性建立四分之一的有限元 模型，在轴向和径向对称面上施加相应的对称约束，在远离裂纹的一端使用约束径 向两个方向上的位移防止有限元计算中的刚性位移，在远离裂纹一段的径向面上根 据薄壁圆筒公式施加与内压相应的轴向力。模型使用三维实体8 节点减缩积分单元 （C3D8R ），对模型使用适当的切分，如图 2 （a ）、（b）、（c ）所示。裂尖附近区域 使用细密网格。为了确定合适的网格尺寸，用30μm、40μm、50μm、60μm、70μm、 80μm 六种网格尺寸进行了网格敏感性测试，结果表明当网格密度小于40μm 时裂 尖正应力及裂纹扩展驱动力 C*就不再随网格密度的变化而变化了，管道的裂尖网 格密度选用40μm。 (a) (b) （c ） 图2 直管轴向裂纹有限元模型（a）整体模型（b）椭圆形裂纹网格（c）裂尖网格 3 结果及分析 图3，蠕变92h 裂尖正应力S 22