Mathematical Model for Production Plan Optimization — A Case Study of Discrete Event Systems †

This paper proposes an optimal scheduling model under production and maintenance constraints for a real case of a discrete event system. The intent was to use the rich mathematical theory and algorithms of optimization in the study of this important class of systems. The current study detailed firstly a new approach for mapping a simulation event relationship graph into a mixed-integer program, with a flexible workshop real case. Several other potential applications of the mathematical model are examined, thanks to the model constraints flexibility characteristics, including a general case of a manufacturing system for optimal resource scheduling, an application on the case of hospital beds’ management. The model extension could be also interesting for other applications like museum systems or the case of big data in complex and social networks.


Introduction
For most companies, maintenance represents a very significant function within the overall production environment. The developments in automation, and the resulting complexity of the systems involved, have made the reliability of the machines even more important. This is especially true in the process industry, characterized by expensive, specialized equipment and stringent environmental considerations. Nowadays, with profit margins decreasing, there is need for a good maintenance planning and control, especially for discrete event systems (DES). The evolution of DES is defined by the instantaneous occurrence of events at possibly unknown and/or irregular time intervals (e.g., communication and transport systems, computer systems, and manufacturing systems). Today, these systems are so complex that intuitive methods are no longer enough to ensure efficient and operational supervision and scheduling. The theory of supervision of DES was developed by Wonham and Ramadge [1,2]. Thanks to formal languages and automata, these authors have introduced concepts like those used in the theory of automatic control of continuous systems. Thus, the synthesized supervisor ensures the desired legal behavior. Indeed, Petri nets (PNs) are shown to be more powerful than automata for modeling DES with synchronization and parallelism. In addition, the PNs are more advantageous for dealing with the control problems defined by specifications on the states of the system [3,4]. They represent a very suitable tool for studying manufacturing systems thanks to their modeling power and mathematical properties. In our previous work [5][6][7] we focused on the development of methods of control and production optimization manufacturing systems thanks to their modeling power and mathematical properties. In our previous work [5][6][7] we focused on the development of methods of control and production optimization of the computation time to calculate controllers that satisfy control specifications and constraints related to manufacturing systems.
The production and maintenance management using new strategies differs greatly from management policies in traditional manufacturing. These strategies are more deterministic and more flexible since they could consider different parameters impacting production and control activities' performances. Production scheduling is considered as a determinant functionality in the production management. Since production activities depend on several parameters, production scheduling may need to consider the manufacturing system status and operational conditions in the production layer.
In a more general context, decision support systems continue to be very popular in industrial engineering thanks to their ability to help decision makers managing their resources in function of several manufacturing constraints. However, manufacturing systems are complex since they are connecting many technical and human resources. These systems are often looking for more efficiency and, so, new strategies are developed, especially in production planning.
The importance of considering production management components linked to manufacturing execution systems (MES) has hugely increased during the last years. The Figure 1 [8] illustrates the MES growth between 2006 and 2016. In literature, many proposals deal with manufacturing system allocation and production planning. Several research works [9][10][11] focused on mathematical modeling to express the production scheduling problem in terms of variables and constraints. In [9], a set of algorithms to minimize the length of production planning was proposed. The proposed algorithms are used to generate all possible optimal schedules grouped into subsets. The authors propose computational algorithms in order to evaluate these performances and select the best ones in function of given performances' parameters. In [10], a hybrid policy to deal with the optimal production scheduling complexity due to various kinds of uncertain factors, such as breakdowns, maintenance actions, and production parameters, was developed. The authors used the genetic algorithm as a resolution tool to minimize the total completion time. In [11], the production planning computation complexity for various types of cost functions, maintenance costs, and production capacity limits was summarized. Many research works deal with several industrial constraints such as subcontracting, warranty, and quality in manufacturing and logistic fields [12]. Several research works [9][10][11] focused on mathematical modeling to express the production scheduling problem in terms of variables and constraints. In [9], a set of algorithms to minimize the length of production planning was proposed. The proposed algorithms are used to generate all possible optimal schedules grouped into subsets. The authors propose computational algorithms in order to evaluate these performances and select the best ones in function of given performances' parameters. In [10], a hybrid policy to deal with the optimal production scheduling complexity due to various kinds of uncertain factors, such as breakdowns, maintenance actions, and production parameters, was developed. The authors used the genetic algorithm as a resolution tool to minimize the total completion time. In [11], the production planning computation complexity for various types of cost functions, maintenance costs, and production capacity limits was summarized. Many research works deal with several industrial constraints such as subcontracting, warranty, and quality in manufacturing and logistic fields [12].
For more details about determination of production and maintenance control policies under different types of constraints, the reader may refer to the work of Rezg et al. [13].
An application example of our methodology such as hospital beds' management is presented in this paper to show the effectiveness of our mathematical model. Then, two illustrative examples are given to highlight, firstly, the importance of the model in a real production issue for a discrete event system and, secondly, to demonstrate the mathematical model efficiency for a typical case of an industrial production unit. The main goal was to find a schedule that satisfies all resource constraints while minimizing the total cost. Other proposals integrate different parameters as subject of optimization conducive to get the best decisions from which we cite tasks' selection and sizing, assignment of tasks to manufacturing systems, and jobs' sequencing.
In other proposals [14], simulation and meta-heuristics tools are used to highlight operating issues and evaluate different decisions of resources' management. In [15], the authors proposed an algorithm dealing with beds' management issues in hospital centers. Maintenance strategies currently used by decision makers are generally classified into two main classes: The corrective maintenance and the preventive one. Corrective maintenance actions are realized when a machine breaks down. The main task of such maintenance activities is to restore the machine to a desired condition. Preventive actions are performed regularly to increase the machine availability and to avoid, as possible, its breaking down. In this paper, we studied preventive maintenance actions that are already planned in the initial production planning.
The junction of maintenance and production scheduling [16] was the subject of many proposals in the literature. In particular, [17] discusses the optimal joint production and maintenance plan in the case of a system consisting of several identical machines. In [18], a new approach for joint production and maintenance scheduling in the case of multiple machines was developed. Machines are subject to failures and an age-dependent preventive maintenance policy was used. The main improvement of this work was to implement a PN supervisor on a flexible manufacturing system (FMS) under production and control constraints using a developed JAVA application. Therefore, we considered merging an efficient mathematical model in order to provide an optimal production planning. This paper is organized thusly: Section 2 presents briefly the Petri nets and the theory of regions concepts and describes the FMS in our laboratory. Then, Section 3 develops the production planning mathematical model. Next, Section 4 provides a numerical example on a real experimental environment. Finally, the Section 5 shows the mathematical model efficiency on a given industrial production unit.

