CN111079275A - Rolled piece temperature obtaining method and device for strip hot rolling production line - Google Patents

Rolled piece temperature obtaining method and device for strip hot rolling production line Download PDF

Info

Publication number
CN111079275A
CN111079275A CN201911236967.9A CN201911236967A CN111079275A CN 111079275 A CN111079275 A CN 111079275A CN 201911236967 A CN201911236967 A CN 201911236967A CN 111079275 A CN111079275 A CN 111079275A
Authority
CN
China
Prior art keywords
rolled piece
temperature
section
area
node
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201911236967.9A
Other languages
Chinese (zh)
Other versions
CN111079275B (en
Inventor
陈全忠
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Johnson Electric SA
Original Assignee
Johnson Electric SA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Johnson Electric SA filed Critical Johnson Electric SA
Priority to CN201911236967.9A priority Critical patent/CN111079275B/en
Publication of CN111079275A publication Critical patent/CN111079275A/en
Application granted granted Critical
Publication of CN111079275B publication Critical patent/CN111079275B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention discloses a rolled piece temperature acquisition method for a plate and strip hot rolling production line, which comprises the following steps: constructing equipment arrangement of a plate strip hot rolling production line, and acquiring related process parameters; dividing a rolled piece heat exchange area according to equipment arrangement on a production line, wherein the heat exchange area comprises an air cooling area, a water cooling area, a roll gap contact heat conduction area and a heat radiation heat preservation area; dividing the grid nodes of the width-thickness section of the rolled piece, and establishing a coordinate system parameter and temperature calculation model of each node; acquiring initial temperature distribution of a rolled piece at an initial calculation position; and from the initial calculation position to the final calculation position, judging the heat exchange area where the section of the rolled piece is positioned when calculating once, and then calculating the temperature distribution of the section of the rolled piece at the position by adopting a corresponding rolled piece temperature model. The invention has the beneficial effects that: the method has clear and definite principle, few assumed and simplified conditions, can accurately forecast the temperature field change of the rolled piece on a given hot rolling production line, and has higher calculation precision than an analytical method.

Description

Rolled piece temperature obtaining method and device for strip hot rolling production line
Technical Field
The invention belongs to the technical field of hot-rolled strips, and particularly relates to a method and a device for acquiring the temperature of a rolled piece in a strip hot-rolling production line.
Background
The accurate prediction of the temperature change of the plate strip in the hot rolling process is an important prerequisite for ensuring the size precision, the plate shape quality and the structure performance of the plate strip, and has important significance for ensuring the stable production of the hot rolling plate strip: on one hand, the plate strip hot rolling process is a very complex heat engine coupling process and comprises comprehensive heat exchange of contact heat conduction, natural convection, forced convection, high-temperature heat radiation, a deformation internal heat source and friction heat; the existing rolled piece temperature model mainly comprises an analytic method model and a finite element method model, wherein the analytic method model does not fully consider the complexity of boundary conditions and field working conditions, and adopts more hypothesis simplifying conditions, so that the calculation precision is low, while the finite element method model belongs to a three-dimensional numerical model, although the calculation precision is higher, the calculation speed is very slow, and the online real-time requirement of engineering calculation cannot be met; on the other hand, the hot rolling production line has complex process flow and multiple changes of equipment arrangement forms, and at present, no general rolled piece temperature calculation model and method suitable for various types of hot rolling production lines exist.
Disclosure of Invention
The invention aims to provide a rolled piece temperature obtaining method and device for a plate and strip hot rolling production line aiming at the defects of the prior art and solving the problems that the temperature model of the prior rolled piece is poor in universality, low in calculation precision, incapable of quickly and conveniently adapting to engineering calculation requirements and the like.
The invention provides a rolled piece temperature obtaining method for a plate and strip hot rolling production line, which is characterized by comprising the following steps of:
step one, constructing equipment arrangement of a plate strip hot rolling production line, and acquiring relevant process parameters;
dividing a rolled piece heat exchange area according to equipment arrangement on the production line, wherein the heat exchange area comprises an air cooling area, a water cooling area, a roll gap contact heat conduction area and a heat radiation heat preservation area, and the water cooling area is divided into a high-pressure water descaling water cooling area and a laminar cooling water cooling area;
dividing rolled piece width direction-thickness direction section grid nodes, and establishing a coordinate system parameter and temperature calculation model of each node;
acquiring initial temperature distribution of the section of the rolled piece at the initial calculation position;
step five, the calculated rolled piece section moves along the rolling direction according to the running speed of the strip steel, and the temperature distribution of the rolled piece section is calculated from the initial calculation position to the termination calculation position; and judging the heat exchange area where the cross section of the rolled piece is positioned when calculating once, and then calculating the temperature distribution of the cross section of the rolled piece at the position by adopting a corresponding rolled piece temperature model.
According to the scheme, the equipment layout of the plate strip hot rolling production line comprises the positions and distances of a rolling mill unit, a high-pressure water descaling device, a heat preservation cover, a heating furnace, a laminar cooling device and a coiling machine, and the process parameters comprise the process parameters of the rolling mill, the rolling parameters of rolled pieces, the process parameters of laminar cooling, thermophysical parameters and medium parameters.
According to the scheme, the technological parameters of the rolling mill comprise the roll body diameter, the roll body length, the roll neck diameter and the roll shifting displacement of the working roll of each frame; the rolling parameters of the rolled piece comprise the speed of the rolled piece, the width of the rolled piece, the inlet thickness and the outlet thickness of a roll gap of each frame and the yield strength of the rolled piece; the laminar cooling process parameters comprise the number of the jet beams, the jet length of the jet beams, the distance between the adjacent jet beams and the opening number and positions of the jet beams; the thermophysical parameters comprise the heat conductivity of the rolled piece, the specific heat capacity, the thermal radiance of the rolled piece, the contact heat exchange coefficient of the roll gap of each frame, the air cooling heat exchange coefficient of the rolled piece, the cooling water heat exchange coefficient of the rolled piece and the high-pressure water descaling heat exchange coefficient of the rolled piece; the medium parameters include cooling water temperature and ambient temperature.
According to the scheme, in the third step, grid nodes of the width direction-thickness direction section of the rolled piece are divided, and a coordinate system position and temperature calculation model of each node is established, and the method specifically comprises the following steps: establishing a y-z rectangular coordinate system, wherein a y axis is positioned at the position of the middle thickness of the rolled piece, and a z axis is positioned at the position of the middle width of the rolled piece; dispersing a half width-thickness cross section of a rolled piece into N multiplied by M grids, wherein the half width of the rolled piece is equally divided into N sections, i is 1, 2 and 3
Figure BDA0002305144310000021
The half thickness of the rolled piece is equally divided into M sections, j is 1, 2 and 3
Figure BDA0002305144310000022
Δzi,j=Δzi,j-1=Δzi,j+1Δ z; wherein, Δ zi,jThe thickness-wise length, Δ z, of the corresponding element of the product node (i, j)i,j-1Rolled piece node (i, j-1)
Corresponding to the thickness-wise length of the cell, Δ zi,j+1-the thickness direction length of the product node (i, j +1) corresponding unit;
the nodes are divided into internal nodes, surface nodes, end nodes, core nodes and corner nodes according to positions;
(1) for an internal node, its sequence number may be represented as (i, j), i ═ 2, 3, 4.. N-1, j ═ 2, 3, 4.. M-1, and the node (i, j) abscissa is:
Figure BDA0002305144310000023
the ordinate is:
Figure BDA0002305144310000024
the expression of the temperature calculation model is as follows:
Figure BDA0002305144310000025
(2) for a surface node, its serial number may be represented as (i, M), i ═ 2, 3, 4.. N-1, and the node (i, M) abscissa is
Figure BDA0002305144310000026
The ordinate is
Figure BDA0002305144310000027
The expression of the temperature calculation model is as follows:
(i) if in the air cooling zone, then
Figure BDA0002305144310000028
(ii) If in the thermal radiation insulation area, then
Figure BDA0002305144310000029
(iii) If in the water cooling area, then
Figure BDA0002305144310000031
(iv) If in the roll gap contact zone, then
Figure BDA0002305144310000032
(3) For the end node, its serial number may be represented as (N, j), j 2, 3, 4
Figure BDA0002305144310000033
The ordinate is
Figure BDA0002305144310000034
The expression of the temperature calculation model is as follows:
(i) if in the air-cooling zone, water-cooling zone or roll gap contact zone
Figure BDA0002305144310000035
(ii) If in the thermal radiation insulation area, then
Figure BDA0002305144310000036
(4) For core nodes, where the serial number of the thickness core end node can be represented as (i,1), i-2, 3, 4.. N-1, with the abscissa being
Figure BDA0002305144310000037
The ordinate is
Figure BDA0002305144310000038
The number of the end face nodes of the width core can be represented by (1, j), j 2, 3, 4
Figure BDA0002305144310000039
The ordinate is
Figure BDA00023051443100000310
The expression of the temperature calculation model is as follows:
(i) if the node is a thick core node
Figure BDA0002305144310000041
(ii) If the core node is a wide core node
Figure BDA0002305144310000042
(5) For corner nodes, the top left corner node has a serial number of (1, M) and the abscissa has
Figure BDA0002305144310000043
The ordinate is
Figure BDA0002305144310000044
The node number at the upper right corner is (N, M), and the abscissa is
Figure BDA0002305144310000045
The ordinate is
Figure BDA0002305144310000046
The node sequence number of the lower left corner is (1,1), and the abscissa is
Figure BDA0002305144310000047
The ordinate is
Figure BDA0002305144310000048
The node number of the lower right corner is (N,1), and the abscissa is
Figure BDA0002305144310000049
The ordinate is
Figure BDA00023051443100000410
The expression of the temperature calculation model is as follows:
(i) if in the air cooling zone, then
For the top left corner node, satisfy
Figure BDA00023051443100000411
For the top right corner node, satisfy
Figure BDA00023051443100000412
For the lower left corner node, satisfy
Figure BDA00023051443100000413
For the lower right corner node, satisfy
Figure BDA0002305144310000051
(ii) If in the heat radiation heat preservation area, then for the upper left corner node, satisfy
Figure BDA0002305144310000052
For the top right corner node, satisfy
Figure BDA0002305144310000053
For the lower left corner node, satisfy
Figure BDA0002305144310000054
For the lower right corner node, satisfy
Figure BDA0002305144310000055
(iii) If the node is in the water cooling area, the requirement for the node at the upper left corner is met
Figure BDA0002305144310000056
For the top right corner node, satisfy
Figure BDA0002305144310000061
For the lower left corner node, satisfy
Figure BDA0002305144310000062
For the lower right corner node, satisfy
Figure BDA0002305144310000063
(iv) If the node is in the roll gap contact area, the condition of the upper left corner node is satisfied
Figure BDA0002305144310000064
For the top right corner node, satisfy
Figure BDA0002305144310000065
For the lower left corner node, satisfy
Figure BDA0002305144310000066
For the lower right corner node, satisfy
Figure BDA0002305144310000067
In the above formulas, hsEquivalent heat transfer coefficient between the rolled piece and the work rolls, in W/(mm)2×℃);hxWater cooling heat transfer coefficient in W/(mm)2X ° c); when descaling and cooling for high-pressure waterx=hhw(ii) a When cooling water for laminar cooling hx=hlwWherein h islwLaminar cooling water convection cooling heat transfer coefficient, hhw-high pressure water convective cooling heat transfer coefficient; h isaAir natural convection cooling heat transfer coefficient, unit W/(mm)2×℃);TrRoll surface temperature in units; t iswCooling water temperature, in units; t isa-ambient temperature in units; t isc-temperature of the heat-retaining cover in units; epsilonrThe thermal emissivity, i.e. blackness, epsilon, of the rolling stockr<1;σ0Radiation coefficient of absolute blackbody, σ0=5.67×10-6W/(mm2×K4) (ii) a c-specific heat capacity of rolled piece, unit J/(kg X DEG C); rho-rolled piece material density, unit kg/mm3(ii) a Lambda-rolled piece thermal conductivity (thermal conductivity), unit W/(mm × ° C); b, width of a rolled piece in unit mm; h-the section thickness of the rolled piece at the current moment, unit mm; Δ t-calculate time increment, unit S;
Figure BDA0002305144310000072
-temperature of node (i, j) at the current time, in units;
Figure BDA0002305144310000073
-temperature of node (i, j) at the last moment in time, in units;
Figure BDA0002305144310000074
-temperature of node (i-1, j) at the last moment in time, in units;
Figure BDA0002305144310000075
-temperature of node (i +1, j) at the last moment in time, in units;
Figure BDA0002305144310000076
-temperature of node (i, j-1) at the last moment in time, in units;
Figure BDA0002305144310000077
-temperature of node (i, j +1) at the last moment in time, in units;
Figure BDA0002305144310000078
heat generated by the heat source per unit volume of time, unit J/(mm)3Xs) when in the nip contact zone
Figure BDA0002305144310000079
When it is other heat exchange area
Figure BDA00023051443100000710
Wherein η -morph thermal equivalent is η -0.9, sigmasAverage yield strength of the rolled stock, hk-the thickness of the rolled stock at the previous moment, hk+1-the rolled piece thickness at the next moment.
According to the scheme, in the fourth step, the method for acquiring the initial temperature distribution of the rolled piece at the initial calculation position specifically comprises the following steps: in the width direction of the rolled piece, the temperature of each node is reduced from the middle width to the edge part in sequence, the temperature difference between adjacent nodes is the same, and the temperature difference between the middle width and the edge part is T0(1,j)-T0(N,j)=ΔT1(ii) a In the thickness direction of a rolled piece, the temperature of each node is increased from the surface to the middle thickness in sequence, the temperature difference between adjacent nodes is in an equal proportion relation, and the proportionality coefficient is gamma (gamma is more than or equal to 1), namely
Figure BDA00023051443100000711
And the temperature difference between the middle thickness position and the surface is T0(i,1)-T0(i,M)=ΔT2. Thus, the initial rolled product temperature profile can be expressed as:
when the gamma is greater than 1, the compound is,
Figure BDA00023051443100000712
when the gamma is 1, the gamma-,
Figure BDA0002305144310000081
according to the scheme, in the step five, the calculated rolled piece section moves along the rolling direction according to the running speed of the strip steel, and the specific method for calculating the temperature distribution of the rolled piece section from the initial calculation position to the final calculation position comprises the following steps:
(1) when the section of the rolled piece is in the thermal radiation heat preservation area, calculating the temperature distribution of the section of the rolled piece by adopting a rolled piece temperature model under the heat preservation condition, then judging whether the section of the rolled piece is still in the current thermal radiation heat preservation area at the next moment, if so, continuously calculating the temperature distribution of the section of the rolled piece under the thermal preservation condition at the next moment, and if not, switching to the next heat exchange area;
(2) when the section of the rolled piece is in a water cooling area, if the section of the rolled piece is in a high-pressure water descaling water cooling area, calculating the temperature distribution of the section of the rolled piece by adopting a rolled piece temperature model under the high-pressure water descaling condition, then judging whether the section of the rolled piece is still in the current high-pressure water descaling water cooling area at the next moment, if the section of the rolled piece is in the current high-pressure water descaling water cooling area, continuously calculating the temperature distribution of the section of the rolled piece under the high-pressure water descaling water cooling condition at the next moment, and if the section of the rolled piece; if the current laminar cooling water cooling area is not the laminar cooling water cooling area, the temperature distribution of the section of the rolled piece under the laminar cooling water cooling condition is continuously calculated, and if the current laminar cooling water cooling area is not the laminar cooling water cooling area, the next heat exchange area is switched to;
(3) when the section of the rolled piece is in a roll gap contact heat conduction area, judging whether the current rack is pressed down, if not, indicating that the roll gap of the rack does not exist, and transferring to the next heat exchange area; if the rolled piece is pressed down, calculating the temperature distribution of the rolled piece section by adopting a rolled piece temperature model of the roll gap contact heat conduction area, then judging whether the rolled piece section at the next moment is still in the current roll gap contact heat conduction area, if so, continuously calculating the temperature distribution of the rolled piece section under the roll gap contact condition at the next moment, and if not, switching to the next heat exchange area;
(4) when the section of the rolled piece is in the air cooling zone, calculating the temperature distribution of the section of the rolled piece by adopting a rolled piece temperature model under the air cooling condition, then judging whether the section of the rolled piece is still in the current air cooling zone at the next moment, if so, continuously calculating the temperature distribution of the section of the rolled piece under the air cooling condition at the next moment, and if not, switching to the next heat exchange zone;
(5) and when the section of the rolled piece is at the position of ending calculation, the calculation is finished.
The invention also provides a rolled piece temperature acquisition device for the plate and strip hot rolling production line, which comprises the following components:
the production line construction unit is used for constructing equipment arrangement of a hot rolling production line and acquiring process parameters;
the heat exchange area dividing unit is used for dividing the heat exchange area of the rolled piece according to the equipment arrangement on the production line;
the grid node dividing unit is used for dividing grid nodes of the width direction-thickness direction section of the rolled piece and establishing a coordinate system parameter and temperature calculation model of each node;
the initial temperature acquisition unit is used for acquiring initial temperature distribution of a rolled piece at an initial calculation position;
the heat exchange area judging and calculating unit is used for judging the temperature distribution of the cross section of the rolled piece from the initial calculating position to the final calculating position; and judging a heat exchange area where the cross section of the rolled piece is positioned when calculating once, and then calculating the temperature distribution of the cross section of the rolled piece at the position by adopting a rolled piece temperature model corresponding to the heat exchange area.
According to the scheme, the heat exchange area judging and calculating unit comprises:
the thermal radiation heat preservation area judgment module is used for judging whether the section of the rolled piece is in the current thermal radiation heat preservation area, if so, the temperature distribution of the section of the rolled piece under the heat preservation condition is calculated by adopting a rolled piece temperature calculation model of the thermal radiation heat preservation area, otherwise, the section of the rolled piece is shown to enter other heat exchange areas;
the laminar cooling water cooling area judging module is used for judging whether the section of the rolled piece is in the current spraying beam spraying area, if so, the rolled piece temperature calculation model of the laminar cooling water cooling area is adopted to calculate the temperature distribution of the section of the rolled piece, and otherwise, the section of the rolled piece is shown to enter other heat exchange areas;
the high-pressure water descaling water cooling area judging module is used for judging whether the section of the rolled piece is in the current high-pressure water descaling water cooling area, if so, the temperature distribution of the section of the rolled piece is calculated by adopting a rolled piece temperature calculation model of the high-pressure water descaling water cooling area, and otherwise, the section of the rolled piece is shown to enter other heat exchange areas;
the roll gap contact heat conduction area judgment module is used for judging whether the section of the rolled piece is in the current frame roll gap contact heat conduction area, if so, the rolled piece temperature distribution of the section of the rolled piece is calculated by adopting a rolled piece temperature calculation model of the roll gap contact heat conduction area, otherwise, the rolled piece section enters other heat exchange areas;
and the air cooling area judging module is used for judging whether the section of the rolled piece is in the current air cooling area, if so, calculating the temperature distribution of the section of the rolled piece by adopting a rolled piece temperature calculation model of the air cooling area, and otherwise, indicating that the section of the rolled piece enters other heat exchange areas.
The invention has the beneficial effects that:
(1) the method establishes a two-dimensional rolled piece temperature field difference equation, divides a rolled piece heat exchange area based on a given hot rolling production line arrangement, and further calculates the temperature of a hot rolled plate strip according to set process parameters (including rolling mill process parameters, rolled piece rolling parameters, laminar cooling process parameters, thermophysical parameters and medium parameters); the method has clear and definite principle, less assumed and simplified conditions, can accurately forecast the temperature field change of a rolled piece on a given hot rolling production line, has higher calculation precision than an analytical method and higher calculation speed than a finite element method, and can meet the online real-time requirements of engineering calculation, thereby ensuring the accurate control of the size precision, the shape quality and the tissue performance of the strip.
(2) The device disclosed by the invention adopts a modular design, can be very conveniently and flexibly suitable for calculating the temperature of rolled pieces in various hot rolling production lines (such as a thin slab continuous casting and rolling production line, a conventional hot continuous rolling production line, a coil rolling production line and the like), and greatly improves the universality of a model method.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
FIG. 2 is a schematic diagram for establishing a width-thickness coordinate system of a rolled piece in embodiment 1.
FIG. 3 is a schematic view of a one-half widthwise-thicknesswise cross-sectional discrete grid of the rolled piece of example 1.
Fig. 4 is a schematic diagram of the main equipment layout of the rolling section of the ESP line of example 2.
FIG. 5 is a rolled piece heat exchange area division schematic diagram of the ESP production line of embodiment 2.
Fig. 6 is a schematic of the laminar cooling strategy of example 2.
FIG. 7 is a graph showing the surface, core and thickness average temperature changes at the mid-width of the rolled stock calculated in example 2.
Fig. 8 is a schematic structural diagram of a rolled piece temperature obtaining device of a strip hot rolling production line provided in embodiment 3.
Wherein: 1. a first heat-insulating cover; 2. a three stand roughing mill train; 2-1, a first roughing mill; 2-2, a second roughing mill; 2-3, a third rough rolling mill; 3. a second heat-insulating cover; 4. an induction heating furnace; 5. a high-pressure water descaling device; 6. a five-stand finishing mill group; 6-1 a first finishing mill; 6-2, a second finishing mill; 6-3, a third finishing mill; 6-4, a fourth finishing mill; 6-5, a fifth finishing mill; 7. a laminar flow cooling device; 8. a coiler; 9. a thermal radiation heat preservation area; 10. an air cooling zone; 11. the roll gap contacts the heat conduction zone; 12. a high-pressure water descaling and water cooling area; 13. a laminar flow cooling zone; 13-1, a laminar cooling water cooling area; 13-2, laminar cooling air cooling area; 14. a spray beam; 14-1, a first spray beam; 14-2, a second spray beam; 14-3, a third spray beam.
Detailed Description
For a better understanding of the present invention, reference is made to the following detailed description taken in conjunction with the accompanying drawings.
The rolled piece temperature obtaining method for the plate and strip hot rolling production line shown in figure 1 comprises the following steps:
step one, constructing equipment arrangement of a plate strip hot rolling production line, and acquiring relevant process parameters; the equipment layout comprises the positions and distances of a rolling mill unit, a high-pressure water descaling device, a heat-insulating cover, a heating furnace, a laminar cooling device and a coiling machine, and the process parameters comprise rolling mill process parameters, rolled piece rolling parameters, laminar cooling process parameters, thermophysical parameters and medium parameters.
The technological parameters of the rolling mill comprise the roll body diameter, the roll body length, the roll neck diameter and the roll shifting displacement of the working roll of each frame; the rolling parameters of the rolled piece comprise the speed of the rolled piece, the width of the rolled piece, the inlet thickness and the outlet thickness of a roll gap of each frame and the yield strength of the rolled piece; the laminar cooling process parameters comprise the number of the jet beams, the jet length of the jet beams, the distance between the adjacent jet beams and the opening number and positions of the jet beams; the thermophysical parameters comprise the heat conductivity of the rolled piece, the specific heat capacity, the thermal radiance of the rolled piece, the contact heat exchange coefficient of the roll gap of each frame, the air cooling heat exchange coefficient of the rolled piece, the cooling water heat exchange coefficient of the rolled piece and the high-pressure water descaling heat exchange coefficient of the rolled piece; the medium parameters include cooling water temperature and ambient temperature.
And secondly, dividing a rolled piece heat exchange area according to equipment arrangement on the production line, wherein the heat exchange area comprises an air cooling area, a water cooling area, a roll gap contact heat conduction area and a heat radiation heat preservation area, and the water cooling area is divided into a high-pressure water descaling water cooling area and a laminar cooling water cooling area. The method specifically comprises the following steps: the area covered by the high-pressure water descaling device is a high-pressure water descaling water cooling area, the area covered by the heat preservation cover and the heating furnace is a heat radiation heat preservation area, the area covered by the laminar flow cooling device is a laminar flow cooling water cooling area, the contact area with a rolling mill roller is a roller gap contact heat conduction area, and the other areas only in contact with air are air cooling areas.
And step three, dividing the rolled piece width direction-thickness direction section grid nodes, and establishing a coordinate system parameter and temperature calculation model of each node.
The method for establishing the grid node division and temperature calculation model specifically comprises the following steps:
as shown in FIG. 2, a y-z rectangular coordinate system is established with the y-axis at the product intermediate thickness position and the z-axis at the product intermediate width position. As shown in fig. 3, the half-width-thickness cross section of the rolled piece is discretized into N × M grids, where the half-width of the rolled piece is equally divided into N segments, i is 1, 2, 3
Figure BDA0002305144310000101
The half thickness of the rolled piece is equally divided into M sections, j is 1, 2 and 3
Figure BDA0002305144310000102
Δzi,j=Δzi,j-1=Δzi,j+1Δ z; wherein, Δ zi,jThe thickness-wise length, Δ z, of the corresponding element of the product node (i, j)i,j-1The thickness-wise length, Δ z, of the corresponding element of the product node (i, j-1)i,j+1-the thickness direction length of the product node (i, j +1) corresponding unit;
the nodes are divided into internal nodes, surface nodes, end nodes, core nodes and corner nodes according to positions;
(1) for an internal node, its sequence number may be represented as (i, j), i ═ 2, 3, 4.. N-1, j ═ 2, 3, 4.. M-1, and the node (i, j) abscissa is:
Figure BDA0002305144310000111
the ordinate is:
Figure BDA0002305144310000112
the expression of the temperature calculation model is as follows:
Figure BDA0002305144310000113
(2) for a surface node, its serial number may be represented as (i, M), i ═ 2, 3, 4.. N-1, and the node (i, M) abscissa is
Figure BDA0002305144310000114
The ordinate is
Figure BDA0002305144310000115
The expression of the temperature calculation model is as follows:
(i) if in the air cooling zone, then
Figure BDA0002305144310000116
(ii) If in the thermal radiation insulation area, then
Figure BDA0002305144310000117
(iii) If in the water cooling area, then
Figure BDA0002305144310000118
(iv) If in the roll gap contact zone, then
Figure BDA0002305144310000119
(3) For the end node, its serial number may be represented as (N, j), j 2, 3, 4
Figure BDA0002305144310000121
The ordinate is
Figure BDA0002305144310000122
The expression of the temperature calculation model is as follows:
(i) if in the air-cooling zone, water-cooling zone or roll gap contact zone
Figure BDA0002305144310000123
(ii) If in the thermal radiation insulation area, then
Figure BDA0002305144310000124
(4) For core nodes, where the serial number of the thickness core end node can be represented as (i,1), i-2, 3, 4.. N-1, with the abscissa being
Figure BDA0002305144310000125
The ordinate is
Figure BDA0002305144310000126
The number of the end face nodes of the width core part canExpressed as (1, j), j 2, 3, 4.. M-1, with abscissa
Figure BDA0002305144310000127
The ordinate is
Figure BDA0002305144310000128
The expression of the temperature calculation model is as follows:
(i) if the node is a thick core node
Figure BDA0002305144310000129
(ii) If the core node is a wide core node
Figure BDA00023051443100001210
(5) For corner nodes, the top left corner node has a serial number of (1, M) and the abscissa has
Figure BDA0002305144310000131
The ordinate is
Figure BDA0002305144310000132
The node number at the upper right corner is (N, M), and the abscissa is
Figure BDA0002305144310000133
The ordinate is
Figure BDA0002305144310000134
The node sequence number of the lower left corner is (1,1), and the abscissa is
Figure BDA0002305144310000135
The ordinate is
Figure BDA0002305144310000136
The node number of the lower right corner is (N,1), and the abscissa is
Figure BDA0002305144310000137
The ordinate is
Figure BDA0002305144310000138
The expression of the temperature calculation model is as follows:
(i) if in the air cooling zone, then
For the top left corner node, satisfy
Figure BDA0002305144310000139
For the top right corner node, satisfy
Figure BDA00023051443100001310
For the lower left corner node, satisfy
Figure BDA00023051443100001311
For the lower right corner node, satisfy
Figure BDA00023051443100001312
(ii) If in the thermal radiation insulation area, then
For the top left corner node, satisfy
Figure BDA0002305144310000141
For the top right corner node, satisfy
Figure BDA0002305144310000142
For the lower left corner node, satisfy
Figure BDA0002305144310000143
For the lower right corner node, satisfy
Figure BDA0002305144310000144
(iii) If the node is in the water cooling area, the requirement for the node at the upper left corner is met
Figure BDA0002305144310000145
For the top right corner node, satisfy
Figure BDA0002305144310000146
For the lower left corner node, satisfy
Figure BDA0002305144310000151
For the lower right corner node, satisfy
Figure BDA0002305144310000152
(iv) If in the roll gap contact zone, then
For the top left corner node, satisfy
Figure BDA0002305144310000153
For the top right corner node, satisfy
Figure BDA0002305144310000154
For the lower left corner node, satisfy
Figure BDA0002305144310000155
For the lower right corner node, satisfy
Figure BDA0002305144310000156
In the above formulas, hsEquivalent heat transfer coefficient between the rolled piece and the work rolls, in W/(mm)2×℃);hxWater cooling heat transfer coefficient in W/(mm)2X ° c); when descaling and cooling for high-pressure waterx=hhw(ii) a When cooling water for laminar cooling hx=hlwWherein h islwLaminar cooling water convection cooling heat transfer coefficient, hhw-high pressure water convective cooling heat transfer coefficient; h isaAir natural convection cooling heat transfer coefficient, unit W/(mm)2×℃);TrRoll surface temperature in units; t iswCooling water temperature, in units; t isa-ambient temperature in units; t isc-temperature of the heat-retaining cover in units; epsilonrThe thermal emissivity, i.e. blackness, epsilon, of the rolling stockr<1;σ0Radiation coefficient of absolute blackbody, σ0=5.67×10-6W/(mm2×K4) (ii) a c-specific heat capacity of rolled piece, unit J/(kg X DEG C); rho-rolled piece material density, unit kg/mm3(ii) a Lambda-rolled piece thermal conductivity (thermal conductivity), unit W/(mm × ° C); b, width of a rolled piece in unit mm; h-the section thickness of the rolled piece at the current moment, unit mm; Δ t-calculate time increment, unit S;
Figure BDA0002305144310000161
-temperature of node (i, j) at the current time, in units;
Figure BDA0002305144310000162
-temperature of node (i, j) at the last moment in time, in units;
Figure BDA0002305144310000163
-temperature of node (i-1, j) at the last moment in time, in units;
Figure BDA0002305144310000164
temperature of node (i +1, j) at the last moment in deg.C;
Figure BDA0002305144310000165
-temperature of node (i, j-1) at the last moment in time, in units;
Figure BDA0002305144310000166
-temperature of node (i, j +1) at the last moment in time, in units;
Figure BDA0002305144310000167
heat generated by the heat source per unit volume of time, unit J/(mm)3Xs) when in the nip contact zone
Figure BDA0002305144310000168
When it is other heat exchange area
Figure BDA0002305144310000169
Wherein η -morph thermal equivalent is η -0.9, sigmasAverage yield strength of the rolled stock, hk-the thickness of the rolled stock at the previous moment, hk+1-the rolled piece thickness at the next moment.
And step four, acquiring initial temperature distribution of the rolled piece section at the initial calculation position. The method specifically comprises the following steps: in the width direction of the rolled piece, the temperature of each node is reduced from the middle width to the edge part in sequence, the temperature difference between adjacent nodes is the same, and the temperature difference between the middle width and the edge part is T0(1,j)-T0(N,j)=ΔT1(ii) a In the thickness direction of a rolled piece, the temperature of each node is increased from the surface to the middle thickness in sequence, the temperature difference between adjacent nodes is in an equal proportion relation, and the proportionality coefficient is gamma (gamma is more than or equal to 1), namely
Figure BDA00023051443100001610
And the temperature difference between the middle thickness position and the surface is T0(i,1)-T0(i,M)=ΔT2Thus, the initial rolled stock temperature profile can be expressed as:
when the gamma is greater than 1, the compound is,
Figure BDA00023051443100001611
when the gamma is 1, the gamma-,
Figure BDA00023051443100001612
step five, the calculated rolled piece section moves along the rolling direction according to the running speed of the strip steel, and the temperature distribution of the rolled piece section is calculated from the initial calculation position to the termination calculation position; and judging the heat exchange area where the cross section of the rolled piece is located when calculating once, and then calculating the temperature distribution of the cross section of the rolled piece at the position by adopting a rolled piece temperature model corresponding to the heat exchange area. The specific method comprises the following steps:
when the section of the rolled piece is in the thermal radiation heat preservation area, calculating the temperature distribution of the section of the rolled piece by adopting a rolled piece temperature model under the heat preservation condition, then judging whether the section of the rolled piece is still in the current thermal radiation heat preservation area at the next moment, if so, continuously calculating the temperature distribution of the section of the rolled piece under the thermal preservation condition at the next moment, and if not, switching to the next heat exchange area.
When the section of the rolled piece is in a water cooling area, if the rolled piece is in a high-pressure water descaling water cooling area, calculating the temperature distribution of the section of the rolled piece by adopting a rolled piece temperature model under the high-pressure water descaling condition, then judging whether the section of the rolled piece is still in the current high-pressure water descaling water cooling area at the next moment, if the section of the rolled piece is in the current high-pressure water descaling water cooling area, continuously calculating the temperature distribution of the section of the rolled piece under the high-pressure water descaling water cooling condition at the next moment, and if the section of the rolled piece is not; and if the temperature distribution is in the laminar cooling water cooling area, calculating the temperature distribution of the section of the rolled piece by using a rolled piece temperature model under the laminar cooling water cooling condition, then judging whether the section of the rolled piece is still in the current laminar cooling water cooling area at the next moment, if so, continuously calculating the temperature distribution of the section of the rolled piece under the laminar cooling water cooling condition at the next moment, and if not, switching to the next heat exchange area.
When the section of the rolled piece is in a roll gap contact heat conduction area, judging whether the current rack is pressed down, if not, indicating that the roll gap of the rack does not exist, and transferring to the next heat exchange area; if the rolled piece is pressed down, calculating the temperature distribution of the rolled piece section by adopting a rolled piece temperature model of the roll gap contact heat conduction area, then judging whether the rolled piece section at the next moment is still in the current roll gap contact heat conduction area, if so, continuously calculating the temperature distribution of the rolled piece section under the roll gap contact condition at the next moment, and if not, switching to the next heat exchange area;
when the section of the rolled piece is in an air cooling zone, calculating the temperature distribution of the section of the rolled piece by adopting a rolled piece temperature model under the air cooling condition, then judging whether the section of the rolled piece is still in the current air cooling zone at the next moment, if so, continuously calculating the temperature distribution of the section of the rolled piece under the air cooling condition at the next moment, and if not, switching to the next heat exchange zone;
and when the rolled piece section is at the calculation termination position, the calculation is finished.
Example two
In this embodiment, an ESP production line is taken as an example, and based on the method described in the first embodiment, the temperature of a rolled piece in the production process of the ESP production line is calculated and compared with a measured temperature value on site, so as to further illustrate the universality and accuracy of the method of the present invention.
The ESP process is a thin slab continuous casting and rolling process, the main equipment arrangement of the rolling section of the ESP process is shown in figure 4, the ESP process sequentially comprises a first heat-preservation cover 1, a three-stand roughing mill group 2, a second heat-preservation cover 3, an induction heating furnace 4, a high-pressure water descaling device 5, a five-stand finishing mill group 6, a laminar cooling device 7 and a coiler 8, and the distance between the main equipment is as follows: l1-3 m, L2-5 m, L3-9.5 m, L4-14 m, L5-26.5 m, L6-37 m, L7-45.5 m, L8-56 m, L9-57 m, L10-58 m, L11-63 m, L12-67.5 m, L13-72 m, L14-76.5 m, L15-81 m, L16-92 m, L17-130 m, L18-150 m.
Table 1, table 2, table 3 and table 4 show the rolling mill process parameters, rolled piece rolling parameters, laminar cooling process parameters and thermophysical parameters, respectively, in addition, the cooling water temperature is 30 ℃, the ambient temperature is 30 ℃, and the induction heating increases the temperature of the rolled piece by 240 ℃ as a whole.
TABLE 1ESP production line Rolling mill technological parameters
Figure BDA0002305144310000181
TABLE 2 rolled stock Rolling parameters
Figure BDA0002305144310000182
TABLE 3 laminar flow Cooling Process parameters
Figure BDA0002305144310000191
TABLE 4 thermophysical parameters
Figure BDA0002305144310000192
FIG. 5 is a schematic diagram of the division of the heat exchange area of the rolled piece in the ESP production line, sequentially including a heat radiation insulation area → an air cooling area → a roller gap contact heat transfer area → an air cooling area → a heat radiation insulation area → an air cooling area → an induction heating temperature of 300 ℃ → an air cooling area → a high pressure descaling water cooling area → an air cooling area → a roller gap contact heat transfer area → an air cooling area → a laminar flow cooling area → an air cooling area, as shown in FIG. 6, the front section cooling strategy is adopted for the middle section flow cooling in this embodiment, and the front section 3 cooling water spray beams (first spray beam 14-1, first spray beam, second spray beam, and third spray beam 14-1, second spray beam, The second injection beam 14-2 and the third injection beam 14-3), and therefore the laminar flow cooling heat exchange area is specifically the laminar flow cooling water cooling area → the air cooling area → the laminar flow cooling water cooling area → the air cooling area in this embodiment.
The half width-thickness cross section of the rolled piece is scattered into 50 multiplied by 40 grids, wherein the half width of the rolled piece is equally divided into 50 sections, i is 1, 2 and 3Has a thickness of
Figure BDA0002305144310000201
h is the cross-sectional thickness of the rolled piece at the current calculation moment, and needs to be calculated in real time in the calculation process, and it can be known from table 2 that h is not less than 2 and not more than 80 mm;
the initial temperature of the rolled piece at the initial calculation position is as follows: surface initial temperature T at intermediate width0(1,40)=1100℃,T0(1,1)=1360℃,T0(50,40) 1070 ℃, and the proportionality coefficient gamma 1.053, so the initial temperature at other nodes of the rolled stock is:
Figure BDA0002305144310000202
from the initial calculation position to the final calculation position, the heat exchange area where the rolled piece section is located is judged at each calculation, then the temperature distribution of the rolled piece section at the position is calculated by adopting a corresponding rolled piece temperature model, and the calculation result is shown in fig. 7. Fig. 7 shows the surface temperature, the core temperature and the thickness direction average temperature variation of the ESP production line at the intermediate width of the rolled piece under the above working conditions, and also shows the measured value of the surface temperature at the intermediate width of the rolled piece at the relevant position of the ESP production line (shown by the black dots in fig. 7).
EXAMPLE III
A rolled piece temperature obtaining apparatus for a strip hot rolling line as shown in fig. 8 comprises:
the production line construction unit is used for constructing equipment arrangement of a hot rolling production line and acquiring relevant process parameters;
the heat exchange area dividing unit is used for dividing the heat exchange area of the rolled piece according to the equipment arrangement on the production line;
the grid node dividing unit is used for dividing grid nodes of the width direction-thickness direction section of the rolled piece and establishing a coordinate system parameter and temperature calculation model of each node;
the initial temperature acquisition unit is used for acquiring initial temperature distribution of a rolled piece at an initial calculation position;
and the heat exchange area judging and calculating unit is used for judging the heat exchange area where the section of the rolled piece is located from the initial calculating position to the final calculating position every time of calculation, and then calculating the temperature distribution of the section of the rolled piece at the position by adopting a corresponding rolled piece temperature model.
In the present invention, the heat exchange area determining and calculating unit includes:
the thermal radiation heat preservation area judgment module is used for judging whether the section of the rolled piece is in the current thermal radiation heat preservation area, if so, the temperature distribution of the section of the rolled piece under the heat preservation condition is calculated by adopting a rolled piece temperature calculation model of the thermal radiation heat preservation area, otherwise, the section of the rolled piece is shown to enter other heat exchange areas;
the laminar cooling water cooling area judging module is used for judging whether the section of the rolled piece is in the current spraying beam spraying area, if so, the rolled piece temperature calculation model of the laminar cooling water cooling area is adopted to calculate the temperature distribution of the section of the rolled piece, and otherwise, the section of the rolled piece is shown to enter other heat exchange areas;
the high-pressure water descaling water cooling area judging module is used for judging whether the section of the rolled piece is in the current high-pressure water descaling water cooling area, if so, the temperature distribution of the section of the rolled piece is calculated by adopting a rolled piece temperature calculation model of the high-pressure water descaling water cooling area, and otherwise, the section of the rolled piece is shown to enter other heat exchange areas;
the roll gap contact heat conduction area judgment module is used for judging whether the section of the rolled piece is in the current frame roll gap contact heat conduction area, if so, the rolled piece temperature distribution of the section of the rolled piece is calculated by adopting a rolled piece temperature calculation model of the roll gap contact heat conduction area, otherwise, the rolled piece section enters other heat exchange areas;
and the air cooling area judging module is used for judging whether the section of the rolled piece is in the current air cooling area, if so, calculating the temperature distribution of the section of the rolled piece by adopting a rolled piece temperature calculation model of the air cooling area, and otherwise, indicating that the section of the rolled piece enters other heat exchange areas.
It should be noted that the above-mentioned embodiments are only preferred embodiments of the present invention, and the present invention is not limited thereto, and although the present invention has been described in detail with reference to the embodiments, it will be apparent to those skilled in the art that modifications can be made to the technical solutions described in the above-mentioned embodiments, or equivalent substitutions of some technical features, but any modifications, equivalents, improvements and the like within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims (8)

