Next Article in Journal
Reliability Assessment of RC Bridges Subjected to Seismic Loadings
Next Article in Special Issue
Prediction of Machine Failure in Industry 4.0: A Hybrid CNN-LSTM Framework
Previous Article in Journal
A Transfer Learning Technique for Inland Chlorophyll-a Concentration Estimation Using Sentinel-3 Imagery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hybrid Salp Swarm Algorithm for Solving the Green Scheduling Problem in a Double-Flexible Job Shop

1
Faculty of Management Engineering, Huaiyin Institute of Technology, Huaian 223200, China
2
College of Economics and Management, Suzhou Chienshiung Institute of Technology, Taicang 215411, China
3
Faculty of Mathematics and Physics, Huaiyin Institute of Technology, Huaian 223003, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(1), 205; https://doi.org/10.3390/app12010205
Submission received: 1 October 2021 / Revised: 10 December 2021 / Accepted: 20 December 2021 / Published: 25 December 2021
(This article belongs to the Special Issue Artificial Intelligence and Optimization in Industry 4.0)

Abstract

:
Green scheduling is not only an effective way to achieve green manufacturing but also an effective way for modern manufacturing enterprises to achieve energy conservation and emission reduction. The double-flexible job shop scheduling problem (DFJSP) considers both machine flexibility and worker flexibility, so it is more suitable for practical production. First, a multi-objective mixed-integer programming model for the DFJSP with the objectives of optimizing the makespan, total worker costs, and total influence of the green production indicators is formulated. Considering the characteristics of the problem, three-layer salp individual encoding and decoding methods are designed for the multi-objective hybrid salp swarm algorithm (MHSSA), which is hybridized with the Lévy flight, the random probability crossover operator, and the mutation operator. In addition, the influence of the parameter setting on the MHSSA in solving the DFJSP is investigated by means of the Taguchi method of design of experiments. The simulation results for benchmark instances show that the MHSSA can effectively solve the proposed problem and is significantly better than the MSSA and the MOPSO algorithm in the diversity, convergence, and dominance of the Pareto frontier.

1. Introduction

Traditional job shop scheduling problems usually do not consider human factors. With increasing research on job shop scheduling problems, researchers believe that human factors should be considered in production scheduling and operation management activities [1,2,3,4]. Based on the flexible job shop scheduling problem (FJSP), Gong et al. [5] proposed a double-flexible job shop scheduling problem (DFJSP) for the first time, which is more in line with the actual production situation. That is, each process of each workpiece can be processed by multiple machines, and each machine can be operated by multiple workers. The DFJSP is an NP hard combinatorial optimization problem that is an extension of FJSP. At present, there is no polynomial time algorithm to solve the problem [5].
In recent years, with the increasingly severe environmental problems, research on green manufacturing has attracted a lot of attention. For example, He and Liu [6] studied the issue of reducing energy consumption for the job-shop scheduling problem, constructed the biobjective model of minimizing both the energy consumption and the makespan, and used the Tabu search algorithm to optimize the model. Mouzon et al. [7] proposed a multi-objective mathematical programming model of minimizing the energy consumption and total completion time. By solving the model, we can get a set of non-dominated solutions to determine the most efficient production sequence. The model only considers the energy saving by shutting down non-bottleneck or idle equipment, and it does not consider the energy saving of running machines and equipment. Dai and Tang et al. [8] modeled an energy-efficient scheduling with makespan and energy consumption. The results showed that the relationship between the two objectives was apparently one of conflict. Similar results were also proved in the study of Liu and Huang [9]. Mansouri et al. [10] developed a mixed integer linear multi-objective optimization model to find the Pareto frontier comprised of makespan and total energy consumption, and they defined lower bounds for the objectives and proposed a heuristic for the problem. In addition, some scholars have also carried out similar studies related to green production (see references [11,12,13,14,15], etc.). The above research results mainly focus on the theoretical model, which is quite different from the actual green manufacturing and production scheduling. Green production indicators need to be used to build scheduling models closer to the actual situation.
On the other hand, green scheduling is an important part of green manufacturing, and some scholars have studied it. For example, in reference [16], green job shop scheduling is defined as efficiency increasing, energy saving, emission reduction, and consumption reduction, which improve economic benefits and lead to a green manufacturing process through resource allocation, operation sequencing, and reasonable optimization of the operation mode. Wang et al. [17] analyzed the complexity of the green job shop scheduling problem and the corresponding processing mechanism and summarized and reviewed the representative research on green job shop scheduling from the problem, method, and application levels. To solve the flexible job shop scheduling problem, Lei et al. [18] proposed a hybrid frog-leaping algorithm that considers the energy consumption index. Wu et al. [19] took completion time and energy consumption as the optimization objectives, studied the flexible job shop problem with multi-speed machines, analyzed the energy consumption distribution of machines, established the energy consumption calculation model, and solved it by the NSGA-II algorithm. Dong and Ye [20] constructed a scheduling model for green manufacturing collaborative optimization of the semiconductor wafer distributed heterogeneous factory and proposed an improved gray wolf algorithms to solve the model. To solve the process scheduling problem, Yao et al. [21] designed a multi-objective salp swarm algorithm to solve the process scheduling problem of a TFT-LCD panel array that considers energy saving. Nanthapodej et al. [22] proposed a variable neighborhood strategy adaptive search method to minimize energy consumption while also considering job priority and makespan control for parallel-machine scheduling problems. In addition, some researchers have also carried out similar studies related to green scheduling (see references [23,24,25], etc.). To sum up, the current research on green scheduling mostly focuses on optimizing green production indicators, such as energy consumption and carbon emissions, while there are few studies considering human factors in green scheduling; Moreover, there are still many indicators related to resource consumption and environmental impact, such as emissions of toxic substances and wastes, noise, workpiece chip recycling, etc. The research on green scheduling considering these factors is also rare. At present, there are still many tasks on green workshop scheduling to deal with, such as constructing new models closer to practical production scheduling using green production indicators and developing more efficient algorithms to solve the green scheduling problem.
In this paper, a double-flexible job shop green scheduling problem model is constructed with the objectives of minimizing the maximum completion time, the total labor cost, and the total green index. A multi-objective hybrid salp swarm algorithm is proposed to solve the complex job shop scheduling problem, such as the DFJSP. By testing a benchmark example and comparing it with two typical algorithms, the feasibility and effectiveness of the algorithm designed in this paper in solving the double flexible job shop green scheduling problem are verified.

2. Green Scheduling Problem Model of Double Flexible Job Shop

2.1. Problem Description

The DFJSP can be described as follows: n workpieces are processed on m parallel machines of different types, which can be operated by l workers with different skill levels. The workpiece i includes r i processes, and the processing sequence of each process is given in advance. The j -th process O i j of the workpiece i can be processed on any machine in the set of optional processing machines. At the same time, each machine can be operated by any worker in the set of optional workers. The solution of the DFJSP includes three sub-problems: machine selection, worker selection, and process scheduling. Its scheduling goal is to select the most suitable machine and matching operators for each process and determine the optimal processing sequence and start-up time of each process on each machine so that some performance indexes of the production system can be optimized.
The double-flexible job shop green scheduling problem constructed in this paper has three optimization objectives: minimizing the maximum completion time, total labor cost, and total green index. Among them, the maximum completion time represents the final completion time of all jobs and all processes, which reflects the production efficiency. The total labor cost is an important part of the production cost of an enterprise. By reducing the total labor cost, the operating cost of an enterprise can be reduced. The total green index in this paper includes energy consumption, noise, workpiece chip recycling, and safety, making it a comprehensive green production index. Minimizing the total green index can help enterprises save energy, reduce emission, and reduce the impact on the environment.
The assumptions and constraints to be met are as follows: (1) Each machine can handle only one process at a time. (2) Each process can be processed only once on one machine. (3) Each worker can operate only one machine at a time. (4) Each process can be processed only once by one worker. (5) There are sequential constraints between the same workpiece processes, and there are no sequential constraints between different workpiece processes. (6) Interruption is not allowed after a process starts. (7) All machines and workpieces are ready at zero hour. (8) The installation and adjustment time of the workpiece is not considered. According to the skill levels, workers are divided into three types. In addition, the processing time, the labor cost, and the energy consumption corresponding to each process, as well as the attribute values of various green indexes of the machine, are known and fixed.

2.2. Mathematical Model of the DFJSP

