Next Article in Journal
A Review of Symmetry-Based Open-Circuit Fault Diagnostic Methods for Power Converters
Previous Article in Journal
A New Test Statistic to Assess the Goodness of Fit of Location-Scale Distribution Based on Progressive Censored Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Adaptive Two-Class Teaching-Learning-Based Optimization for Energy-Efficient Hybrid Flow Shop Scheduling Problems with Additional Resources

School of Automation, Wuhan University of Technology, Wuhan 430070, China
*
Author to whom correspondence should be addressed.
Symmetry 2024, 16(2), 203; https://doi.org/10.3390/sym16020203
Submission received: 11 January 2024 / Revised: 6 February 2024 / Accepted: 7 February 2024 / Published: 8 February 2024

Abstract

:
Energy-efficient scheduling problems with additional resources are seldom studied in hybrid flow shops. In this study, an energy-efficient hybrid flow shop scheduling problem (EHFSP) with additional resources is studied in which there is asymmetry in the machine. An adaptive two-class teaching-learning-based optimization (ATLBO) which has multiple teachers is proposed to simultaneously minimize the makespan and the total energy consumption. After two classes are formed, a teacher phase is first executed, which consists of teacher self-learning and teacher training. Then, an adaptive learner phase is presented, in which the quality of two classes is used to adaptively decide the learner phase or the reinforcement search of the temporary solution set. An adaptive formation of classes is also given. Extensive experiments were conducted and the computational results show that the new strategies are effective and that ATLBO was able to provide better results than comparative algorithms reported in the literature in at least 54 of 68 instances.

1. Introduction

Hybrid flow shop scheduling problems (HFSPs) occur frequently in many manufacturing industries, such as chemical engineering, metallurgy, textiles, petrochemicals, airplane engines, and semiconductors [1]. Compared with a flow shop, a hybrid flow shop has at least two stages, including at least one stage with at least two parallel machines. The redundancy of the machines used may enable dealing with jobs flexibly and avoid bottlenecks [2]. In recent decades, HFSPs have been extensively investigated and a number of results have been obtained [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17].
With the development and application of green manufacturing, companies and society as a whole are paying increased attention to reducing energy consumption in the manufacturing process and improving the rates of energy utilization. In recent years, as a typical meansd of scheduling in green manufacturing, EHFSPs have attracted much attention with a number of investigations being reported. Yan et al. [3] developed a multi-level optimization method to optimize the total energy consumption and the makespan. Zeng et al. [4] investigated the conflicting relationship between total energy consumption and the makespan and employed the non-dominated sorting genetic algorithm-II (NSGA-II). Li et al. [5] presented a two-level imperialist competitive algorithm (TICA). Zuo et al. [6] designed a novel multi-population artificial bee colony (ABC) algorithm. Song [7] constructed a mixed integer linear programming (MILP) model and developed a hybrid multi-objective teaching-learning-based optimization (TLBO) based on decomposition to minimize the makespan and energy consumption. Yue et al. [8] designed a hybrid Pareto spider monkey optimisation algorithm for solving problems, minimizing the total energy consumption and total time taken. Chen et al. [9] reported an improved NSGA-II including optimization strategies. Qin et al. [10] presented a mathematical model of a blocking EHFSP and proposed a modified iterative greedy algorithm based on a swap strategy. Zuo et al. [11] designed a population diversity-based ABC algorithm to minimize the makespan and the total energy consumption simultaneously. Wang et al. [12] considered a fuzzy EHFSP with variable machine speed and proposed use of the NSGA-II algorithm to minimize the fuzzy makespan and total fuzzy energy consumption. Li et al. [13] proposed use of the NSGA-II algorithm combined with Q-learning and general variable neighborhood search. Wang et al. [14] presented an improved multi-objective firefly algorithm.
In scheduling problems, the additional resources involved mainly consist of machine operators, tools, pallets, dies, industrial robots, and automated guided vehicles. The additional resources may be renewable or non-renewable, discrete or continuous, and be used in the processing or setup [18]. When addressing scheduling problems, there are additional resources constraints, that is, the sum of each additional resource used by all machines cannot exceed a given limit at any moment. The existing published studies mainly concern unrelated parallel machine scheduling problems (UPMSPs) with additional resources [19,20,21,22,23,24,25,26,27].
With respect to UPMSPs with additional resources, Zheng and Wang [19,20] presented a two-stage adaptive fruit fly optimization algorithm using a heuristic and knowledge-guided search, and also described use of a collaborative multi-objective fruit fly optimization algorithm. Yepes-Borrero et al. [21] presented three heuristics and a greedy randomized adaptive search procedure for UPMSPs with setup times and additional limited resources in the setup. Fanjul-Peyro et al. [22] constructed two integer linear programming models and proposed three matheuristics. Fanjul-Peyro [23] focused on a UPMSP with a processing resource, setup resource, and a shared resource, and presented a MILP model and a three-phase algorithm. Pinar and Seyda [24] solved a multi-resource-constrained UPMSP with sequence-dependent setup times, a precedence relation, machine eligibility, and release dates using a constraint programming model and two branching strategies. Bitar et al. [25] investigated a UPMSP with auxiliary resources and sequence-dependent and machine-dependent setup times in the photolithography workshop of a semiconductor plant. A mimetic algorithm combined with an integer linear programming model was used to minimize the makespan. Lei and He [26] developed an adaptive ABC approach for addressing a problem with additional resources involving preventive maintenance. An ABC method with adaptive competition was also applied to solve the above problem [27].
Unlike UPMSPs with additional resources, HFSPs with additional resources are seldom considered. For a two-stage HFSP with limited additional resources, Figielska [28,29] proposed a heuristic with several rules for job selection and two heuristic algorithms by combining a column generation technique with a genetic algorithm or column generation with simulated annealing.
As stated above, many works have been published on EHFSPs and UPMSPs with additional resources. Only one paper has considered solving HFSPs with additional resources, while EHFSPs with additional resources in which the additional resource, machine resource, and energy are dealt with simultaneously, have been barely investigated. An additional resource constraint is needed to be handled at each stage and multiple conflicting objectives are optimized simultaneously in EHFSPs with additional resources. As a result, optimization difficulties increase greatly. Moreover, the reasonable usage of additional resources and the optimization of energy consumption will result in EHFSPs with additional resources being close to the realistic situation of a hybrid flow shop. The corresponding optimization results will have high application value, so it is important to consider EHFSPs with additional resources.
TLBO simulates the teaching process of teachers and the learning process of learners in real classes, which consists of a teacher phase and a learner phase [30]. Compared with other algorithms, TLBO has a simple structure, fewer parameters, and is both easy to understand and implement. It has become an effective method for addressing production scheduling problems [7,31,32,33,34,35,36,37,38]. Lei et al. [34] proposed a novel TLBO for EHFSPs. Li et al. [35] applied a discrete TLBO to solve an HFSP with five types of disruption events and makespan minimization. Xu et al. [36] reported an effective TLBO for HFSPs with fuzzy processing time. Buddala and Mahapatra [37] presented a two-stage TLBO to solve an HFSP involving unexpected machine breakdowns. Lei and Xi [38] proposed a Q-learning-based TLBO for two-stage fuzzy HFSPs.
The above studies show that TLBO has strong search capabilities and advantages for solving HFSPs and EHFSPs. EHFSPs with additional resources and EHFSPs have strong similarities. For example, the same coding can be used for two problems. The strong search abilities of TLBO for HFSPs and EHFSPs show that TLBO possesses potential advantages in solving EHFSPs with additional resources, so TLBO is chosen here, and how to design an efficient TLBO for EHFSPs with additional resources is investigated.
In this study, energy-efficient scheduling in a hybrid flow shop is considered. The main contributions are summarized as follows:
  • An EHFSP with additional resources is considered.
  • An adaptive two-class teaching-learning-based optimization (ATLBO) is presented to minimize the makespan and total energy consumption simultaneously. To produce high quality solutions, a teacher phase with teacher self-learning and teacher training are presented and an adaptive learner phase is reported which uses the quality of two classes to decide adaptively the learner phase or the reinforcement search of the temporary solution set. An adaptive formation of classes is also performed.
  • The performance of the ATLBO is tested through a number of experiments. The effectiveness of the new strategies of the ATLBO is validated and the search advantages of the ATLBO for solving the EHFSP with additional resources is also demonstrated.
The rest of the paper is arranged as follows: Section 2 describes the considered EHFSP with additional resources. Section 3 shows the detailed steps of the ATLBO for the EHFSP with additional resources. The computational experiments undertaken are descibed in Section 4. The conclusions and future topics are reported in Section 5.

2. Problem Description

The EHFSP with additional resources is described as follows: There are n jobs J 1 , J 2 , , J n and a hybrid flow shop with m stages. Each stage l has S l unrelated parallel machines M l 1 , M l 2 , , M l S l . Each job J i is processed according to the same production flow: stage 1, stage 2, ⋯, stage m 2 . p i l k indicates the processing time of job J i on M l k . One additional resource is needed for each stage. For job J i processed on M l k , it needs r i l k units of the additional resource. At most, R l units of additional resource can be used at any time on each machine at stage l. Machine M l k has two modes: a processing mode and an idle mode. e l k and i e l k indicate the energy consumption per unit time when M l k is in the processing mode and the idle mode, respectively.
An additional resource constraint is that the sum of each additional resource used by all the machines cannot exceed a given limit at any moment. If there are no additional resources, symmetry exists between any two adjacent jobs on the same machine; after additional resources are considered, asymmetry often exists for adjacent jobs on a machine.
The following constraints on jobs and machines are considered:
Each job can be processed on only one machine at a time.
Each machine handles at most one job at a time.
The operations cannot be interrupted.
All the machines are available at all times.
The goal is to to minimize the makespan and the total energy consumption.
f 1 = min C m a x = max C i i = 1 , 2 , , n
f 2 = T E C = l = 1 m k = 1 S l i = 1 n 0 C max w i l k t × e l k d t + l = 1 m j = 1 S l 0 C max z l k t × i e l k d t ,
where n indicates n jobs, C i indicates the completion time of job J i , C m a x is the maximum completion time of all jobs, m indicates the m stages, S l indicates that there are S l machines at stage l, and e l k and i e l k indicate the energy consumption per unit time when M l k is in the processing mode and the idle mode, respectively. w i l k t is 1 if job J i is processed on M l k at time t and 0 otherwise. z i l k t is 1 if machine M l k is free at time t and 0 otherwise. f 1 is the minimized makespan. f 2 is the total energy consumption. T E C is an acronym for the total energy consumption.
The Pareto domination is a fundamental concept in multi-objective optimization problems [39]. For the above problem, if x and y meet i { 1 , 2 } ,   f i ( x ) f i ( y ) and i { 1 , 2 } , f i ( x ) < f i ( y ) , then x dominates y; x y denotes x dominates y.
An example of an EHFSP with additional resources, two stages and eight jobs is given in Table 1. Every stage has two machines. In Table 1, the number outside the bracket is the processing time and the data in the bracket are r i l k . R l = 10 , ( l = 1 , 2 ) , e 11 = 10 , e 12 = 10 , e 21 = 3 , e 22 = 2 , i e l k = 1 . Figure 1 gives a schedule of the example, where f 1 = 652 , f 2 = 2589 .

3. ATLBO for EHFSP with Additional Resources

TLBO usually has only one class. A multi-class TLBO is considered to improve the performance of TLBO for scheduling problems [32,38]. These algorithms have more than two classes and the number of classes is the algorithm parameter. In this study, a two-class TLBO is presented in which the number of classes is fixed to two and an adaptive learner phase is used.

3.1. Initialization and Formation of Two Classes

A two-string representation is applied. For an EHFSP with n jobs and m stages, its solution is composed of a scheduling string [ q 1 , q 2 , , q n ] and a machine assignment string [ M 1 h 1 , M 1 h 2 , , M 1 h n ] , where M 1 h 1 is the assigned machine for job J i at stage 1 and q 1 { 0 , 1 } is a real number for J i .
The decoding procedure is shown below.
Step 1: For stage l = 1 , assign all jobs on each machine M 1 k , k = 1 , 2 , , S 1 according to the machine assignment string; then, on each machine M 1 k , sort all assigned jobs on M 1 k in the ascending order of their q i ; obtain a permutation, start with the first job of the permutation. For each job J i in the permutation, process J i on machine M 1 k when the additional resource constraint is met.
Step 2: For stage l = 2 , 3 , , m , sort all jobs in the ascending order of their completion time at stage l 1 ; then, obtain a permutation of all jobs [ π 1 , π 2 , , π n ] . Start with the first job of the permutation; for each job π i , choose the machine M l k with the smallest beginning time and process π i on M l k when the additional resource constraint is met.
For job J i , when it is processed on machine M l k at stage l, the handling method [26] for the additional resource is directly applied. J i can be inserted into the idle period means that J i can be started and completed during the idle time and at most R l units of the additional resource can be used at any time on all machines at stage l.
A solution of the example in Table 1 is a machine assignment string [ 0.07 , 0.48 , 0.05 , 0.20 , 0.37 , 0.23 , 0.69 , 0.32 ] and a scheduling string [ 1 , 1 , 2 , 2 , 1 , 2 , 1 , 2 ] . At stage 1, jobs J 1 , J 2 , J 5 , J 7 are assigned on machine M 11 , a permutation of jobs [ 1 , 5 , 2 , 7 ] is obtained on M 1 k in the ascending order of q 1 = 0.07 , q 2 = 0.48 , q 5 = 0.37 , q 7 = 0.69 , J 7 can be processed in the idle period [ 0 , 111 ] , and the additional resource can also be met, J 7 is inserted into idle period [ 0 , 111 ] . At stage 2, obtain a permutation of all jobs [ 7 , 3 , 4 , 1 , 6 , 8 , 5 , 2 ] in the ascending order of their completion time at stage 1, start with the first job J 7 of the permutation, then all the jobs are processed in order. J 5 can be inserted into the idle period [ 77 , 191 ] .
An initial population P with N solutions is randomly generated. Let λ = 0 , λ is an integer parameter. The best b solutions from population P are selected as teachers, which constitute the set of teachers Θ = { x t e a c h e r 1 , x t e a c h e r 2 , , x t e a c h e r b } . Then, the N b solutions as learners are arranged randomly for two classes C l s 1 , C l s 2 .
A formation of two classes is described as follows:
Step 1: Perform non-dominated sorting and compute the crowding distance [40] on all the solutions in P.
Step 2: Compute the fitness c ( x i ) of each solution x i , sort all solutions in the descending order of c ( x i ) , suppose that c ( x 1 ) c ( x 2 ) c ( x N ) , then select the best b solutions as teachers.
Step 3: For solutions x b + 1 , x b + 2 , , x N , let l = b + 1 , τ = 1 ; repeat the following steps until l > N : x l is assigned to C l s τ , l = l + 1 , τ = τ + 1 . If τ = 3 , then let τ = 1 .
c x i = max l = 1 , 2 , , N r a n k l r a n k i + d i s t i / j θ r a n k i d i s t j ,
where r a n k i and d i s t i indicate the r a n k value and the crowding distance of the solution x i .

3.2. Teacher Phase

The teacher phase is depicted below: for each teacher x t e a c h e r r , execute a multiple neighborhood search on it, then randomly assign ( N b ) / b learners, which are included into a set Λ r . For each learner x i Λ r , perform a global search between x i and x t e a c h e r r , obtain a new solution z. If z x i or z and x i are non-dominated over each other, then replace x i with z.
A multiple neighborhood search is executed. For solution x, let g = 1 . Repeat the following steps until g = 4 : produce a new solution z N g ( x ) , if z x or z and x is non-dominated over each other, then replace x with z and let g = 4 ; otherwise, g = g + 1 . Where N 1 and N 2 are s w a p and i n s e r t , respectively, applied on the scheduling string, N 3 is c h a n g e executed on the machine assignment string. The descriptions on s w a p , i n s e r t , c h a n g e can be found in [34]. N i ( x ) is the set of neighborhood solutions generated by N i .
The global search is described below. For two solutions x , y , if a random number r a n d < μ , then order crossover [41] is executed on the scheduling string of x , y ; otherwise, a two-point crossover [42] is applied on the machine assignment string, where μ follows a uniform distribution on [0, 1], and r a n d is a real number.
Generally, the computational complexity of a scheduling sub-problem is notably greater than that of the machine assignment. To obtain the best solution to address the problem efficiently, more computing resource is assigned for scheduling the sub-problem, so we set μ > 0.5 to adapt to the above features of the two sub-problems.
In the teacher phase, each teacher is first improved by self-learning implemented with a multiple neighborhood search. Then, some learners are randomly assigned for each teacher and teaching activities are conducted between each teacher and its assigned learners. The self-learning and teaching take full advantage of high quality solutions. As a result, higher quality solutions can be obtained with high probability; moreover, the global search ability can be improved because of the random assignment of learners.

3.3. Adaptive Learner Phase

The detailed steps of the adaptive learner phase are shown below.
Step 1: Execute a multiple neighborhood search on each learner x i ( C l s 1 C l s 2 ) .
Step 2: Compute r a n k and the crowding distance of learner x i ( C l s 1 C l s 2 ) [40]; then calculate E v o i for each class C l s i .
Step 3: If E v o g < β , E v o 3 g β , g { 1 , 2 } , then let λ = 0 ; for each learner x i C l s g , produce a random number r a n d .
If r a n d E v o 3 g / E v o 1 + E v o 2 , then randomly select y , z C l s 3 g ; execute global search between y and z, produce a new solution z . If z x i or z and x i are non-dominated over each other, then replace x i with z ; otherwise, execute a multiple neighborhood search on x i .
If r a n d < E v o 3 g / E v o 1 + E v o 2 , then randomly select z C l s 3 g . Execute a global search between x i and z, produce a new solution z . If z x i or z and x i are non-dominated over each other, then replace x i with z ; otherwise, carry out a multiple neighborhood search on x i .
Step 4: If E v o 1 < β , E v o 2 < β , then let λ = λ + 1 , a set of solutions Φ , which is composed of the updated solutions in Step 1 and all teachers in the population P. For each x Φ , randomly select y Φ , execute global search between x and y, produce a new solution z. If z x or z and x are non-dominated over each other, then replace x with z; otherwise, perform a multiple neighborhood search on x, produce a new solution z . If z x or z and x are non-dominated over each other, then replace x with z ; otherwise, replace the worst solution of C l s 1 C l s 2 with z . If λ 3 , then return each x Φ to the original class; otherwise, perform an adaptive formation of classes.
Where β is a real number, and E v o i is the evaluation index of class C l s i ( i = 1 , 2 ) and computed by
E v o i = x C l s i c ( x ) / N b 2
The detailed steps of the adaptive formation of classes are shown below.
Step 1: Perform non-dominated sorting and crowding distance [40] on all solutions in P.
Step 2: Compute the fitness c ( x i ) of each solution x i . Sort all the solutions in the descending order of c ( x i ) . Suppose that c ( x 1 ) c ( x 2 ) c ( x N ) ; the best b solutions are selected as teachers.
Step 3: If Φ Θ ( N b ) / 2 , then select the ( N b ) / 2 best solutions in Φ Θ and put them into C l s 1 . All the remaining solutions are put into C l s 2 ; otherwise, assign all solutions x Φ Θ to C l s 1 first. All the remaining solutions are randomly assigned to two classes, ensuring equal numbers in both classes.
Where Φ Θ is a set of all the updated solutions in Step 1 of the adaptive learner phase.
Usually, the learner phase is about mutual learning between learners. If the learners in a class have poor solution quality, then their mutual learning will have low efficiency; as a result, waste of computing resources occurs. To avoid the above case, an adaptive learner phase is used. If the conditions of Step 3 are met, mutual learning is just applied for class with higher E v o i ; if the conditions of Step 4 are met, then no mutual learning occurs and a reinforcement search of the temporary set is executed to update the worst solution of the class. If E v o 1 > β , E v o 2 > β , then the classes evolve well in the teacher phase, so the learner phase is skipped directly.

3.4. Algorithm Description

The detailed steps of ATLBO are shown below.
Randomly produce an initial population P with N solutions, let λ = 0 , compute c ( x i ) for each solution.
While the stopping condition is not met, do
Form two classes.
Execute the teacher phase.
Perform the adaptive learner phase.
End While
Step 6: Output the non-dominated in P.
A flowchart of ATLBO is shown in Figure 2.
Unlike the existing TLBOs [32,38], ATLBO just uses two classes. After two classes are formed, a teacher phase with self-learning and teaching is executed on multiple teachers. Then, an adaptive learner phase is presented in which the quality of the two classes is used to decide adaptively the learner phase or the reinforcement search of the temporary solution set, and an adaptive formation of classes is used. Resources are given to better solutions to produce high quality solutions and to save computing resources so that the search efficiency can be improved.

4. Computational Experiments

Extensive experiments were conducted to test the performance of ATLBO for an EHFSP with additional resources. All experiments were implemented using Microsoft Visual C++ 2019 and run on 16.0 G RAM 2.00 GHz CPU PC.

4.1. Test Instances, Metrics, and Comparative Algorithms

A total of 68 instances are used [26], which are the combinations of n [8, 16, 30, 50, 150, 250, 350] and m [ 2 ,   4 ,   6 ,   8 ] . The detailed descriptions of their data, such as p i l k and r i l k , can be obtained directly from the following link: http://soa.iti.es/ (accessed on 10 January 2024). R l = 5 S l . Each instance is described as n × m × No . There is one way to produce the processing time and three ways to generate the additional resources. There are three combinations of processing time and additional resources. No is defined as a combination of the one way of producing the processing time and the No -th way of generating the additional resources.
The following three metrics are chosen in this study.
The metric C [43] is applied to measure the dominant relationship between the non-dominated solution sets of the two algorithms.
C ( L , B ) = | { b B : h L ,   h b } | | B |
C ( L , B ) measures the fraction of members of B that are dominated by members of L.
The metric ρ [44] is the ratio of { x Ω A x Ω * } to Ω * , where Ω A is the non-dominated set of algorithm A. The reference set Ω * consists of the non-dominated solutions in the union of the non-dominated sets of all the algorithms.
ρ Ω A = x Ω A x Ω * Ω *
The metric D I R [45] measures the non-dominated solution set Ω A of an algorithm with the reference set Ω * .
D I R Ω A = 1 Ω * y Ω * min σ x y x Ω A ,
where σ x y is the distance between a solution x and a reference solution y in the normalized objective space.
Meta-heuristics [46,47,48,49] are often applied to solve complicated optimization problems. In this study, three meta-heuristics are selected, which include an improved multi-objective firefly algorithm (IMOFA) [14], a two-level imperialist competitive algorithm (TICA) [5], and a hybrid multi-objective teaching-learning-based optimization (HMOTLBO) [7]. These algorithms can be used to solve an EHFSP with additional resources after minor revisions are made, so they are chosen as the comparative algorithms.
Related steps on additional resources are added in three comparative algorithms to deal with the considered problem. For IMOFA, the speed determination matrix and its related crossover operators and neighborhood searches are deleted; only the makespan and total energy consumption were chosen as objectives. For TICA, a three-string representation was changed to only a scheduling string and a machine assignment string. Only two crossovers on the scheduling strings and machine assignment string are used. HMOTLBO can be directly used.
For TLBO, the teacher phase and the learner phase are executed by global search. Comparisons between ATLBO and TLBO are performed to show the effect of the new strategies of ATLBO.

