CN109933823A - 塑性瞬态蠕变条件下耦合残余应力和拘束效应的蠕变孕育期预测方法 - Google Patents
塑性瞬态蠕变条件下耦合残余应力和拘束效应的蠕变孕育期预测方法 Download PDFInfo
- Publication number
- CN109933823A CN109933823A CN201711354401.7A CN201711354401A CN109933823A CN 109933823 A CN109933823 A CN 109933823A CN 201711354401 A CN201711354401 A CN 201711354401A CN 109933823 A CN109933823 A CN 109933823A
- Authority
- CN
- China
- Prior art keywords
- stress
- creep
- unit
- mpa
- residual stress
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种塑性瞬态蠕变条件下耦合残余应力和拘束效应的蠕变孕育期预测方法,在Davies工作的基础上,提出了耦合残余应力和拘束效应的蠕变孕育期预测模型。利用韧性耗散损伤模型,引入了与载荷无关拘束参量Q*计算考虑拘束效应的蠕变孕育期。使用紧凑拉伸试样(CT)施加主载荷进行蠕变模拟实验,本发明的有益效果:能够在结构中简洁有效的预测出塑性瞬态蠕变条件下蠕变孕育期。
Description
技术领域
本发明涉及耦合残余应力和拘束效应高温结构在塑性瞬态蠕变条件下的蠕变孕育期工程临界评价,就是确定当结构中存在表面裂纹且处在塑性瞬态蠕变条件下时,评价此高温结构的蠕变裂纹萌生寿命。
背景技术
燃煤为主的能源结构是我国雾霾天气的主要成因之一,而燃煤发电是我国目前最主要的发电方式,该趋势将长期存在。因此,除了改变能源结构,发展高效洁净的超超临界(USC)机组是节能减排的重要途径之一。然而,蒸汽温度和压力等参数的提高导致机组关键高温管道的服役环境非常恶劣,特别是管道中存在裂纹、未焊透、焊接气孔和夹渣等各种缺陷,严重威胁机组的安全运行,需要对其进行科学精确的寿命评估。
几十年来,对于高温下含裂纹构件,国外发展了多种高温蠕变寿命的评定规范和方法。蠕变孕育期是蠕变过程中经历时间最长的阶段,孕育期的准确预测对于高温结构的蠕变寿命预测具有重要意义;Davies等人基于韧性耗散模型提出的孕育期预测模型,考虑了蠕变过程应力变化的完整性,但是结构的残余应力以及拘束效应对孕育期的影响并未得到研究;残余应力、拘束效应广泛的存在于加工制造的高温部件中,并且对部件的服役寿命造成了重大影响。大量针对于高温蠕变情况下残余应力和拘束效应的研究也广泛展开。因此建立耦合残余应力和拘束效应的蠕变孕育期预测模型,可以更加准确完整地评估复合加载结构的蠕变孕育期。
发明内容
本发明的目的是针对现有技术中存在的技术缺陷,而提供一种发明名称。
为实现本发明的目的所采用的技术方案是:
本发明在Davies工作的基础上,提出了耦合残余应力和拘束效应的蠕变孕育期预测模型。利用参考应立法,引入了弹性追随因子Z计算考虑残余应力的蠕变孕育期。使用紧凑拉伸试样(CT)通过预压缩产生残余应力,并施加主载荷进行蠕变实验。
为实现本发明的目的所采用的技术方案是:
本发明的塑性瞬态蠕变条件下耦合残余应力和拘束效应高温结构的蠕变孕育期预测方法,包括以下步骤:
S1:建立模型,包括CT试样本体,CT试样本体的中部前端设有槽,槽的后部设有缺口,CT试样本体上还设有上主载荷销孔、下主载荷销孔,上主载荷销孔、下主载荷销孔上下对称设置,分别设置在槽的上下两端。
S2:首先利用上圆销、下圆销对CT试样本体进行预定大小的压缩加载,然后释放上圆销、下圆销,会在CT试样本体的缺口附近产生残余应力分布;
S3:在含有残余应力的缺口处插入预制裂纹,以进行蠕变试验;
S4:利用销子在上主载荷销孔、下主载荷销孔施加主载荷,进行高温蠕变试验;
S5:通过蠕变有限元模拟可以获得计算含残余应力CT试样孕育期所需要的必要参数,在塑性瞬态蠕变条件下,如图4所示,研究点的初始应力是塑性应力状态,到达转换时间tHRR-RR后进入到瞬态蠕变应力状态,计算孕育期主要包括以下步骤:
(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);是弹性残余应力强度因子,单位为MPa·(m1/2),利用JS计算,JS是残余应力场下断裂参量,单位为MPa·m:
其中:E′是有效弹性模量:E'=E/(1-ν2),E是弹性模量,ν是泊松比,E和ν二者均参见文献:(Zhao L,Jing H,Xu L,Han Y,Xiu J.Evaluation of constraint effects oncreep crack growth by experimental investigation and numericalsimulation.Engng Fract Mech 2012;96:251–66.),和JS都利用有限元模拟结果提取;
(Ⅱ)中:Lr是无量纲参量,描述主载荷幅度:
其中:σy是屈服强度,单位为MPa,参见文献:(Zhao L,Jing H,Xu L,Han Y,XiuJ.Evaluation of constraint effects on creep crack growth by experimentalinvestigation and numerical simulation.Engng Fract Mech 2012;96:251–66.);是主载荷参考应力,单位为MPa,用下式计算:
其中:nL为无量纲裂纹深宽比参数,通过下式计算:
常数
(Ⅱ)中:
其中:是弹性主载荷应力强度因子,单位为MPa·(m1/2),是塑性主载荷应力强度因子,单位为MPa·(m1/2);利用有限元模拟结果计算:
(Ⅱ)中:β描述残余应力的幅度,是无量纲参量;
是二次载荷参考应力,单位为MPa,利用有限元模拟计算;
(Ⅱ)中:Z为无量纲的弹性追随因子,从有限元模拟结果中提取出应力应变关系,取等效蠕变应变增量与等效弹性应变增量的比值:
(2)计算稳态蠕变复合应力场下C*积分数值,其计算公式为:
其中:A是蠕变硬化系数,单位为MPa-n·h-1,A参见文献:(Zhao L,Jing H,Xu L,Han Y,Xiu J.Evaluation of constraint effects on creep crack growth byexperimental investigation and numerical simulation.Engng Fract Mech 2012;96:251–66.),KI是复合应力强度因子,单位为MPa·(m1/2),是初始参考应力,单位为MPa;
(3)然后计算裂尖参量C(t),它是反映瞬态蠕变过程的载荷环路积分值,单位为MPa·mm·(h-1),利用参考应力法计算:
(Ⅴ)中:σref是总参考应力,单位为MPa,利用以下积分式计算:
其中:是总参考应变速率,单位为h-1,是主载荷参考应变速率,单位为h-1,
(Ⅴ)中:εref是总参考应变,利用下式计算:
εref=ε0 ref+A∫σn refdt
其中:ε0 ref是初始参考应变,ε0 ref通过有限元模拟提取。
(4)计算瞬态蠕变条件下的拘束参量Q*RR,其计算公式为:
是利用有限元计算得出的裂纹前沿处的张开应力值,单位是Mpa,σ0是材料的屈服强度,单位是MPa,参见文献:(Zhao L,Xu L,Han Y,Jing H.Two-parametercharacterization of constraint effect induced by specimen size on creep crackgrowth.Engng Fract Mech 2012;96:251–66.),L是标量距离,取1mm;
(Ⅵ)中:σ22是利用HRR应力场计算得出的裂纹前沿的张开应力值,单位是MPa,
其中:r是裂纹后部尖端到裂纹前沿研究点的间距,单位是mm,θ是裂纹尖端角度,是蠕变应变变化率,单位为h-1,与材料高温蠕变属性有关,n为无量纲的蠕变应力硬化指数,n和参见文献:(Zhao L,Jing H,Xu L,Han Y,Xiu J.Evaluation of constrainteffects on creep crack growth by experimental investigation and numericalsimulation.Engng Fract Mech 2012;96:251–66.),In是与n有关的无量纲函数,是与θ和n有关的无量纲函数,In和具体值可以查阅文献得到:Shih,C.F..1983.Tables of Hutchinson-Rice-Rosengren Singular Field Quantities.BrownUniversity Technical Report,MRL E-147.
(5)计算瞬态蠕变等效应力其计算公式为:
其中:是与θ和n有关的无量纲函数,具体值可以查阅文献得到:Shih,C.F..1983.Tables of Hutchinson-Rice-Rosengren Singular Field Quantities.BrownUniversity Technical Report,MRL E-147.
计算塑性等效应力其计算公式为:
σP0是标准化应力,单位为MPa,εP0是标准化应变,单位为1,α为应变硬化系数,N为应变硬化指数,σP0,εP0,α和N参见文献:(Zhao L,Xu L,Han Y,Jing H.Two-parametercharacterization of constraint effectinduced by specimen size on creep crackgrowth.Engng Fract Mech 2012;96:251–66.),IN是与N有关的无量纲函数,是与裂纹尖端角度θ和应变硬化指数N有关的无量纲函数,IN和可查表获得Shih,C.F..1983.Tables of Hutchinson-Rice-Rosengren Singular Field Quantities.BrownUniversity Technical Report,MRL E-147.
(6)利用MATALAB软件计算转换时间tHRR-RR:在此时刻满足:
塑性阶段损伤累计值:
MSFHRR为塑性条件下多轴应力因子,根据Cocks and Ashby关系式计算:
sinh是双曲正弦函数,hHRR为塑性应力三轴度,在塑性应力状态下:
其中:和是与θ和N有关的无量纲函数,具体值可以查阅文献得到:Shih,C.F..1983.Tables of Hutchinson-Rice-Rosengren Singular FieldQuantities.Brown University Technical Report,MRL E-147.。
(7)然后计算瞬态蠕变应力场下孕育期时间ti HRR-RR,其计算公式为:
(Ⅷ)中:d(mm)是判定蠕变萌生发生时裂尖前蠕变损伤达到1所延伸的距离,即蠕变萌生发生的临界距离。
(Ⅷ)中:MSFRR为塑性条件下多轴应力因子,根据Cocks and Ashby关系式计算:
sinh是双曲正弦函数,hRR为瞬态蠕变应力三轴度,在塑性应力状态下:
其中平均应力单位是MPa,其计算公式为:
其中:σ11和σ33是利用RRss应力场计算得出的裂纹前沿的应力值,单位是MPa,
其中:是与θ和n有关的无量纲函数,具体值可以查阅文献得到:Shih,C.F..1983.Tables of Hutchinson-Rice-Rosengren Singular FieldQuantities.Brown University Technical Report,MRL E-147.。
优选的,d取所研究材料的晶粒尺寸。
优选的,Bn=B。
优选的,所述有限元模拟采用ABAQUS6.14进行计算模拟,JS、ε0 ref的提取过程包括以下步骤:
(5)首先建立预压缩加载的CT试样的有限元模型,在材料属性模块里设置弹性塑性参数,在载荷模块里设置压缩载荷,以及拘束条件:包括对称条件和固定条件。在接触模块里设置压缩圆销与试样上下表面的的刚性接触,在分析步模块里设置好输出参量:应力值,在网格模块划分网格;
(6)在作业模块提交任务计算,获得残余应力的计算结果,结果文件中,从场变量中可以直接提取二次载荷参考应力σref S;
(7)建立相同尺寸的试样模型,进行主载荷拉伸试验,在材料属性模块里设置高温下的弹性塑性蠕变参数,在网格模块划分网格,在接触模块里设置拉伸销子与销孔的刚性接触,并在模型中插入预制裂纹,在分析步模块里设置好输出参量:应力应变值、应力强度因子K值、断裂参量J积分值,在载荷模块里设置拉伸载荷,以及拘束条件:包括对称条件和固定条件,在预加载应力场中导入上步算好的残余应力;
(8)在作业模块提交任务计算,获得含有残余应力的蠕变拉伸实验计算结果,结果文件中,在插入裂纹后还未施加拉伸载荷时刻,从场变量中可以获取初始参考应变从历史变量中可以获取弹性残余应力强度因子以及残余应力断裂参量JS,在施加拉伸载荷的初始时刻,可以获取塑性主载荷强度因子从历史变量中可以获取等效应力随总应变增量的变化曲线,从曲线中得到等效蠕变应变增量等效弹性应变增量进而得到弹性追随因子Z计算方法。
与现有技术相比,本发明的有益效果是:
本发明提出了耦合残余应力和拘束效应的弹性瞬态蠕变条件下蠕变孕育期预测模型,与现有模型相比,本设计方法能够将原有的预测模型扩展到含残余应力的模型中,从而提出一种简化的瞬态蠕变条件下蠕变孕育期预测方法,因此能够在结构中简介有效的预测出弹性瞬态蠕变条件下蠕变孕育期。
附图说明
图1紧凑拉伸试样(CT)预压缩示意图;
其中:1-上圆销,2-CT试样本体,3-上主载荷销孔,4-槽,5-缺口,6-预制裂纹,7-下主载荷销孔,8-下圆销。
图2是蠕变裂纹萌生临界条件示意图。
图3是应力应变关系曲线。
图4初始应力向瞬态蠕变转换时间。
具体实施方式
以下结合附图和具体实施例对本发明作进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明的塑性瞬态蠕变条件下耦合残余应力和拘束效应高温结构的蠕变孕育期预测模型,包括CT试样本体2,CT试样本体的中部前端设有槽4,槽的后部设有缺口5,槽4、缺口5在同一平面上,CT试样本体上还设有上主载荷销孔3、下主载荷销孔7,上主载荷销孔、下主载荷销孔上下对称设置,分别设置在槽的上下两端。
选取P92高温耐热钢,以B=20mm,W=40mm,a/W=0.5的CT试样作为研究对象,以预加载为12000N和主载荷P=12000N作为研究载荷。其主要材料属性见下表:
其中E-16为10的-16次方。
本发明的塑性瞬态蠕变条件下耦合残余应力和拘束效应高温结构的蠕变孕育期预测方法,包括以下步骤:
S1:首先利用上圆销1、下圆销8对CT试样本体进行预定大小的压缩加载,然后释放上圆销、下圆销,会在CT试样本体的缺口5附近产生残余应力分布;
S2:在含有残余应力的缺口5处插入预制裂纹6,以进行蠕变试验;
S3:利用销子在上主载荷销孔3、下主载荷销孔7施加主载荷,进行高温蠕变试验;
S4:通过蠕变有限元模拟可以获得计算含残余应力CT试样孕育期所需要的必要参数,在塑性瞬态蠕变条件下,如图4所示,研究点的初始应力是塑性应力状态,到达转换时间tHRR-RR后进入到瞬态蠕变应力状态。
所述有限元模拟采用ABAQUS6.14进行计算模拟,JS、ε0 ref的提取过程包括以下步骤:
(ⅰ)首先建立如图1所示的预压缩加载的CT试样的有限元模型,在材料属性模块里设置弹性塑性参数,在载荷模块里设置压缩载荷,以及拘束条件:包括对称条件和固定条件。在接触模块里设置压缩圆销与试样上下表面的的刚性接触,在分析步模块里设置好输出参量:应力值,在网格模块划分网格;
(ⅱ)在作业模块提交任务计算,获得残余应力的计算结果,结果文件中,从场变量中可以直接提取二次载荷参考应力
(ⅲ)建立如图1所示的试样模型,进行主载荷拉伸试验,在材料属性模块里设置高温下的弹性塑性蠕变参数,在网格模块划分网格,在接触模块里设置拉伸销子与销孔的刚性接触,并在模型中插入预制裂纹,在分析步模块里设置好输出参量:应力应变值、应力强度因子K值、断裂参量J积分值,在载荷模块里设置拉伸载荷,以及拘束条件:包括对称条件和固定条件,在预加载应力场中导入上步算好的残余应力;
(ⅳ)在作业模块提交任务计算,获得含有残余应力的蠕变拉伸实验计算结果,结果文件中,在插入裂纹后还未施加拉伸载荷时刻,在插入裂纹后还未施加拉伸载荷时刻,从场变量中可以获取初始参考应变从历史变量中可以获取弹性残余应力强度因子以及残余应力断裂参量JS=0.013MPa·m,可以计算得塑性残余应力强度因子:在施加拉伸载荷的初始时刻,可以获取塑性主载荷强度因子从历史变量提取出C(t)随时间的变化曲线,如图4所示,从历史变量中可以获取等效应力随总应变增量的变化曲线,如图3所示,从曲线中得到等效蠕变应变增量等效弹性应变增量进而得到弹性追随因子Z计算方法。
计算孕育期主要包括以下步骤:
(1)首先计算各个参量:
(a)弹性主载荷强度因子:
(b)主载荷参考应力:
(c)主载荷幅度:
(d)残余应力参考应力:
残余应力的幅度:
(e)弹性追随因子:
(f)塑性相关项:
(2)所以,复合加载下的应力强度因子
初始参考应力:
稳态蠕变复合应力场下C*积分数值为:
复合加载下的J积分值:
(3)场变量可以获取应力值
查表得:In=4.99,P92钢的材料参数εcrit=0.2;n=5.23,在计算蠕变应力以及拘束时候,我们取裂尖前的距离r=d=0.05mm。
查表得:塑性等效应力
查表得:
瞬态蠕变等效应力
转换时间tHRR-RR:利用以及MATALAB计算得:tHRR-RR=2h
查表得:
塑性阶段损伤累计值:
(4)然后计算瞬态蠕变应力场下发生的萌生:
查表得:
平均应力:
应力三轴度:
多轴应力因子:是判定蠕变萌生发生时裂尖前蠕变损伤达到1所延伸的距离,即蠕变萌生发生的临界距离,一般取所研究材料的晶粒尺寸,如图2所示。
瞬态蠕变条件下孕育期:
利用MATALAB解积分得:ti HRR-RR=1822h。
以上所述仅是本发明的优选实施方式,应当指出的是,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (4)
1.塑性瞬态蠕变条件下耦合残余应力和拘束效应高温结构的蠕变孕育期预测方法,其特征在于,包括以下步骤:
S1:建立模型,包括CT试样本体,CT试样本体的中部前端设有槽,槽的后部设有缺口,CT试样本体上还设有上主载荷销孔、下主载荷销孔,上主载荷销孔、下主载荷销孔上下对称设置,分别设置在槽的上下两端;
S2:首先利用上圆销、下圆销对CT试样本体进行预定大小的压缩加载,然后释放上圆销、下圆销,会在CT试样本体的缺口附近产生残余应力分布;
S3:在含有残余应力的缺口处插入预制裂纹,以进行蠕变试验;
S4:利用销子在上主载荷销孔、下主载荷销孔施加主载荷,进行高温蠕变试验;
S5:通过蠕变有限元模拟可以获得计算含残余应力CT试样孕育期所需要的必要参数,在塑性瞬态蠕变条件下,如图4所示,研究点的初始应力是塑性应力状态,到达转换时间tHRR-RR后进入到瞬态蠕变应力状态,计算孕育期主要包括以下步骤:
(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);是弹性残余应力强度因子,单位为MPa·(m1/2),利用JS计算,JS是残余应力场下断裂参量,单位为MPa·m:
其中:E′是有效弹性模量:E'=E/(1-ν2),E是弹性模量,ν是泊松比,和JS都利用有限元模拟结果提取;
(Ⅱ)中:Lr是无量纲参量,描述主载荷幅度:
其中:σy是屈服强度,单位为MPa;是主载荷参考应力,单位为MPa,用下式计算:
其中:nL为无量纲裂纹深宽比参数,通过下式计算:
常数
(Ⅱ)中:
其中:是弹性主载荷应力强度因子,单位为MPa·(m1/2),是塑性主载荷应力强度因子,单位为MPa·(m1/2);利用有限元模拟结果计算:
(Ⅱ)中:β描述残余应力的幅度,是无量纲参量;
是二次载荷参考应力,单位为MPa,利用有限元模拟计算;
(Ⅱ)中:Z为无量纲的弹性追随因子,从有限元模拟结果中提取出应力应变关系,取等效蠕变应变增量与等效弹性应变增量的比值:
(2)计算稳态蠕变复合应力场下C*积分数值,其计算公式为:
其中:A是蠕变硬化系数,单位为MPa-n·h-1,KI是复合应力强度因子,单位为MPa·(m1 /2),是初始参考应力,单位为MPa;
(3)然后计算裂尖参量C(t),它是反映瞬态蠕变过程的载荷环路积分值,单位为MPa·mm·(h-1),利用参考应力法计算:
(Ⅴ)中:σref是总参考应力,单位为MPa,利用以下积分式计算:
其中:是总参考应变速率,单位为h-1,是主载荷参考应变速率,单位为h-1,
(Ⅴ)中:εref是总参考应变,利用下式计算:
εref=ε0 ref+A∫σn refdt
其中:ε0 ref是初始参考应变,ε0 ref通过有限元模拟提取,
(4)计算瞬态蠕变条件下的拘束参量Q*RR,其计算公式为:
是利用有限元计算得出的裂纹前沿处的张开应力值,单位是Mpa,σ0是材料的屈服强度,单位是MPa,L是标量距离,取1mm;
(Ⅵ)中:σ22是利用HRR应力场计算得出的裂纹前沿的张开应力值,单位是MPa,
其中:r是裂纹后部尖端到裂纹前沿研究点的间距,单位是mm,θ是裂纹尖端角度,是蠕变应变变化率,单位为h-1,与材料高温蠕变属性有关,n为无量纲的蠕变应力硬化指数,In是与n有关的无量纲函数,是与θ和n有关的无量纲函数,
(5)计算瞬态蠕变等效应力其计算公式为:
其中:是与θ和n有关的无量纲函数,
计算塑性等效应力其计算公式为:
σP0是标准化应力,单位为MPa,εP0是标准化应变,单位为1,α为应变硬化系数,N为应变硬化指数,IN是与N有关的无量纲函数,是与裂纹尖端角度θ和应变硬化指数N有关的无量纲函数;
(6)利用MATALAB软件计算转换时间tHRR-RR:在此时刻满足:
塑性阶段损伤累计值:
MSFHRR为塑性条件下多轴应力因子,根据Cocks and Ashby关系式计算:
sinh是双曲正弦函数,hHRR为塑性应力三轴度,在塑性应力状态下:
其中:和是与θ和N有关的无量纲函数;
(7)然后计算瞬态蠕变应力场下孕育期时间ti HRR-RR,其计算公式为:
(Ⅷ)中:d(mm)是判定蠕变萌生发生时裂尖前蠕变损伤达到1所延伸的距离,即蠕变萌生发生的临界距离,
(Ⅷ)中:MSFRR为塑性条件下多轴应力因子,根据Cocks and Ashby关系式计算:
sinh是双曲正弦函数,hRR为瞬态蠕变应力三轴度,在塑性应力状态下:
其中平均应力单位是MPa,其计算公式为:
其中:σ11和σ33是利用RRss应力场计算得出的裂纹前沿的应力值,单位是MPa,
其中:是与θ和n有关的无量纲函数。
2.如权利要求1所述的塑性瞬态蠕变条件下耦合残余应力和拘束效应高温结构的蠕变孕育期预测方法,其特征在于,d取所研究材料的晶粒尺寸。
3.如权利要求1所述的塑性瞬态蠕变条件下耦合残余应力和拘束效应高温结构的蠕变孕育期预测方法,其特征在于,Bn=B。
4.如权利要求1所述的塑性瞬态蠕变条件下耦合残余应力和拘束效应高温结构的蠕变孕育期预测方法,其特征在于,所述有限元模拟采用ABAQUS6.14进行计算模拟,JS、ε0 ref的提取过程包括以下步骤:
(1)首先建立预压缩加载的CT试样的有限元模型,在材料属性模块里设置弹性塑性参数,在载荷模块里设置压缩载荷,以及拘束条件:包括对称条件和固定条件,在接触模块里设置压缩圆销与试样上下表面的的刚性接触,在分析步模块里设置好输出参量:应力值,在网格模块划分网格;
(2)在作业模块提交任务计算,获得残余应力的计算结果,结果文件中,从场变量中可以直接提取二次载荷参考应力
(3)建立相同尺寸的试样模型,进行主载荷拉伸试验,在材料属性模块里设置高温下的弹性塑性蠕变参数,在网格模块划分网格,在接触模块里设置拉伸销子与销孔的刚性接触,并在模型中插入预制裂纹,在分析步模块里设置好输出参量:应力应变值、应力强度因子K值、断裂参量J积分值,在载荷模块里设置拉伸载荷,以及拘束条件:包括对称条件和固定条件,在预加载应力场中导入上步算好的残余应力;
(4)在作业模块提交任务计算,获得含有残余应力的蠕变拉伸实验计算结果,结果文件中,在插入裂纹后还未施加拉伸载荷时刻,从场变量中可以获取初始参考应变ε0 ref,从历史变量中可以获取弹性残余应力强度因子以及残余应力断裂参量JS,在施加拉伸载荷的初始时刻,可以获取塑性主载荷强度因子从历史变量中可以获取等效应力随总应变增量的变化曲线,从曲线中得到等效蠕变应变增量等效弹性应变增量进而得到弹性追随因子Z计算方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711354401.7A CN109933823B (zh) | 2017-12-15 | 2017-12-15 | 塑性瞬态蠕变条件下耦合残余应力和拘束效应的蠕变孕育期预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711354401.7A CN109933823B (zh) | 2017-12-15 | 2017-12-15 | 塑性瞬态蠕变条件下耦合残余应力和拘束效应的蠕变孕育期预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109933823A true CN109933823A (zh) | 2019-06-25 |
CN109933823B CN109933823B (zh) | 2022-11-15 |
Family
ID=66980606
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711354401.7A Active CN109933823B (zh) | 2017-12-15 | 2017-12-15 | 塑性瞬态蠕变条件下耦合残余应力和拘束效应的蠕变孕育期预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109933823B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101046682A (zh) * | 2007-04-28 | 2007-10-03 | 武汉科技大学 | 一种预测热轧含Nb带钢组织及力学性能的方法 |
CN106529017A (zh) * | 2016-11-08 | 2017-03-22 | 华南理工大学 | 一种高铬钢构件高温蠕变变形预测与蠕变损伤分析方法 |
-
2017
- 2017-12-15 CN CN201711354401.7A patent/CN109933823B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101046682A (zh) * | 2007-04-28 | 2007-10-03 | 武汉科技大学 | 一种预测热轧含Nb带钢组织及力学性能的方法 |
CN106529017A (zh) * | 2016-11-08 | 2017-03-22 | 华南理工大学 | 一种高铬钢构件高温蠕变变形预测与蠕变损伤分析方法 |
Non-Patent Citations (3)
Title |
---|
HANNA QUINTERO ETC.: ""Prediction of creep crack initiation behaviour in 316H stainless steel using stress dependent creep ductility"", 《INTERNATIONAL JOURNAL OF SOLIDS AND STRUCTURES》 * |
J.Z.HE ETC.: ""Effect of constraint on creep crack initiation time in test specimens in ASTM-E1457 standard"", 《ENGINEERING FRACTURE MECHANICS》 * |
宋晓梅: ""残余应力对不同几何和材料拘束试样蠕变裂纹起裂和扩展行为的影响"", 《中国优秀硕士学位论文全文数据库》 * |
Also Published As
Publication number | Publication date |
---|---|
CN109933823B (zh) | 2022-11-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108733860A (zh) | 塑性瞬态蠕变条件下考虑拘束效应的蠕变孕育期预测方法 | |
CN108732032A (zh) | 稳态蠕变条件下含残余应力的蠕变孕育期预测方法 | |
CN108732034A (zh) | 弹性瞬态蠕变条件下含残余应力的蠕变孕育期预测方法 | |
CN108733861A (zh) | 塑性条件下含残余应力的蠕变孕育期预测方法 | |
CN109933816A (zh) | 弹性瞬态蠕变条件下耦合残余应力和拘束效应的蠕变孕育期预测方法 | |
CN109933823A (zh) | 塑性瞬态蠕变条件下耦合残余应力和拘束效应的蠕变孕育期预测方法 | |
CN108731989A (zh) | 塑性瞬态蠕变条件下含残余应力的蠕变孕育期预测方法 | |
CN108732029A (zh) | 弹性条件下含残余应力的蠕变孕育期预测方法 | |
CN109933815A (zh) | 稳态蠕变条件下耦合残余应力和拘束效应的蠕变孕育期预测方法 | |
CN109932251A (zh) | 塑性条件下考虑与载荷无关的拘束参量的蠕变孕育期预测方法 | |
CN109932241A (zh) | 塑性条件下耦合残余应力和拘束效应的蠕变孕育期预测方法 | |
CN109932242A (zh) | 稳态蠕变条件下考虑与载荷无关的拘束参量的蠕变孕育期预测方法 | |
CN109933820A (zh) | 弹性条件下考虑与载荷无关的拘束参量的蠕变孕育期预测方法 | |
CN108733862A (zh) | 稳态蠕变条件下考虑拘束效应的蠕变孕育期预测方法 | |
CN109933817A (zh) | 弹性瞬态蠕变条件下考虑与载荷无关的拘束参量的蠕变孕育期预测方法 | |
CN109933822A (zh) | 塑性瞬态蠕变条件下考虑与载荷无关的拘束参量的蠕变孕育期预测方法 | |
CN108732033A (zh) | 弹性瞬态蠕变条件下考虑拘束效应的蠕变孕育期预测方法 | |
Liu et al. | A new energy gradient-based model for LCF life prediction of turbine discs | |
CN108732030A (zh) | 塑性条件下考虑拘束效应的蠕变孕育期预测方法 | |
CN108732031A (zh) | 弹性条件下考虑拘束效应的蠕变孕育期预测方法 | |
Sun et al. | Correlation of crack growth rate and stress ratio for fatigue damage containing very high cycle fatigue regime | |
Buch | Prediction of constant-amplitude fatigue life to failure under pulsating tension by use of the local-strain approach | |
BAGNOLO | Experimental and numerical analysis of the cyclic plastic behavior of a steel for oil and gas application in multiaxial condition | |
Huang et al. | The fatigue crack growth under compressive to compressive fluctuating loading | |
Shlyannikov et al. | Crack growth resistance characteristics determination algorithm for creep-fatigue interaction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |