Uncertain Production Scheduling Based on Fuzzy Theory Considering Utility and Production Rate

Handling uncertainty in an appropriate manner during the real operation of a cyber-physical system (CPS) is critical. Uncertain production scheduling as a part of CPS uncertainty issues should attract more attention. In this paper, a Mixed Integer Nonlinear Programming (MINLP) uncertain model for batch process is formulated based on a unit-specific event-based continuous-time modeling method. Utility uncertainty and uncertain relationship between production rate and utility supply are described by fuzzy theory. The uncertain scheduling model is converted into deterministic model by mathematical method. Through one numerical example, the accuracy and practicability of the proposed model is proved. Fuzzy scheduling model can supply valuable decision options for enterprise managers to make decision more accurate and practical. The impact and selection of some key parameters of fuzzy scheduling model are elaborated.


Introduction
A cyber-physical system (CPS) was identified as a key research area by the US National Science Foundation (NSF) and was listed as the number one research priority by the US President's Council of Advisors on Science and Technology [1].In recent years, CPS is attracting much attention and is being considered as an emerging technology.It has wide potential applications in defense [2,3], transportation [4], industrial automation [5], energy [6,7], health care [1,8] and agriculture [9,10], etc. Cyber-physical systems (CPS) can be viewed as a new generation of systems with integrated control, communication and computational capabilities [1].It combines computation and communication capabilities with the physical world [11].
Uncertainty is intrinsic in most technical systems, including CPS.CPS generally functions in an inherently complex and unpredictable physical environment.A major difficulty with CPS is that it must be designed and operated in the presence of uncertainty [12].Therefore, handling uncertainty in an appropriate manner during the real operation of CPS is critical.In CPS, physical components are inserted by sensors and actuators.Sensors have the abilities of self-calibration, dynamic parameter setting, Ethernet-based communication and self-adaptation [13].So CPS has enough capability to detect uncertainties in itself or from external environment.
For a manufacturing enterprise, the real operation of CPS is just the execution of production scheduling.So the uncertainty problem of production scheduling is a part of the CPS uncertainty issue.Uncertainties in manufacturing are various, which include demand uncertainty, processing time uncertainty, utility uncertainty, due date uncertainty, equipment failure and the absence of employee.For the uncertain scheduling problem of batch process, some research results have been obtained.Guillen, et al. [14] proposed a multistage stochastic optimization model to address the scheduling of supply chains with embedded multipurpose batch chemical plants under demand uncertainty.Zhu, et al. [15] proposed a robust optimization approach for the short-term scheduling of batch plants under demand uncertainty where the uncertain parameters can be described by a normal distribution function.Ierapetritou [16] developed a strategy to quantify scheduling robustness in the face of uncertainty, to increase scheduling flexibility, and to improve system performance when unexpected events occur during scheduling execution.Xu, et al. [17] established a scheduling mathematical model for multi-product batch processes under finite intermediate storage policy with uncertain processing time based on fuzzy programming theory.
Nowadays, most studies mainly deal with uncertainties of demand and processing time in production scheduling.Research on utility uncertainty is rare, and the supply amount of utility is assumed to be constant.In fact, utility disturbance is an important uncertain factor to the production process.In the production process, utilities such as steam, cooling water and electricity, are supplied to areas or product lines to support the normal operation of equipment.Sometimes, these utilities are shared by different areas.However, the utility supply is affected seriously by the energy system, and the energy system is also disturbed by many uncertainty factors.So, utility disturbances usually occur, which will affect the normal running of machines, make the scheduling decision infeasible, even leads to the shutdown of the product line [18][19][20].Consequently, dealing with utility uncertainty well in production scheduling is vital to enterprise competitiveness.Moreover, the production rate is closely related to the utility supply.The relationship between them is roughly assumed to be linear in most studies, which is unable to reflect the actual production situation.So building the uncertain relationship between the production rate and the utility supply is another key point for production scheduling.
For the above two uncertainties, it is time-consuming and material-consuming to establish a mechanism model [21].Existing methods to describe uncertainty include chance constraint programming, fuzzy theory and scenario trees [22,23].Chance constraint programming needs a large amount of historical data to analyze the distribution function of uncertainty [24].The scenario trees method adds many constraints, enables a large programming model scale, and has a computation time too long to meet requirements [25].If the historical data of uncertainty is insufficient or difficult to acquire, then fuzzy theory is an effective means.
In this paper, a MINLP uncertain scheduling model for batch process is formulated based on unit-specific event-based continuous-time (USEBCT) modeling method.Fuzzy method is utilized to describe the uncertainty of utility supply.For the uncertain relationship between the production rate and the utility supply, firstly establish the linear function in deterministic scheduling model, and then use fuzzy method to represent the coefficient of linear function.The uncertain scheduling model is converted into deterministic model by mathematical method.Through one numerical example, the accuracy and practicability of the proposed model are proved.
This paper is organized as follows: Section 2 describes the problem definition of short scheduling for batch process.Section 3 formulates the deterministic scheduling model.Section 4 details the fuzzy scheduling model and the defuzzification method.A numerical example and roles of proposed model are given in Section 5. Conclusions are presented in Section 6.

