Next Article in Journal
Towards the Application of Process Mining in the Mining Industry—An LHD Maintenance Process Optimization Case Study
Next Article in Special Issue
Leveraging Classical Statistical Methods for Sustainable Maintenance in Automotive Assembly Equipment
Previous Article in Journal
Developing a Controlling Model for Analyzing the Subjectivity of Enterprise Sustainability and Expert Group Judgments Using Fuzzy Triangular Membership Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hybrid Flow-Shop Scheduling Problems with Missing and Re-Entrant Operations Considering Process Scheduling and Production of Energy Consumption

1
Institute of Industrial Engineering, College of Mechanical Engineering, Zhejiang University of Technology, Hangzhou 310023, China
2
Jiaxing Sudoku Bridge Technology Co., Ltd., Jiaxing 314599, China
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(10), 7982; https://doi.org/10.3390/su15107982
Submission received: 2 April 2023 / Revised: 6 May 2023 / Accepted: 11 May 2023 / Published: 13 May 2023

Abstract

:
A hybrid flow shop scheduling model with missing and re-entrant operations was designed to minimize the maximum completion time and the reduction in energy consumption. The proposed dual-population genetic algorithm was enhanced with a range of improvements, which include the design of a three-layer gene coding method, hierarchical crossover and mutation techniques, and the development of an adaptive operator that considered gene similarity and chromosome fitness values. The optimal and worst individuals were exchanged between the two subpopulations to improve the exploration ability of the algorithm. An orthogonal experiment was performed to obtain the optimal horizontal parameter set of the algorithm. Furthermore, an experiment was conducted to compare the proposed algorithm with a basic genetic algorithm, particle swarm optimization algorithm, and ant colony optimization, which were all performed on the same scale. The experimental results show that the fitness value of the proposed algorithm is above 15% stronger than the other 4 algorithms on a small scale, and was more than 10% stronger than the other 4 algorithms on a medium and large scale. Under the condition close to the actual scale, the results of ten repeated calculations showed that the proposed algorithm had higher robustness.

1. Introduction

A hybrid flow-shop scheduling problem (HFSP) refers to a production shop with multiple processes and one or more parallel machines per process, arranged in an assembly line, and it is also referred to as a flexible flow-shop [1]. An HFSP is a combination of traditional flow shop scheduling and parallel-machine scheduling problems, and it has significant academic significance and application value. In other words, HFSP is a special kind of flow shop problem [2]. Furthermore, the HFSP can be extended to the re-entrant hybrid flow-shop problem (RHFS) [3] and hybrid flow shop scheduling problem with missing operations (HFSMO) [4], depending on the constraints of the specific process.
The “re-entrant shop” indicates that a workpiece may need to be processed several times on the same machine. In a study on semiconductor manufacturing processes, Kumar [3] proposed a third type of production system that differed from job shops and flow shops, wherein the workpiece needed to be repeatedly processed using certain machines or workstations at different stages of processing. This type of problem can be referred to as an RHFS problem. For the RHFS, most objective functions in the literature have focused on solving the makespan, with some studies examining objectives such as the weighted total completion time, cycle time, maximum delay time, total delay time, on-time completion rate, and early/delay penalties. For example, Chamnanlor et al. [5] aimed to maximize the system throughput by minimizing the maximum completion time based on the requirement of zero loss of work. Zhan et al. [6] discussed the production scheduling problem derived from a real rotor workshop and solved it with a minimum delay. Mousavis et al. [7] considered the setup time position-dependent learning effects and solved the maximum completion time and total delay as objective functions. Marichelvam et al. [8] used a hybrid monkey search algorithm to optimize the total flow time, considering the maximum completion time.
In recent years, with an increasing awareness of sustainability, more and more researchers have linked scheduling optimization to energy consumption and green costs. For example, Geng et al. [9] considered a right-shift operation and adjusted its start-up time to achieve the lowest possible energy cost by avoiding periods of high electricity prices. In addition to the objective function, a method for solving the RHFS is also the focus of this research. The HFSP has been proven to be NP-hard [10] and difficult to solve using traditional methods, whereas the RHFS, as an extension, is considered more difficult; therefore, researchers tend to adopt intelligent optimization algorithms to solve this problem. Geng et al. [11] developed an integrated scheduling model to minimize the maximum completion time, maximize the average agreement index, and design a hybrid NSGA-II based on the characteristics of the problem. Kun et al. [12] proposed a novel shuffle frog-jumping algorithm to handle the total delay and makespan. Wu et al. [13] discussed buffer sizes for batch machines and developed dosing methods, while proposing a greedy selection strategy to select machines for cold draw operations. Eskandadi et al. [14] generalized heuristics based on several basic scheduling rules and proposed a variable neighborhood search (VNS) to solve the problem of sequence-dependent setup times and unrelated parallel machines. Xu et al. [15] developed an improved moth–flame optimization algorithm to solve the green re-entrant hybrid flow shop scheduling problem. Qin et al. [16] proposed a rescheduling-based ant colony algorithm to solve the HFSP with uncertain processing time and introduced the concept of due date deviation to design a rolling-horizon-driven strategy that compressed the path of ant movement and reduced the cycle time in which to find a new solution. Nejad et al. [17] successfully traded off the relationship between process scheduling and production costs.
An HFSMO can be described as follows [18]: The sequence of processing machines in the production system is fixed, and the flow of workpieces to be processed is the same; however, some workpieces can skip process steps that do not need to be processed based on the characteristics of their own processes. This implies that not all workpieces have to undergo every process in a production system. Reading the literature [19] and analyzing it, we can find that “missing operations” are valuable for research. Although this type of problem is similar to traditional HFSP, it cannot simply be considered a unique case of a traditional HFSP due to the phenomenon of process omission. Setting the operation time of the skipped process stage to zero can result in the part being unable to skip the stage that does not require processing and being forced to wait for the process to complete. This not only leads to blockages but also delays in the start of the next process that the workpiece was originally intended for; in addition, this increases the idle time of the machine, which increases the completion time and affects the scheduling results. Therefore, research on HFSMO has important practical implications.
In 1985, Wittrock [4] first proposed an HFS containing missing operations in a practical production context based on an electronics processing plant and suggested that this scheduling model could be optimized using heuristic algorithms. Since then, several works of research have investigated this problem. Some researchers improved the existing intelligent algorithms. For example, Saravanan et al. [20] improved the simulated annealing (SA) algorithm and compared it with particle swarm optimization (PSO), demonstrating that the improved SA algorithm yielded better results and required a shorter computation time than with PSO. Li et al. [21] constructed a mathematical model of the HFSMO and then optimized it by improving the ABC algorithm to achieve certain results. Long et al. [22] proposed an improved GA which they applied to a steelmaking-continuous-casting production shop of a steel company. The final simulation experimental results indicated that the algorithm outperformed the other algorithms. Marichelvam et al. [23] proposed an improved hybrid genetic scatter search (IHGSS) algorithm by combining genetic operators and decentralized search algorithms, and they compared it with GA, a decentralized search algorithm, and the NEH heuristic algorithm to demonstrate the superior performance of the IHGSS algorithm. Saravanan et al. [24] fused GA with the SA algorithm and solved the problem model with an average delay as the scheduling objective using an improved GA.
In addition to the two aforementioned approaches, researchers have also designed new algorithms tailored to the characteristics of the problem. For instance, Lei et al. [25] devised a novel local search using a controlled deterioration algorithm to solve HFSMO. Dios et al. [26] constructed a problem model with completion time as the objective, performed a complexity analysis of HFSMO, and designed an effective heuristic algorithm. Siqueria et al. [27] proposed a new multi-objective VNS algorithm with which to solve problems with a maximum completion time and maximum total weighted delay as the optimization objectives.
After conducting an extensive review of the literature, it was found that research on RHFS is hot and has attracted many researchers. By contrast, the number of works in the literature on HFSMO is much less than the former. Both studies mainly focus on optimizing the maximum completion time, with a smaller number of researchers considering other objectives. Furthermore, the majority of existing studies do not consider important factors such as transport time, processing preparation time, and dynamic process jump constraints. Therefore, research on both needs to be further explored.
Our review of the literature revealed that the study of both the “re-entrant shop” and “workshops with missing operations” is a topic of great relevance and scientific value; however, to the best of the authors’ knowledge, no existing studies have considered them together. However, there are several hybrid flow shops in engineering that contain both re-entrant and missing operations. Therefore, this paper considers a combination of actual workshops, and then studies the new hybrid flow shop scheduling problem.
This study used the production organization of some products of an electric appliance manufacturing company as their research object and considered “re-entry” and “missing operations” together. This is referred to as the hybrid flow-shop scheduling problem with missing and re-entrant operations (HFS-MRO).
The main contributions to this paper are as follows:
  • The integration of RHFS and HFSMO
