CN113627059B - 一种考虑相变热的大规格棒材空冷温度场计算方法 - Google Patents
一种考虑相变热的大规格棒材空冷温度场计算方法 Download PDFInfo
- Publication number
- CN113627059B CN113627059B CN202110888962.5A CN202110888962A CN113627059B CN 113627059 B CN113627059 B CN 113627059B CN 202110888962 A CN202110888962 A CN 202110888962A CN 113627059 B CN113627059 B CN 113627059B
- Authority
- CN
- China
- Prior art keywords
- phase change
- heat
- air cooling
- bar
- temperature field
- 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.)
- Active
Links
- 230000008859 change Effects 0.000 title claims abstract description 62
- 238000001816 cooling Methods 0.000 title claims abstract description 45
- 238000004364 calculation method Methods 0.000 title claims abstract description 13
- 238000000034 method Methods 0.000 claims abstract description 45
- 230000008569 process Effects 0.000 claims abstract description 24
- 238000004458 analytical method Methods 0.000 claims abstract description 15
- 230000001052 transient effect Effects 0.000 claims abstract description 8
- 238000004088 simulation Methods 0.000 claims description 14
- 230000000694 effects Effects 0.000 claims description 13
- 239000000463 material Substances 0.000 claims description 13
- 239000000243 solution Substances 0.000 claims description 12
- 229910000831 Steel Inorganic materials 0.000 claims description 9
- 239000010959 steel Substances 0.000 claims description 9
- 239000007787 solid Substances 0.000 claims description 7
- 230000000295 complement effect Effects 0.000 claims description 3
- 239000003637 basic solution Substances 0.000 claims description 2
- 239000012530 fluid Substances 0.000 claims description 2
- 238000007781 pre-processing Methods 0.000 claims description 2
- 102000015833 Cystatin Human genes 0.000 claims 1
- 108050004038 cystatin Proteins 0.000 claims 1
- 230000007547 defect Effects 0.000 abstract description 3
- 238000009826 distribution Methods 0.000 abstract description 2
- 230000009466 transformation Effects 0.000 description 9
- 230000005855 radiation Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 230000017525 heat dissipation Effects 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 229910052751 metal Inorganic materials 0.000 description 2
- 230000008520 organization Effects 0.000 description 2
- 230000036314 physical performance Effects 0.000 description 2
- 238000002076 thermal analysis method Methods 0.000 description 2
- 230000006032 tissue transformation Effects 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000009851 ferrous metallurgy Methods 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000002436 steel type Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
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)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种考虑相变热的大规格棒材空冷温度场计算方法，首先基于有限元分析软件实现对大规格棒材的精确建模、划分网格、建立边界条件及表面单元及空间节点的建立等工作；然后利用POST26后处理器对大规格棒材空冷过程的温度场进行瞬态分析，计算棒材芯部考虑相变热的温度变化曲线。本发明方法考虑了大规格棒材空冷过程中相变过程产生的热量对棒材空冷温度场的影响，克服了现有研究对棒材空冷温度场计算时忽略相变热的缺陷，不仅提高了空冷温度场的计算精度，而且降低了建模的复杂性和计算量，能够高效预测大规格棒材空冷过程的温度分布。
Description
技术领域
本发明涉及钢铁冶金技术领域，具体涉及热轧大规格棒材空冷温度场计算方法。
背景技术
在热轧大规格棒材(直径超过50mm的棒材)的生产中，一般采取空冷方式进行冷却。空冷过程中往往伴随着组织的转变，轧件温度的变化对组织转变起到了至关重要的作用。在组织转变过程中，相变潜热的释放造成相变过程中金属温度回升并且导致相变的持续时间变长，它使相变在一个波动范围小、相对稳定的温度区间内进行，类似于等温转变。空冷条件下的金属相变的这种等温转变特性，降低了相变过程中的冷却速度，影响组织的形成。温度影响组织，抛开相变潜热研究大规格棒材空冷过程的温度变化，很难准确预测棒材在空冷过程中的组织演变。
发明内容
为了实现对大规格棒材空冷过程中组织演变的调控，掌握相变热对大规格棒材空冷过程温度变化规律的影响是非常重要的。
针对目前对大规格棒材空冷过程中的温度场计算缺陷，本发明采用有限元模拟的方式，基于有限元仿真提出一种考虑相变热的大规格棒材空冷温度场计算方法，其特征在于包括以下步骤：
步骤1：根据棒材几何参数，在有限元分析软件ANSYS的前处理模块中建立棒材的三维几何模型，或在三维建模软件中建立棒材的三维几何模型，并将节点单元数据读入有限元分析软件ANSYS中；
步骤2：建立材料热物性参数模型；
步骤3：在有限元分析软件ANSYS中对棒材的三维几何模型进行网格划分；
步骤4：定义实体单元和表面效应单元类型及表面效应单元实常数；
步骤5：设置初始条件和边界条件；
步骤6：施加载荷及载荷选项；
步骤7：导入相变热模型；
利用有限元分析软件ANSYS计算空冷状态下相变潜热时，利用补热法进行模拟分析，将相变潜热作为一个内热源加入到有限元分析软件ANSYS分析过程中，即在相变区间内，通过内热源释放等效的热量从而替代相变产生的热量，具体相变热模型如下：
即：
式中T代表温度(℃)，x、y分别为中间坯厚度和宽度方向坐标，λ代表试验钢的导热系数(W/m·k)，ρ代表试验钢的密度(kg/m3)，Cp代表试验钢的比热容(J/kg·K)，qv代表内热源强度(W/m3)；
步骤8：进行模拟求解设置；
步骤9：设置计算时间步长和迭代步数，完成求解；
步骤10：进入时间-历程POST26后处理器，对所述相变热模型进行瞬态分析，获取棒材三维几何模型在冷却过程中任意时间点上的温度值，计算棒材芯部考虑相变热的温度变化曲线。
本发明提出了一种考虑相变热的大规格棒材空冷温度场计算方法，首先基于有限元分析软件实现对大规格棒材的精确建模、划分网格、建立边界条件及表面单元及空间节点的建立等工作；然后利用POST26后处理器对大规格棒材空冷过程的温度场进行瞬态分析。本发明方法考虑了大规格棒材空冷过程中相变过程产生的热量对棒材空冷温度场的影响，克服了现有研究对棒材空冷温度场计算时忽略相变热的缺陷，不仅提高了空冷温度场的计算精度，而且降低了建模的复杂性和计算量，能够高效预测大规格棒材空冷过程的温度分布。
附图说明
图1为本发明实施例考虑相变热的大规格棒材空冷温度场计算方法流程图；
图2为本发明实施例基于ANSYS软件构建的棒材三维几何模型示意图；
图3为本发明实施例棒材三维几何模型的网格划分示意图；
图4为本发明实施例三维热实体单元SOLID70节点示意图；
图5为物理试验结果与未考虑相变热的模拟输出结果对比图；
图6为物理试验结果与本发明实施例考虑相变热的模拟输出结果对比图。
具体实施方式
为使本发明要解决的技术问题、技术方案和优点更加清楚，下面将结合具体实施例进行详细描述。应当理解，此处所描述的具体实施实例仅仅用以解释本发明，并不用于限定本发明。
如图1所示，本实施例提供一种考虑相变热的大规格棒材空冷温度场计算方法，包括如下步骤：
步骤1：根据实际棒材几何参数，在有限元分析软件ANSYS的前处理中建立棒材的三维几何模型，设置棒材直径Φ＝120mm，长度300mm。或在三维建模软件中建模，将节点单元数据读入ANSYS中，如图2所示。
步骤2：确定钢种并建立材料热物性参数模型
确定钢种化学成分，C：0.17-0.23％，Mn：0.8-1.1％，Si：0.17-0.37％，Gr：1.00-1.3％，P：≤0.012，S：≤0.035，Alt：0.04％，其余为铁和不可避免的杂质。
材料的物理性能参数随温度的变化而变化，需要考虑不同温度下的材料物理性能参数，包括材料不同温度下的材料密度、导热系数、比热容。如下表1，在温度-相变耦合中定义材料的属性：密度、导热系数、比热容等。
表1材料热物性参数定义
步骤3：在有限元分析软件中对棒材的三维模型进行网格划分，采用mapped网格划分方式，网格采用六面体网格，单元尺寸设为1，如图3所示。
步骤4：定义单元类型与实常数
热实体单元的表面上同时施加热流载荷和对流载荷，又或者施加指定的辐射等，在这种情况，可以用表面效应单元覆盖需要施加载荷的表面并使用它们作为一个管道以施加所需的载荷。针对大棒的空冷问题研究，选用SOLID70实体单元和SURF152表面效应单元。SOLID70实体单元是一个具有导热能力的单元，该单元有八个节点，每个节点只有一个自由度，该单元可用于三维的稳态或瞬态热分析问题，如图4所示。SURF152表面效应单元可用于多种载荷和表面效应的应用。它可以覆盖在任何三维热单元的表面，该单元可用于三维热分析，多种载荷和表面效应可以同时存在。关于实常数的定义就是设置表面效应单元的斯蒂芬·波尔兹曼常数为5.67×10-8，定义任意空间节点。
步骤5：设置初始条件和边界条件
导热微分方程描述的是固体内部导热的一般规律，其解有无数个，在解决工程问题时需要加入边界条件。棒材在空冷过程中与环境的热交换主要是辐射散热和对流散热，因此需要确定对流换热条件和辐射换热条件。
同时存在辐射换热和对流换热的表面边界条件为：
式中TW代表固体的表面温度(℃)，TA代表流体主流温度(℃)，hz代表综合换热系数(W/m2·℃)；ha代表对流换热系数(W/m2·℃)，hr代表辐射换热(W/m2·℃)。
步骤6：施加载荷及载荷选项
将载荷施加在有限元模型的节点和单元上。根据实际情况，对棒材三维几何模型进行加载，包括棒材初始均匀温度场、环境温度，设置棒材初始温度为1150℃，环境温度为25℃。定义分析类型为瞬态分析，载荷步选择Stepped。
步骤7：导入相变热模型
在利用ANSYS计算空冷状态下相变潜热时，利用补热法进行模拟分析。将相变潜热作为一个内热源加入到ANSYS分析过程中，也就是在相变区间内，通过内热源释放等效的热量从而替代相变产生的热量。
具体相变热模型如下：
即：
式中T代表温度(℃)，x、y分别为中间坯厚度和宽度方向坐标，λ代表试验钢的导热系数(W/m·k)，ρ代表试验钢的密度(kg/m3)，Cp代表试验钢的比热容(J/kg·K)，qv代表内热源强度(W/m3)
步骤8：进行模拟求解设置
采用牛顿-拉夫森法进行求解，牛顿-拉夫森法速度快，迭代次数少，求解器选用SPARSE求解器。求解的结果为：(1)节点的自由度值，为基本解；(2)原始解的导出值，为单元解。单元解通常是在单元的公共点上计算出的，ANSYS程序将结果写入数据库和结果文件。
步骤9：根据实际情况设置计算时间步长和迭代步数，完成计算。
步骤10：模拟求解结束后，进入时间—历程POST26后处理器，对棒材三维几何模型进行瞬态分析，分析三维几何模型中心位置在任意时间点上的温度变化，获取温度值，计算棒材芯部考虑相变热的温度变化曲线。
为了验证本发明方法结果的可靠性与准确性，申请人开展了物理模型试验，并将结果与数值模拟的结果进行对比，图5和图6分别为物理试验得到的棒材中心温度与不考虑相变热的数值模拟对比和物理试验得到的棒材中心温度与考虑相变热的数值模拟对比，可以看出，未考虑相变的数值模拟的结果与物理模型试验结果有较大的差异，考虑相变热后的数值模拟结果与实验结果吻合度高，可认为本发明方法构建的数值模型能够较为真实地模拟存在相变的大规格棒材空冷过程中的温度变化。
Claims (7)
1.一种考虑相变热的大规格棒材空冷温度场计算方法，其特征在于，包括以下步骤：
步骤1：根据棒材几何参数，在有限元分析软件ANSYS的前处理模块中建立棒材的三维几何模型，或在三维建模软件中建立棒材的三维几何模型，并将节点单元数据读入有限元分析软件ANSYS中；
步骤2：建立材料热物性参数模型；
步骤3：在有限元分析软件ANSYS中对棒材的三维几何模型进行网格划分；
步骤4：定义实体单元和表面效应单元类型及表面效应单元实常数；
步骤5：设置初始条件和边界条件；
步骤6：施加载荷及载荷选项；
步骤7：导入相变热模型；
利用有限元分析软件ANSYS计算空冷状态下相变潜热时，利用补热法进行模拟分析，将相变潜热作为一个内热源加入到有限元分析软件ANSYS分析过程中，即在相变区间内，通过内热源释放等效的热量从而替代相变产生的热量，具体相变热模型如下：
即：
式中T代表温度(℃)，x、y分别为中间坯厚度和宽度方向坐标，λ代表试验钢的导热系数(W/m·k)，ρ代表试验钢的密度(kg/m3)，Cp代表试验钢的比热容(J/kg·K)，qv代表内热源强度(W/m3)；
步骤8：进行模拟求解设置；
步骤9：设置计算时间步长和迭代步数，完成求解；
步骤10：进入时间-历程POST26后处理器，对所述相变热模型进行瞬态分析，获取棒材三维几何模型在冷却过程中任意时间点上的温度值，计算棒材芯部考虑相变热的温度变化曲线。
2.如权利要求1所述的考虑相变热的大规格棒材空冷温度场计算方法，其特征在于，所述材料热物性参数包括材料不同温度下的材料密度、导热系数、比热容。
3.如权利要求1所述的考虑相变热的大规格棒材空冷温度场计算方法，其特征在于，所述步骤4，选用SOLID70实体单元和SURF152表面效应单元覆盖需要施加载荷的热实体表面，设置SURF152表面效应单元的斯蒂芬·波尔兹曼常数，将SURF152表面效应单元作为施加所需载荷的管道。
4.如权利要求1所述的考虑相变热的大规格棒材空冷温度场计算方法，其特征在于，所述步骤5，设置边界条件为：
式中TW代表固体的表面温度(℃)，TA代表流体主流温度(℃)，hz代表综合换热系数(W/m2·℃)；ha代表对流换热系数(W/m2·℃)，hr代表辐射换热(W/m2·℃)。
5.如权利要求1所述的考虑相变热的大规格棒材空冷温度场计算方法，其特征在于，所述步骤6，将载荷施加在棒材的三维几何模型的节点和单元上，根据实际情况对三维几何模型进行加载，包括棒材初始均匀温度场、环境温度，定义分析类型为瞬态分析，载荷步选择Stepped。
6.如权利要求1所述的考虑相变热的大规格棒材空冷温度场计算方法，其特征在于，所述步骤8，采用牛顿-拉夫森法进行求解，求解器选用SPARSE求解器，设置求解的结果为：(1)节点的自由度值，为基本解；(2)原始解的导出值，为单元解。
7.如权利要求1所述的考虑相变热的大规格棒材空冷温度场计算方法，其特征在于，所述步骤3，采用mapped网格划分方式，网格采用六面体网格，单元尺寸设为1。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110888962.5A CN113627059B (zh) | 2021-08-03 | 2021-08-03 | 一种考虑相变热的大规格棒材空冷温度场计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110888962.5A CN113627059B (zh) | 2021-08-03 | 2021-08-03 | 一种考虑相变热的大规格棒材空冷温度场计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113627059A CN113627059A (zh) | 2021-11-09 |
| CN113627059B true CN113627059B (zh) | 2024-05-14 |
Family
ID=78382487
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110888962.5A Active CN113627059B (zh) | 2021-08-03 | 2021-08-03 | 一种考虑相变热的大规格棒材空冷温度场计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113627059B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118036410B (zh) * | 2024-04-11 | 2024-07-09 | 浙江大学 | 一种原位热固化工艺的温度分布模拟方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102008048479A1 (de) * | 2008-03-05 | 2009-09-10 | Livermore Software Technology Corporation, Livermore | Verfahren und Systeme der technischen Analyse unter Verwendung eines hybriden Ansatzes mit FEM und adaptiver SPH |
| CN105586476A (zh) * | 2016-01-19 | 2016-05-18 | 上海交通大学 | 基于数值模拟技术的高强钢淬火工艺优化方法 |
| CN106055387A (zh) * | 2016-06-22 | 2016-10-26 | 华南理工大学 | 基于ansys apdl与ansys cfx的钢芯铝绞线径向温度仿真方法 |
| CN109858085A (zh) * | 2018-12-26 | 2019-06-07 | 钢铁研究总院 | 一种金属材料热处理过程中的奥氏体化测定方法 |
| CN113130019A (zh) * | 2021-03-31 | 2021-07-16 | 西安理工大学 | 一种tc4钛合金焊接熔池固态相变的模拟方法 |
-
2021
- 2021-08-03 CN CN202110888962.5A patent/CN113627059B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102008048479A1 (de) * | 2008-03-05 | 2009-09-10 | Livermore Software Technology Corporation, Livermore | Verfahren und Systeme der technischen Analyse unter Verwendung eines hybriden Ansatzes mit FEM und adaptiver SPH |
| CN105586476A (zh) * | 2016-01-19 | 2016-05-18 | 上海交通大学 | 基于数值模拟技术的高强钢淬火工艺优化方法 |
| CN106055387A (zh) * | 2016-06-22 | 2016-10-26 | 华南理工大学 | 基于ansys apdl与ansys cfx的钢芯铝绞线径向温度仿真方法 |
| CN109858085A (zh) * | 2018-12-26 | 2019-06-07 | 钢铁研究总院 | 一种金属材料热处理过程中的奥氏体化测定方法 |
| CN113130019A (zh) * | 2021-03-31 | 2021-07-16 | 西安理工大学 | 一种tc4钛合金焊接熔池固态相变的模拟方法 |
Non-Patent Citations (2)
| Title |
|---|
| 基于ANSYS有限元方法对相变材料相变过程的分析;郭志强;《新技术新工艺》;正文第1-4页 * |
| 热连轧棒材温度场及相变模拟研究;管晓光;《中国优秀硕士学位论文全文数据库 (工程科技Ⅰ辑)》;正文第1-70页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113627059A (zh) | 2021-11-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3944125B1 (en) | Method and system for manufacturing a heat exchanger for supercritical pressure fluid | |
| US8437991B2 (en) | Systems and methods for predicting heat transfer coefficients during quenching | |
| Cho et al. | A study on the hot-deformation behavior and dynamic recrystallization of Al–5 wt.% Mg alloy | |
| Yang et al. | Finite element analysis of quenching temperature field, residual stress and distortion in A357 aluminum alloy large complicated thin-wall workpieces | |
| Bobbili et al. | Constitutive modeling of hot deformation behavior of high-strength armor steel | |
| Malinowski et al. | Inverse method implementation to heat transfer coefficient determination over the plate cooled by water spray | |
| He et al. | Inverse determination of temperature-dependent thermophysical parameters using multiobjective optimization methods | |
| Lan et al. | Shape design for heat conduction problems using curvilinear grid generation, conjugate gradient, and redistribution methods | |
| CN113627059B (zh) | 一种考虑相变热的大规格棒材空冷温度场计算方法 | |
| CN112906273A (zh) | 一种板材真空退火过程温度场预测方法 | |
| Li et al. | Quenching residual stress distributions in aluminum alloy plates with different dimensions | |
| Zeng et al. | A new constitutive model for thermal deformation of magnesium alloys | |
| Narayan Prabhu et al. | Inverse modeling of heat transfer with application to solidification and quenching | |
| Liu et al. | Analysis and improvement of material flow during extrusion process using spreading pocket die for large-size, flat-wide, and multi-ribs profile | |
| CN108446414B (zh) | 一种由3d打印多孔结构随机缺陷的逆向预测方法 | |
| Reichl et al. | A numerical simulation of thermodynamic processes for cryogenic metal forming of aluminum sheets and comparison with experimental results | |
| Wu et al. | Coupling a cellular automaton model with a finite element model for simulating deformation and recrystallization of a low-carbon micro-alloyed steel during hot compression | |
| Vu et al. | Comparison of pore-level and volume-averaged computations in highly conductive spherical-void-phase porous materials | |
| Papacharalampopoulos et al. | Towards a surrogate spatiotemporal model of additive manufacturing for digital twin-based process control | |
| Jan et al. | CFD investigation of quench media and orientation effects on structural stress induced in the intense quenching processes for aluminum cylinder heads | |
| Malendowski et al. | Heat transfer model for calculation of the temperature field inside thin-walled sections in fire | |
| CN112084723B (zh) | 光纤预制棒一次拉伸工艺仿真方法及装置 | |
| Huang et al. | Multiscale topology optimization of cellular structures with high thermal conductivity and large convective surface area | |
| Ogorodnikova et al. | Reconstruction of thermo-physical properties to improve material database for casting simulation | |
| Cosentino et al. | Numerical modeling of vacuum heat treatment of nickel-based superalloys |
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 |