Problem State
The short-term scheduling problem for batch processes is stated as follows.Given (i) the production recipe (i.e., the amount of the materials required for the production of each product); (ii) the available units, their capacity limits and utility; (iii) the available storage capacity for each of the materials; and (iv) the time horizon under consideration, then the objective is to determine (i) the optimal sequence of tasks taking place in each unit; (ii) the amount of material being processed at each time in each unit; and (iii) the processing time of each task in each unit.Some assumptions in this paper are given as follows: (i) the switching time between different tasks is neglected [26]; (ii) raw materials are sufficient and supplied at the beginning of the scheduling horizon [27]; (iii) intermediates and final products are only from producing [28][29][30].

Deterministic Scheduling Model Based on USEBCT Modeling Method
The USEBCT model is proposed by Ierapetritou and Floudas [31].The USEBCT model permits event to happen at any moment in the scheduling horizon.The duration time of every event is varying.The scale of mathematical programming is smaller, and the computation time is shorter, so compared with discrete-time model, USEBCT model is more actual and accurate in describing the production situation.

Basic Concepts of USEBCT Model
(1) Event: time measurement, each event has one start time and finish time; (2) State: raw materials, intermediates and final products; (3) Task: action of producing or consuming one state; (4) Unit: machines in which the task can be performed.

Deterministic Scheduling Model
The proposed deterministic scheduling model is formulated based on the model framework proposed by Vooradi and Shaik (V&S model) [26].For the reason of computation time, ∆n is set to be zero.Moreover, material balances, duration constraints, tightening constraints and utility constraints are modified according to the requirements of modeling.
(1) Allocation constraints In every unit, at the most one task can be active at each event as given by constraint (1), in which w(i, n) is binary variable.If task i is active at event n, w(i, n) = 1, else w(i, n) = 0.
(2) Capacity constraints Constraint (2) limits the size of one batch.If task i is active at event n, then the decision variable of processed material b(i, n) is forced in the range of the minimum and maximum batch sizes.Otherwise, b(i, n) is forced to be zero.
(3) Material balances When n > 1, Constraints ( 3) and ( 4) are material balances of raw materials, intermediates and products at initial event.Constraints ( 5)-( 7) are material balances of the above three states at the later events.If n is equal to 1, the amount of states (raw material or intermediate) at event n is equal to the initial value, and decreased by any amounts consumed at event n, as given by constraint (3).The amount of final product at event n is calculated as initial value, as given by constraint (4).
If n > 1, the amount of a state at current event is adjusted by any amounts produced by tasks that ends at the previous event and adjusted by any amounts consumed by the tasks that starts at the current event, as given by constraint (5) to constraint (7).

(4) Duration constraints
In order to formulate the relationship between consumption of utility and production rate, the variable P(i, n) is introduced to represent the production rate of task i at event n.The finish time of task i at event n is equal to the start time plus duration time of the task, which is calculated by b(i, n)/P(i, n).
In constraint (9), production rate P(i, n) is formulated by the summery of minimum production rate, which maintains the equipment running, and linear function of consumption of utility.The parameter c ui is assumed as a constant number in deterministic model.However the relationship between consumption of utility and production rate is not linear simply.In fuzzy model, the parameter c ui is defined as a fuzzy number to describe the uncertain relationship between consumption of utility and production rate as given by constraint (25).
(5) Sequencing constraints 1 Same task in the same unit In constraint (10), the start time of task i at event n + 1 should be greater than or equal to the finish time of the same task at event n.
Different tasks in the same unit The constraint ( 11) is written for two different tasks performed in the same unit at event n and n + 1, respectively.If task i is active at event n + 1, the start time of task i must be greater than or equal to the finish time of previous task i .3 Different tasks in different units Information 2017, 8, 158 5 of 15 According to the production recipe, for different tasks that produce or consume the same state, the start time of the consuming task at the next event is enforced to be later than the finish time of the producing task at the current event as given by constraint (12).Otherwise, the constraint is relaxed and the two times are not related.
(6) Tightening constraints The sum of the durations of all tasks suitable in each unit should be less than the scheduling horizon as given by constraint ( 13).( 7) Utility constraints The consumption of utility is a linear function of amount of processed material, which is given by constraint (14).The consumption of utility is limited by the total supply amount of utility.Constraint (15) expresses that the consumption of utility u at event n should not be greater than the maximum supply.The maximum supply amount of utility U max ui is another uncertain parameter and assumed as a constant number in deterministic model.But the supply amount of utility is affected seriously by energy system, and the energy system is also disturbed by lots of uncertain factors.So, utility disturbances usually occur.In fuzzy model, the parameter U max ui is defined as a fuzzy number to describe the uncertain utility supply as given by constraint (26).
Constraint ( 16) is similar to that for same task in the same unit written for each utility.The start time of consuming utility u at event n + 1 should not be earlier than that at event n.
To account for the utility consumption over different tasks, constraints ( 17) and ( 18) enforce the start time of all suitable tasks that consume the same utility to be equal to T s u (u, n), if the task is active at event n.
In constraint (19), the finish times of all suitable tasks that consume a utility at previous event are enforced to occur before, if the task is active.Otherwise, constraint (19) is relaxed.
Information 2017, 8, 158 6 of 15 Infinite storage and dedicated finite intermediate storage policies are considered in this work.Constraint (20) expresses the limit of storage capacity.For handling dedicated finite intermediate storage case, constraint (21) along with constraint ( 12) is used to avoid real time storage violations without the need for considering storage as a separate task.
If constraints ( 12) and ( 21) are not relaxed, constraint ( 22) is workable.In constraint (22), the finish time of producing task at event n should be greater than or equal to the start time of consuming task at the same event, and less than or equal to the start time of consuming task at event n + 1.
(9) Bounds on variables Constraint ( 23) initials the consumption time of utility at the first event, and limits all the start and finish times of tasks to be within the scheduling horizon.
The objective function is the maximization of profit, and the total amount of the final products produced by the end of last event is calculated by constraint (24).

