Next Article in Journal
Learning Output Reference Model Tracking for Higher-Order Nonlinear Systems with Unknown Dynamics
Previous Article in Journal
The Role of Façade Materials in Blast-Resistant Buildings: An Evaluation Based on Fuzzy Delphi and Fuzzy EDAS
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integration of Production Planning and Scheduling Based on RTN Representation under Uncertainties

1
School of Information and Control Engineering, Liaoning Shihua University, Fushun 113001, China
2
National Experimental Teaching Demonstration Center of Petrochemical Process Control, Liaoning Shihua University, Fushun 113001, China
3
Institute of Intelligence Science and Engineering, Shenzhen Polytechnic, Shenzhen 518055, China
*
Author to whom correspondence should be addressed.
Algorithms 2019, 12(6), 120; https://doi.org/10.3390/a12060120
Submission received: 19 April 2019 / Revised: 2 June 2019 / Accepted: 6 June 2019 / Published: 10 June 2019

Abstract

:
Production planning and scheduling are important bases for production decisions. Concerning the traditional modeling of production planning and scheduling based on Resource-Task Network (RTN) representation, uncertain factors such as utilities are rarely considered as constraints. For the production planning and scheduling problem based on RTN representation in an uncertain environment, this paper formulates the multi-period bi-level integrated model of planning and scheduling, and introduces the uncertainties of demand and utility in planning and scheduling layers respectively. Rolling horizon optimization strategy is utilized to solve the bi-level integrated model iteratively. The simulation results show that the proposed model and algorithm are feasible and effective, can calculate the consumption of utility in every period, decrease the effects of uncertain factors on optimization results, more accurately describe the uncertain factors, and reflect the actual production process.

1. Introduction

Production planning and scheduling are important bases for production decisions, and play an important role in economic benefits. In the traditional control of production processes, production planning and scheduling are optimized, respectively, which has obvious disadvantages such as the infeasibility and sub-optimality of optimization results. With the rapid development of economic globalization and product demand diversification, competition among enterprises is more and more fierce. In order to arrange production reasonably and enhance competitiveness of enterprises, enterprise managers and researchers realize that single production planning or scheduling cannot satisfy the needs of production and management in enterprises. Consequently, the integration problem of production planning and scheduling is receiving more and more attention [1] in the complex and changeable economic situation [2,3].
Shao et al. [4] proposed that production planning and scheduling are complementary, and that integration of production planning and scheduling can enhance the performance of production management greatly. Lasserre [5] formulated a holistic integrated model of production planning and scheduling, and solved the integrated model hierarchically. Xiong et al. [6] explored multi-period multi-workshop production planning and scheduling problems, and formulated an integrated optimization model of the production planning and scheduling based on nonlinear mixed integer programming. An alternant iterative method by hybrid genetic algorithm is employed to solve the integrated model, which promoted the research of the integration problem of production planning and scheduling. Li et al. [7] stated that implementing integrated optimization of planning and scheduling can improve the quality of decision making in the process operations, and also proposed a new decomposition strategy based on two-level optimization. Wang et al. [8] formulated a bi-level integrated model of production planning and scheduling based on State-Task Network (STN) representation, and introduced the rolling horizon optimization strategy to solve the bi-level integrated model iteratively. The feasibility and effectiveness of the proposed model are illustrated through benchmark examples. Vogel et al. [9] developed a hierarchical and integrated model considering both levels, aggregate production planning, and master production scheduling. Computational tests show that it is possible to solve the integrated model and it outperforms the existing hierarchical approach for all instances. Hassani et al. [10] proposed a new model integrating planning and scheduling to make up for the defect that the capacity limitation could not reflect the practical availability of resources, and added the constraint of resource availability. In this model, the optimization total cost of a single-level job shop was taken as the objective function, and the random approximation method was adopted as the genetic algorithm to solve the N-P problem.
A reasonable and elaborate representation for production processes is an important point to establish production planning and scheduling model. In the process industry, general representation methods for production process include STN, State-Sequence Network (SSN) and Resource-Task Network (RTN). Pantelides et al. [11] proposed the RTN representation method, which is a comprehensive and concise process-representation method. This representation method can describe all the resource properties in a unified manner. Any kinds of facilities and other resources are not distinctive, and considered as producible and consumable resources when the tasks are active. Zhang et al. [12] extended the concept of RTN to provide a unified mathematical formulation of the problem of determining the optimal operating conditions of a mixed production facility, comprising a multipurpose plant for both batch and continuous operations. Alternative approaches for solving the general problem were discussed, and some preliminary numerical results were presented. Castro et al. [13] proposed an improved general mathematical programming formulation for optimal scheduling of batch processes based on the RTN representation. The proposed models led to simpler and less degenerate mathematical models. These models can be solved in significantly less CPU time, when compared to other RTN continuous-time formulations. Chen et al. [14] proposed an improved mixed integer linear programming model for short-term scheduling of multipurpose batch plants under maximization of profit based on RTN representation and unit-specific events. To solve the model, a hybrid method based on line-up competition algorithm and linear programming was presented. Simulation results show that the proposed model and hybrid method were effective for short-term scheduling of multipurpose batch plants. Some research achievements on the production planning and scheduling problem based on RTN representation were proposed, however, the majority of these achievements are under determinate condition, and do not involve uncertain factors [15,16]. Chao et al. [17] believed that distributionally robust optimization (DRO) is a new and effective method to deal with the imprecise probability distribution of uncertain parameters in the decision-making process. They proposed an effective DRO framework to solve the production planning and scheduling problem in the uncertain demand environment. Based on this, they considered the multi-stage sequential decision structure in the operation process. Curcio et al. [18] proposed an adaptive strategy applicable to the multi-stage setting of the two-stage planning and robust optimization model, and an approximate heuristic algorithm based on multistage stochastic programming and adjustable robust optimization. The simulation results show that the strategies are promising in solving large-scale problems: the approximate strategy based on adjustable robust optimization has, on average, 6.72% better performance and is 7.9 times faster than the deterministic model. Production planning and scheduling problem based on RTN representation considering uncertain factors as constraints in production process, is rare. This paper formulates the multi-period bi-level integrated model of production planning and scheduling, considers the demand and utility uncertainties for the production planning and scheduling problem based on RTN presentation. A deterministic model with two layers of planning and scheduling is formulated firstly. In the deterministic model, uncertainties of demand and utility are treated as deterministic parameters in the planning layer and scheduling layer respectively. Then the chance constrained programming [19] and fuzzy theory [20,21] are introduced into the planning layer model and scheduling layer model, respectively, to describe the uncertainties of demand and utility [22,23]. The planning layer model is formulated by the discrete-time modeling method [24]. The scheduling layer model is formulated by the unit-specific event-based continuous-time modeling method. Rolling horizon optimization strategy [25] is utilized to solve the bi-level integrated model iteratively. The simulation results show that the model and algorithm can accurately describe the uncertain factors, decrease the effects of uncertain factors on optimization results, and more accurately reflect the actual production processes.
This paper is organized as follows: Section 1 introduces the research status of production planning and scheduling integration. Section 2 establishes the deterministic model and the uncertain model at the planning layer and the scheduling layer respectively. Section 3 introduces the rolling horizon optimization strategy. An example is simulated and analyzed in Section 4. Conclusions are presented in Section 5.

2. Mathematical Model

2.1. Planning Layer Model

Due to the long time scale of the planning layer, it can be divided evenly into several time periods. A discrete-time linear programming model is established at the planning layer, with the maximization of profits as the objective function. Chance constrained programming is introduced to describe the uncertainty of demand. Through a series of mathematical transformations, the uncertain model with chance constraints is transformed into a deterministic model with a specified confidence level. At last, the transformed deterministic model is solved.

2.1.1. Deterministic Model

1. Objective function:
m a x i G p k K [ P r i c e i D ( i , k ) α P r i c e i T d ( i , k ) β F l u ( i , k ) θ I n v ( i , k ) ] .
The maximum of profit is taken as the objective function of the production planning layer model, where the first part represents the product profit, the second part represents the punishment cost caused by the backlog order quantity of the products, the third part represents the cost paid due to the impact of the productivity fluctuation in the production processes, and the fourth part represents the cost paid when the inventory quantity exceeds the safe limit of inventory.
2. Material balance constraints:
S ( i , k ) = S ( i , k 1 ) + P ( i , k ) D ( i , k ) j N i c i j P ( j , k ) ; i G a l l   ,   k > 1   ,   k K .
S ( i , k ) = S i i n i + P ( i , k ) D ( i , k ) j N i c i j P ( j , k )   ; i G a l l   ,   k = 1   ,   k K .
Constraints (2) and (3) respectively represent the material balance in the subsequent period and the initial period. The inventory s ( i , k ) of the material status   i at the end of the planning period k is equal to the inventory at the end of the previous period s ( i , k 1 ) plus the output of current period p ( i , k ) , and finally minus the delivery of current period D ( i , k ) and the corresponding consumption j N i c i j P ( j , k ) , which is consumed by downstream production. Among them, N i is the downstream material set of material state i ; c i j is the conversion coefficient between material state i and material state j .
3. Demand constraints:
T d ( i , k ) = T d ( i , k 1 ) + O r d i k D ( i , k ) ;   i G p ,   k > 1 ,   k K .  
D ( i , k ) O r d i k + T d ( i , k 1 ) + t d i , k ; i G p , k > 1 , k K .
Constraint (4) represents that the backlog order quantity T d ( i , k ) of product i in the period k is equal to the backlog order quantity in the previous period T d ( i , k 1 ) plus the order quantity in the current period O r d i k , minus the final delivery D ( i , k )   . In Constraint (5), t d i , k is a parameter in the model and, as a backlog order quantity generated by the scheduling layer, feedback to the next planning period, where, the delivery quantity D ( i , k ) should be less than or equal to the sum of the order quantity of this period and backlog order quantity of the previous period, as well as the backlog order quantity from the scheduling layer in the previous period.
4. Production Capacity and Productivity Fluctuation Constraints:
P i min P ( i , k ) P i max   ;   i G i p   ,   k K .
P ( i , k ) F l u ( i , k ) P ( i , k 1 ) P ( i , k ) + F l u ( i , k ) ;   i G i p   ,   k > 1 ,   k K .
Constraint (6) represents the production capacity constraint of materials. Material status   i produced in the period k should be between maximum production capacity P i max and the minimum production capacity P i min . In order to avoid excessive fluctuation of productivity in the adjacent periods, auxiliary variable F l u ( i , k ) is introduced into Constraint (7) and corresponding penalty weight β is introduced into the objective function, so as to reduce the loss caused by productivity fluctuation as much as possible.
5. Inventory constraints:
S i min S ( i , k ) S i max   ;   i G i p   ,   k K .
S i L l I n v ( i , k ) S ( i , k ) S i U l + I n v ( i , k ) ; i G i p ,   k K .
Constraint (8) represents the inventory of material status   i at the end of the period k . S ( i , k ) should be between the inventory upper limit S i max and the inventory lower limit S i min . In order to avoid the situation that the inventory storage level is too high or too low, auxiliary variable I n v ( i , k ) is introduced into Constraint (9), and penalty weight θ is introduced into the corresponding objective function. The reference upper limit S i U l and lower limit S i L l of the inventory storage are set, so as to reduce the impact of the inventory storage level on the production processes.

2.1.2. Chance Constrained Uncertain Model

