- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Mathematical Problems in Engineering

Volume 2012 (2012), Article ID 928312, 20 pages

http://dx.doi.org/10.1155/2012/928312

## Variable Neighborhood Search for Parallel Machines Scheduling Problem with Step Deteriorating Jobs

^{1}School of Mechanical Engineering, Southwest Jiaotong University, Chengdu 610031, China^{2}School of Mechanical Engineering and Automation, Xihua University, Chengdu 610039, China

Received 20 February 2012; Revised 27 April 2012; Accepted 8 May 2012

Academic Editor: Jung-Fa Tsai

Copyright © 2012 Wenming Cheng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

In many real scheduling environments, a job processed later needs longer time than the same job when it starts earlier. This phenomenon is known as scheduling with deteriorating jobs to many industrial applications. In this paper, we study a scheduling problem of minimizing the total completion time on identical parallel machines where the processing time of a job is a step function of its starting time and a deteriorating date that is individual to all jobs. Firstly, a mixed integer programming model is presented for the problem. And then, a modified weight-combination search algorithm and a variable neighborhood search are employed to yield optimal or near-optimal schedule. To evaluate the performance of the proposed algorithms, computational experiments are performed on randomly generated test instances. Finally, computational results show that the proposed approaches obtain near-optimal solutions in a reasonable computational time even for large-sized problems.

#### 1. Introduction

Scheduling is a form of decision making that plays a crucial role in manufacturing and service systems. It began to be taken seriously in manufacturing at the beginning of 20th century, and since then has been received the attention of many researchers for years. In the traditional scheduling problems, most research in the literature was usually conducted under the assumption that the processing time of a job is known in advance and remains constant throughout the whole decision process. However, there are many practical situations where the processing times of jobs are not constant but increasing overtime such that the later a job starts, the longer it takes to process. This phenomenon is known as scheduling with deteriorating jobs to many industrial applications, such as equipment maintenance, steel production, medical emergency, firefighting, and other problems. In these cases, the corresponding scheduling problem was first studied by Browne and Yechiali [1] and J. N. D. Gupta and S. K. Gupta [2]. They assumed that the processing time of all jobs is a function of their starting point (), where is the normal processing time, is a deterioration rate, and is the starting time of job . From then onwards, more papers considering scheduling jobs with deterioration effects have been published recently [3–19]. Comprehensive reviews and discussions of different models and problems with time-dependent processing times were given by [20, 21].

Most of the current researches in deteriorating jobs scheduling problems assume that the job processing times are linear functions of their starting times. However, in many practical situations, if some jobs fail to be processed prior to a prespecified deteriorating date, then the jobs will require extra time for successful completion. As a result, the original schedule may become inapplicable under the new environment. This motivates the research of the scheduling problem with piecewise-deteriorating jobs. There are often two types of scheduling problems relevant to the piecewise-deteriorating model in the literature. The first one is step-deterioration scheduling problem [22]. Another one is piecewise-linear-deteriorating scheduling problem [23]. In this paper, we study the first kind of scheduling problem. Considering a set of jobs with deteriorating date for job , the processing time for job is if it is started before the deteriorating date and if the starting time , where is the normal processing time, is a deteriorating penalty. The objective is to schedule the jobs on parallel machines to minimize the sum of job completion times. The scheduling problem with step-deterioration jobs is NP-complete even if in single machine environment [24]. As for such a problem over a certain size, it is handicapped to render an optimal schedule within a reasonable run time, depending on exact methods. Hence, some approximate solution techniques are employed to yield an optimal or near-optimal schedule for such problems.

In this paper, a modified weighted combination search algorithm (MWCSA) is firstly proposed to obtain near optimal solution. Subsequently, variable neighborhood search (VNS) algorithm is employed to yield the better schedule for the problem under consideration. VNS is one of the modern metaheuristic techniques that use systematic changes of the neighborhood structure within a local search to solve optimization problem [25, 26]. Owing to incorporating a lot desirable properties for a metaheuristics such as simplicity, efficiency, effectiveness, generality, and so forth, VNS has been widely used to combinatorial optimization problems in recent years. In addition, to improve the performance of VNS, the proposed heuristic MWCSA is utilized to produce initial solution for the VNS. The effectiveness of the proposed approaches is demonstrated by computational results based on a large set of randomly generated test instances.

This paper is organized as follows: Section 2 reviews the relevant literature in this area of scheduling. Section 3 gives a brief description of the problem under consideration and formulates the corresponding mixed integer programming model. A detailed description of the heuristic MWCSA and VNS algorithm are elaborated in Section 4. The performances of the VNS and the heuristic algorithms are shown in Section 5. Conclusions are given in Section 6.

#### 2. Literature Review

We discuss related work with respect to the scheduling problem with piecewise-deteriorating jobs and with regard to VNS methods in scheduling, especially for parallel machines.

