CN113011068B - Three-dimensional simulation method for walking beam type plate blank heating - Google Patents

Three-dimensional simulation method for walking beam type plate blank heating Download PDF

Info

Publication number
CN113011068B
CN113011068B CN202110319227.2A CN202110319227A CN113011068B CN 113011068 B CN113011068 B CN 113011068B CN 202110319227 A CN202110319227 A CN 202110319227A CN 113011068 B CN113011068 B CN 113011068B
Authority
CN
China
Prior art keywords
dimensional simulation
slab
equation
heating
simulation model
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
Application number
CN202110319227.2A
Other languages
Chinese (zh)
Other versions
CN113011068A (en
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.)
Ganjiang New Area Aobo Particle Technology Research Institute Co ltd
Jiangxi University of Science and Technology
Original Assignee
Ganjiang New Area Aobo Particle Technology Research Institute Co ltd
Jiangxi University of Science and Technology
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 Ganjiang New Area Aobo Particle Technology Research Institute Co ltd, Jiangxi University of Science and Technology filed Critical Ganjiang New Area Aobo Particle Technology Research Institute Co ltd
Priority to CN202110319227.2A priority Critical patent/CN113011068B/en
Publication of CN113011068A publication Critical patent/CN113011068A/en
Application granted granted Critical
Publication of CN113011068B publication Critical patent/CN113011068B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Fluid Mechanics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Control Of Heat Treatment Processes (AREA)

Abstract

The invention provides a walking beam type plate blank heating three-dimensional simulation method, which comprises the following steps: building a three-dimensional simulation model of the heating furnace and the plate blank; generating a grid for the three-dimensional simulation model, and encrypting the combustor and the slab area; acquiring a control equation and a slab motion equation according to the grids; simulating and solving the three-dimensional simulation model based on a control equation and a slab motion equation; and (4) visually solving a result, and obtaining temperature cloud pictures of the plate blank in different positions in the heating furnace. The method can be used for identifying whether the plate blank reaches the tapping temperature standard or not, carrying out fault diagnosis and optimizing the reheating process, and greatly saves the labor cost.

Description

