Abstract

This paper investigates a two-machine flow shop problem with release dates in which the job processing times are variable according to a learning effect. The bicriterion is to minimize the weighted sum of makespan and total completion time subject to release dates. We develop a branch-and-bound (B&B) algorithm to solve the problem by using a dominance property, several lower bounds, and an upper bound to speed up the elimination process of the search tree. We further propose a multiobjective memetic algorithm (MOMA), enhanced by an initialization strategy and a global search strategy, to obtain the Pareto front of the problem. Computational experiments are also carried out to examine the effectiveness and the efficiency of the B&B algorithm and the MOMA algorithm.

1. Introduction

In many industrial systems, logistics, and service settings, the processing times of jobs can decrease due to the learning effects, i.e., firms and employees perform a job (task) over and over and they learn how to perform more efficiently [1, 2]. This is especially observed in seru production system, consisting of some equipment and workers who produce one or more part types. After undertaking cross-training, the learning effects greatly affect their efficiencies [3]. “As there is a significant involvement of humans in scheduling environments, the number of activities subject to learning is high, too. Hence it seems reasonable to consider learning in scheduling environments” [2]. On the other hand, the importance of flow shop scheduling optimization is widely recognized in many manufacturing and assembling processes, e.g., the internet connectivity problem in 3-tiered client-server databases is a two-machine flow shop problem (see [4]); the scheduling of multimedia date objectives for WWW applications reduces to a two-machine flow shop problem (see [57]; see also the reviews given in [6] and [810]).

More recently, Xu et al. [11], Eren [12], Wu and Lee [13], Rudek [14], Kuo et al. [15], Lee and Chung [16], Sun et al. [17, 18], Cheng et al. [19], Li et al. [20], Wang and Zhang [21], J. B. Wang and J. J. Wang [22], Lai et al. [23], Liu and Feng [24], Wang and Zhang [21], Wu et al. [25], Shiau et al. [26], Lu [27], Qin et al. [28], He [29], and Wang et al. [30] considered flow shop scheduling problems with learning effects, but without release dates. Xu et al. [11] considered the flow shop scheduling problems with position-dependent learning effect, i.e., if job is in position of a schedule, the actual processing time of job on machine is and , where denotes the normal processing time of job on and , , and are the learning rates. For several regular objective functions, they presented approximate algorithms. Cheng et al. [19] considered the model . For the maximum lateness minimization, they gave a mathematical programming model. Wang et al. [31] considered the model . For the total completion time minimization, they gave a branch-and-bound algorithm and several well-known heuristics. Lai et al. [23] considered the model . For the total tardiness minimization, they gave a branch-and-bound algorithm and two heuristics. He [29] and Li et al. [20] considered the flow shop scheduling problems with time-dependent learning effect, i.e., if job is in position of a schedule, the actual processing time of job on machine is , where is the learning rate and denotes the job that occupies the th position in a sequence. For several regular objective functions, they presented approximate algorithms. Shiau et al. [26] considered the flow shop scheduling problems with general position-dependent learning effect, i.e., if job is in position of a schedule, the actual processing time of job on machine is , where is a nonincreasing function. Sun et al. [18] considered the total weighted completion time minimization flow shop scheduling problem. For three position-dependent learning effects, they gave heuristic algorithms. Cheng et al. [19], Lai et al. [23], Wang et al. [31], Wu and Lee [13], and Wang et al. [30] considered the flow shop scheduling with truncated learning effects. J.-B. Wang and J.-J. Wang [22] considered the flow shop scheduling with a general exponential learning effect. For five regular objective functions, they presented heuristic algorithms. Wang and Zhang [21] considered the flow shop scheduling with position-weighted learning effects. For the weighted sum of makespan and total completion time minimization, they gave a branch-and-bound algorithm and some heuristic algorithms. Qin et al. [28] considered the flow shop group scheduling with position-based learning effect. For four objectives (i.e., the makespan, total completion time, total weighted completion time, and maximum lateness), they gave several heuristics and metaheuristics. For new trends in flow shop scheduling with learning effects, we refer the reader to Rudek [14], Liu and Feng [24], Shiau et al. [26], Lu [27], and He [29].

