Emulsion flow optimization method for suppressing vibration of cold continuous rolling mill

Abstract

An emulsion flow optimization method suitable for a cold continuous rolling mill that aims to achieve vibration suppression. The method aims to suppress vibrations by an oil film thickness model and a friction coefficient model. An optimum set value of the emulsion flow rate for each rolling stand that aims to achieve vibration suppression is optimized on the basis of an over-lubrication film thickness critical value and an under-lubrication film thickness critical value that are proposed. The method greatly reduces the incidence of rolling mill vibration defects, improves production efficiency and product quality, treats rolling mill vibration defects, and improves the surface quality and rolling process stability of a finished strip of a cold continuous rolling mill.

Claims

1. An emulsion flow optimization method for suppressing vibration of a cold continuous rolling mill, comprising the following steps: (S1) collecting device feature parameters of the cold continuous rolling mill, wherein the device feature parameters comprise: a radius R.sub.i of a working roll of each of a plurality of rolling stands, a surface linear velocity v.sub.ri of a roll of each rolling stand, a original roughness Ra.sub.ir0 of a working roll of each rolling stand, a roughness attenuation coefficient B.sub.L of a working roll, a distance l between rolling stands, and a rolling kilometer L.sub.i after roll change of a working roll of each rolling stand, wherein i is 1, 2, . . . , n, and represents an ordinal number of the rolling stands of the cold continuous rolling mill, and n is the total number of rolling stands; (S2) collecting key rolling process parameters of a strip, wherein the key rolling process parameters comprise: an inlet thickness h.sub.0i of each rolling stand, an outlet thickness h.sub.1i of each rolling stand, a strip width B, an inlet speed v.sub.0i of each rolling stand, an outlet speed v.sub.1i of each rolling stand, an inlet temperature T.sub.1.sup.r, a strip deformation resistance K.sub.i of each rolling stand, a rolling pressure P.sub.i of each rolling stand, a back tension T.sub.0i of each rolling stand, a front tension T.sub.1i of each rolling stand, an emulsion concentration influence coefficient k.sub.c, a pressure-viscosity coefficient of a lubricant, a strip density , a specific heat capacity S of a strip, an emulsion concentration C, an emulsion temperature T.sub.c and a thermal-work equivalent J; (S3) defining process parameters involved in an emulsion flow optimization process, wherein the process parameters comprise an over-lubrication film thickness critical value .sub.i.sup.+ of each rolling stand, a first friction coefficient u.sub.i.sup.+ at this time, an under-lubrication film thickness critical value .sub.i.sup. and a second friction coefficient u.sub.i.sup. at this time, a rolling reduction amount h.sub.i (wherein h.sub.i=h.sub.0ih.sub.1i), a rolling reduction rate .sub.i (wherein ${}_i=\frac{h_i}{h_{0i}}$ ), an inlet temperature T.sub.i.sup.r of each rolling stand, an over-lubrication judgment coefficient A.sup.+, and an under-lubrication judgment coefficient A.sup., and evenly dividing the distance l between the rolling stands into m sections, wherein a temperature in the sections is represented by T.sub.i,j (wherein 1jm, and T.sub.i.sup.r=T.sub.i1,m); (S4) setting an initial set value of an emulsion flow rate comprehensive optimization objective function of the cold continuous rolling mill for achieving a vibration suppression as F.sub.0=1.010.sup.10; wherein an executing order of steps S1-S4 is not limited; (S5) calculating a bite angle .sub.i of each rolling stand according to a rolling theory, wherein a calculation formula is as follows: ${}_i=\sqrt{\frac{h_i}{R_i^{}},}{R}_i^{}$ is a flattening radius of a working roll of an i.sup.th rolling stand, and is a calculation process value of rolling pressure; (S6) calculating a vibration determination index reference value .sub.0i of each rolling stand; (S7) setting an emulsion flow rate w.sub.i of each rolling stand; (S8) calculating a strip outlet temperature T.sub.i of each rolling stand; (S9) calculating an emulsion flow rate comprehensive optimization objective function F(X); $\{\begin{array}{c}F\left(X\right)=\frac{}{n}\underset{i=1}{\overset{n}{.Math.}}\sqrt{{\left({}_i-{}_{0i}\right)}^2}+\left(1-\right)\max .Math.{}_i-{}_{0i}.Math.\\ {}_i^{-}<{}_i<{}_i^{+}\end{array};$ (S10) determining whether an in-equation F(X)<F.sub.0 is established, if yes, enabling w.sub.i.sup.y=w.sub.i, F.sub.0=F(X), and turning to step S11; otherwise, directly turning to step S11; (S11) determining whether an emulsion flow rate w.sub.i exceeds a feasible region range, if yes, turning to step S12, otherwise, turning to step S7, wherein a feasible region of w.sub.i ranges from 0 to a maximum emulsion flow rate value allowed by the rolling mill; and (S12) outputting an optimum emulsion flow rate set value w.sub.i.sup.y, wherein w.sub.i.sup.y is the value of w.sub.i when a calculated value of F(X) in the feasible region is minimum.

2. The emulsion flow optimization method for suppressing vibration of a cold continuous rolling mill according to claim 1, wherein the step S6 comprises the following steps: (S6.1) calculating a neutral angle .sub.i of each rolling stand: ${}_i=\frac{1}{2}\sqrt{\frac{h_i}{R_i^{}}}\left[1+\frac{1}{2u_i}\left(\sqrt{\frac{h_i}{R_i^{}}}+\frac{T_{i0}-T_{i1}}{P_i}\right)\right];$ (S6.2) calculating to obtain $u_i^{+}=\frac{1}{2\left(2A^{+}-1\right)}\left(\sqrt{\frac{h_i}{R_i^{}}}+\frac{T_{i0}-T_{i1}}{P_i}\right)$ from the step S5 and the step S6.1 assuming that when $\frac{{}_i}{{}_i}=A^{+},$ a roll gap is just in an over-lubrication state; (S6.3) calculating an over-lubrication film thickness critical value .sub.i.sup.+ of each rolling stand according to a relation formula between the first friction coefficient u.sub.i.sup.+ and an oil film thickness wherein u.sub.i.sup.+=a.sub.i+b.sub.i.Math.e.sup.B.sup.i.sup..Math..sup.i.sup.+ in the formula, a.sub.i is a liquid friction influence coefficient, b.sub.i is a dry friction influence coefficient, and B.sub.i is a friction coefficient attenuation index, and wherein ${}_i^{+}=\frac{1}{B_i}\ln \frac{u_i^{+}-a_i}{b_i};$ (S6.4) calculating to obtain $u_i^{-}=\frac{1}{2\left(2A^{-}-1\right)}\left(\sqrt{\frac{h_i}{R_i^{}}}+\frac{T_{i0}-T_{i1}}{P_i}\right)$ from the step S5 and the step S6.1 assuming that when $\frac{{}_i}{{}_i}=A^{-},$ a roll gap is just in an under-lubrication state; (S6.5) calculating an under-lubrication film thickness critical value .sub.i.sup. of each rolling stand according to a relation formula between the second friction coefficient u.sub.i.sup. and an oil film thickness wherein u.sub.i.sup.=a.sub.i+b.sub.i.Math.e.sup.B.sup.i.sup..Math..sup.i.sup., and wherein ${}_i^{-}=\frac{1}{B_i}\ln \frac{u_i^{-}-a_i}{b_i};$ and (S6.6) calculating a vibration determination index reference value .sub.0i, wherein ${}_{0i}=\frac{{}_i^{+}+{}_i^{-}}{2}.$

3. The emulsion flow optimization method for suppressing vibration of a cold continuous rolling mill according to claim 2, wherein the step S8 comprises the following steps: (S8.1) calculating an outlet temperature T.sub.1 of a first rolling stand of the plurality of rolling stands, wherein $T_1=T_1^r+\frac{1-\left({}_1/4\right)}{1-\left({}_1/2\right)}.Math.\frac{K_1\ln \left(\frac{1}{1-{}_1}\right)}{SJ};$ (S8.2) enabling i=1; (S8.3) calculating a temperature T.sub.i,1 of a first section of strip behind an outlet of the i.sup.th rolling stand, i.e. T.sub.i,1=T.sub.i; (S8.4) enabling j=2; (S8.5) calculating a temperature T.sub.i,j of a j.sup.th section of strip by a relationship between a temperature of the j.sup.th section and a temperature of a j1.sup.th section shown by the following equation : $T_{i,j}=-\frac{2k_0{w}_i^{0.264}\exp \left(9.45-0.1918C\right)1.163l}{v_{1i}{h}_{1i}\mathrm{Sm}}{T}_{i,j-1}^{0.213}\left(T_{i,j-1}-T_c\right)+T_{i,j-1},$ wherein k.sub.0 is an influence coefficient of nozzle shape and spraying angle; (S8.6) determining whether an in-equation j<m is established, if yes, enabling j=j+1, and then turning to step S8.5; otherwise, turning to step S8.7; (S8.7) obtaining a temperature T.sub.i,m of a m.sup.th section by iterative calculation; (S8.8) calculating an inlet temperature T.sub.i+1.sup.r of an i+1.sup.th rolling stand: T.sub.i+1.sup.r=T.sub.i,m; (58.9) calculating an outlet temperature T.sub.i+1 of the i+1.sup.th rolling stand, wherein $T_{i+1}=T_{i+1}^r+\frac{1-\left({}_{i+1}/4\right)}{1-\left({}_{i+1}/2\right)}.Math.\frac{K_{i+1}\ln \left(\frac{1}{1-{}_{i+1}}\right)}{\mathrm{SJ}};$ (S8.10) determining whether the in-equation i<n is established, if yes, enabling i=i+1, and then turning to step S8.3; otherwise, turning to step S8.11; and (S8.11) obtaining an outlet temperature T.sub.i of each rolling stand.

4. The emulsion flow optimization method for suppressing vibration of a cold continuous rolling mill according to claim 3, wherein the step S9 comprises the following steps: (S9.1) calculating a dynamic viscosity .sub.0i of an emulsion between a roll gap of each of the plurality of rolling stands, wherein .sub.0i=b.Math.exp(a.Math.T.sub.i), and in the formula, a,b are dynamic viscosity parameters of lubricating oil under an atmospheric pressure; (S9.2) calculating an oil film thickness .sub.i between the roll gap of each of the plurality of rolling stands, wherein the calculation formula is as follows: ${}_i=\frac{h_{0i}+h_{1i}}{2h_{0i}}.Math.k_c.Math.\frac{3{}_{0i}\left(v_{\mathrm{ri}}+v_{0i}\right)}{{}_i\left[1-e^{-\left(K\frac{T_{0i}}{h_{0i}.Math.B}\right)}\right]}-k_{rg}.Math.\left(1+K_{rs}\right).Math.{\mathrm{Ra}}_{\mathrm{ir}0}.Math.e^{-B_L.Math.L_i}$ in the formula, k.sub.rg represents a coefficient of the strength of entrainment of lubricant by a longitudinal surface roughness of a work roll and a strip steel and is in a range of 0.09-0.15, K.sub.rs represents an impression rate, wherein a ratio of transferring a surface roughness of the working roll to the strip; and (S9.3) calculating an emulsion flow rate comprehensive optimization objective function, $\{\begin{array}{c}F\left(X\right)=\frac{}{n}\underset{i=1}{\overset{n}{.Math.}}\sqrt{{\left({}_i-{}_{0i}\right)}^2}+\left(1-\right)\max .Math.{}_i-{}_{0i}.Math.\\ {}_i^{-}<{}_i<{}_i^{+}\end{array}$ in the formula, X={w.sub.i} is an optimization variable, and is a distribution coefficient.

5. The emulsion flow optimization method for suppressing vibration of a cold continuous rolling mill according to claim 3, wherein the influence coefficient of nozzle shape and spraying angle is equal to 0.8<k.sub.0<1.2.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) In the present invention, the same reference numerals always represent the same features, wherein:

(2) FIG. 1 is a flowchart of an emulsion flow optimization method of the present invention;

(3) FIG. 2 is a flowchart of calculating the vibration determination index reference value;

(4) FIG. 3 is a flowchart of calculating the strip outlet temperature of each rolling stand; and

(5) FIG. 4 is a flowchart of calculating an emulsion flow comprehensive optimization objective function.

DETAILED DESCRIPTION

(6) The technical solution of the present invention will be further described in combination with the drawings and the embodiments.

(7) Rolling mill vibration defects are very easily caused between roll gaps of each rolling stand of a cold continuous rolling mill, whether in an over-lubrication state or in an under-lubrication state, and the setting of the emulsion flow rate directly affects the lubrication state between the roll gaps of each rolling stand. In order to realize the treatment of the rolling mill vibration defects, starting from the emulsion flow rate, this patent ensures that both the overall lubrication state of the cold continuous rolling mill and the lubrication state of individual rolling stands can be optimum through the comprehensive optimal distribution of the emulsion flow rate of the cold continuous rolling mill, so as to achieve the goal of treating the rolling mill vibration defects, improving the surface quality and rolling process stability of a finished strip of the cold continuous rolling mill.

(8) Referring to FIG. 1, an emulsion flow optimization method for suppressing vibration of a cold continuous rolling mill includes the following steps: S1, collecting device feature parameters of the cold continuous rolling mill, wherein the device feature parameters include: the radius R.sub.i of a working roll of each rolling stand, the surface linear velocity v.sub.ri of a roll of each rolling stand, the original roughness Ra.sub.ir0 of a working roll of each rolling stand, the roughness attenuation coefficient B.sub.L of a working roll, the distance l between rolling stands, and the rolling kilometer L.sub.i after roll change of a working roll of each rolling stand, wherein i is 1, 2, . . . , n, and represents for the ordinal number of rolling stands of the cold continuous rolling mill, and n is the total number of rolling stands; S2, collecting key rolling process parameters of a strip, wherein the key rolling process parameters include: the inlet thickness h.sub.0i of each rolling stand, the outlet thickness h.sub.1i of each rolling stand, strip width B, the inlet speed v.sub.0i of each rolling stand, the outlet speed v.sub.1i of each rolling stand, the inlet temperature T.sub.1.sup.r, strip deformation resistance K.sub.i of each rolling stand, rolling pressure P.sub.i of each rolling stand, back tension T.sub.0i of each rolling stand, front tension T.sub.1i of each rolling stand, emulsion concentration influence coefficient k.sub.c, pressure-viscosity coefficient of a lubricant, strip density , specific heat capacity S of a strip, emulsion concentration C, emulsion temperature T.sub.c and thermal-work equivalent J; S3, defining process parameters involved in the process of emulsion flow optimization, wherein the process parameters include that an over-lubrication film thickness critical value of each rolling stand is .sub.i.sup.+ and the friction coefficient at this time is u.sub.i.sup.+, an under-lubrication film thickness critical value is .sub.i.sup.+ and the friction coefficient at this time is u.sub.i.sup., the rolling reduction amount is h.sub.i=h.sub.0ih.sub.1i, the rolling reduction rate is

(9) ${}_i=\frac{h_i}{h_{0i}},$
the inlet temperature of each rolling stand is T.sub.i.sup.r, the distance l between the rolling stands is evenly divided into m sections, and the temperature in the sections is represented by T.sub.i,j (wherein, 1jm), and T.sub.i.sup.r=T.sub.i1,m, the over-lubrication judgment coefficient is A.sup.+, and the under-lubrication judgment coefficient is A.sup.; S4, setting the initial set value of an emulsion flow rate comprehensive optimization objective function of the cold continuous rolling mill that aims to achieve vibration suppression as F.sub.0=1.010.sup.10;
the executing order of steps S1 to S4 is not limited, and in some cases, steps S1 to S4 can be performed simultaneously. S5, calculating the bite angle .sub.i of each rolling stand according to the rolling theory, wherein the calculation formula is as follows:

(10) ${}_i=\sqrt{\frac{h_i}{R_i^{}}},R_i^{}$
is the flattening radius of the working roll of the i.sup.th rolling stand, and is the calculation process value of rolling pressure; S6, calculating the vibration determination index reference value .sub.0i of each rolling stand, wherein the calculation flowchart is shown in FIG. 2; S6.1, calculating the neutral angle .sub.i of each rolling stand:

(11) ${}_i=\frac{1}{2}\sqrt{\frac{h_i}{R_i^{}}}\left[1+\frac{1}{2u_i}\left(\sqrt{\frac{h_i}{R_i^{}}}+\frac{T_{i0}-T_{i1}}{P_i}\right)\right];$ S6.2, calculating to obtain

(12) 0$u_i^{+}=\frac{1}{2\left(2A^{+}-1\right)}\left(\sqrt{\frac{h_i}{R_i^{}}}+\frac{T_{i0}-T_{i1}}{P_i}\right)$
from steps S5 and S6.1 assuming that when

(13) $\frac{{}_i}{{}_i}=A^{+},$
the roll gap is just in an over-lubrication state; S6.3, calculating an over-lubrication film thickness critical value .sub.i.sup.+ of each rolling stand according to the relationship formula between the friction coefficient and the oil film thickness, namely u.sub.i=a.sub.i+b.sub.i.Math.e.sup.B.sup.i.sup..Math..sup.i (in the formula, a.sub.i is the liquid friction influence coefficient, b.sub.i is the dry friction influence coefficient, and B.sub.i is the friction coefficient attenuation index), wherein

(14) ${}_i^{+}=\frac{1}{B_i}\ln \frac{u_i^{+}-a_i}{b_i};$ S6.4, calculating to obtain

(15) $u_i^{-}=\frac{1}{2\left(2A^{-}-1\right)}\left(\sqrt{\frac{h_i}{R_i^{}}}+\frac{T_{i0}-T_{i1}}{P_i}\right)$
from steps S5 and S6.1 assuming that when

(16) $\frac{{}_i}{{}_i}=A^{-},$
the roll gap is just in an under-lubrication state; S6.5, calculating an under-lubrication film thickness critical value .sub.i.sup. of each rolling stand according to the relationship formula between the friction coefficient and the oil film thickness, namely u.sub.i=a.sub.i+b.sub.i.Math.e.sup.B.sup.i.sup..Math..sup.i, wherein

(17) ${}_i^{-}=\frac{1}{B_i}\ln \frac{u_i^{-}-a_i}{b_i};$
and S6.6, calculating the vibration determination index reference value .sub.0i of each rolling stand, wherein

(18) ${}_{0i}=\frac{{}_i^{+}+{}_i^{-}}{2};$ S7, setting the emulsion flow rate w.sub.i of each rolling stand; S8, calculating the strip outlet temperature T.sub.i of each rolling stand, wherein the calculation flowchart is shown in FIG. 3, S8.1, calculating the outlet temperature T.sub.1 of the first rolling stand, wherein

(19) $T_1=T_1^r+\frac{1-\left({}_1/4\right)}{1-\left({}_1/2\right)}.Math.\frac{K_1\ln \left(\frac{1}{1-{}_1}\right)}{\mathrm{SJ}};$ S8.2, enabling i=1; S8.3, calculating the temperature T.sub.i,1 of the first section of strip behind the outlet of the i.sup.th rolling stand, i.e. T.sub.i,1=T.sub.i; S8.4, enabling j=2; S8.5, showing the relationship between the temperature of the j.sup.th section and the temperature of the j1.sup.th section by the following equation:

(20) ${}_{i,j}=-\frac{2k_0{w}_i^{0.264}\exp \left(9.45-0.1918C\right)1.163l}{v_{1i}{h}_{1i}\mathrm{Sm}}{T}_{i,j-1}^{-0.213}\left(T_{i,j-1}-T_c\right)+T_{i,j-1},$
wherein k.sub.0 is the influence coefficient of the nozzle shape and spraying angle, and 0.80<k.sub.0<1.2; S8.6, determining whether the in-equation j<m is established, if yes, enabling j=j+1, and then turning to step S8.5; otherwise, turning to step S8.7; S8.7, obtaining the temperature T.sub.i,m of the m.sup.th section by iterative calculation; S8.8, calculating the inlet temperature T.sub.i+1.sup.r of the i+1.sup.th rolling stand: T.sub.1+1.sup.r=T.sub.i,m; S8.9, calculating the outlet temperature T.sub.i+1 of the i+1.sup.th rolling stand, wherein

(21) $T_{i+1}=T_{i+1}^r+\frac{1-\left({}_{i+1}/4\right)}{1-\left({}_{i+1}/2\right)}.Math.\frac{K_{i+1}\ln \left(\frac{1}{1-{}_{i+1}}\right)}{\mathrm{SJ}};$ S8.10, determining whether the in-equation i<n is established, if yes, enabling i=i+1, and then turning to step S8.3; otherwise, turning to step S8.11; and S8.11, obtaining the outlet temperature T.sub.i of each rolling stand; S9, calculating an emulsion flow rate comprehensive optimization objective function F(X), wherein the calculation flowchart is shown in FIG. 4, S9.1, calculating the dynamic viscosity .sub.0i of an emulsion between roll gaps of each rolling stand, wherein .sub.0i=b.Math.exp(a.Math.T.sub.i), in the formula, a,b are the dynamic viscosity parameters of lubricating oil under atmospheric pressure; S9.2, calculating the oil film thickness .sub.i between the roll gaps of each rolling stand, wherein the calculation formula is as follows:

(22) 0${}_i=\frac{h_{0i}+h_{1i}}{2h_{0i}}.Math.k_c.Math.\frac{3{}_{0i}\left(v_{\mathrm{ri}}+v_{0i}\right)}{{}_i\left[1-e^{-\left(K-\frac{T_{0i}}{h_{0i}.Math.B}\right)}\right]}-k_{rg}.Math.\left(1+K_{rs}\right).Math.{\mathrm{Ra}}_{\mathrm{ir}0}.Math.e^{-B_L.Math.L_i}$
in the formula, k.sub.rg represents the coefficient of the strength of entrainment of lubricant by the longitudinal surface roughness of the work roll and the strip steel, and is in the range of 0.09-0.15, and K.sub.rs represents the impression rate, that is, the ratio of transferring the surface roughness of the working roll to the strip; and S9.3, calculating an emulsion flow rate comprehensive optimization objective function:

(23) $\{\begin{array}{c}F\left(X\right)=\frac{}{n}\underset{i=1}{\overset{n}{.Math.}}\sqrt{{\left({}_i-{}_{0i}\right)}^2}+\left(1-\right)\max .Math.{}_i-{}_{0i}.Math.\\ {{}_i}^{-}<{}_i<{}_i^{+}\end{array}$
in the formula, X={w.sub.i} is the optimization variable and is the distribution coefficient; S10, determining whether the in-equation F(X)<F.sub.0 is established, if yes, enabling w.sub.i.sup.y=w.sub.i, F.sub.0=F(X), and then turning to step S11; otherwise, turning directly to step S11; S11, determining whether the emulsion flow rate w.sub.i exceeds the a feasible region range, if yes, turning to step S12; otherwise, turning to step S7, wherein the feasible region of w.sub.i ranges from 0 to the maximum emulsion flow rate value allowed by the rolling mill. S12, outputting an optimal emulsion flow rate set value w.sub.i.sup.y, wherein w.sub.i.sup.y is the value of w.sub.i when the calculated value of F(X) in the feasible region is minimum.

Embodiment 1

(24) In order to further explain the application process of the related technology of the present application, the application process of an emulsion flow optimization method for a cold continuous rolling mill that aims to achieve vibration suppression is described by taking a 1730 cold continuous rolling mill in a cold rolling plant as an example.

(25) An emulsion flow optimization method for suppressing vibration of a cold continuous rolling mill includes the following steps:

(26) S1, collecting device feature parameters of the cold continuous rolling mill, wherein the 1730 cold continuous rolling mill in a cold rolling plant has 5 rolling stands in total, and the device feature parameters mainly include: the radius R.sub.i={210,212,230,230,228} mm of a working roll of each rolling stand, the surface linear velocity v.sub.ri={180,320,500,800,1150} m/min of a roll of each rolling stand, the original roughness Ra.sub.ir0={1.0,1.0,0.8,0.8,1.0} um of a working roll of each rolling stand, the roughness attenuation coefficient B.sub.L=0.01 of a working roll, the distance l=2700 mm between rolling stands, and the rolling kilometer L.sub.i={100,110,230,180,90} km after roll change of a working roll of each rolling stand, wherein i is 1, 2, . . . , n, and represents the ordinal number of rolling stands of the cold continuous rolling mill, and n=5 is the total number of rolling stands, the same below; S2, collecting key rolling process parameters of a strip, wherein the key rolling process parameters mainly include: the inlet thickness h.sub.0i={2.0,1.14,0.63,0.43,0.28} mm of each rolling stand, the outlet thickness h.sub.1i={1.14,0.63,0.43,0.28,0.18} mm of each rolling stand, strip width B=966 mm, the inlet speed v.sub.0i={110,190,342,552,848} m/min of each rolling stand, the outlet speed v.sub.1i={190,342,552,848,1214} m/min of each rolling stand, the inlet temperature T.sub.1.sup.r=110 C., strip deformation resistance K.sub.i={360,400,480,590,650} MPa of each rolling stand, rolling pressure P.sub.i={12800,11300,10500,9600,8800} kN of each rolling stand, back tension T.sub.0i={70,145,208,202,229} MPa of each rolling stand, front tension T.sub.1i={145,208,202,229,56} MPa of each rolling stand, emulsion concentration influence coefficient k.sub.c=0.9, pressure-viscosity coefficient =0.034 of a lubricant, strip density =7800 kg/m.sup.3, specific heat capacity S=0.47 kJ/(kg.Math. C.) of a strip, emulsion concentration C=4.2%, emulsion temperature T.sub.c=58 C. and thermal-work equivalent J=1; S3, defining process parameters involved in the process of emulsion flow optimization, wherein the process parameters mainly include that an over-lubrication film thickness critical value of each rolling stand is .sub.i.sup.+ and the friction coefficient at this time is u.sub.i.sup.+, an under-lubrication film thickness critical value is .sub.i.sup. and the friction coefficient at this time is u.sub.i.sup., the rolling reduction amount is h.sub.i=h.sub.0ih.sub.1i, the rolling reduction rate is

(27) ${}_i=\frac{h_i}{h_{0i}},$
the inlet temperature of each rolling stand is T.sub.i.sup.r, and the distance l=2700 mm between the rolling stands is evenly divided into m=30 sections, and the temperature in the sections is represented by T.sub.i,j (wherein, 1jm), and T.sub.i.sup.r=T.sub.i1,m, the over-lubrication judgment coefficient is A.sup.+, and the under-lubrication judgment coefficient is A.sup.; S4, setting the initial set value of an emulsion flow rate comprehensive optimization objective function of the cold continuous rolling mill that aims to achieve vibration suppression as F.sub.0=1.010.sup.10; S5, calculating the bite angle .sub.i of each rolling stand according to the rolling theory, wherein the calculation formula is

(28) ${}_i=\sqrt{\frac{h_i}{R_i^{}}},$
from which it can be obtained that .sub.i={0.0556,0.0427,0.0258,0.0223,0.0184}; S6, calculating the vibration determination index reference value .sub.0i of each rolling stand; S6.1, calculating the neutral angle .sub.i of each rolling stand, wherein the calculation formula is

(29) ${}_i=\frac{1}{2}\sqrt{\frac{h_i}{R_i^{}}}\left[1+\frac{1}{2u_i}\left(\sqrt{\frac{h_i}{R_i^{}}}+\frac{T_{i0}-T_{i1}}{P_i}\right)\right];$ S6.2, calculating to obtain u.sub.i.sup.+={0.0248,0.0186,0.0132,0.0136,0.0191} according to the formula

(30) $u_i^{+}=\frac{1}{2\left(2A^{+}-1\right)}\left(\sqrt{\frac{h_i}{R_i^{}}}+\frac{T_{i0}-T_{i1}}{P_i}\right)$
from steps S5 and S6.1 assuming that when

(31) $\frac{{}_i}{{}_i}=A^{+}=1,$
the roll gap is just in an over-lubrication state; S6.3, calculating an over-lubrication film thickness critical value .sub.i.sup.+ each rolling stand according to the relationship formula between the friction coefficient and the oil film thickness, i.e. u.sub.i=a.sub.i+b.sub.i.Math.e.sup.B.sup.i.sup..Math.*is i (in the formula, a.sub.i is the liquid friction influence coefficient, a.sub.i=0.0126, b.sub.i is the dry friction influence coefficient, b.sub.i=0.1416, and B.sub.i is the friction coefficient attenuation index, B.sub.i=2.4297), wherein the calculation formula is

(32) ${}_i^{+}=\frac{1}{B_i}\ln \frac{u_i^{+}-a_i}{b_i},$
from which it can be obtained that: .sub.i.sup.+={1.009,1.301,2.249,2.039,1.268} um; S6.4, calculating to obtain u.sub.i.sup.={0.1240,0.0930,0.0660,0.0680,0.0955} according to the formula

(33) $u_i^{-}=\frac{1}{2\left(2A^{-}-1\right)}\left(\sqrt{\frac{h_i}{R_i^{}}}+\frac{T_{i0}-T_{i1}}{P_i}\right)$
from steps S5 and S6.1 assuming that when

(34) $\frac{{}_i}{{}_i}=A^{-}=0.6,$
the roll gap is just in an under-lubrication state; S6.5, calculating an under-lubrication film thickness critical value .sub.i.sup. of each rolling stand according to the relationship formula between the friction coefficient and the oil film thickness, i.e. u.sub.i=a.sub.i+b.sub.i.Math.e.sup.B.sup.i.sup..Math.*is i, wherein the calculation formula is

(35) 0${}_i^{-}=\frac{1}{B_i}\ln \frac{u_i^{-}-a_i}{b_i},$
from which it can be obtained that: .sub.i.sup.={0.098,0.233,0.401,0.386,0.220} um; S6.6, calculating the vibration determination index reference value .sub.0i, wherein

(36) ${}_{0i}=\frac{{}_i^{+}+{}_i^{-}}{2},$
from which it can be obtained that: .sub.0i={0.554,0.767,1.325,1.213,0.744}; S7. Setting the emulsion flow rate of each rolling stand to be w.sub.i={900,900,900,900,900} L/min; S8, calculating the strip outlet temperature T.sub.i of each rolling stand, S8.1, calculating the outlet temperature T.sub.1 of the first rolling stand,

(37) $T_1=T_1^r+\frac{1-\left({}_1/4\right)}{1-\left({}_1/2\right)}.Math.\frac{K_1\ln \left(\frac{1}{1-{}_1}\right)}{\mathrm{SJ}}=110+\frac{1-\left(0\text{.43}/4\right)}{1-\left(0\text{.43}/2\right)}.Math.\frac{360\ln \left(\frac{1}{1-0.43}\right)}{7.8.Math.0.47.Math.1}=172.76C.$ S8.2, enabling i=1; S8.3, calculating the temperature T.sub.1,1 of the first section of strip behind the outlet of the first rolling stand, i.e. T.sub.i,1=T.sub.i=172.76 C.; S8.4, enabling j=2; S8.5, showing the relationship formula between the temperature of the j.sup.th section and the temperature of the j1.sup.th section by the following equation:

(38) $T_{i,j}=-\frac{2k_0{w}^{0.264}\exp \left(9.45-0.1918C\right)1.163l}{v_{1i},h_{1i}\mathrm{Sm}}{T}_{i,j-1}^{-0.213}\left(T_{i,j-1}-T_c\right)+T_{i,j-1},$
wherein k.sub.0=1.0; S8.6, determining whether the in-equation j<m is established: if yes, enabling j=j+1. and then turning to step S8.5; otherwise, turning to step S8.7; S8.7, obtaining the temperature T.sub.1,30=103.32 C. of the m=30.sup.th section by iterative calculation finally; S8.8, calculating the inlet temperature T.sub.2.sup.r of the second rolling stand: T.sub.2.sup.r=T.sub.1,m=103.32 C.; S8.9, calculating the outlet temperature T.sub.2 of the second rolling stand:

(39) $T_2=T_2^r+\frac{1-\left({}_2/4\right)}{1-\left({}_2/2\right)}.Math.\frac{K_2\ln \left(\frac{1}{1-{}_2}\right)}{SJ}=103.32+\frac{1-\left(0\text{.45}/4\right)}{1-\left(0\text{.45}/2\right)}.Math.\frac{400\ln \left(\frac{1}{1-0.45}\right)}{7800.Math.0.47.Math.1}=178.02C.;$ S8.10, determining whether the in-equation i<n is established: if yes, enabling i=i+1, and then turning to step S8.3; otherwise, turning to step S8.11; S8.11, obtaining the outlet temperature T.sub.i={172.76,178.02,186.59,194.35,206.33} C. of each rolling stand; S9, calculating an emulsion flow rate comprehensive optimization objective function F(X); S9.1, calculating the dynamic viscosity .sub.0i of an emulsion between roll gaps of each rolling stand, wherein .sub.0i=b.Math.exp(a.Math.T.sub.i), in the formula, a,b are the dynamic viscosity parameters of lubricating oil under atmospheric pressure, and it can be obtained from a=0.05, b=2.5 that .sub.0i={5.39,5.46,5.59,5.69,5.84}; S9.2, calculating the oil film thickness .sub.i between the roll gaps of each rolling stand according to the following formula:

(40) ${}_i=\frac{h_{0i}+h_{1i}}{2h_{0i}}.Math.k_c.Math.\frac{3{}_{0i}\left(v_{\mathrm{ri}}+v_{0i}\right)}{{}_i\left[1-e^{-\left(K-\frac{T_{0i}}{h_{0i}.Math.B}\right)}\right]}-k_{rg}.Math.\left(1+K_{rs}\right).Math.{\mathrm{Ra}}_{\mathrm{ir}0}.Math.e^{-B_{\mathrm{Li}}.Math.L_i}$
wherein in the formula, k.sub.rg represents the coefficient of the strength of entrainment of lubricant by the longitudinal surface roughness of the work roll and the strip steel, k.sub.rg=1.183, and K.sub.rs represents the impression rate, that is, the ratio of transferring the surface roughness of the working roll to the strip, K.sub.rs=0.576, from which it can be obtained that: .sub.i={0.784,0.963,2.101,2.043,1.326} um; S9.3, calculating an emulsion flow rate comprehensive optimization objective function:

(41) $\{\begin{array}{c}F\left(X\right)=\frac{}{n}\underset{i=1}{\overset{n}{.Math.}}\sqrt{{\left(-{}_{0i}\right)}^2}+\left(1-\right)\max .Math.{}_i-{}_{0i}.Math.\\ {}_i^{-}<{}_i<{}_i^{+}\end{array}$
in the formula, X={w.sub.i} is the optimization variable, =0.5 is the distribution coefficient, and thus F(X)=0.94; S10, enabling w.sub.i.sup.y=w.sub.i={900,900,900,900,900} L/min if F(X)=0.94<F.sub.0=110.sup.10 is established, F.sub.0=F(X)=0.94, turning to step S11, wherein in the subsequent x calculation processes, the corresponding F(X) is obtained with the change of w.sub.i, and the x.sup.th F.sub.0 is the x1.sup.th F(X). If the x.sup.th F(X) is smaller than the x1.sup.th F(X), it is judged that F(X)<F.sub.0 is established and turn to step S11; S11, determining whether the emulsion flow rate w.sub.i exceeds the feasible region range. If yes, turning to step S12; otherwise, turning to step S7; and S12, outputting an optimal emulsion flow rate set value w.sub.i.sup.y={1022,1050,1255,1698,1102} L/min.

Embodiment 2

(42) In order to further explain the application process of the related technology of the present application, the application process of an emulsion flow optimization method for a cold continuous rolling mill that aims to achieve vibration suppression is described by taking a 1420 cold continuous rolling mill in a cold rolling plant as an example.

(43) An emulsion flow optimization method for suppressing vibration of a cold continuous rolling mill includes the following steps: S1, collecting device feature parameters of the cold continuous rolling mill, wherein the 1420 cold continuous rolling mill in a cold rolling plant has 5 rolling stands in total, and the device feature parameters mainly include: the radius R.sub.i={211,213,233,233,229} mm of a working roll of each rolling stand, the surface linear velocity v.sub.ri={182,322,504,805,1153} m/min of a roll of each rolling stand, the original roughness Ra.sub.ir0={1.0,1.0,0.9,0.9,1.0} um of a working roll of each rolling stand, the roughness attenuation coefficient B.sub.L=0.015 of a working roll, the distance l=2750 mm between rolling stands, and the rolling kilometer L.sub.i={120,130,230,190,200} km after roll change of a working roll of each rolling stand, wherein i is 1, 2, . . . , n, and represents the ordinal number of rolling stands of the cold continuous rolling mill, and n=5 is the total number of rolling stands, the same below; S2, collecting key rolling process parameters of a strip, wherein the key rolling process parameters mainly include: the inlet thickness h.sub.0i={2.1,1.15,0.65,0.45,0.3} mm of each rolling stand, the outlet thickness h.sub.1i={1.15,0.65,0.45,0.3,0.15} mm of each rolling stand, strip width B=955 mm, the inlet speed v.sub.0i={115,193,346,555,852} m/min of each rolling stand, the outlet speed v.sub.1i={191,344,556,849,1217} m/min of each rolling stand, the inlet temperature T.sub.1.sup.r=115 C. strip deformation resistance K.sub.i={370,410,490,590,660} MPa of each rolling stand, rolling pressure P.sub.i={12820,11330,10510,9630,8820} kN of each rolling stand, back tension T.sub.0i={73,148,210,205,232}MPa of each rolling stand, front tension T.sub.1i={147,212,206,231,60} MPa of each rolling stand, emulsion concentration influence coefficient k.sub.c=0.9, pressure-viscosity coefficient =0.036 of a lubricant, strip density =7800 kg/m.sup.3, specific heat capacity S=0.49 kJ/(kg.Math. C.) of a strip, emulsion concentration C=4.5%, emulsion temperature T.sub.c=59 C. and thermal-work equivalent J=1; S3, defining process parameters involved in the process of emulsion flow optimization, wherein the process parameters mainly include that an over-lubrication film thickness critical value of each rolling stand is .sub.i.sup.+ and the friction coefficient at this time is u.sub.i.sup.+, an under-lubrication film thickness critical value is .sub.i.sup. and the friction coefficient at this time is u.sub.i.sup., the rolling reduction amount is h.sub.i=h.sub.0ih.sub.1i, the rolling reduction rate is

(44) ${}_i=\frac{h_i}{h_{0i}},$
the inlet temperature of each rolling stand is T.sub.i.sup.r, the distance l=2750 mm between the rolling stands is evenly divided into m=30 sections, and the temperature in the sections is represented by T.sub.i,j (wherein, 1jm), and T.sub.i.sup.r=T.sub.i1,m, the over-lubrication judgment coefficient is A.sup.+, and the under-lubrication judgment coefficient is A.sup.; S4, setting the initial set value of an emulsion flow rate comprehensive optimization objective function of a cold continuous rolling mill that aims to achieve vibration suppression as F.sub.0=1.010.sup.10; S5, calculating the bite angle .sub.i of each rolling stand according to the rolling theory, wherein the calculation formula is

(45) ${}_i=\sqrt{\frac{h_i}{R_i^{}}},$
from which it can be obtained that .sub.i={0.0566,0.0431,0.0261,0.0227,0.0188}; S6, calculating the vibration determination index reference value .sub.0i of each rolling stand; S6.1, calculating the neutral angle .sub.i of each rolling stand, wherein the calculation formula is

(46) ${}_i=\frac{1}{2}\sqrt{\frac{h_i}{R_i^{}}}\left[1+\frac{1}{2u_i}\left(\sqrt{\frac{h_i}{R_i^{}}}+\frac{T_{i0}-T_{i1}}{P_i}\right)\right];$ S6.2, calculating to obtain u.sub.i.sup.+={0.0251,0.0187,0.0135,0.0138,0.0193} according to the formula

(47) 0$u_i^{+}=\frac{1}{2\left(2A^{+}-1\right)}\left(\sqrt{\frac{h_i}{R_i^{}}}+\frac{T_{i0}-T_{i1}}{P_i}\right)$
from steps S5 and S6.1 assuming that when

(48) $\frac{{}_i}{{}_i}=A^{+}=1,$
the roll gap is just in an over-lubrication state; S6.3, calculating an over-lubrication film thickness critical value .sub.i.sup.+ of each rolling stand according to the relationship formula between the friction coefficient and the oil film thickness, i.e. u.sub.i=a.sub.i+b.sub.i.Math.e.sup.B.sup.i.sup..Math.*is i (in the formula, a.sub.i is the liquid friction influence coefficient, a.sub.i=0.0128, b.sub.i is the dry friction influence coefficient, b.sub.i=0.1426, and B.sub.i is the friction coefficient attenuation index, B.sub.i2.4307), wherein the calculation formula is

(49) ${{}_i}^{+}=\frac{1}{B_i}\ln \frac{{u_i}^{+}-a_i}{b_i},$
from which it can be obtained that: .sub.i.sup.+={1.011,1.321,2.253,2.041,1.272} um; S6.4, calculating to obtain u.sub.i.sup.={0.1243,0.0936,0.0664,0.0685,0.0955} according to the formula

(50) $u_i^{-}=\frac{1}{2\left(2A^{-}-1\right)}\left(\sqrt{\frac{h_i}{R_i^{}}}+\frac{T_{i0}-T_{i1}}{P_i}\right)$
from steps S5 and S6.1 assuming that when

(51) $\frac{{}_i}{{}_i}=A^{-}=0.6,$
the roll gap is just in an under-lubrication state; S6.5, calculating an under-lubrication film thickness critical value .sub.i.sup. of each rolling stand according to the relationship formula between the friction coefficient and the oil film thickness, i.e. u.sub.i=a.sub.i+b.sub.i.Math.e.sup.B.sup.i.sup..Math.*is i, wherein the calculation formula is

(52) ${}_i^{-}=\frac{1}{B_i}\ln \frac{u_i^{-}-a_i}{b_i},$
from which it can be obtained that: .sub.i.sup.={0.101,0.236,0.411,0.389,0.223} um; S6.6, calculating the vibration determination index reference value .sub.0i, wherein

(53) ${}_{0i}=\frac{{}_i^{+}+{}_i^{-}}{2},$
from which it can be obtained that: .sub.0i={0.557,0.769,1.327,1.215,0.746}; S7, setting the emulsion flow rate of each rolling stand to be w.sub.i={900,900,900,900,900} L/min; S8, calculating the strip outlet temperature T.sub.i of each rolling stand, S8.1, calculating the outlet temperature T.sub.1 of the first rolling stand,

(54) $T_1=T_1^r+\frac{1-\left({}_1/4\right)}{1-\left({}_1/2\right)}.Math.\frac{K_1\ln \left(\frac{1}{1-{}_1}\right)}{\mathrm{SJ}}=110+\frac{1-\left(0.43/4\right)}{1-\left(0.43/2\right)}.Math.\frac{360\ln \left(\frac{1}{1-0.43}\right)}{7.8.Math.0.47.Math.1}=175.81C.$ S8.2, enabling i=1; S8.3, calculating the temperature T.sub.1,1 of the first section of strip behind the outlet of the first rolling stand, i.e. T.sub.i,1=T.sub.i=175.81 C.; S8.4, enabling j=2; S8.5, showing the relationship between the temperature of the j.sup.th section and the temperature of the j1.sup.th section by the following equation:

(55) $T_{i,j}=-\frac{2k_0{w}_i^{0.264}\exp \left(9.45-0.1918C\right)1.163l}{v_{1i}{h}_{1i}Sm}{T}_{i,j-1}^{-0.213}\left(T_{i,j-1}-T_c\right)+T_{i,j-1},\mathrm{wherein}{k}_0=1.;$ S8.6, determining whether the in-equation j<m is established: if yes, enabling j=j+1. and then turning to step S8.5; otherwise, turning to step S8.7; S8.7, obtaining the temperature T.sub.1,30=105.41 C. of the m=30.sup.th section by iterative calculation finally; S8.8, calculating the inlet temperature T.sub.2.sup.r of the second rolling stand: T.sub.2.sup.r=T.sub.1,m=105.41 C.; S8.9, calculating the outlet temperature T.sub.2 of the second rolling stand

(56) $T_2=T_2^r+\frac{1-\left({}_2/4\right)}{1-\left({}_2/2\right)}.Math.\frac{K_2\ln \left(\frac{1}{1-{}_2}\right)}{\mathrm{SJ}}=103.32+\frac{1-\left(0.45/4\right)}{1-\left(0.45/2\right)}.Math.\frac{400\ln \left(\frac{1}{1-0.45}\right)}{7800.Math.0.47.Math.1}=182.52C.$ S8.10, determining whether the in-equation i<n is established: if yes, enabling i=i+1, and then turning to step S8.3; otherwise, turning to step S8.11; S8.11, obtaining the outlet temperature T.sub.i={175.86,179.36,189.77,196.65,207.54} C. of each rolling stand; S9, calculating an emulsion flow rate comprehensive optimization objective function F(X); S9.1, calculating the dynamic viscosity .sub.0i of an emulsion between roll gaps of each rolling stand, wherein .sub.0i=b.Math.exp (a.Math.T.sub.i), in the formula, a,b are the dynamic viscosity parameters of lubricating oil under atmospheric pressure, and it can be obtained from a=0.15, b=3.0 that .sub.0i={5.45,5.78,5.65,5.75,5.89}; S9.2, calculating the oil film thickness .sub.i between the roll gaps of each rolling stand according to the following formula:

(57) 0${}_i=\frac{h_{0i}+h_{1i}}{2h_{0i}}.Math.k_c.Math.\frac{3{}_{0i}\left(v_{\mathrm{ri}}+v_{0i}\right)}{{}_i\left[1-e^{-\left(K-\frac{T_{0i}}{h_{0i}.Math.B}\right)}\right]}-k_{rg}.Math.\left(1+K_{rs}\right).Math.{\mathrm{Ra}}_{\mathrm{ir}0}.Math.e^{-B_{Li}.Math.L_i}$
wherein in the formula, k.sub.rg represents the coefficient of the strength of entrainment of lubricant by the longitudinal surface roughness of the work roll and the strip steel, k.sub.rg=1.196, and K.sub.rs represents the impression rate, that is, the ratio of transferring the surface roughness of the working roll to the strip, K.sub.rs=0.584, from which it can be obtained that: .sub.i={0.795,0.967,2.132,2.056,1.337} um; S9.3, calculating an emulsion flow rate comprehensive optimization objective function:

(58) $\{\begin{array}{c}F\left(X\right)=\frac{}{n}\underset{i=1}{\overset{n}{.Math.}}\sqrt{{\left({}_i-{}_{0i}\right)}^2}+\left(1-\right)\max .Math.{}_i-{}_{0i}.Math.\\ {}_i^{-}<{}_i<{}_i^{+}\end{array}$
in the formula, X={w.sub.i} is the optimization variable, =0.5 is the distribution coefficient, and thus F(X)=0.98; S10, enabling w.sub.i.sup.y=w.sub.i={900,900,900,900,900} L/min if F(X)=0.98<F.sub.0=110.sup.10 is established, F.sub.0=F(X)=0.98, turning to step S11, wherein in the subsequent x calculation processes, the corresponding F(X) is obtained with the change of w.sub.i, and the x.sup.th F.sub.0 is the x1.sup.th F(X). If the x.sup.th F(X) is smaller than the x1.sup.th F(X), it is judged that F(X)<F.sub.0 is established and turn to step S11; S11, determining whether the emulsion flow rate w.sub.i exceeds the feasible region range. If yes, turning to step S12; otherwise, turning to step S7; and S12, outputting an optimal emulsion flow rate set value w.sub.i.sup.y={1029,1055,1261,1703,1109} L/min.

Embodiment 3

(59) In order to further explain the application process of the related technology of the present application, the application process of an emulsion flow optimization method for a cold continuous rolling mill that aims to achieve vibration suppression is described by taking a 1220 cold continuous rolling mill in a cold rolling plant as an example.

(60) An emulsion flow optimization method for suppressing vibration of a cold continuous rolling mill includes the following steps: S1, collecting device feature parameters of the cold continuous rolling mill, wherein the 1220 cold continuous rolling mill in a cold rolling plant has 5 rolling stands in total, and the device feature parameters mainly include: the radius R.sub.i={208,210,227,226,225} mm of a working roll of each rolling stand, the surface linear velocity v.sub.ri={176,317,495,789,1146} m/min of a roll of each rolling stand, the original roughness Ra.sub.ir0={0.9,0.9,0.7,0.7,0.8} um of a working roll of each rolling stand, the roughness attenuation coefficient B.sub.L=0.01 of a working roll, the distance l=2700 mm between rolling stands, and the rolling kilometer L.sub.i={152,102,215,165,70} km after roll change of a working roll of each rolling stand, wherein i is 1, 2, . . . , n, and represents the ordinal number of rolling stands of the cold continuous rolling mill, and n=5 is the total number of rolling stands, the same below; S2, collecting key rolling process parameters of a strip, wherein the key rolling process parameters mainly include: the inlet thickness h.sub.0i={1.8,1.05,0.57,0.39,0.25} mm of each rolling stand, the outlet thickness h.sub.1i={1.05,0.57,0.36,0.22,0.13} mm of each rolling stand, strip width B=876 mm, the inlet speed v.sub.0i={104,185,337,546,844} m/min of each rolling stand, the outlet speed v.sub.1i={188,337,548,845,1201}m/min of each rolling stand, the inlet temperature T.sub.1.sup.r=110 C. strip deformation resistance K.sub.i={355,395,476,580,640} MPa of each rolling stand, rolling pressure P.sub.i={12900,11200,10400,9600,8900} kN of each rolling stand, back tension T.sub.0i={74,141,203,201,219} MPa of each rolling stand, front tension T.sub.1i={140,203,199,224,50} MPa of each rolling stand, emulsion concentration influence coefficient k.sub.c=0.8, pressure-viscosity coefficient =0.035 of a lubricant, strip density =7800 kg/m.sup.3, specific heat capacity S=0.45 kJ/(kg.Math. C.) of a strip, emulsion concentration C=3.7%, emulsion temperature T.sub.c=55 C. and thermal-work equivalent J=1; S3, defining process parameters involved in the process of emulsion flow optimization, wherein the process parameters mainly include that an over-lubrication film thickness critical value of each rolling stand is .sub.i.sup.+ and the friction coefficient at this time is u.sub.i.sup.+, an under-lubrication film thickness critical value is .sub.i.sup. and the friction coefficient at this time is u.sub.i.sup., the rolling reduction amount is h.sub.i=h.sub.0ih.sub.1i, the rolling reduction rate is

(61) ${}_i=\frac{h_i}{h_{0i}},$
the inlet temperature of each rolling stand is T.sub.i.sup.r, the distance l=2700 mm between the rolling stands is evenly divided into m=30 sections, and the temperature in the sections is represented by T.sub.i,j (wherein, 1jm), and T.sub.i.sup.r=T.sub.i+1,m, the over-lubrication judgment coefficient is A.sup.+, and the under-lubrication judgment coefficient is A.sup.; S4, setting the initial set value of an emulsion flow rate comprehensive optimization objective function of a cold continuous rolling mill that aims to achieve vibration suppression as F.sub.0=1.010.sup.10; S5, calculating the bite angle .sub.i of each rolling stand according to the rolling theory, wherein the calculation formula is

(62) ${}_i=\sqrt{\frac{h_i}{R_i^{}}},$
from which it can be obtained that .sub.i={0.0546,0.0406,0.0247,0.0220,0.0179}; S6, calculating the vibration determination index reference value .sub.0i of each rolling stand; S6.1, calculating the neutral angle .sub.i of each rolling stand, wherein the calculation formula is

(63) ${}_i=\frac{1}{2}\sqrt{\frac{h_i}{R_i^{}}}\left[1+\frac{1}{2u_i}\left(\sqrt{\frac{h_i}{R_i^{}}}+\frac{T_{i0}-T_{i1}}{P_i}\right)\right];$ S6.2, calculating to obtain u.sub.i.sup.+={0.0242,0.0179,0.0127,0.0130,0.0185} according to the formula

(64) $u_i^{+}=\frac{1}{2\left(2A^{+}-1\right)}\left(\sqrt{\frac{h_i}{R_i^{}}}+\frac{T_{i0}-T_{i1}}{P_i}\right)$
from steps S5 and S6.1 assuming that when

(65) $\frac{{}_i}{{}_i}=A^{+}=1,$
the roll gap is just in an over-lubrication state; S6.3, calculating an over-lubrication film thickness critical value .sub.i.sup.+ of each rolling stand according to the relationship formula between the friction coefficient and the oil film thickness, i.e. u.sub.i=a.sub.i=b.sub.i.Math.e.sup.B.sup.i.sup..Math..sup.i (in the formula, a.sub.i is the liquid friction influence coefficient, a.sub.i=0.0125, b.sub.i is the dry friction influence coefficient, b.sub.i=0.1414, and B.sub.i is the friction coefficient attenuation index, B.sub.i=2.4280), wherein the calculation formula is

(66) ${}_i^{+}=\frac{1}{B_i}\ln \frac{u_i^{+}-a_i}{b_i},$
from which it can be obtained that: .sub.i.sup.+={1.001,1.289,2.232,2.037,1.268} um; S6.4, calculating to obtain u.sub.i.sup.={0.1241,0.0922,0.0610,0.0630,0.0935} according to the formula

(67) $u_i^{-}=\frac{1}{2\left(2A^{-}-1\right)}\left(\sqrt{\frac{h_i}{R_i^{}}}+\frac{T_{i0}-T_{i1}}{P_i}\right)$
from steps S5 and S6.1 assuming that when

(68) $\frac{{}_i}{{}_i}=A^{-}=0.6,$
the roll gap is just in an under-lubrication state; S6.5, calculating an under-lubrication film thickness critical value .sub.i.sup. of each rolling stand according to the relationship between the friction coefficient and the oil film thickness, i.e. u.sub.i=a.sub.i=b.sub.i.Math.e.sup.B.sup.i.sup..Math..sup.i, wherein the calculation formula is

(69) 0${}_i^{-}=\frac{1}{B_i}\ln \frac{u_i^{-}-a_i}{b_i},$
from which it can be obtained that: .sub.i.sup.={0.097,0.223,0.398,0.385,0.210} um; S6.6, calculating the vibration determination index reference value .sub.0i, wherein

(70) ${}_{0i}=\frac{{}_i^{+}+{}_i^{-}}{2},$
from which it can be obtained that: .sub.0i={0.548,0.762,1.321,1.207,0.736}; S7, setting the emulsion flow rate of each rolling stand to be w.sub.i={900,900,900,900,900} L/min; S8, calculating the strip outlet temperature T.sub.i of each rolling stand, S8.1, calculating the outlet temperature T.sub.1 of the first rolling stand,

(71) $T_1=T_1^r+\frac{1-\left({}_1/4\right)}{1-\left({}_1/2\right)}.Math.\frac{K_1\ln \left(\frac{1}{1-{}_1}\right)}{\mathrm{SJ}}=110+\frac{1-\left(0.43/4\right)}{1-\left(0.43/2\right)}.Math.\frac{360\ln \left(\frac{1}{1-0.43}\right)}{7.8.Math.0.47.Math.1}=169.96C.$ S8.2, enabling i=1; S8.3, calculating the temperature T.sub.1,1 of the first section of strip behind the outlet of the first rolling stand, i.e. T.sub.i,1=T.sub.i=169.96 C.; S8.4, enabling j=2; S8.5, showing the relationship between the temperature of the j.sup.th section and the temperature of the j1.sup.th section by the following equation:

(72) $T_{i,j}=-\frac{2k_0{w}_i^{0.264}\exp \left(9.45-0.1918C\right)1.163l}{v_{1i}{h}_{1i}\mathrm{Sm}}{T}_{i,j-1}^{-0.213}\left(T_{i,j-1}-T_c\right)+T_{i,j-1},$
wherein k.sub.0=1.0; S8.6, determining whether the in-equation j<m is established: if yes, enabling j=j+1. and then turning to step S8.5; otherwise, turning to step S8.7; S8.7, obtaining the temperature T.sub.1,30=101.25 C. of the m=30.sup.th section by iterative calculation finally; S8.8, calculating the inlet temperature T.sub.2.sup.r of the second rolling stand: T.sub.2.sup.r=T.sub.1,m=101.25 C.; S8.9, calculating the outlet temperature T.sub.2 of the second rolling stand:

(73) $T_2=T_2^r+\frac{1-\left({}_2/4\right)}{1-\left({}_2/2\right)}.Math.\frac{K_2\ln \left(\frac{1}{1-{}_2}\right)}{SJ}=103.32+\frac{1-\left(0.45/4\right)}{1-\left(0.45/2\right)}.Math.\frac{400\ln \left(\frac{1}{1-0.45}\right)}{7800.Math.0.47.Math.1}=175.86C.$ S8.10, determining whether the in-equation i<n is established: if yes, enabling i=i+1, and then turning to step S8.3; otherwise, turning to step S8.11; S8.11, obtaining the outlet temperature T.sub.i={177.96,172.78,184.59,191.77,203.33} C. of each rolling stand; S9, calculating an emulsion flow rate comprehensive optimization objective function F(X); S9.1, calculating the dynamic viscosity .sub.0i of an emulsion between roll gaps of each rolling stand, wherein .sub.0i=b.Math.exp(a.Math.T.sub.i), in the formula, a,b are the dynamic viscosity parameter of lubricating oil under atmospheric pressure, and it can be obtained from a=0.15, b=2.0 that .sub.0i={5.45,5.02,5.98,5.45,5.76}; S9.2, calculating the oil film thickness .sub.i between the roll gaps of each rolling stand according to the following formula:

(74) ${}_i=\frac{h_{0i}+h_{1i}}{2h_{0i}}.Math.k_c.Math.\frac{3{}_{0i}\left(v_{\mathrm{ri}}+v_{0i}\right)}{{}_i\left[1-e^{-\left(K-\frac{T_{0i}}{h_{0i}.Math.B}\right)}\right]}-k_{rg}.Math.\left(1+K_{rs}\right).Math.{\mathrm{Ra}}_{\mathrm{ir}0}.Math.e^{-B_{\mathrm{Li}}.Math.L_i}$
wherein in the formula, k.sub.rg represents the coefficient of the strength of entrainment of lubricant by the longitudinal surface roughness of the work roll and the strip steel, k.sub.rg=1.165, and K.sub.rs represents the impression rate, that is, the ratio of transferring the surface roughness of the working roll to the strip, K.sub.rs=0.566, from which it can be obtained that: .sub.i={0.774,0.926,2.088,2.032,1.318} um ;
S9.3, calculating an emulsion flow rate comprehensive optimization objective function:

(75) $\{\begin{array}{c}F\left(X\right)=\frac{}{n}\underset{i=1}{\overset{n}{.Math.}}\sqrt{{\left({}_i-{}_{0i}\right)}^2}+\left(1-\right)\max .Math.{}_i-{}_{0i}.Math.\\ {}_i^{-}<{}_i<{}_i^{+}\end{array}$
In the formula, X={w.sub.i} is the optimization variable, =0.5 is the distribution coefficient, and thus F(X)=0.91; S10, enabling w.sub.i.sup.y=w.sub.i={900,900,900,900,900} L/min if F(X)=0.91<F.sub.0=110.sup.10 is established, F.sub.0=F(X)=0.91, turning to step S11, wherein in the subsequent x calculation processes, the corresponding F(X) is obtained with the change of w.sub.i, and the x.sup.th F.sub.0 is the x1.sup.th F(X). If the x.sup.th F(X) is smaller than the x1.sup.th F(X), it is judged that F(X)<F.sub.0 is established and turn to step S11; S11, determining whether the emulsion flow rate w.sub.i exceeds the feasible region range. If yes, turning to step S12; otherwise, turning to step S7; and S12, outputting an optimal emulsion flow rate set value w.sub.i.sup.y={1016,1040,1266,1681,1111}L /min .

(76) The invention is applied to the five-machine-frame cold continuous rolling mills 1730, 1420 and 1220 in the cold rolling plant. According to the production experience of the cold rolling plant, the solution of the invention is feasible, and the effect is very obvious. The invention can be further applied to other cold continuous rolling mills, and the popularization prospect is relatively broad.

(77) To sum up, the technical solution of the invention is adopted, and the emulsion flow optimization method for suppressing vibration of the cold continuous rolling mill fully combines the device and process features of the cold continuous rolling mill, and aiming at the vibration defect problem, starting from the comprehensive optimization setting of the emulsion flow rate of each rolling stand, the method changes the previous idea of constant emulsion flow control for each rolling stand of the cold continuous rolling mill, and obtains the optimal set value of the emulsion flow rate for each rolling stand that aims to achieve vibration suppression by optimization; and the method greatly reduces the incidence of rolling mill vibration defects, improves production efficiency and product quality, and brings greater economic benefits for enterprises; and achieves the treatment for rolling mill vibration defects, and improves the surface quality and rolling process stability of a finished strip of a cold continuous rolling mill.