Uncertainty Description
Coefficient c ui and maximum supply amount of utility U max ui are fuzzy numbers, as given by constraints (25) and (26).The two fuzzy numbers are described by triangular possibility distributions.This description of fuzzy numbers has been extensively used in literatures due to their various advantages (e.g., Lai and Hwang, 1992 [32], 1993 [33], Liang, 2006 [34] and 2008 [35]).The simplicity in required data acquisition and simplicity in related fuzzy computations are the most important advantages [36][37][38].These advantages are our main motivation in using the triangular fuzzy numbers for modeling the imprecise data in our problem.Figure 1 presents the triangular possibility distribution of fuzzy number A = A P , A m , A o , where A P , A m and A o are the most pessimistic value, the most possible value, and the most optimistic value.The triangular possibility distribution of A is defined as Equation (27).
The universe U of fuzzy number is defined as Equation ( 28) according to Basyigit and Ulu [39], where D min is the minimum value in statistical data of fuzzy number, D max is the maximum value in statistical data of fuzzy number.D 1 and D 2 are defined as Equation ( 29) and (30) according to the experience of engineering.D m f is the value of which the frequency of occurrence is the most in statistical data of fuzzy number.
Typically, the three prominent values, i.e., the most pessimistic, the most possible, and the most optimistic values for fuzzy numbers are estimated by decision maker [32,37,40].In this paper, the three prominent values are defined according to determining method of universe [39] and the experience of engineering.The most pessimistic value A P is equal to the left boundary of universe.The most possible value A m is equal to D m f .The most optimistic value A O is equal to the right boundary of universe.These relationships are given by Equations ( 31) to (33).
Information 2017, 8, 158 7 of 15 The universe U of fuzzy number is defined as Equation ( 28) according to Basyigit and Ulu [39], where min D is the minimum value in statistical data of fuzzy number, max D is the maximum value in statistical data of fuzzy number. 1 D and 2 D are defined as Equation ( 29) and ( 30) according to the experience of engineering.mf D is the value of which the frequency of occurrence is the most in statistical data of fuzzy number.Typically, the three prominent values, i.e., the most pessimistic, the most possible, and the most optimistic values for fuzzy numbers are estimated by decision maker [32,37,40].In this paper, the three prominent values are defined according to determining method of universe [39] and the experience of engineering.The most pessimistic value P A is equal to the left boundary of universe.

Defuzzification of Fuzzy Scheduling Model
Defuzzification is a process whereby the fuzzy scheduling model is converted into a deterministic scheduling model.Constraint (25) can be converted into inequality constraints.All the fuzzy numbers are on one side of inequality constraints, so fuzzy numbers ranking problem is not involved.We apply the well-known weighted average method.This approach seems to be the simplest and most reliable defuzzification method in converting the fuzzy constraints into their deterministic ones [35,37,38,41].In this regard, we also need to determine a minimal acceptable possibility level α , which denotes the minimum acceptable possibility level of occurrence for the corresponding imprecise/fuzzy data.So, the equivalent deterministic constraints can be represented as follows: The deterministic scheduling model is replaced by the following fuzzy scheduling model: Objective function: (24) Constraints: (1)-( 8), ( 10)-( 14), ( 16)-( 23) (25), deterministic constraint ( 9) is replaced by fuzzy constraint (25) (26), deterministic constraint ( 15) is replaced by fuzzy constraint (26).

Defuzzification of Fuzzy Scheduling Model
Defuzzification is a process whereby the fuzzy scheduling model is converted into a deterministic scheduling model.Constraint (25) can be converted into inequality constraints.All the fuzzy numbers are on one side of inequality constraints, so fuzzy numbers ranking problem is not involved.We apply the well-known weighted average method.This approach seems to be the simplest and most reliable defuzzification method in converting the fuzzy constraints into their deterministic ones [35,37,38,41].In this regard, we also need to determine a minimal acceptable possibility level α, which denotes the minimum acceptable possibility level of occurrence for the corresponding imprecise/fuzzy data.So, the equivalent deterministic constraints can be represented as follows: where w a + w b + w c = 1, and w a , w b and w c represent the weights of the most pessimistic, the most possible and the most optimistic value of the related fuzzy numbers, respectively.In practice, the suitable values for these weights as well as α are usually determined subjectively by the experience and knowledge of the decision maker.Based on the most likely value concept proposed by Lai and Hwang [32] and considering several relevant works [42,43], we set these parameters to: w a = 1/6, w b = 4/6, w c = 1/6 in numerical experiments.
α are boundary points of α-cut fuzzy set.Take U max ui as example, boundary points calculated by (37) to (39).The relationships among them are illustrated in Figure 2.
where a b c w w w + + =1 , and a w , b w and c w represent the weights of the most pessimistic, the most possible and the most optimistic value of the related fuzzy numbers, respectively.In practice, the suitable values for these weights as well as α are usually determined subjectively by the experience and knowledge of the decision maker.Based on the most likely value concept proposed by Lai and Hwang [32] and considering several relevant works [42,43], we set these parameters to:  ) Equivalent deterministic model through defuzzification is reformulated as follows: Objective function: (24) Constraints: (1)-( 8), ( 10)-( 14), ( 16)-( 23) and ( 34)- (36).

