938
Views
2
CrossRef citations to date
0
Altmetric
Research Article

Fuzzy bi-criteria scheduling on parallel machines involving weighted flow time and maximum tardiness

, & | (Reviewing Editor)

Abstract

This paper considers a bi-criteria scheduling on parallel machines in fuzzy environment which optimizes the weighted flow time and maximum tardiness. The fuzziness, vagueness or uncertainty in processing time of jobs is represented by triangular fuzzy membership function. The bi-criteria problem with weighted flow time and maximum tardiness as primary and secondary criteria, respectively, for any number of parallel machines is NP-hard. So, a heuristic algorithm is proposed to find the optimal sequence of jobs processing on parallel machines so as to minimize the secondary criterion of maximum tardiness without violating the primary criterion of weighted flow time. A numerical illustration is also given in support of the algorithm proposed.

Public Interest Statement

In this bi-objective parallel machine scheduling problem a system of n-independent jobs are to be processed on any of the available m-identical parallel machines. The processing times of jobs are considered to be uncertain or vague in nature due to the presence of some external sources. Fuzzy set theory is used to handle the uncertainty involved. The delay in processing of jobs is called tardiness or lateness and the delayed job is called tardy job. The flow time of a job is the total time it spends in the system i.e. sum of waiting times and processing times. In case when the jobs have different degree of importance, indicated by the weight of the job, the total weighted flow time is the simplest measure of the quality of service received by the jobs. Therefore, the paper seeks to highlight optimization of maximum tardiness and minimum weighted flow time of jobs.

1. Introduction

Scheduling is concerned with the optimal allocation of scare resources with objective of optimizing one or several criteria’s (Baker, Citation1974). Scheduling of jobs has been a fruitful area of research for many decades, in which order of processing of jobs is decided. If the jobs are scheduled properly, it not only saves time but also increases the efficiency of the system. The classical parallel machine scheduling problem includes the scheduling of n- independent jobs say {1, 2, …, n}, on a set of m-identical parallel machines say {1, 2, …, m}. Every job is considered with a deterministic processing time pi, a release date ri, and a due date di while optimizing the objective function as number of tardy jobs, maximum tardiness, weighted flow time, maximum completion time, etc. However, due to the sources of uncertainty and complex nature of the problems in real-time situations, it is quite erratic to set them as precise values. Some researchers work on the concept of probability distribution to model parallel machine problem that is usually predicted from chronological data. But, whenever production data are unreliable or presented in vague way, stochastic models may not be the best choice. Fuzzy set theory may provide an alternative approach for dealing with the vagueness. Zadeh (Citation1965) was first to introduce fuzzy sets as a mathematical way of representing impreciseness or vagueness in everyday life. Fuzzy set theory has numerous applications in various fields such as engineering, medicine, manufacturing, and others. In real-life situations, the processing times of jobs are not always exact due to incomplete knowledge or uncertain environment which implies the existence of various external sources and types of uncertainty. Fuzzy set theory can be used to handle uncertainty inherent in actual scheduling problems.

On single machine bi-criteria, Smith (Citation1956) was the first to consider minimization of weighted completion time subject to minimal value of Tmax. The other literature on single machine bi-criteria includes (Chand & Schneeberger, Citation1986; Gupta, Hariri, & Potts, Citation1999). Su, Chen, and Chen (Citation2012) dealt with scheduling on parallel machines to minimize the maximum tardiness. For parallel-machine scheduling problems concerning bi-criterion, Parkash (Citation1997) studied the bi-criteria scheduling problems on parallel machines and analyzed that the bi-criteria function is NP-hard. Sarin and Hariharan (Citation2000) proposed a heuristic for a two-machine case with objective as minimizing number of tardy jobs under the constraint of minimizing maximum tardiness. A review of scheduling on parallel machines was given by Mokotoff (Citation2001). Gupta, Ruiztorres, and Webster (Citation2003) considered the scheduling of jobs on identical parallel machines for minimizing maximum tardiness with optimal flow time. In his paper, the type of scheduling problems is presented as: Identical parallel machines; Uniform parallel machines and unrelated parallel machines. Optimizing dual performance measures on parallel machines in fuzzy environment is fairly an open area of research. A survey of literature has revealed little work reported on the bi-criteria scheduling problems on parallel machine in fuzzy environment. Jha, Indumati, Singh, and Gupta (Citation2011) formulated a fuzzy bi-criterion release time problem related to cost minimization and reliability maximization objectives under budgetary constraint. Sunita and Singh (Citation2010a, Citation2010b) studied the optimization of various characteristics on parallel machines in fuzzy environment. Gupta, Gupta, Sharma, and Aggarwal (Citation2012) developed an algorithm for the bi-objective problem which optimizes the number of tardy jobs without violating the value of Tmax with uncertain processing time. Sharma, Gupta, and Seema (Citation2013) developed an algorithm to schedule jobs on parallel identical machines so as to minimize the secondary criteria of weighted flow time without violating the primary criteria of maximum tardiness in fuzzy environment. Sharma, Gupta, and Seema (Citation2012) studied the bi-objective problem with total tardiness and number of tardy jobs as primary and secondary criteria, respectively, for any number of parallel machines in fuzzy environment.

The two approaches can be involved for bi-criteria problems: both the criteria’s are optimized simultaneously by using suitable weights for the criteria’s. Second, the criteria’s are optimized sequentially by first optimizing the primary criterion and then the secondary criterion subject to the value obtained for the primary criterion. In the present paper, we study the bi-objective scheduling on parallel machines in which the problem is divided into two steps: in the primary step, the weighted flow time of jobs is calculated and in the secondary step, the maximum tardiness is minimized with bi-objective function as Tmax/WFT. The concept of fuzziness in processing time of jobs is introduced as in real-life situations the processing time of jobs may vary due to the lack of complete knowledge and the decision-making process gets more complicated with increment in the number of constraints.

The rest of the paper is organized as follows: In Section 2, problem is formulated. Section 3 describes the basics of fuzzy set theory. Section 4 deals with various theorems for optimizing the bi-criteria problem. Section 5 describes the algorithm proposed to find the optimal sequence for bi-criteria problem Tmax/WFT. In Section 6, numerical illustrations are given to support the proposed algorithm. The paper is concluded in Section 7 followed by the references.

2. Problem formulation

The following notations have been used throughout the present paper.

i=

designate the ith job, i = 1, 2, 3, …, n

di=

due date of the ith job

Ci=

completion time of the ith job

wi=

weightage of the ith job

Ti=

tardiness of the ith job = max (Cidi, 0)

Tmax=

maximum tardiness

n=

total number of jobs to be scheduled

m=

total number of parallel machines

WFT(S)=

weighted flow time of jobs for sequence S

The following assumptions are taken into consideration for developing bi-criteria algorithm on parallel machines.

  1. The jobs are available at time zero.

  2. No pre-emption of jobs is allowed.

  3. The machines are identical in all respects.

  4. The machines are available all the time.

  5. The jobs are independent of each other.

  6. Each machine can process at most one job at a time.

Before formulating the bi-criteria problem, the mathematical formulation for the single criterion’s are represented first. These are:

2.1. Criterion: weighted flow time

The weighted flow time is the weighted version of flow time. The formulation to minimize the WFT is as follows:MinimizeWFT(S)=i=1nwiCi

Subject to the condition that all weights for jobs are non-negative,

Here, S is the given sequence and WFT(S) is the corresponding objective function of weighted flow time.

2.2. Criterion: maximum tardiness

Tardiness is given by Ti = max (Ci − di, 0), where Ci and di are the completion time and due date of the ith job, respectively. Maximum tardiness is given by Tmax = max(Ti).

The formulation is as follows:MinimizeZ=Tmax,Subjecttoconstraint:TmaxCi-dii

along with the non-negativity constraints.

The formulation of the bi-criteria problems is similar to that of single criteria problems but with some additional constraints requiring that the optimal value of the primary objective function is not violated. Here, the problem is divided into two steps: one is primary in which weighted flow time of jobs is minimized and in secondary step, the maximum tardiness of jobs is minimized under the objective function value of primary criterion.