However, in fact, the scheduling with release dates is interesting and closer to real problems [3234]. With modeling a realistic production system in mind, Yin et al. [35] considered single-machine scheduling with release dates and position-dependent learning effects. To the best of our knowledge, the flow shop scheduling problems with a learning effect and release dates (i.e., the ready times) simultaneously have barely been investigated. Bai et al. [36] is the only identifiable exception. Hence, in this paper, we consider the two-machine flow shop scheduling problem with position-dependent learning effect and release dates, and the objective is to minimize the weighted sum of makespan and total completion time. Obviously, the problem under study is NP-hard [34]; thus, the branch-and-bound algorithm might be a good way to obtain the optimal schedule. We further design an evolutionary multiobjective optimization algorithm to obtain the Pareto front of high quality. This contribution makes a clear distinction between this paper and Bai et al. [36], which solely focuses on a single aggregated objective.

The remaining of the paper is organized as follows. In Section 2, we formalize the problem. Dominance conditions for the problem are presented in Section 3. Branch-and-bound algorithm is discussed in Section 4. Numerical experiments are conducted in Section 5. The last section concludes the paper.

2. Notations and Assumptions

The problem may be stated as follows. We are given jobs and two machines and . Each machine is available at time zero and the job becomes available at its release date . Let represent the set of jobs which are to be processed on 2 machine permutation flow shop settings. Each job is required to be processed on machine and then . As in Biskup [2], Lee and Chung [16], Lee and Wu [37], and Lee et al. [38], if job is in position of a schedule, then the actual processing times and of job on machines and are respectively, where denotes the normal processing time (i.e., the basic processing time before learning effects happened) of job on and is the learning index for both machines and depending on the learning rate LR. For a given schedule , let be the completion time of job on machine . The objective is to minimize , where is the makespan and is the total completion time, . Using the conventional three-field notation for describing scheduling problems (Graham et al. [39]), the problem under consideration can be described as , where denotes the learning effect [1, 2].

3. Dominance Conditions

In this section, we give some dominance rules for . “Dominance rules are basically necessary conditions for any optimal schedule that can be generated. Applying such rules results in a set of precedence relations between jobs. These precedence relations are then used to reduce the number of branches in a branch-and-bound search tree” (Yin et al. [35]). Let and be obtained from by only interchanging jobs and , where and are partial sequences, and there are jobs in . To show that dominates , it is sufficient to show that , , and . Let denote the last completion time prior to jobs and on machine in ; obviously, and remain unchanged in .

Proposition 1. If the jobs and satisfy the following conditions: (i)Either or or (ii)Either or or (iii)Either or or (iv)Either or (v)Either or or (vi)Either or (vii)Either or or or or or (viii)Either or or or or or (ix)Either or or or or or (x)Either or or or or or (xi)Either or or or or or (xii)Either or or or or or then dominates .

Proof 1. Under schedules and , the completion times of jobs and are If the cases , , , , , , , , , , , and can be satisfied, we have the following: first term in (7) ≤ first term in (3), second term in (7) ≤ first term in (3), third term in (7) ≤ first term in (3), first term in (8) ≤ first term in (4), second term in (8) ≤ first term in (4), third term in (8) ≤ first term in (4), first term in (9) ≤ first term in (5), second term in (9) ≤ first term in (5), third term in (9) ≤ first term in (5), fourth term in (9) ≤ first term in (5), fifth term in (9) ≤ first term in (5), and sixth term in (9) ≤ first term in (5).
Hence, , , , and . Similarly to this case, the other cases can be obtained.

4. Branch-and-Bound (B&B) Algorithm

4.1. Lower Bound

Let be a schedule in which is the scheduled part, and suppose there are jobs in , and is the set of so far unscheduled jobs. Let denote the th smallest normal processing time on machine ( and 2) when jobs in are arranged in an ascending order of the normal processing time on machine . By definition, the completion time of the th job on machine is

Similarly, where .

Hence, where and .

Hence,

It is noticed that the term is a fixed constant and the term is minimized by sequencing the jobs in a nondecreasing order of . Consequently, we obtain the first lower bound which is where , , and .

Similarly, where .

Hence, and the total completion time is

Hence,

Obviously, is a fixed constant, is minimized by choosing and , and can be minimized by sequencing the jobs in a nondecreasing order of and , respectively. Consequently, we have where and is a nondecreasing order of the normal processing times on for the remaining unscheduled jobs (note that and do not necessarily correspond to the same job).

Similarly, where .

Hence, and

Hence,

It is noticed that the term is a fixed constant and the term can be minimized by sequencing the jobs in a nondecreasing order of . Consequently, we have where is a nondecreasing order of the normal processing times on for the remaining unscheduled jobs.

In order to make the lower bound tighter, we choose the maximum values , , as a lower bound for , i.e.,

4.2. The Algorithm for the Upper Bound

Similarly to Framinan and Leisten [40] (they proposed the FL heuristic for solving ), the FL heuristic can be adjusted for solving . Now, we give the modified FL algorithm.

Phase I
 Step 1. Sequence the jobs in nondecreasing order of
 Step 2. Sequence the jobs in nondecreasing order of
 Step 3. Sequence the jobs in nondecreasing order of .
 Step 4. Sequence the jobs in nondecreasing order of
 Step 5. Choose the best solution from Steps 1 to 4
Phase II
 Step 1. Let be the schedule obtained from Phase I
 Step 2. Set . Select the first two jobs from and select the better between the two possible sequences
 Step 3. Increment , . Select the th job from and insert it into possible positions of the best partial sequence obtained so far. Among the sequences, the best -job partial sequence is selected based on minimum . Next, determine all possible sequences by interchanging jobs in positions and of the above partial sequence for all , . Select the best partial sequence among sequences having minimum
 Step 4. If , then STOP; else, go to Step 3

Since Phase I takes time and Phase II takes time, the overall time complexity of the FLTPA algorithm is .

4.3. The B&B Algorithm

A B&B algorithm that uses the depth first search strategy to solve is proposed. Now the B&B algorithm can be described as follows.

Step 1. Initialization: the FL-based two-phase algorithm (FLTPA) is applied to obtain an upper bound.

Step 2. Fathoming: apply Proposition 1 to eliminate the dominated partial sequences from the initial node and their descendants from the tree.

Step 3. Bounding: calculate the lower bound for the node. If the lower bound for an unfathomed partial schedule is larger than the upper bound, eliminate the node and all the nodes following it in the branch. Calculate the objective function value of the completed schedule; if it is less than the upper bound, replace it as the new solution; otherwise, eliminate it.

Step 4. Termination: continue to search all the nodes, and the remaining upper bound is optimal.

4.4. The Improved EMO Algorithm

The B&B algorithm outputs the optimal solution, and the LFTPA heuristic outputs a near-optimal solution to the linearly weighted sum of and . From the perspective of biobjective optimization, evolutionary multiobjective optimization (EMO) algorithm returns at one time a set of nondominated solutions [4143], i.e., the Pareto front, for the decision-maker’s reference. We propose two strategies to improve the multiobjective memetic algorithm (MOMA) which has successfully solved many difficult numerical optimization problems and outperforms NSGA-II (the nondominated sorting genetic algorithm) and SPEA2 (the improved strength Pareto evolutionary algorithm) for the two-objective and three-objective benchmark flow shop scheduling instances [44]. (1)Initialization strategy: the initial population consists of one individual generated by the LFTPA heuristic and some others adjusted from that individual by exchanging operators. The number of the LFTPA heuristic-related individuals is equal to half of the population size. The remaining individuals are randomly generated to guarantee population diversity.(2)Global search strategy: due to the nondominated sorting and crowd distance calculation mechanisms, NSGA-II has a very good performance to solve many multiobjective problems and is used as global search. Its key operators include selection, crossover, and mutation. We apply the tournament selection operator (see the details in [44]). According to Ishibuchi et al. [45], the two-point crossover and insertion mutation operators are adopted.

5. Computational Study

Computational experiments were studied to evaluate the effectiveness of the B&B algorithm and the FLTPA heuristic. We coded the B&B algorithm, the FLTPA heuristic, and the improved MOMA in VC++ 6.0 and ran on CPU Intel® Core (TM) i7-4790 3.6 GHz and 8 GB RAM. For the experiments, the parameters are considered as follows: and the numbers of jobs : 11, 12, 13, 14, and 15 for small size and 500 and 1000 for large size. The normal processing times of jobs ( and ) are randomly generated from the uniform distribution over the integers [1, 100].

