Fuzzy Nonlinear Programming with Applications in Decision Making
View this Special IssueResearch Article  Open Access
HorngRen Tsai, Toly Chen, "A Fuzzy Nonlinear Programming Approach for Optimizing the Performance of a FourObjective Fluctuation Smoothing Rule in a Wafer Fabrication Factory", Journal of Applied Mathematics, vol. 2013, Article ID 720607, 15 pages, 2013. https://doi.org/10.1155/2013/720607
A Fuzzy Nonlinear Programming Approach for Optimizing the Performance of a FourObjective Fluctuation Smoothing Rule in a Wafer Fabrication Factory
Abstract
In theory, a scheduling problem can be formulated as a mathematical programming problem. In practice, dispatching rules are considered to be a more practical method of scheduling. However, the combination of mathematical programming and fuzzy dispatching rule has rarely been discussed in the literature. In this study, a fuzzy nonlinear programming (FNLP) approach is proposed for optimizing the scheduling performance of a fourfactor fluctuation smoothing rule in a wafer fabrication factory. The proposed methodology considers the uncertainty in the remaining cycle time of a job and optimizes a fuzzy fourfactor fluctuationsmoothing rule to sequence the jobs in front of each machine. The fuzzy fourfactor fluctuationsmoothing rule has five adjustable parameters, the optimization of which results in an FNLP problem. The FNLP problem can be converted into an equivalent nonlinear programming (NLP) problem to be solved. The performance of the proposed methodology has been evaluated with a series of production simulation experiments; these experiments provide sufficient evidence to support the advantages of the proposed method over some existing scheduling methods.
1. Introduction
In complex manufacturing systems, such as wafer fabrication factories, job scheduling is subject to many sources of uncertainty or randomness [1]. Such uncertainty or randomness is partly due to manual operations, including the loading and unloading of jobs, the setup or repair of machines, and visual inspections. The other two sources of uncertainty, the unexpected releases of emergency orders and machine breakdowns, are beyond the control of a wafer fabrication factory. The literature provides probabilistic (stochastic) and fuzzy methods that can consider uncertainty or randomness. However, it is difficult to identify the probability distribution of each parameter, which means that a probabilistic (stochastic) method is not easy to use. In addition, fuzzy methods are advantageous because subjective factors can be considered, such as human interpretations of the scheduling performance and the tradeoffs of different scheduling objectives. For example, a job 3 months late and a job 3 days late are both late. However, the first job is difficult to accept, while the second is still acceptable. In other words, there are different degrees of acceptance, even if both jobs have been delayed. Second, to one scheduler, one objective may be much more important than the other, but to another scheduler, the two objectives may be equally important. The concepts of “acceptability” and “relative importance” can both be suitably modeled by a fuzzy method. For example, Murata et al. [2] used trapezoidal fuzzy numbers (TrFNs) to represent the satisfaction levels of due dates. In the literature, due dates, processing times, and precedence relations have been fuzzified [2–4].
Many existing fuzzy scheduling methods take the form of fuzzy inference rules, such as “if the release time is early and the number of operations is large, then the job priority is high” [5, 6]. A fuzzy scheduling system usually uses a number of fuzzy inference rules and can be divided into two types: Mamdani [7] and TakagiSugenoKang (TSK). For example, Xiong et al. [5] scheduled a flexible manufacturing system (FMS) using two fuzzy dispatching rules of the TSK type. Murata et al. [2] used six fuzzy rules to move jobs between different priority classes. Lee et al. [8] established two fuzzy inference rules to select a combination of some existing dispatching rules for scheduling a flexible manufacturing system. Tan and Tang [9] applied Taguchi’s design of experiment (DOE) techniques to improve the design of some fuzzy dispatching rules for a test facility. In Benincasa et al. [6], up to 27 fuzzy inference rules (each with three inputs and one output) were established to schedule automated guided vehicles. Dong and Liu [3] built an adaptive neurofuzzy inference system (ANFIS) to schedule a job shop. For any two jobs, the inputs to the ANFIS were the differences between the two jobs, and the output from the ANFIS determined the sequence of the two jobs. If the output was greater than zero, then the first job was to be processed before the second job. Murata et al. [2] considered a job shop with 10 machines and jobs with different priorities. A fuzzy linear programming (FLP) problem was solved to optimize the total reward. However, the fuzziness came from the satisfaction level of the due date rather than from the parameters. A good review of the literature on fuzzy scheduling methods can be found in Dubois et al. [10].
On the other hand, a scheduling problem can be formulated as a mathematical programming problem. For example, Biggs and Laughton [11] optimized a nonlinear programming (NLP) model for electric power scheduling. A recursive quadratic programming approach was proposed to solve the NLP problem. In Pedro [12], the problem was formulated as a mixed integer programming (MIP) model. The optimal solution of the mathematical programming problem gives the optimal schedule for the manufacturing system. However, sometimes the mathematical programming problem is not easy to solve, and some soft computing methods can be applied to search for the optimal solution of the mathematical programming problem [13–15]. For example, Chiang and Fu [16] minimized the number of tardy jobs for a job shop in which machines have sequencedependent setup times. An NLP problem with a linear objective function and some quadratic constraints was solved by the application of a genetic algorithm (GA). Ishibuchi and Murata [13] applied a similar approach for multiobjective flow shop scheduling. In the NLP model of Connors et al. [17], the inventory level was estimated with a nonlinear equation and then the holding costs were minimized. Chen and Yao applied a deterministic fluid network [18] in order to find an optimal solution. Both Murata et al. [2] and Ishibuchi et al. [19] built FLP models to solve scheduling problems and maximized the satisfaction levels of the due dates. Murata et al. applied GA to solve the FLP problem. For the same purpose, Ishibuchi et al. [19] applied a hybrid GA with neighborhood search. Xia and Wu [15] combined particle swarm optimization (PSO) and simulated annealing (SA) for flexible job shop scheduling.
A summary of existing fuzzy scheduling methods is shown in Table 1. The existing approaches have the following problems.(1)Dubois et al. [10] distinguished two categories of fuzzy scheduling methods: methods that represent preference profiles and methods that model uncertainty distributions. However, in most studies that use fuzzy methods, the fuzziness comes from the fuzzification of the scheduling objective (which belongs to the first category) or from the fuzzy rules that are subjectively chosen by the scheduler rather than from fuzzy parameters (which belong to the second category). In other words, processing time, due date, and precedence relations are all free from uncertainty in these studies.(2)Although an NLP problem is not easy to solve, the soft computing method applied to solve the NLP problem may also pose a considerable challenge.(3)The combination of NLP with fuzziness results in a FNLP approach that considers the uncertainty of parameters and does not need to make simplifying assumptions, and therefore has the potential to solve realistic scheduling problems effectively. However, very few studies applied FNLP methods. Most of them are used for scheduling small manufacturing systems [20].

To tackle these problems, an FNLP approach is proposed in this study to optimize the performance of a job dispatching rule in a wafer fabrication factory. In other words, this study is not going to optimize the performance of a schedule for a wafer fabrication factory, which is known to be an NPhard problem but to optimize the performance of a dispatching rule in a wafer fabrication factory. The use of some special types of fuzzy numbers establishes the corresponding subcategories for these FNLP models, such as the type2 FNLP (with type2 fuzzy numbers), the intervalvalued FNLP (with interval fuzzy numbers), the intuitionistic FNLP (with intuitionistic fuzzy numbers), and the intervalvalued intuitionistic FNLP [21], which has been of interest to researchers in recent years. Fares and Kaminska [22] solved two multipleobjective FNLP problems to find the optimal sets of circuit parameter values for a bipolar emitter follower circuit and an unbuffered twostage complementary metaloxide semiconductor (CMOS) opamp. Chen and Wang [23] defined the yield competitiveness of a semiconductor product, which is uncertain and can be enhanced by allocating more capacity to the product; this functions in a nonlinear way. They therefore constructed an FNLP model to optimize the effects of capacity reallocation on the yield competitiveness of a semiconductor product, which was then converted to an NLP problem to be solved.
In the proposed methodology, the fuzziness comes from the uncertainty of the remaining cycle time, that is, the time still needed to complete a job; this time is highly uncertain [24]. In the proposed methodology, the remaining cycle time of a job is estimated with a triangular fuzzy number (TFN). There are various types of fuzzy numbers with different shapes. Among them, a TFN is easily implemented and has been used for numerous applications (e.g., [23–25]). Subsequently, the remaining cycle time estimate is fed into a fourfactor fluctuation smoothing rule [25] to sequence the jobs in front of a machine. The fourfactor fluctuation smoothing rule is fuzzified in this way, and the slack of a job is also expressed by a TFN. The fuzzy fluctuationsmoothing dispatching rule has five adjustable parameters, the optimization of which constitutes an FNLP problem. To convert the FNLP problem into a more tractable form, cut operations are also applied.
The unique features of the proposed methodology include the following.(1)Combining fuzziness and NLP: scheduling decisions represented in terms of fuzzy sets are flexible in their implementations. NLP relaxes the strict assumptions and constraints of linear programming (LP) and is highly practicable.(2)Considering the uncertainty in the remaining cycle time: dispatching rules that consider dynamic information, such as the remaining cycle time, are more effective for highly complex manufacturing systems [26]. To this end, an effective fuzzy back propagation network (FBPN) approach is applied. (3)Establishing a fuzzy dispatching rule directly from the existing rules: a fuzzy dispatching rule is deduced by fuzzifying the fourobjective fluctuationsmoothing dispatching rule [25] and diversifying the slack. This rule accepts the fuzzy remaining cycle time as an input and uses a fuzzy value to represent the slack of each job.(4)Diversifying the slack by solving an FNLP model: the emergence of ties may lead to incorrect scheduling results. In the proposed methodology, to reduce the number of ties, the slacks of jobs are diversified by maximizing the standard deviation [27], which leads to a FNLP problem. The FNLP problem is not easy to solve; therefore, this study applies cut operations [28].(5)Optimizing four objectives simultaneously: the proposed fuzzy rule fuses four dispatching rules in a nonlinear way. In contrast, most existing methods optimize the weighted sum of multiple objectives (e.g., [13, 15]).
The rest of this paper is organized as follows. Section 2 is divided into four parts: cut operations, effective FBPN, fuzzified dispatching, and FNLP. First, the concepts of cuts and cut operations are introduced. The next part explains the effective FBPN approach that estimates the remaining cycle time of a job with a fuzzy number. The next part shows that the fourobjective fluctuation smoothing rule is fuzzified so that it can accept the remaining cycle time estimate as an input. The final part of Section 2 explains the role of the FNLP. To obtain the best values of the parameters in the fuzzy fourobjective fluctuation smoothing rule, and to diversify the slack, an FNLP model is built. To solve the FNLP problem, cut operations are applied. Section 3 details how a series of production simulation experiments are carried out to assess the advantages and disadvantages of the proposed methodology. Finally, the conclusions of this study are made in Section 4.
2. Methodology
The flow chart of Figure 1 illustrates the steps of the proposed methodology.
Subsequently, the variables and parameters that will be used in the proposed methodology are defined as follows.(1): the cycle time of job .(2): the estimated cycle time of job ; .(3): the due date of job .(4): the release time of job .(5): the remaining cycle time of job from step .(6): the estimated remaining cycle time of job from step ; .(7): the remaining processing time of job from step .(8): the step cycle time of job until step .(9) or : the slack of job at step .(10): the current time.(11): the total processing time of job .(12): mean release rate.(13): the inputs to the threelayer BPN of job .(14): the output from hiddenlayer node .(15): the connection weight between hiddenlayer node and the output node.(16): the connection weight between input node and hiddenlayer node , .(17): the threshold on hiddenlayer node .(18): the threshold on the output node; .
All fuzzy parameters in the proposed methodology are given in TFNs
2.1. Cuts and Cut Operations
The cut operations are applied to solve the FNLP problem. For this reason, the concepts of cuts and cut operations are introduced as follows.
Definition 1 ( cuts). is a fuzzy number. The cut of is an interval number given by
Definition 2 (arithmetic of fuzzy numbers based on cut operations). Given two fuzzy numbers and , and their cuts and , the arithmetic operations of and based on their cuts are as follows: where , , , and denote fuzzy addition, subtraction, multiplication, and division, respectively. Equation (5) is equivalent to if . The min and max function can be replaced by
Theorem 3 (average of fuzzy numbers based on cuts). Given fuzzy numbers , , the average of these fuzzy numbers can be derived as
Proof. Theorem 3 can be directly derived from (2) and (5).
2.2. The Effective FBPN Approach for Estimating the Remaining Cycle Time
Before any job is scheduled, the remaining cycle time of each job needs to be estimated. In this work, the effective FBPN approach is applied and the remaining cycle time is estimated with a fuzzy value.
In the effective FBPN approach, jobs are classified into categories using fuzzy cmeans (FCM). First, in order to facilitate the subsequent calculations and problem solving, all raw data are normalized [29]. Then, we place the (normalized) attributes of job in vector .
FCM classifies jobs by minimizing the following objective function: where is the required number of categories; is the number of jobs; indicates that job belongs to category ; measures the distance from job to the centroid of category ; is a parameter to adjust the fuzziness and is usually set to 2. The procedure of FCM is as follows(1)Produce a preliminary clustering result: the performance of FCM is sensitive to the initial conditions.(2)(Iterations) calculate the centroid of each category as where is the centroid of category . is the membership function that indicates job belongs to category after the th iteration.(3) Remeasure the distance from each job to the centroid of each category, and then recalculate the corresponding membership.(4) Stop if the following condition is met. Otherwise, return to step (2): where is a real number representing the threshold for the convergence of membership.
Finally, the separate distance test ( test) proposed by Xie and Beni [30] can be applied to determine the optimal number of categories : The value that minimizes determines the optimal number of categories.
After clustering, a threelayer FBPN is used to estimate the cycle times of jobs for each category. The configuration of the threelayer FBPN is as follows. First, the inputs are the parameters associated with the th job. Subsequently, there is only a single hidden layer; the hidden layer has twice as many neurons as the input layer. In addition, Chen and Wang [31] and Chen and Lin [32] have described how an NLP model can be constructed to adjust the connection weights and thresholds in an FBPN; this problem is not easy to solve. In the proposed methodology, only the threshold on the output node () will be adjusted. This way is much simpler and can also achieve good results. In other words, only is fuzzy, while the other parameters are crisp. In this way, the fuzzy remaining cycle time estimate is generated with minimal effort. This makes the FBPN approach an effective one. The output from the threelayer FBPN is the (normalized) estimated remaining cycle time (()) of the training examples, where is the normalization function.
The procedure for determining the parameter values is now described. First, to determine the value of each parameter and , the FBPN is treated as a crisp network. Some algorithms are applicable for this purpose, such as gradient descent algorithms, conjugate gradient algorithms, the LevenbergMarquardt algorithm, and others. In this study, the LevenbergMarquardt algorithm is applied. The LevenbergMarquardt algorithm was designed for training with secondorder speed without having to compute the Hessian matrix. It uses approximation and updates the network parameters in a Newtonlike way [33].
Subsequently, is to be determined, so that the actual value will be less than the upper bound of the network output. Assume that the adjustment made to the threshold on the output node is denoted as . The optimal value of should be set as follows: In a similar way, can be determined so that each actual value will be greater than the appropriate lower bound. The optimal value of can be obtained as: This FBPN approach can generate a very precise interval of the remaining cycle time for each job, thereby reducing the risk of misscheduling. An instance has been analyzed in Figure 2 to evaluate the performance of this method. To provide a comparison, a statistical analysis method is also applied to this instance, in which the relationship between the remaining cycle time and job attributes is fitted with a multiple regression equation. The results are shown in Figure 3. Compared with the effective FBPN approach, the statistical analysis method is not only inaccurate but also not precise enough. The remaining cycle time estimated by the statistical analysis method is therefore prone to errors, which may result in incorrect scheduling.
2.3. The Fuzzy FourObjective Dispatching Rule
Lu et al. [26] proposed two fluctuation smoothing rules—the fluctuation smoothing policy for mean cycle time (FSMCT) and the fluctuation smoothing policy for variation of cycle time (FSVCT). FSMCT is aimed at minimizing the mean cycle time, while FSVCT is aimed at minimizing the variance of cycle time:
(FSMCT)
(FSVCT) Jobs with the smallest slack values are given the highest priorities.
If the remaining cycle time is estimated with a TFN, then we have two fuzzy fluctuation smoothing rules as
(fuzzy FSMCT)
(fuzzy FSVCT)
To determine the sequence of jobs, the fuzzy slacks must be compared. To this end, various methods have been proposed in the literature, such as a method based on the probability measure [34], a coefficient of variance (CV) index [35], a method that uses the area between the centroid point and the original point [36], and a method based on the fuzzy mean and standard deviation [37]. For a comparison of these methods, refer to Zhu and Xu [37]. In this study, the method based on the fuzzy mean and standard deviation is applied because it is relatively simple and can yield reasonable comparison results. To put this in context, the following theorem is introduced.
Theorem 4. The fuzzy mean and standard deviation of a triangular fuzzy number = (, , ) can be derived as
Proof. Refer to Zhu and Xu [37]. It is, in fact, the centerofgravity (COG) method.
The following definition details a method based on the fuzzy mean and standard deviation.
Definition 5. For any two fuzzy numbers and , the sequence of and can be determined according to their fuzzy means and standard deviations as follows.(1) if and only if .(2) if and only if .(3)If , then (i) if and only if .(ii) if and only if .(iii) if and only if .
Consider the example in Table 2. The sequencing results by the two fuzzified rules are Fuzzy FSMCT: 7 11 16 6 13 4 19 1 3 14 2 9 10 8 12 18 15 17 5. Fuzzy FSVCT: 7 13 19 16 11 1 6 4 9 15 17 8 3 12 2 14 18 5 10.

