Bi-Objective Re-Entrant Hybrid Flow Shop Scheduling considering Energy Consumption Cost under Time-of-Use Electricity Tariffs

Re-entrant hybrid flow shop scheduling problem (RHFSP) is widely used in industries. However, little attention is paid to energy consumption cost with the raise of green manufacturing concept. )is paper proposes an improved multiobjective ant lion optimization (IMOALO) algorithm to solve the RHFSP with the objectives of minimizing the makespan and energy consumption cost under Time-of-Use (TOU) electricity tariffs. A right-shift operation is then used to adjust the starting time of operations by avoiding the period of high electricity price to reduce the energy consumption cost as far as possible. )e experimental results show that IMOALO algorithm is superior to multiobjective ant lion optimization (MOALO) algorithm, NSGA-II, andMOPSO in terms of the convergence, dominance, and diversity of nondominated solutions. )e proposed model can make enterprises avoid high price period reasonably, transfer power load, and reduce the energy consumption cost effectively. Meanwhile, parameter analysis indicates that the period of TOU electricity tariffs and energy efficiency of machines have great impact on the scheduling results.


Introduction
Re-entrant hybrid flow shop scheduling problem (RHFSP) is a combination of classic hybrid flow shop scheduling problem and re-entrant scheduling problem. It means that all jobs have the same processing path among the stages and go back and forth in the same order many times. It is widely used in some special industries, such as thin-film transistorliquid crystal display (TFT-LCD) panel manufacturing, printed circuit board (PCB), and semiconductor wafer manufacturing. With the rising energy prices and the increasingly serious environmental problems, the issue of reducing energy consumption cost has become a focus for many factory managers. Recently, many power energy suppliers have begun to implement a so-called Time-of-Use (TOU) electricity tariffs, such pricing mechanism represents a huge opportunity to reduce costs for consumers by shifting electricity usage from on-peak period to other period. RHFSP is a typical NP-hard problem [1], which is difficult to solve with exact algorithms, so it will be of academic significance and engineering application value to develop efficient intelligent optimization algorithms for the RHFSP under TOU electricity tariffs.
Since Graves et al. [2] first proposed reentrant scheduling problem (RSP) in 1983, many scholars have studied this field. Bertel and Billaut [3] applied a genetic algorithm for the RHFSP aiming at minimizing the weighted number of delayed jobs. Chen et al. [4] proposed a hybrid tabu search algorithm for the RHFSP with the objective of minimizing the makespan. Choi and Kim [5] proposed several heuristic algorithms for RHFSP, such as improved Nawaz Enscore and Ham (NEH) algorithm, aiming at minimizing the makespan. Kim and Lee [6] studied the RHFSP considering unrelated parallel machines with the objective of minimizing the makespan and total tardiness. In practice, the total tardiness objective is treated as constraints, so the problem is converted into single objective. El-Khouly et al. [7] used Lagrange decomposition to optimize the RSP aiming at minimizing the total tardiness. Wu et al. [8] studied RSP with learning effects, aiming at minimizing the makespan.
Cho et al. [9] proposed a genetic algorithm for the RHFSP in order to minimise the makespan and total tardiness. Ying et al. [10] used an iterative Pareto greedy algorithm for the bi-objective RHFSP based on the research in [9]. Shen et al. [11,12] applied an improved teaching and learning optimization algorithm and a Pareto-based discrete harmony search algorithm to solve the bi-objective RHFSP. Cheng and Lei [13] studied an improved imperialist competition algorithm for the RSP in order to minimise the makespan and total tardiness. In conclusion, production decision makers assume that electricity prices are constant in the production cycle, so they are only interested in time-related objectives, such as maximum tardiness, total tardiness, and makespan.
With the rising international energy costs and the increasingly serious environmental problems in recent years, it is particularly important to achieve green and sustainable development in the manufacturing sector. Energy consumption in the industrial sector is approximately 50% of the total energy consumption [14]. In Germany, the manufacturing enterprises consume around more than 47% of the total energy consumption. In China, the manufacturing industries are responsible for 50% of the national electricity energy each year [15]. erefore, the increasing energy price and the current sustainable development trend bring new pressure to manufacturing enterprises. It must be noted, however, that either adopting advanced production technology or purchasing new equipment needs a lot of investment. Under TOU electricity tariffs, without investing a large amount of money in equipment and technology, it can reduce the energy costs just by adjusting production tasks reasonably, and the benefits could be enjoyed easily by majority of SMEs. Zhang et al. [16] proposed an integer programming model to solve the flow shop scheduling problem under TOU electricity tariffs without affecting the production efficiency in order to minimise the carbon dioxide emissions and electricity cost. Luo et al. [17] addressed the multiobjective ant colony optimization algorithm to solve the hybrid flow shop scheduling problem with unrelated parallel machines under TOU electricity tariffs aiming at minimizing the makespan and energy consumption cost. Wang and Li [18] studied the manufacturing model considering both energy consumption and demand under TOU electricity tariffs. Fang et al. [19] used a mixed integer programming method for the flow shop scheduling problem with peak power constraints. Moon and Park [20] proposed two discrete-time mathematical models for the flexible job shop scheduling problem under TOU electricity tariffs with the objective of weighted sum of maximum completion time and power cost. Che et al. [21] proposed the mixed integer programming model for the single machine scheduling problem with the objective of minimizing the electricity cost under TOU electricity tariffs. Mikhaylidi et al. [22] studied the production and operation scheduling problem of rechargeable batteries under TOU electricity tariffs using a dynamic programming algorithm, aiming at minimizing total power consumption and delaying penalty cost. Shrouf et al. [23] proposed a discretetime integer programming model and genetic algorithm to solve the single machine scheduling problem with power outage mechanism under TOU electricity tariffs. e objective is to minimise the total power cost. Gong et al. [24] modeled the same problem with arbitrary job processing sequence and demonstrated how to reduce power cost under real-time price, TOU electricity price, and critical peak price schemes through a practical case of a surface grinder. In addition, other scholars have also studied the scheduling problems under TOU electricity tariffs in different fields (e.g., Tan et al. [25]; Castro et al. [26]; Sharma et al. [27]; Tan and Liu [28]). Although some achievements have been made in this field, the research on green job shop scheduling under TOU electricity tariffs is still immature, and the study on RHFSP with unrelated parallel machines considering energy consumption cost is even less.
is paper proposes the improved multiobjective ant lion optimization (IMOALO) algorithm with right-shift operation to approximate the Pareto optimal solutions for RHFSP under TOU electricity tariffs with the objective of minimizing the makespan and energy consumption cost. On the premise of not affecting the production efficiency, according to the characteristics of TOU electricity tariffs, the operations are arranged in the period with low electricity price as far as possible. rough this method, enterprises can reduce a large number of energy consumption cost and improve the competitiveness of enterprises, so as to realize the green and sustainable development of economy. e rest of the paper is organized as follows: the RHFSP is described and the mathematical model is established in Section 2. In Section 3, the proposed IMOALO algorithm is introduced in detail. In Section 4, the performance comparisons with other three algorithms are shown. e parameter analysis of the period of TOU electricity tariffs and energy efficiency of machines are studied in Section 5. Finally, some conclusions and future work are given in Section 6.

