CN118052086A - Heating furnace billet three-dimensional heat transfer temperature field prediction model and construction method thereof - Google Patents
Heating furnace billet three-dimensional heat transfer temperature field prediction model and construction method thereof Download PDFInfo
- Publication number
- CN118052086A CN118052086A CN202410455353.4A CN202410455353A CN118052086A CN 118052086 A CN118052086 A CN 118052086A CN 202410455353 A CN202410455353 A CN 202410455353A CN 118052086 A CN118052086 A CN 118052086A
- Authority
- CN
- China
- Prior art keywords
- furnace
- billet
- dimensional
- temperature
- radiation
- 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
- 238000010438 heat treatment Methods 0.000 title claims abstract description 71
- 238000012546 transfer Methods 0.000 title claims abstract description 50
- 238000010276 construction Methods 0.000 title abstract description 7
- 229910000831 Steel Inorganic materials 0.000 claims abstract description 75
- 239000010959 steel Substances 0.000 claims abstract description 75
- 230000005855 radiation Effects 0.000 claims abstract description 66
- 238000000034 method Methods 0.000 claims description 38
- 239000007789 gas Substances 0.000 claims description 26
- 230000001052 transient effect Effects 0.000 claims description 16
- 238000009826 distribution Methods 0.000 claims description 10
- 238000010079 rubber tapping Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 7
- 239000000463 material Substances 0.000 claims description 7
- 230000004907 flux Effects 0.000 claims description 6
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 3
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 3
- 238000010521 absorption reaction Methods 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 239000003546 flue gas Substances 0.000 claims description 3
- 239000002436 steel type Substances 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 2
- 238000010923 batch production Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000002131 composite material Substances 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
- 238000005261 decarburization Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Landscapes
- Control Of Heat Treatment Processes (AREA)
- Heat Treatments In General, Especially Conveying And Cooling (AREA)
Abstract
The invention discloses a heating furnace billet three-dimensional heat transfer temperature field prediction model and a construction method thereof, belonging to the field of heating furnace temperature field prediction, and comprising a full furnace parameter calling model, a furnace section three-dimensional heat transfer model, a billet three-dimensional heat transfer model and a parameter reading model; the whole furnace parameter calling model is used for calling the whole furnace parameter information; the furnace section three-dimensional heat transfer model is used for calculating the radiation heat exchange heat flow density and the convection heat exchange coefficient of each furnace section according to the whole furnace parameter information; the steel billet three-dimensional heat conduction model is used for updating the three-dimensional temperature field of the steel billet, and then transmitting the surface temperature of the steel billet to the furnace section three-dimensional heat conduction model to calculate the updated radiation heat exchange heat flow density and convection heat exchange coefficient. The heating furnace billet three-dimensional heat transfer temperature field prediction model and the construction method thereof can solve a series of limitations of field experiments, and have the advantages of low cost, high speed, high precision and the like.
Description
Technical Field
The invention relates to the technical field of heating furnace temperature field prediction, in particular to a heating furnace billet three-dimensional heat transfer temperature field prediction model and a construction method thereof.
Background
The heating furnaces such as the walking beam heat treatment furnace, the walking bottom, the roller bottom and the like have the advantages of good billet heating quality, high heating efficiency, low product oxidation decarburization rate and the like, and are widely applied to the metallurgical industry. Meanwhile, with the continuous development of the metallurgical automation technology, the requirement on the heating quality of the steel billet is further improved.
In the running process of the heating furnace, if the influence rule of various operation parameters and structural parameters on the heating quality of the steel billet can be obtained, proper parameter setting can be further found, so that the purposes of improving the heating quality of the steel billet, reducing the production cost and the like are achieved. However, the working mode of the walking beam heating furnace determines that the temperature level of a steel billet in the furnace is difficult to detect in real time in the working process, a large-batch production experiment is required to be carried out when proper parameter setting is required to be found, the experiment cost is high, the experiment period is long, and meanwhile, the optimal solution is difficult to find.
Disclosure of Invention
In order to solve the problems, the invention provides a heating furnace billet three-dimensional heat transfer temperature field prediction model and a construction method thereof, which adopt three levels of basic models of a whole furnace, a furnace section, a billet and the like to construct an efficient calculation frame, realize the heating furnace whole furnace three-dimensional number model, meet the requirement of real-time prediction of the temperature field of all billets in the furnace, facilitate the optimization of auxiliary heating process and guide the heat design of the heating furnace.
In order to achieve the purpose, the invention provides a heating furnace billet three-dimensional heat transfer temperature field prediction model, which comprises a full furnace parameter calling model, a furnace section three-dimensional heat transfer model, a billet three-dimensional heat transfer model and a parameter reading model;
The full furnace parameter calling model is used for calling full furnace parameter information, wherein the full furnace parameter comprises furnace type structure information, furnace inner beam position information, furnace temperature information, steel loading information and steel tapping information;
The furnace section three-dimensional heat transfer model is used for calculating the radiation heat exchange heat flow density among the billet, the furnace gas and the furnace lining in each furnace section according to the whole furnace parameter information; calculating the convective heat transfer coefficient of the surface of the steel billet in each furnace section according to an empirical formula;
The steel billet three-dimensional heat conduction model is used for discretizing a three-dimensional transient heat conduction differential equation of each steel billet on the basis of the previous iteration time step by taking the radiation heat exchange heat flow density and the convection heat flow coefficient as steel billet boundary conditions, calculating three-dimensional transient heat conduction, updating the three-dimensional temperature field of the steel billet, and then transmitting the surface temperature of the steel billet to the furnace section three-dimensional heat transfer model to calculate the updated radiation heat flow density and convection heat flow coefficient.
Preferably, the invention further comprises a parameter reading model for reading the whole furnace parameter information contained in the parameter file, carrying out data interpolation according to the current calculation time and each billet position, and transmitting the interpolated whole furnace parameter information to the whole furnace parameter scheduling model in real time.
Preferably, when used in a walking beam furnace, the full furnace parameter information further includes walking beam step information for updating the location of each billet within the furnace and re-associating the emissivity coefficient file for each furnace segment.
Preferably, in the full furnace parameter call model, the furnace type structure information includes: after dividing the furnace sections into a set number, the furnace section size of each furnace section and the number, size and position information of the cross beams and the longitudinal beams;
The furnace temperature information is furnace temperature parameters assigned to each furnace section by using a furnace temperature distribution model;
The steel loading information is a specific position of a billet in the heating furnace determined by steel type, specification and placement position information of the steel loading, and a radiation angle coefficient file associated with the specific position;
the tapping information is stored tapping billet information.
The method for constructing the prediction model of the three-dimensional heat transfer temperature field of the heating furnace billet comprises the following steps:
S1, reading parameter information of a whole furnace by using a parameter reading model, dividing a heating furnace into a plurality of furnace sections by using a whole furnace parameter calling model based on the parameter information of the whole furnace, and setting each furnace section to comprise a plurality of billets;
S2, calculating angle coefficients between surfaces in the heating furnace by using a Monte Carlo method, and generating an angle coefficient file;
S3, solving the radiation heat exchange heat flux density based on the generated angle coefficient file, and estimating the convection heat flux coefficient according to the flow of the high-temperature flue gas in the furnace by using an empirical formula;
S4, taking the radiation heat exchange heat flow density and the convection heat exchange coefficient as steel billet boundary conditions, discretizing a three-dimensional transient heat conduction differential equation of each steel billet, calculating three-dimensional transient heat conduction, and updating a three-dimensional temperature field of the steel billet;
And S5, feeding the calculated three-dimensional temperature field of the steel billet back to the step S3, and updating the radiation heat exchange heat flow density and the convection heat flow coefficient until the iteration condition is met, and outputting the three-dimensional temperature field of the steel billet according to the approval time interval.
Preferably, the step S2 specifically includes the following steps:
Firstly, when solving the radiation heat exchange of the radiation participation gas in the closed space in the heating furnace, the following assumption is made:
Assuming that the enclosed space in the heating furnace is filled with radiation participation gas, and the spatial distribution of the concentration, the temperature and the pressure is uniform, the radiation participation gas is ash body irrespective of the scattering and reflection characteristics of the gas, and the sum of the absorptivity Ag and the transmissivity Dg of the thermal radiation is 1:
(1)
Meanwhile, the N surfaces forming the closed space in the heating furnace are all ash bodies, and the temperature distribution of each surface is uniform;
Then adopting a radiation thermal resistance network method, wherein the temperature of M surfaces in N surfaces in a closed space in a heating furnace is known and is recorded as Ti, i=1-M; the heat flows of the remaining (N-M) surfaces are known and denoted Qi, i=m+1 to N, the form of the radiation heat exchange equation in the furnace is as follows:
(2)
wherein: representing the type of boundary condition,/> When/>When=1, it means that the k surface is a wall temperature boundary condition, when/>When=0, it indicates that the k surface is a constant heat flow boundary condition; δki is a binary function,/>；/>Is the area of surface k; /(I)Emissivity for surface k; /(I)Absorption rate for radiation-participating gases; /(I)Emissivity for radiation-participating gas; /(I)A temperature of the gas involved for radiation; /(I)Is the Stefan Boltzmann constant; /(I)Is the temperature of surface k; /(I)An angular coefficient of k to i;
Finally, solving the formula (2) by using an LU decomposition method to obtain the net radiation heat flow of each surface ：
(3)
Wherein: Is effective radiation for each surface in the heating furnace.
Preferably, the step S4 specifically includes the following steps:
Firstly, setting the following three-dimensional transient heat conduction differential equation:
(4)
Wherein: t is the temperature of the steel billet, and the unit is K; lambda is the heat conductivity of steel billet, and the unit is ；/>The density of the billet material is Kg/m3; /(I)Specific heat capacity of billet material is expressed as/>; T is heating time, and the unit is s;
Setting the initial temperature of the billet heat transfer model to be uniform and the environment temperature, setting the initial temperature in a heating furnace, and setting the heat exchange condition of the billet surface to be a third type boundary condition, wherein the billet boundary condition is as follows:
(5)
wherein: the unit is W/(m2.K) which is the convection heat transfer coefficient; /(I) The radiation heat exchange heat flow density is W/m < 2 >; the furnace temperature is K; /(I) The unit is K, which is the surface temperature of the billet;
and finally discretizing the formula (5) by adopting an implicit differential method, and solving by adopting a Gausserdel iteration method.
Preferably, in step S5, the three-dimensional temperature field of the billet includes a cloud image of the temperature field of the billet and a temperature variation graph.
The invention has the following beneficial effects:
(1) The three-dimensional heat transfer model of the furnace section is utilized to carry out split decoupling on each furnace section of the whole furnace, each furnace section adopts the same mathematical logic to carry out solving, thus achieving the aim of improving the simulation precision, and each furnace section synchronously solves the three-dimensional transient heat conduction differential equation of each steel billet in the furnace section, updates the steel billet temperature field, achieves the aim of improving the simulation precision, and the size of the heating furnace, the number of the furnace sections, the number of the steel billets in each furnace section, the steel loading information, the steel tapping information, the steel material position information, the temperature information of each furnace section and the like can be set arbitrarily;
(2) The generation condition of the black marks of the steel billets in the walking beam heating furnace can be analyzed by considering the shielding effect of the beam in the furnace on the steel billets in radiation heat exchange and the heat conduction effect of the contact part of the beam in the furnace and the steel billets.
(3) The model has high accuracy, adopts convection heat exchange, radiation heat exchange and heat conduction separation solution, wherein the radiation heat exchange is calculated for each furnace section, and the decoupling of the heat transfer process between the furnace sections is realized, so that the prediction of the temperature of the steel billet is more accurate, and the hit rate of the steel billet temperature prediction accuracy within 5% is higher than 80% under typical working conditions; hit rate of prediction accuracy within 7.5% is higher than 90%; meanwhile, the three-dimensional radiation heat exchange in the furnace is analyzed, so that the adjustment of model parameters is not needed, and the number of the model parameter adjustment is greatly reduced.
The technical scheme of the invention is further described in detail through the drawings and the embodiments.
Drawings
FIG. 1 is a flow chart of a method for constructing a predictive model of a three-dimensional heat transfer temperature field of a heating furnace billet;
FIG. 2 is a schematic diagram showing the structure of a walking beam heating furnace and the method for dividing the geometric structure layer in the furnace section according to the experimental example of the invention
FIG. 3 is a three-dimensional temperature field diagram of a steel billet according to an experimental example of the present invention;
FIG. 4 is a black drawing of the bottom of a billet according to the experimental example of the present invention;
FIG. 5 is a graph showing the temperature of a typical position of a billet according to the experimental example of the present invention with time;
FIG. 6 is a three-dimensional temperature distribution diagram of all billets in a walking beam furnace according to an experimental example of the present invention;
fig. 7 is a graph showing the statistical result of the relative error of the temperature prediction of the billet in the walking beam furnace according to the experimental example of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present application more apparent, the embodiments of the present application will be further described in detail below with reference to the accompanying drawings and examples. It should be understood that the detailed description and specific examples, while indicating the embodiment of the application, are intended for purposes of illustration only and are not intended to limit the scope of the application. All other embodiments, which can be made by those skilled in the art based on the embodiments of the application without making any inventive effort, are intended to be within the scope of the application. Examples of the embodiments are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to like or similar elements throughout or elements having like or similar functionality.
It should be noted that the terms "comprises" and "comprising," and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or server that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed or inherent to such process, method, article, or apparatus, but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
Like reference numerals and letters denote like items in the following figures, and thus once an item is defined in one figure, no further definition or explanation thereof is necessary in the following figures.
The heating furnace billet three-dimensional heat transfer temperature field prediction model comprises a full furnace parameter calling model, a furnace section three-dimensional heat transfer model, a billet three-dimensional heat transfer model and a parameter reading model;
The full furnace parameter calling model is used for calling full furnace parameter information, wherein the full furnace parameter comprises furnace type structure information, furnace inner beam position information, furnace temperature information, steel loading information and steel tapping information;
When used in a walking beam heating furnace, the full furnace parameter information also comprises walking beam stepping information, is used for updating the positions of various billets in the heating furnace and re-associating the radiation angle coefficient files of various furnace sections.
In the full furnace parameter call model, furnace type structure information includes: after dividing the furnace sections into a set number, the furnace section size of each furnace section and the number, size and position information of the cross beams and the longitudinal beams;
The furnace temperature information is furnace temperature parameters assigned to each furnace section by using a furnace temperature distribution model;
The steel loading information is a specific position of a billet in the heating furnace determined by steel type, specification and placement position information of the steel loading, and a radiation angle coefficient file associated with the specific position;
the tapping information is stored tapping billet information.
The furnace section three-dimensional heat transfer model is used for calculating the radiation heat exchange heat flow density among the billet, the furnace gas and the furnace lining in each furnace section according to the whole furnace parameter information; calculating the convective heat transfer coefficient of the surface of the steel billet in each furnace section according to an empirical formula;
The steel billet three-dimensional heat conduction model is used for discretizing a three-dimensional transient heat conduction differential equation of each steel billet on the basis of the previous iteration time step by taking the radiation heat exchange heat flow density and the convection heat flow coefficient as steel billet boundary conditions, calculating three-dimensional transient heat conduction, updating the three-dimensional temperature field of the steel billet, and then transmitting the surface temperature of the steel billet to the furnace section three-dimensional heat transfer model to calculate the updated radiation heat flow density and convection heat flow coefficient.
In the process of establishing the steel billet three-dimensional heat conduction model, the steel billet chamfer is ignored to be simplified into a standard hexahedron, and the heat conductivity of the steel billet in all directions is assumed to be consistent. And considering convection between furnace gas and billet in billet boundary conditions, radiation between furnace atmosphere-furnace lining-beam-billet, etc.
The invention also comprises a parameter reading model which is used for reading the whole furnace parameter information contained in the parameter file, carrying out data interpolation according to the current calculation time and each billet position, and transmitting the whole furnace parameter information after interpolation to a whole furnace parameter scheduling model in real time.
As shown in fig. 1, the method for constructing the prediction model of the three-dimensional heat transfer temperature field of the heating furnace billet comprises the following steps:
S1, reading parameter information of a whole furnace by using a parameter reading model, dividing the heating furnace into a plurality of furnace sections by using a whole furnace parameter calling model based on the parameter information of the whole furnace, wherein each furnace section corresponds to an angle coefficient calculation file, has furnace temperature and convection heat transfer coefficients as heat transfer attributes, and is set to comprise a plurality of billets, so that the complexity of radiation heat transfer calculation is simplified;
S2, calculating angle coefficients between surfaces in the heating furnace by using a Monte Carlo method, and generating an angle coefficient file;
The step S2 specifically comprises the following steps:
Because of the complexity of gas radiation, in order to simplify the calculation, firstly, when solving the radiation heat exchange of the radiation participating gas in the closed space in the heating furnace, the following assumption is made:
Assuming that the enclosed space in the heating furnace is filled with radiation participation gas, and the spatial distribution of the concentration, the temperature and the pressure is uniform, the radiation participation gas is ash body irrespective of the scattering and reflection characteristics of the gas, and the sum of the absorptivity Ag and the transmissivity Dg of the thermal radiation is 1:
(1)
Meanwhile, the N surfaces forming the closed space in the heating furnace are all ash bodies, and the temperature distribution of each surface is uniform;
Then adopting a radiation thermal resistance network method, wherein the temperature of M surfaces in N surfaces in a closed space in a heating furnace is known and is recorded as Ti, i=1-M; the heat flows of the remaining (N-M) surfaces are known and denoted Qi, i=m+1 to N, the form of the radiation heat exchange equation in the furnace is as follows:
(2)
wherein: representing the type of boundary condition,/> When/>When=1, it means that the k surface is a wall temperature boundary condition, when/>When=0, it indicates that the k surface is a constant heat flow boundary condition; δki is a binary function,/>；/>Is the area of surface k; /(I)Emissivity for surface k; /(I)Absorption rate for radiation-participating gases; /(I)Emissivity for radiation-participating gas; /(I)A temperature of the gas involved for radiation; /(I)Is the Stefan Boltzmann constant; /(I)Is the temperature of surface k; /(I)The angular coefficient of k to i.
In the embodiment, the nth billet inside the m furnace sections is meshed according to an external node method, imax, jmax, kmax meshes are respectively meshed along the circumferential direction, the radial direction and the height direction, and equation discretization is performed according to external node meshing.
Finally, solving the formula (2) by using an LU decomposition method to obtain the net radiation heat flow of each surface：
(3)
Wherein: Is effective radiation for each surface in the heating furnace.
S3, solving the radiation heat exchange heat flux density based on the generated angle coefficient file, and estimating the convection heat flux coefficient according to the flow of the high-temperature flue gas in the furnace by using an empirical formula;
S4, taking the radiation heat exchange heat flow density and the convection heat exchange coefficient as steel billet boundary conditions, discretizing a three-dimensional transient heat conduction differential equation of each steel billet, calculating three-dimensional transient heat conduction, and updating a three-dimensional temperature field of the steel billet;
the step S4 specifically comprises the following steps:
Because the internal heat conduction of the steel billet accords with the three-dimensional unsteady heat transfer process without internal heat source, the following three-dimensional transient heat conduction differential equation is firstly set:
(4)
Wherein: t is the temperature of the steel billet, and the unit is K; lambda is the heat conductivity of steel billet, and the unit is ；/>The density of the billet material is Kg/m3; /(I)Specific heat capacity of billet material is expressed as/>; T is heating time, and the unit is s;
Setting the initial temperature of the billet heat transfer model to be uniform and the environment temperature, and setting the initial temperature in a heating furnace, wherein the heat exchange condition of the billet surface is a third type of boundary condition, namely a composite boundary condition of convection heat exchange and radiation heat exchange, and the boundary condition of the billet is as follows:
(5)
wherein: the unit is W/(m2.K) which is the convection heat transfer coefficient; /(I) The radiation heat exchange heat flow density is W/m < 2 >; the furnace temperature is K; /(I) The unit is K, which is the surface temperature of the billet;
and finally discretizing the formula (5) by adopting an implicit differential method, and solving by adopting a Gausserdel iteration method.
And S5, feeding the calculated three-dimensional temperature field of the steel billet back to the step S3, and updating the radiation heat exchange heat flow density and the convection heat flow coefficient until the iteration condition is met, and outputting the three-dimensional temperature field of the steel billet according to the approval time interval.
In step S5, the three-dimensional temperature field of the billet includes a billet temperature field cloud image and a temperature variation graph.
Experimental example
In the experimental example, the implementation effect of a heating furnace for a certain walking beam is shown as in fig. 2-7, and the hit rate of the temperature prediction of the steel billet is higher than 90% within 7.5% according to the verification of a black box coupling embedding experiment performed by the walking beam; the hit rate of the prediction accuracy reaching within 5% is higher than 80%.
Therefore, the heating furnace billet three-dimensional heat transfer temperature field prediction model and the construction method thereof can solve a series of limitations of field experiments, and have the advantages of low cost, high speed, high precision and the like.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present invention and not for limiting it, and although the present invention has been described in detail with reference to the preferred embodiments, it will be understood by those skilled in the art that: the technical scheme of the invention can be modified or replaced by the same, and the modified technical scheme cannot deviate from the spirit and scope of the technical scheme of the invention.
Claims (8)
1. The heating furnace billet three-dimensional heat transfer temperature field prediction model is characterized in that: the method comprises a full furnace parameter calling model, a furnace section three-dimensional heat transfer model, a billet three-dimensional heat transfer model and a parameter reading model;
The full furnace parameter calling model is used for calling full furnace parameter information, wherein the full furnace parameter comprises furnace type structure information, furnace inner beam position information, furnace temperature information, steel loading information and steel tapping information;
The furnace section three-dimensional heat transfer model is used for calculating the radiation heat exchange heat flow density among the billet, the furnace gas and the furnace lining in each furnace section according to the whole furnace parameter information; calculating the convective heat transfer coefficient of the surface of the steel billet in each furnace section according to an empirical formula;
The steel billet three-dimensional heat conduction model is used for discretizing a three-dimensional transient heat conduction differential equation of each steel billet on the basis of the previous iteration time step by taking the radiation heat exchange heat flow density and the convection heat flow coefficient as steel billet boundary conditions, calculating three-dimensional transient heat conduction, updating the three-dimensional temperature field of the steel billet, and then transmitting the surface temperature of the steel billet to the furnace section three-dimensional heat transfer model to calculate the updated radiation heat flow density and convection heat flow coefficient.
2. The heating furnace billet three-dimensional heat transfer temperature field prediction model according to claim 1, wherein: the system also comprises a parameter reading model, which is used for reading the whole furnace parameter information contained in the parameter file, carrying out data interpolation according to the current calculation time and each billet position, and transmitting the interpolated whole furnace parameter information to a whole furnace parameter scheduling model in real time.
3. The heating furnace billet three-dimensional heat transfer temperature field prediction model according to claim 2, wherein: when used in a walking beam heating furnace, the full furnace parameter information also comprises walking beam stepping information, is used for updating the positions of various billets in the heating furnace and re-associating the radiation angle coefficient files of various furnace sections.
4. A furnace billet three-dimensional heat transfer temperature field prediction model according to claim 3, wherein: in the full furnace parameter call model, furnace type structure information includes: after dividing the furnace sections into a set number, the furnace section size of each furnace section and the number, size and position information of the cross beams and the longitudinal beams;
The furnace temperature information is furnace temperature parameters assigned to each furnace section by using a furnace temperature distribution model;
The steel loading information is a specific position of a billet in the heating furnace determined by steel type, specification and placement position information of the steel loading, and a radiation angle coefficient file associated with the specific position;
the tapping information is stored tapping billet information.
5. The method for constructing the heating furnace billet three-dimensional heat transfer temperature field prediction model according to claim 4, wherein the method comprises the following steps: the method comprises the following steps:
S1, reading parameter information of a whole furnace by using a parameter reading model, dividing a heating furnace into a plurality of furnace sections by using a whole furnace parameter calling model based on the parameter information of the whole furnace, and setting each furnace section to comprise a plurality of billets;
S2, calculating angle coefficients between surfaces in the heating furnace by using a Monte Carlo method, and generating an angle coefficient file;
S3, solving the radiation heat exchange heat flux density based on the generated angle coefficient file, and estimating the convection heat flux coefficient according to the flow of the high-temperature flue gas in the furnace by using an empirical formula;
S4, taking the radiation heat exchange heat flow density and the convection heat exchange coefficient as steel billet boundary conditions, discretizing a three-dimensional transient heat conduction differential equation of each steel billet, calculating three-dimensional transient heat conduction, and updating a three-dimensional temperature field of the steel billet;
And S5, feeding the calculated three-dimensional temperature field of the steel billet back to the step S3, and updating the radiation heat exchange heat flow density and the convection heat flow coefficient until the iteration condition is met, and outputting the three-dimensional temperature field of the steel billet according to the approval time interval.
6. The method for constructing the heating furnace billet three-dimensional heat transfer temperature field prediction model according to claim 5, wherein the method comprises the following steps of: the step S2 specifically comprises the following steps:
Firstly, when solving the radiation heat exchange of the radiation participation gas in the closed space in the heating furnace, the following assumption is made:
Assuming that the enclosed space in the heating furnace is filled with radiation participation gas, and the spatial distribution of the concentration, the temperature and the pressure is uniform, the radiation participation gas is ash body irrespective of the scattering and reflection characteristics of the gas, and the sum of the absorptivity Ag and the transmissivity Dg of the thermal radiation is 1:
(1)
Meanwhile, the N surfaces forming the closed space in the heating furnace are all ash bodies, and the temperature distribution of each surface is uniform;
Then adopting a radiation thermal resistance network method, wherein the temperature of M surfaces in N surfaces in a closed space in a heating furnace is known and is recorded as Ti, i=1-M; the heat flows of the remaining (N-M) surfaces are known and denoted Qi, i=m+1 to N, the form of the radiation heat exchange equation in the furnace is as follows:
(2)
wherein: representing the type of boundary condition,/> When/>When=1, it means that the k surface is a wall temperature boundary condition, when/>When=0, it indicates that the k surface is a constant heat flow boundary condition; δki is a binary function,/>；/>Is the area of surface k; /(I)Emissivity for surface k; /(I)Absorption rate for radiation-participating gases; /(I)Emissivity for radiation-participating gas; /(I)A temperature of the gas involved for radiation; /(I)Is the Stefan Boltzmann constant; /(I)Is the temperature of surface k; /(I)An angular coefficient of k to i;
Finally, solving the formula (2) by using an LU decomposition method to obtain the net radiation heat flow of each surface ：
(3)
Wherein: Is effective radiation for each surface in the heating furnace.
7. The method for constructing the heating furnace billet three-dimensional heat transfer temperature field prediction model according to claim 6, wherein the method is characterized in that: the step S4 specifically comprises the following steps:
Firstly, setting the following three-dimensional transient heat conduction differential equation:
(4)
Wherein: t is the temperature of the steel billet, and the unit is K; lambda is the heat conductivity of steel billet, and the unit is ；/>The density of the billet material is Kg/m3; /(I)Specific heat capacity of billet material is expressed as/>; T is heating time, and the unit is s;
Setting the initial temperature of the billet heat transfer model to be uniform and the environment temperature, setting the initial temperature in a heating furnace, and setting the heat exchange condition of the billet surface to be a third type boundary condition, wherein the billet boundary condition is as follows:
(5)
wherein: the unit is W/(m2.K) which is the convection heat transfer coefficient; /(I) The radiation heat exchange heat flow density is W/m < 2 >; /(I)The furnace temperature is K; /(I)The unit is K, which is the surface temperature of the billet;
and finally discretizing the formula (5) by adopting an implicit differential method, and solving by adopting a Gausserdel iteration method.
8. The method for constructing the heating furnace billet three-dimensional heat transfer temperature field prediction model according to claim 7, wherein the method comprises the following steps of: in step S5, the three-dimensional temperature field of the billet includes a billet temperature field cloud image and a temperature variation graph.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410455353.4A CN118052086A (en) | 2024-04-16 | 2024-04-16 | Heating furnace billet three-dimensional heat transfer temperature field prediction model and construction method thereof |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410455353.4A CN118052086A (en) | 2024-04-16 | 2024-04-16 | Heating furnace billet three-dimensional heat transfer temperature field prediction model and construction method thereof |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN118052086A true CN118052086A (en) | 2024-05-17 |
Family
ID=91048700
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410455353.4A Pending CN118052086A (en) | 2024-04-16 | 2024-04-16 | Heating furnace billet three-dimensional heat transfer temperature field prediction model and construction method thereof |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118052086A (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102994730A (en) * | 2012-11-30 | 2013-03-27 | 中冶南方(武汉)威仕工业炉有限公司 | Method for tracking temperature of steel billet in heating furnace |
| CN106906350A (en) * | 2017-02-10 | 2017-06-30 | 中冶华天南京工程技术有限公司 | A kind of H profile steel base heating process temperature distribution calculation method |
| CN110348174A (en) * | 2019-08-07 | 2019-10-18 | 中冶赛迪技术研究中心有限公司 | A kind of steel billet temperature calculation method of heating furnace |
| CN113191080A (en) * | 2021-04-26 | 2021-07-30 | 辽宁省交通高等专科学校 | Heating furnace billet temperature field prediction model optimization method based on HMPSO algorithm |
-
2024
- 2024-04-16 CN CN202410455353.4A patent/CN118052086A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102994730A (en) * | 2012-11-30 | 2013-03-27 | 中冶南方(武汉)威仕工业炉有限公司 | Method for tracking temperature of steel billet in heating furnace |
| CN106906350A (en) * | 2017-02-10 | 2017-06-30 | 中冶华天南京工程技术有限公司 | A kind of H profile steel base heating process temperature distribution calculation method |
| CN110348174A (en) * | 2019-08-07 | 2019-10-18 | 中冶赛迪技术研究中心有限公司 | A kind of steel billet temperature calculation method of heating furnace |
| CN113191080A (en) * | 2021-04-26 | 2021-07-30 | 辽宁省交通高等专科学校 | Heating furnace billet temperature field prediction model optimization method based on HMPSO algorithm |
Non-Patent Citations (4)
| Title |
|---|
| 屈成锐, 李义科, 武文斐: "离散坐标法及其在加热炉中的应用", 包头钢铁学院学报, no. 02, 30 June 2004 (2004-06-30) * |
| 沈丙振, 周进, 韩志强, 柳百成, 温亚成, 赵克文, 薛念福: "热轧步进式加热炉内钢坯温度场数值模拟", 冶金能源, no. 04, 30 August 2002 (2002-08-30), pages 1 - 3 * |
| 豆瑞锋 等: "计算机仿真技术在传热传质学 可视化及实践教学中的应用", 高等工程教育研究 2019年增刊I, 31 December 2019 (2019-12-31) * |
| 豆瑞锋 等: "计算机仿真技术在传热传质学可视化及实践教学中的应用", 高等工程教育研究 2019年增刊I, 31 December 2019 (2019-12-31) * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Chugh et al. | A data-driven surrogate-assisted evolutionary algorithm applied to a many-objective blast furnace optimization problem | |
| Huang et al. | A proposed iteration optimization approach integrating backpropagation neural network with genetic algorithm | |
| KR19990074013A (en) | Method and apparatus for surface processing of the outer plate of a ship body | |
| Liao et al. | Hybrid thermal modeling of additive manufacturing processes using physics-informed neural networks for temperature prediction and parameter identification | |
| Pal et al. | Optimization of Phosphorous in Steel Produced by Basic Oxygen Steel Making Process Using Multi‐Objective Evolutionary and Genetic Algorithms | |
| Gomez et al. | Natural gas intermittent kiln for the ceramic industry: A transient thermal analysis | |
| Mojumder et al. | Linking process parameters with lack-of-fusion porosity for laser powder bed fusion metal additive manufacturing | |
| Nikam et al. | A simplified thermal approximation method to include the effects of Marangoni convection in the melt pools of processes that involve moving point heat sources | |
| CN118052086A (en) | Heating furnace billet three-dimensional heat transfer temperature field prediction model and construction method thereof | |
| Wang et al. | A novel time step fusion method with finite volume formulation for accelerated thermal analysis of laser additive manufacturing | |
| Xie et al. | Bayesian inverse uncertainty quantification of a MOOSE-based melt pool model for additive manufacturing using experimental data | |
| CN113515873B (en) | Metal additive manufacturing molten pool shape prediction method based on dimensional analysis | |
| Rincon-Tabares et al. | Sensitivity analysis for transient thermal problems using the complex-variable finite element method | |
| Pereira et al. | A multiobjective optimization parameters applied to additive manufacturing: DOE-based approach to 3D printing | |
| Li et al. | Study on parameter optimization of laser cladding Fe60 based on GA-BP neural network | |
| Kumar et al. | Modeling moving-boundary problems of solidification and melting adopting an arbitrary Lagrangian-Eulerian approach | |
| Schönegger et al. | Computational welding simulation of a plasma wire arc additive manufacturing process for high-strength steel | |
| CN116246724A (en) | Zinc rotary kiln key monitoring variable soft measurement method, device, terminal and medium | |
| Mishra et al. | Development and application of a simplified thermo-evaporative multi-track laser powder bed fusion model for Ti6Al4V | |
| CN118052087B (en) | Online and offline digital twin model synchronous construction method based on three-dimensional heat transfer | |
| CN114548526A (en) | Satellite component layout temperature field prediction method based on physical prior neural network | |
| Shahane | Numerical simulations of die casting with uncertainty quantification and optimization using neural networks | |
| CN113011068A (en) | Three-dimensional simulation method for walking beam type plate blank heating | |
| Elmisaoui et al. | A sparse polynomial surrogate model for the shrinking core model in phosphate ore digestion | |
| Sata | Investment casting defect prediction using neural network and multivariate regression along with principal component analysis |
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 |