CN102489524B - Machine frame load distribution method for decreasing energy consumption of rolling process of hot rolled strip steel - Google Patents

Machine frame load distribution method for decreasing energy consumption of rolling process of hot rolled strip steel Download PDF

Info

Publication number
CN102489524B
CN102489524B CN 201110390910 CN201110390910A CN102489524B CN 102489524 B CN102489524 B CN 102489524B CN 201110390910 CN201110390910 CN 201110390910 CN 201110390910 A CN201110390910 A CN 201110390910A CN 102489524 B CN102489524 B CN 102489524B
Authority
CN
China
Prior art keywords
frame
rolling
formula
energy consumption
roll
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.)
Active
Application number
CN 201110390910
Other languages
Chinese (zh)
Other versions
CN102489524A (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.)
Northeastern University China
Original Assignee
Northeastern University China
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 Northeastern University China filed Critical Northeastern University China
Priority to CN 201110390910 priority Critical patent/CN102489524B/en
Publication of CN102489524A publication Critical patent/CN102489524A/en
Application granted granted Critical
Publication of CN102489524B publication Critical patent/CN102489524B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

A machine frame load distribution method for decreasing energy consumption of a rolling process of hot rolled strip steel includes the following steps: step 1 determining constraint conditions of an initial control scheme, step 2 determining a control target, step 3 determining parameters of a machine and parameters of rolled pieces of the control scheme, step 4 obtaining thickness of an outlet of each machine frame through improved differential evolution algorithm, step 5 determining strip threading speed, temperature, rolling force, rolling power and total energy consumption of each machineframe based on the thickness of the outlet of each machine frame obtained in step 4, and step 6 judging whether rolling force, roll torque and rolling power exceed rated values of the machine or not.If the values exceed the rated values, repeat step 4, and if not, judge whether total energy consumption reaches the minimum value or not. If the total energy consumption reaches the minimum value, output the final value, and if not, repeat step 4 until the minimum value is reached. Through the improved differential evolution algorithm, load of each machine frame is optimally set, rolling reduction of each machine frame is optimally distributed, and actual rolled thickness of each machine frame is determined so that total energy consumption reaches the minimum value, device damage is reduced,and production efficiency and device utilization ratio are improved.

Description

A kind of frame load distribution method that reduces the strip hot rolling process energy consumption
Technical field
The invention belongs to steel plate rolling process control technology field, be specifically related to a kind of frame load distribution method that reduces the strip hot rolling process energy consumption.
Background technology
When hot rolling strip steel, steel billet comes out in heating furnace, through rolling (n is frame number or road number of times) of n frame or passage, roll out the finished strip that meets quality standard, these a series of operations of rolling are to carry out according to the rolling procedure that sets in advance, as shown in Figure 1.The central issue of formulating rolling procedure is the exit thickness that how to distribute the drafts of each frame, determines each frame, just determine the sharing of load of depressing of each frame, in essence, sharing of load has determined the state characteristic of the operation of rolling, and end product quality requirement, equipment adjustment are all had significant effects.Formulate sharing of load, the restrictive condition that needs to consider has the overload of equipment intensity, motor and the restriction of heating, process conditions, plate shape and speed etc., these factors are conflicting often, under different rolling conditions, the limiting factor of each rolling mill also is different, interact again between each factor, concern very complicated.Therefore, can reach production increases with quality up and cost of production down for making the production of band steel, just need to optimize the drafts that distributes each frame.
For a long time, Chinese scholars has been carried out number of research projects to the precision aspect that improves mm finishing mill unit sharing of load model.But still be based on the experiential operating method in essence, rule of thumb distribute earlier each frame load, calculate then reach the finish to gauge target temperature required wear tape speed, determine that rolled piece is in the temperature of each frame, with mathematical model prediction resistance of deformation, roll-force, rolling power, and other parameters, carry out the limit and check; If the rolling power of a certain frame surpasses the main motor current limit, just utilize correction algorithm to redistribute the drafts of each frame, to adjust each frame load, constantly adjust, check, do not transfinite up to the rolling power of institute's organic frame.Check by after calculate roll gap, speed dispatch control system setting value again, finish the setting of parameter.
This method is easy, reliable, can guarantee device security, and little to the automation degree of equipment dependence, but need constantly check and repair the worker, easily causes each frame load uneven, causes load to side arm or to the phenomenon of after-frame accumulation.The computation model of this method generally can not online adaptive, and can not guarantee in each passage in process of production or each frame can both be rolled under the load that allows, let alone energy-conservation and consider belt plate shape.And when production description or rolling condition change, accumulate experience again again, could make sharing of load again.So not only influence the performance of whole unit capacity, and the quality of influence band steel.This method does not consider to carry out the production process operation optimization from the optimal design of production process too much.
In addition, the scholar who has proposes the thickness allocation optimized method based on related algorithm, as waiting methods such as sharing of load, the combination of plurality of target function, these methods are from the angle of balanced load, minimizing energy consumption, because the load (roll-force, power etc.) in this method and drafts are complicated non-linear relation, and roll-force, parameters such as power are to depress the result of generation, are difficult to obtain relatively good solution, can not well satisfy the needs that hot rolling is produced.
In steel rolling was produced, blank was rolling through number frame (passage), produced plastic deformation, finally rolled out the product of requirement up to specification.The operation of rolling that these are a series of is to enter the rules of setting in advance before the milling train according to blank to carry out, and this is the important process of band steel continuous rolling.Wherein sharing of load is to set prerequisite and the basis of calculating, and is to set the key link of calculating.In steel rolling was produced, energy consumption was one of key factor that influences production cost always.Therefore the synergy that lowers consumption also is the target that the producer pursues.
Summary of the invention
At the deficiencies in the prior art, the invention provides a kind of frame load distribution method that reduces the strip hot rolling process energy consumption, optimize the drafts that distributes each frame, determine the actual thickness that shuts out of each frame, make total energy consumption reach minimum.
Technical scheme of the present invention: a kind of frame load distribution method that reduces the strip hot rolling process energy consumption specifically comprises the steps:
Step 1: determine the constraints of initial control scheme, comprising:
(1) bite condition
Bite condition guarantees nipping smoothly of band steel, nips smoothly to guarantee that the band steel can enter mill milling, if the thickness of rolled piece is big, temperature is high, drafts is big, just can not guarantee nipping smoothly with steel.In order to guarantee to guarantee bite condition with the nipping smoothly of steel with drafts, namely rolling each reduction in pass should allow drafts less than maximum, and computing formula as the formula (1)
&Delta; h i < &Delta; h max = D ( 1 - 1 1 + f 2 ) i=1,2,...,n (1)
Δ h wherein iThe drafts of representing i frame; Δ h MaxFor maximum allows drafts; F is the coefficient of friction between roll and rolled piece; D represents the diameter of roll.
(2) roll strength condition
When maximum pressure that draught pressure can bear greater than roll, the roll part may be destroyed.For guaranteeing roll strength, reply draught pressure and roll torque are limited, and guarantee can not surpass maximum roll-force and roll torque in the operation of rolling, and computing formula is suc as formula shown in (2) formula (3):
P i<P imax i=1,2,...,n (2)
M i<M imax i=1,2,...,n (3)
P wherein iIt is the roll-force of i frame; P ImaxIt is the maximum rolling force that i frame allows; M iIt is the roll torque of i frame; M ImaxIt is the maximum rolling force square that i frame allows.
(3) motor ability condition
The motor ability is the restrictive condition of motor overload, and rolling power should satisfy the loading condiction of milling train master motor, requires rolling power must not surpass its rated power, and computing formula as the formula (4)
N i<N imax i=1,2,...,n (4)
N wherein iIt is the rolling power of i frame; N ImaxIt is the maximum rolling power that i frame allows.
(4) strip shape quality restrictive condition
In order to keep strip shape quality to reach customer requirement, optimize the drafts of arranging several of backs, making between its corresponding roll-force has certain ratio, prevents from being with steel the limit wave to occur, and computing formula is as the formula (5)
- 40 ( h i b ) 1.86 < &Delta; CR i < 80 ( h i b ) 1.86 i=n-3,n-2,n-1,n
(5)
Δ CR wherein iEntrance convexity when passing through i frame for slab is poor with the outlet convexity; h iIt is the exit thickness of i frame; B is the width of slab.
Step 2: the target of determining control.The target of this method control is to make rolling total energy consumption reach minimum, and the energy consumption calculation formula is shown in (6) formula:
min &Sigma; i = 1 n N i i=1,2,...,n (6)
Wherein: n-mm finishing mill unit frame number;
N i-Di i frame power of motor, computing formula be as the formula (7):
N i = 2 &pi; * 10 3 60 * 102 c i * M i i=1,2,...,n (7)
C in the formula (7) iBe revolution;
M iBe the moment of the horizontal roller of i frame, computing formula as the formula (8):
i=1,2,...,n (8)
Wherein, It is arm of force coefficient;
P iBe the roll-force computing formula of horizontal roller as the formula (9):
P i = B l c / Q p K i=1,2,...,n (9)
B is the width of blank in the formula (9);
Be contact arc length, computing formula as the formula (10):
l c / = R / &Delta; h i i=1,2,...,n (10)
In the v formula (5), R ' is the roller radius after flattening, computing formula as the formula (11):
R i &prime; = R i &CenterDot; ( 1 + 2.2 &times; 10 - 5 P i B&Delta; h i ) i=1,2,...,n (11)
R iIt is the radius of the working roll of i frame;
Δ h iRepresent the drafts computing formula as the formula (12):
Δh i=Δh i-1-h i i=1,2,...,n (12)
H in the formula I-1It is the inlet thickness of i frame; h iIt is the exit thickness of i frame;
Q PBe external friction stress state coefficient, computing formula as the formula (13):
Q p=0.8206+0.2376l c/h c+0.1006εl c/h c-0.3768ε i=1,2,...,n (13)
In the formula ε be the relative reduction computing formula as the formula (14)
&epsiv; = h i - 1 - h i h i - 1 i=1,2,...,n (14)
l cBe the floor projection length of contact arc, computing formula as the formula (15):
l c = R i &Delta; h i i=1,2,...,n (15)
h cThe mean value of the inlet thickness of rolled piece and exit thickness when being rolling, computing formula as the formula (16)
h c = h i - 1 + h i 2 i=1,2,...,n (16)
Deformation drag under the K plane deformation, computing formula are formula (17) formula (18)
K=1.15σ (17)
&sigma; = &sigma; 0 exp ( a 1 T + a 2 ) ( u m 10 ) ( a 3 T + a 4 ) [ a 6 ( e 0.4 ) a 5 - ( a 6 - 1 ) ( e 0.4 ) ] - - - ( 18 )
The temperature of steel plate when wherein, T is rolling;
E is the actual texturing process degree, and computing formula as the formula (19)
e = ln h i - 1 h i i=1,2,...,n (19)
u mBe average deformation speed, computing formula as the formula (20)
u m = v i l c e i=1,2,...,n (20)
v iBe the mill speed of i frame, computing formula as the formula (21)
v i = v n h n h i i=1,2,...,n (21)
v nRoll linear velocity for last milling train;
h nExit thickness for last milling train;
σ 0, a 1~a 6Be regression coefficient.
T iBe the rolling temperature of i milling train, computing formula as the formula (22)
T i - T W T F 0 - T W = exp ( - K a &Sigma; i = 1 n L i h n v n ) i=1,2,...,n (22)
L 1-finish rolling entrance point for measuring temperature is to F 1Distance;
L i-i-1 frame is to the distance (i=2~7) of i frame; m
L n-finish rolling end frame is to the distance of finish rolling outlet temperature measurer loca; m
T WThe water temperature of spraying water between-frame;
K a-comprehensive convection current cooling ratio;
T F0-mm finishing mill unit porch estimated temperature, computing formula as the formula (23)
T F 0 = 100 ( 6 &epsiv;&sigma; 100 &gamma; c p h 0 &tau; + ( T Rc 100 ) - 3 ) - 1 3 - - - ( 23 )
σ is this graceful constant of Si Difen-bohr;
ε is blackness;
C pBe specific heat capacity;
γ is density;
τ is for exporting to the band steel run duration of finishing mill entrance from roughing mill;
h 0Be roughing unit exit actual measurement thickness;
T RCBe roughing unit exit observed temperature.
Step 3: determine the parameter of machine parameter and the rolled piece of control scheme, machine parameter comprises running time, the rated power of the work roll diameter of mm finishing mill unit, backing roll diameter, motor, the rated speed of motor, maximum rolling force and maximum rolling force square; The parameter of rolled piece comprises the width B of rolled piece, the supplied materials thickness H of rolled piece, the finished product thickness h of rolled piece, the outlet temperature TFC of roughing mill, finish rolling outlet temperature T.
Step 4: utilize improved differential evolution algorithm to obtain each frame exit thickness, concrete steps are as follows:
Step 4.1: the structure initial population, and population scale and iterations are set.Initial population of the present invention adopts real coding, and the step that initial population produces is as follows:
If have N n dimension individual in the population, each individual vectorial X iRepresentation be X i=[x I1, x I2..., x In].
Step 4.1.1: make x 01=500, x 0n=-500,
Step 4.1.2: utilize formula (24) to calculate x 0j, obtain X 0, X 0=[x 01, x 02..., x 0n].
x 0j=x 0j-Slot,j=2,3,…,n-1 (24)
Step 4.1.3: make i=1.
Step 4.1.4: utilize formula (25), with X 0Be template, produce the individual X in the population i, i=1,2 ..., N.x i1=500,x in=-500,
x ij=x 0j+f·Slot,j=2,3,…,n-1 (25)
In the formula, f is random number, f ∈ (1,1).
Step 4.1.5:i=i+1.
Step 4.1.6: if i=N, then the initial population structure finishes; Otherwise, change step 4.1.4 over to.
Step 4.2: the individual vector of improved differential evolution algorithm is mapped to the solution vector space from real number space, and concrete steps are as follows:
Step 4.2.1: find out maximum and the minimum of a value of each component among the individual vectorial X, be designated as x MaxAnd x Min, make v Max=h 0, v Min=h n, v wherein MaxAnd v MinBe respectively x MaxWith x MinThe corresponding component in decoding back; h 0Be mm finishing mill unit strip steel at entry thickness; h nBe mm finishing mill unit outlet belt steel thickness.For other component among the X, be mapped to solution space with formula (26):
v i = v max - ( x max - x i ) ( x max - x min ) * ( v max - v min ) , i=1,2,…,n,
i≠max,min (26)
Step 4.2.2: with the component (v of decoded intermediate vector V 1, v 2..., v n) by arranging from big to small.
Step 4.2.3: with the component v after arranging iUpgrade solution vector U:u 1=v Max..., u n=v Min
Step 4.3: for the solution vector after the decoding, if do not satisfy constraints, carry out and repair strategy, concrete steps are as follows:
Step 4.3.1: if decoded solution vector does not satisfy constraints, then enter step 4.3.2; Otherwise, preserve current solution vector, finish to repair strategy.
Step 4.3.2: the drafts Δ h that calculates each frame of solution vector correspondence i, and these drafts are arranged from big to small: Δ h 1=Δ h Max..., Δ h n=Δ h Min
Step 4.3.3: the drafts of redistributing for each frame, and recomputate the exit thickness of each frame, upgrade solution vector.x 1=h 0+Δh 1,x 2=x 1+Δh 2,...,x n=x n-1+Δh n
Step 4.3.4: if the solution vector after upgrading can satisfy constraints, preserve current solution vector, finish to repair strategy; Otherwise, enter step 4.3.5.
Step 4.3.5: make i=n.
Step 4.3.6: the component x that checks solution vector iWhether satisfy constraints.If x iDo not satisfy constraint, enter step 4.3.7, otherwise, step 4.3.9 entered.
Step 4.3.7:x i=x i-0.1, i=1,2 ..., n.
Step 4.3.8: if x i≤ x I+1(i=1,2 ..., n), change step 4.3.11 over to; If x iSatisfy constraint, then keep x i, enter step 4.3.9; Otherwise, continue step 4.3.7.
Step 4.3.9: if i=1 preserves current solution vector, finish to repair strategy; Otherwise, enter step 4.3.10.
Step 4.3.10:i=i-1 enters step 4.3.6.
Step 4.3.11: repairing failure, add a big penalty value to object function, abandon current solution.
Step 4.4: calculate each the individual fitness value in the population.The present invention will control target as the fitness function that calculates ideal adaptation degree value.
Step 4.5: check the algorithm end condition.See whether fitness value reaches minimum, and satisfy the constraints that the present invention considers.If satisfy the algorithm end condition, then stop and export the individuality of fitness value optimum in the population; Otherwise, continue next step.
Step 4.6: the population to algorithm carries out mutation operation.For individual X I, G: i=1,2 ..., N, a new individual V I, G+1But through type (27) produces:
V i , G + 1 = X r 1 , G + F * ( X r 2 , G - X r 3 , G ) - - - ( 27 )
Wherein, r 1, r 2And r 3Be to go up the mutually different integer of picked at random from interval [1, N], and be different from subscript index i, mutagenic factor F value is between [0,1].
Step 4.7: the population to algorithm carries out interlace operation.Interlace operation as the formula (28).
v i , j G + 1 = u i , j G + 1 , if ( ( randb ( j ) &le; CR ) or ( j = mbr ( i ) ) ) x i , j G , if ( ( randb ( j ) > CR ) and ( j &NotEqual; mbr ( i ) ) ) , - - - ( 28 )
i=1,2,…,N,j=1,2,…,Dv
In the formula, Be J dimension component; Be J dimension component, Individuality for the variation generation; Be J dimension component, Be the parent individuality; Randb (j) is equally distributed probability between [0,1]; An integer that generates at random between mbr (i) expression [1, Dv], CR is crossover probability, generally gets the number between [0,2].Here
CR = CR min + g C ( 1 - CR min &CenterDot; F ( X i ) F ( U i ) + F ( X i ) ) - - - ( 29 )
In the formula, F (X i) expression i parent individuality fitness value; F (U i) fitness value of variation individuality of i parent correspondence of expression; CR MinBe minimum crossover probability; G is current iteration algebraically; G is algorithm greatest iteration algebraically.
Step 4.8: operation is selected in the filial generation in the population.The selection operation of differential evolution algorithm is between parent colony and progeny population, the individual and man-to-man competition of offspring individual of parent.The selection operation of algorithm as the formula (30).
X i G + 1 = V i G + 1 , if ( F ( V i G + 1 ) &le; F ( X i G ) ) X i G , otherwise , i=1,2,…,N (30)
Step 4.9: return step 4.2.
Step 5: what the exit thickness of each frame that obtains according to step 4 was determined each frame wears tape speed, temperature, roll-force, rolling power, total energy consumption.
Step 6: judge whether roll-force, roll torque, rolling power exceed the rated value of machine, exceed, and then repeating step 4, otherwise judge whether total energy consumption reaches minimum, is then to export end value, otherwise repeating step 4, up to reaching minimum of a value.
Beneficial effect: the present invention is optimized setting by improved differential evolution algorithm to the load of each frame of hot rolling production process, reach the purpose that cuts down the consumption of energy, according to frame load distribution method of the present invention, the exit thickness of each frame, its result of rolling power are better than experience sharing of load result, simultaneously because rolling power optimization distribution, reduce equipment damage, improved production efficiency and utilization rate of equipment and installations.
Description of drawings
Schematic diagram is produced in the hot rolling of Fig. 1 iron and steel;
Fig. 2 embodiment of the invention frame load distribution method flow chart;
The rolling power of two kinds of each frames of method of Fig. 3 embodiment of the invention seven frames distributes contrast;
The rolling power of two kinds of each frames of method of Fig. 4 embodiment of the invention six frames distributes contrast.
The specific embodiment
The present invention will be further described below in conjunction with embodiment and accompanying drawing.
Embodiment 1
Be example with seven frame hot fine rolling unit production lines.Steel plate material is Q235, strip width B 0Be 1535mm, belt steel thickness H 0Be 36.7mm, finished product thickness h nBe 5.7mm, roughing outlet temperature t RCBe 1067 ℃, finish rolling outlet temperature t FCBe 891 ℃, be 21s running time, the distance L=5.5m between frame, L 1=14.5m, L 8=8.5m, the coefficient of friction f=0.45 between roll and rolled piece, density γ=7800kg/m 3, this graceful constant σ=5.6662J/m of Si Difen-bohr 2SK 4, last frame is worn tape speed v 7=10m/s, machinery and the technological parameter of other capital equipments see Table 1.
Table 1 mm finishing mill unit major parameter
Carry out the inventive method, specific as follows:
Step 1: determine the constraints of control, constraints is shown below:
&Delta; h i < &Delta; h max = D ( 1 - 1 1 + f 2 ) i=1,...,7
P i<P imax i=1,...,7
M i<M imax i=1,...,7
N i<N imax i=1,...,7
- 40 ( h b ) 1.86 < &Delta; CR i < 80 ( h b ) 1.86 i=4,…,N
Step 2: the target of determining control.The target of the present invention control is to make total energy consumption reach minimum, the energy consumption calculation formula as shown in the formula:
min &Sigma; i = 1 n N i i=1,...,7,n=7
Step 3: the parameter of determining machine parameter and the rolled piece of control scheme.Width B=the 1535mm that comprises rolled piece, the supplied materials thickness H=36.7mm of rolled piece, the finished product thickness h=5.7mm of rolled piece, the outlet temperature T of roughing mill FC=1067 ℃, T=891 ℃ of finish rolling outlet temperature, running time τ=21s, the work roll diameter of mm finishing mill unit, the backing roll diameter, the rated power of motor, the rated speed of motor, maximum rolling force, the maximum rolling force square is as shown in table 1.
Step 4: utilize improved differential evolution algorithm to obtain the exit thickness of each frame, concrete steps are as follows:
Step 4.1: the structure initial population, and population scale and maximum iteration time are set.Population of the present invention is 60, and maximum iteration time is 1000, and initial population adopts real coding, and the step that initial population produces is as follows:
If have 60 7 dimensions individual in the population, each individual vectorial X iRepresentation be X i=[x I1, x I2..., x I7].
Step 4.1.1: make x 1=500, x 7=-500,
Step 4.1.2: utilize formula (24) to calculate x j, obtain X 0: X 0=[x 01, x 02..., x 07].
x j=x j-1-Slot,j=2,3,…,6 (24)
Step 4.1.3: make i=1.
Step 4.1.4: utilize formula (25), with X 0Be template, produce the individual X in the population i, i=1,2 ..., N.x i1=500,x i7-500,x ij=x 0j+f·Slot,j=2,3,…,6 (25)
In the formula, f is random number, f ∈ (1,1).
Step 4.1.5:i=i+1.
Step 4.1.6: if i=N, then the initial population structure finishes; Otherwise, change step 4.1.4 over to.
Step 4.2: the individual vector of adaptive differential descent algorithm is mapped to the solution vector space from real number space.Concrete steps are as follows:
Step 4.2.1: find out maximum and the minimum of a value of each component among the individual vectorial X, be designated as x MaxAnd x Min, make v Max=h 0, v Min=h 7(v MaxAnd v MinBe respectively x MaxWith x MinThe corresponding component in decoding back; h 0Be mm finishing mill unit strip steel at entry thickness; h 7Be mm finishing mill unit outlet belt steel thickness).For other component among the X, be mapped to solution space with formula (26):
v i = v max - ( x max - x i ) ( x max - x min ) * ( v max - v min ) , i=1,2,…,7,i≠max,min (26)
Step 4.2.2: with the component (v of decoded intermediate vector V 1, v 2..., v 7) by arranging from big to small.
Step 4.2.3: with the component v after arranging iUpgrade solution vector U:u 1=v Max..., u 7=v Min
Step 4.3: for the solution vector after the decoding, if do not satisfy constraints, the reparation strategy below adopting, concrete steps are as follows:
Step 4.3.1: if decoded solution vector does not satisfy constraints, then enter step 4.3.2; Otherwise, preserve current solution vector, finish to repair strategy.
Step 4.3.2: the drafts Δ h that calculates each frame of solution vector correspondence i, and these drafts are arranged from big to small: Δ h 1=Δ h Max..., Δ h 7=Δ h Min
Step 4.3.3: the drafts of redistributing for each frame, and recomputate the exit thickness of each frame, upgrade solution vector.x 1=h 0+Δh 1,x 2=x 1+Δh 2,x 7=x 6+Δh 7
Step 4.3.4: if the solution vector after upgrading can satisfy constraints, preserve current solution vector, finish to repair strategy; Otherwise, enter step 4.3.5.
Step 4.3.5: make i=7.
Step 4.3.6: the component x that checks solution vector iWhether satisfy constraints.If x iDo not satisfy constraint, enter step 4.3.7, otherwise, step 4.3.9 entered.
Step 4.3.7:x i=x i-0.1.
Step 4.3.8: if x i≤ x I+1, change step 4.3.11 over to; If x iSatisfy constraint, then keep x i, enter step 4.3.9; Otherwise, continue step 4.3.7.
Step 4.3.9: if i=1 preserves current solution vector, finish to repair strategy; Otherwise, enter step 4.3.10.
Step 4.3.10:i=i-1 enters step 4.3.6.
Step 4.3.11: repairing failure, add a big penalty value to object function, abandon current solution.
Step 4.4: be each the individual fitness value that calculates in the population.The present invention adopts the control target as the fitness value function of adaptive differential descent algorithm, according to formula (6), calculates fitness value
Step 4.5: check the algorithm end condition.See whether fitness value F reaches minimum, and satisfy the constraints that the present invention considers.If satisfy the algorithm end condition, then stop and export the individuality of fitness value optimum in the population; Otherwise, continue next step.
Step 4.6: the population to algorithm carries out mutation operation.For individual X I, G: i=1,2 ..., 60, one new individual V I, G+1But through type (27) produces:
V i , G + 1 = X r 1 , G + F * ( X r 2 , G - X r 3 , G ) - - - ( 27 )
Here r 1, r 2And r 3Be to go up the mutually different integer of picked at random from interval [1, N], and be different from subscript index i, mutagenic factor F value is between [0,1].
Step 4.7: the population to algorithm carries out interlace operation.Interlace operation as the formula (28).
v i , j G + 1 = u i , j G + 1 , if ( ( randb ( j ) &le; 0.4 + g 1000 ( 1 - 0.4 &CenterDot; F ( X i ) F ( U i ) + F ( X i ) ) ) or ( j = mbr ( i ) ) ) x i , j G , if ( ( randb ( j ) > 0.4 + g 1000 ( 1 - 0.4 &CenterDot; F ( X i ) F ( U i ) + F ( X i ) ) ) and ( j &NotEqual; mbr ( i ) ) )
i=1,2,…,60,j=1,2,…,7 (28)
Step 4.8: operation is selected in the filial generation in the population.The selection operation of differential evolution algorithm is between parent colony and progeny population, the individual and man-to-man competition of offspring individual of parent.The selection operation of algorithm as the formula (30).
X i G + 1 = V i G + 1 , if ( F ( V i G + 1 ) &le; F ( X i G ) ) X i G , otherwise , i=1,2,…,60
(30)
Step 4.7: return step 4.2.
What the exit thickness of step 5, each frame of obtaining according to step 4 was determined each frame wears tape speed v, temperature T, roll-force P, rolling power N, total energy consumption F.
Step 4: judge whether roll-force, roll torque, rolling power exceed the rated value of machine, exceed, and then repeating step 4, otherwise judge whether total energy consumption reaches minimum, is then to export end value, otherwise repeating step 4, up to reaching minimum of a value.The sharing of load result of each frame that empirical method and the inventive method obtain is as shown in table 3.
Table 3 empirical method and the inventive method comparative result
From table 3 and Fig. 3, more reasonable than the sharing of load that obtains with empirical method with the inventive method, the general power that consumes also descends to some extent, the object function that obtains with empirical method is 20992.55KW, object function with the inventive method is 17825KW, save energy consumption 3167.55KW, thereby method of the present invention is adopted in explanation, can reach the purpose of saving energy consumption.
Embodiment 2
If the milling train of hot fine rolling unit is not seven frames, the present invention also is suitable for.Be example with six frame hot fine rolling unit production lines.Concrete steps change n into 6 with embodiment 1, each frame exit thickness, rolling power, roll-force, object function such as table 4 that sharing of load obtains, and the rolling power contrast that empirical method and the inventive method obtain is as shown in Figure 4.
Table 4 empirical method and improved difference algorithm comparative result
From table 4 and Fig. 4, more reasonable than the sharing of load that obtains with empirical method with the inventive method, the general power that consumes also descends to some extent, the object function that empirical method obtains is 22457.7KW, object function with the inventive method is 18189.6KW, saves energy consumption 4268.1KW, and energy consumption reduces by 19%, thereby illustrate and adopt method of the present invention, can reach the purpose of saving energy consumption.

