CN104331574A - 基于ansys有限元平台的感应淬火淬硬层深度的预测方法 - Google Patents
基于ansys有限元平台的感应淬火淬硬层深度的预测方法 Download PDFInfo
- Publication number
- CN104331574A CN104331574A CN201410665164.6A CN201410665164A CN104331574A CN 104331574 A CN104331574 A CN 104331574A CN 201410665164 A CN201410665164 A CN 201410665164A CN 104331574 A CN104331574 A CN 104331574A
- Authority
- CN
- China
- Prior art keywords
- workpiece
- finite element
- temperature
- induction
- ansys finite
- 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
Landscapes
- Heat Treatment Of Articles (AREA)
Abstract
本发明公开了一种基于ANSYS有限元平台的感应淬火淬硬层深度的预测方法,该方法的步骤如下:(a)根据工件的技术要求,确定感应淬火时的实际工作环境;(b)根据实际工作环境在ANSYS有限元平台上构建实际感应淬火环境的有限元实体模型;(c)首先定义二维耦合场实体单元PLANE13属性,再定义材料属性,然后划分适当的有限元网格;(d)施加热源载荷和热对流约束以及边界条件至二维耦合场实体单元PLANE13的各节点上,并采用直接法对工件进行电磁热耦合分析,得到电磁热耦合计算结果;(e)对工件温度场进行计算分析。本发明可根据实际工件尺寸和设备参数改变来修改模型、载荷及约束条件,极大地减少了工件的感应淬火试验工作量,提高了分析效率和技术参数测试精度,降低了试验成本。
Description
技术领域
本发明涉及一种感应淬火淬硬层深度预测方法,尤其涉及一种基于ANSYS有限元平台的感应淬火淬硬层深度预测方法。
背景技术
目前,感应淬火热处理具有高质量,再现性及适应性强等特征,是热处理行业中应用最广和发展最快的表面处理工艺之一,对其进行有限元分析可以预测感应淬火工件淬硬层深度,确定最佳工艺参数并指导生产。
为了得到满足要求的表面硬度及淬硬层深度、尺寸、组织等要求,需要对感应淬火过程中的加热、冷却及回火过程的各参数进行严格控制。
传统的材料热处理工艺通过大量试验研究,从中筛选一种比较好的处理工艺,但费时和高成本。
发明内容
本发明所要解决的技术问题是克服现有技术的缺陷,提供一种基于ANSYS有限元平台的感应淬火淬硬层深度预测方法,它可根据实际工件尺寸和设备参数改变来修改模型、载荷及约束条件,极大地减少了工件的感应淬火试验工作量,提高了分析效率和技术参数测试精度,降低了试验成本。
为了解决上述技术问题,本发明的技术方案是:一种基于ANSYS有限元平台的感应淬火淬硬层深度的预测方法,该方法的步骤如下:
(a)根据工件的技术要求,确定感应淬火时的实际工作环境;
(b)根据实际工作环境在ANSYS有限元平台上构建实际感应淬火环境的有限元实体模型;
(c)首先定义二维耦合场实体单元PLANE13属性,再定义材料属性,然后划分适当的有限元网格;
(d)施加热源载荷和热对流约束以及边界条件至二维耦合场实体单元PLANE13的各节点上,并采用直接法对工件进行电磁热耦合分析,得到电磁热耦合计算结果;
(e)定义二维4节点热实体单元PLANE55属性,将经过(d)得到的电磁热耦合计算结果施加到二维4节点热实体单元PLANE55的各节点上,对工件进行温度场分析,得到温度场求解结果;
(f)根据得到的温度场求解结果导出工件的温度分布曲线,根据该温度分布曲线来判断淬硬层深度情况。
进一步,所述的步骤(f)中,根据工件材料选择适当的完全奥氏体化温度和开始奥氏体化温度,结合感应加热升温快的特点在原温度上加上一定的偏移值,再通过温度分布曲线来判断淬硬层深度情况。
进一步,所述的步骤(b)中的实际工作环境包括工件材质、工件和感应器的尺寸、工件的散热条件、工件与感应器之间的距离、空气场。
进一步,所述的步骤(c)中的材料属性包括感应器及工件的相对磁导率、工件的材料密度、工件的比热容和工件的电阻率。
进一步,所述的步骤(d)中的热源载荷和热对流约束以及边界条件包括环境温度、感应加热电流密度、工件的热传导系数及边界换热系数。
采用了上述技术方案后,由于本发明基于ANSYS有限元平台的感应淬火淬硬层深度预测方法,通过模拟计算与对比不同感应加热工艺参数,可确定感应加热工件温度场分布情况及根据工艺要求确定最佳工艺参数,从而减少试验工作量,提高感应淬火工件质量,因此本发明的基于ANSYS有限元平台的感应淬火淬硬层深度预测方法不仅使用成本较低,同时分析效率及精确度较高。
附图说明
图1是建立的轴类零件感应加热几何模型;
图2是工件的有限元分析网格模型;
图3是时间t=4.2s的温度场云图;
图4是感应加热过程中各关键点随时间变化曲线;
图5是时间t=4.2s时工件中部温度的径向温度分布曲线。
具体实施方式
为了使本发明的内容更容易被清楚地理解,下面根据具体实施例并结合附图,对本发明作进一步详细的说明。
如图1~4所示,利用ANSYS有限元分析软件对S45C钢光轴进行感应淬火中的感应加热进行温度场分析,包括以下步骤:
(a)根据工件的技术要求,确定感应淬火时的实际工作环境:其中,需要确定感应加热条件:感应加热频率f=186kHz、加热时间t=4.2s、电流密度Js=700×106A/m2;工件是直径Φ16mm,长度为30mm的光轴;感应线圈内径为Φ20mm,高度为8mm;
(b)对轴对称工件的1/2进行一维建模,感应器以集肤深度进行建模,同时考虑工件和感应圈之间的间隙,以及周围的空气场。工件表面划分较密,心部较粗,纵向划分精度为0.25mm;感应圈划分为20×32个单元格;空气场采用自由划分的方式。模型如图1所示,其中A1区代表工件,A2区为感应圈,A3为空气场;
(c)定义二维耦合场实体单元PLANE13属性及材料属性,材料属性包括材料的密度,在不同温度下的磁导率μ、比热容c、电阻率ρ;
(d)进入电磁热耦合求解过程:定义热源载荷和热对流约束以及边界条件,其包括环境温度、感应加热电流密度、工件的热传导系数及边界换热系数,施加热源载荷和热对流约束以及边界条件至二维耦合场实体单元PLANE13的各节点上,并采用直接法对工件进行电磁热耦合分析,得到电磁热耦合计算结果;
(e)再将耦合计算单元PLANE13转化为温度计算单元PLANE55,读取电磁热耦合计算结果作为载荷施加到各节点上,进行温度场的求解;温度场为不同时刻和不同区域温度,时间t=4.2s的温度场云图如图3所示,
(f)读取温度场求解结果,导出沿工件径向温度分布曲线,根据该温度分布曲线判断淬硬层深度情况,其中,感应加热过程中各关键点随时间变化曲线如图4所示,时间t=4.2s时工件中部温度的径向温度分布曲线如图5所示,具体判断过程如下:材料S45C钢的值分别为778℃、721℃,由于感应加热时间短,需要比普通淬火过热30~50℃以满足完全、开始奥氏体化,所以S45C钢感应淬火取完全奥氏体化温度820℃,超出此温度淬火冷却为100%马氏体组织(M),对应完全淬硬层;开始奥氏体化温度760℃,低于此温度全部为材料处理前的组织,对应淬硬层;完全奥氏体化和开始奥氏体化中间温度为790℃,此处淬火产生50%M组织,对应有效淬硬层,因此可以通过温度分布曲线来判断淬硬层深度,其中完全淬硬层深度为1.62mm,有效淬硬层深度为2.02mm,淬硬层深度为2.40mm。
以上所述的具体实施例,对本发明解决的技术问题、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种基于ANSYS有限元平台的感应淬火淬硬层深度的预测方法,其特征在于该方法的步骤如下:
(a)根据工件的技术要求,确定感应淬火时的实际工作环境;
(b)根据实际工作环境在ANSYS有限元平台上构建实际感应淬火环境的有限元实体模型;
(c)首先定义二维耦合场实体单元PLANE13属性,再定义材料属性,然后划分适当的有限元网格;
(d)施加热源载荷和热对流约束以及边界条件至二维耦合场实体单元PLANE13的各节点上,并采用直接法对工件进行电磁热耦合分析,得到电磁热耦合计算结果;
(e)定义二维4节点热实体单元PLANE55属性,将经过(d)得到的电磁热耦合计算结果施加到二维4节点热实体单元PLANE55的各节点上,对工件进行温度场分析,得到温度场求解结果;
(f)根据得到的温度场求解结果导出工件的温度分布曲线,根据该温度分布曲线来判断淬硬层深度情况。
2.根据权利要求1所述的基于ANSYS有限元平台的感应淬火淬硬层深度的预测方法,其特征在于:所述的步骤(f)中,根据工件材料选择适当的完全奥氏体化温度和开始奥氏体化温度,结合感应加热升温快的特点在原温度上加上一定的偏移值,再通过温度分布曲线来判断淬硬层深度情况。
3.根据权利要求1所述的基于ANSYS有限元平台的感应淬火淬硬层深度的预测方法,其特征在于:所述的步骤(b)中的实际工作环境包括工件材质、工件和感应器的尺寸、工件的散热条件、工件与感应器之间的距离、空气场。
4.根据权利要求1所述的基于ANSYS有限元平台的感应淬火淬硬层深度的预测方法,其特征在于:所述的步骤(c)中的材料属性包括感应器及工件的相对磁导率、工件的材料密度、工件的比热容和工件的电阻率。
5.根据权利要求1所述的基于ANSYS有限元平台的感应淬火淬硬层深度的预测方法,其特征在于:所述的步骤(d)中的热源载荷和热对流约束以及边界条件包括环境温度、感应加热电流密度、工件的热传导系数及边界换热系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410665164.6A CN104331574A (zh) | 2014-11-19 | 2014-11-19 | 基于ansys有限元平台的感应淬火淬硬层深度的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410665164.6A CN104331574A (zh) | 2014-11-19 | 2014-11-19 | 基于ansys有限元平台的感应淬火淬硬层深度的预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104331574A true CN104331574A (zh) | 2015-02-04 |
Family
ID=52406299
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410665164.6A Pending CN104331574A (zh) | 2014-11-19 | 2014-11-19 | 基于ansys有限元平台的感应淬火淬硬层深度的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104331574A (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105586476A (zh) * | 2016-01-19 | 2016-05-18 | 上海交通大学 | 基于数值模拟技术的高强钢淬火工艺优化方法 |
CN105631111A (zh) * | 2015-12-25 | 2016-06-01 | 河海大学常州校区 | 一种台阶轴感应淬火马氏体分布的预测方法 |
CN105653770A (zh) * | 2015-12-25 | 2016-06-08 | 河海大学常州校区 | 一种连续感应淬火淬硬层深度的预测方法 |
CN105930605A (zh) * | 2016-04-29 | 2016-09-07 | 河海大学常州校区 | 感应淬火处理轴零件的淬硬层深度测量方法 |
CN107169166A (zh) * | 2017-04-17 | 2017-09-15 | 燕山大学 | 一种基于ansys仿真平台的大齿圈移动式感应加热的计算方法 |
CN107245552A (zh) * | 2017-06-21 | 2017-10-13 | 河海大学常州校区 | 基于cct曲线预测轴件感应淬火淬硬层组织和硬度方法 |
CN107245557A (zh) * | 2017-06-21 | 2017-10-13 | 河海大学常州校区 | 基于ttt曲线预测轴件感应淬火淬硬层组织和硬度方法 |
CN107907566A (zh) * | 2017-12-11 | 2018-04-13 | 中钢集团邢台机械轧辊有限公司 | 一种预测金属材料激光淬硬深度的试验方法 |
CN112084603A (zh) * | 2020-09-17 | 2020-12-15 | 泰尔重工股份有限公司 | 重载万向轴叉头淬火加热工艺参数的获取方法 |
CN113361159A (zh) * | 2021-05-31 | 2021-09-07 | 西安建筑科技大学 | 一种喷嘴射流冲击淬火的移动平板温度场仿真方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103238064A (zh) * | 2010-10-26 | 2013-08-07 | 高周波热錬株式会社 | 淬火深度测定方法以及淬火深度测定装置 |
-
2014
- 2014-11-19 CN CN201410665164.6A patent/CN104331574A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103238064A (zh) * | 2010-10-26 | 2013-08-07 | 高周波热錬株式会社 | 淬火深度测定方法以及淬火深度测定装置 |
Non-Patent Citations (4)
Title |
---|
吴金富: ""基于ANSYS的感应加热数值模拟分析"", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
周跃庆: ""感应淬火电磁热耦合场的有限元分析"", 《金属热处理》 * |
张根元: ""基于ANSYS 模拟预测S45C钢轴感应加热工艺参数"", 《金属热处理》 * |
许雪峰: ""45钢工件感应淬火温度场有限元模拟分析"", 《第三届中国热处理活动周暨第六次全国热处理生产技术改造会议论文专辑》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105631111A (zh) * | 2015-12-25 | 2016-06-01 | 河海大学常州校区 | 一种台阶轴感应淬火马氏体分布的预测方法 |
CN105653770A (zh) * | 2015-12-25 | 2016-06-08 | 河海大学常州校区 | 一种连续感应淬火淬硬层深度的预测方法 |
CN105653770B (zh) * | 2015-12-25 | 2018-11-06 | 河海大学常州校区 | 一种连续感应淬火淬硬层深度的预测方法 |
CN105586476B (zh) * | 2016-01-19 | 2018-01-30 | 上海交通大学 | 基于数值模拟技术的高强钢淬火工艺优化方法 |
CN105586476A (zh) * | 2016-01-19 | 2016-05-18 | 上海交通大学 | 基于数值模拟技术的高强钢淬火工艺优化方法 |
CN105930605A (zh) * | 2016-04-29 | 2016-09-07 | 河海大学常州校区 | 感应淬火处理轴零件的淬硬层深度测量方法 |
CN107169166A (zh) * | 2017-04-17 | 2017-09-15 | 燕山大学 | 一种基于ansys仿真平台的大齿圈移动式感应加热的计算方法 |
CN107169166B (zh) * | 2017-04-17 | 2020-02-25 | 燕山大学 | 一种基于ansys仿真平台的大齿圈移动式感应加热的计算方法 |
CN107245557A (zh) * | 2017-06-21 | 2017-10-13 | 河海大学常州校区 | 基于ttt曲线预测轴件感应淬火淬硬层组织和硬度方法 |
CN107245552A (zh) * | 2017-06-21 | 2017-10-13 | 河海大学常州校区 | 基于cct曲线预测轴件感应淬火淬硬层组织和硬度方法 |
CN107907566A (zh) * | 2017-12-11 | 2018-04-13 | 中钢集团邢台机械轧辊有限公司 | 一种预测金属材料激光淬硬深度的试验方法 |
CN112084603A (zh) * | 2020-09-17 | 2020-12-15 | 泰尔重工股份有限公司 | 重载万向轴叉头淬火加热工艺参数的获取方法 |
CN113361159A (zh) * | 2021-05-31 | 2021-09-07 | 西安建筑科技大学 | 一种喷嘴射流冲击淬火的移动平板温度场仿真方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104331574A (zh) | 基于ansys有限元平台的感应淬火淬硬层深度的预测方法 | |
Gao et al. | Numerical and experimental analysis of 3D spot induction hardening of AISI 1045 steel | |
Fu et al. | Study on induction heating of workpiece before gear rolling process with different coil structures | |
Steinbeiss et al. | Method for optimizing the cooling design of hot stamping tools | |
Wen et al. | Study on mobile induction heating process of internal gear rings for wind power generation | |
CN105586476B (zh) | 基于数值模拟技术的高强钢淬火工艺优化方法 | |
Zhao et al. | Comparison of the methods for calculating the interfacial heat transfer coefficient in hot stamping | |
CN105930605A (zh) | 感应淬火处理轴零件的淬硬层深度测量方法 | |
CN105653770A (zh) | 一种连续感应淬火淬硬层深度的预测方法 | |
CN112149333B (zh) | 轴承滚道激光-感应复合淬火工艺参数优化方法 | |
CN105414205B (zh) | 一种基于plc的钢板温度在线预测方法 | |
CN105631111A (zh) | 一种台阶轴感应淬火马氏体分布的预测方法 | |
Bae et al. | Derivation of simplified formulas to predict deformations of plate in steel forming process with induction heating | |
Fortunato et al. | A comprehensive model for laser hardening of carbon steels | |
Anusha et al. | Finite element analysis and experimental validation of high-speed laser surface hardening process | |
Sun et al. | Thermo-mechanical coupled analysis of hot ring rolling process | |
Hong et al. | Metallographic method for temperature measurement: Reconstruction of grinding temperature field based on critical austenitizing depth and cyclic feedback algorithm | |
CN105838869B (zh) | 一种钢板淬火炉加热工艺在线调整方法 | |
Zhu et al. | A simplified heat source model for thick plate bending via high-frequency induction line heating | |
Wang et al. | Estimation of heat transfer coefficient and phase transformation latent heat by modified pattern search method | |
Chen et al. | Heat transfer modeling of an annular on-line spray water cooling process for electric-resistance-welded steel pipe | |
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 | |
CN107245557A (zh) | 基于ttt曲线预测轴件感应淬火淬硬层组织和硬度方法 | |
CN102799723B (zh) | 一种马弗炉加热能力计算及预测仿真方法 | |
Kotlan et al. | Combined heat treatment of metal materials |
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 |
Application publication date: 20150204 |
|
RJ01 | Rejection of invention patent application after publication |