4.2. Parameter Settings

ATLBO has the following parameters: N, b, μ , β , λ and a stopping condition. With respect to the stopping condition, ATLBO converges well when 0.15 × n × m seconds CPU time is reached; moreover, when 0.15 × n × m seconds CPU time is applied, IMOFA, TICA, and HMOTLBO also converge well, so the above CPU time is given as the stopping condition for all the algorithms.
The Taguchi method [50] is used to decide the settings for the other parameters by using, for instance, 50 × 2 . Table 2 gives the levels of each parameter. The orthogonal array L 27 ( 3 5 ) is tested. ATLBO with each combination runs 10 times independently, for instance, 50 × 2 .
The results of the D I R and S / N ratio are given in Figure 3, where the S / N ratio is 10 l o g 10 ( D I R 2 ) . As shown in Figure 3, ATLBO with the combination N = 100 , b = 5 , μ = 0.7 , λ = 3 and β = 9 can obtain better results than ATLBO with other combinations, so the above parameter settings are used.
TLBO has two parameters N = 100 , μ = 0.7 . Experiments show that the parameter settings of IMOFA, TICA, and HMOTLBO are still effective, so they continue to be adopted.

4.3. Results and Discussion

The ATLBO, IMOFA, TICA, HMOTLBO, and TLBO algorithms are compared. Each algorithm runs randomly 10 times for each instance. Table 3, Table 4 and Table 5 describe the corresponding results of the five algorithms. A, I, T, H, TL denote ATLBO, IMOFA, TICA, HMOTLBO, and TLBO. Figure 4 shows the convergence curves of ATLBO, IMOFA, TICA, HMOTLBO, and TLBO on 30 × 2 and 50 × 2 instances. Figure 5 shows box plots of the three metrics based on Table 3, Table 4 and Table 5.
As shown in Table 3, Table 4 and Table 5, for eight small-scale instances, the three metrics of ATLBO are greater than TLBO. ATLBO obtains smaller D I R than TLBO and C ( A , T L ) is greater than C ( T L , A ) for 22 of the 24 medium-scale instances. For the large-scale instances, ATLBO obtains smaller D I R than TLBO; moreover, C ( A , T L ) and ρ are 1, that is, all the solutions of TLBO are dominated by non-dominated solutions produced by ATLBO. There are notable performance differences between ATLBO and TLBO. This conclusion can also be seen from Figure 4 and Figure 5. Box plots for all the algorithms illustrate the data in the Table 3, Table 4 and Table 5. The performance differences between ATLBO and TLBO show that the new strategies, such as an adaptive learner phase, have positive impacts on the performances of ATLBO, so the new strategies are effective.
As shown in Table 3, Table 4 and Table 5, ATLBO obtains smaller C ( I , A ) than C ( A , I ) on 68 instances and has C ( A , I ) of 1 on at least 64 instances; that is, all solutions of IMOFA are dominated by non-dominated solutions of ATLBO. ATLBO converges significantly better than IMOFA. It can be found from Table 5 that ρ of ATLBO is better than that of IMOFA on 68 instances and ρ of IMOFA is 0 on 65 instances; that is, IMOFA cannot generate any members for reference set Ω * . Table 4 shows that ATLBO obtains smaller D I R than IMOFA on 64 of 68 instances. Significant performance differences between ATLBO and IMOFA can also be seen from Figure 4 and Figure 5.
When ATLBO is compared with TICA, it is found that ATLBO performs better than TICA for most instances. Table 4 shows that ATLBO produces smaller C ( T , A ) than C ( A , T ) on all instances and obtains C ( A , T ) of 1 on 56 instances. As shown in Table 5, ATLBO produces a larger ρ than TICA on 58 instances and ρ of ATLBO is 1 on 54 instances. ATLBO also performs better than TICA on D I R because ATLBO obtains better D I R than TICA on 65 instances. The above analyses show that ATLBO can provide better results than TICA. This conclusion can also be drawn from Figure 4 and Figure 5. They depict the notable differences between all the algorithms.
It can also be concluded from Table 3, Table 4 and Table 5 that ATLBO outperforms HMOTLBO on the three metrics. ATLBO has smaller C ( H , A ) than C ( A , H ) on all instances, obains larger ρ than HMOTLBO on all instances and better D I R than HMOTLBO on 67 instances. The performance differences between ATLBO and HMOTLBO can also be seen in Figure 4 and Figure 5. Figure 5 shows the significant difference of ATLBO. It can be concluded that ATLBO has promising search abilities for solving EHFSPs with additional resources.
In ATLBO, the teacher phase consists of self-learning and teaching for each teacher, good solutions are used fully, the global ability is intensified, and an adaptive learner phase is applied. As a result, mutual learning on a low quality class is avoided and a reinforcement search of some better solutions is conducted to produce high quality solutions. These new strategies lead to a good balance between exploitation and exploration and better results than for the comparative algorithms; thus, ATLBO is a very competitive algorithm for EHFSPs with additional resources.

5. Conclusions

In this study, an EHFSP with additional resources is studied and a new algorithm called ATLBO, which has multiple teachers, is proposed to minimize the makespan and the total energy consumption. After two classes are formed, a teacher phase is first executed, which consists of the teacher’s self-learning and the teacher’s training. Then, an adaptive learner phase is presented, in which the quality of the two classes is used to decide adaptively the learner phase or the reinforcement search of the temporary solution set, and an adaptive formation of classes is also applied. The computational results obtained demonstrate that new strategies are effective and that ATLBO is a very competitive algorithm for EHFSPs with additional resources.
EHFSPs with various constraints have attracted much attention. In the future, we will consider EHFSPs with a batch processing machine in a casting process by applying new optimization strategies, such as meta-heuristic learning. TLBO with new optimization mechanisms and its application to production scheduling are also topics that will be considered in the future.

Author Contributions

Conceptualization, D.L. and H.L.; methodology, D.L. and J.Z.; software, J.Z.; validation, H.L. and D.L.; formal analysis, D.L. and H.L.; investigation, D.L. and H.L.; resources, J.Z.; data curation, J.Z.; writing—original draft and review, J.Z. and D.L.; visualization, J.Z.; supervision, D.L. and H.L.; project administration, D.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundation of China (No. 61573264).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

S l the number of machines at stage l
M l k the k-th machines at stage l
r i l k the processed job J i needs r i l k units of the additional resource on M l k
R l the total number available of additional resources at stage l
e l k the energy consumption per unit time of processing mode on M l k
i e l k the energy consumption per unit time of idle mode on M l k
C i the completion time of job J i
C m a x the maximum completion time of all jobs
T E C the total energy consumption
w i l k t decision variable, if job J i is processed on M l k at time t, w i l k t = 1;
otherwise w i l k t = 0
z i l k t decision variable, if machine M l k is free at time t, z i l k t = 1;
otherwise z i l k t = 0
E v o i the evaluation index of class C l s i

References

  1. Ruiz, R.; Vázquez-Rodríguez, J.A. The hybrid flow shop scheduling problem. Eur. J. Oper. Res. 2010, 205, 1–18. [Google Scholar] [CrossRef]
  2. Li, J.Q.; Bai, S.C.; Duan, P.Y.; Sang, H.Y.; Han, Y.Y.; Zheng, Z.X. An improved artificial bee colony algorithm for addressing distributed flow shop with distance coefficient in a prefabricated system. Int. J. Prod. Res. 2019, 57, 6922–6942. [Google Scholar] [CrossRef]
  3. Yan, J.H.; Li, L.; Zhao, F.; Zhang, F.Y.; Zhao, Q.L. A multi-level optimization approach for energy-efficient flexible flow shop scheduling. Int. J. Prod. Res. 2016, 137, 1543–1552. [Google Scholar] [CrossRef]
  4. Zeng, Z.Q.; Hong, M.N.; Man, Y.; Li, J.G.; Zhang, Y.Z.; Liu, H.B. Multi-objective optimization of flexible flow shop scheduling with batch process-consideration total electricity consumption and material wastage. J. Clean. Prod. 2018, 183, 925–939. [Google Scholar] [CrossRef]
  5. Li, M.; Lei, D.M.; Cai, J.C. Two-level imperialist competitive algorithm for energy-efficient hybrid flow shop scheduling problem with relative importance of objectives. Swarm. Evol. Comput. 2019, 49, 34–43. [Google Scholar] [CrossRef]
  6. Zuo, Y.D.; Fan, Z.; Zou, T.R.; Wang, P. A novel multi-population artificial bee colony algorithm for energy-efficient hybrid flow shop scheduling problem. Symmetry 2021, 13, 2421. [Google Scholar] [CrossRef]
  7. Song, C.L. A cooperative memetic algorithm with learning-based agent for energy-aware distributed hybrid flow-shop scheduling. IEEE Access 2021, 9, 56822–56835. [Google Scholar] [CrossRef]
  8. Yue, L.; Wang, H.; Mumtaz, J.; Rauf, M.; Li, Z.F. Energy-efficient scheduling of a two-stage flexible printed circuit board flow shop using a hybrid Pareto spider monkey optimisation algorithm. J. Ind. Inf. Integr. 2023, 31, 100412. [Google Scholar] [CrossRef]
  9. Chen, W.D.; Wang, J.N.; Yu, G.Y.; Hu, Y.M. Energy-Efficient Hybrid Flow-Shop Scheduling under Time-of-Use and Ladder Electricity Tariffs. Appl. Sci. 2022, 12, 6456. [Google Scholar] [CrossRef]
  10. Qin, H.X.; Han, Y.Y.; Zhang, B.; Meng, L.L.; Liu, Y.P.; Pan, Q.K.; Gong, D.W. An improved iterated greedy algorithm for the energy-efficient blocking hybrid flow shop scheduling problem. Swarm. Evol. Comput. 2022, 69, 100992. [Google Scholar] [CrossRef]
  11. Zuo, Y.D.; Wang, P.; Li, M. A Population Diversity-Based Artificial Bee Colony Algorithm for Assembly Hybrid Flow Shop Scheduling with Energy Consumption. Appl. Sci. 2023, 13, 10903. [Google Scholar] [CrossRef]
  12. Wang, Y.J.; Wang, G.G.; Tian, F.M.; Gong, D.W.; Pedrycz, W. Solving energy-efficient fuzzy hybrid flow-shop scheduling problem at a variable machine speed using an extended NSGA-II. Eng. Appl. Artif. Intel. 2023, 121, 105977. [Google Scholar] [CrossRef]
  13. Li, P.Z.; Xue, Q.; Zhang, Z.T.; Jian, C.; Zhou, D.Q. Multi-objective energy-efficient hybrid flow shop scheduling using Q-learning and GVNS driven NSGA-II. Comput. Oper. Res. 2023, 159, 106360. [Google Scholar] [CrossRef]
  14. Wang, Z.Y.; Shen, L.S.; Li, X.; Gao, L. An improved multi-objective firefly algorithm for energy-efficient hybrid flowshop rescheduling problem. J. Clean. Prod. 2023, 385, 135738. [Google Scholar] [CrossRef]
  15. Qin, W.; Zhuang, Z.; Liu, Y.; Tang, O. A two-stage ant colony algorithm for hybrid flow shop scheduling with lot sizing and calendar constraints in printed circuit board assembly. Comput. Ind. Eng. 2019, 138, 106115. [Google Scholar] [CrossRef]
  16. Tang, H.T.; Zhang, W.; Li, X.X.; Wei, S.P. A discrete group teaching optimization algorithm for solving many-objective sand casting whole process production scheduling problem. Comput. Oper. Res. 2024, 164, 106563. [Google Scholar] [CrossRef]
  17. Engin, O.; Güçlü, A. A new hybrid ant colony optimization algorithm for solving the no-wait flow shop scheduling problems. Appl. Soft. Comput. 2018, 72, 166–176. [Google Scholar] [CrossRef]
  18. Edis, E.B.; Oguz, C. Integrated production and distribution scheduling in distributed hybrid flow shops. Comput. Ind. Eng. 2012, 63, 433–447. [Google Scholar] [CrossRef]
  19. Zheng, X.L.; Wang, L. A two-stage adaptive fruit fly optimization algorithm for unrelated parallel machine scheduling problem with additional resource constraints. Expert. Syst. Appl. 2016, 65, 28–39. [Google Scholar] [CrossRef]
  20. Zheng, X.L.; Wang, L. A collaborative multiobjective fruit fly optimization algorithm for the resource constrained unrelated parallel machine green scheduling problem. IEEE Trans. Syst. Man Cybern. Syst. 2018, 48, 790–800. [Google Scholar] [CrossRef]
  21. Yepes-Borrero, J.C.; Villa, F.; Perea, F.; Caballero-Villalobos, J.P. GRASP algorithm for the unrelated parallel machine scheduling problem with setup times and additional resources. Expert. Syst. Appl. 2020, 141, 112959. [Google Scholar] [CrossRef]
  22. Fanjul-Peyro, L.; Perea, F.; Ruiz, R. Models and matheuristics for the unrelated parallel machine scheduling problem with additional resources. Eur. J. Oper. Res. 2017, 260, 482–493. [Google Scholar] [CrossRef]
  23. Fanjul-Peyro, L. Models and an exact method for the unrelated parallel machine scheduling problem with setups and resources. Expert Syst. Appl. 2020, 5, 100022. [Google Scholar] [CrossRef]
  24. Pinar, Y.; Seyda, T.Y. Constraint programming approach for multi-resource-constrained unrelated parallel machine scheduling problem with sequence-dependent setup times. Int. J. Prod. Res. 2022, 60, 2212–2229. [Google Scholar]
  25. Bitar, A.; Dauzère-Pérès, S.; Yugma, C. Unrelated parallel machine scheduling with new criteria: Complexity and models. Comput. Oper. Res. 2021, 132, 105291. [Google Scholar] [CrossRef]
  26. Lei, D.M.; He, S.S. An adaptive artificial bee colony for unrelated parallel machine scheduling with additional resource and maintenance. Expert Syst. Appl. 2022, 205, 117577. [Google Scholar] [CrossRef]
  27. Li, M.; Xiong, H.; Lei, D.M. An Artificial Bee Colony with Adaptive Competition for the Unrelated Parallel Machine Scheduling Problem with Additional Resources and Maintenance. Symmetry 2022, 14, 1380. [Google Scholar] [CrossRef]
  28. Figielska, E. A new heuristic for scheduling the two-stage flowshop with additional resources. Comput. Ind. Eng. 2008, 54, 750–763. [Google Scholar] [CrossRef]
  29. Figielska, E. A genetic algorithm and a simulated annealing algorithm combined with column generation technique for solving the problem of scheduling in the hybrid flowshop with additional resources. Comput. Ind. Eng. 2009, 56, 142–151. [Google Scholar] [CrossRef]
  30. Rao, R.V.; Savsani, V.J.; Vakharia, D.P. Teaching learning-based optimization: A novel method for constrained mechanical design optimization problems. Comput. Aided. Design. 2011, 43, 303–315. [Google Scholar] [CrossRef]
  31. Xie, Z.P.; Zhang, C.Y.; Shao, X.Y.; Lin, W.W.; Zhu, H.P. An effective hybrid teaching-learning-based optimization algorithm for permutation flow shop scheduling problem. Adv. Eng. Softw. 2014, 77, 35–47. [Google Scholar] [CrossRef]
  32. Lei, D.M.; Su, B. Cooperated teaching-learning-based optimisation for distributed two-stage assembly flow shop scheduling. Int. J. Prod. Res. 2021, 59, 7232–7245. [Google Scholar] [CrossRef]
  33. Shao, W.S.; Pi, D.C.; Shao, Z.S. A hybrid discrete teaching-learning-based meta-heuristic for solving no-idle flow shop scheduling problem with total tardiness criterion. Comput. Oper. Res. 2018, 94, 89–105. [Google Scholar] [CrossRef]
  34. Lei, D.M.; Gao, L.; Zhen, Y.L. A novel teaching-learning-based optimization algorithm for energy-efficient scheduling in hybrid flow shop. IEEE Trans. Eng. Manag. 2018, 65, 330–340. [Google Scholar] [CrossRef]
  35. Li, J.Q.; Pan, Q.K.; Mao, K. A discrete teachinglearning-based optimisation algorithm for realistic flow shop scheduling problem. Eng. Appl. Artif. Intell. 2015, 37, 279–292. [Google Scholar] [CrossRef]
  36. Xu, Y.; Wang, L.; Wang, S.Y.; Liu, M. An effective teaching-learning-based optimization algorithm for the flexible jobshop scheduling problem with fuzzy processing time. Neurocomputing 2015, 148, 260–268. [Google Scholar] [CrossRef]
  37. Buddala, R.; Mahapatra, S.S. Two-stage teachinglearning-based optimization method for flexible job-shop scheduling under machine breakdown. Int. J. Adv. Manuf. Tech. 2019, 100, 1419–1432. [Google Scholar] [CrossRef]
  38. Lei, D.M.; Xi, B.J. Diversified teaching-learning-based optimization for fuzzy two-stage hybrid flow shop scheduling with setup time. J. Intell. Fuzzy Syst. 2021, 41, 4159–4173. [Google Scholar] [CrossRef]
  39. Chai, R.Q.; Tsourdos, A.; Savvaris, A.; Chai, S.C.; Xia, Y.Q. Multiobjective optimal parking maneuver planning of autonomous wheeled vehicles. IEEE Trans. Ind. Electr. 2020, 67, 10809–10821. [Google Scholar] [CrossRef]
  40. Deb, K.; Pratap, A.; Agarwal, S.; Meyarivan, T. A fast and elitist multiobjective genetic algorithm: NSGA-II. IEEE Transac. Evolut. Comput. 2002, 6, 182–197. [Google Scholar] [CrossRef]
  41. Deng, J.; Wang, L.; Wang, S.Y.; Zheng, X.L. A competitive memetic algorithm for the distributed two-stage assembly flow-shop scheduling problem. Int. J. Prod. Res. 2015, 54, 3561–3577. [Google Scholar] [CrossRef]
  42. Karimi, S.; Ardalan, Z.; Naderi, B.; Mohammadi, M. Scheduling flexible job-shops with transportation times: Mathematical models and a hybrid imperialist competitive algorithm. Appl. Math. Model. 2017, 41, 667–682. [Google Scholar] [CrossRef]
  43. Zitzler, E.; Thiele, L. Multiobjective evolutionary algorithms: A comparative case study and the strength Pareto approach. IEEE Transac. Evolut. Comput. 1999, 3, 257–271. [Google Scholar] [CrossRef]
  44. Lei, D.M. Pareto archive particle swarm optimization for multi-objective fuzzy job shop scheduling problems. Int. J. Adv. Manuf. Tech. 2008, 37, 157–165. [Google Scholar] [CrossRef]
  45. Knowles, J.; Corne, D. On metrics for comparing nondominated sets. In Proceedings of the 2002 Congress on Evolutionary Computation, Honolulu, HI, USA, 12–17 May 2002. [Google Scholar]
  46. Khunkitti, S.; Siritaratiwat, A.; Premrudeepreechacharn, S. A Many-Objective Marine Predators Algorithm for Solving Many-Objective Optimal Power Flow Problem. Appl. Sci. 2022, 12, 11829. [Google Scholar] [CrossRef]
  47. Ezugwu, A.E.; Agushaka, J.O.; Abualigah, L.; Mirjalili, S.; Gandomi, A.H. Prairie Dog Optimization Algorithm. Neural. Comput. Appl. 2022, 34, 20017–20065. [Google Scholar] [CrossRef]
  48. Jia, H.M.; Rao, H.H.; Wen, C.S.; Mirjalili, S. Crayfish optimization algorithm. Artif. Intell. Rev. 2023, 56, 1919–1979. [Google Scholar] [CrossRef]
  49. Degli-Esposti, V.; Fuschini, F.; Bertoni, H.L.; Thomä, R.S.; Kürner, T.; Yin, X.F.; Guan, K. IEEE Access Special Section Editorial: Millimeter-Wave and Terahertz Propagation, Channel Modeling, and Applications. IEEE Access 2021, 9, 67660–67666. [Google Scholar] [CrossRef]
  50. Taguchi, G. Introduction to Quality Engineering; Asian Productivity Organization: Tokyo, Japan, 1986. [Google Scholar]