The mathematical symbols in the model and their meanings are as follows:
m : Number of machines;
n : Number of work pieces;
l : Number of workers;
i , h : Workpiece serial number; i , h = 1 , 2 , , n ;
r i : Total number of operations of workpiece i ;
j , g : Operation sequence number, j , g = 1 , 2 , , r ;
k : Machine serial number, k = 1 , 2 , , m ;
s : Worker serial number, s = 1 , 2 , , l ;
O i j : The j -th process of workpiece i ;
C i j : Completion time of process O i j ;
C k : Completion time of machine k ;
P i j k s : Processing time of the process O i j operated by worker s on machine k ;
L i j k s : Processing cost of the process O i j operated by worker s on machine k ;
E I i j k s : Green index value of the process O i j operated by worker s on machine k ;
E i j k s : Energy consumption of the process O i j operated by worker s on machine k ;
N i j k : Noise generated by the process O i j on machine k ;
R i j k : Chip recovery rate of the process O i j processed on machine k ;
S i j k : Safety of the process O i j processing on machine k ;
x i j k = { 1 , the process   O i j   is   processed   on   machine   k 0 , other
x i j k s = { 1 , the   process   O i j   is   processed   on   machine   k   and   operated   by   the   worker   s 0 , other
Based on the research in [5], this paper constructs the following model where C M represents the maximum completion time target, L C represents the total labor cost, and G represents the total green index.
Target function:
M i n C M = max 1 k m ( C k )
M i n L C = i = 1 n j = 1 r i k = 1 m s = 1 l ( L i j k s × x i j k s )
M i n G = i = 1 n j = 1 r i k = 1 m s = 1 l ( E I i j k s × x i j k s )
E I i j k s = E i j k s + N i j k + ( 1 / R i j k ) + ( 1 / S i j k )
Constraints:
C i j C i ( j 1 ) P i j k s x i j k s , i , j ; i = 1 , 2 , , n ; j = 2 , 3 , , r i
k = 1 m s = 1 l x i j k s = 1 , i , j ; i = 1 , 2 , , n ; j = 1 , 2 , , r i
[ ( C h g C i j t h g k ) x h g k x i j k 0 ] [ ( C i j C h g t i j k ) x h g k x i j k 0 ] , i , j , h , g , k ; i , h = 1 , 2 , , n ; j = 1 , 2 , , r i ; g = 1 , 2 , , r h ; k = 1 , 2 , , m ;
Here, Equations (1)–(3) represent the objective functions of maximum completion time, total labor cost, and total green index, respectively. Equation (4) defines the green index value E I i j k s of the process O i j on machine k operated by worker s . The value of this index is determined by four factors: energy consumption, noise, chip recovery, and safety. The lower the energy consumption and noise value, the better and the higher the chip recovery rate and safety index value, so the reciprocal is taken for calculation. In addition, because the measurement units and meanings of the above four green indicators are different, they cannot be calculated directly. Therefore, it is necessary to normalize the four indicators by the min–max standardization method. The calculation formula of the min–max standardization method is x = x x m i n x m a x x m i n , where x represents the new value after x normalization, x m a x is the maximum value of x indicator, and x m i n is the minimum value of x indicator. Equation (5) is a process constraint to ensure that the workpiece is processed according to the specified process. Equation (6) ensures that each process can be assigned only to a machine in the set of optional processing machines for the process and the machine can be operated only by a worker in the set of optional workers. Equation (7) ensures that each machine can handle only one working procedure at a time.

3. DFJSP Solution Based on the Multi-Objective Hybrid Salp Swarm Algorithm

3.1. Salp Swarm Algorithm

The salp swarm algorithm (SSA) is a population intelligent optimization algorithm proposed by scholars Mirjalili and Gandomi in 2017 [26]. This stochastic optimization algorithm simulates the movement and foraging behavior of the salp swarm group in the ocean, and its effectiveness is verified by solving function optimization problems and engineering design problems. The salp swarm algorithm divides the salp chain into two parts: the first half of the individuals as leaders and the second half of the individuals as followers. It is assumed that there is a food source F in the search space, which is the search target of the group. The position update of the leader is only related to the position of the food source, while the position update of the follower is only related to the position of the salp swimming in front of it. The multi-objective salp swarm algorithm (MSSA) is a multi-objective swarm intelligent optimization algorithm based on the SSA, using the same population division and individual location update methods and referring to the external archive access mechanism of the MOPSO algorithm [27]. The basic idea is as follows: N salp individuals are randomly generated as the initial population; then, the objective function value of these salp individuals is calculated, and all non-dominated individuals in the population are identified according to the Pareto dominance relationship. If the number of individuals in the current external archive does not exceed the specified number, the non-dominated solution is directly saved to the archive. Otherwise, the redundant individuals are deleted from the archive-dense area before a new non-dominated solution is added. When performing the deletion operation, each solution in the archive needs to be sorted and selected by roulette. After updating the archive, according to the sorting value, an individual is selected as the food source F from the sparse area of the archive by the roulette method. Next, the value of c is updated, and the individual position in the population of salp is updated according to the food source F. The above steps are repeated until the termination conditions are met. The key operations involved in the MSSA are the update and maintenance of the archive and the selection of the food source F.
According to the above analysis, the MSSA is easy to implement, because it only has two main control parameters: c and archive size. In addition, the non-dominated solutions obtained so far are stored in the archive, so even if the whole population degenerates, the non-dominated solutions found will not be lost. Moreover, food source F is selected from the list of non-dominated solutions with the least number of neighboring solutions, which leads the search toward the less crowded regions of the Pareto optimal front obtained and improves the coverage of solutions found [26]. As a result of these reasons, the MSSA is logically able to find accurate Pareto optimal solutions with high distribution across all objectives, and it is proved experimentally on benchmark problems.

3.2. Multi-Objective Hybrid Salp Swarm Algorithm for Solving the DFJSP

3.2.1. Algorithm Framework

According to the characteristics of the DFJSP, this paper designs a multi-objective hybrid salp swarm algorithm (MHSSA) to solve the double-flexible job shop green scheduling problem. The flowchart of the MHSSA is shown in Figure 1. The main operators for the studied problems are encoding and decoding, individual position update of the OS part, and crossover and mutation operation of MA and WA parts.

3.2.2. Encoding and Decoding Mechanisms

The double-flexible job shop green scheduling problem needs to solve three sub-problems: machine selection, worker selection, and process scheduling. According to the characteristics of the DFJSP, the multi-objective hybrid salp swarm algorithm designed in this paper needs three layers of salp individuals to represent: operation sequence (OS), machine assignment (MA), and worker assignment (WA). The length of each layer is equal to the total number of processes of all workpieces, as shown in Figure 2. The OS part adopts random key coding based on ascending arrangement, and the corresponding process sequencing is obtained by arranging the element values of each position in ascending order. The specific operation method is shown in the decoding process of the OS part in Figure 2. The MA part adopts machine codes arranged according to processes, which are arranged according to the order of workpieces and workpiece processes, and it is expressed by integers representing the machine number of the processing machine selected in the current process. The WA part adopts worker codes arranged according to processes, which are arranged according to the order of workpieces and workpiece processes, and it is expressed by integers representing the worker type number selected in the current process. For example, in Figure 2, the first numeral, 2, in the MA section indicates that the process O 11 is arranged to be processed on the second machine, and the fourth numeral, 1, indicates that the process O 21 is arranged to be processed on the first machine. The first number, 3, in the WA part indicates that the process O 11 has selected the third type of worker to operate. The other numbers represent the same meaning as above.
To sum up, the following information can be interpreted from Figure 2: in the scheduling scheme represented by the salp individual, process O 11 is processed first, and it is operated by a third type of worker on the second machine. The second process to be processed is O 21 , which will be processed by a third type of worker on the first machine. Then, other processes will be operated.
The processing sequence of each process and the corresponding processing machines and operators can be determined from the individual code of salp. Then, the method in [5] is adopted to decode and generate activity scheduling for the salp individual. The specific steps are as follows:
(1) Determine the machine set required to process each workpiece and the process set processed on each machine.
(2) Determine the allowable starting processing time A S of each process, that is, the processing end time of the previous process.
(3) Adopt a process insertion method to determine the earliest starting processing time S of each process. The method is as follows: Find all interval idle time periods [ T S , T E ] on the machine where the current process is located and judge whether each interval idle time period meets the insertion conditions from left to right. If so, insert it into the current idle time period, S = m a x ( T S , A S ) . Otherwise, S = m a x ( L M , A S ) , where LM represents the end time of the last operation of the current machine.
(4) Calculate the processing completion time C of each operation, where C = S + P and P represents the processing time of the operation.
By the above decoding process, the processing sequence of the workpieces on each machine and the corresponding start time and completion time can be obtained. After calculation, the corresponding objective function value can be obtained.

3.2.3. Individual Location Update Operations in OS Part

In the SSA algorithm, the foraging chain of salp is divided into two parts: leader and follower. Therefore, the individual position update of the population includes the individual position update of leaders and the individual position update of followers. At the same time, combined with the coding method of the OS part, the following methods are used to update the location of the OS part of the salp population.
The position update of the individual leader is related to the position of the food source F, and the position update method is shown in Equation (8).
x j i = { F j + c L r a n d < 0.5 F j c L r a n d 0.5
Suppose there are N individuals in the population and the OS part of each individual has D -dimensional variables (i.e., the total number of processes). x j i represents the j -dimensional variable value of the OS part of the i -th salp individual;
i = 1 , 2 , , N / 2 ; j = 1 , 2 , , D ; and F j represents the j -dimensional variable value of the OS part of the food source. The correlation coefficient c is the most important parameter in the SSA, which is used to balance global exploration and local development. See Equation (9), where t is the current number of iterations, t m a x is the maximum number of iterations of the population, and the value of c decreases from 2 to 0 with the iteration process. L is the Lévy flight step size, which can effectively prevent premature convergence of the algorithm. The calculation method [28] is shown in Equations (10)–(12), where β is a parameter between [1,2] and taken here to be equal to 1.5, u and v obey normal distribution, and r a n d is a random number between [0,1]. When r a n d < 0.5 , x j i changes in the positive direction of the food source F j ; otherwise, it changes in the negative direction.
c = 2 exp ( ( 4 t t max ) 2 )
L = u | ν | β
u ~ N ( 0 , σ u 2 ) , v ~ N ( 0 , σ v 2 )
σ u = Γ ( 1 + β ) sin ( π β / 2 ) Γ [ ( 1 + β ) / 2 ] β 2 ( β 1 ) / 2 β ,   σ v = 1
The individual position update of the follower is only related to the position of the salp individual in front of it, which is updated according to Equation (13), where x j i represents the variable value of the j -th dimension of the OS part of the i -th salp individual, i = ( N / 2 ) + 1 , ( N / 2 ) + 2 , , N , j = 1 , 2 , , D .
x j i = x j i + x j i 1 2

3.2.4. Cross Operations in MA and WA Parts

According to the characteristics of the coding mode of the salp individual, the crossover and mutation operators of the genetic algorithm are embedded in the MSSA to update the position of the salp individual in MA and WA parts. The crossover operation mainly helps obtain better offspring by enabling the exchange of excellent gene fragments between parent chromosomes. To avoid generating infeasible solutions, MA and WA parts are regarded as a whole for the crossover operation. Meanwhile, since the MA and WA parts must ensure that the sequence of each element remains unchanged, the random probability crossover operation method is adopted, as shown in Figure 3. Two individuals are randomly selected in the population, and a random number r a n d on [0,1] is generated from left to right for each position element. If r a n d is less than the crossover probability (CR), the position element remains unchanged; otherwise, the position element is crossed to produce two new individuals. The above steps are repeated until all individuals complete the crossover operation.

3.2.5. Mutation Operations in MA and WA Sections

Usually, mutation operators are used to jump out of the local optimal solution and increase the diversity of the population. In this paper, a random number r a n d is generated for each individual in the salp population. If r a n d is less than the mutation probability (MR), the mutation operation is performed. The specific method of mutation operation for MA and WA parts is as follows: randomly select a process, randomly select a device in the optional processing machine set of the process, and randomly select an operator in the corresponding optional worker set.

4. Simulation Experiment and Analysis

4.1. Experimental Setup and Test Examples

To verify the effectiveness of the MHSSA in solving the double-flexible job shop green scheduling problem, two multi-objective optimization algorithms, the MSSA [26] and the MOPSO [27], are selected for comparison, and each algorithm runs independently 20 times. The simulation platform includes the Win7 operating system, a CPU of 2.60 GHz, and 4G memory; MATLAB R2016a is used for programming.
In this paper, the DFJSP test example constructed in Reference [5] is selected for numerical experiments. The test problem is improved on the basis of the existing Brandimarte benchmark example [29]. Nine groups of data (DFJSP 01-09) are selected for testing the MHSSA in this paper; we can download these test data from https://pan.baidu.com/s/1mhHfv6K (accessed on 15 December 2021) [5] that in this example, with a decrease in the processing time, the labor cost will increase and the energy consumption will decrease. In addition, the example considers complete worker flexibility; that is, each machine can be operated by any worker.

4.2. Performance Evaluation Index

Since the measurement units of the three optimization objectives in this paper (maximum completion time, total labor cost, and total green index) are different, in order to evaluate the performance index of each algorithm, it is necessary to normalize the data of each target value. In this paper, three evaluation indexes SP, IGD, and Ω are selected. Since the optimal Pareto frontier of the tested problem is unknown, the non-dominated solution set obtained from the operation results of the three algorithms is approximated as the optimal Pareto frontier.
The SP index is used to measure the uniformity of the distribution of non-inferior solutions on the Pareto frontier. The calculation process is shown in Equation (14). Where d i = m i n j ( | f 1 i ( x ) f 1 j ( x ) | + | f 2 i ( x ) f 2 j ( x ) | + | f 3 i ( x ) f 3 j ( x ) | ) , i , j = 1 , , | P F | , d ¯ is the mean of all d i , and | P F | is the number of non-inferior solutions on the Pareto frontier. A smaller SP value means better uniformity of the distribution of non-inferior solutions. In particular, when the SP value is 0, the non-inferior solutions are equally uniformly distributed on the Pareto frontier.
S P [ 1 | P F | 1 i = 1 | P F | ( d ¯ d i ) 2 ] 1 / 2
IGD is a comprehensive index for simultaneously evaluating the diversity and frontier convergence of the obtained Pareto non-inferior solutions. The calculation is shown in Equation (15). | P F * | is the number of non-inferior solutions on the optimal Pareto frontier, and d i s t ( x , P F ) is the Euclidean distance between a solution on the optimal Pareto frontier and the nearest solution on the obtained Pareto frontier. The smaller the IGD value, the better the performance of the algorithm.
I G D = x P F d i s t ( x , P F ) | P F |
Ω is the dominant index, which is the percentage of the number of non-inferior solution sets found by the algorithm falling into the optimal Pareto frontier. The calculation process is shown in Equation (16), where H k represents the set of non-dominated solutions obtained by the k-th algorithm, and | P ( S ) | represents the number of non-dominated solutions in the set S . The larger the value of Ω , the better the dominance of the Pareto frontier found by the algorithm.
Ω H k = | P ( i H i ) \ P ( i k H i ) | | P ( i H i ) | × 100 %

4.3. Parameter Setting

The MHSSA proposed in this paper mainly involves three parameters: population size N , crossover probability C R , and mutation probability M R . To obtain the best combination of factor levels, the Taguchi method is used for experimental design to explore the influence of algorithm parameters on the solution effect [30], and DFJSP 01 is used as a test example.
The specific process is to set four factor levels for each parameter. According to the number of experimental groups, the number of experimental factors, and the number of factor levels, the orthogonal table of L 16 ( 4 3 ) is adopted, as shown in Table 1. For each factor-level combination in the orthogonal table, the MHSSA runs independently 20 times and records the obtained non-dominated solutions. Since the index Ω represents the percentage of the number of non-inferior solution sets falling into the optimal Pareto frontier, which can reflect the dominance of the algorithm, the index Ω is used as the response variable (RV). The larger the RV value, the better the performance of the parameter combination. The experimental results are shown in Table 2. The average value of response characteristics (ARV) and the importance ranking of each parameter are shown in Table 3.
As can be seen from Table 3, the population size N is the most significant among the three parameters. The N value is too large, the convergence speed of the algorithm is slow, and the running time is too long. A small N value may cause premature convergence of the algorithm. The variation probability M R is the second most significant, and the factor significance level of the crossover probability C R is the lowest. According to the above parameters, the main parameters of the MHSSA are set to N = 200, C R = 0.7, and M R = 0.3, which are the best combination of factor levels for simulation. Other parameters of the algorithm are set as follows: the maximum number of iterations t m a x = 100, and the external file A r c h i v e _ s i z e = 100.

4.4. Test Results and Analysis

To verify the effect of the proposed HMSSA to solve the DFJSP, nine test problems were selected to simulate the MHSSA. Each algorithm ran 20 times independently, and a group of [SP, IGD, Ω] values was recorded each time. Table 4 lists the mean (Avg) and standard deviation (Std) of each evaluation index of the three algorithms, and the optimal results of each index are marked in bold. It can be seen from the IGD and Ω indexes in Table 4 that the diversity, convergence, and dominance of the non-inferior solutions on the Pareto frontier obtained by the MHSSA are the best in the comparison algorithms and that the solution effects of the other two algorithms are not very different. For the SP index, the difference between the three algorithms is not obvious. For all test problems, the MHSSA designed in this paper has the best uniformity of the distribution of non-inferior solutions on the Pareto frontier.
Furthermore, the significant differences between the algorithms were tested using the Wilcoxon signed-rank test for relevant samples. For the three evaluation indexes of the test problem, the three algorithms are compared in pairs. The test results are shown in Table 5. The MHSSA has significant differences from the comparison algorithm in IGD and Ω indexes, and the effect is the best. For the SP index, the three algorithms have no significant difference. In conclusion, the MHSSA proposed in this paper can effectively solve the dual-flexible job shop green scheduling problem. On the whole, the MHSSA is obviously superior to the other two comparison algorithms in terms of diversity, convergence, and dominance of non-inferior solutions on the Pareto frontier.

5. Conclusions

In this paper, a multi-objective hybrid salp swarm algorithm is proposed to solve the dual-flexible job shop green scheduling problem with the objectives of optimizing the maximum completion time, total labor cost, and total green index. The specific research contents include the following:
(1)
According to the characteristics of the DFJSP, a three-layer individual coding and decoding method of salp is designed, and it is decoded into an active scheduling scheme. In the optimization process of the algorithm, the individual position of the OS part of the individual salp is updated based on the Lévy random walk, and the random probability crossover operation and corresponding mutation operation are performed on the MA and WA parts, which effectively improves the optimization accuracy and efficiency of the algorithm.
(2)
In the part of numerical simulation experiment, the Taguchi method was used to study the influence of the parameter setting on scheduling results and the best factor-level combination of algorithm parameters was identified. Using a simulation test for the benchmark example of the DFJSP, the MHSSA was compared with the MSSA and the MOPSO to verify that the MHSSA can effectively solve the dual-flexible job shop green scheduling problem.
In the next step, the DFJSP will be studied more extensively. In terms of the green scheduling problem, considering that the machine has a variety of speeds and other indicators, green manufacturing will be achieved by an actual production strategy.

Author Contributions

In this article, C.L. and Y.Y. designed the study. Their contributions are as follows: Conceptualization, C.L. and Y.Y.; methodology, C.L.; software, Y.Y.; validation, Y.Y.; writing—original draft preparation, Y.Y.; writing—review and editing, C.L.; visualization, H.Z.; supervision, H.Z.; funding acquisition, C.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The National Social Science Foundation of China, grant number 17BSH040, 20BGL006 and the Natural Science Foundation of the Jiangsu Higher Education Institutions of China, grant number 20KJB520036.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data used to support the findings of this study are included within the article.

Acknowledgments

The authors are thankful to the anonymous reviewers and the editor for their valuable comments and constructive suggestions that have led to an improved version of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Boudreau, J.; Hopp, W.; McClain, J.-O.; Thomas, L.J. On the interface between operations and human resources management. Manu. Serv. Opera. Manag. 2003, 5, 179–202. [Google Scholar] [CrossRef] [Green Version]
  2. Gino, F.; Pisano, G. Toward a theory of behavioral operations. Manu. Serv. Opera. Manag. 2008, 10, 676–691. [Google Scholar] [CrossRef] [Green Version]
  3. Neumann, W.-P.; Dul, J. Human factors: Spanning the gap between OM and HRM. Int. J. Opera. Prod. Manag. 2010, 30, 923–950. [Google Scholar] [CrossRef]
  4. Ye, C. Research on a new model of behavioral production scheduling based on learning effect. Enterp. Econ. 2015, 35, 5–10. [Google Scholar]
  5. Gong, G.; Deng, Q.; Gong, X.; Liu, W.; Ren, Q. A new double flexible job-shop scheduling problem integrating processing time, green production, and human factor indicators. J. Clean. Prod. 2018, 174, 560–576. [Google Scholar] [CrossRef]
  6. He, Y.; Liu, F.; Cao, H.-J.; Li, C.-B. A Bi-objective model for job-shop scheduling problem to minimize both energy consumption and makespan. J. Cent. South Univ. Tech. 2005, 12, 167–171. [Google Scholar] [CrossRef]
  7. Mouzon, G.; Mehmet, B.-Y.; Twomey, J. Operational methods for minimization of energy consumption of manufacturing equipment. Int. J. Prod. Res. 2007, 45, 4247–4271. [Google Scholar] [CrossRef] [Green Version]
  8. Dai, M.; Tang, D.; Giret, A.; Salido, M.-A.; Li, W.-D. Energy-efficient scheduling for a flexible fow shop using an improved genetic-simulated annealing algorithm. Rob. Comput. Integr. Manuf. 2013, 29, 418–429. [Google Scholar] [CrossRef]
  9. Liu, C.-H.; Huang, D.-H. Reduction of power consumption and carbon footprints by applying multi-objective optimisation via genetic algorithms. Int. J. Prod. Res. 2014, 52, 337–352. [Google Scholar] [CrossRef]
  10. Mansouri, S.-A.; Aktas, E.; Besikci, U. Green scheduling of a two-machine fowshop: Trade-off between makespan and energy consumption. Eur. J. Oper. Res. 2016, 248, 772–788. [Google Scholar] [CrossRef] [Green Version]
  11. Subaï, C.; Baptiste, P.; Niel, E. Scheduling issues for environmentally responsible manufacturing: The case of hoist scheduling in an electroplating line. Int. J. Prod. Econ. 2006, 99, 74–87. [Google Scholar] [CrossRef]
  12. Fang, K.; Uhan, N.; Zhao, F.; Sutherland, J.W. A new approach to scheduling in manufacturing for power consumption and carbon footprint reduction. J. Manuf. Syst. 2011, 30, 234–240. [Google Scholar] [CrossRef]
  13. Diaz, N.; Redelsheimer, E.; Dornfeld, D. Energy consumption characterization and reduction strategies for milling machine tool use. In Proceedings of the 18th CIRP International Conference on Life Cycle Engineering, Braunschweig, Germany, 2–4 May 2011; pp. 263–267. [Google Scholar]
  14. Wang, S.; Lu, X.; Li, X.-X.; Li, W.-D. A systematic approach of process planning and scheduling optimization for sustainable machining. J. Clean. Prod. 2015, 87, 914–929. [Google Scholar] [CrossRef]
  15. Zhang, R.; Chiong, R. Solving the energy-efficient job shop scheduling problem: A multi-objective genetic algorithm with enhanced local search for minimizing the total weighted tardiness and total energy consumption. J. Clean. Prod. 2016, 112, 3361–3375. [Google Scholar] [CrossRef]
  16. Gahm, C.; Denz, F.; Dirr, M.; Tuma, A. Energy-efficient scheduling in manufacturing companies: A review and research framework. Eur. J. Opera. Res. 2016, 248, 744–757. [Google Scholar] [CrossRef]
  17. Wang, L.; Wang, J.; Wu, C. Advances in green shop scheduling and optimization. Control Deci. 2018, 33, 385–391. [Google Scholar]
  18. Lei, D.-M.; Zheng, Y.-L.; Guo, X.-P. A shuffled frog-leaping algorithm for flexible job shop scheduling with the consideration of energy consumption. Int. J. Prod. Res. 2017, 55, 3126–3140. [Google Scholar] [CrossRef]
  19. Wu, X.; Sun, Y. Flexible job shop green scheduling problem with multi-speed machine. Comput. Integr. Manu. Sys. 2018, 24, 862–875. [Google Scholar]
  20. Dong, J.; Ye, C. Research on collaborative optimization of green manufacturing in semiconductor wafer distributed heterogeneous factory. Appl. Sci. 2019, 9, 2879–2903. [Google Scholar] [CrossRef] [Green Version]
  21. Yao, Y.-Y.; Ye, C.-M. IMSSA for TFT-LCD panel array process scheduling problem considering energy saving. China Mechan. Eng. 2019, 30, 2994–3003. [Google Scholar]
  22. Nanthapodej, R.; Liu, C.; Nitisiri, K.; Pattanapairoj, S. Variable neighborhood strategy adaptive search to solve parallel-machine scheduling to minimize energy consumption while considering job priority and control makespan. Appl. Sci. 2021, 11, 5311–5332. [Google Scholar] [CrossRef]
  23. He, Y.; Liu, F.; Cao, H.-J.; Li, C.-B. Scheduling optimization of mechanical processing system for green manufacturing model. Chin. J. Mechan. Eng. 2007, 43, 27–33. [Google Scholar] [CrossRef]
  24. Li, X.-X.; Huang, X.-M.; Liu, J.-X.; Liu, F. Optimization simulation for job-shop scheduling for reducing manufacturing energy consumption. J. Sys. Simu. 2016, 28, 114–120. [Google Scholar]
  25. Lu, C.; Gao, L.; Li, X.-Y.; Pan, Q.; Wang, Q. Energy-efficient permutation flow shop scheduling problem using a hybrid multi-objective backtracking search algorithm. J. Clean. Prod. 2017, 144, 228–238. [Google Scholar] [CrossRef]
  26. Mirjalili, S.; Gandomi, A.-H.; Mirjalili, S.-Z.; Saremi, S.; Faris, H.; Mirjalili, S.M. Salp swarm algorithm: A bio-inspired optimizer for engineering design problems. Adv. Eng. Soft. 2017, 114, 163–191. [Google Scholar] [CrossRef]
  27. Coello, C.-A.; Pulido, G.-T.; Lechuga, M.-S. Handling multiple objectives with particle swarm optimization. IEEE Trans. Evolut. Comput. 2004, 8, 256–279. [Google Scholar] [CrossRef]
  28. Liu, C.-P.; Ye, C.-M. Bat algorithm with the characteristics of lévy flight. CAAI Trans. Intell. Sys. 2013, 8, 240–246. [Google Scholar]
  29. Paolo, B. Routing and scheduling in a flexible job shop by tabu search. Ann. Opera. Res. 1993, 41, 157–183. [Google Scholar]
  30. Montgomery, D.-C. Design and Analysis of Experiments, 7th ed; John Wiley & Sons: Hoboken, NJ, USA, 2008; pp. 50–60. [Google Scholar]
Figure 1. Flowchart of a multi-objective hybrid salp swarm algorithm.
Figure 1. Flowchart of a multi-objective hybrid salp swarm algorithm.
Applsci 12 00205 g001
Figure 2. Individual coding and OS part decoding process of the salp agent.
Figure 2. Individual coding and OS part decoding process of the salp agent.
Applsci 12 00205 g002
Figure 3. Random probability crossover operation.
Figure 3. Random probability crossover operation.
Applsci 12 00205 g003
Table 1. Factor level setting of MHSSA parameters.
Table 1. Factor level setting of MHSSA parameters.
FactorFactor Level
1234
N50100150200
CR0.20.50.70.9
MR0.10.20.30.4
Table 2. Orthogonal table and response eigenvalues.
Table 2. Orthogonal table and response eigenvalues.
Experimental Serial NumberFactorRV
NCRMR
11110.0098
21220.0098
31330.0392
41440.0196
52120.0196
62210.0392
72340.0392
82430.0784
93130.0686
103240.0686
113310.0882
123420.0490
134140.0882
144230.1471
154320.1275
164410.1078
Table 3. Response table.
Table 3. Response table.
Factor LevelNCRMR
10.01960.046550.06125
20.04410.0661750.051475
30.06860.0735250.083325
40.117650.06370.0539
Delta0.098050.0269750.03185
Rank132
Table 4. Simulation test results of three algorithms.
Table 4. Simulation test results of three algorithms.
TableEvaluation IndexMSSAMOPSOHMSSA
AvgStdAvgStdAvgStd
DFJSP   01   ( 10   ×   6 ×   3 ) SP0.05620.01540.04880.01500.06910.0223
IGD0.19120.01990.18820.01580.09400.0100
Ω0.00250.00660.00250.00520.04490.0319
DFJSP   02   ( 10   ×   6 ×   3 ) SP0.06500.01790.06470.02390.06090.0164
IGD0.18640.01680.18540.02370.10420.0155
Ω0.00290.00630.00750.01010.03970.0355
DFJSP   03   ( 15 ×   8 ×   3 ) SP0.06770.01690.07610.02970.06920.0150
IGD0.19710.01330.19470.02080.11780.0139
Ω0.00380.00630.00910.01060.03710.0389
DFJSP   04   ( 15 ×   8 ×   3 ) SP0.06190.01890.07080.02590.06490.0148
IGD0.18670.02040.19500.01460.10290.0121
Ω0.00160.00390.00320.00610.04520.0303
DFJSP   05   ( 15 ×   4 ×   3 ) SP0.06380.02020.06760.03160.05990.0121
IGD0.13980.01330.13820.00870.07550.0075
Ω0.00580.00700.00650.00820.03770.0304
DFJSP   06   ( 10 ×   15 ×   3 ) SP0.06030.01210.05700.01380.06780.0151
IGD0.18090.02080.18480.01840.10320.0072
Ω0.00840.01260.00740.00840.03420.0285
DFJSP   07   ( 20 ×   5 ×   3 ) SP0.06260.01740.06720.02750.06900.0144
IGD0.18000.02290.18090.01570.09510.0101
Ω0.00520.00820.00570.00760.03910.0317
DFJSP   08   ( 20 ×   10 ×   3 ) SP0.05900.01950.06180.01730.05660.0078
IGD0.17220.01510.18570.01610.09500.0104
Ω0.00720.01040.00110.00260.04170.0375
DFJSP   09   ( 20 ×   10 ×   3 ) SP0.06040.01100.06900.02260.05740.0148
IGD0.17770.01540.18010.01580.09750.0111
Ω0.00380.00480.00530.00960.04090.0289
TotalSP0.06190.00340.06480.00810.06390.0052
IGD0.17910.01650.18140.01700.09840.0113
Ω0.00460.00230.00540.00260.04010.0036
Table 5. Wilcoxon signed-rank test (the MHSSA vs. the MSSA and the MOPSO).
Table 5. Wilcoxon signed-rank test (the MHSSA vs. the MSSA and the MOPSO).
Evaluation IndexMSSAMOPSO
p ValueSig (p < 0.05)p ValueSig (p < 0.05)
SP0.477N0.515N
IGD0.008Y0.008Y
Ω0.008Y0.008Y
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, C.; Yao, Y.; Zhu, H. Hybrid Salp Swarm Algorithm for Solving the Green Scheduling Problem in a Double-Flexible Job Shop. Appl. Sci. 2022, 12, 205. https://doi.org/10.3390/app12010205

AMA Style

Liu C, Yao Y, Zhu H. Hybrid Salp Swarm Algorithm for Solving the Green Scheduling Problem in a Double-Flexible Job Shop. Applied Sciences. 2022; 12(1):205. https://doi.org/10.3390/app12010205

Chicago/Turabian Style

Liu, Changping, Yuanyuan Yao, and Hongbo Zhu. 2022. "Hybrid Salp Swarm Algorithm for Solving the Green Scheduling Problem in a Double-Flexible Job Shop" Applied Sciences 12, no. 1: 205. https://doi.org/10.3390/app12010205

APA Style

Liu, C., Yao, Y., & Zhu, H. (2022). Hybrid Salp Swarm Algorithm for Solving the Green Scheduling Problem in a Double-Flexible Job Shop. Applied Sciences, 12(1), 205. https://doi.org/10.3390/app12010205

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop