CN106682346A - 一种基于cfd软件对复杂构件气淬系统的优化方法 - Google Patents
一种基于cfd软件对复杂构件气淬系统的优化方法 Download PDFInfo
- Publication number
- CN106682346A CN106682346A CN201710007865.4A CN201710007865A CN106682346A CN 106682346 A CN106682346 A CN 106682346A CN 201710007865 A CN201710007865 A CN 201710007865A CN 106682346 A CN106682346 A CN 106682346A
- Authority
- CN
- China
- Prior art keywords
- quenching system
- air quenching
- cfd software
- method based
- optimization method
- 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
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)
- Control Of Heat Treatment Processes (AREA)
Abstract
本发明公开了一种基于CFD软件对复杂构件气淬系统的优化方法，包括建立构件气淬系统固体计算域和流体计算域几何模型；几何模型网格划分；网格模型导入CFD软件中，设置并求解温度场；通过修改相关参数，重复上述步骤，进行多次数值模拟；将所述多次数值模拟所对应的温度场数据整理并建立数据库，以构件内部温度差最小为目标，通过数据库自动寻优，选择最优工艺参数组合，以保证构件内部温度均匀性。本发明通过计算机模拟技术对各种气淬工况下构件内部的温度场进行预测，以构件内部温度差最小为目标寻优，以保证构件内部温度均匀性，为复杂构件气淬系统优化提供一条快捷途径。
Description
技术领域
本发明涉及一种复杂构件热处理工艺优化方法，具体涉及一种基于CFD软件对复杂构件气淬系统的优化方法。
背景技术
淬火是金属加工的一道重要工序，通过调整和控制冷却介质参数以控制工件的温度场、组织场和应力场以及应变场，使得工件达到相应的组织要求和性能，并尽可能减小不均匀场导致的残余应力及残余应变，以提高金属工件质量。在淬火过程中，由于冷却不均匀导致金属工件温度分布不均匀，从而组织转变不均匀，产生热应力和组织应力，影响工件的机械性能，使得寿命下降，甚至造成工件产生裂纹或开裂而报废，直接影响工件质量。对于淬火的复杂构件，因其几何结构不均匀，工件截面厚度差异大，冷却速率不均匀导致的温度场不均匀问题更加突出。生产实践中，淬火过程是热处理中返修率和废品率最高、最难控制的环节，引起的经济损失不可小觑。因此，如何减少淬火过程中不均匀各场产生的残余应力是相关工作者最为关心的问题，而如何调整和控制工件的冷却速率，以改善温度场均匀性则是这一问题的最基础一步。
淬火过程中，工件与冷却介质间换热情况复杂，受工件形状、表面温度、所处流场及冷却介质的温度、流速等多种因素影响，难以定量分析。目前，冷却介质与工件之间的冷却速率多为定性分析或实验测量，不能满足精确、高效的控制要求。
淬火过程是一个多场耦合的复杂过程，理论求解非常困难，对复杂构件更不可能直接求解；相关实验涉及高温作业和在线测量，参数测量也相当困难；为了解不同冷却方案的效果时，需要进行大量实验，成本浪费太大；同时，由于构件冷却过程中的综合传热系数受构件的表面特性、温度以及冷却介质的流速分布等参数的影响，实验结果的普适性较差。基于上述原因，目前的淬火过程控制多为定性分析、经验调节，理论薄弱，技术落后。随着计算机技术发展，用数值模拟代替实验成为热点。对于复杂构件的气淬过程，将气体流动、传热过程与固体构件的传热过程相耦合，将CFD技术与最优化方法相融合，对气淬装置结构及其操作进行优化，可有效提高复杂构件气淬过程中冷却速率均匀性，减小构件中的残余应力，从而提高复杂构件的质量。
发明内容
本发明的目的在于提供一种基于CFD软件对复杂构件气淬系统的优化方法，利用CFD软件模拟得出复杂构件在各种工况下气体淬火的温度场，将各工况对应的温度场数据整理并建立数据库，以构件内部温度差最小为目标，通过数据库自动寻优，选择最优工艺参数组合，以保证构件内部温度均匀性。
为解决上述技术问题，本发明所提供的技术方案为：
一种基于CFD软件对复杂构件气淬系统的优化方法，包括如下步骤：
1)建立几何模型：根据气淬系统基本尺寸参数，建立气淬系统几何模型，包括定义复杂构件区域为固体计算域及流体区域为流体计算域；
2)划分网格：分别将固体计算域和流体计算域的几何模型导入网格划分软件中进行网格划分；
3)将步骤2)中得到的网格划分后的几何模型导入CFD软件中，定义基本求解器、启动能量方程、指定湍流模型、定义材料的物性参数、边界条件参数及初始化后，开始数值模拟计算；
4)判断是否收敛，如果残差小于设定值，则收敛并进行下一步骤；如果不收敛，对步骤2)中固体计算域和流体计算域的几何模型重新进行网格划分，优化网格后，重复步骤1)～3)直至收敛为止；
5)重复步骤1)～4)，进行多次数值模拟；每一次数值模拟过程中改变步骤1)中基本尺寸参数和步骤3)中边界条件参数；
6)将所述多次数值模拟所对应的温度场数据整理并建立数据库，以构件内部温度差最小为目标，通过数据库自动寻优，选择最优工艺参数组合，以保证构件内部温度均匀性。
所述步骤1)中基本尺寸参数包括构件的几何尺寸、喷嘴的位置、大小及数量。
所述步骤2)中网格划分优先选用结构化网格。
所述步骤3)中基本求解器选用压力基求解器及非稳态。
所述步骤3)中湍流模型根据雷诺数选择K-e模型。
所述步骤3)中材料的物理性质为构件的材质、密度、导热系数、比热容及冷却介质的材质、导热系数、密度、比热容、黏度。
所述步骤5)中改变的基本尺寸参数包括喷嘴的位置、大小及数量，改变的边界条件参数包括冷却介质速度、温度。
所述步骤3)中边界条件参数包括：
速度入口的冷却介质速度、温度、湍流强度、水力直径；
压力出口的压力、湍流强度、水力直径；
构件壁面耦合热边界及喷嘴壁面热边界；
其中，速度入口为冷却介质出口，压力出口为流体计算域的外部空间边界。
步骤4)中残差小于10^-6时收敛。
技术效果：本发明将CFD技术与最优化方法相融合，通过计算机模拟技术对各种气淬工况下构件内部的温度场进行预测，以构件内部温度差最小为目标寻优，并选择最优工艺参数组合，以保证构件内部温度均匀性，降低经验化调整淬火工艺参数造成的损失，为复杂构件气淬系统优化提供一条快捷途径。
附图说明
图1为本发明初始气淬系统示意图；
图2为本发明温度分布示意图；
图3为本发明最优气淬系统网格；
图4为本发明最优气淬系统示意图。
具体实施方式
下面将结合附图和具体实施方式对本发明作进一步的说明。
一种基于CFD软件对复杂构件气淬系统的优化方法，包括如下步骤：
1)建立几何模型：根据气淬系统基本尺寸参数，建立气淬系统的几何模型，包括定义固体计算域和流体计算域；根据构件基本尺寸参数利用建模软件建立固体计算域并选取构件周围220*310mm2空间为流体计算域，得到构件气淬模型；
如图1所示，构件1由三段组成，两端为Φ10×100mm圆柱，中间为Φ20×100mm圆柱，喷嘴2为Φ10mm圆孔，数量为3个，位于构件1的上方，喷嘴口处距离工件100mm.由于圆柱工件为对称结构,为简化计算,建立二维模型即可。固体计算域的几何模型大小为复杂构件实际尺寸，流体计算域的几何模型大小为310mm×220mm。
2)划分网格：将固体计算域和流体计算域的几何模型导入网格划分软件中进行网格划分；网格划分选用结构化Map网格，网格最小尺寸为1*1mm2；
3)将步骤2)中得到的网格划分后的固体计算域和流体计算域几何模型导入CFD软件中，定义基本求解器(压力基求解器及非稳态)、启动能量方程、指定湍流模型(根据雷诺数选择K-e模型)、定义材料的物性参数(设置构件的材质、密度、比热容、导热系数、初始温度及冷却介质的材质、密度、比热容、黏度、导热系数)、边界条件参数(基本工况有三个喷嘴，则分别设置三个速度入口冷却介质的速度、温度、湍流强度和水力直径；设置压力出口的压力、湍流强度和水力直径；构件壁面耦合热边界条件为couple，喷嘴壁面热边界条件为恒温；其中，速度入口为冷却介质出口，压力出口为为流体计算域的外部空间边界，壁面条件为喷嘴壁面及构件表面)及初始化(定义构件初始温度)后，开始数值模拟计算；
4)判断是否收敛，残差小于10^-6则收敛则，可得到构件的温度场云图(图2(a))，从图中的温度场分布能了解到构件中最高温度及最低温度，为优化使得温度场均匀而提供判据，并进行下一步骤；如果不收敛，对步骤2)中固体计算域和流体计算域的几何模型重新进行网格划分，对结构复杂部位进行网格细化后，重复步骤1)～3)直至收敛为止；
5)重复步骤1)～4)，进行多次数值模拟；每一次数值模拟过程中改变步骤1)中喷嘴2的位置、大小及数量和步骤3)中冷却介质的速度和温度，得到多次数值模拟所对应的温度场数据；
6)将所述多次数值模拟所对应的温度场数据整理并建立数据库，以构件内部温度差最小为目标，通过数据库自动寻优，选择最优工艺参数组合，以保证构件内部温度均匀性，得到最优的温度场(图2(b))及其对应的网格(图3)和气淬系统(图4)。
Claims (10)
1.一种基于CFD软件对复杂构件气淬系统的优化方法，其特征在于，包括如下步骤：
1)建立几何模型：根据气淬系统基本尺寸参数，建立气淬系统几何模型，包括定义复杂构件区域为固体计算域及流体区域为流体计算域；
2)划分网格：将固体计算域和流体计算域的几何模型导入网格划分软件中进行网格划分；
3)将步骤2)中得到的网格划分后的几何模型导入CFD软件中，定义基本求解器、启动能量方程、指定湍流模型、定义材料的物性参数、边界条件参数及初始化后，开始数值模拟计算；
4)判断是否收敛，如果残差小于设定值，则收敛并进行下一步骤；如果不收敛，对步骤2)中固体计算域和流体计算域的几何模型重新进行网格划分，优化网格后，重复步骤1)～3)直至收敛为止；
5)重复步骤1)～4)，进行多次数值模拟；每一次数值模拟过程中改变步骤1)中基本尺寸参数和步骤3)中边界条件参数；
6)将所述多次数值模拟所对应的温度场数据整理并建立数据库，以构件内部温度差最小为目标，通过数据库自动寻优，选择最优工艺参数组合，以保证构件内部温度均匀性。
2.根据权利要求1所述的基于CFD软件对复杂构件气淬系统的优化方法，其特征在于：所述步骤1)中基本尺寸参数包括构件的几何尺寸、喷嘴的位置、大小及数量。
3.根据权利要求1所述的基于CFD软件对复杂构件气淬系统的优化方法，其特征在于：所述步骤2)中网格划分优先选用结构化网格。
4.根据权利要求1所述的基于CFD软件对复杂构件气淬系统的优化方法，其特征在于：所述步骤3)中基本求解器选用压力基及非稳态求解器。
5.根据权利要求1所述的基于CFD软件对复杂构件气淬系统的优化方法，其特征在于：步骤3)中湍流模型根据雷诺数选择K-e模型。
6.根据权利要求1所述的基于CFD软件对复杂构件气淬系统的优化方法，其特征在于：所述步骤3)中材料的物性参数包括：构件的材质、密度、比热容、导热系数及冷却介质的种类、密度、比热容、导热系数、黏度。
7.根据权利要求1所述的基于CFD软件对复杂构件气淬系统的优化方法，其特征在于：所述步骤3)中边界条件参数包括：
速度入口的冷却介质速度、温度、湍流强度、水力直径；
压力出口的压力、湍流强度、水力直径；
构件壁面热边界及喷嘴壁面热边界；
其中，速度入口为冷却介质出口，压力出口为流体计算域的外部空间边界。
8.根据权利要求7所述的基于CFD软件对复杂构件气淬系统的优化方法，其特征在于：所述构件壁面热边界条件为couple，喷嘴壁面热边界为恒温热边界。
9.根据权利要求1所述的基于CFD软件对复杂构件气淬系统的优化方法，其特征在于：步骤4)中残差小于10^-6时收敛。
10.根据权利要求1-9中任一项所述的基于CFD软件对复杂构件气淬系统的优化方法，其特征在于：所述步骤5)中，改变的基本尺寸参数包括喷嘴的位置、大小及数量，改变的边界条件参数包括速度入口的冷却介质速度和温度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710007865.4A CN106682346A (zh) | 2017-01-05 | 2017-01-05 | 一种基于cfd软件对复杂构件气淬系统的优化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710007865.4A CN106682346A (zh) | 2017-01-05 | 2017-01-05 | 一种基于cfd软件对复杂构件气淬系统的优化方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106682346A true CN106682346A (zh) | 2017-05-17 |
Family
ID=58849971
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710007865.4A Pending CN106682346A (zh) | 2017-01-05 | 2017-01-05 | 一种基于cfd软件对复杂构件气淬系统的优化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106682346A (zh) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107341291A (zh) * | 2017-06-09 | 2017-11-10 | 中国电子科技集团公司第四十八研究所 | 一种电热辊道炉结构优化方法 |
| CN109033546A (zh) * | 2018-07-03 | 2018-12-18 | 重庆川仪调节阀有限公司 | 一种基于cfd的阀门传热模拟方法 |
| CN109543362A (zh) * | 2019-01-25 | 2019-03-29 | 南京航空航天大学 | 一种具有高效空气减阻功能的仿沙丘微结构及其优化设计方法 |
| CN109583131A (zh) * | 2018-12-18 | 2019-04-05 | 江苏理工学院 | 一种具有减阻效果的表面微结构尺寸参数的优化设计方法 |
| CN109858175A (zh) * | 2019-02-20 | 2019-06-07 | 苏州热工研究院有限公司 | 一种核电厂容控箱氮吹扫效率优化方法 |
| CN109863868A (zh) * | 2019-03-12 | 2019-06-11 | 中国农业科学院农业环境与可持续发展研究所 | 一种基于ai算法与工作需求的文丘里施肥器优化设计方法 |
| CN111027155A (zh) * | 2019-12-13 | 2020-04-17 | 上海市计量测试技术研究院 | 基于cfd技术的飞机空气循环制冷系统的模拟分析方法 |
| CN111143939A (zh) * | 2019-12-13 | 2020-05-12 | 上海市计量测试技术研究院 | 基于cfd技术的飞机空气循环制冷系统的模拟分析系统 |
| CN112069622A (zh) * | 2020-09-08 | 2020-12-11 | 北京航空航天大学 | 一种涡轮导向器叶片冷却结构的智能推荐系统及推荐方法 |
| CN113361159A (zh) * | 2021-05-31 | 2021-09-07 | 西安建筑科技大学 | 一种喷嘴射流冲击淬火的移动平板温度场仿真方法 |
| CN116451341A (zh) * | 2023-03-13 | 2023-07-18 | 中南大学 | 一种直升机减速器螺旋引流飞溅润滑数值模拟方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101220403A (zh) * | 2007-12-20 | 2008-07-16 | 山东建筑大学材料科学研究所 | 一种球墨铸铁贝氏体等温处理方法 |
| CN102043877A (zh) * | 2009-10-22 | 2011-05-04 | 通用汽车环球科技运作公司 | 用于预测淬火期间的传热系数的系统和方法 |
| CN104212969A (zh) * | 2014-09-18 | 2014-12-17 | 上海交通大学 | 一种基于数值模拟的钢管连续淬火过程控制方法 |
-
2017
- 2017-01-05 CN CN201710007865.4A patent/CN106682346A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101220403A (zh) * | 2007-12-20 | 2008-07-16 | 山东建筑大学材料科学研究所 | 一种球墨铸铁贝氏体等温处理方法 |
| CN102043877A (zh) * | 2009-10-22 | 2011-05-04 | 通用汽车环球科技运作公司 | 用于预测淬火期间的传热系数的系统和方法 |
| CN104212969A (zh) * | 2014-09-18 | 2014-12-17 | 上海交通大学 | 一种基于数值模拟的钢管连续淬火过程控制方法 |
Non-Patent Citations (2)
| Title |
|---|
| 王睿智: "变厚度盘件气体淬火系统的设计和仿真", 《中国优秀硕士学位论文全文数据库-工程科技Ⅰ辑》 * |
| 罗筠: "工件气体淬火过程数值模拟", 《热加工工艺》 * |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107341291A (zh) * | 2017-06-09 | 2017-11-10 | 中国电子科技集团公司第四十八研究所 | 一种电热辊道炉结构优化方法 |
| CN109033546A (zh) * | 2018-07-03 | 2018-12-18 | 重庆川仪调节阀有限公司 | 一种基于cfd的阀门传热模拟方法 |
| CN109033546B (zh) * | 2018-07-03 | 2023-10-31 | 重庆川仪调节阀有限公司 | 一种基于cfd的阀门传热模拟方法 |
| CN109583131A (zh) * | 2018-12-18 | 2019-04-05 | 江苏理工学院 | 一种具有减阻效果的表面微结构尺寸参数的优化设计方法 |
| CN109583131B (zh) * | 2018-12-18 | 2023-09-22 | 江苏理工学院 | 一种具有减阻效果的表面微结构尺寸参数的优化设计方法 |
| CN109543362A (zh) * | 2019-01-25 | 2019-03-29 | 南京航空航天大学 | 一种具有高效空气减阻功能的仿沙丘微结构及其优化设计方法 |
| CN109543362B (zh) * | 2019-01-25 | 2020-05-08 | 南京航空航天大学 | 一种具有高效空气减阻功能的仿沙丘微结构及其优化设计方法 |
| CN109858175B (zh) * | 2019-02-20 | 2023-06-09 | 苏州热工研究院有限公司 | 一种核电厂容控箱氮吹扫效率优化方法 |
| CN109858175A (zh) * | 2019-02-20 | 2019-06-07 | 苏州热工研究院有限公司 | 一种核电厂容控箱氮吹扫效率优化方法 |
| CN109863868A (zh) * | 2019-03-12 | 2019-06-11 | 中国农业科学院农业环境与可持续发展研究所 | 一种基于ai算法与工作需求的文丘里施肥器优化设计方法 |
| CN111143939A (zh) * | 2019-12-13 | 2020-05-12 | 上海市计量测试技术研究院 | 基于cfd技术的飞机空气循环制冷系统的模拟分析系统 |
| CN111027155A (zh) * | 2019-12-13 | 2020-04-17 | 上海市计量测试技术研究院 | 基于cfd技术的飞机空气循环制冷系统的模拟分析方法 |
| CN112069622A (zh) * | 2020-09-08 | 2020-12-11 | 北京航空航天大学 | 一种涡轮导向器叶片冷却结构的智能推荐系统及推荐方法 |
| CN113361159A (zh) * | 2021-05-31 | 2021-09-07 | 西安建筑科技大学 | 一种喷嘴射流冲击淬火的移动平板温度场仿真方法 |
| CN116451341A (zh) * | 2023-03-13 | 2023-07-18 | 中南大学 | 一种直升机减速器螺旋引流飞溅润滑数值模拟方法 |
| CN116451341B (zh) * | 2023-03-13 | 2023-11-21 | 中南大学 | 一种直升机减速器螺旋引流飞溅润滑数值模拟方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106682346A (zh) | 一种基于cfd软件对复杂构件气淬系统的优化方法 | |
| Xiong et al. | Finite element analysis and experimental validation of thermal behavior for thin-walled parts in GMAW-based additive manufacturing with various substrate preheating temperatures | |
| Ding et al. | Thermo-mechanical analysis of Wire and Arc Additive Layer Manufacturing process on large multi-layer parts | |
| US10520917B2 (en) | Method of simulating a shaping process | |
| CN104070075B (zh) | 一种热轧带钢层流冷却过程控制装置与方法 | |
| CN107066700B (zh) | 一种电弧增材制造有限元建模方法 | |
| CN109513924B (zh) | 一种激光选区熔化过程中表面粗糙度控制方法 | |
| CN108038266A (zh) | 一种选择性激光修复过程的数值模拟方法 | |
| CN107760830A (zh) | 一种大单重、大断面特厚钢板辊式淬火过程的控制方法 | |
| CN102094528A (zh) | 一种大体积混凝土冷却水管布置方法 | |
| CN102693336A (zh) | 一种预测大型管道焊接热循环参数的方法 | |
| Schreivogel et al. | Study of an optimized trench film cooling configuration using scale adaptive simulation and infrared thermography | |
| Silva et al. | Optimising the energy consumption on pultrusion process | |
| CN109182731A (zh) | 一种基于温度时间控制的高碳铬轴承钢连铸坯加热控制方法 | |
| CN104985298A (zh) | 一种预测旋转电弧低合金结构钢小角度焊接温度场的方法 | |
| Zhang et al. | A physical model and data-driven hybrid prediction method towards quality assurance for composite components | |
| CN114925585B (zh) | 超高强钢热冲压成型模具温度场交互式耦合数值模拟方法 | |
| CN110449587A (zh) | 基于3d打印和msc仿真的义齿热处理工艺方法 | |
| CN212310848U (zh) | 用于基于pbf的三维(3d)打印机的热处理设备 | |
| Jiang et al. | Effects of Heat Input on Morphology of Thin‐Wall Components Fabricated by Wire and Arc Additive Manufacturing | |
| CN106202686B (zh) | 一种涡轮盘等温模锻预成形坯料的多目标设计方法 | |
| CN105799175B (zh) | 一种维持工作缸整体温度均衡的装置及其控制方法 | |
| CN108959793B (zh) | 一种伞齿轮双频分段循环感应加热数值模拟方法 | |
| Hao et al. | 3-D Numerical analysis on heating process of loads within vacuum heat treatment furnace | |
| Yi et al. | Improving the curing cycle time through the numerical modeling of air flow in industrial continuous convection ovens |
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: 20170517 |
|
| RJ01 | Rejection of invention patent application after publication |