Claims (4)

1. a frame load distribution method that reduces the strip hot rolling process energy consumption is characterized in that: specifically comprise the steps:
Step 1: determine the constraints of initial control scheme, comprising:
(1) bite condition
Bite condition guarantees nipping smoothly of band steel, and rolling each reduction in pass should allow drafts less than maximum, and computing formula as the formula (1)
&Delta; h i < &Delta;h max = D ( 1 - 1 1 + f 2 ) , i = 1,2 , . . . , n - - - ( 1 )
△ h wherein iThe drafts of representing i frame; △ h MaxFor maximum allows drafts; F is the coefficient of friction between roll and rolled piece; D represents the diameter of roll;
(2) roll strength condition
For guaranteeing roll strength, draught pressure and roll torque are limited, guarantee can not surpass maximum roll-force and roll torque in the operation of rolling, computing formula is suc as formula shown in (2), the formula (3):
P i<P imax i=1,2,...,n (2)
M i<M imax i=1,2,...,n (3)
P wherein iIt is the roll-force of i frame; P ImaxIt is the maximum rolling force that i frame allows; M iIt is the roll torque of i frame; M ImaxIt is the maximum rolling force square that i frame allows;
(3) motor ability condition
The motor ability is the restrictive condition of motor overload, and rolling power must not surpass its rated power, and computing formula as the formula (4)
N i<N imax i=1,2,...,n (4)
N wherein iIt is the rolling power of i frame; N ImaxIt is the maximum rolling power that i frame allows;
(4) strip shape quality restrictive condition
In order to keep strip shape quality to reach customer requirement, optimize the drafts of arranging several of backs, making between its corresponding roll-force has certain ratio, prevents from being with steel the limit wave to occur, and computing formula is as the formula (5)
- 40 ( h i b ) 1.86 < &Delta;CR i < 80 ( h i b ) 1.86 , i = n - 3 , n - 2 , n - 1 , n - - - ( 5 )
△ CR wherein iEntrance convexity when passing through i frame for slab is poor with the outlet convexity; h iIt is the exit thickness of i frame; B is the width of slab;
Step 2: determining the target of control, is to make rolling total energy consumption reach minimum, and the energy consumption calculation formula is shown in (6) formula:
min &Sigma; i = 1 n N i , i = 1,2 , . . . , n - - - ( 6 )
Wherein: n-mm finishing mill unit frame number;
N i-the i frame power of motor;
Step 3: determine the parameter of machine parameter and the rolled piece of control scheme, machine parameter comprises running time, the rated power of the work roll diameter of mm finishing mill unit, backing roll diameter, motor, the rated speed of motor, maximum rolling force and maximum rolling force square; The parameter of rolled piece comprises the width B of rolled piece, the supplied materials thickness H of rolled piece, the finished product thickness h of rolled piece, the outlet temperature TFC of roughing mill, finish rolling outlet temperature T;
Step 4: utilize improved differential evolution algorithm to obtain each frame exit thickness;
Concrete steps are as follows:
Step 4.1: the structure initial population, and population scale and iterations are set;
Step 4.2: the individual vector of improved differential evolution algorithm is mapped to the solution vector space from real number space;
Step 4.3: for the solution vector after the decoding, if do not satisfy constraints, carry out and repair strategy;
Step 4.4: calculate each the individual fitness value in the population, with the fitness function of control target as calculating ideal adaptation degree value;
Step 4.5: check the algorithm end condition, see whether fitness value reaches minimum, and satisfy constraints; If satisfy the algorithm end condition, then stop and export the individuality of fitness value optimum in the population; Otherwise, continue next step;
Step 4.6: the population to algorithm carries out mutation operation;
For individual X I, G: i=1,2 ..., N, a new individual V I, G+1But through type (27) produces:
V i , G + 1 = X r 1 , G + F * ( X r 2 , G - X r 3 , G ) - - - ( 27 )
Wherein, r 1, r 2And r 3Be to go up the mutually different integer of picked at random from interval [1, N], and be different from subscript index i, mutagenic factor F value is between [0,1];
Step 4.7: the population to algorithm carries out interlace operation; Interlace operation as the formula (28),
v i , j G + 1 = u i , j G + 1 , if ( ( randb ( j ) &le; CR ) or ( j = mbr ( i ) ) ) x i , j G , if ( ( randb ( j ) > CR ) and ( j &NotEqual; mbr ( i ) ) ) , - - - ( 28 )
i=1,2,…,N,j=1,2,…,Dv
In the formula, Be J dimension component; Be J dimension component, Individuality for the variation generation; Be J dimension component, Be the parent individuality; Randb (j) is equally distributed probability between [0,1]; An integer that generates at random between mbr (i) expression [1, Dv], CR is crossover probability, generally gets the number between [0,2]; Here
CR = CR min + g G ( 1 - CR min &CenterDot; F ( X i ) F ( U i ) + F ( X i ) ) - - - ( 29 )
In the formula, F (X i) expression i parent individuality fitness value; F (U i) fitness value of variation individuality of i parent correspondence of expression; CR MinBe minimum crossover probability; G is current iteration algebraically; G is algorithm greatest iteration algebraically;
Step 4.8: operation is selected in the filial generation in the population;
The selection operation of differential evolution algorithm is between parent colony and progeny population, the individual and man-to-man competition of offspring individual of parent, and the selection of algorithm is operated as the formula (30),
X i G + 1 = V i G + 1 , if ( F ( V i G + 1 ) &le; F ( X i G ) ) X i G , otherwise , i = 1,2 , &CenterDot; &CenterDot; &CenterDot; , N - - - ( 30 )
Step 4.9: return step 4.2;
Step 5: what the exit thickness of each frame that obtains according to step 4 was determined each frame wears tape speed, temperature, roll-force, rolling power, total energy consumption;
Step 6: judge whether roll-force, roll torque, rolling power exceed the rated value of machine, exceed, and then repeating step 4, otherwise judge whether total energy consumption reaches minimum, is then to export end value, otherwise repeating step 4, up to reaching minimum of a value.
2. the frame load distribution method of reduction strip hot rolling process energy consumption according to claim 1 is characterized in that:
The step that the described initial population of step 4.1 produces is as follows:
If have N n dimension individual in the population, each individual vectorial X iRepresentation be
Step 4.1.1: make x 01=500, x 0n=-500,
Step 4.1.2: utilize formula (24) to calculate x 0j, obtain X 0, X 0=[x 01, x 02..., x 0n];
x 0j=x 0j-Slot,j=2,3,…,n-1 (24)
Step 4.1.3: make i=1;
Step 4.1.4: utilize formula (25), with X 0Be template, produce the individual X in the population i, i=1,2 ..., N, x I1=500, x In=-500,
x ij=x 0j+f·Slot,j=2,3,…,n-1 (25)
In the formula, f is random number, f ∈ (1,1);
Step 4.1.5:i=i+1;
Step 4.1.6: if i=N, then the initial population structure finishes; Otherwise, change step 4.1.4 over to.
3. the frame load distribution method of reduction strip hot rolling process energy consumption according to claim 1 is characterized in that:
The individual vector of described step 4.2 evolution algorithm is mapped to the solution vector space from real number space, and concrete steps are as follows;
Step 4.2.1: find out maximum and the minimum of a value of each component among the individual vectorial X, be designated as x MaxAnd x Min, make v Max=h 0, v Min=h n, v wherein MaxAnd v MinBe respectively x MaxWith x MinThe corresponding component in decoding back; h 0Be mm finishing mill unit strip steel at entry thickness; h nBe mm finishing mill unit outlet belt steel thickness; For other component among the X, be mapped to solution space with formula (26):
v i = v max - ( x max - x i ) ( x max - x min ) * ( v max - v min ) , i = 1,2 , &CenterDot; &CenterDot; &CenterDot; , n ,
i≠max,min (26)
Step 4.2.2: with the component (v of decoded intermediate vector V 1, v 2..., v n) by arranging from big to small;
Step 4.2.3: with the component v after arranging iUpgrade solution vector U:u 1=v Max..., u n=v Min
4. the frame load distribution method of reduction strip hot rolling process energy consumption according to claim 1 is characterized in that:
The described reparation strategy of step 4.3, concrete steps are as follows:
Step 4.3.1: if decoded solution vector does not satisfy constraints, then enter step 4.3.2; Otherwise, preserve current solution vector, finish to repair strategy;
Step 4.3.2: the drafts △ h that calculates each frame of solution vector correspondence i, and these drafts are arranged from big to small: △ h 1=△ h Max..., △ h n=△ h Min
Step 4.3.3: the drafts of redistributing for each frame, and recomputate the exit thickness of each frame, upgrade solution vector; x 1=h 0+ △ h 1, x 2=x 1+ △ h 2..., x n=x N-1+ △ h n
Step 4.3.4: if the solution vector after upgrading can satisfy constraints, preserve current solution vector, finish to repair strategy; Otherwise, enter step 4.3.5;
Step 4.3.5: make i=n;
Step 4.3.6: the component x that checks solution vector iWhether satisfy constraints, if x iDo not satisfy constraint, enter step 4.3.7, otherwise, step 4.3.9 entered;
Step 4.3.7:x i=x i-0.1, i=1,2 ..., n;
Step 4.3.8: if x i≤ x I+1(i=1,2 ..., n), change step 4.3.11 over to; If x iSatisfy constraint, then keep x i, enter step 4.3.9; Otherwise, continue step 4.3.7;
Step 4.3.9: if i=1 preserves current solution vector, finish to repair strategy; Otherwise, enter step 4.3.10;
Step 4.3.10:i=i-1 enters step 4.3.6;
Step 4.3.11: repairing failure, add a big penalty value to object function, abandon current solution.
CN 201110390910 2011-11-30 2011-11-30 Machine frame load distribution method for decreasing energy consumption of rolling process of hot rolled strip steel Active CN102489524B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110390910 CN102489524B (en) 2011-11-30 2011-11-30 Machine frame load distribution method for decreasing energy consumption of rolling process of hot rolled strip steel

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110390910 CN102489524B (en) 2011-11-30 2011-11-30 Machine frame load distribution method for decreasing energy consumption of rolling process of hot rolled strip steel

Publications (2)

Publication Number Publication Date
CN102489524A CN102489524A (en) 2012-06-13
CN102489524B true CN102489524B (en) 2013-09-04

Family

ID=46181408

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110390910 Active CN102489524B (en) 2011-11-30 2011-11-30 Machine frame load distribution method for decreasing energy consumption of rolling process of hot rolled strip steel

Country Status (1)

Country Link
CN (1) CN102489524B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103272853B (en) * 2013-05-22 2015-01-28 沈阳工业大学 Device and method for setting rolling reduction and rolling speed of each rack in cold continuous rolling
CN105013832A (en) * 2014-04-28 2015-11-04 宝山钢铁股份有限公司 Hot rolled strip steel load distribution method giving consideration to rolling energy consumption and good strip shape
CN106269911B (en) * 2016-08-29 2018-10-23 首钢京唐钢铁联合有限责任公司 A kind of roughing pressure Control System of Load Distribution method and roughing control system
CN109365544B (en) * 2018-09-05 2020-02-21 湖南华菱涟源钢铁有限公司 Load distribution method of reversible single-stand four-roll roughing mill for improving rolling rhythm

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000280015A (en) * 1999-03-31 2000-10-10 Kawasaki Steel Corp Method and equipment to control meandering of hot- rolled sheet steel strip
JP2003001311A (en) * 2001-06-21 2003-01-07 Nisshin Steel Co Ltd Plate width control method in cold tandem rolling
CN101337237A (en) * 2008-08-14 2009-01-07 苏州有色金属研究院有限公司 Thickness automatic allocation process of single-frame double coiling-aluminum hot milling roll
CN101372018A (en) * 2008-10-17 2009-02-25 山西太钢不锈钢股份有限公司 Diversified cross-connection control method for plate shape of hot rolling band steel

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS641208B2 (en) * 1979-12-27 1989-01-10 Mitsubishi Electric Corp

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000280015A (en) * 1999-03-31 2000-10-10 Kawasaki Steel Corp Method and equipment to control meandering of hot- rolled sheet steel strip
JP2003001311A (en) * 2001-06-21 2003-01-07 Nisshin Steel Co Ltd Plate width control method in cold tandem rolling
CN101337237A (en) * 2008-08-14 2009-01-07 苏州有色金属研究院有限公司 Thickness automatic allocation process of single-frame double coiling-aluminum hot milling roll
CN101372018A (en) * 2008-10-17 2009-02-25 山西太钢不锈钢股份有限公司 Diversified cross-connection control method for plate shape of hot rolling band steel

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
双机架热轧中厚板轧机负荷分配的工艺优化研究;路义山等;《青海大学学报(自然科学版)》;20050425(第02期);第1-4页 *
路义山等.双机架热轧中厚板轧机负荷分配的工艺优化研究.《青海大学学报(自然科学版)》.2005,(第02期),第1-4页.

