CN101623708B - Plate-shape control integrated system and executing method thereof - Google Patents

Plate-shape control integrated system and executing method thereof Download PDF

Info

Publication number
CN101623708B
CN101623708B CN2009100750825A CN200910075082A CN101623708B CN 101623708 B CN101623708 B CN 101623708B CN 2009100750825 A CN2009100750825 A CN 2009100750825A CN 200910075082 A CN200910075082 A CN 200910075082A CN 101623708 B CN101623708 B CN 101623708B
Authority
CN
China
Prior art keywords
frame
roll
cross direction
direction profiles
rolling
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
CN2009100750825A
Other languages
Chinese (zh)
Other versions
CN101623708A (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.)
Yanshan University
Original Assignee
Yanshan University
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 Yanshan University filed Critical Yanshan University
Priority to CN2009100750825A priority Critical patent/CN101623708B/en
Publication of CN101623708A publication Critical patent/CN101623708A/en
Application granted granted Critical
Publication of CN101623708B publication Critical patent/CN101623708B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Metal Rolling (AREA)

Abstract

The invention provides an integrated method of a plate-shape control system of a plate-strip cold and hot continuous rolling machine. The plate-shape control system comprises a rolled piece plastic-distortion module, a roller-system elastic-distortion module, a rolled piece and roller temperature field module, a roller-system abrasion module, a flatness well-recognition module, a plate-shaped mode recognition module, a plate-shape standard curve module and a plate-shape control module. The invention integrates all modules according to internal relations, can solve a plurality of plate-shape control problems comprehensively and optimally, adopts a set of calculating flows in actual application and effects saves the labor time.

Description

Plate-shape control integrated system and manner of execution
Technical field
The present invention relates to a kind of analysis and control technology of plate strip rolling process, the manner of execution of particularly a kind of plate-shape control integrated system and this integrated system.
Background technology
Strip is the major product of steel and iron industry, and " elaboration " overwhelming majority in the steel products belongs to the strip product.Strip is than being the important symbol of weighing a national metallurgical industry level.The strip of advanced country in the world reaches 65-70% than at present, exceeds China more than 20%, so the strip of China is produced the space of developing on a large scale very much in addition.Plate shape is the important quality index of strip material, and plate shape control technology is the core technology of strip-mill strip.Existing strip-mill strip tens covers of China, plate shape control technology integral level is compared backward relatively with international most advanced level.The bottleneck that restraining plate shape control technology improves is that with crucial the research of plat control system lags behind the development of technology, lacks advanced and applicable plat control system.
People have carried out a large amount of theoretical researches to the every field that the control of plate shape relates at present.For example: " the stream ready-to-eat type noodle method of simulation board rolling 3 D deformation " (" mechanical engineering journal " 2003,39 the 7th phases of volume: 94-100) plastic deformation is studied to the plate strip rolling process rolled piece; " exploitation of four-high mill roll strain parsing module " (" steel rolling " 2003,20 the 2nd phases of volume: 8-11) roll elastic deformation is studied; " simulation of Temperature in the hot continuous rolling process " (" steel research journal " 2006,18 the 8th phases of volume: 32-34) the rolled piece temperature field is studied; " 1700 continuous hot-rolling mill roller temperature fields and hot convexity research " (" Northeastern University's journal " natural science edition, 2008,29 the 4th phases of volume: 517-520) breaker roll temperature field and thermal deformation are studied; (" iron and steel " 2002,37 the 3rd phases of volume: 24-27) wearing and tearing of pair roller system are studied " researchs of 2050 CVC continuous hot-rolling mill mm finishing mill unit roll wears "; " analysis of cold-rolled strip steel shape discrimination model and discussion " (" iron and steel " 1995,30 the 8th phases of volume: 39-43) good differentiation is studied to glacing flatness; " 6 roller cold continuous rolling types based on plate shape control ability are selected " (" Zhongnan Polytechnic Univ's journal " English edition, 2007,14 the 2nd phases of volume: 278-284) rolled piece plastic deformation, roll elastic deformation have been carried out Coupling Research.
In general, at present the disclosed document individual module that is confined to the control of plate shape is related to or to the Coupling Research of several modules about comprehensive collective method of plat control system, is not also seen relevant report in the document of publishing.Thereby can only analyze and solve simple relatively plate shape problem, can not comprehensively, optimally analyze and solve highly difficult plate shape problem.Set control problem such as plate shape; The general at present method that adopts is to carry out rolled piece plastic deformation module and roll elastic deformation module integrated; Be optimized calculating according to a certain target then, this method generally can be dealt with problems to a certain extent, but can not obtain optimal solution.This is owing to do not consider roll hot-rolling type and roll wear roll shape during computation optimization; The consequence that causes is: machine initial setting control effect is better on the roll; Along with constantly carrying out of the operation of rolling; The comprehensive roll shape of roll constantly changes, and causes production practices and perfect condition to depart from gradually, and the effect of initial plate shape setup parameter is variation gradually just also.Thereby have only each module that the control of plate shape is related to carry out comprehensive integration, and form complete plat control system, could optimally solve the various plate shape control problems that run in the production practices comprehensively.
Summary of the invention
In order to solve the problems of the technologies described above; The present invention provides a kind of plate-shape control integrated system and manner of execution; This system and method can be comprehensively, optimum, solve various common plate shape control problems easily; Form the plurality of plate-shape control technology, specifically comprise: analysis of plate shape control performance and type selection technology, roll shape optimisation technique, plate shape are set control technology, plate shape FEEDFORWARD CONTROL technology, plate shape feedback control technology.
The technical scheme that the present invention adopted is: a kind of plate-shape control integrated system, and this system comprises 8 modules that are mutually related:
Rolled piece plastic deformation module, roll elastic deformation module, rolled piece and roller temperature field module, roller system wearing and tearing module, the good discrimination module of glacing flatness, plate shape pattern recognition module, plate-shape standard curve module, plate shape control module; Wherein,
Said rolled piece plastic deformation module is that roll elastic deformation module and roll wear module provide the draught pressure cross direction profiles, distributes for the good discrimination module of glacing flatness and plate shape pattern recognition module provide flatness and lateral;
Said roll elastic deformation module is that rolled piece plastic deformation module and plate shape pattern recognition module provide band steel exports thickness cross direction profiles, for the roll wear module provides the roll gap pressure cross direction profiles;
Said rolled piece and roller temperature field module are that the calculating of resistance of deformation in the rolled piece plastic deformation module provides the deformed area mean temperature, for the calculating of work roll thermal roll shape in the roll elastic deformation module provides the working roll temperature field;
Said roll wear module is that the calculating of comprehensive roll shape in the roll elastic deformation module provides the roll wear roll shape;
The local integrated formation plate-shape standard curve module of good discrimination module of said glacing flatness and rolled piece plastic deformation module;
Said plate shape control module is calculated plate shape control device regulated quantity according to the output result of plate shape pattern recognition module.
Said rolled piece plastic deformation module adopts bar unit method.According to rolling condition; Consider the lateral flow of rolled piece; Set up deformed area Three-dimensional Flow velocity field, strain field, stress field; Confirm the cross direction profiles of draught pressure and outlet glacing flatness, for roll elastic deformation module, roller system wearing and tearing module, the good discrimination module of glacing flatness and plate shape pattern recognition module provide design conditions.Specifically can be divided into streamline bar unit method, stream ready-to-eat type noodle method, the bar unit calculus of variations again.Wherein streamline bar unit method is applicable to cold rolling calculated off line, and stream ready-to-eat type noodle method is applicable to the hot rolling calculated off line, and bar unit calculus of variations computational speed is very fast, is applicable in line computation.Adopt Mathematical Modeling to represent as follows:
X i=f 11si,μ i,h i-1,h i,B,T i-1,T i) (1a)
X i=f 12si,μ i,h i-1,h i,B,T i-1,T i) (1b)
X i=f 13si,μ i,h i-1,h i,B,T i-1,T i) (1c)
X i=(σ 1i,σ 0i,p i) (1d)
Wherein, i is frame number (i=1~n), f 11Be streamline bar unit method computing formula, f 12Be stream ready-to-eat type noodle method computing formula, f 13Be bar unit calculus of variations computing formula, σ 1i, σ 0iBe respectively the forward and backward tensile stress cross direction profiles of i frame value, p iBe i frame draught pressure cross direction profiles value, σ SiBe i frame band steel average deformation drag, μ iBe i frame coefficient of friction, h I-1, h iBe respectively that i frame band steel is gone into, exit thickness cross direction profiles value, B is band steel supplied materials width, T i, T I-1Be respectively the forward and backward tension force of i frame;
Each frame band steel average deformation drag σ SiBy steel grade, drafts, mill speed and the decision of deformed area band mean temperature, adopt Mathematical Modeling to represent as follows:
σ si=f 14(s s,h i-1,h i,v 0,h 0,T swi) (1e)
Wherein, f 14Be resistance of deformation computing formula, s sBe steel grade, v 0Be the first frame porch window of web velocities, T SwiBe deformed area band ensemble average temperature.
Said roll elastic deformation module adopts influence coefficient method.Calculate amount of deflection, the flattening amount of roll according to Elasticity, finally define and carry a roll gap shape (rolled piece exit thickness cross direction profiles value) and a roll gap pressure cross direction profiles value, for rolled piece plastic deformation module and roller system wearing and tearing module provide design conditions.This method can comprise various types common in the present production practices: normal four-roller milling train, CVC series rolling mill, HC series rolling mill and PC milling train carry out analytical calculation.Adopt Mathematical Modeling to represent as follows:
Y i=f 21(L wi,L mi,L bi,L 1i,L 2i,L 3i,D wi,D mi,D bi,D bni,C wi,C mi,C bi,S 1i,S 2i,RS wi,RS mi,Tilt i,p i,h i)(2a)
Z i=f 22(L wi,L bi,L 1i,L 3i,D wi,D bi,D bni,C wi,C bi,S 1i,RS wi,Tilt i,θ i,p i,h i) (2b)
Y i=(h i,q wmi,q mbi) (2c)
Z i=(h i,q wbi) (2d)
C wi=C wgi+C wti-W wi (2e)
C mi=C mgi-W mi (2f)
C bi=C bgi-W bi (2g)
Wherein, i is frame number (i=1~n), f 21Be six-high cluster mill roll elastic deformation computing formula, f 22Be four-high mill roll elastic deformation computing formula, q WmiBe i frame working roll and intermediate calender rolls roll gap pressure cross direction profiles, q MbiBe i frame intermediate calender rolls and backing roll roll gap pressure cross direction profiles, q WbiBe i frame working roll and backing roll roll gap pressure cross direction profiles, L WiBe i frame working roll barrel length, L MiBe i frame intermediate calender rolls barrel length, L BiBe i frame support roller barrel length, L 1iBe i frame work roll bending cylinder spacing, L 2iBe i frame intermediate calender rolls roller cylinder spacing, L 3iBe i frame depress oil cylinder centre-to-centre spacing, D WiBe i frame work roll diameter, D MiBe i frame intermediate calender rolls diameter, D BiBe i frame support roller diameter, D BniBe i frame support roller roll neck diameter, C WiBe the comprehensive roll shape of i frame working roll, C WgiBe the original grinding roll shape of i frame working roll, C WtiBe i frame work roll thermal roll shape, W WiBe i frame working roller abrasion roll shape, C MiBe the comprehensive roll shape of i frame intermediate calender rolls, C MgiBe the original grinding roll shape of i frame intermediate calender rolls, W MiBe i frame intermediate calender rolls wearing and tearing roll shapes, C BiBe the comprehensive roll shape of i frame support roller, C BgiBe the original grinding roll shape of i frame support roller, W BiBe i frame support roller wearing and tearing roll shapes, S 1iBe i frame work roll bending power, S 2iBe i frame intermediate calender rolls bending roller force, RS WiBe i frame working roll roll shifting amount, RS MiBe i frame intermediate roll shifting amount, Tilt iBe the i frame roller amount of inclining, θ iBe the i breast roller angle of the crossing, h iIt is i frame band steel exports average thickness;
Roll hot-rolling type in the formula (2e) can calculate according to iron Mo Xinke formula, adopts Mathematical Modeling to represent as follows:
C wti j = f 23 ( T wi j , T wi , D wi ) - - - ( 2 h )
Wherein, f 23Be iron Mo Xinke formula, C Wti jI frame work roll thermal roll shape during for rolling j+1 winding steel, T WiBe temperature field during machine on the i frame working roll, T Wi jI frame working roll temperature field during for rolling j+1 winding steel.
Said rolled piece and roller temperature field module adopt finite difference calculus.According to heat transfer principle, because length of rolled piece, is therefore ignored the heat conduction on the length of rolled piece direction, the rolled piece temperature field of adopting the two-dimensional finite calculus of finite differences to calculate rolling district and non-rolling district than rolled piece width and the big several magnitude of rolled piece thickness.Finally confirm each frame rolling deformation district band ensemble average temperature T through integration method Swi, average surface temperature T Ssi, outlet band transverse temperature distribution curve T Δ i, for rolled piece plastic deformation module (1) provides design conditions.Adopt Mathematical Modeling to represent as follows:
T=f 31(T s,v 0,B,h 0,…,h n,μ 1,…,μ n,σ s1,…,σ sn,T ws1,…,T wsn,T c1,…,T cn,T a,U)(x>0) (3a)
T=f 32(T s,v 0,B,h 0,…,h n,μ 1,…,μ n,σ s1,…,σ sn,T ws1,…,T wsn,T c1,…,T cn,T a,V)(x>0) (3b)
U=(L M12,…,L M(n-1)n) (3c)
V=(L M12,…,L M(n-1)n,L T0,L T1,…,L Tn,L C0,L C1,…,L Cn,Bu 1,…,Bu n) (3d)
Wherein, T is the belt steel temperature field, f 31Be divided into cold continuous rolling rolled piece temperature field computing formula, f 32Be continuous hot-rolling mill rolled piece temperature field computing formula, T s(y z) is the two-dimensional temperature field of x=0 place rolled piece, v 0Be supplied materials speed, B is the supplied materials width, h 0~h nBe band steel supplied materials average thickness and each frame outlet average thickness, μ 1~μ nBe each frame coefficient of friction, σ S1~σ SnBe each frame deformation of rolled wire drag, T Ws1~T WsnBe each breast roller average surface temperature, T C1~T CnBe each frame cooling water or emulsion temperature, T aBe air themperature, L M12~L M (n-1) nBe frame spacing, L T0Be De-scaling box center to the first frame central spacing, L T1~L TnThe water-cooled node is to this frame central spacing, L after each frame C0For spraying, De-scaling box covers siding-to-siding block length, L C1~L CnThe water-cooled node sprays and covers siding-to-siding block length, Bu after each frame 1~Bu nFor describing the state variable that whether the water-cooled node sprays water after each frame ( value 0 or 1 is got 0 expression and do not existed or do not spray water, and gets the normal water spray of 1 expression);
When high-speed rolling; Work roll surface alternately contacts with rolled piece, air, cooling water (or emulsion) and intermediate calender rolls (or backing roll); The revolution lap time is very short; Only on working roll superficial layer as thin as a wafer along the circumferential direction temperature have significant change, and roll circumferential direction inside direction variations in temperature is minimum, can it be summed up as a Dynamic Two-dimensional heat conduction problem.For the working roll of every frame, all set up local coordinate system as shown in Figure 4, adopt finite difference calculus to calculate the roller temperature field, represent as follows:
T wi j = f 41 ( T wi j - 1 , T ‾ ssi j , t r j , t p j ) T wi 0 = T wi - - - ( 4 )
Wherein, i is that (i=1~n), j is rolling coil of strip number to the frame number, T Wi jBe the intact back of j winding steel rolling t p jMoment i frame working roll temperature field, t r j, t p jBeing respectively the pure of rolling j winding steel rolls the time and rolls back intermittent time, T WiTemperature during machine on the i frame working roll, T Ssi jEach frame rolling deformation district strip surface mean temperature during for rolling j winding steel.
Said roller is that the wearing and tearing module obtains according to tribology principle.Promptly between two objects wear extent with normal pressure and wearing and tearing distance each other are directly proportional each other.Adopt Mathematical Modeling to represent as follows:
W wi j = f 51 w ( W wi j - 1 , L s j , h 0 j , h i j , p i j , q wmi j ) W wi 0 = 0 - - - ( 5 a )
W mi j = f 51 m ( W mi j - 1 , L s j , h 0 j , h i j , q wmi j , q mbi j ) W mi 0 = 0 - - - ( 5 b )
W bi j = f 51 b ( W bi j - 1 , L s j , h 0 j , h i j , q mbi j ) W bi 0 = 0 - - - ( 5 c )
W wi j = f 52 w ( W wi j - 1 , L s j , h 0 j , h i j , p i j , q wbi j ) W wi 0 = 0 - - - ( 5 d )
W bi j = f 52 b ( W bi j - 1 , L s j , h 0 j , h i j , q wbi j ) W bi 0 = 0 - - - ( 5 e )
Wherein, i is that (i=1~n), j is rolling coil of strip number to the frame number, f 51wBe six-high cluster mill working roll roll wear computing formula, f 51mBe six-high cluster mill intermediate calender rolls roll wear computing formula, f 51bBe six-high cluster mill backing roll roll wear computing formula, f 52wBe four-high mill working roller abrasion computing formula, f 52bBe four-high mill backing roll roll wear computing formula, W Wi jBe the intact back of j winding steel rolling i frame working roller abrasion roll shape, W Mi jBe the intact back of j winding steel rolling i frame intermediate calender rolls wearing and tearing roll shapes, W Bi jBe the intact back of j winding steel rolling i frame support roller wearing and tearing roll shapes, L s jBe j winding steel supplied materials length.
The good discrimination module of said glacing flatness adopts the differentiation factorization method of cutting apart based on bar unit.The band steel exports flatness and lateral that calculates according to rolled piece plastic deformation module distributes, and adopts a plate shape to differentiate the factor and judges plate shape quality, need not the certain unstability pattern of hypothesis before plate shape is differentiated, and has simplified computational methods.Adopt Mathematical Modeling to represent as follows:
ξ i=f 6(B,h i,σ 1i) (6)
Wherein, i is frame number (i=1~n), f 6Be the differentiation factorization method of cutting apart based on bar unit, ξ iI frame band steel exports glacing flatness is differentiated the factor, ξ iThe band steel was straight in>1 o'clock, ξ iThe band steel was in critical instability status, ξ in=1 o'clock i<1 o'clock band steel unstability.
Said plate shape pattern recognition module adopts the least square method based on Legnedre polynomial.This method distributes according to the plate shape of rolling back band steel, it is carried out pattern classification, so that the flatness defect of different mode is taked the Different control measure.Adopt Mathematical Modeling to represent as follows:
A i=f 71i) (7a)
E i=f 7(h i) (7b)
A i=(a 1i,a 2i,a 3i,a 4i) (7c)
E i=(e 1i,e 2i,e 3i,e 4i) (7d)
Wherein, i is frame number (i=1~n), f 7Be the least square method based on Legnedre polynomial, a 1i~a 4iBe respectively the i frame and export strip profile and flatness cross direction profiles 1 time, 2 times, 3 times, 4 sub-eigenvalues, e 1i~e 4iBe respectively the i frame and export belt steel thickness cross direction profiles 1 time, 2 times, 3 times, 4 sub-eigenvalues.
Said plate-shape standard curve module adopts Mathematical Modeling to represent as follows:
Γ ‾ = f 8 ( σ 1 n ob , h n ob , B , h 0 , h ‾ 0 , h ‾ 1 , · · · , h ‾ n )
H ‾ = f 8 ( σ 1 n ob , h n ob , B , h 0 , h ‾ 0 , h ‾ 1 , · · · , h ‾ n )
Γ ‾ = ( σ 11 ob , σ 12 ob , · · · σ 1 n ob )
H ‾ = ( h 1 ob , h 2 ob , · · · h n ob )
Wherein, f 8Be quantitative resolution method, σ 1n ObBe the tandem mill end frame outlet band target glacing flatness curve of industrial requirements decision, h n ObThe tandem mill end frame outlet band target thickness cross direction profiles curve of industrial requirements decision, σ 11 Ob~σ 1n ObBe respectively each frame export goal glacing flatness curve, h 1 Ob~h n ObBe respectively each frame export goal thickness cross direction profiles curve.
Said plate shape control module adopts the influence matrix method.This method connects each time plate shape component and various plate shape control device through plate shape gating matrix; Only need during practical application each component of degree n n of plate shape deviation that detection obtains multiply by plate shape control influence matrix; Can obtain optimum plate shape control device regulated quantity; The conventional searching algorithm of computational accuracy and computational speed all is significantly improved, and adopts Mathematical Modeling to represent as follows:
K n=f 91(Δa 1n,Δa 2n,Δa 3n,Δa 4n) (9a)
K i=f 92(Δe 10,Δe 20,Δe 30,Δe 40)(i=1~n) (9b)
K i=f 93(Δe 10,Δe 20,Δe 30,Δe 40)(i=1~n) (9c)
Wherein, f 91Be the plate shape control device compensation rate computing formula that last frame outlet actual measurement glacing flatness and target glacing flatness deviation cause, f 92For guaranteeing that each frame exports actual glacing flatness curve and equals under the prerequisite of target glacing flatness curve; Front section actual measurement thickness cross direction profiles worked as by supplied materials and the supplied materials head is surveyed each frame plate shape control device compensation rate computing formula that thickness cross direction profiles deviation causes; Be applicable to cold continuous rolling, f 93For guaranteeing that the actual exit thickness cross direction profiles of each frame equals under the prerequisite of target thickness cross direction profiles; Front section actual measurement thickness cross direction profiles worked as by supplied materials and the supplied materials head is surveyed each frame plate shape control device compensation rate computing formula that thickness cross direction profiles deviation causes; Be applicable to continuous hot-rolling mill, Δ a 1n, Δ a 2n, Δ a 3n, Δ a 4nBe respectively last frame outlet poor when front section actual measurement glacing flatness curve 1 time, 2 times, 3 times, 4 sub-eigenvalues and target glacing flatness curve 1 time, 2 times, 3 times, 4 sub-eigenvalues, Δ e 10, Δ e 20, Δ e 30, Δ e 40Be respectively supplied materials poor when front section actual measurement thickness cross direction profiles 1 time, 2 times, 3 times, 4 sub-eigenvalues and band steel toe portion actual measurement thickness cross direction profiles 1 time, 2 times, 3 times, 4 sub-eigenvalues, K iBe all plate shape control device regulated quantity matrixes of i frame.
A kind of manner of execution of plate-shape control integrated system, this method comprise following computer step:
Step 1) is set up the tandem mill coordinate system;
Step 2) confirm the use object of plate-shape control integrated system, the definition integer variable:
{con1,con2 i(i=1,2,…,n),con3,con4,con5},
Wherein, n is a tandem mill frame number, through this n+4 of initialization variable, confirms the use object of plate-shape control integrated system;
Step 3) is collected relevant rolling equipment and technological parameter;
Step 4) is when con4=1 and con1=3; Perhaps when con4=2 and con1=3; Perhaps when con4=3 and con1=3; Calculate the belt steel temperature field according to rolled piece and roller temperature field module, and calculate each stand stretch district average deformation drag according to rolled piece plastic deformation module simultaneously, otherwise execution in step 5);
Step 5) perhaps when con4=2, is coupled to rolled piece plastic deformation module and roll elastic deformation module when con4=1, finds the solution each frame outlet glacing flatness, thickness and draught pressure cross direction profiles, otherwise execution in step 6);
Step 6) is calculated each frame target glacing flatness curve σ according to the plate-shape standard curve module when con4=3 1i Obj(i=1~n) and target thickness cross direction profiles curve h i Obj(i=1~n), otherwise execution in step 7);
Step 7) is carried out optimizing to each frame plate shape control device and is calculated when con4=3, obtains each frame plate shape setup parameter, otherwise execution in step 8);
Step 8) is calculated rolling each the frame work roll thermal roll shape of back that finishes of current volume according to rolled piece and roller temperature field module when con4=3, calculate rolling each the frame working roller abrasion roll shape of back that finishes of current volume according to roller system wearing and tearing module, otherwise execution in step 9);
Step 9) is when con4=1; Perhaps when con4=4; Perhaps when con4=5, glacing flatness or thickness cross direction profiles curve are carried out pattern-recognition, calculate each frame plate shape control device regulated quantity matrix K according to plate shape control module according to plate shape pattern recognition module i j(i=1~n), otherwise execution in step 10);
Step 10) superposes rolling each frame hot-rolling type of back, roll wear roll shape and the original grinding roll shape of roll of finishing of current volume when con4=3.
It is specific as follows to set up the tandem mill coordinate system described in the said method in the step 1):
Set up like Fig. 2 and coordinate system shown in Figure 3 for tandem mill, n represents tandem mill frame sum among the figure.Should be noted that what provide among Fig. 2 and Fig. 3 is the coordinate-system sketch map of four-high mill, but the present invention and be not limited to four-high mill, the coordinate-system of six-high cluster mill is similarly.The length of rolled piece direction is elected to be the x axle, and the rolled piece width is elected to be the y axle, and the rolled piece short transverse is elected to be the z axle, and the origin of coordinates is chosen milling train width center at the y axle before the x axle is chosen tandem mill first rolling mill inlet, choose the rolled piece mid-depth at the z axle.
Step 2 described in the said method) the use object in comprises:
1a) be applied to hot continuous rolling when line computation, make con1=3, con3=2,
Be applied to cold continuous rolling when line computation, make con1=3, con3=1,
When being applied to the hot continuous rolling off-line analysis, make con1=2, con3=2,
When being applied to the cold continuous rolling off-line analysis, make con1=1, con3=1;
When 1b) being applied to six-high cluster mill, make con2 i=1, i=1~n, n are the frame numbers,
When being applied to four-high mill, make con2 i=2, i=1~n, n are the frame numbers;
When 1c) being applied to control performance analysis of plate shape and type selection technology, make con4=1,
When being applied to the roll shape optimisation technique, make con4=2,
When being applied to plate shape setting control technology, make con4=3,
When being applied to plate shape FEEDFORWARD CONTROL technology, make con4=4,
When being applied to plate shape feedback control technology, make con4=5.
Relevant rolling equipment of collection described in the said method in the step 3) and technological parameter may further comprise the steps:
Step 2a) when con4=1, perhaps when con4=2, perhaps when con4=3, collect the rolling equipment relevant parameter, comprising: each frame working roll barrel length L Wi, each frame work roll diameter D Wi, con2 i=1 frame intermediate calender rolls barrel length L Mi, con2 i=1 frame intermediate calender rolls diameter D Mi, each frame support roller barrel length L Bi, each frame support roller diameter D Bi, each frame support roller roll neck diameter D Bni, each frame work roll bending cylinder spacing L 1i, con2 i=1 frame intermediate calender rolls roller cylinder spacing L 2i, each frame depress oil cylinder centre-to-centre spacing L 3i, the original grinding roll shape of each frame working roll C Wgi, con2 iThe original grinding roll shape of=1 frame intermediate calender rolls C Mgi, the original grinding roll shape of each frame support roller C Bgi, otherwise execution in step 2b);
Step 2b) when con4=1 and con1=3, perhaps when con4=2 and con1=3, perhaps when con4=3 and con1=3, collect the rolling equipment relevant parameter, comprising: each frame working roll is initially gone up the machine temperature T Wi, each frame cooling water or emulsion temperature T Ci, air themperature T a, frame spacing L Mi (i+1), otherwise execution in step 2c);
Step 2c) when con4=1 and con1=3 and con3=2; Perhaps when con4=2 and con1=3 and con3=2; Perhaps when con4=3 and con1=3 and con3=2, collect the rolling equipment relevant parameter, comprising: De-scaling box center to the first frame central spacing L T0, the water-cooled node is to this frame central spacing L after each frame Ti, De-scaling box sprays length L between the area of coverage C0, the water-cooled node sprays length L between the area of coverage after each frame Ci, otherwise execution in step 2d);
Step 2d) when con4=1, perhaps when con4=2, collect strip steel related process parameter, comprising: supplied materials width set-point B, supplied materials average thickness set-point h 0, each stand stretch drag set-point σ Si Set, supplied materials thickness cross direction profiles set-point h 0, each frame outlet average thickness set-point h i, first frame inlet given tension value T 0, each frame outlet given tension value T i, each frame coefficient of friction set-point μ i, each frame work roll bending power set-point S 1i, con2 i=1 frame intermediate calender rolls bending roller force set-point S 2i, each frame working roll roll shifting amount set-point RS Wi, con2 i=1 frame intermediate roll shifting amount set-point RS Mi, each frame roller amount set-point Tilt that inclines i, con2 i=2 breast roller angle of the crossing set-point θ i, otherwise execution in step 2e);
Step 2e) when con4=3, collect strip steel related process parameter, comprising: current winding steel steel reel number j, supplied materials steel grade s s j, the width measured value B j, supplied materials average thickness measured value h 0 j, supplied materials length measured value L s j, each stand stretch drag setting value σ Si Setj, supplied materials head thickness cross direction profiles measured value h 0 Hej, supplied materials temperature measured value T s j, supplied materials speed measured value v 0 j, the pure time measured value t that rolls r j, roll the back intermittent time measured value t p j, each frame outlet average thickness setting value h i j, first frame inlet tension force setting value T 0 j, each frame outlet tension force setting value T i j, each frame coefficient of friction setting value μ i j, each frame work roll bending power upper and lower limit S 1i Max, S 1i Min, con2 i=1 frame intermediate calender rolls bending roller force upper and lower limit S 1i Max, S 1i Min, each frame working roll roll shifting amount setting value upper and lower limit RS Wi Max, RS Wi Min, con2 i=1 frame intermediate roll shifting amount upper and lower limit RS Mi Max, RS Mi Min, con2 i=2 breast roller angle of the crossing upper and lower limit θ i Max, θ i Min, the last frame export goal flatness and lateral distribution σ of industrial requirements decision 1n Obj, the last frame export goal thickness cross direction profiles h of industrial requirements decision n Obj, otherwise execution in step 2f);
Step 2f) when con4=3 and con3=2, collect strip steel related process parameter, comprising: the state variable Bu that whether the water-cooled node sprays water after current each frame of winding steel i j, otherwise execution in step 2g);
Step 2g) when con4=3, order Tilt i j = 0 , ( i = 1 ~ n ) , T Wsi=T Wi(i=1~n), C Wi=C Wgi(i=1~n), C Bi=C Bgi(i=1~n), T Wi 0 = T Wi , ( i = 1 ~ n ) , W Wi 0 = 0 , ( i = 1 ~ n ) , W Bi 0 = 0 , ( i = 1 ~ n ) ; As con4=3 and con2 i, make C at=1 o'clock Mi=C Mgi(i=1~n), W Mi 0 = 0 , ( i = 1 ~ n ) , Otherwise execution in step 2h);
Step 2h) when con4=4, collecting belt steel supplied materials head thickness cross direction profiles measured value h 0 Hej, the current section thickness of band steel supplied materials cross direction profiles measured value h 0 Cuj, otherwise execution in step 2i);
Step 2i) when con4=5, collects the band steel end frame export goal flatness and lateral distribution σ of industrial requirements decision 1n Obi, the frame outlet of band steel end is as front section flatness and lateral distribution measured value σ 1n Cuj
Described in the said method in the step 4) when con4=1 and con1=3; Perhaps when con4=2 and con1=3; Perhaps when con4=3 and con1=3; Calculate the belt steel temperature field according to rolled piece and roller temperature field module, and may further comprise the steps according to each stand stretch district average deformation drag of rolled piece plastic deformation module calculating simultaneously:
Step 3a) order σ Si j = σ Si Setj , ( i = 1 ~ n ) ;
Step 3b) when con3=1, with T s j, v 0 j, B j, h 0 j~h n j, μ 1 j~μ n j, σ S1 j~σ Sn j, T Ws1~T Wsn, T C1~T Cn, T a, L M12~L M (n-1) nAs input parameter, can obtain T according to rolled piece and roller temperature field module jWhen con3=2, with T s j, v 0 j, B j, h 0 j~h n j, μ 1 j~μ n j, σ S1 j~σ Sn j, T Ws1~T Wsn, T C1~T Cn, T a, L M12~L M (n-1) n, L T0~L Tn, L C0~L Cn, Bu 1~Bu nAs input parameter, can obtain T according to rolled piece and roller temperature field module jThrough to T jIntegration can obtain each frame rolling deformation district rolled piece mean temperature T Swi j
Step 3c) with s s j, h I-1 j(i=1~n), h i j(i=1~n), v 0 j, h 0, T Swi jAs input parameter, can obtain the new average deformation drag value σ of each frame according to rolled piece plastic deformation module Si' j, calculate ϵ 1 i = | σ Si ′ j - σ Si j | / σ Si j , ( i = 1 ~ n ) ;
Step 3d) judges ε 11~ε 1nWhether all smaller or equal to given in advance little value ε 1If, ε 11~ε 1nAll smaller or equal to ε 1, obtain final σ S1 j~σ Sn j, and according to T at this moment jEach frame rolling deformation district rolled piece average surface temperature T of integral and calculating Ssi j, calculate and finish; If ε 11~ε 1nNot all smaller or equal to ε 1, order σ Si j = σ Si ′ j , ( i = 1 ~ n ) , Repeating step 3b), step 3c), step 3d), up to ε 11~ε 1nAll smaller or equal to ε 1Till.
Described in the said method in the step 5) when con4=1, perhaps when con4=2, rolled piece plastic deformation module and roll elastic deformation module are coupled, find the solution each frame outlet glacing flatness, thickness and draught pressure cross direction profiles and may further comprise the steps:
Step 4a) supposes i frame exit thickness cross direction profiles h i j = h ‾ i j h i - 1 j / h ‾ i - 1 j .
Step 4b) with σ Si j, μ i j, h I-1 j, h i j, B j, T I-1 j, T i jAs input parameter, try to achieve σ according to rolled piece plastic deformation module 1i j, p i j
Step 4c) works as con2 i=1 o'clock, with L Wi, L Mi, L Bi, L 1i, L 2i, L 3i, D Wi, D Mi, D Bi, D Bni, C Wi, C Mi, C Bi, S 1i j, S 2i j, RS Wi j, RS Mi j, Tilt i j, p i j, h i jAs input parameter, try to achieve new h according to the roll elastic deformation module i' jAnd the q corresponding with it Wmi j, q Mbi j, work as con2 i=2 o'clock, with L Wi, L Bi, L 1i, L 3i, D Wi, D Bi, D Bni, C Wi, C Bi, S 1i j, RS Wi j, Tilt i j, θ i, p i j, h i jAs input parameter, try to achieve new h according to the roll elastic deformation module i' jAnd the q corresponding with it Wbi j, calculate ϵ i = ∫ - B j / 2 B j / 2 | h i ′ j - h i j | / ( B h ‾ i j ) Dy .
Step 4d) judges ε iWhether smaller or equal to given in advance little value ε, if ε i≤ε obtains final σ 1i j, h i j, p i jAnd the q corresponding with it Wmi j, q Mbi jOr q Wbi jIf, ε i>ε, order h i j = ( 1 - α ) h i j + α h i ′ j (wherein α is a relaxation factor, 0<α≤0.1), repeating step 4b), step 4c), step 4d), up to ε iTill≤the ε.
Described in the said method in the step 6) when con4=3, calculate each frame target glacing flatness curve σ according to the plate-shape standard curve module 1i Obj(i=1~n) and target thickness cross direction profiles curve h i Obj(i=1~n) may further comprise the steps:
Step 5a) makes i=1;
Step 5b) sets i frame target glacing flatness curve and target thickness cross direction profiles curvilinear equation;
Each frame target glacing flatness curvilinear equation is:
σ 1 i obj ( y ) = a 0 + a 2 y 2 + a 4 y 4
Y in the formula---band steel horizontal relative coordinate, to another side, y from-1 to+1 from one side of band steel; a 0, a 2, a 4---fitting coefficient;
Each frame target thickness cross direction profiles curvilinear equation:
h i obj ( y ) = b 0 + b 2 y 2 + b 4 y 4
B in the formula 0, b 2, b 4---fitting coefficient.
Step 5c) if i is less than or equal to n-1, changes step 5d over to).Otherwise, change step 5f over to);
Step 5d) makes a 4=0, adopt the given unknowm coefficient a of certain optimized Algorithm 2, with B, h i, σ 1i ObjAs input parameter, try to achieve plate shape according to the good discrimination module of glacing flatness and differentiate factor ξ i, the condition of convergence is for making ξ iThe convergence setting value judges whether to satisfy the condition of convergence, as satisfies the condition of convergence as far as possible, obtains i passage target glacing flatness curve standard value, is designated as σ i *(y), as not satisfying the condition of convergence, repeated execution of steps 5d), till satisfying the condition of convergence;
Step 5e) adopts the given h of certain optimized Algorithm i ObjUnknowm coefficient b 2, b 4, with σ Si Set, μ i, h I-1 Obj, h i Obj, B, T I-1, T iAs input parameter, try to achieve σ according to rolled piece plastic deformation module 1i Obj(preceding tensile stress cross direction profiles σ 1iInhomogeneous composition), the condition of convergence is for making σ 1i Obj(preceding tensile stress cross direction profiles σ 1iInhomogeneous composition) convergence i passage target glacing flatness curve standard value σ as far as possible i *(y), judge whether to satisfy the condition of convergence, as satisfy the condition of convergence, obtain i passage target glacing flatness curve σ 1i ObjWith target thickness cross direction profiles curve h i Obj, as not satisfying the condition of convergence, repeated execution of steps 5e), till satisfying the condition of convergence;
Step 5f) according to back step operation two unknowm coefficient a is set in the specific (special) requirements of plate shape 2, a 4, obtain n passage target glacing flatness curve σ 1n Obj
Step 5g) adopts the given h of certain optimized Algorithm n ObjUnknowm coefficient b 2, b 4, with σ Sn Set, μ n, h N-1 Obj, h n Obj, B, T N-1, T nAs input parameter, try to achieve σ according to rolled piece plastic deformation module 1n Obj(preceding tensile stress cross direction profiles σ 1nInhomogeneous composition), the condition of convergence is for making σ 1n Obj(preceding tensile stress cross direction profiles σ 1nInhomogeneous composition) convergence n passage target glacing flatness curve σ as far as possible 1n Obj, judge whether to satisfy the condition of convergence, as satisfy the condition of convergence, obtain n passage target glacing flatness curve σ 1n ObjWith target thickness cross direction profiles curve h n Obj, as not satisfying the condition of convergence, repeated execution of steps 5g), till satisfying the condition of convergence;
Step 5h) make i=i+1, change step 5b over to), up to i=n.
Described in the said method in the step 7) when con4=3, each frame plate shape control device is carried out optimizing calculates, obtain each frame plate shape setup parameter and may further comprise the steps:
Step 6a) makes i=1;
Step 6b) when con4=3 and con3=1; Exporting actual glacing flatness curve with each frame, to equal target glacing flatness curve be target; When con4=3 and con3=2; Equaling target thickness cross direction profiles curve with each frame outlet actual (real) thickness cross direction profiles curve is target, adopts optimized Algorithm given one group of concrete parameter in the upper and lower scope of plate shape control parameter, specifically comprises: each frame work roll bending power S 1i j, each frame working roll roll shifting amount RS Wi j, con2 i=1 frame intermediate calender rolls bending roller force set-point S 2i j, con2 i=1 frame intermediate roll shifting amount RS Mi j, con2 i=2 breast roller angle of the crossing set-point θ i, change step 6c over to).
Step 6c) supposes i frame exit thickness cross direction profiles h i j = h ‾ i j h i - 1 j / h ‾ i - 1 j .
Step 6d) with σ Si j, μ i j, h I-1 j, h i j, B j, T I-1 j, T i jAs input parameter, try to achieve σ according to rolled piece plastic deformation module 1i j, p i j
Step 6e) works as con2 i=1 o'clock, with L Wi, L Mi, L Bi, L 1i, L 2i, L 3i, D Wi, D Mi, D Bi, D Bni, C Wi, C Mi, C Bi, S 1i j, S 2i j, RS Wi j, RS Mi j, Tilt i j, p i j, h i jAs input parameter, try to achieve new h according to the roll elastic deformation module i' jAnd the q corresponding with it Wmi j, q Mbi j, work as con2 i=2 o'clock, with L Wi, L Bi, L 1i, L 3i, D Wi, D Bi, D Bni, C Wi, C Bi, S 1i j, RS Wi j, Tilt i j, θ i, p i j, h i jAs input parameter, try to achieve new h according to the roll elastic deformation module i' jAnd the q corresponding with it Wbi j, calculate ϵ 2 i = ∫ - B j / 2 B j / 2 | h i ′ j - h i j | / ( B h ‾ i j ) Dy .
Step 6f) judges ε 2iWhether smaller or equal to given in advance little value ε 2If, ε 2i≤ε 2, obtain final h i j, p i jAnd the q corresponding with it Wmi j, q Mbi jOr q Wbi j, change step 6g over to), if ε 2i>ε 2, order h i j = ( 1 - α ) h i j + α h i ′ j (wherein α is a relaxation factor, 0<α≤0.1), repeating step 6d), step 6e), step 6f), up to ε 2i≤ε 2Till.
Step 6g) when con4=3, judges whether to satisfy the computation optimization condition of convergence.As satisfy the condition of convergence, write down the optimum plate shape of current winding steel i frame setup parameter, change step 6h over to); As do not satisfy the condition of convergence, repeating step 6b)~step 6g), till satisfying.
Step 6h) judge whether i=n,, change step 6b over to if i ≠ n makes i=i+1), up to i=n.
Described in the said method in the step 8) when con4=3, calculate rolling each the frame work roll thermal roll shape of back that finishes of current volume according to rolled piece and roller temperature field module, calculate rolling each the frame working roller abrasion roll shape of back that finishes of current volume according to roller system wearing and tearing module; May further comprise the steps:
Step 7a) makes i=1;
Step 7b) with T Wi J-1, T Ssi j, t r j, t p j, D WiAs input parameter, try to achieve T by rolled piece and roller temperature field module Wi j
Step 7c) with T Wi j, T Wi, D WiAs input parameter, try to achieve C by the roll elastic deformation module Wti j
Step 7d) judge whether i equals n,, change step 7b over to if i ≠ n makes i=i+1), if i=n changes step 7e over to);
Step 7e) makes i=1;
Step 7f) works as con2 i=1 o'clock, with W Wi J-1, W Mi J-1, W Bi J-1, L s j, h 0 j, h i j, p i j, q Wmi j, q Mbi jAs input parameter, try to achieve W by roller system wearing and tearing module Wi j, W Mi j, W Bi j, work as con2 i=2 o'clock, with W Wi J-1, W Bi J-1, L s j, h 0 j, h i j, p i j, q Wbi jAs input parameter, try to achieve W by roller system wearing and tearing module Wi j, W Bi j
Step 7g) judge whether i equals n,, change step 7f over to if i ≠ n makes i=i+1), up to i=n.
Described in the said method in the step 9) when con4=1; Perhaps when con4=4; Perhaps when con4=5, glacing flatness or thickness cross direction profiles curve are carried out pattern-recognition, calculate each frame plate shape control device regulated quantity matrix K according to plate shape control module according to plate shape pattern recognition module i j(i=1~n) may further comprise the steps:
Step 8a) when con4=1, make i=1, otherwise execution in step 8b);
Step 8b) when con4=1, with σ 1i jAs input parameter, try to achieve a by plate shape pattern recognition module 1i j, a 2i j, a 3i j, a 4i j, with h i jAs input parameter, try to achieve e by plate shape pattern recognition module 1i j, e 2i j, e 3i j, e 4i j, otherwise execution in step 8c);
Step 8c) judge when con4=1 whether i equals n,, change step 8b over to if i ≠ n makes i=i+1), up to i=n, otherwise execution in step 8d);
Step 8d) when con4=4, respectively with h 0 Cuj, h 0 HejAs input parameter, try to achieve e by formula plate shape pattern recognition module 10 Cuj, e 20 Cuj, e 30 Cuj, e 40 CujWith e 10 Hej, e 20 Hej, e 30 Hej, e 40 Hej, and calculate the deviation delta e of corresponding sub-eigenvalue 10 j, Δ e 20 j, Δ e 30 j, Δ e 40 j, otherwise execution in step 8e);
Step 8e) when con4=4 and con3=1, with Δ e 10 j, Δ e 20 j, Δ e 30 j, Δ e 40 j, as input parameter, try to achieve each frame plate shape control device regulated quantity matrix K by plate shape control module i j(i=1~n); When con4=4 and con3=2, with Δ e 10 j, Δ e 20 j, Δ e 30 j, Δ e 40 j, as input parameter, try to achieve each frame plate shape control device regulated quantity matrix K by plate shape control module i j(i=1~n), otherwise execution in step 8f);
Step 8f) when con4=5, respectively with σ 1n Cuj, σ 1n ObjAs input parameter, try to achieve a by plate shape pattern recognition module 1n Cuj, a 2n Cuj, a 3n Cuj, a 4n CujWith a 1n Obj, a 2n Obj, a 3n Obj, a 4n Obj, and calculate the deviation delta a of corresponding sub-eigenvalue 1n j, Δ a 2n j, Δ a 3n j, Δ a 4n j, otherwise execution in step 8g);
Step 8g) when con4=5, with Δ a 1n j, Δ a 2n j, Δ a 3n j, Δ a 4n jAs input parameter, try to achieve last frame plate shape control device regulated quantity matrix K by plate shape control module n j
Described in the said method in the step 10) when con4=3, rolling back each frame hot-rolling type, roll wear roll shape and the original grinding roll shape of roll of finishing of current volume superposeed according to following mode:
Work as con2 i=1 o'clock, order C Wi = C Wgi + C Wti j - W Wi j , C Mi = C Mgi - W Mi j , C Bi = C Bgi - W Bi j ;
Work as con2 i=2 o'clock, order C Wi = C Wgi + C Wti j - W Wi j , C Bi = C Bgi - W Bi j
The invention has the beneficial effects as follows: plat control system is divided into 8 modules that are mutually related; Each module is accomplished an independently function; And for the calculating of correlation module provides input parameter, according to the internal relation of 8 modules carry out integrated after, can form the plurality of plate-shape control technology; Comprise: analysis of plate shape control performance and type selection technology, roll shape optimisation technique, plate shape are set control technology, plate shape FEEDFORWARD CONTROL technology, plate shape feedback control technology; Only need adopt a cover calculation process during Practical Calculation, just can be comprehensively, optimum, solve plurality of plate-shape control problem easily, saved working time effectively.
Description of drawings
Fig. 1 is each module correlation figure of plat control system;
The global coordinate system x-z floor map that Fig. 2 adopts for the present invention;
The global coordinate system x-y floor map that Fig. 3 adopts for the present invention;
Fig. 4 is the local coordinate system sketch map that working roll temperature field module adopts.
Among Fig. 1,1 is rolled piece plastic deformation module, and 2 is the roll elastic deformation module; 3 is rolled piece and roller temperature field module, and 4 are roller system wearing and tearing module, and 5 is the good discrimination module of glacing flatness; 6 is plate shape pattern recognition module, and 7 is the plate-shape standard curve module, and 8 is plate shape control module.
The specific embodiment
First embodiment
The explanation of this example how to adopt the present invention to carry out the control performance analysis of plate shape and type is selected.At first given concrete type clearly need carry out the concrete plate shape control device that plate shape control performance is analyzed.With the complete four roller CVC continuous hot-rolling mills of six frames is example, the plate shape control performance of analytical work roller bending roller force.
The 1st step, execution in step 1), set up the tandem mill coordinate system;
The 2nd step, execution in step 2), the back con1=2 that is finished, con2 i=2 (i=1~6), con3=2, con4=1;
The 3rd step, execution in step 3), collect relevant input parameter;
In the 4th step, because con1=2, step 4) is skipped automatically, execution in step 5), rolled piece plastic deformation module 1 and roll elastic deformation module 2 are coupled, obtain each frame outlet glacing flatness, thickness cross direction profiles σ 1i(i=1~6), h i(i=1~6);
In the 5th step, because con4=1, step 6), step 7), step 8) are skipped automatically, execution in step 9), obtain the characteristic value a that each frame exports glacing flatness, thickness cross direction profiles according to plate shape pattern recognition module 6 1i j, a 2i j, a 3i j, a 4i jWith e 1i j, e 2i j, e 3i j, e 4i j
The 6th step; Repeatedly given as required different operating roller bending roller force; Repeat~the 5 step of the 4th step, can obtain the complete four roller CVC continuous hot-rolling mill work roll bending power control performances of this six frame, analyze the control performance of other plate shape control device like need; Like working roll roll shifting amount, the roller amount of inclining, computational process similarly.
When being used for the type selection; Can import the device parameter of various alternative types; Successively the various plate shape of various types control device is carried out the analysis of plate shape control performance, finally can obtain the overall plate shape control ability of various types, thereby select to provide foundation for type.
Second embodiment
With the full normal four-roller continuous hot-rolling mill of six frames (no work roll bending function) is example, explains that how adopting the present invention that some old milling train that lacks plate shape control device is carried out with plate shape optimum is the backing roll roll shape optimal design of target.Grinding roller efficient is improved for simplifying the grinding roller process management in general production scene, hopes that each frame support roller adopts identical roll shape.Therefore, the unified following formula that adopts of the original grinding roll shape of each frame is represented: C Bgi(y)=c 2(2y/L Bi) 2+ c 4(2y/L Bi) 4+ c 6(2y/L Bi) 6(i=1~6), wherein, C BgiThe original grinding roll shape of each frame support roller, L BiFor each frame support roller body of roll long, c 2, c 4, c 6Backing roll roll shape characteristic parameter.The target of roll shape optimization is decided to be seeks optimum c 2, c 4, c 6, make F 2 = ∫ - B / 2 B / 2 | σ 16 Ob - σ 16 | / B Minimum.Below provide concrete computational process:
In the 1st step, adopt the given one group of backing roll roll shape characteristic parameter c of certain optimized Algorithm 2, c 4, c 6
The 2nd step, execution in step 1), set up the tandem mill coordinate system;
The 3rd step, execution in step 2), the back con1=2 that is finished, con2 i=2 (i=1~6), con3=2, con4=2;
The 4th step, execution in step 3), collect relevant input parameter;
In the 5th step, because con1=2, step 4) is skipped automatically, execution in step 5), rolled piece plastic deformation module (1) and roll elastic deformation module 2 are coupled, obtain each frame outlet glacing flatness, thickness cross direction profiles σ 1i(i=1~6), h i(i=1~6);
In the 6th step, judge whether to satisfy the computation optimization condition of convergence.As satisfy the condition of convergence, and calculate and finish, as not satisfying the condition of convergence, returned for the 1st step to repeat, till satisfying.
The 3rd embodiment
With the full UCM cold continuous rolling of five frames is example, and how the explanation of this example adopts the present invention to carry out the control of plate shape on-line setup.The plate shape that specifically needs to set is controlled the working roll roller bending roller force S that parameter comprises 1~5 frame 11 j~S 15 j, 1~5 frame intermediate calender rolls bending roller force S 21 j~S 25 j, 1~5 frame intermediate roll shifting amount RS M1 j~RS M5 jBelow provide concrete computational process:
The 1st step, execution in step 1), set up the tandem mill coordinate system;
The 2nd step, execution in step 2), the back con1=3 that is finished, con2 i=2 (i=1~6), con3=1, con4=3;
The 3rd step, execution in step 3), collect relevant input parameter;
The 4th step, execution in step 4), calculate the belt steel temperature field according to rolled piece and roller temperature field module 3, and calculate each stand stretch district average deformation drag according to rolled piece plastic deformation module 1 simultaneously;
In the 5th step, because con4=3, step 5) is skipped automatically, execution in step 6), calculate each frame target glacing flatness curve σ according to plate-shape standard curve module 7 1i Obj(i=1~n) and target thickness cross direction profiles curve h i Obj(i=1~n);
The 6th step, execution in step 7), each frame plate shape control device is carried out optimizing calculate, obtain the optimum plate shape of each frame setup parameter S 11 Opj~S 15 Opj, S 21 Opj~S 25 Opj, RS M1 Opj~RS M5 Opj
The 7th step, execution in step 8), calculate rolling each the frame work roll thermal roll shape of back that finishes of current volume according to rolled piece and roller temperature field module 3, calculate rolling each the frame working roller abrasion roll shape of back that finishes of current volume according to roller system wearing and tearing module 4;
In the 8th step, because con4=3, step 9) is skipped automatically, execution in step 10), rolling each frame hot-rolling type of back, roll wear roll shape and the original grinding roll shape of roll of finishing of current volume superposeed, for the back volume calculates the condition that provides;
The 9th step is if rolling down coiled strip changed for the 1st step over to repeat, if not rolling coiled strip down calculates and finishes.
The 4th embodiment
With the full HC six roller cold continuous rollings of five frames is example, explains and how to adopt the present invention to carrying out plate shape FEEDFORWARD CONTROL.Comprise each frame work roll bending power, intermediate roll shifting amount, the roller amount of inclining for this milling train plate shape control parameter.
The 1st step, execution in step 1), set up the cold continuous rolling coordinate system;
The 2nd step, execution in step 2), the back con1=3 that is finished, con2 i=1 (i=1~5), con3=1, con4=4;
The 3rd step, execution in step 3), collect relevant input parameter;
In the 4th step, because con4=4, step 4), step 5), step 6), step 7), step 8) are skipped automatically, execution in step 9), obtain each frame plate shape control device regulated quantity matrix K according to plate shape control module 8 i j(i=1~5), i.e. the regulated quantity of each frame work roll bending power, intermediate roll shifting amount, the roller amount of inclining, calculating finishes.
The 5th embodiment
With last frame is that normal four-roller milling train (the work roll bending function is arranged) hot continuous rolling is an example, explains and how to adopt the present invention to carrying out plate shape FEEDBACK CONTROL.For the normal four-roller milling train that the work roll bending function is arranged, its plate shape control device comprises work roll bending power and inclines the roller amount.
The 1st step, execution in step 1), set up the cold continuous rolling coordinate system;
The 2nd step, execution in step 2), the back con1=3 that is finished, con2 i=1 (i=1~5), con3=1, con4=5;
The 3rd step, execution in step 3), collect relevant input parameter;
In the 4th step, because con4=5, step 4), step 5), step 6), step 7), step 8) are skipped automatically, execution in step 9), obtain last frame plate shape control device regulated quantity matrix K according to plate shape control module 8 n j, the regulated quantity of the promptly last frame work roll bending power and the roller amount of inclining, calculating finishes.

