CN107245557A - 基于ttt曲线预测轴件感应淬火淬硬层组织和硬度方法 - Google Patents
基于ttt曲线预测轴件感应淬火淬硬层组织和硬度方法 Download PDFInfo
- Publication number
- CN107245557A CN107245557A CN201710473878.0A CN201710473878A CN107245557A CN 107245557 A CN107245557 A CN 107245557A CN 201710473878 A CN201710473878 A CN 201710473878A CN 107245557 A CN107245557 A CN 107245557A
- Authority
- CN
- China
- Prior art keywords
- shaft member
- impewdance
- tissue
- impewdance matching
- hardness
- 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
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
- C21D1/00—General methods or devices for heat treatment, e.g. annealing, hardening, quenching or tempering
- C21D1/55—Hardenability tests, e.g. end-quench tests
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Organic Chemistry (AREA)
- Mechanical Engineering (AREA)
- Materials Engineering (AREA)
- Metallurgy (AREA)
- Crystallography & Structural Chemistry (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Thermal Sciences (AREA)
- Heat Treatment Of Articles (AREA)
Abstract
本发明公开了一种基于TTT曲线预测轴件感应淬火淬硬层组织和硬度方法,其特征是,包括以下步骤:轴件感应加热过程有限元分析,获取其感应加热过程轴件径向的奥氏体化分布;轴件感应淬火冷却过程有限元分析,获取感应淬火过程轴件径向任一位置的冷却情况;根据一定奥氏体化温度下的零件材料TTT曲线,将轴件感应淬火的有限元分析结果与TTT曲线结合,计算轴件感应淬火淬硬层的组织和硬度分布。与现有技术相比,本发明具有可在避免损坏工件的情况下进行感应淬火的淬硬层组织硬度预测,也可根据实际工件形状、尺寸和设备参数改变预测模型,快捷方便预测出感应淬火淬硬层的组织和硬度,且预测结果具有较为准确的优点。
Description
技术领域
本发明涉及一种基于TTT曲线预测轴件感应淬火淬硬层组织和硬度方法。
背景技术
传统材料热处理工艺通过大量试验研究筛选一种较好处理工艺,成本既高又费时。具有高质量、重复性及适应性强等特征的感应淬火热处理是热处理行业中应用最广和发展最快的表面处理工艺之一,通过有限元分析预测工件感应淬火的淬硬层组织和硬度,确定其最佳工艺参数并指导生产,避免测定工件的淬硬层深度时损坏试验工件。
发明内容
为解决现有技术的不足,本发明的目的在于提供一种基于TTT曲线预测轴件感应淬火淬硬层的组织和硬度方法,基于ANSYS有限元平台来模拟感应淬火过程,结合材料TTT曲线预测轴件感应淬火淬硬层的组织和硬度,减小感应淬火试验工作量。
为了实现上述目标,本发明采用如下的技术方案:
一种基于TTT曲线预测轴件感应淬火淬硬层组织和硬度方法,按照下列步骤进行:
步骤1)依据感应淬火实际工作环境构建感应淬火有限元实体模型;
步骤2)模型划分网格单元,定义二维耦合场实体单元PLANE13属性及材料属性;
施加热源载荷、热对流约束及边界条件至二维耦合场实体单元PLANE13各节点上,采用直接法对工件进行电磁热耦合分析获取电磁热耦合计算结果;
步骤3)计算轴零件在材料失磁后的升温速率,由轴件材料连续加热升温速率VAc1/Ac3与Ac1和Ac3变化曲线计算得到轴件的奥氏体化温度Ac1和Ac3;
步骤4)由感应淬火的工艺参数计算得到轴件模拟加热温度随径向分布曲线,由步骤3)中计算所得Ac1和Ac3值确定奥氏体在轴径向分布上含量范围;
步骤5)感应加热结束后对轴件表面施加对流换热系数载荷,模拟感应淬火冷却过程温度场变化;得到轴件径向淬硬层深度范围内任一位置的冷却曲线;
步骤6)查找《钢的奥氏体转变曲线—原理、测试与应用》得到上述奥氏体化温度下对应工件TTT图,采用数值处理方法处理TTT图,根据工件材料组织场动力学数学模型编写感应淬火过程中电磁场、温度场和组织场的演变程序;
步骤7)将步骤5)中的感应淬火冷却曲线与步骤6)中的基于TTT曲线程序联合计算,得到感应淬火后的轴径向上各处奥氏体转变组织及其硬度变化;
步骤8)根据步骤4)计算得到轴件感应淬火淬硬层中奥氏体含量分布,结合步骤7)中的计算结果得到轴件感应淬火后的淬硬层组织及硬度;
步骤9)将步骤8)的感应淬火淬硬层的组织及硬度计算结果与最终试验结果对比分析误差,给出模拟预测的可行性。
进一步地,所述步骤2)中采用直接法对工件进行电磁热耦合分析,得到电磁热耦合计算结果。
进一步地,所述步骤1)中感应淬火有限元实体模型包括材料物性参数、模型网格大小和迭代时间步长;所述迭代时间不长依据经验值选取。
进一步地,所述步骤4)中奥氏体在轴径向分布上含量范围的确定采用线性插值方法。
进一步地,所述骤6)中采用数值处理方法对TTT图进行拟合处理;根据扩散型和非扩散型相变动力学模型编写程序描述感应淬火组织场和温度场的关系,并根据图中的数据信息拟合处理得到材料一定奥氏体化后各等温转变获取转变组织与其硬度的关系。
进一步地,所述步骤6)中TTT图数值化处理采用软件GetDate进行信息采集;
感应淬火过程中电磁场、温度场和组织场的演变程序采用MATLAB软件编写;
采用Origin软件对轴件感应淬火轴径向各处冷却速度所获组织与淬硬层硬度的数学关系拟合出相应公式。
进一步地,所述步骤1)中感应淬火实际工作环境包括感应加热频率、加热时间、电流密度、感应线圈内径及高度。
本发明所达到的有益效果:本方法基于ANSYS有限元平台,通过模拟计算轴件感应淬火淬硬层的组织和硬度,调节感应工艺参数来满足工件的技术指标,可在避免损坏工件的情况下进行感应淬火淬硬层的组织硬度预测,也可根据实际工件形状、尺寸和设备参数改变预测模型,快捷方便预测出轴件感应淬火淬硬层的组织和硬度,且预测结果具有较为准确的优点。
附图说明
图1是轴件感应淬火几何模型示意图;
图2是轴件有限元分析网格模型的示意图;
图3是轴件表面加热温度随时间变化曲线示意图;
图4是轴件连续加热升温速率与Ac1和Ac3变化曲线示意图;
图5是轴件径向温度分布及奥氏体含量分布示意图;
图6是轴件径向不同深度的感应淬火冷却曲线示意图;
图7是880℃温度奥氏体化的45钢TTT图;
图8是880℃温度奥氏体化的45钢TTT图数值化处理结果的示意图;
图9是45钢奥氏体化后经等温转变后组织与硬度的关系示意图;
图10是轴件径向距表面2.21mm处的淬火冷却曲线(等时间间隔处理后)与奥氏体化温度为880℃的45钢TTT图结合情况示意图;
图11是轴件径向距表面2.21mm处感应淬火后的金相组织示意图;
图12是轴件径向HV硬度预测分布和实测硬度分布曲线示意图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
如图1~4所示,利用ANSYS有限元分析软件和MATLAB软件对45钢预测轴件感应淬火淬硬层的组织和硬度,包括以下步骤:
(1)依据工件技术要求确定感应淬火时实际加热条件:感应加热频率f=195kHz、加热时间t=3.8s和电流密度Js=420×106A/m2;工件直径Φ15.2mm,长5mm的轴件;感应线圈内径为Φ20mm,高度为5mm。
对轴对称工件的1/2进行一维建模:感应器以集肤深度进行建模,同时考虑工件和感应圈之间间隙及周围空气场如图1,其中A1区代表工件,A2区为感应圈,A3为空气场。仿真模拟模型如图2所示,工件表面至心部采用等间距0.01mm网格大小的方式;感应圈以电流集肤深度0.016mm建立网格,空气场采用自由划分的方式。
(2)定义二维耦合场实体单元PLANE13属性及材料属性,材料属性包括材料的密度,在不同温度下的磁导率μ、比热容c和电阻率ρ。
进入电磁热耦合求解过程:定义热源载荷、热对流约束及边界条件,包括环境温度、感应加热电流密度、工件热传导系数及边界换热系数,施加热源载荷、热对流约束及边界条件至二维耦合场实体单元PLANE13的各节点上,并采用直接法对工件进行电磁热耦合分析计算,迭代耦合时间步长为0.01s,得到电磁热耦合结果;将耦合计算单元PLANE13转化为温度计算单元PLANE55,读取电磁热耦合计算结果作为载荷施加到各节点上,求解得到不同时刻和不同区域的温度场。
(3)由轴件表面部位的加热曲线如图3,以45钢的失磁点为760℃计算工件在失磁后升温速率,根据轴件材料连续加热奥氏体化相变温度-升温速率曲线如图4,得到轴件的奥氏体化开始温度Ac1和结束温度Ac3分别为760℃和820℃。
(4)根据温度场求解得到的时间t=3.8s时轴件径向温度分布曲线如图4所示,判断轴径向上奥氏体含量分布,具体判断过程如下:Ac3对应径向深度δ100%A为完全奥氏体组织;Ac1对应径向深度δ0%A为开始奥氏体组织;完全奥氏体化和开始奥氏体化中间温度790℃处深度δ50%A为50%奥氏体组织如图5。
(5)加热结束后对轴件表面施加对流换热系数载荷,冷却20s,模拟冷却过程轴件的温度场变化;由轴件径向上不同深度处感应淬火温度变化情况如图6。
(6)查找《钢的奥氏体转变曲线—原理、测试与应用》得到上述奥氏体化温度880℃下对应轴件45钢的TTT图如图7,采用数值处理方法得到的TTT图如图8,根据材料组织场动力学数学模型编写感应淬火过程中电磁场、温度场和组织场的演变程序。并根据图中的数据信息拟合处理得到材料奥氏体化后各等温转变温度下获取组织和其硬度的关系如图9。
(7)将步骤(5)中的感应淬火的冷却过程与数值处理的奥氏体化温度880℃的45钢TTT图结合,如轴件径向距离表面2.21mm处的冷却曲线与TTT结合情况如图10,冷却曲线需要处理成等时间间隔冷却曲线,再采用上述程序计算相同时间间隔内组织的转变量。
(8)根据步骤(4)计算得到感应淬火淬硬层中的奥氏体含量分布,结合步骤(7)中的计算结果按百分比计算得到轴件感应淬火后淬硬层的组织及硬度,如轴件径向距表面2.21mm处的组织含量计算结果为49.6%马氏体(M)组织,2.0%贝氏体(B)组织和48.4%铁素体及珠光体(F+P)组织,实测此处的金相组织图片如图11。
(9)将步骤(8)的感应淬火淬硬层的硬度计算结果与最终实验结果对比分析误差如图12,给出模拟预测的可行性。
图8中所标注的数字表示的是洛氏硬度HRC,字母等采用的是国际通用的定义。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (7)
1.一种基于TTT曲线预测轴件感应淬火淬硬层组织和硬度方法,按照下列步骤进行:
步骤1)依据感应淬火实际工作环境构建感应淬火有限元实体模型;
步骤2)模型划分网格单元,定义二维耦合场实体单元PLANE13属性及材料属性;
施加热源载荷、热对流约束及边界条件至二维耦合场实体单元PLANE13各节点上,采用直接法对工件进行电磁热耦合分析获取电磁热耦合计算结果;
步骤3)计算轴件在材料失磁后的升温速率,由轴件材料连续加热升温速率VAc1/Ac3与Ac1和Ac3变化曲线计算得到轴零件的奥氏体化温度Ac1和Ac3;
步骤4)由感应淬火的工艺参数计算得到轴件模拟加热温度随轴径向分布曲线,由步骤3)中计算所得Ac1和Ac3值确定奥氏体在轴件径向分布上含量范围;
步骤5)感应加热结束后对轴件表面施加对流换热系数载荷,模拟感应淬火冷却过程温度场变化;得到轴件径向淬硬层深度范围内任一位置的冷却曲线;
步骤6)查找《钢的奥氏体转变曲线—原理、测试与应用》得到上述奥氏体化温度下对应工件材料的TTT图,采用数值处理方法处理TTT图,根据工件材料组织场动力学数学模型编写感应淬火过程中电磁场、温度场和组织场的演变程序;
步骤7)将步骤5)中感应淬火的轴件径向各处的冷却曲线与步骤6)中的基于材料TTT曲线程序联合计算,得到感应淬火后轴件径向各处的奥氏体转变组织及其硬度变化;
步骤8)根据步骤4)计算得到的感应淬火淬硬层中奥氏体含量分布,结合步骤7)中的计算结果得到轴件感应淬火后淬硬层的组织及硬度;
步骤9)将步骤8)感应淬火淬硬层的组织及硬度计算结果与最终试验结果对比分析误差,给出模拟预测的可行性。
2.根据权利要求1所述的一种基于TTT曲线预测轴件感应淬火淬硬层组织和硬度方法,其特征是,所述步骤2)中采用直接法对工件进行电磁热耦合分析,得到电磁热耦合计算结果。
3.根据权利要求1所述的一种基于TTT曲线预测轴件感应淬火淬硬层组织和硬度方法,其特征是,所述步骤1)中感应淬火有限元实体模型包括材料物性参数、模型网格大小和迭代时间步长;所述迭代时间步长依据经验值选取。
4.根据权利要求1所述的一种基于TTT曲线预测轴件感应淬火淬硬层组织和硬度方法,其特征是,所述步骤4)中奥氏体在轴径向分布上含量范围的确定采用线性插值方法。
5.根据权利要求1所述的一种基于TTT曲线预测轴件感应淬火淬硬层组织和硬度方法,其特征是,所述骤6)中采用数值处理方法对材料TTT图进行拟合处理;根据扩散型和非扩散型相变动力学模型编写程序描述感应淬火组织场和温度场的关系,并根据图中的数据信息拟合处理得到材料一定奥氏体化温度下各等温转变温度下获取的组织与其硬度的关系。
6.根据权利要求1所述的一种基于TTT曲线预测轴件感应淬火淬硬层组织和硬度方法,其特征是,所述步骤6)中TTT图数值化处理采用软件GetDate进行信息采集;
感应淬火过程中电磁场、温度场和组织场的演变程序采用MATLAB软件编写;
采用Origin软件对感应淬火轴件径向各处冷却速度所获组织与淬硬层硬度的数学关系拟合出相应公式。
7.根据权利要求1所述的一种基于TTT曲线预测轴件感应淬火淬硬层组织和硬度方法,其特征是,所述步骤1)中感应淬火实际工作环境包括感应加热频率、加热时间、电流密度、感应线圈内径及高度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710473878.0A CN107245557A (zh) | 2017-06-21 | 2017-06-21 | 基于ttt曲线预测轴件感应淬火淬硬层组织和硬度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710473878.0A CN107245557A (zh) | 2017-06-21 | 2017-06-21 | 基于ttt曲线预测轴件感应淬火淬硬层组织和硬度方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107245557A true CN107245557A (zh) | 2017-10-13 |
Family
ID=60018428
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710473878.0A Pending CN107245557A (zh) | 2017-06-21 | 2017-06-21 | 基于ttt曲线预测轴件感应淬火淬硬层组织和硬度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107245557A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113174467A (zh) * | 2021-03-23 | 2021-07-27 | 中冶南方连铸技术工程有限责任公司 | 铸坯淬冷组织的预测方法及铸坯淬冷工艺的制定方法 |
CN116694919A (zh) * | 2023-06-12 | 2023-09-05 | 无锡市海峰海林精密轴承有限公司 | 一种轴承贝氏体热处理工艺的优化方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103413059A (zh) * | 2013-08-23 | 2013-11-27 | 山东建筑大学 | 一种固体截面温变测定系统 |
CN104331574A (zh) * | 2014-11-19 | 2015-02-04 | 河海大学常州校区 | 基于ansys有限元平台的感应淬火淬硬层深度的预测方法 |
CN105631111A (zh) * | 2015-12-25 | 2016-06-01 | 河海大学常州校区 | 一种台阶轴感应淬火马氏体分布的预测方法 |
CN105653770A (zh) * | 2015-12-25 | 2016-06-08 | 河海大学常州校区 | 一种连续感应淬火淬硬层深度的预测方法 |
CN105930605A (zh) * | 2016-04-29 | 2016-09-07 | 河海大学常州校区 | 感应淬火处理轴零件的淬硬层深度测量方法 |
-
2017
- 2017-06-21 CN CN201710473878.0A patent/CN107245557A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103413059A (zh) * | 2013-08-23 | 2013-11-27 | 山东建筑大学 | 一种固体截面温变测定系统 |
CN104331574A (zh) * | 2014-11-19 | 2015-02-04 | 河海大学常州校区 | 基于ansys有限元平台的感应淬火淬硬层深度的预测方法 |
CN105631111A (zh) * | 2015-12-25 | 2016-06-01 | 河海大学常州校区 | 一种台阶轴感应淬火马氏体分布的预测方法 |
CN105653770A (zh) * | 2015-12-25 | 2016-06-08 | 河海大学常州校区 | 一种连续感应淬火淬硬层深度的预测方法 |
CN105930605A (zh) * | 2016-04-29 | 2016-09-07 | 河海大学常州校区 | 感应淬火处理轴零件的淬硬层深度测量方法 |
Non-Patent Citations (1)
Title |
---|
赵正阳等: "感应淬火层深度的预测方法及其精度探讨", 《金属热处理》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113174467A (zh) * | 2021-03-23 | 2021-07-27 | 中冶南方连铸技术工程有限责任公司 | 铸坯淬冷组织的预测方法及铸坯淬冷工艺的制定方法 |
CN113174467B (zh) * | 2021-03-23 | 2023-02-24 | 中冶南方连铸技术工程有限责任公司 | 铸坯淬冷组织的预测方法及铸坯淬冷工艺的制定方法 |
CN116694919A (zh) * | 2023-06-12 | 2023-09-05 | 无锡市海峰海林精密轴承有限公司 | 一种轴承贝氏体热处理工艺的优化方法及系统 |
CN116694919B (zh) * | 2023-06-12 | 2023-11-21 | 无锡市海峰海林精密轴承有限公司 | 一种轴承贝氏体热处理工艺的优化方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Non-isothermal phase-transformation kinetics model for evaluating the austenization of 55CrMo steel based on Johnson–Mehl–Avrami equation | |
Nie et al. | Experimental study and modeling of H13 steel deposition using laser hot-wire additive manufacturing | |
Inoue et al. | Development and implementation of cae system “hearts” for heat treatment simulation based on metallo-thermo-mechanics | |
Huiping et al. | Inverse heat conduction analysis of quenching process using finite-element and optimization method | |
CN105930605A (zh) | 感应淬火处理轴零件的淬硬层深度测量方法 | |
CN104331574A (zh) | 基于ansys有限元平台的感应淬火淬硬层深度的预测方法 | |
Lee et al. | Optimal die profile design for uniform microstructure in hot extruded product | |
Powar et al. | Modeling of microstructure and mechanical properties of heat treated components by using Artificial Neural Network | |
CN103993138B (zh) | 预测高强塑积钢淬火中马氏体组织演变的方法 | |
Rauch et al. | From high accuracy to high efficiency in simulations of processing of Dual-Phase steels | |
Schwenk et al. | Experimental determination of process parameters and material data for numerical modeling of induction hardening | |
Ahlström | Residual stresses generated by repeated local heating events–Modelling of possible mechanisms for crack initiation | |
Barglik et al. | Influence of the magnetic permeability on modeling of induction surface hardening | |
CN107245557A (zh) | 基于ttt曲线预测轴件感应淬火淬硬层组织和硬度方法 | |
Barglik | Mathematical modeling of induction surface hardening | |
CN105653770A (zh) | 一种连续感应淬火淬硬层深度的预测方法 | |
CN105631111A (zh) | 一种台阶轴感应淬火马氏体分布的预测方法 | |
Schlesselmann et al. | Coupled numerical multiphysics simulation methods in induction surface hardening | |
Spezzapria et al. | Numerical simulation of solid–solid phase transformations during induction hardening process | |
Reti et al. | Prediction of as-quenched hardness after rapid austenitization and cooling of surface hardened steels | |
Barglik | Induction hardening of steel elements with complex shapes | |
Wang et al. | Estimation of heat transfer coefficient and phase transformation latent heat by modified pattern search method | |
Saputro et al. | Mobile Induction Heat Treatment of Large‐Sized Spur Gear—The Effect of Scanning Speed and Air Gap on the Uniformity of Hardened Depth and Mechanical Properties | |
Fuhrmann et al. | Numerical simulation of the surface hardening of steel | |
Liang et al. | Numerical and Experimental Study on Hot Forming by Partition Cooling of 38MnB5Nb |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171013 |
|
RJ01 | Rejection of invention patent application after publication |