# Patent

## Title
Method of emulsion concentration optimization for cold continuous rolling mill set

## Abstract

Disclosed is a method of emulsion concentration optimization for a cold continuous rolling mill set for achieving vibration suppression, the method comprising: defining the process parameters involved in the process of emulsion concentration optimization; setting an initial set value of an emulsion concentration comprehensive optimization target function for a cold continuous rolling mill set for achieving vibration suppression; calculating a bite angle of each stand; calculating a vibration determination index reference value of each stand; setting the emulsion concentration of each stand; calculating the outlet temperature of a strip steel of each stand; calculating the dynamic viscosity of an emulsion in a roll gap of each stand; calculating the oil film thickness in the roll gap of each stand; calculating the emulsion concentration comprehensive optimization target function; determining whether the inequation F(X)<F0 is established; determining whether the concentration of the emulsion exceeds a feasible region range, and outputting the optimal emulsion concentration set value.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
This application is a 371 U.S. National Phase of PCT International Application No. PCT/CN2019/101118 filed on Aug. 16, 2019, which claims benefit and priority to Chinese patent application no. 201811144978.X filed on Sep. 29, 2018, contents of both are incorporated by reference herein in their entireties.
TECHNICAL FIELD
The present invention belongs to the field of cold continuous rolling, and in particular relates to a method of emulsion concentration optimization for a cold continuous rolling mill set for achieving vibration suppression.
BACKGROUND
As an important process parameter in the process lubrication system, emulsion concentration plays a vital role in the lubrication status of the roll gap of each stand of the cold continuous rolling mill set.
At the same time, the lubrication status of the roll gap directly affects the occurrence of vibration defects of the rolling mill.
If the roll gap is in an over-lubricated status, the friction coefficient is too small, which is likely to cause slip in the rolling process and leads to self-excited vibration of the rolling mill. If the roll gap is in an under-lubricated status, the average oil film thickness in the roll gap is less than the required minimum value, which is likely to cause the oil film in the roll gap to crack during the rolling process and leads to the sharp increase of the friction coefficient, thus changing the rolling pressure, causing the periodic fluctuation of the system stiffness and causing the self-excited vibration of the rolling mill as well.
In the past, the site generally relies on the speed of the rolling mill to suppress the occurrence of vibration defects, but this operation restricts the improvement of the production efficiency of the cold continuous rolling mill set and seriously affects the economic benefits of the enterprise.
The Chinese invention patent with the authorized announcement number of CN 103544340 B and the authorized announcement date of Mar. 2, 2016 discloses a “Method for Setting Emulsion Concentration in Extremely Thin Strip Rolling of Five-stand Cold Continuous Rolling Mill Set”. The method for setting emulsion concentration includes the following steps executed by a computer: 1) acquiring the main equipment of the mill set, the characteristics of the strip to be rolled, the main rolling process and the process lubrication system parameters; 2) defining relevant process parameters; 3) calculating the roll bending force and roll shifting amount; 4) assigning values to relevant search process parameters; 5) calculating concentration process parameters; 6) calculating the search process speed of the maximum rolling speed; 7) calculating the friction coefficient of each stand under the current conditions; 8) calculating the rolling force, rolling power, slip factor, thermal slip injury index and vibration coefficient of each stand under the current conditions; 9) calculating the thermal crown of the work roll of each stand; 10) calculating the outlet plate shape and the forced-contact width; and 11) obtaining and outputting the optimum proportion concentration. It can be seen that the patent aims at improving the rolling speed, ensuring the rolling efficiency, and avoiding slipping, thermal slip injury and vibration, so as to ensure the outlet plate shape of the final stand and the minimum forced-contact width at the work roll end.
Through research, it is found that under the premise of determining the process parameters such as rolling schedule, roll process, emulsion flow rate and initial temperature, the setting of emulsion concentration directly determines the lubrication status of the roll gap of each stand of the cold continuous rolling mill set, and can be used as the main process control means to suppress the vibration of the rolling mill.
However, at present, there has not been a report providing the technical solution of suppressing the vibration of the rolling mill by changing the previous mode of constant concentration control of an emulsion in each stand and taking the concentration of the emulsion in each stand as a variable to be optimized.
SUMMARY
The technical problem to be solved by the invention is to provide a method of emulsion concentration optimization for a cold continuous rolling mill set for achieving vibration suppression. The method changes the previous mode of constant concentration control of the emulsion in each stand, takes the concentration of the emulsion in each stand as a variable to be optimized, and carries out comprehensive optimization control on emulsion concentration. The lubrication status of roll gaps in each stand is optimized through the reasonable proportion of emulsion concentration in each stand, thus achieving the purposes of suppressing vibration of the rolling mill, improving product quality and production efficiency, and bringing economic benefits to enterprises.
The technical solution of the invention is to provide a method of emulsion concentration optimization for a cold continuous rolling mill set for achieving vibration suppression, wherein the method includes the following steps:
(a) acquiring apparatus characteristic parameters of a cold continuous rolling mill set;
(b) acquiring critical rolling process parameters of a strip material;
(c) acquiring process parameters involved in the process of emulsion concentration optimization;
(d) setting an initial set value F0=1.0×1010 of an emulsion concentration comprehensive optimization target function for a cold continuous rolling mill set for achieving vibration suppression;
the executing order of steps (a) to (d) is not limited;
(e) calculating a bite angle αi of each stand;
(f) calculating a vibration determination index reference value ξ0i of each stand;
(g) setting an emulsion concentration Ci of each stand;
(h) calculating the outlet temperature Ti of a strip steel of each stand;
(i) calculating the dynamic viscosity η0i of an emulsion in a roll gap of each stand;
(j) calculating an oil film thickness ξi in the roll gap of each stand;
(k) calculating the emulsion concentration comprehensive optimization target function F(X);
(l) determining whether the inequation F(X)<F0, is established, if yes, then setting Ciy=Ci, F0=F(X) and turning to step (m); if no, directly turning to step (m);
(m) determining whether the emulsion concentration Ci exceeds the range of a feasible region, if yes, turning to step (n); if no, turning to the step (g); wherein the feasible region refers to a region from 0 to the maximum emulsion concentration allowed by an apparatus, and wherein the allowed emulsion concentration of the apparatus is usually within 10%, and 0 ̃10% can be set as the feasible range;
(n) outputting the optimal emulsion concentration set value Ciy, wherein Ciy is the value of Ci when the calculated value of F(X) is minimum in the feasible region; and
(o) adjusting and controlling the emulsion concentration of each stand according to the optimal emulsion concentration set value Ciy in the step (n) by a control system of the cold continuous rolling mill set.
wherein calculating the vibration determination index reference value ξ0i of each stand in the step (f) comprising the following steps:
(f1) calculating an over-lubricated oil film thickness critical value ξi+ of each stand as follows:
it is assumed that when



 
  
   
    
     γ
     i
    
    
     α
     i
    
   
   =
   
    A
    +
   
  
  ,
 



the roll gap is just in an over-lubricated status, wherein γi is a neutral angle of each stand, and A+ is an over-lubricated determining coefficient;

calculating to obtain



 
  
   u
   i
   +
  
  =
  
   
    1
    
     2
     ⁢
     
      (
      
       
        2
        ⁢
        
         A
         +
        
       
       -
       1
      
      )
     
    
   
   ⁢
   
    (
    
     
      
       
        Δ
        ⁢
        
         h
         i
        
       
       
        R
        i
        ′
       
      
     
     +
     
      
       
        T
        
         i
         ⁢
         0
        
       
       -
       
        T
        
         i
         ⁢
         1
        
       
      
      
       P
       i
      
     
    
    )
   
  
 



according to




 
  
   
    α
    i
   
   =
   
    
     
      Δ
      ⁢
      
       h
       i
      
     
     
      R
      i
      ′
     
    
   
  
  ,
 



wherein Δhi is the rolling reduction, Δhi=h0i−h1i, h0i is the inlet thickness of each stand, h1i is the outlet thickness of each stand, and Ri′ is the flattening radius of a work roll of the ith stand,

and



 
  
   γ
   i
  
  =
  
   
    1
    2
   
   ⁢
   
    
     
      
       Δ
       ⁢
       
        h
        i
       
      
      
       R
       i
       ′
      
     
    
    [
    
     1
     +
     
      
       1
       
        2
        ⁢
        
         u
         i
        
       
      
      ⁢
      
       (
       
        
         
          
           Δ
           ⁢
           
            h
            i
           
          
          
           R
           i
           ′
          
         
        
        +
        
         
          
           T
           
            i
            ⁢
            0
           
          
          -
          
           T
           
            i
            ⁢
            1
           
          
         
         
          P
          i
         
        
       
       )
      
     
    
    ]
   
  
 



wherein T0i is the back tension of each stand, T1i is the front tension of each stand, and Pi is the rolling pressure of each stand,
according to the relationship between the friction coefficient and the oil film thickness: ui=ai+bi·eHi·ξi, wherein ai is the liquid friction influence coefficient, bi is the dry friction influence coefficient, and Bi is the friction coefficient attenuation index, the over-lubricated oil film thickness critical value ξi+ of each stand is calculated by:



 
  
   
    ξ
    i
    +
   
   =
   
    
     1
     
      B
      i
     
    
    ⁢
    ln
    ⁢
    
     
      
       u
       i
       +
      
      -
      
       a
       i
      
     
     
      b
      i
     
    
   
  
  ;
 



(f2) calculating an under-lubricated oil film thickness critical value ξi− of each stand as follows:
it is assumed that when



 
  
   
    
     γ
     i
    
    
     α
     i
    
   
   =
   
    A
    -
   
  
  ,
 



the roll gap is just in an under-lubricated status, wherein A− is the under-lubricated determining coefficient, the following equation can be obtained:




 
  
   
    u
    i
    -
   
   =
   
    
     1
     
      2
      ⁢
      
       (
       
        
         2
         ⁢
         
          A
          -
         
        
        -
        1
       
       )
      
     
    
    ⁢
    
     (
     
      
       
        
         Δ
         ⁢
         
          h
          i
         
        
        
         R
         i
         ′
        
       
      
      +
      
       
        
         T
         
          i
          ⁢
          0
         
        
        -
        
         T
         
          i
          ⁢
          1
         
        
       
       
        P
        i
       
      
     
     )
    
   
  
  ;
 



and the under-lubricated oil film thickness critical value ξi− of each stand is calculated by:



 
  
   
    ξ
    i
    -
   
   =
   
    
     1
     
      B
      i
     
    
    ⁢
    ln
    ⁢
    
     
      
       u
       i
       -
      
      -
      
       a
       i
      
     
     
      b
      i
     
    
   
  
  ;
 



and

(f3) calculating the vibration determination index reference value ξ0i of each stand as follows:



 
  
   ξ
   
    0
    ⁢
    i
   
  
  =
  
   
    
     
      ξ
      i
      +
     
     +
     
      ξ
      i
      -
     
    
    2
   
   .
  
 



wherein calculating the outlet temperature Ti of the strip steel of each stand in the step (h) comprising the following steps:
(h1) calculating the outlet temperature T1 of the first stand:



 
  
   
    T
    1
   
   =
   
    
     T
     1
     Inlet
    
    +
    
     
      
       1
       -
       
        (
        
         
          ε
          1
         
         /
         4
        
        )
       
      
      
       1
       -
       
        (
        
         
          ε
          1
         
         /
         2
        
        )
       
      
     
     ·
     
      
       
        K
        1
       
       ⁢
       
        ln
        ⁡
        (
        
         1
         
          1
          -
          
           ε
           1
          
         
        
        )
       
      
      
       p
       ⁢
       S
       ⁢
       J
      
     
    
   
  
  ,
 



wherein T1Inlet is the inlet temperature of each stand,



 
  
   
    ε
    i
   
   =
   
    
     Δ
     ⁢
     
      h
      i
     
    
    
     h
     
      0
      ⁢
      i
     
    
   
  
  ,
 



Δhi=h0i−h1i, h0i is the inlet thickness of each stand, h1i is the outlet thickness of each stand, ρ is the density of strip steel, S is the specific heat capacity of strip steel, J is the mechanical equivalent of heat, and Ki is the resistance to deformation of the strip steel of each stand;

(h2) setting i to 1;
(h3) setting the temperature of the first section of strip steel behind the outlet of the ith stand to Ti, 1, i.e., Ti, 1=Ti;
(h4) setting j to 2;
(h5) the relationship between the temperature of the jth section and the temperature of the j−1th section satisfies the following equation:



 
  
   
    T
    
     i
     ,
     j
    
   
   =
   
    
     
      -
      
       
        2
        ⁢
        
         k
         0
        
        ⁢
        
         w
         0.264
        
        ⁢
        
         exp
         ⁡
         (
         
          9.45
          -
          
           0.1918
           
            C
            i
           
          
         
         )
        
        ×
        1.
        1
        ⁢
        6
        ⁢
        3
        ⁢
        l
       
       
        
         v
         
          1
          ⁢
          i
         
        
        ⁢
        
         h
         
          1
          ⁢
          i
         
        
        ⁢
        ρ
        ⁢
        S
        ⁢
        m
       
      
     
     ⁢
     
      
       T
       
        i
        ,
        
         j
         -
         1
        
       
       
        -
        0.213
       
      
      (
      
       
        T
        
         i
         ,
         
          j
          -
          1
         
        
       
       -
       
        T
        c
       
      
      )
     
    
    +
    
     T
     
      i
      ,
      
       j
       -
       1
      
     
    
   
  
  ,
 



wherein k0 is the influence coefficient of nozzle shape and spraying angle, 0.8<k=0<1.2, w is the flow rate of the emulsion, l is the distance between stands, and the distance l between stands is equally divided into m sections, the temperature in the section is represented by Ti, j, v1i is the outlet speed of each stand, h1i is the outlet thickness of each stand, ρ is the density of the strip steel, S is the specific heat capacity of the strip steel, Ti is the outlet temperature of each stand, and Tc is the temperature of the emulsion;
(h6) determining whether the inequation j<m is established, if yes, then setting j=j+1 and turning to the step (h5); if no, turning to step (h7);
(h7) obtaining the temperature Ti,m of the mth section via iterative calculation;
(h8) calculating the inlet temperature Ti+1Inlet of the i+1th stand:

Ti+1Inlet=Ti,m;

(h9) calculating the outlet temperature Ti+1 of the i+1′ stand;



 
  
   
    T
    
     i
     +
     1
    
   
   =
   
    
     T
     
      i
      +
      1
     
     Inlet
    
    +
    
     
      
       1
       -
       
        (
        
         
          ε
          
           i
           +
           1
          
         
         /
         4
        
        )
       
      
      
       1
       -
       
        (
        
         
          ε
          
           i
           +
           1
          
         
         /
         2
        
        )
       
      
     
     ·
     
      
       
        K
        
         i
         +
         1
        
       
       ⁢
       
        ln
        ⁡
        (
        
         1
         
          1
          -
          
           ε
           
            i
            +
            1
           
          
         
        
        )
       
      
      
       ρ
       ⁢
       SJ
      
     
    
   
  
  ;
 



(h10) determining whether the inequation i<n is established, if yes, setting i=i+1 and turning to the step (h3); if no, turning to step (h11); and
(h11) obtaining the outlet temperature Ti of each stand.
Specifically, apparatus characteristic parameters of the cold continuous rolling mill set at least include:
the radius Ri of a work roll of each stand; the surface linear speed vri of a roll of each stand; the original roughness Rair0 of a work roll of each stand; the roughness attenuation coefficient BL, of a work roll; the distance l between stands; and rolling kilometers Li of a work roll of each stand after roll change; wherein i=1, 2, . . . , n; i represents the ordinal number of the stand of the cold continuous rolling mill set, and n is the total number of the stands.
Specifically, the critical rolling process parameters of the strip material at least include:
the inlet thickness h0i of each stand; the outlet thickness h1i of each stand; the width B of the strip steel; the inlet speed v0i of each stand; the outlet speed v1i of each stand; the inlet temperature T1Inlet; the resistance to deformation Ki of the strip steel of each stand; the rolling pressure Pi of each stand; the back tension T0i of each stand; the front tension T11 of each stand; the concentration influence coefficient kc of the emulsion; the viscosity compression coefficient θ of a lubricant; the density ρ of the strip steel; the specific heat capacity S of the strip steel; the flow rate w of the emulsion; the temperature Tc of the emulsion; and the mechanical equivalent of heat J.
Specifically, the process parameters involved in the process of emulsion concentration optimization at least include: the over-lubricated oil film thickness critical value ξi+ and the corresponding friction coefficient ui+ of each stand; the under-lubricated oil film thickness critical value ξi− and the corresponding friction coefficient ui− of each stand; the vibration determination index reference value ξ0i; the rolling reduction Δhi=h0i−h1i; the rolling reduction ratio



 
  
   
    ε
    i
   
   =
   
    
     Δ
     ⁢
     
      h
      i
     
    
    
     h
     
      0
      ⁢
      i
     
    
   
  
  ;
 



the inlet temperature TiInlet of each stand; and the outlet temperature Ti of each stand; the distance l between stands, which is equally divided into m sections, the temperature in the section represented by Ti, j, wherein 1≤j≤m and TiInlet=Ti−1,m; the over-lubricated determining coefficient A+; and the under-lubricated determining coefficient A−.

Further, the calculation formula for calculating the bite angle αi of each stand is as follows:



 
  
   α
   i
  
  =
  
   
    
     Δ
     ⁢
     
      h
      i
     
    
    
     R
     i
     ′
    
   
  
 



wherein Ri′ is the flattening radius of the work roll of the ith stand, and is a process value in rolling pressure calculation.

Further, the calculation formula for calculating the dynamic viscosity η0i of the emulsion in the roll gap of each stand is as follows:

η0i=b·exp(−a·Ti)

wherein a and b are dynamic viscosity parameters of the lubricant under atmospheric pressure.
Further, the calculation formula for calculating the oil film thickness ξi in the roll gap of each stand is as follows:



 
  
   
    ξ
    i
   
   =
   
    
     
      
       
        h
        
         0
         ⁢
         i
        
       
       +
       
        h
        
         1
         ⁢
         i
        
       
      
      
       2
       ⁢
       
        h
        
         0
         ⁢
         i
        
       
      
     
     ·
     
      k
      c
     
     ·
     
      
       3
       ⁢
       θ
       ⁢
       
        
         η
         
          0
          ⁢
          i
         
        
        (
        
         
          v
          ri
         
         +
         
          v
          
           0
           ⁢
           i
          
         
        
        )
       
      
      
       
        α
        i
       
       [
       
        1
        -
        
         e
         
          -
          
           θ
           ⁡
           (
           
            
             K
             i
            
            -
            
             
              T
              
               0
               ⁢
               i
              
             
             
              
               h
               
                0
                ⁢
                i
               
              
              ·
              B
             
            
           
           )
          
         
        
       
       ]
      
     
    
    -
    
     
      k
      rg
     
     ·
     
      (
      
       1
       +
       
        K
        rs
       
      
      )
     
     ·
     
      Ra
      
       ir
       ⁢
       0
      
     
     ·
     
      e
      
       
        -
        
         B
         Li
        
       
       ·
       
        L
        i
       
      
     
    
   
  
  ,
 