In order to study the impact of the parameters, the values of the learning effect are taken to be 90% (i.e., ), 80% (i.e., ), and 70% (i.e., ). The release dates of jobs are randomly generated from the uniform distribution over the integers [1, 10,1, 50] and [1, 100].

As a consequence, for small-size problems, 45 experimental conditions were examined and 50 replications were randomly generated for each condition. A total of 2250 problems were tested. For large-size problems, 18 experimental conditions were examined and 50 replications were randomly generated for each condition. A total of 900 problems were tested. For the B&B algorithm, the average and maximum number of nodes and the average and the maximum time (in milliseconds) are reported. The percentage error of the solution produced by the FLTPA heuristic is calculated as

The results of small-size problems are summarized in Tables 13. From Tables 13, we find that the B&B algorithm can solve a problem of up to 15 jobs in a very short amount of time. The most time-consuming instance took a maximum of 118311 ms for , , and . As for the performance of the FLTPA heuristic, it is seen that the FLTPA heuristic performs very well for the error percentages. The mean error percentage is less than 0.7% for all the tested cases, and the max error percentage is less than 6% for all the tested cases.

For large-size problems, the termination of the proposed MOMA algorithm is controlled by setting a running time limit of 100 seconds. In order to verify the effectiveness and efficiency, the MOMA is compared with several state-of-the-art EMO algorithms [44], including the hybrid quantum-inspired genetic algorithm (HQGA), the hybrid particle swarm optimization algorithm (HPSO), the hybrid simulated annealing algorithm (HSA), the multiobjective evolutionary algorithm based on decomposition (MOEA/D), and NSGA-II. The parameters of the above algorithms are taken from Liu et al. [44]. The performance of an EMO is reflected by the obtained Pareto front’s proximity and diversity. We consider two metrics for performance measurement. The IGD (inverted generational distance) mainly evaluates the proximity of a Pareto front to the optimal reference set, which consists of nondominated solutions in the grouped Pareto fronts. The smaller the IGD metric is, the better proximity a Pareto front has. And the HV (hypervolume) evaluates the maximal area dominated by a Pareto front, which estimates both the convergence and diversity. The larger the value of HV means the better the integrated performance. As is shown in Figure 1 (Pareto front of a problem instance for , , and ), the proposed MOMA clearly dominates the rest of the algorithms. This result is similar for the rest of the problem instances. The dominance relationship is further demonstrated in Tables 4 and 5. MOMA’s smallest IGD values and largest HV values of all instances show that it comprehensively outperforms the other 5 algorithms.

Finally, we examine the impact of choice of parameters, , , and , on algorithm performances. We perform a -way ANOVA using the commercial software SPSS version 17.0. The impacts of , , and on the FLTPA heuristic’s error, MOMA’s IGD metric, and MOMA’s HV metric are shown in Tables 68, respectively. Throughout the results, when the significance level (the last column in the tables) is under 0.01, the impact is deemed as significant. We observe that is significant for the FLTPA heuristic’s performance. The error percentage is much smaller when jobs’ release dates are smaller. This could be explained by smaller release dates that do not tightly restrict the solution space. We further find out that all parameters are not significant for MOMA’s IGD metric, whereas , , and all have a significant impact on MOMA’s HV metric. As and/or change, the hypervolume of a Pareto front would change dramatically as large numerical values are involved.

6. Summary and Future Research

In this study, we considered a bicriterion makespan and total completion time minimization flow shop scheduling with a learning effect subject to release dates. We propose a branch-and-bound algorithm and a heuristic algorithm. Numerical studies showed that the FLTPA heuristic was shown to perform well in obtaining near-optimal solutions. Future research may focus on considering the other bicriteria flow shop scheduling problems with a learning effect subject to release dates or studying more efficient heuristic algorithms.

Data Availability

Anyone who is interested in obtaining the data underlying the findings can send emails to Professor Jian Xu for help.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (no. 71471120, 71872034) and the Support Program for Innovative Talents in Liaoning University (LR2016017). The Liaoning BaiQianWan Talents Program also sponsored this research.