蓝 翔,徐 鸿,李梦阳,倪永中

(华北电力大学电站设备状态监测与控制教育部重点试验室,北京 １０２２０６)

摘 要:采用内压和拉伸组合加载方法,在６５０ ℃下对 P９２钢管状试样进行了不同多轴度的蠕变试验并利用改进的KachanovＧRobatnov蠕变模型对其蠕变行为进行有限元模拟,观察了断口形貌和显微组织,分析了多轴度对 P９２钢应力分布和损伤发展的影响.结果表明:在内压和拉伸组合加载下,当外壁等效应力相同时,多轴度越大,P９２钢的蠕变寿命越短;多轴度对 P９２钢蠕变孔洞的生长具有促进作用;改进的 KachanovＧRobatnov蠕变模型可以准确地描述P９２钢的蠕变行为;在蠕变过程中,多轴度影响应力的分布,进而影响损伤的分布,多轴度大的位置其损伤程度也大.

关键词:P９２钢;多轴度;蠕变行为;有限元模拟

中图分类号:TG１１３．２５ 文献标志码:A 文章编号:１０００Ｇ３７３８(２０１８)０３Ｇ００７４Ｇ０６

０ 引 言

我国是世界上拥有超(超)临界机组最多的国家,随着蒸汽温度和压力参数的提高,在机组效率得到提高的同时,也对电站关键部件材料的性能提出了更高要求,尤其是材料的高温强度、抗高温腐蚀和氧化能力等.P９２铁素体耐热钢由于具有良好的高温导热性能、高强度和低热膨胀系数等,被广泛应用在我国 ６００/６２０ ℃ 等 级 超 超 临 界 机 组 高 温 部 件

中[１],如高温主蒸汽管道、高温再热蒸汽管道及联箱等.这些部件长期处于高温工作状态,其主要失效形式为蠕变失效,同时在复杂应力和部件几何尺寸的影响下,这些部件大多处于多轴应力状态.通过研究多轴应力状态下 P９２钢蠕变损伤规律和机理,可以准确评估部件的剩余寿命,对超超临界机组高温关键部件的寿命管理具有重要意义.国内外学者对多轴应力状态下蠕变的研究主要集中在多轴蠕变试验方法[２],以及基于连续损伤理论并通过有限元方法模拟材料的力学特性和蠕变损伤特征等方面.一些学者采用缺口试样进行多轴蠕变试验,并研究了其蠕变损伤演变[３Ｇ６],这些研究的

试验装置和试样均比较简单,易于操作,但无法实现主应力和静水应力的独立变化,而研究表明利用管状试样可通过不同方式的组合加载来弥补该不足;WANG[７]等利用扭转Ｇ压缩Ｇ拉伸的加载方式,研究了多轴应力状态下 GH３３镍基合金的蠕变疲劳交互作用;SELIGER[８]等在研究多轴度对电站管道所用１４MoV６Ｇ３、１０CrMo９Ｇ１０和 X１CrMoVNb９Ｇ１耐热钢蠕变行 为 的 影 响 时,采 用 了 光 滑 薄 壁 管 和 环向缺口薄壁 管 两 种 试 样,通 过 控 制 轴 向 拉 伸 应 力

和管道内压来调整多轴度;YAO 等[９]对近年来国内外有关蠕 变 行 为 的 研 究 进 展 进 行 了 综 述,通 过引入多轴度 和 参 考 应 力 来 描 述 多 轴 蠕 变,根 据 蠕变过程是由约束孔洞控制还是塑性控制来确定参考应力 系 数;HYDE[１０]等 使 用 KachanovＧRobatnov蠕变模型研 究 光 滑 试 样 和 缺 口 试 样 的 蠕 变 行 为,通过引入不 同 的 等 效 应 力,获 得 与 试 验 结 果 一 致的模拟结果,但 未 考 虑 蠕 变 第 一 阶 段 特 性 和 应 力的多轴性.为研究多轴度对 P９２钢蠕变行为的影响,作者