Problem Statement and Mathematical Model
2.1. Problem Statement. In this paper, the RHFSP with unrelated parallel machine is studied, which can be described as follows: n jobs need to be processed in s serial stages. ere is at least one stage, in which the number of unrelated parallel machines is more than one. Each job can be processed on any machine at the corresponding stage, and part of jobs may visit some stages more than once, as shown in Figure 1. e objective is to allocate all jobs to the machines and fix the starting and ending time for each operation.
Additionally, the following assumptions are made: All machines and jobs are ready at zero time. At any time, each machine can process at most one job, and each job can only be processed by one machine. e number of jobs and the processing time of all operations are given in advance. e number of stages and the number of unrelated parallel machines at each stage are known in advance. e power and speed of each machine are given in advance. All operations of each job have sequence constraints, and all jobs have no influence on each other. e power and speed of unrelated parallel machines at each stage are different. e total energy consumption of each job is independent of the machine 2 Complexity chosen. Preemption is not allowed, and once the job is processed, it cannot be interrupted. e buffer capacity between any two continuous stages is infinite, regardless of machine failure and adjustment time. 1, if the machine q is working at time t, 0, otherwise, ∀q On the basis of literature studies [10,29], a bi-objective RHFSP mathematical model is proposed in this paper. Formulas (1) and (2) represent the two objective functions. Constraint (3) implies that the starting time of the operation O jk+1 is not earlier than the ending time of O jk . Constraint (4) ensures that each operation can only be processed on one machine at the corresponding stage. Constraints (5)- (7) ensure that each machine processes at most one operation simultaneously. Constraints (8) and (9) specify the starting and ending time of O jk . Constraints (10) and (11) describe the maximum completion time.
e first objective of this paper is to minimise the makespan (maximum completion time), as shown in the following formula: Here, C max denotes the makespan. e second objective of this paper is to minimize the total energy consumption cost, which mainly includes the energy consumption cost in machine processing state and in idle state, as shown in the following formula: Here, TEC min denotes the total energy consumption cost.
e two objectives are subject to Complexity S jk ≥ 0, ∀j, k,

