CN107245552A - 基于cct曲线预测轴件感应淬火淬硬层组织和硬度方法 - Google Patents
基于cct曲线预测轴件感应淬火淬硬层组织和硬度方法 Download PDFInfo
- Publication number
- CN107245552A CN107245552A CN201710473892.0A CN201710473892A CN107245552A CN 107245552 A CN107245552 A CN 107245552A CN 201710473892 A CN201710473892 A CN 201710473892A CN 107245552 A CN107245552 A CN 107245552A
- Authority
- CN
- China
- Prior art keywords
- shaft member
- impewdance matching
- hardness
- cct
- tissue
- 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/06—Surface hardening
- C21D1/09—Surface hardening by direct application of electrical or wave energy; by particle radiation
- C21D1/10—Surface hardening by direct application of electrical or wave energy; by particle radiation by electric induction
-
- 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
- C21D11/00—Process control or regulation for heat treatments
- C21D11/005—Process control or regulation for heat treatments for cooling
-
- 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
- C21D9/00—Heat treatment, e.g. annealing, hardening, quenching or tempering, adapted for particular articles; Furnaces therefor
- C21D9/28—Heat treatment, e.g. annealing, hardening, quenching or tempering, adapted for particular articles; Furnaces therefor for plain shafts
-
- 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]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P10/00—Technologies related to metal processing
- Y02P10/25—Process efficiency
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Organic Chemistry (AREA)
- Crystallography & Structural Chemistry (AREA)
- Mechanical Engineering (AREA)
- Materials Engineering (AREA)
- Metallurgy (AREA)
- Thermal Sciences (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Heat Treatment Of Articles (AREA)
Abstract
本发明公开了一种基于CCT曲线预测感应淬火淬硬层组织和硬度的方法,基于ANSYS有限元平台来模拟感应淬火过程,结合材料的CCT曲线预测感应淬火淬硬层组织和硬度,减小感应淬火试验工作量。本发明所达到的有益效果:本方法基于ANSYS有限元平台,通过模拟计算轴件感应淬火时工件淬硬层组织和硬度,调节感应工艺参数来满足工件的技术指标,因此本发明通过改变感应淬火工艺参数预测淬硬层简捷可靠。
Description
技术领域
本发明涉及一种基于CCT曲线预测感应淬火淬硬层组织和硬度的方法。
背景技术
传统材料热处理工艺通过大量试验研究筛选一种较好处理工艺,成本既高又费时。具有高质量、重复性及适应性强等特征的感应淬火热处理是热处理行业中应用最广和发展最快的表面处理工艺之一,通过有限元分析预测工件感应淬火淬硬层的组织和硬度,确定其最佳工艺参数并指导生产,避免测定工件的淬硬层深度时损坏试验工件。
发明内容
为解决现有技术的不足,本发明的目的在于提供一种基于CCT曲线预测轴件感应淬火淬硬层组织和硬度方法,基于ANSYS有限元平台来模拟感应淬火过程,结合材料CCT曲线预测轴件感应淬火淬硬层的组织和硬度,减小感应淬火试验工作量。
为了实现上述目标,本发明采用如下的技术方案:
一种基于CCT曲线和ANSYS有限元分析软件平台预测轴件感应淬火淬硬层组织和硬度的方法,其特征是,包括如下步骤:
步骤1)依据感应淬火实际工作环境构建感应淬火有限元实体模型;
步骤2)将模型划分网格单元,定义二维耦合场实体单元PLANE13属性及材料属性;施加热源载荷、热对流约束及边界条件至二维耦合场实体单元PLANE13各节点上,采用直接法对轴件进行电磁热耦合分析,得到电磁热耦合计算结果;
步骤3)计算轴件在材料失磁后的升温速率,由轴件材料连续加热升温速率VAc1/Ac3与Ac1和Ac3变化曲线计算得到轴件奥氏体化温度Ac1和Ac3;
步骤4)通过感应淬火的工艺参数计算得到轴件模拟加热温度随轴径向分布曲线,由步骤3)中计算所得Ac1和Ac3值确定轴径向分布上奥氏体含量范围;
步骤5)感应加热结束后对轴件表面施加对流换热系数载荷,模拟感应淬火冷却过程温度场变化,得到轴径向淬硬层深度范围内任一位置的冷却曲线;
步骤6)查找《热处理技术数据手册》得到上述奥氏体化温度下对应工件的CCT图,采用数值处理方法得到数值化后的CCT图,并根据图中的数据信息拟合处理得到材料奥氏体化后的冷却速度与其硬度的关系;
步骤7)将步骤5)中感应淬火轴径向各处的冷却曲线与数值处理的CCT图结合,根据感应淬火轴径向各处的冷却曲线得到CCT图上相对应的冷却速度及步骤6)中材料奥氏体化后的冷却速度与其硬度的关系,计算得到感应淬火淬硬层的轴径向相应各处的奥氏体淬火后的组织及硬度;
步骤8)根据步骤4)计算得到感应淬火的淬硬层中轴径向奥氏体含量分布,结合步骤7)中的计算结果得到轴件感应淬火后淬硬层的组织及硬度;
步骤9)将步骤8)感应淬火淬硬层的组织及硬度计算结果与最终实验结果对比分析误差,给出模拟预测的可行性。
进一步地,所述步骤2)中采用直接法对工件进行电磁热耦合分析,得到电磁热耦合计算结果。
进一步地,所述步骤2)中属性包括材料物性参数、模型网格大小和迭代时间步长;所述迭代时间不长依据经验值选取。
进一步地,所示的步骤4)中奥氏体在轴径向分布上含量范围的确定采用线性插值方法。
进一步地,所示的步骤6)中采用数值处理方法对CCT图进行处理时采用软件GetDate进行信息采集和Origin软件对轴件径向各处的冷却速度所获取的组织和其硬度的数学关系拟合出相应公式。
进一步地,所示的步骤1)中感应淬火实际工作环境包括感应加热频率、加热时间、电流密度、感应线圈内径及高度。
本发明所达到的有益效果:本方法基于ANSYS有限元平台,通过模拟计算轴件感应淬火淬硬层的组织和硬度,调节感应工艺参数来满足工件的技术指标,因此本发明通过改变感应淬火工艺参数预测淬硬层简捷可靠。
附图说明
图1是轴件感应淬火几何模型示意图;
图2是轴件有限元分析网格模型的示意图;
图3是轴件表面加热温度随时间分布曲线示意图;
图4是工件材料连续加热升温速率与Ac1和Ac3变化曲线示意图;
图5是轴件径向温度分布及奥氏体含量分布示意图;
图6是轴件径向不同深度的冷却曲线示意图;
图7是880℃温度奥氏体化的45钢CCT图;
图8是880℃温度奥氏体化的45钢CCT图数值化处理结果的示意图;
图9是45钢奥氏体化后的各冷却速度下获取组织与其硬度的关系示意图;
图10是轴件距表面2.21mm处的淬火冷却曲线与奥氏体化温度为880℃的45钢CCT图结合情况示意图;
图11是轴件距表面2.21mm处感应淬火后的金相组织示意图;
图12是轴件径向HV硬度预测分布和实测硬度分布曲线示意图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
如图1~4所示,利本实施例中用ANSYS有限元分析软件对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钢的CCT图如图7,采用数值处理方法得到数值化后的CCT图如图8,并根据图中的数据信息拟合处理得到材料奥氏体化后的各冷却速度下获取组织与其硬度的关系如图9。
(7)将步骤(5)中的感应淬火的冷却过程与数值处理的奥氏体化温度880℃的45钢CCT图结合,根据冷却曲线得到的各冷却速度及步骤(8)中的各冷却速度下的组织与其硬度的关系,计算得到轴径向奥氏体淬火后的组织及硬度。
(8)根据步骤(4)计算得到的轴径向上奥氏体含量分布,结合步骤(7)中的计算结果按百分比计算得到轴件感应淬火后淬硬层的组织及硬度,如轴件距表面2.21mm处的组织含量计算结果为49.6%马氏体(M)组织,2.0%贝氏体(B)组织和48.4%铁素体及珠光体(F+P)组织,实测此处的金相组织图片如图11。
(9)将步骤(8)的感应淬火淬硬层的硬度计算结果与最终试验结果对比分析误差如图12,给出模拟预测的可行性。
图8中所标注的数字,图中上部数字为冷却速度,单位℃/s;下部数字为组织维式硬度(HV),字母等采用的是国际通用的定义。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (6)
1.一种基于CCT曲线预测轴件感应淬火淬硬层组织和硬度方法,其特征是,包括如下步骤:
步骤1)依据感应淬火实际工作环境构建感应淬火有限元实体模型;
步骤2)将模型划分网格单元,定义二维耦合场实体单元PLANE13属性及材料属性;施加热源载荷、热对流约束及边界条件至二维耦合场实体单元PLANE13各节点上,采用直接法对轴件进行电磁热耦合分析,得到电磁热耦合计算结果;
步骤3)计算轴件在材料失磁后的升温速率,由轴件材料连续加热升温速率VAc1/Ac3与Ac1和Ac3变化曲线计算得到轴件奥氏体化温度Ac1和Ac3;
步骤4)通过感应淬火的工艺参数计算得到轴件模拟加热温度随轴径向分布曲线,由步骤3)中计算所得Ac1和Ac3值确定轴径向分布上奥氏体含量范围;
步骤5)感应加热结束后对轴件表面施加对流换热系数载荷,模拟感应淬火冷却过程温度场变化,得到轴径向淬硬层深度范围内任一位置的冷却曲线;
步骤6)查找《热处理技术数据手册》得到上述奥氏体化温度下对应工件的CCT图,采用数值处理方法得到数值化后的CCT图,并根据图中的数据信息拟合处理得到材料奥氏体化后的冷却速度与其硬度的关系;
步骤7)将步骤5)中的感应淬火轴径向各处的冷却曲线与数值处理的CCT图结合,根据感应淬火轴径向各处的冷却曲线得到CCT图上相对应的冷却速度及步骤6)中的材料奥氏体化后的冷却速度与其硬度的关系,计算得到感应淬火淬硬层的轴径向相应各处的奥氏体淬火后的组织及硬度;
步骤8)根据步骤4)计算得到感应淬火的淬硬层中轴径向奥氏体含量分布,结合步骤7)中的计算结果得到轴件感应淬火后淬硬层的组织及硬度;
步骤9)将步骤8)的感应淬火淬硬层的组织及硬度计算结果与最终实验结果对比分析误差,给出模拟预测的可行性。
2.根据权利要求1所述的一种基于CCT曲线预测轴件感应淬火淬硬层组织和硬度方法,其特征是,所述步骤2)中采用直接法对工件进行电磁热耦合分析,得到电磁热耦合计算结果。
3.根据权利要求1所述的一种基于CCT曲线预测轴件感应淬火淬硬层组织和硬度方法,其特征是,所述步骤2)中属性包括材料物性参数、模型网格大小和迭代时间步长;所述迭代时间不长依据经验值选取。
4.根据权利要求1所述的一种基于CCT曲线预测轴件感应淬火淬硬层组织和硬度方法,其特征是,所示的步骤4)中奥氏体在轴径向分布上含量范围的确定采用线性插值方法。
5.根据权利要求1所述的一种基于CCT曲线预测轴件感应淬火淬硬层组织和硬度方法,其特征是,所示的步骤6)中采用数值处理方法对CCT图进行处理时采用软件GetDate进行信息采集和Origin软件对轴件径向各处的冷却速度所获取的组织和其硬度的数学关系拟合出相应公式。
6.根据权利要求1所述的一种基于CCT曲线预测轴件感应淬火淬硬层组织和硬度方法,其特征是,所示的步骤1)中感应淬火实际工作环境包括感应加热频率、加热时间、电流密度、感应线圈内径及高度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710473892.0A CN107245552A (zh) | 2017-06-21 | 2017-06-21 | 基于cct曲线预测轴件感应淬火淬硬层组织和硬度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710473892.0A CN107245552A (zh) | 2017-06-21 | 2017-06-21 | 基于cct曲线预测轴件感应淬火淬硬层组织和硬度方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107245552A true CN107245552A (zh) | 2017-10-13 |
Family
ID=60018554
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710473892.0A Pending CN107245552A (zh) | 2017-06-21 | 2017-06-21 | 基于cct曲线预测轴件感应淬火淬硬层组织和硬度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107245552A (zh) |
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 CN201710473892.0A patent/CN107245552A/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 |
---|
赵正阳等: "感应淬火层深度的预测方法及其精度探讨", 《金属热处理》 * |
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 | |
CN104331574A (zh) | 基于ansys有限元平台的感应淬火淬硬层深度的预测方法 | |
Cajner et al. | Computer simulation of induction hardening | |
CN105930605A (zh) | 感应淬火处理轴零件的淬硬层深度测量方法 | |
Lee et al. | Optimal die profile design for uniform microstructure in hot extruded product | |
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 | |
CN105653770B (zh) | 一种连续感应淬火淬硬层深度的预测方法 | |
CN105631111A (zh) | 一种台阶轴感应淬火马氏体分布的预测方法 | |
Spezzapria et al. | Numerical simulation of solid–solid phase transformations during induction hardening process | |
CN107245557A (zh) | 基于ttt曲线预测轴件感应淬火淬硬层组织和硬度方法 | |
Schlesselmann et al. | Coupled numerical multiphysics simulation methods in induction surface hardening | |
Karban et al. | Continual induction hardening of steel bodies | |
Reti et al. | Prediction of as-quenched hardness after rapid austenitization and cooling of surface hardened steels | |
Wang et al. | Estimation of heat transfer coefficient and phase transformation latent heat by modified pattern search method | |
CN107245552A (zh) | 基于cct曲线预测轴件感应淬火淬硬层组织和硬度方法 | |
Baldan et al. | A parallel multi-fidelity optimization approach in induction hardening | |
Kulawik et al. | Modelling of heat treatment of steel elements with the movement of coolant | |
Milenin et al. | Model of curvature of crankshaft blank during heat treatment, accounting for phase transformations | |
Yermekova et al. | Numerical simulation and automatic optimization of the disk induction heating system | |
Elmaryami et al. | Effect of radius on temperature history of transient industrial quenched chromium Steel-8650H by developing 1-D MM | |
JP4631247B2 (ja) | 鋼材の熱処理方法及びそのプログラム | |
Spezzapria et al. | Multiphysical-multiscale FEM simulation of contour induction hardening on aeronautical gears | |
Barglik | Identification of temperature and hardness distribution during dual frequency induction hardening of gear wheels | |
Kurek et al. | Modeling of induction hardening |
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 |