CN110068507B - 一种对传统再结晶模型进行修正的方法 - Google Patents
一种对传统再结晶模型进行修正的方法 Download PDFInfo
- Publication number
- CN110068507B CN110068507B CN201810059920.9A CN201810059920A CN110068507B CN 110068507 B CN110068507 B CN 110068507B CN 201810059920 A CN201810059920 A CN 201810059920A CN 110068507 B CN110068507 B CN 110068507B
- Authority
- CN
- China
- Prior art keywords
- recrystallization
- strain
- model
- parameters
- correcting
- 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
Images
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/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
- G01N3/18—Performing tests at high or low temperatures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及材料微观组织分析领域,尤其是一种对传统再结晶模型的修正方法。通过物理模拟试验机获得不同变形条件下试样,观察芯部再结晶分数;利用数值模拟软件获得试样芯部真实应变;将所得的数据进行函数拟合,获得模型参数;对传统拐点法求临界应变的不准确性和压缩试样局部应变不均匀引起的误差进行修正;将所得模型参数跟Zener‑Hollomon参数拟合,获得其他条件下的再结晶模型参数。本发明修正后的再结晶模型,不仅适用于单相合金的再结晶预测,也适用于多相合金的再结晶预测,同时预测精度大幅度提高。结合有限元软件后,可以定量的分析材料的动态再结晶体积分数的变化情况,为材料的实际加工工艺的制定提供科学的依据。
Description
技术领域
本发明涉及材料微观组织分析领域,尤其是一种对传统再结晶分数预测模型的修正方法。
背景技术
常规再结晶分数预测模型建立过程如下:通过对高温压缩试验真应力-应变曲线求导得到材料应变硬化率,通过分析该硬化率拐点变化确定再结晶临界应变和最大软化速率应变,再根据经验公式建立临界应变和最大软化速率应变的数学关系模型继而预测动态再结晶体积分数。因此模型需要的参数少,形式简单,故被广泛应用。[S.I.Kim,et.alMater.Sci.Eng.A 311(2001)108.8051.Z.Yang,Materials Science and Engineering A485(2008)487.G.-Z.Quan,et.al Materials Science and Engineering A,528(2011)]。
常规模型的精确度主要依赖于临界应变和最大软化速率应变的准确性,但由于以下原因常常难于得到准确的应变数据:(i)采用物理模拟热压缩实验,由于摩擦力存在,导致试样整体变形不均匀,所测得数据为包括变形不均匀的整体力学性能数据;(ii)并且对于某些多相合金,应变硬化率曲线拐点不一定是由材料再结晶软化引起的,可能由材料强化相软化造成,传统方法局限性未考虑该点。(iii)原始压缩曲线波动常常较大,无法直接从硬化率曲线得到准确的再结晶临界应变和最大软化率曲线,目前常采用的方法是通过各种函数去平滑原始应力应变曲线,但该近似处理使原始曲线失真较大。因此,十分必要对传统方法进行修正以提高模型的预测精度。
发明内容
本发明的目的在于提供一种对传统再结晶模型进行修正的方法,有效地修正传统拐点法和物理模拟压缩变形不均匀引起的误差。
为了实现上述目的,本发明采用的技术方案是:
一种对传统再结晶模型进行修正的方法,通过物理模拟试验机获得不同变形条件下试样,观察并统计芯部再结晶分数,通过有限元软件获得试样芯部对应的真实应变,将再结晶分数与对应的真实应变进行函数拟合,获得模型参数;将模型参数与对应的Zener-Hollomon参数进行函数拟合,获得之间的关系,通过计算不同条件下的Z参数,获得对应条件下的模型参数。
所述的对传统再结晶模型进行修正的方法,还包括:对传统模型拐点法求临界再结晶应变的不准确性和压缩试样应变的不均匀性引起的误差进行修正。
所述修正由于传统拐点法无法准确获得临界再结晶应变引起的误差,具体为:在同一变形条件下,通过物理模拟实验机进行一系列的不同应变量的实验,测量压缩后试样芯部组织的再结晶分数,初步确定其范围,然后与对应的真实应变进行函数拟合,获得该条件下再结晶临界应变。其中,同一变形条件是指:进行两组同一温度、应变速率下不同应变量的压缩变形,获取实验后试样芯部再结晶分数,通过模拟该条件下的变形过程,得到不同应变下的芯部真实应变,应用拟合得到两组εc、ε*和m值,其中εc为再结晶临界应变,ε*为再结晶软化最大速率,s-1;m为材料常数,取平均值。
所述修正由于压缩试样应变不均匀引起的误差,具体为:通过有限元模拟得到压缩后试样形状同实验结果相比,外形吻合较好,模拟结果可以较准确地反映不同位置的真实应变。
所述函数拟合具体为:
(1)进行两组同一温度、应变速率下不同应变量的压缩变形,获取实验后试样芯部再结晶分数,通过模拟该条件下的变形过程,得到不同应变下的芯部真实应变,应用拟合得到两组εc、ε*和m值,其中εc为再结晶临界应变,ε*为再结晶软化最大速率,s-1;m为材料常数,取平均值;
(2)应用εc=A×ZB和ε*=A'×ZB'拟合得到εc、ε*的与Z参数的关系,A、B、A'、B'均为常数。
本发明具有以下有益效果及优点:
1.预测模型精度高。对于真实应力应变曲线波动较大的多相合金,利用传统的拐点法(硬化率曲线法)无法获得准确结果,本发明通过实验数据拟合获得两组不同条件下的εc、ε*,然后根据经验公式得到其他条件的εc、ε*,比传统拐点法求得的εc、ε*更加精确,同时利用有限元的方法对压缩应变的不均匀性进行修正,大大提高模型预测精度。
2.适用范围广。对于某些多相合金,应变硬化率曲线拐点不一定是由材料再结晶软化引起的,可能由材料强化相软化造成,本发明求得的材料再结晶分数预测模型不仅仅适用于单向合金,也适用于多相合金。
附图说明
图1(a)为压缩后试样形貌;图1(b)为波动较大的原始曲线;图1(c)为350℃-0.01s-1试样芯部组织图片。
图2为再结晶模型修正方法全流程图。
图3为压缩试样真实应变模拟结果图。
图4(a)-图4(b)为光学显微镜(OM)数据点拟合曲线;其中,图4(a)为450℃-0.01s-1条件下再结晶分数随应变量变化曲线;图4(b)为450℃-1s-1条件下的再结晶分数随应变量变化曲线。
图5为450℃不同应变速率下,动态再结晶分数随不同应变量的变化曲线。
具体实施方式
如图1(a)所示,两端由于摩擦力存在,压缩后试样出现鼓肚,变形不均匀。
如图1(b)所示,原始真实应力-应变曲线波动较大,从求导所得的应变硬化率曲线中无法准确获得再结晶开始时对应的拐点。
如图1(b)和图1(c)所示,350℃-0.01s-1对应的真实应力-应变曲线发生软化现象,但经过观察并无再结晶发生。
如图4(a)所示,以传统方法求得的再结晶分数预测曲线,误差较大。
如图2所示,在具体实施过程中,本发明再结晶模型修正方法全流程如下:
本发明方法通过物理模拟试验机模拟压缩得到不同条件下的试样,将试样芯部再结晶分数与对应的真实应变进行函数拟合,获得再结晶模型参数;将模型参数与对应的Zener-Hollomon参数进行函数拟合,获得之间的关系,通过计算不同条件下的Z参数,获得对应条件下的模型参数。
本发明方法还包括:对传统拐点法求临界再结晶应变不准确性和压缩试样变形不均匀性引起的误差进行修正;所述修正由于摩擦导致压缩试样应变不均匀性引起的误差具体为:通过有限元模拟得到压缩后试样形状同实验结果相比,外形吻合较好,模拟结果可以较准确地反应不同位置的真实应变,从而用模拟的局部真实应变来代替传统模型使用的整体应变;所述修正由于传统拐点法求临界再结晶应变不准确引起的误差具体为:通过物理模拟试验机获得两组不同温度和应变速率压缩应变量分别为0.2、0.4、0.6、0.8、1.2条件下的模拟压缩试样,金相观测分别获得其芯部组织再结晶分数,结合有限元软件模拟压缩过程获得对应条件下试样芯部的局部真实应变量,应用拟合,得到两组不同条件下的再结晶临界应变εc和再结晶软化最大速率ε*和m值,对求得的两组m值取平均,作为材料常数。
对两组εc、ε*,应用经验公式εc=A×ZB和ε*=A'×ZB'进行拟合,获得常数A、A'、B、B'。根据εc、ε*与Z参数的关系,结合应力应变曲线获得材料的激活能Q,应用计算其他条件下的Z参数,从而通过Z获得不同条件下的再结晶模型参数。
下面,通过实施例对本发明进一步详细阐述。
实施例
本实施例中,以稀土镁合金为例,其化学成分(质量分数)见表1。
表1稀土镁合金化学成分组织
化学成分(wt%) | Gd | Y | Zn | Zr | Mg |
Mg-13Gd-4Y-2Zn-0.6Zr | 10.8 | 3.93 | 1.94 | 0.37 | Bal. |
在Gleeble-3800热模拟试验机上进行压缩试验,试样的尺寸为Φ8×12mm,采用温度为450℃,应变速率1s-1和0.01s-1,压缩应变量为0.2、0.4、0.6、0.8、1.2。使用有限元模拟热压缩变形过程,得到修正后的应变如图3所示,表2为不同变形条件下对应的再结晶分数。
表2不同条件对应的再结晶分数
实验应变 | 0 | 0.2 | 0.4 | 0.6 | 0.8 | 1.2 |
修正应变 | 0 | 0.24 | 0.62 | 1.09 | 1.64 | 2.17 |
450℃-0.01s<sup>-1</sup>的X<sub>DRX</sub>(OM) | 0 | 0 | 0.25 | 0.40 | 0.46 | 0.68 |
450℃-1s<sup>-1</sup>的X<sub>DRX</sub>(OM) | 0 | 0 | 0.1 | 0.22 | 0.36 | 0.53 |
利用经验公式分别对表2中数据进行拟合,如表3所示,R值均大于0.95,说明拟合曲线与实验数据点吻合较好。将m值取平均代入经验公式,再次进行拟合,结果如图4(a)和图4(b)所示,参数见表4,R值均大于0.95,说明曲线拟合较好。
表3不同条件各自拟合JMAK模型参数值
ε<sub>c</sub> | ε<sup>*</sup> | m | R | |
450℃-0.01s<sup>-1</sup>JMAK参数(OM) | 0.278 | 1.200 | 0.750 | 0.965 |
450℃-1s<sup>-1</sup>JMAK参数(OM) | 0.437 | 1.752 | 0.9 | 0.983 |
表4不同条件同一m值拟合JMAK模型参数
ε<sub>c</sub> | ε<sup>*</sup> | m | R | |
450℃-0.01s<sup>-1</sup>JMAK参数(OM) | 0.24 | 1.242 | 0.85 | 0.969 |
450℃-1s<sup>-1</sup>JMAK参数(OM) | 0.46 | 1.75 | 0.85 | 0.983 |
参数εc和ε*与Z参数具有如下关系:
εc=A×ZB
ε*=A’×ZB
将所求的450℃-0.01s-1参数和450℃-1s-1参数代入公式,获得关系如下:
εc=0.00234×Z0.1
ε*=0.157×Z0.0457
表5为450℃-0.1s-1和450℃-10s-1条件下对应的JMAK模型参数值
图5为450℃不同应变速率下,动态再结晶分数随不同应变量的变化曲线,点代表OM数据实验结果,可以看出不同条件下OM数据点与预测曲线吻合较好。
实施例结果表明,本发明修正后的再结晶模型,不仅仅适用于单相合金的再结晶预测,也适用于多相合金的再结晶预测,同时预测精度大幅度提高,结合有限元软件后,可以定量的分析材料的动态再结晶体积分数的变化情况,为材料的实际加工工艺的制定提供科学的依据。
Claims (4)
1.一种对传统再结晶模型进行修正的方法,其特征在于,通过物理模拟试验机获得不同变形条件下试样,观察并统计芯部再结晶分数,通过有限元软件获得试样芯部对应的真实应变,将再结晶分数与对应的真实应变进行函数拟合,获得模型参数;将模型参数与对应的Zener-Hollomon参数进行函数拟合,获得之间的关系,通过计算不同条件下的Z参数,获得对应条件下的模型参数;
所述函数拟合具体为:
(1)进行两组同一温度、应变速率下不同应变量的压缩变形,获取实验后试样芯部再结晶分数,通过模拟该条件下的变形过程,得到不同应变下的芯部真实应变,应用拟合得到两组εc、ε*和m值,其中εc为再结晶临界应变,ε*为再结晶软化最大速率,s-1;m为材料常数,取平均值;
(2)应用εc=A×ZB和ε*=A'×ZB'拟合得到εc、ε*与Z参数的关系,A、B、A'、B'均为常数。
2.根据权利要求1所述的对传统再结晶模型进行修正的方法,其特征在于,还包括:对传统模型拐点法求临界再结晶应变的不准确性和压缩试样应变的不均匀性引起的误差进行修正。
3.根据权利要求2所述的对传统再结晶模型进行修正的方法,其特征在于,所述修正由于压缩试样应变不均匀引起的误差,具体为:通过有限元模拟得到压缩后试样形状同实验结果相比,外形吻合较好,模拟结果可以较准确地反映不同位置的真实应变。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810059920.9A CN110068507B (zh) | 2018-01-22 | 2018-01-22 | 一种对传统再结晶模型进行修正的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810059920.9A CN110068507B (zh) | 2018-01-22 | 2018-01-22 | 一种对传统再结晶模型进行修正的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110068507A CN110068507A (zh) | 2019-07-30 |
CN110068507B true CN110068507B (zh) | 2021-07-23 |
Family
ID=67365019
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810059920.9A Active CN110068507B (zh) | 2018-01-22 | 2018-01-22 | 一种对传统再结晶模型进行修正的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110068507B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110684938A (zh) * | 2019-08-28 | 2020-01-14 | 中南大学 | 一种预测变应变速率工况下金属或合金材料动态再结晶晶粒尺寸的方法 |
CN110751991B (zh) * | 2019-11-20 | 2022-09-23 | 中南大学 | 一种预测时变工况下含Nb镍基合金δ相溶解分数的方法 |
CN110793988B (zh) * | 2019-11-21 | 2022-03-15 | 山东建筑大学 | 一种用高能x射线表征非晶-晶体转变动力学特性的方法 |
CN112129692A (zh) * | 2020-10-09 | 2020-12-25 | 中车青岛四方机车车辆股份有限公司 | 一种补偿材料模拟老化程度偏差的方法及系统 |
CN112210731A (zh) * | 2020-10-16 | 2021-01-12 | 西北工业大学 | 一种高熵合金非均匀组织的调控方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102323981B (zh) * | 2011-08-26 | 2014-04-02 | 首钢总公司 | 一种预测热轧钢材奥氏体静态再结晶组织演变的方法 |
CN102519801B (zh) * | 2011-12-16 | 2015-07-29 | 南京航空航天大学 | 用真应力-真应变曲线建立铝合金动态再结晶模型的方法 |
CN105631183B (zh) * | 2014-11-10 | 2018-12-18 | 沈阳黎明航空发动机(集团)有限责任公司 | 一种建立tc6钛合金锻造成形微观组织预测模型的方法 |
CN104372277B (zh) * | 2014-11-28 | 2016-10-05 | 中国科学院金属研究所 | 一种颗粒状δ相均匀分布的GH4169合金制备方法 |
CN105373683B (zh) * | 2015-12-11 | 2018-09-14 | 武汉理工大学 | 一种20CrMnTiH钢热变形过程微观组织演化规律的预测方法 |
CN106202675B (zh) * | 2016-07-04 | 2019-06-21 | 西北工业大学 | 预测钛合金等温成形与动态再结晶演化耦合响应的方法 |
CN106503397B (zh) * | 2016-11-16 | 2019-06-04 | 中国航空工业集团公司北京航空材料研究院 | 一种获得金属材料再结晶组织演化晶界可动性参数的方法 |
CN107058923B (zh) * | 2017-01-16 | 2018-12-07 | 青岛理工大学 | 一种镁合金热加工工艺优化新方法 |
CN107423469B (zh) * | 2017-04-21 | 2020-06-26 | 太原科技大学 | 一种06Cr19Ni9NbN钢锻透的判定方法 |
-
2018
- 2018-01-22 CN CN201810059920.9A patent/CN110068507B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110068507A (zh) | 2019-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110068507B (zh) | 一种对传统再结晶模型进行修正的方法 | |
CN110261247B (zh) | 金属材料各向异性屈服及硬化本构参数同步表征方法 | |
CN108256179B (zh) | 一种预测材料蠕变曲线的方法 | |
CN103886125B (zh) | 一种钛合金热复合成形数值模拟方法 | |
KR20080002410A (ko) | 인장시험과 유한요소법을 이용한 고 변형률에 대한 진변형률-진응력 곡선의 획득 방법 및 이를 이용한 인장 시험기 | |
CN109933925B (zh) | 一种金属板材的冲压成形性能预测方法 | |
CN108133082B (zh) | 一种基于有限元模拟确定压痕应变法中应力测量常数的方法 | |
CN101852735B (zh) | 用光电直读发射光谱仪测定铝钛硼合金中钛含量的方法 | |
CN111207884B (zh) | 一种残余应力挠度法校准装置的标样的设计制作方法 | |
CN101813620A (zh) | 测定铝钛硼合金中硼元素含量的方法 | |
Spitzer et al. | Improved reliability of pH measurements | |
CN103632013B (zh) | 基于主成分回归分析的晶粒尺寸影响无取向硅钢磁性能的分析方法 | |
CN103678921B (zh) | 一种基于主成分回归分析的织构影响无取向硅钢磁性能的分析方法 | |
CN107818184B (zh) | 一种构建材料变形抗力模型的方法 | |
CN108169019B (zh) | 一种准静态塑性压缩应力应变参量的识别方法 | |
CN104764659A (zh) | 一种半固态坯料多向压缩试验方法 | |
CN109142358B (zh) | 基于神经网络的金属材料力学性能测试方法 | |
CN103115870A (zh) | 一种铁基非晶光谱标样的实验室制作方法 | |
CN108169013B (zh) | 一种适用于多钉连接试件的疲劳寿命校准方法 | |
CN108982564A (zh) | 一种水渣xrf定量分析的方法 | |
CN110567808B (zh) | 超高性能混凝土的抗拉强度和抗折强度的评价方法 | |
CN109870257A (zh) | 一种板材厚度方向淬火残余应力分布预测方法 | |
CN110018046B (zh) | 一种节约型双相不锈钢中trip效应致塑性增量的表征方法 | |
CN109086529B (zh) | 一种基于零压力下应变增量确定压痕应变法中应力计算函数的方法 | |
CN106383946A (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 |