CN116329297B - Plate shape prediction method based on transverse mechanical property difference of rolled piece - Google Patents

Plate shape prediction method based on transverse mechanical property difference of rolled piece Download PDF

Info

Publication number
CN116329297B
CN116329297B CN202310609005.3A CN202310609005A CN116329297B CN 116329297 B CN116329297 B CN 116329297B CN 202310609005 A CN202310609005 A CN 202310609005A CN 116329297 B CN116329297 B CN 116329297B
Authority
CN
China
Prior art keywords
strip steel
rolling
plate shape
strip
plate
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202310609005.3A
Other languages
Chinese (zh)
Other versions
CN116329297A (en
Inventor
李旭
李晓华
韩月娇
Original Assignee
东北大学
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 东北大学 filed Critical 东北大学
Priority to CN202310609005.3A priority Critical patent/CN116329297B/en
Publication of CN116329297A publication Critical patent/CN116329297A/en
Application granted granted Critical
Publication of CN116329297B publication Critical patent/CN116329297B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B21MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21BROLLING OF METAL
    • B21B37/00Control devices or methods specially adapted for metal-rolling mills or the work produced thereby
    • B21B37/28Control of flatness or profile during rolling of strip, sheets or plates
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B21MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21BROLLING OF METAL
    • B21B2261/00Product parameters
    • B21B2261/02Transverse dimensions
    • B21B2261/04Thickness, gauge
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B21MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21BROLLING OF METAL
    • B21B2261/00Product parameters
    • B21B2261/02Transverse dimensions
    • B21B2261/06Width
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B21MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21BROLLING OF METAL
    • B21B2265/00Forming parameters
    • B21B2265/12Rolling load or rolling pressure; roll force
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Abstract

The invention belongs to the technical field of automatic control in a rolling process, and particularly relates to a plate shape prediction method based on transverse mechanical property differences of rolled pieces, which comprises the following steps: cutting part of strip steel before rolling to perform a stretching experiment, obtaining a stretching curve, and calculating the true yield strength and tangential modulus; acquiring roller parameters, rolling process parameters and strip steel parameters before and after rolling; establishing a plate shape simulation model for strip steel-roller deformation coupling analysis; performing a simulation experiment on the strip steel rolling process by using a plate-shaped simulation model; constructing a regulation and control efficacy coefficient calculation model of the plate-shaped executing mechanism, extracting strip steel length data of each simulation experiment in a stable rolling stage, and calculating a plate-shaped value and a regulation and control efficacy coefficient of each plate-shaped executing mechanism; and extracting strip steel width data of each simulation experiment in a stable rolling stage, establishing a strip steel strip shape curve prediction calculation equation based on the strip steel width data and fitting coefficients of the regulation and control efficacy coefficient curve, and inputting regulation and control values of a strip shape executing mechanism to obtain a corresponding strip shape curve.

Description

Plate shape prediction method based on transverse mechanical property difference of rolled piece
Technical Field
The invention belongs to the technical field of automatic control in a rolling process, and relates to a plate shape prediction method based on transverse mechanical property differences of rolled pieces.
Background
In recent years, the demand of cold-rolled sheet strip products keeps a relatively high stable situation for a long time, and along with the increasing of the level of plate shape control precision, downstream enterprises also put forward higher and higher demands on the geometric shape and dimensional precision of the cold-rolled sheet strip products.
Aiming at the problem of controlling the shape of the cold-rolled strip steel, domestic researchers do some related researches. The invention discloses a method for acquiring the regulating and controlling efficacy coefficient of a plate-shaped executing mechanism of a UCM rolling mill, which is disclosed in the Chinese patent publication No. CN 112916624A. The Chinese journal article '1420 six-roller UCM cold continuous rolling mill shape control performance analysis and application' (heavy machinery, 2016, (05): 10-16.) takes a certain 1420 mm six-roller UCM cold rolling mill as a research object, and establishes a three-dimensional nonlinear simulation model of the six-roller UCM rolling mill roll system deformation by utilizing a Marc large-scale finite element software platform. Coupling analysis is carried out on the elastic plastic deformation and the elastic deformation of the roller under the action of various plate-shaped adjustment amounts through numerical simulation, and the regulation and control effects of plate-shaped regulation and control means such as the work roller bending roller, the middle roller bending roller and the middle roller transverse movement of the six-roller UCM rolling mill on the transverse thickness difference of a strip steel outlet, the transverse displacement of strip steel metal and the bending of a rolling mill roller system are deeply researched, so that theoretical basis and reference basis are provided for the optimization design of the roller shape and the online control of the plate shape.
The deficiencies of the above research mainly have three aspects: (1) The rolled piece material has various assumptions, the mechanical properties of the strip steel are not uniformly distributed along the width direction in the actual production process, and the ideal cuboid rolled piece model is adopted to simulate the difference from the actual rolling state; (2) The accuracy of the plate shape control of the actuating mechanism plate shape regulation and control efficacy coefficient obtained by adopting an ideal cuboid rolled piece model is limited; (3) Finite element experiments consume a lot of computation time and do not give a fast method for predicting plate shape based on finite elements.
Disclosure of Invention
In order to solve the technical problems, the invention aims to provide a plate shape prediction method based on the transverse mechanical property difference of rolled pieces.
The invention provides a plate shape prediction method based on transverse mechanical property differences of rolled pieces, which comprises the following steps:
step 1: cutting part of strip steel before rolling to carry out a stretching experiment, obtaining stretching curves at different positions in the width direction, and calculating the true yield strength and tangential modulus at the corresponding positions;
step 2: acquiring roller parameters, rolling process parameters and strip steel parameters before and after rolling, and preprocessing strip steel thickness data before rolling;
step 3: establishing a plate shape simulation model for strip steel-roller deformation coupling analysis according to the calculation data in the step 1 and the parameters obtained in the step 2;
step 4: making experiment plans aiming at different plate-shaped execution mechanisms, and carrying out simulation experiments on the strip steel rolling process by using a plate-shaped simulation model;
step 5: constructing a regulation and control efficacy coefficient calculation model of the plate-shaped executing mechanism, extracting strip steel length data of each simulation experiment in the stable rolling stage in the step 4, and calculating a plate shape value and a regulation and control efficacy coefficient of each plate-shaped executing mechanism;
step 6: and (3) extracting strip steel width data of each simulation experiment stable rolling stage in the step (4), establishing a strip steel plate shape curve prediction calculation equation based on the strip steel width data and fitting coefficients of the regulation and control efficacy coefficient curve, and inputting regulation and control values of the plate shape execution mechanism to obtain a corresponding plate shape curve.
Further, the step 1 specifically includes:
step 1.1: intercepting part of strip steel before rolling, equally dividing the strip steel into N parts along the width direction of the strip steel, and respectively selecting the middle position of each part of strip steel to carry out a stretching experiment to obtain N nominal stress-strain curves;
step 1.2: the true yield strength and tangential modulus of each strip were calculated from the N nominal stress-strain curves.
Further, the step 1.2 specifically includes:
step 1.2.1: the true yield strength was calculated according to the following formula:
wherein, the liquid crystal display device comprises a liquid crystal display device,σ s is true yield strength;σ 0.2 is the nominal yield strength;εis thatσ 0.2 Corresponding true strain;Eis the elastic modulus;
step 1.2.2: the tangential modulus was calculated according to the following formula:
wherein, the liquid crystal display device comprises a liquid crystal display device,Tis the tangential modulus;σ b is true ultimate strength;ε b is true plastic strain;ε s is true elastic strain;TSis the nominal ultimate strength;ELis elongation.
Further, the step 2 specifically includes:
step 2.1: collecting roller parameters, including: roll diameter, roll length, roll density, roll elastic modulus, and roll poisson ratio;
step 2.2: collecting rolling process parameters, including: friction coefficient, rolling speed, front tension, back tension, rolling reduction, work roll bending force, intermediate roll bending force and intermediate roll lateral displacement;
step 2.3: the method for collecting the parameters of the strip steel before and after rolling comprises the following steps: pre-rolling strip width data, pre-rolling strip thickness data, post-rolling strip thickness data and strip density;
step 2.4: selecting a section of strip steel with good plate shape and 500mm length at the inlet of a frame, selecting 770 measuring points at equal intervals in the width direction of the middle position of the section of strip steel, measuring the thickness of the 770 measuring points, and selecting thickness data of 77 measuring points after removing abnormal value points, wherein the abnormal value points are measuring points with thickness deviation exceeding +/-0.1 mm from left and right measuring points, and the thickness data are used as data required for modeling of transverse thickness distribution of a strip steel model.
Further, the step 3 specifically includes:
step 3.1: establishing a three-dimensional rolling mill model according to the acquired roller data;
step 3.2: establishing a strip steel model according to the strip steel width data before rolling, the strip steel thickness data before rolling after pretreatment and the strip steel thickness data after rolling and the strip steel density;
step 3.3: the strip steel model is divided into N fiber strips along the width direction, and the material properties of the N fiber strips are set according to the real yield strength, the tangential modulus and the elastic modulus of the N strip steel obtained in the step 1.2.
Further, the step 4 specifically includes:
step 4.1: according to the acquired rolling process parameters, an experimental plan for different plate-shaped execution mechanisms is formulated, wherein the plate-shaped execution mechanisms comprise: the middle roller transversely moves, the working roller and the middle roller;
step 4.2: and changing the transverse movement amount of the intermediate roll, the bending force of the working roll and the value of the bending force of the intermediate roll in the plate-shaped simulation model according to an experimental plan, and performing simulation.
Further, the step 5 specifically includes:
step 5.1: establishing a three-dimensional coordinate system taking the width, the length and the thickness of the strip steel as coordinate axes, and taking a central position point of the strip steel as a coordinate origin; the strip steel is assumed to be a continuous discretized longitudinal fiber strip,αthe point is set as a transverse position from the center of the strip steelx α Discrete point numbering of the discretized longitudinal fiber strip of (2) is set to the plate-shaped executing mechanism numberingiThe adjustment amount is deltau i
Step 5.2: any position of rolled strip steel in width directionx α The variation delta of the relative length difference of the longitudinal fiber strip after passing through the strip steel in the plate shape value defined by the plate shape unit IUI α,i (x α ) The calculation gives:
in the method, in the process of the invention,L(x α ) For a selected pre-roll length of the longitudinal fiber strip,l(x α ) The length of the selected longitudinal fiber strip after rolling is used;
step 5.3: the control efficiency coefficient is expressed by the change quantity of the relative length difference of the longitudinal fiber strips after strip steel rolling:
in the method, in the process of the invention,E α,i (x α ) Is the firstiPlate-shaped actuating mechanism is arranged onx α A regulatory efficacy coefficient at the site;
step 5.4: extracting the length data of the strip steel in each simulated experiment stable rolling stage in the step 4, and calculating according to the step 5.2 to obtain the positionx α The variation of the relative length difference of the longitudinal fiber strips after rolling the strip steel, namely the positionx α And (3) calculating the plate shape value of each plate shape executing mechanism according to the step (5.3) to obtain the regulating and controlling efficacy coefficient and the distribution curve of each plate shape executing mechanism.
Further, the step 6 specifically includes:
step 6.1: extracting the strip steel width data of each simulation experiment stable rolling stage in the step 4, and carrying out normalization treatment on the width data;
step 6.2: fitting the regulation and control efficacy coefficient curves of the three plate-shaped execution mechanisms obtained in the step 5.4 by using six-degree Legendre orthogonal polynomials to obtain fitting coefficients of the regulation and control efficacy coefficient curves of the three plate-shaped execution mechanismsA 0A 1A 2A 4A 6 The six-degree Legendre orthopolynomial used is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,E(x) Is a strip steel plate-shaped curve;A 0 is a constant term,A 1A 2A 4A 6 The absolute values of the fitting coefficients of the first, second, fourth and sixth terms of the regulation efficacy coefficient curve show the regulation and control division of the first, second, fourth and sixth plate shape defects in the regulation efficacy coefficient curveAn amount of;xis dimensionless coordinate after normalizing the width direction of the strip steel,x∈[-1,1];e(x) Fitting errors;
step 6.3: normalized strip widthxThe values of the axes and the plate shapes areyThe axis establishes a two-dimensional coordinate system, and a prediction calculation equation of the plate shape value of each coordinate point is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,I P for the calculated plate shape value per coordinate point,I FEM the plate shape value of each coordinate point extracted after the simulation experiment;
wherein, the liquid crystal display device comprises a liquid crystal display device,V WRB for the roll bending force of the work rolls,V IRB for the roll bending force of the intermediate roll,V IRS for the amount of lateral movement of the intermediate roll,A W0A I0 andA S0 are all a constant term and are used to determine,A W1A W2A W4 andA W6 respectively fitting coefficients of a first-order, a second-order, a fourth-order and a sixth-order working roll bending regulation efficacy coefficient curve;A I1A I2A I4 andA I6 respectively fitting coefficients of the curve of the regulating efficacy coefficient of the middle roller for one time, two times, four times and six times;A S1A S2A S4 andA S6 respectively the middleFitting coefficients of the roll traversing regulation efficacy coefficient curves for one time, two times, four times and six times;
step 6.4: and (3) carrying out inverse normalization processing on the normalized width data obtained in the step (6.1), substituting the normalized width data into a predictive calculation equation of the step (6.3), calculating to obtain the plate shape value of each coordinate point, and obtaining the integral plate shape curve.
The plate shape prediction method based on the transverse mechanical property difference of the rolled piece has at least the following beneficial effects:
(1) Firstly, the provided plate shape simulation model for the strip steel-roller deformation coupling analysis considers the overall change condition of the thickness, the cross section shape and the plate shape of the strip steel in the cold rolling process, so that the calculation result is more accurate.
(2) Secondly, the mechanical property parameters of the strip steel model built by the invention are built by actually measuring the tensile stress-strain curve of the strip steel at the inlet of the frame on site, so that the analysis of the cold rolling process under the condition of the difference of the transverse mechanical properties of the rolled piece is realized, the convexity, the plate shape and other simulation precision of the strip steel plate at the outlet are further improved, and the strip steel model is more close to reality.
(3) And thirdly, the method provided by the invention can be used for solving the regulation and control efficacy coefficient of the rolling mill plate-shape actuating mechanism in the cold rolling process, and a mathematical regression method is used for providing an outlet plate-shape curve prediction method for the cold-rolled strip steel with different on-site dimensions in combination with the regulation and control mechanism, so that a basis is provided for rapid regulation of the strip steel plate shape in actual production.
(4) Finally, the invention is based on the rolling theory and the finite element display dynamics simulation, which can effectively reduce the equipment and time loss caused by the experiment and reduce the enterprise cost.
Drawings
FIG. 1 is a flow chart of a strip shape prediction method based on the differences in transverse mechanical properties of rolled pieces according to the present invention;
FIG. 2 is a plot of the tensile curve of a strip measured in the field at a location laterally;
FIG. 3 is a graph of transverse thickness distribution data of strip steel at the entrance of a frame;
FIG. 4 is a graph of the control efficiency coefficient of a work roll;
FIG. 5 is a graph of the control efficacy coefficient of an intermediate roll.
FIG. 6 is a graph of the coefficient of efficacy of the regulation of the traversing of the intermediate roll.
Detailed Description
In this example, a 1800mm UCM six-roller cold rolling mill set of a certain factory is taken as an example, the rolling stability is predicted, and the rolling rolls of the rolling mill are all flat rolls.
As shown in FIG. 1, the plate shape prediction method based on the transverse mechanical property difference of the rolled piece comprises the following steps:
step 1: cutting part of strip steel before rolling to carry out a stretching experiment, obtaining stretching curves at different positions in the width direction, and calculating the true yield strength and tangential modulus at the corresponding positions, wherein the step 1 specifically comprises the following steps:
step 1.1: and cutting out part of the strip steel before rolling, equally dividing the strip steel into 32 parts along the width direction of the strip steel, and respectively selecting the middle position of each part of strip steel to carry out a stretching experiment to obtain 32 nominal stress-strain curves. The stretch curve at a certain position obtained in this example is shown in fig. 2.
Step 1.2: calculating the true yield strength and tangential modulus of each strip steel according to 32 nominal stress-strain curves;
taking the obtained single nominal stress-strain curve as an example: the nominal yield strength, elastic modulus, nominal ultimate strength and elongation data can be obtained from the acquired nominal stress-strain curve, and the actual yield strength and tangential modulus can be calculated using the acquired nominal yield strength, elastic modulus, nominal ultimate strength and elongation data, the step 1.2 is specifically:
step 1.2.1: the true yield strength was calculated according to the following formula:
wherein, the liquid crystal display device comprises a liquid crystal display device,σ s is true yield strength;σ 0.2 is the nominal yield strength;εis thatσ 0.2 Corresponding true strain;Eis the elastic modulus;
step 1.2.2: the tangential modulus was calculated according to the following formula:
wherein, the liquid crystal display device comprises a liquid crystal display device,Tis the tangential modulus;σ b is true ultimate strength;ε b is true plastic strain;ε s is true elastic strain;TSis the nominal ultimate strength;ELis elongation.
In this example, the data obtained by the tensile test are shown in table 1, and the true yield strength and tangential modulus of the steel strip at these 32 positions can be finally calculated from the data of the nominal yield strength, elastic modulus, nominal ultimate strength and elongation of the steel strip at these 32 positions.
TABLE 1 tensile test data
Step 2: the method comprises the steps of obtaining roller parameters, rolling process parameters and strip steel parameters before and after rolling, and preprocessing strip steel thickness data before rolling, wherein the step 2 specifically comprises the following steps:
step 2.1: collecting roller parameters, including: roll diameter, roll length, roll density, roll elastic modulus, and roll poisson ratio;
the method specifically comprises the following parameters: work roll neck diameter, work roll neck length, work roll body diameter, work roll body length, work roll density, work roll elastic modulus, work roll poisson ratio; intermediate roll pair diameter, intermediate roll pair length, intermediate roll neck diameter, intermediate roll neck length, intermediate roll body diameter, intermediate roll body length, intermediate roll density, intermediate roll elastic modulus, intermediate roll poisson ratio; the diameter of the supporting roller neck, the length of the supporting roller neck, the diameter of the supporting roller body, the length of the supporting roller body, the density of the supporting roller, the elastic modulus of the supporting roller and the poisson ratio of the supporting roller. The material properties of the working roll, the intermediate roll and the supporting roll are the same, and the density is 7800 kg/m 3 The modulus of elasticity was 205GPa and the Poisson's ratio was 0.3. The dimensional parameters of the rolls are shown in table 2.
Table 2 roll size parameters
Step 2.2: collecting rolling process parameters, including: friction coefficient, rolling speed, front tension, back tension, rolling reduction, work roll bending force, intermediate roll bending force and intermediate roll lateral displacement. The rolling process parameters are shown in table 3.
TABLE 3 Rolling Process parameters
Step 2.3: the method for collecting the parameters of the strip steel before and after rolling specifically comprises the following steps: the width data of the strip steel before rolling, the thickness data of the strip steel before rolling (inlet thickness), the thickness data of the strip steel after rolling (outlet thickness) and the strip steel density are approximately equal. As shown in table 4, which shows some of the parameters collected, only the thickness data at the center of the inlet of the thickness data of the strip before rolling and the thickness data at the center of the outlet of the thickness data of the strip after rolling are shown in table 4.
TABLE 4 parameters of strip steel
Step 2.4: selecting a section of strip steel with good plate shape and 500mm length at the inlet of a frame, selecting 770 measuring points at equal intervals in the width direction of the middle position of the section of strip steel, measuring the thickness of the 770 measuring points, and selecting thickness data of 77 measuring points after removing abnormal value points, wherein the abnormal value points are measuring points with thickness deviation exceeding +/-0.1 mm from left and right measuring points, and the thickness data are used as data required for modeling of transverse thickness distribution of a strip steel model. The obtained transverse thickness distribution data of the strip steel at the inlet of the frame are shown in fig. 3.
Step 3: establishing a plate shape simulation model for strip steel-roller deformation coupling analysis according to the calculation data in the step 1 and the parameters obtained in the step 2, wherein the step 3 specifically comprises the following steps:
step 3.1: establishing a three-dimensional rolling mill model according to the acquired roller data;
step 3.2: establishing a strip steel model according to the strip steel width data before rolling, the strip steel thickness data before rolling after pretreatment and the strip steel thickness data after rolling and the strip steel density;
step 3.3: the strip steel model was divided into 32 fiber strips along the width direction, and the material properties of the 32 fiber strips were set according to the true yield strength, tangential modulus, and elastic modulus of the 32 strip steel obtained in step 1.2.
Step 4: making experiment plans aiming at different plate-shaped execution mechanisms, and carrying out simulation experiments on the strip steel rolling process by using a plate-shaped simulation model, wherein the step 4 specifically comprises the following steps:
step 4.1: according to the acquired rolling process parameters, an experimental plan for different plate-shaped execution mechanisms is formulated, wherein the plate-shaped execution mechanisms comprise: the middle roller transversely moves, the working roller and the middle roller;
in specific implementation, five experimental value points are respectively set for the regulation and control values of the three plate-shaped execution mechanisms, simulation experiments are respectively carried out, and experimental setting conditions are shown in table 5;
TABLE 5 setting values for work roll bending force, intermediate roll bending force and intermediate roll lateral displacement
Step 4.2: and changing the values of the transverse movement amount of the intermediate roll, the bending force of the working roll and the bending force of the intermediate roll in the plate-shaped simulation model according to the comparison experiment plan, and performing simulation.
Step 5: constructing a regulation and control efficacy coefficient calculation model of the plate-shaped executing mechanism, extracting strip steel length data of each simulation experiment stable rolling stage in the step 4, and calculating a plate shape value and a regulation and control efficacy coefficient of each plate-shaped executing mechanism, wherein the step 5 specifically comprises the following steps:
step 5.1: establishing a three-dimensional coordinate system taking the width, the length and the thickness of the strip steel as coordinate axes, and taking a central position point of the strip steel as a coordinate origin; the strip steel is assumed to be a continuous discretized longitudinal fiber strip,αthe point is set as a transverse position from the center of the strip steelx α Discrete point numbering of the discretized longitudinal fiber strip of (2) is set to the plate-shaped executing mechanism numberingiThe adjustment amount is deltau i
Step 5.2: any position of rolled strip steel in width directionx α The variation delta of the relative length difference of the longitudinal fiber strip after passing through the strip steel in the plate shape value defined by the plate shape unit IUI α,i (x α ) The calculation gives:
in the method, in the process of the invention,L(x α ) For a selected pre-roll length of the longitudinal fiber strip,l(x α ) The length of the selected longitudinal fiber strip after rolling is used;
step 5.3: the control efficiency coefficient is expressed by the change quantity of the relative length difference of the longitudinal fiber strips after strip steel rolling:
in the method, in the process of the invention,E α,i (x α ) Is the firstiPlate-shaped actuating mechanism is arranged onx α A regulatory efficacy coefficient at the site;
step 5.4: extracting the length data of the strip steel in each simulated experiment stable rolling stage in the step 4, and calculating according to the step 5.2 to obtain the positionx α The variation of the relative length difference of the longitudinal fiber strips after rolling the strip steel, namely the positionx α And (3) calculating the plate shape value of each plate shape executing mechanism according to the step (5.3) to obtain the regulating and controlling efficacy coefficient and the distribution curve of each plate shape executing mechanism.
Obtaining the regulation and control efficiency coefficients to obtain the regulation and control efficiency coefficient curves of the three plate-shaped execution mechanisms under different regulation amounts, wherein the regulation and control efficiency coefficient curves are shown in figures 4-6. The curve of the regulating efficiency coefficient of the transverse movement of the working roll bending roll, the middle roll bending roll and the middle roll is in inverted V-shaped distribution, which shows that the plate shape value of the edge part of the strip steel can be reduced and the plate shape value of the middle part can be increased along with the increase of the regulating quantity of the transverse movement of the working roll bending roll, the middle roll bending roll and the middle roll. The larger the slope of the regulation and control efficiency coefficient curve is, the larger the difference value between the strip shape of the middle part and the edge part of the strip steel is, and the stronger the strip shape regulation and control capability is. As can be seen from fig. 4 to 6, the intermediate roll traversing has the strongest adjusting capability for the strip steel plate shape, and the working roll bending and the intermediate roll bending are the weakest, the adjusting and controlling effects of the working roll bending and the intermediate roll bending are not changed greatly, the distribution of the adjusting and controlling effect coefficient curves is not changed basically, the adjusting and controlling effect of the intermediate roll traversing is relatively changed greatly, and the distribution of the adjusting and controlling effect coefficient curves is changed obviously along with the change of the traversing amount.
Step 6: extracting strip steel width data of each simulation experiment stable rolling stage in the step 4, establishing a strip steel plate shape curve prediction calculation equation based on the strip steel width data and fitting coefficients of a regulation and control efficacy coefficient curve, and inputting regulation and control values of a plate shape executing mechanism to obtain a corresponding plate shape curve, wherein the step 6 specifically comprises the following steps:
step 6.1: extracting the strip steel width data of each simulation experiment stable rolling stage in the step 4, and carrying out normalization treatment on the width data;
step 6.2: fitting the regulation and control efficacy coefficient curves of the three plate-shaped execution mechanisms obtained in the step 5.4 by using six-degree Legendre orthogonal polynomials to obtain fitting coefficients of the regulation and control efficacy coefficient curves of the three plate-shaped execution mechanismsA 0A 1A 2A 4A 6 The six-degree Legendre orthopolynomial used is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,E(x) Is a strip steel plate-shaped curve;A 0 is a constant term,A 1A 2A 4A 6 The absolute values of the fitting coefficients are respectively the first, second, fourth and sixth fitting coefficients of the regulation efficacy coefficient curve, and the absolute values of the fitting coefficients indicate the regulation components of the regulation efficacy coefficient curve for the first, second, fourth and sixth plate defects;xis dimensionless coordinate after normalizing the width direction of the strip steel,x∈[-1,1];e(x) Fitting errors;
step 6.3: normalized strip widthxThe values of the axes and the plate shapes areyThe axis establishes a two-dimensional coordinate system, and a prediction calculation equation of the plate shape value of each coordinate point is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,I P for the calculated plate shape value per coordinate point,I FEM the plate shape value of each coordinate point extracted after the simulation experiment;
wherein, the liquid crystal display device comprises a liquid crystal display device,V WRB for the roll bending force of the work rolls,V IRB for the roll bending force of the intermediate roll,V IRS for the amount of lateral movement of the intermediate roll,A W0A I0 andA S0 are all a constant term and are used to determine,A W1A W2A W4 andA W6 respectively fitting coefficients of a first-order, a second-order, a fourth-order and a sixth-order working roll bending regulation efficacy coefficient curve;A I1A I2A I4 andA I6 respectively fitting coefficients of the curve of the regulating efficacy coefficient of the middle roller for one time, two times, four times and six times;A S1A S2A S4 andA S6 respectively fitting coefficients of the intermediate roller traversing regulation efficacy coefficient curves for one time, two times, four times and six times;
step 6.4: and (3) carrying out inverse normalization processing on the normalized width data obtained in the step (6.1), substituting the normalized width data into a predictive calculation equation of the step (6.3), calculating to obtain the plate shape value of each coordinate point, and obtaining the integral plate shape curve. The strip steel plate shape curve prediction model is obtained, and the work roll bending force, the intermediate roll bending force and the intermediate roll lateral movement are all used as input variables in the model.
In the specific implementation, the accuracy verification can be carried out on the strip steel plate shape curve prediction model in the step 6.4. The verification steps are as follows:
firstly, inputting the work roll bending force, the intermediate roll bending force and the intermediate roll lateral movement obtained in the step 2.2 into a model as variables to obtain a calculated plate shape value;
and then, comparing the predicted result with the plate shape value of each coordinate point obtained in the step 5.2, wherein the verification result shows that the model calculation result is well matched with the simulation value of the simulation experiment, has higher prediction precision, and can provide data support for quickly adjusting the adjustment quantity of the plate shape executing mechanism in actual production so as to achieve good plate shape.
The foregoing description of the preferred embodiments of the invention is not intended to limit the scope of the invention, but rather to enable any modification, equivalent replacement, improvement or the like to be made without departing from the spirit and principles of the invention.

Claims (5)

1. A plate shape prediction method based on the transverse mechanical property difference of rolled pieces is characterized by comprising the following steps:
step 1: cutting part of strip steel before rolling to carry out a stretching experiment, obtaining stretching curves at different positions in the width direction, and calculating the true yield strength and tangential modulus at the corresponding positions;
step 2: acquiring roller parameters, rolling process parameters and strip steel parameters before and after rolling, and preprocessing strip steel thickness data before rolling;
step 3: establishing a plate shape simulation model for strip steel-roller deformation coupling analysis according to the calculation data in the step 1 and the parameters obtained in the step 2;
step 4: making experiment plans aiming at different plate-shaped execution mechanisms, and carrying out simulation experiments on the strip steel rolling process by using a plate-shaped simulation model;
step 5: constructing a regulation and control efficacy coefficient calculation model of the plate-shaped executing mechanism, extracting strip steel length data of each simulation experiment in the stable rolling stage in the step 4, and calculating a plate shape value and a regulation and control efficacy coefficient of each plate-shaped executing mechanism;
step 6: extracting strip steel width data of each simulation experiment stable rolling stage in the step 4, establishing a strip steel plate shape curve prediction calculation equation based on the strip steel width data and fitting coefficients of a regulation and control efficacy coefficient curve, and inputting regulation and control values of a plate shape executing mechanism to obtain a corresponding plate shape curve;
the step 1 specifically comprises the following steps:
step 1.1: intercepting part of strip steel before rolling, equally dividing the strip steel into N parts along the width direction of the strip steel, and respectively selecting the middle position of each part of strip steel to carry out a stretching experiment to obtain N nominal stress-strain curves;
step 1.2: calculating the true yield strength and tangential modulus of each strip steel according to the N nominal stress-strain curves;
the step 2 specifically comprises the following steps:
step 2.1: collecting roller parameters, including: roll diameter, roll length, roll density, roll elastic modulus, and roll poisson ratio;
the method specifically comprises the following parameters: work roll neck diameter, work roll neck length, work roll body diameter, work roll body length, work roll density, work roll elastic modulus, work roll poisson ratio; intermediate roll pair diameter, intermediate roll pair length, intermediate roll neck diameter, intermediate roll neck length, intermediate roll body diameter, intermediate roll body length, intermediate roll density, intermediate roll elastic modulus, intermediate roll poisson ratio; the diameter of the support roller neck, the length of the support roller neck, the diameter of the support roller body, the length of the support roller body, the density of the support roller, the elastic modulus of the support roller and the poisson ratio of the support roller;
step 2.2: collecting rolling process parameters, including: friction coefficient, rolling speed, front tension, back tension, rolling reduction, work roll bending force, intermediate roll bending force and intermediate roll lateral displacement;
step 2.3: the method for collecting the parameters of the strip steel before and after rolling comprises the following steps: pre-rolling strip width data, pre-rolling strip thickness data, post-rolling strip thickness data and strip density;
step 2.4: selecting a section of strip steel with good plate shape and 500mm length at the inlet of a frame, selecting 770 measuring points at equal intervals in the width direction of the middle position of the section of strip steel, measuring the thickness of the 770 measuring points, and selecting thickness data of 77 measuring points after removing abnormal value points, wherein the abnormal value points are measuring points with thickness deviation exceeding +/-0.1 mm from left and right measuring points, and the thickness data are used as data required by modeling of transverse thickness distribution of a strip steel model;
the step 3 specifically comprises the following steps:
step 3.1: establishing a three-dimensional rolling mill model according to the acquired roller data;
step 3.2: establishing a strip steel model according to the strip steel width data before rolling, the strip steel thickness data before rolling after pretreatment and the strip steel thickness data after rolling and the strip steel density;
step 3.3: the strip steel model is divided into N fiber strips along the width direction, and the material properties of the N fiber strips are set according to the real yield strength, the tangential modulus and the elastic modulus of the N strip steel obtained in the step 1.2.
2. The method for predicting plate shape based on the difference of transverse mechanical properties of rolled pieces according to claim 1, wherein the step 1.2 specifically comprises:
step 1.2.1: the true yield strength was calculated according to the following formula:
σ s =σ 0.2 ×(1+ε);
ε=0.002+σ 0.2 /E;
wherein sigma s Is true yield strength; sigma (sigma) 0.2 Is the nominal yield strength; epsilon is sigma 0.2 Corresponding true strain; e is the elastic modulus;
step 1.2.2: the tangential modulus was calculated according to the following formula:
T=(σ bs )/(ε bs );
σ b =TS×(1+EL);
ε b =ln(l+EL);
ε s =σ s /E;
wherein T is the tangential modulus; sigma (sigma) b Is true ultimate strength; epsilon b Is true plastic strain; epsilon s Is true elastic strain; TS is the nominal ultimate strength; EL is elongation.
3. The method for predicting plate shape based on the difference of transverse mechanical properties of rolled piece according to claim 1, wherein the step 4 specifically comprises:
step 4.1: according to the acquired rolling process parameters, an experimental plan for different plate-shaped execution mechanisms is formulated, wherein the plate-shaped execution mechanisms comprise: the middle roller transversely moves, the working roller and the middle roller;
step 4.2: and changing the transverse movement amount of the intermediate roll, the bending force of the working roll and the value of the bending force of the intermediate roll in the plate-shaped simulation model according to an experimental plan, and performing simulation.
4. The method for predicting plate shape based on the difference of transverse mechanical properties of rolled piece according to claim 1, wherein the step 5 specifically comprises:
step 5.1: establishing a three-dimensional coordinate system taking the width, the length and the thickness of the strip steel as coordinate axes, and taking a central position point of the strip steel as a coordinate origin; assuming the strip steel as a continuous discretized longitudinal fiber strip, the alpha point is set as the transverse position x from the center of the strip steel α The discrete point number of the discretized longitudinal fiber strip is set as i, and the adjustment quantity is deltau i
Step 5.2: any position x in width direction of rolled strip steel α The variation delta I of the relative length difference of the longitudinal fiber strip after passing through the strip steel in the plate shape value defined by the plate shape unit IU α,i (x α ) The calculation gives:
wherein L (x) α ) For the pre-roll length of the selected longitudinal fiber strand, l (x α ) The length of the selected longitudinal fiber strip after rolling is used;
step 5.3: the control efficiency coefficient is expressed by the change quantity of the relative length difference of the longitudinal fiber strips after strip steel rolling:
wherein E is α,i (x α ) For the ith plate-shaped actuator at x α A regulatory efficacy coefficient at the site;
step 5.4: extracting strip steel length data of each simulated experiment stable rolling stage in the step 4, and calculating according to the step 5.2 to obtain a position x α The variation of the relative length difference of the longitudinal fiber strips after rolling the strip steel, namely the position x α And (3) calculating the plate shape value of each plate shape executing mechanism according to the step (5.3) to obtain the regulating and controlling efficacy coefficient and the distribution curve of each plate shape executing mechanism.
5. The method for predicting plate shape based on the difference of transverse mechanical properties of rolled piece according to claim 4, wherein the step 6 specifically comprises:
step 6.1: extracting the strip steel width data of each simulation experiment stable rolling stage in the step 4, and carrying out normalization treatment on the width data;
step 6.2: fitting the regulation and control efficacy coefficient curves of the three plate-shaped execution mechanisms obtained in the step 5.4 by using six-degree Legendre orthogonal polynomials to obtain fitting coefficients A of the regulation and control efficacy coefficient curves of the three plate-shaped execution mechanisms 0 、A 1 、A 2 、A 4 、A 6 The six-degree Legendre orthopolynomial used is as follows:
wherein E (x) is a strip steel plate curve; a is that 0 Is a constant term, A 1 、A 2 、A 4 、A 6 The absolute values of the fitting coefficients are respectively the first, second, fourth and sixth fitting coefficients of the regulation efficacy coefficient curve, and the absolute values of the fitting coefficients indicate the regulation components of the regulation efficacy coefficient curve for the first, second, fourth and sixth plate defects; x is dimensionless coordinate after normalizing the width direction of strip steel, x is [ -1,1 []The method comprises the steps of carrying out a first treatment on the surface of the e (x) is the fitting error;
step 6.3: establishing a two-dimensional coordinate system by taking the normalized strip steel width as an x-axis and the strip shape value as a y-axis, wherein a prediction calculation equation of the strip shape value of each coordinate point is as follows:
I P =I FEM +eff WRB +eff IRB +eff IRS
wherein I is P To calculate the plate shape value of each coordinate point, I FEM The plate shape value of each coordinate point extracted after the simulation experiment;
eff WRB =V WRB ×{[A 0W +A 1W +A 2W ×lgd 2 +A 4W ×lgd 4 +A 6W ×lgd 6 +e(x)]};
eff IRB =V IRB ×{[A 0I +A 1I +A 2I ×lgd 2 +A 4I ×lgd 4 +A 6I ×lgd 6 +e(x)]};
eff IRS =V IRS ×{[A 0S +A 1S +A 2S ×lgd 2 +A 4S ×lgd 4 +A 6S ×lgd 6 +e(x)]};
wherein V is WRB For the work roll bending force, V IRB For the intermediate roll bending force, V IRS For the intermediate roll lateral movement, A 0W 、A 0I And A 0S Are all constant terms, A 1W 、A 2W 、A 4W And A 6W Respectively fitting coefficients of a first-order, a second-order, a fourth-order and a sixth-order working roll bending regulation efficacy coefficient curve; a is that 1I 、A 2I 、A 4I And A 6I Respectively regulating and controlling efficacy coefficient curves of the middle roller, namely one time, two times, four times and six timesA minor term fitting coefficient; a is that 1S 、A 2S 、A 4S And A 6S Respectively fitting coefficients of the intermediate roller traversing regulation efficacy coefficient curves for one time, two times, four times and six times;
step 6.4: and (3) carrying out inverse normalization processing on the normalized width data obtained in the step (6.1), substituting the normalized width data into a predictive calculation equation of the step (6.3), calculating to obtain the plate shape value of each coordinate point, and obtaining the integral plate shape curve.
CN202310609005.3A 2023-05-29 2023-05-29 Plate shape prediction method based on transverse mechanical property difference of rolled piece Active CN116329297B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310609005.3A CN116329297B (en) 2023-05-29 2023-05-29 Plate shape prediction method based on transverse mechanical property difference of rolled piece

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310609005.3A CN116329297B (en) 2023-05-29 2023-05-29 Plate shape prediction method based on transverse mechanical property difference of rolled piece

Publications (2)

Publication Number Publication Date
CN116329297A CN116329297A (en) 2023-06-27
CN116329297B true CN116329297B (en) 2023-10-20

Family

ID=86884395

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310609005.3A Active CN116329297B (en) 2023-05-29 2023-05-29 Plate shape prediction method based on transverse mechanical property difference of rolled piece

Country Status (1)

Country Link
CN (1) CN116329297B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117019884B (en) * 2023-10-08 2023-12-29 东北大学 Cold continuous rolling machine frame outlet plate shape prediction visualization method

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1243047A (en) * 1999-08-27 2000-02-02 冶金工业部钢铁研究总院 Interlinked control method for plate-band rolling course based on coordination law of plate shape and plate thickness
JP2000158025A (en) * 1998-11-20 2000-06-13 Hitachi Ltd Rolling method, rolling mill and rolling equipment
CN101716607A (en) * 2009-12-17 2010-06-02 燕山大学 Method for controlling asymmetric transverses shifting plate shape of asymmetric bending roller of HC rolling mill
CN102601127A (en) * 2012-03-19 2012-07-25 中冶南方工程技术有限公司 High-precision strip shape control prediction method for CVC (continuously variable crown) four-roll cold rolling mill
CN102632085A (en) * 2012-04-23 2012-08-15 中冶南方工程技术有限公司 Cold-rolled strip steel plate shape control system and method
CN102641896A (en) * 2012-04-19 2012-08-22 中冶南方工程技术有限公司 Gauge and flatness comprehensive control system of cold rolled steel sheet
JP2015158443A (en) * 2014-02-25 2015-09-03 Jfeスチール株式会社 Thickness measuring method and thickness measuring apparatus of checkered steel plate
CN110947774A (en) * 2019-12-06 2020-04-03 东北大学 Plate shape prediction method considering rolling width
CN115026136A (en) * 2022-08-11 2022-09-09 东北大学 Method for predicting plate shape of endless rolling wedge-shaped transition process

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000158025A (en) * 1998-11-20 2000-06-13 Hitachi Ltd Rolling method, rolling mill and rolling equipment
CN1243047A (en) * 1999-08-27 2000-02-02 冶金工业部钢铁研究总院 Interlinked control method for plate-band rolling course based on coordination law of plate shape and plate thickness
CN101716607A (en) * 2009-12-17 2010-06-02 燕山大学 Method for controlling asymmetric transverses shifting plate shape of asymmetric bending roller of HC rolling mill
CN102601127A (en) * 2012-03-19 2012-07-25 中冶南方工程技术有限公司 High-precision strip shape control prediction method for CVC (continuously variable crown) four-roll cold rolling mill
CN102641896A (en) * 2012-04-19 2012-08-22 中冶南方工程技术有限公司 Gauge and flatness comprehensive control system of cold rolled steel sheet
CN102632085A (en) * 2012-04-23 2012-08-15 中冶南方工程技术有限公司 Cold-rolled strip steel plate shape control system and method
JP2015158443A (en) * 2014-02-25 2015-09-03 Jfeスチール株式会社 Thickness measuring method and thickness measuring apparatus of checkered steel plate
CN110947774A (en) * 2019-12-06 2020-04-03 东北大学 Plate shape prediction method considering rolling width
CN115026136A (en) * 2022-08-11 2022-09-09 东北大学 Method for predicting plate shape of endless rolling wedge-shaped transition process

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于弹塑性有限元的板形控制机理研究现状与展望;李旭;轧钢;第37卷(第4期);1-11 *

Also Published As

Publication number Publication date
CN116329297A (en) 2023-06-27

Similar Documents

Publication Publication Date Title
CN101623708B (en) Plate-shape control integrated system and executing method thereof
CN116329297B (en) Plate shape prediction method based on transverse mechanical property difference of rolled piece
CN102029294B (en) Control method of lateral thick difference of cold-rolling strip steels
Wang et al. Numerical and experimental analysis of strip cross-directional control and flatness prediction for UCM cold rolling mill
CN110773573B (en) Plate-shaped regulation and control efficiency coefficient actual measurement data processing method
CN102601127A (en) High-precision strip shape control prediction method for CVC (continuously variable crown) four-roll cold rolling mill
Behrens et al. Development of an analytical 3D-simulation model of the levelling process
CN104785543B (en) A kind of hot-strip crown feedback control method based on moving average filter
CN112949108B (en) Hot-rolled high-strength steel plate shape defect full-flow prediction method and graphical user interface device
CN101716604A (en) Forecasting method of rolling pressure and roll gap pressure of DC rolling mill
CN104942019A (en) Automatic control method for width of steel strips during cold rolling
CN104451118B (en) A kind of belt plate shape being suitable for continuous annealing process develops forecasting procedure piecemeal
CN112916624A (en) Method for obtaining regulation and control efficiency coefficient of plate-shaped execution mechanism of UCM rolling mill
CN116371942B (en) Transverse thickness distribution prediction method based on uneven transverse strength of strip steel
CN102179414A (en) Computing method of rigidity of novel rolling mill
Du et al. Optimization of short stroke control preset for automatic width control of hot rolling mill
CN113275387B (en) Method for acquiring transverse and longitudinal stiffness characteristic curve of UCM rolling mill roll system
CN102688896A (en) Mechanism model-based virtual crown gauge setting method for four-roll cold continuous rolling machine set
CN113290061B (en) Zirconium alloy-based multi-rolling-process multi-pass deformation system design method
Liu et al. Axial force analysis and roll contour configuration of four-high CVC mill
CN1179369A (en) Plate shape measuring and controlling method for plate strip rolling process
CN114798755B (en) Method for analyzing influence of hot crown on plate shape in cold continuous rolling process
Kukhta et al. Modelling the barrel body wear of the backup roll: mathematical model and software implementation
CN117019884B (en) Cold continuous rolling machine frame outlet plate shape prediction visualization method
Shin et al. Development and experimental evaluation of strip shape prediction model for Sendzimir rolling mills

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant