Abstract

In recent years, various customer order scheduling (OS) models can be found in numerous manufacturing and service systems in which several designers, who have developed modules independently for several different products, convene as a product development team, and that team completes a product design only after all the modules have been designed. In real-life situations, a customer order can have some requirements including due dates, weights of jobs, and unequal ready times. Once encountering different ready times, waiting for future order or job arrivals to raise the completeness of a batch is an efficient policy. Meanwhile, the literature releases that few studies have taken unequal ready times into consideration for order scheduling problem. Motivated by this limitation, this study addresses an OS scheduling model with unequal order ready times. The objective function is to find a schedule to optimize the total completion time criterion. To solve this problem for exact solutions, two lower bounds and some properties are first derived to raise the searching power of a branch-and-bound method. For approximate solution, four simulated annealing approaches and four heuristic genetic algorithms are then proposed. At last, several experimental tests and their corresponding statistical outcomes are also reported to examine the performance of all the proposed methods.

1. Introduction

In recent years, various researchers have described customer order scheduling (OS) models (Framinan et al. and [1]). These models can be found in numerous manufacturing and service systems in which several designers, who have developed modules independently for several different products, convene as a product development team, and that team completes a product design only after all the modules have been designed. A practical example is that a manufacturer produces semifinished lenses and competes globally provided by Ahmadi et al. [2].

Regarding the literature on order scheduling with order ready times, only Leung et al. [3] discussed and analyzed the order ready times in a customer order scheduling model with regard to minimizing the total weighted completion time. Garg et al. [1] analyzed the complexities of order schedules to optimize the weighted flow time or offline weighted completion time. Therefore, in this paper we investigate a multiple-machine OS model with different order ready times in which the goal is to find a schedule to minimize the total completion time of orders. Regarding the importance of ready times in wafer fabrication, readers may refer to Monch et al. [4]. Moreover, the jobs or orders to be processed may have different priorities such as the due dates or weights or ready times. In a model with different ready times, waiting for future order or job arrivals to raise the completeness of a batch is an efficient policy.

Within the relevant literature on taking the total weighted completion times into consideration, it has been shown as an NP-hard by Sung and Yoon [5]. Some effective heuristics were developed by Wang and Cheng [6] and Leung et al. [3, 7, 8] for approximate solutions, and a branch-and-bound algorithm was proposed by Yoon and Sung [9] for the exact solution. For the total completion time criterion minimization, Ahmadi et al. [2] showed that it is NP-hard and built some heuristics accordingly. More recently, Framinan and Perez-Gonzalez [10] studied an OS model with completion time objective but zero ready times. They proposed a constructive heuristic with the same complexity compared to the ECT. The key element of this heuristic lies in its ability to estimate the expected contribution of the unscheduled orders, therefore mitigating the greedy behavior of the ECT. In addition, they proposed two improvement mechanisms that are further embedded into a greedy search procedure (GSA) for finding high-quality solutions. For the OS models with the due dates, readers may refer to Leung et al. [11] for minimizing the maximum lateness and to Lee [12] for minimizing the total tardiness. Extending the problem by Lee [12], Xu et al. [13] assumed that order processing times can be shortened by a position-based learning factor. Framinan and Perez-Gonzalez [14] studied the same problem as Lee [12] did. As extremely fast (negligible time) solutions are required, Framinan and Perez-Gonzalez [14] derived a new constructive heuristic that incorporates a look-ahead mechanism to estimate the objective function at the time that the solution is being built. For the scenarios where longer decision intervals are allowed, they provided a novel matheuristic strategy to provide extremely good solutions. Wu et al. [15] addressed an OS problem with sum of-processing-time-based learning effect to minimize total tardiness of the orders. For other researches on OS models, readers may refer to Yang and Posner [16], Leung et al. [3, 17], and Leung, Li et al. [18] on identical but parallel machine settings. Lin et al. [19] utilized three variants of a particle swarm optimization (PSO) algorithm including a basic PSO, an opposite-based PSO, and a PSO with a declined inertia weight for a OS model with two agents and ready times. The criterion is to minimize the total completion time for the orders from agent A while the total completion time for the orders from agent B has a given upper bound. Wu et al. [20] considered an OS model to minimize the linear sum of the total flowtime and the maximum tardiness. In particular, they proposed a branch-and-bound (B&B) algorithm along with several dominance relations and a lower bound, three modified heuristics, and a hybrid iterated greedy algorithm to solve the problem. For some substantial works on approximation algorithms for , readers may refer to Chekuri et al. [21], Hochbaum and Shmoys [22], Sevastianov and Woeginger [23], Alon et al. [24], and Schuurnman and Woeginger [25]. Readers may refer to Chen et al. [26] for a general overview of approximation techniques. The main contributions of this work are described as follows: an order scheduling problem of introducing ready times for each orders and minimizing the total completion time criterion is studied. A B&B method with four propositions and two lower bounds is derived for finding the exact solution. Four versions of simulated annealing algorithm (SA) and four versions of genetic algorithm (GA) are provided for searching the good-quality approximate solutions. The instance simulation experiments and some statistical methods are utilized to determine the performances for all proposed algorithms.

