[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN102688896B - Mechanism model-based virtual crown gauge setting method for four-roll cold continuous rolling machine set - Google Patents

Mechanism model-based virtual crown gauge setting method for four-roll cold continuous rolling machine set Download PDF

Info

Publication number
CN102688896B
CN102688896B CN201210150863.8A CN201210150863A CN102688896B CN 102688896 B CN102688896 B CN 102688896B CN 201210150863 A CN201210150863 A CN 201210150863A CN 102688896 B CN102688896 B CN 102688896B
Authority
CN
China
Prior art keywords
inequality
rolling
value
band
calculate
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.)
Expired - Fee Related
Application number
CN201210150863.8A
Other languages
Chinese (zh)
Other versions
CN102688896A (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.)
Tangshan City Delong Iron and Steel Co Ltd
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 CN201210150863.8A priority Critical patent/CN102688896B/en
Publication of CN102688896A publication Critical patent/CN102688896A/en
Application granted granted Critical
Publication of CN102688896B publication Critical patent/CN102688896B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Metal Rolling (AREA)

Abstract

The invention relates to a mechanism model-based virtual crown gauge setting method for four-roll cold continuous rolling machine set. The setting method, executed by a computer, comprises (a) collecting fundamental equipment parameters of each rolling mill of the cold continuous rolling machine set; (b) realizing the outlet section shape display function; and (c) realizing the section characteristic closed-loop feedback control function. The method provided by the invention realizes accurate calculation on the cold rolled finished product section characteristic, and performs feedback regulation to maximally ensure the quality of the cold rolled finished product in respect of plate thickness.

Description