3. Basic fuzzy set theory

A fuzzy set A defined on the universal set of real numbers R, is said to be a fuzzy number if its membership function has the following characteristics:

  1. μA:R[0,1] is continuous.

  2. μA=0 for all x ∈ (−∞, a1) ∪ (a3, ∞).

  3. μA is strictly increasing on [a1, a2] and strictly decreasing on [a2, a3].

  4. μA=1 for x = a2.

3.1. Triangular fuzzy number

Triangular fuzzy number (TFN) is a fuzzy number represented with three points as A=(a1,a2,a3), where a1 and a3 denote the lower and upper limits of support of a fuzzy set A. The membership value of the x denoted by μA(x),xR+, can be calculated according to the following formula (Figure ).μA(x)=0;xa1x-a1a2-a1;a1<x<a2a3-xa3-a2;a2<x<a30;xa3

Figure 1. Triangular fuzzy number.

Figure 1. Triangular fuzzy number.

3.2. Average high ranking

The system characteristics are described by membership function; it preserves the fuzziness of input information. However, the designer would prefer one crisp value for one of the system characteristics rather than fuzzy set. In order to overcome this problem, we defuzzify the fuzzy values of system characteristic by using the Yager’s (Citation1981) formula. If A=(a1,a2,a3), then its crisp value is defined as:(1) crisp(A~)=h(A~)=3a2+a3-a13(1)

3.3. Fuzzy arithmetic operations

If A~1=(mA1,αA1,βA1) and A~2=(mA2,αA2,βA2) be the two TFNs, then

  1. A~1+A~2=(mA1,αA1,βA1)+(mA2,αA2,βA2)=(mA1+mA2,αA1+αA2,βA1+βA1).

  2. A~1-A~2=(mA1,αA1,βA1)-(mA2,αA2,βA2)=(mA1-mA2,αA1-αA2,βA1-βA2), if the following condition is satisfied DPA1DPA2, where DP(A1)=βA1-mA12 and DP(A2)=βA2-mA22. Here, DP denotes difference point of a TFN.

  3. Otherwise; A~1-A~2=(mA1,αA1,βA1)-(mA2,αA2,βA2)=(mA1-βA2,αA1-αA2,βA1-mA2).

  4. kA~1=k(mA1,αA1,βA1)=(kmA1,kαA1,kβA1); if k > 0.

  5. kA~1=k(mA1,αA1,βA1)=(kβA1,kαA1,kmA1); if k < 0.

4. Theorems

The following theorems have been established to optimize the bi-criteria scheduling on parallel machines involving weighted flow time and maximum tardiness.

4.1. Theorem

A sequence S of jobs following weighted smallest processing time (WSPT) rule, in which the jobs are placed to the earliest available location on the machines in WSPT order, minimizes the weighted flow time.

Proof

Let if possible sequence S obtained by using the WSPT rule, i.e. by arranging the jobs in the decreasing order of their weights; break the ties (if any) arbitrary is not optimal. Let there exist a better sequence of jobs S′ (say) in which ith and jth adjacent jobs are interchanged.

Let Ci(S) and Cj(S) be the completion times of ith and jth job, respectively, in schedule S, respectively. Similarly, let Ci(S′) and Cj(S′) be the completion times of ith and jth job, respectively, in schedule S′.

For sequence S: We haveCi(S)=ta+1,Cj(S)=ta+2

For sequence S′: We haveCi(S)=ta+2,Cj(S)=ta+1

Next, the weighted flow time for the sequence S of jobs is(2) WFT(S)=wi×Ci(S)+wj×Cj(S)=wi×(ta+1)+wj×(ta+2)=wita+wi+wjta+2wj=(wi+wj)ta+wi+2wj(2)

Similarly, the weighted flow time for the sequence S′ jobs is(3) WFT(S)=wi×Ci(S)+wj×Cj(S)=wi×(ta+2)+wj×(ta+1)=wita+2wi+wjta+wj=(wi+wj)ta+2wi+wj(3)

Since, the ith and jth jobs are placed by WSPT rule. Therefore, we have wi ≥ wj.

Hence, from results (2) and (3), we have

WFT(S′) ≥ WFT(S), i.e. weighted flow time for the sequence S′ is more as compared to the sequence S.

Hence, the sequence S following the WSPT rule minimizes the weighted flow time.

4.2. Theorem

A sequence S of jobs following early due date (EDD) order is an optimal sequence with maximum tardiness (Tmax).

Proof

Let if possible sequence S following EDD is not an optimal sequence. Let there exist a better sequence of jobs S′ (say) i.e. a late job in S is early in S′. Let S′ be the schedule in which ith and jth adjacent jobs are interchanged, i.e.

S=1-2-3-4--i-j--nS=1-2-3-4--j-i--n

The ith and jth jobs are adjacent on parallel machines if they are placed on adjacent positions as shown in Figure . Further, an adjacent pair of jobs may be located on the same machine, say k, or on different machines, say k and l. The jobs appearing on the machines before an adjacent pair are designated by A and those appearing after the pair are designated by B with an appropriate subscript for the corresponding machine. Let ta be the completion time of the job before ith job on a machine.

Figure 2. Adjacent jobs on parallel machines.

Figure 2. Adjacent jobs on parallel machines.

In schedule S: Irrespective of the fact that the ith and jth adjacent jobs are on same machine or on different machines, we haveTi=maxta+1-di,0,Tj=maxta+2-dj,0

In schedule S′: We haveTi=maxta+2-di,0,Tj=maxta+1-dj,0

Since, di ≤ dj, we have TiTi and TjTj.

Therefore, maxTi,TjmaxTi,Tj.

This implies that the sequence S is better.

Hence, sequence S of jobs following EDD order is the optimal sequence for Tmax.

4.3. Theorem

If all the jobs have distinct weights, then no improvement can be made to the values of any secondary criterion of maximum tardiness.

Proof

Let S be the schedule obtained by arranging the jobs by WSPT order. Let Ci be the completion time of ith job and Cj be the completion time of jth job. Let Si and Sj be the corresponding secondary criterion values. When ith and jth jobs are switched, let S′ be the schedule so obtained.

If Ci < Cj then, Ci=Cj and Ci=Cj. Hence, wiCi+wjCj<wiCi+wjCj, which is a violation to primary criterion as all the other jobs contribute the same. Therefore, it implies that no improvement to the value of secondary criterion can be made with any possible exchange of jobs. The only possibility left is to rearrange the jobs having the same weights so that the value of a secondary criterion can be improved.

4.4. Theorem

If the problem of single criteria, maximum tardiness is NP-hard, the scheduling problem on parallel machines optimizing the bi-objective function maximum tardiness/WFT will also be NP-hard.

Proof

We shall prove the result by the method of contradiction:

Let if possible the bi-objective function maximum tardiness/WFT is not NP-hard. Therefore, there must exist a polynomial algorithm which can solve the problem of optimizing the bi-objective function maximum tardiness/WFT on parallel processing machines.

This implies that single criteria of maximum tardiness can be optimized in polynomial time, i.e. maximum tardiness is not NP-hard. This is a contradiction, as maximum tardiness is NP-hard.

Hence, the scheduling problem optimizing the bi-objective function maximum tardiness/WFT on parallel processing machines will also be NP-hard.

5. Algorithm

The following algorithm is proposed to find the optimal sequence for bi-criteria problem Tmax/WFT:

Step 1: Find average high ranking (AHR) of the fuzzy processing time of various jobs on different machines.

Step 2: Arrange all the jobs in weighted shortest processing time (WSPT) order.

Step 3: If all the jobs have distinct weights then go to Step 5.

Step 4: If all the jobs do not have distinct weights then arrange all the jobs with same weights in EDD order. Break the tie if any arbitrarily. Note the improvement in Tmax subject to the primary set of constraints.

Step 5: Assign the jobs one by one to the earliest available machine. This optimizes the objective Tmax/WFT.

6. Numerical illustration

The following numerical illustrations are given to test the efficiency of algorithm proposed to optimize the bi-criteria Tmax/WFT on parallel machines in fuzzy environment.

Consider an example of jobs with processing time in fuzzy environment, due date and with jobs weightage (being distinct) on three parallel machines in a flow shop. The objective is to obtain a sequence of the jobs processing optimizing the bi-criteria taken as Tmax/WFT (Table ).

Table 1. Jobs with fuzzy processing time

Solution: The AHR of processing time of given jobs as per the Step 1 using Equation 1 is as follows: (Table )

Table 2. Jobs with AHR of processing time

On arranging the jobs as per Step 2 in WSPT order on the parallel available machines M1, M2, and M3, the order of jobs becomes 4–2–1 –3–5. The in–out i.e. the flow of jobs on machines for WSPT order is as in Table .

Table 3. In–out of jobs following WSPT order

Therefore, Tmax = 25/3 units and weighted flow time, WFT=i=1nwiCi=5×233+4×533+3×323+2×553+1×703=201units.

As per the Step 3, all the jobs have distinct weights; therefore, assigning the jobs to the available machines optimizes secondary criteria of Tmax. Hence, the optimal sequence which optimizes the objective Tmax/WFT is 4–2–1–3–5 with WFT = 201 units and minimum Tmax = 25/3.

Consider another example of jobs with processing time in fuzzy environment, due date and with jobs weightage (being not distinct) on three parallel machines in a flow shop is given. The objective is to obtain sequence of the jobs processing optimizing the bi-criteria taken as Tmax/WFT (Table ).

Table 4. Jobs with fuzzy processing time

Solution: The AHR of processing time of given jobs as per Step 1 using Equation 1 is as follows: (Table )

Table 5. Jobs with AHR of processing time

On arranging the jobs as per Step 2 in WSPT order on the parallel available machines M1, M2, and M3, the order of jobs becomes1–6–2–3–4–7–5. The in–out flow of jobs on machines is as in Table .

Table 6. In–out of jobs in WSPT order

Therefore, maximum tardiness, Tmax = 11 units and WFT=i=1nwiCi=5×263+4×383+3×233+3×583+2×553+2×673+1×783=282.3units.

Since the weights corresponding to jobs are not distinct, therefore, as per the Step 3 of algorithm we optimize Tmax with given WFT. Following this we arrange the jobs 2 and 3 with same weightage as 3 in EDD order. Also, arrange the jobs 4 and 7 with same weightage as 2 in EDD order. Since the job 3 has EDD than the job 2, therefore, the job 3 will come before the job 2 in EDD order. The jobs 4 and 7 are already in EDD order. Following this we have the in–out of jobs as in Table .

Table 7. In–out of jobs following EDD order

Therefore, maximum tardiness, Tmax = 34/3 units and WFT = i=1nwiC=5×263+4×383+3×353+3×493+2×643+2×673+1×723=288units which violates the primary criterion of WFT having value 288 units. Therefore, reject the sequence1–6–3–2–4–7–5. Hence, the optimal sequence of jobs processing is 1–6–2–3–4–7–5 with minimum WFT as 282.3 units and minimum Tmax as 11 units.

7. Conclusion

In this paper, an algorithm to optimize the bi-criteria scheduling problem involving maximum tardiness and weighted flow time on parallel machines in fuzzy environment is considered. The concept of fuzzy processing time is introduced in processing of jobs, as in many real-life situations the processing time of jobs may vary due to the lack of complete knowledge, uncertainty, and vagueness. The study may further be extended by introducing the concepts of non-availability constraints, machines processing the jobs with different speeds, and also by considering trapezoidal fuzzy membership functions to represent the fuzziness in processing time.

Additional information

Funding

Funding. The authors received no direct funding for this research.

Notes on contributors

Sameer Sharma

