CN105586476B - 基于数值模拟技术的高强钢淬火工艺优化方法 - Google Patents
基于数值模拟技术的高强钢淬火工艺优化方法 Download PDFInfo
- Publication number
- CN105586476B CN105586476B CN201610032619.XA CN201610032619A CN105586476B CN 105586476 B CN105586476 B CN 105586476B CN 201610032619 A CN201610032619 A CN 201610032619A CN 105586476 B CN105586476 B CN 105586476B
- Authority
- CN
- China
- Prior art keywords
- workpiece
- temperature
- transformation
- numerical simulation
- finite element
- 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.)
- Expired - Fee Related
Links
Classifications
-
- C—CHEMISTRY; METALLURGY
- C21—METALLURGY OF IRON
- C21D—MODIFYING THE PHYSICAL STRUCTURE OF FERROUS METALS; GENERAL DEVICES FOR HEAT TREATMENT OF FERROUS OR NON-FERROUS METALS OR ALLOYS; MAKING METAL MALLEABLE, e.g. BY DECARBURISATION OR TEMPERING
- C21D6/00—Heat treatment of ferrous alloys
-
- C—CHEMISTRY; METALLURGY
- C21—METALLURGY OF IRON
- C21D—MODIFYING THE PHYSICAL STRUCTURE OF FERROUS METALS; GENERAL DEVICES FOR HEAT TREATMENT OF FERROUS OR NON-FERROUS METALS OR ALLOYS; MAKING METAL MALLEABLE, e.g. BY DECARBURISATION OR TEMPERING
- C21D1/00—General methods or devices for heat treatment, e.g. annealing, hardening, quenching or tempering
- C21D1/18—Hardening; Quenching with or without subsequent tempering
Landscapes
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Thermal Sciences (AREA)
- Crystallography & Structural Chemistry (AREA)
- Mechanical Engineering (AREA)
- Materials Engineering (AREA)
- Metallurgy (AREA)
- Organic Chemistry (AREA)
- Control Of Heat Treatment Processes (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Investigating And Analyzing Materials By Characteristic Methods (AREA)
Abstract
一种基于数值模拟技术的高强钢淬火工艺优化方法,根据待处理材料的热物理性能参数、相变动力学参数以及冷却介质的换热系数作为初始条件和边界定义,利用有限元计算方法得出工件的温度场与相变场分布,并根据工件厚度类型确定预冷时间从而实现淬火工艺优化。本发明可以快速有效地确定高强钢淬火过程的预冷时间,提高其淬火工艺的可靠性与稳定性。
Description
技术领域
本发明涉及的是一种钢铁行业热处理领域的技术,具体是一种基于数值模拟技术的高强钢淬火工艺优化方法。
背景技术
传统工艺对预冷时间的确定往往是基于经验公式或实验测定。经验公式往往不针对具体材料与工件外形尺寸,对预冷时间的确定存在较大的偏差;而通过实测实验进行组织与变形验证效率较低,从而增加生产成本。另外对于不同厚度的工件,还没有一个标准的确定预冷时间的准则。
预冷淬火是其预防裂纹和减小变形的方法之一,适用于淬透性好、奥氏体较稳定、临界冷却速率较低的钢材及厚薄悬殊、形状复杂且对变形控制要求较高的工件。对于奥氏体具有较长亚稳定区间的钢材,采用预冷淬火工艺,可以有效的减少淬火畸变的产生。传统工艺对预冷时间的确定往往是基于经验公式或实验测定。经验公式往往不针对具体材料与工件外形尺寸,对预冷时间的确定存在较大的偏差;而通过实测实验进行组织与变形验证效率较低,从而增加生产成本。另外对于不同厚度的工件,还没有一个标准的确定预冷时间的准则。
工件热处理后的性能和变形状况主要取决于零件内部各部位的微观组织构成和表征以及残余应力/应变的分布。这些特征是热处理过程零件内部温度、组织、应力应变状态共同耦合作用和演化的结果,是一个极为复杂的过程。这些特征的变化规律是复杂的,且难于通过有限数量的工艺试验来充分掌握。
热处理过程的数值模拟技术,能够借助少量的试验测试,实现工件热处理性能和变形的预测,已经成为热处理工艺研究中一项重要的辅助手段,近年来得以快速发展和广泛应用。通过数值模拟结合少量实物验证的方式可以较为快速的获取优化的热处理工艺,这样的工艺研究方式已得到广泛的验证。
在零件热处理(淬火冷却)过程中,温度与组织存在相互影响的耦合作用,热处理中淬火过程的温度场可以由经典的热传导方程来描述: 其中:T为温度;λx为x方向上的热导率;Q为内热源;ρ为密度;Cp为比热;τ为时间。
当发生组织相变时,会引发相变潜热。可以将此潜热看作为一种内热源加以处理,即有:其中:ΔH:单位潜热;单位时间内形成的新相体积分数。
发明内容
本发明针对现有技术存在的上述不足,提出一种基于数值模拟技术的高强钢淬火工艺优化方法,利用有限元计算方法(FEM)计算得出工件的温度场与相变场分布,根据工件厚度类型确定预冷时间,通过采用计算机数值模拟技术对热处理工艺过程的温度场、应力应变场与组织场进行准确的预测,结合工件壁厚优化原则,可以快速有效地确定高强钢淬火过程的预冷时间,提高其淬火工艺的可靠性与稳定性。
本发明是通过以下技术方案实现的:
本发明根据待处理材料的热物理性能参数、相变动力学参数以及冷却介质的换热系数作为初始条件和边界定义,利用有限元计算方法得出工件的温度场与相变场分布,并根据工件厚度类型确定预冷时间从而实现淬火工艺优化。
所述的热物理性能参数包括:密度、比热容、热导率。
所述的相变动力学参数是指:等温转变动力学参数。
所述的有限元计算方法,采用但不限于有限元软件MSC.Marc实现。
所述的有限元计算方法,通过对工件淬火过程进行数值模拟分析,获取工件在热处理过程中各个部位的温度场与相变场分布数据,具体包括工件在冷却过程中其组织转变模型中冷却最快与最慢节点的温度与相变量历史曲线;
所述的组织转变模型是通过过冷奥氏体等温转变动力学曲线得出的,对于扩散型相变其转变量计算基于扩展的Avrami方程:f=1-exp(-bτn),其中f表示新相的体积分数,b,n表示相变动力学常数,τ表示等温时间;对于非扩散型马氏体相变,其转变量计算基于经典的K-M方程:fM=1-exp[-αM·(Ms-T)],其中:fM表示马氏体体积分数,αM表示相变常数,Ms表示马氏体相变起始温度,T表示当前温度。
所述的预冷时间是指:当工件属于薄壁件,其淬火时预冷时间根据工件中冷却最慢的点的温度不低于铁素体转变开始温度;当工件属于厚壁件,其淬火时预冷时间根据工件中冷却最快的点的温度不低于铁素体转变开始温度。
技术效果
与现有技术相比,本发明借助有限元软件、热力学计算软件与常规实验,能有效确定复杂结构不同壁厚高强钢工件热处理淬火工艺中的预冷时间;通过本发明方法的应用,通过采用计算机数值模拟技术对热处理工艺过程的温度场、应力应变场与组织场进行准确的预测,结合工件壁厚优化原则,可以快速有效地确定高强钢淬火过程的预冷时间,提高其淬火工艺的可靠性与稳定性。
附图说明
图1为实施例中计算得出的30CrMnSiNi2A的热导率示意图。
图2为实施例中30CrMnSiNi2A先共析铁素体/珠光体等温转变曲线示意图。
图3为实施例中工件模型中冷却最快与最慢节点的温度变化曲线示意图。
具体实施方式
本实施例以30CrMnSiNi2A材质的轴承基座淬火工艺为例,该工件构形复杂,截面尺寸差异大,厚薄不均匀,最厚处达70mm,最薄处仅为4mm,热处理淬火时不同部位受热和冷却程度差异大,容易产生变形,本实施例包括以下步骤:
步骤1)确定材料的热物理性能参数
借助材料热力学计算软件JMatPro对30CrMnSiNi2A材料的各项热物理性能参数进行计算,包括密度、比热容、热导率,由于这些参数是与温度相关的,实际计算出的一条性能参数与温度的曲线,如图1列出了计算得出的热导率与温度的关系图;
步骤2)确定材料的相变动力学参数
通过Gleeble3500热物理模拟试验机对30CrMnSiNi2A材料的等温转变动力学进行测量,对应的相变动力学参数是通过测试得出的试样在不同温度下等温过程中的热膨胀曲线的变化特征得出的,测试得出的30CrMnSiNi2A先共析铁素体/珠光体等温转变曲线如图2所示;
步骤3)确定冷却介质的换热系数
本实施例预冷在空气中进行,空冷的换热系数采用如下经验公式计算得出: 其中:Tw为工件温度,Tc为环境温度,均以K为单位;
步骤4)利用有限元计算方法计算得出工件的温度场与相变场分布
应用有限元软件MSC.Marc以及子程序二次开发对工件淬火过程进行数值模拟分析,获取工件在热处理过程中各个部位的温度场、与相变场分布数据,数据包括工件在冷却过程中其网格模型中冷却最快与最慢节点的温度与相变量历史曲线,计算得出的温度变化曲线如图3所示:
步骤5)根据工件厚度类型确定预冷时间
本实施例的工件属于薄壁件,其预冷时间根据工件中最慢的位置的温度大于铁素体转变开始温度的评判准则进行确定,本实施例材料30CrMnSiNi2A的铁素体转变开始温度约为720℃,根据图3的计算结果,在80s的空气预冷下,工件中冷却最慢的位置的温度,约为850℃,满足评判要求。
上述具体实施可由本领域技术人员在不背离本发明原理和宗旨的前提下以不同的方式对其进行局部调整,本发明的保护范围以权利要求书为准且不由上述具体实施所限,在其范围内的各个实现方案均受本发明之约束。
Claims (3)
1.一种基于数值模拟技术的30CrMnSiNi2A材质的轴承基座的淬火工艺优化方法,其特征在于,根据待处理材料的热物理性能参数、相变动力学参数以及冷却介质的换热系数作为初始条件和边界定义,利用有限元计算方法得出工件的温度场与相变场分布,并根据工件厚度类型确定预冷时间从而实现淬火工艺优化;
所述的热物理性能参数包括:密度、比热容、热导率;
所述的相变动力学参数是指:等温转变动力学参数;
所述的有限元计算方法,通过对工件淬火过程进行数值模拟分析,获取工件在热处理过程中各个部位的温度场与相变场分布数据,具体包括工件在冷却过程中其组织转变模型中冷却最快与最慢节点的温度与相变量历史曲线;
所述的预冷时间是指:当工件属于薄壁件,其淬火时预冷时间根据工件中冷却最慢的点的温度不低于铁素体转变开始温度;当工件属于厚壁件,其淬火时预冷时间根据工件中冷却最快的点的温度不低于铁素体转变开始温度。
2.根据权利要求1所述的基于数值模拟技术的30CrMnSiNi2A材质的轴承基座的淬火工艺优化方法,其特征是,所述的有限元计算方法,采用有限元软件MSC.Marc实现。
3.根据权利要求2所述的基于数值模拟技术的30CrMnSiNi2A材质的轴承基座的淬火工艺优化方法,其特征是,所述的组织转变模型是通过过冷奥氏体等温转变动力学曲线得出的,对于扩散型相变其转变量计算基于扩展的Avrami方程:f=1-exp(-bτn),其中f表示新相的体积分数,b,n表示相变动力学常数,τ表示等温时间;对于非扩散型马氏体相变,其转变量计算基于经典的K-M方程:fM=1-exp[-αM·(Ms-T)],其中:fM表示马氏体体积分数,αM表示相变常数,Ms表示马氏体相变起始温度,T表示当前温度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610032619.XA CN105586476B (zh) | 2016-01-19 | 2016-01-19 | 基于数值模拟技术的高强钢淬火工艺优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610032619.XA CN105586476B (zh) | 2016-01-19 | 2016-01-19 | 基于数值模拟技术的高强钢淬火工艺优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105586476A CN105586476A (zh) | 2016-05-18 |
CN105586476B true CN105586476B (zh) | 2018-01-30 |
Family
ID=55926373
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610032619.XA Expired - Fee Related CN105586476B (zh) | 2016-01-19 | 2016-01-19 | 基于数值模拟技术的高强钢淬火工艺优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105586476B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107480328B (zh) * | 2017-07-04 | 2022-09-20 | 山东建筑大学 | 一种基于q&p工艺的碳配分理论计算方法 |
CN109929969B (zh) * | 2019-03-29 | 2019-12-24 | 上海交通大学 | 一种合金钢水淬工艺的优化方法 |
CN110232214B (zh) * | 2019-05-09 | 2020-07-28 | 北京大学 | 一种通过数值模拟的食品预冷性能评估和优化方法 |
CN110438318B (zh) * | 2019-07-22 | 2020-11-24 | 中南大学 | 一种大型立式淬火炉低能耗最速升温控制方法及系统 |
WO2021092739A1 (zh) * | 2019-11-12 | 2021-05-20 | 常德菲尔美化工技术有限公司 | 一种合金材料淬火确定方法 |
CN111044559B (zh) * | 2019-11-26 | 2021-09-14 | 中南大学 | 一种基于合金端淬有限点温度曲线的相变潜热计算方法 |
CN112084603B (zh) * | 2020-09-17 | 2022-03-15 | 泰尔重工股份有限公司 | 重载万向轴叉头淬火加热工艺参数的获取方法 |
CN113591341A (zh) * | 2021-06-25 | 2021-11-02 | 西部超导材料科技股份有限公司 | 一种基于数值模拟的钛合金锻造工艺优化方法 |
CN113627059B (zh) * | 2021-08-03 | 2024-05-14 | 南京钢铁股份有限公司 | 一种考虑相变热的大规格棒材空冷温度场计算方法 |
CN113642218B (zh) * | 2021-08-18 | 2023-12-15 | 南京钢铁股份有限公司 | 一种钢板淬火处理临界加热速率的确定系统及确定方法 |
CN115029544B (zh) * | 2022-06-22 | 2023-08-01 | 钢铁研究总院有限公司 | 基于模拟预测的薄壁圆环类零件的热处理变形控制方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102005036068A1 (de) * | 2005-08-01 | 2007-02-08 | Siemens Ag | Modellierverfahren für den zeitlichen Verlauf des Zustands eines Stahlvolumens durch einen Rechner und hiermit korrespondierende Gegenstände |
CN101921907A (zh) * | 2009-06-12 | 2010-12-22 | 宝山钢铁股份有限公司 | 一种基于模型控制的钢板淬火时的板形修正方法 |
CN103484646A (zh) * | 2013-09-10 | 2014-01-01 | 中钢集团邢台机械轧辊有限公司 | 轧辊喷雾淬火过程中温度和组织在线模拟控制系统及方法 |
CN103544356A (zh) * | 2013-10-30 | 2014-01-29 | 中冶南方(武汉)信息技术工程有限公司 | 一种基于gpu的热处理炉加热模型计算方法 |
CN104046762A (zh) * | 2014-06-26 | 2014-09-17 | 燕山大学 | 一种大直径厚壁压力气瓶的淬火方法 |
CN104331574A (zh) * | 2014-11-19 | 2015-02-04 | 河海大学常州校区 | 基于ansys有限元平台的感应淬火淬硬层深度的预测方法 |
-
2016
- 2016-01-19 CN CN201610032619.XA patent/CN105586476B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102005036068A1 (de) * | 2005-08-01 | 2007-02-08 | Siemens Ag | Modellierverfahren für den zeitlichen Verlauf des Zustands eines Stahlvolumens durch einen Rechner und hiermit korrespondierende Gegenstände |
CN101921907A (zh) * | 2009-06-12 | 2010-12-22 | 宝山钢铁股份有限公司 | 一种基于模型控制的钢板淬火时的板形修正方法 |
CN103484646A (zh) * | 2013-09-10 | 2014-01-01 | 中钢集团邢台机械轧辊有限公司 | 轧辊喷雾淬火过程中温度和组织在线模拟控制系统及方法 |
CN103544356A (zh) * | 2013-10-30 | 2014-01-29 | 中冶南方(武汉)信息技术工程有限公司 | 一种基于gpu的热处理炉加热模型计算方法 |
CN104046762A (zh) * | 2014-06-26 | 2014-09-17 | 燕山大学 | 一种大直径厚壁压力气瓶的淬火方法 |
CN104331574A (zh) * | 2014-11-19 | 2015-02-04 | 河海大学常州校区 | 基于ansys有限元平台的感应淬火淬硬层深度的预测方法 |
Non-Patent Citations (1)
Title |
---|
低压转子淬火工艺过程数值模拟的研究;王晓燕;《中国优秀硕士学位论文全文数据库•工程科技Ⅰ辑》;20110715;正文第17页第3段、第19页第2段、第20页第2段、第27页最后1段、第28页第1段、第29页第2、5段、第50页第2段、第51页第2段、第52页第3段、第53页最后1段、第55页第2段、第57页-第59页第1段和图4-5、4-6 * |
Also Published As
Publication number | Publication date |
---|---|
CN105586476A (zh) | 2016-05-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105586476B (zh) | 基于数值模拟技术的高强钢淬火工艺优化方法 | |
Zhao et al. | Comparison of the methods for calculating the interfacial heat transfer coefficient in hot stamping | |
Jin et al. | Measurement of boiling heat transfer coefficient in liquid nitrogen bath by inverse heat conduction method | |
Sun et al. | Microstructure evolution modeling and simulation for dynamic recrystallization of Cr12MoV die steel during hot compression based on real metallographic image | |
CN102521439A (zh) | 结合有限元法和反传热法计算淬火介质换热系数的方法 | |
Lopez-Garcia et al. | Effect of immersion routes on the quenching distortion of a long steel component using a finite element model | |
Hao et al. | 3-D Numerical analysis on heating process of loads within vacuum heat treatment furnace | |
Jang et al. | Optimisation of a slab heating pattern with various skid button heights in a walking-beam-type reheating furnace | |
Li et al. | Coupled numerical simulation of hot stamping process and experimental verification | |
Wen et al. | Development and application of cast steel numerical simulation system for heat treatment | |
Wei et al. | Evaluation of elastic-plastic deformation in HSS work roll under coupling of residual stress thermal stress and contact stress during hot rolling | |
Pei et al. | Experimental and CFD investigations on cooling process of end-quench test | |
CN113642212B (zh) | 针对大型超高强钢壳体精细化热处理工艺设计方法及系统 | |
Kim et al. | Inverse characterization method for mechanical properties of strain/strain-rate/temperature/temperature-history dependent steel sheets and its application for hot press forming | |
Goldstein et al. | Modeling of Temperature Gradients During Short Time Dilatometry Testing | |
Liu et al. | Optimization of heating process for bearing rings in a vacuum furnace based on numerical analysis | |
Lin et al. | Deformation analysis of hot stamping tools by thermal–fluid-mechanical coupled approach based on MpCCI | |
Dou et al. | Numerical model and optimization strategy for the annealing process of 3D coil cores | |
Yang et al. | CFD prediction for the performance of a heat treatment furnace | |
Kang et al. | Characterization of cooling heat transfer for various coolant conditions in warm forging process | |
Shang et al. | Optimization Design of Insert Hot Stamping Die’s Cooling System and Research on the Microstructural Uniformity Control of Martensitic Phase Transitions in Synchronous Quenching Process | |
US11306371B1 (en) | Gas quenching system and method for minimizing distortion of heat treated parts | |
Zhang et al. | Analysis of temperature and stress in differential temperature reverse deep drawing of 5A06 aluminum cup | |
Zhao et al. | Calculation of phase transformation latent heat based on cooling curves in end-quench test and its application in nickel-based superalloy | |
Li et al. | Experimental and numerical study on optimization of heating process for small-sized workpieces in vacuum heat treatment furnace |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180130 Termination date: 20220119 |
|
CF01 | Termination of patent right due to non-payment of annual fee |