CN110245382A - 一种确定金属动态再结晶体积分数的Avrami数学模型系数的方法 - Google Patents
一种确定金属动态再结晶体积分数的Avrami数学模型系数的方法 Download PDFInfo
- Publication number
- CN110245382A CN110245382A CN201910387211.8A CN201910387211A CN110245382A CN 110245382 A CN110245382 A CN 110245382A CN 201910387211 A CN201910387211 A CN 201910387211A CN 110245382 A CN110245382 A CN 110245382A
- Authority
- CN
- China
- Prior art keywords
- volume fraction
- curve
- strain
- recrystallization
- avrami
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
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)
- Analysing Materials By The Use Of Radiation (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供一种确定金属动态再结晶体积分数的Avrami数学模型系数的方法,该方法包括如下步骤:步骤1.获得材料的应力应变曲线,试样在变形后,进行水淬处理,观察试样内部的金相组织,测量再结晶体积分数;步骤2.确定动态再结晶的临界应变εc;步骤3.确定回复曲线;步骤4.计算动态再结晶体积分数随应变变化曲线;步骤5.根据步骤4得到的曲线,求解Avrami数学模型系数m和ε*。本发明所述的方法中只需测量变形完成后的试样中的再结晶体积分数,无需变形过程中的试样的再结晶体积分数,直接使用方程即可得到再结晶体积分数随应变变化曲线,无需稳态应力σss参与计算;本发明所述的方法与现有方法相比,计算过程所需参数少、工作量低,试样消耗少,适宜推广。
Description
技术领域
本发明涉及材料微观组织分析技术领域,具体涉及一种确定金属动态再结晶体积分数的Avrami数学模型系数的方法。
背景技术
金属与合金在发生塑性变形后,会在其内部积累大量的畸变能。这时,如果将其加热至适当温度并进行保温处理,其内部的畸变组织会发生再次形核和晶粒长大,生成无畸变的晶粒,这一过程称为再结晶过程;再结晶过程与变形过程同时发生,这样的再结晶称为动态再结晶;动态再结晶体积分数是再结晶晶粒的总体积占试样总体积的百分比,其与加热温度、保温时间和变形程度密切相关,通常可以使用Avrami数学模型描述,
上式中,X代表动态再结晶体积分数,ε为应变,εc为发生动态再结晶的临界应变,m和ε*均为方程系数。
动态再结晶体积分数数学模型的确定,对材料变形组织控制、塑性成形能力改善和成形工艺优化具有重要的作用。
确定方程系数εc、m和ε*的现有方法有两种。
第一种方法,金相试验法
金相试验法通过大量的试验,确定Avrami模型的系数,步骤如下:
1.在特定的变形条件下(温度和应变速率),进行变形试验;
2.高温变形后的试样,立即进行水淬,以保留变形组织;
3.对试样进行打磨、抛光、腐蚀,并观察金相组织,确定金相组织中再结晶晶粒的体积百分数;
4.在多种变形条件下(温度、应变和应变速率),重复前三步骤;
5.在不同温度和应变速率条件下,确定动态再结晶体积分数与应变之间的Avrami数学关系。
第二种方法,应力应变曲线分析法,具体步骤如下:
1.在适当的温度和应变速率范围内,进行变形试验,获得材料的应力应变曲线;
2.确定动态再结晶的临界应变εc:
在此方法中,εc代表发生动态再结晶的临界应变,首先,针对特定的变形条件(温度和应变速率),使用式(1)计算相应的应力应变曲线各处的斜率;
然后使用θ作为纵坐标,σ作为横坐标,绘制曲线,εc即为曲线拐点对应的应变值;
3.确定回复曲线:
假设材料在变形过程中不发生再结晶,仅发生回复,这时的应力应变曲线称为回复曲线,对应力应变曲线最高点之前的数据,使用方程(2)进行拟合,即可得到材料的回复曲线;
拟合时,r值作为曲线的斜率求出,作为截距求出,最终的拟合曲线使用方程(3)计算得到;
4.计算动态再结晶体积分数随应变变化曲线:
在前两步骤的基础之上,对于每一个选定的应变值,使用方程(4)或者(5),计算与之相对应的动态再结晶体积分数,方程中各参数的物理意义如图1所示;
式中,σsat为饱和应力,σss为稳态应力。
5.通过前一步骤获得的动态再结晶体积分数随应变的变化曲线,求解Avrami方程中的其他未知系数m和ε*。
上述的两种方法均可以求得Avrami方程的系数,但是均存在无法避免的缺陷。前一种金相试验法,不仅工作量极为浩大,而且经费消耗量也很巨大。第二种应力应变曲线分析法,需要稳态应力σss参与运算,但是大多数金属和合金的应力应变曲线中,并不存在稳态应力σss,所以其适用性较差。
发明内容
本发明针对现有技术的不足提供一种金属动态再结晶体积分数的Avrami数学模型系数的确定方法。
本发明解决上述技术问题所采用的技术方案是:一种金属动态再结晶体积分数的Avrami数学模型系数的确定方法,其中Avrami数学模型方程为如下所示,包括如下步骤:
步骤1.通过试验,获得材料的应力应变曲线,试样在变形后,立即进行水淬处理,以保留其变形组织,随后观察试样内部的金相组织,测量再结晶体积分数;
步骤2.采用应力应变曲线分析法确定动态再结晶的临界应变εc;
步骤3.采用应力应变曲线分析法确定回复曲线σ回复(ε),及或σsat,σsat为饱和应力;
步骤4.计算动态再结晶体积分数随应变变化曲线X(ε),具体方法为:通过方程(Ⅰ)或(Ⅱ),即可得到再结晶体积分数随应变变化曲线;
式中,X(a)表示真应变为a时,变形后试样中再结晶体积分数;σ再结晶(ε)曲线即为步骤1中获得的曲线,σ再结晶(a)的数值即可随之获得;
步骤5.根据步骤4得到的曲线,求解Avrami数学模型系数m和ε*,方法如下:
首先,将Avrami数学模型方程转换成以下形式:
对方程(Ⅲ)两侧,两次取对数:
ln[-ln(1-X)]=m*ln(ε-εc)-m*lnε*#(Ⅳ)
然后,以ln[-ln(1-X)]作为纵坐标,ln(ε-εc)作为横坐标,绘制曲线,并进行线性拟合,m值即可作为拟合直线的斜率求出,-m*lnε*值为直线的截距,即可求出ε*值。
所述方程(Ⅰ)、(Ⅱ)可通过方程(4)、(5)推导得到,方程(5)的操作方法与方程(4)一致,下面使用方程(4)进行说明,推导过程如下:
首先,假设试样的真应变为a,变形后试样中再结晶体积分数为X(a),根据方程(4),X(a)可以表达为:
对于应力应变曲线上的任意应变ε,动态再结晶体积分数可以表达为:
方程(Ⅴ)与方程(Ⅵ)相除,得到方程(Ⅰ):
方程(Ⅰ)中,X(a)的数值已经在步骤1中测定,方程右侧的四个应力也是已知量,所以,通过方程(Ⅰ),即可得到再结晶体积分数随应变变化曲线,不再需要σss参与计算;
如果使用方程(5),使用上述方法,最终会得到方程(Ⅱ)。
本发明所述的方法中只需测量变形完成后的试样中的再结晶体积分数,无需变形过程中的试样的再结晶体积分数,可直接使用方程(Ⅰ)或(Ⅱ)即可得到再结晶体积分数随应变变化曲线,无需稳态应力σss参与计算;本发明所述的方法与现有方法相比,计算过程简单、工作量低,试样消耗少,适宜推广。
附图说明
图1是方程(4)和(5)各参数的物理意义图示;
图2是实施例材料的真应力真应变曲线图;
图3是实施例材料淬火后的组织图;
图4是实施例材料淬火后的再结晶晶粒图;
图5是实施例中确定动态再结晶的临界点曲线;
图6是实施例中推导的回复曲线与原始试验曲线对比图;
图7是实施例中再结晶体积分数随应变的演化曲线。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
实施例
本发明提供一种金属动态再结晶体积分数的Avrami数学模型系数的确定方法,具体步骤如下:
步骤1.通过试验,获得材料的应力应变曲线,如图2所示;试样在变形后,立即进行水淬处理,以保留其变形组织,淬火后的组织如图3所示,再结晶晶粒如图4所示,随后观察试样内部的金相组织,测量得到的再结晶体积分数为28.1%;
步骤2.确定动态再结晶的临界应变εc,首先,针对特定的变形条件(温度和应变速率),使用公式(1)计算相应的应力应变曲线各处的斜率θ;
然后使用θ作为纵坐标,σ作为横坐标,绘制曲线,εc即为曲线拐点对应的应变值;最终确定的临界应变为0.314,对应的应力为243.6MPa,如图5所示;
步骤3.确定回复曲线σ回复(ε)及或σsat,σsat为饱和应力:
具体方法为:确定方程(2)中的两个参数r和对应力应变曲线最高点之前的数据,使用方程(2)进行拟合:
式中r为常数,拟合时,r值作为曲线的斜率求出,作为截距求出;
使用方程(3)确定回复曲线,式中σ0为屈服强度;
计算得到的回复曲线与原始曲线对比,如图6所示;
步骤4.计算动态再结晶体积分数随应变变化曲线:
σ回复(ε)为图6中的回复曲线,σ再结晶(ε)为图6中的原始曲线,即为步骤(1)中获得的曲线,参考方程(Ⅰ),a代表最大变形量,即为图2中真应变最大值,σ再结晶(a)为与之对应的真应力最大值,X(a)为最大变形量对应的再结晶体积分数,即为步骤1中测得的28.1%;σsat已经在上一步骤中获得;因此,任意应变ε对应的再结晶体积分数X(ε)即为可求,求得的结果如图7所示。
步骤5.以ln[-ln(1-X)]作为纵坐标,ln(ε-εc)作为横坐标,绘制曲线,并进行线性拟合,求得m=2.12,ε*=0.43,该材料的动态再结晶体积分数Avrami模型为:
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (1)
1.一种确定金属动态再结晶体积分数的Avrami数学模型系数的方法,其中Avrami数学模型方程为如下所示,其特征在于,包括如下步骤:
步骤1.通过试验,获得材料的应力应变曲线,试样在变形后,立即进行水淬处理,以保留其变形组织,随后观察试样内部的金相组织,测量再结晶体积分数;
步骤2.采用应力应变曲线分析法确定动态再结晶的临界应变εc:
步骤3.采用应力应变曲线分析法确定回复曲线σ回复(ε)及或σsat,σsat为饱和应力:
步骤4.计算动态再结晶体积分数随应变变化曲线,具体方法为:通过方程(Ⅰ)或(Ⅱ),即可得到再结晶体积分数随应变变化曲线X(ε);
式中,X(a)表示真应变为a时,变形后试样中再结晶体积分数;σ再结晶(ε)曲线即为步骤1中获得的曲线,σ再结晶(a)可根据σ再结晶(ε)曲线获得;
步骤5.根据步骤4得到的曲线,求解Avrami数学模型系数m和ε*,方法如下:
首先,将Avrami数学模型方程转换成以下形式:
对方程(11)两侧,两次取对数:
ln[-ln(1-X)]=m*ln(ε-εc)-m*lnε*#(Ⅳ)
然后,以ln[-ln(1-X)]作为纵坐标,ln(ε-εc)作为横坐标,绘制曲线,并进行线性拟合,m值即可作为拟合直线的斜率求出,-m*lnε*值为直线的截距,即可获得ε*值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910387211.8A CN110245382B (zh) | 2019-05-10 | 2019-05-10 | 一种确定金属动态再结晶体积分数的Avrami数学模型系数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910387211.8A CN110245382B (zh) | 2019-05-10 | 2019-05-10 | 一种确定金属动态再结晶体积分数的Avrami数学模型系数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110245382A true CN110245382A (zh) | 2019-09-17 |
CN110245382B CN110245382B (zh) | 2023-08-22 |
Family
ID=67884146
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910387211.8A Active CN110245382B (zh) | 2019-05-10 | 2019-05-10 | 一种确定金属动态再结晶体积分数的Avrami数学模型系数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110245382B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110702727A (zh) * | 2019-09-30 | 2020-01-17 | 鞍钢股份有限公司 | 一种测定材料静态再结晶体积分数的方法 |
CN110749573A (zh) * | 2019-11-21 | 2020-02-04 | 山东建筑大学 | 一种非晶态合金结晶动力学特征的原位同步辐射表征方法 |
CN110763708A (zh) * | 2019-11-21 | 2020-02-07 | 山东建筑大学 | 一种用exafs表征金属玻璃有序化过程的方法 |
CN110793982A (zh) * | 2019-11-21 | 2020-02-14 | 山东建筑大学 | 一种纳米晶化动力学过程的高能x射线表征方法 |
CN110793989A (zh) * | 2019-11-21 | 2020-02-14 | 山东建筑大学 | 用原位中子衍射表征大块金属玻璃有序化过程的方法 |
CN110793988A (zh) * | 2019-11-21 | 2020-02-14 | 山东建筑大学 | 一种用高能x射线表征非晶-晶体转变动力学特性的方法 |
CN110793990A (zh) * | 2019-11-21 | 2020-02-14 | 山东建筑大学 | 块体非晶合金晶化动力学特性的中子衍射表征方法 |
CN111337356A (zh) * | 2020-04-30 | 2020-06-26 | 西北工业大学 | 一种基于动态再结晶速率的金属热加工性能表征方法 |
CN113901660A (zh) * | 2021-10-11 | 2022-01-07 | 本钢板材股份有限公司 | 一种用于确定材料voce本构方程系数的新方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006337343A (ja) * | 2005-06-06 | 2006-12-14 | Yazaki Corp | 結晶性高分子から成る構造体の真応力−対数ひずみ曲線の推定システム |
CN101046682A (zh) * | 2007-04-28 | 2007-10-03 | 武汉科技大学 | 一种预测热轧含Nb带钢组织及力学性能的方法 |
CN101591729A (zh) * | 2009-06-19 | 2009-12-02 | 东北大学 | 预测板带钢热变形中奥氏体动态再结晶组织演变的方法 |
CN102519801A (zh) * | 2011-12-16 | 2012-06-27 | 南京航空航天大学 | 用真应力-真应变曲线建立铝合金动态再结晶模型的方法 |
CN108169019A (zh) * | 2017-12-26 | 2018-06-15 | 中国兵器工业第五九研究所 | 一种准静态塑性压缩应力应变参量的识别方法 |
CN114864007A (zh) * | 2022-05-25 | 2022-08-05 | 湖北汽车工业学院 | 超高强钢高温成形晶粒尺寸演化模型建立方法 |
-
2019
- 2019-05-10 CN CN201910387211.8A patent/CN110245382B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006337343A (ja) * | 2005-06-06 | 2006-12-14 | Yazaki Corp | 結晶性高分子から成る構造体の真応力−対数ひずみ曲線の推定システム |
CN101046682A (zh) * | 2007-04-28 | 2007-10-03 | 武汉科技大学 | 一种预测热轧含Nb带钢组织及力学性能的方法 |
CN101591729A (zh) * | 2009-06-19 | 2009-12-02 | 东北大学 | 预测板带钢热变形中奥氏体动态再结晶组织演变的方法 |
CN102519801A (zh) * | 2011-12-16 | 2012-06-27 | 南京航空航天大学 | 用真应力-真应变曲线建立铝合金动态再结晶模型的方法 |
CN108169019A (zh) * | 2017-12-26 | 2018-06-15 | 中国兵器工业第五九研究所 | 一种准静态塑性压缩应力应变参量的识别方法 |
CN114864007A (zh) * | 2022-05-25 | 2022-08-05 | 湖北汽车工业学院 | 超高强钢高温成形晶粒尺寸演化模型建立方法 |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110702727B (zh) * | 2019-09-30 | 2022-02-18 | 鞍钢股份有限公司 | 一种测定材料静态再结晶体积分数的方法 |
CN110702727A (zh) * | 2019-09-30 | 2020-01-17 | 鞍钢股份有限公司 | 一种测定材料静态再结晶体积分数的方法 |
CN110749573B (zh) * | 2019-11-21 | 2022-03-04 | 山东建筑大学 | 一种非晶态合金结晶动力学特征的原位同步辐射表征方法 |
CN110763708A (zh) * | 2019-11-21 | 2020-02-07 | 山东建筑大学 | 一种用exafs表征金属玻璃有序化过程的方法 |
CN110793989A (zh) * | 2019-11-21 | 2020-02-14 | 山东建筑大学 | 用原位中子衍射表征大块金属玻璃有序化过程的方法 |
CN110793988A (zh) * | 2019-11-21 | 2020-02-14 | 山东建筑大学 | 一种用高能x射线表征非晶-晶体转变动力学特性的方法 |
CN110793990A (zh) * | 2019-11-21 | 2020-02-14 | 山东建筑大学 | 块体非晶合金晶化动力学特性的中子衍射表征方法 |
CN110793990B (zh) * | 2019-11-21 | 2022-03-18 | 山东建筑大学 | 块体非晶合金晶化动力学特性的中子衍射表征方法 |
CN110793988B (zh) * | 2019-11-21 | 2022-03-15 | 山东建筑大学 | 一种用高能x射线表征非晶-晶体转变动力学特性的方法 |
CN110793982A (zh) * | 2019-11-21 | 2020-02-14 | 山东建筑大学 | 一种纳米晶化动力学过程的高能x射线表征方法 |
CN110763708B (zh) * | 2019-11-21 | 2022-03-04 | 山东建筑大学 | 一种用exafs表征金属玻璃有序化过程的方法 |
CN110793982B (zh) * | 2019-11-21 | 2022-03-04 | 山东建筑大学 | 一种纳米晶化动力学过程的高能x射线表征方法 |
CN110749573A (zh) * | 2019-11-21 | 2020-02-04 | 山东建筑大学 | 一种非晶态合金结晶动力学特征的原位同步辐射表征方法 |
CN111337356A (zh) * | 2020-04-30 | 2020-06-26 | 西北工业大学 | 一种基于动态再结晶速率的金属热加工性能表征方法 |
CN113901660A (zh) * | 2021-10-11 | 2022-01-07 | 本钢板材股份有限公司 | 一种用于确定材料voce本构方程系数的新方法 |
CN113901660B (zh) * | 2021-10-11 | 2024-06-11 | 本钢板材股份有限公司 | 一种用于确定材料voce本构方程系数的新方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110245382B (zh) | 2023-08-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110245382A (zh) | 一种确定金属动态再结晶体积分数的Avrami数学模型系数的方法 | |
CN102519801B (zh) | 用真应力-真应变曲线建立铝合金动态再结晶模型的方法 | |
Yin et al. | Microstructural modeling and simulation for GCr15 steel during elevated temperature deformation | |
Yu et al. | Crack healing in a low-carbon steel under hot plastic deformation | |
Ghazani et al. | Kinetics and critical conditions for initiation of dynamic recrystallization during hot compression deformation of AISI 321 austenitic stainless steel | |
CN1330930C (zh) | 轧制过程钢板内部组织晶粒尺寸的软测量方法 | |
CN101871899A (zh) | 钢铁材料控冷cct曲线测定方法 | |
Massardier et al. | Mn-C interaction in Fe-C-Mn steels: study by thermoelectric power and internal friction | |
Ping et al. | Neural network prediction of flow stress of Ti–15–3 alloy under hot compression | |
CN103940663B (zh) | 不同应力比下材料疲劳门槛值的预测方法 | |
Xu et al. | Determination of the dynamic recrystallization kinetics model for SCM435 steel | |
CN111445960B (zh) | 14Cr17Ni2钢的锻造工艺参数的优化方法 | |
Baracaldo et al. | Charpy impact toughness and transition temperature in ferrite–perlite steel | |
CN108169019B (zh) | 一种准静态塑性压缩应力应变参量的识别方法 | |
Chen et al. | Effects of different process parameters on mechanical properties and microstructures of hot stamping boron steel | |
CN110702727B (zh) | 一种测定材料静态再结晶体积分数的方法 | |
CN110018046B (zh) | 一种节约型双相不锈钢中trip效应致塑性增量的表征方法 | |
Ji et al. | HOT-DEFORMATION BEHAVIOR AND PROCESSING MAP OF 22MnB5 HIGH-STRENGTH STEEL | |
Dong et al. | Evolution of dynamically recrystallized grains in steel SCM435 | |
Ahmad | Modified law of mixture to describe the tensile deformation behavior of thermomechanically processed dual-phase steel | |
CN105642675A (zh) | 防止磷非平衡晶界偏聚热轧工艺控制方法 | |
Matveev et al. | High-temperature plasticity of microalloyed steel | |
Mikó et al. | Determination of the onset of the dynamic recystallization of a 7075 Al alloy | |
Du et al. | Internal friction study of aging hardening and kinetics in low carbon steel | |
Min et al. | Investigation on uniaxial tensile instability of USIBOR1500 steel sheets at elevated temperature |
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 |