To the best of the authors’ knowledge, this is the first study to consider “re-entry” and “missing operations” in a unified manner, bringing scheduling research closer to real production.
  • The design of an improved dual-population genetic algorithm (IDPGA)
According to the processing characteristics, each chromosome of the proposed IDPGA adopts a three-layer gene coding mode and designs a “missing judgment vector” to integrate into the chromosome with which to determine the missing operations. To avoid the generation of non-feasible solutions, different crossover and mutation strategies are used for each coding layer, and adaptive operators based on gene similarity and chromosome fitness values are designed. The superiority and robustness of the algorithm are confirmed through the same-scale comparison experiment.
In this paper, “re-entry” and “missing operations” are considered together for the first time, and transportation and adjustment times are also taken into account, which could facilitate an in-depth study on hybrid flow shops. This could not only solve the actual factory problem but also provide a direction for follow-up researchers. Therefore, the research in this paper is of certain practical significance and scientific value.
The remainder of this paper is organized as follows: Section 2 describes the HFS-MRO and the mathematical model based on the problem characteristics. Section 3 presents the IDPGA, and Section 4 presents the computational experiments. Finally, Section 5 concludes the paper and suggests future research directions.

2. Problem Description and Mathematical Model of HFS-MRO

2.1. Problem Description

The HFS-MRO problem was created for the first time and could be described as having n workpieces that must be processed on s different workstations (processing stages) (Figure 1). Each workstation i comprises a set of m i ( m i 1 ) parallel machines with different machining capabilities. Each job has approximately the same processing flow, with minor differences in individual processes that can be skipped at one workstation or re-entered at another.

2.2. Model Assumptions

The following assumptions were considered based on the above problem description and the characteristics of an actual plant.
  • All workpieces and machines are ready at zero moments.
  • Operations can be ignored and re-entered.
  • At each processing stage, each process is exclusively performed on a single machine, and no machine is capable of processing more than one process simultaneously.
  • All processes for each workpiece are sequential and no misordering is allowed.
  • There is no priority requirement between the different workpiece processes. Once a workpiece has been processed, it cannot be interrupted, and there is no pre-emption.
  • The buffer capacity between two adjacent workstations is unlimited.
  • The workpiece is moved by workers between workstations, and the processing, preparation, and delivery times of each workpiece at each stage, as well as the processing energy consumption and idle energy consumption of each machine, are known and constant.
  • The case of machine breakdown is not considered.

2.3. Definition of Symbols

The parameter symbols and decision variables defined based on the model assumptions, are presented in Table 1.

2.4. Mathematical Planning Models

