CN113805545B - Flexible flow shop combined scheduling rule generation method considering batch processing - Google Patents
Flexible flow shop combined scheduling rule generation method considering batch processing Download PDFInfo
- Publication number
- CN113805545B CN113805545B CN202111074759.0A CN202111074759A CN113805545B CN 113805545 B CN113805545 B CN 113805545B CN 202111074759 A CN202111074759 A CN 202111074759A CN 113805545 B CN113805545 B CN 113805545B
- Authority
- CN
- China
- Prior art keywords
- batch
- workpiece
- processing
- individuals
- scheduling
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000012545 processing Methods 0.000 title claims abstract description 118
- 238000000034 method Methods 0.000 title claims abstract description 95
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 97
- 230000002068 genetic effect Effects 0.000 claims abstract description 53
- 230000008569 process Effects 0.000 claims abstract description 53
- 238000004519 manufacturing process Methods 0.000 claims abstract description 34
- 230000014509 gene expression Effects 0.000 claims abstract description 24
- 238000005457 optimization Methods 0.000 claims abstract description 22
- 230000003044 adaptive effect Effects 0.000 claims abstract description 11
- 230000008030 elimination Effects 0.000 claims abstract description 11
- 238000003379 elimination reaction Methods 0.000 claims abstract description 11
- 230000007547 defect Effects 0.000 claims abstract description 10
- 108090000623 proteins and genes Proteins 0.000 claims description 34
- 238000012163 sequencing technique Methods 0.000 claims description 19
- 238000003780 insertion Methods 0.000 claims description 12
- 230000037431 insertion Effects 0.000 claims description 12
- 210000000349 chromosome Anatomy 0.000 claims description 9
- 239000012634 fragment Substances 0.000 claims description 9
- 230000006872 improvement Effects 0.000 claims description 9
- 238000003754 machining Methods 0.000 claims description 7
- 230000035772 mutation Effects 0.000 claims description 7
- 238000010845 search algorithm Methods 0.000 claims description 4
- 238000011425 standardization method Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 9
- 238000010923 batch production Methods 0.000 description 8
- 238000009827 uniform distribution Methods 0.000 description 8
- 238000000354 decomposition reaction Methods 0.000 description 6
- 230000000052 comparative effect Effects 0.000 description 3
- 230000003252 repetitive effect Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000002028 premature Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 210000001726 chromosome structure Anatomy 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000010353 genetic engineering Methods 0.000 description 1
- 238000002789 length control Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/418—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM]
- G05B19/41885—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM] characterised by modeling, simulation of the manufacturing system
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/32—Operator till task planning
- G05B2219/32339—Object oriented modeling, design, analysis, implementation, simulation language
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- General Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- General Factory Administration (AREA)
Abstract
The invention provides a method for generating a combined scheduling rule of a flexible flow shop, which considers batch processing. The method comprises the following steps: decomposing the FFSP-BPM problem related to the invention, constructing corresponding scheduling rules, and combining the scheduling rules to generate a combined scheduling rule for solving the original problem; the method comprises the steps that a multi-objective optimization model is built by considering the FFSP-BPM problem of incompatible workpiece family batch operation under the environment of workpiece dynamic arrival and uncertain working hours; and generating a combined scheduling rule by adopting an improved gene expression programming algorithm to solve the FFSP-BPM problem, and providing improved strategies of repeated individual elimination, variable neighborhood search and adaptive genetic operators aiming at the defects that the algorithm is easy to generate repeated individuals and fall into local optimum. The method aims at the flexible flow group batch workshop scheduling process considering incompatible workpiece family parallel batch processing, adopts a method of generating a combined scheduling rule, improves the production efficiency, and can effectively guide the scheduling process.
Description
Technical Field
The invention relates to the field of dynamic workshop scheduling, in particular to a method for generating a flexible flow workshop combined scheduling rule considering batch processing.
Background
With the improvement of social and economic levels, the market demand tends to diversify and individually customize products, and the modern manufacturing mode is gradually changed from the original mass production to a flexible multi-variety small-batch production mode. The flexible flow production can flexibly adapt to market changes, realizes the mixed line operation of various products, and is a mainstream mode of discrete manufacturing enterprises for coping with the small-batch production of various products. The flexible flow production considering batch processing is a typical application scene, is widely applied to a plurality of fields such as semiconductor manufacturing and the like, and is mainly characterized in that a batch production mode capable of simultaneously processing a plurality of limited operations is combined on the basis of the original flexible flow production.
The batch production mode can achieve the aims of avoiding frequent production preparation time, simplifying material handling and improving production efficiency. However, the combination of this method and a flexible flow production scenario increases the complexity of production operation management, and the scheduling process needs to select processing equipment, determine the processing sequence of workpieces on the equipment, and consider how to batch workpieces. Meanwhile, the frequent change of the market environment enables the workpiece to have the characteristic of dynamic arrival, and random dynamic events such as working hour fluctuation and the like can occur due to different equipment states in the production process. Due to the dynamic factors, an enterprise cannot acquire scheduled global task information in advance and can only make a decision according to real-time production data, so that full batch operation is difficult to realize in a batch processing stage. Therefore, it is necessary to provide an efficient Scheduling method based on real-time production data for solving the Problem of Flexible Flow-sheet Scheduling with Batch Process Machines (FFSP-BPM) in a Flexible Flow shop with Batch Process Scheduling features in a dynamic environment, so as to adapt to the current production mode of multiple varieties of small batches, and respond to dynamic factors in the production Process in time while improving the production efficiency.
In the current research methods related to real-time scheduling, a scheduling rule method (DR) is widely applied. The method has lower time complexity, can be flexibly suitable for various problems, and in addition, the scheduling decision based on the scheduling rule is given before implementation, so the method can respond to dynamic factors in time and is very suitable for solving the FFSP-BPM problem with higher complexity and real-time scheduling requirement. However, the method for solving the FFSP-BPM problem has two defects: 1) The FFSP-BPM problem is a multi-decision combined optimization problem of multiple stages, multiple equipment types and flexible paths, and the traditional scheduling rule method is single in consideration dimension and cannot be directly solved. 2) The scheduling rule is a method based on local information, the performance of the scheduling rule depends on various factors such as a scheduling target and system configuration, a single scheduling rule cannot adapt to all dynamic environments, and dynamic selection of the scheduling rule needs to be realized by combining different scheduling scenes and performance indexes.
Disclosure of Invention
Aiming at the multi-target FFSP-BPM problem of incompatible workpiece family batch operation under the environment of dynamic arrival of workpieces and uncertain working hours, the method for generating the combined scheduling rule of the flexible flow shop considering batch processing is provided. The invention firstly introduces a scheduling rule method which is widely applied in the field of dynamic scheduling, secondly constructs a combined scheduling rule for solving the original problem by utilizing the thought of problem decomposition, and finally provides an FFSP-BPM combined scheduling rule generation method based on a Gene Expression Programming (GEP) algorithm. The technical means adopted by the invention are as follows:
a method for generating a combined scheduling rule of a flexible flow shop considering batch processing comprises the following steps:
s1, decomposing the FFSP-BPM problem, analyzing problem factors aiming at the decomposed subproblems to construct corresponding scheduling rules, and then combining the scheduling rules corresponding to the subproblems to generate a combined scheduling rule for solving the original problem;
s2, constructing a multi-objective optimization model aiming at FFSP-BPM problems of incompatible workpiece family batch operation under the environment of dynamic arrival and uncertain working hours of the workpiece;
and S3, generating a combined scheduling rule for solving the multi-objective optimization model by adopting a gene expression programming algorithm, and providing improved strategies of repeated individual elimination, variable neighborhood search and adaptive genetic operators aiming at the defects that repeated individuals are easy to appear in the algorithm and the algorithm is trapped in local optimization.
Further, the FFSP-BPM problem is specifically: when different types of workpieces dynamically arrive at a workshop, the production tasks can be completed through K processing stages in sequence, wherein the first stage is a batch processing stage, l is more than or equal to 1 and less than or equal to K, the stage comprises a plurality of batch processors with the load limiting quantity of B, the batch processing tasks of a limited number of workpieces can be completed, the other stages are single-piece processing, and M is used for processing the workpieces (k) A parallel machine component, wherein M (k) More than or equal to 1, workpieces can only be processed by one equipment of the process in each processing stage, the processing time of different types of workpieces in the same processing stage is different, the processing time fluctuation caused by the new and old states of the equipment is considered, the processing time of the same type of workpieces on different equipment in the same stage is also different, batch processing equipment can only process the same type of workpieces at the same time, the operation time depends on the longest processing time in the operation contained in each batch, namely, a parallel batch processing operation mode which is incompatible with a workpiece family is adopted.
Further, based on the FFSP-BPM problem, the following assumptions are made: (1) the buffer area between adjacent stages is an infinite buffer area; (2) batch size is related only to the number of workpieces in a batch; and (3) the batch phase operation has the minimum capacity requirement.
Further, in the S1, decomposing the FFSP-BPM problem specifically means decomposing the original problem into a flexible flow shop scheduling problem and a batch processing scheduling problem;
aiming at the scheduling problem of the flexible flow shop, the job priority is judged according to the job parameters so as to realize job sequencing and the equipment priority is judged according to the equipment parameters so as to realize equipment selection;
aiming at the sub-problem of batch processing scheduling, workpieces are firstly distinguished according to workpiece families, the workpieces in each product type are subjected to job sequencing according to priorities, then feasible batches are generated according to batch requirements, and finally, for all the feasible batches of all the workpiece types, a batch forming process is guided through batch parameters and the batch priorities are determined, so that the batch sequencing of incompatible workpiece families is realized.
Further, the step S2 includes the steps of:
s2-1: and defining parameters and decision variables of the multi-objective optimization model with the maximum completion time and the average delay deviation.
(1) Parameter definition: the workpiece type index is I, wherein I =1,2 … …, I, the workpiece index is j, wherein j =1,2 … …, N, the machining stage index is K, wherein K =1,2 … …, K, the lot set is B, wherein B =1,2 … …, B, the tool index is M, wherein M =1,2 … …, M, d j Is the delivery date of the work j, r i Time of entry of workpiece j into production System, s jk For the time at which the workpiece j begins to be machined at stage k, p jkm For the machining time, p, of the workpiece j on the apparatus m of stage k jk Actual machining time, O, for workpiece j at stage k jk For the workpiece j, the process is completed in stage k, i.e. the k-th process, C, of the workpiece ik Time TO complete processing of workpiece j at stage k, TO jk Is O ijk Type of operation of l ij For correspondence of workpiece j and product type i, TM m Of type m, CM m For the processing capacity of the plant m, MB m Is the minimum batch size of the batch processing equipment m, B m Number of batches, P, for a batch processing apparatus m mb Is the processing time of batch b on batch processing equipment m;
(2) Decision variables:
s2-2: and defining a multi-objective optimization model objective function and constraint conditions of the maximum completion time and the average delay deviation.
The maximum time-out C max Refers to the finish time of the last workpiece; the average delay deviation (AE/L) is used for measuring the difference degree between the completion time and the delivery date of the workpiece and tends to be delivered on time;
the optimization objective is represented as:
f 1 ＝min(max j C jK ) (4)
wherein, f 1 、f 2 Respectively representing a function aimed at minimizing the maximum completion time and the mean delay deviation, C jK Representing the time of completion of the workpiece j in the last processing stage K, d j Indicating the delivery date of the workpiece j;
according to the characteristics of the actual production problem, the problem model should satisfy the following constraints:
constraints (6) indicate that a tool can only process one workpiece at a time, constraints (7) indicate that workpieces can only be processed once at each stage, constraints (8) indicate that parts cannot be processed until they reach the production system, constraints (9) ensure that discrete processing tools can only process one workpiece at a time, constraints (10) indicate that processing of the next process can only begin if the previous stage is complete, constraints (11) define the completion time of a workpiece at stage k, constraints (12) indicate the actual processing time of a workpiece at stage k, constraints (13) indicate that the batch processing time on the batch processing tool is equal to the maximum value of the workpiece processing time in the batch, constraints (14) ensure that a workpiece is allocated to only one batch, constraints (15) specify that the batch size cannot be less than the minimum batch requirement and cannot exceed the batch processing tool capacity, constraints (16) indicate that workpieces allocated to the same batch have the same start time at the batch processing stage, constraints (17) indicate that workpiece j can only belong to one product type i, and constraints (18) indicate that workpieces of one batch type can only come from one product type.
Further, step S3 includes the steps of:
s3-1, constructing an endpoint set and a function set based on the classical scheduling rules, adopting the classical scheduling rules as the endpoint set, and selecting operators { +, -,/, -/() 2 The GEP algorithm represents different scheduling rules through the linear combination of the endpoint set and the function set;
s3-2, based on the codes of the combined dispatching rules, respectively representing the dispatching rules of the equipment selection process, the job sorting process and the batch sorting process according to a three-section real number coding mode;
s3-3, decoding the chromosome to generate a group of dispatching rule function expressions comprising equipment selection rules, job sorting rules and batch sorting rules;
s3-4, defining the effective length of the gene based on the genetic operation of the effective length, designing a genetic operator based on the effective length, and carrying out the genetic operation on elements in the effective length;
s3-5: aiming at the defects that the algorithm is early converged, easily falls into local optimum and is easy to generate repeated individuals, an algorithm improvement strategy is introduced: 1) A repeated individual eliminating algorithm for eliminating the weight of the population after genetic operation to ensure the diversity of the population; 2) Two neighborhood structures for improving the local search capability of the algorithm; 3) An adaptive genetic operator for adjusting the genetic probability.
Further, in step S3-4, the effective length-based genetic operation includes four modes, i.e., a crossover operator, a mutation operator, an RIS insertion operator, and a non-dominated sorting and selecting operator.
Further, the algorithm improvement strategy of step S3-5 includes:
(1) Repeated individual elimination algorithm
Step1: traversing the merged population, and searching for repeated individuals.
Step2: determining an individual to be operated: if the two individuals which are mutually repeated are respectively positioned in the parent population and the offspring population, taking the parent individual corresponding to the repeated individual in the offspring population as an individual to be operated, and if the two individuals which are mutually repeated are both positioned in the offspring population, randomly selecting one individual corresponding to the parent individual as the individual to be operated;
step3: repeated individual removal: repeatedly carrying out variation and RIS string insertion operation on the individuals to be operated until the generated individuals are not repeated with other individuals in the combined population;
after repeated individuals are removed, all the individuals in the combined population after the duplication removal are subjected to non-dominated sorting, and the optimal N individuals are selected as a new generation population P according to the quality degree of the individuals in the population i+1 ；
(2) Variable neighborhood search algorithm
Performing non-domination sorting on the offspring population generated after the repeated individuals are removed, and performing variable neighborhood search on the individuals in the non-domination stage so as to improve the local search capability of the algorithm;
the neighborhood structure definition includes: randomly scrambling and rearranging the neighborhood, namely scrambling and rearranging all end point set elements in the effective length of the current gene segment; forward insertion neighborhood, namely randomly selecting two positions of the end point set of the two elements in the effective length part of the current gene fragment, inserting the element corresponding to the latter position into the former position, and sequentially moving the elements of the other end point sets backwards;
(3) Adaptive genetic operator
Carrying out dimensionless on the performance targets by adopting a max-min standardization method, taking the sum of the standardized targets as the fitness value of the current individual, and endowing the poorer individual with higher genetic operation probability by adopting a self-adaptive genetic operator, wherein the concrete formula is as follows:
in the formula, P max And P min The maximum value and the minimum value of the genetic operation probability are respectively, f is the fitness value of the current individual to be operated, and for the crossover operator, f is the larger value of the fitness of the two individuals to be crossed. f. of max And f avg The maximum value and the average value of the fitness of all individuals in the population to be operated are respectively.
The technical conception of the invention is as follows: the invention aims at the scheduling problem of a flexible flow line batch workshop which is incompatible with parallel batch processing of a workpiece family under the environment of dynamic arrival of workpieces and uncertain working hours, takes the minimized maximum completion time and average delay deviation as optimization targets, introduces a scheduling rule method widely applied to the field of real-time scheduling, constructs a combined scheduling rule for solving the original problem by using the thought of problem decomposition, and further provides a FFSP-BPM combined scheduling rule generation method based on gene expression programming.
The effective effects of the invention are mainly shown in that: in consideration of the defect that the traditional scheduling rule method has single dimension and cannot directly solve the FFSP-BPM combined decision optimization problem related by the invention, the scheduling rule generation algorithm based on gene expression programming and capable of effectively mining functional relation is provided. Aiming at the defects that the traditional gene expression programming algorithm has premature convergence, is easy to fall into local optimum, is easy to generate repeated individuals and the like, the improvement method is provided from the two aspects of problem characteristics and algorithm performance, and the algorithm performance is further improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
FIG. 1 is a diagram of FFSP-BPM scheduling issues in view of incompatible workpiece families in accordance with an embodiment of the present invention;
FIG. 2 is a diagram of a process for generating a problem resolution and combined scheduling rule according to an embodiment of the present invention;
FIG. 3 is a K-expression diagram in an embodiment of the present invention;
FIG. 4 is a diagram illustrating a decoding process according to an embodiment of the present invention;
FIG. 5 is a diagram of a single point forward cross, a two point sequential cross, and a genetic cross in accordance with an embodiment of the present invention;
FIG. 6 is a diagram of single point variations in an embodiment of the present invention;
FIG. 7 is a RIS plug-string diagram in accordance with an embodiment of the present invention;
FIG. 8 is a diagram of repetitive individuals in an embodiment of the present invention;
FIG. 9 is a flow chart of a variable neighborhood search in an embodiment of the present invention;
FIG. 10 is a comparison algorithm coverage boxplot in an embodiment of the present invention;
FIG. 11 is a non-dominated kanji diagram of the s3b4f2n30_0.8 \ u 2 algorithm obtained by improving the GEP algorithm in an embodiment of the present invention;
FIG. 12 is a Gantt chart of an example of the s3b4f2n30_0.8 _2algorithm obtained from the SP-ATC-RMSTB classic schedule rule in an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Aiming at the multi-target FFSP-BPM problem of incompatible workpiece family batch operation under the environment of workpiece dynamic arrival and uncertain working hours, the invention firstly introduces a scheduling rule method which is widely applied in the field of dynamic scheduling, secondly constructs a combined scheduling rule for solving the original problem by using the thought of problem decomposition, and finally provides a FFSP-BPM combined scheduling rule generation method based on gene expression programming. The method specifically comprises the following steps:
s1: constructing a combined scheduling rule based on a problem decomposition idea: the method comprises the steps of decomposing an FFSP-BPM original problem, respectively analyzing problem factors aiming at the decomposed sub-problems to construct corresponding scheduling rule algorithms, and then combining scheduling rules corresponding to the sub-problems to generate a combined scheduling rule for solving the original problem.
S2: and (3) constructing a multi-objective optimization model aiming at the FFSP-BPM problem of incompatible workpiece group batch operation under the environment of dynamic arrival and uncertain working hour of the workpiece.
And S3, generating a combined scheduling rule by adopting an improved gene expression programming algorithm to solve the FFSP-BPM problem, and providing improved strategies of repeated individual elimination, variable neighborhood search and self-adaptive genetic operators aiming at the defects that repeated individuals are easy to occur in the algorithm and the algorithm is trapped into local optimum.
Before step S1, the FFSP-BPM scheduling problem according to the present invention is problem-defined. Specifically, the FFSP-BPM problem is described in FIG. 1, when different types of workpieces dynamically arrive at a workshop, the production task can be completed after K processing stages in sequence, and the FFSP-BPM problem is describedThe middle l (l is more than or equal to 1 and less than or equal to K) stage is a batch processing stage, the stage comprises a plurality of batch processors with the load limiting quantity of B, the batch processing tasks of a limited number of workpieces can be completed, other stages are single-piece processing, and M is used for processing (k) (M (k) Not less than 1), the workpiece can only be processed by one device of the process at each processing stage, the processing time of different types of workpieces at the same processing stage is different, the processing time fluctuation caused by the new and old states of the device is considered, and the processing time of the same type of workpieces at different devices at the same stage is also different. Batch processing facilities can only process the same type of workpiece at the same time, with the job time depending on the longest processing time in the job contained in each batch, i.e., in a parallel batch job with incompatible workpiece families.
The following assumptions are proposed based on the above: (1) the buffer area between adjacent stages is an infinite buffer area; (2) batch size is related only to the number of workpieces in a batch; (3) the batch phase job has a minimum capacity requirement.
The problem decomposition-based combined scheduling construction process in step S1 is shown in fig. 2, and specifically includes the following steps:
s1-1: FFSP-BPM scheduling problem decomposition. The flexible flow shop batch scheduling problem can be decomposed into a flexible flow shop scheduling problem and a batch processing scheduling problem. For the former, firstly, the processing sequence of the workpiece at each stage needs to be decided, namely, the operation priority is judged through the operation parameters so as to realize operation sequencing (the workpiece is selected by the equipment); secondly, the work equipment of the workpiece needs to be determined, namely equipment selection (workpiece selection equipment) is realized by judging equipment priority through equipment parameters. For the sub-problem of batch processing scheduling, the decision process is divided into three steps, firstly, workpieces are distinguished according to workpiece families, job sorting is carried out on the workpieces in each product type according to priorities, then feasible batches are generated according to the batch requirements, and finally, for all the feasible batches of all the workpiece types, the batch forming process is guided through batch parameters and the batch priorities are determined, so that the batch sorting (equipment selection batches) of incompatible workpiece families is realized.
S1-2: the main problem factor of the combined scheduling rule structure influencing the operation sorting subproblem is that workpieces need to be determined according to the operation sorting rule based on the related attributes of the workpieces such as arrival time, delivery date, process number, processing time and the like. The main problem factor affecting the equipment selection subproblem is equipment, and the priority of the equipment needs to be determined according to the equipment selection rule based on information on equipment-related attributes such as work efficiency, earliest available time, number of waiting workpieces, and the like. The problem factor affecting the batch generation subproblem is a workpiece, and firstly, an alternative workpiece range is determined according to a Time Window strategy (TW) based on the attribute of the workpiece arrival Time; and secondly, generating all feasible batches meeting the minimum capacity requirement according to the attribute of the workpiece type. The problem factors affecting the lot ordering sub-problem are the workpieces in all feasible lots, and the candidate lot priorities need to be determined according to the lot ordering rules based on the information about the workpiece family related attributes such as lot, operation time, mean delivery time, lot weight, etc.
And combining the 3 types of heuristic rules analyzed above, and forming a combined dispatching rule algorithm for solving the FFSP-BPM by matching with a time window strategy for solving batch generation subproblems. Combined scheduling rules (DR) HFSP-BPM ) Expressed by equation (1), JS represents job sort rules, MS represents equipment selection rules, and BS represents lot sort rules.
DR HFSP-BPM ＝{MS,JS,BS} (1)
And designing a decision variable from the solution of the sub-problems, providing related parameter definition to construct a mathematical model, and performing joint decision on the three sub-problems to generate an FFSP-BPM real-time scheduling scheme in a dynamic environment. The step S2 includes the steps of:
s2-1: and defining parameters and decision variables of the multi-objective optimization model with the maximum completion time and the average delay deviation.
(1) Decision variables
The final objective of the general flow shop scheduling subproblem is to determine the work-on time of the workpiece at each stage by the job sequencing rule, so the corresponding decision variables are shown in formula (2).
The final purpose of the scheduling subproblem of the parallel machine is to determine the operation equipment of the workpiece at each stage through equipment selection rules, so that the corresponding decision variables are shown in formula (3).
The final purpose of the sub-problem of batch scheduling is to determine the processing batch and the operating equipment of the workpiece in the batch processing stage according to the batch sorting rule, and because the invention considers the batch processing process of incompatible workpiece families, if the operating equipment is determined according to the equipment selection rule and then the workpiece family classification, batch generation and sorting are carried out on the to-be-processed operation of the equipment, the utilization rate of the batch processing equipment is greatly reduced, so the invention does not adopt the equipment selection rule to decide the equipment to be operated by the workpiece in the batch processing stage, but synthesizes the feasible batches of all the equipment for sorting, and the corresponding decision variable is shown in the formula (4).
In summary, the decision variables of the FFSP-BPM scheduling problem are the combined scheduling rules DR including the equipment selection rule MS, the job scheduling rules JS and the batch scheduling rules BS HFSP-BPM = JS, MS, BS. The scheduling objective of the FFSP-BPM scheduling problem researched by the invention is to determine DR HFSP-BPM Therefore, the processing sequence of the workpieces in the batch processing stage, the processing batches and the processing equipment and the optimal processing sequence of the workpieces in the discrete processing stage are determined, and the start-up and completion time of each process of the workpieces is further determined so as to achieve the performance goal of the production system.
(2) Definition of related parameters
Table 1 shows the relevant parameters of the model:
TABLE 1 parameter definitions
S2-2: and defining a multi-objective optimization model objective function and constraint conditions of the maximum completion time and the average delay deviation.
The present invention considers two objective functions: maximum time-out C max Average delay bias (AE/L). The maximum completion time refers to the completion time of the last workpiece; the average delay deviation measures the difference degree between the completion time and the delivery date of the workpiece, and tends to be delivered on time. The optimization objectives of the present invention can be expressed as:
f 1 ＝min(max j C jK ) (5)
wherein f is 1 、f 2 Respectively representing a function aimed at minimizing the maximum completion time and the mean delay deviation, C jK Representing the time of completion of the workpiece j in the last processing stage K, d j Indicating the delivery date of workpiece j.
According to the characteristics of the actual production problem, the problem model should satisfy the following constraints:
constraint (7) indicates that a piece of equipment can only process one workpiece at a time, and constraint (8) indicates that a workpiece can only be processed once at each stage. Constraint (9) indicates that the part cannot be processed until it reaches the production system. The constraints (10) ensure that the discrete processing apparatus can only process one workpiece at a time. Constraint (11) means that only the previous stage is completed and the next process can be started. The constraint (12) defines the completion time of the workpiece at stage k. The constraint (13) represents the actual machining time of the workpiece at the k-th stage. The constraint (14) indicates that the run time of a batch on the batch processing tool is equal to the maximum value of the run times of the workpieces in the batch. Constraints (15) ensure that a workpiece is assigned to only one lot. Constraints (16) specify that the batch size cannot be less than the minimum batch requirement and cannot exceed the batch processing equipment capacity. Constraint (17) indicates that workpieces assigned to the same batch have the same start time in the batch phase. Constraint (18) indicates that a workpiece j can only belong to one product type i, and constraint (19) indicates that a batch of workpieces can only come from one product type.
Aiming at a multi-target model taking minimized maximum completion time and average delay deviation as optimization targets, the characteristic of incompatible workpiece family parallel batch processing is considered, and then an improved multi-target model solving method based on a gene expression programming algorithm in the step S4 is provided, wherein the step S3 comprises the following steps:
and S3-1, constructing an endpoint set and a function set based on the classical scheduling rule. Taking the classic scheduling rules as a set of endpoints, select the operators { +, -,/, √, () 2 The GEP algorithm represents different scheduling rules by linear combination of endpoint set and function set as elements of function set.
And S3-2, coding based on the combined scheduling rule. And according to a three-section real number coding mode, the dispatching rules of the equipment selection process, the job sorting process and the batch sorting process are respectively expressed.
And S3-3, decoding the chromosome to generate a group of dispatching rule function expressions comprising equipment selection rules, job sequencing rules and batch sequencing rules.
S3-4, effective length-based genetic manipulation. Defining effective length of gene, designing genetic operator based on effective length, and carrying out genetic operation on elements in effective length.
S3-5: aiming at the defects that the algorithm is prematurely converged, easily falls into local optimum and is easy to generate repeated individuals, an algorithm improvement strategy is provided: 1) A repeated individual elimination algorithm is provided, and the duplication elimination is carried out on the population after the genetic operation so as to ensure the diversity of the population; 2) Two neighborhood structures are provided to improve the local searching capability of the algorithm; 3) And introducing an adaptive genetic operator to adjust the genetic probability.
The details of the structure relating to the function set and the endpoint set for step S3-1 are as follows:
(1) Function set
The GEP algorithm represents different scheduling rules through linear combination of an endpoint set and a function set, the search space of the GEP algorithm grows exponentially with the increase of the number of set elements, and in order to improve the algorithm efficiency while reducing the search space, the set needs to be designed appropriately. By integrating the results of the related research, the present invention selects the operators { +, -/() 2 As elements of the function set, where the division is a protective division, i.e. returns 1 when divided by 0, otherwise returns the normal quotient, and if the square root parameter is negative, the result is
Non-linear and linear combinations of parameters can be achieved by these operators.
(2) Endpoint set
Different from the traditional research that the use environment and the system parameters are directly used as the endpoint set, the invention uses the widely applied classic scheduling rules in the literature as the elements in the endpoint set and outputs the functional relationship among the classic scheduling rules. In combination with the problem features of the present invention, the following table 2 is used for the device selection process, job sorting process, and batch process, respectively.
TABLE 2 classical scheduling rules
And the step S3-2 adopts a three-section real number coding mode and is respectively used for representing the dispatching rules of the equipment selection process, the job sorting process and the batch sorting process. Based on the fixed-length coding principle of the GEP algorithm, each gene segment is divided into a head part and a tail part, the head region of each gene segment is randomly selected from a function set or an end point set of the function set, and the tail region is only selected from the end point set of the function set.
The gene length is the sum of the head length h and the tail length t, and the two relations satisfy t = h x (n-1) +1, wherein n represents the maximum value of the number of each operator in the function set. The gene generated in this manner is referred to as a K expression. The numbers of corresponding meshes of the operators in the function set are respectively {2,2,2,2,1,1}, so that n =2, the lengths of the heads of the three gene segments are respectively {4,5,4}, the lengths of the tails of the three gene segments are respectively {5,6,5}, the lengths of the genes are respectively {9,11,9}, a chromosome is constructed as shown in fig. 3, the numbers in the chromosome correspond to the elements in the end point set, and Q represents the square.
The chromosome decoding process in step S3-3 is illustrated by FIG. 4.
Step1: the K Expression is converted into a K Expression Tree (ET) using a Depth First Search (DFS) algorithm. And taking the first element of the gene segment as a root node of the tree structure, sequentially traversing from left to right according to the number of the K expression until all the nodes of the last layer are the symbols in the terminator set, and ending the iteration.
Step2: and traversing the K-expression tree by adopting a middle-order traversal algorithm to obtain a function expression of the scheduling rule. The left sub-tree is traversed firstly, then the root node is visited, finally the right sub-tree is traversed, when the left sub-tree and the right sub-tree are traversed, the left sub-tree is still traversed firstly, the root node is visited, the right sub-tree is traversed, and after the traversal is finished, a scheduling rule function expression based on relevant elements in the production process can be generated.
And S3-4, calculating the fitness, namely decoding the chromosome to generate a group of scheduling rule function expressions comprising an equipment selection rule, a job sequencing rule and a batch sequencing rule, applying the mined scheduling rule to an experimental scene, constructing a scheduling rule evaluation simulation model, and calculating the individual fitness value according to the selected 2 individual performance indexes. In order to avoid the problem that the characteristic parameter dimension in the endpoint set is not uniform and the function operation cannot be carried out, the normalization processing is carried out on the relevant parameter values of the workshop according to a formula (20).
In the formula, theta max And theta min The maximum value and the minimum value of the workshop related parameter value theta are respectively, and theta is the value of the parameter theta of the workpiece or equipment. R θ The normalized parameter values are used as the priority of the workpiece determined according to the parameter theta.
The effective length of the defined gene in the step S3-4 is mainly defined by four modes of a crossover operator, a mutation operator, an RIS insertion operator and a non-dominant sequencing and selection operator. The GEP algorithm is easier to realize genetic operation, but a phenomenon that partial codes do not work in a final practical expression occurs, so that invalid elements appear at the tail of a gene sequence. The genetic operation on invalid elements cannot ensure effective iteration of a population, so the effective length of a gene is defined, and a genetic operator based on the effective length is designed to perform the genetic operation on the elements in the effective length.
Definition 1: the effective length L of the gene. And the chromosome is converted into the node number of a tree structure formed after the K expression tree is formed. Such as the three gene segments shown in FIG. 3, which have effective lengths of 6,9,6.
(1) Crossover operator
Three crossover operators are set for two individuals to be crossed: single point crossover, two point crossover, and gene crossover. Single-point forward crossing tends to cross the functional combinatorial relationship between two individual terminal elements, two-point crossing tends to cross a certain substructure in two individual tree structures, and genetic crossing is used to simulate different combinations of three types of rules. And the chromosomes are fully transformed by combining three types of crossing modes to increase the population diversity and enhance the global searching capability of the algorithm.
The single-point crossing and two-point crossing operations are performed for each gene segment of two individuals to be crossed. Firstly, determining the minimum value C of the effective lengths of the current gene segments of two individuals to be crossed min ＝min{L i ,L j And maximum value C max ＝max{L i ,L j }. For single point crossing, randomly in [1,C min ]To select a node for forward single point crossing as shown in fig. 5 (a). For two-point crossing, randomly in [0,C min ]One node is selected as a starting node Cross 1 In [ Cross ] 1 ,C max ]Randomly selects a node as an end node Cross 2 The elements between the two are exchanged as shown in FIG. 5 (b), wherein the yellow segment is the effective length part of the gene segment. The gene crossover operation was performed for a certain gene fragment among the three gene fragments as shown in FIG. 5 (c).
(2) Mutation operator
Setting a single point mutation operator for each gene segment of an individual to be mutated is shown in fig. 6. The smaller value of the effective length and the head length of the gene fragment i is taken as the maximum range of the variation position, namely, the curve M at [0,min i ,h}]And randomly selecting a variation position, if the position is positioned at the head, randomly selecting other elements which are different from the existing elements of the variation position from the function set and the union set of the endpoint set, and if the variation position is positioned at the tail, randomly selecting from the endpoint set.
(3) RIS insert string operator
The RIS (Root Insertion Sequence) Insertion operator is a specific genetic operator of GEP, and can effectively destroy the original structure of individual genotype and phenotype. The RIS insertion string operation performed on each gene fragment of an individual is shown in fig. 7. In [1,h]Randomly selecting one element in the range as the position R of the operator in the function set 1 In [ R ] 1 ,h+t]Randomly selecting a position R within range 2 Will [ R ] 1 ,R 2 ]The fragment is inserted into the head of the fragment, the original element is correspondingly moved backwards, and a gene segment with the same length as the inserted fragment is deleted from the tail of the original head so as to ensure the legality of the chromosome structure.
(4) Non-dominant sorting and selection operator
Introducing a Non-dominant sorting method of an NSGA-II (Non-dominant sorting genetic algorithm-II) belt elite strategy, and evaluating according to fitnessAs a result, all the non-inferior solutions are found from the current population, the crowding distances thereof are calculated, the non-inferior solutions found are removed from the population by assigning the grades thereto in order until the population size is 0. The child selection is performed according to the non-inferior solution grade, and the non-inferior solution of the same grade is selected according to the congestion distance of the non-inferior solution. To avoid premature convergence of the algorithm, each iteration is based on a parent population P of size N i Generating an offspring population Q with the same size as the parent population i A 1 is to P i And Q i And merging the two generations of populations, and generating the next generation of population on the basis.
The improvement strategy for the algorithm described in S3-5 is detailed below:
(1) Repeated individual elimination algorithm
Because the GEP algorithm uses fixed-length codes to represent the characteristics of tree structures with different sizes, although the genetic operation based on the effective length controls the difference of individuals before and after the genetic operation, the situation that two individuals with different genotypes but the same phenotype type and the same generated scheduling rule function expression are generated in a combined population cannot be avoided. The pair of individuals having the same bitsubstance of each gene within the effective length is defined as a repetitive individual. The repeated individuals generated by genetic operation have two sources, including that an offspring individual and a parent individual are mutually repeated individuals, and the inside of an offspring population is mutually repeated individuals. For example, as shown in FIG. 8, the P1 individuals and the P of the parent population are generated after the crossing operation is performed on the P1 and P2 individuals i The individuals are mutually repetitive individuals; to P i C after individual RIS inserting string operation i Individual and pair P j C produced by individual after single point mutation j The individuals are mutually repeated individuals.
In order to avoid that the algorithm is prematurely converged to be in local optimum due to the occurrence of repeated individuals, before non-dominated sorting and filial generation selection, the repeated individuals occurring in the GEP iterative evolution process are removed. The method comprises the following specific steps:
step1: traversing the merged population, and searching for repeated individuals.
Step2: and determining an individual to be operated. And if the two individuals which are mutually repeated are respectively positioned in the parent population and the offspring population, taking the parent individuals corresponding to the repeated individuals in the offspring population as the individuals to be operated. And if the two individuals which are mutually repeated are both positioned in the offspring population, randomly selecting one individual corresponding to the parent individual as the individual to be operated.
Step3: and (5) repeatedly removing individuals. And repeating the variation and RIS string insertion operation on the individuals to be operated until the generated individuals are not repeated with other individuals in the combined population.
After repeated individuals are removed, all the individuals in the combined population after the duplication removal are subjected to non-dominated sorting, and the optimal N individuals are selected as a new generation population p according to the quality degree of the individuals in the population i+1 。
(2) Variable neighborhood search algorithm
The global search capability of the algorithm is enhanced to a great extent through the genetic operation based on the effective length and the selection operation based on repeated individual elimination, and the diversity of the population is ensured. By changing the Neighborhood structure near the current local optimal solution, the Variable Neighborhood Search algorithm (VNS) can continuously Search for a better local optimal solution, and can effectively improve the local Search capability of the GEP algorithm.
Two neighborhood structures are designed according to the structural characteristics of the problem, the offspring population generated after repeated individual elimination is subjected to non-domination sorting, and the individuals in the non-domination level are subjected to variable neighborhood searching, so that the local searching capability of the algorithm is improved. The neighborhood structure definition includes: randomly scrambling and rearranging the neighborhood, namely scrambling and rearranging all end point set elements in the effective length of the current gene segment; and (3) forward insertion neighborhood, namely randomly selecting the positions of two elements belonging to the endpoint set in the effective length part of the current gene segment, inserting the element corresponding to the latter position into the front position, and sequentially moving the elements of the other endpoint sets backwards. The variable neighborhood search process is illustrated in fig. 9.
(3) Adaptive genetic operator
When an evolutionary Algorithm such as Genetic Algorithm (GA) and GEP is used for solving the problem, the higher the Genetic operation probability is, the faster the speed of generating a new individual is, the higher the fitness individual structure is, and the search process tends to random search; if the genetic operation probability is too small, the search process becomes slow even before the search is stopped, the population diversity cannot be maintained, and the local optimum is easily trapped. In order to avoid the disadvantages caused by using the fixed genetic probability, the invention introduces an adaptive strategy, dynamically controls cross operation, mutation operation and RIS string insertion probability, ensures that the population can be quickly converged, simultaneously prevents the diversity of the population from being damaged, improves the adaptability of the search to the space change, and specifically sets a method as shown in a formula (21).
In the formula, P max And P min The maximum value and the minimum value of the genetic operation probability are respectively, f is the fitness value of the current individual to be operated, and for the crossover operator, f is the larger value of the fitness of the two individuals to be crossed. f. of max And f avg The maximum value and the average value of the fitness of all individuals in the population to be operated are respectively. Aiming at the multi-objective optimization problem, the performance targets are subjected to non-dimensionalization by adopting a max-min standardization method, and the sum of the standardized targets is taken as the fitness value of the current individual. And the self-adaptive genetic operator is adopted to endow a poorer individual with a higher genetic operation probability, accelerate the updating and evolution of the individual, and endow a superior individual with a lower genetic operation probability so as to protect the superior individual.
Example 1
An improved gene expression programming algorithm is applied. Design example related parameters, the number of production stages is subject to uniform distribution U (3,7), and the batch processing stage B k Is located at the first
Stage and batch processor capacity B obeys uniform distribution U (2,6), workpiece type number obeys uniform distribution U (2,6), workpiece scale obeys uniform distribution U (300, 1000), and equipment utilization rate obeys uniform distribution U (0.5,0.9); the delivery date tension factor obeys uniform distribution U (4,8), and the invention respectively sets 3 distribution ranges according to the example parametersThe level value.
And setting relevant parameters for each calculation example, and respectively setting operation time distribution for discrete steps of the batch processing procedure. Batch process processing time p b Subject to a uniform distribution U (100, 200), discrete process processing times p d Subject to uniform distribution U (1,30). The workpiece weight setting method employs a 4. 18 orthogonal examples of different scales were generated based on the IBM SPSS Statistics 19 platform.
The algorithm parameters of the present invention relating to the improved GEP algorithm and the comparison algorithm (GEP algorithm, plant parameter GEP algorithm and classical scheduling rules algorithm) are shown in Table 3.
TABLE 3 GEP Algorithm and comparison Algorithm parameter settings
The improved GEP algorithm, the workshop parameter GEP algorithm and the classical scheduling rule algorithm are adopted to solve all the example problems, 5 independent simulation experiments are carried out, the average values of the trans-generation distance and the distributivity index delta in 5 running times are counted and shown in the table 4, wherein the value shown in bold is the optimal value obtained by calculation of the 4 algorithms. As can be seen from the table, for most of the test problems, the Inverse Genetic Distance (IGD) and the distributivity index Δ of the improved GEP algorithm are superior to those of the GEP algorithm, the plant parameter GEP and the classical scheduling rule, and as can be seen from the Total result in the last row in the table, the improved GEP algorithm is superior to the other three comparative algorithms in terms of convergence, diversity and distributivity of non-inferior solutions for all the test problems.
Table 4 comparative algorithm inverse distance IGD mean results
Fig. 10 is a box plot of the relative coverage of the modified GEP algorithm versus the three other comparative algorithms. Each box plot has different examples on the abscissa and relative coverage on the ordinate, representing the relative coverage C (a, B) between row associated algorithm a and column associated algorithm B. The box plot as the top right corner represents the ratio of the non-dominant solution generated by the modified GEP algorithm to cover the solution in the classical scheduling rule. It can be seen from the figure that the solutions generated by the classical scheduling rules are almost all dominated by the non-inferior solutions generated by the improved GEP algorithm, which proves that the improved GEP algorithm is improved in both the maximum completion time and the average delay deviation, and moreover, the non-inferior solutions generated by the improved GEP algorithm can dominate the vast majority of non-column generated by the plant parameter GEP algorithm, while C (improved GEP, GEP) > C (GEP, improved GEP) as a whole. Therefore, the effectiveness of the algorithm improvement strategy provided by the invention, the effectiveness of the strategy that the terminal element adopts the classical scheduling rule and the effectiveness of the strategy that the intelligent algorithm is adopted to generate the scheduling rule can be proved.
Examples adopt S n1 B n2 F n3 N n4 U _ C denotes a flexible line shop with n1 stages and a batch stage capacity of n2, the processing task comprising n4 workpieces of n3 product types, a shop utilization of U and a lead factor of C. In order to show the process of guiding the output of the scheduling scheme by the scheduling rule, a small-scale example of s3b4f2n30_0.8 \ u 2 is additionally designed for explanation, and this example shows that the application scene is a flexible flow shop with 3 stages and a batch processing stage capacity of 4, the processing tasks include 30 workpieces of 2 product types, the shop utilization rate is 0.8, and the workpiece delivery date tension factor is 2. Fig. 11 is a gantt chart corresponding to a non-dominant solution (i.e., a combined dispatching rule including a device selection rule, a job sorting rule, and a lot sorting rule) generated by the GEP algorithm according to this example, and the combined dispatching rule formed after decoding the non-dominant solution is shown in table 5. FIG. 12Gantt chart generated for the s3b4f2n30_0.8 _2algorithm using the SP-ATC-RMSTB classical scheduling rules. Wherein the blue-based workpiece and the orange-based workpiece respectively represent a product type 1 and a product type 2. It can be seen that under the classical scheduling rules approach, the batch processing phase completes 10 batches, each containing an average of 3 workpieces, for a total completion time of 2150. In the scheduling scheme formed by improving the scheduling rule generated by the GEP algorithm, 14 batch jobs are completed in the batch processing stage, each batch contains 2.5 workpieces on average, and the total completion time is shorter than that of the scheduling scheme generated by the classical scheduling rule. It can be shown that the production line can be more balanced by reasonable batch decision. The effectiveness of the improved GEP algorithm proposed by the present invention is also demonstrated.
TABLE 5A non-dominated solution of the s3b4f2n30_0.8 \ u 2 algorithm obtained by refining the GEP algorithm
In summary, aiming at the batch scheduling problem of the flexible flow shop, the performance and the production efficiency of the combined scheduling rule including the equipment selection rule, the operation sequencing rule and the batch sequencing rule generated by the improved GEP algorithm are greatly improved. Compared with a workshop parameter GEP algorithm, the scheduling rule formed by combining functions of the classical scheduling rule is proved to have better performance and can effectively guide the scheduling process compared with the scheduling rule based on the workshop parameter.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.
Claims (6)
1. A method for generating a combined scheduling rule of a flexible flow shop considering batch processing is characterized by comprising the following steps:
s1, decomposing the FFSP-BPM problem, constructing corresponding scheduling rules aiming at the decomposed subproblem analysis problem factors, and combining the scheduling rules corresponding to the subproblems to generate a combined scheduling rule for solving the original problem;
s2, constructing a multi-objective optimization model aiming at FFSP-BPM problems of incompatible workpiece family batch operation under the environment of dynamic arrival and uncertain working hours of the workpiece;
s3, generating a combined scheduling rule for solving the multi-objective optimization model by adopting a gene expression programming algorithm, and providing improved strategies of repeated individual elimination, variable neighborhood search and adaptive genetic operators aiming at the defects that repeated individuals are easy to appear in the algorithm and the algorithm is trapped in local optimization;
in the S1, the FFSP-BPM problem is decomposed into a flexible flow shop scheduling problem and a batch processing scheduling problem;
aiming at the scheduling problem of the flexible flow shop, the job priority is judged according to the job parameters so as to realize job sequencing and the equipment priority is judged according to the equipment parameters so as to realize equipment selection;
aiming at the sub-problem of batch processing scheduling, workpieces are firstly distinguished according to workpiece families, the workpieces in each product type are subjected to job sequencing according to priorities, then feasible batches are generated according to batch requirements, and finally, for all the feasible batches of all the workpiece types, a batch forming process is guided through batch parameters and the batch priorities are determined, so that the batch sequencing of incompatible workpiece families is realized;
the step S2 includes the steps of:
s2-1: defining parameters and decision variables of a multi-objective optimization model with the maximum completion time and the average delay deviation;
(1) Parameter definition: the workpiece type index is I, wherein I =1,2.. Said., I, the workpiece index is j, wherein j =1,2.. Said., N, the machining stage index is K, wherein K =1,2.. Said., K, the batch set is b, and wherein b =1,2.. Said., I, the workpiece index is j, and the workpiece index is j.B, device index is M, where M =1,2 j Is the delivery date of the workpiece j, r i Time of entry of workpiece j into production System, s jk Time for starting processing of workpiece j at stage k, p jkm For the machining time, p, of the workpiece j on the apparatus m of stage k jk For the actual machining time of workpiece j in stage k, O jk For a workpiece j, the process is completed in stage k, i.e. the k-th process, C, of the workpiece ik Time TO complete processing of workpiece j at stage k, TO jk Is O ijk Type of operation of l ij For correspondence of workpiece j and product type i, TM m Of type m, CM m For the processing capacity of the plant m, MB m Is the minimum batch size of the batch processing apparatus m, B m Number of batches, P, for a batch processing apparatus m mb Is the processing time of batch b on batch processing equipment m;
(2) Decision variables:
s2-2: defining a multi-objective optimization model objective function and constraint conditions of maximum completion time and average delay deviation;
said maximum completion time C max Refers to the finish time of the last workpiece; the average delay deviation AE/L is used for measuring the difference degree between the completion time and the delivery date of the workpiece and tends to be delivered on time;
the optimization model objective function is represented as:
f 1 ＝min(max j C jK ) (4)
wherein f is 1 、f 2 Respectively representing a function aimed at minimizing the maximum completion time and the mean delay deviation, C jK Representing the time of completion of the workpiece j in the last processing stage K, d j Indicating a delivery date of the workpiece j;
according to the characteristics of the actual production problem, the problem model should satisfy the following constraints:
constraint (6) indicates that a tool can only process one workpiece at a time, constraint (7) indicates that workpieces can only be processed once at each stage, constraint (8) indicates that parts cannot be processed before they reach the production system, constraint (9) ensures that discrete processing tools can only process one workpiece at a time, constraint (10) indicates that processing of the next process can only begin with completion of the previous stage, constraint (11) defines the completion time of a workpiece at stage k, constraint (12) indicates the actual processing time of a workpiece at stage k, constraint (13) indicates that the batch processing time on the batch processing tool is equal to the maximum of the workpiece processing times in the batch, constraint (14) ensures that a workpiece is assigned to only one batch, constraint (15) specifies that the batch size cannot be less than the minimum batch requirement and cannot exceed the batch processing tool capacity, constraint (16) indicates that workpieces assigned to the same batch have the same start time at the batch processing stage, constraint (17) indicates that workpiece j can only belong to one product type i, and constraint (18) indicates that workpieces of a batch can only come from one product type.
2. The method for generating the batch-considered flexible flow shop combined scheduling rule according to claim 1, wherein the FFSP-BPM problem is specifically: when different types of workpieces dynamically arrive at a workshop, the production tasks can be completed through K processing stages in sequence, wherein the first stage is a batch processing stage, l is more than or equal to 1 and less than or equal to K, the stage comprises a plurality of batch processors with the load limiting quantity of B, the batch processing tasks of a limited number of workpieces can be completed, the other stages are single-piece processing, and M is used for processing the workpieces (k) A parallel machine component, wherein M (k) More than or equal to 1, the workpiece can only be processed by one equipment of the procedure of the processing stage in each processing stage, the processing time of different types of workpieces in the same processing stage is different, the processing time fluctuation caused by the new and old states of the equipment is considered, the processing time of the same type of workpieces on different equipment in the same stage is also different, the batch processing equipment can only process the same type of workpieces at the same time, the operation time depends on the longest processing time in the operation contained in each batch, namely, a parallel batch processing operation mode which is incompatible with a workpiece family is adopted.
3. The batch-considered flexible flow shop combined scheduling rule generating method according to claim 2, wherein based on the FFSP-BPM problem, the following assumptions are made: (1) the buffer area between adjacent stages is an infinite buffer area; (2) batch size is related only to the number of workpieces in a batch; and (3) the batch phase operation has the minimum capacity requirement.
4. The batch-considered flexible flow shop assembly scheduling rule generating method according to claim 1, wherein the step3 comprises the steps of:
s3-1, constructing an endpoint set and a function set based on the classical scheduling rules, adopting the classical scheduling rules as the endpoint set, and selecting operators { +, -, -/() 2 As elements of a function setThe GEP algorithm represents different scheduling rules through the linear combination of an endpoint set and a function set;
s3-2, based on the coding of the combined scheduling rule, respectively representing the scheduling rules of an equipment selection process, an operation sequencing process and a batch sequencing process according to a three-section real number coding mode;
s3-3, decoding the chromosome to generate a group of dispatching rule function expressions comprising equipment selection rules, job sequencing rules and batch sequencing rules;
s3-4, defining the effective length of the gene based on the genetic operation of the effective length, designing a genetic operator based on the effective length, and carrying out the genetic operation on elements in the effective length;
s3-5: aiming at the defects that the algorithm is early converged, easily falls into local optimum and is easy to generate repeated individuals, an algorithm improvement strategy is introduced: 1) A repeated individual eliminating algorithm for eliminating the weight of the population after genetic operation to ensure the diversity of the population; 2) Two neighborhood structures for improving the local search capability of the algorithm; 3) An adaptive genetic operator for adjusting the genetic probability.
5. The batch-considered flexible flow shop combination scheduling rule generating method according to claim 4, wherein in step S3-4, the genetic operation based on the effective length includes four modes of crossover operator, mutation operator, RIS interpolation operator, non-dominated sorting and selection operator.
6. The batch-considered flexible flow shop assembly scheduling rule generating method according to claim 4, wherein the algorithm improvement strategy of the step S3-5 comprises:
(1) Repeated individual elimination algorithm
Step1: traversing the merged population, and searching for repeated individuals;
step2: determining an individual to be operated: if the two individuals which are mutually repeated are respectively located in the parent population and the offspring population, the parent individuals corresponding to the repeated individuals in the offspring population are used as the individuals to be operated, and if the two individuals which are mutually repeated are both located in the offspring population, one individual corresponding to the parent individual is randomly selected as the individual to be operated;
step3: repeated individual removal: repeatedly carrying out variation and RIS string insertion operation on the individuals to be operated until the generated individuals are not repeated with other individuals in the combined population;
after repeated individuals are removed, all the individuals in the combined population after the duplication removal are subjected to non-dominated sorting, and the optimal N individuals are selected as a new generation population P according to the quality degree of the individuals in the population i+1 ；
(2) Variable neighborhood search algorithm
Performing non-domination sorting on the offspring population generated after the repeated individuals are removed, and performing variable neighborhood search on the individuals in the non-domination stage so as to improve the local search capability of the algorithm;
the neighborhood structure definition includes: randomly scrambling and rearranging the neighborhood, namely scrambling and rearranging all end point set elements in the effective length of the current gene segment; forward insertion neighborhood, namely randomly selecting the positions of two elements belonging to an endpoint set in the effective length part of the current gene fragment, inserting the element corresponding to the latter position into the front position, and sequentially moving the elements of the other endpoint sets backwards;
(3) Adaptive genetic operator
Carrying out dimensionless treatment on the performance targets by adopting a max-min standardization method, taking the sum of the standardized targets as the fitness value of the current individual, and endowing a poorer individual with a larger genetic operation probability by adopting an adaptive genetic operator, wherein the specific formula is as follows:
in the formula, P max And P min Respectively representing the maximum value and the minimum value of the genetic operation probability, wherein f is the fitness value of the current individual to be operated, and for the crossover operator, f is the larger value of the fitness of the two individuals to be crossed; f. of max And f avg The maximum value and the average value of the fitness of all individuals in the population to be operated are respectively.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111074759.0A CN113805545B (en) | 2021-09-14 | 2021-09-14 | Flexible flow shop combined scheduling rule generation method considering batch processing |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111074759.0A CN113805545B (en) | 2021-09-14 | 2021-09-14 | Flexible flow shop combined scheduling rule generation method considering batch processing |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113805545A CN113805545A (en) | 2021-12-17 |
| CN113805545B true CN113805545B (en) | 2023-03-21 |
Family
ID=78895201
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111074759.0A Active CN113805545B (en) | 2021-09-14 | 2021-09-14 | Flexible flow shop combined scheduling rule generation method considering batch processing |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113805545B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114415614B (en) * | 2022-01-16 | 2023-09-01 | 西北工业大学 | Service class considered multiprocessor workpiece scheduling method |
| CN114662765A (en) * | 2022-03-25 | 2022-06-24 | 武汉理工大学 | Batch processing scheduling method for discrete manufacturing of irrelevant workpieces |
| CN115081755B (en) * | 2022-08-19 | 2022-12-09 | 合肥工业大学 | Production and maintenance cooperative scheduling method and system based on variable neighborhood search algorithm |
| CN116090788B (en) * | 2023-02-27 | 2023-12-22 | 湘南学院 | Batch scheduling planning method for flexible assembly job shop |
| CN116795054B (en) * | 2023-06-19 | 2024-03-19 | 上海交通大学 | Intermediate product scheduling method in discrete manufacturing mode |
| CN117215275B (en) * | 2023-11-08 | 2024-02-13 | 北京理工大学 | Large-scale dynamic double-effect scheduling method for flexible workshop based on genetic programming |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008047383A2 (en) * | 2006-07-28 | 2008-04-24 | Persistent Systems Private Limited | System and method for network association inference, validation and pruning based on integrated constraints from diverse data |
| CN106250583A (en) * | 2016-07-15 | 2016-12-21 | 西安电子科技大学 | Dynamic job shop scheduling rule optimization methods based on double population gene expression programmings |
| CN106611288A (en) * | 2016-10-12 | 2017-05-03 | 四川用联信息技术有限公司 | Improved solving algorithm for flexible flow shop scheduling problem |
| CN109270904A (en) * | 2018-10-22 | 2019-01-25 | 中车青岛四方机车车辆股份有限公司 | A kind of flexible job shop batch dynamic dispatching optimization method |
| CN109886589A (en) * | 2019-02-28 | 2019-06-14 | 长安大学 | A method of low-carbon Job-Shop is solved based on whale optimization algorithm is improved |
| WO2021036658A1 (en) * | 2019-08-29 | 2021-03-04 | 华中科技大学 | Multi-objective optimization method and system for master production schedule of casting parallel workshops |
-
2021
- 2021-09-14 CN CN202111074759.0A patent/CN113805545B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008047383A2 (en) * | 2006-07-28 | 2008-04-24 | Persistent Systems Private Limited | System and method for network association inference, validation and pruning based on integrated constraints from diverse data |
| CN106250583A (en) * | 2016-07-15 | 2016-12-21 | 西安电子科技大学 | Dynamic job shop scheduling rule optimization methods based on double population gene expression programmings |
| CN106611288A (en) * | 2016-10-12 | 2017-05-03 | 四川用联信息技术有限公司 | Improved solving algorithm for flexible flow shop scheduling problem |
| CN109270904A (en) * | 2018-10-22 | 2019-01-25 | 中车青岛四方机车车辆股份有限公司 | A kind of flexible job shop batch dynamic dispatching optimization method |
| CN109886589A (en) * | 2019-02-28 | 2019-06-14 | 长安大学 | A method of low-carbon Job-Shop is solved based on whale optimization algorithm is improved |
| WO2021036658A1 (en) * | 2019-08-29 | 2021-03-04 | 华中科技大学 | Multi-objective optimization method and system for master production schedule of casting parallel workshops |
Non-Patent Citations (1)
| Title |
|---|
| 刘巍巍 ; 马雪丽 ; 刘晓冰 ; .面向柔性作业车间调度问题的改进变邻域搜索算法.2015,(第04期),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113805545A (en) | 2021-12-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113805545B (en) | Flexible flow shop combined scheduling rule generation method considering batch processing | |
| US20210373888A1 (en) | Multi-objective optimization method and system for master production plan of casting parallel workshops | |
| US7668788B2 (en) | Resource assignment optimization using direct encoding and genetic algorithms | |
| Zolfaghari et al. | Comparative study of simulated annealing, genetic algorithms and tabu search for solving binary and comprehensive machine-grouping problems | |
| CN108460463B (en) | High-end equipment assembly line production scheduling method based on improved genetic algorithm | |
| CN112561225B (en) | Flexible job shop scheduling method based on marker post co-evolution algorithm | |
| CN112580922A (en) | Flexible job shop scheduling method based on multilevel neighborhood structure and hybrid genetic algorithm | |
| Zeng et al. | Robotic disassembly line balancing and sequencing problem considering energy-saving and high-profit for waste household appliances | |
| Yılmaz et al. | Lot streaming in hybrid flowshop scheduling problem by considering equal and consistent sublots under machine capability and limited waiting time constraint | |
| Nikabadi et al. | A hybrid algorithm for unrelated parallel machines scheduling | |
| CN115471133A (en) | Workshop comprehensive scheduling system based on order management and rolling optimization scheduling | |
| CN112381343A (en) | Flexible job shop scheduling method based on genetic-backbone particle swarm hybrid algorithm | |
| CN112990515A (en) | Workshop resource scheduling method based on heuristic optimization algorithm | |
| Teekeng et al. | A combination of shuffled frog leaping and fuzzy logic for flexible job-shop scheduling problems | |
| CN113689066A (en) | Internet of things workshop scheduling method based on NSGA-II algorithm | |
| Ławrynowicz | Genetic algorithms for solving scheduling problems in manufacturing systems | |
| Nidhiry et al. | Scheduling optimization of a flexible manufacturing system using a modified NSGA-II algorithm | |
| CN113139720B (en) | Casting shop production scheduling optimization method with learning effect | |
| Sadaghiani et al. | A Pareto archive floating search procedure for solving multi-objective flexible job shop scheduling problem | |
| CN111144569A (en) | Yield improvement applicable model optimization method based on genetic algorithm | |
| CN116307214A (en) | Automatic cotton distribution method based on NSGA-II algorithm and related device | |
| CN114237166A (en) | Method for solving multi-rotating-speed energy-saving scheduling problem based on improved SPEA2 algorithm | |
| CN114240091A (en) | Flexible job shop scheduling method based on self-adaptive layering strategy | |
| Laoraksakiat et al. | Bi-objective hybrid flow shop scheduling with family setup times using hybrid genetic and migrating birds optimization algorithms | |
| Zhang et al. | Multi-stage hybrid evolutionary algorithm for multiobjective distributed fuzzy flow-shop scheduling problem |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |