Mathematical Problems in Emerging Manufacturing Systems Management
View this Special IssueResearch Article  Open Access
Joint Optimization Approach of Maintenance and Production Planning for a MultipleProduct Manufacturing System
Abstract
This paper deals with the problem of maintenance and production planning for randomly failing multipleproduct manufacturing system. The latter consists of one machine which produces several types of products in order to satisfy random demands corresponding to every type of product. At any given time, the machine can only produce one type of product and then switches to another one. The purpose of this study is to establish sequentially an economical production plan and an optimal maintenance strategy, taking into account the influence of the production rate on the systemâ€™s degradation. Analytical models are developed in order to find the production plan and the preventive maintenance strategy which minimizes sequentially the total production/inventory cost and then the total maintenance cost. Finally, a numerical example is presented to illustrate the usefulness of the proposed approach.
1. Introduction
Manufacturing companies must manage several functional capacities successfully, such as production, maintenance, quality, and marketing. One of the keys to success consists in treating all these services simultaneously. On the other hand, the customer satisfaction is one of the first objectives of a company. In fact, the nonsatisfaction of the customer on time is often due to a random demand or a sudden failure of production system. Therefore, it is necessary to develop maintenance policies relating to production, reducing the total production and maintenance cost. One of the first actions of decisionmaking hierarchy of a company is the development of an economical production plan and an optimal maintenance strategy.
It is necessary to find the best production plan and the best maintenance strategy required by the company to satisfy customers. This is a complex task because there are various uncertainties due to external and internal factors. External factors may be associated with the inability to precisely define the behaviour of the application during periods of production. Internal factors may be associated with the availability of hardware resources of the company. In this context, Filho [1] treated a stochastic scheduling problem in terms of production under the constraints of the inventory.
Establishing an optimal production planning and maintenance strategy has always been the greatest challenge for industrial companies. Moreover, during the last few decades, the integration of production and maintenance policies problem has received much research attention. In this context, Nodem et al. [2] developed a method to find the optimal production, replacement/repair and preventive maintenance policies for a degraded manufacturing system. Gharbi et al. [3] assumed that failure frequencies can be reduced through preventive maintenance, and developed joint production and preventive maintenance policies depending on produced part inventory levels. An analytical model and a numerical procedure which allow determining a joint optimal inventory control and an age based on preventive maintenance policy for a randomly failing production system was presented by Rezg et al. [4].
This work examined a problem of the optimal production planning formulation of a manufacturing system consisting of one machine producing several products in order to meet several random demands. This type of problem was studied by Kenne et al. [5]. They presented an analysis of production control and corrective maintenance problem in a multiplemachine, multipleproduct manufacturing system. They obtained a near optimal control policy of the system through numerical techniques by controlling both production and repair rates. Feng et al. [6] developed a multiproduct manufacturing systems problem with sequence dependent setup times and finite buffers under seven scheduling policies. Sloan and Shanthikumar [7] presented a Markov decision process model that simultaneously determines maintenance and production schedules for a multipleproduct, singlemachine production system, accounting for the fact that equipment condition can affect the yield of different product types differently. Filho [8] developed a stochastic dynamic optimization model to solve a multiproduct, multiperiod production planning problem with constraints on decision variables and finite planning horizon.
Looking at the literature on integrated maintenance policies, we noticed that the influence of the production rate on the degradation system over a finite planning horizon was rarely addressed in depth. Recently, Zied et al. [9â€“11] took into account the influence of production plan on the equipment degradation, in the case of a system composed of single machine producing one type of product under randomly failing and satisfying a random demand over a finite horizon. In the same context, Kenne and Nkeungoue [12] proposed a model, where the failure rate of a machine depends on its age; hence, the corrective and preventive maintenance policies are machineage dependent.
Motivated by the work in the Zied et al. [9â€“11], we treat the production and maintenance problem in another context that we consider a more complex and real industrial system, composed of one machine that produces several products during a finite horizon divided into subperiods. This study displays that it has a novelty and originality relative to this type of problem which considers the influence of several products on the degradation degree of the considered machine and consequently on the average number of failure as well as on the maintenance strategy.
This paper is organized as follows: Section 2 states the problem. Section 3 presents the notations. The production and maintenance models are developed, respectively, in Sections 4 and 5. A numerical example and sensitivity study are presented, respectively, in Sections 6 and 7. Finally, the conclusion is included in Section 8.
2. Statement of the Industrial Problem
This study treated an industrial case. The problem concerns a textile company located in North Africa specialized in clothing manufacturing. The companyâ€™s production system consists of a conversion of three types of fiber into yarn, then fabric, and textiles. These are then fabricated into clothes or other artefacts. The production machine is called the loom, and it uses a jet of air or water to insert the weft. The loom ensures pattern diversity and faultless fabrics by a flexible and gentle material handling process. Fabrics can be in one plain color with or without a simple pattern or they can have decorative designs.
Based on the industrial example described, this study was conducted to deal with the problem of an optimal production and maintenance planning for a manufacturing system. The system is composed of a single machine which produces several products in order to meet corresponding several random demands. The problem is presented in (Figure 1).
The considered equipment is subject to random failures. The degradation of the equipment increases with time and varies according to the production rate. The machine is submitted to a preventive maintenance policy in order to reduce the occurrence of failures. In the literature, the influence of the production rate on the material degradation is rarely studied. In this study, this influence was taken into consideration in order to establish the optimal maintenance strategy.
The model developed in this study is based on the works of Zied et al. [9â€“11]. These studies seek to determine an economical production plan followed by an optimal maintenance policy but for the case of only one product.
Firstly, for a randomly given demand, an optimal production plan was established to minimize the average total storage and production costs while satisfying a service level. Secondly, using the obtained optimal production plan and considering its influence on the manufacturing system failure rate, an optimal maintenance schedule is established to minimize the total maintenance cost.
3. Notations
In this paper, we shall as far as possible use the notation summarized as follows:â€‰: unit production cost of product ,â€‰: holding cost of one unit of product during ,â€‰: setup cost of product ,â€‰: corrective maintenance action cost,â€‰: preventive maintenance action cost,â€‰: total number of periods,â€‰: total number of products,â€‰: total number of subperiods during each period,â€‰: production period duration,â€‰: nominal production quantity of product during ,â€‰: probabilistic index (related to customer satisfaction) of product ,â€‰: demand of product during period ,â€‰: inventory level of product at the end of subperiod of period ,â€‰: the total expected cost of production and inventory over the finite horizon,â€‰: the demand variance of product at period ,â€‰: cumulative Gaussian distribution function,â€‰: inverse distribution function,â€‰: the total cost of maintenance,â€‰: failure rate function at subperiod of the period ,â€‰: nominal failure rate,â€‰: the average number of failures,â€‰: intervention period for preventive maintenance actions.
Decision Variablesâ€‰: production quantity of product during subperiod of period ,â€‰: duration of subperiod at period ,â€‰: a binary variable, which is equal to 1 if product is produced in subperiod of the period , and 0 otherwise,â€‰: number of preventive maintenance actions during the finite horizon.
4. Production Policy
In this section, we developed an analytical model which minimizes the total cost of production and storage. The decision variables are the production quantities , the binary variable , and the duration of subperiods . Our objective consists in determining an economical production plan for a finite time horizon . The production plan must satisfy random demands under the requirement of a given level of service, while minimizing the cost of production and storage. The production of each product will take place at the beginning of subperiods, and delivery to the customer will be at the end of periods.
The state of the stock is determined at the end of each subperiod. Figure 2 shows an example of a production plan.
4.1. Stochastic Model of the Problem
To develop this section, the following assumptions are specifically made:(i)holding and production costs of each product are known and constant;(ii)only a single product can be produced in each subperiod;(iii)as described in (Figure 2), we have divided the period into equal subperiods, with (the total number of products);(iv)the standard deviation of demand and the average demand for each product and each period are known and constant.The model has the following basic structure:under the constraints below:(i)the inventory balance equation,(ii)the service level,(iii)the admissibility of production plan,(iv)the maximum production capacity.
Formally
(i) The Cost Functions. Consider
(ii) The Inventory Balance Equation. The available stock at the end of each subperiod of period for each product is formulated in the form of flow balance constraints (inflow = outflow):where is the initial stock level of product .
This equation shows that the stock of product at the end of each subperiod of each period is determined by the state of the stock of product at the end of the subperiod .
(iii) The Admissibility of Production Plan and Service Level Constraints. The service level of product is determined by the probability constraint on the stock level at the end of each period :We can transform the probabilistic constraint of stock level to a deterministic constraint.
Formally, the function becomeswhere is the estimated demand of product during the period , is the minimum stock level of product required at the end of period , and is the initial stock level of product .
(iv) The Maximum Production Capacity. The production quantity of the machine for each product , , is limited and is presented as follows:The term allows taking into account the influence of duration of subperiods on the maximum quantity of production. If tends to 0, the maximum quantity of production tends also to 0, and if tends to , the maximum quantity of production tends to (with Nominal production quantity of product during ).
However, the term represents the maximum production quantity of product during the subperiod of period .
4.2. Problem Formulation
We recall that, in this study, we assume that the horizon is divided into equal periods and each period is divided into subperiods with different durations. Figure 2 shows the distribution of the production plan for the finite horizon . Each product is produced in a single subperiod in each period . The demand of each product is satisfied at the end of each period .
The mathematical formulation of the proposed problem is based on the extension of the model described by Zied et al. [11] for the one product case study.
Their problem is defined as follows:where is unit production cost and is holding cost of a product unit during the period .
Formally, our stochastic production problem is defined as follows:withwhere is the mathematical expectation.
Under the following constraints:The first constraint stands for the inventory balance equation for each product , during each subperiod , , of period , . Equation (11) refers to the satisfaction level of demand of product in each period . Constraint (12) defines the upper production quantity of the machine for each product . The aim of (13) is to divide each period into different subperiods.
The constraints below should also be taken into account:Equation (14) indicates that only one type of product will be produced in subperiod of period . Constraint (15) states that is a binary variable. We note that is equal to 1 if product is produced in subperiod of the period , and 0 otherwise.
For each subperiod of period , the equation of the stock status is determined by the first constraint. This equation remains random because of the uncertainty of fluctuating demand. Therefore, the variables of production and storage are stochastic. Their statistics depend on a probabilistic distribution function of demand. It is, therefore, necessary to use constraint (11) for decision variables. These constraints can help us to analyse the various production scenarios to improve the performance of the production system.
4.3. The Deterministic Production Model
We admit that a function , represents the cost of storage and production which is relative to the proposed plan and represents the value of the mathematical expectation. The quantity stocked of product at the end of the subperiod of period is stood for by . The production quantity required to satisfy the demand of product at the end of period is , where represents the subperiod during which the product is produced.
Thus, the problem formulation can be presented as follows:The purpose, then, is to determine the decision variables (, and ) required to satisfy economically the various demands under the constraints seen in the previous subsection.
The resolution of the stochastic problem under these assumptions is generally difficult. Thus, its transformation into a deterministic problem facilitates its resolution.
(i) Inventory Balance Equation. The stochastic inventory balance equation iswith being the initial stock level of product .
We suppose that the means and variance of demand are known and constant for each product in each period .
Therefore,The inventory equation is statistically described by its means:We note thatbecause is constant for each interval .
AndThen, the balance equation (10) can be converted into an equivalent inventory balance equation:with being the average initial stock level of product.
(ii) Service Level Constraint. The second step is to convert the service level constraint into a deterministic equivalent constraint by specifying certain minimum cumulative production quantities that depend on the service level requirements.
Lemma 1. Consider the following:
Proof. We know thatNoting that is a Gaussian random variable for demand .
Hence,We recall that represents the probabilistic index (related to customer satisfaction) of product and represents the demand variance of product at period .
The distribution function is invertible because it is an increasing and differentiable function.
Hence,Therefore,
(iii) The Expression of the Total Production and Storage Cost. In this step, we proceed to a simplification of the expected cost of production and storage.
The expression of the total cost of production is presented as follows.
Lemma 2. Consider the following:
Proof. See Appendix A.
(iv) In Summary. The deterministic optimization problem becomes as follows.
(a) The Objective Function. Consider
(b) The Constraints. Consider
5. Maintenance Strategy
5.1. Description of the Maintenance Strategy
The maintenance strategy adopted in this study is known as preventive maintenance with minimal repair. The actions of preventive maintenance are practiced in the period . The replacement rule for this policy is to replace the system with another new system (as good as new) at each period . At each failure between preventive maintenance actions, only one minimal repair is implemented. If we note , the cost of corrective maintenance actions, and , the cost of preventive maintenance actions and degradation of the machine is linear, the total cost of maintenance is expressed as follows:To develop the analytical model, it was assumed that(i)durations of maintenance actions are negligible;(ii) and costs incurred by the preventive and corrective maintenance actions are known and constant, with ;(iii)preventive maintenance actions are always performed at the end of the subperiods of production.The aim of this maintenance strategy is to find the optimal number of preventive maintenance actions minimizing the total cost of maintenance over a given horizon . The existence of an optimal number of partitions and, therefore, the optimal preventive maintenance period is proven in the literature. It has been proven that exists if the failure rate is increasing [13].
Before determining the analytical model minimizing the total cost of maintenance, we need first to develop the expression of the failure rate and then the average number of failures expression , during the finite horizon .
5.2. Expression of Failure Rate
Recall that the key of this study is the influence of the variation of the production rates on the failure rate.
Figure 3 represents the general description of the evolution of the failure rate, which depends on both the production rate and the failure rate of the previous period.
As presented in Figure 3, the failure rate is reset after each , with .
Thus, the expression of the failure rate depending on time and production rate can be written as follows:The term represents the production rate of product during subperiod of period .
The term represents the nominal production rate of product during .
Therefore,The aim of the expression is to reset the failure rate after each with .
Note thatwhere is the integer part of number .
Lemma 3. Consider the following:
Proof. See Appendix B.
5.3. Expression of the Average Number of Failures
In order to reduce the complexity of the generation of the optimal number of preventive maintenance, we assume that interventions are made at the end of subperiods.
Hence, the function of the period of intervention is presented as follows:where is a round number of .
To determine the average number of failures expression during the finite horizon , we will focus on the calculation of the average number of failures during the interval , which we designate . Hence, we have to calculate the three surfaces , , and mentioned in Figure 4.
Therefore, the average number of failures expression during the interval is presented as follows:where is the superior integer part of number .
Thus, the average number of failures expression during the finite horizon is defined byTherefore, we have the following lemma.
Lemma 4. Consider the following: Note that .
5.4. Expression of the Total Cost of Maintenance
We recall that the initial expression of the total cost of maintenance presented in (33) isUsing the average number of failures established in Lemma 4, we can deduce that the analytical expression of the total maintenance cost is expressed as follows: The goal is to find the optimal number of preventive maintenance actions that minimizes the total cost of maintenance . Using this decision variable, we can deduce the optimal period of intervention , knowing that .
5.5. Existence of an Optimal Solution
The following equation determines analytically the optimal solution:Since it is difficult to solve analytically the expression of maintenance cost, we use numerical procedure.
We start by proving the existence of a local minimum.
We have the following.
Limits at the terminals of areNote that is the average number of failures. and represent, respectively, the corrective and the preventive maintenance costs.
Moreover,In addition,In summary, there is an optimal number of partition , which is unique and satisfies the previous relations (46) and (47). The following lemma ensures the existence of a local minimum.
Lemma 5. Consider the following: withTherefore, there exists an optimal number of partition , which satisfies the following expressions:
The resolution of this maintenance policy, using a numerical procedure, is performed by incrementing the number of maintenance intervals until an , satisfying the two first relations in Lemma 5 and minimizing the total cost of maintenance described by (43).
6. Numerical Example
From the industrial example presented in Section 2, we have considered a system producing 3 types of fiber in order to meet three random demands according to every type of product. Using the analytical models developed in previous sections, we start by establishing the optimal production plan and then we determine the optimal maintenance strategy expressed as optimal number of preventive maintenance minimizing the total cost of maintenance over a finite planning horizon: trimesters (two years). We note that the optimal maintenance strategy is obtained while considering of the influence of the production plan on the system degradation. We supposed that the standard deviation of demand of product is the same for all periods. The data required to run this model are given in sequence.
6.1. Numerical Example
(i) The Data Relating to Production. The mean demands (in bobbins) as shown in Table 1:

