CN102527736A - Method for confirming reeling radius of cold-rolled steel strip - Google Patents

Abstract

The invention provides a method for confirming a reeling radius of a cold-rolled steel strip. The method comprises an off-line optimizing calculating stage and an online applying stage, wherein the off-line optimizing calculating stage comprises the following steps of: adjusting a solving sequence of all physical quantities of a reeling coupling model of the steel strip, introducing a grain group optimizing algorithm and solving a coupling calculating problem, thereby avoiding the local optimized point subjected by the calculating process and the generation of diffuse phenomenon. During an executing process of the online applying stage, a reeling radius data of the steel strip obtained in the off-line optimizing calculating stage is stored into a board computer database, and a high-precision instant reeling radius data of the steel strip is supplied to a sequential steel strip reeling compensating module in real time according to a database querying method after a board control system is in an online running process, so that the important reeling compensating function in the board control system is finished. According to the method, the calculating process is prevented from falling into the local optimized point or the generation of diffuse phenomenon is avoided, the increasing of the board control precision of the cold-rolled steel strip is ensured and the significant economic benefit is brought.

Description

Confirm that cold-strip steel batches the method for radius

Technical field

The present invention relates to the cold-strip steel field, relate in particular to the method that a kind of definite cold-strip steel batches radius.

Background technology

Strip shape quality is one of the key technical indexes of weighing cold-rolled steel strip products, through the overall process of band steel production.In order to improve the quality and the precision of the control of plate shape; Advanced cold rolled strip steel production line can be equipped with plate profile instrument usually and carry out the on-line measurement of belt plate shape in the milling train exit; Then measurement feedback is arrived plat control system, accomplish the plate shape closed-loop control function of whole system.In the actual production; Because belt steel product has certain convexity usually after cold rolling, cause strip coiling to form the volume convexity afterwards, cause roll coil of strip outer radius laterally not wait along the band steel; Finally cause being with the steel tightness to change; Produce additional stress and make band steel tensile stress cross direction profiles change, and then influence plate profile instrument, caused the variation of plate shape control effect indirectly the authenticity that plate shape is detected.Therefore, in order further to improve the control accuracy of plat control system, just must give effective compensation to this scroll disturbing factor of plate shape measurement data error that causes.

Can know that by domestic and international pertinent literature can the key that plat control system carries out scroll disturbing factor compensation to the plate shape measurement data be accurately and obtain the strip coiling outer radius values timely.Receive the external tension effect during strip coiling, directly have the interaction of compression layer by layer, and have complicated coupled relation between roll coil of strip outer radius and the compression.Domestic professors Lian Jiachuan of University On The Mountain Of Swallows etc. are through considering the radially characteristic that radially changes of tightening coefficient of coil of strip; Set up a kind of strip coiling model for coupling; In this Mathematical Modeling; Exist the complex mathematical relation that intercouples between strip coiling outer radius, the interlayer radial pressure, this also makes accurately finding the solution of this model become quite difficult.Document " iteration Method coiler mandrel unit pressure successively; " heavy-duty machinery "; connect family's wound etc., 2001, (6): 30-33 " has proposed a kind of iterative calculation method successively: at first batching inner each relative storey displacement through hypothesis is zero with the change in displacement increment; Calculate corresponding interlayer radial compressive stress P1 then; Utilize the interlayer radial pressure to calculate again and batch inner each relative storey displacement and change in displacement increment, utilize the interlayer radial pressure to calculate to batch inner each relative storey displacement and the change in displacement increment recomputates interlayer radial compressive stress P2, with twice calculates before and after the above-mentioned computational process interlayer radial compressive stress P1 and the deviation between the P2 as the measurement index function; If then thinking less than a certain setting value, the measurement index function finds the solution convergence; Otherwise make P1=P2, recomputate, till obtaining satisfied result of calculation according to above-mentioned steps.

It should be noted that above-mentioned traditional iterative method successively has very big dependence for the initial value of separating, when finding the solution the coupling model problem, can run into the thorny problem that result of calculation does not restrain, this also greatly reduces the using value and the realistic meaning of this method.On the other hand, intelligent optimization method is widely used in field of engineering technology having obtained significant progress in the more than ten years recently.Different with conventional calculation; Intelligent optimization method is through simulation and study artificial intelligence; Effectively overcome complicated calculations problem solving result and do not restrained or be absorbed in problem such as local optimum point easily, become the effective tool that solves traditional engineering technology problem.

, the calculating cold-strip steel can run into the technical barrier that result of calculation does not restrain in order to solve traditional alternative manner successively when batching radius; The present invention is through the order of finding the solution between each physical quantity of adjustment strip coiling coupling model; Introduce population (PSO) optimized Algorithm and carry out finding the solution of The Coupling problem, successfully avoided computational process to be absorbed in the generation that local optimum point is perhaps dispersed phenomenon; In addition; The coil of strip that offline optimization is calculated batches radius data and deposits in the plate shape Computer Database;, plat control system can for follow-up band coil of strip shape compensating module high-precision instantaneity coil of strip radius data be provided in real time after dropping into on-line operation through the method for data base querying; Accomplish scroll compensate function comparatively crucial in the plat control system, improve assurance for improving the cold-rolled steel strip products strip shape quality.