Three-dimensional simulation method for walking beam type plate blank heating
Technical Field
The invention relates to the technical field of computer numerical simulation, in particular to a walking beam type plate blank heating three-dimensional simulation method.
Background
Hot rolling of slabs is a key process for producing high quality steel products, and a heating furnace is one of key apparatuses of a steel rolling line and is also a major energy consumption apparatus. Because the heating furnace is located at the upstream of the steel rolling mill, and the defective slab products can damage downstream equipment, unnecessary losses such as product reworking, energy resource waste, shutdown and equipment maintenance/replacement are caused, and the quality of the products output by the heating furnace is particularly important. The purpose of heating the slab using a furnace is to heat the slab to a target rolling temperature to reduce the stresses introduced during rolling. Industrial furnaces generally consist of three zones: a preheating zone, a heating zone, and a soaking zone. During the moving process of the plate blank in the furnace, the heat released by combustion is transferred to the steel plate by radiation, convection, conduction and the like. The heating quality of the slab and the temperature distribution inside the slab are important indicators for determining the heating degree of the slab. However, in the actual production process, the steel mill rolls different steel grades every day, and the optimal balance between slab reheating quality and furnace output can only be achieved if there is a good temperature prediction for the slab.
The traditional slab temperature prediction method generally detects the real-time temperature of the slab through technicians, and then predicts the slab temperature according to the real-time temperature of the slab, thereby being time-consuming, labor-consuming and high in cost.
Disclosure of Invention
Based on the above, in order to solve the problems of time and labor waste and higher cost of the traditional slab temperature prediction method, the invention provides a walking beam type slab heating three-dimensional simulation method, which has the following specific technical scheme:
a walking beam type plate blank heating three-dimensional simulation method comprises the following steps:
building a three-dimensional simulation model of the heating furnace and the plate blank;
generating a grid for the three-dimensional simulation model, and encrypting a combustor and a slab region;
acquiring a control equation and a slab motion equation according to the grids;
simulating and solving the three-dimensional simulation model based on the control equation and the slab motion equation;
and (4) visually solving a result, and obtaining temperature cloud pictures of the plate blank in different positions in the heating furnace.
According to the walking beam type three-dimensional simulation method for heating the plate blank, the heating furnace and the three-dimensional simulation model of the plate blank are constructed, the three-dimensional simulation model is subjected to grid generation processing, then a control equation and a plate blank motion equation are obtained according to the grid, finally the three-dimensional simulation model is subjected to simulation solution and visual solution results based on the control equation and the plate blank motion equation, temperature cloud pictures of the plate blank in different positions in the heating furnace are obtained, the temperature distribution of the plate blank in the heating furnace can be simulated and simulated, compared with the traditional plate blank temperature prediction method, the labor cost is greatly saved, whether the plate blank reaches the temperature standard of discharging the furnace or not can be recognized, fault diagnosis is carried out, the reheating process is optimized, and a steel mill can realize high-efficiency production while the quality of steel products is guaranteed.
Further, a three-dimensional simulation model of the heating furnace and the slab is constructed according to the space structures of the heating furnace and the slab to be simulated.
Further, the specific method for performing the mesh generation processing on the three-dimensional simulation model includes the following steps:
partitioning the three-dimensional simulation model into grids in a partitioning mode to generate unstructured tetrahedral and hexahedral grids;
and carrying out grid refinement treatment on the combustor and the slab region in the three-dimensional simulation model.
Further, the governing equation includes:
conservation of mass equation
Figure BDA0002992125300000031
Equation of conservation of momentum
Figure BDA0002992125300000032
Equation of energy of turbulence
Figure BDA0002992125300000033
Figure BDA0002992125300000034
Component transfer equation of each component
Figure BDA0002992125300000035
Component transfer equation for combustion average component i
Figure BDA0002992125300000036
Equation of gas radiation characteristics
Figure BDA0002992125300000037
Where p is the static pressure, p
Figure BDA0002992125300000038
And
Figure BDA0002992125300000039
gravity and external force, respectively, ρ is the fluid density,
Figure BDA00029921253000000310
is the fluid velocity with x, y and z directional components,
Figure BDA00029921253000000311
is the tensor of the stress(s),
Figure BDA00029921253000000312
C μ for empirical coefficients, k is the kinetic energy of the turbulence, ε is the dissipation ratio of the turbulence, σ k Is the turbulant number of k, G k Is the turbulent kinetic energy, G, caused by the mean velocity gradient b Is the kinetic energy of the turbulence, Y, caused by buoyancy M To account for the contribution of wave expansion in compressible turbulence to the overall dissipation ratio, ε is the turbulent dissipation ratio, S k Is a preset value, R i Is the net rate of chemical reaction of component i,
Figure BDA00029921253000000313
is the diffusion flux of component i,
Figure BDA00029921253000000314
is that component i is reacted at a reaction time τ * After mass fraction, Y i Is the average mass fraction, ξ, of component i * Is the mass fraction, τ, in the structure * Is the residence time of the reaction, a ε,i Emissivity weighting coefficient, k, of the ith gas i Is the absorption coefficient, P is the sum of the partial pressures of all gases, S is the path length, e is the natural constant.
Further, the slab motion equation is as follows:
Figure BDA0002992125300000041
Figure BDA0002992125300000042
where p is the density of the fluid,
Figure BDA0002992125300000043
is a vector of the flow velocity and,
Figure BDA0002992125300000044
is the moving speed of the moving grid, gamma is the diffusion coefficient, V is the control volume, phi is the general scalar,
Figure BDA00029921253000000411
representing the boundaries of the control volume.
Further, the specific method for simulating and solving the three-dimensional simulation model according to the control equation and the slab motion equation comprises the following steps:
determining a reaction mechanism according to the combustion reaction of the three-dimensional simulation model;
setting initial conditions, and acquiring the temperature of an initial plate blank according to the initial conditions;
calculating the furnace temperature under the problem according to the initial slab temperature and taking the furnace temperature as the initial furnace temperature;
and setting a total time step, and solving the three-dimensional simulation model according to the initial furnace temperature and the total time step to obtain the temperature change of the plate blank in the heating process.
Further, the governing equations further include a momentum equation
Figure BDA0002992125300000045
Figure BDA0002992125300000046
Wherein, delta ij Is the stress tensor due to molecular viscosity.
Further, the governing equation also comprises an energy conservation equation
Figure BDA0002992125300000047
Figure BDA0002992125300000048
Wherein k is eff Is the effective electrical conductivity of the conductive material,
Figure BDA0002992125300000049
is the diffusion flux of component j, S h Including the heat of chemical reaction and any other heat of volume, E from
Figure BDA00029921253000000410
By definition, where h is the sensible enthalpy of the incompressible stream.
Furthermore, the furnace wall of the heating furnace is an insulating wall surface, the boundary condition of the plate blank is a coupling wall, and the boundary condition of the slide block in the heating furnace is constant heat flux.
Accordingly, the present invention provides a computer-readable storage medium storing a computer program which, when executed, implements the walking beam type slab heating three-dimensional simulation method as described above.
Drawings
The invention will be further understood from the following description in conjunction with the accompanying drawings. The components in the figures are not necessarily to scale, emphasis instead being placed upon illustrating the principles of the embodiments. Like reference numerals designate corresponding parts throughout the different views.
Fig. 1 is a schematic overall flow chart of a walking beam type slab heating three-dimensional simulation method according to an embodiment of the present invention;
fig. 2 is a temperature cloud diagram of a slab at different positions according to a walking beam type slab heating three-dimensional simulation method in an embodiment of the invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to embodiments thereof. It should be understood that the detailed description and specific examples, while indicating the scope of the invention, are intended for purposes of illustration only and are not intended to limit the scope of the invention.
It will be understood that when an element is referred to as being "secured to" another element, it can be directly on the other element or intervening elements may also be present. When an element is referred to as being "connected" to another element, it can be directly connected to the other element or intervening elements may also be present. The terms "vertical," "horizontal," "left," "right," and the like as used herein are for illustrative purposes only and do not represent the only embodiments.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. The terminology used herein in the description of the invention is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used herein, the term "and/or" includes any and all combinations of one or more of the associated listed items.
The terms "first" and "second" used herein do not denote any particular order or quantity, but rather are used to distinguish one element from another.
As shown in fig. 1, a walking beam type three-dimensional simulation method for slab heating in an embodiment of the present invention includes the following steps:
building a three-dimensional simulation model of the heating furnace and the plate blank;
generating a grid for the three-dimensional simulation model, and encrypting a combustor and a slab region;
acquiring a control equation and a slab motion equation according to the grids;
simulating and solving the three-dimensional simulation model based on the control equation and the slab motion equation;
and (4) visually solving a result, and obtaining temperature cloud pictures of the plate blank in different positions in the heating furnace.
The visual solution result comprises the slab speed on any plane in the heating furnace and temperature clouds of the slabs at different positions, and is shown in fig. 2.
The walking beam type three-dimensional simulation method for heating the plate blank comprises the steps of firstly constructing a heating furnace and a three-dimensional simulation model of the plate blank, generating a grid for the three-dimensional simulation model, then obtaining a control equation and a plate blank motion equation according to the grid, finally carrying out simulation solution and visualization solution on the three-dimensional simulation model based on the control equation and the plate blank motion equation, obtaining temperature cloud charts of the plate blank in different positions in the heating furnace, simulating the temperature distribution of the plate blank in the heating furnace, and accurately predicting the temperature and the temperature distribution of the plate blank in the heating furnace. Compared with the traditional slab temperature prediction method, the walking beam type slab heating three-dimensional simulation method greatly saves labor cost, can be used for identifying whether the slab reaches the tapping temperature standard and carrying out fault diagnosis and optimizing the reheating process, and enables a steel mill to realize high-efficiency production while ensuring the quality of steel products.
In one embodiment, an actual walking beam type slab and a heating furnace drawing are obtained, and then a three-dimensional simulation model of the heating furnace and the slab is constructed according to the space structure of the heating furnace and the slab to be simulated. The heating furnace is installed with a burner, which is divided into a preheating zone, a heating zone, and a soaking zone.
In one embodiment, the specific method for performing the mesh generation processing on the three-dimensional simulation model includes the following steps:
partitioning the three-dimensional simulation model into grids in a partitioning mode to generate unstructured tetrahedral and hexahedral grids;
and carrying out grid refinement treatment on the combustor and the slab region in the three-dimensional simulation model.
Specifically, tetrahedral and hexahedral partition grid division is performed on the three-dimensional simulation model by using an ANSYS Meshing advanced grid division technology, and a three-dimensional simulation grid model with refined grids, namely unstructured tetrahedral and hexahedral grids, is generated. Therefore, the calculation accuracy of the walking beam type plate blank heating three-dimensional simulation method can be improved.
In one embodiment, the governing equation comprises:
conservation of mass equation
Figure BDA0002992125300000071
Equation of conservation of momentum
Figure BDA0002992125300000072
Equation of turbulent energy
Figure BDA0002992125300000073
Figure BDA0002992125300000074
Component transfer equation of each component
Figure BDA0002992125300000075
Component transfer equation for combustion average component i
Figure BDA0002992125300000081
Equation of gas radiation characteristics
Figure BDA0002992125300000082
Where p is the static pressure, p
Figure BDA0002992125300000083
And
Figure BDA0002992125300000084
gravity and external force, respectively, ρ is the fluid density,
Figure BDA0002992125300000085
is the fluid velocity with x, y and z directional components,
Figure BDA0002992125300000086
is the tensor of the stress(s),
Figure BDA0002992125300000087
C μ for empirical coefficients, k is the kinetic energy of the turbulence, ε is the dissipation ratio of the turbulence, σ k Is the turbulant number of k, G k Is the turbulent kinetic energy, G, caused by the mean velocity gradient b Is the kinetic energy of the turbulence, Y, caused by buoyancy M To account for the contribution of wave expansion in compressible turbulence to the overall dissipation ratio, ε is the turbulent dissipation ratio, S k Is a preset value, R i Is a chemical reaction of component iThe net rate of the flow is,
Figure BDA0002992125300000088
is the diffusion flux of component i,
Figure BDA0002992125300000089
is component i at reaction time τ * After mass fraction, Y i Is the average mass fraction, ξ, of component i * Is the mass fraction, τ, in the structure * Is the residence time of the reaction, a ε,i Emissivity weighting coefficient, k, of the ith gas i Is the absorption coefficient, P is the sum of the partial pressures of all gases, S is the path length, e is the natural constant, μ is the kinetic viscosity, u is the fluid mean velocity i Represents the average velocity in the direction of flow I, I is the total number of imaginary gray gases, a epsilon, I (t) is the emissivity weighting factor for the ith imaginary gray gas.
In one embodiment, the moving speed of a specific slab and the three-dimensional simulation grid model are controlled according to a User Definition Function (UDF) code developed under actual conditions, so as to realize the dynamic motion of the slab in the three-dimensional simulation model. The slab motion equation is as follows:
Figure BDA00029921253000000810
where p is the density of the fluid,
Figure BDA00029921253000000811
is a vector of the flow velocity of the fluid,
Figure BDA00029921253000000812
is the moving speed of the moving grid, gamma is the diffusion coefficient, V is the control volume, phi is the general scalar,
Figure BDA00029921253000000813
representing the boundary of the control volume, a being the total surface area of the fluid boundary,
Figure BDA00029921253000000814
is a vector of a.
In one embodiment, the moving speed of the slab is the average moving speed per minute and does not include the moving speed of the slab in the vertical direction, so that the efficiency of calculating the moving speed of the slab by the method is improved.
In one embodiment, the specific method for simulating and solving the three-dimensional simulation model according to the control equation and the slab motion equation comprises the following steps:
firstly, determining a reaction mechanism according to the combustion reaction of the three-dimensional simulation model.
To reduce the complexity of the boundary condition changes, the fuel input is averaged hourly. A two-step reaction mechanism is employed for methane-air to obtain more detailed component concentrations. A one-step reaction mechanism is employed for ethane-air due to the smaller ethane-air combustion ratio.
And secondly, setting initial conditions, and acquiring the temperature of the initial slab according to the initial conditions.
Specifically, the specific slab length, width, initial temperature, steel grade and position of all slabs in the furnace are arranged according to actual production conditions, based on the conditions at a certain point in time of a given rolling mill test as initial conditions. And setting thermal physical parameters (including thermal conductivity, specific heat capacity and the like) of the slab and the slide block according to the material property of the slab.
And thirdly, calculating the furnace temperature under the problem according to the initial slab temperature and taking the furnace temperature as the initial furnace temperature.
And fourthly, setting a total time step length, and solving the three-dimensional simulation model according to the initial furnace temperature and the total time step length to obtain the temperature change of the slab in the heating process.
Specifically, according to the thermophysical parameters of the plate blank, the three-dimensional simulation model is numerically solved by using ANSYS Fluent software.
In one embodiment, the three-dimensional simulation Model comprises a turbulence Model, a combustion Model and a radiation Model, wherein the turbulence Model is a kappa-epsilon Model, the radiation Model is a WSGGM (Weighted Sum of Gray Gases Model), and the combustion Model is an EDC turbulence combustion Model to consider the combustion reaction of natural gas.
In one embodiment, the governing equations further include a momentum equation
Figure BDA0002992125300000101
Figure BDA0002992125300000102
Wherein, delta ij Is the stress tensor caused by molecular viscosity, u' is the pulsation velocity, u i 、u j 、u l Respectively representing the average velocity in the direction of the fluid i, j, l,
Figure BDA0002992125300000103
in one embodiment, the governing equations further comprise an energy conservation equation
Figure BDA0002992125300000104
Figure BDA0002992125300000105
Wherein k is eff Is the effective electrical conductivity of the conductive material,
Figure BDA0002992125300000106
is the diffusion flux of component j, S h Including the heat of chemical reaction and any other heat of volume, E from
Figure BDA0002992125300000107
Where h is the sensible enthalpy of the incompressible flow,
Figure BDA0002992125300000108
which is indicative of a temperature gradient, is,
Figure BDA0002992125300000109
j is the constant pressure specific heat of component j,
Figure BDA00029921253000001010
v is the slab velocity, the effective stress tensor.
In one embodiment, the furnace wall of the heating furnace is an insulating wall surface, the boundary condition of the plate blank is a coupling wall, the plate blank transfers heat in a convection mode and a radiation mode, the boundary condition of a sliding block in the heating furnace is constant heat flux, and the feeding door and the discharging door are both in a closed state.
In one embodiment, the invention provides a computer-readable storage medium storing a computer program which, when executed, implements the walking beam type slab heating three-dimensional simulation method as described above.
The technical features of the embodiments described above may be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the embodiments described above are not described, but should be considered as being within the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present invention, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that various changes and modifications can be made by those skilled in the art without departing from the spirit of the invention, and these changes and modifications are all within the scope of the invention. Therefore, the protection scope of the present patent shall be subject to the appended claims.

