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 PDFInfo
- 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
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
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
the maximum negative bending roller force of working roll
the maximum roller amount of just inclining
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
Wedge shape amount computation model is
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
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,
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:
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 χ, 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
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
fundamental formular is:
x in formula
ifor the relative coordinate of used unit of computational process,
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-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:
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
C18-12) judgement inequality group
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
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-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:
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:
C18-20) calculate target cross section characteristic (being that Target Board convexity distributes) and the deviation that actual convexity distributes, its function expression is
C18-21) judgement inequality
Set up? if inequality is set up, make F
0=F, η
1M=η
1, η
2M=η
2, λ=λ+1, proceed to step c18-22; Otherwise, make λ=λ+1, proceed to step c18-22;
C18-22) judgement inequality group
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
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
Wedge shape amount computation model is
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
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,
Subsequently, in step 14, introduce parameter and reject process variables array { α
k ' r,
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
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 2# milling train outlet
(curve is as shown in Figure 8),
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
kN, bending roller force in 2# frame
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
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
Subsequently, in step 18-10, calculate the actual ratio convexity Distribution Value of the second rolling mill outlet band
(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
Subsequently, in step 18-12, judgement inequality group
Set up? obviously, inequality group
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
Set up? obvious inequality group
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
roller declination amount in 2# frame
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
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
Subsequently, in step 18-19, calculate the actual ratio convexity Distribution Value of the second rolling mill outlet band
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
Subsequently, in step 18-21, judge inequality
Set up? obviously, inequality group
Be false, make λ=λ+1=1, proceed to step 18-22; If but set up, make F
0=F, η
1M=η
1, η
2M=η
2, λ=λ+1, proceed to step 18-22; ;
Subsequently, in step 18-22, judgement inequality group
set up? obvious inequality group
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.
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
Wedge shape amount computation model is
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
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,
Subsequently, in step 14, introduce parameter and reject process variables array { α
k ' r,
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
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
(curve is as shown in figure 16),
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
bending roller force in the second frame
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
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
Subsequently, in step 18-10, calculate the actual ratio convexity Distribution Value of the second rolling mill outlet band
(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
Subsequently, in step 18-12, judgement inequality group
Set up? obviously, inequality group
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
Set up? obvious inequality group
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
roller declination amount in the second frame
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
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
Subsequently, in step 18-19, calculate the actual ratio convexity Distribution Value of the second rolling mill outlet band
(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
Subsequently, in step 18-21, judge inequality
Set up? obviously, inequality group
Be false, make λ=λ+1=1, proceed to step 18-22; If but set up, make F
0=F, η
1M=η
1, η
2M=η
2, λ=λ+1, proceed to step 18-22; ;
Subsequently, in step 18-22, judgement inequality group
set up? obvious inequality group
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
the maximum roller amount of just inclining
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 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
wedge shape amount computation model is
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
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, 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,
C14) introduce parameter and reject process variables array { α
k ' r, and by { α
k ' rwith following formula, represent:
α
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
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
, fundamental formular is:
x in formula
ifor the relative coordinate of used unit of computational process,
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-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:
C18-9) 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
C18-12) judgement inequality group
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
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-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:
C18-18) calculate the critical residual stress σ of current the 2nd rolling mill outlet band generation plate shape unstability
mcr, basic calculating model is:
C18-20) calculate the deviation that Target Board convexity distributes and actual convexity distributes, its function expression is
C18-21) judgement inequality
set up? if inequality is set up, make F
0=F, η
1M=η
1, η
2M=η
2, λ=λ+1, proceed to step c18-22; Otherwise, make λ=λ+1, proceed to step c18-22;
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)
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)
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 |
-
2012
- 2012-05-16 CN CN201210150863.8A patent/CN102688896B/en not_active Expired - Fee Related
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 |