Chen [25] combined four traditional dispatching rules—EDD, critical ratio (CR), the fluctuation smoothing policy for mean cycle time (FSMCT)—and the fluctuation smoothing policy for variation of cycle time (FSVCT), and proposed the fourobjective dispatching rule. In the fourobjective dispatching rule, the slack of job at processing step is defined as where , and and are positive real numbers that satisfy the following constraints: Jobs with the smallest slack values will be given the highest priorities. There are many possible models that can form the combinations of , and . For example, The values of and are within .
If the remaining cycle time is estimated with a triangular fuzzy number, then (22) becomes Job is processed before job if . In Wang et al. [27], to diversify the slack, the standard deviation of the slack was maximized: When the job slack is a fuzzy value, Maximizing is equivalent to maximizing . Finally, the following FNLP problem is to be solved: The proposed FNLP problem is intractable and may need to be converted into an equivalent NLP problem to be solved. First, (27) can be decomposed to The cut of is Subsequently, the objective function is equal to The cut of is Applying (8) to (32) gives where Equation (33) is equivalent to Similarly, (34) can be replaced by In order to facilitate the solving of the problem, the usual practice is to defuzzify the fuzzy objective function , using the centerofgravity defuzzification [28]: where is the defuzzification function. Finally, the following NLP model is optimized instead of the original FNLP problem:
Consider the example in Table 3. Both the proposed methodology and Wang et al.’s method are applied to this example. In the proposed methodology, the optimal objective function value is 187272 when the parameters , and are equal to 0.980, 0.560, , and , respectively. The slacks of the jobs are shown in Figure 4. Please note that in this figure the axis is converted to logarithmic values for clarity. In contrast to this, Wang et al.’s method cannot consider the uncertainty in the remaining cycle time, and therefore only the center value of the remaining cycle time is considered. The slacks obtained by using Wang et al.’s method are shown in Figure 5, in which . The optimal values of the parameters , and are equal to 0.938, 0.587, , and , respectively. Obviously, one has the following.(1)After considering the uncertainty of the remaining cycle time, the best values of the five parameters changed, and the slacks of jobs became different. This might result in different sequencing results. (2)There were 14 ties in Wang et al.’s method. Conversely, after considering the uncertainty in the remaining cycle time, the proposed methodology successfully diversified the slacks of the jobs and reduced the number of ties to 11. In this regard, the advantage of the proposed methodology over Wang et al.’s approach is 21%.(3)In the method of Wang et al., if there is a long tail in the remaining cycle time of a job on the righthand side, then the slack of the job will be underestimated. Conversely, the slack will be overestimated if there is a long tail on the lefthand side.

