Simulation method and device for metal rheology of wide and thick plate rolling

Abstract

The invention relates to a wide and thick plate rolling metal rheology simulation method and device, and belongs to the technical field of wide and thick plate rolling, and the method comprises the steps: determining the physical property parameters of steel and the physical property parameters of a roller according to the steel grade components of a wide and thick plate to be rolled, and storing the relationship between the physical property parameters of each material and the temperature in a database, the rolling pass, rolling reduction and rolling speed of the wide and thick plate are determined according to the field rolling condition of the wide and thick plate; carrying out unsteady three-dimensional thermal / mechanical coupling analysis on an intermediate billet and a roller in the rolling process of the wide and thick plate, and calculating heat transfer and deformation behaviors of the slab in the rolling process of the wide and thick plate; a comprehensive rolling method is adopted to carry out multi-pass rolling on the wide and thick plate, and rolling of the next pass is carried out after calculation of the last pass is completed. According to the method, the rolling process of the wide and thick plate is effectively controlled and optimized; by predicting and simulating various phenomena and behaviors in the wide and thick plate rolling process, the rolling efficiency is improved, the energy consumption is reduced, and the product performance is improved.

Description

Technical Field

[0001] The invention relates to a method and a device for simulating metal rheology during wide and thick plate rolling, belonging to the technical field of wide and thick plate rolling. Background Art

[0002] Plate and strip products are an important support for the development of the national economy. Their preparation process covers large-scale complex industrial processes with multiple processes and multiple control levels. Plate products can be divided into five categories according to their thickness, namely extra-thick plates, thick plates, medium plates, thin plates and ultra-thin strip plates. Plates with a finished steel plate thickness greater than 4.5 mm are collectively referred to as medium and thick plates. Wide and thick plates are one type of medium and thick plates. Products manufactured by rolling mills with a roll width ranging from 2800 mm to 5500 mm are collectively referred to as wide and thick plates. Wide and thick plate products have excellent performance and are widely used in large oil pipelines, marine engineering, bridges, large ships, pressure vessels and other fields.

[0003] Wide and thick plates frequently have large-scale edge crack defects. In order to reduce edge crack defects, in the actual production process, the only way is to increase the width of the steel plate to obtain a wider plate to cut off the edge with edge cracks. Therefore, the yield rate of wide and thick plates is greatly reduced. Therefore, it is necessary to control the edge cracks of wide and thick plates and improve the yield rate of wide and thick plates. The premise of controlling the edge cracks of wide and thick plates is to have a clear understanding of the metal flow at the edge of the wide and thick plates. From the ingot state to the finished product state, the ingot undergoes a rolling process. However, during the rolling process, as the number of rolling passes increases, the metal will undergo rheology. During the rheological process, the flow of the edge metal is the largest. Through the flow of the edge metal, we can have a clear understanding of the edge cracks of the wide and thick plates, while the traditional experimental method is difficult to accurately analyze the flow behavior of the metal. With the development of computer technology, the flow of metal can be simulated by computer simulation technology.

[0004] To this end, the present invention proposes a simulation method for metal rheology during wide and thick plate rolling. Summary of the invention

[0005] In order to solve the above problems, the present invention proposes a method and device for simulating metal rheology during wide and thick plate rolling, which can simulate the metal flow behavior during wide and thick plate rolling.

[0006] The technical solution adopted by the present invention to solve the technical problem is:

[0007] In a first aspect, an embodiment of the present invention provides a method for simulating metal rheology during wide and thick plate rolling, comprising the following steps:

[0008] Step S1, according to the steel composition of the wide and thick plate to be rolled, determine the physical property parameters of the steel and the physical property parameters of the rolling roller, store the relationship between the physical property parameters of each material and the temperature in a database, and determine the pass, reduction and rolling speed of the wide and thick plate according to the on-site rolling conditions of the wide and thick plate; the physical property parameters of the steel include the thermal conductivity, density, enthalpy change, Poisson's ratio, elastic modulus and yield limit of the steel, and the physical property parameters of the rolling roller include the density, elastic modulus and Poisson's ratio of the rolling roller;

[0009] Step S2, performing a non-steady-state three-dimensional thermal / mechanical coupling analysis on the intermediate billet and the rollers during the wide and thick plate rolling process, and calculating the heat transfer and deformation behavior of the slab during the wide and thick plate rolling process;

[0010] Step S3, adopting the comprehensive rolling method to perform multiple passes of rolling on the wide and thick plates, and performing the next pass of rolling after the calculation of the previous pass is completed.

[0011] As a possible implementation of this embodiment, step S2 includes the following steps:

[0012] Step S201, establishing a three-dimensional solid model of the roll and a three-dimensional solid model of the 1 / 4 cast billet according to the sizes of the intermediate billet and the roll in the process of wide and thick plate rolling;

[0013] Step S202, performing eight-node hexahedral mesh division on the three-dimensional solid model of the roll and the three-dimensional solid model of the 1 / 4 ingot respectively;

[0014] Step S203, setting the mesh units of the intermediate billet and roll models after meshing to unit types that can perform thermal / structural analysis, and assigning the physical property parameters of the steel and the physical property parameters of the roll to the intermediate billet and the roll respectively;

[0015] Step S204, defining the roller as a rigid body capable of heat transfer, defining the intermediate blank as a deformable body, and using the Coulomb bilinear analysis method to analyze the contact between the roller and the workpiece;

[0016] Step S205, setting the initial temperatures of the intermediate billet and the rolls based on the type of steel of the thick plate to be rolled, its starting rolling temperature, and the actual temperature of the rolls;

[0017] Step S206, setting the heat transfer conditions at the boundary between the roller and the intermediate billet;

[0018] Step S207, setting the mechanical boundary conditions of the roll and the intermediate billet, setting the friction coefficient between the roll and the workpiece, and using the Coulomb bilinear analysis method to analyze the contact between the roll and the workpiece;

[0019] Step S208, applying three working conditions to calculate the rolling process of the intermediate billet, and obtaining the rolling interval time, pre-contact movement time and rolling time of the intermediate billet;

[0020] Step S209, using transient thermal / structural statics analysis method to analyze the heat transfer and mechanics between the roller and the intermediate billet, obtaining the displacement data of the edge nodes of the intermediate billet, the analysis option of heat transfer and mechanics is large strain, and the analysis task results are stress, equivalent stress, equivalent plastic strain and temperature increment;