The virtual profile gauge establishing method of four roller tandem mills based on mechanism model
Technical field
The invention belongs to the Computer Control Technology field in metallurgy industry, particularly the virtual profile gauge establishing method of a kind of four roller tandem mills based on mechanism model.
Background technology
In cold continuous rolling process, supplied materials strip crown has material impact to cold rolling production board shape and cross section characteristic parameter etc., thereby is directly connected to the quality of cold-rolled products.In the prior art, although normal four-roller cold-rolling mill has work roll bending, the roller equal section shape adjustment of inclining means, but because the initial cross section characteristic of supplied materials can not accurately be predicted in production scene, so in the actual operation of rolling, cannot realize the effective control to strip material plate convexity, wedge shape amount equal section characterisitic parameter.Based on this reason, conventionally the basic configuration that has kept supplied materials section in cold-rolled process, if supplied materials has strip crown or the defect such as wedge shape amount is overproof, cold rolling finished product is " heredity " these defects possibly, and cause thus the product quality objection of lower procedure (as continuous annealing etc.), the even production accident such as sideslip, broken belt.On the one hand, user is also more and more higher to the requirement of Strip Shape and thickness of slab in order.So, how in cold-rolled process, the cross section characteristic of band to be control effectively and just to become the emphasis of site technology tackling key problem.
Be accompanied by the development of the present computer technology, corresponding hot-rolled finished product convexity reproduction and analytical technology are ripe gradually, and some up-to-date scientific and technological achievements can realize transferring at any time and analyzing the cross section characteristic data of hot-rolled finished product supplied materials.Like this, in conjunction with above supplied materials related data, and by mechanism model, calculate final finished product cross section characteristic and just become possibility.Simultaneously, because in front two frames, belt steel thickness is relatively large, resistance of deformation is less, and it is less on the impact of plate shape to depress process, therefore can be by contrast the second frame outlet actual ratio convexity and target proportion convexity, and the rolling technological parameter such as the work roll bending power of these two frames of feedback regulation, the roller amount of inclining, roll shifting amount, thereby realize the FEEDBACK CONTROL to strip section characteristic.
Summary of the invention
The object of the invention is cannot effectively accurately control for cross section characteristic in cold continuous rolling process in the past on-the-spot, thereby easily cause the problem of quality objection or downstream production accident, the virtual profile gauge establishing method of a kind of four roller tandem mills based on mechanism model is provided, this method can realize the accurate calculating to cold rolling finished product cross section characteristic, and carry out feedback regulation, thereby at utmost guarantee the thickness of slab quality of cold rolling finished product.
In order to realize above object, the present invention is through a large amount of field trial and theoretical research, on the basis that can transfer at any time in supplied materials cross section characteristic parameter, take normal four-roller tandem mills as research object, utilize outlet of cold-rolling mill band thickness of slab distribution machine reason model, the virtual profile gauge establishing method of a kind of four roller tandem mills based on mechanism model proposed, can realize the accurate calculating to cold rolling finished product cross section characteristic, and carry out feedback regulation, thereby at utmost guaranteed the thickness of slab quality of cold rolling finished product.By technology of the present invention, can realize following two large functions: (1) Outlet Section shape dynamic display of functional Real time dynamic display is calculated band Outlet Section shape, strip crown, the wedge shape amount of gained based on mechanism model, and corresponding rolling technological parameter etc.; (2) cross section characteristic feedback control function, by the second frame outlet actual ratio convexity and the target proportion convexity of contrast in real time, carries out FEEDBACK CONTROL to the work roll bending power of first, second frame, the roller amount of inclining, roll shifting amount etc.
List of references:
[1] Lian Jiachuan, Liu Hongmin. gauge and shape is controlled [M]. the .1995 of weapon industry publishing house
[2] Wang Guodong. plate shape is controlled and plate shape theoretical [M]. the .1986 of metallurgical industry publishing house
The concrete technical scheme that the present invention adopts is as follows:
The virtual profile gauge establishing method of four roller tandem mills based on mechanism model, comprises following content (the function layout drawing of this virtual profile gauge is shown in accompanying drawing 1):
(a) collect the basic equipment parameter of each rolling mill of tandem mills, comprise the barrel length L of working roll wk, backing roll barrel length L bk, working roll bending cylinder is apart from l wk, backing roll housing screw centre-to-centre spacing l bk, working roll roller footpath D wk, backing roll roller footpath D bk, working roll roll shape D wki, backing roll roll shape D bki, the maximum positive bending roller force of working roll
Figure BDA00001644947300021
the maximum negative bending roller force of working roll
Figure BDA00001644947300022
the maximum roller amount of just inclining
Figure BDA00001644947300023
the maximum negative roller amount of inclining subscript k represents shelf number, k=1, and 2 ..., s, s is total frame number;
(b) realization of Outlet Section shape Presentation Function, comprises the following step that can be carried out by computer (fundamental block diagram is shown in accompanying drawing 2):
B1) collect the sampling period τ of current tandem mills data collecting system s;
B2) provide the period tau that virtual profile gauge section configuration dynamically shows x, this cycle is relevant to dominant frequency, the capacity of the computer of the virtual profile gauge of installation, and must guarantee τ xτ sintegral multiple and within this cycle, can complete twice thickness of slab and distribute and calculate;
B3) provide the interval parameter L to supplied materials sectional sampling 0;
B4) definition length of interval parameter j, supplied materials coil of strip length parameter L j;
B5) definition as-rolled condition parameter ξ, wherein ξ=1 represents to start rolling, ξ=-1 represents to stop rolling, and this instruction is sent according to on-site actual situations by site operation personnel;
B6) given initial steel reel number COILNO1=0, prepares convexity and dynamically shows;
B7) from the three-level system of tandem mills, collect supplied materials coil of strip essential information, comprise reel number COILNO, steel grade code SG, strip width B, supplied materials thickness h 0, last rolling mill exit thickness h s;
B8) judge inequality | does COILNO1-COILNO|>0 set up? if set up, make j=0, L j=0, and read the cross section characteristic data that supplied materials steel reel number is COILNO, proceed to step b9); If inequality | COILNO1-COILNO|>0 is false, and directly proceeds to step b9);
B9) collect current time the first frame entrance and roll band steel length L *;
B10) judgement inequality L *>=L jwhether set up, if set up, turn as step b11; Otherwise, proceed to step b12;
B11) extract and calculate supplied materials and twist in [L j, L j+1] average crown curve function f in length of interval j=(X).L wherein j+1=L j+ L 0; X is illustrated in strip width direction and sets up one-dimensional coordinate system, when the center of width is the origin of coordinates, and the corresponding numerical value of strip width direction each point in coordinate system;
B12) by tandem mills data collecting system, collect the rolling technological parameter actual value ψ of current each frame ky, wherein, y represents to survey the classification of dynamic data, and y=1 represents that draught pressure, y=2 represent that forward pull, y=3 represent that backward pull, y=4 represent that reduction ratio, y=5 represent that roll shifting amount, y=6 represent bending roller force, the y=7 roller amount that represents to incline;
B13) utilizing the described outlet thickness of slab of document [1] distribution machine reason model to calculate actual rolling technological parameter is ψ kyand in supplied materials interval, average crown curve function is f j=(X) time, width is the Outlet Section thickness distribution value h of B, the steel grade code band that is SG i, i is band in horizontal bar unit number, i=1, and 2 ..., n, n is the total bar unit number of band;
B14) utilize and calculate a mouthful section thickness Distribution Value h ithe current section crown curve of matching function f c(X), X is illustrated in strip width direction and sets up one-dimensional coordinate system, when the center of width is the origin of coordinates, and the corresponding numerical value of strip width direction each point in coordinate system;
B15) calculate current section convex value CROWN and wedge shape value WEDGE, convex value computation model is CROWN = f c ( 0 ) - f c ( B / 2 - 40 ) + f c ( - B / 2 + 40 ) 2 , Wedge shape amount computation model is
Figure BDA00001644947300042
f wherein c(0) represent band steel center thickness, f c(B/2-40) represent the thickness apart from band 40mm place, steel one end, f c(B/2+40) represents the thickness apart from band steel other end 40mm place;
B16) utilize the dynamic display of functional of visual software, by current section crown curve function f c(X) dynamically show, show convex value CROWN, wedge shape value WEDGE simultaneously, and corresponding coil of strip information;
B17) do you judge that inequality ξ < 0 sets up? if inequality is set up, finish cross section characteristic and show, virtual plate profile instrument quits work; If inequality is false, COILNO1=COILNO, proceeds to step b7), until inequality ξ < 0 sets up.
(c) cross section characteristic closed loop feedback is controlled the realization of function, comprises the following step that can be carried out by computer (fundamental block diagram is shown in accompanying drawing 3):
C1) collect the target proportion crown curve coefficient a of 2# milling train outlet ε, wherein, ε=0,1,2,3,4;
C2) determine respectively in cross section characteristic feedback procedure the minimum coefficient φ that adjusts of single step of bending roller force in 1# and 2# frame 1, φ 2;
C3) determine respectively the minimum coefficient of adjusting of single step of the roller amount of inclining in 1# and 2# frame in cross section characteristic feedback procedure
Figure BDA00001644947300051
C4) determine the period tau of cross section characteristic feedback f, this cycle is relevant to dominant frequency, the capacity of the computer of the virtual profile gauge of installation, and must guarantee τ fτ sintegral multiple and within this cycle, can complete inferior thickness of slab distributes and calculates;
C5) definition rolling moment procedure parameter t 0, virtual profile gauge cross section characteristic Times of Feedback procedure parameter M;
C6) receive operational order, judge whether to start rolling? ξ=1 represents to start rolling, ξ=-1 represents to stop rolling.If started rolling, record current standard Beijing time t, and make t 0=t, proceeds to step c7); If do not start rolling, enter wait state;
C7) given initial steel reel number COILNO1=0;
C8) from the three-level system of tandem mills, collect supplied materials coil of strip essential information, comprise reel number COILNO, steel grade code SG, strip width B, supplied materials thickness h 0, last rolling mill exit thickness h sdeng;
C9) judge inequality | does COILNO1-COILNO|>0 set up? if set up, make j=0, L j=0, and read the cross section characteristic data that supplied materials steel reel number is COILNO, proceed to step c10; If inequality | COILNO1-COILNO|>0 is false, and directly proceeds to step c10;
C10) collect current time the first frame entrance and roll band steel length L *;
C11) judgement inequality L *>=L jset up? if set up, turn as step c12; Otherwise, proceed to step c13;
C12) extract and calculate supplied materials and twist in [L j, L j+1] average crown curve function f in length of interval mj(X).L wherein j+1=L j+ L 0; X is illustrated in strip width direction and sets up one-dimensional coordinate system, when the center of width is the origin of coordinates, and the corresponding numerical value of strip width direction each point in coordinate system;
C13) by milling train data collecting system, collect 1#, 2# milling train from t 0to t 0+ τ factual rolling technological parameter ψ constantly k ' yMr.The shelf number of front two rolling mills of subscript k ' representative, k '=1,2; Subscript r represents t 0to t 0+ τ fall kinds of actual rolling technological parameter collected is constantly according to the numbering of time order and function order, r=1, and 2 ..., m, m is illustrated in t 0to t 0+ τ fthe number of the rolling technological parameter of collected particular type constantly,
Figure BDA00001644947300061
C14) consider at t 0to t 0+ τ fthere is large jump because of some accidentalia in the actual rolling technological parameter of only a few that may occur constantly, but this parameter can not reflect real as-rolled condition, therefore must reject, and introduces parameter for this reason and reject process variables array { α k ' r, and by { α k ' rwith following formula, represent:
{ &alpha; k &prime; r } = { &Sigma; y = 1 4 ( | &psi; k &prime; yMr - 1 m &Sigma; r = 1 m &psi; k &prime; yMr | 1 m &Sigma; r = 1 m &psi; k &prime; yMr ) + | h 0 Mr - 1 m &Sigma; r = 1 m h 0 Mr | 1 m &Sigma; r = 1 m h 0 Mr y = 1,2,3,4 } , α k ' rrepresent that the k ' rolling mill is at t 0+ r τ srejecting variable constantly, the shelf number of front two rolling mills of subscript k ' representative, k '=1,2;
C15) by { α k ' raccording to the size of numerical value according to being descending sort from big to small, and definition sequence working array { n ll=1,2 ..., m} makes this array meet inequality
Figure BDA00001644947300063
n wherein lrepresentative { α k ' rcorresponding according to the numbering of time sequencing sequence according to l rejecting variable after numerical value descending sort, l represents { α k ' raccording to the big or small sequence number after numerical value descending sort;
C16) given rejecting proportionality coefficient χ, usually, χ=0 ~ 50%;
C17) get rid of the interference of accidentalia, at t 0to t 0+ τ fin the actual rolling technological parameter gathering constantly, weed out the data of χ ratio, then remaining (1-χ) data are averaged, obtain t 0to t 0+ τ fconstantly interior for feeding back the feature rolling technological parameter ψ of calculating k ' yM(y=1,2,3,4), and feature supplied materials thickness h 0M, its expression formula is respectively
Figure BDA00001644947300064
Figure BDA00001644947300065
m wherein 1represent t 0to t 0+ τ fthe number of the actual rolling parameter of rejecting constantly, m 1=int (χ m), m 2represent t 0to t 0+ τ fconstantly, retain the number of the actual rolling parameter of not rejecting, m 2=m-m 1;
C18) according to t 0to t 0+ τ ffeature rolling technological parameter constantly, respectively the value of feedback S of bending roller force in the interior 1# of definite current feedback cycle, 2# frame 1M, S 2M, the roller amount of inclining value of feedback η 1M, η 2M, the value of feedback δ of roll shifting amount 1M, δ 2M, can adopt the following step of being carried out by computer (fundamental block diagram is shown in accompanying drawing 4) to complete:
C18-1) objective definition function initial value F 0, and make F 0=10 10;
C18-2) according to target proportion crown curve coefficient a εcalculate 2# milling train and export corresponding target proportion convexity Distribution Value
Figure BDA00001644947300071
fundamental formular is: x in formula ifor the relative coordinate of used unit of computational process,
Figure BDA00001644947300073
i is band steel bar unit number;
C18-3) make δ k ' M=(L wk '-B-150)/2;
C18-4) make the roller amount η that inclines k '=0;
C18-5) definition bending roller force procedure parameter θ, and make θ=0;
C18-6) make bending roller force in 1# frame
Figure BDA00001644947300074
bending roller force in 2# frame S 2 = - S 2 max - + &theta; &phi; 2 ( S 2 max - + S 2 max + ) ;
C18-7) utilizing the described outlet thickness of slab of document [1] distribution machine reason model to calculate feature rolling technological parameter is ψ k ' yM(y=1,2,3,4), feature supplied materials thickness are h 0M, and roll shifting amount is δ k ' M, the roller amount of inclining is η k ', bending roller force is S k ', average crown curve function is f in supplied materials interval mj(X) time, width is the band that B, steel grade code are SG, at the Outlet Section thickness distribution value h of the 2nd rolling mill miwith band forward pull cross direction profiles value σ 1Mi, i is band in horizontal bar unit number, i=1, and 2 ..., n, n is the total bar unit number of band;
C18-8) calculate the maximum residual stress σ of the 2nd rolling mill outlet band mmax, basic calculating model is: &sigma; M max = max { | &sigma; 1 Mi - &Sigma; i = 1 n &sigma; 1 Mi n | } ;
C18-9) utilize the described method of document [2] to calculate the critical residual stress σ of current the 2nd rolling mill outlet band generation plate shape unstability mcr, basic calculating model is: wherein, k crfor the limit stress coefficient of band generation plate shape unstability, the elastic modelling quantity that E, υ are belt steel material and Poisson's ratio, B is strip width, h *it is the 2nd rolling mill exit thickness of strip;
C18-10) calculate the actual ratio convexity Distribution Value of the 2nd rolling mill outlet band
C18-11) calculate the deviation of actual ratio convexity Distribution Value and target proportion convexity Distribution Value, its function expression is F = &Sigma; i = 1 n ( C i * - C Mi ) 2 ;
C18-12) judgement inequality group &sigma; M max < &sigma; Mcr F < F 0 Set up? if inequality is set up, make F 0=F, S 1M=S 1, S 2M=S 2, θ=θ+1, proceed to step c18-13; Otherwise, make θ=θ+1, proceed to step c18-13;
C18-13) judgement inequality group &theta; &phi; 1 &le; 1.0 &theta; &phi; 2 &le; 1.0 Set up? if inequality is set up, proceed to step c18-6, repeating step c18-6 to c18-13, until no longer set up with upper inequality group; Otherwise, proceed to step c18-14;
C18-14) define the roller amount procedure parameter λ that inclines, and make λ=0, F 0=10 10;
C18-15) make roller declination amount in 1# frame
Figure BDA00001644947300086
roller declination amount in 2# frame
Figure BDA00001644947300087
C18-16) utilizing the described outlet thickness of slab of document [1] distribution machine reason model to calculate feature rolling technological parameter is ψ k ' yM(y=1,2,3,4), feature supplied materials thickness are h 0M, and roll shifting amount is δ k ' M, bending roller force is S k ' M, the roller amount of inclining η k ', average crown curve function is f in supplied materials interval mj(X) time, width is the band that B, steel grade code are SG, at the Outlet Section thickness distribution value C of the 2nd rolling mill miwith band forward pull cross direction profiles value σ 1Mi;
C18-17) calculate the maximum residual stress σ of the 2nd rolling mill outlet band mmax, basic calculating model is: &sigma; M max = max { &sigma; 1 Mi - &Sigma; i = 1 n &sigma; 1 Mi n } ;
C18-18) utilize the critical residual stress σ of current the 2nd rolling mill outlet band generation plate shape unstability of the described method meter of document [2] mcr, basic calculating model is:
Figure BDA00001644947300092
C18-19) calculate the actual ratio convexity Distribution Value of the 2nd rolling mill outlet band
Figure BDA00001644947300093
C18-20) calculate target cross section characteristic (being that Target Board convexity distributes) and the deviation that actual convexity distributes, its function expression is
Figure BDA00001644947300094
C18-21) judgement inequality &sigma; M max < &sigma; Mcr F < F 0 Set up? if inequality is set up, make F 0=F, η 1M1, η 2M2, λ=λ+1, proceed to step c18-22; Otherwise, make λ=λ+1, proceed to step c18-22;
C18-22) judgement inequality group
Figure BDA00001644947300096
set up? if inequality is set up, proceed to step c18-15, repeating step c18-15 to c18-22, until no longer set up with upper inequality group; Otherwise, proceed to step c19;
C19) by t 0to t 0+ τ fthe value of feedback δ of roll shifting amount constantly k ' M, bending roller force value of feedback S k ' M, the roller amount of inclining value of feedback η k ' Msend to a level system of tandem mills, respectively roll shifting amount, the bending roller force of current 1#, 2# milling train, the roller amount of inclining are reset;
C20) do you judge that inequality ξ < 0 sets up? if inequality is set up, finish cross section characteristic feedback, the cross section characteristic reponse system of virtual profile gauge quits work, and proceeds to c6; If inequality is false, make M=M+1, t 0=t 0+ τ f, proceed to step c8, until inequality ξ < 0 sets up.
Advantage of the present invention is: realized the accurate calculating to cold rolling finished product cross section characteristic, and carried out feedback regulation, thereby farthest guaranteed the thickness of slab quality of cold rolling finished product.
Accompanying drawing explanation
Accompanying drawing 1 is the function layout drawing of virtual profile gauge;
Accompanying drawing 2 is FB(flow block)s of Outlet Section shape Presentation Function;
Accompanying drawing 3 is FB(flow block)s that cross section characteristic closed loop feedback is controlled function;
Accompanying drawing 4 is FB(flow block)s of bending roller force in current feedback cycle, the roller amount of inclining, roll shifting amount;
Accompanying drawing 5 is that in embodiment 1, supplied materials twists in the average crown curve function in certain length interval;
Accompanying drawing 6 is in embodiment 1, to calculate a certain section thickness Distribution Value of gained and matched curve function;
Accompanying drawing 7 is software interfaces that in embodiment 1, section crown curve dynamically shows;
Accompanying drawing 8 is the target proportion convexity Distribution Values that calculate gained in embodiment 1;
Accompanying drawing 9 is Outlet Section thickness distribution value and the matched curves thereof that regulate the second rolling mill after bending roller force in embodiment 1;
Accompanying drawing 10 is in embodiment 1, to regulate the actual ratio convexity of the second rolling mill Outlet Section after bending roller force to distribute;
Accompanying drawing 11 is Outlet Section thickness distribution value and the matched curves thereof that regulate the second rolling mill after the roller amount of inclining in embodiment 1;
Accompanying drawing 12 is in embodiment 1, to regulate the actual ratio convexity of the second rolling mill Outlet Section after the roller amount of inclining to distribute;
Accompanying drawing 13 is that in embodiment 2, supplied materials twists in the average crown curve function in certain length interval;
Accompanying drawing 14 is in embodiment 2, to calculate a certain section thickness Distribution Value of gained and matched curve function;
Accompanying drawing 15 is software interfaces that in embodiment 2, section crown curve dynamically shows;
Accompanying drawing 16 is the target proportion convexity Distribution Values that calculate gained in embodiment 2;
Accompanying drawing 17 is Outlet Section thickness distribution value and the matched curves thereof that regulate the second rolling mill after bending roller force in embodiment 2;
Accompanying drawing 18 is in embodiment 2, to regulate the actual ratio convexity of the second rolling mill Outlet Section after bending roller force to distribute;
Accompanying drawing 19 is Outlet Section thickness distribution value and the matched curves thereof that regulate the second rolling mill after the roller amount of inclining in embodiment 2;
Accompanying drawing 20 is in embodiment 2, to regulate the actual ratio convexity of the second rolling mill Outlet Section after the roller amount of inclining to distribute.
The specific embodiment
Embodiment 1
In the present embodiment, utilize accompanying drawing, in conjunction with the service condition of domestic certain steel mill's five Stands Cold Tandem Mill group to this virtual profile gauge, the virtual profile gauge establishing method to four roller tandem mills of the present invention based on mechanism model describes.Concrete steps are as follows:
(1) basic equipment parameter collection.The basic equipment parameter of collecting each rolling mill of tandem mills, partial parameters value is as shown in table 1;
Each rolling mill basic equipment parameter of certain tandem mills of table 1
The first frame The second frame The 3rd frame The 4th frame The 5th frame
Working roll barrel length (mm) 1220 1220 1220 1220 1220
Backing roll barrel length (mm) 1220 1220 1220 1220 1220
Working roll bending cylinder distance (mm) 2100 2100 2100 2100 2100
Housing screw centre-to-centre spacing (mm) 2100 2100 2100 2100 2100
The roller footpath (mm) of working roll 517.5 471.1 474 524.1 510.8
The roller footpath (mm) of backing roll 1290.3 1353.3 1325.3 1245.8 1268.7
Maximum positive bending roller force (kN/chock) 350 350 350 350 350
Maximum negative bending roller force (kN/chock) 280 280 280 280 280
The maximum roller amount (mm) of just inclining 0.35 0.35 0.35 0.35 0.35
The maximum negative roller amount (mm) of inclining 0.35 0.35 0.35 0.35 0.35
(2) realization of Outlet Section shape Presentation Function.Mainly comprise the following step of being carried out by computer:
First, in step 1, collect the sampling period τ of current tandem mills data collecting system s=0.02s;
Subsequently, in step 2, provide the period tau that virtual profile gauge section configuration dynamically shows x=0.5s;
Subsequently, in step 3, provide the interval parameter L to supplied materials sectional sampling 0=50m;
Subsequently, in step 4, definition length of interval parameter j, supplied materials coil of strip length parameter L j;
Subsequently, in step 5, definition as-rolled condition parameter ξ;
Subsequently, in step 6, given initial steel reel number COILNO1=0, prepares convexity and dynamically shows;
Subsequently, in step 7, from the three-level system of tandem mills, collect supplied materials coil of strip essential information, comprise reel number COILNO=1449592000, steel grade code SG=MR2T3, strip width B=816mm, supplied materials thickness h 0=1.82mm, last rolling mill exit thickness h 5=0.211mm;
Subsequently, in step 8, judge inequality | does COILNO1-COILNO|>0 set up? obviously, inequality | 1449592000-0|>0 sets up, and makes j=0, L j=0, and read the cross section characteristic data that supplied materials steel reel number is COILNO, proceed to step 9; If inequality is false, directly proceed to step 9;
Subsequently, in step 9, collect current time the first frame entrance and roll band steel length L *=4.33m;
Subsequently, in step 10, judgement inequality L *>=L jset up? obviously, inequality 4.33>=0 is set up, and proceeds to step 11; If be false, proceed to step 12;
Subsequently, in step 11, extract and calculate supplied materials and twist in the average crown curve function f in [0,50] m interval 0(X)=1.82949+2.54854E-6*X+1.0832E-7*X 2+ 5.89841E-12*X 3-2.56288E-12*X 4(function curve is as shown in Figure 5), wherein X is illustrated in strip width direction and sets up one-dimensional coordinate system, when the center of width is the origin of coordinates, and the corresponding numerical value of strip width direction each point in coordinate system;
Subsequently, in step 12, by tandem mills data collecting system, collect the rolling technological parameter actual value ψ of current each frame ky, concrete numerical value is as shown in table 2;
The rolling technological parameter actual value of each frame of table 2
The first frame The second frame The 3rd frame The 4th frame The 5th frame
Draught pressure (kN) 7462.7 4816.1 5395.979 5115.58 5877.64
Forward pull (kN) 72.54 174.82 97.35999 62.94 40.7
Backward pull (kN) 174.82 97.35999 62.94 40.7 11.75
Reduction ratio (%) 30.99614 39.55313 33.9885 38.1 31.9458
Roll shifting amount (mm) 77 77 77 77 77
Bending roller force (kN) 104.4 230.2 196.96 225.71 240.58
The roller amount of inclining (kN/chock) 0 0 0 0 0
Subsequently, in step 13, utilizing the described outlet thickness of slab of document [1] distribution machine reason model to calculate actual rolling technological parameter is ψ kyand in supplied materials interval, average crown curve function is f j=(X) time, width is the Outlet Section thickness distribution value h of B, the steel grade code band that is SG i=0.2065 0.2084 0.2092 0.20960.2101 0.2103 0.2106 0.2108 0.2109 0.2110 0.2111 0.21120.2112 0.2113 0.2113 0.2114 0.2114 0.2115 0.2115 0.21140.2114 0.2113 0.2113 0.2112 0.2111 0.2110 0.2109 0.21080.2106 0.2103 0.2101 0.2096 0.2092 0.2084 0.2064}(curve maps are as shown in Figure 6);
Subsequently, in step 14, utilize and calculate a mouthful section thickness Distribution Value h ithe current section crown curve of matching function f c(X)=0.21128+1.01431E-7*X+3.95354E-9*X 2-1.04855E-12*X 3-1.77341E-13*X 4(function curve as shown in Figure 6);
Subsequently, in step 15, calculate current section convex value CROWN and wedge shape value WEDGE, convex value is CROWN = f c ( 0 ) - f c ( 368 ) + f c ( - 368 ) 2 = 3.8 &mu;m , Wedge shape amount computation model is WEDGE = | f c ( 368 ) - f c ( - 368 ) 2 | = 0 ;
Subsequently, in step 16, utilize the dynamic display of functional of visual software, by current section crown curve function f c(X) dynamically show, show convex value CROWN, wedge shape value WEDGE simultaneously, and corresponding coil of strip information etc., software interface is as shown in Figure 7;
Do you subsequently, in step 17, judge that inequality ξ < 0 sets up? if inequality is set up, finish cross section characteristic and show, virtual plate profile instrument quits work; If inequality is false, COILNO1=COILNO, proceeds to step 7, until inequality ξ < 0 sets up.
(3) cross section characteristic closed loop feedback is controlled the realization of function.Mainly comprise the following step of being carried out by computer:
First, in step 1, collect the target proportion crown curve coefficient a of 2# milling train outlet 0=1.00374, a 1=0, a 2=-3.2189E-5, a 3=0, a 4=-2.44979E-8;
Subsequently, in step 2, in definite cross section characteristic feedback procedure, in 1# and 2# frame, the single step minimum of bending roller force is adjusted coefficient φ respectively 1=0.05, φ 2=0.048;
Subsequently, in step 3, the minimum coefficient of adjusting of single step of the roller amount of inclining in 1# and 2# frame in definite cross section characteristic feedback procedure respectively
Figure BDA00001644947300141
Figure BDA00001644947300142
Subsequently, in step 4, determine the period tau of cross section characteristic feedback f=2s;
Subsequently, in step 5, definition rolling is procedure parameter t constantly 0, virtual profile gauge cross section characteristic Times of Feedback procedure parameter M;
Subsequently, in step 6, receive operational order, judge whether to start rolling? if started rolling, record current standard Beijing time t, and make t 0=t, proceeds to step 7; If do not start rolling, enter wait state;
Subsequently, in step 7, given initial steel reel number COILNO1=0;
First, in step 8, from the three-level system of tandem mills, collect supplied materials coil of strip essential information, reel number COILNO=1449592000, steel grade code SG=MR2T3, strip width B=816mm, supplied materials thickness h 0=1.82mm, last rolling mill exit thickness h 5=0.211mm etc.;
Subsequently, in step 9, judge inequality | does COILNO1-COILNO|>0 set up? obviously, inequality | 1449592000-0|>0 sets up, and makes j=0, L j=0, and read the cross section characteristic data that supplied materials steel reel number is COILNO, proceed to step 10; If inequality is false, directly proceed to step 10;
Subsequently, in step 10, collect current time the first frame entrance and roll band steel length L *=4.33m;
Subsequently, in step 11, judgement inequality L *>=L jset up? obviously, inequality 4.33>=0 is set up, and proceeds to step 12; If be false, proceed to step 13;
Subsequently, in step 12, extract and calculate supplied materials and twist in the average crown curve function f in [0,50] m interval 0(X)=1.82949+2.54854E-6*X+1.0832E-7*X 2+ 5.89841E-12*X 3-2.56288E-12*X 4(function curve is as shown in Figure 5), wherein X is illustrated in strip width direction and sets up one-dimensional coordinate system, when the center of width is the origin of coordinates, and the corresponding numerical value of strip width direction each point in coordinate system;
Subsequently, in step 13, by milling train data collecting system, collect 1#, 2# milling train from t 0to t 0+ τ fall actual rolling technological parameter ψ constantly k ' yMr, subscript r represents t 0to t 0+ τ fall kinds of actual rolling technological parameter collected is constantly according to the numbering of time order and function order, r=1, and 2 ..., m, m is illustrated in t 0to t 0+ τ fthe number of the rolling technological parameter of collected particular type constantly,
Figure BDA00001644947300151
Subsequently, in step 14, introduce parameter and reject process variables array { α k ' r, { &alpha; k &prime; r } = { &Sigma; y = 1 4 ( | &psi; k &prime; yMr - 1 m &Sigma; r = 1 m &psi; k &prime; yMr | 1 m &Sigma; r = 1 m &psi; k &prime; yMr ) + | h 0 Mr - 1 m &Sigma; r = 1 m h 0 Mr | 1 m &Sigma; r = 1 m h 0 Mr y = 1,2,3,4 } ;
Subsequently, in step 15, by { α k ' raccording to the size of numerical value according to being descending sort from big to small, and definition sequence working array { n ll=1,2 ..., m} makes this array meet inequality &alpha; k &prime; n 1 &GreaterEqual; &alpha; k &prime; n 2 &GreaterEqual; &CenterDot; &CenterDot; &CenterDot; &GreaterEqual; &alpha; k &prime; n l - 1 &GreaterEqual; &alpha; k &prime; n l &GreaterEqual; &alpha; k &prime; n l + 1 &GreaterEqual; &CenterDot; &CenterDot; &CenterDot; &GreaterEqual; &alpha; k &prime; n m ;
Subsequently, in step 16, given rejecting proportionality coefficient χ=25%;
Subsequently, in step 17, get rid of the interference of accidentalia, at t 0to t 0+ τ fin the actual rolling technological parameter gathering constantly, weed out 25% data, then 75% remaining data are averaged, obtain t 0to t 0+ τ fconstantly interior for feeding back the feature rolling technological parameter of calculating
Figure BDA00001644947300161
(y=1,2,3,4), and feature supplied materials thickness
Figure BDA00001644947300162
m wherein 1represent t 0to t 0+ τ fthe number of the actual rolling parameter of rejecting constantly, m 1=int (χ m)=31, m 2represent t 0to t 0+ τ fconstantly, retain the number of the actual rolling parameter of not rejecting, m 2=m-m 1=94;
Subsequently, in step 18, according to t 0to t 0+ τ ffeature rolling technological parameter constantly, respectively the value of feedback S of bending roller force in the interior 1# of definite current feedback cycle, 2# frame 1M, S 2M, the roller amount of inclining value of feedback η 1M, η 2M, the value of feedback δ of roll shifting amount 1M, δ 2M, adopt the following step of being carried out by computer:
First, in step 18-1, objective definition function initial value F 0, and make F 0=10 10;
Subsequently, in step 18-2, according to target proportion crown curve coefficient a εcalculate the target proportion convexity Distribution Value of 2# milling train outlet C i * = &Sigma; &epsiv; = 0 3 a &epsiv; x i &epsiv; = 0.992380952 0.993809524 0.995238095 0.996666667 0.997619048 0.998571429 0.99952381 1.00047619 1.000952381 1.001428571 1.001904762 1.002380952 1.00285714 1.003333333 1.003333333 1.003809524 1.003809524 1.003809524 1.00380952 1.003809524 1.003333333 1.003333333 1.002857143 1.002380952 1.00190476 1.001428571 1.000952381 1.00047619 0.99952381 0.998571429 0.99761905 0.996666667 0.995238095 0.993809524 0.992380952 (curve is as shown in Figure 8), x i = - n 2 + i - 0.5 , i = 0,1,2 , . . . , 35 ;
Subsequently, in step 18-3, make δ k ' M=(L k '-B-150)/2=77mm;
Subsequently, in step 18-4, make the roller amount η that inclines k '=0;
Subsequently, in step 18-5, definition bending roller force procedure parameter θ, and make θ=0;
Subsequently, in step 18-6, make bending roller force in 1# frame
Figure BDA000016449473001611
kN, bending roller force in 2# frame S 2 = - S 2 max - + &theta; &phi; 2 ( S 2 max - + S 2 max + ) = - 280 kN ;
Subsequently, in step 18-7, the Outlet Section thickness distribution value h that utilizes the described outlet thickness of slab of document [1] distribution machine reason model to calculate at the 2nd rolling mill mi={ 0.7322 0.7395 0.7456 0.7507 0.75470.758 0.7606 0.7626 0.764 0.7651 0.7659 0.7664 0.7668 0.7670.7671 0.7672 0.7672 0.7672 0.7672 0.7672 0.7671 0.767 0.76680.7664 0.7659 0.7651 0.764 0.7626 0.7606 0.758 0.7547 0.75070.7456 0.7395 0.7322} and band forward pull cross direction profiles value σ 1Mi=117.4 126.4135.5 144.3 152.9 161.1 168.7 175.9 182.5 188.4 193.7 198.3 202.2205.5 208 209.8 210.9 211.2 210.9 209.8 208 205.5 202.2 198.3193.7 188.4 182.5 175.9 168.7 161.1 152.9 144.3 135.5 126.4 117.4}(curves are as shown in Figure 9);
Subsequently, in step 18-8, calculate the maximum residual stress of the second rolling mill outlet band &sigma; M max = max { | &sigma; 1 Mi - &Sigma; i = 1 n &sigma; 1 Mi n | } = 59 MPa ;
Subsequently, in step 18-9, utilize the described method of document [2] to calculate the critical residual stress of current the second rolling mill outlet band generation plate shape unstability
Figure BDA00001644947300172
Subsequently, in step 18-10, calculate the actual ratio convexity Distribution Value of the second rolling mill outlet band C Mi = h Mi ( &Sigma; i = 1 n h Mi ) / n = 0.964567111 0.9741838 0.982219664 0.988938173 0.994207591 0.998554862 1.001979984 1.004614694 1.00645899 1.00790808 1.008961964 1.009620641 1.010147583 1.010411054 1.01054279 1.010674525 1.010674525 1.010674525 1.010674525 1.010674525 1.01054279 1.010411054 1.010147583 1.009620641 1.008961964 1.00790808 1.00645899 1.004614694 1.001979984 0.998554862 0.994207591 0.988938173 0.982219664 0.9741838 0.964567111 } (curve as shown in Figure 10);
Subsequently, in step 18-11, calculate the deviation of actual ratio convexity Distribution Value and target proportion convexity Distribution Value, its function expression is
Figure BDA00001644947300181
Subsequently, in step 18-12, judgement inequality group &sigma; M max < &sigma; Mcr F < F 0 Set up? obviously, inequality group 59 < 21 0.003741 < 10 10 Be false, make θ=θ+1=1, proceed to step 18-13; If but set up, make F 0=F, S 1M=S 1, S 2M=S 2, θ=θ+1, proceed to step 18-13;
Subsequently, in step 18-13, judgement inequality group &theta; &phi; 1 &le; 1.0 &theta; &phi; 2 &le; 1.0 Set up? obvious inequality group 0 &le; 1.0 0 &le; 1.0 Set up, proceed to step 18-6, repeating step 18-6 to 18-13, until no longer set up with upper inequality group; When inequality group is no longer set up, proceed to step 18-14;
Subsequently, in step 18-14, define the roller amount procedure parameter λ that inclines, and make λ=0, F 0=10 10;
Subsequently, in step 18-15, make roller declination amount in 1# frame
Figure BDA00001644947300186
roller declination amount in 2# frame
Figure BDA00001644947300187
Subsequently, in step 18-16, utilize the described outlet thickness of slab of document [1] distribution machine reason model to calculate the Outlet Section thickness distribution value C of the second rolling mill mi={ 0.7342 0.7416 0.7461 0.7513 0.75480.7582 0.7601 0.7626 0.7627 0.7646 0.7642 0.7636 0.7645 0.7650.7661 0.7651 0.7652 0.7659 0.7652 0.7661 0.7651 0.7653 0.76450.764 0.7642 0.7636 0.7632 0.7616 0.7601 0.7572 0.7558 0.75280.7497 0.74460.7392} and band forward pull cross direction profiles value σ 1Mi=141.8 146.1 150.1154 157.7 161.2 164.5 167.5 170.4 173.1 175.5 177.8 179.9 181.9183.8 185.6 187.3 189 190.4 191.3 191.8 191.9 191.5 190.8 189.8188.5 186.9 185.1 183 180.7 178.3 175.7 173.1 170.3 167.5}(curves are as shown in Figure 11);
Subsequently, in step 18-17, calculate the maximum residual stress of the second rolling mill outlet band &sigma; M max = max { | &sigma; 1 Mi - &Sigma; i = 1 n &sigma; 1 Mi n | } = 34.6 MPa ;
Subsequently, in step 18-18, utilize the described method of document [2] to calculate the critical residual stress of current the second rolling mill outlet band generation plate shape unstability
Figure BDA00001644947300192
Subsequently, in step 18-19, calculate the actual ratio convexity Distribution Value of the second rolling mill outlet band
C Mi = h Mi ( &Sigma; i = 1 n h Mi ) / n = 0.96758039 0.97733263 0.983263047 0.990115973 0.994728519 0.999209278 1.001713231 1.005007907 1.005139694 1.007643648 1.0071165 1.006325778 1.007511861 1.008170796 1.009620453 1.008302583 1.00843437 1.009356879 1.00843437 1.009620453 1.008302583 1.008566157 1.007511861 1.006852926 1.0071165 1.006325778 1.005798629 1.003690037 1.001713231 0.997891407 0.996046389 0.992092778 0.98800738 0.981289241 0.974169742 } (curve as shown in Figure 12);
Subsequently, in step 18-20, calculate target cross section characteristic (being that Target Board convexity distributes) and the deviation that actual convexity distributes, its function expression is
Figure BDA000016449473001911
Subsequently, in step 18-21, judge inequality &sigma; M max < &sigma; Mcr F < F 0 Set up? obviously, inequality group 34.6 < 21 0.002166 < 10 10 Be false, make λ=λ+1=1, proceed to step 18-22; If but set up, make F 0=F, η 1M1, η 2M2, λ=λ+1, proceed to step 18-22; ;
Subsequently, in step 18-22, judgement inequality group
Figure BDA00001644947300201
set up? obvious inequality group 0 &le; 1.0 0 &le; 1.0 Set up, proceed to step 18-15, repeating step c18-15 to c18-22, until no longer set up with upper inequality group; When inequality group is no longer set up, proceed to step 18-19;
Subsequently, in step 19, by t 0to t 0+ τ fthe value of feedback δ of roll shifting amount constantly k ' Mthe value of feedback S of=77mm, bending roller force 1M=115MPa, S 2Mthe value of feedback η of=134.4MPa, the roller amount of inclining 1M=0mm, η 2M=0mm sends to a level system of tandem mills, respectively roll shifting amount, the bending roller force of current 1#, 2# milling train, the roller amount of inclining is reset;
Do you subsequently, in step 20, judge that inequality ξ < 0 sets up? if inequality is set up, finish cross section characteristic feedback, the cross section characteristic reponse system of virtual profile gauge quits work, and proceeds to step 6; If inequality is false, make M=M+1, t 0=t 0+ τ f, proceed to step 8, until inequality ξ < 0 sets up.
Embodiment 2
In the present embodiment, utilize accompanying drawing, in conjunction with the service condition of domestic certain steel mill's five Stands Cold Tandem Mill group to this virtual profile gauge, the virtual profile gauge method for designing to four roller tandem mills of the present invention based on mechanism model describes.Concrete steps are as follows:
(1) basic equipment parameter collection.The basic equipment parameter of collecting each rolling mill of tandem mills, partial parameters value is as shown in table 1;
Each rolling mill basic equipment parameter of certain tandem mills of table 1
The first frame The second frame The 3rd frame The 4th frame The 5th frame
Working roll barrel length (mm) 1220 1220 1220 1220 1220
Backing roll barrel length (mm) 1220 1220 1220 1220 1220
Working roll bending cylinder distance (mm) 2100 2100 2100 2100 2100
Housing screw centre-to-centre spacing (mm) 2100 2100 2100 2100 2100
The roller footpath (mm) of working roll 517.5 471.1 474 524.1 510.8
The roller footpath (mm) of backing roll 1290.3 1353.3 1325.3 1245.8 1268.7
Maximum positive bending roller force (kN/chock) 350 350 350 350 350
Maximum negative bending roller force (kN/chock) 280 280 280 280 280
The maximum roller amount (mm) of just inclining 0.35 0.35 0.35 0.35 0.35
The maximum negative roller amount (mm) of inclining 0.35 0.35 0.35 0.35 0.35
(2) realization of Outlet Section shape Presentation Function.Mainly comprise the following step of being carried out by computer:
First, in step 1, collect the sampling period τ of current tandem mills data collecting system s=0.02s;
Subsequently, in step 2, provide the period tau that virtual profile gauge section configuration dynamically shows x=0.5s;
Subsequently, in step 3, provide the interval parameter L to supplied materials sectional sampling 0=50m;
Subsequently, in step 4, definition length of interval parameter j, supplied materials coil of strip length parameter L j;
Subsequently, in step 5, definition as-rolled condition parameter ξ;
Subsequently, in step 6, given initial steel reel number COILNO1=0, prepares convexity and dynamically shows;
Subsequently, in step 7, from the three-level system of tandem mills, collect supplied materials coil of strip essential information, comprise reel number COILNO=1448630000, steel grade code SG=MRT5CA, strip width B=924mm, supplied materials thickness h 0=2.01mm, last rolling mill exit thickness h 5=0.238mm;
Subsequently, in step 8, judge inequality | does COILNO1-COILNO|>0 set up? obviously, inequality | 1448630000-0|>0 sets up, and makes j=0, L j=0, and read the cross section characteristic data that supplied materials steel reel number is COILNO, proceed to step 9; If inequality is false, directly proceed to step 9;
Subsequently, in step 9, collect current time the first frame entrance and roll band steel length L *=5.28m;
Subsequently, in step 10, judgement inequality L *>=L jset up? obviously, inequality 5.28>=0 is set up, and proceeds to step 11; If be false, proceed to step 12;
Subsequently, in step 11, extract and calculate supplied materials and twist in the average crown curve function f in [0,50] m interval 0(X)=2.02397+4.42582E-6*X-1.08067E-8*X 2-1.14068E-11*X 3-1.25427E-12*X 4(function curve is as shown in figure 13), wherein X is illustrated in strip width direction and sets up one-dimensional coordinate system, when the center of width is the origin of coordinates, and the corresponding numerical value of strip width direction each point in coordinate system;
Subsequently, in step 12, by tandem mills data collecting system, collect the rolling technological parameter actual value ψ of current each frame ky, concrete numerical value is as shown in table 2;
The rolling technological parameter actual value of each frame of table 2
The first frame The second frame The 3rd frame The 4th frame The 5th frame
Draught pressure (kN) 9126.96 8476.561 9718.261 7402.78 7420.12
Forward pull (kN) 91.11 202.64 135.36 81.58001 55.44
Backward pull (kN) 202.64 135.36 81.58001 55.44 14.91
Reduction ratio (%) 26.51067 42.6964 40.24942 33.76749 28.95045
Roll shifting amount (mm) 77 77 77 77 77
Bending roller force (kN) 129.07 192.77 123.02 184.61 230.14
The roller amount of inclining (kN/chock) 0 0 0 0 0
Subsequently, in step 13, utilizing the described outlet thickness of slab of document [1] distribution machine reason model to calculate actual rolling technological parameter is ψ kyand in supplied materials interval, average crown curve function is f j=(X) time, width is the Outlet Section thickness distribution value h of B, the steel grade code band that is SG i=0.2340 0.2352 0.2361 0.23680.2375 0.2380 0.2385 0.2386 0.2389 0.2392 0.2395 0.23960.2396 0.2396 0.2396 0.2398 0.2396 0.2397 0.2397 0.23970.2395 0.2395 0.2392 0.2390 0.2385 0.2380 0.2375 0.23680.2360 0.2351 0.2342}(curves are as shown in figure 14);
Subsequently, in step 14, utilize and calculate a mouthful section thickness Distribution Value h ithe current section crown curve of matching function f c(X)=0.23972+3.94491E-7*X-1.09881E-8*X 2-2.12456E-12*X 3-7.3116E-14*X 4(function curve as shown in Figure 14);
Subsequently, in step 15, calculate current section convex value CROWN and wedge shape value WEDGE, convex value is CROWN = f c ( 0 ) - f c ( 422 ) + f c ( - 422 ) 2 = 4.3 &mu;m , Wedge shape amount computation model is WEDGE = | f c ( 422 ) - f c ( - 422 ) 2 | = 0 ;
Subsequently, in step 16, utilize the dynamic display of functional of visual software, by current section crown curve function f c(X) dynamically show, show convex value CROWN, wedge shape value WEDGE simultaneously, and corresponding coil of strip information etc., software interface is as shown in figure 15;
Do you subsequently, in step 17, judge that inequality ξ < 0 sets up? if inequality is set up, finish cross section characteristic and show, virtual plate profile instrument quits work; If inequality is false, COILNO1=COILNO, proceeds to step 7, until inequality ξ < 0 sets up.
(3) cross section characteristic closed loop feedback is controlled the realization of function.Mainly comprise the following step of being carried out by computer:
First, in step 1, collect the target proportion crown curve coefficient a of the second milling train outlet 0=0.849, a 1=0, a 2=-9.04244E-6, a 3=0, a 4=-1.57174E-7;
Subsequently, in step 2, in definite cross section characteristic feedback procedure, in the first and second frames, the single step minimum of bending roller force is adjusted coefficient φ respectively 1=0.05, φ 2=0.04;
Subsequently, in step 3, the minimum coefficient of adjusting of single step of the roller amount of inclining in the first and second frames in definite cross section characteristic feedback procedure respectively
Figure BDA00001644947300232
Figure BDA00001644947300233
Subsequently, in step 4, determine the period tau of cross section characteristic feedback f=2s;
Subsequently, in step 5, definition rolling is procedure parameter t constantly 0, virtual profile gauge cross section characteristic Times of Feedback procedure parameter M;
Subsequently, in step 6, receive operational order, judge whether to start rolling? if started rolling, record current standard Beijing time t, and make t 0=t, proceeds to step 7; If do not start rolling, enter wait state;
Subsequently, in step 7, given initial steel reel number COILNO1=0;
First, in step 8, from the three-level system of tandem mills, collect supplied materials coil of strip essential information, reel number COILNO=1448630000, steel grade code SG=MRT5CA, strip width B=924mm, supplied materials thickness h 0=2.01mm, last rolling mill exit thickness h 5=0.238mm etc.;
Subsequently, in step 9, judge inequality | does COILNO1-COILNO|>0 set up? obviously, inequality | 1448630000-0|>0 sets up, and makes j=0, L j=0, and read the cross section characteristic data that supplied materials steel reel number is COILNO, proceed to step 10; If inequality is false, directly proceed to step 10;
Subsequently, in step 10, collect current time the first frame entrance and roll band steel length L *=5.28m;
Subsequently, in step 11, judgement inequality L *>=L jset up? obviously, inequality 5.28>=0 is set up, and proceeds to step 12; If be false, proceed to step 13;
Subsequently, in step 12, extract and calculate supplied materials and twist in the average crown curve function f in [0,50] m interval 0(X)=2.02397+4.42582E-6*X-1.08067E-8*X 2-1.14068E-11*X 3-1.25427E-12*X 4(curve is as shown in Figure 13), wherein X is illustrated in strip width direction and sets up one-dimensional coordinate system, when the center of width is the origin of coordinates, the corresponding numerical value of strip width direction each point in coordinate system;
Subsequently, in step 13, by milling train data collecting system, collect first, second milling train from t 0to t 0+ τ fall actual rolling technological parameter ψ constantly k ' yMr, subscript r represents t 0to t 0+ τ fall kinds of actual rolling technological parameter collected is constantly according to the numbering of time order and function order, r=1, and 2 ..., m, m is illustrated in t 0to t 0+ τ fthe number of the rolling technological parameter of collected particular type constantly,
Figure BDA00001644947300241
Subsequently, in step 14, introduce parameter and reject process variables array { α k ' r, { &alpha; k &prime; r } = { &Sigma; y = 1 4 ( | &psi; k &prime; yMr - 1 m &Sigma; r = 1 m &psi; k &prime; yMr | 1 m &Sigma; r = 1 m &psi; k &prime; yMr ) + | h 0 Mr - 1 m &Sigma; r = 1 m h 0 Mr | 1 m &Sigma; r = 1 m h 0 Mr y = 1,2,3,4 } ;
Subsequently, in step 15, by { α k ' raccording to the size of numerical value according to being descending sort from big to small, and definition sequence working array { n ll=1,2 ..., m} makes this array meet inequality &alpha; k &prime; n 1 &GreaterEqual; &alpha; k &prime; n 2 &GreaterEqual; &CenterDot; &CenterDot; &CenterDot; &GreaterEqual; &alpha; k &prime; n l - 1 &GreaterEqual; &alpha; k &prime; n l &GreaterEqual; &alpha; k &prime; n l + 1 &GreaterEqual; &CenterDot; &CenterDot; &CenterDot; &GreaterEqual; &alpha; k &prime; n m ;
Subsequently, in step 16, given rejecting proportionality coefficient χ=25%;
Subsequently, in step 17, get rid of the interference of accidentalia, at t 0to t 0+ τ fin the actual rolling technological parameter gathering constantly, weed out 25% data, then 75% remaining data are averaged, obtain t 0to t 0+ τ fconstantly interior for feeding back the feature rolling technological parameter of calculating (y=1,2,3,4), and feature supplied materials thickness m wherein 1represent t 0to t 0+ τ fthe number of the actual rolling parameter of rejecting constantly, m 1=int (χ m)=31, m 2represent t 0to t 0+ τ fconstantly, retain the number of the actual rolling parameter of not rejecting, m 2=m-m 1=94;
Subsequently, in step 18, according to t 0to t 0+ τ ffeature rolling technological parameter constantly, respectively the value of feedback S of bending roller force in the interior 1# of definite current feedback cycle, 2# frame 1M, S 2M, the roller amount of inclining value of feedback η 1M, η 2M, the value of feedback δ of roll shifting amount 1M, δ 2M, adopt the following step of being carried out by computer:
First, in step 18-1, objective definition function initial value F 0, and make F 0=10 10;
Subsequently, in step 18-2, according to target proportion crown curve coefficient a εcalculate the target proportion convexity Distribution Value of the second milling train outlet C i * = &Sigma; &epsiv; = 0 3 a &epsiv; x i &epsiv; = 0.99118097 0.993780015 0.995906505 0.997560443 0.998978103 1.000041348 1.000868317 1.001577147 1.002049701 1.002404116 1.002640393 1.002758531 1.002876669 1.002876669 1.002994808 1.002994808 1.002994808 1.002876669 1.002876669 1.002758531 1.002640393 1.002404116 1.002049701 1.001577147 1.000868317 1.000041348 0.998978103 0.997560443 0.995906505 1.993780015 0.99118097 (curve is as shown in figure 16), x i = - n 2 + i - 0.5 , i=0,1,2,...,31;
Subsequently, in step 18-3, make δ k ' M=(L k '-B-150)/2=77mm;
Subsequently, in step 18-4, make the roller amount η that inclines k '=0;
Subsequently, in step 18-5, definition bending roller force procedure parameter θ, and make θ=0;
Subsequently, in step 18-6, make bending roller force in the first frame
Figure BDA00001644947300261
Figure BDA00001644947300262
bending roller force in the second frame S 2 = - S 2 max - + &theta; &phi; 2 ( S 2 max - + S 2 max + ) = - 280 kN ;
Subsequently, in step 18-7, the Outlet Section thickness distribution value h that utilizes the described outlet thickness of slab of document [1] distribution machine reason model to calculate at the second rolling mill mi={ 0.8313 0.8359 0.8396 0.8426 0.8450.8468 0.8482 0.8492 0.85 0.8505 0.8508 0.851 0.8511 0.85120.8512 0.8513 0.8512 0.8512 0.8511 0.851 0.8508 0.8505 0.850.8492 0.8482 0.8468 0.845 0.8426 0.8396 0.8359 0.8313} and band forward pull cross direction profiles value σ 1Mi=152.2 157.4 162.5 167.4 172 176.3 180.2 183.6186.6 189.2 191.4 193.2 194.5 195.5 196.1 196.3 196.1 195.5 194.5193.2 191.4 189.2 186.6 183.6 180.2 176.3 172 167.4 162.5 157.4152.2}(curves are as shown in figure 17);
Subsequently, in step 18-8, calculate the maximum residual stress of the second rolling mill outlet band &sigma; M max = max { | &sigma; 1 Mi - &Sigma; i = 1 n &sigma; 1 Mi n | } = 28.2 MPa ;
Subsequently, in step 18-9, utilize the described method of document [2] to calculate the critical residual stress of current the second rolling mill outlet band generation plate shape unstability &sigma; Mcr = k cr &pi; 2 E 12 ( 1 + &upsi; ) ( h * B ) 2 = 19 MPa ;
Subsequently, in step 18-10, calculate the actual ratio convexity Distribution Value of the second rolling mill outlet band C Mi = h Mi ( &Sigma; i = 1 n h Mi ) / n = 0.982095918 0.987530347 0.991901519 0.995445712 0.998281066 1.000407582 1.002061539 1.003242937 1.004188055 1.004778754 1.005133173 1.005369453 1.005487592 1.005605732 1.005605732 1.005723872 1.005605732 1.005605732 1.005487592 1.005369453 1.005133173 1.004778754 1.004188055 1.003242937 1.002061539 1.000407582 0.998281066 0.995445712 0.991901519 0.987530347 0.982095918 (curve as shown in figure 18);
Subsequently, in step 18-11, calculate the deviation of actual ratio convexity Distribution Value and target proportion convexity Distribution Value, its function expression is F = &Sigma; i = 1 n ( C i * - C Mi ) 2 = 0.00038995 ;
Subsequently, in step 18-12, judgement inequality group &sigma; M max < &sigma; Mcr F < F 0 Set up? obviously, inequality group 28.2 < 19 0.0038995 < 10 10 Be false, make θ=θ+1=1, proceed to step 18-13; If but set up, make F 0=F, S 1M=S 1, S 2M=S 2, θ=θ+1, proceed to step 18-13;
Subsequently, in step 18-13, judgement inequality is thin &theta; &phi; 1 &le; 1.0 &theta; &phi; 2 &le; 1.0 Set up? obvious inequality group 0 &le; 1.0 0 &le; 1.0 Set up, proceed to step 18-6, repeating step 18-6 to 18-13, until no longer set up with upper inequality group; When inequality group is no longer set up, proceed to step 18-14;
Subsequently, in step 18-14, define the roller amount procedure parameter λ that inclines, and make λ=0, F 0=10 10;
Subsequently, in step 18-15, make roller declination amount in the first frame
Figure BDA000016449473002710
roller declination amount in the second frame
Figure BDA000016449473002711
Subsequently, in step 18-16, utilize the described outlet thickness of slab of document [1] distribution machine reason model to calculate the Outlet Section thickness distribution value C of the second rolling mill mi={ 0.8208 0.8287 0.8348 0.8397 0.84450.847 0.8481 0.8504 0.8519 0.8518 0.8527 0.8535 0.8529 0.85350.8527 0.8537 0.8536 0.8538 0.8529 0.8532 0.8525 0.8518 0.85190.8504 0.849 0.847 0.8448 0.8411 0.8369 0.8316 0.0253} and band forward pull cross direction profiles value σ 1Mi=161.5 163.5 165.6 167.8 169.9 171.9 173.9 175.6177.2 178.7 179.9 181 181.9 182.6 183.2 183.6 183.8 183.7183.5 182.9 182.2 181.2 180.1 178.7 177.2 175.5 173.7 171.8169.8 167.9 166.1}(curves are as shown in figure 19);
Subsequently, in step 18-17, calculate the maximum residual stress of the second rolling mill outlet band &sigma; M max = max { | &sigma; 1 Mi - &Sigma; i = 1 n &sigma; 1 Mi n | } = 14.5 MPa ;
Subsequently, in step 18-18, utilize the described method of document [2] to calculate the critical residual stress of current the second rolling mill outlet band generation plate shape unstability
Figure BDA00001644947300282
Subsequently, in step 18-19, calculate the actual ratio convexity Distribution Value of the second rolling mill outlet band C Mi = h Mi ( &Sigma; i = 1 n h Mi ) / n = 0.969971993 0.979307737 0.986516349 0.992306874 0.997979225 1.000933574 1.002233488 1.00495149 1.006724099 1.006605925 1.007669491 1.008614883 1.007905839 1.008614883 1.007669491 1.008851231 1.008733057 1.008969405 1.007905839 1.008260361 1.007433143 1.006605925 1.006724099 1.00495149 1.003297054 1.000933574 0.998333747 0.99396131 0.988998003 0.982734782 0.975289822 (curve as shown in figure 20);
Subsequently, in step 18-20, calculate target cross section characteristic (being that Target Board convexity distributes) and the deviation that actual convexity distributes, its function expression is
Figure BDA000016449473002810
Subsequently, in step 18-21, judge inequality &sigma; M max < &sigma; Mcr F < F 0 Set up? obviously, inequality group 14.5 < 19 0.001645 < 10 10 Be false, make λ=λ+1=1, proceed to step 18-22; If but set up, make F 0=F, η 1M1, η 2M2, λ=λ+1, proceed to step 18-22; ;
Subsequently, in step 18-22, judgement inequality group
Figure BDA00001644947300292
set up? obvious inequality group 0 &le; 1.0 0 &le; 1.0 Set up, proceed to step 18-15, repeating step c18-15 to c18-22, until no longer set up with upper inequality group; When inequality group is no longer set up, proceed to step 18-19; ;
Subsequently, in step 19, by t 0to t 0+ τ fthe value of feedback δ of roll shifting amount constantly k ' Mthe value of feedback S of=77mm, bending roller force 1M=129.5MPa, S 2Mthe value of feedback η of=201.95MPa, the roller amount of inclining 1M=0mm, η 2M=0mm sends to a level system of tandem mills, respectively roll shifting amount, the bending roller force of current first, second rolling mill, the roller amount of inclining is reset;
Do you subsequently, in step 20, judge that inequality ξ < 0 sets up? if inequality is set up, finish cross section characteristic feedback, the cross section characteristic reponse system of virtual profile gauge quits work, and proceeds to step 6; If inequality is false, make M=M+1, t 0=t 0+ τ f, proceed to step 8, until inequality ξ < 0 sets up.

