Three-dimensional temperature field dynamic modeling method used in steel plate production process

Abstract

The invention discloses a three-dimensional temperature field dynamic modeling method used in a steel plate production process, and relates to the field of three-dimensional temperature field modeling. The method comprises the following steps: discretizing a cross section area of a steel plate, and dividing each cross section into a plurality of rectangular grid units; calculating the temperature of each grid node on the cross section by using a finite difference model based on the detection position of the existing temperature measurement point and the temperature real-time value corresponding to each temperature measurement point, and constructing a two-dimensional temperature field corresponding to each cross section; based on the two-dimensional temperature field on the cross section, space expansion is carried out in the length direction of the steel plate, and a three-dimensional temperature field initial model is constructed; calculating a temperature correction value corresponding to each grid node on the cross section of the steel plate in the three-dimensional temperature field initial model; according to the temperature correction value, calculating a temperature correction amount used for reflecting time-distance dimension difference, and correcting the grid node temperature calculated based on the finite difference model through the temperature correction amount; according to the method, refined modeling of the three-dimensional temperature field of the steel plate is realized.

Description

Technical Field

[0001] The present invention relates to the field of three-dimensional temperature field modeling, and in particular to a dynamic three-dimensional temperature field modeling method for the steel plate production process. Background Art

[0002] Temperature is a core process parameter in the hot-rolling steel plate production process and has a decisive impact on product quality. The uniformity of the steel plate temperature field is directly related to the flow stress characteristics of the metal during rolling, which in turn affects thickness control and shape quality; while in the cooling stage, the temperature gradient determines the microstructure evolution and final mechanical properties of the material, such as strength and toughness. In addition, temperature changes are also closely related to key process issues such as scale formation, roll thermal fatigue, and wear, significantly affecting production costs and energy consumption control. Therefore, accurately grasping the three-dimensional temperature distribution of the steel plate at each process stage is of great significance for optimizing the production process, improving product dimensional accuracy, and mechanical properties.

[0003] However, the steel plate rolling process is a complex thermo-mechanical coupling process, affected by various process parameters and environmental factors, making temperature field modeling face great challenges. Existing temperature measurement methods have limitations: contact temperature measurement has a response lag and limited installation, making it difficult to meet the dynamic monitoring requirements of high-speed production lines; non-contact infrared temperature measurement has a fast response, but is easily interfered by water mist, scale, and emissivity fluctuations, with large measurement errors, and cannot obtain complete temperature information in the width direction. At the same time, traditional temperature field modeling is mostly based on one-dimensional or two-dimensional simplified models, making it difficult to comprehensively reflect the complex heat transfer behavior of the steel plate in three-dimensional space. Therefore, there is an urgent need to develop a high-precision three-dimensional temperature field modeling method that can integrate measured data and perform dynamic correction to achieve accurate prediction and real-time update of the steel plate temperature distribution, and promote the development of the hot-rolling process towards intelligence, high quality, and low energy consumption. Summary of the Invention

[0004] In order to achieve accurate prediction and real-time update of the steel plate temperature distribution, the present invention proposes a dynamic three-dimensional temperature field modeling method for the steel plate production process, including: Based on the temperature measurement points arranged on the steel plate production line, collect the real-time temperature data of the steel plate at each process stage; the process stages include: heating, rough rolling, finish rolling, and cooling; Preprocess the real-time temperature data to obtain the real-time temperature value; Construct a two-dimensional temperature field: discretize the cross-sectional area of the steel plate, and divide each cross-section into several rectangular grid units; the cross-section is composed of the width direction and thickness direction of the steel plate; based on the detection positions of the existing temperature measurement points and the real-time temperature values corresponding to each temperature measurement point, use the finite difference model to calculate the temperatures of each grid node on the cross-section, and construct the two-dimensional temperature field corresponding to each cross-section accordingly; The method of varying the grid size is adopted to set the distance between adjacent cross-sections at each process stage as the spatial step length in the length direction of the steel plate. Based on the two-dimensional temperature field on the cross-section and the spatial step length corresponding to the current process stage, spatial expansion is carried out along the length direction of the steel plate to construct an initial three-dimensional temperature field model; According to the type of temperature measurement equipment corresponding to the temperature measurement point and the position of the temperature measurement point, calculate the temperature correction value corresponding to each grid node on the cross-section of the steel plate in the initial three-dimensional temperature field model; According to the distance between any two temperature measurement devices, the temperature detection time and the running speed of the steel plate, obtain the moments when the cross-section of the steel plate reaches the temperature measurement points corresponding to the two temperature measurement devices during the production process of the steel plate. According to these moments and the temperature correction values, calculate the temperature correction amount used to reflect the time-distance dimension difference, and correct the temperature of the grid nodes calculated based on the finite difference model through the temperature correction amount to obtain the target three-dimensional temperature field corresponding to the current moment; Based on the latest two-dimensional temperature field, continuously update the target three-dimensional temperature field, thereby constructing a dynamically evolving three-dimensional temperature field model.

[0005] Furthermore, the construction of the two-dimensional temperature field is specifically as follows: Obtain the model parameters of the steel plate; Discretize the cross-sectional area of the steel plate, and divide each cross-section into several rectangular grid units; The cross-section is composed of the width direction and the thickness direction of the steel plate, and determine the spatial step length ΔX in the width direction and the spatial step length ΔY in the thickness direction; According to the spatial step lengths ΔX and ΔY, and the material thermal property parameters of the steel plate, set the stability condition, and based on the stability condition, set the time step length Δt; Use the detection positions of the existing temperature measurement points and the real-time temperature values corresponding to each temperature measurement point at the set moment to perform preliminary interpolation on the grid nodes; Based on the heat transfer control equation, use the model parameters of the steel plate and the real-time temperature values, and adopt the finite difference model to calculate the temperature of each grid node on the cross-section at each moment under the set time step length Δt; Based on the temperature of each grid node, use the interpolation method to calculate the temperature of each point inside the grid unit; Based on the temperature of each grid node and the temperature of each point inside it, construct the two-dimensional temperature field on the cross-section corresponding to each moment.

[0006] Furthermore, the construction of the two-dimensional temperature field is specifically as follows: Discretize the cross-section of the steel plate, divide each cross-section into several rectangular grid cells, and use the detection positions of the existing temperature measurement points and the real-time temperature values corresponding to each temperature measurement point at the set time to perform preliminary interpolation on the grid nodes, calculate the correction coefficients of each grid node in the width direction and thickness direction of the steel plate, adjust the real-time temperature values of the corresponding grid nodes through the correction coefficients, and calculate the temperatures of each grid node on the cross-section using a finite difference model based on the adjusted temperature data, and construct a two-dimensional temperature field on the cross-section accordingly.