Mathematical Model Description
In this part, a planning for manufacturing systems production is presented. Our proposition was to develop a production plan under production constraints. The mathematical model was based on [15]'s proposal.
The proposed model was based on the following assumptions:

•
The production duration was considered as known in advance.

•
Batches could not be operated on more than one manufacturing system.

•
Each batch was associated to two dates: Earliest manufacturing date and latest manufacturing date. Real manufacturing date must be included between these two dates.

•
Two type of workshops: Workshops with two manufacturing systems and ones with only one manufacturing system.

•
Manufacturing of different types of batches in the same workshop at the same time was forbidden.
Hence, the essential objective consisted of assigning the maximum number of batches. Assignment must have respected constraints and did not disturb current planning. Availability (empty and occupied manufacturing systems) was determined and updated after each assignment.

Algorithm for Production Scheduling
The objective function is expressed as following: It aims to minimize the non-operated (rejected) batches.
The model constraints are given as following: NMS l=1 X l,t,i ≤ 1.
Equation (2), a batch cannot be operated by more than one manufacturing system in a time unit t.
Equation (4), non-interrupted manufacturing. The model assumes that manufacturing operation is continuous. The length of stay of the batch in the workshop is strictly equal to its manufacturing duration.
Equations (5) and (6), a batch can occupy at most one manufacturing system.
Equation (9), availability of the manufacturing system to operate the batch i.
Equations (10)- (12), batches of different types cannot be operated in the same workshop at a given unit of time.
Equation (13) expresses the preventive maintenance constraints. Constraint (18) expresses the preventive maintenance constraints. Related to the maintenance planning, we considered the fact that manufacturing systems must be maintained after a given number of operated products. Preventive maintenance is not linked to duration but to a given quantity produced, after which the system must be inspected and maintained. This quantity is mostly fixed based on system characteristics such as the equipment life duration, its mechanic characteristics, and the environmental and operational conditions in which the system operates. It is interpreted to mean that a maintenance action can not interrupt the manufacturing operation, which has to be continuous. This constraint is an extension for the constraint (9) for system availability.

