CN113899471A - Method for acquiring temperature distribution of rolled piece in induction heating system - Google Patents
Method for acquiring temperature distribution of rolled piece in induction heating system Download PDFInfo
- Publication number
- CN113899471A CN113899471A CN202111256027.3A CN202111256027A CN113899471A CN 113899471 A CN113899471 A CN 113899471A CN 202111256027 A CN202111256027 A CN 202111256027A CN 113899471 A CN113899471 A CN 113899471A
- Authority
- CN
- China
- Prior art keywords
- node
- rolled piece
- temperature
- induction heating
- rolled
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K13/00—Thermometers specially adapted for specific purposes
Abstract
The invention discloses a method for acquiring temperature distribution of a rolled piece in an induction heating system, which comprises the following steps: constructing a process arrangement of the induction heating system, and acquiring process parameters; 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; calculating the temperature distribution of the cross section of the rolled piece at the position from the initial calculation position until the final calculation position, judging the heat exchange area where the cross section of the rolled piece is located when calculating once, and calculating the temperature distribution of the cross section of the rolled piece at the position by adopting a corresponding temperature calculation model. The invention has the beneficial effects that: the method has the advantages that the assumed and simplified conditions are few, the temperature field change of the rolled piece in the induction heating process can be accurately forecasted, the calculation precision is higher than that of an analytic method, the calculation speed is higher than that of a finite element method, and the requirements of high precision and high efficiency of online calculation can be met.
Description
Technical Field
The invention relates to the technical field of steel rolling, in particular to a method for acquiring temperature distribution of a rolled piece in an induction heating system.
Background
In the field of strip steel production, induction heating devices are increasingly used with the advantages of small floor area, high heating efficiency, no greenhouse gas emission and the like. However, in the field use process, the temperature prediction error of the rolled piece of the existing induction heating system is large, and the temperature prediction error cannot be used for feedforward control of the temperature of the rolled piece in the induction heating process, so that the control precision of the temperature of the rolled piece is low. The rolled piece temperature model of the existing induction heating system is usually calculated by adopting an analytical formula, and the analytical model does not consider the complexity of boundary conditions and field working conditions, so that more assumptions and simplification conditions are adopted, and the calculation precision is low. In addition, the calculation result of the analytical model is only the average temperature of the rolled piece and cannot be used for calculating the temperature distribution of the rolled piece in the induction heating process, so that the microstructure prediction of the rolled piece in the heating process is not facilitated. Therefore, it is necessary to further develop a calculation method of the temperature distribution of a rolled material in an induction heating system suitable for on-site production.
Disclosure of Invention
The invention aims to provide a method for acquiring the temperature distribution of a rolled piece in an induction heating system aiming at the defects of the prior art, and solves the problems that the temperature forecast error of the existing induction heating rolled piece is large, and the temperature distribution of the rolled piece in the induction heating process cannot be calculated quickly and conveniently.
The technical scheme adopted by the invention is as follows: a method for acquiring temperature distribution of a rolled piece in an induction heating system comprises the following steps:
step one, constructing a process arrangement of an induction heating system, and acquiring process parameters;
dividing rolled piece width direction-thickness direction section grid nodes, and establishing a coordinate system parameter and temperature calculation model of each node;
step three, acquiring initial temperature distribution of a rolled piece at an initial calculation position;
and step four, calculating the temperature distribution of the cross section of the rolled piece at the position from the initial calculation position to the final calculation position, 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 corresponding temperature calculation model.
According to the scheme, in the second step, rolled piece width direction-thickness direction section grid nodes are divided, and a specific method for establishing coordinate system parameters of each node is as follows:
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 one-half of the width-thickness cross section of the 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
The half thickness of the rolled piece is equally divided into M sections, j is 1, 2 and 3
Δzi,j＝Δzi,j-1＝Δzi,j+1Δ z, wherein Δ zi,jFor the thickness-wise length, Δ z, of the corresponding element of the rolled-stock node (i, j)i,j-1The thickness-wise length, Δ z, of the corresponding element for the rolled-stock node (i, j-1)i,j+1The thickness direction length of the corresponding unit of the rolled piece node (i, j +1) is obtained.
According to the scheme, in the step two, the specific method for establishing the temperature calculation model of each node comprises the following steps: according to the positions of the nodes, the nodes are divided into internal nodes, surface nodes, end nodes, core nodes and corner nodes;
(1) the internal node may be denoted by (i, j), i is 2, 3, 4.. N-1, j is 2, 3, 4.. M-1, and the abscissa of the node (i, j) is:
the ordinate is:
the expression of the temperature calculation model of the internal node is as follows:
(2) the serial number of the surface node can be represented as (i, M), i is 2, 3, 4
The ordinate is
The expression of the temperature calculation model of the surface node is as follows:
(3) the end node may be denoted by (N, j), j being 2, 3, 4
The ordinate is
The expression of the temperature calculation model of the end node is as follows:
(4) for the core node, the serial number of the node of the thick end face of the core can be represented as (i,1), i is 2, 3 and 4
The ordinate is
The number of the nodes of the end faces of the widthwise core portion can be represented by (1, j), j being 2, 3, 4
The ordinate is
(i) If the core node is a thick core node, the temperature calculation model expression is as follows:
(ii) if the core node is a wide core node, the temperature calculation model expression is as follows:
(5) and for corner nodes, the serial number of the top left corner node is (1, M), and the abscissa is
The ordinate is
The node number at the upper right corner is (N, M), and the abscissa is
The ordinate is
The node sequence number of the lower left corner is (1,1), and the abscissa is
The ordinate is
The node number of the lower right corner is (N,1), and the abscissa is
The ordinate is
(i) If the node is the upper left corner node, the temperature calculation model expression is as follows:
(ii) if the node is the upper right corner node, the temperature calculation model expression is as follows:
(iii) if the node is the lower left corner node, the temperature calculation model expression is as follows:
(iv) if the node is the lower right corner node, the temperature calculation model expression is as follows:
in the formulae (1) to (9), haIs the natural convection cooling heat exchange coefficient of air, and has the unit W/(mm)2×℃)；TaIs ambient temperature, in units; epsilonrThe thermal emissivity is a rolled piece; sigma0Radiation coefficient of absolute black body, σ0＝5.67×10-6W/(mm2×K4) (ii) a c is the specific heat capacity of the rolled piece, and the unit is J/(kg multiplied by DEG C); rho is the material density of the rolled piece in kg/mm3(ii) a Lambda is the thermal conductivity (heat conductivity) of the rolled piece, and the unit W/(mm X DEG C); b is the width of the rolled piece in unit mm; delta t is the calculation time increment and takes the value of 0.0008S;
is the heat generated by a heat source in a unit time and a unit volume of rolled pieces, and has a unit of W/mm3；
Is the temperature of the node (i, j) at the current time in units;
is the temperature of the node (i, j) at the previous time, in units;
is the temperature of the node (i-1, j) at the last moment in units;
is the temperature of the node (i +1, j) at the last time, in units;
is the temperature of the node (i, j-1) at the last moment in units;
is the temperature of node (i, j +1) at the previous time in units of deg.C.
According to the scheme, the specific method of the step three 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 The temperature of each node is increased in turn from the surface to the middle thickness in the thickness direction of the rolled piece, the temperature difference between adjacent nodes is in equal proportion, and the proportionality coefficient is gamma, namely
And the temperature difference between the middle thickness position and the surface is T0(i,1)-T0(i,M)＝ΔT2The initial rolled piece temperature profile may be expressed as:
(1) when gamma > 1, the respective initial temperatures are:
(2) when γ is 1, each initial temperature is:
according to the scheme, in the fourth step, if the section of the rolled piece is positioned in the induction heating area, the serial number of the induction heating device where the section of the rolled piece is positioned is further judged, and then the corresponding temperature calculation model is adopted for calculation under the heating condition of the induction heating device; and if the section of the rolled piece is positioned in the non-induction heating area, further judging the serial number of an air cooling area where the section of the rolled piece is positioned, and then calculating by adopting a corresponding temperature calculation model under the cooling condition of the air cooling area.
According to the scheme, if the cross section of the rolled piece is positioned in the induction heating area, the induction heating heat absorbed by the rolled piece per unit time unit volume in the ith induction heating device is
In the formula (12), ηiIs the heat transfer coefficient, P, of the i-th induction heating unitiHeating power of the i-th induction heating device in units of W and BiThe heating width of the single side of a rolled piece in the ith induction heating device is unit mm; h is the thickness of a rolled piece in mm; l isiIs the length of the ith induction heating device in mm.
According to the scheme, if the cross section of the rolled piece is positioned in the air cooling area, the induction heating heat absorbed by the rolled piece per unit volume in unit time is
According to the scheme, the heating conversion coefficient of the induction heating device is 0.5-0.9.
According to the scheme, the thermal radiance of the rolled piece in the air cooling area is 0.4-0.9.
The invention has the beneficial effects that: the invention establishes a two-dimensional rolled piece temperature field difference equation, calculates the temperature distribution of the rolled piece based on the set process arrangement and process parameters (including the number, the spacing, the heating length and the heating power of induction heating devices, the thickness, the width, the speed and the thermophysical parameters of the rolled piece, the ambient temperature, the air cooling heat exchange coefficient and the thermal emissivity) of an induction heating system, has less assumption and simplification conditions, can accurately forecast the temperature field change of the rolled piece in the induction heating process, has higher calculation precision than an analytic method and higher calculation speed than a finite element method, can meet the requirements of high precision and high efficiency of online calculation, and improves the control precision of the temperature of the rolled piece in the induction heating process and the forecast precision of the microstructure of the rolled piece.
Drawings
Fig. 1 is a schematic diagram of an induction heating system arrangement according to an embodiment of the present invention.
Fig. 2 is a flowchart of calculating the temperature distribution of the rolled piece in this embodiment.
FIG. 3 is a schematic diagram of a half-width-thickness cross-sectional discrete grid of the rolled piece in this embodiment.
FIG. 4 shows the mean temperature variation of the surface, core and thickness of the rolled piece at the mid-width in the induction heating system calculated in this example.
In the drawing, 1, rolled piece; 2. an induction heating device; 3. and (4) an air cooling area.
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.
A method for acquiring the temperature distribution of a rolled piece in an induction heating system is shown in figure 2, and specifically comprises the following steps:
step one, constructing the process arrangement of the induction heating system, and acquiring process parameters including the thickness, width, speed and thermophysical parameters of the rolled piece 1, the number, the interval, the heating length and the heating power of the induction heating devices 2, the ambient temperature, the air cooling heat exchange coefficient and the thermal radiation coefficient.
In the present embodiment, as shown in fig. 1, the induction heating system includes 12 induction heating devices 2 and 11 air cooling zones 3, and the rolled product 1 is partially located in the induction heating zone of the induction heating device 2 (i.e., the coverage area of the induction heating device 2) and partially located in the air cooling zone 3. The thickness of the rolled piece 1 is 17mm, the width of the rolled piece 1 is 1600mm, the speed of the rolled piece 1 is 28m/min, the thermal conductivity of the rolled piece 1 is 30W/(m x DEG C), the specific heat capacity of the rolled piece 1 is 670J/(kg x DEG C), and the density of the rolled piece is 7800kg/m3The ambient temperature was 30 ℃, and the relevant parameters of the induction heating device 2 and the air cooling zone 3 are shown in table 1 and table 2, respectively.
Table 1 relevant parameters of the induction heating device 2
TABLE 2 relevant parameters of air-cooled zone 3
And step two, dividing the grid nodes of the width-thickness section of the rolled piece 1, and establishing a coordinate system parameter and temperature calculation model of each node. The specific method comprises the following steps: as shown in FIG. 3, a y-z rectangular coordinate system is established, wherein the y-axis is located at the position of the intermediate thickness of the rolled product 1 and the z-axis is located at the position of the intermediate width of the rolled product 1. One half of the width-thickness cross section of the rolled piece 1 is discretized into 20 multiplied by 20 grids, wherein the half width of the rolled piece 1 is equally divided into 20 segments, i is 1, 2, 3
The half thickness of the rolled piece 1 is equally divided into 20 sections, j is 1, 2, 3, 20, and the grid thickness is
Δzi,j＝Δzi,j-1＝Δzi,j+1Δ z, wherein Δ zi,jThe thickness-wise length, Δ z, of the element corresponding to the node (i, j) of the rolled stock 1i,j-1The thickness direction length, Delta z, of the unit corresponding to the node (i, j-1) of the rolled piece 1i,j+1The node (i, j +1) of the rolled piece 1 corresponds to the thickness length of the unit.
According to the positions of the nodes, the nodes are divided into internal nodes, surface nodes, end nodes, core nodes and corner nodes;
(1) the internal node may be denoted by (i, j), i is 2, 3, 4.. N-1, j is 2, 3, 4.. M-1, and the abscissa of the node (i, j) is:
the ordinate is:
the expression of the temperature calculation model of the internal node is as follows:
(2) the serial number of the surface node can be represented as (i, M), i is 2, 3, 4
The ordinate is
The expression of the temperature calculation model of the surface node is as follows:
(3) the end node may be denoted by (N, j), j being 2, 3, 4
The ordinate is
The expression of the temperature calculation model of the end node is as follows:
(4) for the core node, the serial number of the node of the thick end face of the core can be represented as (i,1), i is 2, 3 and 4
The ordinate is
The number of the nodes of the end faces of the widthwise core portion can be represented by (1, j), j being 2, 3, 4
The ordinate is
(i) If the core node is a thick core node, the temperature calculation model expression is as follows:
(ii) if the core node is the widthwise core node, the temperature calculation model expression is as follows:
(5) and for corner nodes, the serial number of the top left corner node is (1, M), and the abscissa is
The ordinate is
The node number at the upper right corner is (N, M), and the abscissa is
The ordinate is
The node sequence number of the lower left corner is (1,1), and the abscissa is
The ordinate is
The node number of the lower right corner is (N,1), and the abscissa is
The ordinate is
(i) If the node is the upper left corner node, the temperature calculation model expression is as follows:
(ii) if the node is the upper right corner node, the temperature calculation model expression is as follows:
(iii) if the node is the lower left corner node, the temperature calculation model expression is as follows:
(iv) if the node is the lower right corner node, the temperature calculation model expression is as follows:
in the formulae (1) to (9), haIs the natural convection cooling heat exchange coefficient of air, and has the unit W/(mm)2×℃)；TaIs ambient temperature, in units; epsilonrThe thermal emissivity is a rolled piece; sigma0Radiation coefficient of absolute black body, σ0＝5.67×10-6W/(mm2×K4) (ii) a c is the specific heat capacity of the rolled piece 1, and the unit J/(kg × ° C); rho is the material density of the rolled piece 1 and the unit kg/mm3(ii) a λ is the thermal conductivity (thermal conductivity) of the rolled piece 1, in units W/(mm × ° c); b isIs the width of the rolled piece 1 and is unit mm; delta t is the calculation time increment and takes the value of 0.0008S;
the heat generated by the internal heat source of the unit time unit volume of the rolled piece 1 (namely the induction heating heat absorbed by the unit time unit volume of the rolled piece 1) is unit of W/mm3；
Is the temperature of the node (i, j) at the current time in units;
is the temperature of the node (i, j) at the previous time, in units;
is the temperature of the node (i-1, j) at the last moment in units;
is the temperature of the node (i +1, j) at the last time, in units;
is the temperature of the node (i, j-1) at the last moment in units;
is the temperature of node (i, j +1) at the previous time in units of deg.C.
And step three, acquiring the initial temperature distribution of the rolled piece 1 at the initial calculation position. The specific method comprises the following steps: in the width direction of the rolled piece 1, the temperature of each node is reduced from the middle width to the edge in sequence, the temperature difference between adjacent nodes is the same, and the temperature difference between the middle width and the edge is T0(1,j)-T0(N,j)＝ΔT1(ii) a The temperature of each node is increased in turn from the surface to the middle thickness in the thickness direction of the rolled piece 1, the temperature difference between the adjacent nodes is in equal proportion, and the proportionality coefficient is gamma, namely
And the temperature difference between the middle thickness position and the surface is T0(i,1)-T0(i,M)＝ΔT2The initial temperature profile of the rolling stock 1 can be expressed as:
(1) when gamma > 1, the respective initial temperatures are:
(2) when γ is 1, each initial temperature is:
in the present embodiment, the initial surface temperature T of the rolled material 1 at the intermediate width0(1,20)＝830℃，T0(1,1)＝850℃，T0The initial temperature at the other nodes of the rolled stock 1 is therefore (20,20) ═ 800 ℃, and the proportionality coefficient γ ═ 1.053:
and step four, calculating the temperature distribution of the section of the rolled piece 1 at the position from the initial calculation position, judging the heat exchange area where the section of the rolled piece 1 is located when calculating once, and then calculating the temperature distribution of the section of the rolled piece 1 at the position by adopting a corresponding temperature model. The specific method comprises the following steps: if the section of the rolled piece 1 is positioned in the induction heating area, further judging the serial number of an induction heating device 2 in which the section of the rolled piece 1 is positioned, and then calculating by adopting a corresponding temperature calculation model under the heating condition of the induction heating device 2; and if the section of the rolled piece 1 is positioned in the air cooling area 3, further judging the serial number of the air cooling area 3 where the section of the rolled piece 1 is positioned, and then adopting a corresponding temperature calculation model to calculate under the cooling condition of the air cooling area 3.
1. If the rolled stock 1 is located in the induction heating zone in cross section, the induction heating heat absorbed per unit time per unit volume of the rolled stock 1 in the i-th induction heating device 2 is
In the formula (12), ηiThe heating conversion coefficient of the ith induction heating device 2 is 0.5-0.9; piIs the heating power of the i-th induction heating apparatus 2 in units of W; b isiThe heating width of the rolled piece 1 in the ith induction heating device 2 on one side is unit mm; h is the thickness of the rolled piece 1 in mm; l isiIs the length of the i-th induction heating unit 2 in mm.
2. If the cross section of the rolled piece 1 is located in the air cooling zone 3, the induction heating heat absorbed by the rolled piece 1 per unit volume per unit time is
And the thermal radiance of the rolled piece in the air cooling area 3 is 0.4-0.9.
Fig. 4 shows the average temperature variation of the surface, the core and the thickness direction at the middle width of the rolled piece 1 in the induction heating system calculated by the embodiment, and also shows the measured surface temperature values (shown by black circles in fig. 4) at the middle width of the rolled piece 1 at two typical positions (respectively in the 3 rd air cooling zone 3, the 6 th air cooling zone 3 and the 9 th air cooling zone 3).
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 (9)
1. A method for acquiring temperature distribution of a rolled piece in an induction heating system is characterized by comprising the following steps:
step one, constructing a process arrangement of an induction heating system, and acquiring process parameters;
dividing rolled piece width direction-thickness direction section grid nodes, and establishing a coordinate system parameter and temperature calculation model of each node;
step three, acquiring initial temperature distribution of a rolled piece at an initial calculation position;
and step four, calculating the temperature distribution of the cross section of the rolled piece at the position from the initial calculation position to the final calculation position, 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 corresponding temperature calculation model.
2. The method for acquiring the temperature distribution of the rolled piece in the induction heating system according to claim 1, wherein in the second step, the width-thickness section grid nodes of the rolled piece are divided, and the specific method for establishing the coordinate system parameters of each node 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 one-half of the width-thickness cross section of the 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
The half thickness of the rolled piece is equally divided into M sections, j is 1, 2 and 3
Δzi,j＝Δzi,j-1＝Δzi,j+1Δ z, wherein Δ zi,jFor the thickness-wise length, Δ z, of the corresponding element of the rolled-stock node (i, j)i,j-1The thickness-wise length, Δ z, of the corresponding element for the rolled-stock node (i, j-1)i,j+1The thickness direction length of the corresponding unit of the rolled piece node (i, j +1) is obtained.
3. The method for acquiring the temperature distribution of the rolled piece in the induction heating system according to claim 2, wherein in the second step, the specific method for establishing the temperature calculation model of each node is as follows: according to the positions of the nodes, the nodes are divided into internal nodes, surface nodes, end nodes, core nodes and corner nodes;
(1) the internal node may be denoted by (i, j), i is 2, 3, 4.. N-1, j is 2, 3, 4.. M-1, and the abscissa of the node (i, j) is:
the ordinate is:
the expression of the temperature calculation model of the internal node is as follows:
(2) the serial number of the surface node can be represented as (i, M), i is 2, 3, 4
The ordinate is
The expression of the temperature calculation model of the surface node is as follows:
(3) the end node may be denoted by (N, j), j being 2, 3, 4
The ordinate is
The expression of the temperature calculation model of the end node is as follows:
(4) for the core node, the serial number of the node of the thick end face of the core can be represented as (i,1), i is 2, 3 and 4
The ordinate is
The number of the nodes of the end faces of the widthwise core portion can be represented by (1, j), j being 2, 3, 4
The ordinate is
(i) If the core node is a thick core node, the temperature calculation model expression is as follows:
(ii) if the core node is a wide core node, the temperature calculation model expression is as follows:
(5) and for corner nodes, the serial number of the top left corner node is (1, M), and the abscissa is
The ordinate is
The node number at the upper right corner is (N, M), and the abscissa is
The ordinate is
The node sequence number of the lower left corner is (1,1), and the abscissa is
The ordinate is
The node number of the lower right corner is (N,1), and the abscissa is
The ordinate is
(i) If the node is the upper left corner node, the temperature calculation model expression is as follows:
(ii) if the node is the upper right corner node, the temperature calculation model expression is as follows:
(iii) if the node is the lower left corner node, the temperature calculation model expression is as follows:
(iv) if the node is the lower right corner node, the temperature calculation model expression is as follows:
in the formulae (1) to (9), haIs the natural convection cooling heat exchange coefficient of air, and has the unit W/(mm)2×℃)；TaIs ambient temperature, in units; epsilonrThe thermal emissivity is a rolled piece; sigma0Radiation coefficient of absolute black body, σ0＝5.67×10-6W/(mm2×K4) (ii) a c is the specific heat capacity of the rolled piece, and the unit is J/(kg multiplied by DEG C); rho is the material density of the rolled piece in kg/mm3(ii) a Lambda is the thermal conductivity (heat conductivity) of the rolled piece, and the unit W/(mm X DEG C); b is the width of the rolled piece in unit mm; delta t is the calculation time increment and takes the value of 0.0008S;
is the heat generated by a heat source in a unit time and a unit volume of rolled pieces, and has a unit of W/mm3；
Is the temperature of the node (i, j) at the current time in units;
is the temperature of the node (i, j) at the previous time, in units;
is the temperature of the node (i-1, j) at the last moment in units;
is the temperature of the node (i +1, j) at the last time, in units;
is the temperature of the node (i, j-1) at the last moment in units;
is the temperature of node (i, j +1) at the previous time in units of deg.C.
4. The method for acquiring the temperature distribution of the rolled piece in the induction heating system as set forth in claim 2, wherein the specific method in the third step is: 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 The temperature of each node is increased in turn from the surface to the middle thickness in the thickness direction of the rolled piece, the temperature difference between adjacent nodes is in equal proportion, and the proportionality coefficient is gamma, namely
And the temperature difference between the middle thickness position and the surface is T0(i,1)-T0(i,M)＝ΔT2The initial rolled piece temperature profile may be expressed as:
(1) when gamma > 1, the respective initial temperatures are:
(2) when γ is 1, each initial temperature is:
5. the method for acquiring the temperature distribution of the rolled piece in the induction heating system according to claim 3, wherein in the fourth step, if the cross section of the rolled piece is located in the induction heating area, the serial number of the induction heating device where the cross section of the rolled piece is located is further judged, and then the calculation is performed by adopting a corresponding temperature calculation model under the heating condition of the induction heating device; and if the section of the rolled piece is positioned in the non-induction heating area, further judging the serial number of an air cooling area where the section of the rolled piece is positioned, and then calculating by adopting a corresponding temperature calculation model under the cooling condition of the air cooling area.
6. A method of obtaining a temperature profile of a rolled product in an induction heating system as defined in claim 5, wherein if the cross-section of the rolled product is in the induction heating zone, the induction heating heat absorbed per unit time per unit volume of the rolled product in the ith induction heating unit is
In the formula (12), ηiIs the heat transfer coefficient, P, of the i-th induction heating unitiHeating power of the i-th induction heating device in units of W and BiThe heating width of the single side of a rolled piece in the ith induction heating device is unit mm; h is the thickness of a rolled piece in mm; l isiIs the length of the ith induction heating device in mm.
7. A method of deriving a temperature profile of a rolled member in an induction heating system as claimed in claim 5, wherein if the cross-section of the rolled member is in the air cooling zone, the amount of induction heating heat absorbed per unit volume of the rolled member per unit time is
8. The method for acquiring the temperature distribution of the rolled piece in the induction heating system according to claim 6, wherein the heating conversion coefficient of the induction heating device is 0.5 to 0.9.
9. The method for obtaining the temperature distribution of the rolled piece in the induction heating system according to claim 6, wherein the thermal emissivity of the rolled piece in the air cooling zone is 0.4-0.9.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111256027.3A CN113899471A (en) | 2021-10-27 | 2021-10-27 | Method for acquiring temperature distribution of rolled piece in induction heating system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111256027.3A CN113899471A (en) | 2021-10-27 | 2021-10-27 | Method for acquiring temperature distribution of rolled piece in induction heating system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113899471A true CN113899471A (en) | 2022-01-07 |
Family
ID=79026619
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111256027.3A Pending CN113899471A (en) | 2021-10-27 | 2021-10-27 | Method for acquiring temperature distribution of rolled piece in induction heating system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113899471A (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1252521A (en) * | 1998-10-22 | 2000-05-10 | 冶金工业部钢铁研究总院 | Temperature measuring method for steel plate and band rolling process |
| CN102564644A (en) * | 2011-12-20 | 2012-07-11 | 东北大学 | Temperature online measuring method for plate blank in production process of heating furnace |
| CN103547385A (en) * | 2011-05-24 | 2014-01-29 | 西门子公司 | Control method for a mill train |
| 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 |
| CN110765671A (en) * | 2019-12-05 | 2020-02-07 | 中冶南方工程技术有限公司 | Method and device for acquiring temperature distribution of hot rolling working roll |
| CN111079275A (en) * | 2019-12-05 | 2020-04-28 | 中冶南方工程技术有限公司 | Rolled piece temperature obtaining method and device for strip hot rolling production line |
-
2021
- 2021-10-27 CN CN202111256027.3A patent/CN113899471A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1252521A (en) * | 1998-10-22 | 2000-05-10 | 冶金工业部钢铁研究总院 | Temperature measuring method for steel plate and band rolling process |
| CN103547385A (en) * | 2011-05-24 | 2014-01-29 | 西门子公司 | Control method for a mill train |
| CN102564644A (en) * | 2011-12-20 | 2012-07-11 | 东北大学 | Temperature online measuring method for plate blank in production process of heating furnace |
| 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 |
| CN110765671A (en) * | 2019-12-05 | 2020-02-07 | 中冶南方工程技术有限公司 | Method and device for acquiring temperature distribution of hot rolling working roll |
| CN111079275A (en) * | 2019-12-05 | 2020-04-28 | 中冶南方工程技术有限公司 | Rolled piece temperature obtaining method and device for strip hot rolling production line |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106636610B (en) | A kind of double dimension walking beam furnace heating curve optimal setting methods based on time and furnace superintendent | |
| CN111079275B (en) | Method and device for acquiring temperature of rolled piece in strip hot rolling production line | |
| CN103028615B (en) | Method for predicting temperature evolution in hot continuous rolling process of strip steel | |
| CN103761370B (en) | A kind of Forecasting Methodology of process of plate belt hot rolling surface film thermal conductance | |
| CN106906350B (en) | A kind of H profile steel base heating process temperature distribution calculation method | |
| CN104073623A (en) | Roller hearth type annealing furnace temperature control method | |
| Hajaliakbari et al. | Analysis of thermal energy performance in continuous annealing furnace | |
| CN114818505B (en) | Method for predicting temperature distribution of steel billet in heating furnace based on particle swarm optimization algorithm | |
| CN114410956B (en) | Online control system and method for intermittent aluminum coil annealing furnace | |
| CN103559334A (en) | Modeling method and system for temperature field in laminar cooling | |
| JP4598586B2 (en) | Cooling control method, apparatus, and computer program | |
| CN113899471A (en) | Method for acquiring temperature distribution of rolled piece in induction heating system | |
| CN115016578B (en) | Strip steel quality regulation and control method based on edge temperature control | |
| CN113930600B (en) | Bell-type furnace annealing process monitoring and control method based on digital twin technology | |
| CN115232955B (en) | Optimization control method for strip steel temperature in dynamic heating process of continuous annealing furnace | |
| JP4408221B2 (en) | Heat transfer coefficient estimation method and cooling control method in water cooling process of steel sheet | |
| CN102994732A (en) | Fuel quantity deciding system for heating furnace and deciding method thereof | |
| CN104673992B (en) | The control method of a kind of production line of bar Controlled cooling process and device | |
| CN102799723A (en) | Muffle furnace heating capacity calculation, prediction and simulation method | |
| JP6981830B2 (en) | Steel piece temperature estimation device and method | |
| CN114178325B (en) | Cooling water flow obtaining method and temperature calculating method for hot-rolled carbon steel laminar cooling jet header | |
| CN117521299B (en) | Prediction method of strip steel three-dimensional temperature field in inter-frame cooling process | |
| CN115828624B (en) | Slab heating temperature accurate forecasting method based on black box experiment tuning | |
| CN114178323B (en) | Cooling water flow obtaining method of hot-rolled dual-phase steel laminar flow cooling spray header | |
| Zheng et al. | Optimization target resetting distributed model predictive control for accelerated cooling 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 |