Next Article in Journal
Bone: An Outstanding Composite Material
Next Article in Special Issue
Special Issue: “The Planning and Scheduling of Manufacturing Systems”
Previous Article in Journal
A Polynomial Inversion-Based Fast Time-Delay Estimation Method for Wideband Large-Scale Antenna Systems
Previous Article in Special Issue
Production Planning Problem of a Two-Level Supply Chain with Production-Time-Dependent Products
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Unreliable Subcontracting on Production and Maintenance Planning Considering Quality Decline

by
Héctor Rivera-Gómez
1,*,
Joselito Medina-Marin
1,
Francisca Santana-Robles
2,
Oscar Montaño-Arango
1,
Irving Barragán-Vite
1 and
Gabriel Cisneros-Flores
1
1
Academic Area of Engineering, Autonomous University of Hidalgo, Pachuca-Tulancingo Road km. 4.5, City of Knowledge, Mineral de la Reforma 42184, Mexico
2
Higher Education School Ciudad Sahagun, Autonomous University of Hidalgo, Sahagun-Otumba Road, Industrial Zone, Sahagun City 43998, Mexico
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(7), 3379; https://doi.org/10.3390/app12073379
Submission received: 3 March 2022 / Revised: 21 March 2022 / Accepted: 22 March 2022 / Published: 26 March 2022
(This article belongs to the Special Issue Planning and Scheduling of Manufacturing Systems)

Abstract

:
Manufacturing systems face several disturbances during production, such as sudden failures, defects, and unreliable subcontractors that reduce their production capacity. Currently, subcontracting represents an efficient alternative to support production decisions. The novelty of the study was the development of a new integrated model that properly coordinates production, subcontracting, and maintenances strategies in the context of stochastic uncertainty, quality deterioration, and random subcontracting availability. Such a set of characteristics has not been addressed before in the literature. A simulation–optimization approach was proposed to address such a stochastic model. A numerical case study was performed as an illustration of the approach and a comprehensive sensitivity analysis was performed to analyze the impact of several costs. Furthermore, the effect of the availability of the subcontractor and the producer was analyzed. The main finding of the study showed that the integrated model led to significant economic cost savings compared to other approaches that address such policies in isolation. The results also indicated that quality deterioration had a strong impact on the subcontracting rate and that the proposed joint control policy adequately coordinated these three key functions. The level of subcontracting participation was directly defined by its availability and the subcontracting cost.

1. Introduction

This introduction section is divided into three subsections. First, we highlight the motivation of the research, we present the reasons that justify the importance of the developed model, and we identify the industrial needs that motivated this study. In the second subsection, an extensive review of the literature is conducted, where we identify the recent research trends in the field of manufacturing systems. Moreover, the research gaps are identified. In the third subsection, the contribution of this study is presented, where we discuss the research gaps that we aimed to fill with the proposed integrated model. Additionally, we present the general resolution approach adopted in this study.

1.1. Motivation

Nowadays, industrial companies face several issues that jeopardize their economic performance. It is common in real production to have random failures, ineffective maintenance strategies, quality decline, etc., which have a strong impact on the general performance of such production systems. Therefore effective production and maintenance strategies are key to mitigating the effects of such random disturbances. Additionally, a third key function defined by product quality has become a critical dimension in recent years and has been gradually incorporated in the field of production systems. Since the determination of effective quality strategies is strongly related to production and maintenance, recently, the development of integrated models incorporating the three key functions of production, quality, and maintenance has provided superior results using traditional models that address these functions separately. However, more research is needed in this domain since there are still several factors that have not been included in several models and such factors have a negative effect on the system’s performance.
For example, gradual deterioration is a very common phenomenon with negative effects in automotive companies, semiconductor industries, manufacturing, etc. The production cost and the performance of production systems depend on the level of system deterioration. In manufacturing, a deterioration process certainly leads to short tool life, frequent setups, and increases in the total cost. Nevertheless, despite the negative effects of the deterioration process, there are still gaps that must be studied in the context of the determination of efficient production, quality, and maintenance strategies. In this uncertain context, companies need effective countermeasures to reach their economic goals; in recent years, subcontracting represents an attractive solution that is extensively used in production. Traditional models, including subcontracting, were developed based on simplistic assumptions that are not representative of the reality of industrial companies. For instance, in many studies, it was assumed that the subcontractor was always available. However, in real production, it is common to observe frequent random delays from the subcontractor. Indeed, stochastic subcontracting availability has a strong effect on the joint control policy and needs to be studied in detail.
From the above paragraphs, it is clear that more research is necessary since there is a need to determine advanced scheduling methods for effective production subcontracting and maintenance planning considering quality deterioration with the aim to foster the profitability of companies.

1.2. Literature Review