Objective function:
f 1 = min C max
f 2 = E C min = a = 1 M total r = 1 R a 1 ( D a r B a r ) P W a + a = 1 M total r = 1 R a 1 ( B a r + 1 D a r ) P I a
Equations (1) and (2) represent the two objective functions of the optimization problem: minimizing the maximum completion time and minimizing the machine energy consumption, respectively.
For the bi-objective problem, this paper converted a multi-objective problem into a single-objective problem by the weighted sum method as:
f = α 1 f 1 + α 2 f 2
where f 1 and f 2 represent the normalized values of the target functions f 1 and f 2 ,respectively, while α 1 and α 2 represent the weights of the target values, respectively. This then satisfies α 1 + α 2 = 1 .
S.t.
l = 1 m i r i j ( k 1 ) l ( S j ( k 1 ) + P j ( k 1 ) + T j ( k 1 ) k a u ) l = 1 m i r i j k l S j k + λ ( 1 Y j k a u ) , i , i , j ; k = { 2 , 3 , , N j } ; l = a , l = u
0 l = 1 m i r i j k l 1 , i , j , O j k U i
λ ( 2 r i j k l r i j k l ) + λ ( 1 Z j k j k ) + ( S j k S j k ) P j k + S T j j a , i , j < j , l = a , O j k U i , O j k U i
λ ( 2 r i j k l r i j k l ) + λ Z j k j k + ( S j k S j k ) P j k + S T j j a , i , j < j , l = a , O j k U i , O j k U i
λ ( 2 r i j k l r i j k l ) + ( S j k S j k ) P j k , i , j < j , l , O j k U i , O j k U i
D a r + S T j j a · x j k a r D a ( r + 1 ) P j k + λ ( 1 x j k a ( r + 1 ) ) j , j , a , k , k ; r = { 1 , 2 , , R a 1 }
Work pieces need to be manually moved between each workstation with transportation time taken into consideration. Therefore, constraint (4) ensures that the processing start time of operation O j k + 1 is not earlier than the sum of the processing completion and transportation times of operation O j k . When a certain operation in the workpiece can be missed, it means that the workpiece can skip the corresponding workstation. Therefore, the cumulative value of constraint (5) is zero. Conversely, it means that the workpiece must be processed on one machine in the workstation, and the cumulative value becomes 1. Constraints (6) and (7) ensure that each machine simultaneously processes one operation at most. Constraint (8) ensures that no wrong order occurs in the first processing and re-entry processing of the workpiece in the same process. There is a setting time before the machine starts processing; therefore, the machine’s starting processing time must be later than the sum of the completion time of the previous process of the current workpiece and the adjustment time of the machine itself. Hence, constraint (9) ensures that the start time of the current operation on a machine is not less than the sum of the completion time of the previous operation on that machine and the set-up time between the two operations.

3. Improved Dual-Population Genetic Algorithm

HFS-MRO, as an extension of HFS, also belongs to the NP-hard family. It was discovered that, as the scale increased, the solution time increased exponentially, and the precise solution method could not solve the problem. Therefore, we used a meta-heuristic algorithm to solve the problem. The genetic algorithm (GA) [28], as an evolutionary algorithm that simulates the law of superiority and inferiority in nature, was a good choice. It is also a neighborhood search algorithm with the core idea of evolving continuously in a solution space, selecting the offspring with high fitness by the selection operator, performing genetic operations on the offspring with the highest fitness, and stopping the algorithm by iterating a certain number of times or when the individuals reach the required fitness value.
GA is commonly used to solve shop scheduling problems because of its simple operation and strong adaptability. However, the limitations of traditional GA, in parent selection, crossover, and variation, cause the algorithm to easily fall into the local optima. Therefore, the IDPGA was designed in this study based on the idea of a dual-population GA to avoid this mentioned problem, as shown in Figure 2. IDPGA used a three-layer gene coding method to better integrate the “re-entrant” and “ignore process”. At the same time, the dual-population model was adopted alongside the hyperbolic tangent function to define the fitness similarity within the population, which successfully gave consideration to the “exploration and development” ability of the algorithm.

3.1. Encoding

There are two sub-problems to the general HFS problem: artifact ordering and machine assignment. The HFS-MRO problem requires “re-entry” and “missing” to be considered in a workpiece sequencing sub-problem. Therefore, this study used a three-level coding approach to represent individuals, including an operation vector, a machine vector, and a missing judgment vector. Each coding level was equal to the sum of the operations of all workpieces. Figure 3 shows a schematic of an individual chromosome with three workpieces, three workstations, and one round of re-entry. The number of machines in the workstation was [2 1 2].

3.1.1. Operation Vector

An integer operation-based coding approach was used, wherein each gene was represented by a workpiece number. The order in which the same workpiece number appeared indicated the different processes of the workpiece. In this study, “re-entry” was considered before “ignore”, as follows: first, a round of re-entry was considered with which to replicate the first three workstations, i.e., the number of operations in the workpiece increased from three to six, after which the ignore judgment vector in the second layer was used to determine which of the six operations should be ignored.

3.1.2. Missing Judgment Vector

Each part of the n vector was displayed as either “0” or “1” to determine if the corresponding operation was a “missing operation”. If “0” was displayed, then the operation corresponding to that part of the workpiece could not be ignored; however, if “1” was displayed, the operation corresponding to that part was an ignorable process. This layer of vector was considered an attribute layer of the first vector and followed the process vector as it changed.

3.1.3. Machine Vector

The machine vector consisted of n parts, each of which indicated the machine corresponding to the operation of a given workpiece. Furthermore, all machines for each workstation were numbered together starting from 1. The total number of machines was the sum of the number of machines in all processing stages.

3.2. Decoding

The decoding steps were as follows.
Step 1: Based on three-level coding, the genes in the operation sequencing section are read from left to right to determine the sequences of all operations;
Step 2: The status of each operation is determined based on the missing judgment vector;
Step 3: The machine corresponding to each operation is determined based on the machine vector;
Step 4: The actual processing and preparation times for each operation are obtained;
Step 5: The start and finish times for each operation are calculated.

3.3. Calculation of Fitness Values

To convert a multi-objective into a single objective, objective functions f 1 and f 2 can be normalized to obtain f 1 and f 2 using Equation (10). Then, using the weighted sum method, the multi-objective problem is converted to a single objective according to Equation (3) to obtain the objective function f .
f t = f t min f t max f t min f t ,   t = 1 , 2
where t = 1 , 2 denotes the objective function, and f t represents the normalized value of the objective f t . Max f t and min f t represent the upper and lower bounds of the objective f t , respectively. The smaller the objective function value of an individual in the population, the better the solution is. Therefore, the fitness function can be expressed as:
F i t = 1 f = 1 α 1 f 1 + α 2 f 2
Makespan and energy-consumption minimization were the two objectives of this study. The primary objective of line optimization was to maximize productivity and efficiency, whereas the secondary objective for green production and reducing energy costs was to minimize energy consumption. Therefore, in this study, the objective weights were set as α 1 = 0.6 and α 2 = 0.4 .

3.4. Initializing the Population

A simultaneous random first-tier permutation order and third-tier machine order were used to generate the initial population to increase population richness. Individual fitness values were calculated, with larger fitness values representing better solutions. Based on the fitness values, the initial population was divided into the top 50% of the outstanding sub-population I and the bottom 50% of the mediocre sub-population II.

3.5. Selection

Both populations were selected using the roulette wheel method.
Step1: Based on the obtained fitness values for each individual { F i t ( 1 ) , F i t ( 2 ) , F i t ( b ) , F i t ( W ) } , where W represents the size of each subpopulation, and b = 1 , 2 , W .
Step 2: The sum of all fitness values is used to obtain v = 1 W F i t ( b ) and to calculate the percentage of each individual fitness value in the sum as p ( b ) and the cumulative probability as q ( b ) .
p ( b ) = F i t ( b ) b = 1 W F i t ( b ) , b = 1 , 2 , W
q ( b ) = b = 1 W p ( b ) , b = 1 , 2 , W
Step 3: A random number is generated γ ( 0 , 1 ) . If γ q ( 1 ) , the first individual can be selected. If q ( b 1 ) γ q ( b ) , the b - t h individual can be selected.
Step 4: This step is then iterated W times.
An elite replacement strategy was used to accelerate the convergence of the outstanding subpopulation. In each iteration of the selection operation, the top 5% of individuals with the highest fitness values were selected to replace the bottom 5% of individuals with the lowest fitness values that completed the overall genetic operation for that round. This facilitated the rapid elimination of poorer individuals from the population.

3.6. Crossover

Considering the characteristics and various constraints of the HFS-MRO problem, a large number of non-feasible solutions were generated if the three vector layers were crossed directly. To avoid this scenario, this study designed different crossover methods for each layer of the vectors.
The first layer of the chromosome was crossed using a job-based crossover. The steps for this were as follows:
Step 1: All workpieces are randomly divided into two subsets j o b 1 and j o b 2 .
Step 2: The operation genes belonging to subset j o b 1 from parents P 1 and P 2 are copied to children O 1 and O 2 at the same location.
Step 3: The operation genes belonging to subset j o b 2 from parents P 1 and P 2 are added to children O 1 and O 2 , respectively, in the original order.
Figure 4 shows the operation vector crossover diagram, where j o b 1 = [ 2 ] and j o b 2 = [ 1 , 3 ] .
The second layer of the chromosome can be considered the attribute layer of the first layer, and, therefore, it follows the changes made by the first layer.
The third layer of the chromosome is crossed using a method similar to that described previously. The steps for this are as follows:
Step 1: Divide all workstations into two subsets s t a t i o n 1 and s t a t i o n 2 .
Step 2: Copy the genes from parents P 1 and P 2 belonging to subset s t a t i o n 1 onto children O 1 and O 2 , respectively.
Step 3: Copy the genes from parents P 1 and P 2 belonging to subset s t a t i o n 2 onto children O 1 and O 2 , respectively.
Figure 5 shows the strategy of the machine vectors for the crossing, where s t a t i o n 1 = [ 3 ] and s t a t i o n 2 = [ 1 , 2 ] .
Crossover probabilities are dynamically adaptive based on the genetic similarity of individuals within a subpopulation. To improve the local search capability of each subpopulation, the crossover probabilities of the outstanding subpopulation were positively correlated with the similarity of individuals within the species when performing crossovers. This increased the efficiency of the subpopulation when searching for optimal solutions. The crossover probability of the mediocre subpopulation was negatively correlated with the similarity of individuals within the species, which expanded the search range to identify new search spaces. A solution for gene similarity has already been reported [29]. In this study, the hyperbolic tangent function was considered to define the relationship between crossover probability and individual similarity, and the value range was taken to be [−1,1] in the domain of its function definition. ϕ ( φ 1 , φ 2 ) denotes the similarity between the individuals of two parents. The similarity between the parental chromosomes must be positive, and, therefore, the first quadrant of the hyperbolic tangent function can be selected. p c 1 denotes the crossover probability of the outstanding subpopulation, and p c 2 denotes the crossover probability of the mediocre subpopulation:
p c 1 = e ϕ e ϕ e ϕ + e ϕ , p c 2 = 1 e ϕ e ϕ e ϕ + e ϕ , ϕ [ 0,1 ]

3.7. Mutation

The mutation operation involves introducing a random change in some genes of the chromosome with a certain probability of generating new chromosomal individuals. This operation can help the algorithm improve its local search ability by increasing the population diversity with a small perturbation. The mutation of each layer of the chromosomal vector remains different, considering the specificity of HFS-MRO and the requirement of avoiding non-feasible solutions while also mimicking the general strategy of the crossover operator.
The first layer of the chromosome uses a two-point swap variant approach, wherein two genes g 1 and g 2 are selected randomly on the operation vector, and their positions are swapped with each other.
The third layer of the chromosomal machines, the machine vector, uses a single-point mutation; that is, a gene g 3 is selected randomly on the machine vector and is changed to another machine in the same workstation.
The variation probability affects the quality of the chromosome. A higher mutation probability can destroy good individual genes, and, therefore, a lower mutation probability is selected when the chromosome fitness value is large. Conversely, the mutation probability can be increased appropriately to widen the search space.
p m = p m m i n , F i t > F i t c v g p m m i n + p m m a x p m m i n F i t c v g F i t F i t c v g F i t m i n , F i t F i t avg  
where p m max and p m min denote the maximum and minimum mutation probabilities, respectively. Here, F i t a v g represents the average fitness value of the current subpopulation, and F i t represents the fitness value of the selected chromosome.

3.8. Elite Exchange

The basic dual-population GA undergoes genetic manipulation to add the best individuals from each population to the other population, facilitating the other population’s evolution and accelerating the convergence of the algorithm.
Based on this concept, the two subpopulations can be set to exchange the two best individuals in the outstanding subpopulation with two of the worst individuals in the mediocre subpopulation in terms of the fitness value after each cycle of the H generations. Simultaneously, the two best individuals in the mediocre subpopulation can be exchanged with the two worst individuals in the outstanding subpopulation to enhance the exploration capability of the algorithm.

4. Experiments and Result Analysis

4.1. Background of the Problem

Given the confidentiality constraints of the business, this section briefly describes the representative product process prototypes. The process flow diagrams for each of the five products are shown in Figure 6.
The process flow diagrams shown in Figure 6 indicate that each product had a consistent flow direction, with “re-entry” and “missing operation”. For example, product “KQ2373” was re-entered three times in operation “c” and never passed through operation “g”, i.e., this operation was ignored. Product “SY6832” was re-entered once in operation “c”, and operation “f” was ignored. A general HFS-MRO can be formulated based on this type of production line.

4.2. Experimental Data

This study used data from a company’s electrical production and processing plants as a case study to verify the feasibility of the algorithm. The data are not disclosed in this study owing to the confidentiality requirements of the actual case data provided by the company. In this study, the features of the data are extracted based on the actual case data of the enterprise; the rules for generating arithmetic cases are listed in Table 2. As shown in Table 2, all quantities are in “pcs”, the times in “min”, and the powers in “kw”. The process route includes a heat treatment process and occupies only one workstation.
This study assumed that the adjustment and transport times were approximately 20–60% and 30–50% of the standard processing times, considering the influence of worker handling. To facilitate the calculation, both processing and transport times were considered as integers with a minimum value of 1.

4.3. Parameter Setting

The parameter settings of the IDPGA algorithm affected the experimental results. The key factors that had to be set included the subpopulation size W , elite exchange H , and adaptive variation probability P m max and p m min . A four-factor, three-level orthogonal test was designed to determine these four parameters. Table 3 lists the orthogonal factor.
To ensure a realistic experiment, a medium-sized workpiece volume was selected for the experiment, i.e., the number of workpieces was 450. Other data were generated randomly from Table 2. Each combination was run ten times, and the average of the optimal objective function values was selected as the response value. When the number of iterations reached 500, the algorithm terminated. Table 4 lists the designed L9 orthogonal test and its corresponding response values. The main effect diagram of the mean value after the analysis in Minitab 20 is presented in Figure 7. The mean responses are listed in Table 5, where “Delta” represents the maximum average minus the minimum average of the orthogonal factors. Figure 7 shows that the IDPGA algorithm performs optimally when the subpopulation size was W = 150 , the elite exchange was H = 5 , p m max = 0.1 , and p m min = 0.01 .
Variance analysis was performed for p m max , which ranked fourth in Table 5. By keeping the other three parameter levels constant, each level of p m max was run 30 times, and the three columns of data were subjected to an ANOVA. The results are listed in Table 6. The test criterion for ANOVA was set at 0.05. The p-value in Table 6 is 0.02084, which is less than the set standard and indicates a significant difference among all levels. Therefore, W = 150 , H = 5 , p m max = 0.1 , and p m min = 0.01 were selected as the initial parameter values.

4.4. Algorithm Test

The effectiveness of the IDPGA was verified by comparing it to the basic GA, ant colony optimization (ACO), and PSO. The convergence of various algorithms was tested by selecting a realistic medium-scale algorithm example. These algorithms were encoded in a similar manner. The population size of the basic GA was set to twice that of W , the crossover probability was 0.9, and the mutation probability was 0.1. For PSO [30], the inertia weight was w = 0.8 and factor was c 1 = 2 and c 2 = 2 . For the basic ACO [31], the heuristic factor was α = 4 , the expectation factor was β = 6 , and the information volatilization coefficient was ρ = 0.5 . The parameters of the mayfly algorithm (MA) are derived from reference [32,33,34]. The algorithm termination condition was set to 500 iterations, MATLAB R2018b was adopted, and the test environment was the Windows 10 system, lntel® CoreTM i7-7700, 3.60 GHz processor with 16 GB of memory.
Figure 8 shows a comparison of the convergence of the four algorithms on a medium scale.
Figure 8 shows that the convergence effects of the five algorithms differed significantly. Although the initial solution of ACO was optimal, the convergence effect was insufficient and easily fell into the local optimal solution. Therefore, the results were not as positive as those of the other three algorithms. The initial solution of the PSO algorithm is insufficient, and the convergence rate was the slowest, even falling into local optima. The initial solution of the mayfly algorithm was similar to the ACO algorithm, but it was clearly seen that the convergence of the mayfly algorithm was stronger and faster, and the final result was only behind that of the IDPGA. The initial GA solution was acceptable; however, the convergence effect was unsatisfactory. In contrast, the proposed IDPGA used an adaptive genetic operator to accelerate the convergence speed and enhance the local search ability. Therefore, even if the initial solution was not satisfactory, it still achieved the highest solution accuracy, fastest solution speed, and optimum stability. Overall, this served to demonstrate the effectiveness of IDPGA.
To further illustrate the effectiveness of the algorithm, ten examples were randomly generated under various order sizes, according to Table 2, for a total of thirty examples. Four algorithms were solved for the thirty cases. To enhance accuracy, each case was run ten times, and the average of the optimal fitness values from the ten replicate trials was obtained. The optimization rates of the IDPGA relative to the various algorithms were also calculated, and the results are listed in Table 7 and Table 8. For example, the optimization rate of IDPGA relative to GA is calculated as follows:
o p t i m i z a t i o n   r a t e = F t I D P G A F i t G A F i t G A
Serial numbers 1–10, 11–20, and 21–30 represent small-scale, medium-scale, and large-scale cases, respectively. For a more intuitive representation, the fitness values of each algorithm listed in Table 7 are selected, and the results of the algorithm are compared as shown in Figure 9.
Table 7 and Table 8 and Figure 9 show that IDPGA outperformed the other four algorithms on all three scales. The average optimal fitness value of the IDPGA was higher than those of the other four algorithms, and a relative optimization rate was the highest when solving small-scale examples. Therefore, compared with other algorithms, IDPGA had significant advantages in solving this problem.
The company aims to improve productivity through algorithmic scheduling, and, therefore, in addition to the superior performance of the algorithms, the robustness of the algorithms was also investigated to solve the problem efficiently and consistently each time. A randomly generated medium-scale example was run ten times for each algorithm, using the objective function value as the evaluation metric. The results obtained were tallied and plotted on a box-line graph, as shown in Figure 10.
The boxplot in Figure 10 shows that IDPGA was more stable and robust than the other four algorithms in solving the HFS-MRO scheduling problem; it could meet the actual production needs of enterprises.

5. Conclusions

An HFS-MRO was extracted and analyzed based on a real shop, and an integer programming model was established. In addition, IDPGA was proposed and applied to solve this problem. An example of a generating table was designed by analyzing and extending the actual data of the enterprise. The optimal parameter set was obtained by an orthogonal experiment, and the IDPGA was compared with four other algorithms on the same scale to confirm the significant advantages of the IDPGA in solving this problem and satisfying the actual production requirements of enterprises. IDPGA was used to solve the practical problems of the factory, including the low efficiency caused by manual scheduling plans and a reduction in production energy consumption. This could serve to enhance the competitiveness of enterprises and conform to the sustainable development policy advocated by the government.
However, there are some limitations to this study. Due to manual handling between each workstation, the transportation time experienced a large fluctuation. In the future, staff training could be considered or AGVs introduced instead of manual handling, aimed at standardizing transportation time. This would not only make the research more consistent with the actual situation but also play a greater role in the enterprise.

Author Contributions

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

Funding

This research was funded by the National Key Research and Development Program of China under Grant 2018YFB1308102.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, Y.; Li, X.; Gao, L. Review on Hybrid Flow Shop Scheduling Problems. Chin. Mech. Eng. 2020, 31, 2798. [Google Scholar]
  2. Marichelvam, M.K.; Geetha, M. A memetic algorithm to solve uncertain energy-efficient flow shop scheduling problems. Int. J. Adv. Manuf. Technol. 2021, 115, 515–530. [Google Scholar] [CrossRef]
  3. Kumar, P.R. Re-entrant lines. Queueing Syst. Theory Appl. 1993, 13, 87–110. [Google Scholar] [CrossRef]
  4. Wittrock, R.J. Scheduling algorithms for flexible flow lines. IBM J. Res. Dev. 1985, 29, 401–412. [Google Scholar] [CrossRef]
  5. Chamnanlor, C.; Sethanan, K.; Chien, C.-F.; Gen, M. Re-entrant flow shop scheduling problem with time windows using hybrid genetic algorithm based on auto-tuning strategy. Int. J. Prod. Res. 2014, 52, 2612–2629. [Google Scholar] [CrossRef]
  6. Zhang, X.Y.; Chen, L. A re-entrant hybrid flow shop scheduling problem with machine eligibility constraints. Int. J. Prod. Res. 2018, 56, 5293–5305. [Google Scholar] [CrossRef]
  7. Mousavi, S.M.; Mahdavi, I.; Rezaeian, J.; Zandieh, M. Bi-objective scheduling for the re-entrant hybrid flow shop with learning effect and setup times. Sci. Iran. 2018, 25, 2233–2253. [Google Scholar] [CrossRef]
  8. Marichelvam, M.K.; Tosun, O.; Geetha, M. Hybrid monkey search algorithm for flow shop scheduling problem under makespan and total flow time. Appl. Soft Comput. 2017, 55, 82–92. [Google Scholar] [CrossRef]
  9. Geng, K.; Ye, C.; Dai, Z.H.; Liu, L. Bi-Objective Re-Entrant Hybrid Flow Shop Scheduling considering Energy Consumption Cost under Time-of-Use Electricity Tariffs. Complexity 2020, 2020, 8565921. [Google Scholar] [CrossRef]
  10. Gupta, J.N.D. Two-stage, hybrid flowshop scheduling problem. J. Oper. Res. Soc. 1988, 38, 359–364. [Google Scholar] [CrossRef]
  11. Geng, K.; Wu, S.; Liu, L. Multi-objective re-entrant hybrid flow shop scheduling problem considering fuzzy processing time and delivery time. J. Intell. Fuzzy Syst. 2022, 43, 7877–7890. [Google Scholar] [CrossRef]
  12. Kun, W.; Lamei, X.; Cong, L. A novel SFLA for re-entrant hybrid flow shop scheduling. J. Phys. Conf. Ser. 2019, 1314, 012126. [Google Scholar]
  13. Wu, X.; Cao, Z. An improved multi-objective evolutionary algorithm based on decomposition for solving re-entrant hybrid flow shop scheduling problem with batch processing machines. Comput. Ind. Eng. 2022, 169, 108236. [Google Scholar] [CrossRef]
  14. Eskandari, H.; Hosseinzadeh, A. A variable neighbourhood search for hybrid flow-shop scheduling problem with rework and set-up times. J. Oper. Res. Soc. 2014, 65, 1221–1231. [Google Scholar] [CrossRef]
  15. Xu, F.; Tang, H.T.; Xun, Q.N.; Lan, H.Y.; Liu, X.; Xing, W.F.; Zhu, T.Y.; Wang, L.; Pang, S.B. Research on Green Reentrant Hybrid Flow Shop Scheduling Problem Based on Improved Moth-Flame Optimization Algorithm. Processes 2022, 10, 2475. [Google Scholar] [CrossRef]
  16. Qin, W.; Zhang, J.; Song, D. An improved ant colony algorithm for dynamic hybrid flow shop scheduling with uncertain processing time. J. Intell. Manuf. 2018, 29, 891–904. [Google Scholar] [CrossRef]
  17. Nejad, M.G.; Shavarani, S.M.; Vizvari, B.; Barenji, R.V. Trade-off between process scheduling and production cost in cyclic flexible robotic cells. Int. J. Adv. Manuf. Technol. 2018, 96, 1081–1091. [Google Scholar] [CrossRef]
  18. Geng, K.-F.; Ye, C.-M. Joint scheduling of machines and AGVs in green hybrid flow shop with missing operations AGV. Kongzhi Juece Control Decis. 2022, 37, 2723–2732. [Google Scholar]
  19. Tseng, C.-T.; Liao, C.-J.; Liao, T.-X. A note on two-stage hybrid flowshop scheduling with missing operations. Comput. Ind. Eng. 2008, 54, 695–704. [Google Scholar] [CrossRef]
  20. Saravanan, M.; Sridhar, S.; Harikannan, N. Optimization of Two-Stage Hybrid Flow Shop Scheduling Problems Using Genetic Algorithm. Appl. Mech. Mater. 2015, 766–767, 962–967. [Google Scholar]
  21. Li, J.Q.; Pan, Q.K.; Duan, P.Y. An Improved Artificial Bee Colony Algorithm for Solving Hybrid Flexible Flowshop With Dynamic Operation Skipping. IEEE Trans. Cybern. 2016, 46, 1311–1324. [Google Scholar] [CrossRef]
  22. Long, J.; Zheng, Z.; Gao, X.; Pardalos, P.M. Scheduling a realistic hybrid flow shop with stage skipping and adjustable processing time in steel plants. Appl. Soft Comput. 2018, 64, 536–549. [Google Scholar] [CrossRef]
  23. Marichelvam, M.K.; Prabaharan, T. Performance evaluation of an improved hybrid genetic scatter search (IHGSS) algorithm for multistage hybrid flow shop scheduling problems with missing operations. Int. J. Ind. Syst. Eng. 2014, 16, 120–141. [Google Scholar] [CrossRef]
  24. Saravanan, M.; Sridhar, S.; Harikannan, N. Minimization of mean tardiness in hybrid flow shop with missing operations using genetic algorithm. J. Adv. Manuf. Syst. 2016, 15, 43–55. [Google Scholar] [CrossRef]
  25. Lei, D.; Guo, X. Hybrid flow shop scheduling with not-all-machines options via local search with controlled deterioration. Comput. Oper. Res. 2016, 65, 76–82. [Google Scholar] [CrossRef]
  26. Dios, M.; Fernandez-Viagas, V.; Framinan, J.M. Efficient heuristics for the hybrid flow shop scheduling problem with missing operations. Comput. Ind. Eng. 2018, 115, 88–99. [Google Scholar] [CrossRef]
  27. Siqueira, E.D.; Souza, M.; Souza, S.D. A Multi-objective Variable Neighborhood Search algorithm for solving the Hybrid Flow Shop Problem. Electron. Notes Discret. Math. 2018, 66, 87–94. [Google Scholar] [CrossRef]
  28. Katoch, S.; Chauhan, S.S.; Kumar, V. A review on genetic algorithm: Past, present, and future. Multimedia Tools Appl. 2021, 80, 8091–8126. [Google Scholar] [CrossRef] [PubMed]
  29. Chen, C.; Zhang, D.; Ren, P. Synthetic similarity measure for case retrieval in case-based reasoning diagnosis system. Chin. J. Mech. Eng. 2004, 40, 48–52. [Google Scholar] [CrossRef]
  30. Shi, Y.; Eberhart, R. A modified particle swarm optimizer. In Proceedings of the 1998 IEEE International Conference on Evolutionary Computation Proceedings, Anchorage, AK, USA, 4–9 May 1998. [Google Scholar]
  31. Alaykýran, K.; Engin, O.; Döyen, A. Using ant colony optimization to solve hybrid flow shop scheduling problems. Int. J. Adv. Manuf. Technol. 2007, 35, 541–550. [Google Scholar] [CrossRef]
  32. Zervoudakis, K.; Tsafarakis, S. A mayfly optimization algorithm. Comput. Ind. Eng. 2020, 145, 106559. [Google Scholar] [CrossRef]
  33. Zhao, C.; Lui, C.F.; Du, S.C.; Wang, D.; Shao, Y.P. An earth mover’s distance based multivariate generalized likelihood ratio control chart for effective monitoring of 3D point cloud surface. Comput. Ind. Eng. 2023, 175, 108911. [Google Scholar] [CrossRef]
  34. Shao, Y.P.; Xu, F.C.; Chen, J.; Lu, J.S.; Du, S.C. Engineering surface topography analysis using an extended discrete modal decomposition. J. Manuf. Process. 2023, 90, 367–390. [Google Scholar] [CrossRef]
Figure 1. Schematic of HFS-MRO.
Figure 1. Schematic of HFS-MRO.
Sustainability 15 07982 g001
Figure 2. Flowchart of the IDPGA.
Figure 2. Flowchart of the IDPGA.
Sustainability 15 07982 g002
Figure 3. Chromosome encoding.
Figure 3. Chromosome encoding.
Sustainability 15 07982 g003
Figure 4. Schematic of the operation vector crossover.
Figure 4. Schematic of the operation vector crossover.
Sustainability 15 07982 g004
Figure 5. Schematic of the machine vector.
Figure 5. Schematic of the machine vector.
Sustainability 15 07982 g005
Figure 6. Process flow diagram.
Figure 6. Process flow diagram.
Sustainability 15 07982 g006
Figure 7. Mean main effects plot.
Figure 7. Mean main effects plot.
Sustainability 15 07982 g007
Figure 8. Comparison of convergence effect.
Figure 8. Comparison of convergence effect.
Sustainability 15 07982 g008
Figure 9. Comparison of algorithmic results: (a) Small scale; (b) Medium scale; (c) Large scale.
Figure 9. Comparison of algorithmic results: (a) Small scale; (b) Medium scale; (c) Large scale.
Sustainability 15 07982 g009
Figure 10. Robustness comparison.
Figure 10. Robustness comparison.
Sustainability 15 07982 g010
Table 1. Definition of symbols.
Table 1. Definition of symbols.
SymbolDefinition
n Number of jobs
j Index for jobs, j = 1 , 2 , , n
S Number of workstations (stages)
i Index of workstations, i = 1 , 2 , , s
m i Number of identical parallel machines in workstation i
l Index for machines in workstation i ; l = 1 , 2 , , m i
N j Number of operations for job j
k Index of operations for job j ; k = 1 , 2 , , N j
O j k The k - t h operation on job j
P j k Processing time of O j k
S j k Starting time of O j k
U i Collection of operations processed on the workstation i
C j Completion time of job j
P W a Machine a ’s energy consumption per unit of time during processing operations
P I a Energy consumption per unit of time when machine a is idle
M t o t a l Total number of machines at each workstation, M t o t a l = m 1 + m 2 + + m s
a Index number of the machine, a = 1 , 2 , , M t o t a l
R a Total number of operations on the machine a
B a r Start time of the r - t h processing on machine a
D a r Completion time of the r - t h processing on machine a
T j ( k 1 ) k a u Transport time of the two consecutive operations O j k 1 and O j k of the job j from the machine a to the machine u
S T j j a Adjustment time between two adjacent workpieces j and j on machine a
λ a sufficiently large positive number
Decision variables:
r i j k l 1, if O j k processed on machine l in workstation i ; 0, other
Z j k j k 1, if O j k prioritizes O j k ; 0, other
Y j k a u 1, if O j k is processed on machine u while O j ( k 1 ) is processed on machine a ; 0, other
x j k a r 1, if O j k corresponds to the r - t h process of machine a ; 0, other
Table 2. Table of arithmetic cases generation rules.
Table 2. Table of arithmetic cases generation rules.
ParametersRange of Values
Number of workpieces{30,450,1200}
Number of workstationsU [5,10]
Number of re-entriesU [1,3]
Number of missing operationsU [1,4]
Number of machines per workstationU [2,5]
Processing time of heat treatmentU [100,200]
Processing time of other stagesU [1,30]
Processing power of heat treatment machinesU [300,450]
Idle power of heat treatment machinesU [100,200]
Processing power of other machinesU [5,10]
Idle power of other machinesU [0.5,3]
Table 3. Orthogonal factor levels.
Table 3. Orthogonal factor levels.
Factor W H P m max p m min
Levels
15050.0010.10
2100100.0050.15
3150150.0100.20
Table 4. Orthogonal test table.
Table 4. Orthogonal test table.
Serial NumberFactorResponse Value
W H P m max p m min
15050.100.0010.123468
250100.150.0050.125570
350150.200.0100.124798
410050.150.0100.120744
5100100.200.0010.126256
6100150.100.0050.123383
715050.200.0050.121117
8150100.100.0100.122076
9150150.150.0010.123383
Table 5. Mean response table.
Table 5. Mean response table.
Levels W H p m max p m min
10.12460.12180.12300.1244
20.12350.12460.12320.1234
30.12220.12390.12410.1225
Delta0.00240.00290.00110.0018
Ranking2143
Table 6. Analysis of variance.
Table 6. Analysis of variance.
Source of DifferenceSSdfMSFp-ValueF Crit
Intergroup difference0.0001426.81 × 10−54.483050.020843.35413
Intra-group variation0.00041271.52 × 10−5
Total0.0005529
Table 7. Comparison of algorithmic results.
Table 7. Comparison of algorithmic results.
Serial NumberIDPGAGAOptimization RatePSOOptimization RateACOOptimization RateMAOptimization Rate
111.019.98510.27%9.84511.83%9.7313.16%10.217.84%
26.2515.05423.68%5.05523.66%4.96625.88%5.1920.44%
37.6156.74612.88%6.49117.32%6.9499.58%6.6314.86%
48.4787.56812.02%7.05220.22%7.31215.95%7.4713.49%
59.7848.9479.36%8.08221.06%8.33617.37%8.0920.94%
65.5214.66218.43%4.60519.89%4.78315.43%4.8314.31%
76.8935.59923.11%5.36428.50%5.65921.81%5.13734.18%
87.4586.39716.59%6.45715.50%6.68711.53%6.3816.90%
99.1428.05213.54%7.66819.22%8.6136.14%8.21211.32%
107.9656.99513.87%6.82816.65%6.68619.13%6.718.88%
1111.78210.6610.53%10.6610.53%10.828.89%10.3314.06%
1213.56312.0812.28%12.429.20%12.3310.00%12.0112.93%
139.6128.71910.24%8.7719.59%8.72110.22%8.7310.10%
1411.2129.79214.50%10.1710.25%10.219.81%10.229.71%
159.4538.52110.94%8.39212.64%8.41512.34%8.4412.00%
168.8977.92712.24%7.87612.96%8.2198.25%7.8812.91%
178.5117.62411.63%7.8278.74%7.5612.58%7.5512.73%
188.1277.36210.39%7.17913.21%7.32710.92%7.212.88%
1911.56511.133.91%11.133.91%11.044.76%115.14%
2010.3789.30611.52%9.27111.94%9.6277.80%9.3311.23%
218.9718.3567.36%8.3287.72%8.01211.97%8.288.35%
227.8867.2169.28%7.113510.86%7.1789.86%7.219.38%
2313.20111.7112.73%12.416.37%12.178.47%12.19.10%
247.896.75116.87%6.68817.97%6.95913.38%6.8115.86%
2512.13410.416.67%10.5115.45%11.277.67%10.5914.58%
2611.15610.0111.45%10.278.63%9.78913.96%10.189.59%
2710.5719.10716.08%9.1415.66%9.46611.67%9.4112.34%
289.458.4112.37%8.9295.83%8.8077.30%8.856.78%
297.2416.33614.28%6.52610.96%6.2715.49%6.4212.79%
306.2515.05423.68%5.05523.66%4.96625.88%5.1920.44%
Table 8. Average relative optimization rate of IDPGA.
Table 8. Average relative optimization rate of IDPGA.
AlgorithmOptimization Rate Relative to GAOptimization Rate Relative to PSOOptimization Rate Relative to ACOOptimization Rate Relative to MA
Scale
Small15.37%19.39%15.60%18.36%
Medium10.82%10.30%9.56%11.37%
Large12.90%11.07%11.43%11.28%
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

