Next Article in Journal
Hydrodynamic and Acoustic Performance Analysis of Marine Propellers by Combination of Panel Method and FW-H Equations
Next Article in Special Issue
Variation Rate to Maintain Diversity in Decision Space within Multi-Objective Evolutionary Algorithms
Previous Article in Journal
Tensors in Newtonian Physics and the Foundations of Classical Continuum Mechanics
Previous Article in Special Issue
Pool-Based Genetic Programming Using Evospace, Local Search and Bloat Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improved Differential Evolution Algorithm for Flexible Job Shop Scheduling Problems

by
Prasert Sriboonchandr
*,
Nuchsara Kriengkorakot
and
Preecha Kriengkorakot
Industrial Engineering, Department, Faculty of Engineering, Ubon Ratchathani University, Ubon Ratchathani 34190, Thailand
*
Author to whom correspondence should be addressed.
Math. Comput. Appl. 2019, 24(3), 80; https://doi.org/10.3390/mca24030080
Submission received: 20 July 2019 / Revised: 19 August 2019 / Accepted: 5 September 2019 / Published: 6 September 2019
(This article belongs to the Special Issue Numerical and Evolutionary Optimization)

Abstract

:
This research project aims to study and develop the differential evolution (DE) for use in solving the flexible job shop scheduling problem (FJSP). The development of algorithms were evaluated to find the solution and the best answer, and this was subsequently compared to the meta-heuristics from the literature review. For FJSP, by comparing the problem group with the makespan and the mean relative errors (MREs), it was found that for small-sized Kacem problems, value adjusting with “DE/rand/1” and exponential crossover at position 2. Moreover, value adjusting with “DE/best/2” and exponential crossover at position 2 gave an MRE of 3.25. For medium-sized Brandimarte problems, value adjusting with “DE/best/2” and exponential crossover at position 2 gave a mean relative error of 7.11. For large-sized Dauzere-Peres and Paulli problems, value adjusting with “DE/best/2” and exponential crossover at position 2 gave an MRE of 4.20. From the comparison of the DE results with other methods, it was found that the MRE was lower than that found by Girish and Jawahar with the particle swarm optimization (PSO) method (7.75), which the improved DE was 7.11. For large-sized problems, it was found that the MRE was lower than that found by Warisa (1ST-DE) method (5.08), for which the improved DE was 4.20. The results further showed that basic DE and improved DE with jump search are effective methods compared to the other meta-heuristic methods. Hence, they can be used to solve the FJSP.

1. Introduction

Nowadays, the goal of businesses and industry is to reduce costs, and this is affected by production scheduling. Efficient production scheduling can reduce production expenses and time, resulting in on-schedule delivery of goods to customers and a competitive advantage for the firm. The issues of production scheduling concern the sequencing and machine assignment for each order. Owing to the requirement of the modern manufacturing processes for greater flexibility, the job shop scheduling problem (JSP) is an important type of production scheduling; furthermore, the flexible job shop scheduling problem (FJSP) was developed from the classical JSP. Job operations are allocated to every given machine in the production process. It is possible to assign any job to more than one machine as per the machine’s capability and, consequently, constructing an environment that is similar to the actual industry [1]. FJSP is an NP-Hard problem of the combinatorial optimization type, which has a complex solution. Applying a mathematical method to determine the exact algorithms for an optimal solution takes a moderate amount of time and wastes money; in particular, there is a big problem with the great number of variables and the limitations of the method. Metaheuristics are developed to solve FJSP to find a near optimal production schedule and to shorten the time required to solve the problem. It includes Tabu Search, the genetic algorithm (GA), particle swarm optimization (PSO), and ant colony optimization (ACO).
Differential evolution was proposed by Storn and Price [2]. It involves the optimization algorithm by using a population of each generation to search for the solution. Differential evolution calls the member of any generation vector and component of the vector point. Furthermore, it calls a number of points for each vector dimension. Each point can be compared to a gene from the GA method. The method of differential evolution is widely popular owing to its simplicity, various types of solution, and appropriate solutions. In some cases, DE results represent the global optimum.
Hence, in this research, we develop the differential evolution algorithm for solving in the flexible job shop scheduling problem with the target of minimizing the makespan.
This paper is structured as follows: In Section 2, the literature review is presented; Section 3 presents the mathematical model of the FJSP; Section 4 describes the general structure of the DE; Section 5 presents the results of the experiment on DE to solve the FJSP; Section 6 presents the comparison of DE with other metaheuristics; and Section 7 presents the conclusions and suggestions.

2. Literature Review

2.1. Flexible Job Shop Scheduling Problems by Using Other Metaheuristic Methods

To summarize the relevant literature and research on the solution to the flexible job shop scheduling problem by using metaheuristics, Xia and Wu [3] studied a hybrid of the PSO and SA methods for multiple purposes. Kanate [4] researched the development of a metaheuristic called the makespan tree for sequencing jobs on machines. Subsequently, two metaheuristics, the genetic algorithm and particle swarm optimization, were developed. Both metaheuristics use the makespan tree as a part of their method to solve the flexible job shop problems with the objective of minimizing the makespan. The findings for the job scheduling problems showed that the makespan tree outperformed the non-delay by 11.80%, improved the earliest finish time by 13.60%, and reduced the shortest processing time by 17.41%. In comparison, PSO had better results than GA by 0.97%. Tang et al. [5] researched the use of a hybrid algorithm for the flexible job-shop scheduling problem by combining chaos particle swarm optimization with the genetic algorithm in order to minimize the makespan. Wannaporn and Arit [6] applied the modified genetic algorithm to the flexible job-shop scheduling problem. This included the following processes: (1) Selecting chromosomes by the fuzzy roulette wheel selection method; (2) operating the crossover by the cluster-crossover operator to calculate the similarity between chromosomes for the crossover; (3) processing mutations using the mutation-local search operator to determine the diversity of the population, resulting in an optimal solution. The objective of this research was to minimize the makespan. Thanyaporn et al. [7] developed and improved the mixed integer programming (MIP) for an advanced planning and scheduling (APS) problem as a technique for production planning and scheduling. Its various considered constraints were the machine capacity, operation sequence, multi-machine due dates, multi-order, and product structure, comprising multiple steps and items. Each item could be processed by any given set of machines, and there was an extension to involve the considered constraints as a preventive maintenance time window. The objective was to minimize the total costs from idle production time, earliness, and tardiness, and furthermore, to optimize production scheduling. Examples of solutions were presented in four models which were (1) APS1, APS without alternative machines, (2) APS1PM, APS without alternative machines and PM, being similar to model 1 but with preventive maintenance included, (3) APS2, APS with alternative machines, and (4) APS2PM, APS with alternative machines and PM, being similar to model 3 but with preventive maintenance included. Hamid et al. [8] studied the machine scheduling in production—a content analysis. Which found that there were 132 surveyed papers regarding machine scheduling problems in production. The results were applied as an approach for proposing future research. It was found that, generally, the objective was to minimize the makespan. For the problem-solving approach, simple heuristics, as the shortest processing time first (SPT), and meta-heuristics, as the genetic algorithm, were employed. Li and Gao [9] studied an effective hybrid genetic algorithm and Tabu Search for flexible job shop scheduling problem. Luan et al. [10] studied improved whale algorithm for solving the flexible job shop scheduling problem and Li et al. [11] studied the hybrid artificial bee colony algorithm with a rescheduling strategy for solving the flexible job shop scheduling problems.

2.2. Flexible Job Shop Scheduling Problems by Using Differential Evolution Algorithm

Based on a study of the literature and research regarding the solutions to flexible job shop scheduling problems by using the differential evolution algorithm, Wisittipanich [12] proposed the application of adapted differential evolution algorithms to minimize the makespan in the FJSP. The modification of algorithms aimed to improve the efficiency of original DE by balancing its exploration and exploitation abilities to avoid the common problem of premature convergence. The first adapted DE was called the DE with subgroup strategy. In this algorithm, the population is divided into groups, and the population in each group employs different strategies to search for the new solution simultaneously in order to extract the strength of any strategy and compensate for the weakness of each strategy. This led to an increase in the overall performance efficiency. The second adapted DE was the DE with the switching strategy. This algorithm allowed the entire population to change searching strategies when there was no improvement to the solution. Thus, the chance of being trapped at a local optimum was decreased. The efficiency of two modified DEs was examined in solving an experimental problem and compared with the result of the original DE. The solutions provided by both modified DE algorithms were comparable or of a higher quality than the solutions from the original DE. Yuan and Xu [13] studied flexible job shop scheduling using hybrid differential evolution algorithms. Hybridization comprised the development of a mechanism to use the discrete differential evolution algorithm to solve the flexible job shop scheduling problem and second, enhancement of the local search ability in the DE framework. The objective was to find minimize the makespan. Bhaskara et al. [14] on the use of the differential evolution algorithm for the flexible job shop scheduling problem, applied the local search algorithm with the objective of minimizing the makespan.

2.3. Differential Evolution Algorithm for Solving Other Problems

The differential evolution algorithm is perceived as a modern method. It has become a favorable method to employ in various areas including operational research problems, such as application in the blocking flow shop scheduling problems to reduce production time [15]. Furthermore, in was applied in job shop production. Zhang et al. [16] adapted this method to deal with the job shop problem in order to minimize the total tardiness. For applications in the supply chain management problem, in particular, for the vehicle routing problem in the supply chain, the relevant studies have been carried out with different objectives and constraints under various conditions. For example, Cao and Lai [17] adapted the method to the open vehicle routing problem with fuzzy demands to reduce the total traveled distance. Lai and Cao [18] applied the method to solve the VRP with pickups and deliveries and time windows aimed at minimizing the total traveled distance. Xu and Wen [19] employed the DE in the unidirectional logistics distribution vehicle routing problem with no time windows and achieved the shortest total distance. By adapting DE with the agricultural management problem, Cruz et al. [20] applied the optimal control problem to determine the optimal control of nitrogen gas in lettuce. Moreover, DE was employed in a proposal on crop planning using the multi-objective differential evolution algorithm. This study had the objective of minimizing irrigation water usage and maximizing the total yield and net profit planting under various conditions [21]. DE was further applied to other problems, for example, in the electric power system, it was used to solve the optimal power flow, resulting in voltage stability enhancement and cost minimization [22]. In the wireless sensor network system, DE was proposed to prolong the lifetime of the system by preventing it from overloading [23]. An adaptation of DE, moreover, was used in the chemical industry in order to determine the optimal criteria for a chemical process [24].

3. Flexible Job Shop Scheduling Problem Pattern and Mathematical Model

We only focused on the FJSP in this study, and the details of the problem are as follows:

3.1. Flexible Job Shop Scheduling Problem (FJSP)

This production system is similar to the job shop scheduling problem but with more flexibility. As an explanation, any job includes a specific operation that can be processed by 1 or more machines owing to the various capabilities of the individual machine. As shown in Table 1, there are 3 jobs, each with a different set of operations. Every operation is allowed to select any machine from a given set, for example, in job 2, operation 1 can be performed on 2 machines, with M2 processing 10 time units and M3 processing 7 time units, while M1 shows “-”, referring to an inability to operate.
Figure 1 shows the general pattern of flexible job shop scheduling problem systems, which consist of a system with C work stations. In each work station, there will be a number of identical parallel machines. Each work station has its own specific route and can choose to perform the tasks assigned to one of the parallel machines and can be at the same work station. Considering the complexity of the model, the model of the flexible job shop scheduling problem systems that allows for recursive operation is the most redundant model.
To measure the production scheduling efficiency, there are various effective evaluations based on the production characteristics such as determining the minimization of the makespan (Cmax), the number of tardy jobs, and the maximum lateness. This paper employed the minimized makespan (Cmax) for evaluation in accordance with Equation (1). While C1, C2, …, Cp are possible solutions to produce scheduling, the solution resulting in the longest processing time was selected. Moreover, Z refers to the objective of production scheduling—to achieve the lowest value—as described in Equation (1). Figure 2 shows a sample of solutions from the total processing time of entire jobs:
Minimize Z = f (C1, C2, …, Cp).
The relative error is the difference between the measured and actual values, and it is generally expressed as a percentage (%). When the measured value is proximate to the actual value, this indicates high correctness or accuracy. The relative error is calculated using Equation (2), as follows;
% Relative   Error :   RE   =   C min BKS BKS × 100
where Cmin is the optimal solution to the algorithm, and BKS is the best known solution.

3.2. Mathematical Model of the Flexible Job Shop Scheduling Problem

The mathematical model of the FJSP, proposed by Kanate [4], includes many relevant binary variables, which are as follows:

3.2.1. Indices

i     Machine
j ,   k   Job
h ,   l   Operation

3.2.2. Parameter

M     Mathematically large real number
P i , j , h  Processing time for operation h of job j on machine i
O i , j , h  Operation h of job j on machine i
a i , j , h  Constant to define if job j at operation h is able to be processed by i —equal to 1 for the ability to process and 0 for the inability to process

3.2.3. Decision Variables