Figure 1. A schedule of the example.
Figure 1. A schedule of the example.
Symmetry 16 00203 g001
Figure 2. Flowchart of ATLBO.
Figure 2. Flowchart of ATLBO.
Symmetry 16 00203 g002
Figure 3. The mean D I R and the mean S / N ratio of D I R .
Figure 3. The mean D I R and the mean S / N ratio of D I R .
Symmetry 16 00203 g003
Figure 4. Distribution of non-dominated solutions of five algorithms.
Figure 4. Distribution of non-dominated solutions of five algorithms.
Symmetry 16 00203 g004
Figure 5. Box plots of five algorithms.
Figure 5. Box plots of five algorithms.
Symmetry 16 00203 g005
Table 1. Processing time of the example.
Table 1. Processing time of the example.
J i J 1 J 2 J 3 J 4 J 5 J 6 J 7 J 8
M 11 21(2)25(5)115(8)20(5)110(9)88(7)16(1)27(2)
M 12 26(6)36(7)111(3)17(4)112(6)89(8)22(5)17(3)
M 21 54(9)51(2)44(3)60(5)41(5)49(5)56(1)47(7)
M 22 61(7)70(5)77(9)64(7)59(4)70(5)76(4)73(4)
Table 2. Parameters and their levels.
Table 2. Parameters and their levels.
ParametersFactor Level
123
N80100120
b357
μ 0.650.70.75
λ 234
β 8.599.5
Table 3. Results of ATLBO, IMOFA, TICA, HMOTLBO, and TLBO on D I R .
Table 3. Results of ATLBO, IMOFA, TICA, HMOTLBO, and TLBO on D I R .
InstanceATLBOIMOFATICAHMOTLBOTLBOInstanceATLBOIMOFATICAHMOTLBOTLBO
8 × 2 × 10.1090.6150.4290.4360.553150 × 2 × 30.0000.9920.6880.6440.732
8 × 2 × 20.1410.2880.6230.2120.359150 × 4 × 10.0000.0730.7220.6940.951
8 × 4 × 10.0000.3660.5860.7321.098150 × 4 × 20.0000.6090.4330.6180.997
8 × 4 × 20.2810.3630.4990.4410.547150 × 4 × 30.0000.5020.3920.5600.705
8 × 6 × 10.0670.3550.5850.2100.381150 × 6 × 10.0000.9420.4410.5150.655
8 × 6 × 20.1440.2430.6030.3100.468150 × 6 × 20.0000.8640.5600.6690.698
8 × 8 × 10.1450.2860.2710.2820.637150 × 6 × 30.0000.5890.4510.6501.078
8 × 8 × 20.1310.0450.3550.5450.533150 × 8 × 10.0000.5860.7140.7021.088
16 × 2 × 10.0000.5000.3780.5400.720150 × 8 × 20.0001.0320.3800.4770.624
16 × 2 × 20.2830.5590.0950.2520.404150 × 8 × 30.0001.0080.6060.6940.733
16 × 4 × 10.0000.4510.6540.9150.801250 × 2 × 10.0000.4920.3310.5721.205
16 × 4 × 20.6140.4100.1520.3360.518250 × 2 × 20.0000.5920.2960.8420.826
16 × 6 × 10.0000.3800.4190.4630.577250 × 2 × 30.0001.3480.4750.4180.580
16 × 6 × 20.0000.5450.4380.5420.725250 × 4 × 10.0001.0190.6430.9460.784
16 × 8 × 10.2220.1500.5620.4160.301250 × 4 × 20.0000.5480.4370.6231.025
16 × 8 × 20.1200.1590.6560.2590.418250 × 4 × 30.0000.5920.5940.7290.972
30 × 2 × 10.0000.6440.3200.4220.676250 × 6 × 10.0000.9600.5490.9170.655
30 × 2 × 20.0000.3600.5490.7290.994250 × 6 × 20.0001.1110.5100.4580.575
30 × 4 × 10.4450.5290.2780.5150.747250 × 6 × 30.0000.6700.2231.0050.980
30 × 4 × 20.7430.4780.4680.5160.587250 × 8 × 10.0000.5010.5970.9721.041
30 × 6 × 10.4600.3690.3350.2880.465250 × 8 × 20.0000.6720.4990.7561.008
30 × 6 × 20.0000.6100.5570.7150.930250 × 8 × 30.0001.0460.4660.6310.853
30 × 8 × 10.0000.4580.5250.6350.881350 × 2 × 10.0000.5780.3590.8480.976
30 × 8 × 20.0000.5850.3580.5700.763350 × 2 × 20.0001.0850.6480.8630.980
50 × 2 × 10.0000.6420.5200.6820.857350 × 2 × 30.0001.0400.6100.6780.907
50 × 2 × 20.0000.6620.5700.8120.880350 × 4 × 10.0000.5900.4930.5200.425
50 × 4 × 10.0000.8690.4320.4630.625350 × 4 × 20.0001.0640.8341.4041.414
50 × 4 × 20.0000.8870.3870.5110.487350 × 4 × 30.0000.6860.5150.8571.030
50 × 6 × 10.0000.5520.4130.6530.640350 × 6 × 10.0000.5090.8320.5240.641
50 × 6 × 20.0000.5530.4940.5530.637350 × 6 × 20.0000.8980.5170.7300.677
50 × 8 × 10.0000.8430.2820.5800.593350 × 6 × 30.0000.7110.3380.8621.013
50 × 8 × 20.0000.2990.4720.5350.788350 × 8 × 10.0000.7360.5890.8841.031
150 × 2 × 10.0000.8000.4700.4310.984350 × 8 × 20.0000.7510.3190.8770.541
150 × 2 × 20.0000.1930.1840.6811.042350 × 8 × 30.0000.8220.4601.0560.822
Table 4. Results of ATLBO, IMOFA, TICA, HMOTLBO, and TLBO on C .
Table 4. Results of ATLBO, IMOFA, TICA, HMOTLBO, and TLBO on C .
Instance C ( A , I ) C ( I , A ) C ( A , T ) C ( T , A ) C ( A , H ) C ( H , A ) C ( A , TL ) C ( TL , A )
8 × 2 × 11.0000.0000.5710.0001.0000.0001.0000.000
8 × 2 × 21.0000.0000.7000.0000.9000.0001.0000.000
8 × 4 × 11.0000.0001.0000.0001.0000.0001.0000.000
8 × 4 × 21.0000.0000.8750.0001.0000.0001.0000.000
8 × 6 × 11.0000.0000.5000.0000.6000.0001.0000.000
8 × 6 × 21.0000.0000.8000.0001.0000.0001.0000.000
8 × 8 × 11.0000.0000.6000.0000.9000.0001.0000.000
8 × 8 × 20.4290.0000.5000.0001.0000.0001.0000.000
16 × 2 × 11.0000.0001.0000.0001.0000.0001.0000.000
16 × 2 × 21.0000.0000.5000.0001.0000.0001.0000.000
16 × 4 × 11.0000.0001.0000.0001.0000.0001.0000.000
16 × 4 × 21.0000.0000.7780.0001.0000.0001.0000.000
16 × 6 × 11.0000.0001.0000.0001.0000.0001.0000.000
16 × 6 × 21.0000.0001.0000.0001.0000.0001.0000.000
16 × 8 × 10.3330.0001.0000.0001.0000.0001.0000.000
16 × 8 × 20.5000.0000.8000.0001.0000.0001.0000.000
30 × 2 × 11.0000.0001.0000.0001.0000.0001.0000.000
30 × 2 × 21.0000.0001.0000.0001.0000.0001.0000.000
30 × 4 × 11.0000.0000.8570.0001.0000.0001.0000.000
30 × 4 × 20.0001.0000.3000.0000.4000.0000.5000.000
30 × 6 × 11.0000.0001.0000.0000.9000.0001.0000.000
30 × 6 × 21.0000.0001.0000.0001.0000.0001.0000.000
30 × 8 × 11.0000.0001.0000.0001.0000.0001.0000.000
30 × 8 × 21.0000.0001.0000.0001.0000.0001.0000.000
50 × 2 × 11.0000.0001.0000.0001.0000.0001.0000.000
50 × 2 × 21.0000.0001.0000.0001.0000.0001.0000.000
50 × 4 × 11.0000.0001.0000.0001.0000.0001.0000.000
50 × 4 × 21.0000.0001.0000.0001.0000.0001.0000.000
50 × 6 × 11.0000.0001.0000.0001.0000.0001.0000.000
50 × 6 × 21.0000.0001.0000.0001.0000.0001.0000.000
50 × 8 × 11.0000.0001.0000.0001.0000.0001.0000.000
50 × 8 × 21.0000.0001.0000.0001.0000.0001.0000.000
150 × 2 × 11.0000.0001.0000.0001.0000.0001.0000.000
150 × 2 × 21.0000.0001.0000.0001.0000.0001.0000.000
150 × 2 × 31.0000.0001.0000.0001.0000.0001.0000.000
150 × 4 × 11.0000.0001.0000.0001.0000.0001.0000.000
150 × 4 × 21.0000.0001.0000.0001.0000.0001.0000.000
150 × 4 × 31.0000.0001.0000.0001.0000.0001.0000.000
150 × 6 × 11.0000.0001.0000.0001.0000.0001.0000.000
150 × 6 × 21.0000.0001.0000.0001.0000.0001.0000.000
150 × 6 × 31.0000.0001.0000.0001.0000.0001.0000.000
150 × 8 × 11.0000.0001.0000.0001.0000.0001.0000.000
150 × 8 × 21.0000.0001.0000.0001.0000.0001.0000.000
150 × 8 × 31.0000.0001.0000.0001.0000.0001.0000.000
250 × 2 × 11.0000.0001.0000.0001.0000.0001.0000.000
250 × 2 × 21.0000.0001.0000.0001.0000.0001.0000.000
250 × 2 × 31.0000.0001.0000.0001.0000.0001.0000.000
250 × 4 × 11.0000.0001.0000.0001.0000.0001.0000.000
250 × 4 × 21.0000.0001.0000.0001.0000.0001.0000.000
250 × 4 × 31.0000.0001.0000.0001.0000.0001.0000.000
250 × 6 × 11.0000.0001.0000.0001.0000.0001.0000.000
250 × 6 × 21.0000.0001.0000.0001.0000.0001.0000.000
250 × 6 × 31.0000.0001.0000.0001.0000.0001.0000.000
250 × 8 × 11.0000.0001.0000.0001.0000.0001.0000.000
250 × 8 × 21.0000.0001.0000.0001.0000.0001.0000.000
250 × 8 × 31.0000.0001.0000.0001.0000.0001.0000.000
350 × 2 × 11.0000.0001.0000.0001.0000.0001.0000.000
350 × 2 × 21.0000.0001.0000.0001.0000.0001.0000.000
350 × 2 × 31.0000.0001.0000.0001.0000.0001.0000.000
350 × 4 × 11.0000.0001.0000.0001.0000.0001.0000.000
350 × 4 × 21.0000.0001.0000.0001.0000.0001.0000.000
350 × 4 × 31.0000.0001.0000.0001.0000.0001.0000.000
350 × 6 × 11.0000.0001.0000.0001.0000.0001.0000.000
350 × 6 × 21.0000.0001.0000.0001.0000.0001.0000.000
350 × 6 × 31.0000.0001.0000.0001.0000.0001.0000.000
350 × 8 × 11.0000.0001.0000.0001.0000.0001.0000.000
350 × 8 × 21.0000.0001.0000.0001.0000.0001.0000.000
350 × 8 × 31.0000.0001.0000.0001.0000.0001.0000.000
Table 5. Results of ATLBO, IMOFA, TICA, HMOTLBO, and TLBO on ρ .
Table 5. Results of ATLBO, IMOFA, TICA, HMOTLBO, and TLBO on ρ .
InstanceATLBOIMOFATICAHMOTLBOTLBOInstanceATLBOIMOFATICAHMOTLBOTLBO
8 × 2 × 10.250.000.750.000.00150 × 2 × 31.000.000.000.000.00
8 × 2 × 20.250.000.750.000.00150 × 4 × 11.000.000.000.000.00
8 × 4 × 11.000.000.000.000.00150 × 4 × 21.000.000.000.000.00
8 × 4 × 20.500.000.500.000.00150 × 4 × 31.000.000.000.000.00
8 × 6 × 10.170.000.830.000.00150 × 6 × 11.000.000.000.000.00
8 × 6 × 20.330.000.670.000.00150 × 6 × 21.000.000.000.000.00
8 × 8 × 10.170.000.670.170.00150 × 6 × 31.000.000.000.000.00
8 × 8 × 20.200.800.000.000.00150 × 8 × 11.000.000.000.000.00
16 × 2 × 11.000.000.000.000.00150 × 8 × 21.000.000.000.000.00
16 × 2 × 20.330.000.670.000.00150 × 8 × 31.000.000.000.000.00
16 × 4 × 11.000.000.000.000.00250 × 2 × 11.000.000.000.000.00
16 × 4 × 20.230.000.770.000.00250 × 2 × 21.000.000.000.000.00
16 × 6 × 11.000.000.000.000.00250 × 2 × 31.000.000.000.000.00
16 × 6 × 21.000.000.000.000.00250 × 4 × 11.000.000.000.000.00
16 × 8 × 10.330.670.000.000.00250 × 4 × 21.000.000.000.000.00
16 × 8 × 20.250.250.500.000.00250 × 4 × 31.000.000.000.000.00
30 × 2 × 11.000.000.000.000.00250 × 6 × 11.000.000.000.000.00
30 × 2 × 21.000.000.000.000.00250 × 6 × 21.000.000.000.000.00
30 × 4 × 10.500.000.500.000.00250 × 6 × 31.000.000.000.000.00
30 × 4 × 20.000.500.500.000.00250 × 8 × 11.000.000.000.000.00
30 × 6 × 10.500.000.000.500.00250 × 8 × 21.000.000.000.000.00
30 × 6 × 21.000.000.000.000.00250 × 8 × 31.000.000.000.000.00
30 × 8 × 11.000.000.000.000.00350 × 2 × 11.000.000.000.000.00
30 × 8 × 21.000.000.000.000.00350 × 2 × 21.000.000.000.000.00
50 × 2 × 11.000.000.000.000.00350 × 2 × 31.000.000.000.000.00
50 × 2 × 21.000.000.000.000.00350 × 4 × 11.000.000.000.000.00
50 × 4 × 11.000.000.000.000.00350 × 4 × 21.000.000.000.000.00
50 × 4 × 21.000.000.000.000.00350 × 4 × 31.000.000.000.000.00
50 × 6 × 11.000.000.000.000.00350 × 6 × 11.000.000.000.000.00
50 × 6 × 21.000.000.000.000.00350 × 6 × 21.000.000.000.000.00
50 × 8 × 11.000.000.000.000.00350 × 6 × 31.000.000.000.000.00
50 × 8 × 21.000.000.000.000.00350 × 8 × 11.000.000.000.000.00
150 × 2 × 11.000.000.000.000.00350 × 8 × 21.000.000.000.000.00
150 × 2 × 21.000.000.000.000.00350 × 8 × 31.000.000.000.000.00
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lei, D.; Zhang, J.; Liu, H. An Adaptive Two-Class Teaching-Learning-Based Optimization for Energy-Efficient Hybrid Flow Shop Scheduling Problems with Additional Resources. Symmetry 2024, 16, 203. https://doi.org/10.3390/sym16020203

AMA Style

Lei D, Zhang J, Liu H. An Adaptive Two-Class Teaching-Learning-Based Optimization for Energy-Efficient Hybrid Flow Shop Scheduling Problems with Additional Resources. Symmetry. 2024; 16(2):203. https://doi.org/10.3390/sym16020203

Chicago/Turabian Style

Lei, Deming, Jiawei Zhang, and Hongli Liu. 2024. "An Adaptive Two-Class Teaching-Learning-Based Optimization for Energy-Efficient Hybrid Flow Shop Scheduling Problems with Additional Resources" Symmetry 16, no. 2: 203. https://doi.org/10.3390/sym16020203

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