CN114936739A - Multi-target flexible job shop scheduling method based on improved cross entropy algorithm - Google Patents

Multi-target flexible job shop scheduling method based on improved cross entropy algorithm Download PDF

Info

Publication number
CN114936739A
CN114936739A CN202210366703.0A CN202210366703A CN114936739A CN 114936739 A CN114936739 A CN 114936739A CN 202210366703 A CN202210366703 A CN 202210366703A CN 114936739 A CN114936739 A CN 114936739A
Authority
CN
China
Prior art keywords
machine
individuals
population
target
search
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.)
Pending
Application number
CN202210366703.0A
Other languages
Chinese (zh)
Inventor
刘鑫宇
张贝贝
徐鹏
朱彤
徐炜翔
全先江
孟凡文
江浩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangsu Jierui Information Technology Co ltd
716th Research Institute of CSIC
Jiangsu Jari Technology Group Co Ltd
CSIC Information Technology Co Ltd
Original Assignee
Jiangsu Jierui Information Technology Co ltd
716th Research Institute of CSIC
Jiangsu Jari Technology Group Co Ltd
CSIC Information Technology Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jiangsu Jierui Information Technology Co ltd, 716th Research Institute of CSIC, Jiangsu Jari Technology Group Co Ltd, CSIC Information Technology Co Ltd filed Critical Jiangsu Jierui Information Technology Co ltd
Priority to CN202210366703.0A priority Critical patent/CN114936739A/en
Publication of CN114936739A publication Critical patent/CN114936739A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06316Sequencing of tasks or work
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/12Computing arrangements based on biological models using genetic models
    • G06N3/126Evolutionary algorithms, e.g. genetic algorithms or genetic programming
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N7/00Computing arrangements based on specific mathematical models
    • G06N7/01Probabilistic graphical models, e.g. probabilistic networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/04Manufacturing
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Biophysics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computing Systems (AREA)
  • Tourism & Hospitality (AREA)
  • Artificial Intelligence (AREA)
  • General Health & Medical Sciences (AREA)
  • Evolutionary Biology (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Pure & Applied Mathematics (AREA)
  • Primary Health Care (AREA)
  • Manufacturing & Machinery (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Physiology (AREA)
  • Genetics & Genomics (AREA)
  • Biomedical Technology (AREA)
  • Computational Linguistics (AREA)
  • Molecular Biology (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)

Abstract

The invention provides a multi-target flexible job shop scheduling method based on an improved cross entropy algorithm, aiming at solving the problems existing in the multi-target flexible job shop scheduling problem at the present stage, and the method comprises the following steps: initializing iteration times M and a strategy switching threshold value n _ switch, initializing a population P, a process probability distribution matrix, a machine probability distribution matrix and an external memory library E, obtaining a new population according to a process search operator and a machine search operator, removing duplication of individuals with the same target vector by using a mutation operator, screening excellent individuals by adopting an environment selection operator based on SPEA2, and updating the external memory library E by adopting the population P; performing search optimization on the population P and an external memory bank E by hierarchical multi-target neighborhood search, and iterating for M times; and obtaining a final non-dominated solution set by using non-dominated sorting for the external memory bank E. The method can effectively solve the relevant non-dominated solution set, and has better diversity and convergence.

Description

Multi-target flexible job shop scheduling method based on improved cross entropy algorithm
Technical Field
The invention belongs to the technical field of multi-target flexible job shop scheduling, and particularly relates to a multi-target flexible job shop scheduling method based on an improved cross entropy algorithm.
Background
The ship manufacturing industry is one of the reflected national comprehensive strength, concerns the national civilization, the national safety and the national revival, is also a representative of the technology-intensive and labor-intensive enterprises, and belongs to the large-scale equipment manufacturing industry. With the gradual emphasis of the country on ship manufacturing, a strong and comprehensive ship and marine equipment manufacturing system is established so far, and the development of the country in the aspects of ocean development, ocean transportation and the like is promoted.
The shop manufacturing is one of the bases of the ship enterprise production, and the related development and intelligent digitization level thereof closely influence the overall cost and the promotion progress of the project. However, the current domestic workshop intelligent manufacturing management layer is relatively short of the whole scheduling planning of the operation workshop, and is relatively lagged behind in international level. In addition, the relevant top design of the shipyard workshop level is not perfect, and the bottom logic between workshops and projects is not opened. Therefore, under the condition that the current global epidemic situation is repeated to cause the ship market to be in a low state, the ship manufacturing enterprises urgently need to accelerate the intelligent manufacturing of workshops, further promote the production and trust integration of the enterprises and improve the global competitiveness.
The main domestic ship manufacturing plants are investigated and analyzed, so that a large gap exists between the domestic ship manufacturing enterprises and the ship manufacturing industry of developed countries, the degree of informatization of the inter-vehicle level is not high, a subsequent production plan is determined temporarily according to the regular production schedule of every day, every week and the like, and the workshop receives the schedule arrangement of the upper level in a file mode, so that the problem that the production plan is not received timely exists; the resources such as the workshop site, the equipment, the team and the like cannot be effectively utilized, and the production and processing plan is difficult to arrange reasonably.
The flexible workshop scheduling problem is an abstract description of a ship manufacturing process, a single-target flexible workshop scheduling problem is mainly focused on innovation and improvement of an algorithm to obtain a better result, most of researches are focused on three targets of latest completion time, total machine load and maximum machine load for a multi-target flexible workshop scheduling problem, production cost and energy loss are also considered in the research, finally, for the dynamic flexible workshop scheduling problem, the scheduling is a current research hotspot in predicting reaction time, two aspects of production efficiency and stability are mainly considered, the production cost and energy consumption are less considered, the problems of poor diversity and convergence and inaccurate solution exist.
Disclosure of Invention
The invention aims to provide a multi-target flexible job shop scheduling method based on an improved cross entropy algorithm, which can accurately solve the problem of multi-target flexible job shop scheduling and has better convergence, population diversity and robustness.
The technical scheme for realizing the purpose of the invention is as follows: a multi-target flexible job shop scheduling method based on an improved cross entropy algorithm comprises the following steps:
step 1, initializing iteration times M and a strategy switching threshold value n _ switch, initializing a population P, a process probability distribution matrix S, a machine probability distribution matrix Q and an external memory library E, searching by adopting a cross entropy algorithm, obtaining a new population H according to a process search operator and a machine search operator, removing the weight of individuals with the same target vector by using a mutation operator, screening good individuals by adopting an environment selection operator based on SPEA2, and updating the external memory library E by adopting the population P;
step 2, performing search optimization on the population P and an external memory bank E by adopting hierarchical multi-target neighborhood search, and iterating for M times;
and 3, obtaining a final non-dominated solution set by using non-dominated sorting for the external memory library E.
Compared with the prior art, the invention has the following remarkable effects: the cross entropy algorithm is improved by using the coevolution strategy, the problem that the local searching capability of a naive cross entropy algorithm is insufficient is solved, a hierarchical multi-target neighborhood searching strategy is adopted, a random weight mode is used as a replacement condition for a solution during neighborhood searching, and the diversity of a population is enhanced on the premise of ensuring certain convergence of the algorithm; the invention has better superiority and robustness, and can effectively obtain the high-quality solution of the corresponding problem.
Drawings
The accompanying drawings are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the principles of the invention and not to limit the invention.
Fig. 1 is a flowchart of an improved cross entropy algorithm in the present embodiment.
FIG. 2 is a high-quality calamine-solving diagram of example Mk07 in the present embodiment.
Fig. 3 is a high-quality kanji diagram using Mk10 in the present embodiment.
Fig. 4 is a gantt chart corresponding to the extracted graph model in the present embodiment.
Detailed Description
For solving the scheduling problem of the multi-target flexible job shop, a genetic algorithm based on operators such as crossover operators and mutation operators is generally used for solving, due to the blindness of search of the operators, the invention adopts a crossover entropy algorithm for solving, the crossover entropy algorithm belongs to a distribution estimation calculation method, a probability model is maintained, and then the new solution is obtained by sampling the probability model. However, the naive cross entropy algorithm uses a probability model to perform global search, the corresponding local search capability is weak, and the finally obtained result has a further optimization space, so that the invention provides the multi-target flexible job shop scheduling method based on the improved cross entropy algorithm, the method not only can effectively solve the relevant non-dominated solution set, but also has diversity and convergence which are obviously superior to some latest algorithms under certain conditions. The hierarchical multi-target neighborhood search provided by the invention mainly solves two key problems: firstly, the local searching capability of a naive cross entropy algorithm is insufficient; second, consideration of multiple targets in neighborhood searching. The invention mainly comprises the following steps: according to the problem that the local search capability of the cross entropy algorithm is insufficient, a coevolution strategy is used for improving the cross entropy algorithm. The method comprises the following steps of (1) considering the optimization problem of multiple targets for the multi-target flexible job shop scheduling problem: meanwhile, a hierarchical multi-target neighborhood searching strategy is provided, a random weight mode is used as a replacement condition of a solution during neighborhood searching, and the diversity of the population is enhanced on the premise of ensuring certain convergence of the algorithm.
The invention adopts a coevolution strategy to improve a naive cross entropy algorithm, when the algorithm can not obtain a better solution for a plurality of generations, a procedure search operator and a machine search operator are adopted to carry out local search, then hierarchical multi-target neighborhood search is used for further optimization, and when the algorithm can not obtain a better solution for a plurality of generations, a probability model is reused for global search.
The hierarchical multi-target neighborhood search is used for solving the problem that optimization of multiple targets is required to be considered during local search under multi-target flexible job shop scheduling, and is different from single-target flexible job shop scheduling, only the latest completion time is required to be considered during optimization, and the maximum machine load and the total machine load are required to be additionally considered in the multi-target problem. The hierarchical multi-target neighborhood search mainly aims at a multi-target search strategy and a hierarchical search strategy respectively, wherein the multi-target search strategy mainly solves the problem that the local search considers a plurality of targets, the hierarchical search strategy optimizes two sub-problems of the flexible workshop scheduling problem, namely a machine allocation sub-problem and a process sequence sub-problem, and meanwhile, the neighborhood search based on resource switching and the neighborhood search based on idle time are used and have a complementary relationship, wherein the former mainly considers the idle time of other machines, and the latter mainly considers the idle time of the current machine.
In order to achieve the objective of the present invention, the improved cross-entropy algorithm (Co-CEM) based on hierarchical multi-target neighborhood search of the present invention, described with reference to fig. 1 based on the above problem, includes the following steps:
and 001, initializing iteration times M, a strategy switching threshold value n _ switch, a population P, a process probability distribution matrix S, a machine probability distribution matrix Q and an external memory library E.
And 002, carrying out search optimization on the population P and the external memory bank E, and carrying out iteration for M times in total.
And step 003, obtaining a final non-dominated solution set by using non-dominated sorting to the external memory bank E.
The step 001 specifically comprises the following steps:
step 00101, initializing a population P: the method comprises a machine allocation part and a process sequence part for initializing solution vectors in the population P. The machine allocation portion of the solution vector is initialized using 4 machine allocation rules, including 2 global initialization rules, 1 local allocation rule, and 1 random allocation rule. The detailed description is as follows:
(1) global minimum load: selecting the machine with the minimum load from all the machines for distribution each time;
(2) random arrangement: firstly, randomly sequencing process vectors, sequentially traversing the process vectors from left to right, and selecting a machine with the minimum load for each process;
(3) minimum processing time: the working procedure is distributed to the machine with the shortest processing time;
(4) random distribution: each procedure randomly distributes corresponding processing machines;
for the process part in the population P initialization solution vector, the invention uses 3 sort order rules, which are described in detail as follows:
(1) at most remaining man-hours: and performing process distribution according to the remaining working hours of each task. Selecting a task with the minimum remaining working hours from all tasks to distribute the working procedure of the task, and deleting the task from the task set to be distributed if the next working procedure does not exist; randomly selecting one if there are a plurality of tasks with the smallest remaining man-hours;
(2) the most remaining processes: and performing process distribution according to the residual number of processes of each task. Selecting a task with the least number of processes from all the tasks to distribute the process, and if the next process does not exist, deleting the task from the task set to be distributed; randomly selecting one task if a plurality of tasks with the minimum number of processes exist;
(3) random sequencing: randomly distributing the processing sequence of all the working procedures;
00102, initializing a probability distribution matrix (the probability distribution matrix comprises a process probability distribution matrix and a machine probability distribution matrix). Let S be an element of R N×N Is a process probability distribution matrix for generating process vectors, Q ∈ R N×M For a machine probability distribution matrix used to generate machine vectors, where N is the total number of processes and M is the total number of machines, then S [ x, y]The probability of assigning process y to the x-th position of the process vector is shown, and S [ x, y ] is given in the initial state to ensure that the whole solution space can be uniformly sampled]=1/N。Q[x,u]Probability of assigning machine u to the x-th position of machine vector, where m y The number of machines that can process step y is initially as shown in the following equation:
Figure BDA0003587442920000041
because the precedence relationship between the working procedures of the same task is restricted, the working procedures are towardsThe x-th position of the quantity selectable process cannot be taken arbitrarily from 1, 2.., N, but is related to the process selected by the process vector subscript 0 to x-1. Also, due to the inherent relationship between machine selection and process selection, the machine vector is constrained by the same constraint that the x-th location-selectable machine cannot be taken from all machines at will, but is instead associated with the set of machines that correspond to the process selectable at that location. Aiming at the constraint condition, a procedure vector order _ mask and a machine vector machine _ mask are respectively introduced. Taking an order _ mask vector as an example, defining the order _ mask as an element R N All the processes with selectable current positions are represented, and the values are shown in the following formula:
Figure BDA0003587442920000051
when a process is selected for the position x of a process vector, an order _ mask is used for screening to obtain a probability distribution vector sample _ vec, then the probability sum is ensured to be 1 through processing, so that the generated process is a feasible solution, and the calculation step is shown as follows:
sample_vec=P[x,:]*order_mask
sample_vec=sample_vec/sum(sample_vec)
sum () represents summation, and once the process vector is determined to the process of the x-th position, the order _ mask is updated to ensure that the subsequently generated process is valid.
Similar processing is adopted for the generation step of the machine vector. Useless solutions in the sampling process are removed through the machine _ mask, and the generated solutions are guaranteed to be effective.
After a new population is obtained through sampling, the cross entropy algorithm selects n _ elite individuals according to the fitness to update a process probability distribution matrix S and a machine probability distribution matrix Q, as shown in the following formula:
Figure BDA0003587442920000052
Figure BDA0003587442920000053
wherein alpha and beta are the learning rates of S and Q, respectively,
Figure BDA0003587442920000054
and
Figure BDA0003587442920000055
the conditional function for the kth solution in the population is shown in the following formula:
Figure BDA0003587442920000056
Figure BDA0003587442920000057
00103 initializing external memory bank
Figure BDA0003587442920000058
Indicating that the external memory pool is an empty set.
00104, obtaining a new population H according to the procedure search operator and the machine search operator. For the process search operator, if parent individuals m and n exist, the detailed steps are as follows:
step 1: firstly, judging whether m is the same as n, if so, randomly exchanging two working procedures which do not belong to the same task to obtain a new individual and directly returning; otherwise, executing step 2;
step 2: setting the current task number as N, numbering the N tasks, randomly selecting a plurality of tasks from the N tasks to a set u, defining v as a complementary set of the set u, and reselecting if the selected task number is 0 or equal to N;
and step 3: the processes belonging to the set u in the offspring p individuals are sequentially copied to corresponding positions from the parent individuals m, and the remaining processes belonging to the set v are sequentially copied to corresponding positions from the parent individuals n;
and 4, step 4: conversely, the processes belonging to the set u in the child q are sequentially copied from the parent n to the corresponding positions, and the remaining processes belonging to the set v are sequentially copied from the parent m to the corresponding positions.
When the process searching is carried out, in order to avoid obtaining equivalent offspring individuals from the same parent individuals, similarity checking is firstly carried out on the parent individuals, and if the parent individuals and the parent individuals are completely the same, two processes are directly and randomly exchanged to obtain new individuals. The detailed steps of the corresponding machine search operator are as follows:
step 1: and setting the current process number as N, randomly selecting k channels from the N processes to form a set u, and if the number of elements in the set u is 0 or equal to N, reselecting.
Step 2: and (4) exchanging the assigned machines of the parent individuals m and n belonging to the process in the set u to obtain the offspring individuals p and q.
Step 00105, let the population P ← P ≧ H, use mutation operator to deduplicate to the individual with the same target vector, wherein the mutation operator or local search operator used is as follows:
(1) maximum load based machine reallocation operator: from the process of machining on the machine with the largest load, one is randomly selected and redistributed to another machine with a relatively smaller load for machining.
(2) The machine reassignment operator based on the latest completion time: randomly selecting a process on the machine with the latest completion time, and redistributing the process to another machine with a smaller latest completion time for processing.
(3) The key process moves: randomly acquiring a key path for a current individual, randomly finding a key block on the current key path, and randomly selecting a working procedure to move forwards or backwards to the head or tail of the key block.
(4) Key procedure reinsertion: a key process is randomly selected and moved to another processing facility.
00106, screening excellent individuals by using an environment selection operator based on SPEA2, and assigning the excellent individuals to a population P. Since SPEA2 uses density information in the population to better preserve population diversity, it is used as an environment selection operator in the present invention. SPEA2 incorporates some means of ensuring population diversity, such as a specific fitness allocation strategy, density estimation algorithms and modified truncation methods, which first calculate the dominance between each pair of individuals in the population, and then define the associated intensity values str (x), as shown below:
Str(x i )=|{j|j∈P∧x i <x j }|
wherein P represents a population, x i And x j Respectively, the individuals in the population P, Str (x) i ) Representing the number of other individuals for whom the intensity value of an individual dominates. Then, an original fitness value raw (x) is defined, which is defined as follows:
Figure BDA0003587442920000061
representing the sum of intensity values that dominate the current individual. In addition, an additional density information value Den (x) is defined for identifying individuals having the same Raw value, as shown in the following formula:
Figure BDA0003587442920000071
wherein k is the square root of the population size,
Figure BDA0003587442920000072
denotes x i And the Euclidean distance to other k-th near individuals in the population, wherein the final fitness value is expressed as the sum of the original fitness value and the density information value, and is shown as the following formula:
Fitness(x i )=Raw(x i )+Den(x i )
when the environment selection operator based on SPEA2 is used, firstly, individuals with the fitness value smaller than 1 are preferably selected from candidate solutions, namely, non-dominant solutions, and if the number of the non-dominant solutions is not enough, the individuals with the proper number are sequentially selected according to the fitness value; otherwise, corresponding stage steps are executed to eliminate redundant individuals, namely individuals with smaller Euclidean distance.
And 00107, updating the external memory bank E by using the population P. The updating process of the memory bank is as follows:
step 1: comparing the objective function vector of each good individual with each individual in an external memory bank;
step 2: replace if the individual in the external memory bank is dominant;
and step 3: if the objective function vectors of the two are the same, calculating the Hamming distance of the machine vector, if the Hamming distance is 0, using a variation strategy until different objective function values are contained, and checking whether an alternative individual exists; otherwise, the backward comparison is continued.
The step 002 specifically includes the following steps:
00201, selecting an evolution strategy according to a strategy identifier flag: by sampling the process matrix S and the machine matrix Q or using dynamic crossover operators. The diversity of the population can be enhanced by generating offspring individuals based on a dynamic crossover probability mode, and the corresponding crossover probability gradually changes along with the increase of the iteration times. If the current Iteration number is Iteration and the Total Iteration number is Total _ Iteration, the corresponding crossover probability is P c 1-Iteration/Total _ Iteration, which represents the probability P of selecting one individual from the memory bank and the population to intersect at present c All of them select individuals from the current population to cross with a probability of 1-P c
00202, carrying out local search on Q by using a hierarchical multi-target neighborhood search strategy to obtain a new individual K. First, an extraction graph model and related symbols are defined. A disjunctive graph model was originally proposed and applied to job shop scheduling, which is a directed acyclic graph. First, a description FJSP is defined in which (V, U, E) denotes a set of nodes formed by all the processes, and a virtual start node s and an end node E are determined. U denotes a set of all joint arcs (joint edges) that determine the machining priority between the processes. E represents the set of all disjunctive arcs (disjunctive edges) and satisfies
Figure BDA0003587442920000081
Figure BDA0003587442920000082
The set of all the extracted arcs on the kth machine is shown, the corresponding processing time is shown below the corresponding node, and the processing time of the virtual start node and the processing time of the virtual end node are both 0.
Fig. 4 shows an extracted graph obtained after determining the allocation machine for a feasible solution and a gantt chart corresponding to the extracted graph. The method is characterized in that a plurality of symbols are defined based on the disjunctive graph, so that the hierarchical multi-target neighborhood searching strategy can be conveniently described later. Defining the node on the analysis graph G to correspond to the process j under the task i, and using the symbol O i,j Denotes that the machining machine is k and the corresponding machining time is t i,j,k
Definition of
Figure BDA0003587442920000083
To represent a process O i,j At the earliest start time on the machine k,
Figure BDA0003587442920000084
represents the step O i,j Without delaying C on machine k max The latest start time allowed under the conditions of (1). Accordingly, define
Figure BDA0003587442920000085
Is a process O i,j At the earliest end time on the machine k,
Figure BDA0003587442920000086
is a process O i,j Not postponing C on machine k max And the latest end time allowed under the condition (1) and satisfies the following formula:
Figure BDA0003587442920000087
Figure BDA0003587442920000088
order to
Figure BDA0003587442920000089
For working on machine k i,j The pre-processing procedure of (2) is carried out,
Figure BDA00035874429200000810
for working on machine k i,j And (4) post-processing. Let PJ i,j (G) Is a process O i,j A precursor step of (4), SJ i,j (G) Is a process O i,j The subsequent processes of (2). Process O i,j Is a key process and only
Figure BDA00035874429200000811
In the directed graph G, the critical path is a path formed by adjacent critical processes. As shown in FIG. 4, there is a critical path (s → O) 2,1 →O 1,1 →O 1,2 →O 1,3 →O 3,3 → e). The key process comprises the following steps: { O 2,1 ,O 1,1 ,O 1,2 ,O 1,3 ,O 3,3 }, key block: { (O) 2,1 ,O 1,1 ),(O 1,2 ),(O 1,3 ,O 3,3 )}. The makespan for this schedule is 16.
A resource switch based neighborhood search is then defined. Defining process ω, the movable time period is the predecessor process PJ ω The earliest completion time and the subsequent process SJ ω The latest start time of (c):
Figure BDA00035874429200000812
then releasing and switching process ω from the current resource to the appropriate resource may result in an improved solution.
Let k be the current processing machine of process ω and t be the corresponding processing time ω,k The set of machinable machines defining process ω is Ψ ω To, for
Figure BDA00035874429200000813
Has working procedures p and q processed on a machine m, and
Figure BDA00035874429200000814
the process ω can reduce the latest completion time after the resource switch if and only if the earliest end time of p and the latest start time of q satisfy:
Figure BDA0003587442920000091
if the above equation is not satisfied, the process ω cannot be shifted after the resource switching because the latest completion time cannot be reduced.
Then the critical process ω is defined as well. There is an optimum insertion position if and only if the following:
Figure BDA0003587442920000092
wherein the content of the first and second substances,
Figure BDA0003587442920000093
g is obtained by deleting key processes of claim movement from the graph G - Procedure (2)
Figure BDA0003587442920000094
The earliest time of completion of the time,
Figure BDA0003587442920000095
the same is true.
Figure BDA0003587442920000096
Shows the step v at G - In the middle, the original latest completion time C is not delayed max (G) The latest start time under the conditions of (a),
Figure BDA0003587442920000097
the same is true. t is t ω,k Representing the processing time of the critical process omega on the machine k.
Finally, define neighborhood based on idle timeAnd (6) searching. The processes p and q are adjacent to each other in FIG. G, and the processing machines are m. If the critical process ω is a machining machine m, then ω has an exchangeable time period of
Figure BDA0003587442920000098
The idle time period of the process p, q is
Figure BDA0003587442920000099
The critical process w can then be inserted into the free time period of the processes p, q by a process transformation if and only if the following equation is satisfied:
Figure BDA00035874429200000910
that is, if the intersection of the two time periods is not 0, then shifting the critical process ω to a new location may reduce the latest completion time.
The hierarchical multi-target neighborhood search strategy comprises two parts: a multi-target oriented search strategy and a hierarchical search strategy.
(1) Multi-objective oriented search strategy
Because multiple indexes are always considered at the same time in the multi-objective FJSP, and the latest completion time is the most difficult to optimize, a multi-objective search strategy is proposed in the patent. Firstly, for the scheduling solution G, psi (G) ═ co is defined 1 ,co 2 ,...,co n Is the set of all key process steps, and defines pi (G) ═ co i →M k I 1, 2.. n is a set of all neighborhood actions, called action set, wherein each neighborhood action is composed of two elements of a process and a target machine, which represents the redistribution of a key process to another machine, wherein n represents the number of key processes, M k Representing the target machine. To pair
Figure BDA00035874429200000911
The total machine load variation Δ t and the maximum machine load Δ c can be expressed as follows:
Figure BDA0003587442920000101
Figure BDA0003587442920000102
wherein m is * Representing target machine, m representing source machine, t co,m* Indicating machine m * Time required for working the process, t co,m Indicates the processing time required by the machine m,
Figure BDA0003587442920000103
represents m * The total load of (c).
The multi-target-oriented search strategy mainly comprises two stages, namely neighborhood search based on resource switching and neighborhood search based on key procedures. First, since the neighborhood search for resource switching guarantees that the resulting solution may be improved at the latest completion time, it is preferentially used as the first-stage neighborhood search. Then, according to a grading strategy, all actions corresponding to the current scheduling solution G are sorted in a non-descending order according to the Δ t and the Δ c, so that the smaller Δ t is considered preferentially and the Δ c is considered secondly under the condition that the latest completion time is possibly improved. Then, sequentially using a neighborhood search based on resource switching in a first stage for the action set pi (G) obtained after sorting, and if a feasible action is found, ending the search; otherwise, using the neighborhood search based on key process in the second stage to finish the search after finding a feasible action.
(2) Search strategy based on hierarchy
The hierarchical search strategy is mainly used for respectively solving two sub-problems of FJSP machine allocation and process sequence, wherein the first layer uses a multi-target-oriented search strategy to solve the sub-problems of machine allocation, and the second layer uses a neighborhood search strategy based on idle time to solve the sub-problems of process sequence. In addition, the replacement of the old solution by the new solution is accomplished using an aggregation function approach that randomly generates weight vectors: randomly generating a weight vector lambda according to a specific method, and replacing if the dot product of the target vector of the new solution and lambda is smaller than the dot product of the old solution and lambda; otherwise, the replacement is not performed.
The first layer is optimized aiming at idle time of a cross machine, the second layer uses neighborhood search based on the idle time, and the processing sequence of the current working procedure is further optimized on the basis of the first layer of results. Since only the processing sequence between the processes is changed and the machines to which the processes are assigned are not changed, only the latest completion time is affected. The former layer aims at the utilization of larger idle time of a machine in the current scheduling, while the current layer aims at the utilization of smaller idle time of the machine, and a complementary relationship exists between the former layer and the latter layer, so that the local optimization capability of the algorithm can be effectively enhanced.
Step 00203, let P ← P ^ H ^ K, and use mutation operator to remove repeated individual.
Step 00204.P ← screening good individuals for P using SPEA2 based context selection operator.
Step 00205.E '. ae' uses a hierarchical multi-target neighborhood search strategy to perform local search on E.
Step 00206.E ← E'. E, and removing repeated individuals using a mutation operator.
Step 00207.E ← screen for E good individuals using SPEA2 based context selection operator.
Step 00208.E ← updates the external memory bank E with P.
Step 00209.S, Q ← updating the probability matrix using the external memory bank in accordance with flag.
Step 00210. strategy identification flag ← switching evolution strategy if the continuous n _ switch generation does not update the external memory base.
The use of non-dominated sorting to the external memory bank E in step 003 to obtain the final non-dominated solution set is a method well known in the art and will not be described in detail herein.
To verify the effectiveness of the Co-CEM method of the invention, a public data set was used: the Kacem data set comprises 5 use cases and the BRdata use case data set comprises 10 use cases. The number of machines and the number of tasks and specific information of the working procedure under each task are determined by each use case. And respectively comparing Co-CEM with hDPSO, DABC, BEG-NSGA-II and INSBBO algorithms by using the contrast indexes, namely the hyper-volume HV and the reverse generation distance IGD, wherein the bold shows that the contrast algorithms are optimal, and the result is shown in the following table.
TABLE 1 Co-CEM comparison with other algorithm results
Figure BDA0003587442920000111
First, by comparing HV indicators, it can be seen that Co-CEM achieves relatively optimal results in 5 cases, in total, exceeding hDPSOA in 8 cases, BEG-NSGA-II in 9 cases, INSBBO in 4 cases, and DABC in 1 case. It is noteworthy that in the largest data size Mk10 case, Co-CEM is significantly better than all algorithms compared
Then, comparing IGD indexes, it can be seen that hDPSO is optimal in 2 cases, DABC and INSBBO are superior to other algorithms in 3 cases, and Co-CEM obtains results superior to other algorithms in 6 cases. Meanwhile, Co-CEM exceeds hDPSO in 8 cases, BEG-NSGA-II in all cases, DABC in 2 cases and INSBBO in 5 cases.
In conclusion, the experiment proves that the Co-CEM has certain superiority compared with other algorithms, and can effectively obtain a high-quality solution of the corresponding problem.
The optimization algorithm solves the traditional scheduling problem, thereby ensuring the stability and effectiveness of the production plan of the domestic ship manufacturing enterprise and reducing the blindness of the plan. Except for ship manufacturing, other fields such as power systems and medical resource distribution systems have the above requirements, so that the traditional scheduling problem is solved based on the flexible workshop scheduling model, and the method has important practical significance for the actual production and manufacturing such as national industry 4.0, China manufacturing 2025 and other plans and the people life aspects such as medical systems and power systems.