[0007] Further, the calculation of the correction coefficients of each grid node in the width direction and thickness direction of the steel plate, adjusting the real-time temperature values of the corresponding grid nodes through the correction coefficients, calculating the temperatures of each grid node on the cross-section using a finite difference model based on the adjusted temperature data, and constructing a two-dimensional temperature field on the cross-section accordingly are specifically as follows: Based on the temperature non-uniformity of the grid nodes in the width direction, construct a transverse correction function; Based on the temperature gradient of the grid nodes in the thickness direction, construct a vertical correction function; Use the transverse correction function and the vertical correction function to calculate the transverse correction coefficient of the grid nodes in the width direction of the steel plate and the vertical correction coefficient in the thickness direction of the steel plate; Correct the real-time temperature values of the grid nodes through the transverse correction coefficient and the vertical correction coefficient to obtain the temperature adjustment value; Based on the heat transfer control equation, use the model parameters of the steel plate and the temperature adjustment value, and calculate the temperatures of each grid node on the cross-section using a finite difference model; Based on the temperature adjustment value and the calculated temperatures of the grid nodes, use the interpolation algorithm to calculate the temperature values of the interpolation points inside the grid cells; Combine the temperature adjustment values of the grid nodes, the calculated temperatures of the grid nodes and the temperature values of the internal interpolation points to construct a two-dimensional temperature field on the cross-section.

[0008] Further, the model parameters of the steel plate include: preset geometric dimension information, material thermophysical properties parameters and convective heat transfer coefficients under various boundary conditions; the boundary conditions include: Dirichlet boundary condition and Neumann boundary condition.

[0009] Further, the types of temperature measurement devices include single-point pyrometers and / or transverse scanning pyrometers; according to the types of temperature measurement devices corresponding to the temperature measurement points and the positions of the temperature measurement points, calculate the temperature correction values corresponding to each grid node on the cross-section of the steel plate in the initial model of the three-dimensional temperature field, specifically as follows: When there are only temperature measurement points arranged by single-point pyrometers on the upper surface of the steel plate, calculate the temperature correction values corresponding to each grid node in the width direction of the cross-section of the steel plate in the initial model of the three-dimensional temperature field through the first single-point temperature correction model; When temperature measurement points are arranged on both the upper surface and the lower surface of the steel plate through single-point pyrometers, the temperature correction values corresponding to each grid node in the thickness direction of the cross-section of the steel plate in the initial three-dimensional temperature field model are calculated through the second single-point temperature correction model; When temperature measurement points are arranged on the surface of the steel plate through a transverse scanning pyrometer, the temperature correction values corresponding to each grid node in the width direction and the thickness direction of the cross-section of the steel plate in the initial three-dimensional temperature field model are calculated through the transverse temperature correction model.

[0010] Further, the formula expression of the first single-point temperature correction model is: ; where: ; In the formula, represents the temperature correction amount at the center point of the upper surface of the steel plate, which is used to compensate for the deviation between the calculated value of the finite difference model and the measured temperature value; represents the grid node temperature of the grid node at the i-th column in the width direction and the j-th row in the thickness direction of the cross-section of the steel plate at the k-th moment calculated based on the finite difference model; represents the temperature correction value of the grid node at the i-th column in the width direction and the j-th row in the thickness direction of the cross-section of the steel plate at the k-th moment; represents the measured temperature value at the center point of the upper surface of the steel plate, where represents measured; n represents the total number of grid nodes in the width direction of the cross-section of the steel plate, and m represents the total number of grid nodes in the thickness direction of the cross-section of the steel plate; represents the grid node temperature of the grid node at the center point of the upper surface of the steel plate, i.e., the -th column and the m-th row at the k-th moment calculated through the finite difference model.

[0011] Further, the formula expression of the second single-point temperature correction model is: ; where: ; ; In the formula, represents the temperature correction value of the grid node at the i-th column in the width direction and the j-th row in the thickness direction of the cross-section of the steel plate at the k-th moment; represents the grid node temperature of the grid node at the i-th column in the width direction and the j-th row in the thickness direction of the cross-section of the steel plate at the k-th moment calculated based on the finite difference model; represents the temperature correction amount at the center point of the upper surface of the steel plate; Represents the temperature correction amount at the center point of the lower surface of the steel plate; Represents the distance between adjacent grid nodes along the thickness direction of the steel plate in the cross-section of the steel plate; Represents the coordinate position of the current grid node in the thickness direction; Represents the total height of the steel plate, which is equal to the sum of the distances between all grid nodes in the thickness direction; Represents the measured temperature value at the center point of the upper surface of the steel plate, where, Represents measurement; Represents the grid node temperature at the center point of the upper surface of the steel plate, i.e., the grid node at the nth column and the mth row, calculated by the finite difference model at the kth moment; n represents the total number of grid nodes in the width direction of the cross-section of the steel plate, and m represents the total number of grid nodes in the thickness direction of the cross-section of the steel plate; Represents the grid node temperature at the center point of the upper surface of the steel plate, i.e., the grid node at the nth column and the mth row, calculated by the finite difference model at the kth moment; n represents the total number of grid nodes in the width direction of the cross-section of the steel plate, and m represents the total number of grid nodes in the thickness direction of the cross-section of the steel plate; Represents the measured temperature value at the center point of the lower surface of the steel plate; Represents the grid node temperature at the center point of the lower surface of the steel plate, i.e., the grid node at the nth column and the 0th row, calculated by the finite difference model at the kth moment. Represents the grid node temperature at the center point of the lower surface of the steel plate, i.e., the grid node at the nth column and the 0th row, calculated by the finite difference model at the kth moment.

[0012] Furthermore, the formula expression of the transverse temperature correction model is: ; where: ; In the formula, Represents the temperature correction value of the grid node at the ith column in the width direction and the jth row in the thickness direction of the cross-section of the steel plate at the kth moment; Represents the grid node temperature of the grid node at the ith column in the width direction and the jth row in the thickness direction of the cross-section of the steel plate calculated based on the finite difference model at the kth moment; Represents the transverse overall correction amount corresponding to the ith column in the width direction; Represents the measured temperature value of the grid node at the ith column in the width direction and the topmost layer, i.e., the mth row, in the thickness direction of the upper surface of the steel plate; Represents the grid node temperature of the grid node at the ith column in the width direction and the mth row in the thickness direction of the cross-section of the steel plate calculated based on the finite difference model at the kth moment.

[0013] Furthermore, the calculation formula used to calculate the temperature correction amount reflecting the time-distance dimension difference according to these moments and temperature correction values is: ; In the formula, Represents the moment when the cross-section of the steel plate reaches the first temperature measurement point S, Indicates the moment when the cross-section of the steel plate reaches the second temperature measurement point E. Indicates the moment when the cross-section of the steel plate reaches the set position between the two temperature measurement points ; Indicates the difference between the temperature correction value of the grid node at the position on the cross-section of the steel plate where the width direction is the i-th column and the thickness direction is the j-th row and the temperature of the grid node calculated by the corresponding finite difference model; Indicates the temperature correction amount corresponding to the grid node at the position on the cross-section of the steel plate where the width direction is the i-th column and the thickness direction is the j-th row; The temperature of the grid node calculated by the finite difference model is corrected by the temperature correction amount, and the calculation formula used is: ; In the formula, Indicates the grid node temperature of the grid node at the position on the cross-section of the steel plate where the width direction is the i-th column and the thickness direction is the j-th row at the k-th moment calculated by the finite difference model; Indicates the final temperature value after dimension correction of the grid node at the position on the cross-section of the steel plate where the width direction is the i-th column and the thickness direction is the j-th row at the k-th moment.

[0014] Compared with the prior art, the present invention has at least the following beneficial effects:

[0015] (1) In the present invention, the cross-sectional area of the steel plate is discretized, and each cross-section is divided into a number of rectangular grid units; based on the detection positions of the existing temperature measurement points and the real-time temperature values corresponding to each temperature measurement point, the finite difference model is used to calculate the temperature of each grid node on the cross-section, and the two-dimensional temperature field corresponding to each cross-section is constructed based on this; based on the two-dimensional temperature field on the cross-section and the space step corresponding to the current process stage, spatial expansion is carried out along the length direction of the steel plate to construct an initial three-dimensional temperature field model; according to the type of temperature measurement equipment corresponding to the temperature measurement point and the position of the temperature measurement point, the temperature correction value corresponding to each grid node on the cross-section of the steel plate in the initial three-dimensional temperature field model is calculated; the moments when the cross-section of the steel plate reaches the temperature measurement points corresponding to the two temperature measurement devices during the production process of the steel plate are obtained, and the temperature correction amount reflecting the time-distance dimension difference is calculated based on these moments and the temperature correction value, and the temperature of the grid node calculated by the finite difference model is corrected by the temperature correction amount; this method realizes the refined modeling of the three-dimensional temperature field of the steel plate, solves the problems that traditional one-dimensional or two-dimensional models are difficult to accurately reflect complex heat transfer behaviors, have large temperature measurement errors, and poor model adaptability, significantly improves the modeling accuracy of the temperature field and the dynamic response ability, and provides reliable data support for subsequent process control and product quality improvement.

[0016] (2) The present invention calculates the temperature correction values of each grid node according to the type and position information of the temperature measurement device corresponding to the temperature measurement point. This process not only takes into account the characteristics of different temperature measurement devices (such as single-point pyrometers or transverse scanning pyrometers), but also performs precise correction calculations for the specific position of each grid node, significantly improving the accuracy of the target three-dimensional temperature field construction.

[0017] (3) Based on the distance between any two temperature measurement devices, the temperature detection time, and the running speed of the steel plate, the present invention obtains the moments when the cross-section of the steel plate reaches the temperature measurement points corresponding to the two temperature measurement devices, and calculates the temperature correction amount for reflecting the time-distance dimension difference according to these moments and the temperature correction values. This correction amount dynamically adjusts the calculation results of the finite difference model through a physically-driven time ratio interpolation method, making the final temperature field model closer to the temperature distribution under actual production conditions. This multi-step and multi-level temperature correction mechanism ensures that the model can not only accurately reflect the initial temperature field, but also continuously update and optimize with the changes in the production line, achieving high-precision simulation of the complex heat transfer behavior during the hot rolling process. Therefore, the present invention solves the problem of inaccurate temperature field modeling caused by large measurement errors and poor model adaptability in traditional methods, and greatly improves the product quality control ability and production process optimization level.

[0018] (4) By introducing a temperature correction amount based on the time-distance relationship, the present invention realizes a reasonable estimation and compensation of temperature deviations at different positions. This time ratio interpolation method based on the physical process has more theoretical basis and engineering applicability than simple linear or empirical corrections.

[0019] (5) By superimposing the temperature correction amount obtained based on the measured data on the calculation results of the finite difference model to form the corrected target temperature value, the present invention effectively integrates the continuity of the finite difference model and the authenticity of the measured data, improving the overall reliability of the target three-dimensional temperature field. Description of the Drawings

[0020] Figure 1 It is a flowchart of a three-dimensional temperature field dynamic modeling method for the steel plate production process according to an embodiment of the present invention. Detailed Embodiment

[0021] The following are specific embodiments of the present invention in combination with the drawings to further describe the technical solutions of the present invention, but the present invention is not limited to these embodiments.

[0022] In order to achieve accurate prediction and real-time update of the temperature distribution of the steel plate, as Figure 1 shown, the present invention proposes a three-dimensional temperature field dynamic modeling method for the steel plate production process, including: Based on the temperature measurement points arranged on the steel plate production line, collect the real-time temperature data of the steel plate in each process stage; the process stages include: heating, rough rolling, finish rolling, and cooling; The temperature measurement equipment arranged on the steel plate production line in this embodiment is shown in Table 1 below: Table 1:

[0023] Among them: R1 represents the first rough rolling mill, and R2 represents the second rough rolling mill.

[0024] Preprocess the real-time temperature data to obtain the real-time temperature value; Construct a two-dimensional temperature field: discretize the cross-sectional area of the steel plate, and divide each cross-section into several rectangular grid units; the cross-section is composed of the width direction and the thickness direction of the steel plate; based on the detection positions of the existing temperature measurement points and the real-time temperature values corresponding to each temperature measurement point, use the finite difference model to calculate the temperatures of each grid node on the cross-section, and construct the two-dimensional temperature field corresponding to each cross-section accordingly; It should be noted that the construction of the two-dimensional temperature field in this embodiment includes two implementation methods, among which: The first implementation method is: Specifically, the construction of the two-dimensional temperature field is as follows: Obtain the model parameters of the steel plate; the model parameters of the steel plate include: preset geometric dimension information, material thermophysical parameters, and convective heat transfer coefficients under various boundary conditions; the boundary conditions include: Dirichlet boundary conditions and Neumann boundary conditions.

[0025] Discretize the cross-sectional area of the steel plate, and divide each cross-section into several rectangular grid units; the cross-section is composed of the width direction and the thickness direction of the steel plate, and determine the spatial step size ΔX in the width direction and the spatial step size ΔY in the thickness direction; According to the spatial step sizes ΔX and ΔY, and the material thermophysical parameters of the steel plate, set the stability condition, and based on the stability condition, set the time step size Δt; The formula expression of the stability condition is: ; In the formula, represents the material thermal diffusivity.

[0026] Before performing the finite difference calculation, set the time step size Δt according to the set stability condition, so as to meet the stability requirements of numerical calculation and avoid numerical instability phenomena during the temperature field simulation process.

[0027] Use the detection positions of the existing temperature measurement points and the real-time temperature values corresponding to each temperature measurement point at the set moment to perform preliminary interpolation on the grid nodes; Based on the heat transfer control equation, using the model parameters of the steel plate and the real-time temperature values, the finite difference model is used to calculate the temperature of each grid node on the cross-section at each moment under the set time step Δt; Based on the temperature of each grid node, the interpolation method is used to calculate the temperature of each point inside the grid element; Based on the temperature of each grid node and the temperature of each point inside it, a two-dimensional temperature field corresponding to the cross-section at each moment is constructed.

[0028] The first implementation method first obtains the geometric dimension information of the steel plate, the thermophysical properties of the material, and the convective heat transfer coefficient under the boundary conditions, ensuring that the model can accurately reflect the actual physical characteristics of the steel plate. By determining the spatial step sizes in the width direction and the thickness direction and setting the time step size that satisfies the stability condition, the finite difference model has good stability and accuracy during the calculation process. This method not only provides a solid foundation for subsequent temperature field calculations but also avoids numerical instability or calculation errors caused by inappropriate step size settings. Then, using preliminary interpolation and the heat transfer control equation, the temperature of each grid node on the cross-section is calculated at each moment, and then the two-dimensional temperature field corresponding to each moment is constructed. This method based on actual physical parameters and a strict mathematical model can accurately capture the temperature changes of the steel plate in different process stages, improving the accuracy and reliability of temperature field modeling. At the same time, the calculation process at each moment ensures the continuity and consistency of the dynamic change of the temperature field over time, further enhancing the practicality and prediction ability of the model.

[0029] The second implementation method is as follows: The specific construction of the two-dimensional temperature field is as follows: The cross-section of the steel plate is discretized, and each cross-section is divided into several rectangular grid cells. Using the detection positions of the existing temperature measurement points and the real-time temperature values corresponding to each temperature measurement point at the set moment, preliminary interpolation is performed on the grid nodes, the correction coefficients in the width direction and the thickness direction of the steel plate for each grid node are calculated, the real-time temperature values of the corresponding grid nodes are adjusted through the correction coefficients, and based on the adjusted temperature data, the finite difference model is used to calculate the temperature of each grid node on the cross-section, and a two-dimensional temperature field on the cross-section is constructed accordingly.

[0030] The calculation of the correction coefficients in the width direction and the thickness direction of the steel plate for each grid node, the adjustment of the real-time temperature values of the corresponding grid nodes through the correction coefficients, the use of the finite difference model to calculate the temperature of each grid node on the cross-section based on the adjusted temperature data, and the construction of a two-dimensional temperature field on the cross-section accordingly are specifically as follows: Based on the temperature non-uniformity of the grid nodes in the width direction, a transverse correction function is constructed; The formula expression of the transverse correction function is: ; In the formula, represents the lateral correction coefficient of the i-th column grid node; represents the column index of the grid node in the width direction; represents the total number of grid nodes of the steel plate cross-section in the width direction; is the temperature non-uniformity correction factor in the width direction (set according to the actual situation, usually an empirical parameter).

[0031] Based on the temperature gradient of the grid node in the thickness direction, a vertical correction function is constructed; The formula expression of the vertical correction function is: ; In the formula, represents the vertical correction coefficient of the j-th row grid node; j represents the row index of the grid node in the thickness direction; m represents the total number of grid nodes of the steel plate cross-section in the thickness direction; represents the temperature gradient correction factor in the thickness direction (set according to the actual situation, usually an empirical parameter).

[0032] Using the lateral correction function and the vertical correction function, calculate the lateral correction coefficient of the grid node in the width direction of the steel plate and the vertical correction coefficient in the thickness direction of the steel plate; Correct the real-time temperature value of the grid node through the lateral correction coefficient and the vertical correction coefficient to obtain the temperature adjustment value; In this embodiment, the calculation formula of the temperature adjustment value is: ; In the formula, represents the real-time temperature value of the grid node where the width direction is the i-th column and the thickness direction is the j-th row; represents the temperature adjustment value corresponding to the grid node where the width direction is the i-th column and the thickness direction is the j-th row.

[0033] In the second embodiment, first, the cross-section of the steel plate is discretized, divided into several rectangular grid cells, and the grid nodes are preliminarily interpolated using the detection positions of the existing temperature measurement points and their corresponding real-time temperature values at the set moments. To accurately reflect the temperature non-uniformity in the width direction and the temperature gradient change in the thickness direction of the steel plate, this method further introduces a transverse correction function and a vertical correction function. Through these correction functions, the transverse correction coefficients in the width direction and the vertical correction coefficients in the thickness direction of each grid node are calculated, and based on this, the temperature values after preliminary interpolation are corrected to obtain temperature adjustment values closer to the actual temperature distribution. Specifically, the transverse correction function is used to adjust the temperature non-uniformity caused by the edge effect in the width direction, while the vertical correction function is used to compensate for the temperature gradient difference in the thickness direction. Based on these corrected temperature data, a finite difference model is used to calculate the temperatures of each grid node on the cross-section, thereby constructing an accurate two-dimensional temperature field. This method ensures that the initial temperature field can accurately reflect the actual temperature distribution characteristics of the steel plate at this moment, providing a reliable basis for the subsequent dynamic evolution of the three-dimensional temperature field modeling.

[0034] Based on the heat transfer control equation, using the model parameters of the steel plate and the temperature adjustment values, a finite difference model is used to calculate the temperatures of each grid node on the cross-section; Based on the temperature adjustment values and the calculated temperatures of the grid nodes, an interpolation algorithm is used to calculate the temperature values of the interpolation points inside the grid cells; Combining the temperature adjustment values of the grid nodes, the calculated temperatures of the grid nodes, and the temperature values of the internal interpolation points, a two-dimensional temperature field on the cross-section is constructed.

[0035] The method of variable grid size is used to set the distance between adjacent cross-sections at each process stage as the spatial step ΔZ in the length direction of the steel plate, and based on the two-dimensional temperature field on the cross-section and the spatial step ΔZ corresponding to the current process stage, spatial expansion is carried out along the length direction of the steel plate to construct an initial model of the three-dimensional temperature field; To simulate the temperature change situation within the entire length range of the steel plate during the entire hot rolling process, it is necessary to expand the above two-dimensional temperature field along the length direction of the steel plate. Specifically, the steps of constructing the two-dimensional temperature field in this embodiment are to repeatedly apply the same two-dimensional temperature field modeling process to different length positions of the steel plate. By splicing multiple such two-dimensional temperature fields along the length direction of the steel plate, a three-dimensional temperature field that can accurately reflect the actual temperature distribution is finally formed.

[0036] In this embodiment, the method of varying the grid size is adopted to set the distance between adjacent cross-sections as the space step in the length direction of the steel plate. Specifically, according to the process characteristics of each process stage and the requirements of the temperature field change, the size of the grid unit in the length direction is dynamically adjusted. For example, in areas where the temperature changes drastically (such as the heating section or the cooling section), a smaller space step is used to improve the model resolution and ensure the accurate capture of the temperature distribution; while in areas where the temperature changes relatively gently (such as the finish rolling section), a larger space step can be used to reduce the computational complexity and improve the computational efficiency. In this way, not only can the spatial resolution of the model be flexibly adjusted in different process stages, but also the requirements of computational accuracy and computational resources can be effectively balanced, thereby constructing a three-dimensional temperature field initial model that more conforms to the actual production situation.

[0037] According to the type of temperature measurement device corresponding to the temperature measurement point and the position of the temperature measurement point, calculate the temperature correction values corresponding to each grid node on the cross-section of the steel plate in the three-dimensional temperature field initial model; The present invention calculates the temperature correction values of each grid node according to the type of temperature measurement device corresponding to the temperature measurement point and the position information. This process not only considers the characteristics of different temperature measurement devices (such as single-point pyrometers or transverse scanning pyrometers), but also performs accurate correction calculations for the specific position of each grid node, significantly improving the accuracy of constructing the target three-dimensional temperature field.

[0038] The type of temperature measurement device includes a single-point pyrometer and / or a transverse scanning pyrometer; According to the type of temperature measurement device corresponding to the temperature measurement point and the position of the temperature measurement point, calculate the temperature correction values corresponding to each grid node on the cross-section of the steel plate in the three-dimensional temperature field initial model.

[0039] According to Table 1, this embodiment sets multiple types of pyrometers, and the second single-point temperature correction model or the transverse temperature correction model can be preferentially selected for calculating the temperature correction values.

[0040] When there are only temperature measurement points arranged by single-point pyrometers on the upper surface of the steel plate, through the first single-point temperature correction model, calculate the temperature correction values corresponding to each grid node in the width direction of the cross-section of the steel plate in the three-dimensional temperature field initial model; The formula expression of the first single-point temperature correction model is: ; where: ; In the formula, represents the temperature correction amount at the center point of the upper surface of the steel plate, which is used to compensate for the deviation between the calculated value of the finite difference model and the measured temperature value; represents the grid node temperature at the k-th moment of the grid node in the i-th column in the width direction and the j-th row in the thickness direction of the steel plate cross-section calculated based on the finite difference model; represents the temperature correction value at the k-th moment of the grid node in the i-th column in the width direction and the j-th row in the thickness direction of the steel plate cross-section; represents the measured temperature value at the center point of the upper surface of the steel plate, where represents measured; n represents the total number of grid nodes in the width direction of the steel plate cross-section, and m represents the total number of grid nodes in the thickness direction of the steel plate cross-section; represents the grid node temperature at the k-th moment of the grid node at the center point of the upper surface of the steel plate, that is, the i-th column and the m-th row grid node calculated by the finite difference model.

[0041] When temperature measurement points are arranged through single-point pyrometers on both the upper and lower surfaces of the steel plate, the temperature correction values corresponding to each grid node in the thickness direction of the steel plate cross-section in the initial three-dimensional temperature field model are calculated through the second single-point temperature correction model; The formula expression of the second single-point temperature correction model is: ; where: ; ; In the formula, represents the temperature correction value at the k-th moment of the grid node in the i-th column in the width direction and the j-th row in the thickness direction of the steel plate cross-section; represents the grid node temperature at the k-th moment of the grid node in the i-th column in the width direction and the j-th row in the thickness direction of the steel plate cross-section calculated based on the finite difference model; represents the temperature correction amount at the center point of the upper surface of the steel plate; represents the temperature correction amount at the center point of the lower surface of the steel plate; represents the distance between adjacent grid nodes in the thickness direction of the steel plate cross-section; represents the coordinate position of the current grid node in the thickness direction; represents the total height of the steel plate, which is equal to the sum of the distances of all grid nodes in the thickness direction; represents the measured temperature value at the center point of the upper surface of the steel plate, where represents measured; represents the grid node temperature at the center point of the upper surface of the steel plate, i.e., the grid node at the nth column and the mth row, at the kth moment, calculated by the finite difference model; n represents the total number of grid nodes in the width direction of the cross-section of the steel plate, and m represents the total number of grid nodes in the thickness direction of the cross-section of the steel plate; represents the measured temperature value at the center point of the lower surface of the steel plate; represents the grid node temperature at the center point of the lower surface of the steel plate, i.e., the grid node at the nth column and the 0th row, at the kth moment, calculated by the finite difference model.

[0042] When there are temperature measurement points arranged on the steel plate surface through a transverse scanning pyrometer, the temperature correction values corresponding to each grid node in the width direction and thickness direction of the cross-section of the steel plate in the initial model of the three-dimensional temperature field are calculated through the transverse temperature correction model.

[0043] The formula expression of the transverse temperature correction model is: ; where: ; In the formula, represents the temperature correction value of the grid node at the ith column in the width direction and the jth row in the thickness direction of the cross-section of the steel plate at the kth moment; represents the grid node temperature of the grid node at the ith column in the width direction and the jth row in the thickness direction of the cross-section of the steel plate at the kth moment, calculated based on the finite difference model; represents the overall transverse correction amount corresponding to the ith column in the width direction; represents the measured temperature value of the grid node at the ith column in the width direction and the topmost layer, i.e., the mth row, in the thickness direction of the upper surface of the steel plate; represents the grid node temperature of the grid node at the ith column in the width direction and the mth row in the thickness direction of the cross-section of the steel plate at the kth moment, calculated based on the finite difference model.

[0044] According to the distance between any two temperature measurement devices, the temperature detection time, and the running speed of the steel plate, the moments when the cross-section of the steel plate reaches the corresponding temperature measurement points of the two temperature measurement devices during the production process of the steel plate are obtained. Based on these moments and the temperature correction values, the temperature correction amount reflecting the time-distance dimension difference is calculated. The grid node temperature calculated based on the finite difference model is corrected by the temperature correction amount to obtain the target three-dimensional temperature field corresponding to the current moment; based on the latest two-dimensional temperature field, the target three-dimensional temperature field is continuously updated, thereby constructing a dynamically evolving three-dimensional temperature field model.

