CN110411864A - 一种基于蠕变激活能的高温蠕变寿命预测解析计算方法 - Google Patents
一种基于蠕变激活能的高温蠕变寿命预测解析计算方法 Download PDFInfo
- Publication number
- CN110411864A CN110411864A CN201810388386.6A CN201810388386A CN110411864A CN 110411864 A CN110411864 A CN 110411864A CN 201810388386 A CN201810388386 A CN 201810388386A CN 110411864 A CN110411864 A CN 110411864A
- Authority
- CN
- China
- Prior art keywords
- creep
- activation energy
- temperature
- life
- span
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/28—Investigating ductility, e.g. suitability of sheet metal for deep-drawing or spinning
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0001—Type of application of the stress
- G01N2203/0003—Steady
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/006—Crack, flaws, fracture or rupture
- G01N2203/0067—Fracture or rupture
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
- G01N2203/0071—Creep
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/022—Environment of the test
- G01N2203/0222—Temperature
- G01N2203/0226—High temperature; Heating means
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
Abstract
本发明提供了一种基于蠕变激活能的高温蠕变寿命预测解析计算方法,基于蠕变激活能的应力相关性,建立一种蠕变寿命预测方法,将推导出的参数k1,μ,D,m代入以下寿命预测模型:即可精准有效的预测出材料的高温蠕变寿命,解决不同应力水平下由于蠕变机制变化引起的蠕变寿命预测难题,提高长时寿命预测精度。
Description
技术领域
本发明属于材料科学与工程应用技术,具体的说一种基于蠕变激活能的高温蠕变寿命预测解析计算方法。
背景技术
为了响应国家政策、提高燃料利用效率、降低生产成本、保护生态环境,火电行业、核电行业、航空航天等领域不断提高构件的工作温度。但是在提高工作温度的同时,高温构件所处环境发生了变化,蠕变性能也随之改变。这就导致高温构件在实际使用过程中发生蠕变失效的概率大大提高,从而影响生命财产安全和正常的生产生活。因此,如何精确对材料蠕变过程的寿命进行预测,使高温构件更加安全稳定的运行,是一项非常重要的任务,对其的使用和设计具有十分深远的现实意义。
目前,国内外研究人员针对高温材料蠕变寿命预测进行了大量的研究和分析,主要从微观组织演变,宏观断裂模式等方面着手进行,并开展了一系列相关试验进行解释说明。对于蠕变寿命的预测,国内外通常采用的有基于持久强度试验的数据外推得到的以时间温度参数法和Robinson断裂法的参数唯象模型,基于微观损伤机理的蠕变损伤力学模型和基于应变的连续损伤力学模型。通过不同的模型,研究人员提出了多种不同的预测寿命的方法,建立了多种蠕变损伤本构模型。但是由于不同模型所涉及的反应机理的差别,不同的蠕变损伤本构模型需要通过不同的方法进行参数拟合,且考虑因素越多,所拟合的参数就越多,进而导致拟合过程十分繁琐、材料参数对材料性质和结构等有所依赖的问题较为突出。这在很大程度上限制了实际生产生活过程中对高温构建蠕变寿命预测的发展。因此,研究和提出新的预测蠕变断裂寿命方法,对于拓展寿命预测的研究领域、寻找更加适合某类型钢材的寿命预测理论,具有十分重要的现实意义。近几年,从自然界事物普遍适用的热力学基本定律推导而出的蠕变断裂寿命解析计算方法越来越受到研究学者们的重视,此类模型具有材料参数拟合方法简单,且一般情况下预测精度较高的特点,为高温构件损伤评估和寿命预测提供了一个新的研究方向。
发明内容
本发明的目的在于通过考虑蠕变激活能与应力相关性,解决由于蠕变机制转变引起的长期寿命预测精度降低的问题,为应用在高温高压环境中的大型重要设备上的承压构件的蠕变断裂寿命预测提供一种新的解析计算方法。为此,本发明提供一种通过计算材料的在各个应力下的蠕变激活能,结合中高应力范围内蠕变数据,揭示蠕变激活能与蠕变断裂时间的定量关系。
一种基于蠕变激活能理论的高温蠕变寿命预测解析计算方法,通过如下步骤实现:
步骤1,获取材料不同温度、不同应力水平下蠕变性能的数据,每个试验点包括材料的应力σ(单位为MPa)、温度T(单位为℃)、材料屈服强度σys(单位为MPa)、断裂时间tf(单位为h)、最小蠕变速率(单位为h-1)以及气体常数R(单位为J/(mol·K));
步骤2,把试验数据按照式利用数学分析软件,按最小二乘法回归,求出待定系数α和M;
步骤3,根据最小蠕变速率与应力、温度和蠕变激活能的关系:
两边取对数:
利用试验数据,绘制曲线和并采用最小二乘数法的斜率求到材料常数n和蠕变激活能Q*(单位为kJ/mol)的值。
步骤4,利用步骤3求得不同应力水平σ下的蠕变激活能Q*,采用最小二乘数法获得Q*=f(σ)=D*σ+m的D和m参数。
步骤5,将步骤1得到的数据和步骤3得到的蠕变激活能Q*按照式3
进行最小二乘数法进行拟合试验数据,获得材料系数k1与μ的值。
步骤6,将步骤1-5得到的参数k1,μ,D,m代入寿命预测模型,如式4所示:
考虑蠕变激活能和应力相关性:
在上述技术方案中,所述高温蠕变寿命预测解析计算方法适用应力水平为0.2σys-σys,其中σys为材料屈服强度。
在上述技术方案中,所述高温蠕变寿命预测解析计算方法适用作业温度为400-1200℃。
本发明的优势有以下几点:
1.考虑蠕变应力变化引起的蠕变激活能的变化;
2.建立一种基于蠕变激活能的高温蠕变寿命预测模型;
3.解决短时试验数据预测长时寿命的预测;
4.提高了蠕变寿命预测的精度,扩大高温金属材料的适用范围;
5.预测方法简单,所需要的数据均可由常规的材料蠕变性能的测试获得。
发明内容
图1是拟合Monkman-Grant模型中的参数图。
图2是计算的斜率,从而求得应力指数n。
图3是计算的斜率,从而求得蠕变激活能Q*。
图4是蠕变激活能Q*与应力σ的关系。
图5是线性拟合得k1与μ的值。
图6是基于蠕变激活能理论的寿命预测曲线与试验值的比较。
具体实施方式
下面结合具体实施例和附图,进一步阐述本发明。
本发明提供一种更加精确地高温蠕变寿命预测解析计算方法,具体步骤如下:
第一步,基于700℃、725℃、750℃三个温度下不同应力水平下的材料单轴蠕变试验:
试验按照GB/T2039-2012《金属拉伸蠕变持久试验方法》进行。试样尺寸:直径为5mm的标准圆棒试样,标距为50mm。试验设备为高温蠕变持久强度试验机。该试验机的组成如下:主机;加热炉;温度测控系统;变形测量系统。其负荷范围为0.3-30KN,载荷误差小于等于±1%。蠕变自动记录仪的量程为:0-10mm,测量误差不超过±0.1%。首先将试样安装在试验机上,安装好引伸计,检查试样的同轴度在规定的范围内,若超出应当按要求调整。试样安装好之后,先施加预载荷200N,并开始升温加热到预定温度然后保持60min。最后施加总载荷,记录断裂时间。在此试验中,试验温度为700℃和750℃,应力水平为:87-240Mpa。试验结束,通过数据整理,获得材料不同温度、不同应力水平下的蠕变性能参数,如应力(σ)、温度(T)、材料屈服强度(σys)、断裂时间(tf)、最小蠕变速率
表1 材料的最小蠕变应变速率数据表。
表1中,E-06是指×10-6,E-05是指×10-5,E-04是指×10-4
第二步,将所得到的蠕变性能参数,按照Monkman-Grant模型,如式(1)所示,通过最小二乘法进行回归拟合:
700℃、725℃、750℃三个温度下的参数的拟合过程图如图1所示。由表1所示的试验数据,通过最小二乘法拟合(可采用Matlab,Origin等软件)得到待定系数α和M,如表2所示。
表2 Monkman-Grant模型拟合的参数
第三步,根据最小蠕变速率与温度、蠕变激活能和应力的关系:
式中:为最小蠕变应变速率,A为与材料有关的常数,n为应力指数,Q*为蠕变激活能,R为气体常数(R=8.314,单位为J/(mol·K)),T为温度。
将式(2)两边均取对数,可得:
利用试验数据,绘制式(3)中的和曲线,采用最小二乘数法拟合获得的斜率(可采用Matlab,Origin等软件),即可得到材料常数n和蠕变激活能Q*的值。如图2和图3所示。
700℃时n值为7.7465,725℃时n值为6.043,750℃时n值为5.009。各个应力值对应的激活能Q*,如表3所示。
表3 计算得到的各个应力值对应的蠕变激活能Q*
第四步:利用第三步求得不同应力水平σ下的蠕变激活能Q*,采用最小二乘数法(可采用Matlab,Origin等软件)获得Q*=f(σ)=D*σ+m的D和m参数,如图4所示。
两者之间的线性方程Q*=f(σ),如式(4)所示:
Q*=-2.9712σ+1140.89
第五步:基于蠕变破断时间与应力、蠕变激活能之间的关系式,如式(5)所示。
式中:Q*是蠕变激活能,tf是断裂时间,R是气体常数,T为温度值,k1、μ是材料常数。
对式(5)进行取对数变换,可以得到下式:
ln[-ln(σ/σys)]=lnk1+μln[tf·exp(-Q*/RT)] (6)
做出700℃、725℃、750℃三个温度下ln[tf·exp(-Q*/RT)]与ln[-ln(σ/σTS)]的关系点图,线性拟合即可得出k1与μ的值,如图5所示。
至此可获得式(4)与(5)中的所有需要的参数值,如表4所示。
表4 基于激活能法蠕变寿命预测方法的所需参数值。
第六步:为寻找蠕变激活能的应力相关性,建立一种蠕变寿命预测方法,将式(5)进行变换,可以得出断裂时间与应力、蠕变激活能、温度等参量的公式,如式(7)所示。
将式(4)以及表4中的参数值代入式(7)中,即可得到断裂时间与应力的关系式,如式(8)所示。代入各个温度下的屈服强度σys、k1、μ的值,即可得到不同温度下不同应力下的蠕变断裂寿命预测解析模型。
利用以上得到的基于蠕变激活能理论推导得到的蠕变断裂寿命解析模型预测的寿命预测曲线,如图6所示,图中散点是文献(Chai G,Hernblom J,Peltola T,et al.Creepbehavior in a newly developed heat resistant austenitic stainless steel[J].BHM Berg-undMonatshefte,2015,160(9):400-405.)报道的蠕变寿命数值,曲线是利用本发明的预测方法模拟出的曲线。可以发现本发明的基于蠕变激活能理论的高温蠕变寿命预测解析计算方法,考虑蠕变激活能与应力相关性,解决蠕变机制变化引起长期寿命预测精度偏低的难题,能够简便并且较高精度地计算得到蠕变寿命,使应用具有更强的操作性和说服力。
Claims (3)
1.一种基于蠕变激活能理论的高温蠕变寿命预测解析计算方法,其特征在于,包括以下步骤:
步骤1,获取材料不同温度、不同应力水平下蠕变性能的数据,每个试验点包括材料的应力σ、温度T、材料屈服强度σys、断裂时间tf、最小蠕变速率以及气体常数R;
步骤2,把试验数据按照式按最小二乘法回归,求出待定系数α和M;
步骤3,根据最小蠕变速率与应力、温度和蠕变激活能的关系:
两边取对数:
利用试验数据,绘制曲线和并采用最小二乘数法的斜率求到材料常数n和蠕变激活能Q*的值;
步骤4,利用步骤3求得不同应力水平σ下的蠕变激活能Q*,采用最小二乘数法获得Q*=f(σ)=D*σ+m的D和m参数;
步骤5,将步骤1得到的数据和步骤3得到的蠕变激活能Q*按照式3
进行最小二乘数法进行拟合试验数据,获得材料系数k1与μ的值;
步骤6,将步骤1-5得到的参数k1,μ,D,m代入寿命预测模型,如式4所示:
2.如权利要求1所述的一种基于蠕变激活能理论的高温蠕变寿命预测解析计算方法,其特征在于,所述高温蠕变寿命预测解析计算方法适用应力水平为0.2σys-σys,其中σys为材料屈服强度。
3.如权利要求1所述的一种基于蠕变激活能理论的高温蠕变寿命预测解析计算方法,其特征在于,所述高温蠕变寿命预测解析计算方法适用作业温度为400-1200℃。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810388386.6A CN110411864B (zh) | 2018-04-26 | 2018-04-26 | 一种基于蠕变激活能的高温蠕变寿命预测解析计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810388386.6A CN110411864B (zh) | 2018-04-26 | 2018-04-26 | 一种基于蠕变激活能的高温蠕变寿命预测解析计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110411864A true CN110411864A (zh) | 2019-11-05 |
CN110411864B CN110411864B (zh) | 2022-02-25 |
Family
ID=68346019
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810388386.6A Active CN110411864B (zh) | 2018-04-26 | 2018-04-26 | 一种基于蠕变激活能的高温蠕变寿命预测解析计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110411864B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113252465A (zh) * | 2021-05-20 | 2021-08-13 | 天津理工大学 | 一种基于m-h法的耐热钢蠕变寿命预测方法 |
CN113866008A (zh) * | 2021-09-07 | 2021-12-31 | 南京航空航天大学 | 一种基于阈值应力和抗拉强度的蠕变寿命预测方法 |
CN114088517A (zh) * | 2021-09-24 | 2022-02-25 | 核工业理化工程研究院 | 一种材料蠕变寿命试验的加速条件的评价方法 |
CN114295491A (zh) * | 2021-12-14 | 2022-04-08 | 南京航空航天大学 | 一种蠕变损伤及变形随时间演化行为的预测方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5682042A (en) * | 1991-06-28 | 1997-10-28 | International Business Machines Corporation | Nonbolometric superconductive photoresponsive |
CN101710053A (zh) * | 2009-11-06 | 2010-05-19 | 上海师范大学 | 高温材料的蠕变寿命预测方法 |
CN104165811A (zh) * | 2014-08-01 | 2014-11-26 | 天津大学 | 一种环形焊接试样蠕变试验装置 |
CN105158084A (zh) * | 2015-09-15 | 2015-12-16 | 华东理工大学 | 一种材料的蠕变-疲劳寿命预测方法 |
CN105158080A (zh) * | 2015-05-12 | 2015-12-16 | 上海发电设备成套设计研究院 | 一种预测高温材料蠕变寿命的加速试验法 |
CN106446390A (zh) * | 2016-09-19 | 2017-02-22 | 核工业理化工程研究院 | 金属材料稳态蠕变速率拟合方程的计算方法 |
CN106557630A (zh) * | 2016-11-21 | 2017-04-05 | 中国石油大学(华东) | 一种材料在多轴应力状态下的蠕变‑损伤寿命预测方法 |
CN107895088A (zh) * | 2017-11-30 | 2018-04-10 | 天津大学 | 一种航空发动机燃烧室寿命预测方法 |
-
2018
- 2018-04-26 CN CN201810388386.6A patent/CN110411864B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5682042A (en) * | 1991-06-28 | 1997-10-28 | International Business Machines Corporation | Nonbolometric superconductive photoresponsive |
CN101710053A (zh) * | 2009-11-06 | 2010-05-19 | 上海师范大学 | 高温材料的蠕变寿命预测方法 |
CN104165811A (zh) * | 2014-08-01 | 2014-11-26 | 天津大学 | 一种环形焊接试样蠕变试验装置 |
CN105158080A (zh) * | 2015-05-12 | 2015-12-16 | 上海发电设备成套设计研究院 | 一种预测高温材料蠕变寿命的加速试验法 |
CN105158084A (zh) * | 2015-09-15 | 2015-12-16 | 华东理工大学 | 一种材料的蠕变-疲劳寿命预测方法 |
CN106446390A (zh) * | 2016-09-19 | 2017-02-22 | 核工业理化工程研究院 | 金属材料稳态蠕变速率拟合方程的计算方法 |
CN106557630A (zh) * | 2016-11-21 | 2017-04-05 | 中国石油大学(华东) | 一种材料在多轴应力状态下的蠕变‑损伤寿命预测方法 |
CN107895088A (zh) * | 2017-11-30 | 2018-04-10 | 天津大学 | 一种航空发动机燃烧室寿命预测方法 |
Non-Patent Citations (5)
Title |
---|
DONGQUAN WU: "Theoretical and numerical analysis of creep crack initiation combined with primary and secondary stresses", 《THEORETICAL AND APPLIED FRACTURE MECHANICS》 * |
M.T. WHITTAKER 等: "Creep and creep fracture of 2.25Cr-1.6W steels", 《MATERIALS SCIENCE AND ENGINEERING A》 * |
刘贤翠 等: "3003铝合金蠕变行为与本构方程", 《材料导报 B:研究篇》 * |
张国尚 等: "80Au-20Sn钎料焊点可靠性研究现状与展望", 《机械工程材料》 * |
王磊 等: "《材料强韧学基础》", 30 September 2012, 上海交通大学出版社 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113252465A (zh) * | 2021-05-20 | 2021-08-13 | 天津理工大学 | 一种基于m-h法的耐热钢蠕变寿命预测方法 |
CN113866008A (zh) * | 2021-09-07 | 2021-12-31 | 南京航空航天大学 | 一种基于阈值应力和抗拉强度的蠕变寿命预测方法 |
CN114088517A (zh) * | 2021-09-24 | 2022-02-25 | 核工业理化工程研究院 | 一种材料蠕变寿命试验的加速条件的评价方法 |
CN114295491A (zh) * | 2021-12-14 | 2022-04-08 | 南京航空航天大学 | 一种蠕变损伤及变形随时间演化行为的预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110411864B (zh) | 2022-02-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110411864A (zh) | 一种基于蠕变激活能的高温蠕变寿命预测解析计算方法 | |
CN105117536B (zh) | 一种rpv含裂纹类缺陷的简化弹塑性断裂力学分析方法 | |
Arioka | 2014 WR Whitney Award Lecture: Change in Bonding Strength at Grain Boundaries Before Long-Term SCC Initiation | |
CN105527311B (zh) | 大体积混凝土细观破裂试验系统 | |
CN104458400A (zh) | 一种高温高压水汽环境结构材料试验装置 | |
Sun et al. | Constitutive relationship of IN690 superalloy by using uniaxial compression tests | |
Kim et al. | A direct assessment of creep life based on small punch creep test | |
CN110907475A (zh) | 一种马氏体耐热钢的剩余寿命评估方法 | |
Wang et al. | Temperature induced creep in low-alloy structural Q345 steel | |
Zheng et al. | Creep life prediction of small punch creep testing specimens for service-exposed Cr5Mo using the theta-projection method | |
CN110411863B (zh) | 一种基于蠕变延性的高温蠕变寿命预测方法 | |
Sawada et al. | Analysis of long-term creep curves by constitutive equations | |
Perez-Feró et al. | Experimental database of E110 claddings exposed to accident conditions | |
Tu et al. | Effect of stress regime‐dependent creep behaviour on measurement of creep strain rate based on small specimen techniques | |
CN103278525A (zh) | 一种承压设备火灾后的安全评估方法 | |
Yan et al. | Creep-fatigue tests on full scale directionally solidified turbine blades | |
CN112730112B (zh) | 适用于反应堆结构部件材料长期服役后的环境断裂评估方法 | |
Tu et al. | Effect of large deformation on creep property evaluation by small-specimen bending tests | |
Brett et al. | Impression creep workshop agreement activity | |
Han et al. | Research on creep behavior of three-point bending specimen with fixed constraints at the large deformation stage | |
PODRO et al. | Small-punch testing of a weld’s heat-affected zones | |
Zhang et al. | Assessment of prior fatigue damage and a new approach to predict remanent creep rupture of P92 steel | |
Tan et al. | Environmentally Assisted Fatigue Evaluation Model of Alloy 690 Steam Generator Tube in High-Temperature Water | |
Sanrutsadakorn et al. | Constitutive modeling of flow behaviour of AISI 4340 steel under hot working conditions | |
CN114894647B (zh) | 一种合金预生膜弹性模量测试方法及其应用 |
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 | ||
CP02 | Change in the address of a patent holder | ||
CP02 | Change in the address of a patent holder |
Address after: 300452 Binhai Industrial Research Institute Campus of Tianjin University, No. 48 Jialingjiang Road, Binhai New Area, Tianjin Patentee after: Tianjin University Address before: 300072 Tianjin City, Nankai District Wei Jin Road No. 92 Patentee before: Tianjin University |