The rest of this study is organized as follows: Section 2 introduces the notation and formulates the problem. Section 3 provides some dominance propositions and two lower bounds used in a branch-and-bound (B&B) algorithm. Section 4 proposes four variants of a simulated annealing (SA) approach and four heuristic genetic algorithms (GA) for finding approximate solutions. Section 5 presents computational results to evaluate the performance of the eight proposed algorithms. The conclusion and suggestions for future study are provided in Section 6.

2. Problem Statement

This section first introduces some notation that is applied in this study.n: total order numberm: total machine number: machine k, and : schedules of orders where and denote partial sequences of the orders: processing time for order on machine k, ; : ready time for order i, : completion time for the last order in on machine k, k = 1, 2, …, m: completion time of order in S, i.e., is the order which is arranged in the lth position in a given schedule.

The problem is formulated as follows. Suppose that orders come in from different clients at different times (ready times), and each order comprises components. Let a facility with different machines be operated in parallel form. A machine can produce its own item; that is, an item can be done on a dedicated machine. If each order has its own ready time, then unforced idleness is also allowed. Preemptions are not forbidden. Let denote the processing time for component of order to be processed on machine k and let denote order to be placed in the rth position in a given sequence. The objective is to seek a sequence to let the total completion time criterion be minimized. To find that the optimal sequence is equivalent to finding a schedule that optimizes the objective function , an illustrative example follows.

Example 1. Given an instance problem with n = 2, m = 2. The processing times of each component of order , i = 1, 2, on machine , and the ready times of order , i = 1, 2, are provided in Table 1.
Given an order schedule S = (O1, O2), according to the definition of , one have the following computation, and Figure 1 shows the Gantt diagram.

3. Propositions and Lower Bound

Branch-and-bound (B&B) has received a considerable amount of attention for solving numerous combinatorial problems (Pinedo 2008; [2729]). On the basis of the proof by Lenstra et al. [30] that the problem is strongly NP-hard when each order has one component only, we know our proposed problem is also an NP-hard problem. Garg et al. [1] also demonstrated that even the unweighted completion time problem with a single ready time is NP-hard. Thus, two lower bounds and certain dominant properties can be applied in a B&B algorithm to find an optimal solution. Following that, four variants of a simulated annealing and four heuristic genetic algorithms are proposed for finding near-optimal solutions.

Let and denote two schedules, in which and are partial sequences. Assume that and are the completion times of orders and in S, whereas and denote the completion times of orders and in . To demonstrate that dominates , it suffices to show that and .

Proposition 2. If , , and , then dominates .

Proof. From , we haveFrom and , we have andBecause and , we haveTherefore, dominates .

Proposition 3. If and , then dominates .

Proposition 4. If and , then dominates .

Proposition 5. If and for , then dominates .

The proofs of Propositions 35 are omitted because they are similar to that of Proposition 2.

Next, a lower bound is another element to raise the efficiency of a B&B algorithm; the algorithm can be expedited by eliminating the partial nodes. In the following, two lower bounds will be proposed for the B&B method. Let be a partial sequence, in which the order of the first jobs is determined and let be the unscheduled part consisting of = (n-v) jobs. Furthermore, let be the starting time of the first order in US for machine k, and . Following the definition of (1), we have

From the above discussions, we haveRegarding the second right-hand term of (8), according to Hardy et al. [31], the value of may be optimized when the orders and are ordered in opposite ways. Therefore, the first lower bound (LB1) can be obtained as follows.where for i = 1,2,.., n1, denotes the nondecreasing form of , and .

Adopting similar steps, another lower bound is easily derived as follows.where for the remaining orders in US.

To tighten the lower bound, the maximum value from (9) and (10) can be set as the lower bound of AS; i.e.,

4. Four Variants of SA and Four GA Algorithms

Metaheuristics have been widely utilized to numerous discrete or continuous NP-hard dynamic combinatorial models. For more examples, readers may refer to Kirkpatrick et al. [32] and Holland [33], among other works. Because the problem under study is considerably difficult, to minimize the solution time cost increased by the solving process, we developed four simulated annealing (SA) approaches and four genetic algorithms (GAs) for near-optimal solutions. Notably, both SA and GA are not simple heuristics; they are metaheuristics. This is crucial because heuristics are usually simple algorithms, but they are complex metaheuristics and they are also difficult to build and carry when comparing with simple heuristics.

The SA is one of the most common metaheuristics to apply for many discrete or continuous optimization problems (see [32, 34]). The main procedures of SA are described in the following.

The main steps of SA procedure:Do while (iteration time < 500n)(i)Input three initial sequences in SA.(ii)For each initial sequence, improve it with the pairwise interchange neighborhood generation method.(iii)Determine if a new sequence is preferable to the original sequence; if the new sequence is preferable, then replace the original sequence with the new sequence; however, if the new sequence is not preferable to the original sequence, the original sequence is accepted if , where 0<q<1 and , , and is an experimental constant.End do

GA widely used another metaheuristic using an intelligent random searching policy to solve numerous difficult NP-hard problems. The important steps are including a structure, a population, a fitness function, reproduction mechanism, and a crossover operator to generate offspring. For more details, readers may refer to Holland [33] and Essafi et al. [35]. The procedure of the proposed GA is detailed as follows.

The steps of GA procedure:Do while (generation sizes < 500)(i)Input three initial sequences to generate N = 20 populations in GA.(ii)Compute the total completion time for each population and then calculate fitness function values of the strings as , where is the ith string chromosome in the vth generation; determine the selection probability of a schedule as .(iii)Use a linear order crossover in the GA (see Falkenauer and Bouffouix, [36]), and set the crossover rate at 1.(iv)Set the value of mutation () at 0.35.(v)Excluding the ideal schedule with a minimum total completion time of orders, create the other offspring from the parent sequences by using the roulette wheel method.End do

In the present work, to obtain the final solution as quickly as possible, three initial sequences are adopted in SA or GA. In SA_min or GA_min, orders are arranged in nondecreasing sequence of , whereas, in SA_max or GA_max, orders are arranged in nondecreasing sequence of . In SA_mean or GA_mean, orders are arranged in nondecreasing sequence of .

The weighted earliest completion time (WECT) heuristic from Leung et al. [3] is incorporated into SA and GA; they are recoded as SA_wect and GA_wect, respectively. Because the objective of this study is to minimize the total completion function of all orders, instead of the total weighted completion time, we set all . Let denote the set of orders. Initially, place the order returned from in the first position and find its completion time, denoted as . Next, delete the scheduled order from N, and place the order returned from in the second position. Repeat this process similarly for all orders. This is called the WECT method. For more details, readers may refer to Leung et al. [3].

5. Computational Experiments

Several experimental tests were conducted to determine the performance of the B&B, the four SAs, and four GAs. Those proposed algorithms are written in FORTRAN.

According to the scheme of Lee [12] and Leung et al. [3, 8, 17], the order processing times are created from the uniform distribution U (1, 100) in every case. According to the scheme of Reeves [37], the order ready times are created from another uniform distribution U (0, 100n) where is the order-size and was a control variable. The value of was set at 0.25, 0.5, and 0.75.

To determine the B&B performance, we report the mean and maximum numbers of nodes and the mean and maximum CPU times (in seconds), respectively. For the variants of SA or GA algorithms based on small numbers of orders (n = 12, 14, or 16), we report the mean and maximum error percentages in which the average error percentage (AEP) was defined as , where was the value of the total completion time yielded from SA or GA and was the value of total completion time yielded from the B&B algorithm.

For the variants of SA or GA algorithms based on large numbers of orders (n = 50 or 100), we recorded the mean and maximum relative percentages of deviation (RPDs) in which the RPD was defined as , where was the value of total completion time yielded from SA or GA and was the smallest objective value selected from among the four SAs and the four GAs.

The two parts of the computational experiment involved small numbers of orders and large numbers of orders.

For the small numbers of orders, we set the numbers of machines at m = 2, 3, and 4. The number of orders was tested at n = 12, 14, and 16. A set of 100 instances was randomly generated for each case.

In the experiment, the B&B was set to skip to the next data set if the number of nodes exceeded 108. The computational results of the B&B method are presented in Table 2, where the last column FS reports the number of instances for which the B&B achieved the optimal solution. The results from Table 2 can be examined in detail; the changes of number of nodes obtained by the B&B method with variations in the parameters n, m, and λ are illustrated in Tables 3 and 4. From Table 3, it can be found that the number of nodes tends to decrease for each fixed value of as the value of λ increases. On the contrary, it appears from Table 4 that the number of nodes tends to increase for each fixed value of as the value of increases except for n=16.