Claims (2)

1. the virtual profile gauge establishing method of four roller tandem mills based on mechanism model, is characterized in that: described establishing method comprises the following step of being carried out by computer:
(a) collect the basic equipment parameter of each rolling mill of tandem mills, comprise the barrel length L of working roll wk, backing roll barrel length L bk, working roll bending cylinder is apart from l wk, backing roll housing screw centre-to-centre spacing l bk, working roll roller footpath D wk, backing roll roller footpath D bk, working roll roll shape D wki, backing roll roll shape D bki, the maximum positive bending roller force of working roll the maximum negative bending roller force of working roll
Figure FDA0000455937850000012
the maximum roller amount of just inclining
Figure FDA0000455937850000013
the maximum negative roller amount of inclining
Figure FDA0000455937850000014
subscript k represents shelf number, k=1, and 2 ..., s, s is total frame number;
(b) realization of Outlet Section shape Presentation Function, comprises the following step of being carried out by computer:
B1) collect the sampling period τ of current tandem mills data collecting system s;
B2) provide the period tau that virtual profile gauge section configuration dynamically shows x, this cycle is relevant to dominant frequency, the capacity of the computer of the virtual profile gauge of installation, and must guarantee τ xτ sintegral multiple and within this cycle, can complete twice thickness of slab and distribute and calculate;
B3) provide the interval parameter L to supplied materials sectional sampling 0;
B4) definition length of interval parameter j, supplied materials coil of strip length parameter L j;
B5) definition as-rolled condition parameter ξ, wherein ξ=1 represents to start rolling, ξ=-1 represents to stop rolling, and this instruction is sent according to on-site actual situations by site operation personnel;
B6) given initial steel reel number COILNO1=0, prepares convexity and dynamically shows;
B7) from the three-level system of tandem mills, collect supplied materials coil of strip essential information, comprise reel number COILNO, steel grade code SG, strip width B, supplied materials thickness h 0, last rolling mill exit thickness h s;
B8) judge inequality | does COILNO1-COILNO|>0 set up? if set up, make j=0, L j=0, and read the cross section characteristic data that supplied materials steel reel number is COILNO, proceed to step b9); If inequality | COILNO1-COILNO|>0 is false, and directly proceeds to step b9);
B9) collect current time the first frame entrance and roll band steel length L *;
B10) judgement inequality L *>=L jwhether set up, if set up, proceed to step b11; Otherwise, proceed to step b12;
B11) extract and calculate supplied materials and twist in [L j, L j+1] average crown curve function f in length of interval mj(X), L wherein j+1=L j+ L 0; X is illustrated in strip width direction and sets up one-dimensional coordinate system, when the center of width is the origin of coordinates, and the corresponding numerical value of strip width direction each point in coordinate system;
B12) by tandem mills data collecting system, collect the rolling technological parameter actual value ψ of current each frame ky, wherein, y represents to survey the classification of dynamic data, and y=1 represents that draught pressure, y=2 represent that forward pull, y=3 represent that backward pull, y=4 represent that reduction ratio, y=5 represent that roll shifting amount, y=6 represent bending roller force, the y=7 roller amount that represents to incline;
B13) utilize outlet thickness of slab distribution machine reason model to calculate actual rolling technological parameter for ψ kyand in supplied materials interval, average crown curve function is f mj(X), time, width is the Outlet Section thickness distribution value h of B, the steel grade code band that is SG i, i is band in horizontal bar unit number, i=1, and 2 ..., n, n is the total bar unit number of band;
B14) utilize and calculate a mouthful section thickness Distribution Value h ithe current section crown curve of matching function f c(X), X is illustrated in strip width direction and sets up one-dimensional coordinate system, when the center of width is the origin of coordinates, and the corresponding numerical value of strip width direction each point in coordinate system;
B15) calculate current section convex value CROWN and wedge shape value WEDGE, convex value computation model is
Figure FDA0000455937850000021
wedge shape amount computation model is
Figure FDA0000455937850000022
f wherein c(0) represent band steel center thickness, f c(B/2-40) represent the thickness apart from band 40mm place, steel one end, f c(B/2+40) represents the thickness apart from band steel other end 40mm place;
B16) utilize the dynamic display of functional of visual software, by current section crown curve function f c(X) dynamically show, show convex value CROWN, wedge shape value WEDGE simultaneously, and corresponding coil of strip information;
B17) do you judge that inequality ξ < 0 sets up? if inequality is set up, finish cross section characteristic and show, virtual plate profile instrument quits work; If inequality is false, COILNO1=COILNO, proceeds to step b7), until inequality ξ < 0 sets up;
(c) cross section characteristic closed loop feedback is controlled the realization of function, comprises the following step that can be carried out by computer:
C1) collect the target proportion crown curve coefficient a of 2# milling train outlet ε, wherein, ε=0,1,2,3,4;
C2) determine respectively in cross section characteristic feedback procedure the minimum coefficient φ that adjusts of single step of bending roller force in 1# and 2# frame 1, φ 2;
C3) determine respectively the minimum coefficient of adjusting of single step of the roller amount of inclining in 1# and 2# frame in cross section characteristic feedback procedure
Figure FDA0000455937850000032
C4) determine the period tau of cross section characteristic feedback f, this cycle is relevant to dominant frequency, the capacity of the computer of the virtual profile gauge of installation, and must guarantee τ fτ sintegral multiple and within this cycle, can complete
Figure FDA0000455937850000031
inferior thickness of slab distributes and calculates;
C5) definition rolling moment procedure parameter t 0, virtual profile gauge cross section characteristic Times of Feedback procedure parameter M;
C6) receive operational order, judge whether to start rolling? ξ=1 represents to start rolling, ξ=-1 represents to stop rolling, if started rolling, records current standard Beijing time t, and makes t 0=t, proceeds to step c7); If do not start rolling, enter wait state;
C7) given initial steel reel number COILNO1=0;
C8) from the three-level system of tandem mills, collect supplied materials coil of strip essential information, comprise reel number COILNO, steel grade code SG, strip width B, supplied materials thickness h 0, last rolling mill exit thickness h s;
C9) judge inequality | does COILNO1-COILNO|>0 set up? if set up, make j=0, L j=0, and read the cross section characteristic data that supplied materials steel reel number is COILNO, proceed to step c10; If inequality | COILNO1-COILNO|>0 is false, and directly proceeds to step c10;
C10) collect current time the first frame entrance and roll band steel length L *;
C11) judgement inequality L *>=L jset up? if set up, proceed to step c12; Otherwise, proceed to step c13;
C12) extract and calculate supplied materials and twist in [L j, L j+1] average crown curve function f in length of interval mj(X), L wherein j+1=L j+ L 0; X is illustrated in strip width direction and sets up one-dimensional coordinate system, when the center of width is the origin of coordinates, and the corresponding numerical value of strip width direction each point in coordinate system;
C13) by milling train data collecting system, collect 1#, 2# milling train from t 0to t 0+ τ factual rolling technological parameter ψ constantly k ' yMr; The shelf number of front two rolling mills of subscript k ' representative, k '=1,2; Subscript r represents t 0to t 0+ τ fall kinds of actual rolling technological parameter collected is constantly according to the numbering of time order and function order, r=1, and 2 ..., m, m is illustrated in t 0to t 0+ τ fthe number of the rolling technological parameter of collected particular type constantly,
Figure FDA0000455937850000042
C14) introduce parameter and reject process variables array { α k ' r, and by { α k ' rwith following formula, represent:
Figure FDA0000455937850000043
α k ' rrepresent that the k ' rolling mill is at t 0+ r τ srejecting variable constantly, the shelf number of front two rolling mills of subscript k ' representative, k '=1,2;
C15) by { α k ' raccording to the size of numerical value according to being descending sort from big to small, and definition sequence working array { n ll=1,2 ..., m} makes this array meet inequality n wherein lrepresentative { α k ' rcorresponding according to the numbering of time sequencing sequence according to l rejecting variable after numerical value descending sort, l represents { α k ' raccording to the big or small sequence number after numerical value descending sort;
C16) given rejecting proportionality coefficient χ, χ=0~50%;
C17) at t 0to t 0+ τ fin the actual rolling technological parameter gathering constantly, weed out the data of χ ratio, then remaining (1-χ) data are averaged, obtain t 0to t 0+ τ fconstantly interior for feeding back the feature rolling technological parameter ψ of calculating k ' yM(y=1,2,3,4), and feature supplied materials thickness h 0M, its expression formula is respectively
Figure FDA0000455937850000045
m wherein 1represent t 0to t 0+ τ fthe number of the actual rolling parameter of rejecting constantly, m 1=int (χ m), m 2represent t 0to t 0+ τ fconstantly, retain the number of the actual rolling parameter of not rejecting, m 2=m-m 1;
C18) according to t 0to t 0+ τ ffeature rolling technological parameter constantly, respectively the value of feedback S of bending roller force in the interior 1# of definite current feedback cycle, 2# frame 1M, S 2M, the roller amount of inclining value of feedback η 1M, η 2M, the value of feedback δ of roll shifting amount 1M, δ 2M;
C19) by t 0to t 0+ τ fthe value of feedback δ of roll shifting amount constantly k ' M, bending roller force value of feedback S k ' M, the roller amount of inclining value of feedback η k ' Msend to a level system of tandem mills, respectively roll shifting amount, the bending roller force of current 1#, 2# milling train, the roller amount of inclining are reset;
C20) do you judge that inequality ξ < 0 sets up? if inequality is set up, finish cross section characteristic feedback, the cross section characteristic reponse system of virtual profile gauge quits work, and proceeds to c6; If inequality is false, make M=M+1, t 0=t 0+ τ f, proceed to step c8, until inequality ξ < 0 sets up.
2. the virtual profile gauge establishing method of four roller tandem mills according to claim 1 based on mechanism model, is characterized in that: step c18) adopt the following step of being carried out by computer to complete:
C18-1) objective definition function initial value F 0, and make F 0=10 10;
C18-2) according to target proportion crown curve coefficient a εcalculate 2# milling train and export corresponding target proportion convexity Distribution Value
Figure FDA0000455937850000055
, fundamental formular is:
Figure FDA0000455937850000051
x in formula ifor the relative coordinate of used unit of computational process,
Figure FDA0000455937850000052
i is band steel bar unit number;
C18-3) make δ k ' M=(L wk '-B-150)/2;
C18-4) make the roller amount η that inclines k '=0;
C18-5) definition bending roller force procedure parameter θ, and make θ=0;
C18-6) make bending roller force in 1# frame
Figure FDA0000455937850000053
bending roller force in 2# frame
Figure FDA0000455937850000054
C18-7) utilize outlet thickness of slab distribution machine reason model to calculate feature rolling technological parameter for ψ k ' yM(y=1,2,3,4), feature supplied materials thickness are h 0M, and roll shifting amount is δ k ' M, the roller amount of inclining is η k ', bending roller force is S k ', average crown curve function is f in supplied materials interval mj(X) time, width is the band that B, steel grade code are SG, at the Outlet Section thickness distribution value h of the 2nd rolling mill miwith band forward pull cross direction profiles value σ 1Mi, i is band in horizontal bar unit number, i=1, and 2 ..., n, n is the total bar unit number of band;
C18-8) calculate the maximum residual stress σ of the 2nd rolling mill outlet band mmax, basic calculating model is:
Figure FDA0000455937850000061
C18-9) calculate the critical residual stress σ of current the 2nd rolling mill outlet band generation plate shape unstability mcr, basic calculating model is:
Figure FDA0000455937850000062
wherein, k crfor the limit stress coefficient of band generation plate shape unstability, the elastic modelling quantity that E, υ are belt steel material and Poisson's ratio, B is strip width, h *it is the 2nd rolling mill exit thickness of strip;
C18-10) calculate the actual ratio convexity Distribution Value of the 2nd rolling mill outlet band
C18-11) calculate the deviation of actual ratio convexity Distribution Value and target proportion convexity Distribution Value, its function expression is
Figure FDA0000455937850000064
C18-12) judgement inequality group
Figure FDA0000455937850000065
set up? if inequality is set up, make F 0=F, S 1M=S 1, S 2M=S 2, θ=θ+1, proceed to step c18-13; Otherwise, make θ=θ+1, proceed to step c18-13;
C18-13) judgement inequality group
Figure FDA0000455937850000066
set up? if inequality is set up, proceed to step c18-6, repeating step c18-6 to c18-13, until no longer set up with upper inequality group; Otherwise, proceed to step c18-14;
C18-14) define the roller amount procedure parameter λ that inclines, and make λ=0, F 0=10 10;
C18-15) make roller declination amount in 1# frame
Figure FDA0000455937850000071
roller declination amount in 2# frame
Figure FDA0000455937850000072
C18-16) utilize outlet thickness of slab distribution machine reason model to calculate feature rolling technological parameter for ψ k ' yM(y=1,2,3,4), feature supplied materials thickness are h 0M, and roll shifting amount is δ k ' M, bending roller force is S k ' M, the roller amount of inclining η k ', average crown curve function is f in supplied materials interval mj(X) time, width is the band that B, steel grade code are SG, at the Outlet Section thickness distribution value C of the second rolling mill miwith band forward pull cross direction profiles value σ 1Mi;
C18-17) calculate the maximum residual stress σ of the 2nd rolling mill outlet band mmax, basic calculating model is:
Figure FDA0000455937850000073
C18-18) calculate the critical residual stress σ of current the 2nd rolling mill outlet band generation plate shape unstability mcr, basic calculating model is:
Figure FDA0000455937850000074
C18-19) calculate the actual ratio convexity Distribution Value of the 2nd rolling mill outlet band
Figure FDA0000455937850000075
C18-20) calculate the deviation that Target Board convexity distributes and actual convexity distributes, its function expression is
Figure FDA0000455937850000076
C18-21) judgement inequality
Figure FDA0000455937850000077
set up? if inequality is set up, make F 0=F, η 1M1, η 2M2, λ=λ+1, proceed to step c18-22; Otherwise, make λ=λ+1, proceed to step c18-22;
C18-22) judgement inequality group
Figure FDA0000455937850000078
set up? if inequality is set up, proceed to step c18-15, repeating step c18-15 to c18-22, until no longer set up with upper inequality group; Otherwise, proceed to step c19.
CN201210150863.8A 2012-05-16 2012-05-16 Mechanism model-based virtual crown gauge setting method for four-roll cold continuous rolling machine set Expired - Fee Related CN102688896B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210150863.8A CN102688896B (en) 2012-05-16 2012-05-16 Mechanism model-based virtual crown gauge setting method for four-roll cold continuous rolling machine set

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210150863.8A CN102688896B (en) 2012-05-16 2012-05-16 Mechanism model-based virtual crown gauge setting method for four-roll cold continuous rolling machine set

Publications (2)

Publication Number Publication Date
CN102688896A CN102688896A (en) 2012-09-26
CN102688896B true CN102688896B (en) 2014-04-16

Family

ID=46854726

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210150863.8A Expired - Fee Related CN102688896B (en) 2012-05-16 2012-05-16 Mechanism model-based virtual crown gauge setting method for four-roll cold continuous rolling machine set

Country Status (1)

Country Link
CN (1) CN102688896B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104289528A (en) * 2013-07-18 2015-01-21 上海宝钢钢材贸易有限公司 Rolling tension control method of double-rack four-roller mill
CN104550267B (en) * 2014-12-30 2017-01-18 燕山大学 Strip steel thickness lateral distribution characteristic parameter extraction method
CN105203731B (en) * 2015-09-10 2017-03-01 北京首钢股份有限公司 A kind of strip steel cross-sectional outling defect local high spot quantization method and device
CN105251781B (en) * 2015-11-26 2017-11-28 首钢京唐钢铁联合有限责任公司 Control method of cold continuous rolling mill
CN106944483B (en) * 2017-03-22 2018-11-09 宁波宝新不锈钢有限公司 A kind of measurement control method of reversing cold mill band transverse gage Distribution Value
CN112122358B (en) * 2020-09-11 2021-05-28 燕山大学 Plate shape target curve setting method considering plate shape requirements of subsequent process

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE69121789T2 (en) * 1990-06-04 1997-04-03 Hitachi Ltd Control device for controlling a controlled system and control method therefor
DE19622825B4 (en) * 1996-06-07 2005-03-31 Betriebsforschungsinstitut VDEh - Institut für angewandte Forschung GmbH Presetting for cold rolling reversing stand
JP3319356B2 (en) * 1997-08-29 2002-08-26 住友金属工業株式会社 Method and apparatus for controlling thickness of material to be rolled by rolling mill
JP2000167611A (en) * 1998-12-03 2000-06-20 Toshiba Corp Method and device for collecting rolling data
JP4268582B2 (en) * 2004-10-14 2009-05-27 新日本製鐵株式会社 Plate thickness control method and plate thickness / shape non-interference control method
CN101618402B (en) * 2008-06-30 2011-04-13 宝山钢铁股份有限公司 Method for controlling planeness of cold-rolling strip steel
CN101602068B (en) * 2009-07-07 2011-08-17 东北大学 Control method and control system of tension in the process of rolling periodic variable-thickness strips
CN102029292B (en) * 2009-09-28 2014-04-30 宝山钢铁股份有限公司 Strip steel thickness feedforward control method based on mechanical property detection
CN102029294B (en) * 2009-09-28 2013-06-19 宝山钢铁股份有限公司 Control method of lateral thick difference of cold-rolling strip steels
CN102189117B (en) * 2010-03-16 2013-05-29 宝山钢铁股份有限公司 Cold rolled steel strip straightness feedforward control method based on transverse performance detection

Also Published As

Publication number Publication date
CN102688896A (en) 2012-09-26

Similar Documents

Publication Publication Date Title
CN102688896B (en) Mechanism model-based virtual crown gauge setting method for four-roll cold continuous rolling machine set
Nakasaki et al. Application of hydrostatic integration parameter for free-forging and rolling
CN101623708B (en) Plate-shape control integrated system and executing method thereof
CN102303050B (en) Rough-rolling width self-learning method
CN103586286B (en) Rolling schedule comprehensive optimization method for cold continuous rolling unit taking scratch prevention as objective
CN101513647B (en) Method for leveling strip produced by secondary cold rolling unit
CN103962391B (en) Rolling load optimization method for hot continuous finishing mill group
CN112100745B (en) Automobile girder steel mechanical property prediction method based on LDA theory
CN113434994B (en) Method for predicting cold rolling deformation resistance based on hot rolling process parameters of hot rolled original plate
CN101602067A (en) Five frame UCM tandem mills plate shapes and the online integrated control method of strip crown
CN101716607A (en) Method for controlling asymmetric transverses shifting plate shape of asymmetric bending roller of HC rolling mill
CN101648208A (en) Rolling procedure making method of aluminium single-stand cold-rolling machine
CN103567228B (en) Method for forecasting strip shape and leaning during abnormal rolling of ultrathin strip of six-roller mill
CN102671959B (en) Method for plate shape closed-loop control by using virtual plate gauge for six-roller flattening unit
CN107931329A (en) A kind of improvement CSP double fluids change the control method of specification rolling force model precision
CN110434172B (en) Load distribution calculation method for continuous rolling of furnace coil and finishing mill group
CN103357670A (en) Reduction schedule optimizing method applicable to five-stand UCM (universal crown mill) type cold continuous rolling unit
CN102553945A (en) Abnormal shape forecasting method suitable for four-high rolling mill
CN105032945B (en) A kind of hot tandem plate shape and plate convexity Comprehensive Control merit rating method
CN107716560B (en) A kind of Hot Strip Rolling load distribution method
CN101898202B (en) Method for forecasting edge reduction in rolling strips of SMS-EDC rolling mill
CN102161054A (en) Plate shape closed-loop control method based on influencing self learning of matrix
Ma et al. Influence of profile indicators of hot-rolled strip on transverse thickness difference of cold-rolled silicon steel
CN1641356B (en) Rigid line product microscopic structure and mechanical property forecasting system
Cho et al. Analysis of the cogging process for heavy ingots by finite element method and physical modelling method

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
TR01 Transfer of patent right

Effective date of registration: 20191028

Address after: 063600 No.1, Lingang Road, Xiangyun island forest farm, leting County, Tangshan City, Hebei Province

Patentee after: Tangshan City Delong Iron & Steel Co., Ltd.

Address before: Hebei Street West Harbor area, 066004 Hebei city of Qinhuangdao province No. 438

Patentee before: Yanshan University

TR01 Transfer of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140416

Termination date: 20200516

CF01 Termination of patent right due to non-payment of annual fee