CN108959793A - 一种伞齿轮双频分段循环感应加热数值模拟方法 - Google Patents

一种伞齿轮双频分段循环感应加热数值模拟方法 Download PDF

Info

Publication number
CN108959793A
CN108959793A CN201810767541.5A CN201810767541A CN108959793A CN 108959793 A CN108959793 A CN 108959793A CN 201810767541 A CN201810767541 A CN 201810767541A CN 108959793 A CN108959793 A CN 108959793A
Authority
CN
China
Prior art keywords
induction heating
frequency induction
frequency
bevel gear
cycle
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.)
Granted
Application number
CN201810767541.5A
Other languages
English (en)
Other versions
CN108959793B (zh
Inventor
韩毅
肖瑶
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changsha Hanghua Electronic Technology Co ltd
Dragon Totem Technology Hefei Co ltd
Original Assignee
Yanshan University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Yanshan University filed Critical Yanshan University
Priority to CN201810767541.5A priority Critical patent/CN108959793B/zh
Publication of CN108959793A publication Critical patent/CN108959793A/zh
Application granted granted Critical
Publication of CN108959793B publication Critical patent/CN108959793B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • 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)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Induction Heating (AREA)
  • Heat Treatment Of Articles (AREA)

Abstract

一种伞齿轮双频分段循环感应加热数值模拟方法，具体步骤如下：步骤1、模型参数确定；步骤2、几何模型创建；步骤3、计算中频感应加热；步骤4、计算高频感应加热；步骤5、设置双频感应加热过程。本发明对中、高频双频感应加热进行仿真模拟，三次嵌套使用DO循环编写APDL程序，实现双频率的分段循环加热，循环次数的增加提高了加热均匀性；可根据实际加热情况修改物性参数和工艺参数，匹配不同尺寸和材料的伞齿轮的加热循环次数和时间比，实现更好加热效果。

Description