Claims (10)

1. A multi-target flexible job shop scheduling method based on an improved cross entropy algorithm is characterized by comprising the following steps:
step 1, initializing iteration times M and a strategy switching threshold value n _ switch, initializing a population P, a process probability distribution matrix S, a machine probability distribution matrix Q and an external memory library E, searching by adopting a cross entropy algorithm, obtaining a new population H according to a process search operator and a machine search operator, removing the weight of individuals with the same target vector by using a mutation operator, screening good individuals by adopting an environment selection operator based on SPEA2, and updating the external memory library E by adopting the population P;
step 2, performing search optimization on the population P and an external memory bank E by adopting hierarchical multi-target neighborhood search, and iterating for M times;
and 3, obtaining a final non-dominated solution set by using non-dominated sorting to the external memory library E.
2. The multi-target flexible job shop scheduling method according to claim 1, wherein the initializing the population P in step 1 specifically includes initializing a machine allocation vector in a population P initialization solution vector and initializing a process vector in the population P initialization solution vector;
machine allocation vector initialization in the population P initialization solution vector is initialized by adopting 4 machine allocation rules, wherein the 4 machine allocation rules specifically comprise:
(1) global minimum load distribution rule: selecting the machine with the minimum load from all the machines for distribution each time;
(2) random permutation assignment rule: firstly, randomly sequencing process vectors, sequentially traversing the process vectors from left to right, and selecting a machine with the minimum load for each process;
(3) minimum processing time allocation rule: the working procedure is distributed to the machine with the shortest processing time;
(4) random allocation rule: each procedure randomly distributes corresponding processing machines;
the initialization of the process vector in the population P initialization solution vector adopts 3 sort rules for initialization, wherein the 3 sort rules specifically are as follows:
(1) most remaining man-hour ordering rules: performing procedure distribution according to the remaining working hours of each task, selecting a task with the minimum remaining working hours from all the tasks to distribute the working procedure of the task, and deleting the task from the task set to be distributed if the next working procedure does not exist; randomly selecting one task if there are a plurality of tasks with the smallest remaining man-hours;
(2) at most, the remaining process ordering rules: performing process distribution according to the remaining process number of each task, selecting one task with the least process number from all the tasks to distribute the process, and deleting the task from the task set to be distributed if the next process does not exist; if a plurality of tasks with the minimum number of processing steps exist, one task is randomly selected;
(3) random ordering rule: and randomly distributing the processing sequence of all the procedures.
3. The multi-target flexible job shop scheduling method according to claim 1, wherein the process probability distribution matrix S and the machine probability distribution matrix Q are initialized as follows:
and S is formed as R N×N For the process probability distribution matrix used to generate the process vector, Q ∈ R N×M For a machine probability distribution matrix used to generate machine vectors, where N is the total number of processes and M is the total number of machines, then S [ x, y]The probability of assigning the process y to the x-th position of the process vector is represented by S [ x, y ] in the initial state]=1/N,Q[x,u]Probability of assigning machine u to the x-th position of machine vector, where m y The number of machines that can process the process y is shown, and the initial state is:
Figure FDA0003587442910000021
the external memory bank E is initialized to an empty set.
4. The multi-target flexible job shop scheduling method according to claim 1, wherein the obtaining of the new population according to the process search operator specifically comprises the steps of:
judging whether m is the same as n, if so, randomly exchanging two working procedures which do not belong to the same task to obtain a new individual and directly returning;
the current task number is set as N, the N tasks are numbered, a plurality of tasks are randomly selected from the N tasks to a set u, v is defined as a complementary set of the set u, and if the selected task number is 0 or equal to N, the selection is carried out again;
the processes belonging to the set u in the offspring p individuals are sequentially copied to corresponding positions from the parent individuals m, and the remaining processes belonging to the set v are sequentially copied to corresponding positions from the parent individuals n;
on the contrary, the processes belonging to the set u in the offspring q individuals are sequentially copied to the corresponding positions from the parent individuals n, and the remaining processes belonging to the set v are sequentially copied to the corresponding positions from the parent individuals m;
the step of obtaining a new population according to the machine search operator specifically comprises the following steps:
setting the current process number as N, randomly selecting k channels from the N processes to form a set u, and if the number of elements in the set u is 0 or equal to N, reselecting;
and (4) exchanging the assigned machines of the parent individuals m and n belonging to the process in the set u to obtain the offspring individuals p and q.
5. The multi-target flexible job shop scheduling method according to claim 1, wherein the de-duplication using mutation operators for individuals having the same target vector specifically comprises:
maximum load based machine reassignment operator: randomly selecting one of the processing procedures from the machine with the maximum load, and redistributing the selected processing procedure to the other machine with relatively smaller load for processing;
the machine reassignment operator based on the latest completion time: randomly selecting a process on a machine with the latest completion time, and redistributing the process to another machine with smaller latest completion time for processing;
the key process moves: randomly acquiring a key path for a current individual, randomly finding a key block on the current key path, and randomly selecting a working procedure to move forwards or backwards to the head or tail of the key block;
the key process reinsertion: a key process is randomly selected and moved to another processing facility.
6. The multi-target flexible job shop scheduling method according to claim 1, wherein the screening of good individuals using an environment selection operator based on SPEA2 specifically comprises: preferably selecting individuals with the fitness value smaller than 1 from the candidate solutions, namely non-dominated solutions, and if the number of the non-dominated solutions is not enough, sequentially selecting the individuals with the proper number according to the fitness value; otherwise, the redundant individuals are removed, namely the individuals with smaller Euclidean distance; the fitness is as follows:
Fitness(x i )=Raw(x i )+Den(x i )
Figure FDA0003587442910000031
Figure FDA0003587442910000032
wherein P represents a population, x i And x j Respectively, the individuals in the population P, Str (x) i ) Representing the number of other individuals for which the intensity value of an individual dominates, raw (x) the original fitness value, den (x) the additional density information value, k the square root of the population size,
Figure FDA0003587442910000033
denotes x i Euclidean distance to other kth near individuals in the population, Fitness (x) i ) Is the fitness.
7. The multi-target flexible job shop scheduling method according to claim 1, wherein the updating the external memory library E with the population P specifically comprises:
comparing the target function vector of each good individual of the population P with the target function vector of each individual in the external memory bank;
replace if the individual in the external memory bank is dominant;
if the objective function vectors of the two are the same, calculating the Hamming distance of the machine vector, if the Hamming distance is 0, using a variation strategy until different objective function values are contained, and checking whether an alternative individual exists; otherwise, the backward comparison is continued.
8. The multi-target flexible job shop scheduling method according to claim 1, wherein the step 2 specifically comprises:
step 2-1, selecting an evolution strategy according to a variable flag: obtaining a new population H by sampling the process matrix S and the machine matrix Q or using a dynamic crossover operator; generating offspring individuals based on a dynamic cross probability mode, setting the current Iteration number as Iteration and the Total Iteration number as Total _ Iteration, and then setting the corresponding cross probability as P c 1-Iteration/Total _ Iteration, which represents the probability P of selecting one individual from the memory bank and the population to intersect at present c All of them select individuals from the current population to cross with a probability of 1-P c
Step 2-2, carrying out local search on the H by using a hierarchical multi-target neighborhood search strategy to obtain a new individual K;
step 2-3, updating P ← P £ H $, and removing repeated individuals by using a mutation operator;
step 2-4, updating P ← screening good individuals for P by using an environment selection operator based on SPEA 2;
step 2-5, updating E' ← and using a hierarchical multi-target neighborhood search strategy to perform local search on E;
step 2-6, updating E ← E' U E, and removing repeated individuals by using a mutation operator;
step 2-7, updating E ← screening good individuals for E by using an environment selection operator based on SPEA 2;
step 2-8, updating E ← updating external memory bank E by using P;
step 2-9, updating S, Q ← updating probability matrix using external memory base according to flag;
and 2-10, updating a variable flag ← switching the evolution strategy if the continuous strategy switching threshold n _ switch generation does not update the external memory base.
9. The multi-target flexible job shop scheduling method according to claim 1, wherein the hierarchical multi-target neighborhood search strategy comprises a multi-target-oriented search strategy and a hierarchical search strategy, and the multi-target-oriented search strategy comprises:
firstly, for the scheduling solution G, psi (G) ═ co is defined 1 ,co 2 ,...,co n Is the set of all key process steps, and defines pi (G) ═ co i →M k I 1, 2.. n is a set of all neighborhood actions, called action set, wherein each neighborhood action is composed of two elements of a process and a target machine, which represents the redistribution of a key process to another machine, wherein n represents the number of key processes, M k Represents a target machine, pair
Figure FDA0003587442910000041
The total machine load variation Δ t and the maximum machine load Δ c are respectively expressed by the following formulas:
Figure FDA0003587442910000042
Figure FDA0003587442910000043
wherein m is * Representing the target machine, m representing the source machine,
Figure FDA0003587442910000044
showing the machine m * Time required for working the process, t co,m Indicates the processing time required by the machine m,
Figure FDA0003587442910000045
represents m * The total load of (c);
performing neighborhood search based on resource switching;
respectively sequencing the delta t and the delta c in a non-descending order for all actions corresponding to the current scheduling solution G, so that the actions with smaller delta t are considered preferentially and the delta c is considered secondarily under the condition that the latest completion time is possibly improved;
sequentially using neighborhood search based on resource switching to the action set pi (G) obtained after sorting, and if a feasible action is found, finishing the search; otherwise, using neighborhood search based on key process, and ending the search after finding a feasible action;
the hierarchical search strategy specifically comprises the following steps: optimizing the idle time of the cross machine by adopting a multi-objective-oriented search strategy, optimizing the processing sequence of the current procedure on the basis of the idle time optimization result on the basis of the neighborhood search of the idle time, and acquiring a polymerization function mode of randomly generating a weight vector to complete the replacement of the old solution by the new solution: randomly generating a weight vector lambda, and replacing if the dot product of the target vector of the new solution and the lambda is smaller than the dot product of the old solution and the lambda; otherwise, the replacement is not carried out.
10. The multi-target flexible job shop scheduling method according to claim 9, wherein the neighborhood search based on resource switching is specifically:
defining the process ω such that its movable time period is the predecessor process PJ ω Earliest completion time and subsequent process SJ ω The latest start time of (c):
Figure FDA0003587442910000051
let k be the current processing machine of process ω and t be the corresponding processing time ω,k Machinable machine defining a process ωSet to Ψ ω To, for
Figure FDA0003587442910000052
Has working procedures p and q processed on a machine m, and
Figure FDA0003587442910000053
Figure FDA0003587442910000054
if and only if the earliest end time of p and the latest start time of q satisfy:
Figure FDA0003587442910000055
when the above formula does not hold, the process ω is not moved;
the neighborhood search based on idle time specifically comprises:
if the processes p and q are adjacent processes in (V, U, E) and the processing machines are m, and if the processing machine is m, the exchangeable time period of ω is m
Figure FDA0003587442910000056
The idle time period of the process p, q is
Figure FDA0003587442910000057
The critical process w can then be inserted into the free time period of the processes p, q by a process transformation if and only if the following equation is satisfied:
Figure FDA0003587442910000058
if the intersection of the two time periods is not 0, then the key procedure omega is transformed to a new position; wherein V represents a set of nodes formed by all the processes, and a virtual start node s and a virtual end node e are determined at the same time, and U represents a set formed by all conjunctive arcsDetermining the processing priority between the processes; e represents the set of all disjunct arcs and satisfies
Figure FDA0003587442910000059
Figure FDA00035874429100000510
Representing a set of all disjunct arcs on a kth machine, representing corresponding processing time below corresponding nodes, wherein the processing time of the virtual start node and the virtual end node is 0;
the neighborhood search based on the key process is specifically as follows:
defining the critical process step ω, there is an optimal insertion position if and only if the following equation is satisfied:
Figure FDA0003587442910000061
wherein the content of the first and second substances,
Figure FDA0003587442910000062
g is obtained by deleting the key process of the requested movement from the graph G - Procedure (2)
Figure FDA0003587442910000063
The earliest time of completion of the time,
Figure FDA0003587442910000064
in the same way, the method has the advantages of,
Figure FDA0003587442910000065
shows the step v at G - In the middle, the original latest completion time C is not delayed max (G) The latest start time under the conditions of (a),
Figure FDA0003587442910000066
the same process is carried out; t is t ω,k Representing the processing time of the critical process omega on the machine k.
CN202210366703.0A 2022-04-08 2022-04-08 Multi-target flexible job shop scheduling method based on improved cross entropy algorithm Pending CN114936739A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210366703.0A CN114936739A (en) 2022-04-08 2022-04-08 Multi-target flexible job shop scheduling method based on improved cross entropy algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210366703.0A CN114936739A (en) 2022-04-08 2022-04-08 Multi-target flexible job shop scheduling method based on improved cross entropy algorithm