wherein h0i is the inlet thickness of each stand, h1i is the outlet thickness of each stand, kc is the emulsion concentration influence coefficient, θ is the viscosity compression coefficient of the lubricant, Ki is the resistance to deformation of the strip steel of each stand, η0i is the dynamic viscosity of the emulsion in the roll gap of each stand, v0i is the inlet speed of each stand, vri is the surface linear speed of a roll of each stand, T0i is the back tension of each stand, B is the width of the strip steel, krg represents the coefficient of the strength of the surface longitudinal roughness of the work roll and the strip steel to entrain the lubricant, the value of which is from 0.09 to 0.15, Krs represents the impression rate, i.e., the ratio of transmitting surface roughness of the work roll to the strip steel, the value of which is from 0.2 to 0.6, Rair0 is the original roughness of a work roll of each stand, BL is the roughness attenuation coefficient of the work roll, and Li is the rolling kilometers of a work roll of each stand after roll change.
Further, the emulsion concentration comprehensive optimization target function is calculated according to the following formula:



 
  {
  
   
    
     
      
       
        F
        ⁡
        (
        X
        )
       
       =
       
        
         
          λ
          n
         
         ⁢
         
          
           ∑
           
            i
            =
            1
           
           n
          
          
           
            
             (
             
              
               ξ
               i
              
              -
              
               ξ
               
                0
                ⁢
                i
               
              
             
             )
            
            2
           
          
         
        
        +
        
         
          (
          
           1
           -
           λ
          
          )
         
         ⁢
         max
         ⁢
         
          
           ❘
           "\[LeftBracketingBar]"
          
          
           
            ζ
            i
           
           -
           
            ξ
            
             0
             ⁢
             i
            
           
          
          
           ❘
           "\[RightBracketingBar]"
          
         
        
       
      
     
    
    
     
      
       
        ξ
        i
        -
       
       <
       
        ξ
        i
       
       <
       
        ξ
        i
        +
       
      
     
    
   
   ,
  
 



wherein x={Ci} is an optimized variable, and λ is a distribution coefficient.
In the disclosure, as long as the execution of the next step is not conditional on the result of the previous step, it is not necessary to proceed according to the above order, unless the execution of the next step depends on the previous step.
Compared with the prior art, the invention has the following advantages:
1. The lubricated status of roll gaps in each stand is optimized through reasonable proportion of the emulsion concentration of each stand, thus achieving the purposes of suppressing vibration of the rolling mill, and improving product quality and production efficiency.
2. As a result of a large number of field tests and theoretical research, based on the apparatus characteristics and rolling process features of the cold continuous rolling mill set, a method of emulsion concentration optimization for a cold continuous rolling mill set for achieving vibration suppression is put forward. The method realizes the optimal ratio of emulsion concentration for each stand of the cold continuous rolling mill set, achieves suppressing vibration of the rolling mill and improving product quality and production efficiency, and brings great economic benefits to enterprises.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a schematic flow diagram of the general technical solutions according to the present invention;
FIG. 2 is a schematic diagram of the calculation process of the reference value of the vibration determination index according to the present invention; and
FIG. 3 is a schematic diagram of the calculation process of the strip steel outlet temperature of each stand according to the present invention.



DETAILED DESCRIPTION
The invention will be further described below in conjunction with the accompany drawings and embodiments.
In order to further explain the application process of the related technologies mentioned in the invention, a 1730 cold continuous rolling mill set of a cold rolling plant is taken as an example. The application process of the method of emulsion concentration optimization for a cold continuous rolling mill set for achieving vibration suppression is described in detail.
At first, according to the various steps shown in FIG. 1, the relevant parameters are determined in turn. Then the parameters are substituted into corresponding formulas for calculation, and the desired optimal emulsion concentration set value Ciy is determined or obtained. Finally, the emulsion concentration of each stand is controlled according to the determined optimal emulsion concentration set value, and the comprehensive optimization control is carried out to achieve suppressing the vibration of the rolling mill.
Specifically, in the step (a), the acquired apparatus characteristic parameters of a cold continuous rolling mill set mainly include:
the radius of a work roll of each stand: Ri={210, 212, 230, 230, 228} mm;
the surface linear speed of a roll of each stand: vri={180, 320, 500, 800, 1150} m/min;
the original roughness of a work roll of each stand: Rair0={1.0, 1.0, 0.8, 0.8, 1.0} um;
the roughness attenuation coefficient of the work roll: BL=0.01;
the distance between stands: l=2700 mm; and
the rolling kilometers of a work roll of each stand after roll change: Li={100, 110, 230, 180, 90} km; wherein i=1, 2, . . . , n; i represents the ordinal number of the stand of the cold continuous rolling mill set; and n=5, which is the total number of the stands, the same below;
subsequently, in the step (b), the acquired critical rolling process parameters of the strip material mainly include:
the inlet thickness of each stand: h0i={2.0, 1.14, 0.63, 0.43, 0.28} mm;
the outlet thickness of each stand: h1i={1.14, 0.63, 0.43, 0.28, 0.18} mm;
the width of the strip steel: B=966 mm;
the inlet speed of each stand: v0i={110, 190, 342, 552, 848} m/min;
the outlet speed of each stand: v1i={190, 342, 552, 848, 1214} m/min;
the inlet temperature: T1Inlet=10° C.;
the resistance to deformation of the strip steel of each stand:
Ki={36, 40, 480, 590, 650} MPa;
the rolling pressure of each stand: Pi={12800, 11300, 10500, 9600, 8800} kN;
the back tension of each stand: T0i={70, 145, 208, 202, 229} MPa;
the front tension of each stand: T1i={145, 208, 202, 229, 56} MPa;
the concentration influence coefficient of the emulsion: kc=0.9;
the viscosity compression coefficient of the lubricant: θ=0.034;
the density of the strip steel: ρ=7800 kg/m3;
the specific heat capacity of the strip steel: S=0.47 kJ/(kg·° C.);
the flow rate of the emulsion: w=900 m/min;
the temperature of the emulsion: Tc=58° C.;
the mechanical equivalent of heat: J=1; and
the flattening radius of the work roll of the it stand: Ri′={278.2, 279.7, 300.5, 301.6, 295.4};
subsequently, in the step (c), the acquired process parameters involved in the process of emulsion concentration optimization mainly include: the over-lubricated oil film thickness critical value ξi+ and the corresponding friction coefficient ui+ of each stand; the under-lubricated oil film thickness critical value ξi− and the corresponding friction coefficient ui− of each stand; the vibration determination index reference value ξ0i, the rolling reduction: Δhi=h0i−h1i={0.86, 0.51, 0.2, 0.15, 0.1}; the rolling reduction ratio:



 
  
   
    ε
    i
   
   =
   
    
     
      Δ
      ⁢
      
       h
       i
      
     
     
      h
      
       0
       ⁢
       i
      
     
    
    =
    
     {
     
      
       
        0
        .
        4
       
       ⁢
       3
      
      ,
      
       
        0
        .
        4
       
       ⁢
       5
      
      ,
      
       
        0
        .
        3
       
       ⁢
       2
      
      ,
      
       
        0
        .
        3
       
       ⁢
       5
      
      ,
      
       
        0
        .
        3
       
       ⁢
       6
      
     
     }
    
   
  
  ;
 



the inlet temperature TiInlet each stand; and the outlet temperature Ti of each stand; the distance l between stands is 2700 mm, which is equally divided into m=30 sections, the temperature in the section is represented by Ti, j (wherein 1≤j≤m), and TiInlet=Ti−1,m; the over-lubricated determining coefficient A+; and the under-lubricated determining coefficient A−.

subsequently, in the step (d), an initial set value F0=1.0×1010 of an emulsion concentration comprehensive optimization target function for a cold continuous rolling mill set for achieving vibration suppression is set;
subsequently, in the step (e), the bite angle αi of each stand is calculated according to the rolling theory using the calculation formula:



 
  
   
    α
    i
   
   =
   
    
     
      Δ
      ⁢
      
       h
       i
      
     
     
      R
      i
      ′
     
    
   
  
  ,
 



obtaining αi={0.0556, 0.0427, 0.0258, 0.0223, 0.0184};

subsequently, in the step (f), the vibration determination index reference value ξ0i of each stand is calculated according to the sub-steps shown in the FIG. 2:
(f1) calculating the over-lubricated oil film thickness critical value ξi+ of each stand:
it is assumed that when



 
  
   
    
     γ
     i
    
    
     α
     i
    
   
   =
   
    
     A
     +
    
    =
    1
   
  
  ,
 



the roll gap is just in an over-lubricated status; according to




 
  
   α
   i
  
  =
  
   
    
     Δ
     ⁢
     
      h
      i
     
    
    
     R
     i
     ′
    
   
  
 



and the calculation formula of the neutral angle γi of each stand,




 
  
   u
   i
   +
  
  =
  
   
    1
    
     2
     ⁢
     
      (
      
       
        2
        ⁢
        
         A
         +
        
       
       -
       1
      
      )
     
    
   
   ⁢
   
    (
    
     
      
       
        Δ
        ⁢
        
         h
         i
        
       
       
        R
        i
        ′
       
      
     
     +
     
      
       
        T
        
         i
         ⁢
         0
        
       
       -
       
        T
        
         i
         ⁢
         1
        
       
      
      
       P
       i
      
     
    
    )
   
  
 



is obtained, obtaining ui+={00248, 0.0186, 0.0132, 0.0136, 0.0191};

according to the relationship between the friction coefficient and the oil film thickness: ui=ai+bi·eBi·ξi (in the embodiment, ai=0.0126, bi=0.1416, and Bi=−2.4297), the over-lubricated oil film thickness critical value ξi+ of each stand is calculated by the calculation formula:



 
  
   
    ξ
    i
    +
   
   =
   
    
     1
     
      B
      i
     
    
    ⁢
    ln
    ⁢
    
     
      
       u
       i
       +
      
      -
      
       a
       i
      
     
     
      b
      i
     
    
   
  
  ,
 



obtaining ξi+={1.009, 1.301, 2249, 2.039, 1.268} um;

(f2) calculating the under-lubricated oil film thickness critical value ξi+ of each stand is calculated:
it is assumed that when



 
  
   
    
     γ
     i
    
    
     α
     i
    
   
   =
   
    
     A
     -
    
    =
    0.6
   
  
  ,
 



the roll gap is just in an under-lubricated status,




 
  
   
    u
    i
    -
   
   =
   
    
     1
     
      2
      ⁢
      
       (
       
        
         2
         ⁢
         
          A
          -
         
        
        -
        1
       
       )
      
     
    
    ⁢
    
     (
     
      
       
        
         Δ
         ⁢
         
          h
          i
         
        
        
         R
         i
         ′
        
       
      
      +
      
       
        
         T
         
          i
          ⁢
          0
         
        
        -
        
         T
         
          i
          ⁢
          1
         
        
       
       
        P
        i
       
      
     
     )
    
   
  
  ,
 



ui−={0.1240, 0.0930, 0.0660, 0.0680, 0.0955};

according to the relationship between the friction coefficient and the oil film thickness: ui=ai+bi·eBi·ξi, the under-lubricated oil film thickness critical value ξi− of each stand is calculated by the calculation formula



 
  
   
    ξ
    i
    -
   
   =
   
    
     1
     
      B
      i
     
    
    ⁢
    ln
    ⁢
    
     
      
       u
       i
       -
      
      -
      
       a
       i
      
     
     
      b
      i
     
    
   
  
  ,
 



obtaining ξi−={0.098, 0.233, 0.401, 0.386, 0.220} um;

subsequently, in the step (f3), the vibration determination index reference value ξ0i is calculated, wherein



 
  
   
    ξ
    
     0
     ⁢
     i
    
   
   =
   
    
     
      ξ
      i
      +
     
     +
     
      ξ
      i
      -
     
    
    2
   
  
  ,
 



obtaining ξ0i={0.554, 0.767, 1.325, 1.213, 0.744};

subsequently, in the step (g), the emulsion concentration of each stand is set as Ci={4.2, 4.2, 4.2, 4.2, 4.2} %.
Then, in the step (h), the outlet temperature Ti of the strip steel of each stand is calculated according to the sub-steps shown in the FIG. 3.
Subsequently, in the step (h1), the outlet temperature T1 of the first stand is calculated as follows:



 
  
   T
   1
  
  =
  
   
    
     T
     1
     Inlet
    
    +
    
     
      
       1
       -
       
        (
        
         
          ε
          1
         
         /
         4
        
        )
       
      
      
       1
       -
       
        (
        
         
          ε
          1
         
         /
         2
        
        )
       
      
     
     ·
     
      
       
        K
        1
       
       ⁢
       
        ln
        ⁡
        (
        
         1
         
          1
          -
          
           ε
           1
          
         
        
        )
       
      
      
       ρ
       ⁢
       SJ
      
     
    
   
   =
   
    
     
      1
      ⁢
      1
      ⁢
      0
     
     +
     
      
       
        1
        -
        
         (
         
          0.
          4
          ⁢
          
           3
           /
           4
          
         
         )
        
       
       
        1
        -
        
         (
         
          0.43
          /
          2
         
         )
        
       
      
      ·
      
       
        360
        ⁢
        
         ln
         ⁡
         (
         
          1
          
           1
           -
           0.13
          
         
         )
        
       
       
        7.8
        ·
        0.47
        ·
        1
       
      
     
    
    =
    
     
      172.76
      ∘
     
       
     
      C
      .
     
    
   
  
 



Subsequently, in the step (h2), i is set to 1.
Subsequently, in the step (h3), the temperature of the first section of strip steel behind the outlet of the f0 stand is set to Ti, 1, i.e., Ti, 1=Ti=172.76° C.
Subsequently, in the step (h4), j is set to 2.
Subsequently, in the step (h5), the relationship between the temperature of the jth section and the temperature of the j−1th section satisfies the following equation:



 
  
   
    T
    
     i
     ,
     j
    
   
   =
   
    
     
      -
      
       
        2
        ⁢
        
         k
         0
        
        ⁢
        
         w
         0.264
        
        ⁢
        
         exp
         ⁡
         (
         
          9.45
          -
          
           
            0
            .
            1
           
           ⁢
           9
           ⁢
           1
           ⁢
           8
           ⁢
           
            C
            i
           
          
         
         )
        
        ×
        1.163
        l
       
       
        
         v
         
          1
          ⁢
          i
         
        
        ⁢
        
         h
         
          1
          ⁢
          i
         
        
        ⁢
        ρ
        ⁢
        S
        ⁢
        m
       
      
     
     ⁢
     
      
       T
       
        i
        ,
        
         j
         -
         1
        
       
       
        
         -
         
          0
          .
          2
         
        
        ⁢
        1
        ⁢
        3
       
      
      (
      
       
        T
        
         i
         ,
         
          j
          -
          1
         
        
       
       -
       
        T
        c
       
      
      )
     
    
    +
    
     T
     
      i
      ,
      
       j
       -
       1
      
     
    
   
  
  ,
 



wherein k0 is the influence coefficient of the nozzle shape and the spraying angle, and k0=1.
Subsequently, in the step (h6), whether the inequation j<m is established is determined, if yes, then setting j=j+1 and turning to step (h5); if no, turning to step (h7).
Subsequently, in the step (h7), the temperature of the 30th section (when m=30) is obtained via iterative calculation: T1.30=103.32° C.
Subsequently, in the step (h8), the inlet temperature T2Inlet of the second stand is calculated as follows: T2Inlet=T1,m=103.32° C.
Subsequently, in the step (h9), the outlet temperature T2 of the second stand is calculated as follows:



 
  
   T
   2
  
  =
  
   
    
     T
     2
     Inlet
    
    +
    
     
      
       1
       -
       
        (
        
         
          ε
          2
         
         /
         4
        
        )
       
      
      
       1
       -
       
        (
        
         
          ε
          2
         
         /
         2
        
        )
       
      
     
     ·
     
      
       
        K
        2
       
       ⁢
       
        ln
        ⁡
        (
        
         1
         
          1
          -
          
           ε
           2
          
         
        
        )
       
      
      
       p
       ⁢
       S
       ⁢
       J
      
     
    
   
   =
   
    
     103.32
     +
     
      
       
        1
        -
        
         (
         
          0.45
          /
          4
         
         )
        
       
       
        1
        -
        
         (
         
          0.45
          /
          2
         
         )
        
       
      
      ·
      
       
        400
        ⁢
        
         ln
         ⁡
         (
         
          1
          
           1
           -
           0.45
          
         
         )
        
       
       
        7800
        ·
        0.47
        ·
        1
       
      
     
    
    =
    
     1
     ⁢
     7
     ⁢
     
      8
      .
      0
     
     ⁢
     
      2
      ∘
     
     ⁢
        
     
      C
      .
     
    
   
  
 



Subsequently, in the step (h10), whether the inequation i<n is established is determined; if yes, setting i=i+1 and turning to the step (h3); if no, turning to step (h11).
Subsequently, in the step (h11), the outlet temperature of each stand is obtained: Ti={17276, 17802, 186.59, 194.35, 206.33}° C.
Subsequently, in the step (i), the dynamic viscosity η0i of the emulsion in the roll gap of each stand is calculated by equation: η0i=b·exp(−a·Ti) (in the equation, a and b are dynamic viscosity parameters of the lubricant under the atmospheric pressure, wherein a=0.05, b=2.5), η0i={5.39, 5.46, 5.59, 5.69, 5.84} is obtained.
Subsequently, in the step (j), the thickness ξi of the oil film in the roll gap of each stand is calculated by the following calculation formula:



 
  
   ξ
   i
  
  =
  
   
    
     
      
       h
       
        0
        ⁢
        i
       
      
      +
      
       h
       
        1
        ⁢
        i
       
      
     
     
      2
      ⁢
      
       h
       
        0
        ⁢
        i
       
      
     
    
    ·
    
     k
     c
    
    ·
    
     
      3
      ⁢
      θ
      ⁢
      
       
        η
        
         0
         ⁢
         i
        
       
       (
       
        
         v
         ri
        
        +
        
         v
         
          0
          ⁢
          i
         
        
       
       )
      
     
     
      
       α
       i
      
      [
      
       1
       -
       
        e
        
         -
         
          θ
          (
          
           
            K
            i
           
           ⁢
           
            
             T
             
              0
              ⁢
              i
             
            
            
             
              h
              
               0
               ⁢
               i
              
             
             ·
             B
            
           
          
         
        
       
      
      ]
     
    
   
   -
   
    
     k
     rg
    
    ·
    
     (
     
      1
      +
      
       K
       rs
      
     
     )
    
    ·
    
     Ra
     
      ir
      ⁢
      0
     
    
    ·
    
     e
     
      
       
        -
        
         B
         Li
        
       
       ·
       
        L
        i
       
      
      ,
     
    
   
  
 



wherein, krg represents the coefficient of the strength of the surface longitudinal roughness of the work roll and the strip steel to entrain the lubricant, wherein krg=1.183, Krs represents the impression rate, i.e., the ratio of transmitting surface roughness of the work roll to the strip steel, wherein Krs=0.576, and S={0.784, 0.963, 2.101, 2.043, 1.326} um is obtained.
Subsequently, in the step (k), the emulsion concentration comprehensive optimization target function is calculated as follows:



 
  {
  
   
    
     
      
       
        F
        ⁡
        (
        X
        )
       
       =
       
        
         
          λ
          n
         
         ⁢
         
          
           ∑
           
            i
            =
            1
           
           n
          
          
           
            
             (
             
              
               ξ
               i
              
              -
              
               ξ
               
                0
                ⁢
                i
               
              
             
             )
            
            2
           
          
         
        
        +
        
         
          (
          
           1
           -
           λ
          
          )
         
         ⁢
         max
         ⁢
         
          
           ❘
           "\[LeftBracketingBar]"
          
          
           
            ξ
            i
           
           -
           
            ξ
            
             0
             ⁢
             i
            
           
          
          
           ❘
           "\[RightBracketingBar]"
          
         
        
       
      
     
    
    
     
      
       
        ξ
        i
        -
       
       <
       
        ξ
        i
       
       <
       
        ξ
        i
        +
       
      
     
    
   
   ;
  
 



wherein X={Ci} is an optimized variable, and λ=0.5 is a distribution coefficient, and F(X)=0.94 is obtained.
Subsequently, in the step (1), if the equation, F(X)=0.94<F0=1×1010, is established, then Ciy=Ci={4.2, 4.2, 4.2, 4.2, 4.2} %, F0=F(X)=0.94, and step (m) is performed.
Subsequently, in the step (m), whether the emulsion concentration Ci exceeds the range of the feasible region is determined, if yes, the step (n) is performed; if no, the step (g) is performed.
Subsequently, in the step (n), the optimal emulsion concentration set value Ciy={4.2, 4.2, 4.5, 4.6, 4.3} % is outputted, wherein Ciy is the value of Ci when the calculated value of F(X) is minimum in the feasible region.
Finally, during the whole rolling process, the control system of the cold continuous rolling mill set adjusts and controls the emulsion concentration of each stand respectively according to the optimal emulsion concentration set value obtained in the step (n).
To sum up, the technical solution of the invention changes the mode in the prior art that the emulsion in each stand adopts constant concentration control, takes the concentration of the emulsion in each stand as a variable to be optimized, and carries out comprehensive optimization control on the emulsion concentration, thus achieving suppressing the vibration of the rolling mill.
The method of the invention can be widely used in the field of controlling emulsion concentration of the cold continuous rolling mill set.

## Claims

1. A method of emulsion concentration optimization for a cold continuous rolling mill set, comprising the following steps:
(a) acquiring apparatus characteristic parameters of a cold continuous rolling mill set;
(b) acquiring critical rolling process parameters of a strip material;
(c) acquiring process parameters involved in the process of emulsion concentration optimization;
(d) setting an initial set value of an emulsion concentration comprehensive optimization target function for a cold continuous rolling mill set for achieving vibration suppression: F0=1.0×1010;
the executing order of steps (a) to (d) is not limited;
(e) calculating a bite angle α, of each stand;
(f) calculating a vibration determination index reference value ξ0i of each stand;
(g) setting an emulsion concentration Ci of each stand;
(h) calculating the outlet temperature Ti of a strip steel of each stand;
(i) calculating the dynamic viscosity η0i of an emulsion in a roll gap of each stand;
(j) calculating an oil film thickness ξi in the roll gap of each stand;
(k) calculating the emulsion concentration comprehensive optimization target function F(X);
(l) determining whether the inequation F(X)<F0 is established, if yes, then setting Ciy=Ci, F0=F(X) and turning to step (m); if no, directly turning to step (m);
(m) determining whether the emulsion concentration Ci exceeds the range of a feasible region, if yes, turning to step (n); if no, turning to step (g), wherein the feasible region refers to a region from 0 to the maximum emulsion concentration allowed by the cold continuous rolling mill set;
(n) outputting the optimal emulsion concentration set value Ciy, wherein Ciy is the value of Ci when the calculated value of F(X) is minimum in the feasible region; and
(o) adjusting and controlling the emulsion concentration of each stand according to the optimal emulsion concentration set value Ciy in the step (n) by a control system of the cold continuous rolling mill set; and
in each formula, i represents the stand ordinal number of the cold continuous rolling mill set.

2. The method of claim 1, wherein calculating the vibration determination index reference value ξ0i of each stand in the step (f) comprising the following steps:
(f1) calculating an over-lubricated oil film thickness critical value ξi+ of each stand as follows:

it is assumed that when



 
  
   
    
     γ
     i
    
    
     α
     i
    
   
   =
   
    A
    +
   
  
  ,
 



the roll gap is just in an over-lubricated status, wherein γi is a neutral angle of each stand, and A+ is an over-lubricated determining coefficient;
calculating to obtain




 
  
   u
   i
   +
  
  =
  
   
    1
    
     2
     ⁢
     
      (
      
       
        2
        ⁢
        
         A
         +
        
       
       -
       1
      
      )
     
    
   
   ⁢
   
    (
    
     
      
       
        Δ
        ⁢
        
         h
         i
        
       
       
        R
        i
        ′
       
      
     
     +
     
      
       
        T
        
         i
         ⁢
         0
        
       
       -
       
        T
        
         i
         ⁢
         1
        
       
      
      
       P
       i
      
     
    
    )
   
  
 



according to



 
  
   
    α
    i
   
   =
   
    
     
      
       Δ
       ⁢
       h
      
      i
     
     
      R
      i
      ′
     
    
   
  
  ,
 



wherein Δhi is the rolling reduction, Δhi=h0i−h1i, h0i is the inlet thickness of each stand, h1i is the outlet thickness of each stand, and Ri′ is the flattening radius of a work roll of the ith stand, and



 
  
   
    γ
    i
   
   =
   
    
     1
     2
    
    ⁢
    
     
      
       
        Δ
        ⁢
        
         h
         i
        
       
       
        R
        i
        ′
       
      
     
     [
     
      1
      +
      
       
        1
        
         2
         ⁢
         
          u
          i
         
        
       
       ⁢
       
        (
        
         
          
           
            Δ
            ⁢
            
             h
             i
            
           
           
            R
            i
            ′
           
          
         
         +
         
          
           
            T
            
             i
             ⁢
             0
            
           
           -
           
            T
            
             i
             ⁢
             1
            
           
          
          
           P
           i
          
         
        
        )
       
      
     
     ]
    
   
  
  ,
 


wherein T0i is the back tension of each stand, T1i is the front tension of each stand, and Pi is the rolling pressure of each stand;

then calculating to obtain the over-lubricated oil film thickness critical value of each stand ξi+ by



 
  
   ξ
   1
   +
  
  =
  
   
    1
    
     B
     i
    
   
   ⁢
   ln
   ⁢
   
    
     
      u
      i
      +
     
     -
     
      a
      i
     
    
    
     b
     i
    
   
  
 



according to the relationship formula between the friction coefficient and the oil film thickness: ui=ai+bi·eBi·ξi, wherein ai is the liquid friction influence coefficient, bi is the dry friction influence coefficient, and Bi is the friction coefficient attenuation index;
(f2) calculating an under-lubricated oil film thickness critical value ξi− of each stand as follows:
calculating to obtain




 
  
   u
   i
   -
  
  =
  
   
    1
    
     2
     ⁢
     
      (
      
       
        2
        ⁢
        
         A
         -
        
       
       -
       1
      
      )
     
    
   
   ⁢
   
    (
    
     
      
       
        Δ
        ⁢
        
         h
         i
        
       
       
        R
        i
        ′
       
      
     
     +
     
      
       
        T
        
         i
         ⁢
         0
        
       
       -
       
        T
        
         i
         ⁢
         1
        
       
      
      
       P
       i
      
     
    
    )
   
  
 



assuming that when



 
  
   
    
     γ
     i
    
    
     α
     i
    
   
   =
   
    A
    -
   
  
  ,
 



the roll gap is just in an under-lubricated status, wherein A− is the under-lubricated determining coefficient;
then calculating to obtain the under-lubricated oil film thickness critical value ξi− of each stand by




 
  
   
    ξ
    i
    -
   
   =
   
    
     1
     
      B
      i
     
    
    ⁢
    ln
    ⁢
    
     
      
       u
       i
       -
      
      -
      
       a
       i
      
     
     
      b
      i
     
    
   
  
  ;
 



and
(f3) calculating the vibration determination index reference value ξ0i by




 
  
   ξ
   
    0
    ⁢
    i
   
  
  =
  
   
    
     
      ξ
      i
      +
     
     +
     
      ξ
      i
      -
     
    
    2
   
   .

3. The method of claim 1, wherein calculating the outlet temperature Ti of the strip steel of each stand in the step (h) comprising the following steps:
(h1) calculating the outlet temperature Ti of the first stand:




 
  
   
    T
    1
   
   =
   
    
     T
     1
     Inlct
    
    +
    
     
      
       1
       -
       
        (
        
         
          ε
          1
         
         /
         4
        
        )
       
      
      
       1
       -
       
        (
        
         
          ε
          1
         
         /
         2
        
        )
       
      
     
     ·
     
      
       
        K
        1
       
       ⁢
       
        ln
        ⁡
        (
        
         1
         
          1
          -
          
           ε
           1
          
         
        
        )
       
      
      
       ρ
       ⁢
       S
       ⁢
       J
      
     
    
   
  
  ,
 


wherein T1Inlet is the inlet temperature of each stand,




 
  
   
    ε
    i
   
   =
   
    
     Δ
     ⁢
     
      h
      i
     
    
    
     h
     
      0
      ⁢
      i
     
    
   
  
  ,
 



Δhi=h0i−h1ih0i is the inlet thickness of each stand, h1i is the outlet thickness of each stand, ρ is the density of strip steel, S is the specific heat capacity of strip steel, J is the mechanical equivalent of heat, and Ki is the resistance to deformation of the strip steel of each stand;
(h2) setting i to 1;
(h3) setting Ti, 1 the temperature of the first section of strip steel behind the outlet of the ith stand, i.e., Ti, 1=Ti;
(h4) setting j to 2;
(h5) calculating the temperature Ti,j of the jth section of strip steel by the relationship between the temperature of the jth section and the temperature of the j−1th section shown by the following equation:




 
  
   
    T
    
     i
     ,
     j
    
   
   =
   
    
     
      -
      
       
        2
        ⁢
        
         k
         0
        
        ⁢
        
         w
         0.264
        
        ⁢
        
         exp
         ⁡
         (
         
          9.45
          -
          
           0.1918
           
            C
            i
           
          
         
         )
        
        ×
        1.163
        l
       
       
        
         v
         
          1
          ⁢
          i
         
        
        ⁢
        
         h
         
          1
          ⁢
          i
         
        
        ⁢
        ρ
        ⁢
        Sm
       
      
     
     ⁢
     
      
       T
       
        i
        ,
        
         j
         -
         1
        
       
       
        -
        0.213
       
      
      (
      
       
        T
        
         i
         ,
         
          j
          -
          1
         
        
       
       -
       
        T
        c
       
      
      )
     
    
    +
    
     T
     
      i
      ,
      
       j
       -
       1
      
     
    
   
  
  ,
 


wherein k0 is the influence coefficient of nozzle shape and spraying angle, w is the flow rate of the emulsion, l is a distance between stands, and the distance l between stands is equally divided into m section, the temperature in the section is represented by Ti, j, v1i is the outlet speed of each stand, h1i is the outlet thickness of each stand, ρ is the density of the strip steel, S is the specific heat capacity of the strip steel, Ti is the outlet temperature of each stand, and Tc is the temperature of the emulsion;
(h6) determining whether the inequation j<m is established, if yes, then setting j=j+1 and turning to the step (h5); if no, turning to step (h7);
(h7) obtaining the temperature Ti,m of the mth section via iterative calculation;
(h8) calculating the inlet temperature Ti+1Inlet of the i+1th stand:


Ti+1Inlet=Ti,m;


(h9) calculating the outlet temperature Ti+1Inlet of the i+1th stand:




 
  
   
    T
    
     i
     +
     1
    
   
   =
   
    
     T
     
      i
      +
      1
     
     Inlet
    
    +
    
     
      
       1
       -
       
        (
        
         
          ε
          
           i
           +
           1
          
         
         /
         4
        
        )
       
      
      
       1
       -
       
        (
        
         
          ε
          
           i
           +
           1
          
         
         /
         2
        
        )
       
      
     
     ·
     
      
       
        K
        
         i
         +
         1
        
       
       ⁢
       ln
       ⁢
       
        (
        
         1
         
          1
          -
          
           ε
           
            i
            +
            1
           
          
         
        
        )
       
      
      
       ρ
       ⁢
       SJ
      
     
    
   
  
  ;
 


(h10) determining whether the inequation i<n is established, if yes, setting i=i+1 and turning to the step (h3); if no, turning to step (h11); and
(h11) obtaining the outlet temperature Ti of each stand.

4. The method of claim 1, wherein the apparatus characteristic parameters of the cold continuous rolling mill set at least comprise:
the radius Ri of a work roll of each stand; the surface linear speed vri of a roll of each stand; the original roughness Rair0 of a work roll of each stand; the roughness attenuation coefficient BL of a work roll; the distance l between stands; and rolling kilometers Li of a work roll of each stand after roll change; wherein i=1, 2, . . . , n; i represents the ordinal number of the stand of the cold continuous rolling mill set, and n is the total number of the stands.

5. The method of emulsion concentration optimization for the cold continuous rolling mill set as claimed in claim 1, wherein the critical rolling process parameters of the strip material at least comprise:
the inlet thickness h0i of each stand; the outlet thickness h1i of each stand; the width B of the strip steel; the inlet speed v0i of each stand; the outlet speed v1i of each stand; the inlet temperature T1Inlet; the resistance to deformation Ki of the strip steel of each stand; the rolling pressure Pi of each stand; the back tension T0i of each stand; the front tension T1i of each stand; the concentration influence coefficient kc of the emulsion; the viscosity compression coefficient θ of a lubricant; the density ρ of the strip steel; the specific heat capacity S of the strip steel; the flow rate w of the emulsion; the temperature Tc of the emulsion; and the mechanical equivalent of heat J.

6. The method of claim 1, wherein the process parameters involved in the process of emulsion concentration optimization at least comprise: the over-lubricated oil film thickness critical value ξi+ and the corresponding friction coefficient ui+ of each stand; the under-lubricated oil film thickness critical value ξi− and the corresponding friction coefficient ui− of each stand; and the vibration determination index reference value ξ0i;
the rolling reduction Δhi=h0i−h1i; the rolling reduction ratio




 
  
   
    ε
    i
   
   =
   
    
     Δ
     ⁢
     
      h
      i
     
    
    
     h
     
      0
      ⁢
      i
     
    
   
  
  ;
 



the inlet temperature TiInlet and the outlet temperature Ti of each stand;
the distance l between stands, which is equally divided into m section, the temperature in the section represented by Ti, j, wherein 1≤j≤m and TiInlet=Ti−1,m;
the over-lubricated determining coefficient A+; and the under-lubricated determining coefficient A−.

7. The method of claim 1, wherein the calculation formula for calculating the bite angle αi of each stand is as follows:



 
  
   α
   i
  
  =
  
   
    
     Δ
     ⁢
     
      h
      i
     
    
    
     R
     i
     ′
    
   
  
 


wherein Ri′ is the flattening radius of the work roll of the ith stand, and is a process value in rolling pressure calculation.

8. The method of claim 1, wherein the calculation formula for calculating the dynamic viscosity η0i of the emulsion in the roll gap of each stand is as follows:


η0i=b·exp(−a·Ti)

wherein a and b are dynamic viscosity parameters of the lubricant under atmospheric pressure.

9. The method of claim 1, wherein the calculation formula for calculating the oil film thickness in the roll gap of each stand is as follows:



 
  
   
    ξ
    i
   
   =
   
    
     
      
       
        h
        
         0
         ⁢
         i
        
       
       +
       
        h
        
         1
         ⁢
         i
        
       
      
      
       2
       ⁢
       
        h
        
         0
         ⁢
         i
        
       
      
     
     ·
     
      k
      c
     
     ·
     
      
       3
       ⁢
       θ
       ⁢
       
        
         η
         
          0
          ⁢
          i
         
        
        (
        
         
          v
          ri
         
         +
         
          v
          
           0
           ⁢
           i
          
         
        
        )
       
      
      
       
        α
        i
       
       [
       
        1
        -
        
         e
         
          -
          
           θ
           ⁡
           (
           
            
             K
             i
            
            ⁢
            
             
              T
              
               0
               ⁢
               i
              
             
             
              
               h
               
                0
                ⁢
                i
               
              
              ·
              B
             
            
           
           )
          
         
        
       
       ]
      
     
    
    -
    
     
      k
      rg
     
     ·
     
      (
      
       1
       +
       
        K
        rs
       
      
      )
     
     ·
     
      Ra
      
       ir
       ⁢
       0
      
     
     ·
     
      e
      
       
        -
        
         B
         Li
        
       
       ·
       
        L
        i
       
      
     
    
   
  
  ,
 


wherein h0i is the inlet thickness of each stand, h1i is the outlet thickness of each stand, kc is the emulsion concentration influence coefficient, θ is the viscosity compression coefficient of the lubricant, Ki is the resistance to deformation of the strip steel of each stand, η0i is the dynamic viscosity of the emulsion in the roll gap of each stand, v0i is the inlet speed of each stand, vri is the surface linear speed of a roll of each stand, T0i is the back tension of each stand, B is the width of the strip steel, krg represents the coefficient of the strength of entrainment of lubricant by the longitudinal surface roughness of the work roll and the strip steel, the value of which is from 0.09 to 0.15, Krs represents the impression rate, i.e., the ratio of transmitting surface roughness of the work roll to the strip steel, Rair0 is the original roughness of a work roll of each stand, BL is the roughness attenuation coefficient of the work roll, and Li is the rolling kilometers of a work roll of each stand after roll change.

10. The method of in claim 1, wherein the emulsion concentration comprehensive optimization target function is calculated according to the following formula:



 
  {
  
   
    
     
      
       
        F
        ⁡
        (
        X
        )
       
       =
       
        
         
          λ
          n
         
         ⁢
         
          
           ∑
           
            i
            =
            1
           
           n
          
          
           
            
             (
             
              
               ξ
               i
              
              -
              
               ξ
               
                0
                ⁢
                i
               
              
             
             )
            
            2
           
          
         
        
        +
        
         
          (
          
           1
           -
           λ
          
          )
         
         ⁢
         max
         ⁢
         
          
           ❘
           "\[LeftBracketingBar]"
          
          
           
            ξ
            i
           
           -
           
            
             ζ
             
              0
              ⁢
              i
             
            
            
             ❘
             "\[LeftBracketingBar]"
            
           
          
         
        
       
      
     
    
    
     
      
       
        ξ
        i
        -
       
       <
       
        ξ
        i
       
       <
       
        ξ
        i
        +
       
      
     
    
   
   ,
  
 


wherein X={Ci} is an optimization variable, and λ is a distribution coefficient.