[0021] Step S210, obtaining a relative displacement cloud diagram of edge metal relative to the metal at the center position of the edge of the intermediate billet according to the node displacement data of the edge of the intermediate billet.

[0022] As a possible implementation method of this embodiment, step S202 is specifically as follows: according to the specifications of the ingot, a non-uniform grid division method is used to encrypt the intermediate ingot between 0 and 60 mm close to the surface, and the center of the wide surface and the center of the narrow surface of the ingot are gradually encrypted toward the corner. The internal grid can be divided loosely, and the grid size is controlled to be 2 to 15 mm; the rolling roller is divided into an O-shaped grid.

[0023] As a possible implementation of this embodiment, the heat transfer conditions include: convection heat transfer between cooling water and high-temperature rolled pieces, heat generation by plastic deformation during the deformation of the intermediate billet, contact conduction heat transfer between the rolled pieces and the rollers, and radiation heat transfer between the rolled pieces and the air;

[0024] The convective heat transfer between the cooling water and the high-temperature rolled piece is:

[0025] h w (x) = αW 0.55 (1-0.075T w ) (1)

[0026] In the formula, h w (x) is the heat transfer coefficient between the spray water and the rolled piece; α is the correction coefficient; W is the water flow density; Tw is the cooling water temperature;

[0027] The contact conduction heat transfer between the rolled piece and the roller is:

[0028] q c =h c (T s -T r ) (2)

[0029] In the formula, q c is the contact heat flow; h c is the heat transfer coefficient on the contact surface; T s is the surface temperature of the rolled piece; T r is the surface temperature of the roller;

[0030] Heat transfer coefficient h c for:

[0031]

[0032] In the formula, k glu is the thermal conductivity of the interface between the workpiece and the roller; λ is the roughness of the contact surface; P is the pressure on the contact interface; σ rup t ure is the ultimate strength of the material;

[0033] The radiation heat transfer between the rolled product and the air is:

[0034] h rad =εσ((T surf +273)+(T env +273))((T surf +273) 2 +(T env +273) 2 ) (4) In the formula, h rad is the radiation heat transfer coefficient; ε is the radiation coefficient; σ is the Boltzmann constant; T surf is the surface temperature of the casting billet; T env is the ambient temperature;

[0035] The heat generated by plastic deformation during the intermediate billet deformation process is:

[0036] Q p =ηW p , 0<η<1 (6)

[0037] Where W p is the plastic deformation work;

[0038] The plastic deformation work per unit volume is:

[0039]

[0040] In the formula, σ f is the deformation resistance; f is the strain produced by the rolled piece.

[0041] As a possible implementation method of this embodiment, in step S206, the contact heat transfer between the workpiece and the roller, the radiation heat transfer between the workpiece and the air, and the heat transfer boundary conditions of the contact surface between the roller and the intermediate billet are set to convective heat transfer. The heat transfer boundary is imposed by the convective heat transfer coefficient, and the heat generated by the plastic deformation during the deformation of the intermediate billet is applied to all unit nodes of the intermediate billet model.

[0042] As a possible implementation method of this embodiment, in step S207, the workpiece moves in the rolling direction with the rolling speed as the speed boundary condition, the symmetric surface of the wide face of the ingot is set to have a displacement of 0 along the narrow face of the ingot, and the symmetric surface of the narrow face of the ingot is set to have a displacement of 0 along the wide face of the ingot. The rolling is carried out under the condition that the rollers are in contact with the workpiece, the metal slides along the surface of the rollers, and the rollers generate friction to prevent the flow of metal.

[0043] As a possible implementation method of this embodiment, in step S208, the calculation time of working condition 1 is the interval time. In working condition 1, the workpiece is not in motion and has no contact with the roller, and the workpiece is only subjected to radiation heat transfer; the calculation time of working condition 2 is the movement time before contact, which is the bite process of the workpiece. The driving surface drives the workpiece to contact the roller according to the rolling speed and the distance between the workpiece and the roller. The loads in working condition 2 include radiation heat transfer and convection heat transfer of the nozzle; the calculation time of working condition 3 is the rolling time, which is a stable rolling process. In this stage, the driving surface retreats after the workpiece contacts the roller, and the movement of the workpiece depends entirely on the friction between the workpiece and the roller. The loads in working condition 3 include radiation heat transfer, convection heat transfer of the nozzle, and contact heat transfer between the workpiece and the roller.

[0044] As a possible implementation of this embodiment, in step S3, the model needs to be updated before each pass calculation to achieve the continuity of metal flow of the rolled piece during multi-pass rolling deformation and to be more in line with the actual multi-pass hot rolling process.

[0045] As a possible implementation of this embodiment, step S3 includes the following steps:

[0046] The first pass of rolling simulation is modeled according to the rolling process and material parameters, and the initial load and boundary conditions are applied to complete the first pass of rolling simulation. The model calculation results are then saved using the pre-state analysis.

[0047] The stress, strain and temperature field distribution of the rolled piece in the first pass are used as the initial load for the second pass. The position and speed of the roll are adjusted according to the reduction and speed, and the initial conditions and boundary conditions of the second pass are applied to perform the second pass rolling calculation. The above operation is repeated until all rolling passes are completed.

[0048] In a second aspect, an embodiment of the present invention provides a simulation device for metal rheology during wide and thick plate rolling, comprising:

[0049] The physical property parameter determination module is used to determine the physical property parameters of the steel and the physical property parameters of the rolling roller according to the steel composition of the wide and thick plate to be rolled, store the relationship between the physical property parameters of each material and the temperature in the database, and determine the pass, reduction and rolling speed of the wide and thick plate according to the on-site rolling conditions of the wide and thick plate; the physical property parameters of the steel include the thermal conductivity, density, enthalpy change, Poisson's ratio, elastic modulus and yield limit of the steel, and the physical property parameters of the rolling roller include the density, elastic modulus and Poisson's ratio of the rolling roller;

[0050] Thermal / mechanical coupling analysis module, used to perform non-steady-state three-dimensional thermal / mechanical coupling analysis on the intermediate billet and the rolls during the wide and thick plate rolling process, and calculate the heat transfer and deformation behavior of the slab during the wide and thick plate rolling process;

[0051] The multi-pass rolling module is used to perform multi-pass rolling on wide and thick plates using a comprehensive rolling method, and the next pass of rolling is performed after the calculation of the previous pass is completed.