3. Simulation Experiment
A real wafer fabrication factory mainly used for the production of dynamic random access memory (DRAM) was simulated. The wafer fabrication factory is located in Taiwan’s Taichung Science Park and has a monthly capacity of about 25,000 wafers. However, the following assumptions were made to generate data that are less noisy than realworld data.(1)The distributions of the times between machine breakdowns are exponential.(2)The distribution of the time required to repair a machine is uniform.(3)The percentages of jobs with different priorities released into the wafer fabrication factory are controlled.(4)A job has equal chances to be processed on each alternative machine or head that is available at a step.(5)A job cannot proceed to the next step until the processing of every wafer in the job has been finished.(6)No preemption is allowed.
In the simulated wafer fabrication factory, there are more than 10 types of memory products and more than 500 workstations for performing singlewafer or batch operations using 58 nm~110 nm technologies. Jobs released into the fabrication factory are assigned three types of priorities, that is, “normal,” “hot,” and “super hot.” Usually, a job will only be “super hot” if it is part of an emergency order; “super hot” jobs will be processed first. The large scale and the reentrant process flows of this wafer fabrication factory exacerbate the difficulties of job dispatching. Currently, the longest average cycle time exceeds three months with a variation of more than 300 hours. The managers of this wafer fabrication factory are therefore seeking better dispatching rules to replace FIFO and EDD, in order to shorten the average cycle times and ensure ontime delivery to customers.
One hundred replications of the simulation were successively run. The simulation horizon of each replication was twentyfour months. The warmup period was the first four months. The time required for each simulation replication was about 45 minutes using a PC with Intel Dual E2200 2.2 GHz CPUs and 1.99 G RAM.
To make comparisons with some existing approaches, eight methods were tested. FIFO, EDD, shortest remaining processing time (SRPT), CR, FSVCT, FSMCT, the nonlinear fluctuation smoothing rule (NFS), and the fourobjective slackdiversifying rule (4oSDR) [25] were applied to schedule the simulated wafer fabrication factory. The data of 1000 jobs were collected and separated by product types and priorities.
For FIFO, jobs were sequenced on each machine first by their priorities, then by their arrival times at the machine. For EDD, jobs were also sequenced first by their priorities, then by their due dates. The performance of EDD depends on how jobs’ due dates are determined. In the experiment, the due date of each job was determined as follows: where indicates the cycle time multiplier.
FSVCT and FSMCT consisted of two stages. First, jobs were scheduled based on FIFO, in which the remaining cycle times of all jobs were recorded and averaged at each step. Then, FSVCT/FSMCT policy was applied to schedule the jobs based on the average remaining cycle times obtained earlier. In other words, jobs were sequenced on each machine first by their priorities, and then by their slack values, which were determined by (17) and (18). With SRPT, the remaining processing time of each job was calculated. Then, jobs were sequenced first by their priorities, then by their remaining processing times. With CR, jobs were sequenced first by their priorities, then by their critical ratios. NFS is a nonlinear fusion of FSMCT and FSVCT. In the simulation experiment, a weight of 0.8 was given to FSMCT. With 4oSDR, the remaining cycle time of a job was estimated using the fuzzy means and back propagation network (FCMBPN) approach [31]; it was a crisp value. The five adjustable parameters were set to after initial scenarios had been examined.
In the proposed methodology, the remaining cycle time of a job was estimated using the effective FBPN approach; it was a fuzzy value. After the fuzzy remaining cycle time estimate had been fed into the fuzzy fourobjective fluctuation smoothing rule, an FNLP problem was solved to determine the values of the five parameters in the rule, so as to optimize the scheduling performance.
The average cycle time, cycle time standard deviation, the number of tardy jobs, and the maximum lateness of all cases were calculated to assess the scheduling performance. The results are summarized in Tables 4, 5, 6, and 7.