Tang, H.; Zhou, J.; Shao, Y.; Yang, Z. Hybrid Flow-Shop Scheduling Problems with Missing and Re-Entrant Operations Considering Process Scheduling and Production of Energy Consumption. Sustainability 2023, 15, 7982. https://doi.org/10.3390/su15107982

AMA Style

Tang H, Zhou J, Shao Y, Yang Z. Hybrid Flow-Shop Scheduling Problems with Missing and Re-Entrant Operations Considering Process Scheduling and Production of Energy Consumption. Sustainability. 2023; 15(10):7982. https://doi.org/10.3390/su15107982

Chicago/Turabian Style

Tang, Hongtao, Jiahao Zhou, Yiping Shao, and Zhixiong Yang. 2023. "Hybrid Flow-Shop Scheduling Problems with Missing and Re-Entrant Operations Considering Process Scheduling and Production of Energy Consumption" Sustainability 15, no. 10: 7982. https://doi.org/10.3390/su15107982

APA Style

Tang, H., Zhou, J., Shao, Y., & Yang, Z. (2023). Hybrid Flow-Shop Scheduling Problems with Missing and Re-Entrant Operations Considering Process Scheduling and Production of Energy Consumption. Sustainability, 15(10), 7982. https://doi.org/10.3390/su15107982

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