Since many different problems are encountered in production and service environments where the processing time of jobs or tasks can be modeled as the piecewise-deteriorating function, many active researchers in the field of scheduling have perceived the promising advantages of considering the piecewise-deteriorating jobs in scheduling problems. Kunnathur and Gupta [27] firstly assumed that job processing times are piecewise-linear function of their starting times with two pieces and addressed five heuristics and two optimizing algorithms based on dynamic programming and branch and bound techniques. However, Sundararaghavan and Kunnathur [22] considered that the processing time can be modeled by step function and gave some solvable cases for minimizing the sum of the weighted completion times in single machine scheduling. In addition, Mosheiov [28] studied the problem of stepwise deteriorating jobs for minimizing makespan, proved that the case of single machine step-deterioration is NP-complete, and suggested heuristic methods in single and multimachine scheduling environments.

Cheng and Ding [24] studied a piecewise model where the job processing time deteriorates as a step function if its starting time exceeds a given threshold, and presented NP-complete proofs for minimizing makespan, flow time, and weighted flow time in single machine scheduling. Jeng and Lin [29] presented a pseudopolynomial time dynamic programming algorithm for a single-machine scheduling problem. Then, they described an efficient branch and bound algorithm based on two dominance rules and a lower bound for deriving optimal solution. Subsequently, a single machine scheduling problem of minimizing the total completion time was solved by branch and bound algorithms [30]. Owing to the complexity of the problem, He et al. [31] proposed a weight combination search algorithm to yield a near-optimal solution. They also developed a branch and bound algorithm incorporated with several properties and a lower bound to obtain the optimal solution. Afterwards, Layegh et al. [32] applied a memetic algorithm based on three dominance properties to minimize the total weighted completion time on a single machine under step-deterioration.

With regards to the piecewise-linear-deterioration model, Kovalyov and Kubiak [33] investigated that the job processing time can be expressed by a piecewise-linear-increasing function with three pieces and presented a fully polynomial approximation scheme for minimizing makespan of single machine scheduling problem. Then Kubiak and van de Velde [23] proved that the problem is NP-hard, and proposed a branch and bound algorithm that solves instances with up to 100 jobs in a reasonable amount of time. Alternatively, the processing job times was expressed in a piecewise-linear-decreasing function of their start times. The single machine scheduling problem was researched by Cheng et al. [34]. Moslehi and Jafari [35] dealt with the same problem proposed by [33] with the minimization of the number of tardy jobs. Owing to the complexity of the studied problem, they proposed a heuristic algorithm with and a branch and bound algorithm.

As far as reviewed, there are some literature focused on single machine scheduling problem with piecewise-deteriorating jobs. However, to the best of our knowledge, parallel machines scheduling problem with piecewise deteriorating jobs has not been considered in the existing literature. And it still lacks some effective approaches to solve the intractable problem, especially metaheuristics. Therefore, we intend to propose a VNS approach for parallel machines scheduling problem with step-deteriorating jobs.

VNS approaches have successfully applied to solve several scheduling problems. Gupta and Smith [36] described a new hybrid of VNS for single machine total tardiness scheduling with sequence-dependence setups. Paula et al. [37] applied a VNS algorithm to solve large instances of parallel machines scheduling problems with sequence-dependent times. Anghinolfi and Paolucci [38] contributed to the design of a hybrid meta-heuristic approach which integrates several features from tabu search (TS), simulated annealing (SA), and VNS for a generalized parallel machine total tardiness scheduling problem. Moreover, Driessel and Mönch [39] proposed several variants of VNS schemes for the same problem considered precedence constraints to minimize the total weighted tardiness of all jobs. Behnamian et al. [40] proposed a hybrid meta-heuristic method which combines some advantages of ant colony optimization, SA and VNS for the minimization of makespan in parallel machine scheduling problem with sequence-dependent times. C.-L. Chen and C.-L. Chen [41] proposed several hybrid metaheuristics for unrelated parallel machine scheduling with sequence dependent setup times by integrating VNS and TS principles. Furthermore, some job shop scheduling problems were solved by VNS [42–45]. A latest systemic survey of state-of-the-art development of VNS can be found in [46].

#### 3. Problem Formulation

The problem under consideration is to schedule independent jobs with step-deterioration effects, noncommon deadlines and varying processing times on identical parallel machines. The jobs are available for processing at time zero and the machines are available in the whole process. No job preemption is permitted. For each job , there is a normal processing time and associated with a deteriorating date . If the starting point of its process is less than or equal to its deteriorating date, then it only requires a normal processing time . Otherwise, it requires an extra processing time , which is called the deteriorating penalty. Thus, the actual processing time of job depends on its starting time and deteriorating date , and can be defined as a step-function: if ; + , otherwise. Without loss of generality, it is assumed that parameters , , and are all integers. Let denote the completion time of job . The goal of the problem is to find a schedule such that the total completion time, or the sum of completion times, of all jobs is minimized. This problem is denoted as or by adopting the standard three-field notation.