一种伞齿轮双频分段循环感应加热数值模拟方法
技术领域
本发明涉及热处理领域，尤其是一种伞齿轮双频分段循环感应加热数值模拟方法。
背景技术
伞齿轮作为机械传动重要零部件，其表面接触疲劳强度和耐磨性影响正常服役以及工作寿命。利用电磁感应技术对伞齿轮表面加热，具有快速、高效、节能的特点，但同时由于伞齿轮齿面复杂，感应加热时存在尖端效应，齿根和齿顶难以同时实现硬化。利用双频率电流对伞齿轮加热可以有效改善热处理后伞齿轮表面硬化层分布均匀性。
双频感应加热涉及复杂多场耦合分析过程，采用实验的方法对其研究可以得到真实可靠的加热结果，但是耗费时间及物力较大，不便于探究参数改变时对双频加热的量化影响规律。随着计算机技术的飞速发展，数值模拟为这类影响因素复杂、实际加工耗能、成本高的研究提供了有效的分析手段。目前关于双频感应加热的数值模拟分析研究较少，已有的同步双频感应加热模拟方法，是将简谐载荷指定为有时间历程的载荷函数，从而进行瞬态动力学分析；已有的异步双频感应加热模拟方法，是先进行中频感应加热，结束后再进行高频感应加热。
发明内容
本发明的目的是提供一种利用“绕齿形”感应线圈加热的伞齿轮双频分段循环感应加热数值模拟方法。
本发明基于ANSYS软件进行模拟，利用顺序耦合物理环境法模拟双频感应加热，中、高频加热循环实施。本发明通过以下技术方案实现：
伞齿轮双频分段循环感应加热数值模拟方法的具体步骤如下：
步骤1、模型参数确定：根据伞齿轮双频感应加热实际生产过程确定模型建立的几何参数，提取实际加热时感应线圈内的电流载荷，确定工艺参数；
步骤2、几何模型创建：根据步骤1确定的几何参数建立几何模型以及空气模型，为所建立模型赋予材料物理性能参数；划分网格，对模型进行单元离散化；
步骤3、计算中频感应加热：在线圈上施加中频感应加热条件下的载荷，限定边界条件；读入中频电磁场物理文件，创建电磁分析物理环境；创建热分析物理环境；执行计算程序，采用DO循环计算中频感应加热过程中电磁-热耦合；存储结果数据；
步骤4、计算高频感应加热：在线圈上施加高频感应加热条件下的载荷，限定边界条件；以中频感应加热后的温度场分布作为高频感应加热的初始热源，读入高频电磁场物理文件，创建电磁分析物理环境；创建热分析物理环境；执行计算程序，采用DO循环计算高频感应加热过程中电磁-热耦合；存储结果数据；
步骤5、设置双频感应加热过程：采用DO循环计算中频感应加热和高频感应加热过程，设置循环次数，中、高频感应加热时间比，总加热时间；完成加热，将最终计算结果存储到数据库。
在步骤1中，所述模型建立的几何参数，是指伞齿轮尺寸参数和“绕齿形”感应线圈尺寸参数。
在步骤3中，所述中频感应加热条件下的载荷是指电流频率为6kHZ，电流密度为1.5×108A/m2。
在步骤4中，所述高频感应加热条件下的载荷是指电流频率为100kHZ，电流密度为6.0×107A/m2。
在步骤5中，所述循环次数为30次，中、高频感应加热时间比为2:1，总加热时间为45s；所述循环计算，是指每次循环包括中频感应加热1s和高频感应加热0.5s，并且每次循环都是基于上一次计算结果进行的。
本发明与现有技术相比具有如下优点：
(1)对中、高频双频感应加热进行仿真模拟，三次嵌套使用DO循环编写APDL程序，实现双频率的分段循环加热，循环次数的增加提高了加热均匀性；
(2)可根据实际加热情况修改物性参数和工艺参数，匹配不同尺寸和材料的伞齿轮的加热循环次数和时间比，实现更好加热效果。
附图说明
图1为本发明的双频感应加热模拟流程图；
图2为本发明一个实例的装配示意图；
图3为本发明的几何建模示意图；
图4为本发明的网格划分示意图；
图5为本发明的双频感应加热温度模拟云图。
具体实施方式
在图1至图4所示的本发明的示意简图中，本发明针对伞齿轮双频感应加热过程进行仿真，在保证对研究目标影响作用不大的前提下，对模拟过程作出如下假设及处理：(1)伞齿轮密度不变；(2)在温度场分析过程中，不考虑感应线圈、空气等的温度变化：感应线圈采用紫铜材料制作，电阻小，因此电阻热影响很小，并且感应加热时间极短，线圈中通入冷却水，线圈上热量变化对工件加热影响很小，可以忽略；(3)利用“绕齿形”感应线圈加热伞齿轮时，主要是分析双频感应加热对两齿齿间的加热作用效果，因此建立模型时为减少计算量，提高计算效率，只选取伞齿轮两齿以及齿间部分感应线圈进行建模。
伞齿轮双频分段循环感应加热数值模拟方法的具体步骤如下：
步骤1、模型参数确定，根据伞齿轮双频感应加热实际生产过程确定模型建立的几何参数，根据实际加热时感应线圈内的电流载荷确定工艺参数，如图2所示为加工实例的装配示意图；
步骤2、几何模型创建，根据步骤1确定的几何参数建立几何模型以及空气模型，如图3所示为几何模型示意图，为所建立模型赋予材料物理性能参数；划分网格，如图4所示为网格划分示意图；
步骤3、计算中频感应加热，在线圈上施加中频感应加热条件下电流频率为6kHZ，电流密度为1.5×108A/m2，限定边界条件；读入中频电磁场物理文件，创建电磁分析物理环境；创建热分析物理环境；执行计算程序，采用DO循环计算中频感应加热过程中电磁-热耦合，设置中频感应加热时间为1s；计算结束，存储结果数据；
步骤4、计算高频感应加热，在线圈上施加高频感应加热条件下100kHZ，电流密度为6.0×107A/m2，限定边界条件；以中频感应加热后的温度场分布作为高频感应加热的初始热源，读入高频电磁场物理文件，创建电磁分析物理环境；创建热分析物理环境；执行计算程序，采用DO循环计算高频感应加热过程中电磁-热耦合，设置高频感应加热时间为0.5s；计算结束，存储结果数据；
步骤5、设置双频感应加热过程，具体包括：采用DO循环计算中频感应加热和高频感应加热过程，设置循环次数为30，中、高频感应加热时间比为2:1，总加热时间为45s，其中每次循环包括中频感应加热1s和高频感应加热0.5s，并且每次循环都是基于上一次计算结果进行的；完成加热，将最终计算结果存储到数据库。

Claims (5)

