Simulation calculation method for roller moving induction heating temperature field

Abstract

The invention discloses a simulation calculation method for a roller movement induction heating temperature field. The simulation calculation method comprises the following steps: (1) establishing a finite element simulation model of the roller movement induction heating temperature field; (2) setting finite element simulation thermophysical parameters of a roller moving induction heating temperature field; (3) importing a result of a finite element simulation electromagnetic module of a roller movement induction heating temperature field; (4) setting finite element simulation thermal radiation and thermal convection parameters of a roller moving induction heating temperature field; and (5) calculating a roller movement induction heating temperature field simulation thermal module. According to the method, the overall time is decomposed into continuous time points, the surface area of the roller is decomposed into different areas, the roller at different moments and in different areas is subjected to thermal parameter processing endowing, the thermal physical property parameter in the highest temperature time interval of the decomposed temperature measuring point is used as a reference value, and the thermal parameter of the temperature measuring point is calculated according to the moving speed and the temperature translation principle. And three-phase coupling simulation of the thermal physical parameter transient magnetic thermal velocity in the temperature field in the actual roller moving induction heating process is realized.

Description

Technical Field

[0001] The invention belongs to the technical field of temperature field simulation calculation, and in particular relates to a simulation calculation method for the temperature field of mobile induction heating of a rolling roller. Background Art

[0002] With the development of industry, the requirements for rolled steel are becoming increasingly stringent, and therefore the quality requirements for rolls are also gradually increasing. In order to optimize the relevant parameters in the mobile induction heating of rolls, save costs, guide the actual roll induction heating process, and improve roll quality, it is necessary to conduct finite element simulation analysis of mobile induction heating of rolls.

[0003] Due to the continuous relative motion between the rollers and the induction coils, the temperature field exhibits dynamic, non-steady-state characteristics. This involves complex factors such as magnetic, thermal, and multi-field coupling, nonlinear material changes, and dynamic boundary conditions. Traditional modeling often uses static assumptions or simplified boundary conditions, resulting in large prediction errors, making it difficult to support real-time control. The computational overhead is also long, failing to meet the real-time demands of industry. Existing debugging relies on manual trial and error, which is inefficient and increases energy consumption and roller wear.

[0004] Therefore, there is an urgent need to build an efficient dynamic temperature field model and a rapid debugging method to resolve the contradiction between accuracy and efficiency and realize the optimal control of the process during the mobile induction heating of the rolling mill rollers. This is a technical problem that needs to be solved urgently in this field. Summary of the Invention

[0005] The purpose of the present invention is to solve the above problems and provide a simulation calculation method for the temperature field of roller movement induction heating to fill the gap in this field.

[0006] The technical solution for achieving the purpose of the present invention is: a simulation calculation method for the temperature field of roller movement induction heating, comprising the following steps: ① Establishment of finite element simulation model of temperature field of roller moving induction heating; ② Setting of thermophysical parameters for finite element simulation of temperature field of roller moving induction heating; ③ Import the results of electromagnetic module of finite element simulation of temperature field of roller moving induction heating; ④ Finite element simulation of heat radiation and heat convection settings for roller moving induction heating temperature field; ⑤ Calculation of thermal module for roller moving induction heating temperature field simulation.

[0007] The specific method of step ① above is: combining the structural characteristics of the roller and the heating process, simplifying the model into a two-dimensional axisymmetric model, and performing block processing on the surface of the roller, performing encrypted division near the temperature measurement point, selecting the solver type as the transient thermal solver, importing the model into the geometry of the transient thermal module in ANSYS, and meshing it.

[0008] In the above step ②, the thermophysical parameters include density, thermal conductivity and specific heat capacity.

[0009] The specific method of step ③ above is: import all the calculation results in the Maxwell module into the transient thermal module, divide the analysis steps in the analysis step settings, keep the analysis time of each analysis step corresponding to the final value of the time interval in Maxwell, and call the imported ohmic loss value in the corresponding Maxwell time interval in the time interval of each analysis step in the thermal module.

[0010] The specific setting described in the above step ④ is: setting the outer contour of the roller body as its heat radiation surface and convection heat transfer surface at different times and in different areas, and setting the radiation coefficient and convection heat transfer coefficient.

[0011] The above step ⑤ calculates and obtains the temperature rise curve of the roll under the proposed process parameter conditions and the temperature field distribution and changes during the heating process.

[0012] The present invention has the positive effects: (1) Based on the actual structure of the factory roller and the mobile induction heater, the present invention rationally simplifies them into a two-dimensional axisymmetric model to simulate the entire mobile heating process of the roller. The model is simple and the number of grids is small, which greatly saves calculation time and improves calculation efficiency.

[0013] (2) The present invention decomposes the overall time into continuous time points and decomposes the surface area of ​​the roller into different areas, and processes and assigns thermal parameters to the rollers in different areas at different times. The thermal physical parameters in the highest temperature time interval of the temperature measuring point after decomposition are used as the reference value. According to the moving speed and the temperature translation principle, the transient magnetic thermal velocity three-phase coupling simulation of the thermal physical parameters in the temperature field of the actual roller moving induction heating process is realized. BRIEF DESCRIPTION OF THE DRAWINGS

[0014] Figure 1 The figure is a flow chart of the simulation calculation method of the temperature field of the mobile induction heating of the rolling mill roll according to the present invention.

[0015] Figure 2 A two-dimensional axisymmetric finite element model of the roller established for the present invention.

[0016] Figure 3 Schematic diagram of mesh division of the finite element model of the present invention.

[0017] Figure 4 This is a diagram of the roll density setting in the transient thermal simulation of the present invention.

[0018] Figure 5 This is a diagram showing the thermal conductivity settings of the rolls in the transient thermal simulation of the present invention.

[0019] Figure 6The setting diagram for specific heat capacity of the roll in transient thermal simulation of the application.

[0020] Figure 7 The import diagram for the calculation result of the application.

[0021] Figure 8 The setting diagram for solving analysis step in transient thermal simulation of the application.

[0022] Figure 9 The import diagram for heat loss in transient thermal simulation of the application.

[0023] Figure 10 The setting for thermal radiation and thermal convection of the application.

[0024] Figure 11 The setting for convective heat transfer surface of the application.

[0025] Figure 12 The setting for thermal radiation surface of the application.

[0026] Figure 13 The highest temperature point convective heat transfer coefficient curve of the application.

[0027] Figure 14 The temperature cloud chart of the roll at 3460s of the application.

[0028] Figure 15 The temperature cloud chart of the roll at 4880s of the application.

[0029] Figure 16 The temperature cloud chart of the roll at 6700s of the application.

[0030] Figure 17 The comparison chart of the measured temperature rise and the simulation calculation temperature rise of the roll at the radial direction of 12mm of the application.

[0031] Figure 18 The comparison chart of the measured temperature rise and the simulation calculation temperature rise of the roll at the radial direction of 20mm of the application.

[0032] Figure 19 The comparison chart of the measured temperature rise and the simulation calculation temperature rise of the roll at the radial direction of 30mm of the application.

[0033] Figure 20 The comparison chart of the measured temperature rise and the simulation calculation temperature rise of the roll at the radial direction of 40mm of the application. DETAILED DESCRIPTION

[0034] (Example 1) This example takes the GP850x175mm / 200mm type inductor heating 20# steel material Ф800 roll as an example, and combines the attached drawings to specifically illustrate the simulation calculation method of the temperature field of the roll moving inductive heating. Figure 1 , the simulation calculation method of the temperature field of the roll moving inductive heating.

[0035] ① Establishment of finite element simulation model of temperature field of roller moving induction heating.

[0036] Combined with the structural characteristics of the roller and the heating process, a two-dimensional axisymmetric finite element model of the roller was established in SolidWorks. The surface of the roller was divided into blocks, and the area near the temperature measurement point was divided into denser blocks. Figure 2 Since the temperature measurement point is at the center of the roller, the block width at the center is set to 40mm, and increases to 80mm, 100mm, and 150mm at the edge.

[0037] Select the transient thermal solver as the solver type, import the model into the geometry of the transient thermal module in ANSYS, and mesh it. Figure 3 , the module with a width of 40mm is divided into 2mm grids, the modules with other widths are divided into 8mm grids, the remaining main body of the roller is divided into 20mm grids, and the rest are divided into automatic grids according to the system.

[0038] ② Setting of thermophysical parameters for finite element simulation of temperature field of rolling roller moving induction heating.

[0039] Set the density, thermal conductivity and specific heat capacity of the roller respectively. Figures 4-6 .

[0040] ③ Import the results of the electromagnetic module of the finite element simulation of the temperature field of the moving induction heating of the rolling mill roller.

[0041] Import all calculation results from Maxwell modules into transient thermal module, see Figure 7 .

[0042] The analysis steps are divided in the analysis step settings. The analysis time of each analysis step corresponds to the final value of the time interval in Maxwell. Figure 8 .

[0043] In each analysis step of the thermal module, the ohmic loss value in the corresponding Maxwell time interval is retrieved. Figure 9 .

[0044] ④ Finite element simulation of heat radiation and heat convection settings of the temperature field of rolling mill moving induction heating.

[0045] See also Figures 10 to 13 , set the heat radiation surface, radiation coefficient, convection heat transfer surface and convection heat transfer coefficient of the roller.

[0046] During the entire heating process, since the convection heat transfer coefficient between the roller body surface and the air can be regarded as basically unchanged, the convection heat transfer coefficient of the roller body surface is set to the same value.

[0047] In the heat radiation setting, since its coefficient is greatly affected by the movement of the roll and the temperature measuring point is at the center of the roll, the roll surface is divided into 7 areas for setting, and the width of the surface area of ​​the roll is the smallest in the center and the largest at both ends.

[0048] ⑤ Calculation of thermal module for roller moving induction heating temperature field simulation.

[0049] After all settings are completed, the roller transient thermal module is calculated to calculate the temperature rise of the roller. The temperature cloud diagrams at different times (3460s, 4880s, 6700s) are shown in the figure. Figures 14-16 .

[0050] After calculating the temperature results through simulation, the accuracy of the simulation results is verified through actual temperature measurement experiments. The comparison chart of the measured temperature rise and the simulated temperature rise at different radial depths (12mm, 20mm, 30mm, 40mm) of the roller center point is shown in Figures 17-20 .

[0051] The results show that at the end of heating, the temperature difference between the calculated result and the actual maximum temperature result is within ±35℃. Combined with the temperature curve of the roller, the maximum temperature error range is 0-3.36%, and the minimum temperature error range is 0-7.09%. Since the indicator during the heating process is the highest temperature, the accuracy of the entire simulation process reaches more than 96.64%, indicating that the simulation calculation results are basically consistent with the experimental measured results. The simulation model can be used to simulate other heating parameters to optimize the relevant heating parameters.

Claims

1. A simulation calculation method for the temperature field of roller moving induction heating, characterized in that The following steps are involved: ① Establishment of finite element simulation model of temperature field of roller moving induction heating; ② Setting of thermophysical parameters for finite element simulation of temperature field of roller moving induction heating; ③ Import the results of electromagnetic module of finite element simulation of temperature field of roller moving induction heating; ④ Setting of thermal radiation and convection parameters for finite element simulation of temperature field of roller moving induction heating; ⑤ Calculation of thermal module for roller moving induction heating temperature field simulation.

2. The simulation calculation method of the roller moving induction heating temperature field according to claim 1 is characterized in that The specific method of step ① above is: combining the structural characteristics of the roller and the heating process, simplifying the model into a two-dimensional axisymmetric model, and performing block processing on the surface of the roller, performing encrypted division near the temperature measurement point, selecting the solver type as the transient thermal solver, importing the model into the geometry of the transient thermal module in ANSYS, and meshing it.

3. The simulation calculation method for the roller moving induction heating temperature field according to claim 1 is characterized in that: In the above step ②, the thermophysical parameters include density, thermal conductivity and specific heat capacity.

4. The method for simulating and calculating the temperature field of roller moving induction heating according to claim 1 is characterized in that: The specific method of step ③ above is: import all the calculation results in the Maxwell module into the transient thermal module, divide the analysis steps in the analysis step settings, keep the analysis time of each analysis step corresponding to the final value of the time interval in Maxwell, and call the imported ohmic loss value in the corresponding Maxwell time interval in the time interval of each analysis step in the thermal module.

5. The simulation calculation method for the temperature field of roller moving induction heating according to claim 1 is characterized in that: The specific setting described in the above step ④ is: setting the outer contour of the roller body as its heat radiation surface and convection heat transfer surface at different times and in different areas, and setting the radiation coefficient and convection heat transfer coefficient.

6. The method for simulating and calculating the temperature field of roller moving induction heating according to claim 1, characterized in that: The above step ⑤ calculates and obtains the temperature rise curve of the roll under the proposed process parameter conditions and the temperature field distribution and changes during the heating process.

Images