The uncertain factors in production planning have attracted more and more attention from scholars and researchers, as well as many methods and theories have been proposed to deal with the uncertain factors in production planning, including stochastic programming [26], fuzzy programming [27], robust optimization [28,29], recoverable robustness [30], sensitivity analysis [31] etc. At present, many research results have been obtained by using stochastic programming to deal with the problem of demand uncertainty [32,33]. Chance constrained programming is an important branch of stochastic programming. It allows decisions, to some extent, to not meet the constraints, but these decisions should make the probability of the constraints set up no less than a certain confidence level. The chance constrained method can more truly and effectively respond to demand uncertainty in the production planning.
The deterministic planning model is replaced by the following uncertainty planning model:
Objective function:(1)
Constraints: (2), (3), (6)–(9)
Constraint (10), deterministic constraint (4) is replaced by uncertainty constraint (10)
Constraint (11), deterministic constraint (5) is replaced by uncertainty constraint (11).
T d ( i , k ) = T d ( i , k 1 ) + O ˜ r d i k D ( i , k ) ;   i G p ,   k > 1 ,   k K .  
D ( i , k ) O ˜ r d i k + T d ( i , k 1 ) + t d i , k ; i G p , k > 1 , k K .
In Constraints (10) and (11), the parameter order demand quantity O ˜ r d i k is a random variable, assuming that the parameter obeys normal distribution, namely N ( E { O ˜ r d i k } , v a r { O ˜ r d i k } ) . Under a certain confidence level, Constraints (10) and (11) can be rewritten into Constraints (12) and (13).
P r { O ˜ r d i k D ( i , k ) + T d ( i , k ) T d ( i , k 1 ) } λ ; i G p .
P r { O ˜ r d i k D ( i , k ) T d ( i , k 1 ) t d i , k 1 } β ; i G p   k K .
By using the knowledge of Theorem 1 and Theorem 2 the probability density function can be used to further transform Constraints (12) and (13) into Constraints (14) and (15), and finally transform the model with demand uncertainty into a deterministic model with a confidence level. This transformation of the chance constraint has been solved in detail in many related studies [34,35], and it will not be repeated here in this paper.
[ D ( i , k ) + T d ( i , k ) T d ( i , k 1 ) ] E { O r d i k } v a r { O r d i k } Φ 1 ( i λ ) ; i G p .
[ D ( i , k ) T d ( i , k 1 ) t d i , k 1 ] E { O r d i k } v a r { O r d i k } Φ 1 ( i β ) ; i G p   k K .
Theorem 1.
For the random linear programming model:
m i n   Z = j = 1 n h j x j   ; s . t . { j = 1 n g i j x j r i ( ξ )   ( i = 1 , 2 , , m ) , x j 0   ( j = 1 , 2 , , n ) .
r i ( ξ )   ( i = 1 , 2 , , m ) is a single random variable, and the deterministic programming model with the formula at the confidence level ν is as follows:
m i n   Z = j = 1 n h j x j   ; s . t . { j = 1 n g i j x j K ν i   ( i = 1 , 2 , , m ) , x j 0   ( j = 1 , 2 , , n ) .
In the formula, K ν i =   s u p { K i | K i = F i 1 ( 1 ν ) } , F i ( ) is r i ( ξ )   distribution function.
Theorem 2.
For the random linear programming model:
m i n   Z = j = 1 n h j x j   ; s . t . { j = 1 n g i j x j r i ( ξ )   ( i = 1 , 2 , , m ) , x j 0   ( j = 1 , 2 , , n ) .
r i ( ξ )   ( i = 1 , 2 , , m ) is a single random variable, and the deterministic programming model with the formula at the confidence level ν is as follows:
m i n   Z = j = 1 n h j x j   ; s . t . { j = 1 n g i j x j K ν i   ( i = 1 , 2 , , m ) , x j 0   ( j = 1 , 2 , , n ) .
In the formula, K ν i =   i n f { K i | K i = F i 1 ( ν ) } , F i ( ) is r i ( ξ )   distribution function.

2.2. Scheduling Layer Model

The continuous-time modeling method was adopted to establish the scheduling layer model. The optimization objective of the scheduling layer is to complete the production tasks assigned from the planning layer as much as possible. Therefore, the minimization of the value of the difference between the task amount assigned by the planning layer and the task amount completed by the scheduling layer was taken as the objective function. Resource balance constraint, capacity constraint, demand constraint, time constraint, sequence constraint, utility supply and utility time constraint were also considered. In scheduling layer, fuzzy theory is utilized to describe the uncertainty of the utility. Finally, the uncertainty model of utility is transformed into a deterministic model with membership degree by corresponding mathematical method.

2.2.1. Deterministic Model

1. Objective function:
m i n D e f r , k n N d ( r , n )   ;   k K   ,   r R s .
The objective function is the minimization of the value of the difference between the task quantity assigned by the planning layer and the task quantity completed by the scheduling layer. The first part represents the reference task quantity assigned at the end of the planning period, and the second part represents the total final delivery of products completed by the scheduling layer within the planning period.
2. Resource balance constraints:
E ( r , n ) = E ( r , n 1 ) + i I r [ μ r i p w ( i , n 1 ) + ρ r i p b ( i , n 1 ) ] + i I r [ μ r i c w ( i , n ) + ρ r i c b ( i , n ) ] ;   r R   ,   n > 1 ,   n N .
E ( r , n ) = E i n i + i I r [ μ r i c w ( i , n ) + ρ r i c b ( i , n ) ] ;   r R , n = 1 .
Constraints (21) and (22) respectively represent resource balance constraints of subsequent event points and the initial event point. The amount of available resources at event n E ( r , n ) is equal to the amount of available resources at the previous event point E ( r , n 1 ) plus the production amount of the previous event point, and the consumption of the current point. Moreover, in Constraints (21) and (22), μ r i c 0 , ρ r i c 0 .
3. Capacity constraints:
E r min E ( r , n ) E r max   ;   r R ,   n N .
w ( i , n ) B i min b ( i , n ) w ( i , n ) B i max ;   n N , i I r .
Constraint (23) represents that the availability of the resource status r at the event point n is limited between the resource upper limit E r max and the resource lower limit E r min . For device resources, this constraint limits the appropriate number of device units; for material resources, this constraint limits the storage of material resources. In Constraint (24), B i max is the upper limit of material handling amount and B i min is the lower limit of material handling amount. When w ( i , n ) = 1 , material handing amount b ( i , n ) of task i should be between the upper limit of material handling amount B i max and the lower limit of material handling amount B i min ; otherwise, when w ( i , n ) = 0 , this constraint forces b ( i , n ) = 0 .
4. Demand constraints:
n N d ( r , n ) D e f r , k   ;   k K   ,   r R s .
t d ( r , k ) = D e f r , k n N d ( r , n )   ;   k K   ,   r R s .
Constraint (25) represents that the reference task quantity D e f r , k assigned by the planning layer within the specified period   k shall be greater than or equal to the total product delivery quantity. In Constraint (26), the backlog order quantity of product t d ( r , k ) of the scheduling layer should be equal to the difference between the reference task quantity D e f r , k of the planning layer and the final delivery quantity of the product, the backlog order quantity t d ( r , k ) in the obtained optimization results of the scheduling layer is fed back to the next period of planning model as t d i , k + 1 in Constraint (5). When we first solve the planning layer, we treat t d i , k as 0 and get the first period of optimization result of the planning layer. By inputting this result into the scheduling layer, the optimization result of scheduling layer, optimized for the first time, will be obtained. The backlog order quantity t d ( r , k ) in the obtained optimization results of the scheduling layer is fed back to the next period of planning model as t d i , k + 1 . The final optimization results of planning layer and scheduling layer are obtained by iterative solution.
t d ( r , k ) = t d i , k + 1 .
5. Time constraints:
T s ( i , n ) H   ;   i I r   ,   n N .
T f ( i , n ) = T s ( i , n ) + α i w ( i , n ) + β i b ( i , n ) ; i I r   ,   n N .
Constraint (28) represents that the start time T s ( i , n ) of tasks i at the event point n should be less than or equal to the scheduling horizon H . In Constraint (29), α i w ( i , n ) and β i b ( i , n ) respectively represent the constant term and variable term of task duration time, namely the end time T f ( i , n ) of tasks i at the event point n should be equal to the duration time of the task plus the start time of the task.
6. Sequence constraints:
a. Same task in the same unit:
T s ( i , n + 1 ) T f ( i , n )   ;   i I r   ,   n < N   ,   n N .
Constraint (30) represents that the start time of tasks i at the event point n + 1 should be greater than or equal to the end time of tasks i executing on the same device at the event point n .
b. Different tasks in the same unit:
T s ( i , n + 1 ) T f ( i , n )   ;   i I r ,   i I r ,   i i ,   r R J ,   n < N   ,   n N .
Constraint (31) represents that the start time of tasks i at the event point n + 1 should be greater than or equal to the end time of different tasks i executing on the same device at the event point n .
c. Different tasks in different units:
T s ( i , n + 1 ) T f ( i , n ) S l a ( 1 w ( i , n ) )   ;   i I r   ,   i I r   ,   i i   ,   r R S   ,   n < N   ,   n N .
In Constraint (32), if w ( i , n ) = 1 , the start time of tasks i at the event point n + 1 should be greater than or equal to the end time of different tasks   i on different devices at the event point n ; If w ( i , n ) = 0 , Constraint (32) is relaxed, where S l a is a positive big number.
7. Utility supply and utility time constraints:
i I u ( λ i u w ( i , n ) ) + i I u ( φ i u b ( i , n ) ) U u max ; u U ,   n N .
T u s ( u , n ) H   ;   u U   ,   n N .
T u s ( u , n + 1 ) T u s ( u , n )   ;   u U ,   n < N   ,   n N .  
Constraint (33) represents that the maximum supply U u max of utility u should be greater than or equal to the total consumption of the same utility u at the same event point n . The consumption of utility is determined by constant term λ i u w ( i , n ) and variable term φ i u b ( i , n ) . Constraint (34) represents that the start time consuming utility u should be less than or equal to the scheduling time horizon H . Constraint (35) represents that the start time of the utility u at the event point   n + 1 should be greater than or equal to the start time of the utility u at the event point n .

2.2.2. Fuzzy Uncertain Model

It is an effective method to describe the uncertainty of utility with fuzzy theory, because it is much easier to determine the membership function of a fuzzy number than the distribution function of a random variable. Different from the order demand, utilities are greatly influenced by the energy system, and the energy system is disturbed by many factors, so it is difficult to determine the specific distribution function to express the uncertainty of utilities. In this chapter, the fuzzy theory [36] is used to describe the uncertainty of utilities.
The deterministic scheduling model is replaced by the following fuzzy scheduling model:
Objective function: (20)
Constraints: (21)–(26), (28)–(32), (34), (35)
Constraint (36), deterministic constraint (33) is replaced by fuzzy constraint (36).
i I u ( λ i u w ( i , n ) ) + i I u ( φ i u b ( i , n ) ) U ^ u max ; u U ,   n N .
The utility supply U ^ u max in Constraint (36) is an uncertain parameter. It is assumed that the triangular membership function of U ^ u max is T r i ( u ) , which is expressed as follows:
T r i ( u ) = {   u u 1 u 2 u 1   ,   u 1 u u 2   ;   u 3 u u 3 u 2 ,   u 2 < u u 3   ;   0   ,   other .
Using the theory of fuzzy knowledge, Constraint (36) can be converted to Constraint (38). Finally, the uncertain model with utility constraints is transformed into a deterministic model with membership degree.
i I u ( λ i u w ( i , n ) ) + i I u ( φ i u b ( i , n ) ) w 1 U 1 , δ max + w 2 U 2 , δ max + w 3 U 3 , δ max ;   u U ,   n N .
w 1 , w 2 , w 3 are the weights, and w 1 + w 2 + w 3 = 1 ; the fuzzy parameters of High-pressure Steam (HS) and Cooling Water (CW) include the most pessimistic value, the most possible value, and the most optimistic value. The three prominent values, i.e., the most pessimistic, the most possible, and the most optimistic values for each fuzzy number are usually estimated by decision maker [37,38]. In this paper, the three prominent values are defined according to determining method of universe and the experience of engineering. The universe U of fuzzy number is defined as an equation according to Basyigit and Ulu [39], where D m i n is the minimum value in statistical data of fuzzy number, D m a x is the maximum value in statistical data of fuzzy number. D 1 and D 2 are defined as an equation and according to the experience of engineering. D m f is the value of which the frequency of occurrence is the highest in statistical data of fuzzy numbers.
U = [ D m i n D 1 , D m a x + D 2 ] ,
D 1 = 0.2 | D m i n D m f | ,
D 2 = 0.2 | D m a x D m f | .
In order to describe fuzzy numbers more generally, in this paper we define u 1 , u 2 and u 3 as the most pessimistic value, the most possible value, and the most optimistic value. The most pessimistic value u 1 is equal to D m i n D 1 (the left boundary of universe). The most possible value u 2 is equal to D m f . The most optimistic value u 3 is equal to D m a x + D 2 (the right boundary of universe).
u 1 = D m i n D 1 ,
u 2 = D m f ,
u 3 = D m a x + D 2 .
The triangular membership degrees of U 1 , δ max , U 2 , δ max , U 3 , δ max are the boundary points of U u max under the cut set δ , the corresponding calculation formula of the boundary point is as follows, and the relationship between them is shown in Figure 1.
U 1 , δ max = u 1 + δ ( u 2 u 1 ) ,
U 2 , δ max = u 2 ,
U 3 , δ max = u 3 δ ( u 3 u 2 ) .

3. Integrated Model and Solving Strategy

In this paper, the integrated model of production planning and scheduling considering demand uncertainty and utility uncertainty is solved based on the integrated optimization method and rolling horizon optimization strategy. The solving method of integrated optimization is to construct a bi-level integrated model of planning and scheduling to deal with the demand and utility uncertainties. As the uncertainty of demand and utility are different in time scales, it is necessary to deal with the uncertainty of demand and utility in production planning layer and scheduling layer respectively. The discrete-time model is established in the production planning layer, and chance constrained programming is used to describe demand uncertainty. The continuous-time model is established in the production scheduling layer, and fuzzy theory is used to describe utility uncertainty. When the planning layer completes the Kth optimization, the optimization result D ( i , k ) corresponding to the first planning period is equal to D e f r , k . Parameter D e f r , k is assigned as a task to the scheduling layer. The double-layer structure of production planning and scheduling is shown in Figure 2.
Rolling horizon optimization strategy is a general time decomposition method. The planning horizon of the discrete-time model in the planning layer is composed of m planning periods with equal time durations. The scheduling horizon of the continuous-time model in the scheduling layer is equal to the length of the planning period. The rolling horizon optimization strategy is used to iteratively solve the integrated model of production planning and scheduling. The specific steps of rolling horizon optimization strategy are as follows:
a. Solve the production planning model and assign the delivery quantity of the first planning period to the scheduling model.
b. Solve the scheduling model of the current period and judge whether the scheduling layer has completed the delivery quantity assigned by the planning layer. If it can be completed, roll the planning horizon and repeat step a; if it cannot be completed, the unfinished tasks will be passed to the next planning period as the backlog order quantity, then roll the planning horizon and repeat the step a.
c. The calculation is terminated after obtaining the planning and scheduling schemes for all periods.
The flow chart of iterative solving strategy of integrated model is shown in Figure 3, and the schematic diagram of rolling horizon optimization strategy is shown in Figure 4.

4. Case Study

The classical example [40] was used to verify the feasibility and effectiveness of the two-layer integrated model of planning and scheduling in dealing with the uncertain problems of demand and utility based on the RTN representation. The example uses RTN representation to describe the production process, as shown in Figure 5. The three raw materials S1, S2 and S3 produce two products S8 and S9 through a series of endothermic and exothermic reactions, in which the endothermic reaction requires high-pressure steam (HS) and the exothermic reaction requires condensed water (CW). Due to the limitations of the energy system, the average supply of HS and CW is 64 and 69 kg/min, respectively. Heating can be performed on heater J1. Reaction 1 can be performed on reactors J2, J3. Reaction 2 can be performed on reactors J2, J3. Reaction 3 can be performed on reactors J2, J3. Separation can be performed on separator J4. The study involved a total of eight different production processes. Unlimited storage policy is adopted for raw materials and products, and limited storage policy is adopted for intermediate materials [41,42]. The relevant data of the calculation example are shown in Table 1, Table 2 and Table 3.

4.1. The Influence of the Planning Layer Confidence Factor on the Optimization Results

The stochastic programming with chance constraints is introduced into the production planning layer to describe the uncertainty of demand. Through a series of mathematical transformation, the uncertain model with chance constraints can be transformed into a deterministic model with specified confidence level. Different confidence levels have a very important impact on the optimization results. The objective functions under different confidence factors are shown in Table 4. The value of the confidence factor in the part of integrated model simulation analysis is 0.9. The variation trend of maximum profit under different confidence levels is shown in Figure 6.

4.2. The Influence of Scheduling Layer Fuzzy Weight on the Optimization Results

In the production scheduling layer, fuzzy theory is introduced to describe the uncertainty of utilities. It is assumed that triangular membership function is adopted to express the uncertainty of utilities. Using the knowledge of fuzzy theory, the uncertainty model with utility is transformed into a certainty model with fuzzy weights. The fuzzy parameter values of HS is u 1 = 62   ,   u 2 = 64   ,   u 3 = 68 ; and the fuzzy parameter values of CW is u 1 = 67   ,   u 2 = 69   ,   u 3 = 73 . A visualized sensitivity analysis of fuzzy weights is illustrated in Figure 7. The optimization results of scheduling layer under different weights are shown in Table 5. It can be seen from Figure 7 that when the weight of the maximum possible value increases from 0.1 to 0.8, the optimal solution of the uncertain model is closer to the optimal solution of the deterministic model. This is because the weights w 1 , w 2 , and w 3 represent the degree of inclination of the model to the most optimistic value, the maximum possible value, and the most pessimistic value in the membership function. In other words, if the weight of the boundary is large, the influence of the boundary on the optimization result will be greater. For the fuzzy weights, w 1 , w 2 , w 3 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 [38], and considering several relevant works [43,44], we set these parameters to: w 1 = 0.1   ,   w 2 = 0.5   ,   w 3 = 0.4 in this paper.

4.3. Integrated Model Simulation Analysis

The order data of two products in nine planning periods and the corresponding average supply of two utilities are given. Based on the RTN representation, integrated optimization solving method and rolling horizon optimization strategy were used to iteratively solve the integrated model to complete the planning and scheduling tasks of the first five planning periods. The software Lingo 11.0 was used for simulation calculation on a 2.20 GHz computer, where the planning horizon is 5 days and the planning period and scheduling horizon are both 1 day. The number of event points is 12. Figure 8 and Figure 9 respectively show the gantt chart and utility consumption level of the deterministic method and the uncertainty method. From Figure 8 and Figure 9, we can intuitively get the respective task number and the corresponding specific utility consumption under the deterministic method and the uncertain method. Under the deterministic method, the solution time of the rolling horizon algorithm is 28 s. Under the uncertain method, the solution time of the rolling horizon algorithm is 34 s.
In the comparative analysis between the deterministic method and the uncertain method shown in Table 6 and Table 7, S 8 p l a n and S 9 p l a n respectively represent the planned production of S 8 and S 9 ; S 8 s c h e d u l e and S 9 s c h e d u l e represent the scheduled production of S 8 and S 9 respectively. The objective function of the planning layer of the deterministic method is $99,688.52, and the output of products S 8 and S 9 in the planning layer are, respectively, 1330 and 1550 kg. In the uncertain method with the confidence factor of 0.9, the objective function of the planning layer is $97,285.11, and the output of products S 8 and S 9 are respectively 1296.998 and 1514.416 kg. The objective function in the scheduling layer of the deterministic method is 128.9 kg, and the output of products S 8 and S 9 are, respectively, 1201.1 and 1550 kg. The objective function in the scheduling layer of the uncertain method is 87.277 kg, and the output of products S 8 and S 9 are, respectively, 1209.721 and 1514.416 kg. The final profit of the enterprise under the deterministic method is $86,792.96, and the final profit of the enterprise under the uncertain method is $88,557.392. The comparison of the final profit shows that the uncertain method reflects the actual production process more realistically due to considering the uncertainty of the demand and the utility, which brings greater profits to the enterprise. Table 8 and Table 9 respectively represent the consumption amount of utility in each period under the deterministic method and uncertain method, where u s e i   ( i = 1 5 ) represents the consumption amount of utility in the planning period i and u s e a l l represents the total usage amount of utility. The total amounts of utility HS used under the deterministic method and the uncertain method are 1804.94419 and 1735.6687 kg, respectively. The total amounts of utility CW used under the deterministic method and the uncertain method are 1685.06007 and 1732.3723 kg, respectively.
Backlog order quantity of 87.277 kg under the uncertain method is significantly less than that of 128.9 kg under the deterministic method. These comparative results show that introducing the uncertain programming method into the planning and scheduling layers, respectively, can greatly reduce the generation of backlog order quantity, and the scheduling layer can more effectively complete the tasks assigned by the planning layer. Based on the RTN representation, the uncertain method adopted in the scheduling layer can effectively reduce the impact of uncertainty factors on the optimization results, and can clearly count the consumption amount of utilities in each period. Therefore, the proposed integrated model can more realistically reflect the actual production situation, save energy, and reduce production costs, so that the final profit of the enterprise obtained by the uncertain method is $88,557.392, which is obviously better than the final profit of the enterprise under the deterministic method of $86,792.96.
Considering the uncertainty of demand and utility, the proposed method is slightly more complex than the deterministic model in modeling and solving, but it belongs to the category of pre-scheduling, and the solving time is 34 s, which can fully meet the real-time requirements of actual production. In addition, this method has more application value for chemical enterprises whose demand is affected by seasons and which need utilities, especially for small businesses. In order to reduce the cost, small businesses will be gathered and purchase utilities from other companies to support their own production. Compared with large enterprises which have their own energy systems, the utilities supply for small businesses experiences disturbance more easily. Consequently, for this kind of small business, considering the uncertain utility problem for effective enterprise production process is of great application value and practical significance.

5. Conclusions

This paper introduces and analyzes uncertain factors involved in production planning and scheduling. According to the property of time scale, these uncertain factors are dealt with distinctively in the planning and scheduling layers. The planning layer and scheduling layer models are formulated by discrete-time and continuous-time modeling methods, respectively. Chance constrained programming and fuzzy theory are introduced into planning layer model and scheduling layer model, respectively, to describe the uncertainties of demand and utility. Rolling horizon optimization strategy is utilized to solve the bi-level integrated model iteratively. Through comparison, simulation results show that the proposed model and algorithm are feasible and effective, can calculate the consumption of utility in every period, decrease the effects of uncertain factors on optimization results, and more accurately reflect the actual production process.

Author Contributions

T.Z. conceived, formulated and verified the scheduling model. Y.W., X.J. and S.L. analyzed the results, discussed the manuscript and corrected the manuscript. All authors have read and approved the final manuscript.

Funding

This work was supported by Scientific Research Fund of Liaoning Provincial Education Department (No. L2017LQN030, No. L2017LQN032), supported by Talent Scientific Research Fund of LSHU (No. 2016XJJ-101 and No. 2016XJJ-102), and supported by Youth Innovation Talent Program by Department of Education of Guangdong Province, China (601821K42050).

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

(1) Production planning layer
a. Indices:
iMaterial status
kPlanning period
b. Sets:
G a l l Raw material, intermediate product, final product material status set
G i p Intermediate product and final product material status set
G p Final product material status set
N i The downstream material status set of material status i
c. Parameters:
P r i c e i The market price of product i , $ k g 1
α The penalty weight of backlog order item, which is 100 in this paper
β The penalty weight of productivity fluctuation item, which is 0.0001 in this paper
θ The penalty weight of inventory limit item, which is 0.0001 in this paper
S i i n i Initial inventory of material status i
c i j Conversion coefficient between material status i and material status j
O r d i k Order quantity for product i in period k
p i min Minimum production capacity of material status i
p i max Maximum production capacity of material status i
s i min Inventory lower limit of material status i
s i max Inventory upper limit of material status i
s i L l Inventory reference lower limit of material status i
s i U l Inventory reference upper limit of material status i
t d i , k 1 Backlog order quantity from material status i of the scheduling layer
P r { } Probability computing operator
λ Confidence factor
d. Variables:
S ( i , k ) Inventory of material status i in planning period k
P ( i , k ) Amount of material status i produced in planning period k
D ( i , k ) Delivery of material status i in planning period k
T d ( i , k ) Backlog order quantity of material status i in planning period k
F l u ( i , k ) Auxiliary variable of productivity fluctuation
I n v ( i , k ) Auxiliary variable of inventory limit
(2) Production scheduling layer
a. Indices:
i , i Task
r Resource
n Event point
u Utility
b. Sets:
I r A set of tasks that consume material resources
I u A set of tasks that consume utilities
R   , R S   , R J Resource set, Material resource set, Device resource set
N Event point set
H Scheduling horizon
c. Parameters:
D e f r , k The amount of reference tasks assigned at the end of the planning period k
E i n i Initial resource quantity
u s e i Consumption amount of utility in the planning period i ( i = 1 5 )
u s e a l l Total usage amount of utility
μ r i p   , μ r i c Production and consumption coefficients of device resources involved in task i ; μ r i p 0   μ r i c 0
ρ r i p   , ρ r i c Production and consumption coefficients of material resources involved in task i ; ρ r i p 0   ρ r i c 0
E r min The lower limit of resource r
E r max The upper limit of resource r
B i min Lower limit of material handing amount of task i
B i max Upper limit of material handing amount of task i
S l a Large positive number in relaxation constraints.
w 1   , w 2   , w 3 Weights
U 1 δ max   , U 2 δ max   , U 3 δ max The boundary points of the triangular membership function under the cut set δ
T r i ( u ) Triangle membership function of U u max
d. Variables:
d ( r , n ) The delivery of material resources r at the event point n
E ( r , n ) The amount of available resources for resource r at the event point n
w ( i , n ) Binary variable for task i active at event n
b ( i , n ) The amount of resources processed by task i at event point n
t d ( r , k ) The backlog order quantity of product in the period k
T s ( i , n ) Start time of task i at event point n
T f ( i , n ) End time of task i at event point n
T u s ( u , n ) The start time of the utility u consumed at event point n
α i w ( i , n ) Constant term of duration time
β i b ( i , n ) Variable term of duration time
λ i u w ( i , n ) Constant term of utility consumption
φ i u b ( i , n ) Variable term of utility consumption

References

  1. Chu, Y.; You, F.; Wassick, J.M.; Agarwal, A. Integrated planning and scheduling under production uncertainties: Bi-level model formulation and hybrid solution method. Comput. Chem. Eng. 2015, 72, 255–272. [Google Scholar] [CrossRef]
  2. An, Y.W.; Yan, H.S. Solution Strategy of Integrated Optimization of Production Planning and Scheduling in a Flexible Job-shop. Acta Autom. Sin. 2013, 39, 1476–1491. [Google Scholar] [CrossRef]
  3. Li, W.; Ong, S.K.; Nee, A.Y.C. A Unified Model-based Integration of Process Planning and Scheduling. In Process Planning & Scheduling for Distributed Manufacturing; Springer: London, UK, 2007; pp. 295–309. [Google Scholar]
  4. Shao, X.; Li, X.; Liang, G.; Zhang, C. Integration of process planning and scheduling—A modified genetic algorithm-based approach. Comput. Oper. Res. 2009, 36, 2082–2096. [Google Scholar] [CrossRef]
  5. Lasserre, J.B. An Integrated Model for Job-Shop Planning and Scheduling. Manag. Sci. 1992, 38, 1201–1211. [Google Scholar] [CrossRef]
  6. Xiong, F.; Yan, H.S. Integrated production planning and scheduling of multi-stage workshop based on alternant iterative genetic algorithm. J. Southeast Univ. (Nat. Sci. Ed.) 2012, 42, 183–187. [Google Scholar]
  7. Li, Z.K.; Ierapetritou, M.G. Production planning and scheduling integration through augmented Lagrangian optimization. Comput. Chem. Eng. 2010, 34, 996–1006. [Google Scholar] [CrossRef]
  8. Wang, Y.; Su, H.Y.; Shao, H.S.; Lu, S.; Xie, L. Integration of Production Planning and Scheduling under Demand and Utility Uncertainties. J. Zhejiang Univ. (Eng. Sci.) 2017, 51, 57–67. [Google Scholar]
  9. Vogel, T.; Almada-Lobo, B.; Almeder, C. Integrated versus hierarchical approach to aggregate production planning and master production scheduling. OR Spectr. 2017, 39, 193–229. [Google Scholar] [CrossRef]
  10. Hassani, Z.I.M.; El Barkany, A.; Jabri, A.; El Abbassi, I.; Darcherif, A.M. New Approach to Integrate Planning and Scheduling of Production System: Heuristic Resolution. Int. J. Eng. Res. Afr. 2018, 39, 156–169. [Google Scholar] [CrossRef]
  11. Pantelides, C.C. Unified Frameworks for the Optimal Process Planning and Scheduling. In Proceedings of the Second Conference on the Foundations of Computer Aided Process Operations, Crested Butte, CO, USA, 18–23 July 1993; CACHE Publications: New York, NY, USA, 1994; pp. 253–274. [Google Scholar]
  12. Zhang, X.; Sargent, R.W.H. The optimal operation of mixed production facilities—A general formulation and some approaches for the solution. Comput. Chem. Eng. 1996, 20, 897–904. [Google Scholar] [CrossRef]
  13. Castro, P.; Barbosa-Póvoa, A.P.F.D.; Matos, H. An Improved RTN Continuous-Time Formulation for the Short-term Scheduling of Multipurpose Batch Plants. Ind. Eng. Chem. Res. 2001, 40, 2059–2068. [Google Scholar] [CrossRef]
  14. Chen, G. Modeling and Optimization for Short-term Scheduling of Multipurpose Batch Plants. Chin. J. Chem. Eng. 2014, 22, 682–689. [Google Scholar] [CrossRef]
  15. Schilling, G.; Pantelides, C.C. A simple continuous-time process scheduling formulation and a novel solution algorithm. Comput. Chem. Eng. 1996, 20, S1221–S1226. [Google Scholar] [CrossRef]
  16. Dimitriadis, A.D.; Shah, N.; Pantelides, C.C. RTN-based rolling horizon algorithms for medium term scheduling of multipurpose plants. Comput. Chem. Eng. 1997, 21, S1061–S1066. [Google Scholar] [CrossRef]
  17. Chao, S.; You, F. Distributionally robust optimization for planning and scheduling under uncertainty. Comput. Chem. Eng. 2017, 110, 53–68. [Google Scholar]
  18. Curcio, E.; Amorim, P.; Qi, Z.; Almada-Lobo, B. Adaptation and approximate strategies for solving the lot-sizing and scheduling problem under multistage demand uncertainty. Int. J. Prod. Econ. 2018, 202, 81–96. [Google Scholar] [CrossRef]
  19. Li, P.; Wendt, M.; Wozny, G. Optimal production planning under uncertain market conditions. Comput. Aided Chem. Eng. 2003, 15, 511–516. [Google Scholar]
  20. Wang, Y.; Jin, X.; Xie, L.; Zhang, Y.; Lu, S. Uncertain Production Scheduling Based on Fuzzy Theory Considering Utility and Production Rate. Information 2017, 8, 158. [Google Scholar] [CrossRef]
  21. Chen, Q.X.; Jiao, B.; Yan, S.B. Flow Shop Production Scheduling under Uncertainty within Infinite Intermediate Storage. Adv. Mater. Res. 2011, 204–210, 777–783. [Google Scholar] [CrossRef]
  22. Lamba, N.; Karimi, I.A. Scheduling Parallel Production Lines with Resource Constraints. 1. Model Formulation. Ind. Eng. Chem. Res. 2002, 41, 790–800. [Google Scholar] [CrossRef]
  23. Méndez, C.A.; Cerdá, J. An MILP framework for short-term scheduling of single-stage batch plants with limited discrete resources. Comput. Aided Chem. Eng. 2002, 10, 721–726. [Google Scholar]
  24. Tian, Y. Chemical Production Planning and Scheduling Integration under Uncertainty. J. Chem. Ind. 2014, 65, 3552–3558. [Google Scholar]
  25. Chand, S.; Traub, R.; Uzsoy, R. Rolling horizon procedures for the single machine deterministic total completion time scheduling problem with release dates. Ann. Oper. Res. 1997, 70, 115–125. [Google Scholar] [CrossRef]
  26. Bonfill, A.; Bagajewicz, M.; Espuña, A.; Puigjaner, L. Risk management in the scheduling of batch plants under uncertain market demand. Ind. Eng. Chem. Res. 2004, 43, 741–750. [Google Scholar] [CrossRef]
  27. Chanas, S.; Kasperski, A. Minimizing maximum lateness in a single machine scheduling problem with fuzzy processing times and fuzzy due dates. Eng. Appl. Artif. Intell. 2001, 14, 377–386. [Google Scholar] [CrossRef]
  28. Ning, C.; You, F. A data-driven multistage adaptive robust optimization framework for planning and scheduling under uncertainty. AIChE J. 2017, 63, 4343–4369. [Google Scholar] [CrossRef]
  29. Leung, S.C.H.; Tsang, S.O.S.; Ng, W.L.; Wu, Y. A robust optimization model for multi-site production planning problem in an uncertain environment. Eur. J. Oper. Res. 2007, 181, 224–238. [Google Scholar] [CrossRef]
  30. Iris, Ç.; Lam, J.S.L. Recoverable robustness in weekly berth and quay crane planning. Transp. Res. Part B Methodol. 2019, 122, 365–389. [Google Scholar] [CrossRef]
  31. Cheng, H.N.; Xiang, S.G.; Yang, X.; Han, F.Y. Supply chain planning model under uncertainty. Comput. Appl. Chem. 2004, 21, 97–102. [Google Scholar]
  32. Safaei, A.S.; Farsad, S.; Paydar, M.M. Emergency logistics planning under supply risk and demand uncertainty. Oper. Res. 2018. Available online: https://doi.org/10.1007/s12351-018-0376-3 (accessed on 30 January 2018).
  33. Englberger, J.; Herrmann, F.; Manitz, M. Two-stage stochastic master production scheduling under demand uncertainty in a rolling planning environment. Int. J. Prod. Res. 2016, 54, 6192–6215. [Google Scholar] [CrossRef]
  34. Fajie Wei, R.Z. Activity-Oriented Production Theory and Cost Theory of Enterprises; Beijing University of Aeronautics and Astronautics Press: Beijing, China, 2007. [Google Scholar]
  35. Wang, Y. Research on Production Planning and Scheduling under Uncertainty. Ph.D. Thesis, Zhejiang University, Hangzhou, China, 2016. [Google Scholar]
  36. Iris, C.; Cevikcan, E. A fuzzy linear programming approach for aggregate production planning. In Supply Chain Management under Fuzziness; Springer: Berlin/Heidelberg, Germany, 2014; pp. 355–374. [Google Scholar]
  37. Torabi, S.A.; Ebadian, M.; Tanha, R. Fuzzy hierarchical production planning (with a case study). Fuzzy Sets Syst. 2010, 161, 1511–1529. [Google Scholar] [CrossRef]
  38. Lai, Y.J.; Hwang, C.L. A new approach to some possibilistic linear programming problems. Fuzzy Sets Syst. 1992, 49, 121–133. [Google Scholar] [CrossRef]
  39. Basyigit, A.; Ulu, C.; Guzelkaya, M. A New Fuzzy Time Series Model Using Triangular and Trapezoidal Membership Functions. In Proceedings of the International Work-Conference on Time Series, Granada, Spain, 25–27 June 2014; pp. 25–27. [Google Scholar]
  40. Li, Z.; Ierapetritou, M.G. Rolling horizon based planning and scheduling integration with production capacity consideration. Chem. Eng. Sci. 2010, 65, 5887–5900. [Google Scholar] [CrossRef]
  41. And, M.G.I.; Floudas, C.A. Effective Continuous-Time Formulation for Short-Term Scheduling. 1. Multipurpose Batch Processes. Ind. Eng. Chem. Res. 1998, 37, 4341–4359. [Google Scholar]
  42. Mokhtari, H.; Abadi, I.N.K.; Zegordi, S.H. Production capacity planning and scheduling in a no-wait environment with controllable processing times: An integrated modeling approach. Expert Syst. Appl. 2011, 38, 12630–12642. [Google Scholar] [CrossRef]
  43. Garciasabater, J.P. Capacity and material requirement planning modelling by comparing deterministic and fuzzy models. Int. J. Prod. Res. 2008, 46, 5589–5606. [Google Scholar] [Green Version]
  44. Christou, I.; Lagodimos, A.; Lycopoulou, D. Hierarchical production planning for multi-product lines in the beverage industry. Prod. Plan. Control 2007, 18, 367–376. [Google Scholar] [CrossRef]
Figure 1. The triangular membership function of parameter U u max .
Figure 1. The triangular membership function of parameter U u max .
Algorithms 12 00120 g001
Figure 2. The double-layer structure of production planning and scheduling.
Figure 2. The double-layer structure of production planning and scheduling.
Algorithms 12 00120 g002
Figure 3. Flow chart of integrated model iterative solving strategy.
Figure 3. Flow chart of integrated model iterative solving strategy.
Algorithms 12 00120 g003
Figure 4. Rolling horizon optimization strategy.
Figure 4. Rolling horizon optimization strategy.
Algorithms 12 00120 g004
Figure 5. Resource-Task Network diagram.
Figure 5. Resource-Task Network diagram.
Algorithms 12 00120 g005
Figure 6. The variation trend of maximum profit under different confidence levels.
Figure 6. The variation trend of maximum profit under different confidence levels.
Algorithms 12 00120 g006
Figure 7. Scheduling layer optimization results under different weights.
Figure 7. Scheduling layer optimization results under different weights.
Algorithms 12 00120 g007
Figure 8. Gantt chart and utility consumption of the deterministic method.
Figure 8. Gantt chart and utility consumption of the deterministic method.
Algorithms 12 00120 g008
Figure 9. Gantt chart and utility consumption of the uncertain method.
Figure 9. Gantt chart and utility consumption of the uncertain method.
Algorithms 12 00120 g009
Table 1. Task processing time coefficient and capacity data.
Table 1. Task processing time coefficient and capacity data.
Task (i)Unit (j) α i β i B i min B i max
Heating (i = 1)Heater0.6670.00667-100
Reaction 1 (i = 2)Reactor 11.3340.02664-50
Reaction 1 (i = 3)Reactor 21.3340.01665-80
Reaction 2 (i = 4)Reactor 11.3340.02664-50
Reaction 2 (i = 5)Reactor 21.3340.01665-80
Reaction 3 (i = 6)Reactor 10.6670.01332-50
Reaction 3 (i = 7)Reactor 20.6670.008325-80
Separation (i = 8)Separator1.33420.00666-200
Table 2. Material and device resource related data.
Table 2. Material and device resource related data.
ResourceCapacity (kg)Initial (kg)Price ($/kg)
S1ULAA0
S2ULAA0
S3ULAA0
S410000
S520000
S615000
S720000
S8UL040
S9UL030
Heater11-
Reactor111-
Reactor211-
Separator11-
UL—Unlimited; AA—available as and when required.
Table 3. Utility consumption coefficient.
Table 3. Utility consumption coefficient.
Task (i) λ i H S
( k g min 1 )
φ i H S
( k g min 1 )
λ i C W
( k g min 1 )
ϕ i C W
( k g min 1 )
i = 160.25--
i = 2--40.25
i = 350.25--
i = 4--40.3
i = 540.5--
i = 6--30.3
i = 740.2--
i = 8--60.35
Table 4. The objective function under different confidence factors.
Table 4. The objective function under different confidence factors.
Confidence FactorInverse FunctionObjective Function ($)
0.6−0.2599,206.36
0.65−0.3998,945.22
0.7−0.5298,702.73
0.75−0.6798,422.94
0.8−0.8498,105.84
0.85−1.0497,732.78
0.9−1.2897,285.11
0.95−1.6596,594.95
Table 5. Scheduling layer optimization result data under different weights.
Table 5. Scheduling layer optimization result data under different weights.
w 1 w 2 w 3 Objective Function (kg)
0.10.10.867.59
0.10.20.772.51
0.10.30.677.43
0.10.40.582.35
0.10.50.487.27
0.10.60.392.19
0.10.70.297.11
0.10.80.1102.03
Table 6. Comparative analysis of the planning layer.
Table 6. Comparative analysis of the planning layer.
Modeling MethodObjective Function ($) S 8 p l a n ( k g ) S 9 p l a n ( k g )
Deterministic method99,688.5213301550
Uncertain method97,285.111296.9981514.416
Table 7. Comparative analysis of the scheduling layer.
Table 7. Comparative analysis of the scheduling layer.
Modeling MethodObjective Function ($) S 8 p l a n ( k g ) S 9 p l a n ( k g )
Deterministic method128.91201.11550
Uncertain method87.2771209.7211514.416
Table 8. The consumption amount of each period of utility under the deterministic method.
Table 8. The consumption amount of each period of utility under the deterministic method.
u t i l i t y u s e 1 ( k g ) u s e 2 ( k g ) u s e 3 ( k g ) u s e 4 ( k g ) u s e 5 ( k g ) u s e a l l ( k g )
HS352.31419345.83433383.03704370.31418353.444451804.94419
CW337.26051339.90991357.33333322.22298328.333341685.06007
Table 9. The consumption amount of each period of utility under the uncertain method.
Table 9. The consumption amount of each period of utility under the uncertain method.
u t i l i t y u s e 1 ( k g ) u s e 2 ( k g ) u s e 3 ( k g ) u s e 4 ( k g ) u s e 5 ( k g ) u s e a l l ( k g )
HS339.07849371.35303340.01582343.52171341.697821735.66687
CW340.84803359.11669372.16621330.04896330.192411732.3723

Share and Cite

MDPI and ACS Style

Zhang, T.; Wang, Y.; Jin, X.; Lu, S. Integration of Production Planning and Scheduling Based on RTN Representation under Uncertainties. Algorithms 2019, 12, 120. https://doi.org/10.3390/a12060120

AMA Style

Zhang T, Wang Y, Jin X, Lu S. Integration of Production Planning and Scheduling Based on RTN Representation under Uncertainties. Algorithms. 2019; 12(6):120. https://doi.org/10.3390/a12060120

Chicago/Turabian Style

Zhang, Tao, Yue Wang, Xin Jin, and Shan Lu. 2019. "Integration of Production Planning and Scheduling Based on RTN Representation under Uncertainties" Algorithms 12, no. 6: 120. https://doi.org/10.3390/a12060120

APA Style

Zhang, T., Wang, Y., Jin, X., & Lu, S. (2019). Integration of Production Planning and Scheduling Based on RTN Representation under Uncertainties. Algorithms, 12(6), 120. https://doi.org/10.3390/a12060120

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop