首页 / 专利库 / 物理 / 应力 / 拉伸应力 / 塑性条件下含残余应力的蠕变孕育期预测方法

塑性条件下含残余应的蠕变孕育期预测方法

阅读:1019发布:2020-07-17

专利汇可以提供塑性条件下含残余应的蠕变孕育期预测方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种塑性条件下含残余应 力 的蠕变孕育期预测方法,在Davies工作的 基础 上,提出了考虑残余 应力 的蠕变孕育期 预测模型 。利用参考应立法,引入了弹性追随因子Z计算考虑残余应力的蠕变孕育期。使用紧凑拉伸试样(CT)通过预压缩产生残余应力,并施加主 载荷 进行蠕变实验。本发明的有益效果是:本设计方法能够将原有的预测模型扩展到含残余应力的模型中,从而提出一种简化的塑性条件下蠕变孕育期预测方法,因此能够在结构中简洁有效的预测出塑性条件下蠕变孕育期。,下面是塑性条件下含残余应的蠕变孕育期预测方法专利的具体信息内容。

1.塑性条件下含残余应的蠕变孕育期预测方法,包括以下步骤:
S1:建立模型:所述模型包括CT试样本体,所述CT试样本体的中部前端设有槽,槽的后部设有缺口,CT试样本体上还设有上主载荷销孔、下主载荷销孔,上主载荷销孔、下主载荷销孔上下对应设置,分别设置在槽的上下两端;
S2:首先利用上圆销、下圆销对CT试样本体的上下两端进行预定大小的压缩加载,然后释放上圆销、下圆销,会在CT试样本体的缺口附近产生残余应力分布;
S3:在含有残余应力的缺口处插入预制裂纹,以进行蠕变试验;
S4:利用销子在上主载荷销孔、下主载荷销孔施加主载荷,进行高温蠕变试验;
S5:通过蠕变有限元模拟可以获得计算含残余应力CT试样孕育期所需要的必要参数,在塑性条件下,计算孕育期主要包括以下步骤:
(1)首先计算复合加载下的应力强度因子,其计算公式为:
(Ⅰ)中:
其中: 是模拟计算的只含有残余应力下的应力强度因子,单位为MPa·(m1/2); 是主载荷应力强度因子,单位为MPa·(m1/2);P是主载荷,单位为N;B是试样厚度,单位为mm,Bn是试样净厚度,单位为mm;a/W是预制裂纹长度比率,a是预制裂纹长度,采用上主载荷销孔圆心到预制裂纹后端的平直线距离,单位为mm;W是名义试样宽度,采用上主载荷销孔圆心到CT试样本体后端的水平直线距离,单位为mm;f(a/W)是CT试样几何系数,只与a/W有关;
V是无量纲的塑性相关项,计算如下:
(Ⅱ)中:V0是无量纲参量,
其中: 是塑性残余应力强度因子,单位为MPa·(m1/2); 是弹性残余应力强度因子,
1/2
单位为MPa·(m ), 利用JS计算,JS是残余应力场下断裂参量,单位为MPa·m;
其中:E′是有效弹性模量:E'=E/(1-ν2),E是弹性模量,ν是泊松比, 和JS都利用有限元模拟结果提取;
(Ⅱ)中:Lr是无量纲参量,描述主载荷幅度:
其中:σy是屈服强度,与材料属性有关,单位为MPa, 是主载荷参考应力,单位为MPa,用下式计算:
其中:nL为无量纲裂纹深宽比参数,通过下式计算:
常数
(Ⅱ)中:
是弹性主载荷应力强度因子, 是塑性主载荷应力强度因子,单位为MPa·(m1/2);
利用有限元模拟结果计算:
(Ⅱ)中:β描述残余应力的幅度,是无量纲参量;
是二次载荷参考应力,利用有限元模拟计算,
(Ⅱ)中:Z为无量纲的弹性追随因子,从有限元模拟结果中提取出应力应变关系,取等效蠕变应变增量 与等效弹性应变增量 的比值:
(2)计算塑性复合应力场下J积分数值,其计算公式为:
其中:KI是复合应力强度因子,E′是有效弹性模量:E'=E/(1-ν2);
(3)然后计算塑性应力场下孕育期时间tiHRR,其计算公式为:
(Ⅲ)中:εcrit是单轴蠕变韧性,与材料属性有关,单位为1, 是蠕变应变变化率,与材料高温蠕变属性有关,σP0是标准化应力,单位为MPa,εP0是标准化应变,单位为1,α为应变硬化系数,N为应变硬化指数,IN是与N有关的无量纲函数, 是与θ和N有关的无量纲函数;
(Ⅲ)中:d(mm)是判定蠕变萌生发生时裂尖前蠕变损伤达到1所延伸的距离,即蠕变萌生发生的临界距离;
(Ⅲ)中:MSFHRR为塑性条件下多轴应力因子,根据Cocks and Ashby关系式计算:
n为无量纲的蠕变应力硬化指数,sinh是双曲正弦函数,hHRR为塑性应力三轴度,在塑性应力状态下:
其中:θ是裂纹尖端度,N是无量纲的应力硬化指数, 和 是与θ和N有关
的无量纲函数。
2.如权利要求1所述的塑性条件下含残余应力的蠕变孕育期预测方法,其特征在于:d取所研究材料的晶粒尺寸。
3.如权利要求1所述的塑性条件下含残余应力的蠕变孕育期预测方法,其特征在于:Bn=B。
4.如权利要求1所述的塑性条件下含残余应力的蠕变孕育期预测方法,其特征在于:所述有限元模拟采用ABAQUS6.14进行计算模拟, JS、 Z的提取过程包括以下步骤:
(1)首先建立预压缩加载的CT试样的有限元模型,在材料属性模里设置弹性塑性参数,在载荷模块里设置压缩载荷,以及拘束条件:包括对称条件和固定条件,在接触模块里设置压缩圆销与试样上下表面的的刚性接触,在分析步模块里设置好输出参量:应力值,在网格模块划分网格;
(2)在作业模块提交任务计算,获得残余应力的计算结果,结果文件中,从场变量中可以直接提取二次载荷参考应力σrefS;
(3)建立相同尺寸的试样模型,进行主载荷拉伸试验,在材料属性模块里设置高温下的弹性塑性蠕变参数,在网格模块划分网格,在接触模块里设置拉伸销子与销孔的刚性接触,并在模型中插入预制裂纹,在分析步模块里设置好输出参量:应力应变值、应力强度因子K值、断裂参量J积分值,在载荷模块里设置拉伸载荷,以及拘束条件:包括对称条件和固定条件,在预加载应力场中导入上步算好的残余应力;
(4)在作业模块提交任务计算,获得含有残余应力的蠕变拉伸实验计算结果,结果文件中,在插入裂纹后还未施加拉伸载荷时刻,从历史变量中可以获取弹性残余应力强度因子以及残余应力断裂参量JS,在施加拉伸载荷的初始时刻,可以获取塑性主载荷强度因子从历史变量中可以获取等效应力随总应变增量的变化曲线,从曲线中得到等效蠕变应变增量 等效弹性应变增量 进而得到弹性追随因子Z。

说明书全文

塑性条件下含残余应的蠕变孕育期预测方法

技术领域

[0001] 本发明涉及含残余应力高温结构在塑性条件下的蠕变孕育期工程临界评价,就是确定当结构中存在表面裂纹且处在塑性应力条件下时,评价此高温结构的蠕变裂纹萌生寿命。

背景技术

[0002] 燃为主的能源结构是我国雾霾天气的主要成因之一,而燃煤发电是我国目前最主要的发电方式,该趋势将长期存在。因此,除了改变能源结构,发展高效洁净的超超临界(USC)机组是节能减排的重要途径之一。然而,蒸汽温度和压力等参数的提高导致机组关键高温管道的服役环境非常恶劣,特别是管道中存在裂纹、未焊透、焊接气孔和夹渣等各种缺陷,严重威胁机组的安全运行,需要对其进行科学精确的寿命评估。
[0003] 几十年来,对于高温下含裂纹构件,国外发展了多种高温蠕变寿命的评定规范和方法。蠕变孕育期是蠕变过程中经历时间最长的阶段,孕育期的准确预测对于高温结构的蠕变寿命预测具有重要意义;Davies等人基于韧性耗散模型提出的孕育期预测模型,考虑了蠕变过程应力变化的完整性,但是结构的残余应力对孕育期的影响并未得到研究;残余应力广泛的存在于加工制造的高温部件中,并且对部件的服役寿命造成了重大影响。大量针对于高温蠕变情况下残余应力(残余应力)的研究也广泛展开。因此建立考虑残余应力的蠕变孕育期预测模型,可以更加准确完整地评估复合加载结构的蠕变孕育期。

发明内容

[0004] 本发明在Davies工作的基础上,提出了考虑残余应力的蠕变孕育期预测模型。利用参考应立法,引入了弹性追随因子Z计算考虑残余应力的蠕变孕育期。使用紧凑拉伸试样(CT)通过预压缩产生残余应力,并施加主载荷进行蠕变实验。
[0005] 为实现本发明的目的所采用的技术方案是:
[0006] 本发明的塑性条件下含残余应力的蠕变孕育期预测方法,包括以下步骤:
[0007] S1:建立模型:所述模型包括CT试样本体,所述CT试样本体的中部前端设有槽,槽的后部设有缺口,CT试样本体上还设有上主载荷销孔、下主载荷销孔,上主载荷销孔、下主载荷销孔上下对应设置,分别设置在槽的上下两端;
[0008] S2:首先利用上圆销、下圆销对CT试样本体的上下两端进行预定大小的压缩加载,然后释放上圆销、下圆销,会在CT试样本体的缺口附近产生残余应力分布
[0009] S3:在含有残余应力的缺口处插入预制裂纹,以进行蠕变试验;
[0010] S4:利用销子在上主载荷销孔、下主载荷销孔施加主载荷,进行高温蠕变试验;
[0011] S5:通过蠕变有限元模拟可以获得计算含残余应力CT试样孕育期所需要的必要参数,在塑性条件下,计算孕育期主要包括以下步骤:
[0012] (1)首先计算复合加载下的应力强度因子,其计算公式为:
[0013]
[0014] (Ⅰ)中:
[0015]
[0016] 其中: 是模拟计算的只含有残余应力下的应力强度因子,单位为MPa·(m1/2);是主载荷应力强度因子,单位为MPa·(m1/2);P是主载荷,单位为N;B是试样厚度,单位为mm,Bn是试样净厚度,单位为mm;a/W是预制裂纹长度比率,a是预制裂纹长度,采用上主载荷销孔圆心到预制裂纹后端的平直线距离,单位为mm;W是名义试样宽度,采用上主载荷销孔圆心到CT试样本体后端的水平直线距离,单位为mm;f(a/W)是CT试样几何系数,只与a/W有关;V是无量纲的塑性相关项,计算如下:
[0017]
[0018] (Ⅱ)中:V0是无量纲参量,
[0019]
[0020] 其中: 是塑性残余应力强度因子,单位为MPa·(m1/2); 是弹性残余应力强度因子,单位为MPa·(m1/2), 利用JS计算,JS是残余应力场下断裂参量,单位为MPa·m;
[0021]
[0022] 其中:E′是有效弹性模量:E'=E/(1-ν2),E是弹性模量,ν是泊松比,E、ν参见文献:(Zhao L,Jing H,Xu L,Han Y,Xiu J.Evaluation of constraint effects on creep crack growth by experimental investigation and numerical simulation.Engng Fract Mech 2012;96:251–66.), 和JS都利用有限元模拟结果提取;
[0023] (Ⅱ)中:Lr是无量纲参量,描述主载荷幅度:
[0024]
[0025] 其中:σy是屈服强度,与材料属性有关,单位为MPa,参见文献:(Zhao L,Jing H,Xu L,Han Y,Xiu J.Evaluation of constraint effects on creep crack growth by experimental investigation and numerical simulation.Engng Fract Mech 2012;96:251–66.);σrefP是主载荷参考应力,单位为MPa,用下式计算:
[0026]
[0027] 其中:nL为无量纲裂纹深宽比参数,通过下式计算:
[0028]
[0029] 常数
[0030] (Ⅱ)中:
[0031] 是弹性主载荷应力强度因子, 是塑性主载荷应力强度因子,单位为MPa·(m1/2); 利用有限元模拟结果计算:
[0032] (Ⅱ)中:β描述残余应力的幅度,是无量纲参量;
[0033]
[0034] σrefS是二次载荷参考应力,利用有限元模拟计算,
[0035] (Ⅱ)中:Z为无量纲的弹性追随因子,从有限元模拟结果中提取出应力应变关系,取等效蠕变应变增量 与等效弹性应变增量 的比值:
[0036]
[0037] (2)计算塑性复合应力场下J积分数值,其计算公式为:
[0038]
[0039] 其中:KI是复合应力强度因子,E′是有效弹性模量:E'=E/(1-ν2),ν是泊松比,E是弹性模量,ν是泊松比,E和ν二者均参见文献:(Zhao L,Jing  H,Xu L,Han  Y,Xiu J.Evaluation of constraint effects on creep crack growth by experimental investigation and numerical simulation.Engng Fract Mech 2012;96:251–66.);
[0040] (3)然后计算塑性应力场下孕育期时间tiHRR,其计算公式为:
[0041]
[0042] (Ⅲ)中:εcrit是单轴蠕变韧性,与材料属性有关,单位为1, 是蠕变应变变化率,与材料高温蠕变属性有关,εcrit和 参见文献:(Zhao L,Jing H,Xu L,Han Y,Xiu J.Evaluation of constraint effects on creep crack growth by experimental investigation and numerical simulation.Engng Fract Mech 2012;96:251–66.),σP0是标准化应力,单位为MPa,εP0是标准化应变,单位为1,α为应变硬化系数,N为应变硬化指数,σP0、εP0、α和N参见文献:(Zhao L,Xu L,Han Y,Jing H.Two-parameter characterization of constraint effect induced by specimen size on creep crack growth.Engng Fract Mech 2012;96:251–66.)。IN是与N有关的无量纲函数, 是与θ和N有关的无量纲函数,具体值可以查阅文献得到:(Shih,C.F..1983.Tables of Hutchinson-Rice-Rosengren Singular Field Quantities.Brown University Technical Report,MRL E-147.)
[0043] (Ⅲ)中:d(mm)是判定蠕变萌生发生时裂尖前蠕变损伤达到1所延伸的距离,即蠕变萌生发生的临界距离。
[0044] (Ⅲ)中:MSFHRR为塑性条件下多轴应力因子,根据Cocks and Ashby关系式计算:
[0045]
[0046] n为无量纲的蠕变应力硬化指数,sinh是双曲正弦函数,hHRR为塑性应力三轴度,在塑性应力状态下:
[0047]
[0048] 其中:θ是裂纹尖端度,N是无量纲的应力硬化指数, 和 是与θ和N有关的无量纲函数,具体值可以查阅文献得到:(Shih,C.F..1983.Tables of Hutchinson-Rice-Rosengren Singular Field Quantities.Brown University Technical Report,MRL E-147.)。
[0049] 优选的,d取所研究材料的晶粒尺寸。
[0050] 优选的,Bn=B。
[0051] 优选的,所述有限元模拟采用ABAQUS6.14进行计算模拟,σrefS、 JS、 Z的提取过程包括以下步骤:
[0052] (5)首先建立预压缩加载的CT试样的有限元模型,在材料属性模里设置弹性塑性参数,在载荷模块里设置压缩载荷,以及拘束条件:包括对称条件和固定条件。在接触模块里设置压缩圆销与试样上下表面的的刚性接触,在分析步模块里设置好输出参量:应力值,在网格模块划分网格;
[0053] (6)在作业模块提交任务计算,获得残余应力的计算结果,结果文件中,从场变量S中可以直接提取二次载荷参考应力σref;
[0054] (7)建立相同尺寸的试样模型,进行主载荷拉伸试验,在材料属性模块里设置高温下的弹性塑性蠕变参数,在网格模块划分网格,在接触模块里设置拉伸销子与销孔的刚性接触,并在模型中插入预制裂纹,在分析步模块里设置好输出参量:应力应变值、应力强度因子K值、断裂参量J积分值,在载荷模块里设置拉伸载荷,以及拘束条件:包括对称条件和固定条件,在预加载应力场中导入上步算好的残余应力;
[0055] (8)在作业模块提交任务计算,获得含有残余应力的蠕变拉伸实验计算结果,结果文件中,在插入裂纹后还未施加拉伸载荷时刻,从历史变量中可以获取弹性残余应力强度因子 以及残余应力断裂参量JS,在施加拉伸载荷的初始时刻,可以获取塑性主载荷强度因子 从历史变量中可以获取等效应力随总应变增量的变化曲线,从曲线中得到等效蠕变应变增量, 等效弹性应变增量 进而得到弹性追随因子Z。
[0056] 与现有技术相比,本发明的有益效果是:
[0057] 本发明提出了修正的塑性条件下蠕变孕育期预测模型,与现有模型相比,本设计方法能够将原有的预测模型扩展到含残余应力的模型中,从而提出一种简化的塑性条件下蠕变孕育期预测方法,因此能够在结构中简洁有效的预测出塑性条件下蠕变孕育期。附图说明
[0058] 图1紧凑拉伸试样(CT)预压缩示意图;
[0059] 图2蠕变裂纹萌生临界条件示意图;
[0060] 图3应力应变关系曲线。

