A Fuzzy Nonlinear Programming Approach for Optimizing the Performance of a Four-Objective Fluctuation Smoothing Rule in a Wafer Fabrication Factory

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 four-factor 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 four-factor fluctuation-smoothing rule to sequence the jobs in front of each machine. The fuzzy four-factor fluctuation-smoothing 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 proposedmethodology 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.


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][3][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 Takagi-Sugeno-Kang (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   || max 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][14][15].For example, Chiang and Fu [16] minimized the number of tardy jobs for a job shop in which machines have sequence-dependent 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 NP-hard 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 type-2 FNLP (with type-2 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 multiple-objective FNLP problems to find the optimal sets of circuit parameter values for a bipolar emitter follower circuit and an unbuffered two-stage complementary metaloxide semiconductor (CMOS) op-amp.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][24][25]).Subsequently, the remaining cycle time estimate is fed into a four-factor fluctuation smoothing rule [25] to sequence the jobs in front of a machine.The four-factor 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 four-objective fluctuation-smoothing 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 four-objective 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 four-objective 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.

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) CT  : the cycle time of job .
(3) DD  : the due date of job .
(5) RCT  : the remaining cycle time of job  from step .
(7) RPT  : the remaining processing time of job  from step .
(8) SCT  : the step cycle time of job  until step .
(9) SK  or SK  : the slack of job  at step .
(11) TPT  : the total processing time of job .
(15)    : the connection weight between hidden-layer node  and 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 B, and their  cuts () and (), the arithmetic operations of Ã and B based on their  cuts are as follows: where (+), (−), (×), and (/) denote fuzzy addition, subtraction, multiplication, and division, respectively.Equation ( 5) is equivalent to if   (),   () ̸ = 0.The min( ) and max( ) function can be replaced by Theorem 3 (average of fuzzy numbers based on  cuts).

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 c-means (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 ;  ∈ [1, ∞) 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 three-layer FBPN is used to estimate the cycle times of jobs for each category.The configuration of the three-layer 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 three-layer FBPN is the (normalized) estimated remaining cycle time (( RCTE  )) 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   2 , the FBPN is treated as a crisp network.Some algorithms are applicable for this purpose, such as gradient descent algorithms, conjugate gradient algorithms, the Levenberg-Marquardt algorithm, and others.In this study, the Levenberg-Marquardt algorithm is applied.The Levenberg-Marquardt algorithm was designed for training with second-order speed without having to compute the Hessian matrix.It uses approximation and updates the network parameters in a Newton-like way [33].
Subsequently,   3 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,   1 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 Four-Objective 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 Ã = ( 0 − ,  0 ,  0 + ) can be derived as Proof.Refer to Zhu and Xu [37].It is, in fact, the center-ofgravity (COG) method.
The following definition details a method based on the fuzzy mean and standard deviation.Definition 5.For any two fuzzy numbers Ã and B ∈ (), the sequence of Ã and B can be determined according to their fuzzy means and standard deviations as follows. (1) Consider the example in Table 2.The sequencing results by the two fuzzified rules are 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 four-objective dispatching rule.In the four-objective dispatching rule, the slack of job  at processing step  is defined as where , , , and  and are positive real numbers that satisfy the following constraints: ( The values of  and  are within [0 1].If the remaining cycle time is estimated with a triangular fuzzy number, then (22) becomes Job  is processed before job  if SK  < SK  .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  SK  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 S  is Subsequently, the objective function is equal to The -cut of Z1 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 Z1 , using the center-of-gravity 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  * 2 is 187272 when the parameters , , , , and  are equal to 0.980, 0.560, −1.100, −1.100, and −2.082, 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  * SK  = 93.The optimal values of the parameters , , , , and  are equal to 0.938, 0.587, −1.112, −1.112, and −2.050, 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 right-hand side, then the slack of the job will be underestimated.Conversely, the slack will be overestimated if there is a long tail on the left-hand side.

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 real-world 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 single-wafer 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 on-time delivery to customers.One hundred replications of the simulation were successively run.The simulation horizon of each replication was twenty-four months.The warm-up 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 four-objective slack-diversifying rule (4o-SDR) [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 4o-SDR, the remaining cycle time of a job was estimated using the fuzzy cmeans and back propagation network (FCM-BPN) approach [31]; it was a crisp value.The five adjustable parameters were set to (, , , , ) = (0.6, 0.2, 0, −0.6, 0) 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 four-objective 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 four-objective 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.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. 0 : When shortening the average cycle time, the scheduling performance of the proposed methodology is the same as that of the existing approach being compared. 1 : When shortening the average cycle time, the scheduling performance of the proposed methodology is better than that of the existing approach being compared. 0 : When reducing the maximum lateness, the scheduling performance of the proposed methodology is the same as that of the existing approach being compared.
1 : When reducing the maximum lateness, the scheduling performance of the proposed methodology is better than that of the existing approach being compared.
0 : 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.
1 : When reducing the cycle time standard deviation, the scheduling performance of the proposed methodology is better than that of the existing approach being compared.
0 : 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.
1 : 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  0 was rejected at  = 0.025, which showed that the fuzzy slack-diversifying fluctuation-smoothing rule was superior to seven existing approaches at reducing the average cycle time.With regard to maximum lateness, the advantage of the fuzzy slack-diversifying fluctuation-smoothing rule over FIFO, SRPT, and FSVCT was significant.Similar results could be observed with cycle time standard deviation.However, the advantage of the fuzzy slack-diversifying fluctuationsmoothing rule was not statistically significant for the number of tardy jobs.

Conclusions and Directions for
Future Research 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 four-objective 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.Takagi-Sugeno-Kang.

FuzzifyFigure 1 :
Figure 1: The flowchart of the proposed methodology.

Figure 2 :
Figure 2: The estimation results by the effective FBPN approach.

Figure 3 :
Figure 3: The estimation results by the statistical analysis approach.

Figure 4 :
Figure 4: The fuzzy slacks obtained by the proposed methodology.

Figure 5 :
Figure 5: The fuzzy slacks obtained by using Wang et al. 's method.

Table 1 :
Some fuzzy scheduling methods.

Table 4 :
The performance levels of various approaches for average cycle time.

Table 5 :
The performance levels of various approaches for maximum lateness.

Table 6 :
The performance levels of various approaches for cycle time standard deviation.

Table 7 :
The performance levels of various approaches for number of tardy jobs.