Claims (9)

1. A walking beam type plate blank heating three-dimensional simulation method is characterized by comprising the following steps:
building a three-dimensional simulation model of the heating furnace and the plate blank;
generating a grid for the three-dimensional simulation model, and encrypting a combustor and a slab region;
acquiring a control equation and a slab motion equation according to the grids;
simulating and solving the three-dimensional simulation model based on the control equation and the slab motion equation;
visualizing a solution result to obtain temperature cloud charts of the plate blank in different positions in the heating furnace;
wherein the governing equation comprises:
conservation of mass equation
Figure FDA0003786842120000011
Equation of conservation of momentum
Figure FDA0003786842120000012
Equation of energy of turbulence
Figure FDA0003786842120000013
Figure FDA0003786842120000014
Component transfer equation of each component
Figure FDA0003786842120000015
Component transfer equation for combustion average component i
Figure FDA0003786842120000016
Equation of gas radiation characteristics
Figure FDA0003786842120000017
Wherein, p is a static pressure,
Figure FDA0003786842120000018
and
Figure FDA0003786842120000019
respectively gravity and outerThe force, p, is the density of the fluid,
Figure FDA00037868421200000110
is the fluid velocity with x, y and z directional components,
Figure FDA00037868421200000111
is the tensor of the stress(s),
Figure FDA00037868421200000112
C μ for empirical coefficients, k is the kinetic energy of the turbulence, ε is the dissipation ratio of the turbulence, σ k Is the turbulant number of k, G k Is the turbulent kinetic energy, G, caused by the mean velocity gradient b Is the kinetic energy of the turbulence, Y, caused by buoyancy M Contribution of wave expansion in compressible turbulence to the overall dissipation ratio, S k Is a preset value, R i Is the net rate of chemical reaction of component i,
Figure FDA0003786842120000021
is the diffusion flux of component i,
Figure FDA0003786842120000022
is that component i is reacted at a reaction time τ * After mass fraction, Y i Is the average mass fraction, ξ, of component i * Is the mass fraction, τ, in the structure * Is the residence time of the reaction, a ε,i Emissivity weighting coefficient, k, of the ith gas i Is the absorption coefficient, P is the sum of the partial pressures of all gases, S is the path length, e is the natural constant.
2. The method as claimed in claim 1, wherein the three-dimensional simulation model of the heating furnace and the slab is constructed according to the space structure of the heating furnace and the slab to be simulated.
3. The walking beam type slab heating three-dimensional simulation method as claimed in claim 1, wherein the specific method for generating the grid of the three-dimensional simulation model comprises the following steps:
partitioning the three-dimensional simulation model into grids in a partitioning mode to generate unstructured tetrahedral and hexahedral grids;
and carrying out grid refinement treatment on the combustor and the slab region in the three-dimensional simulation model.
4. The walking beam type slab heating three-dimensional simulation method as claimed in claim 1, wherein the slab motion equation is as follows:
Figure FDA0003786842120000023
Figure FDA0003786842120000024
where p is the density of the fluid,
Figure FDA0003786842120000025
is a vector of the flow velocity of the fluid,
Figure FDA0003786842120000026
is the moving speed of the moving grid, gamma is the diffusion coefficient, V is the control volume, phi is the general scalar,
Figure FDA0003786842120000027
representing the boundaries of the control volume.
5. The walking beam type slab heating three-dimensional simulation method as claimed in claim 1, wherein the specific method for simulating and solving the three-dimensional simulation model according to the control equation and the slab motion equation comprises the following steps:
determining a reaction mechanism according to the combustion reaction of the three-dimensional simulation model;
setting initial conditions, and acquiring the temperature of an initial plate blank according to the initial conditions;
calculating the furnace temperature under the problem according to the initial slab temperature, and taking the furnace temperature as the initial furnace temperature;
and setting a total time step, and solving the three-dimensional simulation model according to the initial furnace temperature and the total time step to obtain the temperature change of the plate blank in the heating process.
6. The walking beam type slab heating three-dimensional simulation method of claim 1, wherein the governing equation further comprises a momentum equation
Figure FDA0003786842120000031
Figure FDA0003786842120000032
Wherein, delta ij Is the stress tensor caused by molecular viscosity.
7. The walking beam type slab heating three-dimensional simulation method of claim 4, wherein the governing equation further comprises an energy conservation equation
Figure FDA0003786842120000033
Figure FDA0003786842120000034
Wherein k is eff Is the effective electrical conductivity of the conductive material,
Figure FDA0003786842120000035
is the diffusion flux of component j, S h Including the heat of chemical reaction and any other heat of volume, E from
Figure FDA0003786842120000036
By definition, where h is the sensible enthalpy of the incompressible stream.
8. The walking beam type slab heating three-dimensional simulation method according to claim 1, wherein the wall of the heating furnace is a heat insulation wall surface, the boundary condition of the slab is a coupling wall, and the boundary condition of the slide block in the heating furnace is a constant heat flux.
9. A computer-readable storage medium, characterized in that it stores a computer program which, when executed, implements a walking beam type slab heating three-dimensional simulation method as claimed in any one of claims 1 to 8.
CN202110319227.2A 2021-03-25 2021-03-25 Three-dimensional simulation method for walking beam type plate blank heating Active CN113011068B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110319227.2A CN113011068B (en) 2021-03-25 2021-03-25 Three-dimensional simulation method for walking beam type plate blank heating

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110319227.2A CN113011068B (en) 2021-03-25 2021-03-25 Three-dimensional simulation method for walking beam type plate blank heating