1. A rolled piece temperature obtaining method for a plate and strip hot rolling production line is characterized by comprising the following steps:
step one, constructing equipment arrangement of a plate strip hot rolling production line, and acquiring relevant process parameters;
dividing a rolled piece heat exchange area according to equipment arrangement on the production line, wherein the heat exchange area comprises an air cooling area, a water cooling area, a roll gap contact heat conduction area and a heat radiation heat preservation area, and the water cooling area is divided into a high-pressure water descaling water cooling area and a laminar cooling water cooling area;
dividing rolled piece width direction-thickness direction section grid nodes, and establishing a coordinate system parameter and temperature calculation model of each node;
acquiring initial temperature distribution of the section of the rolled piece at the initial calculation position;
and step five, the calculated rolled piece section moves along the rolling direction according to the running speed of the strip steel, the heat exchange area where the rolled piece section is located is judged from the initial calculation position to the final calculation position every time of calculation, and then the temperature distribution of the rolled piece section at the position is calculated by adopting a corresponding rolled piece temperature model.
2. The rolled piece temperature acquisition method for the strip hot rolling production line according to claim 1, wherein the equipment layout of the strip hot rolling production line comprises positions and distances of a rolling mill unit, a high-pressure water descaling device, a heat preservation cover, a heating furnace, a laminar cooling device and a coiler, and the process parameters comprise rolling mill process parameters, rolled piece rolling parameters, laminar cooling process parameters, thermophysical parameters and medium parameters.
3. The rolled piece temperature acquisition method for the plate and strip hot rolling production line as claimed in claim 2, wherein the rolling mill process parameters include the roll body diameter, the roll body length, the roll neck diameter, the work roll shifting displacement of each stand work roll; the rolling parameters of the rolled piece comprise the speed of the rolled piece, the width of the rolled piece, the inlet thickness and the outlet thickness of a roll gap of each frame and the yield strength of the rolled piece; the laminar cooling process parameters comprise the number of the jet beams, the jet length of the jet beams, the distance between the adjacent jet beams and the opening number and positions of the jet beams; the thermophysical parameters comprise the heat conductivity of the rolled piece, the specific heat capacity, the thermal radiance of the rolled piece, the contact heat exchange coefficient of the roll gap of each frame, the air cooling heat exchange coefficient of the rolled piece, the cooling water heat exchange coefficient of the rolled piece and the high-pressure water descaling heat exchange coefficient of the rolled piece; the medium parameters include cooling water temperature and ambient temperature.
4. The rolled piece temperature obtaining method for the plate and strip hot rolling production line according to claim 1, characterized in that in step three, grid nodes of a rolled piece width direction-thickness direction section are divided, and a coordinate system position and temperature calculation model of each node is established, specifically: establishing a y-z rectangular coordinate system, wherein a y axis is positioned at the position of the middle thickness of the rolled piece, and a z axis is positioned at the position of the middle width of the rolled piece; dispersing a half width-thickness cross section of a rolled piece into N multiplied by M grids, wherein the half width of the rolled piece is equally divided into N sections, i is 1, 2 and 3
Figure FDA0002305144300000011
The half thickness of the rolled piece is equally divided into M sections, j is 1, 2 and 3
Figure FDA0002305144300000012
Δzi,j=Δzi,j-1=Δzi,j+1Δ z; wherein, Δ zi,jThe thickness-wise length, Δ z, of the corresponding element of the product node (i, j)i,j-1The thickness-wise length, Δ z, of the corresponding element of the product node (i, j-1)i,j+1-the thickness direction length of the product node (i, j +1) corresponding unit;
the nodes are divided into internal nodes, surface nodes, end nodes, core nodes and corner nodes according to positions;
(1) for internal nodes, the order thereofThe number may be represented as (i, j), i 2, 3, 4.. N-1, j 2, 3, 4.. M-1, with the node (i, j) having the abscissa:
Figure FDA0002305144300000013
the ordinate is:
Figure FDA0002305144300000014
the expression of the temperature calculation model is as follows:
Figure FDA0002305144300000021
(2) for a surface node, its serial number may be represented as (i, M), i ═ 2, 3, 4.. N-1, and the node (i, M) abscissa is
Figure FDA0002305144300000022
The ordinate is
Figure FDA0002305144300000023
The expression of the temperature calculation model is as follows:
(i) if in the air cooling zone, then
Figure FDA0002305144300000024
(ii) If in the thermal radiation insulation area, then
Figure FDA0002305144300000025
(iii) If in the water cooling area, then
Figure FDA0002305144300000026
(iv) If in the roll gap contact zone, then
Figure FDA0002305144300000027
(3) For the end node, its serial number may be represented as (N, j), j 2, 3, 4
Figure FDA0002305144300000028
The ordinate is
Figure FDA0002305144300000029
The expression of the temperature calculation model is as follows:
(i) if in the air-cooling zone, water-cooling zone or roll gap contact zone
Figure FDA0002305144300000031
(ii) If in the thermal radiation insulation area, then
Figure FDA0002305144300000032
(4) For core nodes, where the serial number of the thickness core end node can be represented as (i,1), i-2, 3, 4.. N-1, with the abscissa being
Figure FDA0002305144300000033
The ordinate is
Figure FDA0002305144300000034
The number of the end face nodes of the width core can be represented by (1, j), j 2, 3, 4
Figure FDA0002305144300000035
The ordinate is
Figure FDA0002305144300000036
The expression of the temperature calculation model is as follows:
(i) if the node is a thick core node
Figure FDA0002305144300000037
(ii) If the core node is a wide core node
Figure FDA0002305144300000038
(5) For corner nodes, the top left corner node has a serial number of (1, M) and the abscissa has
Figure FDA0002305144300000039
The ordinate is
Figure FDA00023051443000000310
The node number at the upper right corner is (N, M), and the abscissa is
Figure FDA00023051443000000311
The ordinate is
Figure FDA00023051443000000312
The node sequence number of the lower left corner is (1,1), and the abscissa is
Figure FDA00023051443000000313
The ordinate is
Figure FDA00023051443000000314
The node number of the lower right corner is (N,1), and the abscissa is
Figure FDA0002305144300000041
The ordinate is
Figure FDA0002305144300000042
The expression of the temperature calculation model is as follows:
(i) if in the air cooling zone, then
For the top left corner node, satisfy
Figure FDA0002305144300000043
For the top right corner node, satisfy
Figure FDA0002305144300000044
For the lower left corner node, satisfy
Figure FDA0002305144300000045
For the lower right corner node, satisfy
Figure FDA0002305144300000046
(ii) If in the thermal radiation insulation area, then
For the top left corner node, satisfy
Figure FDA0002305144300000047
For the top right corner node, satisfy
Figure FDA0002305144300000051
For the lower left corner node, satisfy
Figure FDA0002305144300000052
For the lower right corner node, satisfy
Figure FDA0002305144300000053
(iii) If in the water cooling area, then
For the top left corner node, satisfy
Figure FDA0002305144300000054
For the top right corner node, satisfy
Figure FDA0002305144300000055
For the lower left corner node, satisfy
Figure FDA0002305144300000056
For the lower right corner node, satisfy
Figure FDA0002305144300000057
(iv) If in the roll gap contact zone, then
For the top left corner node, satisfy
Figure FDA0002305144300000062
For the top right corner node, satisfy
Figure FDA0002305144300000063
For the lower left corner node, satisfy
Figure FDA0002305144300000064
For the lower right corner node, satisfy
Figure FDA0002305144300000065
In the above formulas, hsEquivalent heat transfer coefficient between the rolled piece and the work rolls, in W/(mm)2×℃);hxWater cooling heat transfer coefficient in W/(mm)2X ° c); when descaling and cooling for high-pressure waterx=hhw(ii) a When cooling water for laminar cooling hx=hlwWherein h islwLaminar cooling water convection cooling heat transfer coefficient, hhw-high pressure water convective cooling heat transfer coefficient; h isaAir natural convection cooling heat transfer coefficient, unit W/(mm)2×℃);TrRoll surface temperature in units; t iswCooling water temperature, in units; t isa-ambient temperature in units; t isc-temperature of the heat-retaining cover in units; epsilonrThe thermal emissivity, i.e. blackness, epsilon, of the rolling stockr<1;σ0Radiation coefficient of absolute blackbody, σ0=5.67×10-6W/(mm2×K4) (ii) a c-specific heat capacity of rolled piece, unit J/(kg X DEG C); rho-rolled piece material density, unit kg/mm3(ii) a Lambda-rolled piece thermal conductivity (thermal conductivity), unit W/(mm × ° C); b, width of a rolled piece in unit mm; h-the section thickness of the rolled piece at the current moment, unit mm; Δ t-calculate time increment, unit S;
Figure FDA0002305144300000066
-temperature of node (i, j) at the current time, in units;
Figure FDA0002305144300000067
-temperature of node (i, j) at the last moment in time, in units;
Figure FDA0002305144300000068
-temperature of node (i-1, j) at the last moment in time, in units;
Figure FDA0002305144300000069
-temperature of node (i +1, j) at the last moment in time, in units;
Figure FDA0002305144300000071
-temperature of node (i, j-1) at the last moment in time, in units;
Figure FDA0002305144300000072
-temperature of node (i, j +1) at the last moment in time, in units;
Figure FDA0002305144300000073
heat generated by the heat source per unit volume of time, unit J/(mm)3Xs) when in the nip contact zone
Figure FDA0002305144300000074
When it is other heat exchange area
Figure FDA0002305144300000075
Wherein η -morph thermal equivalent is η -0.9, sigmasAverage yield strength of the rolled stock, hk-the thickness of the rolled stock at the previous moment, hk+1-the rolled piece thickness at the next moment.
5. The rolled piece temperature obtaining method for the plate and strip hot rolling production line according to claim 1, wherein in the fourth step, the method for obtaining the initial temperature distribution of the rolled piece at the initial calculation position specifically comprises the following steps: in the width direction of the rolled piece, the temperature of each node is reduced from the middle width to the edge part in sequence, the temperature difference between adjacent nodes is the same, and the temperature difference between the middle width and the edge part is T0(1,j)-T0(N,j)=ΔT1(ii) a In the thickness direction of a rolled piece, the temperature of each node is increased from the surface to the middle thickness in sequence, the temperature difference between adjacent nodes is in an equal proportion relation, and the proportionality coefficient is gamma (gamma is more than or equal to 1), namely
Figure FDA0002305144300000076
And the temperature difference between the middle thickness position and the surface is T0(i,1)-T0(i,M)=ΔT2. Thus, the initial rolled product temperature profile can be expressed as:
when the gamma is greater than 1, the compound is,
Figure FDA0002305144300000077
when the gamma is 1, the gamma-,
Figure FDA0002305144300000078
6. the method for acquiring the temperature of the rolled piece in the plate and strip hot rolling production line according to claim 1, wherein in the fifth step, the calculated rolled piece section moves along the rolling direction according to the running speed of the strip steel, and the specific method for calculating the temperature distribution of the rolled piece section from the initial calculation position to the final calculation position comprises the following steps:
(1) when the section of the rolled piece is in the thermal radiation heat preservation area, calculating the temperature distribution of the section of the rolled piece by adopting a rolled piece temperature model under the heat preservation condition, then judging whether the section of the rolled piece is still in the current thermal radiation heat preservation area at the next moment, if so, continuously calculating the temperature distribution of the section of the rolled piece under the thermal preservation condition at the next moment, and if not, switching to the next heat exchange area;
(2) when the section of the rolled piece is in a water cooling area, if the section of the rolled piece is in a high-pressure water descaling water cooling area, calculating the temperature distribution of the section of the rolled piece by adopting a rolled piece temperature model under the high-pressure water descaling condition, then judging whether the section of the rolled piece is still in the current high-pressure water descaling water cooling area at the next moment, if the section of the rolled piece is in the current high-pressure water descaling water cooling area, continuously calculating the temperature distribution of the section of the rolled piece under the high-pressure water descaling water cooling condition at the next moment, and if the section of the rolled piece; if the current laminar cooling water cooling area is not the laminar cooling water cooling area, the temperature distribution of the section of the rolled piece under the laminar cooling water cooling condition is continuously calculated, and if the current laminar cooling water cooling area is not the laminar cooling water cooling area, the next heat exchange area is switched to;
(3) when the section of the rolled piece is in a roll gap contact heat conduction area, judging whether the current rack is pressed down, if not, indicating that the roll gap of the rack does not exist, and transferring to the next heat exchange area; if the rolled piece is pressed down, calculating the temperature distribution of the rolled piece section by adopting a rolled piece temperature model of the roll gap contact heat conduction area, then judging whether the rolled piece section at the next moment is still in the current roll gap contact heat conduction area, if so, continuously calculating the temperature distribution of the rolled piece section under the roll gap contact condition at the next moment, and if not, switching to the next heat exchange area;
(4) when the section of the rolled piece is in the air cooling zone, calculating the temperature distribution of the section of the rolled piece by adopting a rolled piece temperature model under the air cooling condition, then judging whether the section of the rolled piece is still in the current air cooling zone at the next moment, if so, continuously calculating the temperature distribution of the section of the rolled piece under the air cooling condition at the next moment, and if not, switching to the next heat exchange zone;
(5) and when the section of the rolled piece is at the position of ending calculation, the calculation is finished.
7. The utility model provides a strip hot rolling production line rolled piece temperature acquisition device which characterized in that includes:
the production line construction unit is used for constructing equipment arrangement of a hot rolling production line and acquiring process parameters;
the heat exchange area dividing unit is used for dividing the heat exchange area of the rolled piece according to the equipment arrangement on the production line;
the grid node dividing unit is used for dividing grid nodes of the width direction-thickness direction section of the rolled piece and establishing a coordinate system parameter and temperature calculation model of each node;
the initial temperature acquisition unit is used for acquiring initial temperature distribution of a rolled piece at an initial calculation position;
the heat exchange area judging and calculating unit is used for judging the temperature distribution of the cross section of the rolled piece from the initial calculating position to the final calculating position; and judging a heat exchange area where the cross section of the rolled piece is positioned when calculating once, and then calculating the temperature distribution of the cross section of the rolled piece at the position by adopting a rolled piece temperature model corresponding to the heat exchange area.
8. The rolled piece temperature acquisition device for the plate and strip hot rolling production line of claim 7, wherein the heat exchange area judgment and calculation unit comprises:
the thermal radiation heat preservation area judgment module is used for judging whether the section of the rolled piece is in the current thermal radiation heat preservation area, if so, the temperature distribution of the section of the rolled piece under the heat preservation condition is calculated by adopting a rolled piece temperature calculation model of the thermal radiation heat preservation area, otherwise, the section of the rolled piece is shown to enter other heat exchange areas;
the laminar cooling water cooling area judging module is used for judging whether the section of the rolled piece is in the current spraying beam spraying area, if so, the rolled piece temperature calculation model of the laminar cooling water cooling area is adopted to calculate the temperature distribution of the section of the rolled piece, and otherwise, the section of the rolled piece is shown to enter other heat exchange areas;
the high-pressure water descaling water cooling area judging module is used for judging whether the section of the rolled piece is in the current high-pressure water descaling water cooling area, if so, the temperature distribution of the section of the rolled piece is calculated by adopting a rolled piece temperature calculation model of the high-pressure water descaling water cooling area, and otherwise, the section of the rolled piece is shown to enter other heat exchange areas;
the roll gap contact heat conduction area judgment module is used for judging whether the section of the rolled piece is in the current frame roll gap contact heat conduction area, if so, the rolled piece temperature distribution of the section of the rolled piece is calculated by adopting a rolled piece temperature calculation model of the roll gap contact heat conduction area, otherwise, the rolled piece section enters other heat exchange areas;
and the air cooling area judging module is used for judging whether the section of the rolled piece is in the current air cooling area, if so, calculating the temperature distribution of the section of the rolled piece by adopting a rolled piece temperature calculation model of the air cooling area, and otherwise, indicating that the section of the rolled piece enters other heat exchange areas.
CN201911236967.9A 2019-12-05 2019-12-05 Method and device for acquiring temperature of rolled piece in strip hot rolling production line Active CN111079275B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911236967.9A CN111079275B (en) 2019-12-05 2019-12-05 Method and device for acquiring temperature of rolled piece in strip hot rolling production line

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911236967.9A CN111079275B (en) 2019-12-05 2019-12-05 Method and device for acquiring temperature of rolled piece in strip hot rolling production line