Petri Net and Theory of Regions' Concepts
A Petri net (PN) is a graphical modeling and mathematical tool for divers' systems, especially for discrete event systems. From an informal standpoint, a PN is a directed graph with two types of nodes (transitions and places) and has a dynamic behavior. The transitions and places are connected by oriented arcs. An arc links either a transition to a place/a place to a transition but never links a place to a place/a transition to a transition.
Regularly, a PN is a bipartite graph: PN = P, T, Pre, Post where: [19] • P is a set of places.

•
T is a set of controllable and uncontrollable transitions: Pre : P × T in N represent a pre-incidence function that specifies weighted arcs from P to T. N is a set of non-negative integers.

•
Post : P × T in N is the post-incidence function that specifies weighted arcs from T to P.
Let the set p (t) (respectively, (t) p) be the output transitions (respectively, input transitions) of a place p. Naturally, a PN can be formalized by an indexed matrix C as: i , and 0 elsewhere. The w : F → N is a valuation function of arcs (the finite set of arcs F ⊆ (P * T) ∪ (T * P). Moreover, let t (p) (respectively, (p) t) be the set of output places (respectively, input places) of a transition t. The reachability graph constructed from the initial marking The theory of regions in control synthesis was proposed by [19][20][21][22] using the properties of PN tools for adding a PN controller Pc to the initial Petri net model with initial marking M 0 (P c ) and its incidence vectors C(P c , ·). This theory is represented by a linear system composed of three types of equations: The reachability equations, the marking/transition separation instance (MTSI) equations, and the cycle equations. Its resolution leads to the design of PN supervisor. This control place P c has to fulfill the reachability conditions, i.e., where → Γ M is the path between M 0 and M. In addition, each PN controller to add has to solve at least one prohibited event (M, t) in the set Ω of prohibited state transitions. The MTSI condition relative to the couple (M * , t) is: where M * is a forbidden state. Finally, P c must satisfy cycle equations: is the algebraic sum of occurrences of t in σ. Considerably, many MTSI conditions may obtain the same solutions by resolving the linear system composed of Equations (15)- (17). Afterwards, the number of Petri net controllers will be much smaller than the set of MTSI conditions. The following example of Figure 2 is an application of the theory of regions for control synthesis. where is the path between and . In addition, each PN controller to add has to solve at least one prohibited event ( , ) in the set Ω of prohibited state transitions. The MTSI condition relative to the couple ( * , ) is: where * is a forbidden state. Finally, must satisfy cycle equations: where [ is the algebraic sum of occurrences of in . Considerably, many MTSI conditions may obtain the same solutions by resolving the linear system composed of Equations (15)-(17). Afterwards, the number of Petri net controllers will be much smaller than the set of MTSI conditions. The following example of Figure 2 is an application of the theory of regions for control synthesis.

Real Case Presentation
A flexible manufacturing system (FMS), composed of two production stations installed in the university UFR_MIM in Metz-FRANCE and its Petri net model, are presented in Figures 3 and 4. These stations, modeled by the places P2 and P3, pack engraved blue and red pieces in black boxes. The capacity of each station is modeled by P4 and P5. The transitions (t 1 , t 2 ) and (t 4 , t 5 ) represent, respectively, the events of entry and exit of the pallets through their annex's conveyors. Otherwise, the pallets can continue their way to the main conveyor by crossing t 3 . Note that only the transitions t 1 and t 2 are controllable.

Real Case Presentation
A flexible manufacturing system (FMS), composed of two production stations installed in the university UFR_MIM in Metz-FRANCE and its Petri net model, are presented in Figures 3 and 4. These stations, modeled by the places P2 and P3, pack engraved blue and red pieces in black boxes. The capacity of each station is modeled by P4 and P5. The transitions (t1, t2) and (t4, t5) represent, respectively, the events of entry and exit of the pallets through their annex's conveyors. Otherwise, the pallets can continue their way to the main conveyor by crossing t3. Note that only the transitions t1 and t2 are controllable.
The control part of the system makes it imperative not to pack two pieces of different colors for the same unit of time. Our objective in this example was to sort a number N of blue and red engraved pieces in order to package them in black boxes and ship them to customers. In fact, in our previous works, the control constraint was feasible on a number N of pieces in order to dispatch the engraved parts, but the sorting according to the color of the piece remains an insoluble problem in our field. The application of the theory of regions to design the controller was not suitable until now for colored Petri nets. In addition, this work shed new light on the integration of the scheduling in our previous work. In the next section, we develop a mathematical model that will be used to define the optimal production and maintenance planning for the described FMS.  In this section, the algorithm was simulated to highlight its importance in some modeling complex cases. Related to the maintenance planning, we considered the fact that the system must be

Real Case Presentation
A flexible manufacturing system (FMS), composed of two production stations installed in the university UFR_MIM in Metz-FRANCE and its Petri net model, are presented in Figures 3 and 4. These stations, modeled by the places P2 and P3, pack engraved blue and red pieces in black boxes. The capacity of each station is modeled by P4 and P5. The transitions (t1, t2) and (t4, t5) represent, respectively, the events of entry and exit of the pallets through their annex's conveyors. Otherwise, the pallets can continue their way to the main conveyor by crossing t3. Note that only the transitions t1 and t2 are controllable.
The control part of the system makes it imperative not to pack two pieces of different colors for the same unit of time. Our objective in this example was to sort a number N of blue and red engraved pieces in order to package them in black boxes and ship them to customers. In fact, in our previous works, the control constraint was feasible on a number N of pieces in order to dispatch the engraved parts, but the sorting according to the color of the piece remains an insoluble problem in our field. The application of the theory of regions to design the controller was not suitable until now for colored Petri nets. In addition, this work shed new light on the integration of the scheduling in our previous work. In the next section, we develop a mathematical model that will be used to define the optimal production and maintenance planning for the described FMS.  In this section, the algorithm was simulated to highlight its importance in some modeling complex cases. Related to the maintenance planning, we considered the fact that the system must be The control part of the system makes it imperative not to pack two pieces of different colors for the same unit of time. Our objective in this example was to sort a number N of blue and red engraved pieces in order to package them in black boxes and ship them to customers. In fact, in our previous works, the control constraint was feasible on a number N of pieces in order to dispatch the engraved parts, but the sorting according to the color of the piece remains an insoluble problem in our field.
The application of the theory of regions to design the controller was not suitable until now for colored Petri nets. In addition, this work shed new light on the integration of the scheduling in our previous work. In the next section, we develop a mathematical model that will be used to define the optimal production and maintenance planning for the described FMS.
In this section, the algorithm was simulated to highlight its importance in some modeling complex cases. Related to the maintenance planning, we considered the fact that the system must be maintained after a given number of operated batches. Preventive maintenance is not linked to duration but to a given quantity produced, after which the system must be inspected and maintained. This quantity is mostly fixed based on system characteristics such as the equipment life duration, its mechanic characteristics, and the environmental and operational conditions in which the system operates.
The data that will be used in the numerical example are expressed as following: • The number of machines is two, as mentioned in the previous subsection.

•
The two machines are installed in the same workshop.

•
Two types of batches are considered: Blue and red.

•
The number of production batches is 8 (N = 8).

•
The time horizon is expressed in hours and is equal to seven.

•
Production time of a simple piece is 10 min (1/6 h).

•
Each production batch has an execution time proportional to its size, which must be a multiple of 6. For example, a batch having 12 pieces, takes 12*10min = 2 h of production.

•
The batches' characteristics (type, earliest and latest manufacturing dates) are given in the Table 1.  Based on this data, the production program given by the mathematical model application is illustrated as following: The simulation results show a deterministic optimal solution for production and maintenance scheduling. Machines' unavailability is expressed by three letters: B (respectively, R) for pre-programmed production for a blue batch (respectively, red batch) and M for preventive maintenance.

Batch i EM (i) LM (i) T (i) Batch Size (i) Production Time
The results shown in the Table 2 led us to conclude that the model constraints were respected (machines' unavailability and non-mixed batches in each time unit. For example, batch 2 was operated in the station 2, during the third and fourth time units. From a total of 11 batches, one batch was not scheduled for operation. The batch number 6 was rejected. This rejection was due to the production constraints (batch size, earliest and latest manufacturing dates, etc.).
Conforming to our new control methodology, 2 PN monitors, Pc1 and Pc2, can then be located by eliminating redundant PN monitors when the linear system composed of reachability MTSI conditions and cycle equations is solved. The controlled system net is obtained by applying the theory of regions presented in the previous section and shown in the Table 3 and the Figure 5.  Table 3. The controlled Petri net model of FMS.

PN Controllers
The Initial Marking The Incidence Vector In fact, the JAVA application was developed to model any PN controller and then generate the corresponding LIST code that will be sent to the target in a second time. The interface is functional from the creation of the project until the code generation. It is possible to save a project, to model the controlled PN as in any editor, and to associate the inputs to their addresses in the PLC. Then, the generated code (Figure 7) will be tested using the SIEMENS STEP7 compiler (Figure 8). The proposed policy was implemented in the flexible manufacturing system installed in our laboratory in France and the computation of PN controllers was done by using CPLEX software. The results analysis showed that our control approach using PN was more effective than previous methodologies based on the theory of regions. Indeed, for the same number of controllers, one can supervise the system under less computational burden. Our proposed control approach was merged with mathematical modeling tools to face one of the real constraints in the industrial world.
The proposed control policy was implemented in the FMS through a JAVA application ( Figure 6), (Rezig et al. 2018) [23], wherein the controlled PN model was edited and then transformed into a comprehensible code supported by the Programmable Logic Controller (PLC). Once the internal input/output signals of the PLC are detected, one can transform the PN model into a language supported by the PLC. The LIST language of IEC norm 61131-3 was chosen.
In fact, the JAVA application was developed to model any PN controller and then generate the corresponding LIST code that will be sent to the target in a second time. The interface is functional from the creation of the project until the code generation. It is possible to save a project, to model the controlled PN as in any editor, and to associate the inputs to their addresses in the PLC. Then, the generated code ( Figure 7) will be tested using the SIEMENS STEP7 compiler (Figure 8).

Mathematical Model Extensions
The main advantage of the current proposed model is its mathematical flexibility that permits adapting it to many production applications by modifying the model constraints.
In this section, we present some different applications where the model has shown its efficiency to validate the flexibility characteristic.

Case of Manufacturing System
This application was discussed in an existing work in the literature [24] with an algorithm for manufacturing planning under unavailability constraints. In this section, we extended the database used for the numerical example to show the model efficiency dealing with a bigger size problem. Table 4 gives the initial unavailability dates of a manufacturing system in the time horizon that will be the subject of scheduling. As mentioned before, the unavailability, in this paper, corresponds to preprogrammed production activities: Based on Table 4, as an example, preprogrammed production actions in the studied time horizon for the system 1 are the manufacturing of batch type B for fifth and sixth time units.

Mathematical Model Extensions
The main advantage of the current proposed model is its mathematical flexibility that permits adapting it to many production applications by modifying the model constraints.
In this section, we present some different applications where the model has shown its efficiency to validate the flexibility characteristic.

Case of Manufacturing System
This application was discussed in an existing work in the literature [24] with an algorithm for manufacturing planning under unavailability constraints. In this section, we extended the database used for the numerical example to show the model efficiency dealing with a bigger size problem. Table 4 gives the initial unavailability dates of a manufacturing system in the time horizon that will be the subject of scheduling. As mentioned before, the unavailability, in this paper, corresponds to preprogrammed production activities: Based on Table 4, as an example, preprogrammed production actions in the studied time horizon for the system 1 are the manufacturing of batch type B for fifth and sixth time units.
For the numerical example, we considered five workshops having eight manufacturing systems (three double workshops and two single). The following database was considered (Table 5):

•
The horizon time was fixed at 14 units of time.

•
The manufacturing duration (MD) was fixed at three units of time for batch type A and two units of time for batch type B. • Initial unavailability was expressed by the preventive maintenance action plan by the letter M in the Tables 4 and 6 and preprogrammed production actions expressed by Letter B for batch type B and Letter A for batch type A.    Table 6 illustrates the obtained results. The simulation results highlight the optimal production scheduling in the given time horizon and in function of model constraints. Results led us to conclude that production constraints, that two batches of different types could not be in the same workshop at a given time t, were respected. For instance, batch 2 was operated in the first workshop, on the first machine during the second, third, and the fourth time units. Batch 5 was operated on the third machine of the Workshop 3 on the 10th, the 11th, and 12th time units.
The optimal solution was obtained with the solver CPLEX, after 2 min and 50 s of calculation, which corresponds to 190,134 iterations. From a total of 11 batches, one batch was not scheduled for operation. The batch was rejected due to the earliest and latest manufacturing dates' constraint.
We proposed a performance comparison between the solvers used, which were, in this case, CPLEX and Xpress. We considered five experiments to schedule 15,20,25,30, and 50 products. The results obtained are illustrated in Table 7. We note that the results obtained with CPLEX and Xpress were nearly the same for the case of 15 products. For the case of 20 and 25 products, we note that CPLEX was faster with the same given optimal solution. Moreover, in the case of 50 products, CPLEX reached the optimal solution within a few seconds, whereas Xpress spent much more time before getting a solution, which was not the optimal one. We can conclude then that, in the case of this algorithm resolution, it would be better to use the CPLEX solver.

Case of Hospital Beds' Management
The proposed mathematical model was also used for beds' management planning and management in healthcare establishments [15].
Beds' management is a complex problem that lots of hospital centers decide to integrate, to well manage the financial impact on the budget and to maximize the acceptance rate of patients.
Many proposed models are based on reducing hospitalization cost and optimizing the use of resources. In our proposal [15], the model considered the capacity of beds for two types of patients: Scheduled and non-scheduled patients. Otherwise, the model managed also preventive maintenance actions during the planning horizon.
The idea of our work was to apply the proposed mathematical algorithm with some modifications to be in good agreement with real specific constraints of this application from which we can cite:

•
The need of considering some cases where some patients need to be isolated, so assigned to a room with a single bed.

•
The need of replanning the refused patients given from the first planning proposal, by modifying, if possible, the earliest and latest acceptance dates.
The model was used in an analogical way with the example presented for manufacturing system. The workshops designed the rooms, the manufacturing systems designed the beds, and the products were assimilated to patients. For more details and numerical illustrations, the reader may refer to our previous work [15].
The actual proposed mathematical model could also be extended to several other applications with adaptations on constraints, objective function, and parameters characteristics such as museum system or the study of big data in complex and social networks [25][26][27][28].
For instance, the problem of visitors' museum touring could be addressed, while considering physical parameters, as well as visitors' interests and social interactions, towards maximizing their perceived quality of experience (QoE). In this case, it would be interesting to develop this proposal by establishing an analogy between hospital environment (hospital, patients, beds, and characteristics) and museum case (museum, visitors, resources, and exhibits' characteristics).

Conclusions
This work aimed to propose a new approach in mathematical modeling optimization. The contribution of this paper was, firstly, to propose a linear optimization model for a joint production and maintenance planning and, secondly, to show the importance of such models in different industrial issues and particularly in the case of flexible manufacturing system where the Petri net modeling is more commonly used. Furthermore, an application example of our methodology, such as hospital beds management, was presented in this paper to show the effectiveness of our mathematical model. Thus, we backed up our proposals by another illustrative example and a real case of study. Therefore, the algorithm for solving an optimal production planning under production and unavailability constraints was discussed.
Thanks to the modeling tools, decision makers could optimize the production planning which depends on several constraints. The mathematical model was then applied on a discrete event system installed in our laboratory (University of Lorraine, Metz, France). Even though the main used tool for the flexible manufacturing system was the PN modeling, we showed in this article that it could be controlled with mathematical modeling, in the case of complex production issues and, so, there is no need for more complex tools in PN, such as colored PN tool. Therefore, thanks to the mathematical model solution output, the PN supervisor was implemented via STEP7 software.
Future works aim to integrate other production constraints in the manufacturing planning algorithm to more generalize the model usefulness for other industrial applications.
Author Contributions: The authors of this paper are specialists in the industrial engineering field and more precisely in production and maintenance planning for optimizing for FMS scheduling. Our teamwork tested the control model on a (FMS) flexible manufacturing system installed in our laboratory (LGIPM) by using all the appropriate software (SCHNEIDER Unity Pro XL and Siemens STEP7) for control software, simulation, and modeling software, CPN Tools, and, finally, supervision software, WINCC FLEXIBLE. The decision variables of this model are as follows: X i,l,t : 1 if batch i is assigned in manufacturing system l during the period t, 0 otherwise. F i : 0 if batch i is operated, 0 otherwise. A i,l : 1 if batch i is assigned in manufacturing system l, 0 otherwise. J i : Indicates the real manufacturing date End i : Indicates the end manufacturing date Variables i, l, and t have to be included respectively in [1,NP], [1, L], and [1, T].