1.一种伞齿轮双频分段循环感应加热数值模拟方法，其特征在于：具体步骤如下：
步骤1、模型参数确定：根据伞齿轮双频感应加热实际生产过程确定模型建立的几何参数，提取实际加热时感应线圈内的电流载荷，确定工艺参数；
步骤2、几何模型创建：根据步骤1确定的几何参数建立几何模型以及空气模型，为所建立模型赋予材料物理性能参数；划分网格，对模型进行单元离散化；
步骤3、计算中频感应加热：在线圈上施加中频感应加热条件下的载荷，限定边界条件；读入中频电磁场物理文件，创建电磁分析物理环境；创建热分析物理环境；执行计算程序，采用DO循环计算中频感应加热过程中电磁-热耦合；存储结果数据；
步骤4、计算高频感应加热：在线圈上施加高频感应加热条件下的载荷，限定边界条件；以中频感应加热后的温度场分布作为高频感应加热的初始热源，读入高频电磁场物理文件，创建电磁分析物理环境；创建热分析物理环境；执行计算程序，采用DO循环计算高频感应加热过程中电磁-热耦合；存储结果数据；
步骤5、设置双频感应加热过程：采用DO循环计算中频感应加热和高频感应加热过程，设置循环次数，中、高频感应加热时间比，总加热时间；完成加热，将最终计算结果存储到数据库。
2.根据权利要求1所述的一种伞齿轮双频分段循环感应加热数值模拟方法，其特征在于：在步骤1中，模型建立的几何参数，是指伞齿轮尺寸参数和“绕齿形”感应线圈尺寸参数。
3.根据权利要求1所述的一种伞齿轮双频分段循环感应加热数值模拟方法，其特征在于：在步骤3中，中频感应加热条件下的载荷是指电流频率为6kHZ，电流密度为1.5×108A/m2。
4.根据权利要求1所述的一种伞齿轮双频分段循环感应加热数值模拟方法，其特征在于：在步骤4中，所述高频感应加热条件下的载荷是指电流频率为100kHZ，电流密度为6.0×107A/m2。
5.根据权利要求1所述的一种伞齿轮双频分段循环感应加热数值模拟方法，其特征在于：在步骤5中，循环次数为30次，中、高频感应加热时间比为2:1，总加热时间为45s；循环计算，是指每次循环包括中频感应加热1s和高频感应加热0.5s，并且每次循环都是基于上一次计算结果进行的。
CN201810767541.5A 2018-07-13 2018-07-13 一种伞齿轮双频分段循环感应加热数值模拟方法 Active CN108959793B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810767541.5A CN108959793B (zh) 2018-07-13 2018-07-13 一种伞齿轮双频分段循环感应加热数值模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810767541.5A CN108959793B (zh) 2018-07-13 2018-07-13 一种伞齿轮双频分段循环感应加热数值模拟方法

Publications (2)

Publication Number Publication Date
CN108959793A true CN108959793A (zh) 2018-12-07
CN108959793B CN108959793B (zh) 2023-01-17

Family

ID=64483149

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810767541.5A Active CN108959793B (zh) 2018-07-13 2018-07-13 一种伞齿轮双频分段循环感应加热数值模拟方法

Country Status (1)

Country Link
CN (1) CN108959793B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111400944A (zh) * 2020-03-05 2020-07-10 燕山大学 一种基于动态坐标系的斜齿轮螺旋式电流加载方法
CN112877508A (zh) * 2021-01-12 2021-06-01 东北大学 一种基于异形感应器的异步双频感应加热数值模拟方法
CN113151658A (zh) * 2021-04-19 2021-07-23 东北大学 基于异步双频感应加热均匀自控制齿轮热处理装置及方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030222656A1 (en) * 2001-12-20 2003-12-04 Phillips Alan D. On-line oil condition sensor system for rotating and reciprocating machinery
CN1973578A (zh) * 2004-04-21 2007-05-30 感应加热有限公司 通过感应加热对工件进行多频加热处理
US20100088033A1 (en) * 2008-10-03 2010-04-08 Baker Hughes Incorporated Data acquisition and processing for invasion profile and gas zone analysis with nmr dual or multiple interecho spacing time logs
CN105740536A (zh) * 2016-01-28 2016-07-06 燕山大学 对直缝焊管中频热处理过程进行动态仿真的模拟方法
US20160339494A1 (en) * 2014-02-04 2016-11-24 Toshiba Mitsubishi-Electric Industrial Systems Corporation Temperature control apparatus of hot-rolling mill
CN107769603A (zh) * 2017-10-30 2018-03-06 西安理工大学 二极管箝位型逆变器双频感应加热电源的解耦控制方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030222656A1 (en) * 2001-12-20 2003-12-04 Phillips Alan D. On-line oil condition sensor system for rotating and reciprocating machinery
CN1973578A (zh) * 2004-04-21 2007-05-30 感应加热有限公司 通过感应加热对工件进行多频加热处理
US20100088033A1 (en) * 2008-10-03 2010-04-08 Baker Hughes Incorporated Data acquisition and processing for invasion profile and gas zone analysis with nmr dual or multiple interecho spacing time logs
US20160339494A1 (en) * 2014-02-04 2016-11-24 Toshiba Mitsubishi-Electric Industrial Systems Corporation Temperature control apparatus of hot-rolling mill
CN105740536A (zh) * 2016-01-28 2016-07-06 燕山大学 对直缝焊管中频热处理过程进行动态仿真的模拟方法
CN107769603A (zh) * 2017-10-30 2018-03-06 西安理工大学 二极管箝位型逆变器双频感应加热电源的解耦控制方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
SCHWENK,M等: "Experimentally validated residual stresses and distortion prediction for dualfrequency induction hardening", 《INTERNATIONAL JOURNAL OF APPLIED ELECTROMAGENTICS AND MECHANICS》 *
易小开等: "输送管道中频感应加热的磁-热耦合仿真", 《金属热处理》 *
赵天旭: "高频直缝焊管动态中频感应加热过程数值分析及实验研究", 《中国优秀硕士学位论文全文数据库 (工程科技Ⅰ辑)》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111400944A (zh) * 2020-03-05 2020-07-10 燕山大学 一种基于动态坐标系的斜齿轮螺旋式电流加载方法
CN111400944B (zh) * 2020-03-05 2022-03-25 燕山大学 一种基于动态坐标系的斜齿轮螺旋式电流加载方法
CN112877508A (zh) * 2021-01-12 2021-06-01 东北大学 一种基于异形感应器的异步双频感应加热数值模拟方法
CN112877508B (zh) * 2021-01-12 2022-02-01 东北大学 一种基于异形感应器的异步双频感应加热数值模拟方法
CN113151658A (zh) * 2021-04-19 2021-07-23 东北大学 基于异步双频感应加热均匀自控制齿轮热处理装置及方法
CN113151658B (zh) * 2021-04-19 2022-07-22 东北大学 基于异步双频感应加热均匀自控制齿轮热处理装置及方法