Publications (2)

Publication Number Publication Date
CN111079275A true CN111079275A (en) 2020-04-28
CN111079275B CN111079275B (en) 2023-04-28

Family

ID=70313235

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911236967.9A Active CN111079275B (en) 2019-12-05 2019-12-05 Method and device for acquiring temperature of rolled piece in strip hot rolling production line

Country Status (1)

Country Link
CN (1) CN111079275B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113899471A (en) * 2021-10-27 2022-01-07 中冶南方工程技术有限公司 Method for acquiring temperature distribution of rolled piece in induction heating system
CN114178325A (en) * 2021-10-29 2022-03-15 中冶南方工程技术有限公司 Method for obtaining cooling water flow of hot-rolled carbon steel laminar cooling injection header and temperature calculation method
CN114178323A (en) * 2021-10-29 2022-03-15 中冶南方工程技术有限公司 Method for obtaining cooling water flow of hot-rolled dual-phase steel laminar flow cooling injection header
CN114178324A (en) * 2021-10-29 2022-03-15 中冶南方工程技术有限公司 Method for obtaining cooling water flow of hot-rolled alloy steel laminar cooling injection header and correction method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002282919A (en) * 2001-03-26 2002-10-02 Nippon Steel Corp Learning calculation method for roll profile in plate rolling
CN101221416A (en) * 2007-12-28 2008-07-16 东北大学 Finite element method for on-line board temperature calculation in course of hot rolling
CN107066737A (en) * 2017-04-14 2017-08-18 北京科技大学 A kind of two-dimentional staggered difference method for predicting hot rolling process plate belt temperature field

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002282919A (en) * 2001-03-26 2002-10-02 Nippon Steel Corp Learning calculation method for roll profile in plate rolling
CN101221416A (en) * 2007-12-28 2008-07-16 东北大学 Finite element method for on-line board temperature calculation in course of hot rolling
CN107066737A (en) * 2017-04-14 2017-08-18 北京科技大学 A kind of two-dimentional staggered difference method for predicting hot rolling process plate belt temperature field

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
杨金光;孙丽荣;刘华强;杨荃;王连生;: "基于PSO和DE优化算法的热轧工作辊热辊形的研究" *
汪娟;曹燕;郝振宇;黄贞益;: "热轧工作辊温度场有限元分析" *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113899471A (en) * 2021-10-27 2022-01-07 中冶南方工程技术有限公司 Method for acquiring temperature distribution of rolled piece in induction heating system
CN114178325A (en) * 2021-10-29 2022-03-15 中冶南方工程技术有限公司 Method for obtaining cooling water flow of hot-rolled carbon steel laminar cooling injection header and temperature calculation method
CN114178323A (en) * 2021-10-29 2022-03-15 中冶南方工程技术有限公司 Method for obtaining cooling water flow of hot-rolled dual-phase steel laminar flow cooling injection header
CN114178324A (en) * 2021-10-29 2022-03-15 中冶南方工程技术有限公司 Method for obtaining cooling water flow of hot-rolled alloy steel laminar cooling injection header and correction method
CN114178324B (en) * 2021-10-29 2023-06-23 中冶南方工程技术有限公司 Cooling water flow obtaining method and correction method for hot-rolled alloy steel laminar cooling jet header

