CN106446390B - 金属材料稳态蠕变速率拟合方程的计算方法 - Google Patents
金属材料稳态蠕变速率拟合方程的计算方法 Download PDFInfo
- Publication number
- CN106446390B CN106446390B CN201610828973.3A CN201610828973A CN106446390B CN 106446390 B CN106446390 B CN 106446390B CN 201610828973 A CN201610828973 A CN 201610828973A CN 106446390 B CN106446390 B CN 106446390B
- Authority
- CN
- China
- Prior art keywords
- creep rates
- creep
- secondary creep
- fit equation
- parameter
- 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.)
- Active
Links
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
本发明公开了一种金属材料稳态蠕变速率拟合方程的计算方法包括以下步骤:(ⅰ)验证材料蠕变持久性能的Monkman‑Grant关系;(ⅱ)将时间温度参数修正为速率温度参数;(ⅲ)选择主曲线P(σ)的形式;(ⅳ)拟合速率温度参数模型中的未知常数项;(ⅴ)建立材料的稳态蠕变速率拟合方程。本发明解决了传统时间温度参数无法计算稳态蠕变速率的问题,在设备构件可靠性设计、优化方面的应用更加广泛;可以拟合较宽实验条件范围的稳态蠕变速率,降低了蠕变机理改变对蠕变速率的影响,提高了材料稳态蠕变速率拟合方程的适用性和准确度;具有通用的研究和应用价值。
Description
技术领域
本发明属于一种金属材料稳态蠕变速率拟合方程的计算方法,具体涉及一种基于速率温度参数模型金属材料稳态蠕变速率拟合方程的计算方法。
背景技术
在能源、航空航天、石油化工等工业领域,金属材料的蠕变性能及寿命评估技术是确保材料构件安全运行又避免资源浪费的重要环节。尤其是当设计条件要求金属材料在更宽范围的实验条件下长时间服役,材料蠕变性能、寿命预测模型的准确性显得十分重要。材料稳态蠕变速率拟合方程的计算对于研究材料的蠕变性能及相关构件的寿命可靠性有着至关重要的作用。
目前,蠕变领域的研究人员已经建立了很多关于材料蠕变性能、寿命的预测方法,其中包括基于蠕变机理的双曲正弦模型,等温直线外推法、基于经验的时间温度参数法等。现有预测方法主要存在以下三个问题:(1)由于稳态蠕变速率受温度和应力载荷的影响较大,一旦材料的蠕变机制发生改变,基于蠕变机理的双曲正弦模型可能无法精确拟合较宽范围实验条件内的蠕变速率,从而影响稳态蠕变速率拟合方程的结果准确性。(2)受主观因素影响,等温直线外推法精度较低,具有较大的局限性。(3)虽然时间温度参数法的适用范围较广,但是该方法仅能计算材料的持久性能,无法计算材料的稳态蠕变速率。
因此,亟需研究适用于较宽实验范围蠕变速率拟合方程的计算方法。
发明内容
本发明是为了克服现有技术中存在的缺点而提出的,其目的是提供一种金属材料稳态蠕变速率拟合方程的计算方法。
本发明的技术方案是:
一种金属材料稳态蠕变速率拟合方程的计算方法,包括以下步骤:
(ⅰ)验证材料蠕变持久性能的Monkman-Grant关系
将材料在不同实验条件下测量得到的稳态蠕变速率、持久断裂时间标示在双对数坐标图上,采用线性拟合的方法,验证材料的蠕变持久性能是否满足Monkman-Grant关系;
(ⅱ)将时间温度参数修正为速率温度参数
根据Monkman-Grant的关系式,利用稳态蠕变速率项替换传统时间温度参数模型中的持久时间项,将时间温度参数修正为对应的速率温度参数;
(ⅲ)选择主曲线P(σ)的形式
采用多阶次lgσ的多项式拟合主曲线P(σ)函数,
;
式中:a0、a1、a2、a3、…、an为多项式的待定系数,σ为试验应力水平,单位为MPa;
(ⅳ)拟合速率温度参数模型中的未知常数项
将材料的蠕变实验应力、实验温度、对应条件下的稳态蠕变速率代入到公式中,拟合公式中所有的未知常数项;
(ⅴ)建立材料的稳态蠕变速率拟合方程
将步骤(ⅳ)计算得到的未知常数项拟合值代入到公式中,并整理该方程,建立材料的稳态蠕变速率拟合方程。
所述Monkman-Grant的关系式如下:
;
其中:代表材料在一定实验条件下的稳态蠕变速率,单位为1/h;代表相同条件下的持久断裂时间,单位为h;m为拟合直线的斜率,b为常数。
所述拟合直线的斜率m的取值范围在-0.7~-1之间。
所述步骤(ⅱ)中时间温度参数,包含Larson-Miller参数、Orr-Sherby-Dorn参数、Manson-Haferd参数、Manson-Succop参数和Restrained-Manson-Brown参数。
所述步骤(ⅲ)中用于拟合主曲线P(σ)函数的lgσ多项式阶次n不大于7。
所述步骤(ⅳ)中拟合未知常数项的方法采用最小二乘法、遗传算法、退火算法或者粒子群优化算法中的任意一种。
所述蠕变实验应在至少3个不同的温度条件下进行,每个温度水平下试样数量不少于9个。
本发明的有益效果是:
本发明从金属材料蠕变持久性能的Monkman-Grant关系出发,利用时间温度参数模型应用广泛、拟合程度良好的优点,推导出对应的速率温度参数模型,并根据速率温度参数模型计算金属材料的稳态蠕变速率拟合方程。该计算方法解决了传统时间温度参数无法计算稳态蠕变速率的问题。基于速率温度参数模型计算得到的稳态蠕变速率拟合方程在设备构件可靠性设计、优化方面的应用更加广泛。结合时间温度参数的优点,基于速率温度参数的计算方法可以拟合较宽实验条件范围的稳态蠕变速率,降低了蠕变机理改变对蠕变速率的影响,提高了材料稳态蠕变速率拟合方程的适用性和准确度。
本发明提出的方法具有通用的研究和应用价值。本发明不仅适用于普通金属材料的蠕变性能研究,也适用于计算其他金属基复合材料及非金属基复合材料的稳态蠕变速率拟合方程。
附图说明
图1 是本发明的方法流程图;
图2 是本发明实施例1、2使用的高强铝合金7075蠕变数据分布情况;
图3 是本发明实施例1、2使用的高强铝合金7075稳态蠕变速率与持久时间的拟合关系;
图4 是本发明实施例1中速率温度参数MH’模型下等温预测线与实测数据的对比;
图5 是本发明实施例2中速率温度参数LM’模型下等温预测线与实测数据的对比。
具体实施方式
下面结合说明书附图及实施例对本发明金属材料稳态蠕变速率拟合方程的计算方法进行详细说明:
如图1所示,金属材料稳态蠕变速率拟合方程的计算方法,包括以下步骤:
(ⅰ)验证材料蠕变持久性能的Monkman-Grant关系
将材料在不同实验条件下测量得到的稳态蠕变速率、持久断裂时间标示在双对数坐标图上,采用线性拟合的方法,验证材料的蠕变持久性能是否满足Monkman-Grant关系;所述Monkman-Grant的关系式如下:
其中:代表材料在一定实验条件下的稳态蠕变速率,单位为1/h;代表相同条件下的持久断裂时间,单位为h;m为拟合直线的斜率,b为常数;
(ⅱ)将时间温度参数修正为速率温度参数
根据Monkman-Grant的关系式,利用稳态蠕变速率项替换传统时间温度参数模型中的持久时间项,将时间温度参数修正为对应的速率温度参数;
(ⅲ)选择主曲线P(σ)的形式
采用多阶次lgσ的多项式拟合主曲线P(σ)函数,
;
式中:a0、a1、a2、a3、…、an为多项式的待定系数,σ为试验应力水平,单位为MPa;
(ⅳ)拟合速率温度参数模型中的未知常数项
将材料的蠕变实验应力、实验温度、对应条件下的稳态蠕变速率代入到公式中,拟合公式中所有的未知常数项;
(ⅴ)建立材料的稳态蠕变速率拟合方程
将步骤(ⅳ)计算得到的未知常数项拟合值代入到公式中,并整理该方程,建立材料的稳态蠕变速率拟合方程。
所述拟合直线的斜率m的取值范围在-0.7~-1之间。
所述步骤(ⅱ)中时间温度参数,包含Larson-Miller参数、Orr-Sherby-Dorn参数、Manson-Haferd参数、Manson-Succop参数和Restrained-Manson-Brown参数。
所述步骤(ⅲ)中用于拟合主曲线P(σ)函数的lgσ多项式阶次n不大于7。
所述步骤(ⅳ)中拟合未知常数项的方法采用最小二乘法、遗传算法、退火算法或者粒子群优化算法中的任意一种。
所述蠕变实验应在至少3个不同的温度条件下进行,每个温度水平下试样数量不少于9个。
实施例1
选取高强铝合金7075铝合金板材(T6态)作为研究对象,蠕变试样的取样方向为板内短横向(垂直于挤压板的挤压方向),试样尺寸参照GB/T 2039-1997“金属拉伸蠕变及持久试验方法”。实验温度水平分别为34.4℃、99.4℃、148.9℃,共获得31个有效蠕变数据,数据分布见图2。
(ⅰ)将不同温度下的7075铝合金试样的稳态蠕变速率、持久断裂时间标示在双对数坐标图上,如图3所示。采用线性拟合的方法计算得到的拟合直线表达式,其中线性拟合度R2=0.941,拟合情况如图3所示;
(ⅱ)根据Monkman-Grant的关系式,利用稳态蠕变速率项替换时间温度参数Manson-Haferd(MH)模型中的持久时间项,修正后的速率温度参数MH’模型的表达式为:
;
式中:T—绝对温度;lgta、Ta为待定的常数项;—稳态蠕变速率;
(ⅲ)采用四阶次lgσ的多项式拟合主曲线P(σ)函数,即:
式中:a0、a1、a2、a3、a4为多项式的待定系数,σ为试验应力水平,单位为MPa;
(ⅳ)将7075铝合金试样的蠕变数据代入到基于MH’参数模型的表达式:
;
采用最小二乘法,使与之间的残差平方和最小,拟合得到参数项lgta、Ta、a0、a1、a2、a3、a4,从图4中可以看出,在主要实验温度水平下,根据以上参数项,拟合得到的预测等温线与蠕变速率实测值吻合程度良好,对应的与之间的残差平方和最小值为0.3487;
(ⅴ)将步骤(ⅳ)计算得到的未知常数项拟合值代入到MH’参数模型的表达式中,整理该方程,得到基于MH’参数模型的高强铝合金7075稳态蠕变速率拟合方程:
。
实施例2
选取高强铝合金7075铝合金板材(T6态)作为研究对象,蠕变试样的取样方向为板内短横向(垂直于挤压板的挤压方向),试样尺寸参照GB/T 2039-1997“金属拉伸蠕变及持久试验方法”。实验温度水平分别为34.4℃、99.4℃、148.9℃,共获得31个有效蠕变数据,数据分布见图2。
(ⅰ)将不同温度下的7075铝合金试样的稳态蠕变速率、持久断裂时间标示在双对数坐标图上,如图3所示。采用线性拟合的方法计算得到的拟合直线表达式,其中线性拟合度R2=0.941,拟合情况如图3所示。
(ⅱ)根据Monkman-Grant的关系式,利用稳态蠕变速率项替换时间温度参数Larson-Miller(LM)模型中的持久时间项,修正后的速率温度参数LM’模型的表达式为:
;
式中:T—绝对温度;C为待定的常数项;—稳态蠕变速率。
(ⅲ)采用四阶次lgσ的多项式拟合主曲线P(σ)函数,即:
;
式中:a0、a1、a2、a3、a4为多项式的待定系数,σ为试验应力水平;
(ⅳ)将7075铝合金试样的蠕变数据代入到基于MH’参数模型的表达式:
;
采用遗传算法,使与之间的残差平方和最小,交叉概率设置为0.95,变异概率为0.08,经100次迭代得到参数项lgta、Ta、a0、a1、a2、a3、a4,从图5中可以看出,在主要实验温度水平下,根据以上参数项,拟合得到的预测等温线与蠕变速率实测值吻合程度良好,对应的与之间的残差平方和最小值为0.8285。
(ⅴ)将步骤(ⅳ)计算得到的未知常数项拟合值代入到LM’参数模型的表达式中,整理该方程,得到基于LM’参数模型的高强铝合金7075稳态蠕变速率拟合方程:
。
本发明从金属材料蠕变持久性能的Monkman-Grant关系出发,利用时间温度参数模型应用广泛、拟合程度良好的优点,推导出对应的速率温度参数模型,并根据速率温度参数模型计算金属材料的稳态蠕变速率拟合方程。该计算方法解决了传统时间温度参数无法计算稳态蠕变速率的问题。基于速率温度参数模型计算得到的稳态蠕变速率拟合方程在设备构件可靠性设计、优化方面的应用更加广泛。结合时间温度参数的优点,基于速率温度参数的计算方法可以拟合较宽实验条件范围的稳态蠕变速率,降低了蠕变机理改变对蠕变速率的影响,提高了材料稳态蠕变速率拟合方程的适用性和准确度。
本发明提出的方法具有通用的研究和应用价值。本发明不仅适用于普通金属材料的蠕变性能研究,也适用于计算其他金属基复合材料及非金属基复合材料的稳态蠕变速率拟合方程。
Claims (5)
1.一种金属材料稳态蠕变速率拟合方程的计算方法,其特征在于:包括以下步骤:
(ⅰ)验证材料蠕变持久性能的Monkman-Grant关系
将材料在不同实验条件下测量得到的稳态蠕变速率持久断裂时间tr标示在双对数坐标图上,采用线性拟合的方法,验证材料的蠕变持久性能是否满足Monkman-Grant关系;
(ⅱ)将时间温度参数修正为速率温度参数
根据Monkman-Grant的关系式,利用稳态蠕变速率项替换传统时间温度参数模型中的持久时间项lgtr,将时间温度参数TTP(tr,T)修正为对应的速率温度参数
(ⅲ)选择主曲线P(σ)的形式
采用多阶次lgσ的多项式拟合主曲线P(σ)函数,
P(σ)=a0+a1·lgσ+a2·lg2σ+a3·lg3σ+...+an·lgnσ
式中:a0、a1、a2、a3、…、an为多项式的待定系数,σ为试验应力水平,单位为MPa;
(ⅳ)拟合速率温度参数模型中的未知常数项
将材料的蠕变实验应力、实验温度、对应条件下的稳态蠕变速率代入到公式中,拟合公式中所有的未知常数项;
(ⅴ)建立材料的稳态蠕变速率拟合方程
将步骤(ⅳ)计算得到的未知常数项拟合值代入到公式中,并整理该方程,建立材料的稳态蠕变速率拟合方程;
所述Monkman-Grant的关系式如下:
其中:代表材料在一定实验条件下的稳态蠕变速率,单位为1/h;tr代表相同条件下的持久断裂时间,单位为h;m为拟合直线的斜率,b为常数;
所述拟合直线的斜率m的取值范围在-0.7~-1之间。
2.根据权利要求1所述的一种金属材料稳态蠕变速率拟合方程的计算方法,其特征在于:所述步骤(ⅱ)中时间温度参数,包含Larson-Miller参数、Orr-Sherby-Dorn参数、Manson-Haferd参数、Manson-Succop参数和Restrained-Manson-Brown参数。
3.根据权利要求1所述的一种金属材料稳态蠕变速率拟合方程的计算方法,其特征在于:所述步骤(ⅲ)中用于拟合主曲线P(σ)函数的lgσ多项式阶次n不大于7。
4.根据权利要求1所述的一种金属材料稳态蠕变速率拟合方程的计算方法,其特征在于:所述步骤(ⅳ)中拟合未知常数项的方法采用最小二乘法、遗传算法、退火算法或者粒子群优化算法中的任意一种。
5.根据权利要求1所述的一种金属材料稳态蠕变速率拟合方程的计算方法,其特征在于:所述蠕变实验应在至少3个不同的温度条件下进行,每个温度条件下试样数量不少于9个。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610828973.3A CN106446390B (zh) | 2016-09-19 | 2016-09-19 | 金属材料稳态蠕变速率拟合方程的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610828973.3A CN106446390B (zh) | 2016-09-19 | 2016-09-19 | 金属材料稳态蠕变速率拟合方程的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106446390A CN106446390A (zh) | 2017-02-22 |
CN106446390B true CN106446390B (zh) | 2019-03-29 |
Family
ID=58168589
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610828973.3A Active CN106446390B (zh) | 2016-09-19 | 2016-09-19 | 金属材料稳态蠕变速率拟合方程的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106446390B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110411864B (zh) * | 2018-04-26 | 2022-02-25 | 天津大学 | 一种基于蠕变激活能的高温蠕变寿命预测解析计算方法 |
CN109142083B (zh) * | 2018-09-05 | 2020-04-24 | 南京航空航天大学 | 一种变载荷历程下的蠕变损伤计算方法 |
CN112926230A (zh) * | 2019-12-06 | 2021-06-08 | 核工业理化工程研究院 | 结构件的材料蠕变性能指标确立方法 |
CN113299351B (zh) * | 2020-02-21 | 2022-08-16 | 核工业理化工程研究院 | 基于统计检验及状态检查的蠕变数据筛选方法 |
CN113008677B (zh) * | 2021-03-02 | 2022-06-10 | 中南大学 | 一种镍基高温合金的蠕变持久性能预测方法 |
CN113252465B (zh) * | 2021-05-20 | 2022-07-01 | 天津理工大学 | 一种基于m-h法的耐热钢蠕变寿命预测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101308127A (zh) * | 2008-07-08 | 2008-11-19 | 华东理工大学 | 一种材料的高温蠕变寿命加速试验方法 |
JP4367845B2 (ja) * | 2004-06-03 | 2009-11-18 | 九州電力株式会社 | クリープひずみ速度を利用した金属材料の余寿命評価方法 |
CN101625351A (zh) * | 2008-07-08 | 2010-01-13 | 华东理工大学 | 一种蠕变数据转换为材料高温应力松弛数据的方法 |
CN104931336A (zh) * | 2015-06-15 | 2015-09-23 | 天津理工大学 | 一种测试金属材料在低应力下的蠕变试验装置 |
CN105004617A (zh) * | 2015-07-20 | 2015-10-28 | 沈阳工业大学 | 一种描述金属材料蠕变曲线的方法 |
-
2016
- 2016-09-19 CN CN201610828973.3A patent/CN106446390B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4367845B2 (ja) * | 2004-06-03 | 2009-11-18 | 九州電力株式会社 | クリープひずみ速度を利用した金属材料の余寿命評価方法 |
CN101308127A (zh) * | 2008-07-08 | 2008-11-19 | 华东理工大学 | 一种材料的高温蠕变寿命加速试验方法 |
CN101625351A (zh) * | 2008-07-08 | 2010-01-13 | 华东理工大学 | 一种蠕变数据转换为材料高温应力松弛数据的方法 |
CN104931336A (zh) * | 2015-06-15 | 2015-09-23 | 天津理工大学 | 一种测试金属材料在低应力下的蠕变试验装置 |
CN105004617A (zh) * | 2015-07-20 | 2015-10-28 | 沈阳工业大学 | 一种描述金属材料蠕变曲线的方法 |
Non-Patent Citations (2)
Title |
---|
Evaluation of creep behavior of near-eutectic Sn-Ag solders containing small amount of alloy additions;F.Guo;《MATERIALS SCIENCE AND ENGINEERING A-STRUCTURAL MATERIALS PROPERTIES MICROSTRUCTURE AND PROCESSING》;20031231;第351卷(第1-2期);全文 |
铸态稀土镁合金热变形应变补偿型本构模型;张龙等;《热加工工艺》;20150827;第44卷(第16期);全文 |
Also Published As
Publication number | Publication date |
---|---|
CN106446390A (zh) | 2017-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106446390B (zh) | 金属材料稳态蠕变速率拟合方程的计算方法 | |
Stenzel et al. | Predicting effective conductivities based on geometric microstructure characteristics | |
CN104236889B (zh) | 橡胶密封圈的贮存寿命优化设计方法 | |
CN103617816B (zh) | 反应堆堆芯功率分布的测量方法 | |
CN103530449B (zh) | 一种弹上寿命件的多变量加速贮存试验优化设计方法 | |
Wu et al. | Structural and dynamical properties of liquid Cu 80 Si 20 alloy studied experimentally and by ab initio molecular dynamics simulations | |
Li et al. | Hot deformation characteristics—Constitutive equation and processing maps—Of 21-4N heat-resistant steel | |
Gomez et al. | Natural gas intermittent kiln for the ceramic industry: A transient thermal analysis | |
Maheshwari et al. | Modified Johnson–Cook material flow model for hot deformation processing | |
Manzolaro et al. | Thermal and mechanical characterization of carbides for high temperature nuclear applications | |
CN106202725B (zh) | 基于多源数据的成败型产品可靠性试验设计和决策方法 | |
Gombár et al. | Experimental analysis of the influence of factors acting on the layer thickness formed by anodic oxidation of aluminium | |
Kugalur Palanisamy et al. | Identification of the parameter values of the constitutive and friction models in machining using EGO algorithm: application to Ti6Al4V | |
Ren et al. | A comprehensive review of key technologies for enhancing the reliability of lithium-ion power batteries | |
Beardsley et al. | A deformation mechanism map for Incoloy 800H optimized using the genetic algorithm | |
Prataviera et al. | A new bimodal Maxwell regression model with engineering applications | |
Li et al. | A review of the latest developments in the field of additive manufacturing techniques for nuclear reactors | |
Xiong et al. | A new method to determine isothermal flow curves for integrated process and microstructural simulation in metal forming | |
CN116150830A (zh) | 一种SiC复合包壳的多物理场耦合计算方法 | |
Park et al. | Flow Stress optimization of Inconel 718 based on a coupled simulation of material-forming analysis and joule heating analysis | |
Do Carmo et al. | Thermal Simulation Studies for the Characterization of a Cyclotron Liquid Target with Thick Niobium Target Windows | |
CN113722961A (zh) | 一种基于广义协同高斯过程模型的结构不确定性量化解析方法 | |
CN102360430A (zh) | 一种基于支持向量机和误差估计的供热负荷区间预报方法 | |
Shahani et al. | Investigation on the mechanical properties and fracture toughness of graphite | |
Gołębiowski et al. | Transient Thermal Field Analysis in ACCC Power Lines by the Green’s Function Method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |