A method for predicting temperature field of stelmor air cooling line

Abstract

The present application belongs to the field of metallurgy, and provides a Stelmor air cooling line temperature field prediction method, which comprises the following steps: obtaining basic parameters, fan parameters, heat preservation cover parameters and roller parameters, calculating the wire rod position according to the basic parameters, fan parameters, heat preservation cover parameters and roller parameters, obtaining the cooling process information of the wire rod according to the wire rod position, obtaining the thermal physical property parameters of the wire rod according to the initial temperature, performing heat exchange model calculation according to the cooling process information and the thermal physical property parameters, obtaining the heat exchange coefficient, bringing the heat exchange coefficient into the difference equation for iterative calculation, and obtaining the temperature of the wire rod cross section at the next moment. The problems that the Stelmor air cooling line in the prior art can only detect the surface temperature of the rolled piece, cannot detect the core temperature of the rolled piece, is difficult to obtain the temperature difference of the whole cross section of the rolled piece, and cannot optimize the cooling process are solved. The cooling process is adjusted and optimized according to the temperature difference of the cross section of the rolled piece, and the product performance, wire passing property and qualified rate are improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of metallurgy, in particular to a Stelmor air cooling line temperature field prediction method. BACKGROUND

[0002] In the production process of threaded steel wire rod, the key to organization and performance control lies in the development of controlled rolling and controlled cooling process and post-rolling air cooling process, and temperature is one of the important factors affecting the rolling mill load distribution, rolling system, water cooling system, air cooling system and ensuring the quality of the final product in the production process of threaded steel wire rod. In recent years, with the digitalization and intelligentization transformation and upgrading of the steel industry, the threaded steel wire rod whole-process temperature prediction and control technology has developed rapidly. The phase transition process of high-speed wire rod from austenite to ferrite and pearlite mainly occurs in the Stelmor air cooling line, which has an important influence on the internal organization, mechanical properties and through-bar properties of the wire rod. Due to the different cooling speeds, there is usually a certain deviation in the temperature and organization performance of the lap joint points and non-lap joint points of the threaded steel rod. The existing inspection means in industrial production can only measure the surface temperature of the rolled piece, and cannot measure the core temperature of the rolled piece. Therefore, it is difficult to know the temperature difference of the whole cross section of the rolled piece in the industrial field, and it is of great significance to develop a system that can simulate the temperature evolution history of the cross section of the rolled piece in the Stelmor air cooling line to optimize the post-rolling cooling process and control the final organization and performance. SUMMARY

[0003] In view of the above-mentioned shortcomings of the prior art, the purpose of the present application is to provide a Stelmor air cooling line temperature field prediction method, which solves the technical problems in the prior art that the Stelmor air cooling line can only detect the surface temperature of the rolled piece, cannot detect the core temperature of the rolled piece, and cannot obtain the temperature difference of the whole cross section of the rolled piece, and cannot optimize the control of the cooling process.

[0004] To achieve the above-mentioned purposes and other related purposes, the present application provides a Stelmor air cooling line temperature field prediction method, which comprises:

[0005] obtaining basic parameters, fan parameters, heat preservation cover parameters and roller parameters, the basic parameters at least including steel grade, rolling product diameter and wire feeding temperature, the fan parameters at least including opening and closing state, air volume and opening degree of the fan, the heat preservation cover parameters at least including opening and closing state, position and length of the heat preservation cover, and the roller parameters at least including roller speed and length;

[0006] calculating the position of the rod according to the basic parameters, the fan parameters, the heat preservation cover parameters and the roller parameters;

[0007] obtaining cooling process information of the rod according to the position of the rod, the cooling process information at least including initial temperature and cooling schedule of the rod;

[0008] According to the initial temperature, a thermal physical parameter of the wire rod is obtained, the thermal physical parameter at least including a thermal conductivity coefficient, a density and a specific heat capacity of the wire rod;

[0009] According to the cooling process information and the thermal physical parameter, a heat exchange model is calculated to obtain a heat exchange coefficient;

[0010] The heat exchange coefficient is brought into a difference equation to perform iterative calculation to obtain a temperature of a wire rod section at a next time.

[0011] Optionally, the heat exchange model includes an air cooling heat exchange model and a wind cooling heat exchange model, the air cooling heat exchange model is calculated when the fan is in a closed state, and the wind cooling heat exchange model is calculated when the fan is in an opened state.

[0012] Optionally, the air cooling heat exchange model is:

[0013] H a =εk B ((T+273.15) 2 +(T f +273.15) 2 )((T+273.15)+(T f +273.15))+(T-T f ) 0.25 ;

[0014] In the formula, H a is an air cooling heat exchange coefficient, ε is a thermal emissivity, k B is a Boltzmann constant, T is a rolling piece surface temperature, and T f is an environmental temperature.

[0015] Optionally, the wind cooling heat exchange model is:

[0016] H w =εk B ((T+273.15) 2 +(T f +273.15) 2 )((T+273.15)+( T f+273.15))+λ a *N u / D;

[0017] In the formula, H w is a wind cooling heat exchange coefficient, ε is a thermal emissivity, k B is a Boltzmann constant, T is a rolling piece surface temperature, T f is an environmental temperature, λ a is an air thermal conductivity coefficient, and N uNu is the Nusselt number; D is the diameter of the rod.

[0018] Optionally, when calculating the non-lap joint point of the rod, the Nusselt number calculation model is:

[0019] N u = C * (G r * P r ) n ;

[0020] In the formula, G r is the Grashof number, P r is the Prandtl number, C is a constant, the value range is 0.0165-0.48, and n is an exponential coefficient, the value range is 0.25-0.42.

[0021] Optionally, the Grashof number calculation model is:

[0022]

[0023] In the formula, g is the acceleration of gravity, a a is the air thermal expansion coefficient, T is the rod surface temperature, T f is the air temperature, D is the rod diameter, and v a is the air kinematic viscosity.

[0024] Optionally, when calculating the lap joint point of the rod, the Nusselt number calculation model is:

[0025]

[0026] In the formula, R ef is the Reynolds number at the interface temperature; P rf is the Prandtl number at the interface temperature; P r is the Prandtl number at the air temperature; and n is a constant, the value range is 0.4-0.6.

[0027] Optionally, the Prandtl number calculation model is:

[0028]

[0029] In the formula, v a is the air kinematic viscosity, and a is the air thermal diffusion coefficient.

[0030] The Reynolds number calculation model is:

[0031] R e = u w *D / v a ;

[0032] In the formula, u w is the wind speed, D is the rod diameter, and v aFor air kinematic viscosity.

[0033] Optionally, the wind speed calculation model is:

[0034]

[0035] In the formula, Q is the rated air volume of the fan, P is the rated power of the fan, K is the fan opening degree, and n is a coefficient, and the value range is 10-20.

[0036] Optionally, the initial condition of the wire rod is obtained according to the basic parameters.

[0037] The boundary condition of the wire rod is calculated according to the thermal physical parameters and the cooling process information.

[0038] The difference equation is obtained through the initial condition and the boundary condition.

[0039] As described above, the present application provides a Stelmor air cooling line temperature field prediction method, which has at least the following beneficial effects:

[0040] By obtaining basic parameters, fan parameters, heat preservation cover parameters and roller parameters, the basic parameters at least include steel grade, rolled product diameter and wire drawing temperature, the wire rod position is calculated according to the basic parameters, fan parameters, heat preservation cover parameters and roller parameters, the cooling process information of the wire rod is obtained according to the wire rod position, the thermal physical parameters of the wire rod are obtained according to the initial temperature, the heat exchange coefficient is obtained through the heat exchange model calculation according to the cooling process information and the thermal physical parameters, the heat exchange coefficient is brought into the difference equation for iterative calculation to obtain the temperature of the wire rod cross section at the next moment, the temperature field prediction of the Stelmor air cooling line is realized, the problem that the Stelmor air cooling line can only detect the surface temperature of the rolled piece in the prior art, cannot detect the core temperature of the rolled piece, is difficult to obtain the temperature difference of the whole cross section of the rolled piece, and cannot control and optimize the cooling process is solved, the temperature evolution history of the wire rod at the joint point and the non-joint point in the Stelmor air cooling line interval can be completely predicted, the adjustment and optimization of the high-speed wire rod cooling process after rolling of the deformed steel bar is carried out according to the temperature difference of the rolled piece cross section, and the product performance, the wire passing property and the qualified rate are improved.

[0041] It should be understood that the foregoing general description and the following detailed description are only exemplary and explanatory, and cannot limit the present application. BRIEF DESCRIPTION OF DRAWINGS

[0042] The drawings herein are incorporated into the specification and form part of the specification, show embodiments consistent with the present application, and together with the specification serve to explain the principles of the present application. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor. In the drawings:

[0043] Figure 1 is a flow chart of a Stelmor air cooling line temperature field prediction method according to an example embodiment of the present application. DETAILED DESCRIPTION

[0044] Other advantages and effects of the present application can be easily understood by those skilled in the art from the above description of the embodiments of the present application. The present application can also be implemented or applied in other different specific embodiments, and various modifications or changes can be made to the details of the present application based on different views and applications without departing from the spirit of the present application. It should be understood that the preferred embodiments are only for illustrating the present application, but not for limiting the protection scope of the present application.

[0045] It should be noted that the diagrams provided in the embodiments only schematically illustrate the basic concepts of the present application, and only the components related to the present application are shown in the diagrams, but not the number, shape and size of the components when actually implemented. The shapes, number and proportions of the components when actually implemented can be arbitrarily changed, and the layout pattern of the components can also be more complex. The structures, proportions, sizes, etc. shown in the diagrams attached to the present specification are only used to cooperate with the content disclosed in the specification, so as to be understood and read by those skilled in the art, and do not define the limiting conditions for implementing the present application, and therefore do not have technical significance. Any modification of the structure, change of the proportional relationship or adjustment of the size, without affecting the effects and purposes that can be achieved by the present application, should still fall within the scope of the technical content disclosed by the present application. Meanwhile, the terms such as "upper", "lower", "left", "right", "middle" and "one" used in the present specification are only for the convenience of clear understanding of the description, and not for limiting the scope of implementation of the present application. The change or adjustment of the relative relationship, without substantially changing the technical content, is also considered as the implementation scope of the present application.

[0046] In the following description, a large number of details are discussed to provide a more thorough explanation of the embodiments of the present application, however, it is obvious to those skilled in the art that the embodiments of the present application can be implemented without these specific details, and in other embodiments, the well-known structures and devices are shown in the form of block diagrams instead of details, so as to avoid making the embodiments of the present application difficult to understand.

[0047] In an example embodiment, the present application exemplarily provides a Stelmor air cooling line temperature field prediction method, please refer to Figure 1 , Figure 1 is a flow chart of a Stelmor air cooling line temperature field prediction method according to an example embodiment of the present application, which comprises at least steps S110 to S160, which are described in detail as follows:

[0048] In step S110, the basic parameters, the fan parameters, the heat preservation cover parameters and the roller parameters are obtained, the basic parameters at least including the steel grade, the rolled product diameter and the wire bar temperature, the fan parameters at least including the opening and closing state, the air volume and the opening degree of the fan, the heat preservation cover parameters at least including the opening and closing state, the position and the length of the heat preservation cover, and the roller parameters at least including the roller speed and the length.

[0049] In step S120, the wire rod position is calculated according to the basic parameters, the fan parameters, the heat preservation cover parameters and the roller parameters.

[0050] In step S130, the cooling process information of the wire rod is obtained according to the wire rod position, the cooling process information at least including the initial temperature and the cooling schedule of the wire rod.

[0051] In step S140, the thermal physical parameters of the wire rod are obtained according to the initial temperature, the thermal physical parameters at least including the thermal conductivity, the density and the specific heat capacity of the wire rod.

[0052] In step S150, the heat exchange model calculation is performed according to the cooling process information and the thermal physical parameters, and the heat exchange coefficient is obtained.

[0053] In step S160, the heat exchange coefficient is brought into the difference equation for iterative calculation, and the temperature of the wire rod cross section at the next time is obtained.

[0054] The steps S110 to S160 are described in detail as follows:

[0055] In step S110, the basic parameters, the fan parameters, the heat preservation cover parameters and the roller parameters are obtained, the basic parameters at least including the steel grade, the rolled product diameter and the wire bar temperature, the fan parameters at least including the opening and closing state, the air volume and the opening degree of the fan, the heat preservation cover parameters at least including the opening and closing state, the position and the length of the heat preservation cover, and the roller parameters at least including the roller speed and the length, and the obtained basic parameters, the fan parameters, the heat preservation cover parameters and the roller parameters are input into the system front end, the system front end can be programmed by JavaScript language, mainly responsible for data input and result display, wherein the data input can be automatically input by the system or input by the operator, which is not limited herein.

[0056] In step S120, the wire rod position is calculated according to the basic parameters, the fan parameters, the heat preservation cover parameters and the roller parameters, the basic parameters, the fan parameters, the heat preservation cover parameters and the roller parameters received by the system front end are processed by the data processing back end and sent to the calculation algorithm, the data processing back end can be programmed by Java language, mainly responsible for the communication and data transmission between the system front end and the calculation algorithm, processing, and managing and maintaining the business process.

[0057] In step S130, the calculation algorithm calculates the coil position according to the transmitted information, matches the cooling process information according to the coil position, and the algorithm system can be written in Python language, mainly responsible for the calculation of the Stelmor air cooling line temperature field, taking the heat exchange model based on heat transfer and fluid mechanics and combined with the actual temperature of the high-speed wire rod workshop as the core, using the finite difference method to simulate the temperature evolution process of the deformed steel high-speed wire rod on the Stelmor line.

[0058] In step S140, the calculation algorithm obtains the thermal physical parameters of the coil rod according to the initial temperature, and the thermal physical parameters at least include the thermal conductivity, density and specific heat capacity of the coil rod.

[0059] In step S150, the calculation algorithm calculates the heat exchange coefficient according to the cooling process information and the thermal physical parameters, wherein the heat exchange model includes an air cooling heat exchange model and an air cooling heat exchange model, when the fan is in the closed state, the air cooling heat exchange model is calculated, and when the fan is in the open state, the air cooling heat exchange model is calculated.

[0060] The air cooling heat exchange model is:

[0061] H a =εk B ((T+273.15) 2 +(T f +273.15) 2 )((T+273.15)+(T f +273.15))+(T-T f ) 0.25 ;

[0062] In the formula, H a is the air cooling heat exchange coefficient, ε is the thermal radiation rate, the value range is 0.7-0.9, k B is the Boltzmann constant, T is the surface temperature of the rolled piece; T f is the ambient temperature.

[0063] The air cooling heat exchange model is:

[0064] H w =εk B ((T+273.15) 2 +(T f +273.15) 2 )((T+273.15)+(T f +273.15))+λ a *N u / D;

[0065] In the formula, H wε is the heat radiation rate, and its value range is 0.7-0.9, k B is the Boltzmann constant, T is the surface temperature of the rolled piece, T f is the ambient temperature, λ a is the air thermal conductivity, N u is the Nusselt number; D is the wire diameter.

[0066] When the non-lap joint of the wire is calculated, the Nusselt number calculation model is:

[0067] N u =C*(G r *P r ) n ;

[0068] In the formula, G r is the Grashof number, P r is the Prandtl number, C is a constant, and its value range is 0.0165-0.48, and n is an index coefficient, and its value range is 0.25-0.42.

[0069] Wherein, the Grashof number calculation model is:

[0070]

[0071] In the formula, g is the acceleration of gravity, α a is the air thermal expansion coefficient, T is the surface temperature of the wire, T f is the air temperature, D is the wire diameter, and v a is the air kinematic viscosity.

[0072] When the lap joint of the wire is calculated, the Nusselt number calculation model is:

[0073]

[0074] In the formula, R ef is the Reynolds number at the interface temperature; P rf is the Prandtl number at the interface temperature; P r is the Prandtl number at the air temperature; and n is a constant, and its value range is 0.4-0.6.

[0075] Wherein, the Prandtl number calculation model is:

[0076]

[0077] In the formula, v a is the air kinematic viscosity, and a is the air thermal diffusion coefficient;

[0078] The Reynolds number calculation model is:

[0079] R e =uw D / v a ;

[0080] wherein u w is the wind speed, D is the diameter of the rod, v a is the kinematic viscosity of air.

[0081] wherein the wind speed calculation model is:

[0082]

[0083] wherein Q is the rated air volume of the fan, P is the rated power of the fan, K is the opening degree of the fan, and n is a coefficient, and the value range is 10-20.

[0084] The system selects the corresponding heat exchange model for calculation according to the opening and closing state of the fan. When the fan is in the closed state, the air cooling heat exchange model is called for calculation. When the fan is in the open state, the air cooling heat exchange model is called for calculation. Since the cooling speeds of the lap joint points and the non-lap joint points on the rod are different, there is a certain deviation in the temperature and the organizational performance. Different Nusselt number calculation models are called for the lap joint points and the non-lap joint points, which can completely predict the temperature evolution history of the lap joint points and the non-lap joint points of the rod in the Stelmor air cooling line section, so as to optimize the post-rolling cooling process parameters, regulate and control the organizational performance of the final rolled product, and further reduce the organizational performance difference between the lap joint points and the non-lap joint points.

[0085] Table 1 is the comparison result of the predicted value and the measured value of the lap joint point temperature of the rod in the embodiment. The predicted value and the actual value are both the surface temperature of the blank.

[0086] Table 1 is the comparison result of the predicted value and the measured value of the lap joint point temperature of the rod in the embodiment. The predicted value and the actual value are both the surface temperature of the blank.

[0087] Table 1 is the comparison result of the predicted value and the measured value of the lap joint point temperature of the rod in the embodiment. The predicted value and the actual value are both the surface temperature of the blank.

[0088]

[0089] In step S160, the heat exchange coefficient is brought into the difference equation for iterative calculation to obtain the temperature of the wire rod section at the next moment. The calculation algorithm solves the high-speed wire rod section temperature at the next moment according to the heat exchange coefficient and the difference equation. According to the calculated wire rod section temperature at the next moment, the thermal physical property parameters and the cooling process information are reacquired, the boundary conditions of the wire rod are calculated, the updated initial conditions and boundary conditions are combined, the calculation heat exchange coefficient is obtained through the air cooling heat exchange model and the air cooling heat exchange model, the heat exchange coefficient is brought into the difference equation to obtain the temperature value of the wire rod section at the next moment, and the iteration is repeated in this way until the complete Stelmor air cooling line temperature field data are obtained. The final obtained temperature field data are stored in the database. The data processing backend obtains the data from the database according to the rules and then transmits the data to the system front end for result display. The final display result includes the temperature evolution curve and the rolling piece section temperature value. The temperature evolution curve includes the lap joint temperature-distance curve and the non-lap joint temperature-distance curve, which facilitates subsequent data analysis.

[0090] It can be seen that the technical scheme provided by the embodiment couples the heat exchange model by using the difference idea, builds a Stelmor air cooling line temperature field calculation model of high-speed wire rod, simulates the whole section temperature evolution process from the wire laying machine to the center of the coiling machine, acquires the basic parameters, fan parameters, heat preservation cover parameters and roller parameters, the basic parameters at least including the steel grade, the rolling product diameter and the wire laying temperature, the fan parameters at least including the opening and closing state, the air volume and the opening degree of the fan, the heat preservation cover parameters at least including the opening and closing state, the position and the length of the heat preservation cover, and the roller parameters at least including the roller speed and the length, calculates the wire rod position according to the basic parameters, the fan parameters, the heat preservation cover parameters and the roller parameters, obtains the cooling process information of the wire rod according to the wire rod position, the cooling process information at least including the initial temperature and the cooling schedule of the wire rod, obtains the thermal physical property parameters of the wire rod according to the initial temperature, the thermal physical property parameters at least including the thermal conductivity, the density and the specific heat capacity of the wire rod, performs heat exchange model calculation according to the cooling process information and the thermal physical property parameters to obtain the heat exchange coefficient, brings the heat exchange coefficient into the difference equation for iterative calculation to obtain the temperature of the wire rod section at the next moment, realizes the temperature field prediction of the Stelmor air cooling line, solves the problem that the Stelmor air cooling line can only detect the surface temperature of the rolling piece in the prior art, cannot detect the core temperature of the rolling piece, is difficult to obtain the temperature difference of the whole rolling piece section, and cannot optimize the cooling process, can completely predict the temperature evolution process of the wire rod at the lap joint and the non-lap joint in the Stelmor air cooling line interval, adjusts and optimizes the post-rolling cooling process of the high-speed wire rod of the screw steel according to the temperature difference of the rolling piece section, and improves the product performance, the wire passing property and the qualified rate.

[0091] The above embodiments merely illustrate the principles of the present application and its efficacy, and are not intended to limit the present application. Any modification or change made by those skilled in the art without departing from the spirit and scope of the present application shall be covered by the claims of the present application.

Claims

1. A method of predicting the temperature field of a Stelmor cooling line, characterized by, The method comprises: acquiring basic parameters, fan parameters, heat preservation cover parameters and roller table parameters, the basic parameters at least including steel grade, rolled product diameter and wire rod temperature, the fan parameters at least including fan opening and closing state, air volume and opening degree, the heat preservation cover parameters at least including heat preservation cover opening and closing state, position and length, and the roller table parameters at least including roller table speed and length; calculating wire rod position according to the basic parameters, the fan parameters, the heat preservation cover parameters and the roller table parameters; obtaining wire rod cooling process information according to the wire rod position, the cooling process information at least including wire rod initial temperature and cooling schedule; obtaining wire rod thermal physical property parameters according to the initial temperature, the thermal physical property parameters at least including wire rod thermal conductivity, density and specific heat capacity; carrying out heat exchange model calculation according to the cooling process information and the thermal physical property parameters to obtain heat exchange coefficient; bringing the heat exchange coefficient into difference equation to carry out iterative calculation to obtain wire rod section temperature at next time.

2. The method of claim 1, wherein, The step of carrying out heat exchange model calculation according to the cooling process information and the thermal physical property parameters to obtain heat exchange coefficient further comprises: The heat exchange model comprises air cooling heat exchange model and air cooling heat exchange model, when the fan is in the closed state, the air cooling heat exchange model is calculated, and when the fan is in the open state, the air cooling heat exchange model is calculated.

3. The method for predicting the temperature field of the Stellmore wind-cooled line according to claim 2, characterized in that, In the step of carrying out heat exchange model calculation according to the cooling process information and the thermal physical property parameters to obtain heat exchange coefficient, The air cooling heat exchange model is: H a = εk B ((T+273.15) 2 +(T f +273.15) 2 )((T+273.15)+(T f +273.15))+(T-T f ) 0.25 ; In the formula, H a is the heat exchange coefficient of air cooling, ε is the thermal emissivity, the value range is 0.7-0.9, k B is the Boltzmann constant, T is the surface temperature of the rolled piece; T f is the ambient temperature.

4. The method for predicting the temperature field of the Stellmore wind-cooled line according to claim 2, characterized in that, In the step of carrying out heat exchange model calculation according to the cooling process information and the thermal physical property parameters to obtain heat exchange coefficient, The air cooling heat exchange model is: H w = εk B ((T+273.15) 2 +(T f +273.15) 2 )((T+273.15)+(T f +273.15))+λ a *N u / D; where H w is the heat transfer coefficient of air cooling, ε is the thermal emissivity, and its value range is 0.7-0.9, k B is the Boltzmann constant, T is the surface temperature of the rolled piece, T f is the ambient temperature, λ a is the air thermal conductivity, N u is the Nusselt number; D is the wire rod diameter.

5. The method for predicting the temperature field of the Stellmore wind-cooled line according to claim 4, characterized in that, In the step of carrying out heat exchange model calculation according to the cooling process information and the thermal physical property parameters to obtain heat exchange coefficient, When the wire rod non-overlapping point is calculated, the Nusselt number calculation model is: N u = C * (G r * P r ) n ; In the formula, G r is the Grashof number, P r is the Prandtl number, C is a constant, the value range is 0.0165-0.48, and n is an index coefficient, the value range is 0.25-0.

42.

6. The method of claim 5, wherein the temperature field is predicted by solving the following equation: ###0001### where T is the temperature field, V is the velocity field, k is the thermal conductivity, and p is the density of the Stirling air-cooled generator. 5 In the step of carrying out heat exchange model calculation according to the cooling process information and the thermal physical property parameters to obtain heat exchange coefficient, The Grashof number calculation model is: where g is the acceleration due to gravity, a a is the air expansion coefficient, T is the wire rod surface temperature, T f is the air temperature, D is the wire rod diameter, v a is the air kinematic viscosity.

7. The method for predicting the temperature field of the Stellmore wind-cooled line according to claim 6, characterized in that, In the step of carrying out heat exchange model calculation according to the cooling process information and the thermal physical property parameters to obtain heat exchange coefficient, When the wire rod overlapping point is calculated, the Nusselt number calculation model is: wherein R ef is the Reynolds number at the interface temperature; P rf is the Prandtl number at the interface temperature; P r is the Prandtl number at the air temperature; n is a constant having a value in the range of 0.4-0.

6.

8. The method for predicting the temperature field of the Stellmore wind-cooled line according to claim 7, characterized in that, In the step of carrying out heat exchange model calculation according to the cooling process information and the thermal physical property parameters to obtain heat exchange coefficient, The Prandtl number calculation model is: where v a is the kinematic viscosity of air, a is the thermal diffusivity of air; The Reynolds number calculation model is: R e = u w D / v a ; where u is the wind speed, D is the diameter of the rod, v is the air kinematic viscosity. w a where u is the wind speed, D is the diameter of the rod, v is the air kinematic viscosity.​ 9. The method for predicting the temperature field of the Stellmore wind-cooled line according to claim 8, characterized in that, In the heat exchange model includes air cooling heat exchange model and air cooling heat exchange model, when the fan is in the closed state, the air cooling heat exchange model is calculated, and when the fan is in the open state, the air cooling heat exchange model is calculated. The wind speed calculation model is: In the formula, Q is the rated air volume of the fan, P is the rated power of the fan, K is the opening degree of the fan, n is a coefficient, and the value range is 10-20.

10. The method of claim 1, wherein the temperature field is predicted for a Stelmor cooling line. The step of bringing the heat exchange coefficient into the difference equation for iterative calculation to obtain the temperature of the wire rod section at the next time point comprises: According to the basic parameters, the initial condition of the wire rod is obtained. According to the thermal physical parameters and the cooling process information, the boundary condition of the wire rod is calculated. The difference equation is obtained through the initial condition and the boundary condition.

Images