The group of authors involved in this critically researched paper consists in Deepak Gupta, Sameer Sharma, and Kewal Krishan Nailwal. The team under the competent supervision of Deepak Gupta works in the area of scheduling problems. Some research work on optimizing scheduling problems in industry with various parameters along with the linkage models has been contributed. This paper particularly optimizes a fuzzy bi-criteria problem involving weighted flow time and tardiness. This paper will go in a long way to solve the problems faced in production management where a quality service means product is to be processed with minimum weighted flow time and meeting customer’s required due date or satisfaction of the customer. Besides this, a fuzzy environment of the problem brings the problem closer to real-life application area as fuzzy theory tackles well with vagueness and uncertainty in processing times of jobs.

References

  • Baker, K. R. (1974). Introduction to sequencing and scheduling. New York, NY: Wiley.
  • Chand, S., & Schneeberger, H. (1986). A note on the single-machine scheduling problem with minimum weighted completion time and maximum allowable tardiness. Naval Research Logistics Quarterly, 33, 551–557.10.1002/(ISSN)1931-9193
  • Gupta, D., Sharma, S., & Aggarwal, S. (2012). Bi-objective scheduling on parallel machines with uncertain processing time. Advances in Applied Science Research, 3, 1020–1026.
  • Gupta, J. N. D., Hariri, A. M. A., & Potts, C. N. (1999). Single-machine scheduling to minimize maximum tardiness with minimum number of tardy jobs. Annals of Operations Research, 92, 107–123.10.1023/A:1018974428912
  • Gupta, J. N. D., Ruiz-torres, A. J., & Webster, S. (2003). Minimizing maximum tardiness and number of tardy jobs on parallel machines subject to minimum flow-time. Journal of the Operational Research Society, 54, 1263–1274.10.1057/palgrave.jors.2601638
  • Jha, P. C., Indumati, N. A., Singh, O., & Gupta, D. (2011). Bi-criterion release time problem for a discrete SRGM under fuzzy environment. International Journal of Mathematics in Operational Research, 3, 680–696.10.1504/IJMOR.2011.043016
  • Mokotoff, E. (2001). Parallel machine scheduling problems: A survey. Asia-Pacific Journal of Operational Research, 18, 193–242.
  • Parkash, D. (1997). Bi-criteria scheduling problems on parallel machines ( PhD thesis). University of Birekshurg, Virginia, VA.
  • Sarin, S. C., & Hariharan, R. (2000). A two machine bicriteria scheduling problem. International Journal of Production Economics, 65, 125–139.10.1016/S0925-5273(99)00050-X
  • Sharma, S., Gupta, D., & Seema., S. (2012). Bi-objective scheduling on parallel machines in fuzzy environment. Proceedings of the Second International Conference on Soft Computing for Problem Solving (SocProS 2012), 1, 365–372.
  • Sharma, S., Gupta, D., & Seema, S. (2013). Bicriteria scheduling on parallel machines: Total tardiness and weighted flowtime in fuzzy environment. International Journal of Mathematics in Operational Research, 5, 492–507.10.1504/IJMOR.2013.054733
  • Smith, W. E. (1956). Various optimizers for single stage production. Naval Research Logistics Quarterly, 3, 59–66.10.1002/(ISSN)1931-9193
  • Su, L.-H., Chen, P.-S., & Chen, S.-Y. (2012). Scheduling on parallel machines to minimize maximum tardiness for the customer order problem. International Journal of System Sciences, 44, 926–936. doi:10.1080/00207721.2011.649366
  • Sunita, & Singh, T. P. (2010a). Bi-objective in fuzzy scheduling on parallel machines. Arya Bhatta Journal of Mathematics and Informatics, 2, 149–152.
  • Sunita, & Singh, T. P. (2010b). Single criteria scheduling with fuzzy processing time on parallel machines. In Proceedings of International Conference in IISN 2010, ISTK (pp. 444–446). Yamuna Nagar, Haryana.
  • Yager, R. R. (1981). A procedure for ordering fuzzy subsets of the unit interval. Information Sciences, 24, 143–161.10.1016/0020-0255(81)90017-7
  • Zadeh, L. A. (1965). Fuzzy sets. Information and Control, 8, 338–353.10.1016/S0019-9958(65)90241-X