For convenience, a job is called early if its starting time is before or at its deteriorating date; tardy, otherwise. For schedule *S*, the objective value is denoted by *Z*(*S*). In addition, let be a sufficiently large positive number. Before formulating the proposed mixed integer programming model, the following variables have to be defined. : Binary, set to 1 if job is immediately followed by job in sequence on machine *k*; 0, otherwise () : Binary, set to 1 if job is assigned to machine *; 0, *otherwise ().

Based on the definitions and notation described above, the considered problem can now be formulated as 0-1 mixed integer programming model, as shown below. Subject to:

In the above mathematical model, the objective (3.1) minimizes the total completion time. Constraints (3.2) depict the processing time of the step-deteriorating jobs. Constraints (3.3) and (3.4) ensure that only one job can be processed at the first and the last position on each machine. Constraints (3.5) and (3.6) guarantee that each job is scheduled only once and processed by at most one machine. Constraints (3.7) define the completion time of the first job assigned to a machine. Constraints (3.8) represent that the completion time of a job in sequence on each machine will be at least equal to the sum of the completion time of the preceding job and the processing time of the present job, if the job is immediately schedule after the previous job. Constraints (3.9) define the starting time of the first job assigned to a machine. Constraints (3.10) state that the starting time of a job in sequence on each machine is greater than or equal to the completion time of the preceding job. Constraints (3.11) confirm that each job is only processed to exactly one machine. Constraints (3.12) specify that the decision variable and is binary over all domains.

This described problem is NP-complete because the single machine problem or has proven to be NP-complete by [24]. There exists no polynomial time algorithm for the exact solution of the considered problem. Therefore, we have to look for some efficient heuristic or meta-heuristic algorithms to yield near optimal solutions of large-sized problems in a reasonable computational time.

#### 4. Heuristic Solution Approaches

In this section, we start with presenting one property. Then we modify a heuristic proposed in a literature. Subsequently, some detailed descriptions of VNS meta-heuristic are given.

Since the identical parallel machines scheduling problem is solvable in *O*(*N*log*N*) time using the generalized shortest processing time (GSPT) rule [47]. Therefore, we can also use the same method to determine the sequence of early jobs and tardy jobs assigned to a given machine for minimizing the total completion time. Based on the GSPT, we obtain the following property directly.

*Property 1. * On each machine of an optimal solution for the considered problem, early jobs and tardy jobs are sequenced in the nondecreasing order of and , respectively.

*Proof. *The proof is straightforward from the SPT rule on and , respectively.

##### 4.1. Modified Weight Combination Search Approach

He et al. [31] proposed weight combination search approach (WCSA) to solve single machine total completion time scheduling problem with step-deteriorating jobs. The performance of the WCSA is relatively excellent compared to the proposed branch and bound algorithm within 24 jobs. Based on the good performance of the algorithm, it is modified to adapt the parallel machine model by incorporating Property 1. But the ranges of the weights in the literature are not suitable for parallel machine scheduling. For specifying the suitable ranges of the weights we carried out the preliminary tests of 30 10-job instances and found that the linear combination of , , and will acquire better results if and . The procedure of the modified WCSA (MWCSA) is designed as show in Algorithm 1.

From Algorithm 1, all jobs are sorted in non-decreasing order of their normal processing times . The first jobs are assigned to machines, respectively. Then the last jobs are arranged in a non-decreasing order of the weight combination of processing time, deteriorating date and deterioration penalty. Of course, if the completion time of the last scheduled job is greater than the maximal deteriorating date of the unscheduled jobs, then the remaining jobs are indexed in the order of sum of their normal processing times and deterioration penalties.

##### 4.2. Variable Neighborhood Search Metaheuristic

VNS, a new local search technique, attempts to escape from local optimum by exploring more than one type of neighborhood search structure (NSS) during the course of the algorithm. In practice, VNS is very similar to iterated local search (ILS). Owing to iterating over one constant type of neighborhood structure, ILS is easy to stick in local optima: the single move required to improve the solution cannot escape from the found local optima, and even lead to a deterioration of the solution quality. The disadvantage is also existed in some other metaheuristics, such as simulated annealing, and tabu search. However, VNS makes full use of some NSSs until some stopping criterion is met. By exploring various NSSs, we expect VNS to enjoy a systematic diversification mechanism. Besides this diversification mechanism, the other reasons to high acceptability and popularity of VNS among researchers are due to conceptual simplicity to understand and implement and the high flexibility and brilliant adaptability to different problems. Consequently, we intend to use it to solve our problem.

Since the VNS was proposed in 1997, some variants of VNS are developed by Hansen and Mladenović [25]. Variable neighborhood descent (VND), as one of these variants of VNS, performs change of neighborhoods in a deterministic way and adopts the first improvement rather than the best improvement. In the subsequent subsection, we are going to employ a VNS with five powerful NSSs based on insertion, swap, and inversion under the framework of VND.

###### 4.2.1. Encoding, Decoding, and Initialization Schemes

Encoding scheme is a procedure to make a solution recognizable for an algorithm. Permutation list is a common used encoding scheme for combinatorial optimization problems. Since the parallel machines scheduling problem is a combination of machine assignment and operation sequencing decisions, we construct five neighborhood structures that expediently use insertion, swap, and inversion operations based on permutation list. A good initial solution can reduce considerably the computational time. According to observing the processing sequence of jobs, we find that a job with the smallest normal processing time should be processed at first, and a job with the largest penalty should be processed at last. Therefore, we can obtain the initial solution by arranging jobs in the non-decreasing order of the ratio of /. The sample sequencing method is called smallest rate first (SRF). A permutation with all jobs sequence can express the order in which the jobs are processed. To decode a permutation list into a schedule, we use LIST method. For one solution sequence, whenever a machine is freed the job that is chosen from the sequence in succession is put on the machine.

The procedures of encoding and decoding a candidate solution are illustrated according to an example. Consider a problem with 8 jobs and 2 machines. The parameters of all jobs are given in Table 1. The initial solution is obtained by using SRF, as shown in Figure 1. The job that is selected from the solution sequence one by one is processed by the machine with the smallest completion time. The schedule is also depicted in Figure 1. We calculated that the total completion processing time is 1113.

###### 4.2.2. Neighborhood Search Structures

The main purpose of applying a neighborhood structure is to produce a neighboring solution from the current solution via making some changes in it. A variety of neighborhood structures have been applied to scheduling problems. Most of these neighborhood structures are based on insertion or swap operations. In addition, the inversion operation has been reported to be significantly surpassing all the modified forms of crossover of genetic algorithm especially tailored to deal with combinational problems [48]. Therefore, we define five types of neighborhood search structures based on insertion, swap, and inversion.

To generate a neighboring solution, NSS_{1} makes some slight changes in the candidate solution by swapping the positions of all jobs in the sequence with one selected job at random and without repetition. All jobs are selected one after another without repetition in a random order. If we observe the first improvement, the associated sequence is accepted and the procedure restarts. According to the NSS_{1}, some available neighboring searches are found and used to improve the candidate solution. The whole procedure repeats so long as no improvement is obtained through swap all jobs with another randomly selected job. The procedure of NSS_{1} is described in Algorithm 2 in detail. After implementing all swap operations, we believe that there is little hope for further improvement just by swapping two jobs on their positions. Hence, it necessitates consider another neighborhood structure to escape from this local optimum of the NSS_{1}. Subsequently, we need to introduce another NSS to search potential improvement based on insertion moves.

In our NSS_{2}, the procedure is similar to the NSS_{1}. The only difference between them is their move pattern that NSS_{2} changes the position of one job based on insertion neighborhood. A job is removed from the sequence at random and without repetition, and then relocated another random selected position. The other procedure is the same as the NSS_{1}. To avoid duplicated description, it is unnecessary to go into these details here. After relocating all jobs, we feel that the algorithm cannot rely solely on single move to improve the quality of the current solution. Hence, we need to introduce other neighborhood search structure to generate more complex neighbors than the above two structures.

In the NSS_{3}, the number of randomly selected jobs is 2. The manner of choosing these 2 jobs is all the combinations of two-out-of-*n* jobs. Swap the selected jobs on their positions and observe if the solution is improved. Once observing the first improvement, the associated solution is accepted and the procedure restarts. If not, this search strategy repeats for the subsequent combinations. Algorithm 3 illustrated the whole procedure of the NSS_{3}. By the same way of switching from NSS_{1} to NSS_{2}, the double insertions are employed in the fourth neighborhood structure NSS_{4}. In the NSS_{4}, the number of removed jobs is set equal to 2. To relocate them, the two jobs are reinserted into two new randomly selected positions. The other procedure is the same to the NSS_{3}.

To improve further the search performance of the proposed VNS, an inversion operator is embedded into the NSS_{5} when the foregoing neighborhood search methods terminate. In the NSS_{5}, it randomly selects two positions, known as the points of inversion, and inverts the sequence between these positions. The inversion procedure is repeated * φ* times for producing a new job sequence. Owing to the impact of the parameter on neighborhood search, it must be tuned. The inversion operation is unidirectional, that is, only the inverted sequence is improved, the corresponding solution is accepted to replace the incumbent one. The procedure of is described in Algorithm 4.

The general outline of the proposed VNS is shown in Algorithm 5. According to using the above neighborhood search structures, the proposed VNS can be quick to yield the near optimal solution. If the neighborhood search structure continues to increase the number of selected jobs, the quality of the algorithm will be not improved according to our primary experiments.

#### 5. Computational Analysis

To test the performance of our approaches, computational experiments were carried out for the * = * or , problem. More specifically, we also carried out comprehensive computational and statistical tests to assess the performance of MWCSA, and VNS. Firstly, we gave a detailed description of the different instance sets that we have employed. After reporting the results, we carried out comprehensive statistical analyses in order to soundly test the significance of the reported results.

Because there is not a widely available set of benchmark instances for the problem under study, data for the test instances are generated randomly. The normal processing times () are randomly generated from an integer uniform distribution on . The deterioration penalties () are randomly picked from an integer uniform distribution on , where . Let , 0 ≤ *f *≤ 1. The deteriorating dates are randomly selected from three intervals , and . In this paper, we have divided the benchmark instances into 2 different types according to the number of jobs. For small-sized instances, we test all combinations of and . For large-sized instances, we use the combinations of and to evaluate the proposed algorithms. Since the deteriorating date may be generated from three different intervals, we totally have instances.

In order to validate the performance of the proposed approaches, attempts are made to solve the MIP model presented in Section 3 using the software IBM-ILOG CPLEX 12.3 solvers. The software is run on a PC with Intel Pentium Dual-Core 2.60 GHz CPU and 2 GB RAM. Because the problem under consideration is NP-complete, it is impossible to obtain optimal solutions by using some polynomial time algorithms. Therefore, only small-sized instances can be solved to optimality using CPLEX. For the small-sized instances, we use relative percentage deviation (RPD) as a performance measure to compare those results of the CPLEX software, the MWCSA and the VNS. RPD is obtained by given formula below:
where *Z*(Alg) is the total completion time of the solution obtained by a given algorithm and instance and *Z*(OPL) is the total completion time of the optimal schedule given by CPLEX. The SRF (small rate first) of Section 4 is also brought into the large-sized instances for comparison. Our purpose to utilize SRF is to use its result as upper bound for a given instance to evaluate the VNS. In addition, the MWCSA is used to obtain the initial solution of the VNS for improving the quality of solutions. The hybrid algorithm, denoted by VNS + MWCSA, is employed for the large-sized instances.

All proposed algorithms in the previous sections were coded in MATLAB 7.11 and implement on the same PC. According to the preliminary tests, the stopping criterion used when testing all instances with the proposed algorithms is set to the maximum iterations time fixed to 200. Thus, there is only one parameter () is tuned in our proposed VNS. The considered candidate levels are 10, 30, 50, and 70. A set of 15 large-sized instances are randomly generated when the number of jobs is fixed to 100. All the 15 instances are solved by these VNSs with different parameter candidate levels. The results were analyzed by the one-way ANOVA test. The means plot and Fisher LSD intervals for the 4 levels of parameter are shown in Figure 2. As we can see, of 50 and 70 provide the better results among all levels. However, it needs more computational time when is 70. Therefore, the most reasonable value of parameter is 50.

CPLEX, SRF, and MWCSA are deterministic and only one run is necessary. The VNS and VNS + MWCSA are stochastic and we need to run some replicates in order to better assess the results. We run each algorithm five times. The listed objective value and computational time are the means of five reported results.

Table 2 provides a comparison of the solutions to the mixed integer programming model generated by the CPLEX and those solutions provided by execution of the proposed MWCSA and VNS algorithms. All small-sized instances were solved to optimality by CPLEX as shown in Table 1. The computational time of the CPLEX increases rapidly as the instance become larger. When the number of jobs is equal to 10, solutions to MIP models require far longer than 1 hour. From Table 1, the longest computational time of CPLEX is 144761.59 seconds (40.21 hours). It is impractical in scheduling for such a long time. In addition, we found that CPLEX is not greatly affected by the distribution of the deteriorating dates almost. While the deteriorating dates are randomly picked from the interval , all corresponding instances are relatively easy. This is because the number of deteriorating jobs is less. It is worthwhile to note that the MWCSA behaves well for all small-sized instances and gives an average result of 0.62 deviation from the optimal solution. Specially, eleven of these instances can be solved optimally by the proposed heuristic. Its execution time is far less than 1 second so that its log is not necessary. It can also be noticed that the solutions from the VNS algorithm are optimal for these small-sized instances. Meanwhile, the computational time of these instances given by VNS are not more than 1 second. This means that the VNS algorithm is effective in solving the scheduling problem under study and has significantly better optimization performance than the others.