The other data are presented as shown in Table 2.

(ii) The Data Relating to System Reliability. System reliability and costs related to maintenance actions are defined by the following data:(1)the law of failure characterizing the nominal conditions is Weibull. It is defined by(a)scale parameter : 12 months,(b)shape parameter : 2,(c)position parameter : 0,(2)the initial failure rate: .These parameters provide information on the evolution of the failure rate in time.
This failure rate is increasing and linear over time. Thus, the function of the nominal failure rate is expressed byThe preventive and corrective maintenance costs are, respectively, â€‰mu and â€‰mu.
6.2. Determination of the Economic Production Plan
The economic production plan obtained is presented in Table 3.

6.3. Determination of the Optimal Maintenance Plan
As described in Figure 5, the optimal maintenance strategy is obtained based on the optimal production plan given in the previous section.
Figure 6 shows the curve of the total cost of maintenance according to (number of preventive maintenance actions). We conclude that the optimal number of preventive maintenance actions that minimizes the total cost of maintenance during the finite horizon (two years) is times. Hence, the optimal period to intervene for the preventive maintenance is months, and the minimal total cost of maintenance â€‰mu.
7. The Economical Profit of the Study
We recall that the specificity of this study is that it considered the impact of the production rate variation on the system degradation and consequently on the optimal maintenance strategy adopted in the case of multiple product. In order to show the significance of our study we will consider, in this section, the case of not considering the influence of the production rate variation on the systemâ€™s degradation. That is to say, we assume that the manufacturing system is exploited at its maximal production rate every time. Analytically, we will consider the nominal failure rate which depends only on time. The results of this study are presented in Table 4.