Publications (2)

Publication Number Publication Date
CN113011068A CN113011068A (en) 2021-06-22
CN113011068B true CN113011068B (en) 2022-09-16

Family

ID=76407014

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110319227.2A Active CN113011068B (en) 2021-03-25 2021-03-25 Three-dimensional simulation method for walking beam type plate blank heating

Country Status (1)

Country Link
CN (1) CN113011068B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113935149A (en) * 2021-09-09 2022-01-14 国网天津市电力公司电力科学研究院 High-voltage cable buffer layer repairing effect simulation method based on graphite powder injection
CN114048696B (en) * 2021-11-18 2024-06-07 长春工业大学 Pneumatic simulation analysis method for flapping-rowing-torsion flapping wing aircraft

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009057190A1 (en) * 2007-10-29 2009-05-07 Japan Agency For Marine-Earth Science And Technology Meteorological phenomena simulation device and method
CN107346353A (en) * 2017-06-05 2017-11-14 民政部零研究所 A kind of solid burning article combustion process emulation mode and server
CN108363857A (en) * 2018-02-05 2018-08-03 燕山大学 Recuperative heater flow field and inside workpiece temperature, thermal-stress analysis method
CN111737891A (en) * 2020-05-11 2020-10-02 自然资源部第二海洋研究所 Three-dimensional simulation method for submarine hydrothermal fluid circulation

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040210429A1 (en) * 2003-02-05 2004-10-21 Moldflow Ireland Ltd. Apparatus and methods for performing process simulation using a hybrid model
US10956625B2 (en) * 2015-09-01 2021-03-23 Siemens Industry Software Inc. Mesh generation system and method
CN107885955B (en) * 2017-11-29 2021-03-09 中昇创举(天津)科技有限公司 Simulation method and device for motor temperature field
CN210984686U (en) * 2020-01-20 2020-07-10 江西理工大学 Rectangular frame-double-rod slow wave structure
CN111859697A (en) * 2020-07-28 2020-10-30 西安建筑科技大学 Simulation method for heating simulation of plate blank of radiant tube heating furnace
CN111961776B (en) * 2020-07-31 2022-04-01 中南大学 Thermocouple position mapping method for corner area of blast furnace hearth lining

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009057190A1 (en) * 2007-10-29 2009-05-07 Japan Agency For Marine-Earth Science And Technology Meteorological phenomena simulation device and method
CN107346353A (en) * 2017-06-05 2017-11-14 民政部零研究所 A kind of solid burning article combustion process emulation mode and server
CN108363857A (en) * 2018-02-05 2018-08-03 燕山大学 Recuperative heater flow field and inside workpiece temperature, thermal-stress analysis method
CN111737891A (en) * 2020-05-11 2020-10-02 自然资源部第二海洋研究所 Three-dimensional simulation method for submarine hydrothermal fluid circulation

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
板坯加热炉三维温度场数学模型系统的开发及应用;商桂梅等;《冶金自动化》;20090331;第33卷(第2期);38-42 *

Also Published As

Publication number Publication date
CN113011068A (en) 2021-06-22

Similar Documents

Publication Publication Date Title
CN113011068B (en) Three-dimensional simulation method for walking beam type plate blank heating
Shaikh et al. Effect of multiple phase change materials (PCMs) slab configurations on thermal energy storage
Hamman et al. A CFD simulation process for fast reactor fuel assemblies
CN109992900A (en) A kind of more real-time online cooperative intelligent emulation modes of mass concrete and system
Khatir et al. Multi-objective Computational Fluid Dynamics (CFD) design optimisation in commercial bread-baking
Bai et al. Thermal stratification in a cylindrical tank due to heat losses while in standby mode
António et al. Air temperature fields inside refrigeration cabins: A comparison of results from CFD and ANN modelling
Zheng et al. Numerical study on impact of non-heating surface temperature on the heat output of radiant floor heating system
CN106906351A (en) A kind of board briquette forecasting model and optimum furnace method
CN105087882B (en) Partitioning method for heat treatment stages of vertical quenching furnace
Kumlutaş et al. Investigation of design parameters of a domestic refrigerator by artificial neural networks and numerical simulations
Zhang et al. A physical model and data-driven hybrid prediction method towards quality assurance for composite components
Karakoyun et al. A comprehensive approach to analyze the discrepancies in heat transfer characteristics pertaining to radiant ceiling heating system
Lee et al. Optimum residence time for steel productivity and energy saving in a hot rolled reheating furnace
Shen et al. Soft sensor modeling of blast furnace wall temperature based on temporal–spatial dimensional finite-element extrapolation
CN114417664B (en) Online simulation and visualization method for steel hot rolling microstructure evolution based on cellular automaton
Oumarou et al. Transient process model of open anode baking furnace
Arkhazloo et al. Determination of temperature distribution during heat treatment of forgings: Simulation and experiment
Arreola-Villa et al. Numerical simulation of a metallic load in continuous movement during heating in an experimental vertical-type furnace
Ben-Nakhi et al. Inter-model comparison of CFD and neural network analysis of natural convection heat transfer in a partitioned enclosure
Hu et al. Zone modelling coupled with dynamic flow pattern for the prediction of transient performance of metal reheating
Arkhazloo et al. A numerical thermal analysis of the heating process of large size forged ingots
Dou et al. Numerical model and optimization strategy for the annealing process of 3D coil cores
Zhang et al. The study of heat-mass transfer characteristics and optimization on electric arc furnace with the novel fence electrode
Oguma et al. Numerical Simulations of a Ground Source Heat Pump System With Pile Heat Exchangers

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