[0045] The temperature correction amount used to reflect the time - distance dimension difference is calculated based on these moments and temperature correction values, and the calculation formula adopted is: ; In the formula, represents the moment when the cross - section of the steel plate reaches the first temperature measurement point S, represents the moment when the cross - section of the steel plate reaches the second temperature measurement point E, represents the moment when the cross - section of the steel plate reaches the set position between the two temperature measurement points; represents the difference between the temperature correction value of the grid node at the position on the cross - section of the steel plate, where the width direction is the i - th column and the thickness direction is the j - th row, and the temperature of the grid node calculated by the corresponding finite - difference model; represents the temperature correction amount corresponding to the grid node at the i - th column in the width direction and the j - th row in the thickness direction on the cross - section of the steel plate; Based on the distance between any two temperature - measuring devices, the temperature - detection time, and the running speed of the steel plate, the present invention obtains the moments when the cross - section of the steel plate reaches the corresponding temperature - measurement points of the two temperature - measuring devices respectively, and calculates the temperature correction amount used to reflect the time - distance dimension difference according to these moments and temperature correction values. This correction amount dynamically adjusts the calculation result of the finite - difference model through a physically - driven time - proportional interpolation method, making the final temperature - field model closer to the temperature distribution under actual production conditions. This multi - step and multi - level temperature correction mechanism ensures that the model can not only accurately reflect the initial temperature field, but also continuously update and optimize with the changes in the production line, achieving a high - precision simulation of the complex heat - transfer behavior during the hot - rolling process. Therefore, the present invention solves the problem of inaccurate temperature - field modeling caused by large measurement errors and poor model adaptability in traditional methods, and greatly improves the product - quality control ability and the production - process optimization level.

[0046] The grid - node temperature calculated based on the finite - difference model is corrected by the temperature correction amount, and the calculation formula adopted is: ; In the formula, represents the grid - node temperature of the grid node at the i - th column in the width direction and the j - th row in the thickness direction on the cross - section of the steel plate at the k - th moment calculated based on the finite - difference model; represents the final temperature value after dimension correction of the grid node at the i - th column in the width direction and the j - th row in the thickness direction on the cross - section of the steel plate at the k - th moment.

[0047] The present invention realizes a reasonable estimation and compensation of temperature deviations at different positions by introducing a temperature correction amount based on the time-distance relationship. This time-ratio interpolation method based on physical processes has more theoretical basis and engineering applicability compared to simple linear or empirical corrections.

[0048] Meanwhile, the present invention forms a corrected target temperature value, i.e., the final temperature value, by superimposing the temperature correction amount obtained based on measured data onto the calculation result of the finite difference model, effectively integrating the continuity of the finite difference model and the authenticity of the measured data, and improving the overall reliability of the target three-dimensional temperature field.

[0049] In this embodiment, the method for obtaining the moment when the steel plate cross-section reaches a set position between two temperature measurement points is as follows:

[0050] It is known that the two temperature measurement points are respectively located at fixed positions on the production line, and the distance between them is ;

[0051] It is known that the moving speed of the steel plate on the production line is v (the unit is usually meters per second);

[0052] It is known that the moment when the steel plate cross-section reaches the first temperature measurement point S is ;

[0053] It is known that the moment when the steel plate cross-section reaches the second temperature measurement point E is ;

[0054] The set position is located between S and E, and the distance from the first temperature measurement point S is . Then it can be calculated through the following method :

[0055] Determine the distance between P and S: Assume that the distance from the set position to the first temperature measurement point S is ;

[0056] Calculate the time difference: The time required for the steel plate to move from S to P: ;

[0057] Calculate : .

[0058] In the present invention, based on the latest two-dimensional temperature field, the target three-dimensional temperature field is continuously updated to construct a dynamically evolving three-dimensional temperature field model. Specifically, in each process stage (including heating, rough rolling, finish rolling, and cooling), first, real-time temperature data is collected by temperature measurement points arranged on the steel plate production line and preprocessed to obtain the real-time temperature values. Then, based on the positions of the existing temperature measurement points and their corresponding real-time temperature values, the finite difference model is used to calculate the temperatures of each grid node on the cross-section, and the two-dimensional temperature field at the current moment is constructed. Next, based on the two-dimensional temperature field on the current cross-section and the spatial step corresponding to the current process stage, spatial expansion is carried out along the length direction of the steel plate, and the initial three-dimensional temperature field model is gradually constructed. According to the type of temperature measurement equipment corresponding to the temperature measurement points and the position information, the temperature correction values of each grid node in the initial three-dimensional temperature field model are calculated; and through the distance between any two temperature measurement devices, the temperature detection time, and the running speed of the steel plate, the moments when the steel plate cross-section reaches different temperature measurement points are obtained, and the temperature correction amount reflecting the time-distance dimension difference is calculated. The temperature of the grid node calculated based on the finite difference model is corrected by this temperature correction amount to obtain the target three-dimensional temperature field corresponding to the current moment. This process is continuously repeated, and each time it is updated based on the latest two-dimensional temperature field and the corresponding temperature correction values, ensuring that the three-dimensional temperature field model can reflect the actual temperature distribution characteristics of the steel plate in each production stage in real time, realizing dynamic evolution and accurate prediction.

[0059] The present invention discretizes the cross-sectional area of the steel plate and divides each cross-section into several rectangular grid units; based on the detection positions of the existing temperature measurement points and the real-time temperature values corresponding to each temperature measurement point, the finite difference model is used to calculate the temperatures of each grid node on the cross-section, and the two-dimensional temperature field corresponding to each cross-section is constructed based on this; based on the two-dimensional temperature field on the cross-section and the spatial step corresponding to the current process stage, spatial expansion is carried out along the length direction of the steel plate to construct the initial three-dimensional temperature field model; according to the type of temperature measurement equipment corresponding to the temperature measurement points and the positions of the temperature measurement points, the temperature correction values corresponding to each grid node on the steel plate cross-section in the initial three-dimensional temperature field model are calculated; the moments when the steel plate cross-section reaches the temperature measurement points corresponding to two temperature measurement devices during the steel plate production process are obtained, and the temperature correction amount reflecting the time-distance dimension difference is calculated based on these moments and the temperature correction values, and the temperature of the grid node calculated based on the finite difference model is corrected by the temperature correction amount; this method realizes the refined modeling of the three-dimensional temperature field of the steel plate, solves the problems that traditional one-dimensional or two-dimensional models are difficult to accurately reflect complex heat transfer behaviors, have large temperature measurement errors, and poor model adaptability, significantly improves the temperature field modeling accuracy and dynamic response ability, and provides reliable data support for subsequent process control and product quality improvement.

[0060] It should be noted that all directional indications (such as up, down, left, right, front, back...) in the embodiments of the present invention are only used to explain the relative positional relationship and movement conditions between components in a specific posture (as shown in the accompanying drawings). If the specific posture changes, the directional indications will also change accordingly.

[0061] In addition, in the present invention, descriptions such as "first", "second", "one", etc. are only for descriptive purposes and should not be construed as indicating or implying their relative importance or implicitly indicating the quantity of the indicated technical features. Thus, features defined with "first" and "second" may explicitly or implicitly include at least one such feature. In the description of the present invention, the meaning of "a plurality" is at least two, such as two, three, etc., unless otherwise specifically defined.

[0062] In the present invention, unless otherwise clearly specified and limited, the terms "connection", "fixation", etc. shall be understood in a broad sense. For example, "fixation" may be a fixed connection, a detachable connection, or integrated; it may be a mechanical connection or an electrical connection; it may be directly connected or indirectly connected through an intermediate medium, and it may be the internal communication of two components or the interaction relationship between two components, unless otherwise clearly limited. For those of ordinary skill in the art, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0063] In addition, the technical solutions between various embodiments of the present invention can be combined with each other, but it must be based on the ability of those of ordinary skill in the art to implement. When the combination of technical solutions results in contradictions or cannot be implemented, it should be considered that such a combination of technical solutions does not exist and is not within the protection scope required by the present invention.

Claims

1. A dynamic modeling method for three-dimensional temperature field in the steel plate production process, characterized in that Including: Collecting real-time temperature data of the steel plate in each process stage based on the temperature measurement points arranged on the steel plate production line; The process stages include: heating, rough rolling, finish rolling, and cooling; Preprocessing the real-time temperature data to obtain the real-time temperature value; Constructing a two-dimensional temperature field: discretizing the cross-sectional area of the steel plate and dividing each cross-section into several rectangular grid cells; the cross-section is composed of the width direction and the thickness direction of the steel plate; based on the detection positions of the existing temperature measurement points and the real-time temperature values corresponding to each temperature measurement point, using the finite difference model to calculate the temperatures of each grid node on the cross-section, and constructing the two-dimensional temperature field corresponding to each cross-section accordingly; Using the method of variable grid size to set the distance between adjacent cross-sections in each process stage as the space step in the length direction of the steel plate, and based on the two-dimensional temperature field on the cross-section and the space step corresponding to the current process stage, performing spatial expansion along the length direction of the steel plate to construct an initial three-dimensional temperature field model; Calculating the temperature correction value corresponding to each grid node on the cross-section of the steel plate in the initial three-dimensional temperature field model according to the type of temperature measurement equipment corresponding to the temperature measurement point and the position of the temperature measurement point; According to the distance between any two temperature measurement devices, the temperature detection time, and the running speed of the steel plate, obtaining the moments when the cross-section of the steel plate reaches the temperature measurement points corresponding to the two temperature measurement devices during the production process of the steel plate, calculating the temperature correction amount for reflecting the time-distance dimension difference according to these moments and the temperature correction value, and correcting the temperature of the grid node calculated based on the finite difference model through the temperature correction amount to obtain the target three-dimensional temperature field corresponding to the current moment; Based on the latest two-dimensional temperature field, continuously updating the target three-dimensional temperature field, thereby constructing a dynamically evolving three-dimensional temperature field model.

2. The three-dimensional temperature field dynamic modeling method for the steel plate production process according to claim 1, characterized in that, The specific construction of the two-dimensional temperature field is as follows: Obtaining the model parameters of the steel plate; Discretizing the cross-sectional area of the steel plate and dividing each cross-section into several rectangular grid cells; the cross-section is composed of the width direction and the thickness direction of the steel plate, and determining the space step ΔX in the width direction and the space step ΔY in the thickness direction; Setting the stability condition according to the space steps ΔX and ΔY and the material thermal property parameters of the steel plate, and setting the time step Δt based on the stability condition; Performing preliminary interpolation on the grid nodes using the detection positions of the existing temperature measurement points and the real-time temperature values corresponding to each temperature measurement point at the set moment; Based on the heat transfer control equation, using the model parameters of the steel plate and the real-time temperature value, and using the finite difference model to calculate the temperatures of each grid node on the cross-section at each moment under the set time step Δt; Calculating the temperatures of each point inside the grid cell using the interpolation method based on the temperatures of each grid node; Constructing the two-dimensional temperature field corresponding to each cross-section at each moment based on the temperatures of each grid node and the temperatures of each point inside it.

3. A three-dimensional temperature field dynamic modeling method for the steel plate production process according to claim 1, characterized in that The specific construction of the two-dimensional temperature field is as follows: Discretize the cross-section of the steel plate, divide each cross-section into several rectangular grid cells, and use the detection positions of the existing temperature measurement points and the real-time temperature values corresponding to each temperature measurement point at the set time to perform preliminary interpolation on the grid nodes, calculate the correction coefficients of each grid node in the width direction and thickness direction of the steel plate, adjust the real-time temperature values of the corresponding grid nodes through the correction coefficients, and calculate the temperatures of each grid node on the cross-section using the finite difference model based on the adjusted temperature data, and construct a two-dimensional temperature field on the cross-section with this.

4. A three-dimensional temperature field dynamic modeling method for the steel plate production process according to claim 3, characterized in that, The calculation of the correction coefficients of each grid node in the width direction and thickness direction of the steel plate, adjusting the real-time temperature values of the corresponding grid nodes through the correction coefficients, calculating the temperatures of each grid node on the cross-section using the finite difference model based on the adjusted temperature data, and constructing a two-dimensional temperature field on the cross-section with this is specifically as follows: Based on the temperature non-uniformity of the grid nodes in the width direction, construct a transverse correction function; Based on the temperature gradient of the grid nodes in the thickness direction, construct a vertical correction function; Use the transverse correction function and the vertical correction function to calculate the transverse correction coefficient of the grid nodes in the width direction of the steel plate and the vertical correction coefficient in the thickness direction of the steel plate; Correct the real-time temperature values of the grid nodes through the transverse correction coefficient and the vertical correction coefficient to obtain the temperature adjustment value; Based on the heat transfer control equation, use the model parameters of the steel plate and the temperature adjustment value, and calculate the temperatures of each grid node on the cross-section using the finite difference model; Based on the temperature adjustment value and the calculated temperatures of the grid nodes, use the interpolation algorithm to calculate the temperature values of the interpolation points inside the grid cells; Combine the temperature adjustment values of the grid nodes, the calculated temperatures of the grid nodes, and the temperature values of the internal interpolation points to construct a two-dimensional temperature field on the cross-section.

5. A dynamic three-dimensional temperature field modeling method for the steel plate production process according to any one of claims 2 to 4, characterized in that The model parameters of the steel plate include: preset geometric dimension information, material thermophysical properties parameters, and convective heat transfer coefficients under various boundary conditions; the boundary conditions include: Dirichlet boundary conditions and Neumann boundary conditions.

6. A three-dimensional temperature field dynamic modeling method for the steel plate production process according to claim 4, characterized in that The types of temperature measurement devices include single-point pyrometers and / or transverse scanning pyrometers; according to the type of temperature measurement device corresponding to the temperature measurement point and the position of the temperature measurement point, calculate the temperature correction values corresponding to each grid node on the cross-section of the steel plate in the initial model of the three-dimensional temperature field, specifically as follows: When there are only temperature measurement points arranged by single-point pyrometers on the upper surface of the steel plate, calculate the temperature correction values corresponding to each grid node in the width direction of the cross-section of the steel plate in the initial model of the three-dimensional temperature field through the first single-point temperature correction model; When there are temperature measurement points arranged by single-point pyrometers on both the upper surface and the lower surface of the steel plate, calculate the temperature correction values corresponding to each grid node in the thickness direction of the cross-section of the steel plate in the initial model of the three-dimensional temperature field through the second single-point temperature correction model; When there are temperature measurement points arranged by transverse scanning pyrometers on the steel plate surface, calculate the temperature correction values corresponding to each grid node in the width direction and thickness direction of the cross-section of the steel plate in the initial model of the three-dimensional temperature field through the transverse temperature correction model.

7. A dynamic three-dimensional temperature field modeling method for the steel plate production process according to claim 6, characterized in that The formula expression of the first single-point temperature correction model is: ; wherein: ； In the formula, represents the temperature correction amount at the center point of the upper surface of the steel plate, which is used to compensate for the deviation between the calculated value of the finite difference model and the measured temperature value; represents the grid node temperature of the grid node at the i-th column in the width direction and the j-th row in the thickness direction on the cross-section of the steel plate at the k-th moment calculated based on the finite difference model; represents the temperature correction value of the grid node at the i-th column in the width direction and the j-th row in the thickness direction on the cross-section of the steel plate at the k-th moment; represents the measured temperature value at the center point of the upper surface of the steel plate, where represents measurement; n represents the total number of grid nodes in the width direction of the cross-section of the steel plate, and m represents the total number of grid nodes in the thickness direction of the cross-section of the steel plate; represents the grid node temperature at the center point of the upper surface of the steel plate, i.e., the nth column and the mth row grid node at the kth moment, obtained by calculation using the finite difference model.

8. A three-dimensional temperature field dynamic modeling method for the steel plate production process according to claim 6, characterized in that, The formula expression of the second single-point temperature correction model is: ; wherein: ； ； In the formula, represents the temperature correction value of the grid node at the i-th column in the width direction and the j-th row in the thickness direction of the steel plate cross-section at the k-th moment; represents the grid node temperature of the grid node at the i-th column in the width direction and the j-th row in the thickness direction of the steel plate cross-section at the k-th moment calculated based on the finite difference model; represents the temperature correction amount at the center point of the upper surface of the steel plate; represents the temperature correction amount at the center point of the lower surface of the steel plate; represents the distance between adjacent grid nodes along the thickness direction of the steel plate cross-section; represents the coordinate position of the current grid node in the thickness direction; represents the total height of the steel plate, which is equal to the sum of the spacings of all grid nodes in the thickness direction; represents the measured temperature value at the center point of the upper surface of the steel plate, where represents measurement; represents the grid node temperature of the grid node at the center point of the upper surface of the steel plate, that is, the -th column and the m-th row, calculated by the finite difference model at the k-th moment; n represents the total number of grid nodes in the width direction of the steel plate cross-section, and m represents the total number of grid nodes in the thickness direction of the steel plate cross-section; represents the measured temperature value at the center point of the lower surface of the steel plate; represents the grid node temperature of the grid node at the center point of the lower surface of the steel plate, that is, the -th column and the 0-th row, calculated by the finite difference model at the k-th moment.

9. A three-dimensional temperature field dynamic modeling method for the steel plate production process according to claim 6, characterized in that, The formula expression of the horizontal temperature correction model is as follows: ; wherein: ； In the formula, represents the temperature correction value of the grid node at the i-th column in the width direction and the j-th row in the thickness direction of the steel plate cross-section at the k-th moment; represents the grid node temperature of the grid node at the i-th column in the width direction and the j-th row in the thickness direction of the steel plate cross-section at the k-th moment calculated based on the finite difference model; represents the overall lateral correction amount corresponding to the i-th column in the width direction; represents the measured temperature value of the grid node at the i-th column in the width direction and the topmost layer, i.e., the m-th row in the thickness direction, of the steel plate upper surface; represents the grid node temperature of the grid node at the i-th column in the width direction and the m-th row in the thickness direction of the steel plate cross-section at the k-th moment calculated based on the finite difference model.

10. A three-dimensional temperature field dynamic modeling method for the steel plate production process according to any one of claims 6 to 9, characterized in that, The calculation formula used to calculate the temperature correction amount reflecting the time-distance dimension difference based on these moments and the temperature correction value is as follows: ； Wherein, represents the moment when the cross-section of the steel plate reaches the first temperature measurement point S; represents the moment when the cross-section of the steel plate reaches the second temperature measurement point E; represents the moment when the cross-section of the steel plate reaches the set position between the two temperature measurement points ; represents the difference between the temperature correction value of the grid node at the position where the width direction is the i-th column and the thickness direction is the j-th row on the cross-section of the steel plate and the temperature of the grid node calculated by the corresponding finite difference model; represents the temperature correction amount corresponding to the grid node at the i-th column in the width direction and the j-th row in the thickness direction on the cross-section of the steel plate; The calculation formula used to correct the grid node temperature calculated based on the finite difference model by the temperature correction amount is as follows: ； In the formula, represents the grid node temperature of the grid node at the i-th column in the width direction and the j-th row in the thickness direction on the cross-section of the steel plate at the k-th moment calculated based on the finite difference model; represents the final temperature value after dimension correction of the grid node at the i-th column in the width direction and the j-th row in the thickness direction on the cross-section of the steel plate at the k-th moment.

Images