采用内压和拉伸组合加载的蠕变试验方法,通过控制内压和拉伸载荷获得不同的多轴度,然后在６５０℃下对 P９２钢管状试样进行了３组不同多轴度的蠕变试验,观察了断口形貌和显微组织,分析多轴应力对P９２钢蠕变行为的影响机制;利用能够描述蠕变第一阶段特性的改进的 KachanovＧRobatnov蠕变模型对内压和拉伸组合加载下的 P９２钢蠕变行为进

行有限元模拟,并进行试验验证,分析了多轴度对P９２钢应力分布和损伤发展的影响.

１ 试样制备与试验方法

试验原料为日本JFE 钢铁公司生产的 ASTMA３３５GRADEP９２钢,热处理工艺为１０５０ ℃保温２０min正火,７８５ ℃保温６０min回火,化学成分如表１所示.

多轴蠕变试样的形状和尺寸如图１所示:该试样由管上端及封头、中部、下端及封头、填充芯棒等４部分组成;上端与中部、下端与中部均通过焊接方式连接;下端有与内压加载系统连接的通气孔以及与中部连接的导管,试样上下端通过机械螺纹与蠕变试验机连接.

多轴蠕变试验系统包括压力加载系统和拉力加载系统两个部分,如图２所示.压力加载系统通过增压泵将氮气加压到指定压力后注入试样内部,并保持试验过程中压力波动范围不超过±１％.拉力加载系统由力学加载部分和加热炉部分组成,试验过程中拉力波动范围不超过±１％,试样中部的温度变化范围不超过±３ ℃.试验过程中采用引伸计测得试样轴向方向的变形情况.材料的应力状态可以用多轴度来衡量,多轴度是一个无量纲量,一般定义为静水应力与等效应力之比[１１].为研究多轴应力对 P９２钢蠕变行为的影响,作者选取如表２所示的３组载荷组合,在６５０℃下进行多轴蠕变试验.

试验结束后测量试样断口的尺寸.在试样断口上取厚度约１０mm 的圆环,切割成１/４圆环试样,对其轴向和环向截面进行粗磨、细磨、抛光,用体积分 数 ４％ 的 硝 酸 酒 精 溶 液 腐 蚀 后,利 用 JSMＧ６４９０LV 型扫描电镜(SEM)观察其断口形貌和显微组织.

２ 试验结果

多轴蠕变试验结束后,不同载荷组合下的试样均在标距间有效管段中间偏上或偏下位置断裂.由图３和表３可知,蠕变试验后试样断口处的外径减小,内径增大,管壁发生颈缩导致壁厚减薄,最终发生断裂.由表３还可以看出:在相同等效应力条件

下,内压越大试样管壁的颈缩量越大,在相同内压条件下,等效应力越大管壁的颈缩量越大.由图４可知:试样３的外壁等效应力最大,最先发生断裂;试样１和试样２的外壁等效应力相同,试样２先于试样１断裂,这说明多轴度的增加会加速P９２钢的断裂.

由图５可知:试样１断口中韧窝的尺寸最大、数 量最少;试样３断口中韧窝的数量最多、尺寸最小; 试样２断口中韧窝的数量和尺寸介于试样１和试样 ３之间.

由图６可以看出:在试样１的轴向与环向截面中孔洞的数量和尺寸相当,孔洞基本分布于晶界和马氏体板条边界上;试样２轴向截面中孔洞的数量比环向截面中的多,孔洞尺寸也更大,这说明增加内压后,垂直于轴向的应力增大,导致孔洞沿环向方向长大,使具有较大多轴度的试样２比试样１更早发生断裂.综上可知,多轴度对蠕变孔洞的生长具有促进作用,从而影响 P９２钢的蠕变寿命.

３ 有限元分析及讨论

KachanovＧRobatnov(KＧR)蠕变模型[１２]常用于描述蠕变第二、第三阶段特性,但未考虑蠕变第一阶段特性和应力的多轴性.在 KＧR 模型的基础上,加入描述蠕变第一阶段特性的方程,并考虑应力的多轴性及局部损伤效应[１３],改进后的 KＧR 蠕变模型为