When analyzing the literature on the integration of production and maintenance strategies, we observed that several works were conducted, which makes sense given that current production systems have become more complex. For instance, Kouedeu et al. [1] studied a hierarchical decision-making process, where, in the first level, they determined the mean time to failure, and in the second level, they determined the production and maintenance rates. They assumed imperfect maintenance actions. Later, Dellagi et al. [2] studied production and maintenance planning for a system that satisfied random demand by considering a required service level. They considered a periodic preventive maintenance option, where, at failure, a minimal amount of repair was conducted. A subsequent work conducted by Fitouhi et al. [3] treated a production system that could degrade into several discrete states, which indicated different performance states. They analyzed the trade-off of a preventive maintenance policy by considering their impact on the production rate and the buffer. Polotski et al. [4] considered a hybrid system, where the machine was capable of manufacturing activities, and remanufacturing, where returned products were used in production. In their study, the system was prone to deterioration and they determined production, remanufacturing, and preventive maintenance rates. Rivera-Gómez et al. [5] studied production and repair strategies for a system subject to deterioration. They modeled quality deterioration with an aging process and several operational states. They also determined the repair efficiency to partially restore the machine after failure. More recently, Hajej and Rezg [6] developed a lot-sizing model in terms of energy, service level, and capacity constraints. They determined maintenance actions by taking into account the impact of the production rates on the level of deterioration and the increase in the failure rate. Despite the relevance of the discussed studies, more research is needed since the complex trade-offs between quality issues and the determination of production maintenance strategies were disregarded. However, companies must manage these three key functions simultaneously for economic success.
There is a vast literature emerging that investigated the interaction of production, quality, and maintenance strategies. For example, Bouslah et al. [7] determined the size of production lots and the values of the parameters of an acceptance sampling plan. They also defined overhaul rates for a deteriorating system considering an outgoing quality constraint. Fakher et al. [8] integrated production planning and imperfect maintenance, where process inspections were used to determine the system conditions. Their system comprised multiple products and multiple machines that deteriorated over time, producing more defects. Cheng et al. [9] developed a joint strategy of production and quality control for a degrading system, where preventive maintenance was conducted after inspection if the deterioration level surpassed a threshold. The strategy of quality control was based on inspecting all the units. Abubakar et al. [10] suggested a plan for production and maintenance for a system prone to deterioration, where the quality of the products was monitored with statistical process control. Their model determined the sample size, the interval for sampling, and the limits of quality control. The deterioration degree of their system was influenced by the production rate and time. Ait-El-Cadi et al. [11] proposed a model for the simultaneous control of production planning and maintenance scheduling considering a dynamic sampling inspection. Their quality policy adjusted the sample size according to the degree of deterioration. They assumed that the deterioration process had effects on the failure rate and the failure intensity. Hajej et al. [12] highlighted the advantage of using a dynamic sampling strategy for quality control for cases of progressive deterioration. In their model, the rates of production and maintenance were determined for the case of a system subject to quality deterioration. They also considered a service level and an outgoing quality constraint. As can be noted from the above-mentioned papers, the emphasis in this literature is that integrating production, maintenance, and quality control is critical in modern production systems. However, there is a limited amount of literature that considered subcontracting strategies in the integration of these three key functions. Subcontracting was demonstrated to be an effective alternative for deteriorating systems that face progressive capacity reduction.
In the literature, several researchers developed integrated models to study the strong link between production and subcontracting planning, as in Assid et al. [13], who addressed production and subcontracting strategies for a multiple-production facility with different capacities. In the optimization, they included a customer satisfaction constraint. Ben-Salem et al. [14] analyzed production and subcontracting plans for a manufacturing system that generated harmful emissions to the environment. In their model, items were available from a subcontractor at a higher cost. Such a subcontractor possessed more efficient technology but it had random availability. Haoues et al. [15] studied a two-echelon supply chain comprising multiple outsourcers and multiple subcontractors. Their model served to determine optimal production rates, periods for maintenance, and the determination of appropriate outsourcing options and optimal outsourcing quantities. Ayed et al. [16] treated an industrial problem of a production system that satisfied random demand. In their model, the subcontractor was available with a stochastic service level. Their system was prone to degradation and such a process had an impact on its availability. Rivera-Gómez et al. [17] investigated the problem of a deteriorating system with effects on quality and reliability with the aim to determine production and subcontracting strategies; in this study, preventive maintenance could be conducted to mitigate the effects of the deterioration process. In Kammoun et al. [18], subcontracting was allowed in a dynamic lot-sizing problem for a system that satisfied random demand and had requirements in the service level and consumption of energy. In their maintenance strategy, they regrouped the maintenance plans of several machines into only one plan. From these studies, it is evident that in the pursuit of industrial competitiveness, subcontracting represents an alternative to increasing the service level. However, the discussed models have several drawbacks, for instance, the effect of progressive degradation on the strategies of production and subcontracting was rarely considered. Indeed, the influence of such degradation on the control policy must be assessed and that is the object of this study.
In the real-life industrial field, production systems experience progressive deterioration that reduces their production capacity. Great interest has been dedicated to the field of deteriorating systems, for instance, Boudhar et al. [19] analyzed the case of a system subject to deterioration. They determined a maintenance policy that defined the inspection dates and the quality of the spare parts to be installed in the machine during maintenance. Martinod et al. [20] sought to compare several maintenance strategies, such as the age-based and periodic type for a production system with multiple components, taking into account the fact that such components are subject to deterioration. Ouaret et al. [21] investigated a production system prone to an aging process. The effect of such deterioration was mainly reflected in the rate of failures and also in the defect rate. The system faced random demand and replacement was available to partially reduce the impact of deterioration. In another work, Polotski et al. [22] developed a production policy considering two machines, the first one uses raw materials for manufacturing and the second machine uses return products for remanufacturing. They considered that product demand and return rates varies in function of time due to seasonal market behavior. An analytical model was presented in Dellagi et al. [23] to determine a maintenance strategy and the management of spare parts. Their system was subject to an increasing deterioration rate. They considered the influence of varying the production rate and the degree of degradation of the system and integrated the optimal quantity of spare parts to order for conducting maintenance actions. Magnanini and Tolio [24] proposed the integration of production and maintenance strategies, where the control policy used buffer thresholds to regulate the production rate and they considered the dependency of the thresholds on the deterioration condition. Their policy also included switching points that activated preventive maintenance. In the industrial environment, deterioration has a strong relationship with production and maintenance policies, as highlighted in the presented studies. Nevertheless, currently, there is a lack of research regarding the interaction of deterioration with subcontracting and quality issues since the focus has moved to a system perspective, where high levels of coordination are needed to achieve a required balance between such strategies and thus provide economic success. Indeed, such interactions and trade-offs are complex issues to be managed.
Mathematical techniques and simulation expertise are useful alternatives for analyzing the arising complexity of manufacturing models. For instance, Hosseini and Tan [25] presented an analytical simulation approach to analyze the performance of a continuous production system; their method provided solutions considerably faster than solely evaluating the system dynamics at critical times. They determined time instances of the trajectory of the buffer, stock level dynamics, and changing flow rates. Guiras et al. [26] examined different optimization algorithms and simulations to determine optimal plans for production and maintenance. They also analyzed the integration of imperfect maintenance with the consideration of returned products. Abdolmaleki et al. [27] used a simulation model to manage the production and maintenance strategies of a transported material network system. They considered that the machines were subject to deterioration with failures and their machine increased its deterioration level with each repair conducted. Rivera-Gómez et al. [28] studied production and maintenance strategies in the context of deterioration. They included a dynamic sampling inspection policy to ensure the satisfaction of an outgoing quality constraint. Further simulation research was proposed by Assid et al. [29], who determined production and remanufacturing rates. They assumed that returns were categorized into two categories based on their quality condition and time for reprocessing. They also defined the switching policy for the remanufacturing modes of the facility. Ait-El-Cadi et al. [30] proposed a simulation model to analyze the effects of production and preventive maintenance rates of a production system that was prone to aging on reliability and quality. They included the case of imperfect inspection activities, where a non-negligible duration was assumed for inspection and rectification activities. Despite the relevance of the presented studies, more research is needed since these studies were based on simple and unrealistic assumptions that did not consider the influence of degradation on the law of control. A similar argument shows that such models also disregarded the interdependence between random subcontracting availability and control actions.
Additionally, recent works in the area of production management and maintenance planning applied innovative methods based on metaheuristics and artificial intelligence (AI). Fathollahi-Fard et al. [31] studied a production-scheduling system that integrated air transportation. They considered a capacitated transportation system with a time window and they did not permit idle time. Their model was solved through nature-inspired metaheuristics. Jian et al. [32] formulated an optimization model to design the assembly system configuration, where their model defined the subassembly planning and the task assignments for uncertain product evolution, achieving workload balancing. Villalonga et al. [33] proposed an automated decision-making process using a fuzzy inference system to predict the condition of the asset and appropriate re-scheduling the production of a cyber-physical system. Chen et al. [34] studied the economic dependence of equipment cost and modeled the relationship between the effective service age and the reliability of components through a Weibull distribution. They established a selective maintenance model of a complex system with multiple states by considering economic dependence. Fathollahi-Fard et al. [35] proposed a multi-objective mixed-integer linear model to minimize the total energy consumption due to production. They considered social factors related to job opportunities and lost working days. They developed a sustainable flow shop scheduling problem with the assumption of different production centers and technologies on machines that had a strong impact on environmental and social criteria. Gholizadeh et al. [36] presented an optimization model for the problem of flexible flowshop scheduling for a waste-to-energy system. They proposed a preventive maintenance policy to determine an optimal sequence for processing tasks and minimizing delays. They considered the work processing time as being uncertain. Villalonga et al. [37] introduced a cloud-to-edges-based approach for a cyber-physical system to increase its smartness and the autonomy for monitoring and controlling the behavior of such a system at the shop-flow level. They evaluated their proposed solution with a pilot line. As can be noted from the presented papers, the authors only focused on the production or the production–maintenance relationship. The interconnections between production–quality–subcontracting–maintenance were not considered. Furthermore, a decision maker must take into account the hidden costs of AI methods due to IT infrastructures investments, the salary of data science professionals, etc., given that, in some cases, it is more convenient in terms of cost to use a traditional production–maintenance approach rather than expensive AI solutions, as indicated in Florian et al. [38].
To summarize, Table 1 serves to identify the research gaps that the present study aimed to fill. In this table, the rows classify the studies presented in the literature review in five main areas, and the columns present the key factors investigated in these studies.
The discussion of Table 1 is as follows: we analyzed in detail the studies of this table to identify their main contribution and define the key parameters presented in the columns. We note that the key parameters presented in the columns of Table 1 follow the evolution of the field of manufacturing systems. For instance, from the presented studies, it is clear that the first research trend focused on the development of production–maintenance strategies, as observed in the studies of Kouedeu et al. [1], Dellagi et al. [2], Fitouhi et al. [3], etc. Then, since 2013, the quality trend appeared, and thus, several studies included quality issues in their formulation, as in Bouslah et al. [7], Fakher et al. [8], Cheng et al. [9], etc. They aimed to jointly study production–quality–maintenance strategies. We note from the literature review that subcontracting strategies were only considered with production–maintenance models, as in Assid et al. [13], Ben-Salem et al. [14], Haoues et al. [15], etc. It is apparent to see that quality control was disregarded in such subcontracting models. Another common industrial phenomenon that has attracted the attention of several researchers is deterioration, as in the studies of Boudhar et al. [19], Martinod et al. [20], Ouaret et al. [21], etc. However, we observed that deterioration, mainly quality deterioration, was only studied with production–maintenance strategies. The influence of quality control or subcontracting has not been addressed in deterioration models. The same observation applied to the studies that adopted a simulation approach, such as Hosseini and Tan [25], Guiras et al. [26], Abdolmaleki et al. [27], etc. Therefore, based on this discussion, we concluded that there is a current need to develop an integrated model that considers the key functions of production, subcontracting, quality, and maintenance in a stochastic and quality deteriorating context.

1.3. Contribution

To summarize, the contribution of this research is the deep analysis of an integrated model that took into account the main functions of production planning, quality control, and maintenance, where unreliable subcontracting strategies were available to support demand satisfaction. The proposed model defined appropriate countermeasures by considering the level of inventory, the evolution of defectives, and the dynamics of failures and repairs, and a continuous-time stochastic process was used to model the underlying dynamics of the production system. This research aimed to extend previous results of the conventional system presented in Assid et al. [13] and mainly in the studies of Rivera-Gómez et al. [17,18,19,20,21,22,23,24,25,26,27,28], which investigated the effective coordination of production and subcontracting policies in the case of progressive quality decline. The mathematical complexity and the stochastic behavior of the proposed production system imposed several challenges, which were overcome with the use of a methodology based on several techniques, such as a semi-Markov process, stochastic optimal control, and simulation-based optimization. In our view, this research has three major contributions that make a distinction concerning the existing literature:
(i)
The influence of quality decline on the production, subcontracting, and maintenance control rule in a dynamic and stochastic context, which remain important issues in the field of manufacturing systems, was studied.
(ii)
Several production and maintenance models were developed by considering subcontracting in the context of deterioration, but they were limited to the case where subcontracting was always reliable. The key observation about the present model was that we aimed to extend such an assumption to study the impact of unreliable subcontracting.
(iii)
We proposed a novel stochastic optimal control model that is characterized by the use of the history of the deterioration process, leading to a semi-Markov model that effectively coordinated production and unreliable subcontracting.
The proposed model was extensively examined with a performance analysis of a numerical example and a comprehensive sensitivity analysis. The economic advantage of the proposed integrated model was highlighted in a comparative study. The next sections present the assumptions, the problem description, the control, and the simulation model used to reach the objectives proposed by the three contributions.
The rest of this paper is structured as follows: The assumptions; the problem description; and the definitions of the production, subcontracting, and maintenance policies are presented in Section 2. The resolution approach is detailed in Section 3. The simulation model used in the optimization phase is introduced and validated in Section 4. A numerical case study is presented in Section 5. An extensive sensitivity analysis and a comparative study are discussed in Section 6. An indifference curve that was used to analyze the trade-off between internal production and subcontracting is also presented in Section 6. Finally, conclusions are given in Section 7.

2. Problem Description and Assumptions

In order to clarify the control model under investigation, this section is divided into three subsections. First, the problem description is presented, where we describe in detail the set of features of the manufacturing system through a conceptual diagram. In the second subsection, the list of assumptions that support the proposed model is presented. In the third subsection, the control model is formulated analytically, where its dynamics, decision variables, control policies, and objective costs are identified.

2.1. Problem Description

Before beginning the description of the conceptual diagram presented in Figure 1, we first state that the proposed production system comprised an unreliable production unit that was prone to random failures and repairs and produced a single product. Assuming that the machine was subject to random failures and progressive deterioration had a strong impact on the rate of defects, defining a quality deterioration phenomenon as shown in Figure 1. Minimal repairs were available to the decision maker to overcome sudden failures. However, the level of deterioration of the machine was only rejuvenated with major maintenance, which returned the rate of defects to the initial conditions. The manufacturing system used inventory as protection against the shortages caused by the numerous breakdowns of the machine. This inventory level was a decision variable of the model that must be optimized. Another available strategy in the case of high levels of deterioration was subcontracting, which satisfied a fraction of the demand. It was assumed that the manufacturing system was not capable of satisfying product demand at high levels of deterioration. Furthermore, it was conjectured that the subcontractor was unreliable because it had several companies to serve simultaneously, and also because it was subject to technical problems that limited its production responsiveness. Nevertheless, the minimum amount of availability that the subcontractor must comply with was stipulated in a contract. The advantage of subcontracting is that it provides units free of defects. In this stochastic context, the main objective of the study was to develop an innovative integral model that simultaneously determined the optimal production–subcontracting plan and major maintenance strategies that minimized the total cost in the case of quality deterioration and unreliable subcontracting. In such an optimization, several costs were considered, such as inventory, backlog, defects, major maintenance, minimal repair, and subcontracting costs. The configuration of the whole production system is presented in Figure 1.

2.2. Assumptions

The stochastic system analyzed in this research was based on the following assumptions:
  • A1. In the beginning, a new production unit was operational and the influence of deterioration was minimal, as in Bouslak et al. [7].
  • A2. The demand was known and constant during the considered production horizon, as in Kouedeu et al. [1].
  • A3. The production unit experienced increasing wear with usage and it deteriorated. The deterioration was modeled with the use of an increasing function, which varied with the production rate, as in Ait-El-Cadi [11].
  • A4. The deterioration of the production system had a negative effect on the number of defects generated, as in Rivera-Gomez et al. [5].
  • A5. At failure, minimal maintenance was performed, where the condition of the machine was left in as-bad-as-old (ABAO) conditions, as in Assid et al. [29].
  • A6. Major maintenance was available and was a perfect repair that rejuvenated the production unit to as-good-as-new (AGAN) conditions, as in Rivera-Gomez et al. [28].
  • A7. Shortages may be observed if the time required to perform repairs and major maintenance took longer than the time to deplete the stock level used as protection, as in Polotski et al. [4].
  • A8. Subcontracting was available to support internal production and was an alternative to satisfy product demand, as in Ben-Salem et al. [14].
  • A9. The subcontractor was unreliable due to several disturbances caused by random failures, periods where it was busy serving other customers, other random events, etc.; as such, the subcontractor exhibited random availability. This assumption has not been used in the context of production–quality–maintenance joint optimization.
The notations used in the model are presented in Appendix A.

2.3. Control Model Formulation

To formulate the problem, for any specific time t , the production unit and the subcontractor modes were described using a stochastic process { α 1 ( t ) ,   t > 0 } Ω 1 = { 1 , 2 , 3 } and { α 2 ( t ) ,   t > 0 } Ω 2 = { 1 , 2 } , respectively. In particular, we note that when α 1 ( t ) = 1 , the production system was operational and generating defects at a rate of β . Further, when α 1 ( t ) = 2 , the production system was in failure mode, where minimal maintenance was performed, which left the machine in an ABAO condition before the occurrence of the failure. When α 1 ( t ) = 3 , major maintenance was undertaken, which completely mitigated the effects of the deterioration process since it restored the machine to an AGAN condition. Regarding the unavailability of the subcontractor, contrary to the assumption made in the literature, we assumed that the contractor was not always available to supply units due to technical disruptions; therefore, to make things simple, we proposed Markov dynamics for α 2 ( t ) . More precisely, when α 2 ( t ) = 1 , the subcontractor was available and they send flawless units to the production system, and when α 2 ( t ) = 2 , subcontracting was not available, and thus, it stopped its operation due to the conduction of maintenance activities.
The dynamics of the level of inventory was given by the next differential equation:
x ( t ) ( t ) = u 1 ( t ) 1 β   ( a ) + u 2 ( t ) d ,     x ( 0 ) = x 0      
where u 1 ( t ) is the production rate of the machine at time t , u 2 ( t ) is the subcontracting rate, β ( a ) is the rate of defects at age a , and d is the demand rate. The production rate at any time t satisfied 0 u 1 ( t ) u m a x , where u m a x is the maximum production rate for internal production. We assumed that the subcontractor satisfied only a fraction k 1 of the customer demand to be profitable for the company, as presented in the following equation:
u 2 ( t ) = k 1 · d  
For the current problem, we conjectured that the production unit experienced progressive deterioration. We modeled this condition using the age of the production unit as a key metric of its amount of deterioration. Furthermore, by assuming that the age of the machine at time t implied an increasing function calculated from its production rate, then the cumulative age was defined as follows:
a ( t ) ( t ) = η o · u 1 ( t )    
a ( T ) = 0
where η o is a given constant and T indicates the last restart time of the machine. Equation (3) uses the production rate to determine the age of the production unit. Furthermore, of significant importance is the fact that the model incorporated the influence of a degradation process in its formulation. In particular, we state that the defect rate was increased by the pace of deterioration of the machine, as in Colledani and Tolio [39], as defined in the next equation:
β ( a ) = b 0 + b 1 [ e η 1 a ( t ) η 2 ]  
where b 0 is the defect rate in an AGAN condition, b 1 defines the upper limit for the rate of defects, and η 1 and η 2 are positive constants. The production unit was unreliable and exhibited Markov dynamics due to random failures and repairs; hence, it was not always operational. The availability of the production unit could be calculated with the next equation:
π i p r o d · Q ( · ) = 0   y   i = 1 3 π i p r o d = 1  
where Q ( · ) refers to the transition rate matrix that is defined by the transition rates q i j between mode i and mode j of the production unit. Additionally, π i p r o d is the fraction of time that the production unit is in mode i . The resolution of Expression (6) defined the availability of the production unit in the operational mode; in this case, such a resolution led to
π 1 p r o d = 1 1 + q 12 q 21 + q 13 q 31
The major uncertainty in the model was that the subcontractor was unreliable, and thus, exhibited random unavailability due to sudden failures and maintenance activities. Hence, its availability in the operational mode could be defined with an expression similar to Equation (6); in this case, we used
π 1 s u b = 1 1 + q 12 s u b q 21 s u b  
where q i j s u b refers to the transition rate between mode i and mode j of the subcontractor. Due to the influence of the degradation process, the decision maker must design effective countermeasures to guarantee that the whole system was capable to meet the long-term demand, mainly in the scenario of extensive levels of degradation. In this case, when considering unreliable subcontracting, the total capacity constraint for the whole system was determined using
u m a x 1 β   ( a ) · π 1 p r o d + u 2 ( · ) · π 1 s u b d
In modeling terms, Equation (9) indicates that the level of production of the machine tended to decrease due to the effect of its deterioration, and thus, unreliable subcontracting was needed as an alternative to satisfy product demand.

2.3.1. Production and Subcontracting Strategies

