CN108256179A - 一种预测材料蠕变曲线的方法 - Google Patents
一种预测材料蠕变曲线的方法 Download PDFInfo
- Publication number
- CN108256179A CN108256179A CN201711483312.2A CN201711483312A CN108256179A CN 108256179 A CN108256179 A CN 108256179A CN 201711483312 A CN201711483312 A CN 201711483312A CN 108256179 A CN108256179 A CN 108256179A
- Authority
- CN
- China
- Prior art keywords
- creep
- curve
- creep curve
- parameter
- prediction
- 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
- 238000000034 method Methods 0.000 title claims abstract description 67
- 239000000463 material Substances 0.000 title claims abstract description 21
- 239000007769 metal material Substances 0.000 claims abstract description 15
- 238000010998 test method Methods 0.000 claims abstract description 15
- 238000012360 testing method Methods 0.000 claims abstract description 7
- 229910045601 alloy Inorganic materials 0.000 claims description 15
- 239000000956 alloy Substances 0.000 claims description 15
- 239000002131 composite material Substances 0.000 claims description 7
- 229910052751 metal Inorganic materials 0.000 claims description 5
- 239000002184 metal Substances 0.000 claims description 5
- 150000002739 metals Chemical class 0.000 claims description 3
- 239000000203 mixture Substances 0.000 abstract description 2
- 229910000831 Steel Inorganic materials 0.000 description 20
- 239000010959 steel Substances 0.000 description 20
- FYYHWMGAXLPEAU-UHFFFAOYSA-N Magnesium Chemical compound [Mg] FYYHWMGAXLPEAU-UHFFFAOYSA-N 0.000 description 5
- RTAQQCXQSZGOHL-UHFFFAOYSA-N Titanium Chemical compound [Ti] RTAQQCXQSZGOHL-UHFFFAOYSA-N 0.000 description 5
- 239000004411 aluminium Substances 0.000 description 5
- 229910052782 aluminium Inorganic materials 0.000 description 5
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 5
- 229910052749 magnesium Inorganic materials 0.000 description 5
- 239000011777 magnesium Substances 0.000 description 5
- 239000010936 titanium Substances 0.000 description 5
- 229910052719 titanium Inorganic materials 0.000 description 5
- PXHVJJICTQNCMI-UHFFFAOYSA-N Nickel Chemical compound [Ni] PXHVJJICTQNCMI-UHFFFAOYSA-N 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 4
- 238000013213 extrapolation Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 239000003779 heat-resistant material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000005088 metallography Methods 0.000 description 1
- 229910052759 nickel Inorganic materials 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000005482 strain hardening Methods 0.000 description 1
- 229910000601 superalloy Inorganic materials 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
一种预测材料蠕变曲线的方法，该方法采用如下两种方法之一：第一种方法包括以下步骤：根据《GB/T2039‑2012金属材料单轴拉伸蠕变试验方法》，在拟预测温度及应力下测定拉伸蠕变曲线，第二种方法：本方法同样适用于通过测试一定应力、不同温度下的拉伸蠕变曲线预测相同应力其它温度下的拉伸蠕变曲线，本发明揭示了蠕变曲线模型中，主要描述蠕变第一阶段组成项参数、在蠕变进入稳态阶段以后基本不变的事实，采用多点加权方式确定参数，进一步提高蠕变曲线的预测精度。
Description
技术领域
本发明属于材料高温力学性能研究领域，涉及一种材料蠕变曲线的预测方法。
背景技术
蠕变曲线反映材料的蠕变性能特征和蠕变失效过程，准确预测材料的蠕变曲线对于研究材料的蠕变行为、蠕变机制、预测材料蠕变寿命，加快研究开发高性能耐热材料步伐意义重大。多年来，许多学者对表征材料蠕变曲线的模型进行了研究，并努力探索蠕变曲线的预测方法。
上世纪八十年代，Evans RW和Wilshire B等人基于时间硬化和应变硬化原理提出了描述蠕变曲线的θ影射法，其模型为：
式中：ε为应变，t为时间，θj(j＝1，2，3，4)为参数，可通过对相关实测数据分析获得。表达了蠕变硬化过程，主要描述蠕变的第一阶段；表达了蠕变软化过程，主要描述蠕变的第三阶段。用θ影射法进行蠕变曲线预测时，可采用相同温度、不同应力或相同应力、不同温度蠕变曲线拟合θ方程参数外推，并通过以下公式实现。
相同温度、不同应力θ参数经验公式：
lgθj＝aj+bjσ(j＝1，2，3，4) (2)
相同应力、不同温度θ参数经验公式：
lgθj＝aj+bjT(j＝1，2，3，4) (3)
运用θ方程能较好地反映蠕变第一阶段和第二阶段，而与蠕变第三阶段偏差很大。日本Maruyama K和Oikawa H等人对θ方程进行了改进，提出修正θ方程：
修正θ方程与θ方程相比，参数少，能相对精确地反映蠕变第三阶段，但蠕变第一阶段误差较大，以至对一些蠕变曲线难以表达。因此，应用θ方程和修正θ方程预测蠕变曲线以及预测蠕变寿命很难达到理想效果。日本学者伊津野仁史等人针对单晶镍基合金，提出复函数形式，包含9个参数在内，由三项式组成的蠕变曲线模型，见“Ni基超合金のためのクリープ构成式の因子分析”.日本金属学会志第71卷第2号(2007)223-225。该模型结构复杂，参数多，预测蠕变曲线难度大，且需要大量的实验。申请号为201510239105.7的专利“一种预测高温材料蠕变寿命的加速试验法”提出：通过对快速蠕变试验数据分析处理，外推高温材料蠕变寿命的方法，该方法简单快捷，但操作中主要靠人工判断蠕变第二、第三阶段起始点，精度有限，同时，无法了解蠕变过程。
近年来，人们又提出一些表达蠕变曲线的新模型，如幂指数模型式中：ε为应变，t为时间，δj(j＝1，2，3，4)为待定参数。该模型能很好地描述各种材料的蠕变曲线，但目前尚无确定参数δ1、δ2的方法，而拟合蠕变曲线方程参数对试验条件敏感，简单利用参数经验公式外推基本误差无法控制，难以保证预测精度。
发明内容：
发明目的：
本发明提供一种预测材料蠕变曲线的方法，其目的是解决以往所存在的问题，本发明基于材料蠕变曲线的一般特征以及模型的特点，经过对多种材料不同条件蠕变曲线模拟试验研究，提出采用实测与外推相结合确定预测蠕变曲线方程参数的方法进行蠕变曲线预测。上述模型中：ε为应变，t为时间，δj(j＝1，2，3，4)为参数。实践证明，本发明提出的方法简捷，精度高。在试验条件下，对403Nb钢600℃、260MPa及260MPa、580℃蠕变曲线预测精度分别为94.4％和97.6％。
此外，进一步明确本方法还可用于通过其它二项式模型对各种耐热钢、耐热合金、高温合金、钛基、镁基、铝基等合金、部分陶瓷以及一些复合材料蠕变曲线包括蠕变寿命的预测。其中：ε为应变，t为时间，α系列及β1、β2代表蠕变参数
技术方案：
一种预测材料蠕变曲线的方法，其特征在于：该方法采用如下两种方法之一：
第一种方法包括以下步骤：
(1)、根据《GB/T2039-2012金属材料单轴拉伸蠕变试验方法》，在拟预测温度Tm及应力σm下测定拉伸蠕变曲线，要求蠕变进入稳态阶段；在蠕变曲线的稳态蠕变区间选取3个以上时间点，通过蠕变曲线模型对上述不同时间的蠕变数据进行拟合，求出各拟合蠕变曲线方程参数δ1，δ2的平均值δm1，δm2；式中：ε为应变，t为时间，δj(j＝1，2，3， 4)为待定参数；
(2)根据《GB/T2039-2012金属材料单轴拉伸蠕变试验方法》，在Tm温度下测定3个以上不同应力σi(i＝1，2，3，……)的拉伸蠕变曲线，直至样品断裂；通过模型对不同应力的蠕变数据进行拟合，确定各拟合蠕变曲线方程参数δ3，δ4的值δij(i＝1，2，3，……；j＝3，4)；
(3)通过公式lgδ＝a+bσ分别对由步骤(2)得到不同应力σi(i＝1，2，3，……)的拟合蠕变曲线方程参数值δij(i＝1，2，3，……；j＝3，4)进行拟合，确定常数a，b的值am3，bm3及am4，bm4；得到参数δ3，δ4与应力σ的关系式：
lgδ3＝am3+bm3σ (1)
lgδ4＝am4+bm4σ (2)
将应力σm代入公式(1)、(2)，求出δ3及δ4的估值δ'm3，δ'm4；
(4)将由步骤(1)得到的参数值δm1，δm2及由步骤(3)得到的参数值δ'm3，δ'm4代入蠕变曲线模型得到温度Tm应力σm下的预测蠕变曲线方程：
根据预测蠕变曲线方程(3)绘制预测蠕变曲线；
第二种方法：本方法同样适用于通过测试一定应力、不同温度下的拉伸蠕变曲线预测相同应力其它温度下的拉伸蠕变曲线，包括以下步骤：
(1)根据《GB/T2039-2012金属材料单轴拉伸蠕变试验方法》，在拟预测应力σr及温度Tr下测定拉伸蠕变曲线，要求蠕变进入稳态阶段；在蠕变曲线的稳态蠕变区间选取3个以上时间点，通过蠕变曲线模型对上述不同时间的蠕变数据进行拟合，求出各拟合蠕变曲线方程参数δ1，δ2的平均值δr1，δr2；式中：ε为应变，t为时间，δj(j＝1，2，3，4)为 待定参数；
(2)根据《GB/T2039-2012金属材料单轴拉伸蠕变试验方法》，在拟预测应力σr下，测定3个以上不同温度Tk(k＝1，2，3，……)的拉伸蠕变曲线，直至样品断裂；通过模型对不同温度的蠕变数据进行拟合，确定各拟合蠕变曲线方程参数δ3，δ4的值δkj(k＝1，2，3，……；j＝3，4)；
(3)通过公式lgδ＝a+bT分别对由步骤(2)得到不同温度Tk(k＝1，2，3，……)的拟合蠕变曲线方程参数值δkj(k＝1，2，3，……；j＝3，4)进行拟合，确定常数a，b的值ar3，br3及ar4，br4；得到参数δ3，δ4与温度T的关系式：
lgδ3＝ar3+br3T (4)
lgδ4＝ar4+br4T (5)
将温度Tr代入公式(4)、(5)，求出δ3及δ4的值δ'r3，δ'r4；
(4)将由步骤(1)得到的参数值δr1，δr2及由步骤(3)得到的参数值δ'r3，δ'r4代入蠕变曲线模型得到应力σr温度Tr下的预测蠕变曲线方程：
根据预测蠕变曲线方程(6)绘制预测蠕变曲线。
本方法适合用于通过其它二项式模型进行蠕变曲线以及蠕变寿命的预测，其中：fα(t)主要描述前期的蠕变过程，ε为应变，t为时间，α系列及β1、β2代表蠕变参数；
本方法应用于各种金属、金属基合金及其复合材料蠕变曲线包括蠕变寿命的预测。
试验证明，模型几乎能完整地描述耐热钢、耐热合金、镍基、铜基、铝基、镁基、钛基等合金及一些金属基复合材料蠕变的全过程。式中：ε为应变，t为蠕变时间，δj(j＝1，2，3，4)为待定参数。表征蠕变硬化过程，主要描述蠕变第一阶段，当蠕变进入稳态阶段以后，δ1与δ2的值比较稳定；表征蠕变软化过程，主要描述蠕变后期的变形行为，δ3与δ4随着蠕变时间的延长而变化。因此，在预测蠕变曲线方程时，δ1和与δ2的值可通过拟合进入稳态蠕变阶段较短时间的待预测蠕变曲线直接获得。而δ3和δ4的值则可通过相同温度、不同应力蠕变参数经验公式lgδ＝a+bσ或相同应力、不同温度蠕变参数经验公式lgδ＝a+bT外推得到，上述式中：a、b为常数。于是，预测蠕变曲线方程中待定参数由原来的4个减少到2个，简化了试验和数据分析过程，并且显著降低预测参数的分散性，从而使蠕变曲线的预测精度将大幅度提高。
本发明基于以上原理，先通过蠕变曲线模型对拟预测温度和应力下进入稳态的短时间蠕变曲线进行拟合分析，确定参数δ1和δ2的值，再根据模型对预测温度、不同应力或预测应力、不同温度下的蠕变曲线进行拟合，借助参数经验公式lgδ＝a+bσ或lgδ＝a+bT外推，确定预测蠕变曲线方程参数δ3和δ4的值，由此得到预测蠕变曲线方程，根据预测蠕变曲线方程获得预测蠕变曲线。
由以上分析可以发现，本方法还可用于通过其它二项式模型实现对各种金属、耐热钢、耐热合金、高温合金、钛基、镁基、铝基等合金及金属基复合材料蠕变曲线包括蠕变寿命的预测。其中：fα(t)主要描述前期的蠕变过程，ε为应变，t为时间，α系列及β1、β2代表蠕变参数。
本发明的有益效果：
(1)本发明揭示了蠕变曲线模型中，主要描述蠕变第一阶段组成项参数δ1、δ2在蠕变进入稳态阶段以后基本不变的事实，据此，通过对进入稳态蠕变阶段的短时间预测蠕变曲线拟合分析直接获得δ1与δ2的值，使预测蠕变曲线方程中待定参数由原来的4个减少到2个，从根本上降低了蠕变曲线的预测误差，提高了预测效率。上述模型中：ε为应变，t为蠕变时间，δj(j＝1，2，3，4)为参数。
(2)采用多点加权方式确定参数δ1、δ2，进一步提高蠕变曲线的预测精度。
(3)本方法还可用于通过其它二项式模型实现对各种耐热钢、耐热合金、高温合金、钛基、镁基、铝基等合金及金属基复合材料蠕变曲线包括蠕变寿命的预测。其中：fα(t)主要描述前期的蠕变过程，ε为应变，t为时间，α系列及β1、β2代表蠕变参数。因此，具有一定的理论意义及广泛的应用价值。
附图说明：
图1是实施例1中403Nb钢600℃，260MPa、270MPa、280MPa下实测蠕变曲线及其拟合蠕变曲线。
图2是实施例1中403Nb钢600℃、240MPa预测蠕变曲线与实测蠕变曲线的对比。
图3是实施例2中403Nb钢260MPa，600℃、610℃、620℃下实测蠕变曲线及其拟合蠕变曲线。
图4是实施例2中403Nb钢260MPa、580℃预测蠕变曲线与实测蠕变曲线的对比。
具体实施方式：
一种预测材料蠕变曲线的方法实施方式分两种情况：
第一种情况：根据一定温度，不同应力的蠕变曲线预测相同温度其它应力的蠕变曲线，步骤如下：
(1)根据《GB/T2039-2012金属材料单轴拉伸蠕变试验方法》，在拟预测温度Tm及应力σm下测定拉伸蠕变曲线，要求蠕变进入稳态阶段；在蠕变曲线的稳态蠕变区间选取3个以上时间点，通过蠕变曲线模型对上述不同时间的蠕变数据进行拟合，求出各拟合蠕变曲线方程参数δ1，δ2的平均值δm1，δm2；式中：ε为应变，t为时间，δj(j＝1，2，3，4)为待定参数；
(2)根据《GB/T2039-2012金属材料单轴拉伸蠕变试验方法》，在Tm温度下测定3个以上不同应力σi(i＝1，2，3，……)的拉伸蠕变曲线，直至样品断裂；通过蠕变曲线模型对不同应力的蠕变数据进行拟合，确定各拟合蠕变曲线方程参数δ3，δ4的值δij(i＝1，2，3，……；j＝3，4)；
(3)通过公式lgδ＝a+bσ分别对由步骤(2)得到不同应力σi(i＝1，2，3，……)的拟合蠕变曲线方程参数值δij(i＝1，2，3，……；j＝3，4)进行拟合，确定常数a，b的值am3，bm3及am4，bm4，得到参数δ3，δ4与应力σ的关系式：
lgδ3＝am3+bm3σ (1)
lgδ4＝am4+bm4σ (2)
将应力σm代入公式(1)、(2)，求出δ3及δ4的值δ'm3，δ'm4；
(4)将由步骤(1)得到的参数值δm1，δm2及由步骤(3)得到的参数值δ'm3，δ'm4代入模型得到温度Tm应力σm下的预测蠕变曲线方程：
根据预测蠕变曲线方程(3)绘制预测蠕变曲线；
第二种情况：根据一定应力、不同温度的蠕变曲线预测相同应力、其它温度的蠕变曲线，步骤如下：
(1)根据《GB/T2039-2012金属材料单轴拉伸蠕变试验方法》，在拟预测应力σr及温度Tr下测定拉伸蠕变曲线，要求蠕变进入稳态阶段；在蠕变曲线的稳态蠕变区间选取3个以上时间点，通过蠕变曲线模型对上述不同时间的蠕变数据进行拟合，求出各拟合蠕变曲线方程参数δ1，δ2的平均值δr1，δr2；式中：ε为应变，t为时间，δj(j＝1，2，3，4)为待定参数；
(2)根据《GB/T2039-2012金属材料单轴拉伸蠕变试验方法》，在拟预测应力σr下，测定3个以上不同温度Tk(k＝1，2，3，……)的拉伸蠕变曲线，直至样品断裂；通过模型对不同温度的蠕变数据进行拟合，确定各拟合蠕变曲线方程参数δ3，δ4的值δkj(k＝1，2，3，……；j＝3，4)；
(3)通过公式lgδ＝a+bT分别对由步骤(2)得到不同温度Tk(k＝1，2，3，……)的拟合蠕变曲线方程参数值δkj(k＝1，2，3，……；j＝3，4)进行拟合，确定常数a、b的值ar3，br3及ar4，br4；得到参数δ3，δ4与温度T的关系式：
lgδ3＝ar3+br3T (4)
lgδ4＝ar4+br4T (5)
将温度Tr代入公式(4)、(5)，求出δ3及δ4的值δ'r3，δ'r4；
(4)将由步骤(1)得到的参数值δr1，δr2及由步骤(3)得到的参数值δ'r3，δ'r4代入模型得到应力σr温度Tr下的预测蠕变曲线方程：
根据预测蠕变曲线方程(6)绘制预测蠕变曲线；
本方法可同样用于通过其它二项式模型对各种耐热钢、耐热合金、高温合金、钛基、镁基、铝基等合金及金属基复合材料蠕变曲线包括蠕变寿命的预测。其中：fα(t)主要描述前期的蠕变过程，ε为应变，t为时间，α系列及β1、β2代表蠕变参数。
实施例：
实施例1403Nb钢600℃、240MPa下拉伸蠕变曲线的预测
步骤如下：
(1)根据GB/T2039-2012金属材料单轴拉伸蠕变试验方法》，采用截面2.5mm×4.5mm，标距14mm的片状“工”形试样，测定403Nb钢600℃、240MPa下的拉伸蠕变曲线，当时间达到25h时，蠕变进入稳态阶段。通过蠕变曲线模型采用最小二乘法分别对27h，35h，45h及60h的蠕变数据进行拟合，确定各蠕变曲线拟合方程参数δ1，δ2的值，求出其平均值分别为：0.6354和0.2496，如表1所示。可见，该值与实际245h拉断时的数值非常接近。
表1403Nb钢600℃/240MPa不同时间蠕变曲线方程参数δ1及δ2的值
| t/h | δ1 | δ2 |
| 27 | 0.6253 | 0.2496 |
| 35 | 0.6348 | 0.2495 |
| 45 | 0.6375 | 0.2495 |
| 60 | 0.6439 | 0.2496 |
| 平均 | 0.6354 | 0.2496 |
| 实际245 | 0.6346 | 0.2496 |
(2)按照步骤(1)测定蠕变曲线的方法测定403Nb钢600℃，260MPa、270MPa、280MPa下的拉伸蠕变曲线，直至样品拉断。通过模型对上述蠕变数据进行拟合，如图1所示，确定各应力拟合蠕变曲线方程参数δσj(σ＝260，270，280；j＝3，4)的值，如表2所示；通过公式lgδj＝aj+bjσ对δσj(σ＝260，270，280；j＝3，4)值进行拟合，确定常数a3，b3及a4，b4的值分别为-9.3512，0.0185及-3.0655，0.0073。得到参数计算公式：
lgδ3＝-9.3512+0.0185σ (1)
lgδ4＝-3.0655+0.0073σ (2)
根据公式(1)、(2)求出600℃、240MPa预测蠕变曲线方程参数δ3及δ4的值分别为1.223×10-5，0.0486。
表2 403Nb钢600℃不同应力蠕变曲线拟合方程参数δ3及δ4的值
| σ/MPa | 260 | 270 | 280 |
| δσ3 | 3.8828×10-5 | 2.4039×10-5 | 9.1010×10-5 |
| δσ4 | 0.0667 | 0.0890 | 0.0935 |
(3)将由步骤(1)得到的实测蠕变曲线拟合方程参δ1，δ2的值0.6354，0.2496及由步骤(2)求出预测蠕变曲线方程参数δ3，δ4的值1.223×10-5，0.0486代入蠕变曲线模型，得到600℃、240MPa预测蠕变曲线方程：
ε＝0.6354t 0.2496+1.223×10-5(e0.0486t-1) (3)
(4)根据方程(3)绘制预测蠕变曲线。图2为600℃、240MPa下预测蠕变曲线与实测蠕变曲线的对比。其预测蠕变寿命262h与实测蠕变寿命248h相比精度为94.4％。
实施例2 403Nb钢260MPa、580℃下拉伸蠕变曲线的预测
步骤如下：
(1)根据GB/T2039-2012金属材料单轴拉伸蠕变试验方法》采用截面2.5mm×4.5mm，标距14mm的片状“工”形试样，测定403Nb钢260MPa、580℃下的拉伸蠕变曲线，当时间达到38h时，蠕变进入稳态阶段。通过蠕变曲线模型采用最小二乘法对40h、55h、70h蠕变数据进行拟合，确定各蠕变曲线拟合方程参数δ1，δ2的值，求出其平均值分别为：0.4653、0.2203。如表3所示。
表3 403Nb钢260MPa/580℃不同时间蠕变曲线方程参数δ1，δ2的值
| t/h | δ1 | δ2 |
| 40 | 0.4644 | 0.22 |
| 55 | 0.4657 | 0.2201 |
| 70 | 0.4658 | 0.221 |
| 平均 | 0.4653 | 0.2203 |
(2)按照步骤(1)测定蠕变曲线的方法测定403Nb钢260MPa，600℃、610℃、620℃下的拉伸蠕变曲线，直至样品拉断。通过模型对上述蠕变数据进行拟合，如图3所示，确定各温度拟合蠕变曲线方程参数δTj(T＝600，610，620；j＝3，4)的值，如表4所示；通过公式lgδj＝aj+bjT对δTj(T＝600，610，620；j＝3，4)的值进行拟合，确定常数a3，b3及a4，b4的值分别为80.026，-0.0959及-24.1671，0.026。得到参数计算公式：
lgδ3＝80.026-0.0959T (4)
lgδ4＝-24.0471+0.026T (5)
根据公式(4)、(5)求出预测260MPa、580℃蠕变曲线方程参数δ3及δ4的值分别为0.021，0.01074。
表4 403Nb钢260MPa不同温度蠕变曲线拟合方程参数δ3及δ4的值
| T/℃ | 600 | 610 | 620 |
| δT3 | 0.0011 | 5.0797×10-5 | 1.3343×10-5 |
| δT4 | 0.0475 | 0.0886 | 0.1577 |
(3)将由步骤(1)得到的实测蠕变曲线拟合方程参数δ1，δ2的值0.4653，0.2203及由步骤(2)求出预测蠕变曲线方程参数δ3，δ4的值0.021、0.01074代入蠕变曲线模型，得到预测260MPa、580℃蠕变曲线方程：
ε＝0.4653t0.2203+0.021(e0.01074-1) (6)
(4)根据预测蠕变曲线方程(6)绘制预测蠕变曲线。图4为260MPa、580℃预测蠕变曲线与实测蠕变曲线的对比。其预测蠕变寿命559h与实测蠕变寿命571h相比，预测精度为97.6％。
Claims (3)
1.一种预测材料蠕变曲线的方法，其特征在于：该方法采用如下两种方法之一：
第一种方法包括以下步骤：
(1)、根据《GB/T2039-2012金属材料单轴拉伸蠕变试验方法》，在拟预测温度Tm及应力σm下测定拉伸蠕变曲线，要求蠕变进入稳态阶段；在蠕变曲线的稳态蠕变区间选取3个以上时间点，通过蠕变曲线模型对上述不同时间的蠕变数据进行拟合，求出各拟合蠕变曲线方程参数δ1，δ2的平均值δm1，δm2；式中：ε为应变，t为时间，δj(j＝1，2，3，4)为待 定参数；
(2)根据《GB/T2039-2012金属材料单轴拉伸蠕变试验方法》，在Tm温度下测定3个以上不同应力σi(i＝1，2，3，……)的拉伸蠕变曲线，直至样品断裂；通过模型对不同应力的蠕变数据进行拟合，确定各拟合蠕变曲线方程参数δ3，δ4的值δij(i＝1，2，3，……；j＝3，4)；
(3)通过公式lgδ＝a+bσ分别对由步骤(2)得到不同应力σi(i＝1，2，3，……)的拟合蠕变曲线方程参数值δij(i＝1，2，3，……；j＝3，4)进行拟合，确定常数a，b的值am3，bm3及am4，bm4；得到参数δ3，δ4与应力σ的关系式：
lgδ3＝am3+bm3σ (1)
lgδ4＝am4+bm4σ (2)
将应力σm代入公式(1)、(2)，求出δ3及δ4的估值δ'm3，δ'm4；
(4)将由步骤(1)得到的参数值δm1，δm2及由步骤(3)得到的参数值δ'm3，δ'm4代入蠕变曲线模型得到温度Tm应力σm下的预测蠕变曲线方程：
根据预测蠕变曲线方程(3)绘制预测蠕变曲线；
第二种方法：本方法同样适用于通过测试一定应力、不同温度下的拉伸蠕变曲线预测相同应力其它温度下的拉伸蠕变曲线，包括以下步骤：
(1)根据《GB/T2039-2012金属材料单轴拉伸蠕变试验方法》，在拟预测应力σr及温度Tr下测定拉伸蠕变曲线，要求蠕变进入稳态阶段；在蠕变曲线的稳态蠕变区间选取3个以上时间点，通过蠕变曲线模型对上述不同时间的蠕变数据进行拟合，求出各拟合蠕变曲线方程参数δ1，δ2的平均值δr1，δr2；式中：ε为应变，t为时间，δj(j＝1，2，3，4)为待 定参数；
(2)根据《GB/T2039-2012金属材料单轴拉伸蠕变试验方法》，在拟预测应力σr下，测定3个以上不同温度Tk(k＝1，2，3，……)的拉伸蠕变曲线，直至样品断裂；通过模型对不同温度的蠕变数据进行拟合，确定各拟合蠕变曲线方程参数δ3，δ4的值δkj(k＝1，2，3，……；j＝3，4)；
(3)通过公式lgδ＝a+bT分别对由步骤(2)得到不同温度Tk(k＝1，2，3，……)的拟合蠕变曲线方程参数值δkj(k＝1，2，3，……；j＝3，4)进行拟合，确定常数a，b的值ar3，br3及ar4，br4；得到参数δ3，δ4与温度T的关系式：
lgδ3＝ar3+br3T (4)
lgδ4＝ar4+br4T (5)
将温度Tr代入公式(4)、(5)，求出δ3及δ4的值δ'r3，δ'r4；
(4)将由步骤(1)得到的参数值δr1，δr2及由步骤(3)得到的参数值δ'r3，δ'r4代入蠕变曲线模型得到应力σr温度Tr下的预测蠕变曲线方程：
根据预测蠕变曲线方程(6)绘制预测蠕变曲线。
2.根据权利要求1所述的一种预测材料蠕变曲线的方法，其特征在于：本方法适合用于通过其它二项式模型进行蠕变曲线以及蠕变寿命的预测，其中：fα(t)主要描述前期的蠕变过程，ε为应变，t为时间，α系列及β1、β2代表蠕变参数。
3.根据权利要求1或2所述的一种预测材料蠕变曲线的方法，其特征在于：本方法应用于各种金属、金属基合金及其复合材料蠕变曲线包括蠕变寿命的预测。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711483312.2A CN108256179B (zh) | 2017-12-29 | 2017-12-29 | 一种预测材料蠕变曲线的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711483312.2A CN108256179B (zh) | 2017-12-29 | 2017-12-29 | 一种预测材料蠕变曲线的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108256179A true CN108256179A (zh) | 2018-07-06 |
| CN108256179B CN108256179B (zh) | 2021-06-15 |
Family
ID=62725382
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711483312.2A Active CN108256179B (zh) | 2017-12-29 | 2017-12-29 | 一种预测材料蠕变曲线的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108256179B (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109142083A (zh) * | 2018-09-05 | 2019-01-04 | 南京航空航天大学 | 一种变载荷历程下的蠕变损伤计算方法及模型 |
| CN109871649A (zh) * | 2019-03-11 | 2019-06-11 | 中国石油大学(华东) | 一种Norton-Bailey模型参数标定方法及系统 |
| CN110940572A (zh) * | 2019-12-11 | 2020-03-31 | 北京科技大学 | 一种高Cr铁素体耐热钢的蠕变寿命预测方法 |
| CN113008677A (zh) * | 2021-03-02 | 2021-06-22 | 中南大学 | 一种镍基高温合金的蠕变持久性能预测方法 |
| CN113125275A (zh) * | 2021-04-06 | 2021-07-16 | 西北工业大学 | 镍基单晶高温合金蠕变模型参数确定与蠕变寿命预测方法 |
| CN113118458A (zh) * | 2021-04-20 | 2021-07-16 | 江西省科学院应用物理研究所 | 一种激光选区熔化成形金属构件拉伸性能的预测方法 |
| CN114088517A (zh) * | 2021-09-24 | 2022-02-25 | 核工业理化工程研究院 | 一种材料蠕变寿命试验的加速条件的评价方法 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101308127A (zh) * | 2008-07-08 | 2008-11-19 | 华东理工大学 | 一种材料的高温蠕变寿命加速试验方法 |
| CN104597757A (zh) * | 2014-12-31 | 2015-05-06 | 北京理工大学 | 一种波纹管驱动定位平台蠕变建模与补偿控制方法 |
| CN105004617A (zh) * | 2015-07-20 | 2015-10-28 | 沈阳工业大学 | 一种描述金属材料蠕变曲线的方法 |
| CN105910921A (zh) * | 2016-04-11 | 2016-08-31 | 沈阳工业大学 | 一种预测dz125合金蠕变曲线的方法 |
| JP2017049066A (ja) * | 2015-08-31 | 2017-03-09 | 新日鐵住金株式会社 | 耐熱鋼材の寿命予測方法 |
| WO2017051739A1 (ja) * | 2015-09-25 | 2017-03-30 | 株式会社デンソー | ヘッドアップディスプレイ装置 |
| CN106568655A (zh) * | 2016-10-28 | 2017-04-19 | 沈阳工业大学 | 一种预测耐热合金蠕变寿命的方法 |
| CN107391903A (zh) * | 2017-06-07 | 2017-11-24 | 上海电力学院 | 用于马氏体钢的蠕变‑疲劳寿命预测及其可靠性分析方法 |
-
2017
- 2017-12-29 CN CN201711483312.2A patent/CN108256179B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101308127A (zh) * | 2008-07-08 | 2008-11-19 | 华东理工大学 | 一种材料的高温蠕变寿命加速试验方法 |
| CN104597757A (zh) * | 2014-12-31 | 2015-05-06 | 北京理工大学 | 一种波纹管驱动定位平台蠕变建模与补偿控制方法 |
| CN105004617A (zh) * | 2015-07-20 | 2015-10-28 | 沈阳工业大学 | 一种描述金属材料蠕变曲线的方法 |
| JP2017049066A (ja) * | 2015-08-31 | 2017-03-09 | 新日鐵住金株式会社 | 耐熱鋼材の寿命予測方法 |
| WO2017051739A1 (ja) * | 2015-09-25 | 2017-03-30 | 株式会社デンソー | ヘッドアップディスプレイ装置 |
| CN105910921A (zh) * | 2016-04-11 | 2016-08-31 | 沈阳工业大学 | 一种预测dz125合金蠕变曲线的方法 |
| CN106568655A (zh) * | 2016-10-28 | 2017-04-19 | 沈阳工业大学 | 一种预测耐热合金蠕变寿命的方法 |
| CN107391903A (zh) * | 2017-06-07 | 2017-11-24 | 上海电力学院 | 用于马氏体钢的蠕变‑疲劳寿命预测及其可靠性分析方法 |
Non-Patent Citations (2)
| Title |
|---|
| M.T.WHITTAKER 等: "Long-term creep data prediction for type 316H stainless steel", 《MATERIALS SCIENCE AND ENGINEERING:A》 * |
| 于兴福 等: "无铼二代镍基单晶高温合金中温高应力蠕变机制", 《中国有色金属学报》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109142083A (zh) * | 2018-09-05 | 2019-01-04 | 南京航空航天大学 | 一种变载荷历程下的蠕变损伤计算方法及模型 |
| CN109142083B (zh) * | 2018-09-05 | 2020-04-24 | 南京航空航天大学 | 一种变载荷历程下的蠕变损伤计算方法 |
| CN109871649A (zh) * | 2019-03-11 | 2019-06-11 | 中国石油大学(华东) | 一种Norton-Bailey模型参数标定方法及系统 |
| CN110940572A (zh) * | 2019-12-11 | 2020-03-31 | 北京科技大学 | 一种高Cr铁素体耐热钢的蠕变寿命预测方法 |
| CN113008677A (zh) * | 2021-03-02 | 2021-06-22 | 中南大学 | 一种镍基高温合金的蠕变持久性能预测方法 |
| CN113008677B (zh) * | 2021-03-02 | 2022-06-10 | 中南大学 | 一种镍基高温合金的蠕变持久性能预测方法 |
| CN113125275A (zh) * | 2021-04-06 | 2021-07-16 | 西北工业大学 | 镍基单晶高温合金蠕变模型参数确定与蠕变寿命预测方法 |
| CN113125275B (zh) * | 2021-04-06 | 2024-05-03 | 西北工业大学 | 镍基单晶高温合金蠕变模型参数确定与蠕变寿命预测方法 |
| CN113118458A (zh) * | 2021-04-20 | 2021-07-16 | 江西省科学院应用物理研究所 | 一种激光选区熔化成形金属构件拉伸性能的预测方法 |
| CN113118458B (zh) * | 2021-04-20 | 2023-04-07 | 江西省科学院应用物理研究所 | 一种激光选区熔化成形金属构件拉伸性能的预测方法 |
| CN114088517A (zh) * | 2021-09-24 | 2022-02-25 | 核工业理化工程研究院 | 一种材料蠕变寿命试验的加速条件的评价方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108256179B (zh) | 2021-06-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108256179A (zh) | 一种预测材料蠕变曲线的方法 | |
| Shang et al. | A ductile fracture model considering stress state and Zener–Hollomon parameter for hot deformation of metallic materials | |
| Yuan et al. | Effect of mean stress and ratcheting strain on the low cycle fatigue behavior of a wrought 316LN stainless steel | |
| Lin et al. | High-temperature creep behavior of Al–Cu–Mg alloy | |
| Majzoobi et al. | Determination of the constants of damage models | |
| CN109255202A (zh) | 一种用于机械构件疲劳裂纹萌生寿命的预估方法 | |
| Yao et al. | A new fracture criterion for ductile materials based on a finite element aided testing method | |
| CN1207550C (zh) | 高温部件蠕变寿命的测试方法 | |
| Tang et al. | Behavior and modeling of microstructure evolution during metadynamic recrystallization of a Ni-based superalloy | |
| CN104504185A (zh) | 一种基于开裂准则和有限元优化钛合金开坯锻造工艺方法 | |
| Huo et al. | Multi-axis low-cycle creep/fatigue life prediction of high-pressure turbine blades based on a new critical plane damage parameter | |
| CN110826267B (zh) | 发动机缸盖蠕变疲劳分析方法 | |
| CN108893689A (zh) | Inconel718合金盘锻件均匀化制造方法 | |
| Fu et al. | A modified θ projection model for constant load creep curves-I. Introduction of the model | |
| CN107449894A (zh) | 氮化硅结合碳化硅耐火陶瓷的疲劳寿命预测方法 | |
| Almroth et al. | Out-of-phase thermo-mechanical fatigue crack growth and the effect of the compressive minimum load level on crack closure at notches | |
| Zhu et al. | Creep rupture life prediction of nickel-based superalloys based on data fusion | |
| Lan et al. | Artificial neural network approach for mechanical properties prediction of as-cast A380 aluminum alloy | |
| Celentano | Thermomechanical simulation and experimental validation of wire drawing processes | |
| Wang et al. | Cyclic cohesive zone model damage parameter acquisition for fatigue crack growth considering crack closure effect | |
| Liu et al. | Crack growth analysis and fatigue life estimation in the piston rod of a Kaplan hydro turbine | |
| CN106202937B (zh) | M50钢锻造组织中碳化物尺寸预测方法 | |
| Chen et al. | A nonlinear fatigue damage accumulation model under variable amplitude loading considering the loading sequence effect | |
| Chen et al. | An enhanced lemaitre damage model and verification method for X12 alloy steel in hot forming process | |
| Wu et al. | Probabilistic fatigue modeling of notched components under size effect using modified energy field intensity approach |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240619 Address after: B1-a-F202-11, No. 158 Puhe Road, Shenbei New District, Shenyang City, Liaoning Province, 110000 Patentee after: Liaoning Shenghong Qinghang Technology Co.,Ltd. Country or region after: China Address before: 110870 No. 111 Shenyang West Road, Shenyang economic and Technological Development Zone, Liaoning Patentee before: SHENYANG University OF TECHNOLOGY Country or region before: China |
|
| TR01 | Transfer of patent right |