C m a x   Makespan
t j , h    Start time of the processing operation h of job j on any machine
f j , h    Finish time of the processing operation h of job j on any machine
y i , j , h    Binary variable equal to 1 while processing operation h of job j on machine i
X i , j , h , k , l  Binary decision variable equal to 1 while processing operation h of job j on machine ( O i , j , h ) , which comes before operation l of job k on machine i ( O i , k , l ) .
The mathematical model of the flexible job shop scheduling problem can be written as:
Minimize   C m a x   =   max   { C i } , i   =   1 ,   2 ,   3 ,   , n st .
t j , h + y i , j , h × P i , j , h f i , h ;           ( i , j , h )
f j , h t j , h + 1 ;           ( j , h )
f j , h C m a x ;           ( j , h )
y i , j , h a i , j , h ;         ( i , j , h )
t j , h + P i , j , h t k . l + ( 1 x i , j , h , k , l ) M ;           ( i , j , h , k , l )
i y i , j , h = 1 ;           ( h , j )
j h x i , j , h , k , l = y i , k , l ;           ( i , k , l )
k j x i , j , h , l , k = y i , j , h ;           ( i , j , h )
x i , j , h , j , h = 0 ;             ( i , j , h )
t j , h 0 ;         ( j , h )
f j , h 0 ;         ( j , h )
y i , j , h { 0 , 1 }
x i , j , h , k , l { 0 , 1 }
Equation (3) is a target equation to produce the minimum makespan, in order to reduce the production time, contributing to lower costs and product delivery time. Next, Equations (4) and (5) restrict the order based on the priority order in which each job is processed (precedence constraint). Equation (6) imposes a constraint on the makespan, whereby all job operations must be finished in a time less than or equal to the makespan. For Equation (7), only machines available for processing can be selected. Equation (8) is used to ensure that each machine can process, at most, one job at a time. Constrained by Equation (9), any job operation can be assigned to process on one machine only. Equations (10) and (11) create the sequence of job priority on machine i —Equation (10) selects the predecessor job and Equation (11) selects the next job. Equation (12) is the constraint that ensures that each job operation cannot be processed before its release time on machine i . Subsequently, Equations (13) and (14) restrict the start and finish times of every job to positive real numbers. Furthermore, Equations (15) and (16) specify y i , j , h and x i , j , h , k , l as binary variables.

4. General Differential Evolution Algorithm