The computational results obtained by the four SAs and four GAs are given in Table 5, and the AEPs are obtained by these algorithms as the parameters n, m, and λ vary. From Table 5, it can be seen that the AEP values obtained by all algorithms tend to linearly decrease as λ increases for n = 12 and n = 14. For n = 16, the AEP values obtained by these algorithms first increase but then decrease as λ increases. Table 5 also indicates that the variation range of AEP obtained by each algorithm is relatively small and the results reveal no trace concerning the impact of on the AEP.

Regarding the performance of these algorithms, the box plots in Figure 6 display the AEP distributions of the eight algorithms. To compare the statistically significant differences in solution quality levels for solutions from the four SAs and the four GAs, we applied analysis of variance (ANOVA); the p value was less than 0.05, as seen in Table 6. We then employed SAS package software to do Fisher’s least significant difference (LSD) tests; the outcomes are shown in Table 7. At the 0.05 level of significance, the AEP means were divided into three groups. The AEPs of SA-mean and SA-max were in the smallest (optimal) group, and AEP of the GA-wect was the largest (worst), whereas the AEPs of the other five algorithms were in the middle group. The SA-max and SA-mean algorithms can be recommended for small numbers of orders.

In experiments with large numbers of orders, we set the numbers of orders at n = 50 and 100 and the numbers of machines at m = 5, 10, and 20. One hundred instances were randomly created for each case.

The performance results of the four SAs and the four GAs are presented in Tables 10 and 11, and the RPDs and corresponding CPU times obtained by these algorithms as the parameters n, m, and λ varied are illustrated in Figures 25. Concerning the RPD metric, no pattern is evident in Figure 2 regarding the impact of λ on the performance levels of these SAs and GAs. The reason is that for n = 50 the RPD slightly increases (except for SA-mean) as λ increases, whereas for n = 100 the RPD slightly decreases (GA-max first increases but then decreases) as λ increases. From Figure 3, it is found that for n = 50 the RPD obtained by each algorithm first increases but then decreases quickly as increases. For n = 100, the RPD values obtained by SAs remain unchanged but the RPD values obtained by GAs tend to decrease linearly as increases. In addition, from Figures 2 and 3, it can be concluded that the performance levels of SAs are superior to those of GAs. Regarding the CPU time metric, Figures 4 and 5 indicate that the values of λ have little impact on the CPU times used by SAs and GAs; however, as increases, the CPU times used by these algorithms quickly increase.

As shown in Figure 7, the box plots of RPDs illustrate the performance levels of the eight algorithms for large numbers of orders. It is obvious that the four SA algorithms outperformed the GAs, and the GA-wect performed the worst. However, another ANOVA, shown in Table 8, and another set of Fisher’s least significant difference tests, shown in Table 9, indicate that, at the 0.05 significance level, there are no differences among the algorithms except for GA-wect. Furthermore, the box plots in Figure 7 also indicate that the RPDs of the four SAs have less dispersion; thus, the SA algorithms solve these problems with both accuracy and stability when the numbers of orders are large.

6. Conclusions

Multifacility order scheduling model has become a vital challenge in numerous practical service and manufacturing situations. Unlike previous studies, which have assumed that all orders are available at time zero, this study investigates orders with different ready times and provides several contributions. First, two lower bounds and several dominant properties are developed, based on which a B&B algorithm is proposed to solve the problem with a small number of orders. Second, four SA variants and four GA variants are presented to obtain near-optimal solutions for problems with large numbers of orders. Computational experiments based on extensive test problems were conducted to evaluate the performance levels of the proposed B&B algorithm and the various SAs and GAs, as well as the impacts of parameters on their performance levels. The experimental outcomes reported that, with the help of the lower bounds and properties, the B&B can tackle problem instances up to n = 16 within a considerable time. In addition, the tested outcomes also illustrate that the four SA variants and the four GA variants deliver satisfactory performance in terms of efficacy and efficiency for problem instances with both small and large numbers of orders. Based on the outcomes obtained by the algorithms, it is recommended that SA-mean or SA-max be utilized to solve the proposed problem because they both outperform GAs and can quickly achieve favorable solutions with a very small value of AEPs. It can be also seen that the problem scale that can be solved by the proposed B&B algorithm only up to n = 16 is limited. The dominance relations or lower bounds of the B&B method appear to be not strong or tight enough. A future study may develop more powerful dominance relations or lower bounds.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This paper was supported in part by Taiwan’s Ministry of Science and Technology under Grant no. MOST105-2221-E-035-053-MY3.