Publications (1)

Publication Number Publication Date
CN114936739A true CN114936739A (en) 2022-08-23

Family

ID=82863399

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210366703.0A Pending CN114936739A (en) 2022-04-08 2022-04-08 Multi-target flexible job shop scheduling method based on improved cross entropy algorithm

Country Status (1)

Country Link
CN (1) CN114936739A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115826537A (en) * 2023-01-29 2023-03-21 广东省科学院智能制造研究所 Flexible scheduling method for multi-robot production line

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115826537A (en) * 2023-01-29 2023-03-21 广东省科学院智能制造研究所 Flexible scheduling method for multi-robot production line
CN115826537B (en) * 2023-01-29 2023-05-02 广东省科学院智能制造研究所 Flexible scheduling method for multi-robot production line

Similar Documents

Publication Publication Date Title
CN107590603B (en) Based on the dispatching method and system for improving change neighborhood search and differential evolution algorithm
WO2021036658A1 (en) Multi-objective optimization method and system for master production schedule of casting parallel workshops
Segura et al. A novel diversity-based replacement strategy for evolutionary algorithms
CN110738365B (en) Flexible job shop production scheduling method based on particle swarm algorithm
CN110288185B (en) Distributed flexible pipeline scheduling method
Chen et al. A fast memetic multi-objective differential evolution for multi-tasking optimization
Bouabda et al. A branch and bound enhanced genetic algorithm for scheduling a flowline manufacturing cell with sequence dependent family setup times
CN113569483A (en) Method for solving multi-target flexible job shop scheduling based on artificial bee colony algorithm
CN114936739A (en) Multi-target flexible job shop scheduling method based on improved cross entropy algorithm
He et al. Differential evolution with adaptive mutation and parameter control using Lévy probability distribution
Kumar et al. Proposed selection technique of evolutionary algorithm and its implementation for combinatorial problems
CN115271163A (en) Strategy gradient-based vehicle path optimization method based on super heuristic algorithm
CN116985146B (en) Robot parallel disassembly planning method for retired electronic products
Tavakkoli-Moghaddam et al. A genetic algorithm and memetic algorithm to sequencing and scheduling of cellular manufacturing systems
Yu et al. A route planning strategy for the automatic garment cutter based on genetic algorithm
CN116070761A (en) Project scheduling rule mining method and system based on gene expression programming
Rifai et al. Multi-operator hybrid genetic algorithm-simulated annealing for reentrant permutation flow-shop scheduling
CN111665799A (en) Time constraint type parallel machine energy-saving scheduling method based on collaborative algorithm
Wang et al. A tailored NSGA-III for multi-objective flexible job shop scheduling
Zaman et al. Evolutionary algorithm for project scheduling under irregular resource changes
Dominik Solving multi-objective job shop problem using nature-based algorithms: new Pareto approximation features
CN112884370B (en) Multi-project rescheduling method and system for high-end equipment development process considering order insertion list
Wang et al. An estimation of distribution algorithm for the flexible job-shop scheduling problem
Weise et al. Frequency fitness assignment: Making optimization algorithms invariant under bijective transformations of the objective function value
Luo et al. Solving Flexible Job-Shop Problem with Sequence-Dependent Setup Times by Using Selection Hyper-Heuristics

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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 222061 No.18, Shenghu Road, Haizhou District, Lianyungang City, Jiangsu Province

Applicant after: JIANGSU JARI TECHNOLOGY GROUP Co.,Ltd.

Applicant after: Jiangsu Jierui Information Technology Co.,Ltd.

Applicant after: The 716th Research Institute of China Shipbuilding Corp.

Applicant after: CSIC Information Technology Co.,Ltd.

Address before: 222061 No.18, Shenghu Road, Haizhou District, Lianyungang City, Jiangsu Province

Applicant before: JIANGSU JARI TECHNOLOGY GROUP Co.,Ltd.

Applicant before: Jiangsu Jierui Information Technology Co.,Ltd.

Applicant before: 716TH RESEARCH INSTITUTE OF CHINA SHIPBUILDING INDUSTRY Corp.

Applicant before: CSIC Information Technology Co.,Ltd.