Summary of the invention

Technical problem to be solved by this invention is: provide a kind of definite cold-strip steel to batch the method for radius, often run into the technical problem that result of calculation does not restrain to solve traditional iterative calculation method successively.

The technical scheme that the present invention solves its technical problem employing is: comprise that offline optimization calculates link and two parts of online applied links.Wherein: it is through the order of finding the solution between each physical quantity of adjustment strip coiling coupling model that offline optimization calculates link, introduces particle swarm optimization algorithm and carries out finding the solution of The Coupling problem, is absorbed in the generation that local optimum point is perhaps dispersed phenomenon to avoid computational process; In the implementation process of online applied links; Be offline optimization to be calculated the resulting coil of strip of link batch radius data and deposit in the plate shape Computer Database; After plat control system drops into on-line operation, can for follow-up band coil of strip shape compensating module high-precision instantaneity coil of strip radius data be provided in real time, accomplish scroll compensate function comparatively crucial in the plat control system through the method for data base querying.

Described offline optimization calculates link and can may further comprise the steps:

(1) collect the related process parameter that cold-strip steel batches radius calculation:

This related process parameter comprises: each measuring section size of the elasticity modulus of materials of cold-strip steel, cold-strip steel Poisson's ratio, cold-strip steel contacting surface asperity maximum, cold-strip steel width, plate profile instrument measuring roller, cold-strip steel transverse gage distribute, the cold-strip steel coiling tension distributes;

(2) confirm bar unit number:

Confirm shaping unit number according to each measuring section size of strip width and plate profile instrument measuring roller; Specifically: will be with the steel center line to align with plate profile instrument measuring roller center line; The plate profile instrument measuring section that each is covered with steel fully is one and calculates bar unit; Both sides are not are not then accepted or rejected according to the coverage rate size by the measuring section that covers fully, accept or reject rule and be if coverage rate during greater than critical coverage rate this measuring section first as a calculating bar, otherwise cast out this measuring section;

(3) obtain cold-strip steel and batch radius data:

Accomplish offline optimization and calculate link, calculate the cold-strip steel that produces in the dynamic coiling process successively and batch radius data.

In above-mentioned steps (2), should set critical coverage rate by operating personnel, this critical coverage rate is >=50%.

In above-mentioned steps (2), distribute and tensile stress distributes all along band steel center line when symmetrical at belt steel thickness, only choose and be with the calculating bar unit of steel center line one side to calculate, result of calculation and this side result of calculation of an other side are symmetry equivalent relation at this moment.

In above-mentioned steps (3), can adopt following method to obtain cold-strip steel and batch radius data:

Step 1: make j=1, promptly batch ground floor and batch dynamically successively calculating of radius since i bar unit.

Step 2: according to roll coil of strip interlayer radial compressive stress roughly the regularity of distribution generate n group roll coil of strip interlayer radial compressive stress p at random _{I, k, j}, the p of each group wherein _{I, k, j}Create-rule be:

p _i，k，j＝p _i，k，j-1+c×rand

In the formula: p _{I, j, j-1}=0; If during j=1 k=1,2, Λ, j are arranged; C is the random number penalty coefficient; Rand is that the one-dimensional random on [0,1] interval is counted generating function.

Step 3: according to the n group p that obtains _{I, k, j}, utilize following mathematical formulae to calculate n group Δ u respectively _{I, k, j}And r _{I, k, j}:

$m_i=1+k\frac{{\mathrm{\δ}}_{\max }}{h_i}{\left(\frac{E}{p_{i,k,j}}\right)}^{0.8},$

In the formula: m _iIt is the coil of strip tightening coefficient of i bar unit; K is for having considered a correction factor theoretical and that actual value exists deviation to do, k=1,2, Λ, j; δ _MaxBe band steel contacting surface asperity maximum, it comprises microcosmic surface unevenness and meso-and mega-relief; h _iBe thickness before the i bar unit strip coiling; E is the elastic modelling quantity of belt steel material; p _{I, k, j}Be i bar unit inboard interlayer radial compressive stress of k layer when batching the j layer;

${\mathrm{\ϵ}}_{\mathrm{ri},k,j-1}=\frac{m_i}{E}[(1-\frac{v}{m_i})p_{i,k,j-1}-\frac{v}{m_i}({\mathrm{\σ}}_{\mathrm{Ti},j-1}+r_{i,k,j-1}\frac{p_{i,k+1,j-1}-p_{i,k,j-1}}{h_{i,k,j-1}})],$

${\mathrm{\ϵ}}_{\mathrm{ri},k,j}=\frac{m_i}{E}[(1-\frac{v}{m_i})p_{i,k,j}-\frac{v}{m_i}({\mathrm{\σ}}_{\mathrm{Ti},j}+r_{i,k,j}\frac{p_{i,k+1,j}-p_{i,k,j}}{h_{i,k,j}})],$

In the formula: v is the material Poisson's ratio; σ _{Ti, j}Be the coiling tension of i bar unit when batching the j layer, unit is MPa; r _{I, k, j}Be i bar unit k layer radius corresponding when batching the j layer; h _{I, k, j}Be that i bar unit k layer when batching the j layer batches back thickness; ε _{Ri, k, j}It is the radial deformation of i bar unit k layer band steel when batching the j layer;

Δε _ri，k，j＝ε _ri，k，j-ε _ri，k，j-1，

In the formula: Δ ε _{Ri, k, j}It is the k layer band steel of i bar unit radially deformation variable quantity when batching the j-1 layer when batching the j layer;

Δh _i，k，j＝h _iΔε _ri，k，j，

In the formula: Δ h _{I, k, j}It is the k layer band steel of i bar unit amounts of thickness variation when batching the j-1 layer when batching the j layer;

$u_{0i}=\frac{r_{0i}}{E_c}[\frac{r_{0i}^2+r_{\mathrm{ci}}^2}{r_{0i}^2-r_{\mathrm{ci}}^2}-v_c](p_{i,1,j}-p_{i,1,j-1}),$

In the formula: u _0iBe the displacement that the corresponding reel outer surface of i bar unit takes place when batching the j-1 layer when strip coiling j layer, r _0iBe i the corresponding reel outer radius of bar unit, r _CiBe i the corresponding reel inside radius of bar unit, E _cElasticity modulus of materials for reel;

Δu _i，1，j＝u _0i+Δh _i，1，j，

In the formula: Δ u _{I, 1, j}Be that i bar unit ground floor band steel batching the displacement that takes place behind the j layer;

Δu _i，k，j＝Δu _i，k-1，j+Δh _i，k，j，

In the formula: Δ u _{I, k, j}Be that i the bar k of unit layer band steel batching the displacement that takes place behind the j layer;

r _i，k，j＝r _i，k-1，j+h _i，k，j-Δu _i，k，j，

In the formula: r _{I, k, j}Be that i bar k of unit layer band steel batches radius after batching the j layer.

Step 4: with the n group Δ u that obtains in the step 3 _{I, k, j}And r _{I, k, j}The following strip coiling radial stress of substitution computing formula is asked for the n group respectively

$p_{i,k,m}^{*}=\frac{1}{(1-v)\frac{r_{i,k,j}}{E}-\frac{r_{i,k,\mathrm{<figure-callout\; id="2"\; label="j"\; filenames="HDA0000129010920000021.png,HDA0000129010920000022.png"\; state="\{\{state\}\}">j</figure-callout>}}^2}{{\mathrm{Eh}}_{i,k,j}}}[{\mathrm{\&Delta;u}}_{i,k,j}+u_{i,k,j-1}-\frac{r_{i,k,j}}{E}({\mathrm{\&sigma;}}_i+r_{i,k,j}\frac{p_{i,k+1,j}}{h_{i,k,j}})].$

Step 5: with

Judge for the measurement index function and respectively to organize p _{I, k, j}With

Between error size; If there is one group of p _{I, k, j}With

Make J≤ε set up, then go to step 6; Then produce new n group p if be false according to PSO algorithm particle create-rule _{I, k, j}, the p that each group is new _{I, k, j}Create-rule following:

V _i，k，j＝V _i，k，j+c ₁×rand×(pb _i，k，j-p _i，k，j)+c ₂×rand×(gb _i，k，j-p _i，k，j)

，

p _i，k，j＝p _i，k，j+V _i，k，j

In the formula: V _{I, k, j}Be PSO algorithm computation p _{I, k, j}Velocity component, the initial velocity component is one and generates number, k=1,2, Λ, j at random; c ₁And c ₂Acceleration weight for algorithm; Pb _{I, k, j}Corresponding p when obtaining minimum J value when this cycle calculations for this group particle _{I, k, j}Gb _{I, k, j}The p of correspondence when obtaining minimum J value when this cycle calculations for n group particle _{I, k, j},

Generating new n group p _{I, k, j}After, go to step 3;

Step 6: if j=N calculates and finishes, N batches the number of plies for the roll coil of strip maximum of confirming according to practical condition here.If J is not equal to N, then makes j=j+1, and forward step 2 continuation calculating to.

The present invention compared with prior art has following main beneficial effect:

1. offline optimization calculates link through the order of finding the solution between each physical quantity of adjustment strip coiling coupling model; Introduce population (PSO) optimized Algorithm and carry out finding the solution of The Coupling problem, successfully avoided computational process to be absorbed in the generation that local optimum point is perhaps dispersed phenomenon;

2. online applied links is calculated coil of strip that link obtains with offline optimization and is batched radius data and deposit in the plate shape Computer Database; After plat control system drops into on-line operation, can for follow-up band coil of strip shape compensating module high-precision instantaneity coil of strip radius data be provided in real time, accomplish scroll compensate function comparatively crucial in the plat control system through the method for data base querying.

3. improved assurance for improving the cold-rolled strip steel shape control accuracy, can considerable economic benefit.

Description of drawings

Fig. 1 is the strip coiling radius calculation flow chart of the inventive method offline optimization calculating section.

Fig. 2 is a cold-strip steel transverse gage distribution map in this instance.

Fig. 3 is a cold-strip steel coiling tension distribution map in this instance.

Internal compressive stress three-dimensional distribution map when the coil of strip that Fig. 4 obtains for the inventive method offline optimization calculating link batches.

Fig. 5 batches radius distribution figure for each bar unit of coil of strip that the inventive method offline optimization calculating link is got.

Fig. 6 calculates the error of calculation tendency chart of link for the inventive method offline optimization.

Fig. 7 uses front and rear panel shape control effect comparison diagram for the inventive method input.

The specific embodiment

The core calculations of the inventive method partly is divided into offline optimization and calculates link and two parts of online applied links; Wherein offline optimization calculates link through the order of finding the solution between each physical quantity of adjustment strip coiling coupling model; Introduce population (PSO) optimized Algorithm and carry out finding the solution of The Coupling problem, avoided computational process to be absorbed in the generation that local optimum point is perhaps dispersed phenomenon; The main effect of online applied links is offline optimization to be calculated the resulting coil of strip of link batch radius data and deposit in the plate shape Computer Database; After plat control system drops into on-line operation, can for follow-up band coil of strip shape compensating module high-precision instantaneity coil of strip radius data be provided in real time, accomplish scroll compensate function comparatively crucial in the plat control system through the method for data base querying.The inventive method simple and feasible satisfies the real-time requirement that existing cold-rolled strip steel shape is controlled system fully simultaneously, for further improving the cold-rolled steel strip products strip shape quality strong guarantee is provided.

Below in conjunction with embodiment and accompanying drawing the present invention is described further, but does not limit the present invention.

The method of batching radius based on definite cold-strip steel of the present invention can be used for four rollers, six roller single chassis or multimachine frame tandem mills.Present embodiment is disclosed to be the method that the definite cold-strip steel in certain single chassis six roller reversable cold-rolling machine plat control system batches radius.This six-high cluster mill can comprise common plate, high-strength steel, part stainless steel and silicon steel etc. by rolling product.What present embodiment was rolling is middle high grade silicon steel, and type is the UCM milling train, and plate shape control device comprises roller declination, the positive and negative roller of working roll, the positive roller of intermediate calender rolls, intermediate roll shifting and emulsion section cooling etc.The basic mechanical design feature index and the device parameter of this unit are:

Mill speed: Max 900m/min, draught pressure: Max 18000KN, maximum rolling force square: 140.3KN * m, coiling tension: Max 220KN, main motor current: 5500KW;

Supplied materials thickness range: 1.8～2.5mm, supplied materials width range: 850～1280mm, outgoing gauge scope: 0.3mm～1.0mm;

Work roll diameter: 290～340mm, working roll height: 1400mm, intermediate calender rolls diameter: 440～500mm, intermediate calender rolls height: 1640mm, backing roll diameter: 1150～1250mm, backing roll height: 1400mm;

Every side work roll bending power :-280～350KN, every side intermediate calender rolls bending roller force: 0～500KN, the axial traversing amount of intermediate calender rolls :-120～120mm, auxiliary hydraulic system pressure: 14MPa, balance roller system pressure: 28MPa, press down system pressure: 28MPa;

The plate shape measurement device adopts the ABB plate shape roller of Sweden; This plate shape roller roller footpath 313mm; Form by the single solid steel axle; Broad ways is every to be divided into a measured zone at a distance from 52mm, around measuring roller, is uniform-distribution with four grooves in each measured zone vertically to place magnetoelasticity power sensor, and the outside of sensor is wrapped up by steel loop.Product specification (thickness * width): 0.80mm * 840mm.Reel outer radius initial value is 300mm, and reel inside radius initial value is 100mm, and the elasticity modulus of materials of reel is taken as E _c=2.06 * 10 ⁵

Because the band steel finished product after rolling has certain convexity usually; Form the volume convexity after causing strip coiling; Cause roll coil of strip outer radius laterally not wait, finally cause being with the steel tightness to change, produce additional stress and make band steel tensile stress cross direction profiles change along the band steel; And then influence plate profile instrument to the authenticity that plate shape is detected, finally caused the variation of plate shape control effect.Production practices show, if do not give effective compensation to this scroll disturbing factor of plate shape measurement data error that causes, will increase the defect rate of product greatly, have a strong impact on economic benefit of enterprises.From the control technology angle, can the key that plat control system carries out scroll disturbing factor compensation to the plate shape measurement data be accurately and obtain the strip coiling outer radius values timely; But because the many factors of the strip coiling outer radius calculation Design operation of rolling; And has the complicacy contact that intercouples between each factor; Utilize the routine techniques method to be difficult to obtain satisfied result, this also becomes the technical barrier that obstruction cold-rolled strip steel shape control system realizes effective scroll compensate function.

The radius method of batching definite cold-strip steel that present embodiment provides adopts and comprises the steps:

1. collect the related process parameter that cold-strip steel batches radius calculation:

The related process coefficient that cold-strip steel batches radius calculation is the reference frame of calculating strip coiling half w footpath, and the collected technological parameter of this instance mainly contains: elasticity modulus of materials E=2.06 * 10 of cold-strip steel ⁵, cold-strip steel Poisson's ratio μ=0.3, cold-strip steel contacting surface asperity maximum δ _Max=1.52 * 10 ^-5M, cold-strip steel width B=840mm, each measuring section width dimensions of plate profile instrument measuring roller are 52mm; Distribution of cold-strip steel transverse gage and coiling tension distribute respectively as shown in Figures 2 and 3 in addition.

2. determine according to each measuring section size of strip width and plate profile instrument measuring roller and calculate bar unit number:

The main application that cold-strip steel batches radius calculation is the instant outside dimension that the horizontal zones of different of coil of strip is provided for scroll compensating module comparatively crucial in the plat control system; Therefore calculating bar unit number establishes rules really and then is: will be with the steel center line to align with plate profile instrument measuring roller center line; The plate profile instrument measuring section that each is covered with steel fully is one and calculates bar unit; Both sides are not are not then accepted or rejected according to the coverage rate size by the measuring section that covers fully; Accept or reject rule and be if coverage rate during greater than critical coverage rate this measuring section first as a calculating bar, otherwise cast out this measuring section.Here critical coverage rate can be come artificial the setting according to technological requirement, is not less than 50% usually.In particular cases; Distribute and tensile stress distributes all when being with the steel center line symmetrical like belt steel thickness; The calculating bar unit that can only choose band steel center line one side is calculated, and this moment, result of calculation and this side result of calculation broad ways center line of an other side were symmetry equivalent relation.

Because strip width is 840mm in this instance; Each measuring section width dimensions of plate profile instrument measuring roller is 52mm; Therefore be with steel to cover 840/52=16.1538 measuring section altogether; Because therefore two measuring section coverage rates very little (less than 8%) of milling train fore side and transmission side outermost end can cast out milling train fore side and two measuring sections of transmission side outermost end, that is to say only needs to calculate 16 bar units; In addition; Can know by Fig. 2 and Fig. 3; Belt steel thickness distributes and tensile stress distributes all along band steel center line symmetry; Calculate so can only choose 8 bar units of band steel center line one side this moment, result of calculation and this side result of calculation broad ways center line of 8 bar units of opposite side are symmetry equivalent relation.

3. accomplish offline optimization and calculate link, calculate the cold-strip steel that produces in the dynamic coiling process successively and batch radius data.

Fig. 1 provided this instance offline optimization calculate link i bar unit (8 bar units altogether) batch the radius calculation flow chart, the calculation procedure of this part is divided into:

Step 1: make j=1, promptly batch ground floor and batch dynamically successively calculating of radius since i bar unit;

Step 2: according to roll coil of strip interlayer radial compressive stress roughly the regularity of distribution generate n group roll coil of strip interlayer radial compressive stress p at random _{I, k, j}(k=1,2, Λ, j), wherein each organizes roll coil of strip interlayer radial compressive stress p _{I, k, j}(k=1,2, Λ, create-rule j) is:

p _i，k，j＝p _i，k，j-1+c×rand，

P is wherein arranged _{I, k, j-1}=0; If p is arranged during j=1 _{I, k, j-1}=0; C is the random number penalty coefficient; Rand is that the one-dimensional random on [0,1] interval is counted generating function.

Step 3: according to the n group p that obtains _{I, k, j}(k=1,2, Λ, j), utilization is batched Mathematical Modeling like next group cold-strip steel and is calculated n group Λ u respectively _{I, k, j}And r _{I, k, j}:

$m_i=1+k\frac{{\mathrm{\&delta;}}_{\max }}{h_i}{\left(\frac{E}{p_{i,k,j}}\right)}^{0.8},$

In the formula: m _iIt is the coil of strip tightening coefficient of i bar unit; K is for having considered a correction factor theoretical and that actual value exists deviation to do, and value is k=0.45 in this instance; δ _MaxBe band steel contacting surface asperity maximum, it comprises microcosmic surface unevenness and meso-and mega-relief, is δ according to correlation technique data value _Max=1.52 * 10 ^-5h _iBe thickness before the i bar unit strip coiling, its value is as shown in Figure 2; E=2.06 * 10 ⁵Elastic modelling quantity for belt steel material; p _{I, k, j}Be i bar unit inboard interlayer radial compressive stress of k layer when batching the j layer;

${\mathrm{\&epsiv;}}_{\mathrm{ri},k,j-1}=\frac{m_i}{E}[(1-\frac{v}{m_i})p_{i,k,j-1}-\frac{v}{m_i}({\mathrm{\&sigma;}}_{\mathrm{Ti},j-1}+r_{i,k,j-1}\frac{p_{i,k+1,j-1}-p_{i,k,j-1}}{h_{i,k,j-1}})],$

${\mathrm{\&epsiv;}}_{\mathrm{ri},k,j}=\frac{m_i}{E}[(1-\frac{v}{m_i})p_{i,k,j}-\frac{v}{m_i}({\mathrm{\&sigma;}}_{\mathrm{Ti},j}+r_{i,k,j}\frac{p_{i,k+1,j}-p_{i,k,j}}{h_{i,k,j}})],$

In the formula: v=0.3 is the material Poisson's ratio; σ _{Ti, j}Be the coiling tension (unit be MPa) of i bar unit when batching the j layer, its value is as shown in Figure 3; r _{I, k, j}Be i bar unit k layer radius corresponding when batching the j layer, h _{I, k, j}Be that i bar unit k layer when batching the j layer batches back thickness, ε _{Ri, k, j}It is the radial deformation of i bar unit k layer band steel when batching the j layer;

Δε _ri，k，j＝ε _ri，k，j-ε _ri，k，j-1，

In the formula: Δ ε _{Ri, k, j}It is the k layer band steel of i bar unit radially deformation variable quantity when batching the j-1 layer when batching the j layer;

Δh _i，k，j＝h _iΔε _ri，k，j，

In the formula: Δ h _{I, k, j}It is the k layer band steel of i bar unit amounts of thickness variation when batching the j-1 layer when batching the j layer;

$u_{0i}=\frac{r_{0i}}{E_c}[\frac{r_{0\mathrm{<figure-callout\; id="2"\; label="i"\; filenames="HDA0000129010920000021.png,HDA0000129010920000022.png"\; state="\{\{state\}\}">i</figure-callout>}}^2+r_{\mathrm{ci}}^2}{r_{0i}^2-r_{\mathrm{ci}}^2}-v_c](p_{i,1,j}-p_{i,1,j-1}),$

In the formula: u _0iBe the displacement that the corresponding reel outer surface of i bar unit takes place when batching the j-1 layer when strip coiling j layer, r _0iBe i the corresponding reel outer radius of bar unit, r _CiBe i the corresponding reel inside radius of bar unit, E _cElasticity modulus of materials for reel;

Δu _i，1，j＝u _0i+Δh _i，1，j，

In the formula: Δ u _{I, l, j}Be that i bar unit ground floor band steel batching the displacement that takes place behind the j layer;

Δu _i，k，j＝Δu _i，k-1，j+Δh _i，k，j，

In the formula, Δ u _{I, k, j}Be that i the bar k of unit layer band steel batching the displacement that takes place behind the j layer;

r _i，k，j＝r _i，k-1，j+h _i，k，j-Δu _i，k，j，

In the formula, r _{I, k, j}Be that i bar k of unit layer band steel batches radius after batching the j layer.

Step 4: the n group Δ u that will in step 3, obtain _{I, k, j}And r _{I, k, j}The following strip coiling radial stress of substitution computing formula is asked for the n group respectively

$p_{i,k,m}^{*}=\frac{1}{(1-v)\frac{r_{i,k,j}}{E}-\frac{r_{i,k,\mathrm{<figure-callout\; id="2"\; label="j"\; filenames="HDA0000129010920000021.png,HDA0000129010920000022.png"\; state="\{\{state\}\}">j</figure-callout>}}^2}{{\mathrm{Eh}}_{i,k,j}}}[{\mathrm{\&Delta;u}}_{i,k,j}+u_{i,k,j-1}-\frac{r_{i,k,j}}{E}({\mathrm{\&sigma;}}_i+r_{i,k,j}\frac{p_{i,k+1,j}}{h_{i,k,j}})].$

Step 5: with

Judge for the measurement index function and respectively to organize p _{I, k, j}With

Between error size.If there is one group of p _{I, k, j}With

Make J≤ε set up, then go to step 6; Then produce new n group p if be false according to PSO algorithm particle create-rule _{I, k, j}(k=1,2, Λ, j), the p that each group is new _{I, k, j}(k=1,2, Λ, create-rule j) is following:

V _i，k，j＝V _i，k，j+c ₁×rand×(pb _i，k，j-p _i，k，j)+c ₂×rand×(gb _i，k，j-p _i，k，j)

p _i，k，j＝p _i，k，j+V _i，k，j

V wherein _{I, k, j}Be PSO algorithm computation p _{I, k, j}Velocity component, the initial velocity component is one and generates number, c at random ₁And c ₂Be the acceleration weight of algorithm, pb _{I, k, j}Corresponding p when obtaining minimum J value when the cycle calculations for this group particle _{I, k, j}, gb _{I, k, j}The p of correspondence when obtaining minimum J value when the cycle calculations for n group particle _{I, k, j}

Generating new n group p _{I, k, j}(k=1,2, Λ goes to step 3 after j).

Step 6:, forward step 2 to and continue to calculate if j=N (N batches the number of plies for the roll coil of strip maximum of confirming according to practical condition here) calculates and finishes, otherwise makes j=j+1.

For this instance, it is N=60 that the setting maximum is batched the number of plies, and this moment, resulting result dispersed, and that is to say that conventional method can not be used for the calculating of this instance if adopt traditional iterative calculation method successively.We adopt the offline optimization computational methods that the present invention proposes in this example, set-up and calculated resultant error index need less than ε (here the ε value be satisfy commercial Application 0.01), this method can obtain satisfied result of calculation.Internal compressive stress three-dimensional distribution map when Fig. 4 has provided the inventive method offline optimization and calculates coil of strip that link obtains and batch; Fig. 5 has provided the inventive method offline optimization and has calculated coil of strip each bar unit that link gets and batch radius distribution figure; Result of calculation and field data are identical basically, have bigger using value.Fig. 6 has provided the error of calculation tendency chart of the inventive method offline optimization calculating link, and as can be seen from Figure 6, the error criterion function J of gained satisfies the computational accuracy requirement all less than 0.01.In sum, calculated off line link of the present invention has successfully been avoided computational process to be absorbed in local optimum point or has been dispersed the generation of phenomenon, for obtaining high-precision strip coiling radius data assurance is provided.

4. the online applied links of this method: at first will calculate cold-strip steel that link obtains and batch radius data and deposit plate shape Computer Database at offline optimization; After plat control system puts into operation; Cold-strip steel has whenever batched one deck; The control system accomplishes the scroll compensate function in the plat control system quickly and easily through the query manipulation of database is immediately batched radius value for band coil of strip shape compensating module provides each bar unit after batching this layer.

Fig. 7 has provided in this instance and has used the plate shape control effect comparison diagram that the same size belt steel products are rolled up in front and back two in the inventive method input.Can find out that the inventive method can improve the cold-rolled strip steel shape control accuracy, significantly reduce the flatness defect of cold-rolled steel strip products, improve the quality and the class of product, can be enterprise's considerable economic benefit.

Above embodiment only is used to explain calculating thought of the present invention and characteristics; Its purpose is to make those skilled in the art can understand content of the present invention and implements according to this; Protection scope of the present invention is not limited to the foregoing description; The disclosed principle of all foundations, equivalent variations or the modification that mentality of designing is done are all within protection scope of the present invention.

Claims (6)

1. method that definite cold-strip steel batches radius; It is characterized in that comprising that offline optimization calculates link and two parts of online applied links; Wherein: it is through the order of finding the solution between each physical quantity of adjustment strip coiling coupling model that offline optimization calculates link; Introduce particle swarm optimization algorithm and carry out finding the solution of The Coupling problem, be absorbed in the generation that local optimum point is perhaps dispersed phenomenon to avoid computational process; In the implementation process of online applied links; Be offline optimization to be calculated the resulting coil of strip of link batch radius data and deposit in the plate shape Computer Database; After plat control system drops into on-line operation, can for follow-up band coil of strip shape compensating module high-precision instantaneity coil of strip radius data be provided in real time, accomplish scroll compensate function comparatively crucial in the plat control system through the method for data base querying.

2. definite cold-strip steel according to claim 1 batches the method for radius, it is characterized in that offline optimization calculates link and may further comprise the steps:

(1) collect the related process parameter that cold-strip steel batches radius calculation:

This related process parameter comprises: each measuring section size of the elasticity modulus of materials of cold-strip steel, cold-strip steel Poisson's ratio, cold-strip steel contacting surface asperity maximum, cold-strip steel width, plate profile instrument measuring roller, cold-strip steel transverse gage distribute, the cold-strip steel coiling tension distributes;

(2) confirm bar unit number:

Confirm shaping unit number according to each measuring section size of strip width and plate profile instrument measuring roller; Specifically: will be with the steel center line to align with plate profile instrument measuring roller center line; The plate profile instrument measuring section that each is covered with steel fully is one and calculates bar unit; Both sides are not are not then accepted or rejected according to the coverage rate size by the measuring section that covers fully, accept or reject rule and be if coverage rate during greater than critical coverage rate this measuring section first as a calculating bar, otherwise cast out this measuring section;

(3) obtain cold-strip steel and batch radius data:

Accomplish offline optimization and calculate link, calculate the cold-strip steel that produces in the dynamic coiling process successively and batch radius data.

3. definite cold-strip steel according to claim 2 batches the method for radius, it is characterized in that in step (2), should setting critical coverage rate by operating personnel, and this critical coverage rate is >=50%.

4. definite cold-strip steel according to claim 2 batches the method for radius; It is characterized in that in the step (2); Distribute and tensile stress distributes all when being with the steel center line symmetrical at belt steel thickness; The calculating bar unit that only chooses band steel center line one side is calculated, and this moment, result of calculation and this side result of calculation of an other side were symmetry equivalent relation.

5. definite cold-strip steel according to claim 2 batches the method for radius, it is characterized in that adopting following method to obtain cold-strip steel and batching radius data in the step (3):

Step 1: make j=1, promptly batch ground floor and batch dynamically successively calculating of radius since i bar unit;

Step 2: according to roll coil of strip interlayer radial compressive stress roughly the regularity of distribution generate n group roll coil of strip interlayer radial compressive stress p at random _{I, k, j}, the p of each group wherein _{I, k, j}Create-rule be:

p _i，k，j＝p _i，k，j-1+c×rand

In the formula: p _{I, j, j-1}=0; If during j=1 k=1,2, Λ, j are arranged; C is the random number penalty coefficient; Rand is that the one-dimensional random on [0,1] interval is counted generating function;

Step 3: according to the n group p that obtains _{I, k, j}, utilize following mathematical formulae to calculate n group Δ u respectively _{I, k, j}And r _{I, k, j}:

$m_i=1+k\frac{{\mathrm{\&delta;}}_{\max }}{h_i}{\left(\frac{E}{p_{i,k,j}}\right)}^{0.8},$

In the formula: m _iIt is the coil of strip tightening coefficient of i bar unit; K is for having considered a correction factor theoretical and that actual value exists deviation to do, k=1,2, Λ, j; δ _MaxBe band steel contacting surface asperity maximum, it comprises microcosmic surface unevenness and meso-and mega-relief; h _iBe thickness before the i bar unit strip coiling; E is the elastic modelling quantity of belt steel material; p _{I, k, j}Be i bar unit inboard interlayer radial compressive stress of k layer when batching the j layer;

${\mathrm{\&epsiv;}}_{\mathrm{ri},k,j-1}=\frac{m_i}{E}[(1-\frac{v}{m_i})p_{i,k,j-1}-\frac{v}{m_i}({\mathrm{\&sigma;}}_{\mathrm{Ti},j-1}+r_{i,k,j-1}\frac{p_{i,k+1,j-1}-p_{i,k,j-1}}{h_{i,k,j-1}})],$

${\mathrm{\&epsiv;}}_{\mathrm{ri},k,j}=\frac{m_i}{E}[(1-\frac{v}{m_i})p_{i,k,j}-\frac{v}{m_i}({\mathrm{\&sigma;}}_{\mathrm{Ti},j}+r_{i,k,j}\frac{p_{i,k+1,j}-p_{i,k,j}}{h_{i,k,j}})],$

In the formula: v is the material Poisson's ratio; σ _{Ti, j}Be the coiling tension of i bar unit when batching the j layer, unit is MPa; r _{I, k, j}Be i bar unit k layer radius corresponding when batching the j layer; h _{I, k, j}Be that i bar unit k layer when batching the j layer batches back thickness; ε _{Ri, k, j}It is the radial deformation of i bar unit k layer band steel when batching the j layer;

Δε _ri，k，j＝ε _ri，k，j-ε _ri，k，j-1，

In the formula: Δ ε _{Ri, k, j}It is the k layer band steel of i bar unit radially deformation variable quantity when batching the j-1 layer when batching the j layer;

Δh _i，k，j＝h _iΔε _ri，k，j，

In the formula: Δ h _{I, k, j}It is the k layer band steel of i bar unit amounts of thickness variation when batching the j-1 layer when batching the j layer;

$u_{0i}=\frac{r_{0i}}{E_c}[\frac{r_{0i}^2+r_{\mathrm{ci}}^2}{r_{0i}^2-r_{\mathrm{ci}}^2}-v_c](p_{i,1,j}-p_{i,1,j-1}),$

In the formula: u _0iBe the displacement that the corresponding reel outer surface of i bar unit takes place when batching the j-1 layer when strip coiling j layer, r _0iBe i the corresponding reel outer radius of bar unit, r _CiBe i the corresponding reel inside radius of bar unit, E _cElasticity modulus of materials for reel;

Δu _i，1，j＝u _0i+Δh _i，1，j，

In the formula: Δ u _{I, l, j}Be that i bar unit ground floor band steel batching the displacement that takes place behind the j layer;

Δu _i，k，j＝Δu _i，k-1，j+Δh _i，k，j，

In the formula: Δ u _{I, k, j}Be that i the bar k of unit layer band steel batching the displacement that takes place behind the j layer;

r _i，k，j＝r _i，k-1，j+h _i，k，j-Δu _i，k，j，

In the formula: r _{I, k, j}Be that i bar k of unit layer band steel batches radius after batching the j layer;

Step 4: with the n group Δ u that obtains in the step 3 _{I, k, j}And r _{I, k, j}The following strip coiling radial stress of substitution computing formula is asked for the n group respectively

$p_{i,k,m}^{*}=\frac{1}{(1-v)\frac{r_{i,k,j}}{E}-\frac{r_{i,k,j}^2}{{\mathrm{Eh}}_{i,k,j}}}[{\mathrm{\&Delta;u}}_{i,k,j}+u_{i,k,j-1}-\frac{r_{i,k,j}}{E}({\mathrm{\&sigma;}}_i+r_{i,k,j}\frac{p_{i,k+1,j}}{h_{i,k,j}})];$

Step 5: with

Judge for the measurement index function and respectively to organize p _{I, k, j}With Between error size; If there is one group of p _{I, k, j}With Make J≤ε set up, then go to step 6; Then produce new n group p if be false according to PSO algorithm particle create-rule _{I, k, j}, the p that each group is new _{I, k, j}Create-rule be:

V _i，k，j＝V _i，k，j+c ₁×rand×(pb _i，k，j-p _i，k，j)+c ₂×rand×(gb _i，k，j-p _i，k，j)，

p _i，k，j＝p _i，k，j+V _i，k，j

In the formula: V _{I, k, j}Be PSO algorithm computation p _{I, k, j}Velocity component, the initial velocity component is one and generates number, k=1,2, Λ, j at random; c ₁And c ₂Acceleration weight for algorithm; Pb _{I, k, j}Corresponding p when obtaining minimum J value when this cycle calculations for this group particle _{I, k, j}Gb _{I, k, j}The p of correspondence when obtaining minimum J value when this cycle calculations for n group particle _{I, k, j},

Generating new n group p _{I, k, j}After go to step 3;

Step 6: if j=N calculates and finishes, N batches the number of plies for the roll coil of strip maximum of confirming according to practical condition here.

6. definite cold-strip steel according to claim 5 batches the method for radius, it is characterized in that in the step 6, if J is not equal to N, then makes j=j+1, and forwards step 2 continuation calculating to.

Images