Tables 3, 4, and 5 show the results of large-sized instances where the deteriorating dates are generated from three different intervals. In order to compare the performance of the proposed algorithms, the ratios of the results of a given algorithm to the values obtained by SRF were calculated out. For the instances with the deteriorating dates generated from interval , VNS behaves slightly better than MWCSA. The difference in average rate between two algorithms is only 2.26%. For other two types of instances with intervals and , the VNS and VNS + MWCSA strikingly outperformed the MWCSA heuristic because of its larger search space. As a rule of thumb, a given algorithm takes possession of the better initial solution and should yield better result. However, the hybrid algorithm, with longer run times, does not perform noticeably better than the VNS. It is chiefly because the VNS carries out enough neighborhood searches before the stopping criterion meets. The rates of three different types of instances were plotted into smooth curves shown in Figures 3, 4, and 5. It can be seen from those figures that the performance of MWCSA heuristic is not particularly robust as regards the distribution of the deteriorating date. Alternatively, VNS and VNS + MWCSA are statistically better than MWCSA for three types of instances.

#### 6. Conclusions

This paper considers the identical parallel machines scheduling problem with step-deteriorating jobs. The processing time of each job is a step function of its starting time and a deteriorating date that is individual to all jobs. The problem is to determine the allocation of jobs to machines as well as the sequence of the jobs assigned to each machine for the criteria of minimizing the total completion time. A mathematical model for this problem has been formulated. Since the problem under study is NP-complete, it is impossible to solve large-sized instances to optimality. To solve the tackled problem, a heuristic MWCSA and a VNS are proposed to obtain the near optimal solutions. In order to further improve the quality of solution, the heuristic MWCSA has been hybridized with the VNS algorithm and implemented to provide a good initial solution. Numerical experiments are conducted on small- and large-sized instances. Computational results show that MWCSA produces some good solutions compared to CPLEX, but the performance is greatly affected by the distribution of the deteriorating date. In contrast, VNS and VNS + WMCSA are robust as regards three types of instances. Therefore, fairly good solutions can be obtained by the proposed methods within reasonable amount of time.

For further study, it is worth considering the setup times between jobs for the problem of scheduling jobs with piecewise-deterioration on multimachine. In addition, other efficient constructive heuristics and neighborhood properties are worthwhile to investigate for improving the quality of solutions.

#### Acknowledgments

This work is partially supported by the National Natural Science Foundation of China under Grant no. 51175442, the Youth Foundation for Humanities and Social Sciences of Ministry of Education of China under Grant no. 12YJCZH296, the Ph.D. Programs Foundation of Ministry of Education of China under Grant no. 200806131014, and the Fundamental Research Funds for the Central Universities under Grant nos. SWJTU09CX022, 2010ZT03.

#### References

- S. Browne and U. Yechiali, “Scheduling deteriorating jobs on a single processor,”
*Operations Research*, vol. 38, no. 3, pp. 495–498, 1990. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - J. N. D. Gupta and S. K. Gupta, “Single facility scheduling with nonlinear processing times,”
*Computers and Industrial Engineering*, vol. 14, no. 4, pp. 387–393, 1988. View at Publisher · View at Google Scholar · View at Scopus - A. Bachman and A. Janiak, “Minimizing maximum lateness under linear deterioration,”
*European Journal of Operational Research*, vol. 126, no. 3, pp. 557–566, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - A. Bachman, A. Janiak, and M. Y. Kovalyov, “Minimizing the total weighted completion time of deteriorating jobs,”
*Information Processing Letters*, vol. 81, no. 2, pp. 81–84, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - C.-C. Wu, W.-C. Lee, and Y.-R. Shiau, “Minimizing the total weighted completion time on a single machine under linear deterioration,”
*International Journal of Advanced Manufacturing Technology*, vol. 33, no. 11-12, pp. 1237–1243, 2007. View at Publisher · View at Google Scholar · View at Scopus - J.-B. Wang, L. Lin, and F. Shan, “Single-machine group scheduling problems with deteriorating jobs,”
*International Journal of Advanced Manufacturing Technology*, vol. 39, no. 7-8, pp. 808–812, 2008. View at Publisher · View at Google Scholar · View at Scopus - J.-B. Wang, C. T. Ng, and T. C. E. Cheng, “Single-machine scheduling with deteriorating jobs under a series-parallel graph constraint,”
*Computers & Operations Research*, vol. 35, no. 8, pp. 2684–2693, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - M. M. Mazdeh, F. Zaerpour, A. Zareei, and A. Hajinezhad, “Parallel machines scheduling to minimize job tardiness and machine deteriorating cost with deteriorating jobs,”
*Applied Mathematical Modelling*, vol. 34, no. 6, pp. 1498–1510, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - L. Sun, L. Sun, K. Cui, and J. B. Wang, “A note on flow shop scheduling problems with deteriorating jobs on no-idle dominant machines,”
*European Journal of Operational Research*, vol. 200, no. 1, pp. 309–311, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - D. Wang and J. B. Wang, “Single-machine scheduling with simple linear deterioration to minimize earliness penalties,”
*International Journal of Advanced Manufacturing Technology*, vol. 46, no. 1–4, pp. 285–290, 2010. View at Publisher · View at Google Scholar · View at Scopus - J.-B. Wang, “Flow shop scheduling with deteriorating jobs under dominating machines to minimize makespan,”
*International Journal of Advanced Manufacturing Technology*, vol. 48, no. 5–8, pp. 719–723, 2010. View at Publisher · View at Google Scholar · View at Scopus - C.-M. Wei and J.-B. Wang, “Single machine quadratic penalty function scheduling with deteriorating jobs and group technology,”
*Applied Mathematical Modelling*, vol. 34, no. 11, pp. 3642–3647, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - X. Huang and M.-Z. Wang, “Parallel identical machines scheduling with deteriorating jobs and total absolute differences penalties,”
*Applied Mathematical Modelling*, vol. 35, no. 3, pp. 1349–1353, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J.-B. Wang and C. Wang, “Single-machine due-window assignment problem with learning effect and deteriorating jobs,”
*Applied Mathematical Modelling*, vol. 35, no. 8, pp. 4017–4022, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J.-B. Wang, J.-J. Wang, and P. Ji, “Scheduling jobs with chain precedence constraints and deteriorating jobs,”
*Journal of the Operational Research Society*, vol. 62, no. 9, pp. 1765–1770, 2011. View at Publisher · View at Google Scholar · View at Scopus - J.-B. Wang and C.-M. Wei, “Parallel machine scheduling with a deteriorating maintenance activity and total absolute differences penalties,”
*Applied Mathematics and Computation*, vol. 217, no. 20, pp. 8093–8099, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - S.-H. Yang and J.-B. Wang, “Minimizing total weighted completion time in a two-machine flow shop scheduling under simple linear deterioration,”
*Applied Mathematics and Computation*, vol. 217, no. 9, pp. 4819–4826, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J.-B. Wang and M.-Z. Wang, “Single-machine scheduling with nonlinear deterioration,”
*Optimization Letters*, vol. 6, no. 1, pp. 87–98, 2012. View at Publisher · View at Google Scholar - J.-B. Wang, M.-Z. Wang, and P. Ji, “Single machine total completion time minimization scheduling with a time-dependent learning effect and deteriorating jobs,”
*International Journal of Systems Science*, vol. 43, no. 5, pp. 861–868, 2012. View at Publisher · View at Google Scholar · View at Scopus - B. Alidaee and N. K. Womer, “Scheduling with time dependent processing times: review and extensions,”
*Journal of the Operational Research Society*, vol. 50, no. 7, pp. 711–720, 1999. View at Google Scholar · View at Zentralblatt MATH · View at Scopus - T. C. E. Cheng, Q. Ding, and B. M. T. Lin, “A concise survey of scheduling with time-dependent processing times,”
*European Journal of Operational Research*, vol. 152, no. 1, pp. 1–13, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - P. S. Sundararaghavan and A. S. Kunnathur, “Single machine scheduling with start time dependent processing times: some solvable cases,”
*European Journal of Operational Research*, vol. 78, no. 3, pp. 394–403, 1994. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - W. Kubiak and S. van de Velde, “Scheduling deteriorating jobs to minimize makespan,”
*Naval Research Logistics*, vol. 45, no. 5, pp. 511–523, 1998. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - T. C. E. Cheng and Q. Ding, “Single machine scheduling with step-deteriorating processing times,”
*European Journal of Operational Research*, vol. 134, no. 3, pp. 623–630, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - P. Hansen and N. Mladenović, “Variable neighborhood search: principles and applications,”
*European Journal of Operational Research*, vol. 130, no. 3, pp. 449–467, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - N. Mladenović and P. Hansen, “Variable neighborhood search,”
*Computers & Operations Research*, vol. 24, no. 11, pp. 1097–1100, 1997. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - A. S. Kunnathur and S. K. Gupta, “Minimizing the makespan with late start penalties added to processing times in a single facility scheduling problem,”
*European Journal of Operational Research*, vol. 47, no. 1, pp. 56–64, 1990. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - G. Mosheiov, “Scheduling jobs with step-deterioration; Minimizing makespan on a single- and multi-machine,”
*Computers and Industrial Engineering*, vol. 28, no. 4, pp. 869–879, 1995. View at Publisher · View at Google Scholar · View at Scopus - A. A. K. Jeng and B. M. T. Lin, “Makespan minimization in single-machine scheduling with step-deterioration of processing times,”
*Journal of the Operational Research Society*, vol. 55, no. 3, pp. 247–256, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - A. A. K. Jeng and B. M. T. Lin, “Minimizing the total completion time in single-machine scheduling with step-deteriorating jobs,”
*Computers and Operations Research*, vol. 32, no. 3, pp. 521–536, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - C.-C. He, C.-C. Wu, and W.-C. Lee, “Branch-and-bound and weight-combination search algorithms for the total completion time problem with step-deteriorating jobs,”
*Journal of the Operational Research Society*, vol. 60, no. 12, pp. 1759–1766, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - J. Layegh, F. Jolai, and M. S. Amalnik, “A memetic algorithm for minimizing the total weighted completion time on a single machine under step-deterioration,”
*Advances in Engineering Software*, vol. 40, no. 10, pp. 1074–1077, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - M. Y. Kovalyov and W. Kubiak, “A fully polynomial approximation scheme for minimizing makespan of deteriorating jobs,”
*Journal of Heuristics*, vol. 3, no. 4, pp. 287–297, 1998. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - T. C. E. Cheng, Q. Ding, M. Y. Kovalyov, A. Bachman, and A. Janiak, “Scheduling jobs with piecewise linear decreasing processing times,”
*Naval Research Logistics*, vol. 50, no. 6, pp. 531–554, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - G. Moslehi and A. Jafari, “Minimizing the number of tardy jobs under piecewise-linear deterioration,”
*Computers and Industrial Engineering*, vol. 59, no. 4, pp. 573–584, 2010. View at Publisher · View at Google Scholar · View at Scopus - S. R. Gupta and J. S. Smith, “Algorithms for single machine total tardiness scheduling with sequence dependent setups,”
*European Journal of Operational Research*, vol. 175, no. 2, pp. 722–739, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - M. R. D. Paula, M. G. Ravetti, G. R. Mateus, and P. M. Pardalos, “Solving parallel machines scheduling problems with sequence-dependent setup times using variable neighbourhood search,”
*IMA Journal of Management Mathematics*, vol. 18, no. 2, pp. 101–115, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - D. Anghinolfi and M. Paolucci, “Parallel machine total tardiness scheduling with a new hybrid metaheuristic approach,”
*Computers & Operations Research*, vol. 34, no. 11, pp. 3471–3490, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - R. Driessel and L. Mönch, “Variable neighborhood search approaches for scheduling jobs on parallel machines with sequence-dependent setup times, precedence constraints, and ready times,”
*Computers and Industrial Engineering*, vol. 61, no. 2, pp. 336–345, 2011. View at Publisher · View at Google Scholar · View at Scopus - J. Behnamian, M. Zandieh, and S. M. T. F. Ghomi, “Parallel-machine scheduling problems with sequence-dependent setup times using an ACO, SA and VNS hybrid algorithm,”
*Expert Systems with Applications*, vol. 36, no. 6, pp. 9637–9644, 2009. View at Publisher · View at Google Scholar · View at Scopus - C.-L. Chen and C.-L. Chen, “Hybrid metaheuristics for unrelated parallel machine scheduling with sequence-dependent setup times,”
*International Journal of Advanced Manufacturing Technology*, vol. 43, no. 1-2, pp. 161–169, 2009. View at Publisher · View at Google Scholar · View at Scopus - J. Gao, L. Sun, and M. Gen, “A hybrid genetic and variable neighborhood descent algorithm for flexible job shop scheduling problems,”
*Computers & Operations Research*, vol. 35, no. 9, pp. 2892–2907, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - V. Roshanaei, B. Naderi, F. Jolaib, and M. Khalili, “A variable neighborhood search for job shop scheduling with set-up times to minimize makespan,”
*Future Generation Computer Systems*, vol. 25, no. 6, pp. 654–661, 2009. View at Publisher · View at Google Scholar · View at Scopus - M. A. Adibi, M. Zandieh, and M. Amiri, “Multi-objective scheduling of dynamic job shop using variable neighborhood search,”
*Expert Systems with Applications*, vol. 37, no. 1, pp. 282–287, 2010. View at Publisher · View at Google Scholar · View at Scopus - M. Amiri, M. Zandieh, M. Yazdani, and A. Bagheri, “A variable neighbourhood search algorithm for the flexible job-shop scheduling problem,”
*International Journal of Production Research*, vol. 48, no. 19, pp. 5671–5689, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - P. Hansen, N. Mladenović, and J. A. M. Pérez, “Variable neighbourhood search: methods and applications,”
*Annals of Operations Research*, vol. 175, pp. 367–407, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - M. L. Pinedo,
*Scheduling: Theory, Algorithms and Systems*, Springer, New York, NY, USA, 2008. - C. Ferreira, “Combinatorial optimization by gene expression programming: inversion revisited,” in
*Proceedings of the Argentine Symposium on Artificial Intelligence*, pp. 160–174, Santa Fe, NM, USA, 2002.