CN105930605A - 感应淬火处理轴零件的淬硬层深度测量方法 - Google Patents
感应淬火处理轴零件的淬硬层深度测量方法 Download PDFInfo
- Publication number
- CN105930605A CN105930605A CN201610280580.3A CN201610280580A CN105930605A CN 105930605 A CN105930605 A CN 105930605A CN 201610280580 A CN201610280580 A CN 201610280580A CN 105930605 A CN105930605 A CN 105930605A
- Authority
- CN
- China
- Prior art keywords
- workpiece
- depth
- temperature
- heating
- 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
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (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
本发明提供一种感应淬火处理轴零件的淬硬层深度测量方法,该方法基于ANSYS有限元平台,通过模拟计算轴零件感应淬火时工件淬硬层深,可以通过调节感应工艺参数计算不同工件的淬硬层深,因此本发明通过改变感应淬火工艺参数预测淬硬层深度简捷可靠。
Description
技术领域
本发明涉及一种感应淬火处理轴零件的淬硬层深度预测方法。
背景技术
传统材料热处理工艺通过大量试验研究筛选一种较好处理工艺,成本既高又费时。具有高质量、重复性及适应性强等特征的感应淬火热处理是热处理行业中应用最广和发展最快的表面处理工艺之一,通过有限元分析预测工件感应淬火的淬硬层深度,确定其最佳工艺参数并指导生产,避免测定工件的淬硬层深度时损坏试验工件。
发明内容
本发明为了解决现有技术中存在的上述缺陷和不足,提供了一种基于ANSYS有限元平台来模拟感应淬火加热和冷却两个阶段分别预测感应淬火淬硬层深度。
为解决上述技术问题,本发明提供一种感应淬火处理轴零件的淬硬层深度测量方法,利用ANSYS有限元分析平台测量感应淬火工件的淬硬层深度,具体包括以下步骤:
步骤一,建立模型;根据感应淬火实际工作环境构建工件的感应淬火有限元实体模型;将有限元实体模块划分为若干个二维耦合场实体单元,确定二维耦合场实体单元的属性及工件的材料属性,即工艺参数;
步骤二,模拟加热过程;施加热源载荷、热对流约束条件及边界条件至二维耦合场实体单元的各节点上,模拟感应淬火加热过程温度场变化,并采用直接法对工件进行电磁热耦合分析计算,得到工件模拟加热温度随径向分布曲线图;计算工件在失磁后的升温速率,并根据工件连续加热奥氏体化相变温度-升温速率曲线图,计算工件奥氏体化的开始温度Ac1和结束温度Ac3;利用工件奥氏体化的开始温度Ac1和结束温度Ac3在工件模拟加热温度随径向分布曲线图上确定工件奥氏体在径向上的分布范围;
步骤三,模拟冷却过程;模拟加热处理完毕后,加载对流换热系数,模拟感应淬火冷却过程温度场变化;由工件径向上不同深度处的冷却曲线图得到工件在径向不同深度处的700℃冷却速度分布图;采用Maynier组织预测模型与HV硬度计算出工件径向不同深度处的硬度值,得到工件HV硬度径向分布曲线图;根据完全马氏体、半马氏体和无马氏体层的硬度结合工件HV硬度径向分布曲线图确定完全马氏体、半马氏体和无马氏体层的深度;
步骤四,将步骤二和步骤三相互校验。
其中,所述步骤一中,建模过程考虑工件和感应圈之间的间隙和空气场;其中,工件表面划分为密集区,中心部划分为粗糙区,纵向划分精度为0.125mm,感应圈划分为20个单元格;空气场采用自由划分方式。
进一步,工件的材料属性包括感应器及工件相对磁导率、工件的材料密度、工件比热容和工件电阻率。
进一步,边界条件包括环境温度、感应加热电流密度、工件热传导系数及边界换热系数。
进一步,利用工件奥氏体化的开始温度Ac1和结束温度Ac3在工件模拟加热温度随径向分布曲线图上确定工件奥氏体在径向上的分布范围,具体过程为:结束温度Ac3对应径向深度δ100%A为完全奥氏体组织;开始温度Ac1对应径向深度δ0%A为开始奥氏体组织;完全奥氏体组织和开始奥氏体组织的中间温度791℃处径向深度δ50%A为50%奥氏体组织。
进一步,HV硬度计算公式为:
HVM=569.95+21lgV
HVB=-41.8+82.95lgV
HVF+P=175.63+7.6lgV
其中,HVM为马氏体M的维氏硬度;HVB为贝氏体B的维氏硬度;HVF+P为铁素体与珠光体混合F+P的维氏硬度;V为轴径向节点在700℃瞬间冷却速度。
本发明所达到的有益技术效果:本发明提供一种感应淬火处理轴零件的淬硬层深度测量方法,该方法基于ANSYS有限元平台,通过模拟计算轴零件感应淬火时工件淬硬层深,可以通过调节感应工艺参数计算不同工件的淬硬层深,因此本发明通过改变感应淬火工艺参数预测淬硬层深度简捷可靠。
附图说明
图1本发明之实施例工件感应淬火几何模型;
图2本发明之实施例工件有限元实体模型;
图3本发明之实施例工件模拟加热温度随径向分布曲线图;
图4本发明之实施例工件连续加热奥氏体化相变温度-升温速率曲线图;
图5本发明之实施例工件径向-温度分布及奥氏体含量分布图;
图6本发明之实施例工件径向不同深度的冷却曲线图;
图7本发明之实施例工件径向不同深度处的700℃冷却速度分布图;
图8本发明之实施例Maynier组织预测模型;
图9本发明之实施例工件径向硬度分布和实测硬度分布对比曲线图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
本发明提供一种感应淬火处理轴零件的淬硬层深度测量方法,利用ANSYS有限元分析平台测量感应淬火工件的淬硬层深度,具体包括以下步骤:
步骤一,建立模型;根据感应淬火实际工作环境构建工件的感应淬火有限元实体模型;将有限元实体模块划分为若干个二维耦合场实体单元,确定二维耦合场实体单元的属性及工件的材料属性,即工艺参数;
建模过程考虑工件和感应圈之间的间隙和空气场;其中,工件表面划分为密集区,中心部划分为粗糙区,纵向划分精度为0.125mm,感应圈划分为20个单元格;空气场采用自由划分方式。工件的材料属性包括感应器及工件相对磁导率、工件的材料密度、工件比热容和工件电阻率。
步骤二,模拟加热过程;施加热源载荷、热对流约束条件及边界条件至二维耦合场实体单元的各节点上,边界条件包括环境温度、感应加热电流密度、工件热传导系数及边界换热系数;模拟感应淬火加热过程温度场变化,并采用直接法对工件进行电磁热耦合分析计算,得到工件模拟加热温度随径向分布曲线图;计算工件在失磁后的升温速率,并根据工件连续加热奥氏体化相变温度-升温速率曲线图,计算工件奥氏体化的开始温度Ac1和结束温度Ac3;利用工件奥氏体化的开始温度Ac1和结束温度Ac3在工件模拟加热温度随径向分布曲线图上确定工件奥氏体在径向上的分布范围,具体过程为:结束温度Ac3对应径向深度δ100%A为完全奥氏体组织;开始温度Ac1对应径向深度δ0%A为开始奥氏体组织;完全奥氏体组织和开始奥氏体组织的中间温度791℃处径向深度δ50%A为50%奥氏体组织。
步骤三,模拟冷却过程;模拟加热处理完毕后,加载对流换热系数,模拟感应淬火冷却过程温度场变化;由工件径向上不同深度处的冷却曲线图得到工件在径向不同深度处的700℃冷却速度分布图;采用Maynier组织预测模型与HV硬度公式计算出工件径向不同深度处的硬度值,得到工件HV硬度径向分布曲线图;根据完全马氏体、半马氏体和无马氏体层的硬度结合工件HV硬度径向分布曲线图确定完全马氏体、半马氏体和无马氏体层的深度;
HV硬度计算公式为:
HVM=569.95+21lgV
HVB=-41.8+82.95lgV
HVF+P=175.63+7.6lgV
其中,HVM为马氏体M的维氏硬度;HVB为贝氏体B的维氏硬度;HVF+P为铁素体与珠光体混合F+P的维氏硬度;V为轴径向节点在700℃瞬间冷却速度。
步骤四,将步骤二和步骤三相互校验。
实施例
为了更好的说明本发明的具体过程,现利用ANSYS有限元分析平台对45#钢轴零件进行感应淬火进行数值模拟分析,包括以下步骤:
(1)依据45#钢轴零件的技术要求确定感应淬火时实际加热工艺条件:感应加热频率f=195kHz、加热时间t=3.8s、电流密度Js=700×106A/m2;工件的直径Φ16mm,长5mm;感应圈内径为Φ20mm,高为5mm;
(2)对轴对称工件的1/2进行建模:感应器以集肤深度进行建模,同时考虑工件和感应圈之间的间隙及空气场;工件表面划分密集区,中心部划分为粗糙区,纵向划分精度为0.125mm;感应圈划分为20个单元格;空气场采用自由划分的方式,如图1所示,其中A1区代表轴零件,A2区为感应圈,A3为空气场;
(3)确定二维耦合场实体单元PLANE13属性及工件材料属性,材料属性包括材料的密度,在不同温度下的磁导率μ、比热容c、电阻率ρ;
(4)进入电磁热耦合求解过程:确定热源载荷、热对流约束及边界条件,边界条件包括环境温度、感应加热电流密度、工件热传导系数及边界换热系数;施加热源载荷、热对流约束及边界条件至二维耦合场实体单元PLANE13的各节点上,并采用直接法对边界条件进行电磁热耦合分析计算;将二维耦合场实体单元PLANE13转化为温度计算单元PLANE55施加到各节点上,计算得到不同时刻和不同区域的温度场;
(5)由工件模拟加热温度随径向分布曲线图如图3所示,计算工件在失磁后升温速率,根据工件材料连续加热奥氏体化相变温度-升温速率曲线如图4,得到工件的奥氏体化开始温度Ac1和结束温度Ac3,分别为761℃,821℃;
(6)根据温度场计算得到的模拟加热时间t=3.8s时工件径向-温度分布及奥氏体含量分布图如图5所示,奥氏体含量分布具体判断过程如下:结束温度Ac3对应径向深度δ100%A为完全奥氏体组织;开始温度Ac1对应径向深度δ0%A为开始奥氏体组织;完全奥氏体组织和开始奥氏体组织的中间温度791℃处径向深度δ50%A为50%奥氏体组织;
(7)加热结束后对工件表面施加对流换热系数载荷,冷却20s,模拟冷却过程温度场变化;由工件径向上不同深度处冷却曲线,如图6所示,计算工件径向上不同深度处的700℃时冷却速度,得到工件径向不同深度处的700℃冷却速度分布图,如图7所示。
(8)采用Maynier组织预测模型,如图8所示,与HV硬度,即维氏硬度计算公式HVM=569.95+21lgV,HVB=-41.8+82.95lgV,HVF+P=175.63+7.6lgV计算出工件径向上不同深度处的硬度值,得到工件径向硬度分布曲线图,与实际测量硬度曲线相比,采用本方法测得的硬度与实际测量得到的硬度值一致,如图9所示,根据完全马氏体、半马氏体和无马氏体层的硬度结合工件HV硬度径向分布曲线图确定完全马氏体、半马氏体和无马氏体层的深度。
(9)将步骤(6)感应加热过程中预测淬硬层深度和步骤(8)感应淬火冷却过程预测淬硬层深度进行相互校验。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (6)
1.感应淬火处理轴零件的淬硬层深度测量方法,其特征在于:利用ANSYS有限元分析平台测量感应淬火工件的淬硬层深度,具体包括以下步骤:
步骤一,建立模型;根据感应淬火实际工作环境构建工件的感应淬火有限元实体模型;将有限元实体模块划分为若干个二维耦合场实体单元,确定二维耦合场实体单元的属性及工件的材料属性,即工艺参数;
步骤二,模拟加热过程;施加热源载荷、热对流约束条件及边界条件至二维耦合场实体单元的各节点上,模拟感应淬火加热过程温度场变化,并采用直接法对工件进行电磁热耦合分析计算,得到工件模拟加热温度随径向分布曲线图;计算工件在失磁后的升温速率,并根据工件连续加热奥氏体化相变温度-升温速率曲线图,计算工件奥氏体化的开始温度Ac1和结束温度Ac3;利用工件奥氏体化的开始温度Ac1和结束温度Ac3在工件模拟加热温度随径向分布曲线图上确定工件奥氏体在径向上的分布范围;
步骤三,模拟冷却过程;模拟加热处理完毕后,加载对流换热系数,模拟感应淬火冷却过程温度场变化;由工件径向上不同深度处的冷却曲线图得到工件在径向不同深度处的700℃冷却速度分布图;采用Maynier组织预测模型与HV硬度计算出工件径向不同深度处的硬度值,得到工件HV硬度径向分布曲线图;根据完全马氏体、半马氏体和无马氏体层的硬度结合工件HV硬度径向分布曲线图确定完全马氏体、半马氏体和无马氏体层的深度;
步骤四,将步骤二和步骤三相互校验。
2.根据权利要求1所述的感应淬火处理轴零件的淬硬层深度测量方法,其特征在于:所述步骤一中,建模过程考虑工件和感应圈之间的间隙和空气场;其中,工件表面划分为密集区,中心部划分为粗糙区,纵向划分精度为0.125mm,感应圈划分为20个单元格;空气场采用自由划分方式。
3.根据权利要求1所述的感应淬火处理轴零件的淬硬层深度测量方法,其特征在于:工件的材料属性包括感应器及工件相对磁导率、工件的材料密度、工件比热容和工件电阻率。
4.根据权利要求1所述的感应淬火处理轴零件的淬硬层深度测量方法,其特征在于:边界条件包括环境温度、感应加热电流密度、工件热传导系数及边界换热系数。
5.根据权利要求1所述的感应淬火处理轴零件的淬硬层深度测量方法,其特征在于:利用工件奥氏体化的开始温度Ac1和结束温度Ac3在工件模拟加热温度随径向分布曲线图上确定工件奥氏体在径向上的分布范围,具体过程为:结束温度Ac3对应径向深度δ100%A为完全奥氏体组织;开始温度Ac1对应径向深度δ0%A为开始奥氏体组织;完全奥氏体组织和开始奥氏体组织的中间温度791℃处径向深度δ50%A为50%奥氏体组织。
6.根据权利要求1所述的感应淬火处理轴零件的淬硬层深度测量方法,其特征在于:HV硬度计算公式为:
HVM=569.95+21lgV
HVB=-41.8+82.95lgV
HVF+P=175.63+7.6lgV
其中,HVM为马氏体M的维氏硬度;HVB为贝氏体B的维氏硬度;HVF+P为铁素体与珠光体混合F+P的维氏硬度;V为轴径向节点在700℃瞬间冷却速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610280580.3A CN105930605A (zh) | 2016-04-29 | 2016-04-29 | 感应淬火处理轴零件的淬硬层深度测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610280580.3A CN105930605A (zh) | 2016-04-29 | 2016-04-29 | 感应淬火处理轴零件的淬硬层深度测量方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105930605A true CN105930605A (zh) | 2016-09-07 |
Family
ID=56837740
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610280580.3A Pending CN105930605A (zh) | 2016-04-29 | 2016-04-29 | 感应淬火处理轴零件的淬硬层深度测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105930605A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107245552A (zh) * | 2017-06-21 | 2017-10-13 | 河海大学常州校区 | 基于cct曲线预测轴件感应淬火淬硬层组织和硬度方法 |
CN107245557A (zh) * | 2017-06-21 | 2017-10-13 | 河海大学常州校区 | 基于ttt曲线预测轴件感应淬火淬硬层组织和硬度方法 |
CN107904393A (zh) * | 2017-12-08 | 2018-04-13 | 徐工集团工程机械有限公司 | 机械零件热处理强化工艺要求的确定方法 |
CN107907566A (zh) * | 2017-12-11 | 2018-04-13 | 中钢集团邢台机械轧辊有限公司 | 一种预测金属材料激光淬硬深度的试验方法 |
CN112084603A (zh) * | 2020-09-17 | 2020-12-15 | 泰尔重工股份有限公司 | 重载万向轴叉头淬火加热工艺参数的获取方法 |
CN114250343A (zh) * | 2021-12-10 | 2022-03-29 | 浙江欧迪恩传动科技股份有限公司 | 一种降低芯轴强度差异的生产、计算、加热及校验方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101379203A (zh) * | 2006-02-01 | 2009-03-04 | 本田技研工业株式会社 | 感应淬火方法 |
EP2586548A1 (en) * | 2011-10-31 | 2013-05-01 | Alstom Technology Ltd | Component or coupon for being used under high thermal and stress load and method for manufacturing such component or coupon |
CN104212969A (zh) * | 2014-09-18 | 2014-12-17 | 上海交通大学 | 一种基于数值模拟的钢管连续淬火过程控制方法 |
CN104331574A (zh) * | 2014-11-19 | 2015-02-04 | 河海大学常州校区 | 基于ansys有限元平台的感应淬火淬硬层深度的预测方法 |
-
2016
- 2016-04-29 CN CN201610280580.3A patent/CN105930605A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101379203A (zh) * | 2006-02-01 | 2009-03-04 | 本田技研工业株式会社 | 感应淬火方法 |
EP2586548A1 (en) * | 2011-10-31 | 2013-05-01 | Alstom Technology Ltd | Component or coupon for being used under high thermal and stress load and method for manufacturing such component or coupon |
CN104212969A (zh) * | 2014-09-18 | 2014-12-17 | 上海交通大学 | 一种基于数值模拟的钢管连续淬火过程控制方法 |
CN104331574A (zh) * | 2014-11-19 | 2015-02-04 | 河海大学常州校区 | 基于ansys有限元平台的感应淬火淬硬层深度的预测方法 |
Non-Patent Citations (2)
Title |
---|
张根元 等: "感应淬火工艺参数优化和组织硬度分布预测", 《材料热处理学报》 * |
陈珺 等: "45钢光轴连续感应淬火过程的数值模拟", 《金属热处理》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107245552A (zh) * | 2017-06-21 | 2017-10-13 | 河海大学常州校区 | 基于cct曲线预测轴件感应淬火淬硬层组织和硬度方法 |
CN107245557A (zh) * | 2017-06-21 | 2017-10-13 | 河海大学常州校区 | 基于ttt曲线预测轴件感应淬火淬硬层组织和硬度方法 |
CN107904393A (zh) * | 2017-12-08 | 2018-04-13 | 徐工集团工程机械有限公司 | 机械零件热处理强化工艺要求的确定方法 |
CN107907566A (zh) * | 2017-12-11 | 2018-04-13 | 中钢集团邢台机械轧辊有限公司 | 一种预测金属材料激光淬硬深度的试验方法 |
CN112084603A (zh) * | 2020-09-17 | 2020-12-15 | 泰尔重工股份有限公司 | 重载万向轴叉头淬火加热工艺参数的获取方法 |
CN114250343A (zh) * | 2021-12-10 | 2022-03-29 | 浙江欧迪恩传动科技股份有限公司 | 一种降低芯轴强度差异的生产、计算、加热及校验方法 |
CN114250343B (zh) * | 2021-12-10 | 2023-10-31 | 浙江欧迪恩传动科技股份有限公司 | 一种降低芯轴强度差异的生产、计算、加热及校验方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105930605A (zh) | 感应淬火处理轴零件的淬硬层深度测量方法 | |
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有限元平台的感应淬火淬硬层深度的预测方法 | |
Choi et al. | Prediction of high-frequency induction hardening depth of an AISI 1045 specimen by finite element analysis and experiments | |
Barglik | Mathematical modeling of induction surface hardening | |
Eastwood et al. | An induction hardening process model to assist sustainability assessment of a steel bevel gear | |
CN105631111A (zh) | 一种台阶轴感应淬火马氏体分布的预测方法 | |
CN105653770A (zh) | 一种连续感应淬火淬硬层深度的预测方法 | |
Barglik | Induction hardening of steel elements with complex shapes | |
Spezzapria et al. | Numerical simulation of solid–solid phase transformations during induction hardening process | |
Schlesselmann et al. | Coupled numerical multiphysics simulation methods in induction surface hardening | |
Wang et al. | Estimation of heat transfer coefficient and phase transformation latent heat by modified pattern search method | |
CN107245557A (zh) | 基于ttt曲线预测轴件感应淬火淬硬层组织和硬度方法 | |
Wang et al. | Finite-element simulation of moving induction heat treatment | |
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 | |
Barglik | Identification of temperature and hardness distribution during dual frequency induction hardening of gear wheels | |
CN105838869A (zh) | 一种钢板淬火炉加热工艺在线调整方法 | |
Barglik | Induction contour hardening of gear wheels made of steel 300M | |
Bukanin et al. | Simulation of induction heat treatment of steel articles with the help of ELTA 6.0 and 2DELTA software | |
Li et al. | Controlling gear distortion and residual stresses during induction hardening | |
Zhu et al. | Theoretical and experimental analysis of two-pass spot continual induction hardening of AISI 1045 steel | |
Barglik et al. | Hardness and microstructure distributions in gear wheels made of steel AISI 4340 after consecutive dual frequency induction hardening | |
Rudnev | A common misassumption in induction hardening | |
Barglik et al. | Numerical modeling of induction hardening of gear wheels made of steel AMS 6419 | |
Suwanpinij et al. | Phase transformation modelling and parameter identification from dilatometric investigations |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160907 |