Also Published As

Publication number Publication date
CN102489524A (en) 2012-06-13

Similar Documents

Publication Publication Date Title
CN102513351B (en) Rolling method and device for strip steel tandem cold rolling
CN103586286B (en) Rolling schedule comprehensive optimization method for cold continuous rolling unit taking scratch prevention as objective
CN102489524B (en) Machine frame load distribution method for decreasing energy consumption of rolling process of hot rolled strip steel
CN104858241B (en) A kind of emulsion flow comprehensive optimization method of tandem mills paper-thin strip rolling
CN1301810C (en) Comprehensive optimized control method of rolling stardard for cold band-steel continuous rolling mill
CN103544340B (en) The establishing method of concentration of emulsion used in five Stands Cold Tandem Mill group strip in razor-thin rollings
CN105312321A (en) Method for optimizing technological lubrication system of cold continuous rolling unit
CN103611732B (en) The technological lubrication system optimization method that tandem mills is target with plucking control
CN105512804B (en) Cold continuous rolling process is using cost comprehensively control as the emulsion flow set method of target
CN103357670B (en) Reduction schedule optimizing method applicable to five-stand UCM (universal crown mill) type cold continuous rolling unit
CN102266865A (en) Hot/cold rolling load distribution method
CN103372575B (en) Pickling and rolling mill load distribution method
CN107716560B (en) A kind of Hot Strip Rolling load distribution method
CN102921743B (en) Method for determining five stand tandem cold mill depressing distribution
CN105013832A (en) Hot rolled strip steel load distribution method giving consideration to rolling energy consumption and good strip shape
CN104942019B (en) A kind of cold rolling of strip steel process Automatic control method of width
CN101507978B (en) Worker roller abrasion forecasting method of cold continuous rolling device based on mechanism and operation condition
CN101658871A (en) Optimization method of rolling schedule of non-reversible aluminum strip cold rolling mill
CN101648208A (en) Rolling procedure making method of aluminium single-stand cold-rolling machine
CN104785540B (en) A kind of rolling efficiency method for improving for being suitable for five Stands Cold Tandem Mill groups
CN105234188B (en) Mill speed optimization method in cold continuous rolling process with Result coutrolling as target
CN1091008C (en) Interlinked control method for plate-band rolling course based on coordination law of plate shape and plate thickness
Hu et al. Multi-parameter deep-perception and many-objective autonomous-control of rolling schedule on high speed cold tandem mill
CN104923571B (en) A kind of hot-rolling plate belt steel rolling code dynamic setting method
CN107442577A (en) A kind of fine-rolling strip steel sharing of load establishing method

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
C14 Grant of patent or utility model