Abstract:
The present invention provides a coordinated coil scheduling method of multiple production lines in cold rolling area in steel plant, comprising: acquiring the information about the candidate coils to be scheduled in each of the multiple production lines in cold rolling area; building a coordinated coil scheduling model of multiple lines; obtaining the initial coordinated coil schedule of multiple lines by using a heuristic algorithm; improving the initial coordinated coil schedule in real time; and sending the improved coordinated coil schedule to the ACS of each production line. In this method, besides different production technological requirements of each coil, the cohesion relationship of different production lines are taken into account to make material supply to each equipment reasonable, and the smooth production process of coils in each of the lines is guaranteed with the changeover minimized, such that the whole efficiency of the manufacturing process is improved.

Description:
TECHNICAL FIELD 
       [0001]    The present invention belongs to the field of information technology for manufacturing in metallurgy and relates to the field of automation technology, and more specifically, to a coordinated coil scheduling method of multiple production lines in cold rolling area in steel plant. 
       BACKGROUND ART 
       [0002]    Coordination scheduling of metal materials for multiple lines in cold rolling area is an important part in the production operations for metallurgic enterprises. The effectiveness of the coordinated schedule has directly relationship on the stability of equipments and the operation efficiency in the multiple production lines in cold rolling area, and, on the quality of the final product and the production cost. 
         [0003]    Taking the steel plant as an example, the cold rolling area comprises multiple production stages, such as, acid pickling, acid rolling, continuous annealing, hot galvanizing, hot galvanizing and aluminizing, finishing, electro-tinning, crosscutting, rewinding and etc. The hot rolled coils must pass the continuous acid pickling line to clean the iron oxide on the surface of the strip steel so as to perform cold rolling operation and the consequent surface processes smoothly. The hot rolled coils will be processed to be cold rolled coils in acid-rolling line, which can be directly sold, however, cold hardening caused by continuous cold deformation will lead to the increasing the strength and hardness of the coils and reducing the toughness, consequently, deteriorating the stamping characteristic. In order to improve the mechanical characteristic of the coils, the coils need to be annealed to softening the steel, smooth the surface, reduce the hardness, get a better resistance against deformation and crack, then a uniform distribution of structures and composites and a better performance of the material will achieve. In addition, the annealed coils also pass the electro-tinning or galvanizing line to form a firm, bright, compact, uniform and continuous alloy coating, which will realize good anti-corrosion and tarnish resistance characteristics. The tin-plated coils can be sold after being cut crossly according to the orders.  FIG. 1  shows the arrangement of the respective production lines in the cold rolling area. The coordinated coil scheduling among multiple lines in cold rolling area has the features including multiple production lines arranged in parallel and in series concurrently, complex logistics, existence of both direct flows and cross flows simultaneously. Further, the coils are processed continuously in each production line, the processing of the coil cannot be started in the downstream production line until it is finished in the upstream production line, thus, the coil will be reprocessed or become wasted material, if its terrible defect incurred in the upstream production line is not be repaired in the downstream production line. Therefore, it must be guaranteed that the production in each production line is stably, and every product must be processed precisely in each line, so as to satisfy the customers&#39; requirements. Reasonable coordinated coil schedule of multiple lines can equilibrate the production rhythm of the upstream and downstream of production lines in order to keep the production in each production line continuously and unblocked, reduce the changeover times, and keep the reasonable utilization of the production capacity and stably production in each production line. 
         [0004]    Currently, the coil schedules in acid rolling, continuous annealing, hot galvanizing, hot galvanizing and aluminizing are made by the planner in each production line and the coordinated planner, in practice. The coordinated planner determines the production volume in each line per day, based on information such as the orders&#39; requirements, the production capacity, the current state of lines and the coils in inventory etc. The planner in each production line determines the coil schedule based on the determined production volume in the lines each day. However, the practical production conditions are various, and often involve several decision makers, so that the current manual planning method could not effectively coordinate the production rhythms among the lines and determining the schedule with good performance, thus causing conflicts among the lines and increasing the changeover times, which couldn&#39;t guarantee the production stability and the quality of the product. On the other hand, the manual plan is strongly dependent on the personal experience of the planner, and the planner in each production line makes the coil schedule according to the scheduling requirements in the respective line. Each line may have different scheduling requirements to be taken into account thus causing conflicts between the coil schedules in different lines. If a coil schedule is made with only consideration of the scheduling requirements of a line, other lines may be negatively influenced and suffer from frequent changeover and material shortage, resulting in a high changeover costs. Therefore, it is important to determine the coil schedule for each production line in cold rolling area in a view of whole supply chain to coordinate coil schedules comprehensively, so as to optimize the overall production effect of all the lines in the cold rolling area. 
       SUMMARY OF INVENTION 
       [0005]    In order to overcome the above mentioned drawbacks in the prior art, the present invention provides a coordinated coil scheduling method of multiple lines in cold rolling area in steel plant. 
         [0006]    The technical solution of the present invention is: 
         [0007]    a coordinated coil scheduling method of multiple lines in cold rolling area in steel plant, comprising: 
         [0008]    Step 1: acquiring the information about the candidate coils to be scheduled for each of multiple production lines in the cold rolling area; 
         [0009]    Step 2: building a coordinated coil scheduling model of multiple lines in cold rolling section in steel plant; 
         [0010]    Step 2-1: according to production technological requirements of the production lines, calculating abrasion penalties and changeover penalties incurred by equipment adjustment due to different production requirements, comprising: abrasion penalties of equipment incurred by equipment adjustment due to different widths of two adjacent coils in a production line, abrasion penalties of equipment incurred by equipment adjustment due to different thicknesses of two adjacent coils in a production line, abrasion penalties of equipment incurred by equipment adjustment due to different surface requirements of two adjacent coils in a production line, abrasion penalties of equipment incurred by equipment adjustment due to different annealing temperatures of two adjacent coils in a production line, abrasion penalties of equipment incurred by equipment adjustment due to different weights of the surface plates of two coils in a production line, and abrasion penalties of equipment incurred by equipment adjustment due to different post-processing requirements of two coils in a production line; 
         [0011]    Step 2-2: among the coil schedules corresponding to the production lines in the cold rolling supply chain composed by multiple production lines, an objective function of coordinated coil scheduling of multiple lines in the cold rolling area is built as minimizing the abrasion penalties and changeover penalties incurred by equipment adjustment due to different production requirements, 
         [0000]    
       
         
           
             
               i 
               . 
               e 
               . 
             
             , 
             
               
 
             
              
             
               
                 minimize 
                  
                 
                   
                     ∑ 
                     
                       k 
                       ∈ 
                       M 
                     
                     
                         
                     
                   
                    
                   
                     
                       ∑ 
                       
                         j 
                         ∈ 
                         
                           N 
                           k 
                         
                       
                     
                      
                     
                       
                         ∑ 
                         
                           i 
                           ∈ 
                           
                             N 
                             k 
                           
                         
                       
                        
                       
                         
                           F 
                           ijk 
                         
                          
                         
                           x 
                           ijk 
                         
                       
                     
                   
                 
               
               ; 
             
           
         
       
     
         [0012]    where a decision variable of the objective function is represented as x ijk , which indicates whether the coil i is processed immediately before the coil j, in the kth production line, if yes, the decision variable x ijk  is assigned as 1, and if not, it is assigned as 0; 
         [0013]    Step 2-3: constraints comprise the constraints for each of the multiple lines and the cohesion constraints among the production lines; 
         [0014]    where the constraints for each of the multiple lines include: (1) for each of all actual coils, there must be exactly one coil being processed immediately before it; (2) for each of all the actual coils, there must be exactly one coil being processed immediately after it; (3) the changeover in width, thickness, surface requirement, and annealing temperature of the lines should not exceed the given upper limits; (4) it should be guaranteed that the determination to the processing time of the coils in each line will not incur idle time between any two adjacent processed coils; (5) for each of the coils, it should be processed exactly once in a production line along its required processing path; 
         [0015]    where the cohesion constraints among the production lines imply the cohesion relationship between the production lines that the coil cannot be processed in a downstream production line until it is finished in the upstream production line; 
         [0016]    Step 3: an initial coordinated coil schedule of multiple production lines in the cold rolling area is obtained by using heuristic algorithm, which includes an initial coil schedule in acid rolling line, an initial coil schedule in a continuous annealing line, an initial coil schedule in a hot galvanizing line, and an initial coil schedule in a hot aluminizing and galvanizing line; 
         [0017]    Step 3-1: Determine the initial coil schedule in acid rolling line: an initial coil schedule in the acid rolling line is obtained by sequencing the coils to be processed in a descending order of the surface requirement code, and then sequencing the coils with the same surface requirement in a descending order of their widths, and then sequencing the coils with the same surface requirement and the same width in a descending order of their post-processing codes, and then sequencing the coils with the same surface requirement and the same width and the same post-processing code in an descending order of their annealing temperatures; 
         [0018]    Step 3-2: calculating the starting time and the finishing time of each coil in the initial coil schedule in the acid rolling line; 
         [0019]    Step 3-3: determine an initial coil schedule in the continuous annealing line: an initial annealing schedule of coils in the material inventory of continuous annealing line is obtained by sequencing the coils to be annealed in a descending order of their widths, and then sequencing the coils with the same width in a descending order of their annealing temperatures, and then sequencing the coils with the same width and the same annealing temperature in a descending order of their thicknesses, a coil schedule in continuous annealing line is obtained by selecting the coils to be annealed from the coil schedule in the acid rolling line, and sequentially inserting them after the initial annealing schedule of coils in the material inventory of continuous annealing line; 
         [0020]    Step 3-4: determine an initial coil schedule in the hot galvanizing line: an initial schedule of the coils in the material inventory of hot galvanizing line is obtained, by sequencing the coils in the material inventory of hot galvanizing line in a descending order of their post-processing codes; and then sequencing the coils with the same post-processing code in a descending order of their widths; and then sequencing the coils with the same width and the same post-processing code in a descending order of their thicknesses; and then sequencing the coils with the same post-processing code, and the same width and the same thickness in a descending order of their annealing temperatures; an initial coil schedule in the hot galvanizing line is obtained by selecting the coils to be processed in hot galvanizing line from the initial coil schedule in acid rolling line, sequentially inserting them after the initial schedule of the coils in the material inventory of hot galvanizing line; 
         [0021]    Step 3-5: determine an initial coil schedule in the hot aluminizing and galvanizing line: an initial schedule of coils in material inventory of the hot aluminizing and galvanizing line is obtained, by sequencing the coils in material inventory of the hot aluminizing and line galvanizing in a descending order of their post-processing codes; and then sequencing the coils with the same post-processing code in a descending order of their widths; and then sequencing the coils with the same post-processing code and the same width in a descending order of their thicknesses; and then sequencing the coils with the same post-processing code, the same width and the same thickness in a descending order of their annealing temperatures; an initial coil schedule in the hot aluminizing and galvanizing line is obtained, by selecting the coils to be processed in hot aluminizing and galvanizing line from the initial coil schedule in acid rolling line, sequentially inserting them after the initial schedule of coils in material inventory of the hot aluminizing and galvanizing line; 
         [0022]    Step 3-6: calculating the starting time of each coil in the initial coil schedule in the continuous annealing line, the starting time of each coil in the initial coil schedule in the hot galvanizing line, and the starting time of each coil in the initial coil schedule in the hot galvanizing and aluminizing line; 
         [0023]    Step 3-7: recording an initial coordinated coil schedule of multiple production lines in the cold rolling area, and calculating the objective function value corresponding to the initial coordinated coil schedule of multiple lines in the cold rolling area; 
         [0024]    Step 4: improving in real time the initial coordinated coil schedule of multiple lines in the cold rolling area; 
         [0025]    Step 5: sending the improved coordinated coil schedule of multiple lines in the cold rolling area to an automatic control systems (ACS) of each production lines in the cold rolling area, and realizing the coordinated coil scheduling of multiple lines in the cold rolling area. 
         [0026]    the step 4 improves in real time the initial coordinated coil schedule of multiple lines in the cold rolling area, comprises the following steps: 
         [0027]    Step 4-1: initializing population of the coordinated coil schedule of multiple lines in the cold rolling area, and setting a given maximum iteration number and a given maximum iteration number without improvement, population of the coordinated coil schedule of multiple lines in the cold rolling area comprises the initial coordinated coil schedule of multiple lines in the cold rolling area, and other coil schedules in the population of the coordinated coil schedule of multiple lines in the cold rolling area are generated randomly; and during the improvement, an adaptive value of each coordinated coil schedule is calculated according to the way for calculating the objective function of each coordinated coil schedule; 
         [0028]    Step 4-2: judging whether the given maximum iteration number is reached, and if it is reached, then taking the schedule with the best adaptive value in the current population as the optimal schedule, that is, the improved coordinated coil schedule of multiple lines in the cold rolling area, and if not, go to step 4-3; 
         [0029]    Step 4-3: judging whether the given maximum iteration number without improvement is reached, and if it is reached, then taking the schedule with the best adaptive value in the current population as the optimal schedule, that is, the improved coordinated coil schedule of multiple lines in the cold rolling area, and if not, go to step 4-4; 
         [0030]    Step 4-4: if the individual with the best adaptive value in the population is not updated, the maximum iteration number without improvement is added by 1; 
         [0031]    Step 4-5: performing a mutation operation to the population; 
         [0032]    Step 4-6: performing a crossover operation to the population; 
         [0033]    Step 4-7: performing a selecting operation on the son generation in the population; 
         [0034]    Step 4-8: calculating the adaptive value for each individual in the population, that is, calculating the objective function value of coordinated coil schedule of multiple lines in the cold rolling area corresponding to the current population; 
         [0035]    Step 4-9: performing an increment operation to the current iteration number and updating the population, and then returning to step 4-2. 
         [0036]    The step 4-5 performs a mutation operation to the population, comprises the following: 
         [0000]    
       
         
           
             
               V 
               i 
               t 
             
             = 
             
               
                 X 
                 a 
                 
                   t 
                   - 
                   1 
                 
               
               ⊕ 
               
                 F 
                 ⊗ 
                 
                   ( 
                   
                     
                       X 
                       
                         b 
                          
                         
                             
                         
                       
                       
                         t 
                         - 
                         1 
                       
                     
                     - 
                     
                       X 
                       c 
                       
                         t 
                         - 
                         1 
                       
                     
                   
                   ) 
                 
               
             
           
         
       
       
         
           
             
               Δ 
               i 
               t 
             
             = 
             
               
                 
                   F 
                   ⊗ 
                   
                     ( 
                     
                       
                         X 
                         b 
                         
                           t 
                           - 
                           1 
                         
                       
                       - 
                       
                         X 
                         c 
                         
                           t 
                           - 
                           1 
                         
                       
                     
                     ) 
                   
                 
                 ⇔ 
                 
                   δ 
                   ij 
                   t 
                 
               
               = 
               
                 { 
                 
                   
                     
                       
                         
                           
                             
                               x 
                               bj 
                               
                                 t 
                                 - 
                                 1 
                               
                             
                             - 
                             
                               x 
                               cj 
                               
                                 t 
                                 - 
                                 1 
                               
                             
                           
                         
                         
                           
                             
                               if 
                                
                               
                                   
                               
                                
                               
                                 rand 
                                 ( 
                                 ) 
                               
                             
                             &lt; 
                             F 
                           
                         
                       
                       
                         
                           0 
                         
                         
                           otherwise 
                         
                       
                     
                      
                     
                       
 
                     
                      
                     
                       Δ 
                       i 
                       t 
                     
                   
                   = 
                   
                     
                       
                         ( 
                         
                           
                             δ 
                             
                               i 
                               , 
                               0 
                             
                             t 
                           
                           , 
                           
                             δ 
                             
                               i 
                               , 
                               1 
                             
                             t 
                           
                           , 
                           … 
                            
                           
                               
                           
                           , 
                           
                             δ 
                             
                               i 
                               , 
                               n 
                             
                             t 
                           
                         
                         ) 
                       
                        
                       
                         
 
                       
                        
                       
                         V 
                         i 
                         t 
                       
                     
                     = 
                     
                       
                         
                           
                             X 
                             a 
                             
                               t 
                               - 
                               1 
                             
                           
                           ⊕ 
                           
                             Δ 
                             i 
                             t 
                           
                         
                         ⇔ 
                         
                           v 
                           ij 
                           t 
                         
                       
                       = 
                       
                         
                           
                             x 
                             aj 
                             
                               t 
                               - 
                               1 
                             
                           
                           ⊕ 
                           
                             δ 
                             ij 
                             t 
                           
                         
                         = 
                         
                           mod 
                            
                           
                             ( 
                             
                               
                                 ( 
                                 
                                   
                                     x 
                                     aj 
                                     
                                       t 
                                       - 
                                       1 
                                     
                                   
                                   + 
                                   
                                     δ 
                                     ij 
                                     t 
                                   
                                   + 
                                   n 
                                 
                                 ) 
                               
                               , 
                               n 
                             
                             ) 
                           
                         
                       
                     
                   
                 
               
             
           
         
       
     
         [0037]    where V i   t  is the ith variant chromosome in the tth generation, and v ij   t  is the jth individual in V i   t , and X a   t−1 , X b   t−1  and X c   t−1  are three different coordinated coil schedules of multiple lines in the cold rolling area in the population of a t−1th generation, x ij   t−1  is a jth individual in X i   t−1 , F is a scaling factor, Δ i   t  is the ith intermediate variant chromosome in the mutation process of the tth generation, δ ij   t  is a jth individual in the intermediate variant chromosome Δ i   t  in the mutation process of the tth generation, and mod means modulus. 
         [0038]    Step 4-6 performs the crossover operation to the population, comprises the following: 
         [0039]    generating a trial chromosome U i   t , i=1, 2, . . . , PS, and u ij   t  is an individual in the chromosome U i   t ; 
         [0000]    
       
         
           
             
               u 
               ij 
               t 
             
             = 
             
               { 
               
                 
                   
                     
                       v 
                       ij 
                       t 
                     
                   
                   
                     
                       
                         
                           if 
                            
                           
                               
                           
                            
                           
                             rand 
                             ( 
                             ) 
                           
                         
                         &lt; 
                         
                           CR 
                            
                           
                               
                           
                            
                           or 
                            
                           
                               
                           
                            
                           j 
                         
                       
                       = 
                       
                         D 
                         j 
                       
                     
                   
                 
                 
                   
                     
                       x 
                       ij 
                       
                         t 
                         - 
                         1 
                       
                     
                   
                   
                     otherwise 
                   
                 
               
             
           
         
       
     
         [0040]    where CR is in [0,1] which is the crossover parameter to ensure the variety of the s of the population, rand( ) generates a group of uniformly distributed randomly numbers between 0 and 1, D j  is a number randomly selected from the set{1, 2, . . . , n}, so as to guarantee that at least one dimension of U i   t  is different from the ith coordinated coil schedule of multiple lines in the cold rolling area in the t−1th generation X i   i−1 ; 
         [0041]    CR is a self-updating number as the iteration number increases, and is updated as follows: 
         [0000]    
       
         
           
             CR 
             = 
             
               
                 CR 
                 min 
               
               + 
               
                 t 
                 * 
                 
                   
                     ( 
                     
                       
                         CR 
                         max 
                       
                       - 
                       
                         CR 
                         min 
                       
                     
                     ) 
                   
                   T 
                 
               
             
           
         
       
     
         [0042]    where CR min  is a minimum crossover probability, CR max  is a maximum crossover probability, t is the current iteration number, and T is a maximum iteration number. 
       Advantageous Effects 
       [0043]    (1) Not only different production technological requirements but also the cohesion relationship between different production lines in cold rolling area are taken in to account in the method of the present invention, so as to equilibrate material supply of the lines, and to guarantee a steady and smoothly production of coils in the lines and finally to improve the overall efficiency of the production in cold rolling area. 
         [0044]    (2) Compared to manual production schedule, the total changeover times are reduced, and the frequency of equipment adjustment is reduced, thereby improving the quality of the product and production efficiency. 
         [0045]    (3) Since there are similarities between the process in steel plant and that of production of other metals, the method of the present invention can be easily applied to production of other kinds of metals. 
     
    
     
       BRIEF DESCRIPTION OF DRAWINGS 
         [0046]      FIG. 1  illustrates the arrangement of lines in cold rolling area in steel plant in the embodiments of the present invention; 
           [0047]      FIG. 2  illustrates the flowchart of coordinated coil scheduling method of multiple lines in cold rolling area in steel plant in the embodiments of the present invention; and 
           [0048]      FIG. 3  illustrates the flowchart of real time improvement method for the initial coordinated coil schedule of multiple lines in cold rolling area in steel plant in the embodiments of the present invention. 
       
    
    
     DETAILED DESCRIPTION OF THE EMBODIMENTS 
       [0049]    The detailed description of the present invention is given here after with reference to the accompanied drawings. 
         [0050]    The method of the present invention is implemented in a coordinated coil scheduling system of multiple lines in cold rolling area in steel plant, comprises at least one personal computer (PC), at least one electric cable interface/optic fiber interface/telephone line interface, and at least one router. 
         [0051]    The PC is provided therein with a condition-setting module, an initial coordinated scheduling module and a real time improving module. 
         [0052]    The condition-setting module is configured to set processing parameters of the processing environment in the lines in cold rolling area, and describe the details of the processes and determine the objective and constraints for adjustment of the processes. 
         [0053]    The initial coordinated scheduling module is configured to sequence the coils based on the current supply of coils and obtain the initial coordinated coil schedule of multiple lines in cold rolling area. 
         [0054]    The real time improving module improves the current coordinated coil schedule in real time and produces a better schedule. 
         [0055]    The PC is configured to be connected to the front end of each of the lines in cold rolling area through a network and an internal server. 
         [0056]    The electric cable interface/optical fiber interface/telephone line interface is configured for communications between the information management platform within the enterprise, the coordinated coil scheduling system of multiple lines in cold rolling area and the automatic control system in work site, for the purpose of optimizing the production process and improving the quality of product. 
         [0057]    The software embedded in the PC includes a Windows operating system as a supporting platform, a Microsoft SQL Server 2000 database system supporting data management and configuring the information communicating ports. 
         [0058]    The arrangement of lines in cold rolling area in steel plant in the embodiment is illustrated as in  FIG. 1 , including an acid rolling line, a continuous annealing line, a hot galvanizing line and a hot aluminizing and galvanizing line. 
         [0059]    A coordinated coil scheduling method of multiple lines in cold rolling area is illustrated as in  FIG. 2 , including the following steps: 
         [0060]    Step 1: acquiring information about the coils to be scheduled for the lines in cold rolling area; 
         [0061]    downloading data information about the coils from the ERP data server of the enterprise (involved fields including the coil number, the status of the material, width, thickness, actual weight, metallurgic specifications, specifications of product, the order number, the repairing type, the downstream lines, theoretical weight, the mark of the coil, stock, original stock, the steel grade, mother coil, production date, blocking reason, reason for surplus material, processing ways, due date, etc), so as to obtain the operating objects for the system and the scheduling method, at the same time, the information about the coil is added/deleted one-by-one according to current coil consumption in inventory and coil supply from the upstream lines, and the coil information is amended in order to confirm the coil information to be taken into account in the production scheduling; then proceeding into an automatic optimized scheduling stage. In the embodiment, the automatic optimized scheduling may be used for scheduling 16 coils in inventory. 
         [0000]    
       
         
               
             
               
               
               
               
               
               
               
               
             
               
               
               
               
               
               
               
               
             
           
               
                 TABLE 1 
               
             
             
               
                   
               
               
                 Information About Coils To Be Scheduled 
               
             
          
           
               
                   
                 Surface 
                   
                 Thick- 
                 Annealing 
                   
                   
                 Located 
               
               
                 Coil 
                 require- 
                 Width 
                 ness 
                 Temper- 
                 Post- 
                 processing 
                 produc- 
               
               
                 No. 
                 ment 
                 (mm) 
                 (mm) 
                 ature 
                 processing 
                 route 
                 tion line 
               
               
                   
               
             
          
           
               
                 1 
                 1 
                 926 
                 2 
                 580 
                 — 
                 Acid rolling −&gt; 
                 Acid 
               
               
                   
                   
                   
                   
                   
                   
                 Continuous 
                 rolling 
               
               
                   
                   
                   
                   
                   
                   
                 annealing 
               
               
                 2 
                 1 
                 860 
                 2 
                 580 
                 — 
                 Acid rolling −&gt; 
                 Acid 
               
               
                   
                   
                   
                   
                   
                   
                 Continuous 
                 rolling 
               
               
                   
                   
                   
                   
                   
                   
                 annealing 
               
               
                 3 
                 1 
                 860 
                 2 
                 610 
                 — 
                 Acid rolling −&gt; 
                 Acid 
               
               
                   
                   
                   
                   
                   
                   
                 Continuous 
                 rolling 
               
               
                   
                   
                   
                   
                   
                   
                 annealing 
               
               
                 4 
                 1 
                 937 
                 2.5 
                 620 
                 — 
                 Acid rolling −&gt; 
                 Acid 
               
               
                   
                   
                   
                   
                   
                   
                 Continuous 
                 rolling 
               
               
                   
                   
                   
                   
                   
                   
                 annealing 
               
               
                 5 
                 1 
                 800 
                 0.18 
                 680 
                 — 
                 Acid rolling −&gt; 
                 Acid 
               
               
                   
                   
                   
                   
                   
                   
                 Continuous 
                 rolling 
               
               
                   
                   
                   
                   
                   
                   
                 annealing 
               
               
                 6 
                 1.2 
                 800 
                 0.18 
                 620 
                 Anti-finger 
                 Acid rolling −&gt; 
                 Acid 
               
               
                   
                   
                   
                   
                   
                 plating 
                 hot galvanizing 
                 rolling 
               
               
                 7 
                 1.2 
                 800 
                 0.18 
                 650 
                 Anti-finger 
                 Acid rolling −&gt; 
                 Acid 
               
               
                   
                   
                   
                   
                   
                 plating 
                 hot galvanizing 
                 rolling 
               
               
                 8 
                 1.2 
                 896 
                 0.18 
                 580 
                 Phospha- 
                 Acid rolling −&gt; 
                 Acid 
               
               
                   
                   
                   
                   
                   
                 tization 
                 hot aluminizing 
                 rolling 
               
               
                   
                   
                   
                   
                   
                   
                 and galvanizing 
               
               
                 9 
                 1.2 
                 896 
                 0.18 
                 620 
                 Phospha- 
                 Acid rolling −&gt; 
                 Acid 
               
               
                   
                   
                   
                   
                   
                 tization 
                 hot aluminizing 
                 rolling 
               
               
                   
                   
                   
                   
                   
                   
                 and galvanizing 
               
               
                 10 
                 0.5 
                 896 
                 0.18 
                 620 
                 — 
                 Acid rolling −&gt; 
                 Continuous 
               
               
                   
                   
                   
                   
                   
                   
                 Continuous 
                 annealing 
               
               
                   
                   
                   
                   
                   
                   
                 annealing 
               
               
                 11 
                 0.5 
                 896 
                 0.18 
                 680 
                 — 
                 Acid rolling −&gt; 
                 Continuous 
               
               
                   
                   
                   
                   
                   
                   
                 Continuous 
                 annealing 
               
               
                   
                   
                   
                   
                   
                   
                 annealing 
               
               
                 12 
                 0.5 
                 896 
                 0.18 
                 620 
                 — 
                 Acid rolling −&gt; 
                 Continuous 
               
               
                   
                   
                   
                   
                   
                   
                 Continuous 
                 annealing 
               
               
                   
                   
                   
                   
                   
                   
                 annealing 
               
               
                 13 
                 0.5 
                 1024 
                 2.3 
                 — 
                 — 
                 Acid rolling −&gt; 
                 Acid 
               
               
                   
                   
                   
                   
                   
                   
                 Directly sale 
                 rolling 
               
               
                 14 
                 0.5 
                 1024 
                 2.3 
                 — 
                 — 
                 Acid rolling −&gt; 
                 Acid 
               
               
                   
                   
                   
                   
                   
                   
                 Directly sale 
                 rolling 
               
               
                 15 
                 0.5 
                 1225 
                 2.3 
                 — 
                 — 
                 Acid rolling −&gt; 
                 Acid 
               
               
                   
                   
                   
                   
                   
                   
                 Directly sale 
                 rolling 
               
               
                 16 
                 0.5 
                 1225 
                 2.3 
                 — 
                 — 
                 Acid rolling −&gt; 
                 Acid 
               
               
                   
                   
                   
                   
                   
                   
                 Directly sale 
                 rolling 
               
               
                   
               
             
          
         
       
     
         [0062]    Step 2: building a model for coordinated coil scheduling of multiple lines in cold rolling area in steel plant; 
         [0063]    Step 2-1: according to production technological requirements of the production lines, calculating abrasion penalties and changeover penalties incurred by equipment adjustment due to different production requirements, comprising: abrasion penalties of equipment incurred by equipment adjustment due to different widths of two adjacent coils in a production line, abrasion penalties of equipment incurred by equipment adjustment due to different thicknesses of two adjacent coils in a production line, abrasion penalties of equipment incurred by equipment adjustment due to different surface requirements of two adjacent coils in a production line, abrasion penalties of equipment incurred by equipment adjustment due to different annealing temperatures of two adjacent coils in a production line, abrasion penalties of equipment incurred by equipment adjustment due to different weights of the surface plates of two coils in a production line, and abrasion penalties of equipment incurred by equipment adjustment due to different post-processing requirements of two coils in a production line; 
         [0064]    (1) abrasion penalty F ijk   w  of equipment incurred by equipment adjustment due to width changeover of coils i and j in line k: 
         [0000]        F   ijk   w   =K   k   w   ·|w   i   −w   j   |, i,j∈N, k∈M    (1)
 
         [0065]    where N represents the set of all coils to be processed, N={0, 1, . . . , n, n+1}, and 0 and n+1 represent dummy coils, n represents the number of actual coils to be scheduled, M represents the set of all production lines, M={0, . . . , m−1}, m represents the number of the production lines, K k   w  represents the unit weight of penalty incurred by width changeover on the kth line, and w i  and w j  represent the width of the coil i and the width of the coil j, respectively. 
         [0066]    (2) abrasion penalty F ijk   h  of equipment incurred by equipment adjustment due to thickness changeover of coils i and j in line k: 
         [0000]        F   ijk   h   =K   k   h   ·|h   i   −h   j   |, i,j∈N, k∈M    (2)
 
         [0067]    where K h   k  represents the unit weight of penalty incurred by thickness changeover in the kth line, and h i  and h j  represent the thickness of the coil i and the thickness of the coil j, respectively. 
         [0068]    (3) abrasion penalty F ijk   R  of equipment incurred by equipment adjustment due to surface requirement changeover of coils i and j in line k: 
         [0000]    
       
         
           
             
               
                 
                   
                     F 
                     ijk 
                     R 
                   
                   = 
                   
                     { 
                     
                       
                         
                           
                             
                               K 
                               k 
                               
                                 R 
                                  
                                 
                                     
                                 
                                  
                                 1 
                               
                             
                           
                           
                             
                               
                                 R 
                                 i 
                               
                               ≤ 
                               
                                 1 
                                  
                                 
                                     
                                 
                                  
                                 and 
                                  
                                 
                                     
                                 
                                  
                                 
                                   R 
                                   j 
                                 
                               
                               &gt; 
                               1 
                             
                           
                         
                         
                           
                             
                               K 
                               k 
                               
                                 R 
                                  
                                 
                                     
                                 
                                  
                                 2 
                               
                             
                           
                           
                             
                               
                                 R 
                                 j 
                               
                               ≤ 
                               
                                 1 
                                  
                                 
                                     
                                 
                                  
                                 and 
                                  
                                 
                                     
                                 
                                  
                                 
                                   R 
                                   i 
                                 
                               
                               &gt; 
                               1 
                             
                           
                         
                       
                       , 
                       
                           
                       
                        
                       i 
                       , 
                       
                         j 
                         ∈ 
                         N 
                       
                       , 
                       
                         k 
                         ∈ 
                         M 
                       
                     
                   
                 
               
               
                 
                   ( 
                   3 
                   ) 
                 
               
             
           
         
       
     
         [0069]    where KR k   R1  and K k   R2  represent the unit weight of penalty incurred by surface requirement changeover in the kth line, and R i  and R j  represent the surface requirement of the coil i and the surface requirement of the coil j, respectively. 
         [0070]    (4) abrasion penalty F ijk   i  of equipment incurred by equipment adjustment due to annealing temperature changeover of coils i and j in line k: 
         [0000]        F   ijk   t   =K   k   t   ·|t   i   −t   j   |, i,j∈N, k∈M    (4)
 
         [0071]    where K k   t  represents the unit weight of penalty incurred by annealing temperature changeover in the kth line, and t i  and t j  represent the annealing temperature of the coil i and the annealing temperature of the coil j, respectively. 
         [0072]    (5) abrasion penalty F ijk   wt  of equipment incurred by equipment adjustment due to weight changeover of the surface plates of coils i and j in line k 
         [0000]        F   ijk   wt   =K   k   wt   ·|wt   i   −wt   j   |, i,j∈N, k∈M    (5)
 
         [0073]    where K k   wt  represents the unit weight of penalty incurred by weight changeover of the surface plates in the kth line, and wt i  and wt j  represent the weight changeover of the surface plates of coil i and coil j, respectively. 
         [0074]    (6) abrasion penalty F ijk   G  of equipment incurred by equipment adjustment due to post-processing changeover of coils i and j in line k 
         [0000]    
       
         
           
             
               
                 
                   
                     F 
                     ijk 
                     G 
                   
                   = 
                   
                     { 
                     
                       
                         
                           
                             
                               K 
                               k 
                               G 
                             
                           
                           
                             
                               
                                 G 
                                 i 
                               
                               ≠ 
                               
                                 G 
                                 j 
                               
                             
                           
                         
                         
                           
                             0 
                           
                           
                             
                               
                                 G 
                                 i 
                               
                               = 
                               
                                 G 
                                 j 
                               
                             
                           
                         
                       
                       , 
                       
                           
                       
                        
                       i 
                       , 
                       
                         j 
                         ∈ 
                         N 
                       
                       , 
                       
                         k 
                         ∈ 
                         M 
                       
                     
                   
                 
               
               
                 
                   ( 
                   6 
                   ) 
                 
               
             
           
         
       
     
         [0075]    where K k   G  represents the unit weight of penalty incurred by post-processing changeover in the kth line, and G i  and G j  represent the post-processing code of the coil i and the post-processing code of the coil j respectively. 
         [0076]    (7) abrasion penalties and changeover penalties incurred by equipment adjustment due to different production requirements of coils i and j in each line, which is calculated in detail as follows: 
         [0077]    The abrasion penalties and changeover penalties F ij0  incurred by equipment adjustment due to different production requirements of coils i and j in the acid rolling line: 
         [0000]        F   ij0   =F   ij0   w   +F   ij0   h   +F   ij0   R    i, j∈N    (7)
 
         [0078]    The abrasion penalties and changeover penalties F ij1  incurred by equipment adjustment due to different production requirements of coils i and j in the continuous annealing line: 
         [0000]        F   ij1   =F   ij1   w   +F   ij1   h   +F   ij1   t    i, j∈N    (8)
 
         [0079]    The abrasion penalties and changeover penalties F ij2  incurred by equipment adjustment due to different production requirements of coils i and j in the hot galvanizing line: 
         [0000]        F   ij2   =F   ij2   w   +F   ij2   h   +F   ij2   t   +F   ij2   wt   +F   ij2   G    i, j∈N    (9)
 
         [0080]    The abrasion penalties and changeover penalties F ij3  incurred by equipment adjustment due to different production requirements of coils i and j on the hot aluminizing and galvanizing line: 
         [0000]        F   ij3   =F   ij3   w   +F   ij3   h   +F   ij3   t   +F   ij3   wt   +F   ij3   G    i, j∈N    (10)
 
         [0081]    Step 2-2: among the coil schedules corresponding to the production lines in the cold rolling supply chain composed by multiple production lines, an objective function of coordinated coil scheduling of multiple lines in the cold rolling area is built as minimizing the abrasion penalties and changeover penalties incurred by equipment adjustment due to different production requirements: 
         [0000]    
       
         
           
             
               
                 
                   Minimize 
                    
                   
                       
                   
                    
                   
                     
                       ∑ 
                       
                         k 
                         ∈ 
                         M 
                       
                     
                      
                     
                         
                     
                      
                     
                       
                         ∑ 
                         
                           j 
                           ∈ 
                           
                             N 
                             k 
                           
                         
                       
                        
                       
                           
                       
                        
                       
                         
                           ∑ 
                           
                             i 
                             ∈ 
                             
                               N 
                               k 
                             
                           
                         
                          
                         
                             
                         
                          
                         
                           
                             F 
                             ijk 
                           
                            
                           
                             x 
                             ijk 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   11 
                   ) 
                 
               
             
           
         
       
     
         [0082]    where the decision variable of the objective function is represented as x ijk , which indicates whether the coil i is processed immediately before the coil j, in the kth production line, if yes, the decision variable x ijk  is assigned as 1, and if not, it is assigned as 0; 
         [0083]    Step 2-3: the constraints for the objective function comprise the constraints for each of the multiple lines and the cohesion constraints among the production lines. 
         [0084]    The constraints of each of the multiple respective lines include: (1) for each of all the actual coils, there must be exactly one coil being processed immediately before it; (2) for each of all the actual coils, there must be exactly one coil being processed immediately after it; (3) the changeover in width, thickness, surface requirement, and annealing temperatures of the lines should not exceed the predetermined upper limits; (4) it should be guaranteed that the determination to the processing time of the coils in each line will not incur idle time between any two adjacent processed coils; (5) for each of the coils, it is only processed exactly once in a production line along its processing path; 
         [0085]    (1) the constraints of each of the multiple production lines include: 
         [0086]    1) for each of all the actual coils, there must be exactly one coil being processed immediately before it; 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         ∑ 
                         
                           i 
                           ∈ 
                           
                             
                               N 
                               k 
                             
                             ⋃ 
                             
                               { 
                               0 
                               } 
                             
                           
                         
                       
                        
                       
                           
                       
                        
                       
                         x 
                         ijk 
                       
                     
                     = 
                     1 
                   
                   , 
                   
                       
                   
                    
                   
                     ∀ 
                     
                       j 
                       ∈ 
                       
                         N 
                         k 
                       
                     
                   
                   , 
                   
                     k 
                     ∈ 
                     
                       M 
                       . 
                     
                   
                 
               
               
                 
                   ( 
                   12 
                   ) 
                 
               
             
           
         
       
     
         [0087]    where N k  represents the set of coils to be processed in line k. 
         [0088]    2) for each of all the actual coils, there must be exactly one coil being processed immediately after it; 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         ∑ 
                         
                           j 
                           ∈ 
                           
                             
                               N 
                               k 
                             
                             ⋃ 
                             
                               { 
                               
                                 N 
                                 + 
                                 1 
                               
                               } 
                             
                           
                         
                       
                        
                       
                           
                       
                        
                       
                         x 
                         ijk 
                       
                     
                     = 
                     1 
                   
                   , 
                   
                       
                   
                    
                   
                     ∀ 
                     
                       i 
                       ∈ 
                       
                         N 
                         k 
                       
                     
                   
                   , 
                   
                     k 
                     ∈ 
                     
                       M 
                       . 
                     
                   
                 
               
               
                 
                   ( 
                   13 
                   ) 
                 
               
             
           
         
       
     
         [0089]    3) constraints for sequencing dummy coils, constraint (14) means that, there must be exactly one actual coil being processed immediately after the head dummy coil, and constraint (15) means that, there must be exactly one actual coil being processed immediately before the last dummy coil, 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         ∑ 
                         
                           j 
                           ∈ 
                           
                             N 
                             k 
                           
                         
                       
                        
                       
                           
                       
                        
                       
                         x 
                         
                           0 
                            
                           
                               
                           
                            
                           jk 
                         
                       
                     
                     = 
                     1 
                   
                   , 
                   
                     ∀ 
                     
                       k 
                       ∈ 
                       
                         M 
                         . 
                       
                     
                   
                 
               
               
                 
                   ( 
                   14 
                   ) 
                 
               
             
             
               
                 
                   
                     
                       
                         ∑ 
                         
                           i 
                           ∈ 
                           
                             N 
                             k 
                           
                         
                       
                        
                       
                           
                       
                        
                       
                         x 
                         
                           i 
                           , 
                           
                               
                           
                            
                           
                             n 
                             + 
                             1 
                           
                           , 
                           k 
                         
                       
                     
                     = 
                     1 
                   
                   , 
                   
                     ∀ 
                     
                       k 
                       ∈ 
                       
                         M 
                         . 
                       
                     
                   
                 
               
               
                 
                   ( 
                   15 
                   ) 
                 
               
             
           
         
       
     
         [0090]    where the dummy coils are used to indicate the beginning and finishing of each line. 
         [0091]    4) the changeover in width, thickness, thicknesses of the surface plates and annealing temperatures of the lines should not exceed the predetermined upper limits: 
         [0000]        |w   i   −w   j   |×x   ijk   ≦w   0   , ∀i, j∈N, k∈M    (16)
 
         [0000]        |h   i   −h   j   |×x   ijk   ≦h   0   , ∀i, j∈N, k∈M    (17)
 
         [0000]        |t   i   −t   j   |×x   ijk   ≦t   0   , ∀i, j∈N, k∈M    (18)
 
         [0000]        |wt   i   −wt   j   |×x   ijk   ≦wt   0   , ∀i, j∈N, k∈M    (19)
 
         [0092]    where w 0  is the predetermined upper limit of the width changeover, h 0  is the predetermined upper limit of the thickness changeover, t 0  is the predetermined upper limit of the annealing temperature changeover, and wt 0  is the predetermined upper limit of the thickness changeover of the surface plates; 
         [0093]    5) it should be guaranteed that the determination to the processing time of the coils in each line will not incur idle time between any two adjacent processed coils; 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       C 
                       ik 
                     
                     = 
                     
                       
                         
                           x 
                           
                             0 
                              
                             
                                 
                             
                              
                             ik 
                           
                         
                          
                         
                           p 
                           ik 
                         
                       
                       + 
                       
                         
                           ∑ 
                           
                             j 
                             ∈ 
                             
                               N 
                               k 
                             
                           
                         
                          
                         
                             
                         
                          
                         
                           
                             x 
                             jik 
                           
                            
                           
                             ( 
                             
                               
                                 C 
                                 jk 
                               
                               + 
                               
                                 p 
                                 ik 
                               
                             
                             ) 
                           
                         
                       
                     
                   
                   , 
                   
                       
                   
                    
                   
                     ∀ 
                     
                       i 
                       ∈ 
                       
                         N 
                         k 
                       
                     
                   
                   , 
                   
                     k 
                     ∈ 
                     
                       M 
                       . 
                     
                   
                 
               
               
                 
                   ( 
                   20 
                   ) 
                 
               
             
           
         
       
     
         [0094]    where C ik  represents the completion time of the coil i in the kth line, and p ik  represents the processing time of coil i in line k. 
         [0095]    6) for each of the coils, it should be processed exactly once in a production line along its required processing path: 
         [0000]        d   jk ≧( d   ik +1)−(1− x   ijk )·( n+ 1), ∀ i∈N   k ∪{0},  j∈N   k   ∪{n+ 1},  k∈M,    (21)
 
         [0000]        d   ik   ≦n+ 1, ∀ i∈N   k ∪{0,  n+ 1},  k∈M.    (22)
 
         [0096]    where d ik  represents the processing position number of coil i in the coil schedule of line k. 
         [0097]    7) constraints (19)˜(23) are constraints for range of the variables. 
         [0000]        x   iik =0, ∀ i∈N   k ∪{0,  n+ 1},  k∈M.    (23)
 
         [0000]        x   i0k =0, ∀ i∈N   k   ∪{n+ 1},  k∈M.    (24)
 
         [0000]        x   n+1, i,k =0, ∀ i∈N   k ∪{0},  k∈M.    (25)
 
         [0000]        x   ijk ={0,1}, ∀ i∈N   k ∪{0},  j∈N   k   ∪{n+ 1},  k∈M.    (26)
 
         [0000]        C   ik ≧0, ∀ i∈N   k ∪{0,  n+ 1},  k∈M.    (27)
 
         [0098]    The cohesion constraints among the production lines imply the cohesion relationship between the production lines for guaranteeing that the coil cannot be processed in a downstream production line until it is finished in the upstream production line; 
         [0000]        C   ik   ≦C   ik′   −p   ik′   , ∀i∈N   k   , k∈M, k′∈L   k . 
         [0099]    where L k  represents the set of the downstream lines immediately following the line k. 
         [0100]    Step 3: an initial coordinated coil schedule of multiple production lines in the cold rolling area is obtained by using heuristic algorithm, which includes an initial coil schedule in acid rolling line, an initial coil schedule in a continuous annealing line, an initial coil schedule in a hot galvanizing line, and an initial coil schedule in a hot aluminizing and galvanizing line; 
         [0101]    Step 3-1: Determine the initial coil schedule in acid rolling line: an initial coil schedule in the acid rolling line is obtained by sequencing the coils to be processed in a descending order of the surface requirement code, and then sequencing the coils with the same surface requirement in a descending order of their widths, and then sequencing the coils with the same surface requirement and the same width in a descending order of their post-processing codes, and then sequencing the coils with the same surface requirement and the same width and the same post-processing code in an descending order of their annealing temperatures; 
         [0102]    Step 3-2: calculating the starting time and the finishing time of each coil in the initial coil schedule in the acid rolling line; 
         [0103]    Step 3-3: determine an initial coil schedule in the continuous annealing line: an initial annealing schedule of coils in the material inventory of continuous annealing line is obtained by sequencing the coils to be annealed in a descending order of their widths, and then sequencing the coils with the same width in a descending order of their annealing temperatures, and then sequencing the coils with the same width and the same annealing temperature in a descending order of their thicknesses, a coil schedule in continuous annealing line is obtained by selecting the coils to be annealed from the coil schedule in the acid rolling line, and sequentially inserting them after the initial annealing schedule of coils in the material inventory of continuous annealing line; 
         [0104]    Step 3-4: determine an initial coil schedule in the hot galvanizing line: an initial schedule of the coils in the material inventory of hot galvanizing line is obtained, by sequencing the coils in the material inventory of hot galvanizing line in a descending order of their post-processing codes; and then sequencing the coils with the same post-processing code in a descending order of their widths; and then sequencing the coils with the same width and the same post-processing code in a descending order of their thicknesses; and then sequencing the coils with the same post-processing code, and the same width and the same thickness in a descending order of their annealing temperatures; an initial coil schedule in the hot galvanizing line is obtained by selecting the coils to be processed in hot galvanizing line from the initial coil schedule in acid rolling line, sequentially inserting them after the initial schedule of the coils in the material inventory of hot galvanizing line; 
         [0105]    Step 3-5: determine an initial coil schedule in the hot aluminizing and galvanizing line: an initial schedule of coils in material inventory of the hot aluminizing and galvanizing line is obtained, by sequencing the coils in material inventory of the hot aluminizing and galvanizing line in a descending order of their post-processing codes; and then sequencing the coils with the same post-processing code in a descending order of their widths; and then sequencing the coils with the same post-processing code and the same width in a descending order of their thicknesses; and then sequencing the coils with the same post-processing code, the same width and the same thickness in a descending order of their annealing temperatures; an initial coil schedule in the hot aluminizing and galvanizing line is obtained, by selecting the coils to be processed in hot aluminizing and galvanizing line from the initial coil schedule in acid rolling line, sequentially inserting them after the initial schedule of coils in material inventory of the hot aluminizing and galvanizing line; 
         [0106]    Step 3-6: calculating the starting time of each coil in the initial coil schedule in the continuous annealing line, the starting time of each coil in the initial coil schedule in the hot galvanizing line, and the starting time of each coil in the initial coil schedule in the hot aluminizing and galvanizing line; 
         [0107]    Step 3-7: recording an initial coordinated coil schedule of multiple production lines in the cold rolling area, and calculating the objective function value corresponding to the initial coordinated coil schedule of multiple lines in the cold rolling area; 
         [0108]    An initial coordinated coil schedule of multiple lines in cold rolling area is obtained as follows: 
         [0109]    the coil schedule in acid rolling line: 8→9→6→7→4→1→2→3→5→15→16→13→14; 
         [0110]    the coil schedule in continuous annealing line: 11→10→12→4→1→2→3→5; 
         [0111]    the coil schedule in hot galvanizing line: 6→7 
         [0112]    the coil schedule in hot aluminizing and galvanizing line: 8→9 
         [0113]    Step 4: improving in real time the initial coordinated coil schedule of multiple lines in cold rolling area, as shown in  FIG. 3 ; 
         [0114]    Step 4 specifically comprises the following steps: 
         [0115]    Step 4-1: initializing the population of the coordinated coil schedule of multiple lines in cold rolling area, and setting the given maximum iteration number UN iter  and the given maximum iteration number without improvement UN iternoimpro , population of the coordinated coil schedule of multiple lines in the cold rolling area comprises the initial coordinated coil schedule of multiple lines in the cold rolling area, and other coil schedules in the population of the coordinated coil schedule of multiple lines in the cold rolling area are generated randomly; and during the improvement, an adaptive value of each coordinated coil schedule is calculated according to the way for calculating the objective function of each coordinated coil schedule, and the calculation is in accordance with the equation (11); 
         [0116]    Setting the crossover factor CR in a range of [0.1, 0.5], the scaling factor F is 0.7, the population size PS is 80, the given maximum iteration number is 300, and the given maximum iteration number without improvement is 50. 
         [0117]    Step 4-2: judging whether the given maximum iteration number is reached, and if it is reached, then taking the schedule with the best adaptive value in the current population as the optimal schedule, that is, the improved coordinated coil schedule of multiple lines in the cold rolling area, and if not, go to step 4-3; 
         [0118]    Step 4-3: judging whether the given maximum iteration number without improvement is reached, and if it is reached, then taking the schedule with the best adaptive value in the current population as the optimal schedule, that is, the improved coordinated coil schedule of multiple lines in the cold rolling area, and if not, go to step 4-4; 
         [0119]    Step 4-4: if the individual with the best adaptive value in the population is not updated, then the maximum iteration number without improvement N iternoimpro =N iternoimpro +1. 
         [0120]    Step 4-5: performing a mutation operation to the population, specifically as follows: 
         [0000]    
       
         
           
             
               
                 
                   
                     V 
                     i 
                     t 
                   
                   = 
                   
                     
                       X 
                       a 
                       
                         t 
                         - 
                         1 
                       
                     
                     ⊕ 
                     
                       F 
                       ⊗ 
                       
                         ( 
                         
                           
                             X 
                             b 
                             
                               t 
                               - 
                               1 
                             
                           
                           - 
                           
                             X 
                             c 
                             
                               t 
                               - 
                               1 
                             
                           
                         
                         ) 
                       
                     
                   
                 
               
               
                 
                   ( 
                   29 
                   ) 
                 
               
             
             
               
                 
                   
                     Δ 
                     i 
                     t 
                   
                   = 
                   
                     
                       
                         F 
                         ⊗ 
                         
                           ( 
                           
                             
                               X 
                               b 
                               
                                 t 
                                 - 
                                 1 
                               
                             
                             - 
                             
                               X 
                               c 
                               
                                 t 
                                 - 
                                 1 
                               
                             
                           
                           ) 
                         
                       
                       ⇔ 
                       
                         δ 
                         ij 
                         t 
                       
                     
                     = 
                     
                       { 
                       
                         
                           
                             
                               
                                 x 
                                 bj 
                                 
                                   t 
                                   - 
                                   1 
                                 
                               
                               - 
                               
                                 x 
                                 cj 
                                 
                                   t 
                                   - 
                                   1 
                                 
                               
                             
                           
                           
                             
                               
                                 if 
                                  
                                 
                                     
                                 
                                  
                                 
                                   rand 
                                   ( 
                                   ) 
                                 
                               
                               &lt; 
                               F 
                             
                           
                         
                         
                           
                             0 
                           
                           
                             otherwise 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   30 
                   ) 
                 
               
             
             
               
                 
                   
                     Δ 
                     i 
                     t 
                   
                   = 
                   
                     ( 
                     
                       
                         δ 
                         
                           i 
                           , 
                           0 
                         
                         t 
                       
                       , 
                       
                         δ 
                         
                           i 
                           , 
                           1 
                         
                         t 
                       
                       , 
                       … 
                        
                       
                           
                       
                       , 
                       
                         δ 
                         
                           i 
                           , 
                           n 
                         
                         t 
                       
                     
                     ) 
                   
                 
               
               
                 
                   ( 
                   31 
                   ) 
                 
               
             
             
               
                 
                   
                     V 
                     i 
                     t 
                   
                   = 
                   
                     
                       
                         
                           X 
                           a 
                           
                             t 
                             - 
                             1 
                           
                         
                         ⊕ 
                         
                           Δ 
                           i 
                           t 
                         
                       
                       ⇔ 
                       
                         v 
                         ij 
                         t 
                       
                     
                     = 
                     
                       
                         
                           x 
                           aj 
                           
                             t 
                             - 
                             1 
                           
                         
                         ⊕ 
                         
                           δ 
                           ij 
                           t 
                         
                       
                       = 
                       
                         mod 
                          
                         
                           ( 
                           
                             
                               ( 
                               
                                 
                                   x 
                                   aj 
                                   
                                     t 
                                     - 
                                     1 
                                   
                                 
                                 + 
                                 
                                   δ 
                                   ij 
                                   t 
                                 
                                 + 
                                 n 
                               
                               ) 
                             
                             , 
                             n 
                           
                           ) 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   32 
                   ) 
                 
               
             
           
         
       
     
         [0121]    where V i   t  is the ith variant chromosome in the tth generation, and v ij   t  is the jth individual in the V i   t , and X a   t−1 , X b   t−1  and X c   t−1  are three different coordinated coil schedules of multiple lines in cold rolling area in the population of the t−1th generation, x ij   t−1  means the jth individual in the X i   t−1 , F means the scaling factor, Δ i   t  means the ith intermediate variant chromosome in the mutation process of the tth generation, δ ij   t  means the jth individual in the intermediate variant chromosome Δ i   t  in the mutation process of the tth generation, and mod means modulus. 
         [0122]    Step 4-6: performing a crossover operation to the population, specifically as follows: 
         [0123]    generating a trial chromosome U i   t , i=1, 2, . . . , PS, and u ij   t  is an individual in the chromosome U i   t ; 
         [0000]    
       
         
           
             
               
                 
                   
                     u 
                     ij 
                     t 
                   
                   = 
                   
                     { 
                     
                       
                         
                           
                             v 
                             ij 
                             t 
                           
                         
                         
                           
                             
                               
                                 
                                   if 
                                    
                                   
                                       
                                   
                                    
                                   
                                     rand 
                                     ( 
                                     ) 
                                   
                                 
                                 &lt; 
                                 
                                   CR 
                                    
                                   
                                       
                                   
                                    
                                   or 
                                    
                                   
                                       
                                   
                                    
                                   j 
                                 
                               
                               = 
                               
                                 D 
                                 j 
                               
                             
                              
                             
                                 
                             
                           
                         
                       
                       
                         
                           
                             x 
                             ij 
                             
                               t 
                               - 
                               1 
                             
                           
                         
                         
                           otherwise 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   33 
                   ) 
                 
               
             
           
         
       
     
         [0124]    where PS represents the number of chromosomes of the initial population, CR is in [0,1] and in order to ensure the variety of the population, rand( ) generates a group of uniformly distributed random numbers between 0 and 1, D j  is a number randomly selected from the set{1, 2, . . . , n}, so as to guarantee that at least one dimension of the U i   t  is different from the ith coordinated coil schedule for the multiple lines in the cold rolling area in the t−1th generation X i   t−1 ; 
         [0125]    CR is an self-updating number as the iteration number increases, and is updated as follows: 
         [0000]    
       
         
           
             
               
                 
                   CR 
                   = 
                   
                     
                       CR 
                       min 
                     
                     + 
                     
                       t 
                       * 
                       
                         
                           ( 
                           
                             
                               CR 
                               max 
                             
                             - 
                             
                               CR 
                               min 
                             
                           
                           ) 
                         
                         T 
                       
                     
                   
                 
               
               
                 
                   ( 
                   34 
                   ) 
                 
               
             
           
         
       
     
         [0126]    where CR min  is the minimum crossover probability, CR max  is the maximum crossover probability, t means the current iteration number, and T is the maximum iteration number. 
         [0127]    Step 4-7: performing a selecting operation on the son generation in the population; 
         [0000]    
       
         
           
             
               
                 
                   
                     X 
                     i 
                     t 
                   
                   = 
                   
                     { 
                     
                       
                         
                           
                             U 
                             i 
                             t 
                           
                         
                         
                           
                             
                               
                                 if 
                                  
                                 
                                     
                                 
                                  
                                 
                                   f 
                                   ( 
                                   
                                     U 
                                     i 
                                     t 
                                   
                                   ) 
                                 
                               
                               ≤ 
                               
                                 f 
                                  
                                 
                                   ( 
                                   
                                     X 
                                     i 
                                     
                                       t 
                                       - 
                                       1 
                                     
                                   
                                   ) 
                                 
                               
                             
                              
                             
                                 
                             
                           
                         
                       
                       
                         
                           
                             X 
                             i 
                             
                               t 
                               - 
                               1 
                             
                           
                         
                         
                           otherwise 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   35 
                   ) 
                 
               
             
           
         
       
     
         [0128]    where f( ) is the calculating equation for adaptive value based on the objective function (11), and the equation (31) generates a new object chromosome finally; 
         [0129]    Step 4-8: calculating an adaptive value for each individual in the population, that is, calculating the objective function value of coordinated coil schedule of multiple lines in cold rolling area corresponding to the current population; 
         [0130]    Step 4-9: the current iteration number is increment N iter =N iter +1 and the population is updated, and then return to step 4-2. 
         [0131]    the improved coil schedule in acid rolling line: 8→9→6→7→4→2→3→1→5→15→16→14→13; 
         [0132]    the improved coil schedule in continuous annealing line: 11→10→12→2→3→4→1→5 
         [0133]    the improved coil schedule in hot galvanizing line: 6→7 
         [0134]    the improved coil schedule in hot aluminizing and galvanizing line: 8→9 
         [0135]    Step 5: sending the improved coordinated schedule to the ACS (automatic control system) of the respective lines in cold rolling area, and the coordinated coil scheduling of multiple lines in cold rolling area is completed.