The general differential evolution algorithm has general procedures [25], as shown in Figure 3.
The procedures of the differential evolution algorithm consist of (1) the initial population, and (2) mutation to generate mutant vector by differentiating a dimension of the vector. The calculation for mutating a vector’s dimensions is shown in Equation (17):
V m , n , G = X r 1 , n , G + F ( X r 2 , n , G X r 3 , n , G )
where m is the number of vectors of each generation, n is vector dimension, G is the iteration round, such as round 1, 2, or 3, and r1, r2, and r3 are 3 random vectors. F is the scaling factor, where V m , n , G is a mutant vector for vector m at vector dimension n in iteration round G. Moreover, X r 1 , n , G is firstly the random target vector at dimension n in iteration round G. X r 2 , n , G and X r 3 , n , G are, respectively, the second and third random target vectors. Thus, the mutant vector of vector m at dimension n in round G equals the value of target vector r1 at dimension n round G plus scaling factor F times the difference between target vectors r2 and r3. (3) Recombination of the trial vector is generated by exchanging vector dimensions. The equation employed for generating a trial vector is Equation (18):
U m , n , G = { V m , n , G   i f   r a n d m n C R   o r   D m = D m r a n d X m , n , G ,     e l s e .
where U m , n , G is the trial vector of vector m at dimension n in round G, r a n d m n is a random real number between [0, 1] of vector m at dimension n, CR refers to the crossover rate, Dm is the dimension of vector m. Furthermore, D m r a n d denotes a random integer number of vector m in the range [1, N], where N is the vector size. Consequently, the value of trial vector m in dimension n in iteration round G equals the mutant vector when there is a random number of vector m at dimension n that is less than CR; in other words, D m r a n d equals Dm. (4) The formula for selecting a target vector for the next round described in Equation (19). To solve the equation for the minimum target value (minimization), the less than or equal to sign is used:
X m , n , G + 1 = { U m , n , G   i f   f ( U m , n , G ) f ( X m , n , G ) X m , n , G ,     e l s e .
where, X m , n , G + 1 is the target vector m at dimension n in round G + 1, where the vector with a better fitness function value is selected in comparison with the value of the target vector and trial vector in round G.

4.1. Procedure of FJSP by Using Differential Evolution Algorithm

The procedure of FJSP uses the general differential evolution algorithm. The values used in the calculation must be set. The variables are as follows: Iterate = round, NP = number in population, F = scaling factor, and CR = crossover rate. In this problem’s calculation, these suitable variables were set from the experiment as Iterate = 1, NP = 5, F = 0.8, and CR = 0.8. The general differential evolution algorithm “DE/rand/1” and binomial crossover were used in the calculation.

4.1.1. Calculation Using the General Differential Evolution Algorithm DE/rand/1 and Binomial Crossover

The flexible job shop scheduling problem in Table 2 provides the details of the processing time on each machine. By explanation,
Row 1 indicates the numbers of jobs and machines, including 4 jobs and 5 machines.
Row 2 shows the data on job 1, comprising 3 operations. In operation 1, there are 5 available machines, which are machine 1 with a processing time of 2, machine 2 with a processing time of 5, machine 3 with a processing time of 4, machine 4 with a processing time of 1, and machine 5 with a processing time of 2. Operation 2 includes machine 1 with a processing time of 5, machine 2 with a processing time of 4, machine 3 with a processing time of 5, machine 4 with a processing time of 7, and machine 5 with a processing time of 5. Furthermore, operation 3 has machine 1 with a processing time of 4, machine 2 with a processing time of 5, machine 3 with a processing time of 5, machine 4 with a processing time of 4, and machine 5 with a processing time of 5.
Rows 3–5 show the data on jobs 2–4. The processing time of each operation can be likewise described, as shown in row 2, and put into categories, as shown in Table 3.
To improve the solution to round 1, the following procedure is completed:
Step 1: Generate the initial population
The initial population is randomized from a number between [0, 1], where the dimension or position (D) equals the number of operations (12), and the number of populations (P; 5) results in the target vector of a sample, as shown in Table 4.
Step 2: Differentiation of a dimension or mutation
Three vectors, r 1 , r 2 , and   r 3 , are randomly picked from the total population equal to 5 in order to generate a mutant vector of the entire 5 target vectors. As indicated by the sample problems, the randomized vectors ( r 1 , r 2 , r 3 ) from target vector 1 consist of vectors 2, 5, and 3. The other concerned target vectors are demonstrated in Table 5.
Calculation of the mutant vector ( V m , n , G ) can be performed by substituting the randomized vector into Equation (20). Thus, the mutant vectors of sample problems are presented in Table 6.
V m , n , G = X r 1 , n , G + F ( X r 2 , n , G X r 3 , n , G )
Step 3: Crossover
The number in the range [0, 1] at the target vector position is randomly picked and is comparative to the crossover rate (CR) predefined for the crossover. Table 7 shows the random numbers for the crossover.
Subsequently, the trial vector is calculated with Equation (21). In comparison, if a random number is less than or equal to CR = 0.8, the mutant vector in the same position will be selected as the obtained trial vector of that position. For differential cases, the target vector in the same position is the answer to the trial vector of the concerned position. Table 8 shows the obtained trial vectors.
U ji , G + 1   =   { V ji , G + 1   i f ( r a n d b ( j )     C R )   o r   j   =   r n b r   ( i )   X ji , G           i f ( r a n d b ( j )   >   C R )   o r   j     r n b r   ( i )  
Step 4: Fitness Evaluation
To decode the value of each dimension, the values are sorted based on the rank order value (ROV) in ascending order. Furthermore, the values are positioned in accordance with the vector order, as described in Table 9.
As a consequence of the sequencing machine operation, the machine with the earliest completion time is selected. If the machine operation consumes time equally, selection becomes random. This is described in Table 10 and illustrated by the Gantt chart in Figure 4.
Figure 4 presents a Gantt chart for solving the problem of a sample with 4 jobs, 5 machines, 12 operations, and a makespan of 17. The machines on the critical part are machines 1, 5, and 2 which consequentially interfere with solutions using the local search method.
Step 5: Selection
By comparing the fitness factor values of target vectors (Table 11) with the fitness factor values of the trial vectors (Table 12), the better vectors are selected as the target values for the next round. This is shown in Table 13.

4.1.2. Procedure of FJSP by Using the Improved Differential Evolution Algorithm

(1) The improved DE was developed by applying four mutation equations [28], DE/rand-to-best/1/bin, DE/rand/2/bin, DE/rand/1/exp 2 position, and DE/best/2/exp 2 position, as seen in Equations (22), (23), (24), and (25), as follows:
V m , n , G = X r 1 , n , G + F 1 ( X r 2 , n , G X r 3 , n , G ) + F 2 ( X b e s t , n , G X r 1 , n , G )
V m , n , G = X r 1 , n , G + F ( X r 2 , n , G X r 3 , n , G + X r 4 , n , G X r 5 , n , G )
V m , n , G = X r 1 , n , G + F ( X r 2 , n , G X r 3 , n , G + X r 4 , n , G X r 5 , n , G )
V m , n , G = X b e s t , n , G + F ( X r 2 , n , G X r 3 , n , G + X r 4 , n , G X r 5 , n , G )
Let r1, r2, r3, r4, and r5 denote the vectors which are randomly selected from a set of target vectors j, which represents the best vector found so far in the algorithm. F is a predefined integer parameter (scaling factor). In the proposed heuristics, F is set to 2; i is the vector number which ranges from 1 to NP, and j is the position of a vector which runs from 1 to D.
(2) The improved DE was developed by applying one crossover or recombination equation at exponential crossover position 2, as seen in Equation (26) [29], as follows:
U i , j , G = { V i , j , G   w h e n   j   r a n d i , 1   a n d   j     r a n d i , 2 X i , j , G   w h e n   r a n d i , 1   <   j   <   r a n d i , 2
As the predefined parameters in the proposed heuristics, let r a n d b i be a random number between 0 and 1 and let CR be the recombination probability. r a n d b i , r a n d b i , 1 , and r a n d b i , 2 are random integer numbers which represent the position of a vector; these random numbers range from 1 to D.
On the basis of the explanations in steps 1–4, the improved DE is shown in Algorithm 1.
Algorithm 1. Pseudo-code of the improved DE for the FJSP
Setup the initial DE parameter
Do while from first iteration to final iteration
Do while from first DE to final DE
  Setup the initial parameters: job, operation, machine, processing time,
  operation sequence, machine assignment.
  Do while from first task to final task
   Find the start/following task where the fitness is the makespan of the data instances
   Input the scaling factor, crossover rate, NP, job assignment, machine assignment, and
   local search to data list
   Produce the four mutation equations:
V m , n , G = X r 1 , n , G + F 1 ( X r 2 , n , G X r 3 , n , G ) + F 2 ( X best , n , G X r 1 , n , G )
V m , n , G = X r 1 , n , G + F ( X r 2 , n , G X r 3 , n , G + X r 4 , n , G X r 5 , n , G )
V m , n , G = X r 1 , n , G + F ( X r 2 , n , G X r 3 , n , G + X r 4 , n , G X r 5 , n , G )
V m , n , G = X best , n , G + F ( X r 2 , n , G X r 3 , n , G + X r 4 , n , G X r 5 , n , G )
   Developed by applying the two crossover or recombination equations:
U i , j , G = { V i , j , G   when   j   rand i , 1   and   j     rand i , 2 X i , j , G   when   rand i , 1   <   j   <   rand i , 2 .
   Produce the new target vector (selection/process):
X m , n , G + 1 = { U m , n , G   i f   f ( U m , n , G ) f ( X m , n , G ) X m , n , G , else .
  End do
End do
 Select the best solution from all DEs in the iteration
End do
Show/select the best solution from all DEs in all iterations

4.1.3. Procedure of FJSP by the Using Local Search with the Jump Search

The flexible job shop scheduling problem points to the optimal target of the minimum makespan. Considering the pathway of any latest complete operation, it denotes the critical pathway for flexible job shop scheduling. As shown in Figure 5 and Figure 6, the critical pathway is S → O2,1 → O4,1 → O4,2 → O4,3 → O2,2 → O2,3 → O2,4 → T.
Figure 6 reveal the critical pathway for flexible job shop scheduling. This pathway is able to determine solutions using the local search and jump search methods through an algorithm that identifies target values from the sorted critical pathway by thoroughly checking any operation on the critical pathway. In Figure 6, operation O2,1 is checked for possible intervention by a predecessor under these circumstances and priorities with a lower processing time and compatibility with an operating machine. Accordingly, operations are checked in the following consecutive order O4,1 → O4,2 → O4,3 → O2,2 → O2,3 → O2,4. Then, the fitness value is calculated in each round until completion of a set number of iterations.

5. Analysis of the Results from the Experiment on DE for Solving FJSP

To solve the flexible job shop scheduling problem, the Matlab program running on a personal computer (Core i5, 2.5 GHz, 8.00 GB RAM, Windows 7 operating system) was applied. It was developed by metaheuristic algorithms for the solutions focusing on the makespan.
Calculation factors were derived from the experiment based on optimization and relevant research. The findings were as follows: NP = 150, Iterate = 200, F = 2, CR = 0.8, number of iterated local search = 500 rounds. Moreover, the operation sequence that gave priority to operating the most remaining operations as well as machine assignment by choosing a machine with the minimum workload (MWL) were determined.

5.1. Results of Solving the Flexible Job Shop Scheduling Problem with Sample Problems from Kacem et al.

The results of solving the flexible job shop scheduling problem with sample problems from Kacem et al. [26,27] are shown in Table 14, and a Gantt chart of solutions to problem K01 is shown in Figure 7.
In Table 14, the result of solving the flexible job shop scheduling problem with small size problems [26,27] is revealed. The solutions optimizing differential evolution algorithms that provide the most optimal solutions are a combination of DE/rand/1 and exponential position 2 as well as a combination of DE/best/2 and exponential position 2, resulting in a makespan of 11. These solutions are illustrated in the Gantt chart in Figure 7 for problem K01. A further finding is that the lowest MRE is 3.25.

5.2. Results of Solving the Flexible Job Shop Scheduling Problem with Sample Problems of Brandimarte

The results of solving the flexible job shop scheduling problem with the sample problems of Brandimarte [30] can be found in Table 15 and the Gantt chart of solutions for problem Mk1 (Figure 8).
Table 15 shows the results of solving the flexible job shop scheduling problem for medium size problems using the sample problems of Brandimarte [30]. The solution optimizing the differential evolution algorithm that provides the most optimal solution is DE/best/2 combined with exponential position 2 and DE/rand/1 combined with exponential position 2, which results in a makespan of 40 and 204. This value is equal to the BKS value of the data set from sample Mk1 and Mk3. Furthermore, the lowest MRE value obtained is 7.11, as shown in the Gantt chart in Figure 8.

5.3. Results of Solving the Flexible Job Shop Scheduling Problem with Sample Problems of Dauzere-Peres and Paulli

The results of solving the flexible job shop scheduling problem with the sample problems of Dauzere-Peres and Paulli [31] can be found in Table 16.
Table 16 shows the results of solving the flexible job shop scheduling problem for large size problems using the sample problems of Dauzere-Peres and Paulli [31]. The solution optimizing the differential evolution algorithm that provides the best solution is DE/best/2 combined with exponential position 2, which results in a makespan of 2610. When compared with the other DE algorithm, it obtained the lowest value that is 4.20.

6. The Results of the Comparison of the DE Algorithm with Other Metaheuristic Methods

6.1. Results of Solving the Flexible Job Shop Scheduling Problem with Sample Problems of Brandimarte

A comparison of the differential evolution algorithm with other metaheuristic algorithms, GA and PSO, for the 10 comparative sample problems of Brandimarte [30] is given in Table 17.
Table 17 shows the comparative results for solving the flexible job shop scheduling problem using the differential evolution algorithm versus various dimension-optimizing algorithms. The value of mean of relative error (MRE) is lower than the value obtained by Girish and Jawahar [33] with PSO (7.75), while the improved DE has a value of 7.11. The work by Chen et al. [32] showed a value of 2.08 with the, while DE provided a greater value of 7.11. However, some problems, including Mk1 and Mk3, resulted in an equally good makespan in comparison with the BKS value of the data set.

6.2. Results of Solving the Flexible Job Shop Scheduling Problem with Sample Problems of Dauzere-Peres and Paulli

A comparison of the differential evolution algorithm with other metaheuristic algorithms, 1ST-DE, for the 5 comparative sample problems of Dauzere-Peres and Paulli [31] is given in Table 18.
Table 18 shows the comparative results for solving the flexible job shop scheduling problem using the differential evolution algorithm versus various dimension-optimizing algorithms. The value of mean of relative error (MRE) is lower than the value obtained by Wisittipanich [12] with 1ST-DE (5.08), while the improved DE has a value of 4.20. From the computational results, we can see that the improved DE algorithms with jump search are effective methods when compare with the basic DE and some meta-heuristic method.

7. Conclusions and Suggestions

This section presents the conclusions of this study. The differential evolution algorithm was used to solve the flexible job shop scheduling problem and optimize the dimensions. Among sample problems, the makespan and the BKS value of the data set were compared and the mean of relative error (MRE) was calculated. The sample problems of Kacem et al. were used as examples of small-sized problems. The dimensions were optimized by using “DE/rand/1” combined with exponential crossover position 2, as well as “DE/best/2” combined with exponential crossover position 2, which resulted in the minimum MRE value of 3.25. The sample problems of Brandimarte were used as examples of medium-sized problems. The dimensions were optimized with the combination of “DE/best/2” and exponential crossover position 2, providing a minimized MRE value of 7.11. Furthermore, the sample problems of Dauzere-Peres and Paulli were used as examples of large-sized problems. The dimensions were optimized with the combination of “DE/best/2” and exponential crossover position 2, providing a minimized MRE value of 4.20. Hence, the improved differential evolution in this research was able to solve the flexible job shop scheduling problem.
The DE algorithm proposed in this study can be applied to solve problems in the manufacturing industry in Thailand, such as mold and die manufacturing, the flexible job shop scheduling problem of the work center (FJSSPWC), or flexible job shop scheduling problem with preventive maintenance of the machine.

Author Contributions

P.S. designed the algorithm, summarized the computation, wrote the conclusions, and validated the algorithm; N.K. and P.K. gathered data and validated the algorithm.

Acknowledgments

We would like to thank the Industrial Engineering Department, Faculty of Engineering, Ubon Ratchathani University for funding this research.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Karthikeyan, S.; Asokan, P.; Nickolas, S. A hybrid discrete firefly algorithm for multi-objective flexible job shop scheduling problem with limited resource constraints. Int. J. Adv. Manuf. Technol. 2014, 72, 1567–1579. [Google Scholar] [CrossRef]
  2. Storn, R.; Price, K. Differential evolution—A simple and efficient heuristic for global optimization over continuous space. J. Glob. Optim. 1997, 11, 341–359. [Google Scholar] [CrossRef]
  3. Xia, W.; Wu, Z. An effective hybrid optimization approach for multi-objective flexible job-shop scheduling problems. Comput. Ind. Eng. 2005, 48, 409–425. [Google Scholar] [CrossRef]
  4. Kanate, P. Algorithm Development for Solving Flexible Job Shop Scheduling Problem. Ph.D. Thesis, Kasetsart University, Krung Thep Maha Nakhon, Thailand, 2011. [Google Scholar]
  5. Tang, J.; Zhang, G.; Lin, B.; Zhang, B. A Hybrid Algorithm for Flexible Job-shop Scheduling Problem. Procedia Eng. 2011, 15, 3678–3683. [Google Scholar] [CrossRef]
  6. Wannaporn, T.; Arit, T. Modified Genetic Algorithm for Flexible Job-Shop Scheduling Problems. Procedia Comput. Sci. 2012, 12, 122–128. [Google Scholar] [Green Version]
  7. Thanyapon, U.; Pupong, P.; Kwanniti, K. Solving an Advanced Planning and Scheduling Problem with Preventive Maintenance Time Window Constraints by Mixed Integer Programming Models. Thai J. Oper. Res. 2016, 4, 1–15. Available online: https://www.tci-thaijo.org/index.php/TJOR/article/view/60994 (accessed on 22 July 2019).
  8. Hamid, A.; Christoph, H.; Michael, D. Machine scheduling in production: A content analysis. Appl. Math. Model. 2017, 50, 279–299. [Google Scholar]
  9. Li, X.; Gao, L. An Effective Hybrid Genetic Algorithm and Tabu Search for Flexible Job Shop Scheduling Problem. Int. J. Prod. Econ. 2016, 174, 93–110. [Google Scholar] [CrossRef]
  10. Luan, F.; Cai, Z.; Wu, S.; Jiang, T.; Li, F.; Yang, J. Improved Whale Algorithm for Solving the Flexible Job Shop Scheduling Problem. Mathematics 2019, 7, 384. [Google Scholar] [CrossRef]
  11. Li, X.; Peng, Z.; Du, B.; Guo, J.; Xu, W.; Zhuang, K. Hybrid artificial bee colony algorithm with a rescheduling strategy for solving flexible job shop scheduling problems. Comput. Ind. Eng. 2017, 113, 10–26. [Google Scholar] [CrossRef]
  12. Wisittipanich, W. Minimizing Makespan in Flexible Job Shop Problems by Adapting the Differential Evolution. Thai J. Oper. Res. 2015, 3, 40–50. [Google Scholar]
  13. Yuan, Y.; Xu, H. Flexible job shop scheduling using hybrid differential evolution algorithms. Comput. Ind. Eng. 2013, 65, 246–260. [Google Scholar] [CrossRef]
  14. Bhaskara, P.; Padmanabhan, G.; Satheesh, K.B. Differential Evolution Algorithm for Flexible Job Shop Scheduling Problem. Int. J. Adv. Prod. Mech. Eng. 2015, 5, 71–79. [Google Scholar]
  15. Wang, L.; Pan, Q.K.; Tasgetiren, M.F. Minimizing the total flow time in a flow shop with blocking by using hybrid harmony search algorithms. Expert Syst. Appl. 2010, 37, 7929–7936. [Google Scholar] [CrossRef]
  16. Zhang, R.; Song, S.; Wu, C. A hybrid differential evolution algorithm for job shop scheduling problems with expected total tardiness criterion. Appl. Soft Comput. 2013, 13, 1448–1458. [Google Scholar] [CrossRef]
  17. Cao, E.; Lai, M. The open vehicle routing problem with fuzzy demands. Expert Syst. Appl. 2010, 37, 15. [Google Scholar] [CrossRef]
  18. Lai, M.; Cao, E. An improved differential evolution algorithm for vehicle routing problem with simultaneous pickups and deliveries and time windows. Eng. Appl. Artif. Intell. 2010, 23, 188–195. [Google Scholar]
  19. Xu, H.; Wen, J. Differential Evolution Algorithm for the Optimization of the Vehicle Routing Problem in Logistics. In Proceedings of the 2012 Eighth International Conference on Computational Intelligence and Security, Guangzhou, China, 17–18 November 2012; pp. 48–51. [Google Scholar]
  20. Cruz, I.L.; van Willigenburg, L.G.; van Straten, G. Efficient differential evolution algorithms for multimodal optimal control problems. Appl. Soft Comput. 2003, 3, 97–122. [Google Scholar] [CrossRef]
  21. Adeyemo, J.; Otieno, F. Differential evolution algorithm for solving multi-objective crop planning model. Agric. Water Manag. 2010, 97, 848–856. [Google Scholar] [CrossRef]
  22. Abou El Ela, A.A.; Abido, M.A.; Spea, S.R. Optimal power flow using differential Evolution algorithm. Electr. Power Syst. Res. 2010, 80, 878–885. [Google Scholar] [CrossRef]
  23. Kuila, P.; Jana, P.K. A novel differential evolution based clustering Algorithm for wireless sensor networks. Appl. Soft Comput. 2014, 25, 414–425. [Google Scholar] [CrossRef]
  24. Sharma, S.; Rangaiah, G.P. An improved multi-objective differential evolution with a termination criterion for optimizing chemical processes. Comput. Chem. Eng. 2013, 56, 155–173. [Google Scholar] [CrossRef]
  25. Pitakaso, R.; Parawech, P.; Jirasirierd, G. Comparisons of Different Mutation and Recombination Processes of the DEA for SALB-1. In Proceedings of the Institute of Industrial Engineers Asian Conference, Taipei, Taiwan, 18–20 July 2013; pp. 1571–1579. [Google Scholar]
  26. Kacem, I.; Hammadi, S.; Borne, P. Approach by localization and multiobjective evolutionary optimization for flexible job-shop scheduling problems. IEEE Trans. Syst. Man Cybern. 2002, 32, 1–13. [Google Scholar] [CrossRef]
  27. Kacem, I.; Hammadi, S.; Borne, P. Pareto-optimality approach for flexible job-shop scheduling problems: hybridization of evolutionary algorithms and fuzzy logic. Math. Comput. Simul. 2002, 60, 245–276. [Google Scholar] [CrossRef]
  28. Nguyen, S.; Kitchitvichynukul, V.; Wisittipanich, W. ET-Lib User’s Guide. Volume 2. Differential Evolution; Asian Institute of Technology: Khlong Luang, Thailand, 2013. [Google Scholar]
  29. Poontana, S.; Nuchsara, K.; Preecha, K.; Krit, C. U-Shaped Assembly Line Balancing by Using Differential Evolution Algorithm. Math. Comput. Appl. 2018, 23, 79. [Google Scholar] [Green Version]
  30. Brandimarte, P. Routing and Scheduling in a Flexible Job Shop by Tabu Search. Ann. Oper. Res. 1993, 41, 157–183. [Google Scholar] [CrossRef]
  31. Dauzere-Peres, S.; Paulli, J. An integrated approach for modeling and solving the general multiprocessor job-shop scheduling problem using tabu search. Ann. Oper. Res. 1997, 70, 281–306. [Google Scholar] [CrossRef]
  32. Chen, H.; Ihlow, J.; Lehmann, C. A genetic algorithm for flexible job-shop scheduling. In Proceedings of the IEEE International Conference on Robotics and Automation, Detroit, MI, USA, 10–15 May 1999; pp. 1120–1125. [Google Scholar]
  33. Girish, B.S.; Jawahar, N. A particle swarm optimization algorithm for flexible job shop scheduling problem. In Proceedings of the IEEE International Conference on Automation Science and Engineering, Bangalore, India, 22–25 August 2009; pp. 298–303. [Google Scholar]
Figure 1. Flexible job shop scheduling problem pattern.
Figure 1. Flexible job shop scheduling problem pattern.
Mca 24 00080 g001
Figure 2. Gantt chart of solutions.
Figure 2. Gantt chart of solutions.
Mca 24 00080 g002
Figure 3. General differential evolution algorithm procedure.
Figure 3. General differential evolution algorithm procedure.
Mca 24 00080 g003
Figure 4. Gantt chart of sample problems.
Figure 4. Gantt chart of sample problems.
Mca 24 00080 g004
Figure 5. Disjunction graph of flexible job shop scheduling.
Figure 5. Disjunction graph of flexible job shop scheduling.
Mca 24 00080 g005
Figure 6. Critical pathway for flexible job shop scheduling.
Figure 6. Critical pathway for flexible job shop scheduling.
Mca 24 00080 g006
Figure 7. Gantt chart of solutions to flexible job shop scheduling problem K01 presented by Kacem et al. [26,27].
Figure 7. Gantt chart of solutions to flexible job shop scheduling problem K01 presented by Kacem et al. [26,27].
Mca 24 00080 g007
Figure 8. Gantt chart of solutions for flexible job shop scheduling problem Mk1 presented by Brandimarte [30].
Figure 8. Gantt chart of solutions for flexible job shop scheduling problem Mk1 presented by Brandimarte [30].
Mca 24 00080 g008
Table 1. Example of flexible job shop scheduling problem pattern.
Table 1. Example of flexible job shop scheduling problem pattern.
JobOperationMachines
M1M2M3
J1O1,15-3
O1,2-510
O1,359-
J2O2,1-107
O2,2206-
O3,32-11
J3O3,1254
O3,32510
Table 2. Sample problems of Kacem et al. [26,27].
Table 2. Sample problems of Kacem et al. [26,27].
45.
35. 1 2 2 5 3 4 4 1 5 2 5 1 5 2 4 3 5 4 7 5 5 5 1 4 2 5 3 5 4 4 5 5
35. 1 2 2 5 3 4 4 7 5 8 5 1 5 2 6 3 9 4 8 5 5 5 1 4 2 5 3 4 4 54 5 5
45. 1 9 2 8 3 6 4 7 5 9 5 1 6 2 1 3 2 4 5 5 4 5 1 2 2 5 3 4 4 2 5 4 5 1 4 2 5 3 2 4 1 5 5
25. 1 1 2 5 3 2 4 4 5 12 5 1 5 2 1 3 2 4 1 5 2
Table 3. Sample problem with four jobs and five machines.
Table 3. Sample problem with four jobs and five machines.
JobsOperationsMachines
M1M2M3M4M5
J1O1,125412
O1,254575
O1,345545
J2O2,125478
O2,256985
O2,3454545
J3O3,198679
O3,261254
O3,325424
O3,445215
J4O4,1152412
O4,251212
Table 4. Target vectors of sample problems.
Table 4. Target vectors of sample problems.
NPDimensions, D
123456789101112
10.550.320.700.120.640.890.960.810.380.550.270.71
20.170.800.940.930.440.360.770.350.130.420.170.11
30.420.350.150.610.100.340.930.510.080.590.630.50
40.650.720.300.580.020.740.590.170.140.070.730.31
50.720.320.040.200.890.280.420.670.150.490.090.81
Table 5. Randomized vectors r 1 , r 2 , and   r 3 .
Table 5. Randomized vectors r 1 , r 2 , and   r 3 .
Random Vectorr1r2r3
1253
2335
3422
4511
5144
Table 6. Mutant vectors of sample problems in round 1.
Table 6. Mutant vectors of sample problems in round 1.
Mutation123456789101112
10.770.740.420.112.020.24−0.250.670.270.22−0.910.73
20.800.130.740.241.780.50−0.470.041.810.38−0.430.37
30.160.700.560.611.450.43−0.740.241.30−0.080.020.40
40.810.430.24−0.360.030.090.970.69−0.701.041.191.04
50.871.300.650.691.051.390.23−0.861.310.600.310.91
Table 7. Random numbers for the crossover.
Table 7. Random numbers for the crossover.
Vector123456789101112
10.830.050.750.800.950.390.290.600.300.440.590.65
20.680.360.480.470.700.960.040.760.640.420.160.44
30.320.400.970.380.630.690.710.920.650.830.920.49
40.560.180.060.380.470.230.110.850.800.300.650.02
50.810.350.700.500.890.890.840.290.010.210.410.83
Table 8. Trial vectors from round 1.
Table 8. Trial vectors from round 1.
Trial Vector123456789101112
10.550.740.420.120.640.24−0.250.670.270.22−0.910.73
20.800.130.740.241.780.50−0.470.041.810.38−0.430.37
30.160.700.560.611.450.43−0.740.241.30−0.080.020.40
40.650.720.300.580.020.740.590.170.140.070.730.31
50.720.320.040.200.890.280.420.670.150.490.090.81
Table 9. Results of decoding.
Table 9. Results of decoding.
Vector117410693158122
1−0.91−0.250.120.220.240.270.420.550.640.670.730.74
Oi,j1, 11, 21, 32, 12, 22, 33, 13, 23, 33, 44, 14, 2
Table 10. Results of the sequencing machine operation.
Table 10. Results of the sequencing machine operation.
Vector123456789101112
10.550.740.420.120.640.24−0.250.670.270.22−0.910.73
Oi,j3,24,23,11,33,32,21,23,42,32,11,14,1
M224115243141
PT116425414211
Table 11. Target vectors of round 1.
Table 11. Target vectors of round 1.
Target Vector123456789101112Target
10.550.320.700.120.640.890.960.810.380.550.270.7121
20.170.800.940.930.440.360.770.350.130.420.170.1120
30.420.350.150.610.100.340.930.510.080.590.630.5022
40.650.720.300.580.020.740.590.170.140.070.730.3124
50.720.320.040.200.890.280.420.670.150.490.090.8119
Table 12. Trial vectors of round 1.
Table 12. Trial vectors of round 1.
Trial Vector123456789101112Target
10.550.740.420.120.640.24−0.250.670.270.22−0.910.7318
20.800.130.740.241.780.50−0.470.041.810.38−0.430.3725
30.160.700.560.611.450.43−0.740.241.30−0.080.020.4016
40.650.720.300.580.020.740.590.170.140.070.730.3117
50.720.320.040.200.890.280.420.670.150.490.090.8126
Table 13. Target vectors for the next round.
Table 13. Target vectors for the next round.
Vector123456789101112Target
10.550.740.420.120.640.24−0.250.670.270.22−0.910.7318
20.170.800.940.930.440.360.770.350.130.420.170.1120
30.160.700.560.611.450.43−0.740.241.30−0.080.020.4016
40.650.720.300.580.020.740.590.170.140.070.730.3117
50.720.320.040.200.890.280.420.670.150.490.090.8119
Table 14. Summary of solving the flexible job shop scheduling problem with sample problems from Kacem et al. [26,27].
Table 14. Summary of solving the flexible job shop scheduling problem with sample problems from Kacem et al. [26,27].
ProblemBKSMutation Strategy
DE *DE **DE ***DE ****
K011112
(9.09)
12
(9.09)
11
(0.00)
11
(0.00)
K021415
(7.14)
15
(7.14)
15
(7.14)
15
(7.14)
K031111
(0.00)
11
(0.00)
11
(0.00)
11
(0.00)
K0477
(0.00)
7
(0.00)
7
(0.00)
7
(0.00)
K051112
(9.09)
12
(9.09)
12
(9.09)
12
(9.09)
MRE5.065.063.253.25
* DE/rand to best/1/Bin; ** DE/rand/2/Bin; *** DE/rand/1/Exp Crossover Position 2/Local Search with Jump Search; **** DE/best/2/Exp Crossover Position 2/Local Search with Jump Search; Mean relative error (MRE); Best known solution (BKS).
Table 15. Summary of solving the flexible job shop scheduling problem with sample problems Brandimarte [30].
Table 15. Summary of solving the flexible job shop scheduling problem with sample problems Brandimarte [30].
ProblemBKSMutation Strategy
DE *DE **DE ***DE ****
Mk14043
(7.50)
43
(7.50)
40
(0.00)
40
(0.00)
Mk22728
(7.69)
28
(7.69)
28
(7.69)
28
(7.69)
Mk3204204
(0.00)
204
(0.00)
204
(0.00)
204
(0.00)
Mk46071
(18.33)
71
(18.33)
71
(18.33)
71
(18.33)
Mk5174178
(2.30)
178
(2.30)
179
(2.87)
179
(2.87)
Mk65973
(23.73)
73
(23.73)
73
(23.73)
73
(23.73)
Mk7143149
(4.20)
149
(4.20)
148
(3.50)
146
(2.10)
Mk8523528
(0.96)
528
(0.96)
528
(0.96)
528
(0.96)
Mk9307324
(5.54)
321
(4.56)
323
(5.21)
321
(4.56)
Mk10212234
(10.38)
233
(9.90)
236
(11.32)
235
(10.85)
MRE8.067.927.367.11
* DE/rand to best/1/Bin; ** DE/rand/2/Bin; *** DE/rand/1/Exp Crossover Position 2/Local Search with Jump Search; **** DE/best/2/Exp Crossover Position 2/Local Search with Jump Search; Mean relative error (MRE).
Table 16. Summary of solving the flexible job shop scheduling problem with sample problems Dauzere-Peres and Paulli.
Table 16. Summary of solving the flexible job shop scheduling problem with sample problems Dauzere-Peres and Paulli.
ProblemBKSMutation Strategy
DE *DE **DE ***DE ****
01a25302895
(14.42)
2750
(8.70)
2615
(3.36)
2645
(4.55)
04a25552859
(11.90)
2770
(8.41)
2650
(3.72)
2610
(2.15)
07a23962759
(15.15)
2650
(10.60)
2650
(10.60)
2510
(4.76)
09a20742281
(9.98)
2269
(9.40)
2210
(6.56)
2150
(3.66)
11a20782378
(14.44)
2366
(13.86)
2221
(6.88)
2200
(5.87)
MRE13.1810.196.224.20
* DE/rand to best/1/Bin; ** DE/rand/2/Bin; *** DE/rand/1/Exp Crossover Position 2/Local Search with Jump Search; **** DE/best/2/Exp Crossover Position 2/Local Search with Jump Search; Mean relative error (MRE).
Table 17. Summary of comparing the differential evolution algorithm with other metaheuristic algorithms.
Table 17. Summary of comparing the differential evolution algorithm with other metaheuristic algorithms.
Problemn × m × k *BKS **Chen et al. (GA) [32]Girish and Jawahar (PSO) [33]DE-FJSP
CmaxCmaxCmaxCmax
Mk0110 × 6 × 554040
(0.00)
40
(0.00)
40
(0.00)
Mk0210 × 6 × 582729
(6.89)
27
(0.00)
28
(7.69)
Mk0315 × 8 × 150204204
(0.00)
204
(0.00)
204
(0.00)
Mk0415 × 8 × 906063
(4.76)
62
(3.22)
71
(18.33)
Mk0515 × 4 × 106174181
(3.86)
178
(2.24)
179
(2.87)
Mk0610 × 15 × 1505960
(1.66)
78
(24.35)
73
(23.73)
Mk0720 × 5 × 100143148
(3.38)
147
(2.72)
146
(2.10)
Mk0820 × 10 × 225523523
(0.00)
523
(0.00)
528
(0.96)
Mk0920 × 10 × 240307308
(0.32)
341
(9.97)
321
(4.56)
Mk1020 × 15 × 240212212
(0.00)
252
(15.07)
235
(10.85)
MRE2.087.757.11
* n = Job, m = Machine, k = Operation; ** Best known solution (BKS) ; FJSP = Flexible job shop scheduling problem; GA = Genetic algorithm; PSO = Particle swarm optimization.
Table 18. Summary of comparing the differential evolution algorithm with other metaheuristic algorithms.
Table 18. Summary of comparing the differential evolution algorithm with other metaheuristic algorithms.
Problemn × m × k *BKS **Wisittipanich (1ST-DE) [12]DE-FJSP
CmaxCmaxCmax
01a10 × 5 × 19625302645
(4.55)
2645
(4.55)
04a10 × 5 × 19625552616
(2.39)
2610
(2.15)
07a15 × 8 × 29323962582
(7.76)
2510
(4.76)
09a15 × 8 × 29320742153
(3.81)
2150
(3.66)
11a15 × 8 × 29320782221
(6.88)
2200
(5.87)
MRE5.084.20
* n = Job, m = Machine, k = Operation; ** Best known solution (BKS); FJSP = Flexible job shop scheduling problem.

Share and Cite

MDPI and ACS Style

Sriboonchandr, P.; Kriengkorakot, N.; Kriengkorakot, P. Improved Differential Evolution Algorithm for Flexible Job Shop Scheduling Problems. Math. Comput. Appl. 2019, 24, 80. https://doi.org/10.3390/mca24030080

AMA Style

Sriboonchandr P, Kriengkorakot N, Kriengkorakot P. Improved Differential Evolution Algorithm for Flexible Job Shop Scheduling Problems. Mathematical and Computational Applications. 2019; 24(3):80. https://doi.org/10.3390/mca24030080

Chicago/Turabian Style

Sriboonchandr, Prasert, Nuchsara Kriengkorakot, and Preecha Kriengkorakot. 2019. "Improved Differential Evolution Algorithm for Flexible Job Shop Scheduling Problems" Mathematical and Computational Applications 24, no. 3: 80. https://doi.org/10.3390/mca24030080

APA Style

Sriboonchandr, P., Kriengkorakot, N., & Kriengkorakot, P. (2019). Improved Differential Evolution Algorithm for Flexible Job Shop Scheduling Problems. Mathematical and Computational Applications, 24(3), 80. https://doi.org/10.3390/mca24030080

Article Metrics

Back to TopTop