According to the experimental results, the following points can be made.(1)In various respects, the proposed methodology was obviously superior to the existing dispatching rules. For example, the fuzzy fourobjective fluctuation smoothing rule was the best at reducing the average cycle time for all cases. Its advantage over the current rule FIFO was 26% on average. The average cycle time is one of the most important scheduling goals of a wafer fabrication factory; such experimental results are very valuable. Of the traditional scheduling rules, SRPT performed well for reducing the average cycle times but posed the risk of high cycle time variation.(2)Ontime delivery is another important scheduling objective. The maximum lateness is often used to assess this. The proposed methodology can effectively reduce the maximum lateness and so can enhance ontime delivery. This tends to improve the firm’s customer relationships. The performance levels of the two traditional methods in this field, EDD and CR, were not as good as expected.(3)Through reduced cycle time variability, we can more accurately estimate the cycle time and promise our customer a more reliable due date. The fuzzy fourobjective fluctuation smoothing rule has very good performance in this respect, with an average advantage of 28% over the existing scheduling rules.(4)The number of tardy jobs is another indicator by which one can assess ontime delivery. The proposed methodology outperformed the existing methods in most cases.
In order to confirm the advantages of the proposed methodology over the existing methods, a Wilcoxon signedrank test [38] was used to test the following hypotheses.: When shortening the average cycle time, the scheduling performance of the proposed methodology is the same as that of the existing approach being compared. : When shortening the average cycle time, the scheduling performance of the proposed methodology is better than that of the existing approach being compared. : When reducing the maximum lateness, the scheduling performance of the proposed methodology is the same as that of the existing approach being compared. : When reducing the maximum lateness, the scheduling performance of the proposed methodology is better than that of the existing approach being compared. : When reducing the cycle time standard deviation, the scheduling performance of the proposed methodology is the same as that of the existing approach being compared. : When reducing the cycle time standard deviation, the scheduling performance of the proposed methodology is better than that of the existing approach being compared. : When reducing the number of tardy jobs, the scheduling performance of the proposed methodology is the same as that of the existing approach being compared. : When reducing the number of tardy jobs, the scheduling performance of the proposed methodology is better than that of the existing approach being compared.
The results are summarized in Table 8. The null hypothesis was rejected at , which showed that the fuzzy slackdiversifying fluctuationsmoothing rule was superior to seven existing approaches at reducing the average cycle time. With regard to maximum lateness, the advantage of the fuzzy slackdiversifying fluctuationsmoothing rule over FIFO, SRPT, and FSVCT was significant. Similar results could be observed with cycle time standard deviation. However, the advantage of the fuzzy slackdiversifying fluctuationsmoothing rule was not statistically significant for the number of tardy jobs.
 
P < 0.05. ** P < 0.025. *** P < 0.01. 
4. Conclusions and Directions for Future Research
Multiobjective scheduling is an important task in a wafer fabrication factory. It is also a difficult task owing to the uncertainty and complexity of the wafer fabrication system. This study demonstrates that an FNLP approach can consider such uncertainties and optimize the performance of multiobjective scheduling in a wafer fabrication factory.
The proposed methodology starts from the estimation of the remaining cycle time for each job. To this end, an effective FBPN approach has been proposed. The estimated remaining cycle time from the FBPN is a fuzzy number. After the fuzzy remaining cycle time estimate is fed to the fourfactor fluctuation smoothing rule, the rule is fuzzified, and the slack of each job is expressed by a fuzzy number. To reduce the number of ties, the slacks of jobs need to be diversified, which results in an FNLP problem. Since the FNLP problem is not easy to solve, some cut operations are applied to convert it into an equivalent NLP problem.
After a simulation study, the following points are concluded(1)The simulation experiment results showed that the proposed method indeed enhanced the scheduling performance of the wafer fabrication factory in four respects—average cycle time, maximum lateness, cycle time standard deviation, and number of tardy jobs.(2)Consideration of the uncertainty in the remaining cycle time has a considerable degree of influence on scheduling performance. The proposed FBPN approach not only enhances the accuracy of estimating the remaining cycle time but can also generate a precise range for the remaining cycle time. This effectively avoids incorrect scheduling.(3)The experimental results in this study also confirmed that the scheduling performance of a complex production system can be significantly improved by optimizing the existing dispatching rule with an FNLP model.
However, the proposed methodology only optimizes the performance of a dispatching rule in a wafer fabrication factory; this does not optimize every aspect of scheduling for a wafer fabrication factory. Although slack diversification optimizes the fourobjective fluctuation smoothing rule, there are other approaches that can achieve the same effect. Further, some soft computing techniques can be applied to solve the FNLP problem. All of these issues constitute directions for future research.
Abbreviations
4oSDR:  Fourobjective slackdiversifying 
ANFIS:  Adaptive neurofuzzy inference system 
COG:  Center of gravity 
CR:  Critical ratio 
CV:  Coefficient of variance 
DOE:  Design of experiment 
EDD:  Earliest due date 
FBPN:  Fuzzy back propagation network 
FCM:  Fuzzy cmeans 
FLP:  Fuzzy linear programming 
FMS:  Flexible manufacturing system 
FNLP:  Fuzzy nonlinear programming 
FSMCT:  The fluctuation smoothing policy for mean cycle time 
FSVCT:  The fluctuation smoothing policy for cycle time variation 
GA:  Genetic algorithm 
LP:  Linear programming 
MIP:  Mixed integer programming 
NFS:  Nonlinear fluctuation smoothing 
NLP:  Nonlinear programming 
PSO:  Particle swarm optimization 
SA:  Simulated annealing 
SRPT:  Shortest remaining processing time 
TFN:  Triangular fuzzy number 
TrFN:  Trapezoidal fuzzy number 
TSK:  TakagiSugenoKang. 
Acknowledgments
The authors thank Richard Romanowski for his valuable contribution to this study. This work was supported by the National Science Council of Taiwan.
References
 C. Grimme and J. Lepping, “Combining basic heuristics for solving multiobjective scheduling problems,” in Proceedings of IEEE Symposium on Computational Intelligence in Scheduling, pp. 9–16, 2011. View at: Google Scholar
 T. Murata, H. Oshida, and M. Gen, “Rulebased weight definition for multiobjective fuzzy scheduling with the OWA operator,” in Proceedings of the 26th Annual Conference of the IEEE Electronics Society (IECON '00), vol. 4, pp. 2756–2761, October 2000. View at: Google Scholar
 M. Dong and M. Liu, “An ANFISbased dispatching rule for complex fuzzy job shop scheduling problem,” in Proceedings of the International Conference on Information Science and Technology (ICIST '11), pp. 263–266, March 2011. View at: Publisher Site  Google Scholar
 T. Murata, G. Mitsuo, and H. Ishibuchi, “Multiobjective scheduling with fuzzy duedate,” Computers and Industrial Engineering, vol. 35, no. 34, pp. 439–442, 1998. View at: Publisher Site  Google Scholar
 H. H. Xiong, M. Zhou, and C. N. Manikopoulos, “Scheduling flexible manufacturing systems based on timed petri nets and fuzzy dispatching rules,” in Proceedings of the INRIA/IEEE Symposium on Emerging Technologies and Factory Automation (ETFA '95), vol. 3, pp. 309–315, October 1995. View at: Google Scholar
 A. X. Benincasa, O. Morandin, and E. R. R. Kato, “Reactive fuzzy dispatching rule for automated guided vehicles,” in Proceedings of IEEE International Conference on Systems, Man and Cybernetics, vol. 5, pp. 4375–4380, October 2003. View at: Google Scholar
 E. H. Mamdani, “Application of fuzzy logic to approximate reasoning using linguistic synthesis,” IEEE Transactions on Computers, vol. 26, no. 12, pp. 1182–1191, 1977. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 K. K. Lee, W. C. Yoon, and D. H. Baek, “Generating interpretable fuzzy rules for adaptive job dispatching,” International Journal of Production Research, vol. 39, no. 5, pp. 1011–1030, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 K. K. Tan and K. Z. Tang, “Vehicle dispatching system based on Taguchituned fuzzy rules,” European Journal of Operational Research, vol. 128, no. 3, pp. 545–557, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 D. Dubois, H. Fargier, and P. Fortemps, “Fuzzy scheduling: modelling flexible constraints vs. coping with incomplete knowledge,” European Journal of Operational Research, vol. 147, no. 2, pp. 231–252, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. C. Biggs and M. A. Laughton, “Optimal electric power scheduling: a large nonlinear programming test problem solved by recursive quadratic programming,” Mathematical Programming, vol. 13, no. 2, pp. 167–182, 1977. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. L. Pedro, Scheduling Theory, Algorithms, and Systems, Springer, New York, NY, USA, 2008.
 H. Ishibuchi and T. Murata, “A multiobjective genetic local search algorithm and its application to flowshop scheduling,” IEEE Transactions on Systems, Man and Cybernetics C, vol. 28, no. 3, pp. 392–403, 1998. View at: Publisher Site  Google Scholar
 V. Kumar, N. Mishra, F. T. S. Chan, and A. Verma, “Managing warehousing in agile supply chain environment: an FAIS algorithm based approach,” International Journal of Production Research, vol. 49, no. 21, pp. 6407–6426, 2011. View at: Publisher Site  Google Scholar
 W. Xia and Z. Wu, “An effective hybrid optimization approach for multiobjective flexible jobshop scheduling problems,” Computers and Industrial Engineering, vol. 48, no. 2, pp. 409–425, 2005. View at: Publisher Site  Google Scholar
 T. C. Chiang and L. C. Fu, “Using a family of critical ratiobased approaches to minimize the number of tardy jobs in the job shop with sequence dependent setup times,” European Journal of Operational Research, vol. 196, no. 1, pp. 78–92, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 D. Connors, G. Feigin, and D. Yao, “Scheduling semiconductor lines using a fluid network model,” IEEE Transactions on Robotics and Automation, vol. 10, no. 2, pp. 88–98, 1994. View at: Publisher Site  Google Scholar
 H. Chen and D. D. Yao, “Dynamic scheduling of a multiclass fluid network,” Operations Research, vol. 41, no. 6, pp. 1104–1115, 1993. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Ishibuchi, N. Yamamoto, T. Murata, and H. Tanaka, “Genetic algorithms and neighborhood search algorithms for fuzzy flowshop scheduling problems,” Fuzzy Sets and Systems, vol. 67, no. 1, pp. 81–100, 1994. View at: Publisher Site  Google Scholar  MathSciNet
 M. Saidi Mehrabad and A. Pahlavani, “A fuzzy multiobjective programming for scheduling of weighted jobs on a single machine,” International Journal of Advanced Manufacturing Technology, vol. 45, no. 12, pp. 122–139, 2009. View at: Publisher Site  Google Scholar
 D. F. Li, “Closeness coefficient based nonlinear programming method for intervalvalued intuitionistic fuzzy multiattribute decision making with incomplete preference information,” Applied Soft Computing Journal, vol. 11, no. 4, pp. 3402–3418, 2011. View at: Publisher Site  Google Scholar
 M. Fares and B. Kaminska, “FPAD: a fuzzy nonlinear programming approach to analog circuit design,” IEEE Transactions on ComputerAided Design of Integrated Circuits and Systems, vol. 14, no. 7, pp. 785–793, 1995. View at: Publisher Site  Google Scholar
 T. Chen and Y.C. Wang, “A fuzzy set approach for evaluating and enhancing the midterm competitiveness of a semiconductor factory,” Fuzzy Sets and Systems, vol. 160, no. 5, pp. 569–585, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 T. Chen, “Job remaining cycle time estimation with a postclassifying fuzzyneural approach in a wafer fabrication plant: a simulation study,” Proceedings of the Institution of Mechanical Engineers B, vol. 223, no. 8, pp. 1021–1031, 2009. View at: Publisher Site  Google Scholar
 T. Chen, “The optimizedrulefusion and certainrulefirst approach for multiobjective job scheduling in a wafer fabrication factory,” International Journal of Innovative Computing, Information and Control, vol. 9, no. 6, pp. 2283–2302, 2013. View at: Google Scholar
 S. C. H. Lu, D. Ramaswamy, and P. R. Kumar, “Efficient scheduling policies to reduce mean and variance of cycletime in semiconductor manufacturing plants,” IEEE Transactions on Semiconductor Manufacturing, vol. 7, no. 3, pp. 374–388, 1994. View at: Publisher Site  Google Scholar
 Y. C. Wang, T. Chen, and C. W. Lin, “A slackdiversifying nonlinear fluctuation smoothing rule for job dispatching in a wafer fabrication factory,” Robotics and ComputerIntegrated Manufacturing, vol. 29, no. 3, pp. 41–47, 2013. View at: Google Scholar
 T. Chen, “A fuzzy midterm singlefab production planning model,” Journal of Intelligent Manufacturing, vol. 14, no. 34, pp. 273–285, 2003. View at: Publisher Site  Google Scholar
 T. Chen, Y. C. Wang, and H. R. Tsai, “Lot cycle time prediction in a rampingup semiconductor manufacturing factory with a SOMFBPNensemble approach with multiple buckets and partial normalization,” International Journal of Advanced Manufacturing Technology, vol. 42, no. 1112, pp. 1206–1216, 2009. View at: Publisher Site  Google Scholar
 X. L. Xie and G. Beni, “A validity measure for fuzzy clustering,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 13, no. 8, pp. 841–847, 1991. View at: Publisher Site  Google Scholar
 T. Chen and Y. C. Wang, “Incorporating the FCMBPN approach with nonlinear programming for internal due date assignment in a wafer fabrication plant,” Robotics and ComputerIntegrated Manufacturing, vol. 26, no. 1, pp. 83–91, 2010. View at: Publisher Site  Google Scholar
 T. Chen and Y. C. Lin, “A collaborative fuzzyneural approach for internal due date assignment in a wafer fabrication plant,” International Journal of Innovative Computing, Information and Control, vol. 7, no. 9, pp. 5193–5210, 2011. View at: Google Scholar
 J. Nocedal and S. J. Wright, Numerical Optimization, Springer, Berlin, Germany, 2006. View at: MathSciNet
 E. S. Lee and R.J. Li, “Comparison of fuzzy numbers based on the probability measure of fuzzy events,” Computers & Mathematics with Applications, vol. 15, no. 10, pp. 887–896, 1988. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C.H. Cheng, “A new approach for ranking fuzzy numbers by distance method,” Fuzzy Sets and Systems, vol. 95, no. 3, pp. 307–317, 1998. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 T.C. Chu and C.T. Tsao, “Ranking fuzzy numbers with an area between the centroid point and original point,” Computers & Mathematics with Applications, vol. 43, no. 12, pp. 111–117, 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L. Zhu and R. Xu, “Ranking fuzzy numbers based on fuzzy mean and standard deviation,” in Proceedings of the 8th International Conference on Fuzzy Systems and Knowledge Discovery, pp. 854–857, 2011. View at: Google Scholar
 F. Wilcoxon, “Individual comparisons by ranking methods,” Biometrics Bulletin, vol. 1, no. 6, pp. 80–83, 1945. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 HorngRen Tsai and Toly Chen. 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.