Claims (20)

1. plate-shape control integrated system, it is characterized in that: this system comprises 8 modules that are mutually related: rolled piece plastic deformation module (1), roll elastic deformation module (2); Rolled piece and roller temperature field module (3); Roller system wearing and tearing module (4), the good discrimination module of glacing flatness (5), plate shape pattern recognition module (6); Plate-shape standard curve module (7), plate shape control module (8); Wherein,
Said rolled piece plastic deformation module (1) provides the draught pressure cross direction profiles for roll elastic deformation module (2) and roll wear module (4), distributes for good discrimination module of glacing flatness (5) and plate shape pattern recognition module (6) provide flatness and lateral;
Said roll elastic deformation module (2) provides band steel exports thickness cross direction profiles for rolled piece plastic deformation module (1) and plate shape pattern recognition module (6), for roll wear module (4) provides the roll gap pressure cross direction profiles;
Said rolled piece and roller temperature field module (3) provide the deformed area mean temperature for the calculating of resistance of deformation in the rolled piece plastic deformation module (1), for the calculating of work roll thermal roll shape in the roll elastic deformation module (2) provides the working roll temperature field;
Said roll wear module (4) provides the roll wear roll shape for the calculating of comprehensive roll shape in the roll elastic deformation module (2);
Good discrimination module of said glacing flatness (5) and the local integrated formation plate-shape standard curve module of rolled piece plastic deformation module (1) (7);
Said plate shape control module (8) is calculated plate shape control device regulated quantity according to the output result of plate shape pattern recognition module (6).
2. plate-shape control integrated system according to claim 1; It is characterized in that: said rolled piece plastic deformation module (1) is according to rolling condition; Consider the lateral flow of rolled piece; Set up deformed area Three-dimensional Flow velocity field, strain field, stress field, confirm the cross direction profiles of draught pressure and outlet glacing flatness, adopt Mathematical Modeling to represent as follows:
Figure FSB00000344127600011
Figure FSB00000344127600012
Figure FSB00000344127600013
Figure FSB00000344127600014
Wherein, i is the frame number,
Figure FSB00000344127600015
Be a row vector, comprised 3 element: σ 1i, σ 0i, p i, f 11Be streamline bar unit method computing formula, f 12Be stream ready-to-eat type noodle method computing formula, f 13Be bar unit calculus of variations computing formula, σ 1i, σ 0iBe respectively the forward and backward tensile stress cross direction profiles of i frame value, p iBe i frame draught pressure cross direction profiles value, σ SiBe i frame band steel average deformation drag, μ iBe i frame coefficient of friction, h I-1, h iBe respectively that i frame band steel is gone into, exit thickness cross direction profiles value, B is band steel supplied materials width, T i, T I-1Be respectively the forward and backward tension force of i frame;
Each frame band steel average deformation drag σ SiBy steel grade, drafts, mill speed and the decision of deformed area band mean temperature, adopt Mathematical Modeling to represent as follows:
Figure FSB00000344127600021
Wherein, f 14Be resistance of deformation computing formula, s sBe steel grade, v 0Be the first frame porch window of web velocities,
Figure FSB00000344127600022
Be deformed area band ensemble average temperature.
3. plate-shape control integrated system according to claim 1; It is characterized in that; Said roll elastic deformation module (2) is calculated amount of deflection, the flattening amount of roll according to Elasticity; Final definite rolled piece exit thickness cross direction profiles value and roll gap pressure cross direction profiles value, adopt Mathematical Modeling to represent as follows:
Figure FSB00000344127600023
Figure FSB00000344127600024
Figure FSB00000344127600025
Figure FSB00000344127600026
C wi=C wgi+C wti-W wi
C mi=C mgi-W mi
C bi=C bgi-W bi
Wherein, i is the frame number,
Figure FSB00000344127600027
Be a row vector, comprised 3 element: h i, q Wmi, q Mbi,
Figure FSB00000344127600028
Be a row vector, comprised 2 element: h i, q Wbi, f 21Be six-high cluster mill roll elastic deformation computing formula, f 22Be four-high mill roll elastic deformation computing formula, h iBe i frame band steel exports thickness cross direction profiles value, p iBe i frame draught pressure cross direction profiles value, q WmiBe i frame working roll and intermediate calender rolls roll gap pressure cross direction profiles, q MbiBe i frame intermediate calender rolls and backing roll roll gap pressure cross direction profiles, q WbiBe i frame working roll and backing roll roll gap pressure cross direction profiles, L WiBe i frame working roll barrel length, L MiBe i frame intermediate calender rolls barrel length, L BiBe i frame support roller barrel length, L 1iBe i frame work roll bending cylinder spacing, L 2iBe i frame intermediate calender rolls roller cylinder spacing, L 3iBe i frame depress oil cylinder centre-to-centre spacing, D WiBe i frame work roll diameter, D MiBe i frame intermediate calender rolls diameter, D BiBe i frame support roller diameter, D BniBe i frame support roller roll neck diameter, C WiBe the comprehensive roll shape of i frame working roll, C WgiBe the original grinding roll shape of i frame working roll, C W1iBe i frame work roll thermal roll shape, W WiBe i frame working roller abrasion roll shape, C MiBe the comprehensive roll shape of i frame intermediate calender rolls, C MgiBe the original grinding roll shape of i frame intermediate calender rolls, W MiBe i frame intermediate calender rolls wearing and tearing roll shapes, C BiBe the comprehensive roll shape of i frame support roller, C BgiBe the original grinding roll shape of i frame support roller, W BiBe i frame support roller wearing and tearing roll shapes, S 1iBe i frame work roll bending power, S 2iBe i frame intermediate calender rolls bending roller force, RS WiBe i frame working roll roll shifting amount, RS MiBe i frame intermediate roll shifting amount, Tilt iBe the i frame roller amount of inclining, θ iBe the i breast roller angle of the crossing,
Figure FSB00000344127600031
It is i frame band steel exports average thickness;
Formula C Wi=C Wgi+ C Wti-W WiMiddle roll hot-rolling type can calculate according to iron Mo Xinke formula, adopts Mathematical Modeling to represent as follows:
Figure FSB00000344127600032
Wherein, f 23Be iron Mo Xinke formula,
Figure FSB00000344127600033
I frame work roll thermal roll shape during for rolling j+1 winding steel, T WiBe temperature field during machine on the i frame working roll,
Figure FSB00000344127600034
I frame working roll temperature field during for rolling j+1 winding steel.
4. plate-shape control integrated system according to claim 1; It is characterized in that: the rolled piece temperature field that the rolled piece temperature field in said rolled piece and the roller temperature field module (3) adopts the two-dimensional finite calculus of finite differences to calculate rolling district and non-rolling district, confirm each frame rolling deformation district band ensemble average temperature through integration method
Figure FSB00000344127600035
, average surface temperature
Figure FSB00000344127600036
, outlet band transverse temperature distribution curve T Δ i, adopt Mathematical Modeling to represent as follows:
Figure FSB00000344127600037
(x>0)
Figure FSB00000344127600038
(x>0)
Figure FSB00000344127600039
Figure FSB000003441276000310
Wherein, T is a two dimensional field amount, and expression belt steel temperature field is the function of y, z, and coordinate system is defined as: x is that belt steel rolling direction, y are that strip width direction, z are short transverse, and the origin of coordinates is finish rolling inlet temperature measurer place, f 31Be divided into cold continuous rolling rolled piece temperature field computing formula, f 32Be continuous hot-rolling mill rolled piece temperature field computing formula, n is the total frame number of tandem mill, Be a row vector, comprised n-1 element: L M12~L M (n-1) n,
Figure FSB000003441276000312
Be a row vector, comprised 4n+1 element: L M12~L M (n-1) n, L T0, L T1~L Tn, L C0, L C1~L Cn, Bu 1~Bu n, T s(y z) is the two-dimensional temperature field of x=0 place rolled piece, v 0Be supplied materials speed, B is the supplied materials width,
Figure FSB00000344127600041
Be band steel supplied materials average thickness and each frame outlet average thickness, μ 1~μ nBe each frame coefficient of friction, σ S1~σ SnBe each frame deformation of rolled wire drag,
Figure FSB00000344127600042
Be each breast roller average surface temperature, T C1~T CnBe each frame cooling water or emulsion temperature, T aBe air themperature, L M12~L M (n-1) nBe frame spacing, L T0Be De-scaling box center to the first frame central spacing, L T1~L TnThe water-cooled node is to this frame central spacing, L after each frame C0For spraying, De-scaling box covers siding-to-siding block length, L C1~L CnThe water-cooled node sprays and covers siding-to-siding block length, Bu after each frame 1~Bu nFor describing the state variable that whether the water-cooled node sprays water after each frame, value 0 or 1 is got 0 expression and is not existed or do not spray water, and gets the normal water spray of 1 expression;
Calculate through finite difference calculus roller temperature field in said rolled piece and the roller temperature field module (3), adopts Mathematical Modeling to represent as follows:
Figure FSB00000344127600043
Wherein, i is the frame number, and j is rolling coil of strip number,
Figure FSB00000344127600044
Be after j winding steel rolling finishes
Figure FSB00000344127600045
Moment i frame working roll temperature field,
Figure FSB00000344127600046
Being respectively the pure of rolling j winding steel rolls the time and rolls back intermittent time, T WiTemperature during machine on the i frame working roll,
Figure FSB00000344127600047
Each frame rolling deformation district strip surface mean temperature during for rolling j winding steel.
5. plate-shape control integrated system according to claim 1 is characterized in that: said roller system's wearing and tearing modules (4) adopt Mathematical Modeling to represent as follows:
Figure FSB00000344127600048
Figure FSB00000344127600049
Figure FSB000003441276000411
Figure FSB000003441276000412
Wherein, i is the frame number, and j is rolling coil of strip number, f 51wBe six-high cluster mill working roll roll wear computing formula, f 51mBe six-high cluster mill intermediate calender rolls roll wear computing formula, f 51bBe six-high cluster mill backing roll roll wear computing formula, f 52wBe four-high mill working roller abrasion computing formula, f 52bBe four-high mill backing roll roll wear computing formula,
Figure FSB00000344127600051
Be j winding steel inlet thickness cross direction profiles value rolling behind the machine on the roll,
Figure FSB00000344127600052
Be rolling j winding steel exports thickness cross direction profiles value behind the machine on the i breast roller,
Figure FSB00000344127600053
Be rolling j winding steel rolling pressure cross direction profiles value behind the machine on the i breast roller,
Figure FSB00000344127600054
Be rolling j winding steel working roll and intermediate calender rolls roll gap pressure cross direction profiles value behind the machine on the i breast roller,
Figure FSB00000344127600055
Be rolling j winding steel intermediate calender rolls and backing roll roll gap pressure cross direction profiles value behind the machine on the i breast roller,
Figure FSB00000344127600056
Be rolling j winding steel working roll and backing roll roll gap pressure cross direction profiles value behind the machine on the i breast roller,
Figure FSB00000344127600057
Be the intact back of j-1 winding steel rolling i frame working roller abrasion roll shape,
Figure FSB00000344127600058
Be the intact back of j-1 winding steel rolling i frame intermediate calender rolls wearing and tearing roll shapes,
Figure FSB00000344127600059
Be the intact back of j-1 winding steel rolling i frame support roller wearing and tearing roll shapes,
Figure FSB000003441276000510
Be machine time wearing and tearing roll shape on the i frame working roll, Be machine time wearing and tearing roll shape on the i frame intermediate calender rolls,
Figure FSB000003441276000512
Be machine time wearing and tearing roll shape on the i frame support roller,
Figure FSB000003441276000513
Be the intact back of j winding steel rolling i frame working roller abrasion roll shape,
Figure FSB000003441276000514
Be the intact back of j winding steel rolling i frame intermediate calender rolls wearing and tearing roll shapes,
Figure FSB000003441276000515
Be the intact back of j winding steel rolling i frame support roller wearing and tearing roll shapes,
Figure FSB000003441276000516
Be j winding steel supplied materials length.
6. plate-shape control integrated system according to claim 1; It is characterized in that: the good discrimination module of said glacing flatness (5) distributes according to the band steel exports flatness and lateral that rolled piece plastic deformation module (1) calculates; Differentiate the factor through plate shape and judge plate shape quality, adopt Mathematical Modeling to represent as follows:
Figure FSB000003441276000517
Wherein, i is the frame number, f 6Be the differentiation factorization method of cutting apart based on bar unit, B is a strip width, I frame exit band steel average thickness, σ 1iBe tensile stress cross direction profiles value before the i frame, ξ iI frame band steel exports glacing flatness is differentiated the factor, ξ iThe band steel was straight in>1 o'clock, ξ iThe band steel was in critical instability status, ξ in=1 o'clock i<1 o'clock band steel unstability.
7. plate-shape control integrated system according to claim 1; It is characterized in that: said plate shape pattern recognition module (6) adopts the least square method based on Legnedre polynomial; Pattern classification is carried out in plate shape distribution to rolling back band steel, adopts Mathematical Modeling to represent as follows:
Figure FSB000003441276000519
Figure FSB000003441276000520
Figure FSB000003441276000522
Wherein, i is the frame number, σ 1iBe tensile stress cross direction profiles value before the i frame, h iBe i frame band steel exports thickness cross direction profiles value,
Figure DEST_PATH_FSB00000613083600011
Be a row vector, comprised 4 element: a 1i, a 2i, a 3i, a 4i,
Figure DEST_PATH_FSB00000613083600012
Be a row vector, comprised 4 element: e 1i, e 2i, e 3i, e 4i, f 7Be the least square method based on Legnedre polynomial, a 1i~a 4iBe respectively the i frame and export strip profile and flatness cross direction profiles 1 time, 2 times, 3 times, 4 sub-eigenvalues, e 1i~e 4iBe respectively the i frame and export belt steel thickness cross direction profiles 1 time, 2 times, 3 times, 4 sub-eigenvalues.
8. plate-shape control integrated system according to claim 1 is characterized in that: said plate-shape standard curve module (7) adopts Mathematical Modeling to represent as follows:
Figure RE-FSB00000613083600014
Figure RE-FSB00000613083600015
Figure RE-FSB00000613083600016
Wherein, f 8Be tandem mill plate-shape standard curve computing formula, n is the total frame number of tandem mill, and B is a strip width, h 0Be finishing mill incoming band steel thickness cross direction profiles,
Figure RE-FSB00000613083600017
Be finishing mill incoming band steel average thickness,
Figure RE-FSB00000613083600018
Be the 1st~n frame band steel exports average thickness, Be a row vector, comprised n element:
Figure RE-FSB000006130836000111
Be a row vector, comprised n element:
Figure RE-FSB000006130836000112
Figure RE-FSB000006130836000113
Be the tandem mill end frame outlet band target glacing flatness curve of industrial requirements decision,
Figure RE-FSB000006130836000114
The tandem mill end frame outlet band target thickness cross direction profiles curve of industrial requirements decision,
Figure RE-FSB000006130836000115
Be respectively each frame export goal glacing flatness curve,
Figure RE-FSB000006130836000116
Be respectively each frame export goal thickness cross direction profiles curve.
9. plate-shape control integrated system according to claim 1 is characterized in that: said plate shape control module (8) adopts Mathematical Modeling to represent as follows:
Figure RE-FSB000006130836000117
Figure RE-FSB000006130836000119
Wherein, f 91Be the plate shape control device compensation rate computing formula that last frame outlet actual measurement glacing flatness and target glacing flatness deviation cause, f 92For guaranteeing that each frame exports actual glacing flatness curve and equals under the prerequisite of target glacing flatness curve; Front section actual measurement thickness cross direction profiles worked as by supplied materials and the supplied materials head is surveyed each frame plate shape control device compensation rate computing formula that thickness cross direction profiles deviation causes; Be applicable to cold continuous rolling, f 93For guaranteeing that the actual exit thickness cross direction profiles of each frame equals under the prerequisite of target thickness cross direction profiles; Front section actual measurement thickness cross direction profiles worked as by supplied materials and the supplied materials head is surveyed each frame plate shape control device compensation rate computing formula that thickness cross direction profiles deviation causes; Be applicable to continuous hot-rolling mill, Δ a 1n, Δ a 2n, Δ a 3n, Δ a 4nBe respectively last frame outlet poor when front section actual measurement glacing flatness curve 1 time, 2 times, 3 times, 4 sub-eigenvalues and target glacing flatness curve 1 time, 2 times, 3 times, 4 sub-eigenvalues, Δ e 10, Δ e 20, Δ e 30, Δ e 40Be respectively supplied materials poor when front section actual measurement thickness cross direction profiles 1 time, 2 times, 3 times, 4 sub-eigenvalues and band steel toe portion actual measurement thickness cross direction profiles 1 time, 2 times, 3 times, 4 sub-eigenvalues,
Figure RE-FSB00000613083600021
Be all plate shape control device regulated quantity matrixes of i frame, n is the total frame number of tandem mill, Be all plate shape control device regulated quantity matrixes of n frame.
10. the manner of execution of a plate-shape control integrated system, it is characterized in that: this method comprises following computer step:
Step 1) is set up the tandem mill coordinate system;
Step 2) confirm the use object of plate-shape control integrated system, the definition integer variable:
{con1,con2 1,…,con2 n,con3,con4,con5},
Wherein, n is a tandem mill frame number, through this n+4 of initialization variable, confirms the use object of plate-shape control integrated system;
Step 3) is collected relevant rolling equipment and technological parameter;
Step 4) is when con4=1 and con1=3; Perhaps when con4=2 and con1=3; Perhaps when con4=3 and con1=3; Calculate the belt steel temperature field according to rolled piece and roller temperature field module (3), and calculate each stand stretch district average deformation drag according to rolled piece plastic deformation module (1) simultaneously;
Step 5) perhaps when con4=2, is coupled to rolled piece plastic deformation module (1) and roll elastic deformation module (2) when con4=1, finds the solution each frame outlet glacing flatness, thickness and draught pressure cross direction profiles;
Step 6) is calculated each frame target glacing flatness curve
Figure RE-FSB00000613083600023
and target thickness cross direction profiles curve according to plate-shape standard curve module (7) when con4=3
Step 7) is carried out optimizing to each frame plate shape control device and is calculated when con4=3, obtains each frame plate shape setup parameter;
Step 8) is calculated rolling each the frame work roll thermal roll shape of back that finishes of current volume according to rolled piece and roller temperature field module (3) when con4=3, calculate rolling each the frame working roller abrasion roll shape of back that finishes of current volume according to roller system's wearing and tearing modules (4);
Step 9) is when con4=1; Perhaps when con4=4; Perhaps when con4=5; According to plate shape pattern recognition module (6) glacing flatness or thickness cross direction profiles curve are carried out pattern-recognition, calculate each frame plate shape control device regulated quantity matrix according to plate shape control module (8)
Step 10) superposes rolling each frame hot-rolling type of back, roll wear roll shape and the original grinding roll shape of roll of finishing of current volume when con4=3.
11. the manner of execution of plate-shape control integrated system according to claim 10 is characterized in that: it is specific as follows to set up the tandem mill coordinate system in the said step 1):
The length of rolled piece direction is elected to be the x axle, and the rolled piece width is elected to be the y axle, and the rolled piece short transverse is elected to be the z axle, and the origin of coordinates is chosen milling train width center at the y axle before the x axle is chosen tandem mill first rolling mill inlet, choose the rolled piece mid-depth at the z axle.
12. the manner of execution of plate-shape control integrated system according to claim 10 is characterized in that: the use object said step 2) comprises:
1a) be applied to hot continuous rolling when line computation, make con1=3, con3=2,
Be applied to cold continuous rolling when line computation, make con1=3, con3=1,
When being applied to the hot continuous rolling off-line analysis, make con1=2, con3=2,
When being applied to the cold continuous rolling off-line analysis, make con1=1, con3=1;
When 1b) being applied to six-high cluster mill, make con2 i=1, i=1~n, n are the frame numbers,
When being applied to four-high mill, make con2 i=2, i=1~n, n are the frame numbers;
When 1c) being applied to control performance analysis of plate shape and type selection technology, make con4=1,
When being applied to the roll shape optimisation technique, make con4=2,
When being applied to plate shape setting control technology, make con4=3,
When being applied to plate shape FEEDFORWARD CONTROL technology, make con4=4,
When being applied to plate shape feedback control technology, make con4=5.
13. the manner of execution of plate-shape control integrated system according to claim 10 is characterized in that: relevant rolling equipment of the collection in the said step 3) and technological parameter may further comprise the steps:
Step 2a) when con4=1, perhaps when con4=2, perhaps when con4=3, collect the rolling equipment relevant parameter, comprising: each frame working roll barrel length L Wi, each frame work roll diameter D Wi, con2 i=1 o'clock frame intermediate calender rolls barrel length L Mi, con2 i=1 o'clock frame intermediate calender rolls diameter D Mi, each frame support roller barrel length L Bi, each frame support roller diameter D Bi, each frame support roller roll neck diameter D Bni, each frame work roll bending cylinder spacing L 1i, con2 i=1 o'clock frame intermediate calender rolls roller cylinder spacing L 2i, each frame depress oil cylinder centre-to-centre spacing L 3i, the original grinding roll shape of each frame working roll C Wgi, con2 i=1 o'clock the original grinding roll shape of frame intermediate calender rolls C Mgi, the original grinding roll shape of each frame support roller C Bgi
Step 2b) when con4=1 and con1=3, perhaps when con4=2 and con1=3, perhaps when con4=3 and con1=3, collect the rolling equipment relevant parameter, comprising: each frame working roll is initially gone up the machine temperature T Wi, each frame cooling water or emulsion temperature T Ci, air themperature T a, frame spacing L Mi (i+1)
Step 2c) when con4=1 and con1=3 and con3=2; Perhaps when con4=2 and con1=3 and con3=2; Perhaps when con4=3 and con1=3 and con3=2, collect the rolling equipment relevant parameter, comprising: De-scaling box center to the first frame central spacing L T0, the water-cooled node is to this frame central spacing L after each frame Ti, De-scaling box sprays length L between the area of coverage C0, the water-cooled node sprays length L between the area of coverage after each frame Ci
Step 2d) when con4=1, perhaps when con4=2, collect strip steel related process parameter, comprising: supplied materials width set-point B, supplied materials average thickness set-point
Figure RE-FSB00000613083600041
Each stand stretch drag set-point Supplied materials thickness cross direction profiles set-point h 0, each frame outlet average thickness set-point
Figure RE-FSB00000613083600043
First frame inlet given tension value T 0, each frame outlet given tension value T i, each frame coefficient of friction set-point μ i, each frame work roll bending power set-point S 1i, con2 i=1 o'clock frame intermediate calender rolls bending roller force set-point S 2i, each frame working roll roll shifting amount set-point RS Wi, con2 i=1 o'clock frame intermediate roll shifting amount set-point RS Mi, each frame roller amount set-point Tilt that inclines i, con2 i=2 o'clock breast roller angle of the crossing set-point θ i
Step 2e) when con4=3, collect strip steel related process parameter, comprising: current winding steel steel reel number j, supplied materials steel grade The width measured value B j, supplied materials average thickness measured value Supplied materials length measured value
Figure RE-FSB00000613083600046
Each stand stretch drag setting value
Figure RE-FSB00000613083600047
Supplied materials head thickness cross direction profiles measured value
Figure RE-FSB00000613083600048
Supplied materials temperature measured value
Figure RE-FSB00000613083600049
Supplied materials speed measured value
Figure RE-FSB000006130836000410
The pure time measured value that rolls
Figure RE-FSB000006130836000411
Roll back intermittent time measured value Each frame outlet average thickness setting value First frame inlet tension force setting value
Figure RE-FSB000006130836000414
Each frame outlet tension force setting value Each frame coefficient of friction setting value
Figure RE-FSB000006130836000416
Each frame work roll bending power upper and lower limit
Figure RE-FSB000006130836000417
Con2 i=1 o'clock frame intermediate calender rolls bending roller force upper and lower limit
Figure RE-FSB000006130836000418
Each frame working roll roll shifting amount setting value upper and lower limit
Figure RE-FSB000006130836000419
Con2 i=1 frame intermediate roll shifting amount upper and lower limit
Figure RE-FSB000006130836000420
Con2 i=2 breast roller angle of the crossing upper and lower limits
Figure RE-FSB000006130836000421
The last frame export goal flatness and lateral of industrial requirements decision distributes
Figure RE-FSB000006130836000422
The last frame export goal thickness cross direction profiles of industrial requirements decision
Figure RE-FSB000006130836000423
Step 2f) when con4=3 and con3=2; Collect strip steel related process parameter, comprising: the state variable that whether the water-cooled node sprays water after current each frame of winding steel
Figure RE-FSB00000613083600051
Step 2g) when con4=3, makes each frame roller amount of inclining
Figure RE-FSB00000613083600052
Each frame work roll surface mean temperature
Figure RE-FSB00000613083600053
The comprehensive roll shape C of each frame working roll roller Wi=C Wgi, the comprehensive roll shape C of each frame support roller Bi=C Bgi, temperature during machine on each frame working roll
Figure RE-FSB00000613083600054
Wear extent during machine on each frame working roll
Figure RE-FSB00000613083600055
Wear extent during machine on each frame support roller
Figure RE-FSB00000613083600056
As con4=3 and con2 i, make the comprehensive roll shape C of each frame intermediate calender rolls at=1 o'clock Mi=C Mgi, wear extent during machine on each frame intermediate calender rolls
Figure RE-FSB00000613083600057
Step 2h) when con4=4, collecting belt steel supplied materials head thickness cross direction profiles measured value
Figure RE-FSB00000613083600058
band steel supplied materials current section thickness cross direction profiles measured value
Figure RE-FSB00000613083600059
Step 2i) when con4=5, band steel end frame export goal flatness and lateral distribution
Figure RE-FSB000006130836000510
the band steel end frame outlet of collecting the industrial requirements decision is when front section flatness and lateral distribution measured value
Figure RE-FSB000006130836000511
14. the manner of execution of plate-shape control integrated system according to claim 10; It is characterized in that; In the said step 4) when con4=1 and con1=3, perhaps when con4=2 and con1=3, perhaps when con4=3 and con1=3; Calculate the belt steel temperature field according to rolled piece and roller temperature field module (3), and calculate each stand stretch district average deformation drag according to rolled piece plastic deformation module (1) simultaneously and may further comprise the steps:
Step 3a) so that each rack strip deformation resistance of the actual value?
Figure RE-FSB000006130836000512
is equal to the current volume of each rack strip deformation resistance settings?
Figure RE-FSB000006130836000513
Step 3b) when con3=1, with the supplied materials temperature
Figure RE-FSB000006130836000514
Supplied materials speed
Figure RE-FSB000006130836000515
The supplied materials width B j, the supplied materials average thickness With each frame outlet average thickness
Figure RE-FSB000006130836000517
Each frame coefficient of friction
Figure RE-FSB000006130836000518
Each stand stretch drag actual value of band steel
Figure RE-FSB000006130836000519
Each frame work roll surface mean temperature
Figure RE-FSB000006130836000520
Each frame emulsion temperature T C1~T Cn, air themperature T a, frame spacing L M12~L M (n-1) nAs input parameter, can obtain strip temperature field T according to rolled piece and roller temperature field module (3) jWhen con3=2, with the supplied materials temperature
Figure RE-FSB000006130836000521
Supplied materials speed The supplied materials width B j, the supplied materials average thickness
Figure RE-FSB000006130836000523
With each frame outlet average thickness
Figure RE-FSB000006130836000524
Each frame coefficient of friction
Figure RE-FSB000006130836000525
Each stand stretch drag actual value of band steel
Figure RE-FSB000006130836000526
Each frame work roll surface mean temperature Each frame cooling water temperature T C1~T Cn, air themperature T a, frame spacing L M12~L M (n-1) n, De-scaling box center to the first frame central spacing L T0, the water-cooled node is to this frame central spacing L after each frame T1~L Tn, De-scaling box sprays length L between the area of coverage C0, the water-cooled node sprays length L between the area of coverage after each frame C1~L Cn, whether the water-cooled node sprays water after each frame state variable Bu 1~Bu nAs input parameter, can obtain current winding steel temperature field T according to rolled piece and roller temperature field module (3) j, through to T jIntegration can obtain each frame rolling deformation district rolled piece mean temperature
Figure RE-FSB00000613083600061
Step 3c) with the supplied materials steel grade
Figure RE-FSB00000613083600062
I frame inlet thickness cross direction profiles
Figure RE-FSB00000613083600063
I frame exit thickness cross direction profiles Supplied materials speed
Figure RE-FSB00000613083600065
Supplied materials thickness cross direction profiles h 0, i frame rolling deformation district rolled piece mean temperature
Figure RE-FSB00000613083600066
As input parameter, can obtain the new average deformation drag value of each frame according to rolled piece plastic deformation module (1)
Figure RE-FSB00000613083600067
Calculate the convergence judgment variable
Figure RE-FSB00000613083600068
Step 3d) judges ε 11~ε 1nWhether all smaller or equal to given in advance little value ε 1If, ε 11~ε 1nAll smaller or equal to ε 1, obtain final
Figure RE-FSB00000613083600069
And according to T at this moment jEach frame rolling deformation district rolled piece average surface temperature of integral and calculating
Figure RE-FSB000006130836000610
Calculate and finish; If ε 11~ε 1nNot all smaller or equal to ε 1, order
Figure RE-FSB000006130836000611
Repeating step 3b), step 3c), step 3d), up to ε 11~ε 1nAll smaller or equal to ε 1Till.
15. the manner of execution of plate-shape control integrated system according to claim 10; It is characterized in that; In the said step 5) when con4=1; Perhaps when con4=2, rolled piece plastic deformation module (1) and roll elastic deformation module (2) are coupled, find the solution each frame outlet glacing flatness, thickness and draught pressure cross direction profiles and may further comprise the steps:
Step 4a) suppose i frame exit thickness cross direction profiles wherein
Figure RE-FSB000006130836000613
Figure RE-FSB000006130836000614
be respectively i frame inlet, outlet average thickness,
Figure RE-FSB000006130836000615
is i frame inlet thickness cross direction profiles;
Step 4b) with i stand stretch drag actual value
Figure RE-FSB000006130836000616
I frame coefficient of friction
Figure RE-FSB000006130836000617
I frame inlet thickness cross direction profiles
Figure RE-FSB000006130836000618
I frame exit thickness cross direction profiles
Figure RE-FSB000006130836000619
The supplied materials width B j, i frame inlet tension force
Figure RE-FSB000006130836000620
I frame outlet tension force
Figure RE-FSB000006130836000621
As input parameter, try to achieve the preceding tensile stress cross direction profiles of current winding steel i frame according to rolled piece plastic deformation module (1)
Figure RE-FSB000006130836000622
Current winding steel i frame draught pressure cross direction profiles
Figure RE-FSB000006130836000623
Step 4c) works as con2 i=1 o'clock, with i frame working roll barrel length L Wi, i frame intermediate calender rolls barrel length L Mi, i frame support roller barrel length L Bi, i frame work roll bending cylinder spacing L 1i, i frame intermediate calender rolls roller cylinder spacing L 2i, i frame depress oil cylinder centre-to-centre spacing L 3i, i frame work roll diameter D Wi, i frame intermediate calender rolls diameter D Mi, i frame support roller diameter D Bi, i frame support roller roll neck diameter D BniThe comprehensive roll shape C of i frame working roll Wi, the comprehensive roll shape C of i frame intermediate calender rolls Mi, the comprehensive roll shape C of i frame support roller Bi, i frame work roll bending power I frame intermediate calender rolls bending roller force
Figure RE-FSB000006130836000625
I frame working roll roll shifting amount
Figure RE-FSB000006130836000626
I frame intermediate roll shifting amount
Figure RE-FSB000006130836000627
The i frame roller amount of inclining I frame draught pressure cross direction profiles I frame outlet average thickness
Figure RE-FSB00000613083600072
As input parameter, try to achieve the new exit thickness cross direction profiles of i frame according to roll elastic deformation module (2)
Figure RE-FSB00000613083600073
And i frame working roll corresponding and the distribution of intermediate calender rolls roll gap pressure with it
Figure RE-FSB00000613083600074
I frame intermediate calender rolls and backing roll roll gap pressure distribute
Figure RE-FSB00000613083600075
Work as con2 i=2 o'clock, with i frame working roll barrel length L Wi, i frame support roller barrel length L Bi, i frame work roll bending cylinder spacing L 1i, i frame depress oil cylinder centre-to-centre spacing L 3i, i frame work roll diameter D Wi, i frame support roller diameter D Bi, i frame support roller roll neck diameter D Bni, the comprehensive roll shape C of i frame working roll Wi, the comprehensive roll shape C of i frame support roller Bi, i frame work roll bending power
Figure RE-FSB00000613083600076
I frame working roll roll shifting amount
Figure RE-FSB00000613083600077
The i frame roller amount of inclining
Figure RE-FSB00000613083600078
I breast roller angle of the crossing θ i, i frame draught pressure cross direction profiles
Figure RE-FSB00000613083600079
I frame outlet average thickness
Figure RE-FSB000006130836000710
As input parameter, try to achieve the new exit thickness cross direction profiles of i frame according to roll elastic deformation module (2)
Figure RE-FSB000006130836000711
And i frame working roll corresponding and the distribution of backing roll roll gap pressure with it
Figure RE-FSB000006130836000712
Calculate the convergence judgment variable
Figure RE-FSB000006130836000713
Step 4d) judges ε iWhether smaller or equal to given in advance little value ε, if ε i≤ε obtains the preceding tensile stress cross direction profiles of final i frame
Figure RE-FSB000006130836000714
I frame exit thickness cross direction profiles I frame draught pressure cross direction profiles
Figure RE-FSB000006130836000716
And i frame working roll corresponding and the distribution of intermediate calender rolls roll gap pressure with it
Figure RE-FSB000006130836000717
I frame intermediate calender rolls and backing roll roll gap pressure distribute Or i frame working roll and the distribution of backing roll roll gap pressure If ε i>ε, order
Figure RE-FSB000006130836000720
α is a relaxation factor, repeating step 4b), step 4c), step 4d), up to ε iTill≤the ε.
16. the manner of execution of plate-shape control integrated system according to claim 10; It is characterized in that; In the said step 6) when con4=3, calculate each frame target glacing flatness curve and target thickness cross direction profiles curve
Figure RE-FSB000006130836000722
according to plate-shape standard curve module (7) and may further comprise the steps:
Step 5a) makes i=1;
Step 5b) sets i frame target glacing flatness curve and target thickness cross direction profiles curvilinear equation;
Each frame target glacing flatness curvilinear equation is:
Figure RE-FSB000006130836000723
Y in the formula---band steel horizontal relative coordinate, to another side, y from-1 to+1 from one side of band steel;
a 0, a 2, a 4---fitting coefficient;
Each frame target thickness cross direction profiles curvilinear equation:
Figure RE-FSB00000613083600081
B in the formula 0, b 2, b 4---fitting coefficient;
Step 5c), change step 5d over to if i is less than or equal to n-1), otherwise, change step 5f over to);
Step 5d) makes a 4=0, adopt the given unknowm coefficient a of certain optimized Algorithm 2, with supplied materials width B, i frame outlet average thickness
Figure RE-FSB00000613083600082
Tensile stress cross direction profiles σ before the i frame 1iInhomogeneous composition
Figure RE-FSB00000613083600083
As input parameter, try to achieve plate shape according to the good discrimination module of glacing flatness (5) and differentiate factor ξ i, the condition of convergence is for making ξ iThe convergence setting value judges whether to satisfy the condition of convergence, as satisfies the condition of convergence as far as possible, obtains i passage target glacing flatness curve standard value, is designated as
Figure RE-FSB00000613083600084
As do not satisfy the condition of convergence, repeated execution of steps 5d), till satisfying the condition of convergence;
Step 5e) adopt certain optimized Algorithm given
Figure RE-FSB00000613083600085
Unknowm coefficient b 2, b 4, with i stand stretch drag setting value
Figure RE-FSB00000613083600086
I frame coefficientoffriction i, i frame inlet thickness cross direction profiles
Figure RE-FSB00000613083600087
I frame exit thickness cross direction profiles
Figure RE-FSB00000613083600088
Supplied materials width B, i frame inlet tension force T I-1, i frame outlet tension force T iAs input parameter, try to achieve preceding tensile stress cross direction profiles σ according to rolled piece plastic deformation module (1) 1iInhomogeneous composition
Figure RE-FSB00000613083600089
The condition of convergence is tensile stress cross direction profiles σ before making 1iInhomogeneous composition
Figure RE-FSB000006130836000810
Convergence i passage target glacing flatness curve standard value as far as possible
Figure RE-FSB000006130836000811
Judge whether to satisfy the condition of convergence, as satisfy the condition of convergence, obtain i passage target glacing flatness curve
Figure RE-FSB000006130836000812
With target thickness cross direction profiles curve
Figure RE-FSB000006130836000813
As do not satisfy the condition of convergence, repeated execution of steps 5e), till satisfying the condition of convergence;
Step 5f) according to back step operation two unknowm coefficient a is set in the specific (special) requirements of plate shape 2, a 4, obtain n passage target glacing flatness curve
Figure RE-FSB000006130836000814
Step 5g) adopt certain optimized Algorithm given
Figure RE-FSB000006130836000815
Unknowm coefficient b 2, b 4, with n stand stretch drag setting value
Figure RE-FSB000006130836000816
N frame coefficientoffriction n, n frame inlet thickness cross direction profiles
Figure RE-FSB000006130836000817
N frame exit thickness cross direction profiles
Figure RE-FSB000006130836000818
Supplied materials width B, n frame inlet tension force T N-1, n frame outlet tension force T nAs input parameter, try to achieve preceding tensile stress cross direction profiles σ according to rolled piece plastic deformation module (1) 1nInhomogeneous composition
Figure RE-FSB000006130836000819
The condition of convergence is tensile stress cross direction profiles σ before making 1nInhomogeneous composition
Figure RE-FSB000006130836000820
Convergence n passage target glacing flatness curve as far as possible Judge whether to satisfy the condition of convergence, as satisfy the condition of convergence, obtain n passage target glacing flatness curve
Figure RE-FSB000006130836000822
With target thickness cross direction profiles curve
Figure RE-FSB000006130836000823
As do not satisfy the condition of convergence, repeated execution of steps 5g), till satisfying the condition of convergence;
Step 5h) make i=i+1, change step 5b over to), up to i=n.
17. the manner of execution of plate-shape control integrated system according to claim 10 is characterized in that, in the said step 7) when con4=3, each frame plate shape control device is carried out optimizing calculates, obtain each frame plate shape setup parameter and may further comprise the steps:
Step 6a) makes i=1;
Step 6b) when con4=3 and con3=1; Exporting actual glacing flatness curve with each frame, to equal target glacing flatness curve be target; When con4=3 and con3=2; Equaling target thickness cross direction profiles curve with each frame outlet actual (real) thickness cross direction profiles curve is target, adopts optimized Algorithm given one group of concrete parameter in the upper and lower scope of plate shape control parameter, specifically comprises: each frame work roll bending power
Figure RE-FSB00000613083600091
Each frame working roll roll shifting amount
Figure RE-FSB00000613083600092
Con2 i=1 o'clock frame intermediate calender rolls bending roller force set-point
Figure RE-FSB00000613083600093
Con2 i=1 o'clock frame intermediate roll shifting amount
Figure RE-FSB00000613083600094
Con2 i=2 o'clock breast roller angle of the crossing set-point θ i, change step 6c over to);
Step 6c) suppose i frame exit thickness cross direction profiles
Figure RE-FSB00000613083600095
wherein
Figure RE-FSB00000613083600096
be respectively i frame inlet, outlet average thickness,
Figure RE-FSB00000613083600097
is i frame inlet thickness cross direction profiles;
Step 6d) with i stand stretch drag actual value
Figure RE-FSB00000613083600098
I frame coefficient of friction
Figure RE-FSB00000613083600099
I frame inlet thickness cross direction profiles
Figure RE-FSB000006130836000910
I frame exit thickness cross direction profiles
Figure RE-FSB000006130836000911
The supplied materials width B j, i frame inlet tension force
Figure RE-FSB000006130836000912
I frame outlet tension force T i jAs input parameter, try to achieve the preceding tensile stress cross direction profiles of current winding steel i frame according to rolled piece plastic deformation module (1)
Figure RE-FSB000006130836000913
Current winding steel i frame draught pressure cross direction profiles
Figure RE-FSB000006130836000914
Step 6e) works as con2 i=1 o'clock, with i frame working roll barrel length L Wi, i frame intermediate calender rolls barrel length L Mi, i frame support roller barrel length L Bi, i frame work roll bending cylinder spacing L 1i, i frame intermediate calender rolls roller cylinder spacing L 2i, i frame depress oil cylinder centre-to-centre spacing L 3i, i frame work roll diameter D Wi, i frame intermediate calender rolls diameter D Mi, i frame support roller diameter D Bi, i frame support roller roll neck diameter D Bni, the comprehensive roll shape C of i frame working roll WiThe comprehensive roll shape C of i frame intermediate calender rolls Mi, the comprehensive roll shape C of i frame support roller Bi, i frame work roll bending power
Figure RE-FSB000006130836000915
I frame intermediate calender rolls bending roller force
Figure RE-FSB000006130836000916
I frame working roll roll shifting amount
Figure RE-FSB000006130836000917
I frame intermediate roll shifting amount The i frame roller amount of inclining
Figure RE-FSB000006130836000919
I frame draught pressure cross direction profiles
Figure RE-FSB00000613083600101
I frame outlet average thickness
Figure RE-FSB00000613083600102
As input parameter, try to achieve the new exit thickness cross direction profiles of i frame according to roll elastic deformation module (2)
Figure RE-FSB00000613083600103
And i frame working roll corresponding and the distribution of intermediate calender rolls roll gap pressure with it
Figure RE-FSB00000613083600104
I frame intermediate calender rolls and backing roll roll gap pressure distribute
Figure RE-FSB00000613083600105
Work as con2 i=2 o'clock, with i frame working roll barrel length L Wi, i frame support roller barrel length L Bi, i frame work roll bending cylinder spacing L 1i, i frame depress oil cylinder centre-to-centre spacing L 3i, i frame work roll diameter D Wi, i frame support roller diameter D Bi, i frame support roller roll neck diameter D Bni, the comprehensive roll shape C of i frame working roll Wi, the comprehensive roll shape C of i frame support roller Bi, i frame work roll bending power
Figure RE-FSB00000613083600106
I frame working roll roll shifting amount
Figure RE-FSB00000613083600107
The i frame roller amount of inclining I breast roller angle of the crossing θ i, i frame draught pressure cross direction profiles
Figure RE-FSB00000613083600109
I frame outlet average thickness
Figure RE-FSB000006130836001010
As input parameter, try to achieve the new exit thickness cross direction profiles of i frame according to roll elastic deformation module (2)
Figure RE-FSB000006130836001011
And i frame working roll corresponding and the distribution of backing roll roll gap pressure with it
Figure RE-FSB000006130836001012
Calculate the convergence judgment variable
Figure RE-FSB000006130836001013
Step 6f) judges ε 2iWhether smaller or equal to given in advance little value ε 2If, ε 2i≤ε 2, obtain final
Figure RE-FSB000006130836001014
And it is corresponding with it
Figure RE-FSB000006130836001015
Or
Figure RE-FSB000006130836001016
Change step 6g over to), if ε 2i>ε 2, order α is a relaxation factor, repeating step 6d), step 6e), step 6f), up to ε 2i≤ε 2Till;
Step 6g) when con4=3, judge whether to satisfy the computation optimization condition of convergence, as satisfy the condition of convergence, write down the optimum plate shape of current winding steel i frame setup parameter, change step 6h over to); As do not satisfy the condition of convergence, repeating step 6b)~step 6g), till satisfying;
Step 6h) judge whether i=n,, change step 6b over to if i ≠ n makes i=i+1), up to i=n.
18. the manner of execution of plate-shape control integrated system according to claim 10; It is characterized in that; In the said step 8) when con4=3; Calculate rolling each the frame work roll thermal roll shape of back that finishes of current volume according to rolled piece and roller temperature field module (3), calculate rolling each the frame working roller abrasion roll shape of back that finishes of current volume according to roller system's wearing and tearing modules (4); May further comprise the steps:
Step 7a) makes i=1;
Step 7b) with i frame working roll temperature field behind the rolling j-1 winding steel J winding steel i frame rolling deformation district rolled piece average surface temperature
Figure RE-FSB00000613083600111
Pure the rolling the time of j winding steel J winding steel rolls the back intermittent time
Figure RE-FSB00000613083600113
I frame work roll diameter D WiAs input parameter, try to achieve i frame working roll temperature field behind the rolling j winding steel by rolled piece and roller temperature field module (3)
Figure RE-FSB00000613083600114
Step 7c) with i frame working roll temperature field
Figure RE-FSB00000613083600115
I frame working roll is initially gone up the machine temperature T Wi, i frame work roll diameter D WiAs input parameter, try to achieve i frame work roll thermal roll shape behind the rolling j winding steel by roll elastic deformation module (2)
Figure RE-FSB00000613083600116
Step 7d) judge whether i equals n,, change step 7b over to if i ≠ n makes i=i+1), if i=n changes step 7e over to);
Step 7e) makes i=1;
Step 7f) works as con2 i=1 o'clock, with i frame working roller abrasion roll shape behind the rolling j-1 winding steel
Figure RE-FSB00000613083600117
I frame intermediate calender rolls wearing and tearing roll shape behind the rolling j-1 winding steel
Figure RE-FSB00000613083600118
I frame support roller wearing and tearing roll shape behind the rolling j-1 winding steel J winding steel supplied materials length J winding steel supplied materials thickness cross direction profiles
Figure RE-FSB000006130836001111
J winding steel i frame exit thickness cross direction profiles
Figure RE-FSB000006130836001112
J winding steel i frame draught pressure cross direction profiles
Figure RE-FSB000006130836001113
J winding steel i frame working roll and intermediate calender rolls roll gap pressure distribute
Figure RE-FSB000006130836001114
I frame intermediate calender rolls and backing roll roll gap pressure distribute As input parameter, try to achieve i frame working roller abrasion roll shape behind the rolling j winding steel by roller system wearing and tearing modules (4)
Figure RE-FSB000006130836001116
I frame intermediate calender rolls wearing and tearing roll shape behind the rolling j winding steel
Figure RE-FSB000006130836001117
I frame support roller wearing and tearing roll shape behind the rolling j winding steel
Figure RE-FSB000006130836001118
Work as con2 i=2 o'clock, with i frame working roller abrasion roll shape behind the rolling j-1 winding steel
Figure RE-FSB000006130836001119
I frame support roller wearing and tearing roll shape behind the rolling j-1 winding steel J winding steel supplied materials length
Figure RE-FSB000006130836001121
J winding steel supplied materials thickness cross direction profiles
Figure RE-FSB000006130836001122
J winding steel i frame exit thickness cross direction profiles J winding steel i frame draught pressure cross direction profiles
Figure RE-FSB000006130836001124
I frame working roll and backing roll roll gap pressure distribute As input parameter, try to achieve i frame working roller abrasion roll shape behind the rolling j winding steel by roller system wearing and tearing modules (4)
Figure RE-FSB000006130836001126
I frame support roller wearing and tearing roll shape behind the rolling j winding steel
Figure RE-FSB000006130836001127
Step 7g) judge whether i equals n,, change step 7f over to if i ≠ n makes i=i+1), up to i=n.
19. the manner of execution of plate-shape control integrated system according to claim 10; It is characterized in that; In the said step 9) when con4=1; Perhaps when con4=4; Perhaps when con4=5; According to plate shape pattern recognition module (6) glacing flatness or thickness cross direction profiles curve are carried out pattern-recognition, calculate each frame plate shape control device regulated quantity matrix
Figure RE-FSB000006130836001128
according to plate shape control module (8) and may further comprise the steps:
Step 8a) when con4=1, makes i=1;
Step 8b) when con4=1; The i frame is exported preceding tensile stress cross direction profiles
Figure RE-FSB00000613083600121
as input parameter; By plate shape pattern recognition module (6) try to achieve rolling j winding steel i frame outlet strip profile and flatness cross direction profiles 1 time, 2 times, 3 times, 4 sub-eigenvalues
Figure RE-FSB00000613083600122
with i frame exit thickness cross direction profiles as input parameter, by plate shape pattern recognition module (6) try to achieve rolling j winding steel i frame outlet belt steel thickness cross direction profiles 1 time, 2 times, 3 times, 4 sub-eigenvalues
Figure RE-FSB00000613083600124
Step 8c) judge when con4=1 whether i equals the total frame of tandem mill and count n, if i ≠ n makes i=i+1, changes step 8b over to), up to i=n;
Step 8d) when con4=4; To be with the current section thickness of steel supplied materials cross direction profiles measured value
Figure RE-FSB00000613083600125
band steel supplied materials head thickness cross direction profiles measured value
Figure RE-FSB00000613083600126
as input parameter respectively, by formula plate shape pattern recognition module (6) try to achieve the current section thickness of rolling j winding steel supplied materials cross direction profiles 1 time, 2 times, 3 times, 4 sub-eigenvalues
Figure RE-FSB00000613083600127
and rolling j winding steel supplied materials head thickness cross direction profiles 1 time, 2 times, 3 times, 4 sub-eigenvalues
Figure RE-FSB00000613083600128
and calculate the deviation
Figure RE-FSB00000613083600129
of corresponding sub-eigenvalue
Step 8e) when con4=4 and con3=1; With
Figure RE-FSB000006130836001210
as input parameter; Try to achieve each frame plate shape control device regulated quantity matrix
Figure RE-FSB000006130836001211
when con4=4 and the con3=2 by plate shape control module (8);
Figure RE-FSB000006130836001212
as input parameter, tried to achieve each frame plate shape control device regulated quantity matrix
Figure RE-FSB000006130836001213
by plate shape control module (8)
Step 8f) when con4=5; Flatness and lateral distributes
Figure RE-FSB000006130836001215
as input parameter will to export the band steel end frame export goal that determines when front section flatness and lateral distribution measured value
Figure RE-FSB000006130836001214
industrial requirements with steel end frame respectively, tries to achieve the deviation
Figure RE-FSB000006130836001218
that the last frame outlet of rolling j winding steel is worked as j winding steel end frame export goal flatness and lateral distribution 1 time, 2 times, 3 times, 4 sub-eigenvalues
Figure RE-FSB000006130836001217
of front section flatness and lateral distribution measured value 1 time, 2 times, 3 times, 4 sub-eigenvalues
Figure RE-FSB000006130836001216
and industrial requirements decision and calculated corresponding sub-eigenvalue by plate shape pattern recognition module (6)
Step 8g) when con4=5;
Figure RE-FSB000006130836001219
as input parameter, tried to achieve last frame plate shape control device regulated quantity matrix by plate shape control module (8)
20. the manner of execution of plate-shape control integrated system according to claim 10; It is characterized in that: in the said step 10) when con4=3, rolling back each frame hot-rolling type, roll wear roll shape and the original grinding roll shape of roll of finishing of current volume superposeed according to following mode:
Work as con2 i=1 o'clock, order
Figure RE-FSB00000613083600131
Figure RE-FSB00000613083600132
Figure RE-FSB00000613083600133
C wherein WiBe the comprehensive roll shape of i frame working roll, C MiBe the comprehensive roll shape of i frame intermediate calender rolls, C BiBe the comprehensive roll shape of i frame support roller, C WgiBe the original grinding roll shape of i frame working roll, C MgiBe the original grinding roll shape of i frame intermediate calender rolls, C BgiBe the original grinding roll shape of i frame support roller,
Figure RE-FSB00000613083600134
Be i frame working roller abrasion roll shape behind the rolling j winding steel,
Figure RE-FSB00000613083600135
Be i frame intermediate calender rolls wearing and tearing roll shape behind the rolling j winding steel,
Figure RE-FSB00000613083600136
Be i frame support roller wearing and tearing roll shape behind the rolling j winding steel,
Figure RE-FSB00000613083600137
I frame work roll thermal roll shape behind the rolling j winding steel;
Work as con2 i=2 o'clock, order
Figure RE-FSB00000613083600138
C wherein WiBe the comprehensive roll shape of i frame working roll, C BiBe the comprehensive roll shape of i frame support roller, C WgiBe the original grinding roll shape of i frame working roll, C BgiBe the original grinding roll shape of i frame support roller,
Figure RE-FSB000006130836001310
Be i frame working roller abrasion roll shape behind the rolling j winding steel,
Figure RE-FSB000006130836001311
Be i frame support roller wearing and tearing roll shape behind the rolling j winding steel,
Figure RE-FSB000006130836001312
I frame work roll thermal roll shape behind the rolling j winding steel.
CN2009100750825A 2009-08-05 2009-08-05 Plate-shape control integrated system and executing method thereof Active CN101623708B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100750825A CN101623708B (en) 2009-08-05 2009-08-05 Plate-shape control integrated system and executing method thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100750825A CN101623708B (en) 2009-08-05 2009-08-05 Plate-shape control integrated system and executing method thereof

Publications (2)

Publication Number Publication Date
CN101623708A CN101623708A (en) 2010-01-13
CN101623708B true CN101623708B (en) 2012-02-29

Family

ID=41519786

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100750825A Active CN101623708B (en) 2009-08-05 2009-08-05 Plate-shape control integrated system and executing method thereof

Country Status (1)

Country Link
CN (1) CN101623708B (en)

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102161054B (en) * 2010-12-24 2012-12-12 燕山大学 Plate shape closed-loop control method based on influencing self learning of matrix
CN102896459B (en) * 2011-07-26 2014-12-17 上海宝钢工业技术服务有限公司 Online repair method for window rolling stress surfaces of rack archways of hot continuous rolling units
CN102581029B (en) * 2011-08-18 2014-05-07 苏州有色金属研究院有限公司 Method for determining effect of strip shape execution mechanism of metal strips
CN103191932B (en) * 2012-01-09 2015-03-04 宝山钢铁股份有限公司 Detection method for strip steel tension fluctuation between racks of hot continuous rolling
CN102581027B (en) * 2012-01-18 2014-02-12 中冶南方工程技术有限公司 Control system for transverse integrative optimization of shape of cold rolled steel strip
CN102989787B (en) * 2012-12-21 2014-09-03 鞍钢股份有限公司 Hot-rolled strip steel head and tail plate shape control method
CN104096714B (en) * 2013-04-11 2016-06-29 宝山钢铁股份有限公司 A kind of hot-strip convexity autocontrol method
CN103341505B (en) * 2013-06-25 2015-01-28 中冶南方工程技术有限公司 Cold-rolled sheet shape signal online mode recognition method
CN105013835B (en) * 2014-04-23 2017-06-06 宝山钢铁股份有限公司 Original roll gap establishing method based on roll thermal crown in the rolling of tandem mills strip in razor-thin
CN105032945B (en) * 2015-07-08 2017-03-29 燕山大学 A kind of hot tandem plate shape and plate convexity Comprehensive Control merit rating method
CN106886191A (en) * 2015-12-16 2017-06-23 鞍钢股份有限公司 High-speed data acquisition method for hot galvanizing production process of strip steel
US10994317B2 (en) 2016-03-08 2021-05-04 Novelis Inc. Method and apparatus for controlling metal strip profile during rolling with direct measurement of process parameters
CN105867312B (en) * 2016-04-12 2018-07-06 燕山大学 A kind of Strip Shape intelligent remote cloud TT&C system
CN106944483B (en) * 2017-03-22 2018-11-09 宁波宝新不锈钢有限公司 A kind of measurement control method of reversing cold mill band transverse gage Distribution Value
CN108500067B (en) * 2018-04-11 2020-02-04 攀钢集团攀枝花钢钒有限公司 Strip shape adjusting method for hot continuous rolling strip steel
CN112934978A (en) * 2019-12-10 2021-06-11 宝山钢铁股份有限公司 Technological method for controlling thermal convexity of working roll of front frame of secondary cold rolling unit
CN111856932B (en) * 2020-07-02 2022-05-27 唐山钢铁集团有限责任公司 Cold strip mill plate shape closed-loop control method based on influence matrix recursive identification
CN111983983B (en) * 2020-08-07 2023-10-10 东北大学 Intelligent control method and system for heat treatment plate shape of hot rolled steel plate
CN112893480B (en) * 2021-01-18 2021-11-23 燕山大学 Optimization method of non-square plate-shaped control system

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000343112A (en) * 1999-03-31 2000-12-12 Kawasaki Steel Corp Shape controlling method of plate rolling and device therefor
CN1485156A (en) * 2002-09-28 2004-03-31 鞍钢集团新钢铁有限责任公司 Hot-rolled strip steel moderate sea profile shape controlling method
CN1923392A (en) * 2005-08-30 2007-03-07 宝山钢铁股份有限公司 Planisher draught pressure setting method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000343112A (en) * 1999-03-31 2000-12-12 Kawasaki Steel Corp Shape controlling method of plate rolling and device therefor
CN1485156A (en) * 2002-09-28 2004-03-31 鞍钢集团新钢铁有限责任公司 Hot-rolled strip steel moderate sea profile shape controlling method
CN1923392A (en) * 2005-08-30 2007-03-07 宝山钢铁股份有限公司 Planisher draught pressure setting method

Also Published As

Publication number Publication date
CN101623708A (en) 2010-01-13

Similar Documents

Publication Publication Date Title
CN101623708B (en) Plate-shape control integrated system and executing method thereof
Cao et al. Integrated design of roll contours for strip edge drop and crown control in tandem cold rolling mills
CN107377634B (en) A kind of hot-strip outlet Crown Prediction of Media method
CN101648216B (en) Method for setting plate shape and plate convexity off-line prediction of PC rolling mill
CN101653786B (en) Method for improving control accuracy of hot rolling width
CN102029294B (en) Control method of lateral thick difference of cold-rolling strip steels
CN102513351B (en) Rolling method and device for strip steel tandem cold rolling
Peng et al. Strip shape modeling and its setup strategy in hot strip mill process
CN102236322A (en) Roller wear optimization control method for improving quality of band steel
CN103567229B (en) Roll bending force combination plate shape control method for six-roll rolling mill
CN104998913A (en) Prediction method of motor power in cold rolling process
CN1091008C (en) Interlinked control method for plate-band rolling course based on coordination law of plate shape and plate thickness
CN102161054B (en) Plate shape closed-loop control method based on influencing self learning of matrix
CN102688896A (en) Mechanism model-based virtual crown gauge setting method for four-roll cold continuous rolling machine set
CN106345817A (en) Plate shape prediction method of six-roller mill comprising working roller with large length-diameter ratio and small diameter absolute value
CN105013835A (en) Original roller seam setting method based on thermal crown in ultra-thin strip rolling conducted by cold continuous rolling unit
CN116329297B (en) Plate shape prediction method based on transverse mechanical property difference of rolled piece
CN105290122B (en) A kind of AGC system measurer for thickness of RBF network supports
CN101898202B (en) Method for forecasting edge reduction in rolling strips of SMS-EDC rolling mill
Zhu et al. A fuzzy algorithm for flatness control in hot strip mill
CN1179369A (en) Plate shape measuring and controlling method for plate strip rolling process
CN106825068A (en) A kind of Forecasting Methodology of operation of rolling belt steel surface roughness
Ma et al. A new model for thermo-mechanical coupled analysis of hot rolling
CN114798755B (en) Method for analyzing influence of hot crown on plate shape in cold continuous rolling process
Han et al. Prediction and control of profile for silicon steel strip in the whole tandem cold rolling based on PSO-BP algorithm

Legal Events

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