[0052] In a third aspect, an embodiment of the present invention provides an electronic device, comprising a processor, a memory and a bus, wherein the memory stores machine-readable instructions executable by the processor, and when the electronic device is running, the processor communicates with the memory via the bus, and the processor executes the machine-readable instructions to perform the steps of any of the above-mentioned methods for simulating the rheology of metal rolling of wide and thick plates.

[0053] In a fourth aspect, an embodiment of the present invention provides a storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the steps of the method for simulating metal rheology during rolling of any thick and wide plate described above are executed.

[0054] The beneficial effects of the technical solution of the embodiment of the present invention are as follows:

[0055] The present invention analyzes the physical properties of steel and rollers and the rolling process, simulates the metal flow behavior during the rolling of wide and thick plates, and realizes effective control and optimization of the wide and thick plate rolling process; by combining thermodynamics and non-steady-state three-dimensional simulation technology, various phenomena and behaviors in the wide and thick plate rolling process can be accurately predicted and simulated, providing new possibilities for improving rolling efficiency, reducing energy consumption and improving product performance.

[0056] The present invention can clarify the formation mechanism of edge cracks of wide and thick plates by simulating the metal rheological behavior, stress distribution and metal rollover behavior of the edge of the intermediate billet in multi-pass rolling of wide and thick plates. In this way, the real wide and thick plate rolling process is simulated and optimized, which can not only help understand and improve the existing rolling process, but also provide a basis for developing new rolling processes. BRIEF DESCRIPTION OF THE DRAWINGS

[0057] Figure 1 is a flow chart of a simulation method for metal rheology during wide and thick plate rolling according to an exemplary embodiment;

[0058] Figure 2 is a schematic structural diagram of a device for simulating metal rheology during wide and thick plate rolling according to an exemplary embodiment;

[0059] Figure 3 A specific flow chart of using the simulation device of the present invention to simulate the metal rheology of wide and thick plate rolling;

[0060] Figure 4 It is a motion state diagram of various working conditions for wide and thick plate rolling simulation;

[0061] Figure 5It is a flow chart of pre-state analysis of wide and thick plate rolling simulation. DETAILED DESCRIPTION

[0062] In order to more clearly illustrate the technical features of the solution of the present invention, the present invention is described in detail below through specific implementation methods and in conjunction with the accompanying drawings.

[0063] like Figure 1 As shown, a method for simulating metal rheology during wide and thick plate rolling provided by an embodiment of the present invention comprises the following steps:

[0064] Step S1, according to the steel composition of the wide and thick plate to be rolled, determine the physical property parameters of the steel and the physical property parameters of the rolling roller, store the relationship between the physical property parameters of each material and the temperature in a database, and determine the pass, reduction and rolling speed of the wide and thick plate according to the on-site rolling conditions of the wide and thick plate; the physical property parameters of the steel include the thermal conductivity, density, enthalpy change, Poisson's ratio, elastic modulus and yield limit of the steel, and the physical property parameters of the rolling roller include the density, elastic modulus and Poisson's ratio of the rolling roller;

[0065] Step S2, performing a non-steady-state three-dimensional thermal / mechanical coupling analysis on the intermediate billet and the rollers during the wide and thick plate rolling process, and calculating the heat transfer and deformation behavior of the slab during the wide and thick plate rolling process;

[0066] Step S3, adopting the comprehensive rolling method to perform multiple passes of rolling on the wide and thick plates, and performing the next pass of rolling after the calculation of the previous pass is completed.

[0067] As a possible implementation of this embodiment, step S2 includes the following steps:

[0068] Step S201, establishing a three-dimensional solid model of the roll and a three-dimensional solid model of the 1 / 4 cast billet according to the sizes of the intermediate billet and the roll in the wide and thick plate rolling process;

[0069] Step S202, performing eight-node hexahedral mesh division on the three-dimensional solid model of the roll and the three-dimensional solid model of the 1 / 4 ingot respectively;

[0070] Step S203, setting the mesh units of the intermediate billet and roll models after meshing to unit types that can perform thermal / structural analysis, and assigning the physical property parameters of the steel and the physical property parameters of the roll to the intermediate billet and the roll respectively;

[0071] Step S204, defining the roller as a rigid body capable of heat transfer, defining the intermediate blank as a deformable body, and using the Coulomb bilinear analysis method to analyze the contact between the roller and the workpiece;

[0072] Step S205, setting the initial temperatures of the intermediate billet and the rolls based on the type of steel of the thick plate to be rolled, its starting rolling temperature, and the actual temperature of the rolls;

[0073] Step S206, setting the heat transfer conditions at the boundary between the roller and the intermediate billet;

[0074] Step S207, setting the mechanical boundary conditions of the roll and the intermediate billet, setting the friction coefficient between the roll and the workpiece, and using the Coulomb bilinear analysis method to analyze the contact between the roll and the workpiece;

[0075] Step S208, applying three working conditions to calculate the rolling process of the intermediate billet, and obtaining the rolling interval time, pre-contact movement time and rolling time of the intermediate billet;

[0076] Step S209, using transient thermal / structural statics analysis method to analyze the heat transfer and mechanics between the roller and the intermediate billet, obtaining the displacement data of the edge nodes of the intermediate billet, the analysis option of heat transfer and mechanics is large strain, and the analysis task results are stress, equivalent stress, equivalent plastic strain and temperature increment;

[0077] Step S210, obtaining a relative displacement cloud diagram of edge metal relative to the metal at the center position of the edge of the intermediate billet according to the node displacement data of the edge of the intermediate billet.

[0078] As a possible implementation method of this embodiment, step S202 is specifically as follows: according to the specifications of the ingot, a non-uniform grid division method is used to encrypt the intermediate ingot between 0 and 60 mm close to the surface, and the center of the wide surface and the center of the narrow surface of the ingot are gradually encrypted toward the corner. The internal grid can be divided loosely, and the grid size is controlled to be 2 to 15 mm; the rolling roller is divided into an O-shaped grid.

[0079] As a possible implementation of this embodiment, the heat transfer conditions include: convection heat transfer between cooling water and high-temperature rolled pieces, heat generation by plastic deformation during the deformation of the intermediate billet, contact conduction heat transfer between the rolled pieces and the rollers, and radiation heat transfer between the rolled pieces and the air;

[0080] The convective heat transfer between the cooling water and the high-temperature rolled piece is:

[0081] h w (x) = αW 0.55 (1-0.075T w ) (1)

[0082] In the formula, h w (x) is the heat transfer coefficient between the spray water and the rolled piece; α is the correction coefficient; W is the water flow density; Tw is the cooling water temperature;

[0083] The contact conduction heat transfer between the rolled piece and the roller is:

[0084] q c =h c (T s -T r) (2)

[0085] In the formula, q c is the contact heat flow; h c is the heat transfer coefficient on the contact surface; T s is the surface temperature of the rolled piece; T r is the surface temperature of the roller;

[0086] Heat transfer coefficient h c for:

[0087]

[0088] In the formula, k glu is the thermal conductivity of the interface between the workpiece and the roller; λ is the roughness of the contact surface; P is the pressure on the contact interface; σ rupture is the ultimate strength of the material;

[0089] The radiation heat transfer between the rolled product and the air is:

[0090] h rad =εσ((T surf +273)+(T env +273))((T surf +273) 2 +(T env +273) 2 ) (4)

[0091] In the formula, h rad is the radiation heat transfer coefficient; ε is the radiation coefficient; σ is the Boltzmann constant; T surf is the surface temperature of the casting billet; T env is the ambient temperature;

[0092] The heat generated by plastic deformation during the intermediate billet deformation process is:

[0093] Q p =ηW p , 0<η<1 (6)

[0094] Where W p is the plastic deformation work;

[0095] The plastic deformation work per unit volume is:

[0096]

[0097] In the formula, σ f is the deformation resistance; f is the strain produced by the rolled piece.

[0098] As a possible implementation method of this embodiment, in step S206, the contact heat transfer between the workpiece and the roller, the radiation heat transfer between the workpiece and the air, and the heat transfer boundary conditions of the contact surface between the roller and the intermediate billet are set to convective heat transfer. The heat transfer boundary is imposed by the convective heat transfer coefficient, and the heat generated by the plastic deformation during the deformation of the intermediate billet is applied to all unit nodes of the intermediate billet model.

[0099] As a possible implementation method of this embodiment, in step S207, the workpiece moves in the rolling direction with the rolling speed as the speed boundary condition, the symmetric surface of the wide face of the ingot is set to have a displacement of 0 along the narrow face of the ingot, and the symmetric surface of the narrow face of the ingot is set to have a displacement of 0 along the wide face of the ingot. The rolling is carried out under the condition that the rollers are in contact with the workpiece, the metal slides along the surface of the rollers, and the rollers generate friction to prevent the flow of metal.

[0100] As a possible implementation method of this embodiment, in step S208, the calculation time of working condition 1 is the interval time. In working condition 1, the workpiece is not in motion and has no contact with the roller, and the workpiece is only subjected to radiation heat transfer; the calculation time of working condition 2 is the movement time before contact, which is the bite process of the workpiece. The driving surface drives the workpiece to contact the roller according to the rolling speed and the distance between the workpiece and the roller. The loads in working condition 2 include radiation heat transfer and convection heat transfer of the nozzle; the calculation time of working condition 3 is the rolling time, which is a stable rolling process. In this stage, the driving surface retreats after the workpiece contacts the roller, and the movement of the workpiece depends entirely on the friction between the workpiece and the roller. The loads in working condition 3 include radiation heat transfer, convection heat transfer of the nozzle, and contact heat transfer between the workpiece and the roller.

[0101] As a possible implementation of this embodiment, in step S3, the model needs to be updated before each pass calculation to achieve the continuity of metal flow of the rolled piece during multi-pass rolling deformation and to be more in line with the actual multi-pass hot rolling process.

[0102] As a possible implementation of this embodiment, step S3 includes the following steps:

[0103] The first pass of rolling simulation is modeled according to the rolling process and material parameters, and the initial load and boundary conditions are applied to complete the first pass of rolling simulation. The model calculation results are then saved using the pre-state analysis.

[0104] The stress, strain and temperature field distribution of the rolled piece in the first pass are used as the initial load for the second pass. The position and speed of the roll are adjusted according to the reduction and speed, and the initial conditions and boundary conditions of the second pass are applied to perform the second pass rolling calculation. The above operation is repeated until all rolling passes are completed.

[0105] like Figure 2 As shown, a simulation device for metal rheology during wide and thick plate rolling provided by an embodiment of the present invention comprises:

[0106] The physical property parameter determination module is used to determine the physical property parameters of the steel and the physical property parameters of the rolling roller according to the steel composition of the wide and thick plate to be rolled, store the relationship between the physical property parameters of each material and the temperature in the database, and determine the pass, reduction and rolling speed of the wide and thick plate according to the on-site rolling conditions of the wide and thick plate; the physical property parameters of the steel include the thermal conductivity, density, enthalpy change, Poisson's ratio, elastic modulus and yield limit of the steel, and the physical property parameters of the rolling roller include the density, elastic modulus and Poisson's ratio of the rolling roller;

[0107] Thermal / mechanical coupling analysis module, used to perform non-steady-state three-dimensional thermal / mechanical coupling analysis on the intermediate billet and the rolls during the wide and thick plate rolling process, and calculate the heat transfer and deformation behavior of the slab during the wide and thick plate rolling process;

[0108] The multi-pass rolling module is used to perform multi-pass rolling on wide and thick plates using a comprehensive rolling method, and the next pass of rolling is performed after the calculation of the previous pass is completed.

[0109] Taking the wide and thick plates actually produced by a steel plant as an example, the specific process of simulating the metal rheology of the wide and thick plates by using the simulation device of the present invention is as follows: Figure 3 shown.

[0110] Step 1. According to the composition of the wide and thick plate steel actually produced by a steel plant as shown in Table 1, the JMatPro software is used to calculate the thermal physical properties of the steel, including density, Poisson's ratio, expansion coefficient, thermal conductivity, and elastic modulus. On this basis, the density, elastic modulus, and Poisson's ratio of the rollers shown in Table 2 are obtained by consulting the standard tool manual, and the relationship between the physical properties of each material and the temperature is imported into the corresponding commercial software database. Further, according to the actual production of the enterprise, the number of passes, reduction, and rolling speed of the wide and thick plate rolling as shown in Table 3 are determined, including 3 passes of rough rolling and horizontal rolling, 4 passes of rough rolling and longitudinal rolling, and 8 passes of finishing rolling. The 1800×250mm slab is rolled into a 2935×30mm thick plate with a width ratio of 1.63.

[0111] Table 1: Composition of wide and thick plate steel grades

[0112]

[0113]

[0114] Table 2: Roller physical parameters

[0115] Material Poisson's ratio Elastic modulus, Pa Density, kg / m3 rigid body 0.3 2×1011 7850

[0116] Table 3: Rolling pass information

[0117] Rolling pass Pressing amount (mm) Rolling speed (rad / s) R1 22.56 1 R2 19.74 1 R3 15.99 1 R4 18.85 2 R5 29.86 2 R6 27.34 2 R7 27.78 2 F1 10.76 3.4 F2 10.57 3.7 F3 10.05 4 F4 9.06 4.3 F5 6.83 4.7 F6 5.37 5 F7 4.15 5.3 F8 3.22 5.6

[0118] Step 2: Use MSC Marc commercial software to perform a three-dimensional thermal / mechanical coupling analysis of the intermediate billet and the rolls to calculate the heat transfer and deformation behavior of the slab during the wide and thick plate rolling process.

[0119] Step 2-1, according to the actual intermediate billet and roller size, establish a three-dimensional solid model of the roller and 1 / 4 cast billet. In this embodiment, the roller diameter of a steel plant is 1120mm, the length is 1300mm, the length of the rolled piece is selected to be 1000mm, the thickness is 250mm, and the width is 1800mm.

[0120] Step 2-2, import the solid models of the roll and the billet into the meshing software ICEM CFD for eight-node hexahedral meshing. According to the billet specifications, the non-uniform meshing method is adopted. Since the main research is on the metal rheological behavior of the edge of the intermediate billet, the intermediate billet between 0 and 60 mm close to the surface is encrypted, and the center of the wide face and the center of the narrow face of the billet are gradually encrypted towards the corner. The internal grid can be divided loosely, and the grid size is controlled to be 2 to 15 mm. Due to the special shape of the roll, O-shaped meshing is adopted, and the grid size is 45 mm.

[0121] Step 2-3: Import the meshed roll and intermediate billet models into the MSC Marc software, set the mesh unit of the model to unit type 7 for performing thermal / structural analysis, and assign the physical property parameters imported into the material database in step 1 to the roll and intermediate billet respectively.

[0122] Step 2-4: Define the roller as a heat-conducting rigid body, the intermediate billet as a deformable body, and the contact between the roller and the workpiece is set using a rigid-flexible contact analysis algorithm. The contact between the roller and the workpiece is analyzed using Coulomb bilinear (displacement) analysis.

[0123] Step 2-5: Set the initial temperatures of the intermediate slab and the rolls based on the actual mainstream steel type produced on site, its rolling temperature, and the actual roll temperature. Assuming that the temperature of the ingot is uniform after being heated in the heating furnace, the initial temperature of the slab is set to 1200°C and the initial temperature of the roll is set to 25°C.

[0124] Step 2-6, set the heat transfer conditions of the roller and the intermediate billet boundary. During the whole hot rolling simulation process, the heat transfer phenomenon is different at different positions of the intermediate billet, and the heat transfer boundary conditions are also different, mainly in the following three aspects: convection heat transfer between cooling water and high-temperature rolled pieces, heat generated by plastic deformation during the deformation of the intermediate billet, contact conduction heat transfer between the rolled pieces and the rollers, and radiation heat transfer between the rolled pieces and the air. Set the contact heat transfer between the rolled pieces and the rollers, the radiation heat transfer between the rolled pieces and the air, and the heat transfer boundary conditions of the contact surface between the rollers and the intermediate billet as convection heat transfer. The heat transfer boundary is imposed by the convection heat transfer coefficient, and the heat generated by plastic deformation during the deformation of the intermediate billet is applied to all unit nodes of the intermediate billet model.

[0125] The convective heat transfer between the cooling water and the high temperature rolled piece is shown in the relationship (1):

[0126] h w (x) = αW 0.55 (1-0.075T w ) (1)

[0127] In the formula, h w (x) is the heat transfer coefficient between the spray water and the rolled piece, W / (m 2 ·℃); α is the correction factor; W is the water flow density, L / (s·m 2 );T w is the cooling water temperature, °C,

[0128] The contact conduction heat transfer between the workpiece and the roller is obtained according to the Fourier heat conduction equation as shown in formula (2):

[0129] q c =h c (T s -T r ) (2)

[0130] In the formula, q c is the contact heat flow; h c is the heat transfer coefficient on the contact surface; T s is the surface temperature of the rolled piece; T r is the surface temperature of the roller. Heat transfer coefficient h c The function related to rolling state, interface medium and temperature is shown in formula (3) using the heat transfer coefficient:

[0131]

[0132] In the formula, k glu is the thermal conductivity of the interface between the workpiece and the roller; λ is the roughness of the contact surface; P is the pressure on the contact interface; σ rup t ure The ultimate strength of the material.

[0133] Although the radiation heat transfer between the rolled product and the environment accounts for a small proportion of the total heat transfer, it exists in the entire rolling process and cannot be ignored. The radiation heat transfer coefficient is determined by formula (4):

[0134] h rad =εσ((T surf +273)+(T env +273))((T surf +273) 2 +(T env +273) 2 ) (4)

[0135] In the formula, h rad is the radiation heat transfer coefficient, W / (m 2 ·℃); ε is the radiation coefficient, which is a constant of 0.8; σ is the Boltzmann constant, which is 5.67×10-8(W / m 2 / ℃4);T surf is the surface temperature of the casting, °C; T env is the ambient temperature, ℃.

[0136] During the hot rolling process, the rolled piece undergoes a large plastic deformation, so it is necessary to consider the heat generated by the rolled piece during plastic deformation. The plastic deformation work is related to the instantaneous yield stress and strain of the rolled piece. The plastic deformation work per unit volume Wp is expressed by formula (5):

[0137]

[0138] In the formula, σ f is the deformation resistance; f is the strain generated by the rolled piece. In this way, the heat generated by plastic work Qp is expressed by formula (6):

[0139] Q p =ηW p , 0<η<1 (6).

[0140] Step 2-7, set the mechanical boundary conditions of the roll and the intermediate billet. The present invention selects the calculation model as the 1 / 4 model. The rolled piece moves in the rolling direction with the rolling speed as the speed boundary condition. The symmetric surface of the wide surface of the billet is set to have a displacement of 0 along the narrow surface of the billet, and the symmetric surface of the narrow surface of the billet is set to have a displacement of 0 along the wide surface of the billet. The rolling is carried out under the condition that the roll and the rolled piece are in contact. The metal slides along the surface of the roll, and the roll generates friction to prevent the metal from flowing. The friction coefficient between the roll and the rolled piece is set. In this embodiment, the friction coefficient is 0.3. The contact between the roll and the rolled piece adopts Coulomb bilinear (displacement) analysis.

[0141] The basic laws of plasticity used in rolling include the conservation of mass, momentum, and energy, which are shown in (7) to (13):

[0142] (1) Conservation of mass

[0143] If the mass m(B) remains constant for any selected configuration, it is called the law of conservation of mass. In region R, the rate of increase of mass entering object B is equal to the rate of decrease of mass leaving object B, that is:

[0144]

[0145] In the formula, is the growth rate of ρ at a fixed point in R. Applying the divergence theorem to the surface integral, the above formula can be rewritten as:

[0146]

[0147] Region R is an arbitrary region, and the integrand is zero everywhere in the above formula, and the continuity equation commonly used in engineering calculations can be obtained:

[0148]

[0149] (2) Conservation of momentum

[0150] At a certain moment, the rate of change of linear momentum of each particle of the continuum B in R is proportional to the resultant force acting in R, which is called the law of conservation of momentum. The resultant force is composed of the body force b and the surface traction t (surface tractions), which can be expressed as formula (10):

[0151]

[0152] The integrand is zero, let a j is the acceleration, then a j =dv j / dt, we can get the equation of motion of the continuum (11):

[0153]

[0154] (3) Conservation of Energy

[0155] In R, the time derivatives of the kinetic energy and internal energy of material B are equal to the sum of the mechanical power exerted by the volume force and surface force acting on R and other powers entering R, which is called the law of conservation of energy, as shown in formula (12):

[0156]

[0157] In the formula, K is kinetic energy; E is internal energy; W1 is the work done by volume force; W2 is the work done by surface force; Q1 is the thermal energy rate; Q2 is the energy rate of internal chemical changes; Q3 is the radiation energy rate; Q4 is other energy rates.

[0158] Combining the energy conservation equation with Fourier's law, we can get the unsteady three-dimensional heat conduction equation:

[0159]

[0160] Where λ is the thermal conductivity, W / (m·℃); ρ is the density of the object, kg / m 3 ; Q is the internal heat source, J; c is the specific heat, J / (kg·℃); T is the temperature, ℃; t is the time, s; x, y, z are the coordinates in three directions in three-dimensional space, m.

[0161] During the hot rolling process, after the workpiece and the roller bite into each other, the movement of the workpiece and the flow of the metal are achieved by friction. The present invention uses the Coulomb friction model, as shown in formula (14):

[0162] σ fr ≤-μσ n t (14)

[0163] In the formula, σ n is the normal stress at the contact node; σ fr is the tangential (friction) stress; μ is the friction coefficient; t is the tangential unit vector in the direction of relative sliding velocity.

[0164] The form of the node resultant force of the Coulomb friction model is shown in formula (15):

[0165] f t ≤-μf n t (15)

[0166] In the formula, f t is the tangential force; f n is the normal reaction force.

[0167] The friction coefficient is related to the contact surface temperature and rolling speed. The friction coefficient between the workpiece and the roller is expressed by formula (16):

[0168] μ=1.05-0.0005T-0.056v (16)

[0169] Where T is the workpiece temperature, °C; v is the speed of the roller relative to the workpiece, m / s.

[0170] Step 2-8, in order to reflect the movement process of the intermediate billet in different states during the rolling process, the present invention calculates the rolling process of the intermediate billet by applying three working conditions, such as Figure 4As shown, the calculation time of working condition 1 is the interval time. In this working condition, the workpiece is not in motion and has no contact with the roller. The workpiece is only subjected to radiation heat transfer. The calculation time of working condition 2 is the movement time before contact, which is the bite process of the workpiece. The driving surface drives the workpiece to contact the roller according to the rolling speed and the distance between the workpiece and the roller. The loads in this working condition include radiation heat transfer and convection heat transfer of the nozzle. In this embodiment, the calculation time of working condition 1 and working condition 2 for each pass is 0.1s. The calculation time of working condition 3 is the rolling time, which is a stable rolling process. In this stage, the driving surface retreats after the workpiece contacts the roller, and the movement of the workpiece depends entirely on the friction between the workpiece and the roller. The loads in this working condition include radiation heat transfer, convection heat transfer of the nozzle, and contact heat transfer between the workpiece and the roller. As the rolling process proceeds, the length or thickness of the intermediate billet changes, and the length of the rolled piece in each pass gradually increases as the rolling process proceeds. Therefore, the calculation time of working condition 3 in each pass also increases accordingly. In this embodiment, the calculation time range of working condition 3 is 3.6s to 6.2s.

[0171] Step 2-9: Set the heat transfer and mechanical analysis between the roll and the intermediate billet to transient thermal / structural static analysis. The analysis option is large strain, and the analysis task results are stress, equivalent stress, equivalent plastic strain, and temperature increment.

[0172] Step 2-10, the node displacement data of the intermediate billet edge obtains a relative displacement cloud diagram of the edge metal relative to the metal at the center position of the intermediate billet edge.

[0173] This is achieved by using basic thermodynamic equations, such as conservation of energy and momentum, combined with finite element or boundary element methods. This simulation can help understand the temperature and stress distribution of thick plates during rolling, thereby helping to optimize rolling conditions and improve rolling efficiency.

[0174] Step 3: In actual enterprise production, wide and thick plates are rolled in multiple passes using the integrated rolling method. The second pass is performed after the first pass calculation is completed. Therefore, the model needs to be updated before each calculation to achieve the continuity of the metal flow of the rolled piece during the multi-pass rolling deformation and to be more in line with the actual multi-pass hot rolling process. Figure 5 As shown, in the present invention, the first pass is firstly modeled according to the rolling process and the rolled material parameters, and the initial load and boundary conditions are applied to complete the first pass rolling simulation. Then, the model calculation results are saved using the pre-state analysis in the Marc software, and the stress, strain and temperature field distribution of the first pass rolled piece are used as the initial load for the second pass rolling. The position and speed of the rollers are adjusted according to the reduction and rotation speed, and the initial conditions and boundary conditions of the second pass are applied to perform the second pass rolling calculation, and so on until all rolling passes are completed.

[0175] An electronic device provided by an embodiment of the present invention includes a processor, a memory and a bus, wherein the memory stores machine-readable instructions executable by the processor. When the device is running, the processor communicates with the memory through the bus, and the processor executes the machine-readable instructions to perform the steps of any of the above-mentioned methods for simulating the rheology of metal rolling of wide and thick plates.

