1. Introduction
Manufacturing companies aim to transform raw materials or components received from their suppliers and to assemble them into finished products to be delivered to their customers. As well as market volatility and the competitive context, the management of this concept differs from one company to another. Supply chain management draws the attention of several researchers in the field of industrial engineering. A growing body of recent books has been published whose researchers present supply chain management with different concepts and techniques in diverse entities related to the logistics chain [
1,
2,
3]. Stadtler [
2] also investigated the essence of supply chain management and advanced planning. He presented the latest research results for the resolution described in his paper. Management approaches in the field of industrial engineering are diverse. On the other side, with the rise of consumption, concepts of supply chain management are extended. Indeed, more and more of production waste, resources reduction, and even environmental degradation encourage us to develop new approaches that manage these problems. Several firms are consequently increasingly responsible for their post-consumer products owing to stricter regulations for managing industrial sustainability [
4,
5]. Accordingly, today, the supply chain is driven by sustainable supply chain management [
6,
7]. Sustainable supply chain management can be defined as the integration of environmental and financial approaches into the entire supply chain lifecycle. All supply chains can be optimized using sustainable approaches. To achieve profitability, supply chain sustainability approaches must deliver improved environmental performance within a financially viable operating construct. Thus, we consider the issue of planning management regarding assembly and disassembly systems in a supply chain under different uncertainties.
The production of durable goods has exploded, resulting in an increased flow of returns from customers to suppliers. This is called “reverse logistics”. Govindan et al. [
8] proposed a review paper. They reviewed and analyzed recently published papers in reverse logistics and closed-loop supply chains to clarify and suggest future research opportunities. In fact, they mentioned that reverse logistics start from end users. The used products are collected from customers and then they attempt to manage them through different decisions (recycling, remanufacturing, repairing, and disposing of). In our study, we consider the remanufacturing process—in other words, the remanufacturing of components derived from good parts of used/end-of-life product disassembly in order to recover them in the assembly process and have a saleable product. The remanufacturing of used products is an area that is very much addressed in the literature on reverse logistics [
9,
10]. These remanufacturing processes are significantly more environmentally friendly than recycling or first-time manufacturing. In fact, these operations use fewer materials and less energy. Most of the concerned firms are specialized in the remanufacture of automobile, truck, and other vehicle components such as starters, alternators, transmissions, etc. [
11]. Many researchers tried to find solutions within their countries by studying problems and roadblocks that prevent the start of the process of remanufacturing [
12]. Therefore, compared to previous studies that deal with disassembly-assembly issues, this work considers the remanufacturing of used product with disassembly and assembly activities.
Disassembly operations are among the main process in remanufacturing process. The objective of these operations is to obtain the valuable components of used/end-of-life products while separating those which affect the environment. In the literature, many researchers were interested in the problems of the disassembly system such as disassembly planning and scheduling, disassembly shop scheduling and control, etc. Priyono et al. [
13] presented a literature review paper where they analyzed what had been done the first time and what needs to be done the second time in the field of disassembly in remanufacturing. Therefore, Bentaha et al. [
14] studied the problem of disassembly process planning considering the states or quality of the products to be disassembled. Meanwhile, the problem of disassembly process planning and the definition of the disassembly sequences have been studied in numerous works; we cite Yu et al. [
15], Sanchez et al. [
16], Ren et al. [
17] and Peeters et al. [
18]. In other studies [
19,
20], authors also considered the problem of disassembly scheduling under a single period of disassembly. In the same context, Ferrer et al. [
11] treated the problem of determining the quantity and timing of disassembling each parent item for a given disassembly structure in order to satisfy the demand of leaf items over a planning horizon with discrete time periods. However, these above-mentioned studies have treated the disassembly process under different uncertainties independently of the manufacturing process. Indeed, in real-life systems, the disassembly operation depends on the demand for manufacturing of finished products.
Our paper is based on Material Requirements Planning (MRP) taking into account remanufacturing components. Within recent decades, MRP systems were among the most popular methods of production and inventory planning for companies, and of interest to many researchers [
21,
22,
23,
24,
25]. However, MRP tends to be an assembly-oriented scheduling system [
17]. Some studies noted that most MRP does not provide capabilities to plan disassembly [
18]. Inderfurth and Jensen [
26] proposed a mathematical analysis of remanufacturing within the MRP framework to develop control rules for undertaking production of new components, refurbishing returned cores, and disposing of excess cores. Their model is also limited to remanufacturing processes that maintain a finished good identity. In the work of Raupp et al. [
27], they presented an MRP optimization model where they treated a particular production system. In addition to manufacturing a final product, its main assembling component can be recovered from parts of the returned product. Gotzel and Inderfurth [
28] developed an extended MRP model for a hybrid single-stage multi-period production and remanufacturing system. Compared to those studies cited above, the MRP examined in our system considering not only finding release dates of new components, but also the order date of used product from the suppliers.
To solve this kind of problem, Taleb and Gupta [
29] proposed an algorithm whose objective it is to determine the quantity and the disassembly schedule for the item to fulfill that demand on time. In another paper [
30], they proposed two companion algorithms to obtain the quantity and disassembly operations schedule; they treated a more complex problem with the existence of multiple product structures and common parts and/or materials. In addition, Langella [
19] dealt with the same problem and proposed a modified version of the algorithm of Taleb and Gupta [
30]. Indeed, Kim et al. [
20] proposed a cost objective function to minimize the sum of setup, inventory holding costs, and disassembly operation. They implemented a heuristic algorithm based on linear programming relaxation to solve the proposed cost function. Kim et al. [
31] treated a production planning problem in disassembly systems with a branch and bound algorithm. Prakash et al. [
32] presented a new form of constraint-based simulated annealing (CBSA) in his paper to determine the ordering and disassembly schedule, to minimize inventory level for products with general assembly product structure.
This paper aims to bridge the gap in research on the disassembly-assembly system by developing an optimal disassembly-assembly plan. Indeed, we consider stochastic lead-times of intermediate disassembled and assembled components. Moreover, we propose an integrated maintenance plan that considers deterioration of the system, which also depends on the disassembly-assembly activities. The integration of a maintenance plan aims to study the risk of disruption of the disassembly-assembly plan caused by machine breakdowns.
This paper is organized as follows: in
Section 2, the global assembly-disassembly system is described as well as the different assumptions considered. An analytical model is presented in
Section 3. A maintenance plan is implemented for the proposed assembly-disassembly system in
Section 4. The last section is denoted by numerical examples as well as sensitivity studies of some parameters that influence the system.
3. Description of the MRP Mathematical Model with Remanufacturing Parts
In this section, an optimization MRP model with remanufacturing parts is considered. First, a mathematical model is developed to determine the optimal ordered date for used/end-of-life product as well as the optimal vector of release dates for components at level 2 (or new raw materials) to satisfy customer demand over a finite horizon. In the second step, from the optimal dates founded, we built a disassembly-assembly plan for the system. The proposed system is disposed to random failures that lead to an unavailability of the system and therefore a delay in the delivery of the customer demand. To reduce breakdown frequency, a mathematical model is implemented to determine the optimal preventive maintenance plan (
Figure 4).
3.1. Disassembly-Assembly Plan
The mathematical model proposed below leads to minimize these costs: inventory holding costs of different components (disassembled components, leaf-items, new components ordered from external resources, assembled components at level 1, finished product), remanufacturing costs, disposal costs, purchasing costs of new components and of used/end-of-life product, and backlogging cost of finished product. This mathematical function optimizes the order date of used/end-of-life product as well as the release dates of new components at level 2 over a finite horizon.
The mathematical expressions mentioned below are used to clarify the analytical proposed model:
The disassembly date
for each component
is issued from the disassembly of the used/end-of-life product. This date represents the initiation of the disassembly process of
. The disassembly process of each component begins after a random lead-time that includes the component transmission, component preparation, etc. If the component is ready to be disassembled and the system is occupied by another operation, the component is stocked until the system is released.
We note by the stock period of each component is which is calculated as follows:
We suppose that:
: If the machine is occupied in .
: If the machine is available in
.
1: | |
2: | |
3: | |
4: | |
5: | |
6: | |
7: | |
8: | |
The assembly date
for each component is
, which begins when all the necessary components required for its assembly are ready. The assembly of the component
Ci,1 consists of a set of new components ordered from external resources and a set of remanufacturing components resulting from the disassembly process of used/end-of-life component. Components arriving early will be stored until all other components are received. At the assembly, if the system is busy, all the components must be stored until the system is released. The expression of the assembly date is as follows:
We note by the stock period of each component is which is calculated as follows:
Where in this expression, the stock is calculated as follow:
: If the machine is occupied in .
: If the machine is available in .
1: | |
2: | |
3: | |
4: | |
5: | |
6: | |
7: | |
The assembly date for the finished product is
: the assembly of the finished product is started when all components of level 1 are ready:
The maximum between the end of the finished product assembly and the planned delivery time is noted by
:
The minimum between the end of the finished product assembly and the planned delivery time is noted by
:
The objective function includes:
The inventory costs of disassembled components:
The inventory costs of components at level 2:
The inventory costs of components at level 1:
Backlogging costs for the finished product:
Inventory holding costs for the finished product:
Therefore, the total function cost denoted by
which is composed of inventory holding costs of the different components (disassembled components issued from the disassembly of the used/end-of-life product, components at level 1 and 2, and the finished product), backlogging costs, purchasing costs (purchasing of the used/end-of-life product and ordered components of level 2 from external resources), remanufacturing, and disposal costs.
With , where , and .
The procurement lead times of used/end-of-life product , disassembled components issued from this latter and components of level 1 and 2 are random variables. The determination of the ordered dates ( and ) and the maintenance plan when minimizing the cost function is done sequentially. First, we determine the ordered dates; the cost function does not include maintenance costs. Following these dates, we determine the disassembly-assembly system plan. This plan will be used to determine the optimal preventive maintenance plan.
To have an optimal disassembly-assembly system plan, we need to resolve the expected total cost function
, which is given by:
3.2. Maintenance Plan
According to optimal dates founded as
and
, we determine the optimal disassembly-assembly plan. This plan allows determination of the system’s solicitation rate, which has an influence on the effectiveness of the system; likewise, it leads to failures. The objective is to find the optimal number of preventive maintenance
N*. We assume that, after each corrective maintenance, the system is in the state just before breakdown. In addition, preventive and corrective actions have negligible times. Between two preventive maintenance actions, the failure rate depends on the disassembly or assembly processes. Each disassembly or assembly process solicits the machine at a solicitation rate, which depends on the process type and relative to the highest level of system solicitation. We propose in this paper an example of solicitation levels
u for disassembly types in
Table 1 and assembly types in
Table 2. We take the example of five types of components which can be assembled together or disassembled from each other; each disassembly or assembly operation of a specific type gives a machine solicitation rate. We assume that the maximum level of machine solicitation is
U = 2.
The following analytical expression of the total maintenance cost includes preventive and corrective maintenance actions:
where
and
which represents the expected function of failures number during the horizon
depending on the disassembly-assembly plan. The failure rate
is cumulative at the beginning of each interval
k and is relative to the level of machine solicitation
is determined from the optimal disassembly-assembly plan. We give the failure rate in the interval
k as follows:
with:
where
is the failure rate for nominal conditions is equivalent to the failure rate of the maximum level of machine solicitation during all horizon
.Our work is based on Hajej et al. [
33] where they established an optimal integrated maintenance-production policy under subcontractor constraint, considering the influence of the production rate variation on the manufacturing system degradation and consequently on the maintenance strategy adopted.
Based on the work of Hajej et al. [
33], the average failure number expressed as follows:
We replace
by its expression and
T with
:
4. Numerical Results
In this section, we try to prove the profitability and the effectiveness of the proposed MRP system when considering remanufacturing parts. At first, we test an MRP model without taking the recovery of used/end-of-life products into account, which means that all raw materials (components of level 2) are ordered and new. Second, we interpret the results given as compared with our proposed model. Based on the calculation of the expected cost in Ben-Ammar et al. [
24], we calculate
to minimize the expected total cost (details in the
Appendix). Then, we determine the disassembly-assembly plan during the system’s work horizon. From this disassembly-assembly plan, we find the vector of the system solicitation levels for each period and type of disassembly or assembly process. We calculate the optimal preventive maintenance number
We consider a finished product with the two-level assembly structure described in
Figure 5.
The data noted below are used for the computation:
Table 3 and
Table 4 illustrate the different values of the unit inventory holding costs for components
and
in monetary units.
The other values of our example are:
,
.
Table 5 and
Table 6 represent respectively the remanufacturing and purchasing costs of components
.
The procurement time of the different components (used/end-of-life component, disassembled component, components at level 2 and assembled component at level 1) follow a uniform probability law.
For the calculation of the expected average number of failures, we adopt for a Weibull distribution with two parameters . This distribution is generally used to describe the system’s failure rate in industrial areas.
4.1. Optimal Ordered Date of Used/End-of-Life Product and Release Dates Plan of Components at Level 2
To calculate the optimal ordered date of used/end-of-life product and the vector of release dates for components at level 2, we implement a genetic algorithm with a simulation method. These components at level 2 must be ordered following the destruction of certain components after the disassembly process. The algorithm is implemented within JAVA language to resolve
.
Figure 6 illustrates the proposed algorithm. Indeed, the choice of the genetic algorithm is decided upon from its widespread use to solve discrete optimization problems, which proves that it is among the most used algorithms [
34] considering that is a relevant, accurate and fast algorithm. The objective of using the classical simulation method is to evaluate all possible values of the decision variables one by one. Indeed, our search space for the decision variable
is very reasonable, which helps us to calculate the total expected cost for each value, whereas the decision vector
depends on
and other deterministic and stochastic parameters of the system, which seems impossible to list all the possible combinations of the vector. This vector depends too on the number of components to be ordered. Therefore, as indicated in the flowchart of
Figure 6, we initialize at first the set of optimal solutions founded for each value of
. Then, we generate a random vector that represents the state of each component of the disassembled product. As mentioned above in
Figure 5, the proposed product consists of 10 parts at level 2, so among its components, we have remanufacturable and unremanufacturable parts. Then, we apply the genetic algorithm steps (selection, crossing, mutation) to find the best vector
The parameters of the genetic algorithms are noted in
Figure 6. In the second step, to calculate the optimal number of preventive maintenance actions, we implement the mathematical function
CM within the “Wolfram Mathematica” tool.
The results found are shown in tables below (
Table 7 and
Table 8). In
Table 7, the disassembly process is not considered; all components of level 2 are ordered from suppliers (external resources). On the other hand,
Table 8 shows results of the proposed model. The number of remanufactured components is generated randomly. As shown, the vector of dates of new components ordered from external suppliers consists of five dates, which proves that the disassembly of the used/end-of-life product has given us five remanufacturable components. Indeed, the optimal ordered date of the used/end-of-life product founded with the proposed algorithm is given in the table. We also observe the difference between costs, which proves the advantage of integrating the disassembly and remanufacturing process.
4.2. Optimal Maintenance Plan
In this sub-section, we try to find the optimal preventive maintenance plan for the disassembly-assembly system. We adopt the sequential strategy whose objective is to integrate the system plan obtained from the optimal dates calculated above (
Table 8). This preventive maintenance plan is determined by the disassembly-assembly plan. Based on dates founded in
Table 8, we built the disassembly-assembly system plan as shown in
Figure 7. According to the different levels of the system solicitation presented in
Table 1 and
Table 2 for different types, we suppose that our disassembly and assembly operations types vary between type 1 and type 5. In addition, depending on the different disassembly and assembly process and types of our example,
Figure 7 present levels of system solicitation on a time scale with the disassembly-assembly schedule planning. We suppose that the used/end-of-life disassembly operation solicits the machine at
and the solicitation rate of the finished product assembly at
.
The mathematical function
takes into consideration the influence of the variation of system solicitation on the disassembly-assembly system. These solicitation rates presented in
Figure 7 represents the vector
. This vector is used to calculate the average number of failures
.
To find the optimal value of preventive maintenance number, we implement the maintenance cost function
within the “Wolfram Mathematica” tool. We try to vary the variable
N from 0 to
H-1 and we calculate, at each value of
N, the total maintenance cost using
. As shown in
Figure 8, the curve has a convex form. The value
represents the number of preventive maintenance with minimum cost 7109.4 mu.
4.3. The Influence of Ordering Cost of the Used/End-of-Life Product on the System Profitability
In this subsection, we discuss the impact of ordering cost of the used/end-of-life product on system profitability, at a time when we are talking about environmental management, ecological responsibility on the future of end-of-life manufactured products, environmental certification, and finally, when many companies are looking after techniques and processes for re-evaluating end-of-life products. Besides the interest on the environment, entrepreneurs are also interested in the profitability of these processes and techniques. In this subsection, we are looking for the profitability of disassembly and remanufacturing processes compared to manufacturing process from new items. In fact, generally, the purchase of used/end-of-life products are cheaper than a new product. We are looking for the influence of the purchase cost of the used/end-of-life product.
Figure 9 shows a sensitivity study of the expected total cost
as a function of the purchase cost of used/end-of-life product. The chosen values of costs reflect the real case study (Assumption 7). We tried to find the ordering cost threshold of the used/end-of-life product to avoid falling into the profit-loss risk.
As shown in
Figure 9, using the values of the chosen example, we try to vary the ordering cost of the used/end-of-life product from 0 mu to 250 mu. If we exceed 230 mu, the expected total cost of our model exceeds the expected total cost without manufactured components of level 2.
In that case, the proposed system with remanufacturing parts becomes unprofitable. Consequently, by integrating the remanufacturing process, to have an effective cost, even as the cost of remanufacturing is very low compared to the purchase cost of the new products, it is inevitable to take into consideration the cost of used end-of-life product compared to total cost of new components.
4.4. The Influence of the Used/end-of-life Quality on the System Profitability
In this sub-section, we study the influence of the used/end-of-life quality on system profitability; in other words, the influence of the percentage of remanufactured parts. The number of remanufactured items is usually uncertain, and it will not be known until the end of the disassembly process. In fact, we cannot detect the quality of disassembled items from the arrival of the returned product. This uncertainty may have caused profit-loss risk compared to the expected cost of manufacturing a product from new items. In
Figure 10, a sensitivity study is illustrated by the expected total cost
as a function of the percentage of remanufactured items. The percentage of remanufactured items is adjusted to see their influence on the total expected cost. To a percentage found, a profit-loss risk will be discussed.
As shown in
Figure 10, increasing the percentage of remanufacturing items from 0 to 90%
under 10% of remanufactured parts, the expected total cost exceeds the expected cost with new items. Therefore, the proposed system with remanufacturing parts becomes unprofitable compared to the manufactured system with the purchase of new items in level 2. As a result, it is crucial to put requirements and standards for the quality of the used end-of-life product. Indeed, the remanufacturing process requires us to add used purchase cost of an end-of-life product, remanufacturing cost of disassembled parts, cost of storage of disassembled components, system maintenance cost, etc. Therefore, it is the handling of the used/end-of-life product with poor quality that causes a very low percentage of remanufactured items. The integration of a remanufacturing process becomes unprofitable and can even engender profit losses.
4.5. The Influence of the Corrective Maintenance Periods on the Lead-Time of the Finished Product
Many uncertainties can disturb the manufacturing process. In the implementation of the maintenance plan realized above, repair periods are negligible. This has no influence on disassembly-assembly plan. In a real industrial case, some failures take important periods, which disturb the manufacturing process and cause a delay in the delivery. Thus, decision makers must consider this repair time in planning before starting the order processing. In our work, we discuss the period that can take corrective maintenance actions according to the optimal plan implemented by the manager.
We suppose a threshold period noted
(
Figure 11). This period represents the time that corrective maintenance actions can take it.
As described above, we propose to study the system in a sequential way. From
Figure 11, we first implement the plan of disassembly and assembly which is determined from the optimal ordered date of used/end-of-life product and the vector of release dates for components at level 2. However, this optimal plan does not take into consideration the system repair time that causes the machine breakdowns, therefore the delay in disassembly and assembly processes. In fact, according to the representation of the optimal disassembly-assembly plan (
Figure 11), the product can be delivered to the customer at
.
Using this plan, we calculated the average number of random failures during the working horizon.
In the example illustrated above, we found that the average number of failures . Therefore, if we consider the same plan already calculated, we can note that .
Therefore, if we can quantify each repair time, we can say that each failure should not exceed 1/6 of a period: .
Beyond this period, the decision makers must re-implement the finished product manufacturing plan. Otherwise, a profit-loss risk must be discussed.
5. Conclusions
In a volatile industrial market, decision-makers find it difficult to choose the right approach that suits them to satisfy customers on time with a good-quality product in a first instance and with the lowest cost ahead of their competitors in the second instance. In this context, an optimization of disassembly, remanufacturing and assembly planning system is considered. The proposed system is composed of a single machine for disassembly and assembly operations, inventory of the disassembled components issued from the disassembly process of used/end-of-life product, inventory of components at level 2 and 1 and inventory of the finished product. The disassembly-assembly system is subject to random failures during the working horizon, which leads us to integrate a maintenance policy to deal with the unavailability of the system.
A mathematical model is proposed to formulate the problem. An optimization algorithm based on genetic algorithm and simulation is developed to find the optimal ordered date of used/end-of-life product and the vector of release dates for components at level 2. Using the optimization algorithm, numerical results reveal that the proposed model is more cost effective, under certain costs, than a simple MRP with new items. The costs that can influence the system are the ordering cost of the used/end-of-life product, and purchase costs of new components and remanufacturing costs. For this reason, the influence of ordering cost of the used end-of-life product on system profitability as well as the influence of the quality of the used/end-of-life product are studied. Finally, the risk of considering periods of corrective maintenance is investigated. This evaluation shows that the repair periods influence the disassembly-assembly plan. Beyond a certain threshold time, corrective maintenance actions must be integrated into the manufacturing planning. Otherwise, the profit-loss risk due to delay caused by the repair periods must be evaluated.
This research can serve as a decision support to manage certain product manufacturers. This approach can be used in many examples of product manufacturing. In fact, diverse products consist of expensive components as well as limited natural resources. In that case, remanufactured components become more cost effective than buying new ones. A lot of companies in the automotive and aeronautic sectors are looking to improve their productivity by integrating recovery and remanufacturing concepts in their process. Indeed, incorporating such an approach helps a manager reduce costs if components are very expensive compared to buying and processing a used one. In addition, this proposed method helps decision makers to know when to order a used product to exploit its different remanufacturable parts and when ordering new components that are missing in order to begin customer request processing. This method has been proved by the experiments proposed in this paper to encourage concerned researchers to develop it.
For future research, we will consider a multi-level disassembly-assembly product and a new maintenance policy. We will look to take other uncertainties, which can disrupt the supply chain.