Also Published As

Publication number Publication date
CN111079275B (en) 2023-04-28

Similar Documents

Publication Publication Date Title
CN111079275A (en) Rolled piece temperature obtaining method and device for strip hot rolling production line
CN107066737B (en) A kind of two-dimentional staggered difference method for predicting hot rolling process plate belt temperature field
CN101456038B (en) Plate-belt temperature monitoring method during hot-rolled strip steel stream cooling process
CN102941232B (en) Method for controlling hot continuous rolling finish rolling process
CN103225017B (en) Rod and wire billet heating furnace model control method and apparatus
CN101780478B (en) Method and device for controlling strip shape and precision of hot rolling steel strips and plates
CN103286141A (en) Automatic control method of width of hot continuous rolling and finish rolling strip steel
CN103028615A (en) Method for predicting temperature evolution in hot continuous rolling process of strip steel
CN106984652A (en) The method that finishing stand sideslip is controlled according to breakdown bar camber
CN109013717B (en) A kind of hot continuous rolling centre base center portion temperature computation method
CN104298884A (en) Finite element and finite difference coupling method for fast calculating rolled piece section temperature
CN105414205B (en) PLC-based online predication method for temperatures of steel plates
CN104324949B (en) A kind of acquisition methods of roughing vertical roll passage edger roll opening degree
CN102601128B (en) Temperature difference control method for rolling pieces for steckel mill
CN108941196A (en) A method of with 200 big cross section square billet continuous rolling small-sized bars
CN103559334A (en) Modeling method and system for temperature field in laminar cooling
CN105855297B (en) A kind of control method for improving first piece of non-orientation silicon steel head thickness precision of hot rolling
CN103990653B (en) Entry temperature at finishing accuracy at target method of assuring
CN106825060B (en) The coding method of Heat Working Rolls chilling spray amount control system and spray flux record
CN101811143B (en) Temperature control method for heat-preservation casting blank to be rolled
CN201033333Y (en) Strip steel coiling temperature control device
KR100510841B1 (en) Method for designing the second optimum cooling pattern of continuous slab casting
CN105149361A (en) Self-learning method for heat transfer coefficients of belt materials of hot continuous rolling mill
KR100568346B1 (en) A method for controlling the bar temperature by predicting the rougher delivery temperature
Peng et al. Mathematical modeling and simulated analysis of metal flow behavior during the FGC of ESP rolling process

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant