P91钢高温管道蠕变裂纹扩展寿命评价.pdf

2023.04 设备监理60检验与技术I n s p e c t i o n a n d T e c h n o l o g yP91 钢高温管道蠕变裂纹扩展寿命评价谢知伟（福建省锅炉压力容器检验研究院 福州 350000）摘 要：含有缺陷的设备在高温、高压条件下容易发生蠕变损伤破坏，进而需要准确评价高温结构的蠕变裂纹扩展寿命。以某电厂主蒸汽管道用 P91 钢为研究对象，结合理论分析、实验及数值模拟的方法对其展开研究，采用 ABAQUS 软件，建立管道半椭圆形裂纹有限元模型，研究了 2 种压力工况下不同初始裂纹尺寸的蠕变裂纹扩展寿命，并与理论计算结果进行了对比，结果表明有限元计算结果更加符合实际要求。关键词：蠕变损伤 有限元模型 裂纹扩展寿命Life Evaluation of Creep Crack Growth in P91 Steel High Temperature PipelineXie Zhiwei(Fujian Boiler and Pressure Vessel Inspection and Research Institute Fuzhou 350000)Abstract There are some initial defects which make the equipment prone to creep damage under high temperature and high pressure.Therefore,it is necessary to accurately evaluate the creep crack growth life of high temperature structure.The P91 steel used for main steam pipeline of a power plant is taken as the research object,and the theoretical analysis,experiment and numerical simulation methods are used to study the P91 steel,the finite element model of the semi elliptical crack in the pipeline is established by using ABAQUS software,and the creep crack growth life of different initial crack size under two kinds of pressure conditions is studied,and compared with the theoretical calculation results.The results show that the finite element results are more reasonable.Keywords Creep damage The finite element model Crack propagation life中图分类号：TB497 文献标志码：B文章编号：2095-2465(2023)08-0060-05 DOI:10.19919/j.issn.2095-2465.2023.08.015作者简介：谢知伟（1987 )，男，硕士，质量工程师，从事承压类特种设备检验检测工作。通讯作者：谢知伟，E-mail：。（收稿日期：2023-04-04）随着科学技术的不断发展，环保要求的进一步提高，为了提高设备的效率，高温设备朝着高压力、高温度的方式逐步发展，然而高参数的发展不可避免地对设备的安全提出了更高的要求1-3。由于制造工艺水平的限制使得工程中的高温管道、容器在制造过程中存在一些缺陷，因而在苛刻环境下长时间运行的设备经常发生蠕变损伤和断裂4-5，为了更好地评价高温结构寿命则需要准确评价高温结构的蠕变裂纹扩展寿命，而能准确描述材料蠕变变形的本构模型是高温结构寿命评价的保障。P91 材料可承受 600 的高温条件，并能在高温条件下保持良好的性能，因此被广泛应用在高温场合，如主蒸汽管道、高温再热管道。然而材料在高温下的性能相比常温下有很大差别。材料在高温条件下的蠕变曲线如图 1 所示，图 1 蠕变曲线?IIIIII?PLANT ENGINEERING CONSULTANTS 2023.0461检验与技术I n s p e c t i o n a n d T e c h n o l o g y从图 1 中可知，该曲线由 3 部分组成，分别为蠕变第一阶段、蠕变第二阶段和蠕变第三阶段，其中第二阶段也称为稳态蠕变，其对应的蠕变速率称为稳态蠕变速率或者最小蠕变速率。此外材料的蠕变曲线与应力水平和温度相关6-7，提高温度或者应力中的其中一个因素，蠕变应变速率都会提高，因此研究不同温度和应力条件下材料的蠕变性能很有必要。1 实验方案1.1 试样加工实验所采用的材料为某电厂高温蒸汽管道材料P91，其成分见表 1。材料CMnSiNCrMoNbCuV母材0.110.380.0230.0468.810.970.130.0760.21焊缝0.0881.020.120.048.731.020.180.030.22表 1 P91 材料成分%试 样 的 加 工 根 据 GB/T 20392012 金 属 材料 单轴拉伸蠕变试验方法8进行，试样示意图如图 2 所示，本次所有试样采用圆棒试样，其标距为 100 mm，两端螺纹与夹具连接起来进行实验加载。1.2 实验过程 1.2.1 试验机参数本次实验采用的实验设备为 CTM-105B1 型电子高温蠕变持久试验机，见图 3。图 3 电子高温蠕变持久试验机该试验机的主要参数如下：1）试验力为 50 kN；2）试验力准确度 0.5%；3）试验力分辨力为试验力的 1/200 000(全程分辨力不变）；4）计时误差为0.2%；5）主机功率为 2 kW；6）主机电源为 380 V；7）电炉结构采用三段电炉丝加热，分段控制；8）工作温度范围为 200 1 000；9）有效均温区长度为 150 mm；10）电炉升降方式为电动升降机构提升和下降电炉。1.2.2 实验流程1）打开实验机电源，点击控制面板上的炉子上升、下降按钮，控制炉子的位置；2）通过调节上下拉杆的位置，将加工好的圆棒试样安装在上下拉杆上，然后在试样上、中、下 3 个位置绑上热电偶；3）给试样初始加一个小载荷，稳定好试样，然后固定位移引伸计，并清零位移记录点；4）启动电脑，通过电脑控制界面进行实验参数的设置，包括温度、载荷以及时间；5）启动实验，观察电脑面板及实验操作面板各参数的变化情况；6）在实验操作过程中时刻关注实验进程，如有特殊情况，可按紧急停止按钮；7）实验结束后，点击控制面板上的停止实验，然后对实验数据进行保存，并退出实验系统；8）手动控制实验机面板来调节上下拉杆的位置进行卸载并关掉加温装置，让实验温度降下来；9）当炉温降到室温时关闭电源，手动拆掉试样和引伸计。1.3 实验结果与分析本次实验所采用的材料为 P91 材料，实验温度为图 2 试样示意图1851003A13.525245R0.50.02A0.031.6?mmGB 145?850.8117R12IIM16?6h15132023.04 设备监理62检验与技术I n s p e c t i o n a n d T e c h n o l o g y600，试验应力水平根据材料的屈服应力的 80%左右来进行估算，该材料在 600 下的屈服强度在 240 MPa左右，因此本次蠕变实验的应力水平设定了5 组，分别为110 MPa、125 MPa、160 MPa、180 MPa、200 MPa。图 4 为该材料在 600 下，不同应力水平下的蠕变应变与时间的蠕变曲线。从图 4 中可以看出，随着应力水平的增大，材料的蠕变应变也是增大的。图 4 不同应力水平下蠕变拉伸曲线0.000 10.1110100?h1 00010 000100 0000.001?0.010.1110 MPa 125 MPa160 MPa180 MPa200 MPa11.4 蠕变本构方程为了准确评价高温结构蠕变裂纹扩展寿命，需要准备蠕变本构方程来表征材料的蠕变变形情况。经过几十年的发展，蠕变本构方程有了很大的改进，包括Norton 本构、Norton-Bailey 本构、K-R 本构9以及一些多项式本构方程。分别采用上述本构方程来拟合试验曲线获取本构方程参数，进而为后续结构的寿命评价提供参考。2 有限元分析首先构建含轴向表面裂纹管道的数值模型，用有限元法分析裂纹扩展形状、裂纹周边变形及应力分布。根据 GB/T 196242019在用含缺陷压力容器安全评定10，管道中的表面裂纹一般可以简化为半椭圆形裂纹，管道承受内压作用时，裂纹会沿着裂纹尖端 扩展，其中包括沿轴向扩展和沿厚度方向的径向扩展。2.1 材料及参数工作温度为 600，材料的蠕变性能及本构模型见式（1）：(1)?m?An式中：A与材料特性和温度有关的常数；n稳态蠕变速率的应力指数；应力。材料参数可以通过对上述方程两边取对数求得。由于式（1）中蠕变速率是应力的幂函数，所以符合式(1)的蠕变称为幂律蠕变。一般幂律蠕变材料的应力指数为 5，颗粒强化材料的应力指数相对高点。当然一些材料随着应力升高,log log?曲线的斜率开始增大且偏离线性关系，这种现象称为幂律失效。一般而言幂律失效情况下，材料蠕变性能进入一个高应力区，蠕变机制发生了改变。这种情况可以用其他的应力本构方程表示。通过对式(1)两边取对数可以得到式(2)：(2)lglglg?m?nA然后将不同应力水平条件下对应力的蠕变应变速率放到 origin 图中可以得到图 5。图 5 应力与最小蠕变速率关系10?810?710?610?510?410?310?2101102103?MPa?h通 过 数 据 拟 合 可 以 获 得 材 料 参 数 A=3.368 45E-28，n=10.464 7，R2=0.944 83。2.2 有限元建模本次有限元模型以半椭圆形裂纹进行建模，并进行有关分析，裂纹的主要尺寸为：a/c=0.3，0.6，1.0；裂纹深/厚比=a/t=0.1，0.2。其中 a 为裂纹沿径向的深度尺寸，c 为裂纹沿轴向方向的尺寸，t 为管道厚度。本次模拟了3种裂纹尺寸下的管道裂纹扩展情况，其中内压为 10 MPa，具体尺寸如图 6 所示，有限元模型如图 7 所示。裂纹初始尺寸为：t=38 mm,a0/t=0.1 (a0=4 mm),a0/c0=0.8(c0=5 mm)，由于具有对称性，因此在进行建模时利用对称性只需要建立管道的 1/8 模型即可，模型的边界条件为：在管道对称面上施加对称约束，管道所受的载荷工况为内压及轴向力。采用PLANT ENGINEERING CONSULTANTS 2023.0463检验与技术I n s p e c t i o n a n d T e c h n o l o g yABAQUS 软件进行有限元建模，对于管道模型本次选用的单元为八节点完全积分实体单元(C3D8)，模型包含 49 888 个单元。（a）半椭圆形裂纹管道示意图（b）半椭圆形裂纹管道尺寸图图 6 半椭圆形裂纹管道示意图2LRiRotPi?2catRil图 7 有限元模型示意图2.3 损伤模型所采用的损伤模型为基于延性耗竭的损伤模型，其表达式见式(3)：(3)?=cf*该模型被大多数学者用来模拟蠕变裂纹的扩展，其中损伤速率定义为蠕变应变率c与多轴蠕变延性f*的比值。在式(3)中损伤的累积可以通过对式(3)进行积分获得，见式(4)：(4)?=c00tttt?ddf*在式(4)中，损伤变量 来表征模型的损伤量，根据其定义，该值介于01之间。当损伤值达到1时，则认为该单元失去承载能力，进而裂纹进行扩展。2.4