式中:εcij 为应变张量;Sij 为偏应力张量;σ１ 为第一主应力;σe 为等效应力;D 为损伤变量;α 为表征多轴断裂准则的材料常数;A,B,m,n０,n,M 和ν均为材料常数;φ 为损伤常数;ρ 为考虑损伤非均匀性的材料常数;t为蠕变时间.在６５０ ℃、内压和拉伸组合加载下 P９２钢的蠕变模型参数如表４所示.考虑多轴蠕变试样的对称性,选取其１/４结构作为计算模型,采用 SOLID１８６平面单元进行网格划分,有限元模型及加载方向如图７所示.

将改进的 KＧR 蠕变模型嵌入到有限元分析软件 ANSYS中,根据表４中的模型参数计算在内压和拉伸组合加载下试样的蠕变曲线,并与试验结果进行对比.由图８可知,蠕变曲线的模拟结果与试验结果基本吻合.因此,该模型可用于模拟在内压

和拉伸组合加载下 P９２钢的蠕变行为.从有限元模拟结果中提取试样外壁等效应力、第一主应力和多轴度随时间的变化曲线,并将时间做归一化处理,如图９所示,图中tr 为断裂时间.由图９可知:所有试样外壁等效应力和第一主应力均在加载初始时刻瞬间增加到一定值,之后随着时间的延长而逐渐减小;试样２的外壁多轴度始终大于试样１和试样３的,试样３的外壁多轴度介于试样１和试样２的之间.

在轴向拉伸应力１００ MPa、内压１５ MPa条件下 ,不同时刻的多轴度、等效应力、第一主应力和静水应力沿壁厚的变化曲线如图１０所示,图中 Ri 为内径,Ro 为外径,r 为管壁不同厚度处的半径.由图１０可知:外壁多轴度大于内壁的,在蠕变

第一阶段和第二阶段,多轴度基本不随试样壁厚的变化而变化,在蠕变第三阶段,损伤累积导致应力重新分布,多轴度也随之变化;在蠕变初始阶段,等效应力、第一主应力和静水应力沿壁厚分布不均匀,外壁的第一主应力和静水压力均大于内壁的,外壁的等效应力小于内壁的;在蠕变过程中,内壁的这３种应力均先减小,之后维持在一个稳定范围,并且在蠕变断裂时刻突然增大;外壁的这３种应力均先增大,之 后维持在一个稳定范围,并且在蠕变断裂时刻突然减小,这说明此时试样外壁已经变形失效.在蠕变过程中,试样在半径约为１１．２５mm 处的多轴度、等效应力、第一主应力和静水应力基本不随时间而变化,即为该试验条件下试样的骨点[１４Ｇ１５].

图１１为试样３有效管段损伤随时间的分布,其中左侧为试样内壁,右侧为外壁.由图１１可知:试样有效管段的损伤沿壁厚分布不均匀,由于外壁的多轴度大于内壁的,使得外壁损伤的增长速率大于内壁处的.随着损伤的进行,外壁产生裂纹,试样的承载能力下降,应力重新分布,内壁的等效应力和第一应力增加,外壁的减小,最终导致试样断裂.综上所述,多轴度影响应力的分布,进而影响损伤的分布,多轴度大的位置其损伤程度也大,从而导致试样在该位置失效.

４ 结 论

(１)在内压和拉伸组合加载下,当外壁等效应力相同时,多轴度越大,P９２钢的蠕变寿命越短;多轴度对 P９２钢蠕变孔洞的生长具有促进作用.

(２)采用改进的 KＧR 蠕变模型对６５０ ℃、内压和拉伸组合加载下 P９２钢的蠕变行为进行有限元模拟,模拟结果与试验结果基本吻合,该模型可以准确描述 P９２钢的蠕变行为.

(３)在蠕变过程中,多轴度影响应力的分布,进而影响损 伤 的 分 布,多 轴 度 大 的 位 置 其 损 伤 程度也大.

来源：材料与测试网