Case Study
Numerical example is introduced firstly by Maravelias and Grossmann [44] and has been studied in several studies [31,45].As shown in Figure 3, this example involves four tasks in three units processing seven states.The corresponding data is given in Tables 1 and 2. There are two types of reactors available for the process (type I and II) with different number of corresponding units available: two reactors (RI1 and RI2) of type I, while one reactor (RII) of type II.Reactions R1 and R2 require type I reactor, while reactions R3 and R4 require type II reactor.Furthermore, reactions R1 and R3 require heat, provided by steam (HS) produced in limited amounts in the plant, while reactions R2 and R4 are exothermic and require cooling water (CW), also available in limited amounts.Due to safety reasons and temperature restrictions, the heat integration of the process is not possible.

Case Study
Numerical example is introduced firstly by Maravelias and Grossmann [44] and has been studied in several studies [31,45].As shown in Figure 3, this example involves four tasks in three units processing seven states.The corresponding data is given in Tables 1 and 2. There are two types of reactors available for the process (type I and II) with different number of corresponding units available: two reactors (RI1 and RI2) of type I, while one reactor (RII) of type II.Reactions R1 and R2 require type I reactor, while reactions R3 and R4 require type II reactor.Furthermore, reactions R1 and R3 require heat, provided by steam (HS) produced in limited amounts in the plant, while reactions R2 and R4 are exothermic and require cooling water (CW), also available in limited amounts.Due to safety reasons and temperature restrictions, the heat integration of the process is not possible.The minimum batch-size for all tasks is the half of maximum capacities of units.Finite storage policy for intermediates and unlimited storage policy for the raw materials, final products, are applied.Numerical example is solved using Lingo 11.0 (Lindo Systems, Chicago, IL, USA)with a 2.20 GHz computer.Scheduling horizon is seven hours, and event number is set to 5. A better optimization solution cannot be found by larger event numbers.The data of fuzzy numbers are given in Tables 3 and 4. Three prominent values (most pessimistic value, most possible value, and most optimistic value) of fuzzy numbers are in parentheses.
In order to demonstrate the advantages of proposed fuzzy mode, we input a set of actual data into the deterministic model.The results based on the actual data are viewed as a criterion for comparing the accuracy of deterministic model and fuzzy model.Actual data of fuzzy number  ui c and  max ui U are given in Tables 3 and 4.  The minimum batch-size for all tasks is the half of maximum capacities of units.Finite storage policy for intermediates and unlimited storage policy for the raw materials, final products, are applied.Numerical example is solved using Lingo 11.0 (Lindo Systems, Chicago, IL, USA) with a 2.20 GHz computer.Scheduling horizon is seven hours, and event number is set to 5. A better optimization solution cannot be found by larger event numbers.The data of fuzzy numbers are given in Tables 3  and 4. Three prominent values (most pessimistic value, most possible value, and most optimistic value) of fuzzy numbers are in parentheses.
In order to demonstrate the advantages of proposed fuzzy mode, we input a set of actual data into the deterministic model.The results based on the actual data are viewed as a criterion for comparing the accuracy of deterministic model and fuzzy model.Actual data of fuzzy number c ui and U max ui are given in Tables 3 and 4. Results based on actual data, of deterministic model and fuzzy model are compared in Table 5.In this numerical example, w a = 1/6, w b = 4/6, w c = 1/6.The objective function (maximization of profit) based on actual data is 5448$.The corresponding completion time is 11.20 h, which is the sum of processing time of three machine units (RI1, RI2 and RII).The objective function of deterministic model is 5400$.The completion time of deterministic model is 12.30 h.For fuzzy model, the objective function ranges from 5426.25$ to 5410.50$, correspondingly, α varies from 0.5 to 0.8.Completion time also varies depending on different α.In order to compare the optimization results better, deviations of deterministic and fuzzy models as well as average values of deviations for fuzzy model are given in Table 6.Based on data in Table 5, taking results of actual data as criterion, the deviations of deterministic model and fuzzy model are calculated by formula |x−a| a , where x is the measured data and a is truth-value.Results of deterministic and fuzzy models are used as measured data.Results based on actual data are used as truth-value.In the deterministic scheduling model, some uncertain parameters like U max ui and c ui are assumed as constant numbers.In the fuzzy scheduling model, these uncertain parameters are defined as fuzzy numbers, which considers the uncertainties in scheduling problem.Actual data of these uncertain parameters are fluctuant, but not constant.Results based on actual data are optimized through the collected data from actual operating condition, which can reflect the practical situation of production line.So the results closer to that based on actual data should be more practical.For accuracy, in this paper, we view the results based on actual data as the truth-value to evaluate the accuracy of deterministic model and fuzzy model, because it is fairer than just comparing the results of deterministic model and fuzzy model.Consequently, the results closer to that based on actual data should be more accurate.As shown in Table 6, the average values of deviations in terms of objective function and completion time about fuzzy model are 0.545% and 5.09%, which are all smaller than that of deterministic model.It means that the results of fuzzy model are closer to results based on actual data.So the proposed fuzzy model is more practical and accurate than deterministic model.
Fuzzy scheduling model can supply valuable decision options for enterprise managers to make decision more accurate and practical.Decision maker can adjust the acceptable possibility level α, according to risks taken by enterprises, and then obtain more accurate and practical scheduling results.The smaller α decider chooses, the higher risks are taken by enterprise, but the profit will be greater and the scheduling results will be more far away from that of deterministic model.In contrast, the larger α decider chooses, the lower risks are taken by enterprise, but the profit will be less and the scheduling results will be closer to that of deterministic model.In addition to objective function and completion time, Gantt charts for deterministic model and fuzzy model are illustrated in Figures 4 and 5.
more practical.For accuracy, in this paper, we view the results based on actual data as the truth-value to evaluate the accuracy of deterministic model and fuzzy model, because it is fairer than just comparing the results of deterministic model and fuzzy model.Consequently, the results closer to that based on actual data should be more accurate.As shown in Table 6, the average values of deviations in terms of objective function and completion time about fuzzy model are 0.545% and 5.09%, which are all smaller than that of deterministic model.It means that the results of fuzzy model are closer to results based on actual data.So the proposed fuzzy model is more practical and accurate than deterministic model.
Fuzzy scheduling model can supply valuable decision options for enterprise managers to make decision more accurate and practical.Decision maker can adjust the acceptable possibility level α , according to risks taken by enterprises, and then obtain more accurate and practical scheduling results.The smaller α decider chooses, the higher risks are taken by enterprise, but the profit will be greater and the scheduling results will be more far away from that of deterministic model.In contrast, the larger α decider chooses, the lower risks are taken by enterprise, but the profit will be less and the scheduling results will be closer to that of deterministic model.In addition to objective function and completion time, Gantt charts for deterministic model and fuzzy model are illustrated in Figures 4 and 5.It is noteworthy that the membership function of fuzzy number is vital for the quality of fuzzy scheduling results.If membership function is not symmetry and tends to the superior side, the  Fuzzy scheduling model can supply valuable decision options for enterprise managers to make decision more accurate and practical.Decision maker can adjust the acceptable possibility level α , according to risks taken by enterprises, and then obtain more accurate and practical scheduling results.The smaller α decider chooses, the higher risks are taken by enterprise, but the profit will be greater and the scheduling results will be more far away from that of deterministic model.In contrast, the larger α decider chooses, the lower risks are taken by enterprise, but the profit will be less and the scheduling results will be closer to that of deterministic model.In addition to objective function and completion time, Gantt charts for deterministic model and fuzzy model are illustrated in Figures 4 and 5.It is noteworthy that the membership function of fuzzy number is vital for the quality of fuzzy scheduling results.If membership function is not symmetry and tends to the superior side, the It is noteworthy that the membership function of fuzzy number is vital for the quality of fuzzy scheduling results.If membership function is not symmetry and tends to the superior side, the objective function of fuzzy model will be better than that of deterministic model, and vice versa.So one precise membership function analyzed by historical data, is an important support for the accuracy of optimization results of fuzzy scheduling model.Some parameters, for example acceptable possibility level α and weights w a , w b and w c , also have important impact on optimization result of fuzzy model.However these parameters only affect optimization result on quantitative value but not on the quality as membership function.With the increasing of acceptable possibility level α, the optimization result gets close to that of deterministic model.That is because the greater the acceptable possibility level is, the more the fuzzy number is close to the most possible value.The selection of α depends on decision maker.If decision maker doesn't want to take too much risk, the acceptable possibility level should be greater, and vice versa.Weights w a , w b and w c , represent the inclination extent to most pessimistic value, the most possible value, and the most optimistic value in membership function.The three weights are usually selected by experience.