Also Published As

Publication number Publication date
CN108959793B (zh) 2023-01-17

Similar Documents

Publication Publication Date Title
CN105740536B (zh) 对直缝焊管中频热处理过程进行动态仿真的模拟方法
CN108959793A (zh) 一种伞齿轮双频分段循环感应加热数值模拟方法
CN105608273B (zh) 一种基于cfd软件优化动力电池包温度分布的系统
CN103699744B (zh) 一种基于有限元法的风电主控柜体热管理分析及优化方法
CN103970591B (zh) 一种Matlab‑ATP集成自动化调用方法
CN104020763B (zh) 输入输出接口的自动测试方法及系统
CN103309727B (zh) 一种基于ansys-cfx软件的叶片泵空化数值模拟自动运行方法
CN107169166A (zh) 一种基于ansys仿真平台的大齿圈移动式感应加热的计算方法
CN112231945A (zh) 一种基于star CCM+与Amesim的动力电池系统热扩散联合仿真方法
Park et al. Optimization of the in-line induction heating process for hot forging in terms of saving operating energy
CN109920036A (zh) 一种基于bim技术的弧形钢结构施工方法
CN110096812B (zh) 一种基于cfd平台的离心泵汽蚀性能自动模拟方法
CN105260572A (zh) 一种风机叶片建模软件调用系统及其调用方法
CN113221422A (zh) 基于非线性流体数据信息处理技术的流体仿真模拟方法
CN106557597A (zh) 一种热虚拟试验方法
CN114580222B (zh) 用于电子器件的高速瞬态热仿真方法及系统
CN103020337A (zh) 利用参数化建模控制电加热炉熔融矿石的方法
CN107330204A (zh) 闭式叶轮焊接工艺的有限元模拟分析方法及装置
CN114117869B (zh) 一种工件运动状态下多频感应淬火加热的有限元模拟方法
CN110955999B (zh) 轮胎三维温度场仿真方法
CN110162897B (zh) 一种大口径弯管加热导磁体优化方法
Zhang et al. Research on digital twin model of energy equipment
CN110009746B (zh) 反应堆燃料组件的带有边界层的六面体网格自动生成方法
CN107729699B (zh) 一种基于matlab的增压器涡轮叶轮设计优化方法
CN104570737A (zh) 基于随机算法的pid控制器参数选定方法

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
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230630

Address after: 230000 floor 1, building 2, phase I, e-commerce Park, Jinggang Road, Shushan Economic Development Zone, Hefei City, Anhui Province

Patentee after: Dragon totem Technology (Hefei) Co.,Ltd.

Address before: 066004 No. 438 west section of Hebei Avenue, seaport District, Hebei, Qinhuangdao

Patentee before: Yanshan University

Effective date of registration: 20230630

Address after: Plant 7, No. 27, Leifeng Avenue, Wangcheng Economic and Technological Development Zone, Changsha, Hunan 410000

Patentee after: Changsha HANGHUA Electronic Technology Co.,Ltd.

Address before: 230000 floor 1, building 2, phase I, e-commerce Park, Jinggang Road, Shushan Economic Development Zone, Hefei City, Anhui Province

Patentee before: Dragon totem Technology (Hefei) Co.,Ltd.