具体实施方式

[0061] 以下结合附图和具体实施例对本发明作进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
[0062] 本实施例选取P92高温耐热,以B=20mm,W=40mm,a/W=0.5的CT试样作为研究对象,以预加载为12000N和主载荷P=12000N作为研究载荷。其主要材料属性见下表:
[0063]
[0064] 其中:E-16为10的-16次方。
[0065] 本发明的塑性条件下含残余应力的蠕变孕育期预测方法,包括以下步骤:
[0066] S1:建立如图1所示的模型:所述模型包括CT试样本体2,所述CT试样本体2的中部前端设有槽4,槽4的后部设有缺口5,CT试样本体1上还设有上主载荷销孔3、下主载荷销孔7,上主载荷销孔3、下主载荷销孔7上下对应设置,分别设置在槽4的上下两端;
[0067] S2:首先利用上圆销1、下圆销8对CT试样本体2进行预定大小的压缩加载,然后释放上圆销1、下圆销8,会在CT试样本体2的缺口5附近产生一定的残余应力分布;
[0068] S3:在含有残余应力的缺口处插入预制裂纹6,槽4、缺口5、预制裂纹6在同一平面上,以进行蠕变试验;
[0069] S4:利用销子在上主载荷销孔3、下主载荷销孔7施加主载荷,进行高温蠕变试验;
[0070] S5:通过蠕变有限元模拟可以获得计算含残余应力CT试样本体2孕育期所需要的必要参数。在塑性条件下,计算孕育期主要包括以下步骤:
[0071] (1)首先计算各个参量:
[0072] (a)弹性主载荷强度因子:
[0073]
[0074]
[0075] 通过有限元结果提取以下几个数据:
[0076] ⅰ)首先按照尺寸,建立预压缩加载的CT试样的有限元模型。在材料属性模块里设置弹性塑性参数。在载荷模块里设置压缩载荷,以及拘束条件,所述拘束条件包括对称条件和固定条件,在接触模块里设置压缩圆销与试样上下表面的的刚性接触,在分析步模块里设置好输出参量:应力值,在网格模块划分网格;
[0077] ⅱ)在作业模块提交任务计算,获得残余应力的计算结果。结果文件中,从场变量中可以直接提取二次载荷参考应力
[0078] ⅲ)建立相同尺寸的试样模型,进行主载荷拉伸试验,可参见图1。在材料属性模块里设置高温下的弹性塑性蠕变参数,在网格模块划分网格,在接触模块里设置拉伸销子与销孔的刚性接触,并在模型中插入预制裂纹,在分析步模块里设置好输出参量:应力值,应力强度因子K值,断裂参量J积分值,在载荷模块里设置拉伸载荷,以及拘束条件:包括对称条件和固定条件,在预加载应力场中导入上步算好的残余应力;
[0079] ⅳ)在作业模块提交任务计算,获得含有残余应力的蠕变拉伸实验计算结果,结果文件中,在插入裂纹后还未施加拉伸载荷时刻,从历史变量中可以获取弹性残余应力强度因子 以及残余应力断裂参量JS=0.013MPa·m,可以计算得塑性残余应力强度因子: 在施加拉伸载荷的初始时刻,可以获取塑性
主载荷强度因子 从历史变量中可以获取等效应力随总应变增量
的变化曲线,如图3所示,从曲线中得到等效蠕变应变增量 等效弹性应变增量 进而得到弹性追随因子Z计算方法。
[0080] (b)主载荷参考应力:
[0081]
[0082] (c)主载荷幅度:
[0083]
[0084]
[0085] (d)残余应力参考应力:
[0086] 残余应力的幅度:
[0087]
[0088] (e)弹性追随因子:
[0089]
[0090]
[0091]
[0092] (f)塑性相关项:
[0093]
[0094]
[0095] (2)所以,复合加载下的应力强度因子
[0096] 复合加载下的J积分值:
[0097] (3)然后计算塑性应力场下发生的萌生:
[0098] (a)查表得: IN=4.49,P92钢的材料参数n=5.23,εcrit=0.2;N=11,
[0099] 应力三轴度:
[0100] 多轴应力因子:
[0101] d(mm)是判定蠕变萌生发生时裂尖前蠕变损伤达到1所延伸的距离,即蠕变萌生发生的临界距离,一般取所研究材料的晶粒尺寸,如图2所示。
[0102] 塑性条件下孕育期:
[0103] 以上所述仅是本发明的优选实施方式,应当指出的是,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