[0176] Specifically, the above-mentioned memory and processor can be general-purpose memory and processor, which are not specifically limited here. When the processor runs the computer program stored in the memory, the above-mentioned simulation method of metal rheology in wide and thick plate rolling can be executed.

[0177] Those skilled in the art will appreciate that the structure of the electronic device does not limit the electronic device and may include more or fewer components than shown in the figure, or combine or split certain components, or arrange the components differently.

[0178] In some embodiments, the electronic device may also include a touch screen that can be used to display a graphical user interface (e.g., a startup interface of an application) and receive user operations on the graphical user interface (e.g., startup operations for an application). The specific touch screen may include a display panel and a touch panel. The display panel may be configured in the form of LCD (Liquid Crystal Display), OLED (Organic Light-Emitting Diode), etc. The touch panel may collect the user's contact or non-contact operations on or near it, and generate pre-set operation instructions, for example, the user uses any suitable object or accessory such as a finger, stylus, etc. on the touch panel or near the touch panel. In addition, the touch panel may include two parts: a touch detection device and a touch controller. Among them, the touch detection device detects the user's touch position and posture, detects the signal brought by the touch operation, and transmits the signal to the touch controller; the touch controller receives the touch information from the touch detection device, converts it into information that the processor can process, and then sends it to the processor, and can receive and execute commands sent by the processor. In addition, the touch panel can be implemented by various types such as resistive, capacitive, infrared and surface acoustic wave, and any technology developed in the future can also be used to implement the touch panel. Further, the touch panel can cover the display panel, and the user can operate on or near the touch panel covered on the display panel according to the graphical user interface displayed on the display panel. After the touch panel detects the operation on or near it, it is transmitted to the processor to determine the user input, and then the processor provides corresponding visual output on the display panel in response to the user input. In addition, the touch panel and the display panel can be implemented as two independent components or integrated.

[0179] Corresponding to the startup method of the above-mentioned application, an embodiment of the present invention further provides a storage medium on which a computer program is stored. When the computer program is executed by a processor, the steps of the simulation method of metal rheology of any wide and thick plate rolling are executed.

[0180] The application startup device provided in the embodiment of the present application can be specific hardware on the device or software or firmware installed on the device. The device provided in the embodiment of the present application, its implementation principle and the technical effect produced are the same as those in the aforementioned method embodiment. For the sake of brief description, the parts not mentioned in the device embodiment can refer to the corresponding contents in the aforementioned method embodiment. Those skilled in the art can clearly understand that for the convenience and simplicity of description, the specific working processes of the systems, devices and units described above can all refer to the corresponding processes in the aforementioned method embodiment, and will not be repeated here.

[0181] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented in one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that include computer-usable program code.

[0182] In the embodiments provided in the present application, it should be understood that the disclosed devices and methods can be implemented in other ways. The device embodiments described above are merely schematic. For example, the division of modules is only a logical function division. There may be other division methods in actual implementation. For example, multiple modules or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the mutual coupling or direct coupling or communication connection shown or discussed can be through some communication interfaces, indirect coupling or communication connection of devices or modules, which can be electrical, mechanical or other forms.

[0183] The modules described as separate components may or may not be physically separated, and the components shown as modules may or may not be physical modules, that is, they may be located in one place or distributed on multiple network modules. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0184] In addition, each functional module in the embodiments provided in the present application may be integrated into one processing module, or each module may exist physically separately, or two or more modules may be integrated into one module.

[0185] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0186] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0187] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0188] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the above embodiments, ordinary technicians in the relevant field should understand that the specific implementation methods of the present invention can still be modified or replaced by equivalents. Any modification or equivalent replacement that does not depart from the spirit and scope of the present invention should be covered within the scope of protection of the claims of the present invention.

Claims

1. A method for simulating metal rheology during wide and thick plate rolling, characterized in that: The steps include: Step S1, according to the steel composition of the wide and thick plate to be rolled, determine the physical property parameters of the steel and the physical property parameters of the rolling roller, store the relationship between the physical property parameters of each material and the temperature in a database, and determine the pass, reduction and rolling speed of the wide and thick plate according to the on-site rolling conditions of the wide and thick plate; the physical property parameters of the steel include the thermal conductivity, density, enthalpy change, Poisson's ratio, elastic modulus and yield limit of the steel, and the physical property parameters of the rolling roller include the density, elastic modulus and Poisson's ratio of the rolling roller; Step S2, performing a non-steady-state three-dimensional thermal / mechanical coupling analysis on the intermediate billet and the rollers during the wide and thick plate rolling process, and calculating the heat transfer and deformation behavior of the slab during the wide and thick plate rolling process; Step S3, adopting the comprehensive rolling method to perform multiple passes of rolling on the wide and thick plates, and performing the next pass of rolling after the calculation of the previous pass is completed.

2. The method for simulating metal rheology during thick plate rolling according to claim 1, characterized in that: The step S2 comprises the following steps: Step S201, establishing a three-dimensional solid model of the roll and a three-dimensional solid model of the 1 / 4 cast billet according to the sizes of the intermediate billet and the roll in the wide and thick plate rolling process; Step S202, performing eight-node hexahedral mesh division on the three-dimensional solid model of the roll and the three-dimensional solid model of the 1 / 4 ingot respectively; Step S203, setting the mesh units of the intermediate billet and roll models after meshing to unit types that can perform thermal / structural analysis, and assigning the physical property parameters of the steel and the physical property parameters of the roll to the intermediate billet and the roll respectively; Step S204, defining the roller as a rigid body capable of heat transfer, defining the intermediate blank as a deformable body, and using the Coulomb bilinear analysis method to analyze the contact between the roller and the workpiece; Step S205, setting the initial temperatures of the intermediate billet and the rolls based on the type of steel of the thick plate to be rolled, its starting rolling temperature, and the actual temperature of the rolls; Step S206, setting the heat transfer conditions at the boundary between the roller and the intermediate billet; Step S207, setting the mechanical boundary conditions of the roll and the intermediate billet, setting the friction coefficient between the roll and the workpiece, and using the Coulomb bilinear analysis method to analyze the contact between the roll and the workpiece; Step S208, applying three working conditions to calculate the rolling process of the intermediate billet, and obtaining the rolling interval time, pre-contact movement time and rolling time of the intermediate billet; Step S209, using transient thermal / structural statics analysis method to analyze the heat transfer and mechanics between the roller and the intermediate billet, obtaining the displacement data of the edge nodes of the intermediate billet, the analysis option of heat transfer and mechanics is large strain, and the analysis task results are stress, equivalent stress, equivalent plastic strain and temperature increment; Step S210, obtaining a relative displacement cloud diagram of edge metal relative to the metal at the center position of the edge of the intermediate billet according to the node displacement data of the edge of the intermediate billet.

3. The method for simulating metal rheology during thick plate rolling according to claim 2, characterized in that: The step S202 is specifically as follows: according to the specifications of the ingot, a non-uniform grid division method is used to encrypt the intermediate ingot between 0 and 60 mm close to the surface, and the center of the wide surface and the center of the narrow surface of the ingot are gradually encrypted toward the corner. The internal grid can be divided loosely, and the grid size is controlled to be 2 to 15 mm; the rolling roller is divided into an O-shaped grid.

4. The method for simulating metal rheology during wide and thick plate rolling according to claim 2, characterized in that: The heat transfer conditions include: convection heat transfer between cooling water and high-temperature rolled pieces, heat generation by plastic deformation during the deformation of the intermediate blank, contact conduction heat transfer between the rolled pieces and the rolls, and radiation heat transfer between the rolled pieces and the air; The convective heat transfer between the cooling water and the high-temperature rolled piece is: h w (x)=αW 0.55 (1-0.075T w ) (1) In the formula, h w (x) is the heat transfer coefficient between the spray water and the rolled piece; α is the correction coefficient; W is the water flow density; Tw is the cooling water temperature; The contact conduction heat transfer between the rolled piece and the roller is: q c =h c (T s -T r ) (2) In the formula, q c is the contact heat flow; h c is the heat transfer coefficient on the contact surface; T s is the surface temperature of the rolled piece; T r is the surface temperature of the roller; Heat transfer coefficient h c for: In the formula, k glu is the thermal conductivity of the interface between the workpiece and the roller; λ is the roughness of the contact surface; P is the pressure on the contact interface; σ rup t ure is the ultimate strength of the material; The radiation heat transfer between the rolled piece and the air is: h rad =εσ((T surf +273)+(T env +273))((T surf +273) 2 +(T env +273) 2 ) (4) In the formula, h rad is the radiation heat transfer coefficient; ε is the radiation coefficient; σ is the Boltzmann constant; T surf is the surface temperature of the casting billet; T env is the ambient temperature; The heat generated by plastic deformation during the intermediate billet deformation process is: Q p =ηW p ,0<η<1 (6) Where W p is the plastic deformation work; The plastic deformation work per unit volume is: In the formula, σ f is the deformation resistance; f is the strain produced by the rolled piece.

5. The method for simulating metal rheology during wide and thick plate rolling according to claim 2, characterized in that: In step S207, the workpiece moves in the rolling direction with the rolling speed as the speed boundary condition, the symmetric plane of the wide face of the ingot is set to have a displacement of 0 along the narrow face of the ingot, and the symmetric plane of the narrow face of the ingot is set to have a displacement of 0 along the wide face of the ingot. Rolling is carried out under the condition that the rollers are in contact with the workpiece, the metal slides along the surface of the rollers, and the rollers generate friction to prevent the flow of metal.

6. The method for simulating metal rheology during wide and thick plate rolling according to claim 2, characterized in that: In step S208, the calculation time of working condition 1 is the interval time. In working condition 1, the workpiece is not in motion and has no contact with the roller. The workpiece is only subjected to radiation heat transfer. The calculation time of working condition 2 is the pre-contact movement time, which is the bite process of the workpiece. The driving surface drives the workpiece to contact the roller according to the rolling speed and the distance between the workpiece and the roller. The load in working condition 2 includes radiation heat transfer and convection heat transfer of the nozzle. The calculation time of working condition 3 is the rolling time. In order to stabilize the rolling process, the driving surface retreats after the workpiece contacts the roller. The movement of the workpiece depends entirely on the friction between the workpiece and the roller. The loads in working condition 3 include radiation heat transfer, convection heat transfer of the nozzle, and contact heat transfer between the workpiece and the roller.

7. The method for simulating metal rheology during thick plate rolling according to any one of claims 1 to 6, characterized in that: The step S3 comprises the following steps: The first pass of rolling simulation is modeled according to the rolling process and material parameters, and the initial load and boundary conditions are applied to complete the first pass of rolling simulation. The model calculation results are then saved using the pre-state analysis. The stress, strain and temperature field distribution of the rolled piece in the first pass are used as the initial load for the second pass. The position and speed of the roll are adjusted according to the reduction and speed, and the initial conditions and boundary conditions of the second pass are applied to perform the second pass rolling calculation. The above operation is repeated until all rolling passes are completed.

8. A simulation device for metal rheology during wide and thick plate rolling, characterized in that: include: The physical property parameter determination module is used to determine the physical property parameters of the steel and the physical property parameters of the rolling roller according to the steel composition of the wide and thick plate to be rolled, store the relationship between the physical property parameters of each material and the temperature in the database, and determine the pass, reduction and rolling speed of the wide and thick plate according to the on-site rolling conditions of the wide and thick plate; the physical property parameters of the steel include the thermal conductivity, density, enthalpy change, Poisson's ratio, elastic modulus and yield limit of the steel, and the physical property parameters of the rolling roller include the density, elastic modulus and Poisson's ratio of the rolling roller; Thermal / mechanical coupling analysis module, used to perform non-steady-state three-dimensional thermal / mechanical coupling analysis on the intermediate billet and the rolls during the wide and thick plate rolling process, and calculate the heat transfer and deformation behavior of the slab during the wide and thick plate rolling process; The multi-pass rolling module is used to perform multi-pass rolling on wide and thick plates using a comprehensive rolling method, and the next pass of rolling is performed after the calculation of the previous pass is completed.

9. An electronic device, characterized in that: It includes a processor, a memory and a bus, the memory stores machine-readable instructions executable by the processor, when the electronic device is running, the processor and the memory communicate through the bus, and the processor executes the machine-readable instructions to perform the steps of the simulation method of metal rheology of wide and thick plate rolling as described in any one of claims 1-7.

10. A storage medium, characterized in that: The storage medium stores a computer program, which, when executed by a processor, executes the steps of the method for simulating metal rheology during wide and thick plate rolling as described in any one of claims 1-7.

Images