Conclusions
In this paper, a MINLP uncertain model for batch process is formulated based on the unit-specific event-based continuous-time modeling method.Utility uncertainty and uncertain relationship between production rate and utility supply are described by fuzzy theory.The uncertain scheduling model is converted into deterministic model by mathematical method.The accuracy and practicability of the proposed model is proved through one numerical example.Fuzzy scheduling model can supply valuable decision options for enterprise managers to make decision more accurate and practical.The impact and selection of some key parameters of fuzzy scheduling model are elaborated.

Figure 1 .
Figure 1.Triangular possibility distribution of fuzzy number.

Figure 1 .
Figure 1.Triangular possibility distribution of fuzzy number.

Figure 2 .
Figure 2. Triangular possibility distribution of fuzzy number.

Figure 2 .
Figure 2. Triangular possibility distribution of fuzzy number.

B
in kg , δ iHS in kg per kg of batch.

Figure 4 .
Figure 4. Gantt chart for the deterministic model.

Figure 5 .
Figure 5. Gantt chart for the fuzzy model.

Figure 4 .
Figure 4. Gantt chart for the deterministic model.

Figure 5 .
Figure 5. Gantt chart for the fuzzy model.

Table 1 .
Data of production capability and coefficients of utility consumption.

Table 2 .
Data related States.

Table 3 .
Data of fuzzy number  ui c .

Table 1 .
Data of production capability and coefficients of utility consumption.

ffi iHS ffi iCW ffi iHS ffi iCW
B min , B max in kg, δ iHS in kg per kg of batch.

Table 2 .
Data related States.

Table 3 .
Data of fuzzy number c ui .

Table 4 .
Data of fuzzy number U max ui .

Table 5 .
Optimization results of deterministic model and fuzzy model.
For accuracy, in this paper, we view the results based on actual data as the truth-value to evaluate the accuracy of deterministic model and fuzzy model, because it is fairer than just comparing the results of deterministic model and fuzzy model.Consequently, the results closer to that based on actual data should be more accurate.As shown in Table6, the average values of deviations in terms of objective function and completion time about fuzzy model are 0.545% and 5.09%, which are all smaller than that of deterministic model.It means that the results of fuzzy model are closer to results based on actual data.So the proposed fuzzy model is more practical and accurate than deterministic model.