CN102234716A - 纯铜热变形过程中动态再结晶组织演变的预测方法 - Google Patents
纯铜热变形过程中动态再结晶组织演变的预测方法 Download PDFInfo
- Publication number
- CN102234716A CN102234716A CN2011102201340A CN201110220134A CN102234716A CN 102234716 A CN102234716 A CN 102234716A CN 2011102201340 A CN2011102201340 A CN 2011102201340A CN 201110220134 A CN201110220134 A CN 201110220134A CN 102234716 A CN102234716 A CN 102234716A
- Authority
- CN
- China
- Prior art keywords
- grain
- calculate
- dynamic
- dynamic recrystallization
- formula
- 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.)
- Pending
Links
- JAPMJSVZDUYFKL-LWOQYNTDSA-N C1C2[C@@H]1CCC2 Chemical compound C1C2[C@@H]1CCC2 JAPMJSVZDUYFKL-LWOQYNTDSA-N 0.000 description 1
Images
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Description
技术领域
本发明涉及一种纯铜热变形过程中动态再结晶组织演变的预测方法,属于金属轧制技术领域。
背景技术
金属材料在热轧过程中存在着复杂的组织演变过程,常常伴随着动态回复以及动态再结晶过程。在多阶段轧制过程中,发生静态回复现象,对最终的晶粒尺寸产生重要影响。
目前研究金属材料热轧过程主要依靠实验手段,不仅耗费大量的人力、物力、财力,而且无法实时观测组织的变化过程。随着计算机技术的迅速发展,数值模拟已成为材料研究中的一种重要手段,通过建立适当物理冶金模型,并结合相关计算机模拟技术,可以实时、定量的再现热加工过程中的显微组织演化过程,大大节约研发成本,缩短新产品的开发周期,对于合理制定加工工艺和优化制品的组织和性能具有重要意义。采用物理冶金理论与多相场法结合的方式对热变形过程中的动态再结晶过程进行预测,能实现对热变形过程中晶粒组织的实时观测,并得到应力应变曲线,从而得到流变应力的变化规律,以及晶粒尺寸随应变变化的转变规律等重要信息,对进一步分析微观组织的演变规律有重要意义。因此开发出用来预测纯铜热变形中动态再结晶组织演变的方法具有重要的意义。
有关专利报道(授权号101591729B),根据元胞自动机法对板带钢再结晶组织演变的预测,所构造的模型简单,选取参数缺乏一定的物理意义,实施例无具体参数,可信度不高。
发明内容
本发明的目的在于提供一种纯铜热变形过程中动态再结晶组织演变的预测方法,该方法具有过程简单、准确和可信度高特点。
本发明是通过以下技术方案加以实现的,一种纯铜热变形过程中动态再结晶组织演变的预测方法,该方法以相场理论为基础,应用相场变量演化方程:
通过计算机对上述模型模拟计算,对纯铜在动态再结晶过程中的晶粒形态、晶粒尺寸、平均晶粒半径、再结晶动力学以及流变应力进行预测,其特征在于包括以下过程:
(1)通过程序编制,建立二维相场空间,输入初始参数:
将计算模型的模拟区域划分为m×m1个二维空间,每个格点ΔX=ΔY=0.5μm,每个模拟区域的试样尺寸为n×n1mm,计算相场演变时用到的参数:界面能γ=0.208J/m2、系数M0=0.139m4K/Js、系数μ=42.1GPa、激活能Qb=110KJ/mol、温度T=750~950K、应变速率计算形核率时用到的参数:激活能Qa=275KJ/mol、系数c=5.0×1025、系数d=1,将它们输入计算程序中;
(2)采用周期性边界条件,通过对计算模拟模型式(1)的初步计算,得到母相基体初始晶粒,初始晶粒设置为规则的正六边形,用灰色表示;新形成的无规则的再结晶晶粒用白色表示;
(3)给定初始位错密度:
计算模型式(1)赋予每个晶粒2个变量:其中一个为位错密度变量,它的初始位错密度为ρini=1.0×109/m2;另一个为晶粒相场变量,该晶粒相场变量用两个数组P、Q表示,其中P数组存放每个晶粒的相场值,Q数组存放每一个晶粒的编号;
(4)计算位错密度随热压缩模拟区域变形增加的变化关系:
根据Mecking和Kocks提出的位错密度模型式(3),计算位错密度随热压缩模拟区域变形增加的变化关系:
式(3)中等式右端第一项表示加工硬化项,k1为硬化系数,k1=4×108,第二项表示软化项,代表回复过程,k2为软化系数,是温度和应变速率的函数,k2=32;
(5)根据位错密度对每个再结晶晶粒判断形核条件:
(6)计算位错密度随动态再结晶的变化的流变应力:
根据式(4)计算位错密度随动态再结晶的变化的流变应力,纯铜应力与平均位错密度满足:
其中:α为0.5,μ是剪切模量,b是柏格斯矢量,是平均位错密度;
(7)输出动态再结晶晶粒组织的动态演化图像、再结晶过程中平均晶粒尺寸以及动态再结晶分数随应变的变化关系:
经S时间段计算,输出动态再结晶晶粒组织的动态演化图像、再结晶过程中平均晶粒尺寸变化关系,再依据动态再结晶分数计算式(5),计算输出动态再结晶分数随应变的变化关系:
X=Nnew/N (5)
式中,Nnew为新生成的再结晶晶粒所占的格点数目,N为原始母相基体所占的格点数目,根据输出结果与设计精度进行比较,调整S时间段,重新开始计算,直至达到设计精度要求为止。
本发明具有以下优点:本发明实现对纯铜热变形过程中动态再结晶组织的演化过程的预测,得到动态再结晶过程中应力应变曲线变化形式、平均晶粒尺寸的变化关系、动态再结晶动力学的转变关系。实现了计算机模拟对纯铜热变形过程中组织演变的再现,不仅能够节约实验成本,而且有助于热加工新工艺的制定。所开发的用来预测纯铜热变形过程中动态再结晶组织演变的方法,实现了在动态再结晶过程中,晶粒形态,晶粒尺寸、平均晶粒半径,再结晶转变分数以及流变应力的预测,对于合理制定加工工艺和优化制品的组织和性能具有重要意义。
附图说明
图1为本发明预测纯铜热变形过程中动态再结晶组织演变方法具体的计算模拟流程框图。
图2~图5为本发明实施例一得到的动态再结晶组织演变过程结果输出图。
图6为本发明实施例二得到的应力应变曲线输出图。
图7为本发明实施例三得到的动态再结晶动力学曲线输出图。
图8为本发明实施例四得到的平均晶粒尺寸变化曲线的输出图。
图9为本发明实施例五得到的平均晶粒尺寸变化曲线的输出图。
具体实施方式
实施例一
将式(1)模型的模拟区域划分为1045×1101的二维空间,整个模拟区域代表0.5225mm×0.5505mm的试样尺寸。选取参数界面能γ=0.208J/m2、M0=0.139m4K/Js、μ=42.1GPa、激活能Qa=275KJ/mol、激活能Qb=110KJ/mol、c=5.0×1025、d=1、温度T=800k、应变速率然后通过程序计算得到平均半径为13.56μm的规则的正六边形的初始母相基体晶粒组织。对每一个晶粒给定初始位错密度ρini=1.0×109/m2,利用一个数组P存储每一个晶粒相场值,利用另一个数组Q存储每一个晶粒的编号。初始状态给定后,开始施加应力变形,计算总步数为19200步,利用式(3)计算位错密度的变化,当位错密度达到临界值ρc=5.5×1013/m2时,开始在晶界上随机形核,利用式(1)计算再结晶晶粒的相场变化,当应变ε分别为0.0702、0.0832、0.0988和0.5000时输出计算结果,如图2~图5所示,当计算步数为第19200步时停止计算。
实施例二
将式(1)模型的模拟区域划分为913×769的二维空间,整个模拟区域代表0.4565mm×0.3845mm的试样尺寸。选取参数界面能γ=0.208J/m2、M0=0.139m4K/Js、μ=42.1GPa、激活能Qa=275KJ/mol、激活能Qb=110KJ/mol、c=5.0×1025、d=1、应变速率变形温度分别选择为750K、800K、850K和950K。然后通过程序计算得到平均半径116.62μm的规则的正六边形的初始母相基体晶粒组织。对每一个晶粒给定初始位错密度ρini=1.0×109/m2,利用一个数组P存储每一个晶粒相场值,利用另一个数组Q存储每一个晶粒的编号。初始状态给定后,开始施加应力变形,计算总步数为19200步,利用式(3)计算位错密度的变化,当位错密度达到临界值ρc=5.5×1013/m2时,开始在晶界上随机形核,利用式(1)计算再结晶晶粒的相场变化,利用式(4)计算流变应力随应变的变化,每计算100步统计流变应力随应变比的变化量,当计算步数为第19200步时停止计算,输出应力应变曲线,如图6所示。
实施例三
将式(1)模型的模拟区域划分别划分为913×769、1065×1025、1065×897和1045×1101二维空间,整个模拟区域代表0.4565mm×0.3845mm、0.5325×0.5125mm、0.5325mm×0.4485mm和0.5225mm×0.5505mm试样尺寸。选取参数界面能γ=0.208J/m2、M0=0.139m4K/Js、μ=42.1GPa、激活能Qa=275KJ/mol、激活能Qb=110KJ/mol、c=5.0×1025、d=1、应变速率温度T=800K。然后通过程序计算得到平均半径为13.56μm、37.76μm、67.34μm和116.62μm的规则的正六边形的初始母相基体晶粒组织。对每一个晶粒给定初始位错密度ρini=1.0×109/m2,利用一个数组P存储每一个晶粒相场值,利用另一个数组Q存储每一个晶粒的编号。初始状态给定后,开始施加应力变形,计算总步数为19200步,利用式(3)计算位错密度的变化,当位错密度达到临界值ρc=5.5×1013/m2时,开始在晶界上随机形核,利用式(1)计算再结晶晶粒的相场变化,利用式(5)计算动态再结晶分数随应变的变化,每计算100步统计再结晶分数随应变的变化量,当计算步数为第19200步时停止计算,输出动态再结晶动力学曲线,如图7所示。
实施例四
将式(1)模型的模拟区域划分别划分为1065×897二维空间,整个模拟区域代表0.5325mm×0.4485mm试样尺寸。选取参数界面能γ=0.208J/m2、M0=0.139m4K/Js、μ=42.1Gpa、激活能Qa=275KJ/mol、激活能Qb=110KJ/mol、c=5.0×1025、d=1、应变速率温度T=800K。然后通过程序计算得到平均半径为67.34μm的规则的正六边形的初始母相基体晶粒组织。对每一个晶粒给定初始位错密度ρini=1.0×109/m2,利用一个数组P存储每一个晶粒相场值,利用另一个数组Q存储每一个晶粒的编号。初始状态给定后,开始施加应力变形,计算总步数为19200步,利用式(3)计算位错密度的变化,当位错密度达到临界值ρc=5.5×1013/m2时,开始在晶界上随机形核,利用式(1)计算再结晶晶粒的相场变化,每计算3000步统计平均晶粒尺寸随应变的变化量,当计算步数为第19200步时停止计算,输出平均晶粒半径随应变的变化曲线,如图8所示,由于本实施例选取S=3000远远大于其他实施例中的输出步数S=100,因此输出结果精度较差,不符合实际规律,为了提高输出结果精度,将输出步数S值逐步减小,从而提高输出精度,直至达到设计要求。
实施例五
将式(1)模型的模拟区域划分别划分为1065×897二维空间,整个模拟区域代表0.5325mm×0.4485mm试样尺寸。选取参数界面能γ=0.208J/m2、M0=0.139m4K/Js、μ=42.1Gpa、激活能Qa=275KJ/mol、激活能Qb=110KJ/mol、c=5.0×1025、d=1、应变速率温度T=800K。然后通过程序计算得到平均半径为67.34μm的规则的正六边形的初始母相基体晶粒组织。对每一个晶粒给定初始位错密度ρini=1.0×109/m2,利用一个数组P存储每一个晶粒相场值,利用另一个数组Q存储每一个晶粒的编号。初始状态给定后,开始施加应力变形,计算总步数为19200步,利用式(3)计算位错密度的变化,当位错密度达到临界值ρc=5.5×1013/m2时,开始在晶界上随机形核,利用式(1)计算再结晶晶粒的相场变化,每计算100步统计平均晶粒尺寸随应变的变化量,当计算步数为第19200步时停止计算,输出平均晶粒半径随应变的变化曲线,如图9所示。
Claims (1)
1.一种纯铜热变形过程中动态再结晶组织演变的预测方法,该方法以相场理论为基础,应用相场变量演化方程:
通过计算机对上述模型模拟计算,对纯铜在动态再结晶过程中的晶粒形态、晶粒尺寸、平均晶粒半径、再结晶动力学以及流变应力进行预测,其特征在于包括以下过程:
(1)通过程序编制,建立二维相场空间,输入初始参数:
将计算模型的模拟区域划分为m×m1个二维空间,每个格点ΔX=ΔY=0.5μm,每个模拟区域的试样尺寸为n×n1mm,计算相场演变时用到的参数:界面能γ=0.208J/m2、系数M0=0.139m4K/Js、系数μ=42.1GPa、激活能Qb=110KJ/mol、温度T=750~950K、应变速率计算形核率时用到的参数:激活能Qa=275KJ/mol、系数c=5.0×1025、系数d=1,将它们输入计算程序中;
(2)采用周期性边界条件,通过对计算模拟模型式(1)的初步计算,得到母相基体初始晶粒,初始晶粒设置为规则的正六边形,用灰色表示;新形成的无规则的再结晶晶粒用白色表示;
(3)给定初始位错密度:
计算模型式(1)赋予每个晶粒2个变量:其中一个为位错密度变量,它的初始位错密度为ρini=1.0×109/m2;另一个为晶粒相场变量,该晶粒相场变量用两个数组P、Q表示,其中P数组存放每个晶粒的相场值,Q数组存放每一个晶粒的编号;
(4)计算位错密度随热压缩模拟区域变形增加的变化关系:
根据Mecking和Kocks提出的位错密度模型式(3),计算位错密度随热压缩模拟区域变形增加的变化关系:
式(3)中等式右端第一项表示加工硬化项,k1为硬化系数,k1=4×108,第二项表示软化项,代表回复过程,k2为软化系数,是温度和应变速率的函数,k2=32;
(5)根据位错密度对每个再结晶晶粒判断形核条件:
(6)计算位错密度随动态再结晶的变化的流变应力:
根据式(4)计算位错密度随动态再结晶的变化的流变应力,纯铜应力与平均位错密度满足:
(7)输出动态再结晶晶粒组织的动态演化图像、再结晶过程中平均晶粒尺寸以及动态再结晶分数随应变的变化关系:
经S时间段计算,输出动态再结晶晶粒组织的动态演化图像、再结晶过程中平均晶粒尺寸变化关系,再依据动态再结晶分数计算式(5),计算输出动态再结晶分数随应变的变化关系:
X=Nnew/N (5)
式中,Nnew为新生成的再结晶晶粒所占的格点数目,N为原始母相基体所占的格点数目,根据输出结果与设计精度进行比较,调整S时间段,重新开始计算,直至达到设计精度要求为止。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011102201340A CN102234716A (zh) | 2011-08-03 | 2011-08-03 | 纯铜热变形过程中动态再结晶组织演变的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011102201340A CN102234716A (zh) | 2011-08-03 | 2011-08-03 | 纯铜热变形过程中动态再结晶组织演变的预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102234716A true CN102234716A (zh) | 2011-11-09 |
Family
ID=44885835
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011102201340A Pending CN102234716A (zh) | 2011-08-03 | 2011-08-03 | 纯铜热变形过程中动态再结晶组织演变的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102234716A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104602830A (zh) * | 2012-08-31 | 2015-05-06 | 东芝三菱电机产业系统株式会社 | 材料组织预测装置、产品制造方法及材料组织预测方法 |
CN106053754A (zh) * | 2016-07-06 | 2016-10-26 | 中南大学 | 一种预测时变工况下高合金化材料动态再结晶分数的方法 |
CN108681629A (zh) * | 2018-05-03 | 2018-10-19 | 中国科学院金属研究所 | 一种冷拉拔工业纯铝导线的屈服强度预测方法 |
CN107058923B (zh) * | 2017-01-16 | 2018-12-07 | 青岛理工大学 | 一种镁合金热加工工艺优化新方法 |
CN109753691A (zh) * | 2018-12-11 | 2019-05-14 | 西安交通大学 | 针对钠冷快堆单排组件耦合热变形的模拟方法 |
-
2011
- 2011-08-03 CN CN2011102201340A patent/CN102234716A/zh active Pending
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104602830A (zh) * | 2012-08-31 | 2015-05-06 | 东芝三菱电机产业系统株式会社 | 材料组织预测装置、产品制造方法及材料组织预测方法 |
CN104602830B (zh) * | 2012-08-31 | 2017-03-01 | 东芝三菱电机产业系统株式会社 | 材料组织预测装置、产品制造方法及材料组织预测方法 |
CN106053754A (zh) * | 2016-07-06 | 2016-10-26 | 中南大学 | 一种预测时变工况下高合金化材料动态再结晶分数的方法 |
CN106053754B (zh) * | 2016-07-06 | 2017-11-14 | 中南大学 | 一种预测时变工况下高合金化材料动态再结晶分数的方法 |
CN107058923B (zh) * | 2017-01-16 | 2018-12-07 | 青岛理工大学 | 一种镁合金热加工工艺优化新方法 |
CN108681629A (zh) * | 2018-05-03 | 2018-10-19 | 中国科学院金属研究所 | 一种冷拉拔工业纯铝导线的屈服强度预测方法 |
CN108681629B (zh) * | 2018-05-03 | 2023-12-15 | 中国科学院金属研究所 | 一种冷拉拔工业纯铝导线的屈服强度预测方法 |
CN109753691A (zh) * | 2018-12-11 | 2019-05-14 | 西安交通大学 | 针对钠冷快堆单排组件耦合热变形的模拟方法 |
CN109753691B (zh) * | 2018-12-11 | 2019-10-29 | 西安交通大学 | 针对钠冷快堆单排组件耦合热变形的模拟方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101591729B (zh) | 预测板带钢热变形中奥氏体动态再结晶组织演变的方法 | |
CN102234716A (zh) | 纯铜热变形过程中动态再结晶组织演变的预测方法 | |
Takaki et al. | Multiscale modeling of hot-working with dynamic recrystallization by coupling microstructure evolution and macroscopic mechanical behavior | |
Chen et al. | Modeling the dynamic recrystallization in austenitic stainless steel using cellular automaton method | |
Kalidindi et al. | On the prediction of yield surfaces by the crystal plasticity models for fcc polycrystals | |
Jin et al. | Investigation on dynamic recrystallization using a modified cellular automaton | |
Derazkola et al. | Review on modeling and simulation of dynamic recrystallization of martensitic stainless steels during bulk hot deformation | |
CN102254057A (zh) | 薄板轧制离线力学性能预报方法 | |
CN112131633B (zh) | 一种基于粗粒化计算理论的流固耦合仿真方法及系统 | |
Kim et al. | Microstructural modelling of equal channel angular pressing for producing ultrafine grained materials | |
CN102339344B (zh) | 一种动态再结晶模型参数反分析识别方法 | |
JIN et al. | Identification of nucleation parameter for cellular automaton model of dynamic recrystallization | |
Zhang et al. | Three-dimensional virtual grain structure generation with grain size control | |
Frydrych et al. | Microstructure evolution in cold-rolled pure titanium: modeling by the three-scale crystal plasticity approach accounting for twinning | |
Guan-feng et al. | Constitutive model of 25CrMo4 steel based on IPSO-SVR and its application in finite element simulation | |
Jin et al. | Modelling the effect of initial grain size on dynamic recrystallization using a modified cellular automata and a adaptive response surface method | |
Sitko et al. | The role of the cellular automata cell size and time step length in the microstructure evolution model—The static recrystallization case study | |
Shi et al. | Non-parametric shape design optimization of elastic-plastic shear panel dampers under cyclic loading | |
Su et al. | Comparative investigation of shear-band evolution using discrete and continuum-based particle methods | |
Yang et al. | Inhomogeneous Deformation and Texture Evolution of 30SiMn2MoVA Steel Gun Barrel Processed by Radial Forging Based on Cross-scale Crystal Plasticity Finite Element Method | |
He et al. | Robust parameters control methodology of microstructure for heavy forgings based on Taguchi method | |
CN105740513A (zh) | 一种gh4169合金热变形动态再结晶模拟方法 | |
Bugat et al. | The toughness module of the PERFECT platform: A predictive tool for the fracture toughness of RPV steels | |
Esche | Monte Carlo simulations of grain growth in metals | |
Smolyakov et al. | Three-dimensional numerical simulations of multi-pass equal-channel angular pressing by a variation difference method and comparison with experiment |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20111109 |