The optimal number of preventive maintenance obtained in the case when we did not consider the variation of production rate is times and it corresponds to a total cost of maintenance during the finite horizon (two years), â€‰mu. We recall that in our case study when we consider the variation of production rate we have obtained and â€‰mu. We can easily note that we have reduced the optimal number of preventive maintenance with performing an economical gain estimated at 10%.
Several studies have addressed issues related to production and maintenance problem. But, the consideration of the materiel degradation according to the production rate in the case of multipleproduct has been rarely studied.
This study was conducted to deal with the problem of an optimal production and maintenance planning for a manufacturing system. The significance of the present study is that we took into account the influence of the production plan on the system degradation in order to establish an optimal maintenance strategy. The considered system is composed of a single machine which produces several products in order to meet corresponding several random demands.
8. Conclusion
In this paper, we have discussed the problem of integrated maintenance to production for a manufacturing system consisting of a single machine which produces several types of products to satisfy several random demands. As the machine is subject to random failures, preventive maintenance actions are considered in order to improve its reliability. At failure, a minimal repair is carried out to restore the system into the operating state without changing its failure rate.
At first we have formulated a stochastic production problem. To solve this problem, we have used a production policy to achieve a level of economic output. This policy is characterized by the transformation of the problem to a deterministic equivalent problem in order to obtain the economic production plan. In the second step, taking into account the economic production plan obtained, we have studied and optimized the maintenance policy. This policy is defined by preventive actions carried out at constant time intervals. The objective of this policy is to determine the optimal number of preventive maintenance and the optimal intervention periods over a finite horizon. This policy is characterized by a failure rate for a linear degradation of the equipment considering the influence of production rate variation on the system degradation and on the optimal maintenance plan in the case of multiple products represents.
The promising results obtained in this thesis can lead to interesting perspectives. A perspective that we are looking for at the short term, is to consider maintenance durations. We recall that, throughout our study, we neglected the durations of actions of preventive and corrective maintenance. It is clear that the consideration of these durations impacts the optimal maintenance plan and the established production plan. In the medium term, it is interesting to concretely consider the impact of logistics service on the study. It is clear that the inmaintenance logistics are absent in most researches. The combination of maintenance logistics and production represents a motivating perspective in this field of study.
Another interesting perspective specifying the manufactured product can be explored.
Appendices
A. Expression of the Total Production and Storage Cost
We haveAlso,Therefore, and are independent random variables, so we can deduceOn the other hand, we note thatTherefore,We know thatTherefore,Finally,