In this section, we present a proposed extension to the hedging point policy, which included the contribution of unreliable subcontracting. The proposed policy used the age a ( t ) of the machine at time t to indicate the adequate time to trigger subcontracting participation. Equations (10) and (11) are convenient representations of the proposed production–subcontracting control policy.
If a ( t ) <   A s u b :
{ u 1 ( t ) = { u m a x 1 β ( a ) i f   x ( t ) < Z *   y     α 1 ( t ) = 1 d 1 β ( a ) i f   x ( t ) = Z *   y   α 1 ( t ) = 1 0 i f   x ( t ) > Z *   o r   α 1 ( t ) = 0   u 2   ( t ) = 0   x ( t ) ,   α 1 ( t ) ,   α 2 ( t )
If a ( t )   A s u b :
{ u 1 ( t ) = { u m a x i f   x ( t ) < Z *   a n d   α 1 ( t ) = 1   a n d   α 2 ( t ) = 0 M a x ( u m a x 1 β ( a ) u 2 * )   i f   x ( t ) = Z *   a n d   α 1 ( t ) = 1   a n d   α 2 ( t ) = 1 d 1 β ( a ) i f   x ( t ) = Z *   a n d   α 1 ( t ) = 1   a n d   α 2 ( t ) = 0 M a x ( ( d 1 β ( a ) u 2 *   ) , 0 )   i f   x ( t ) = Z *   a n d   α 1 ( t ) = 1   a n d   α 2 ( t ) = 1 0 i f   x ( t )   Z *   a n d   α 1 ( t ) = 0   a n d   α 2 ( t ) = 0 0 i f   x ( t ) > Z * u 2   ( t ) = {       u 2   *         i f   x ( t )   Z *   a n d   α 2 ( t ) = 1   a n d   α 1 ( t )       0       i f   x ( t )   Z *   a n d   α 2 ( t ) = 0   a n d   α 1 ( t )         0       i f   x ( t ) > Z *  
For practical purposes, A s u b is the critical age where subcontracting is triggered and Z * is the optimal production threshold. To an extent, Equation (10) represents a classical HPP strategy in the production unit, and it indicates that subcontracting is not requested when the age of the machine is less than A s u b since, in such a period, the machine is operational in acceptable conditions and the effects of deterioration are not so severe. By applying Equation (11), one could control the participation of the production unit and subcontracting when the age of the machine surpasses A s u b . In such an equation, internal production and subcontracting depend on their stochastic modes ( α 1 ( t ) , α 1 ( t ) ) , the inventory level x ( t ) , and the age of the machine a ( t ) .

2.3.2. Major Maintenance Strategy

The major maintenance strategy defined that this activity could be conducted once the machine reached the critical age A o . With the performance of this maintenance option, the machine was rejuvenated to an AGAN condition, eliminating the impairments caused by degradation. In this case, the major maintenance strategy was determined using
w * ( 1 , x , a ) = { 1 i f     a ( t ) A o 0 o t h e r w i s e          
To be more precise, A o was the critical age that triggered major maintenance activities. Summing up, the control parameters of the proposed policy were ( Z , A s u b , A o ) , which served to define the appropriate rate of production, subcontracting, and major maintenance ( u 1 * , u 2 * , w * ) . Such rates were defined taking into account the state variables of the system ( α 1 ,   α 2 , x , a ) . The strong interactions between these three strategies required jointly optimizing the control parameters to minimize the total cost since they were interlinked. Classical approaches that dissociate decisions turn out to be less efficient than integrated models, as is discussed in the next sections.

2.3.3. Optimization Problem

The main feature of the proposed model was the joint optimization of the control parameters ( Z , A s u b , A o ) , which were necessary to optimally control the subcontracting contribution, regulate the inventory level, and define the critical age to conduct major maintenance. These control parameters must be jointly optimized to minimize the total incurred cost T C . In our formulation, the total cost consisted of the sum of the inventory cost and the backlog cost I B ( t ) , the quality cost Q C ( t ) , and the maintenance cost M C ( t ) . Regarding the inventory–backlog cost, this component also comprised the production cost C p r o and the subcontracting cost C s u b , defined as
I B ( t ) = 1 T · 0 T ( C + x + ( t ) + C x ( t ) + C p r o u 1 ( t ) + C s u b u 2 ( t ) ) d t  
where x + = max ( 0 , x ) and x = max ( x , 0 ) and the constants C + and C serve to penalize inventory and shortages, respectively. The indicator Q C ( t ) represents the average cost of quality per unit of time in the period [ 0 , T ] ; it is determined by the cost of defects C d e f and is given by the following equation:
Q C ( t ) = 1 T · ( C d e f 0 T ( β ( t ) · d ) d t )
Further, M C ( t ) denotes the average maintenance cost per unit of time in the period [ 0 , T ] and it takes into account the cost of minimal repair and the cost of major maintenance; it is defined as:
M C ( t ) = 1 T · ( C R · N R ( t ) + C M · N M ( t ) )
where N R ( t ) and N M ( t ) define the number of minimal repairs and the number of major maintenances performed in the period [ 0 , T ] , respectively. Summing up, Equations (1)–(15) led to defining the next stochastic optimization model:
{ M i n   T C ( Z 0 , A s u b , A o ) = lim t ( I B ( t ) + Q C ( t ) + M C ( t ) ) s u b j e c t   t o E q u a t i o n s   ( 1 ) ( 9 ) ( d y n a m i c s   o f   i n v e n t o r y   a n d   q u a l i t y ) ( Z 0 , A s u b , A o ) 0 ( n o n n e g a t i v e   c o n s t r a i n t s )
Another important aspect was the fact that this optimization model implied a nonlinear and stochastic problem; for this reason, traditional mathematical methods could not be applied to solve this kind of problem since it was difficult to determine an analytical solution for the total cost. The difficulty was due to the random transitions, the Markov dynamics for internal production and subcontracting, and mainly because of the gradual growth in the number of defects caused by deterioration. One way to overcome these difficulties was with the implementation of a simulation–optimization approach. In the next section, we detail how we determined a solution to this model and determined the value of the parameters of the control rule in a reasonable amount of time and with acceptable accuracy.

3. Resolution Approach

In line with the previous section, we now introduce the proposed simulation–optimization approach, which combined simulation techniques with optimization methods with the aim to accurately reproduce the stochastic dynamics of the system under investigation. The optimization strategy and applied procedure began with a mathematical formulation of the manufacturing system under investigation; then, a simulation model was developed for the purpose of simulating and analyzing the system’s behavior. Statistical techniques were used for this analysis and optimization methods based on the response surface method were applied. The applied procedure provided remarkable results in previous studies, as in Rivera-Gómez et al. [17]. The approach comprised the following steps:
i.
Mathematical formulation: In this phase, an analytical model was formulated for the proposed system, where the decision variables were also identified, as well as the objective function to be minimized. Further in this step, the control parameters ( Z , A s u b , A o ) of the proposed control strategy were determined. Equations (1)–(16) defined the mathematical formulation of the production system under study.
ii.
Simulation model: In this step, a simulation model was built that served to reproduce the stochastic behavior of the production system. The inputs of this simulation model were the control parameters ( Z , A s u b , A o ) that were defined in the mathematical formulation step. The output of the simulation model was the total incurred cost. In the following section, the simulation model developed is presented in detail.
iii.
Statistical analysis: In this step, the design of experiments technique was used to determine the significant variables and interactions that must be included in the optimization phase. Moreover, using the total cost reported from several simulation runs, a second-order regression model was determined for the total incurred cost. A design of experiments technique that included three independent variables ( Z , A s u b , A o ) was applied for the determination of this regression model.
iv.
Optimization: In this step, a domain for the three independent variables ( Z , A s u b , A o ) was defined to perform the optimization of the regression model. The optimization strategy applied in this study was based on the response surface technique. Then, simulation runs were conducted to have enough data regarding the total cost and to perform the optimization. The response surface defined the optimal values ( Z * , A s u b * , A o * ) and minimized the total cost.
v.
Verification and validation: Several activities are conducted to ensure that the applied procedure provided accurate and credible results. Hence, an operational validation was performed, several key performance measures of the simulation model were displayed graphically, and a parameter variability analysis was conducted in the sensitivity analysis section to assess the effect of the variation of several cost and system parameters. Further, statistical analysis was carried out through a confidence interval for the total cost. Finally, a comparative study was done to highlight the economic advantage of the proposed model.
Figure 2 presents the resolution approach used in the research.

4. Simulation Model Description

As mentioned earlier, a simulation model was developed with a mixture of discrete and continuous components, where the simulation software ARENA was used in this phase. The model had the objective of reproducing the flow of material and the logic defined in the production–subcontracting–maintenance control strategy of Equations (10)–(12). The inputs for the simulation were the control parameters ( Z , A s u b , A o ) . Furthermore, the differential Equations (1) and (3) for the inventory and age dynamics were constantly integrated through the Runge–Kutta–Fehlberg technique. In the simulation model, the random duration of the minimal repairs and the duration of the major maintenance were defined by the discrete elements. Furthermore, the model defined the periods of inactivity of the unreliable subcontractor. The continuous part of the model was incorporated with C++ subroutines that actualized the defect rate β ( · ) with Equation (5), and it updated the system age a ( · ) . This C++ file defined the level of inventory and backlog. In this simulation model, all the components worked simultaneously to replicate the interactions and the random events of the whole system. Furthermore, in this model, the production and subcontracting planning and the maintenance strategies were included with the use of Equations (10)–(12). Once the simulation ended, the model displayed indices to estimate the inventory cost I B ( t ) , the quality cost Q C ( t ) , and the maintenance cost M C ( t ) with the use of Equations (13)–(15), respectively. Figure 3 presents the flow diagram of the simulation model.

Validation of the Simulation Model

With the aim to guarantee the accuracy of the simulation results, validation was carried out using several relevant indicators of the system’s behavior. The dynamics of these key indicators are presented in Figure 4. As an illustration, this validation was performed during a time interval where the parameters of the control rule were set at the following values: Z   = 20, A s u b   = 130, and A o = 175. As presented in Figure 4, it is clear that at time t = 0, the production unit was in perfect condition and the inventory level increased steadily until it reached its optimal level (see arrow 1 in Figure 4g). Subsequently, the unit experienced some random failures, and at time t = 37, a decrease in the inventory level was observed because the machine suffered a failure (see arrow 2 in Figure 4g). The occurrence of random failures and the periods of inactivity of the machine are identified in Figure 4d (see arrow 3). At this point, the machine experienced progressive deterioration, and then, at time t = 227, the unit could no longer satisfy the product demand by itself; therefore, the subcontractor started to supply units to satisfy customer demand (see arrow 4 in Figure 4c). At time t = 230, the inventory level was maintained at its optimal level, where internal production and subcontractor were both needed (see arrow 5 in Figure 4g). Then, at time t = 250, the unit operated at the rate d / ( ( 1 β ( a ) ) to mitigate the increase in the defect rate (as indicated by arrow 6 in Figure 4e) because the level of deterioration was significantly high, as shown by arrow 7 in Figure 4b. Later, at time t = 287, the unit age surpassed the critical value A o = 175, which initiated the conduction of major maintenance activities (as in arrow 8 in Figure 4f). Before the conduction of major maintenance, the maximum level of deterioration of the machine was reached (see arrow 9 in Figure 4a). Moreover, once the conduction of major maintenance was terminated, arrow 10 in Figure 4g shows a decrease in the inventory level due to the performance of such maintenance. It is evident from the results of Figure 4 that after the end of the conduction of major maintenance, there was a significant decrease in the rate of defects (see arrow 11 in Figure 4a) since the major maintenance rejuvenated the system to its initial conditions and completely mitigated the harmful effects of deterioration. At this point, the inventory level reached its optimal level of Z = 20, indicating the beginning of a new deterioration cycle. Based on the analysis of the dynamics of Figure 4, the proposed control strategy was validated since it worked appropriately in the simulation model and it was observed that the role of the subcontractor was very important once the machine age surpassed the critical value A s u b and major maintenance was triggered if the age of the machine reached A o . Summing up, the operational validity of the simulation model is presented in Figure 4.

5. Statistical Analysis and Optimization

In the rest of this section, a numerical case study of the system under investigation was conducted, where the data used in the analysis satisfied the feasibility condition (Equation (9)) and such data were based on previous studies of Rivera-Gómez et al. [17]. Furthermore, by taking into consideration the validation of the simulation model from the previous section, the data used ensured that the obtained results were representative of the productive system. In this numerical case study, three independent variables were considered ( Z , A s u b , A o ) and the dependent variable was defined by the total incurred cost. Several simulation runs were conducted based on a factorial 33 design. For each combination of the independent variables ( Z , A s u b , A o ) , the factorial design 33 was replicated two times; this meant that a total of 33 × 2 = 54 simulation runs were needed for the analysis. The simulation time used in the analysis was defined in 100,000 time units to guarantee steady-state conditions. The input data for the parameters of the simulation model are presented in Table 2.
Table 3 shows the values of several cost categories used in the simulation.
To simplify the calculations, it should be noted that to ensure that A s u b < A o , this age was defined as a fraction of age A o , i.e., A s u b = k · A o , where k   ϵ   [ 0 ,   1 ] . The values of the independent variables ( Z ,   k , A o ) shown in Table 4 were selected via the observation of several off-line simulations.
In the following, taking into account the information in Table 2, Table 3 and Table 4, the simulation model generated the necessary output data to determine a second-order regression model for the total cost. Such a regression model was analyzed through the use of statistical software called STATGRAPHICS, where an ANOVA was performed. Then, the regression model for the total cost was the following:
T C ( Z ,   k , A o ) =      809.977 + 4.74576   ·   Z 3.70433   · A o 1705.16   ·   k + 0.043567   ·   Z 2      0.0161327 · Z · A o   7.5515 · Z · k + 0.00580246   · A o 2    + 3.67873   ·   A o ·   k + 1036.04   ·   k 2
The analysis of variance of the regression model is presented in Table 5, where it is evident that all the main factors, most of the quadratic effects, and most of the interactions were significant since the p-value indicators were less than 5%.
The contour plots of the total cost reported by the regression model are presented in Figure 5 in two-dimensional space.
The regression model presented in Equation (17) was optimized with a response surface; in this case, the optimal solution that minimized the total cost is presented in Table 6. The value of k * = 0.6707 implied an age to trigger subcontracting of A s u b * = 100.81. Table 5 presents the optimal values of the control parameters that must be applied in the production system to simultaneously control the production rate, subcontracting, and major maintenance with minimal cost.
In closing this section, fifty extra simulation replications were performed using the optimal values reported in Table 6, whereby the accuracy of the proposed approach was corroborated since the total cost of 34.33 was within the confidence interval [25.78, 42.79]; this interval was obtained with the 50 simulation replications.

6. Sensitivity Analysis

So far, we have discussed the proposed simulation–optimization approach to the production system under investigation. Now, this section presents the detailed results of the robustness analysis of the proposed model and contrasts the value of the control parameters and the total cost under several system scenarios. In Table 7, fourteen different configurations are presented; the shaded values in Table 7 were selected by decreasing and increasing the parameters of the base case shown in Table 3 by a certain percentage. Such variations were selected to facilitate the observation of their impact on the control parameters and the total cost.
Table 8 presents the sensitivities for the variations devised in Table 7.
The sensitivity analysis of the variation of the cost parameters produced the following observations:
Inventory cost variation: When the inventory cost C + increased (case ii), the proposed strategy reduced the inventory level Z * since the accumulation of inventory was more sanctioned. However, by increasing C + , it was observed that the inventory was reduced to almost half of the level of the base case and this had the consequence that the subcontractor supplied parts earlier to avoid shortages. Thus, the system decreased the age A s u b * . Moreover, the system opted to perform more major maintenance to limit the generation of more defects. For this reason, the critical parameter of major maintenance A o * decreased. By decreasing the inventory cost (case i), the opposite effects were observed as compared to the case when the inventory cost increased.
Backlog cost variation: As the cost parameter C increased (case iv), shortages had a stronger economic penalty, hence the inventory level Z * need to increase with the aim to protect the system against the risk of not being able to meet the demand. It was also observed that by increasing C , the age to perform major maintenance A o * increased to operate the machine for longer and avoid product shortages caused by the performance of major maintenance. In addition, as the production unit was operative for longer, more inventory accumulated; therefore, the subcontractor was required less and this increased the age A s u b * . Decreasing the cost of backlog (case iii) had the opposite effects to increasing C .
Defective cost variation: The increase in the defect cost C d e f (case vi) mainly caused major maintenance to be performed more frequently to renew the unit as quickly as possible and significantly reduce the quantity of defective items; this measure decreased the age A o * . Additionally, it was observed that when increasing C d e f , the system reduced the inventory level since the presence of defects was more heavily penalized. This decision was complemented by decreasing the critical age A s u b * , which implied that the subcontractor must supply units at an earlier age to mitigate the presence of defects. The opposite effects were observed when the cost of defects decreased (case v).
Production cost variation: As the production cost C p r o increased (case viii), it was observed that major maintenance was performed at an earlier age as an attempt to ensure that the machine produced only flawless units; hence, the unit was functional for a shorter period, decreasing A o * . Moreover, as C p r o increased, it was evident that the inventory level Z * decreased because internal production was more costly. Moreover, as the production cost increased, it led to reducing the critical age of utilizing the subcontractor A s u b * because, in this context, it was more convenient for the company to have the subcontractor supply parts at an earlier age. As the production cost decreased, the opposite effects were observed (case vii).
Variation in the subcontracting cost: It was noted that an increase in the subcontracting cost C s u b (case x) had the logical effect of the major maintenance being performed with more frequency to rejuvenate the system and limit the presence of shortages. Likewise, the increase in the subcontractor’s cost C s u b caused that the system to increase the inventory level Z * because it was more profitable to opt for internal production than to use the subcontractor. Furthermore, it was obvious that by increasing C s u b , the critical age of subcontracting A s u b * increased due to the high cost of external supplies. Reducing the cost of subcontracting had the opposite effect of increasing C s u b (case ix).
Variation in the minimum repair cost: The reduction in the minimum repair cost C R (case xii) triggered the conduction of more major maintenance since this action kept the machine in a better condition, decreasing the presence of defects; this is the reason that explained the decrease in the critical age for major maintenance A o * . Furthermore, it is shown in Table 7 that with the decrease of A o * , the unit was functional for less time; therefore, less inventory accumulated, reducing Z * . On the other hand, the critical age of subcontracting A s u b * was reduced to guarantee the product availability due to the reduction in the inventory level Z * . By decreasing the minimum repair cost, the opposite effects were observed (case xi).
Variation in the cost of major maintenance: From the results, it was noted that the increase in the major maintenance cost C M (case xiv) led to observing that this action was performed less frequently; therefore, the age A o * increased. As C M increased, the system reached greater levels of deterioration to account for the high cost of major maintenance. The objective of postponing this maintenance action was to let the system operate for longer as an attempt to increase the inventory level, and thus, increase the capacity to respond to customer demand; this increased the inventory level Z * . Moreover, as the machine was operational for longer periods, it was obvious that less subcontracting was required; therefore, the age A s u b * increased because it was not economical for the company to have the subcontractor supply units in advance. Reducing the cost of major maintenance had the opposite effects (case xiii).
Based on the results discussed in Table 8, it was validated that the proposed control policy worked efficiently and that the logic of the control strategy was consistent throughout the analysis.

6.1. Effect of the Variation of the Producer Availability and the Subcontractor Availability

As shown in the following paragraphs, the main intention of this section was to present the analysis results regarding the influence of varying the subcontractor and producer availability. To make things simple, we defined the availability of the production system D p r o as the fraction of time when the manufacturing unit was available for production, as denoted previously in Equation (7). Additionally, to develop a realistic production model, it was assumed that the subcontractor was unreliable due to random disturbances. Then, the fraction of time where the subcontractor was available was denoted by D s u b , as noted earlier in Equation (8). Table 9 presents the effects of the availability of the producer and the availability of the subcontractor on the parameters of the control strategy and also on the total cost.
From the results presented in Table 9, it was observed that when the availability of the subcontractor D s u b was low (case xv), more protection was needed against shortages, hence the inventory level Z * increased. Moreover, in this context, it was required to operate the system for longer to compensate for the low availability of the subcontractor, thus major maintenance must be postponed; this measure led to an increased age A o * . With the low availability of the subcontractor, it was normal that its participation was delayed, increasing the age A s u b * . When the availability of the subcontractor increased (case xvi), the opposite effects were observed.
In the case where the availability of the producer D p r o reduced (case xvii), the system required more protection against random failures; thus, the inventory level Z * increased. Major maintenance was postponed to allow the unit to operate for longer without interruptions and to increase the capacity to satisfy customer demand, thus increasing the age A o * . Since the subcontractor’s cost was more expensive than in-house production, when the availability of the producer decreased, then it was preferable to continue with in-house production, thus increasing the age A s u b * . When the availability D p r o increased, the opposite effects were observed (case xviii).

6.2. Comparative Study

To ensure that the proposed joint control policy provided economic benefits over other existing approaches, in this section, a comparison of the total cost was conducted between the proposed policy and an alternative policy derived from typical assumptions used in the literature. The alternative policy was denoted as Policy-B and it did not consider a subcontracting strategy. Policy-B did not include the parameter A s u b in the decision variables. The control parameters used in Policy-B were therefore only the decisions regarding production and maintenance planning ( Z , A o ) . A similar optimization approach to that used in the previous section was applied with the aim to define the optimal control factors ( Z , A o ) for Policy-B. The definitions of the control strategies included in the study were as follows:
Policy-A referred to the control policy ( Z , A s u b , A o ) developed in this study and consisted of the simultaneous optimization of the inventory level, the determination of the most adequate age to use subcontracting, and the determination of the age to conduct major maintenance activities.
Policy-B was a generalization of our policy and was intended to examine the impact of subcontracting in the optimization. In this policy, the control parameters were reduced to only the production and maintenance decisions ( Z , A o ) , disregarding the subcontracting decisions in the optimization. In Policy-B, subcontracting was not taken into account in the control policy and the optimization.
Table 10 presents the results of the comparative study conducted.
An important consideration regarding Table 10 is that the cost difference between the two policies was 32.14%, where this gap was explained by the fact that Policy-B used larger inventory levels to guarantee that product demand could be satisfied and this countermeasure considerably increased the total cost. In addition, the production unit remained operational for longer periods in Policy-B before major maintenance was performed; this action led to an increase in the total cost since if the machine operated for longer without adequate maintenance, then more defects were generated due to deterioration, increasing the total cost. This comparison gave us a measure of the utility of the proposed joint control policy.

6.3. Indifference Curve to Determine the Subcontracting Participation

Based on the results discussed in the last section, it was observed that subcontracting was an effective countermeasure to improve the level of demand satisfaction and minimize costs. The key idea of the proposed model was to optimize in-house production and compensate with the use of the subcontractor if needed when the deterioration process reduced the production capacity of the production unit. In this section, we present further analysis results that served to properly select the participation of subcontracting based on two key characteristics, namely, the cost of subcontracting C s u b and its availability D s u b . The optimization procedure described in the previous section was used to determine under which conditions it was more advisable to use subcontracting than in-house production. The approach consisted of considering the cost of subcontracting C s u b as an independent variable that was optimized with the rest of the control parameters ( Z , A s u b , A o ) of Policy-A. In addition, the total cost was defined as the dependent variable of the optimization. In this way, different values of the subcontractor’s availability were analyzed to define until what amount of cost C s u b it was still profitable to use subcontracting. By analyzing different subcontracting costs and availabilities, the inference curve of Figure 6 was obtained.
The interpretation of Figure 6 is in terms of two main zones. The first zone (see indicator 1 in Figure 6) is where the system recommended using a combination of subcontracting and in-house production to reduce the total incurred cost. Furthermore, zone 2 (see indicator 2 in Figure 6) is where in-house production was the most cost-effective and profitable option. For example, when the subcontractor had availability of D s u b = 91 % , the manager could manage to pay a cost of C s u b = 16 to the subcontractor. Otherwise, there was no economic incentive in negotiating with the subcontractor. However, the company could pay more to the subcontractor, for instance, a cost of C s u b = 17.84 , only if the subcontractor guaranteed a higher availability of D s u b = 99 % . In this sense, it could be stated that if the availability of the subcontractor was high, a higher cost could be paid for its services.

7. Conclusions

The present research was devoted to the joint analysis of production, unreliable subcontracting, and perfect maintenance strategies for a manufacturing system that experienced a quality decline process. The studied production system was prone to random failures and increased its defect rate as it deteriorated. Thus, to ensure demand satisfaction, subcontracting was available at a higher cost as a countermeasure to cope with the effect of such deterioration. However, as observed in real production, the subcontractor was unreliable since it was not always available to supply products to the company due to random disturbances caused by, e.g., sudden failures or because it was occupied serving other customers.
The major outcome of this research was the consideration of the four key aspects of production, subcontracting, maintenance, and quality in an integrated model under a dynamic and stochastic context, where the subcontractor was unreliable and the system underwent deterioration. In this domain, no prior work in the literature has treated this set of characteristics. A control policy derived from the well-known HPP was proposed to properly coordinate internal production and subcontracting. From the results, it was observed that the proposed policy appropriately achieved the coordination of internal production and subcontracting for any possible case. Decisions on these strategies and maintenance were defined as a function of the level of inventory, the availability of the subcontractor, the availability of the producer, and the amount of degradation of the machine. A simulation–optimization approach was adopted to solve such a stochastic problem. The proposed approach involved developing a meta-model comprising a quantitative study based on the design of experiments and simulation analysis, where the control parameters were optimized through a response surface. The proposed approach was validated through a graphical analysis and the assessment of the variation of several costs; in all the analyzed cases, the proposed model managed to adequately control the production system. Moreover, the impact of the variation of the availability of the subcontractor and the producer were analyzed; it was found that more inventory was needed to cope with shortages when the subcontractor was less available and more maintenance was required to restore the production unit; in this context, the participation of the subcontractor was delayed. Through a comparative study, it was observed that the proposed policy reported considerable cost savings of 32.14% compared to other classical policies that disregarded subcontracting. Such a cost difference could even be more considerable if the cost C s u b or the availability of the subcontractor D s u b increased. Furthermore, an indifference curve was analyzed and it served to define the cases when only internal production was more profitable to the company and the cases when a combination of internal production and subcontracting was more cost-effective.
Regarding future research, we note that the maintenance strategy used to overhaul the system in the proposed model was a perfect repair. However, imperfect maintenance strategies are more realistic regarding model current production systems. The consideration of imperfect maintenance in the context of unreliable subcontracting will be the subject of future research. Additionally, the effect of random demand in the context of unreliable subcontracting can be included in the models.

Author Contributions

H.R.-G. and J.M.-M. developed the mathematical model, F.S.-R. and O.M.-A. performed the numerical approach, I.B.-V. conducted the sensitivity analysis, and G.C.-F. performed the comparative study. All authors read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

This section describes the major notations and symbols contained in this paper.
State variables:
x ( t ) Level of inventory at time t
a ( t ) Production unit age at time t
α 1 ( t ) Production unit mode
α 2 ( t ) Unreliable subcontractor mode
Model parameters:
d Constant demand rate
u m a x Maximal production rate of the unit
β ( a ) Defect rate at age a
k 1 The ratio of demand satisfied by the subcontractor
η o Positive constant used to model aging
b 0 Value of defectives at initial conditions
b 1 The upper limit for the increase in the rate of defects
η 1 Positive constant used to model the increase in defects
η 2 Positive constant used to model the increase in defects
q i j Transition rates from state i to state j of the production unit
q i j s u b Transition rates from state i to state j of the unreliable subcontractor
Cost parameters:
T C Total incurred cost
I B ( t ) Average inventory–backlog cost per time unit
Q C ( t ) Average quality cost per time unit
M C ( t ) Average maintenance cost per time unit
N R ( t ) Number of minimal repairs conducted in the simulation time interval
N M ( t ) Number of major maintenance conducted in the simulation time interval
C + Cost of goods in stock
C Cost of backlogs
C p r o Cost of production
C s u b Cost of subcontracting
C d e f Cost of defects
C R Cost of a minimal repair
C M Cost of major maintenance
Decision variables:
u 1 ( t ) Internal production rate at time t
u 2 ( t ) The unreliable subcontracting rate at time t
w ( t ) Major maintenance rate at time t
Control parameters:
Z Optimal inventory level
A s u b Critical age that controls the participation of subcontracting
A o Critical age to trigger major maintenance activities

References

  1. Kouedeu, A.F.; Kenné, J.P.; Songmene, V. Production, preventive and corrective maintenance planning in manufacturing systems under imperfect repairs. In Proceedings of the 3rd International Workshop on Dependable Control of Discrete Systems, Saarbruecken, Germany, 15–17 June 2011. [Google Scholar] [CrossRef] [Green Version]
  2. Dellagui, S.; Chelbi, A.; Trabelsi, W. Joint integrated production-maintenance policy with production plan smoothing through production rate control. J. Manuf. Syst. 2017, 42, 262–270. [Google Scholar] [CrossRef]
  3. Fitouhi, M.C.; Nourelfath, M.; Gershwin, S.B. Performance evaluation of a two-machine line with a finite buffer and condition-based maintenance. Reliab. Eng. Syst. Saf. 2017, 166, 61–72. [Google Scholar] [CrossRef]
  4. Polostki, V.; Kenné, J.P.; Gharbi, A. Joint production and maintenance optimization in flexible hybrid manufacturing-remanufacturing systems under age-dependent deterioration. Int. J. Prod. Econ. 2019, 216, 239–254. [Google Scholar] [CrossRef]
  5. Rivera-Gómez, H.; Lara, J.; Montaño-Arango, O.; Hernández-Gress, E.S.; Corona-Armenta, J.R.; Santana-Robles, F. Joint production and repair efficiency planning of a multiple deterioration system. Flex. Serv. Manuf. J. 2019, 31, 446–471. [Google Scholar] [CrossRef]
  6. Hajej, Z.; Rezg, N. An optimal integrated lot sizing and maintenance strategy for multi-machines system with energy consumption. Int. J. Prod. Res. 2020, 58, 4450–4470. [Google Scholar] [CrossRef]
  7. Bouslah, B.; Gharbi, A.; Pellerin, R. Integrated production, sampling quality control and maintenance of deteriorating production systems with AOQL constraint. Omega 2011, 61, 110–126. [Google Scholar] [CrossRef]
  8. Fakher, H.B.; Nourelfath, M.; Gendreau, M. A cost minimization model for joint production and maintenance planning under quality constraints. Int. J. Prod. Res. 2017, 55, 2163–2176. [Google Scholar] [CrossRef]
  9. Chen, G.Q.; Zhou, B.H.; Li, L. Integrated production, quality control and condition-based maintenance for imperfect production systems. Reliab. Eng. Syst. Saf. 2018, 175, 251–264. [Google Scholar] [CrossRef]
  10. Abubakar, A.S.; Nyoungue, A.; Hajej, Z. Integrated production, maintenance and control chart of supply chain management under quality constraint. In Proceedings of the 2020 IEEE International Conference on Industrial Engineering and Engineering Management (IEEM), Singapore, 14–17 December 2020. [Google Scholar] [CrossRef]
  11. Ait-El-Cadi, A.; Gharbi, A.; Dhouib, K.; Artiba, A. Integrated production, maintenance and quality control policy for unreliable manufacturing systems under dynamic inspection. Int. J. Prod. Econ. 2021, 236, 108140. [Google Scholar] [CrossRef]
  12. Hajej, Z.; Rezg, N.; Gharbi, A. Joint production preventive maintenance and dynamic inspection for a degrading manufacturing system. Int. J. Adv. Manuf. Technol. 2021, 112, 221–239. [Google Scholar] [CrossRef]
  13. Assid, M.; Gharbi, A.; Dhouib, K. Joint production and subcontracting planning of unreliable multi-facility multi-product production systems. Omega 2015, 50, 54–69. [Google Scholar] [CrossRef]
  14. Ben-Salem, A.; Gharbi, A.; Hajji, A. Production and uncertain green subcontracting control for an unreliable manufacturing system facing emissions. Int. J. Adv. Manuf. Technol. 2016, 83, 1787–1799. [Google Scholar] [CrossRef] [Green Version]
  15. Haoues, M.; Dahane, M.; Mouss, N.K. Outsourcing optimization in two-echelon supply chain network under integrated production-maintenance constraints. J. Intell. Manuf. 2019, 30, 701–725. [Google Scholar] [CrossRef]
  16. Ayed, S.; Hajej, Z.; Turki, S.; Rezg, N. FPA method for optimal production planning under availability/degradation machine and subcontracting constraint. Int. J. Prod. Res. 2017, 55, 2135–2148. [Google Scholar] [CrossRef]
  17. Rivera-Gómez, H.; Gharbi, A.; Kenné, J.P.; Montaño-Arango, O.; Hernández-Gress, E.S. Subcontracting strategies with production and maintenance policies for a manufacturing system subject to progressive deterioration. Int. J. Prod. Econ. 2018, 200, 103–118. [Google Scholar] [CrossRef]
  18. Kammoun, M.A.; Hajej, Z.; Rezg, N. Risk analysis for lot-sizing and maintenance optimization problem under energy constraint with subcontractor solution. Int. J. Qual. Reliab. 2020, 37, 873–904. [Google Scholar] [CrossRef]
  19. Boudhar, H.; Dahane, M.; Rezg, N. New dynamic heuristic for the optimization of opportunities to use new and remanufactured spare part in stochastic degradation context. J. Intell. Manuf. 2017, 28, 437–454. [Google Scholar] [CrossRef]
  20. Martinod, R.M.; Bistorin, O.; Castañeda, L.F.; Rezg, N. Maintenance policy optimization for multi-component systems considering degradation of components and imperfect maintenance actions. Comput. Ind. Eng. 2018, 124, 100–112. [Google Scholar] [CrossRef] [Green Version]
  21. Ouaret, S.; Kenné, J.P.; Gharbi, A. Production and replacement policies for a deteriorating manufacturing system under random demand and quality. Eur. J. Oper. Res. 2018, 264, 623–636. [Google Scholar] [CrossRef]
  22. Polotski, V.; Kenné, J.P.; Gharbi, A. Production control of hybrid manufacturing-remanufacturing systems under demand and return variations. Int. J. Prod. Res. 2019, 57, 100–123. [Google Scholar] [CrossRef]
  23. Dellagi, S.; Trabelsi, W.; Hajej, Z.; Rezg, N. Integrated maintenance/spare parts management for manufacturing system according to variable production rate impacting the system degradation. Concurr. Eng. 2020, 28, 72–84. [Google Scholar] [CrossRef]
  24. Magnanini, M.C.; Tolio, T. Switching-and hedging-point policy for preventive maintenance with degrading machines: Application to a two-machine line. Flex. Serv. Manuf. J. 2020, 32, 241–271. [Google Scholar] [CrossRef] [Green Version]
  25. Hosseini, B.; Tan, B. Simulation and optimization of continuous-flow production systems with a finite buffer by using mathematical programming. IISE Trans. 2017, 49, 255–267. [Google Scholar] [CrossRef]
  26. Guiras, Z.; Hajej, Z.; Rezg, N.; Dolgui, A. Comparative analysis of heuristic algorithms used for solving a production and maintenance planning problem (PMPP). Appl. Sci. 2018, 8, 1088. [Google Scholar] [CrossRef] [Green Version]
  27. Abdolmaleki, S.; Shokouhyar, S.; Afshar-Bakeshloo, M. A simulation based optimization approach to combine two hedging control policies in a transported degrading failure-prone manufacturing system. Transp. Manag. 2018, 1, 832. [Google Scholar] [CrossRef]
  28. Rivera-Gómez, H.; Gharbi, A.; Kenné, J.P.; Montaño-Arango, O. Joint optimization of production and maintenance strategies considering a dynamic sampling strategy for a deteriorating system. Comput. Ind. Eng. 2020, 140, 106273. [Google Scholar] [CrossRef]
  29. Assid, M.; Gharbi, A.; Hajji, A. Production planning and control of unreliable hybrid manufacturing-remanufacturing systems with quality-based categorization of returns. J. Clean. Prod. 2021, 312, 127800. [Google Scholar] [CrossRef]
  30. Ait-El-Cadi, A.; Gharbi, A.; Dhouib, K.; Artiba, A. Joint production and preventive maintenance controls for unreliable and imperfect manufacturing systems. J. Manuf. Syst. 2021, 58, 263–279. [Google Scholar] [CrossRef]
  31. Fathollahi-Fard, A.M.; Hajiaghaei-Keshteli, M. Integrated capacitated transportation and production scheduling problem in a fuzzy environment. Int. J. Ind. Eng. Prod. Res. 2018, 29, 197–211. [Google Scholar] [CrossRef]
  32. Jian, Z.; Wang, H.; Dulebenets, M.A.; Pasha, J. Assembly system configuration design for reconfigurability under uncertain production evolution. J. Manuf. Sci. Eng. 2019, 141, 071001. [Google Scholar] [CrossRef]
  33. Villalonga, A.; Negri, E.; Biscardo, G.; Castano, F.; Haber, R.E.; Fumagalli, L.; Macchi, M. A decision-making framework for dynamic scheduling of cyber-physical production system based on digital twins. Annu. Rev. Control 2021, 51, 357–373. [Google Scholar] [CrossRef]
  34. Chen, Z.; Zhang, L.; Tian, G.; Nasr, E.A. Economic maintenance planning of complex systems based on discrete artificial bee colony algorithm. IEEE Access 2020, 8, 108062–108071. [Google Scholar] [CrossRef]
  35. Fathollahi-Fard, A.M.; Woodward, L.; Akhrif, O. Sustainable distributed permutation flow-shop scheduling model based on a triple bottom line concept. J. Ind. Inf. Integr. 2021, 24, 100233. [Google Scholar] [CrossRef]
  36. Gholizadeh, H.; Fazlollahtabar, H.; Fathollahi-Fard, A.M.; Dulebenets, M.A. Preventive maintenance for the flexible flowshop scheduling under uncertainty: A waste-to-energy system. Environ. Sci. Pollut. Res. 2021. [Google Scholar] [CrossRef]
  37. Villalonga, A.; Beruvides, G.; Castaño, F.; Haber, R. Cloud-based industrial cyber-physical system for data-driven reasoning. A review and use case on an industry 4.0 pilot line. IEEE Trans. Ind. Inform. 2020, 16, 5975–5984. [Google Scholar] [CrossRef]
  38. Florian, E.; Sgarbossa, F.; Zennaro, I. Machine learning-based predictive maintenance: A cost-oriented model for implementation. Int. J. Prod. Econ. 2021, 236, 108114. [Google Scholar] [CrossRef]
  39. Colledani, M.; Tolio, T. Integrated analysis of quality and production logistics performance in manufacturing lines. Int. J. Prod. Res. 2011, 49, 485–518. [Google Scholar] [CrossRef]
Figure 1. Block diagram of the production system.
Figure 1. Block diagram of the production system.
Applsci 12 03379 g001
Figure 2. Resolution approach.
Figure 2. Resolution approach.
Applsci 12 03379 g002
Figure 3. Logic for the simulation.
Figure 3. Logic for the simulation.
Applsci 12 03379 g003
Figure 4. Simulation model validity.
Figure 4. Simulation model validity.
Applsci 12 03379 g004
Figure 5. Total cost contour plots.
Figure 5. Total cost contour plots.
Applsci 12 03379 g005
Figure 6. Indifference curve to determine the subcontractor participation.
Figure 6. Indifference curve to determine the subcontractor participation.
Applsci 12 03379 g006
Table 1. Literature review.
Table 1. Literature review.
Production PlanningQuality IssuesMaintenance PlanningSubcontracting StrategyUnreliable SubcontractingDeterioration Quality DeteriorationDevelopment of a New PolicyPerformance Metrics Analysis
(a) Production–maintenance strategies
Kouedeu et al. [1]
Dellagi et al. [2]
Fitouhi et al. [3]
Polotski et al. [4]
Rivera-Gómez et al. [5]
Hajej and Rezg [6]
(b) Production–quality–maintenance strategies
Bouslah et al. [7]
Fakher et al. [8]
Cheng et al. [9]
Abubakar et al. [10]
Ait-El-Cadi et al. [11]
Hajej et al. [12]
(c) Production and subcontracting planning
Assid et al. [13]
Ben-Salem et al. [14]
Haoues et al. [15]
Ayed et al. [16]
Rivera-Gómez et al. [17]
Kammoun et al. [18]
(d) Deteriorating systems
Boudhar et al. [19]
Martinod et al. [20]
Ouaret et al. [21]
Polotski et al. [22]
Dellagi et al. [23]
Magnanini and Tolio [24]
(e) Simulation approach
Hosseini and Tan [25]
Guiras et al. [26]
Abdolmaleki et al. [27]
Rivera-Gómez et al. [28]
Assid et al. [29]
Ait-El-Cadi et al. [30]
(f) Additional issues
Fathollahi-Fard et al. [31]
Jian et al. [32]
Villalonga et al. [33]
Chen et al. [34]
Fathollahi-Fard et al. [35]
Gholizadeh et al. [36]
Villalonga et al. [37]
The proposed model
Table 2. Data for the numerical case study.
Table 2. Data for the numerical case study.
Parameter: q 12 q 21 q 13 q 31 q 12 s u b q 21 s u b
Value:0.11.50.1550.250.3
Parameter: u m a x d k 1 η o b 0 b 1
Value:950.5 0.10.010.49
Parameter: η 1 η 2 k 1
Value:15 × 10−6.22.40.5
Table 3. Cost parameters.
Table 3. Cost parameters.
Parameter: C + C C R C M C d e f C p r o C s u b
Value:140100300020515
Table 4. Domains of the independent variables.
Table 4. Domains of the independent variables.
FactorLow LevelHigh LevelDescription
Z 550Inventory level
k 0.50.99Critical parameter to trigger subcontracting
A o 50300Critical age to trigger major maintenance
Table 5. ANOVA Table.
Table 5. ANOVA Table.
SourceSum of SquaresGlMedium SquareF-Ratiop-Value
A:Zo31,139.5131,139.54.090.0493
B:Ao218,6991218,69928.750.0000
C:k163,0111163,01121.430.0000
AA5837.515837.50.770.3859
AB49,409.4149,409.46.50.0145
AC41,588.6141,588.65.470.0241
BB98,638.4198,638.412.970.0008
BC304,6201304,62040.050.0000
CC46,408.7146,408.76.10.0176
Blocks115.4451115.4450.020.9025
Total error327,073 7606.36
Total (corr.)1.29 × 106
Table 6. Optimal control parameters.
Table 6. Optimal control parameters.
Z * k * A o * Total CostCross-Check Confidence Interval (95%)
Value31.50.6707150.3234.33[25.78, 42.79]
* show the optimal control.
Table 7. Cost variation.
Table 7. Cost variation.
Cost Variation
CasesPar. C + C C d e f C p r o C s u b C R C M
Base case-140205151003000
Case i C + 0.540205151003000
Case ii 240205151003000
Case iii C 110205151003000
Case iv 150205151003000
Case v C d e f 14055151003000
Case vi 140305151003000
Case vii C p 140201151003000
Case viii 1402010151003000
Case ix C s 14020551003000
Case x 140205301003000
Case xi C R 14020515303000
Case xii 140205152003000
Case xiii C M 140205151001000
Case xiv 140205151004500
Table 8. Sensitivities for the cost variations.
Table 8. Sensitivities for the cost variations.
Total
Optimal Control Parameters Cost
CasesPar. Z * k * A o * A s u b * TIC*Effect
Base Case-31.50.6707150.32100.819634.33Base of comparison
Case i C + 39.830.692154.84107.149317.05 Z * , Ao*, Asub*↑
Case ii 14.930.627140.3688.005756.86 Z * , Ao*, Asub*↓
Case iii C 16.580.6604132.7787.681366.10 Z * , Ao*, Asub*↓
Case iv 34.60.6718154.47103.772916.28 Z * , Ao*, Asub*↑
Case v C d e f 32.210.6471165.17106.881523.42 Z * , Ao*, Asub*↑
Case vi 310.6875139.8596.146940.97 Z * , Ao*, Asub*↓
Case vii C p r o 32.190.6711154.27103.530614.10 Z * , Ao*, Asub*↑
Case viii 30.670.6705145.4697.530959.43 Z * , Ao*, Asub*↓
Case ix C s u b 30.660.6441154.8899.758224.68 Z * , Ao*, Asub*↓
Case x 33.610.7095143.82102.040347.66 Z * , Ao*, Asub*↑
Case xi C R 31.530.6702150.7100.999127.96 Z * , Ao*, Asub*↑
Case xii 31.470.6715149.87100.637743.43 Z * , Ao*, Asub*↓
Case xiii C M 30.780.6985133.5993.312625.46 Z * , Ao*, Asub*↓
Case xiv 31.930.6543160.29104.877740.26 Z * , Ao*, Asub*↑
↑, and ↓ shows the effect of the variation of the cost parameter.
Table 9. Effect of the variation of the producer and subcontractor availability.
Table 9. Effect of the variation of the producer and subcontractor availability.
Case AvailabilityControl ParametersEffect
Par. D s u b   (%) D p r o
(%)
Z * k * A o * A s u b * Cost*
Base case 92.3093.7531.500.6707150.32100.819634.33Base Case
Case xv D s u b 89.0293.7532.490.6733151.58102.0634.53 Z * ↑, Ao*↑, Asub*↑
Case xvi 95.2393.7529.980.6710146.6298.3833.72 Z * ↓, Ao*↓, Asub*↓
Case xvii D p r o 92.3091.4633.920.6611153.68101.6014.59 Z * ↑, Ao*↑, Asub*↑
Case xviii 92.3095.5428.310.6682149.69100.0250.61 Z * ↓, Ao*↓, Asub*↓
↑, and ↓ shows the effect of the variation of the availability of the subcontractor and the availabi-lity of the producer.
Table 10. Policy Comparison.
Table 10. Policy Comparison.
DescriptionControl ParametersTotal CostObserved Difference
Δ-Cost (%)
Z *   k * A o *
Policy-A31.50.6707150.3234.33
Policy-B43.26-189.9350.59+32.14
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rivera-Gómez, H.; Medina-Marin, J.; Santana-Robles, F.; Montaño-Arango, O.; Barragán-Vite, I.; Cisneros-Flores, G. Impact of Unreliable Subcontracting on Production and Maintenance Planning Considering Quality Decline. Appl. Sci. 2022, 12, 3379. https://doi.org/10.3390/app12073379

AMA Style

Rivera-Gómez H, Medina-Marin J, Santana-Robles F, Montaño-Arango O, Barragán-Vite I, Cisneros-Flores G. Impact of Unreliable Subcontracting on Production and Maintenance Planning Considering Quality Decline. Applied Sciences. 2022; 12(7):3379. https://doi.org/10.3390/app12073379

Chicago/Turabian Style

Rivera-Gómez, Héctor, Joselito Medina-Marin, Francisca Santana-Robles, Oscar Montaño-Arango, Irving Barragán-Vite, and Gabriel Cisneros-Flores. 2022. "Impact of Unreliable Subcontracting on Production and Maintenance Planning Considering Quality Decline" Applied Sciences 12, no. 7: 3379. https://doi.org/10.3390/app12073379

APA Style

Rivera-Gómez, H., Medina-Marin, J., Santana-Robles, F., Montaño-Arango, O., Barragán-Vite, I., & Cisneros-Flores, G. (2022). Impact of Unreliable Subcontracting on Production and Maintenance Planning Considering Quality Decline. Applied Sciences, 12(7), 3379. https://doi.org/10.3390/app12073379

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

Article Metrics

Back to TopTop