The Proposed Improved Multiobjective Ant Lion Optimization Algorithm
"No Free Lunch eorem" [30] points out that there is no algorithm that can solve all optimization problems. In this paper, we use the ant lion optimization ( [32], of which the search mechanism is very similar to ALO. At present, the ALO algorithm has been widely used in power grid and power optimization [33][34][35], optimal power flow optimization [36], link state routing protocol optimization [37], feature selection [38], and integrated process planning and scheduling [39]. To get better optimization results, the MOALO algorithm is improved from three aspects: (1) If the initial population is generated by a completely random method, its solutions may be concentrated in a local range, which is not conducive to convergence to the global optimal solution. is paper adopted the Latin hypercube sampling technology to initialize the population. (2) e adaptive elite ant lion updating strategy is used to improve the exploring ability in the early stage and the exploiting ability in the later stage of MOALO. (3) e local search ability of MOALO is improved by two neighbourhood structures. e flowchart for the proposed algorithm is shown in Figure 2.

Encoding and Decoding.
ALO algorithm is mainly used to solve continuous optimization problems, but rarely to deal with combination optimization problems. Random key ascending sequence coding is used in this paper to construct the mapping from individual position to job sequencing and then using the decoding method in [9,10] to decode it. Finally, scheduling schemes are generated according to various constraints. e mathematical model contains a large number of 0-1 variables. If the minimum time unit is too small, it will seriously affect the calculation speed of the model. If the minimum time unit is too large, the error will be relatively large. In this paper, the minimum time unit is set to one tenth of an hour. e decoding process is detailed as follows: Step 1: select the first gene i (i.e., job i) from the chromosome, arrange all the operations of job i on the machine that can finish it as early as possible, and record the starting time and ending time of each operation Step 2: select the next gene i′ and arrange machines for it and obtain the available machine set MS for operation O i′j Step 3: select machine q from MS, obtain all idle time periods [MStart, MEnd] of machine q, and traverse all idle time periods of machine q, and the earliest starting time t i′j of operation O i′j is shown in the following formula: Step 4: find a suitable insertion point for O i′j according to formula (13). If no idle time period satisfying the Step 5: traverse all machines in MS, repeat step 3 and step 4 to select the minimum t i′j : Step 6: calculate the ending time of operation O i′j , as shown in the following formula: Step 7: repeat step 2-step 6 until all genes in the chromosome have been processed.
Taking an example, in which there are 3 stages and 4 jobs. e number of unrelated identical parallel machines is 3, 2, and 2 at each stage, respectively. e standard processing time of each operation on the corresponding stage is shown in Table 1. In addition, the power and speed of machines at each stage are different, as shown in Table 2. For example, at stage 1, the processing speed and power of M 1,1 , M 2,1 , M 3,1 are 1, 2, and 3 and 10, 20, and 30, respectively. e idle power of all machines is 2.
In this paper, according to [14], TOU electricity tariffs function (yuan/Kwh) is as follows: If the processing time P jk � 0, it means that the job is not processed in a certain pass. Using random key coding, the permutation of all jobs represents each individual, and each element in the individual is selected in [0, 1] arbitrarily. e corresponding permutation can be obtained by ascending the ranking of each element. To get a feasible schedule, the decoding method is used to determine the processing order of all jobs at each stage and a suitable machine is assigned for each operation. en, the starting and ending time of each operation can be obtained. At last, the objective function values can be calculated. Taking the individual [0.8147, 0.1270, 0.9058, 0.9134] as an example, the corresponding permutation is 2-1-3-4 according to the ascending order. e detailed process of decoding is as follows. Firstly, all the operations of job 2 are arranged on the machine that can finish it at the earliest time. en, each operation of job 1 will be arranged on the machine that can finish it as early as possible. If the ending time of a certain operation of the job 2 is smaller than the starting time of the arranged operation of job 4 on a certain machine, the operation of the job 2 will be arranged before it, otherwise, it will be arranged behind it. By analogy, all operations of job 3 and 4 are arranged on the right machines, the Gantt charts are shown in Figure 3.
en, C max � 7.8333 hour and TEC min � 184.1634 yuan are obtained.

Population Initialization.
Generally, the region of the optimal solution is difficult to predict, and the solution space characteristics of the initial population can represent the information of all individuals to a certain extent, so the distribution performance of the initial population directly affects the convergence of the algorithm. In this paper, Latin hypercube sampling technique is used to initialize the population.
Assuming that m samples need to be extracted in n-dimensional vector space, the specific steps of Latin hypercube sampling are as follows: (1) Each dimension vector is divided into m intervals and they do not overlap with each other (2) e point is randomly selected in all intervals of each dimension (3) e selected points in (2) are randomly extracted from each dimension to form a new vector

Adaptive Elite Ant Lion Updating Strategy.
According to the elite strategy, the ALO algorithm updates the ant position according to the following formula: Among them, Ant t j is the current position of ant and R t A and R t E are ants selected by the roulette wheel selection and ants around the elite ant lion in generation t, respectively. However, in general, the initial stage of the algorithm should have a strong exploring ability, the update step is larger, and the later stage of the algorithm should have a strong exploiting ability, the update step is smaller, so the adaptive elite ant lion update strategy is introduced to update the ant position in this paper, as shown in the following formula: where f(t) � (1 − (t/T)) 2 is the adaptive coefficient and T and t represent the maximum number and current number of iterations. At the beginning of iterations, f(t) is about 1 and the value of Ant t j is mainly determined by R t A , which can ensure that the algorithm has strong exploratory ability and global search ability in the early stage. As the number of iterations increases, f(t) gets smaller and smaller. At the later stage of iterations, the value of Ant t j is mainly determined by R t E , and ants mainly walk around elite ant lions, which ensures that the algorithm has strong exploiting ability in the later stage. In summary, ants can choose different walking modes adaptively, which help to improve the exploring and exploiting ability of the algorithm.

Neighbourhood Structure.
In this paper, insertion neighbourhood and exchange neighbourhood are used to improve the local search ability of MOALO. Each individual has a 50% probability to perform insertion or exchange operations, respectively.

Complexity
Insertion neighbourhood: first of all, a new sequence π with n − 1 jobs is obtained by deleting one job from the individual randomly. en, the job j is inserted into n − 1 intervals of sequence π, and the objective function values fitness(i) of n − 1 new individuals are calculated, i ∈ [1, n − 1]. e corresponding objective function values of the original sequence is fitness(n). en, the dominance is judged and the optimal solution best is assumed to fitness(1). If fitness(i) ≺ best, then best � fitness(i), and so on; the individual corresponding to best is taken as the optimal insertion sequence, as shown in Figure 4(a).
Exchange neighbourhood: firstly, a job j is randomly selected from the individual, and the job j is exchanged with the remaining n − 1 jobs, respectively. en, the objective function values fitness(i) of the individual after exchange is calculated, i ∈ [1, n − 1], and the corresponding objective function values of the original sequence is fitness(n). Finally, the dominance is judged and the optimal solution best is assumed to fitness(1). If fitness(i) ≺ best, then best � fitness(i), and so on; the individual corresponding to best is taken as the optimal exchange sequence, as shown in Figure 4(b).

Right-Shift Operation Based on TOU Electricity Tariffs.
Under TOU electricity tariffs, the RHFSP considering energy consumption cost should not only select the appropriate machines for all jobs but also determine the starting time and ending time for each operation reasonably. Due to the complexity of RHFSP, it is inevitable that there will be waiting time for machines and jobs, especially in the case of large number of jobs, so the waiting time can be fully utilized to adjust the processing time of each operation. A right-shift operation is then added in the decoding process, and the time to minimise the energy consumption cost of the operation is selected from the possible starting time as the appropriate starting time. If there are more than one qualified starting time, the maximum value is chosen as the starting time to ensure that there is enough right-shift space. e right-shift operation does not change the makespan and job sequence but helps to reduce the energy consumption cost. Due to the limitation of the decoding method, each operation is decoded according to the earliest ending time, so there is no space to move left. It can be seen that the adjustment of the latter operation will affect the adjustment of the former one, so the right-shift operation should follow the rules from back to front in turn. Firstly, all operations are arranged in a nonincremental order according to the ending time, and then the right-shift operation is carried out in this order. Taking the operation O jk as an example, the detailed adjustment process is as follows.
Assuming that the operation O jk is processed on machine M il , the index of machine M il is q and the adjustment range of S jk is t ∈ [t min , t max ]. So without affecting the production efficiency, the job j can choose any time in this range as the starting time.
(1) Suppose that O jk is the last operation on machine M il and k � N j . If E jk = C max , then the operation does not need to move right. If E jk < C max , then t min � S jk , t max � C max − P jk /V il , and a new idle period will be added in this case; the adjustment range is shown in Figure 5. erefore, the additional idle energy consumption cost should be considered. e starting time of operation O jk after adjustment is (2) Suppose that O jk is the last operation of job j(k � N j ), but not the last operation of machine M il . If the next operation on machine M il is O j′k′ , then t min � S jk and t max � (S j′k′ − P jk /V il ); the adjustment range is shown in Figure 6. e starting time of operation O jk after adjustment is (3) Suppose that O jk is the last operation on machine M il , job j has the next operation O jk+1 and O jk+1 is processed on machine M i′l′ . e adjustment range is shown in Figure 7; meanwhile, t min � S jk and t max � (S jk+1 − P jk )/V il . e starting time of operation O jk after adjustment is (4) In other cases, t min � S jk and t max � min S j′k′ , S jk+1 − P jk /V il . If S j′k′ ≤ S jk+1 , the adjustment range is shown in Figure 8. If S j′k′ > S jk+1 , the adjustment range is shown in Figure 9. Adjustment of starting time of process O jk can be divided into two situations: (a) If O jk is the first operation on M il , part of idle time will be removed, so the removal idle energy consumption cost should be considered, then

Computational Results
In order to evaluate the effectiveness of the designed IMOALO algorithm in solving RHFSP, three other multiobjective optimization algorithms, MOALO, MOPSO, and NSGA-II, are selected for comparative study. e experimental environment is Windows 8, Intel Core i7-7700 CPU@3.60 GHz, 8G memory, and the computer simulation program is written with Matlab R2017a. In addition, to objectively compare the performance of the algorithms, the common parameters of various algorithms are taken the same value.

Description of Test Data.
is paper chooses the smallsized test problems by Cho et al. [9] as benchmark instances to verify the performance of the proposed IMOALO algorithm for solving bi-objective RHFSP. Since there are no parameters about the speed and power of machines in the instances, it is assumed that the faster the machine processing speed is, the more the power will be. In addition, the energy consumption of the same job is fixed and independent of the selected machine. is paper selects six small-sized test problems randomly for testing, and Table 3 shows the range of parameter values. For example, there are two unrelated parallel machines at a stage, numbered 1 and 2, corresponding speed is 1 and 2, and corresponding power is 10 and 20, respectively.
In this paper, a three-period TOU electricity tariffs function with a total period of 24 hours is used. e expression of the function is shown in formula (23), and the graph of the piecewise function is shown in Figure 10. As can be seen from the figure, the on-peak price (1.15 yuan/kwh) is nearly three times as high as that of off-peak (0.42 yuan/ kwh). Such a large price gap means that there is a huge potential to save energy consumption cost for the manufacturing industry.

Performance Measures.
ree performance measures [40] are used for performance comparison, including convergence measure IGD, dominance measure Ω, and diversity measure Δ. In addition, because the real optimal Pareto fronts of the tested problem are unknown, this paper approximates the union set of nondominant solutions of the four algorithms as the optimal Pareto solutions.
IGD evaluates the convergence of the algorithm and the distribution of nondominated solutions by calculating the minimum distance between each point on the optimal Pareto front and the Pareto front obtained by using a certain algorithm. It is a comprehensive performance measure. Because the units of two optimization objectives are different, it is necessary to normalize the objective values before participating in the calculation of IGD. e calculation formula is shown in (23), where |N * | is the number of nondominated solutions on the Pareto optimal front and N is the set of nondominated solutions obtained by using a certain algorithm. Obviously, the smaller the IGD, the better the convergence and distribution quality:  Complexity 9 e measure Ω denotes the percentage of nondominated solutions obtained by a certain algorithm in the optimal Pareto set, which is given by where |P(S)| is the number of nondominated solutions in the Pareto set S and H k is the set of nondominated solutions obtained by the k − th algorithm. e larger the value of Ω, the better the dominance performance of the algorithm. If the optimal Pareto set is obtained by using a certain algorithm H k , then Ω k = 100%. e measure Δ represents the diversity of the nondominated solutions, and it can be calculated by where N is the number of nondominated solutions in the obtained Pareto set. d f and d l are the Euclidean distances between two boundary solutions in the nondominated Pareto set and two extreme solutions on the optimal Pareto front. d i is the Euclidean distance between i and i + 1 consecutive solutions in the Pareto set. d is the mean Euclidean distance of all nondominated solutions. e smaller the value of Δ, the better the diversity performance.

Results and Discussion.
For each instance, each of the four algorithms runs 20 times independently, and each run gets a set of [IGD, Ω, Δ]. e minimum (Min), average (Avg), and maximum (Max) of each performance measure after 20 runs for the four algorithms are reported in the Table 4. e bold font represents the optimal results of each performance measure.
Minimum value, mean value, and maximum value can only show the solution effect of each algorithm from the macro level. Table 5 lists the statistical difference between the IMOALO and other algorithms based on the T-test. e bold font represents that it is significantly different at the significant level 0.05. It can be seen that IGD, Δ, and Ω measures of IMOALO have significant differences on almost all the problems with other three algorithms from Tables 4  and 5. erefore, at 95% confidence level, the proposed IMOALO algorithm has significant advantages over other algorithms. e box charts of three performance measures are shown in Figure 11. It can be seen from the figure that Ω measures of IMOALO are greater than in other algorithms; meanwhile, Δ and Ω measures of IMOALO are smaller than in other algorithms significantly, which further verifies the conclusion.
Take Sproblem-04-02 as an example to analyse the effectiveness of adding right-shift operation on the scheduling results. ere are 16 jobs, 8 machines, 1 reentrance, and 6 stages in the problem. Besides, the number of unrelated parallel machines at each stage is 1, 2, 1, 2, 1, and 1, respectively. Taking a nondominated solution 2-8-13-12-5-11-9-6-1-4-10-7-16-15-14-3 of IMOALO as an example, its objective function values C max � 151.5 hour and TEC min � 4991 yuan. In Figure 12, the Gantt charts with right-shift operation under TOU electricity tariffs are illustrated, and Figure 13 shows the Gantt charts (C max �151.5 hour and TEC min � 5133 yuan) without right-shift operation under TOU electricity tariffs. Comparisons show that the processing times of some operations have been shifted without affecting the production efficiency, and energy consumption cost is reduced by 2.76% to 142 yuan. In the processing cycle, the trend charts of total energy consumption of all machines are shown in Figure 14. It can be seen from Figure 14 that the electric load shifts obviously in some period. For example, in the vicinity where 45 hours and 75 hours are the periods with low electricity price, the total energy consumption of all machines increases after the right-shift procedure; on the contrary, in the vicinity where 55 hours and 115 hours are the periods with high electricity price, the total energy consumption of all machines decreases after the right-shift procedure. Obviously, the energy consumption cost can be saved by avoiding the periods with high electric price. To further prove the effectiveness of the right-shift operation, the above six smallscale problems run 20 times independently without considering right-shift operation under TOU electricity tariffs. Figure 15 shows the relative change of the average C max and TEC min of each problem (Sproblem is abbreviated as S), and the negative sign denotes relative decrease. From Figure 15, we can see that the energy consumption cost of each case has decreased. Although the makespan of most

Parameters
Range Remarks Number of jobs [10,20] Discrete integer values in [10,20] Number of stages [5,10] Discrete integer values in [5,10]   instances has increased slightly, it is acceptable that sacrificing a small amount of production efficiency for a large reduction in production cost. In summary, the RHFSP model considering energy consumption cost under the TOU electricity tariffs can reduce the energy consumption cost and enhance the competitiveness of enterprises; meanwhile, there is no significant impact on the makespan.  Figure 11: Box charts of three performance measures.

Analysis of Machines with Different Power and Speed.
Each stage contains several unrelated parallel machines with different power and speed in the RHFSP. For the purposes of cost saving, many enterprises mix new machines with old ones in practice. Although the new machine has a high processing capacity, it also consumes a large amount of energy per unit time. Taking the scheduling of 10 jobs, 2 reentrances, 2 stages, 2 parallel machines at each stage, and processing time ranges of [10,40] as an example, the effects of different parameter combinations on the scheduling results are studied. In this case, the speed and power of machines are divided into three levels. Assuming that the idle power of all machines is SW = 1. e specific parameters are shown in Table 6. To ensure the validity of the comparison, all machines at one stage have the same energy efficiency ratio (PW/V) and each stage has the same production capacity. For example, the PW/V at the stage 1 is 2 and that at the stage 2 is 3. In addition, the total speed at every stage is 6. e Pareto fronts obtained by running each parameter combination once randomly based on IMOALO algorithm are shown in Figure 16. As shown in Figure 16, there is Complexity 13 no significant difference on the makespan. Combining the high-power and low-power machines together is more energy efficient than using two identical machines with middle power. Although the faster machines consume more power, when assigning the jobs to the faster machines, the processing time of each operation is reduced. is helps to transfer the operations to the off-peak periods more flexibly.

Analysis of TOU Electricity Tariffs with Different Periods.
In practice, the period of TOU electricity tariffs varies with season. In general, the period of on-peak in summer is longer than that in winter, and there are differences between holidays and working days. us, it is important to study how the period of TOU electricity tariffs impacts the scheduling results. Taking Sproblem-04-02 as an example, we just change the total period without changing the corresponding electricity price for each period in the following five tests. Table 7 shows the specific parameters. e Pareto fronts obtained by each parameter combination based on IMOALO algorithm are shown in Figure 17.
As can be seen from Table 7 and Figure 17, maintaining the same electricity price in each period, energy    consumption cost is affected obviously by the length of total periods. However, the change has little effect on the makespan, and the lower bound of the makespan remains near 148 hours. e cure of Pareto front shows a downward trend with the increase of the periods. Tracing it to its cause, with increasing period of TOU electricity tariffs, there is more space for the operations transferring to off-peak periods after adding the right-shift operation. erefore, keeping the same electricity price in each period and the total number of TOU electricity tariffs periods, increasing the length of total periods appropriately can reduce energy consumption cost effectively and maintain the makespan unchanged.

Conclusions and Future Work
is paper mainly studies the RHFSP under TOU electricity tariffs with the objective of minimizing the makespan and energy consumption cost. e effectiveness of the proposed algorithm is evaluated by considering the benchmark problems. From the results, we can conclude that the proposed IMOALO algorithm can solve the RHFSP effectively and it is significantly superior to other algorithms in terms of convergence, dominance, and diversity measures of nondominated solutions. In addition, parameter analysis shows that the energy efficiency of machines and period of TOU electricity tariffs have great impact on the scheduling results. Under the TOU electricity tariffs, although the total electricity consumption has not decreased, the energy consumption cost can be saved by avoiding the periods of high electric price. Meanwhile, it helps to promote the use of green energy to generate electricity in low valleys, such as wind power, thus saving fossil energy and reducing pollutant emissions. In the future, we will further study the RHFSP, such as designing better coding and decoding methods and designing better algorithms and joint optimization of production scheduling and maintenance.

Data Availability
e data used to support the findings of this study were supplied by Professor Cho HM. ese datasets are cited at relevant places within the text as references [9]. Authors are not authorized to make them public.

Conflicts of Interest
e authors declare that they have no conflicts of interest.