Next Article in Journal
Linear Diophantine Fuzzy Fairly Averaging Operator for Suitable Biomedical Material Selection
Next Article in Special Issue
MRL-Based Model for Diverse Bidding Decision-Makings of Power Retail Company in the Wholesale Electricity Market of China
Previous Article in Journal
Orbit Tracing Properties on Hyperspaces and Fuzzy Dynamical Systems
Previous Article in Special Issue
Corn Disease Recognition Based on Attention Mechanism Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bayesian Statistical Method Enhance the Decision-Making for Imperfect Preventive Maintenance with a Hybrid Competing Failure Mode

1
School of Computer Science and Software, Zhaoqing University, Zhaoqing 526061, China
2
Department of International Business, Ming Chuan University, Taipei 33300, Taiwan
3
College of Management, Shu-Te University, Kaohsiung 82445, Taiwan
*
Authors to whom correspondence should be addressed.
Axioms 2022, 11(12), 734; https://doi.org/10.3390/axioms11120734
Submission received: 18 October 2022 / Revised: 6 December 2022 / Accepted: 9 December 2022 / Published: 15 December 2022
(This article belongs to the Special Issue A Hybrid Analysis of Information Technology and Decision Making)

Abstract

:
The study aims to provide a Bayesian statistical method with natural conjugate for facilities’ preventive maintenance scheduling related to the hybrid competing failure mode. An effective preventive maintenance strategy not only can improve a system’s health condition but also can increase a system’s efficiency, and therefore a firm needs to make an appropriate strategy for increasing the utilization of a system with reasonable costs. In the last decades, preventive maintenance issues of deteriorating systems have been studied in the related literature, and hundreds of maintenance/replacement models have been created. However, few studies focused on the issue of hybrid deteriorating systems which are composed of maintainable and non-maintainable failure modes. Moreover, due to the situations of the scarcity of historical failure data, the related analyses of preventive maintenance would be difficult to perform. Based on the above two reasons, this study proposed a Bayesian statistical method to deal with such preventive maintenance problems. Non-homogeneous Poisson processes (NHPP) with power law failure intensity functions are employed to describe the system’s deterioration behavior. Accordingly, the study can provide useful ways to help managers to make effective decisions for preventive maintenance. To apply the proposed models in actual cases, the study provides solution algorithms and a computerized architecture design for decision-makers to realize the computerization of decision-making.

1. Introduction

Preventive maintenance (PM) is in charge of maintaining equipment or facilities (repairable systems) in good condition. However, they might cause catastrophic damage and consequential losses when such equipment or facilities in production lines fail. As a result, it is critical to proceed with adequately preventative maintenance to avoid any damages and losses in order to keep these equipment or facilities in a healthy state. Preventive maintenance is capable of delaying system deterioration and returning systems to better condition, lowering failure rates and extending system lifetime. In other words, preventative maintenance has a significant influence on quality and cost, and therefore a firm should be concerned with developing appropriate maintenance programs in order to boost its competitiveness. There has been a lot of attention previous research works in the field of preventive maintenance modeling and optimization during the last several years.
Usually, preventive maintenance policies are mostly based on time intervals. There are two PM policies found in the literature: periodic and sequential (non-periodic). Periodic preventive maintenance policy is the manufacturer provides its maintenance work with equal time intervals. A sequential preventive maintenance policy is characterized by the search for the optimal number of maintenance actions at the optimal intervals. Therefore, the sequential preventive policy may provide an unequal sequence of intervals for minimizing the related costs. Park et al. [1] calculated the ideal duration and amount of PM activities based on a periodic PM program with minimum repair service after breakdowns. Yeh and Lo [2] demonstrated that the ideal PM interval between two consecutive PM activities is equal to the degree of PM and that providing an equivalent degree of PM is the best approach to decrease predicted warranty costs. Jung and Park [3] proposed an optimum post-warranty PM policy by reducing predicted long-run PM expenses. Seo and Bai [4] illustrated a periodic PM strategy for two scenarios where the operating time of PM might be disregarded or not. Yeh and Chang [5] determined the best failure rate and maintenance strategy for the lifetime of the equipment. Das and Sarmah [6] provided an overview of optimization models for preventative replacement with the related constraints in heavy-process industries. Yeh et al. [7] evaluated the impact of various PM cost functions on a leased product with a Weibull lifetime distribution’s periodic PM policy. Bouguerra et al. [8] proposed a mathematical model for various PM plans when customers choose to purchase an extended warranty. They discovered a viable compromise to establish a win-win situation between producers and customers in terms of warranty costs. Chang and Lin [9] developed an ideal PM policy for repairable items with extended warranties. They considered that manufacturers could give a slight discount to consumers to incent their intention to purchase extended warranty contracts. Under the context of reliability-based optimization, Beaurepaire et al. [10] proposed the best model of mechanical component maintenance scheduling. Their model is different from the traditional approaches based on linear fracture mechanics. Schutz and Rezg [11] provided a methodology for determining an optimal product maintenance program to guarantee that the minimal reliability meets the consumers’ requirements. Kim and Ozturkoglu [12] reduced the classic preventive maintenance issue to an integers programming problem. Khojandi et al. [13] investigated the optimal lifetime reward maintenance strategies for perfect and imperfect maintenance situations. They showed the tradeoff between the virtual age of the systems and the incentive rate for decision-makers. Yuan and Lu [14] suggested an effective way to solve a reliability-based optimization issue that combines the weighted approach with sequential approximation optimization. Lu et al. [15] proposed a joint model of sequential PM and quality improvement for deteriorating systems in manufacturing industries. The study is superior and more appropriate for maintaining machines in a production system. Wang and Djurdjanovic [16] also proposed a joint model for PM scheduling with consideration of stocks and logistics issues. They proposed an integrated policy to trigger PM for working parts. Zhou et al. [17] proposed a sequential PM model with a reducing failure factor. The model can be applied to urban bus systems’ maintenance works. García and Salgado [18] presented a case study to analyze the selection of PM strategies in multistage industrial facilities and equipment. Their study utilized some individual indicators to evaluate which PM strategies are better. Diatte et al. [19] proposed a methodology for improving brake systems in automobile industries, and the methodology can make the integration of the machine’s reliability, availability, and maintainability into system engineering and dependability analyses for reducing related costs and increasing system’s reliability.
Some related studies regarding the issues of competing failures or risks models have been proposed or developed in the past. Generally, maintainable and non-maintainable failure modes compete to cause the system to fail. Furthermore, the three maintenance operations should be taken into account: imperfect preventative maintenance, minor repairs, and replacements. Some competing risk models were also proposed to reflect the complex systems’ deterioration [20,21,22]. Salinas-Torres et al. [23] proposed a competing-risks model with a Bayesian statistical method to estimate a system component’s survival time. They used Dirichlet multivariate processes to deduce the parameters’ estimator. Yousef et al. [24] applied Bayesian and Non-Bayesian analyses to the reliability of the stress-strength system. The proposed Bayesian estimators can be achieved by the Markov chain Monte Carlo method. Wang and Miao [25] applied a semi-Markov model to optimize firms’ preventive maintenance policy. Their model would be useful to any two symmetric components to avoid system unbalance. Alotaibi et al. [26] also utilized a Bayesian analysis and Monte-Carlo simulations to estimate the parameters of a mixture bivariate exponential model. They applied the proposed method to motor data analysis. Yousef et al. [27] applied a Bayesian estimation method and Monte Carlo simulation to evaluate a system reliability. They also compared the performance between maximum likelihood estimators and Bayesian estimators. Liu et al. [28] applied a Bayesian estimation method to evaluate products’ reliability. This method can achieve both point and confidence interval estimation for the critical parameters. Zequeira and Berenguer [29] proposed a hybrid model in which the two failure modes (maintainable and non-maintainable) were believed to be dependent. It is intriguing to investigate the impact of the two failure models. El-Ferik and Ben-Daya [30] also proposed a hybrid model that includes adjustment variables in both the hazard rate and the effective age. In El-Ferik and Ben-Daya’s model, the effect of PM is assumed to be imperfect. Kahrobaee and Asgarpoor [31] proposed a hybrid analytical-simulation approach to solve deteriorating equipment’s PM works under some systems’ constraints. They apply their approach to the wind turbines industry, and expected rewards and penalties are also taken into consideration. Rafiee et al. [32] proposed a condition-based PM policy considering competing risks of internal deterioration and external shocks. In their model, the external shocks arrive at random times and can be divided into two categories based on their impacts on the system: (1) fatal shocks that can cause the system to fail immediately (2) non-fatal shocks that can only damage the system by randomly. Their model can apply in micro-electro-mechanical systems. Zhou et al. [33] proposed a hybrid PM model to reflect the reliability status of leased equipment. Furthermore, it also can clearly discriminate between the impacts from external shocks and from internal deterioration. Yang et al. [34] modeled complex industrial systems involving a hybrid failure mode, degradation-based failure and sudden failure. Their proposed condition-based PM strategy can be applied in oil pipeline industries. Cao [35] also proposed a condition-based PM strategy, but he took competing failure processes into consideration. He successfully utilized the condition-based PM strategy and a genetic algorithm to make periodic inspection and maintenance plans for wind-driven generators. Liu et al. [36] proposed hybrid maintenance models to deal with the issue of competing failures. In their model, all the maintenance actions are perfect, and all the actions are instant, it implies that the time cost can be negligible. However, their model didn’t take imperfect repair issues into consideration. Basílio et al. [37] proposed a review study regarding the applications of multi-criteria decision aid methods. They provided a complete overview of multi-criteria methods to comprehend the current and future development patterns of multi-criteria decision-making studies. The study can give useful research directions for applying multi-criteria techniques in industries.
Based on the above considerations, this work provided Bayesian decision models to cope with preventative maintenance with the hybrid deteriorating system. The deteriorating behavior of the repairable system can be described by a non-homogeneous Poisson process (NHPP) with a hybrid of power law failure intensity functions. Furthermore, the study investigates how different levels of preventive maintenance influence the related costs and the frequency of the system’s breakdowns. Accordingly, we consider the study has the following advantages: (1) Few studies focused on the issue of hybrid deteriorating systems, which are composed of maintainable and non-maintainable failure modes. Furthermore, due to some situations of the scarcity of historical failure data, the related analyses of preventive maintenance will be difficult to perform. For the above two reasons, this study proposed Bayesian decision models to deal with such a preventive maintenance problem with the hybrid deteriorating system. (2) The study’s Bayesian decision process can be divided into two phases. In the first phase, the decision maker can use the domain experts’ knowledge and judgment to evaluate the deterioration using four prior distribution statistical characteristics. In the second phase, the decision maker can collect failure data to re-estimate the deterioration by using the NHPP likelihood functions and the properties of the natural conjugate. It can effectively extend the related applications of the study. (3) In order to apply the proposed models in actual cases, the study provides solution algorithms and a computerized architecture design for decision-makers to realize the computerization of decision-making. Accordingly, the study can easily be applied to various manufacturing industries that want to effectively manage their durable and high-priced equipment or facilities in the factory.
The rest of this paper is organized as follows: Section 2 presents the proposed hybrid competing mode, the estimation of a deteriorating system’s failures under periodic preventive maintenance, and the estimation of the related costs. Section 3 provides a Bayesian decision process by using domain experts’ knowledge with collected information. The design of a computerized information system is also presented in the section. Section 4 present the numerical application and sensitivity analysis. Finally, Section 5 presents the concluding remarks and the future study.

2. Preventive Maintenance for Deteriorating Systems with a Hybrid Deterioration

Both non-maintainable and maintainable failure modes may exist in a multi-component system, and therefore they will compete to cause failures. In the maintainable failure mode, a component can be replaced with a new one to restore the system to its original status. However, in the non-maintainable failure mode, a component cannot be replaced to alleviate the system’s aging. From the whole system’s perspective, the effect of preventive maintenance can only present imperfect recovery. Therefore, a hybrid model which is related to competing failure processes would be more adequate for measuring such deterioration model in practice. The following subsections will introduce the hybrid deterioration model with mathematical analysis.

2.1. Hybrid Competing Failure Mode

In considering system reliability and operational stability, scheduled PM works can be performed for reducing equipment or facilities’ breakdowns and also to avoid potential disasters. Either periodic PM or non-periodic PM policy can be adopted for enhancing the equipment or facilities’ reliability. However, a scheduled periodic PM policy is taken into consideration in this study since it may be more practically manageable to the manager. Suppose that a deteriorating system (equipment) will undergo N 1 PM works during its lifetime T L , where the intervals of periodic PM are equal and designated to be x by the reliability engineering department, and the whole system is replaced at the timing of the Nth PM work. Any breakdown which occurs within an interval two PM actions would cause minimal repairs, and minimal repairs cannot reduce the system aging. Accordingly, a NHPP with a power law intensity function λ ( t | α , β ) = α β t β 1 is for describing the process of system’s deterioration, where α and β denote the scale factor and the shape factor, respectively. The domains of α and β are within the range [0, ) according to the property of the power law intensity function. Therefore, the power law form would be more flexibility and manageability than other intensity functions.
Since PM works can partially improve the system reliability and can also extend the original system lifetime. However, it is still unable to stop the system’s aging process. Ultimately, the whole system will need to be replaced after the Nth PM action (the system lifetime would be T L = N x ) in the consideration of cost- effectiveness. In this study, C ( N ) denotes the expected cost per unit time with respect to the PM number N . The expected cost includes the minimal repair cost C m r , the penalty cost C p l , the ith PM cost C p m i , and the replacement cost C r p . In general, the cost of performing a PM action should be relatively higher than the cost of performing repairs in the initial stage since the machine’s abrasion status might not be serious when the machine is in its early age. However, when the machine’ abrasion status becomes more serious over time, the repair cost will more than the PM cost if the age reduction factor of PM δ is relatively low and cannot effectively reduce the exponential increase of a system’s breakdowns. Figure 1 illustrates of the preventive maintenance model of an equipment or facility.
Furthermore, a facility or equipment may include non-maintainable and maintainable components in a system, and therefore the two modes (non-maintainable mode and maintainable mode) are needed to be integrated into one for presenting the phenomenon of a system’s imperfect recovery. Therefore, the two intensity functions of the system deterioration will be devised as λ o ( t | α o , β o ) = α o β o t β o 1 (non-maintainable mode) and λ p ( t | α p , β p ) =   α p β p t β p 1 (maintainable mode) respectively. Besides, please note that the values of scale and the shape parameters may be unknown and they will be discussed as uncertain states in the Bayesian analysis.
The assumptions of the proposed model are stated as follows:
(1)
The system’s deterioration behaves as a non-homogeneous Poisson process (NHPP).
(2)
The system’s deterioration is composed of maintainable and non-maintainable failure modes.
(3)
A PM cannot restore the whole system to a brand-new state; instead, it can restore the whole system to some state as better-than-now.
(4)
Any breakdowns occurring within the interval between two PM actions cause a minimal repair.
The following notations are used in the analysis throughout this study (Table 1):

2.2. Estimation of a System’s Failures under Preventive Maintenances

Since the interval time of PM is x which was set by maintenance engineers, and t 1 + denotes the effective age of a system after the time point of the first PM. It can be deduced by t 1 + = x δ x = ( 1 δ ) x because t 1 + is influenced by the age reduction factor δ under the maintainable mode. Therefore, the effective age of a system before and after the time point of the kth PM can be presented as
t k = k x ( k 1 ) δ x = ( k 1 ) ( 1 δ ) x + x = ( ( k 1 ) ( 1 δ ) + 1 ) x ,
and
t k + = t k δ x = k x ( k 1 ) δ x δ x = k ( 1 δ ) x
respectively.
Due to the system’s deterioration which includes non-maintainable and maintainable modes, the hybrid intensity function can be written as follows:
λ h ( T L | α o , β o , α p , β p ) = λ o ( k x ) + λ p ( ( ( k 1 ) ( 1 δ ) + 1 ) x ) = α o β o ( k x ) β o 1 + α p β p ( ( ( k 1 ) ( 1 δ ) + 1 ) x ) β p 1 ,
Moreover, it means that the system belongs to perfect maintenance if the age reduction factor δ is equal to one. Thus, the expected number of failures of the equipment or facility is given by
N m r ( N , x , δ | α o , β o , α p , β p ) = k = 1 N ( k 1 ) x k x λ h ( t , δ | α o , β o , α p , β p ) d t = k = 1 N ( k 1 ) x k x λ o ( t | α o , β o ) d t + k = 1 N t k 1 + t k λ p ( t | α p , β p ) d t = k = 1 N ( k 1 ) x k x λ o ( t ) d t + k = 1 N ( k 1 ) ( 1 δ ) x ( ( k 1 ) ( 1 δ ) + 1 ) x λ p ( t ) d t = α o ( N x ) β o + k = 1 N α p ( ( ( ( k 1 ) ( 1 δ ) + 1 ) x ) β p ( ( k 1 ) ( 1 δ ) x ) β p ) .
However, if the age reduction factor δ is equal to one, the expected number of the failures of the equipment or facility can be rewritten as follows:
N m r ( N , x | α o , β o , α p , β p ) = α o ( N x ) β o + N α p x β p       ( δ = 1 ) .
Since the breakdown process of a deterioration system can be modeled as a Non-Homogenous Poisson Process, the probability of the number of breakdowns N 0 in the interval ( N x , ( N + 1 ) x ) is thus given by
P r { N m r ( N + 1 , x , δ | α o , β o , α p , β p ) N m r ( N , x , δ | α o , β o , α p , β p ) = N 0 } = ( N m r ( N + 1 , x , δ | α o , β o , α p , β p ) N m r ( N , x , δ | α o , β o , α p , β p ) ) N 0 × e ( N m r ( N + 1 , x , δ | α o , β o , α p , β p ) N m r ( N , x , δ | α o , β o , α p , β p ) ) N 0 ! .
The reliability of a product R ( T L = N x ) will decline with time, and therefore it can denote as
R ( T L = N x ) = P r { N m r ( N , x , δ | α o , β o , α p , β p ) = 0 } = e N m r ( N , x , δ | α o , β o , α p , β p ) .
Figure 2 illustrates preventive maintenance between perfect recovery and imperfect recovery under a hybrid deterioration. The difference between perfect and imperfect recoveries is that the critical component of a perfect recovery system can be maintainable or can be replaced with a new one to restore the system to its original status. As can be seen in the middle-left side of Figure 2, the maintainable component can be fully restored to its original status after preventive maintenance.

2.3. Evaluation of Repair and Maintenance Costs of a Facility

Some expenditures will be incurred as long as equipment or facility operates during the system lifetime. The incurred expenditures are mainly from repair, penalty, replacement, and PM costs. The repair cost ( C m r ) means the expected cost to perform a minimal repair. The penalty cost ( C p l ) means that the cost was incurred from the loss of a production line shutdown if the actual repair time exceed the time limit ( φ ). Since any breakdown of the facility will be rectified by minimal repairs, the time of a minimal repair ( t r ) will need to be measured. Therefore, the repair time can be regarded as a random variable, and it is assumed to follow a Gamma probability distribution. The expected repair time over the tolerable waiting time limit φ can be expressed as
E [ t r | ω , η ] = φ t r ( η ω t r ω 1 Γ ( ω ) e η t r ) d t r = Γ ( 1 + ω ) ω Γ ( ω ) + Γ ( 1 + ω , η φ ) η Γ ( ω )
The parameters ω and η can be estimated by ω = ( E ( t r ) σ ( t r ) ) 2 and η = E ( t r ) σ ( t r ) 2 under engineers’ judgment or historical data of repairs. If the repair time is over the time limit φ , the penalty cost will be incurred by the owner of the equipment or facility. Γ ( z 0 ) denotes a Gamma function with the parameter z 0 and Γ ( z 1 , z 2 ) denotes an upper incomplete Gamma function with the parameters z 1 and z 2 .
The estimation of the equipment or facility’s deterioration is critical to the owner, and therefore the manager needs to accurately evaluate the expected failure number of the equipment or facility during its lifetime. Supposed that the failure process follows an NHPP with a power-law intensity function λ ( t ) . Therefore, the expected number of failures during the lifetime [ 0 ,   T L ] under the age reduction factor in effective age δ p m q and the interval of PM x is N m r ( N , x , δ p m q | α o α o , β o , α p , β p ) . Accordingly, the total expected repair cost during the lifetime can be given as
( C m r + C p l E [ t r | ω , η ] ) N m r ( N , x , δ p m q | α o α o , β o , α p , β p ) .
Moreover, the PM cost will get higher and higher for sequential PM activities during the lifetime due to mechanical aging of a deteriorating system. Therefore, the PM cost should be related to the ith number of PM actions with the age reduction factor. Based on this, the PM cost is defined as
C p m k = C F ( 1 + τ ( k 1 ) x )
where τ denotes the periodically increasing rate of the PM cost, and C F is the base cost of performing a PM work. Furthermore, different PM alternatives bring a different degree of the system’s recovery but it also influences different PM costs to the firm. Suppose that a series of PM alternatives M P = { M P 1 ,   M P 2 ,   M P q , , M P Q } can be selected, and the corresponding PM cost and the expected failure numbers can be redefined as follows:
C p m q ( C F q , τ q , x , T L ) = k = 1 T L / x 1 C p m k q = k = 1 T L / x 1 C F q ( 1 + τ q ( k 1 ) x ) .
It is important to realize the process of deterioration of the product in order to evaluate the costs of repair during the system’s lifetime. Based on an assumption that the failure times can be drawn from an NHPP with a specific intensity function λ h ( · ) , an estimate of the number of expected failures N m r ( N , x , δ p m q | · ) under the time interval of PM x as well as the age reduction factor δ p m q in the effective age of the system, the total repair cost during the system’s lifetime [ 0 ,   T L ] can be calculated as follows:
C m r N m r ( N , x , δ p m q | α o , β o , α p , β p ) = C m r ( α o ( N x ) β 1 + k = 0 N 1 α p ( ( 1 + k ( 1 δ p m q ) ) x ) β 2 ( k ( 1 δ p m q ) x ) β 2 ) )
If the age reduction factor δ p m q is equal to one, the total repair cost can be rewritten as follows:
C m r N m r ( N , x , α o , β o , α p , β p ) = C m r ( α o ( N x ) β o + N α p x β p ) .

2.4. Optimal Preventive Maintenance Schedule with Consideration of Multiple PM Alternatives

Generally, managers need to consider how to determine the optimal maintenance schedule to minimize the total expected costs associated with the project requirements. Therefore, in consideration of all the candidate PM alternatives, the total expected cost per unit time under the system lifetime T L can be given as follows:
C ( N | M P q ) = k = 1 N 1 C p m k + ( C m r + C p l E [ t r | ω , η ] ) N m r ( N , x , δ p m q | α o , β o , α p , β p ) + C r p T L = C p m q ( C F q , τ p , x , T L ) + ( C m r + C p l E [ t r | ω , η ] ) N m r ( N , x , δ p m q | α o , β o , α p , β p ) + C r p N x .
The convexity of the cost function with respect to N under a specific PM alternative M P q can be justified if the two inequalities C ( N + 1 | M P q ) C ( N | M P q ) and C ( N | M P q ) < C ( N 1 | M P q ) are both held, and the optimal N * can therefore be obtained. Proposition 1 give the proof of the convexity of C ( N | M P q ) .
Proposition 1. 
Given the intensity function N m r ( N , x , α 1 , β 1 , α 2 , β 2 , δ ) is strictly increasing. As long as the two inequalities C ( N + 1 | M P q ) C ( N | M P q ) and C ( N | M P q ) < C ( N 1 | M P q ) can be both held under a specific number N, the convexity of the cost function C ( N | M P q ) with respect to N can be assured.
Proof. 
For C ( N + 1 | M P q ) C ( N | M P q ) , we have C ( N + 1 | M P q ) C ( N | M P q ) 0
k = 1 N C p m k + ( C m r + C p l E [ t r | ω , η ] ) N m r ( N + 1 , x , δ p m q | α o , β o , α p , β p ) + C r p ( N + 1 ) x k = 1 N 1 C p m k + ( C m r + C p l E [ t r | ω , η ] ) N m r ( N , x , δ p m q | α o , β o , α p , β p ) + C r p N x 0 k = 1 N C p m k + ( C m r + C p l E [ t r | ω , η ] ) N m r ( N + 1 , x , δ p m q | α o , β o , α p , β p ) + C r p ( N + 1 ) x ( N + 1 N ) ( k = 1 N 1 C p m k + ( C m r + C p l E [ t r | ω , η ] ) ( N , x , δ p m q | α o , β o , α p , β p ) + C r p ) ( N + 1 ) x 0 ( C m r + C p l E [ t r | ω , η ] ) × { N m r ( N + 1 , x , δ p m q | α o , β o , α p , β p ) ( 1 + 1 N ) N m r ( N , x , δ p m q | α o , β o , α p , β p ) } ( N + 1 ) x ( 1 N ) C r p C p m N + ( 1 N ) k = 1 N 1 C p m k ( N + 1 ) x ( C m r + C p l E [ t r | ω , η ] ) { N m r ( N + 1 , x , δ p m q , α o , β o , α p , β p ) ( 1 + 1 N ) N m r ( N , x , δ p m q | α o , β o , α p , β p ) } ( 1 N ) C r p C p m N + ( 1 N ) k = 1 N 1 C p m k ( C m r + C p l E [ t r | ω , η ] ) { ( N ) N m r ( N + 1 , x , δ p m q | α o , β o , α p , β p ) ( N + 1 ) N m r ( N , x , δ p m q , α o , β o , α p , β p ) } k = 1 N 1 C p m k N C p m N + C r p
( N ) N m r ( N + 1 , x , δ p m q | α o , β o , α p , β p ) ( N + 1 ) N m r ( N , x , δ p m q | α o , β o , α p , β p ) k = 1 N 1 C p m k N C p m N + C r p C m r + C p l E [ t r | ω , η ] .
For C ( N | M P q ) < C ( N 1 | M P q ) , we have C ( N | M P q ) C ( N 1 | M P q ) < 0
k = 1 N 1 C p m k + ( C m r + C p l E [ t r | ω , η ] ) N m r ( N , x , δ p m q | α o , β o , α p , β p ) + C r p N x k = 1 N 2 C p m k + ( C m r + C p l E [ t r | ω , η ] ) N m r ( N 1 , x , δ p m q | α o , β o , α p , β p ) + C r p ( N 1 ) x < 0 ( 1 1 N ) k = 1 N 1 C p m k + ( 1 1 N ) ( C m r + C p l E [ t r | ω , η ] ) N m r ( N , x , δ p m q | α o , β o , α p , β p ) + ( 1 1 N ) C r p ( k = 1 N 2 C p m k + ( C m r + C p l E [ t r | ω , η ] ) N m r ( N 1 , x , δ p m q | α o , β o , α p , β p ) + C r p ) < 0 ( C m r + C p l E [ t r | ω , η ] ) { ( N 1 N ) N m r ( N , x , δ p m q | α o , β o , α p , β p ) N m r ( N + 1 , x , δ p m q | α o , β o , α p , β p ) } + ( N 1 N ) C p m N 1 1 N k = 1 N 2 C p m k ( 1 N ) C r p < 0 ( C m r + C p l E [ t r | ω , η ] ) { ( N 1 ) N m r ( N , x , δ p m q | α o , β o , α p , β p ) ( N ) N m r ( N 1 , x , δ p m q | α o , β o , α p , β p ) } + ( N 1 ) C p m N 1 k = 1 N 2 C p m k C r p < 0 ( C m r + C p l E [ t r | ω , η ] ) { ( N 1 ) N m r ( N , x , δ p m q | α o , β o , α p , β p ) ( N ) N m r ( N 1 , x , δ p m q , α o , β o , α p , β p ) } < k = 1 N 2 C p m k ( N 1 ) C p m N 1 + C r p
( N 1 ) N m r ( N , x , δ p m q | α o , β o , α p , β p ) ( N ) N m r ( N 1 , x , δ p m q | α o , β o , α p , β p ) < k = 1 N 2 C p m k ( N 1 ) C p m N 1 + C r p C m r + C p l E [ t r | ω , η ] .
Let L ( N ) = ( N ) N m r ( N + 1 , x , δ p m q | α o , β o , α p , β p ) ( N + 1 ) N m r ( N , x , δ p m q | α o , β o , α p , β p ) for N =1, 2,…, and L ( N ) = 0 for N = 0 . The condition that L ( N ) is strictly increasing with N is supported and the two inequalities (15) and (16) can hold simultaneously, and there would exist the PM number N can minimize the total expected cost per unit time under the system lifetime. Due to the fact that the failure intensity function is increasing with time for deteriorating systems, i.e., λ h ( ( N + 1 ) x ) > λ h ( N x ) > > λ h ( 0 ) for x , 2 x , , N x , ( N + 1 ) x , , we then have as follows:
L ( N ) L ( N 1 ) = { ( N ) N m r ( N + 1 , x , δ p m q | α o , β o , α p , β p ) ( N + 1 ) N m r ( N , x , δ p m q | α o , β o , α p , β p ) } { ( N 1 ) N m r ( N , x , δ p m q | α o , β o , α p , β p ) ( N ) N m r ( N 1 , x , δ p m q | α o , β o , α p , β p ) } > 0 = ( N ) { N m r ( N + 1 , x , δ p m q | α o , β o , α p , β p ) + N m r ( N 1 , x , δ p m q | α o , β o , α p , β p ) } ( 2 N ) N m r ( N , x , δ p m q | α o , β o , α p , β p ) > 0 = { N m r ( N + 1 , x , δ p m q | α o , β o , α p , β p ) + N m r ( N 1 , x , δ p m q | α o , β o , α p , β p ) } 2 N m r ( N , x , δ p m q | α o , β o , α p , β p ) > 0
Therefore, the convexity of the cost function C ( N | M P q ) with respect to N is thus assured according to Jensen’s inequality. □
The heuristic solution algorithm for obtaining the minimal cost by setting N * and M P q * can be described in Figure 3.

3. Bayesian Decision Process by Using Domain Experts’ Judgment and Collected Information

3.1. Analysis by the Natural Conjugate Probability Distribution

It might be not easy to perform Bayesian decision analysis due to the fact that numerical integration is needed to derive the prior and posterior distributions. As the state space in our case contains multiple random variables, the previous analysis would have been much more complicated. Huang and Bier [38] proposed a natural conjugate prior distribution for the power law deteriorating model for repairable systems, and the form is as follows:
f ( α , β ) = K α κ 1 β κ 1 ( e ϖ υ κ ) β 1 e α ψ υ β .
In order to make sure the distribution sums up to one, K is used to be a normalizing factor. The main advantage is to using the natural conjugate prior distribution to proceed with a straightforward analysis instead of the complicated traditional calculations. In Equation (17), the joint probability distribution with the desired prior marginal moments are composed of the four parameters ( ψ , ϖ , κ , and υ ). By applying Equation (17) to the prior probability distribution of α and β , it will be easy to deduce the expected failure number based on the prior distribution, and its form is given by
E P r i [ N m r ( T L | α , β ) ] = 0 0 0 T L α β t β 1 f ( α , β ) d t d β d α κ ϖ κ ψ ( ϖ + l n [ υ ] l n [ T L ] ) κ .
Here, T L denotes the actual age of the equipment or facility. The prior analysis of the expected failure number can be performed by straightforwardly calculating E P r i [ N m r ( T L | α , β ) ] in Equation (18) with the four parameters ( ψ , ϖ , κ , and υ ) which are specified by the reliability or domain experts with their prior knowledge and judgment about the deteriorating system ( μ α , σ α , μ β , and σ β ). μ α , σ α , μ β , and σ β denote the mean values and standard deviations of α and β , respectively. Furthermore, in order to calculate E P r i [ N m r ( T L | α , β ) ] , the values of the four parameters need to be obtained first. Equations (19)–(22) can be used to obtain the four parameters values as follows:
ϖ = μ β σ β 2
κ = ( μ β σ β ) 2
υ = e ϖ ( ξ 1 / κ 1 + ξ 1 / κ ( ξ 1 / κ 1 ) )   ( where   ξ = μ α 2 σ α 2 + 1 μ β 2 σ β 2 + 1 ) ,
and
ψ = ( κ μ α ) ( ϖ ϖ + l n [ υ ] ) κ
However, if the decision maker may not be convinced by the result of the prior analysis, he/she may collect failure datasets in practice to adjust the prior analysis. It is called the posterior analysis in a Bayesian decision process. Once the decision maker wants to proceed with the posterior analysis, he/she needs to prepare the extra budget for accelerated deterioration experiments to increase the accuracy of the prediction. If a further investigation is undertaken, the sample size should be carefully examined in consideration of the budget of the experiments. If the n breakdown times are collected from accelerated deterioration experiments ( x 1 , x 2 , , x n ) , the property of natural conjugate families can be used to obtain the posterior distribution of α and β without further computation. Its form is as follows:
g ( α , β , D ( n ) ) L ( D ( n ) | α , β ) f ( α , β ) = K α κ + n 1 β κ + n 1 ( e ϖ υ κ i = 1 n x i ) β 1 e α ( ψ υ β + x n β ) .
L ( D ( n ) | α , β ) = α n β n ( i = 1 n x i ) β 1 e α x n β is the likelihood function of NHPP with the power-law intensity function, and K denotes a normalizing factor to ensure the distribution sums up to one. However, due to the math complexity of the expected failure number of the posterior analysis, the closed form expression cannot be obtained. Fortunately, numerical integration methods (e.g., Monte Carlo numerical integration) can be applied to calculate the prediction of the expected failure number based on the posterior distribution, and it can be given by
E P o s [ N m r ( T L | α , β , D ( n ) ) ] = 0 0 0 T L α β t β 1 g ( α , β , D ( n ) ) d t d β d α 0 0 α T L β g ( α , β , D ( n ) ) d β d α .
Therefore, both experts’ knowledge and sampling information can be considered in the posterior analysis. Moreover, in consideration of the hybrid deterioration and the number of PM with the age reduction factor, the state space can be presented as Θ : { α o , β o , α p , β p | α o ( 0 , ) ,   β o ( 0 , ) ,   α p ( 0 , ) , β p ( 0 , ) } , and the expected failure numbers of the prior analysis and the posterior analysis can be given as follows:
E P r i [ N m r ( N , x , δ p m q | α o , β o , α p , β p ) ] = κ o ϖ o κ o ψ o ( ϖ o + l n [ υ o ] l n [ N x ] ) κ o + ( κ p ϖ p κ p ψ p ) k = 0 N 1 { ( ϖ p + l n [ υ p ] l n [ ( 1 + k ( 1 δ p m q ) ) x ] ) κ p ( ϖ p + l n [ υ p ] l n [ k ( 1 δ p m q ) x ] ) κ p } ,
and
E P o s [ N m r ( N , x , δ p m q | α o , β o , α p , β p ) ] = 0 0 α o T L β o g ( α o , β o ) d β o d α o + k = 0 N 1 { 0 0 α p ( ( 1 + k ( 1 δ p m q ) ) x ) β p g ( α p , β p ) d β p d α p 0 0 α p ( k ( 1 δ p m q ) x ) β p g ( α p , β p ) d β p d α p } .
It should be noted that the parameters of the hybrid deteriorating function ( ϖ o , κ o , υ o , ψ o , ϖ p , κ p , υ p , and ψ p ) can be calculated by applying Equations (19)–(22). After getting the parameters’ value, the expected failure numbers, E P r i [ N m r ( · ) ] and E P o s [ N m r ( · ) ] ), can be obtained to calculate the expected costs per unit time E P r i [ C ( N | M P q ) ] and E P o s [ C ( N | M P q ) ] from Equations (27) and (28), and their form are as follows:
E P r i [ C ( N | M P q ) ] = C p m q ( C F q , τ q , x , T L ) + ( C m r + C p l E [ t r | ω , η ] ) E P r i [ N m r ( N , x , δ p m q | α o , β o , α p , β p ) ] + C r p N x
and
E P o s [ C ( N | M P q ) ] = C p m q ( C F q , τ q , x , T L ) + ( C m r + C p l E [ t r | ω , η ] ) E P o s [ N m r ( N , x , δ p m q | α o , β o , α p , β p ) ] + C r p N x .
By applying the solution algorithm shown in Figure 3, the minimal expected cost of the prior analysis or the posterior analysis can be determined with the settings of the Bayesian parameters.

3.2. The Bayesian Decision Process

The decision maker should confirm whether the relevant assumptions are satisfied or not in the case before deciding to proceed with the Bayesian analysis. The reliability or domain experts can provide the judgment of the parameters’ value ( μ α o , σ α o , μ β o , σ β o , μ α p , σ α p , μ β p , and σ β p ) according to their prior knowledge and experience. After this step, the decision makers can proceed with the prior analysis to minimize the expected cost. Suppose the decision makers consider that the prior analysis is reliable. In that case, it is not necessary to collect extra experimental data to adjust their prior judgment and just proceed with the choice of the candidate PM alternatives.
However, if the decision makers feel that the prior analysis might not be reliable or convinced, they would request more information to justify or amend the prior analysis. However, how many experimental data will be enough? It depends on the firm’s budget and possible benefit for amending the original PM alternative. Once the decision makers decide to proceed with the posterior analysis, the cost of collecting this additional information should also be considered in the decision process. The posterior analysis can utilize additional information from accelerated deterioration experiments along with domain experts’ judgment to make the final decision. Nevertheless, the calculation of the posterior analysis is not easy to performance since numerical methods, Monte Carlo integration and computation engines would be needed to get the non-closed form solution. Figure 4 illustrates the analysis of the Bayesian decision process.

3.3. Computerized Information System Design

In general, the optimal decision would not be easily obtained without a computerized information system when dealing with such complicated mathematical models. The entire system can be split into two subsystems to improve manageability. Engineers and domain experts can use the model management system to update the model base and database. A decision support system is also implemented to give decision makers the knowledge they need to make informed choices. The engineers should inspect the expenses, failure intensity functions, likelihood of repair time, penalty cost, replacement cost, and related probabilities, etc. characteristics before operating the model management system. When the engineers have these data, they can use the model management system to store them in the database. Additionally, the lack of the deterioration information would make it difficult to assess how the new system is deteriorating and what the frequency of failures is within a time interval. Therefore, based on their expertise, the statistical characteristics and associated factors can be assessed by the reliability engineering domain experts. In addition, the engineers must gather failure data for the facilities from various engineering experiments if a posterior analysis is desired by the decision maker. Moreover, we can store or access data more effectively through the use of a data formalizing mechanism. This mechanism converts inconsistent data to more consistent data for storage and access to the database and model base. Additionally, computing engines could be required to handle the complexity of finding the best answer. By utilizing an application programming interface (API), system developers can move forward with all of the mathematical analyses for the decision support system. The design of the computerized information system is shown in Figure 5.

4. Application and Sensitivity Analyses

4.1. Application of Prior and Posterior Analyses

Suppose that a firm plans to purchase a batch of new industrial equipment for manufacturing its products. After the purchase and setting of the industrial equipment up in its production lines, the firm has to make an effective preventive maintenance plan for factor management. However, due to the lack of the equipment’s deterioration information, the decision makers of the firm would be hard to make an effective preventive maintenance plan to reduce unexpected equipment breakdowns for saving the related cost. Moreover, the industrial equipment’s deterioration can be categorized into maintainable and non-maintainable failure modes. It means that some components can be repaired or replaced to restore the equipment to a younger status but some components cannot be restored to the system’s status by repair or replacement works. Therefore, the firm needs to estimate the two modes’ deterioration to handle the equipment’s failures occurring. Since the new equipment did not take any accelerating deterioration experiments for some reasons, the equipment’s deterioration cannot be estimated by traditional statistical methods. Due to the fact that to proceed with a complete accelerating deterioration experiment will need a lot of time and expenditure, the firm may not proceed with such experiment in time. Accordingly, in order to solve the issue of insufficient data, the decision makers would try to apply a Bayesian decision process to estimate the deterioration. The Bayesian decision process can be separated into two phases: the prior analysis, which is assessed by reliability domain experts, and the posterior analysis, which requires failure data from engineering tests. In the first phase, the decision makers will ask domain experts to evaluate the maintainable and non-maintainable parameters’ statistical characteristics for the prior analysis. After the experts’ evaluation, the eight prior parameters are set as μ α o = 0.6, μ β o = 0.15, σ α o = 1.25, σ β o = 0.3125; μ α p = 0.85, μ β p = 0.2125, σ α p = 1.75, and σ β p = 0.4375. Besides, the firm’s engineering department proposed five candidate PM alternatives for decision makers. This information can be referred to Table 2. Besides, since different PM alternatives can bring different degrees of system’s recovery, the corresponding PM costs are also different. PM alternatives 4 & 5’s age reduction ( δ p m 4 = 0.95, δ p m 5 = 1.0) are higher than the others ( δ p m 1 = 0.80, δ p m 2 = 0.85, δ p m 3 = 0.90). The higher age reduction PM alternatives can effectively reduce the increase of the possible repair cost and penalty cost but they also increases the related PM cost. Therefore, it is hard to judge which PM alternative is best for the firm before the evaluation of the all PM alternatives. Moreover, the experts also evaluate the repair time’s statistical characteristics and consider the time is a random variable and follows a Gamma probability distribution. The statistical characteristics of the repair time are estimated as E ( t r ) = 5 h and σ ( t r ) = 3 h, respectively. However, the tolerable waiting time is only 4.5 h in practice. If the repair time is over the tolerable waiting time, the production line will be halted and bring the related loss. The loss can be evaluated as the penalty cost. The parameters ω and η of the Gamma distribution can be calculated by the simultaneous equations ω = E ( t r ) 2 / σ ( t r ) 2 and η = E ( t r ) / σ ( t r ) 2 . Based on the above mentioned, the detailed information of the five candidate PM alternatives with the firm’s domain experts evaluation and the cost parameters is given as Table 2.
After performing the prior analysis of the Bayesian decision process based on the proposed solution algorithm, the trend of the average cost of PM alternatives 1–5 are presented in Figure 6. Table 3 provides the related results of PM alternatives 2–4 in detail. The best PM alternative is no.3, and the annual average cost, PM cost, and replacement cost are estimated to be $6212, $2956, and $2105 at 9.5 years. According to the prior analysis results in Figure 6, the optimal lifetime of equipment’s replacement of PM alternatives 1–5 should be set to 9, 9, 9.5, 8.5, and 9 years respectively, and the annual average cost will be $6329, $6268, $6212, $6634, and $6585. It can be observed that the highly intensive PM alternative may not have a beneficial influence on decreasing the cost. In general, highly intensive PM alternatives (4 & 5) can lower equipment failure times to save expenditure on repairs but they need more the related PM costs. On the contrary, the lower intensive PM alternatives (1 & 2) are able to decrease the related PM costs but increase expenditure on repairs. However, although lower intensive PM alternatives may bring serious equipment’ failures at the post-phase, the firm may adopt the strategy of shortening equipment’s lifetime to prevent this disadvantage. According to the above mentioned, we understand that the firm cannot judge and decide the best PM alternative by qualitative analysis or without cost calculation. Moreover, it can be seen that the average cost of PM alternative 3 is lower than the other PM alternatives and the optimal equipment’s lifetime is 9.5 years (marked with * in Table 3 and Table 4). However, the average cost of the lower intensive PM alternatives are very close to PM alternative 3 before 8 years. It indicates that the firm may choose a medium to lower intensive PM alternative with a shorter lifetime if the firm encounter a budget or financial issue. Accordingly, the firm may adjust the present PM alternative to adapt different scenarios according to its financial health in practice.
However, if the firm didn’t have enough confidence to believe the outcomes of the prior analysis results, the decision makers will request more evidences to verify the system’s deterioration. Therefore, in order to adjust the previous analysis, it is necessary to gather the extra failure data by proceeding with accelerated deterioration experiments. In consideration of the limits budget and time, the firm can only gather few experimental data to revise the previous analysis result. After getting the extra experimental data to proceed with the posterior analysis, the data can integrate with the previous experts’ opinions of the prior analysis, and the complete analysis results are reported in Table 4. Figure 7 shows that the overall cost per unit and year of the posterior analysis is always smaller than that of the prior analysis. The best lifetime of an equipment will be extended to 10 years (marked with ** in Table 4) and the overall cost per unit and year will be $6056. After reviewing the posterior analysis’s results, it is clear that the prior analysis’s conclusions may be overly pessimistic since the re-estimating deterioration in the posterior analysis is less severe than that in the prior study. In other words, the firm should moderately extend the planned lifetime of equipment to save the related costs.

4.2. Sensitivity Analyses

Misjudging the parameters values of the maintainable and non-maintainable modes μ α p , μ β p , μ α o , and μ β o may have an influence on the predictions of the overall cost per unit and year, thus the firm should be aware of potential changes in the projections. As a result, sensitivity analysis may be undertaken to evaluate differences in the equipment’s lifetime and the overall cost. It is logical to assume that if the firm underestimates these parameters’ value, the related cost are also underestimated, and it will lead to poor judgments such as incorrectly extending equipment’s lifetime. According to Figure 8, the firm will prolong the lifetime of equipment by taking advantage of the lower growing degradation if these parameters values are decreasing. Moreover, since the variation effect of the non-maintainable parameters is greater than that of the maintainable parameters, the firm needs carefully to evaluate the non-maintainable parameters to avoid inappropriate decisions especially for shape factors. It can be seen in that the variation of μ β o between the range (−30%, +30%) causes the change of decision about lifetime from 7 years to 10 years. The variation of the overall cost is also huge. However, the effect of misjudging the shape factor of the maintainable parameter μ β p will be smaller because the system’s deterioration can be restored to a younger status after a PM action. It means that the growth of the effect age can be slowed down with time by periodic PM actions. Besides, misjudging the parameters of standard deviation σ α p , σ β p , σ α o , and σ β o will influence the range of confidence intervals of the expected overall cost but it may not result in risky decisions because the optimal decision of the equipment’s lifetime will be less changed with these parameters of standard deviation.
Besides, the related parameters regarding PM, repair, penalty would also influence the estimation of the overall cost. The base cost C F q seems the most important influencing factor in this case because the portion of this factor is almost 30–50% of the overall cost. Increasing the base cost will cause the optimal decision of shortening equipment’s lifetime. It can be seen that the variation of C F q between the range (−30%, +30%) leads to the change of decision about lifetime from 10.5 years to 8.5 years, and the overall cost at the individual optimal lifetime will go up from $5306 to $7061. Moreover, the increaing rate τ q is an important parameter for measuring the expected overall cost. Inflation, wage increases, or the cost of supplying components may be contributing factors to the increasing rate τ q . As a result, if the company believes or predicts that the vital components for the replacement of PM actives will rise in the future, the company should plan to reduce the equipment’s lifetime since the average replacement cost of the total equipment cannot cover the cumulative PM cost. It can be seen in Figure 9. The variation of τ q between the range (−30%, +30%) lead to the change of the optimal lifetime T L * from 10.5 years to 8.5 years, and the estimated overall cost will go up from $5763 to $6609. Besides, the variations of repair cost C m r and the penalty cost C p l between the range (−30%, +30%) also lead to the change of the optimal lifetime from 9.5 years to 9 years. However, the effects are less than PM’s parameters in this case because the equipment can adopt some high intensive PM strategies to lower the occurances of breakdowns. However, in some cases, the equipment may cannot adopt high intensive PM strategies due to the equipment’s characteristics. In such cases, the repair and penalty costs will be the major portion of the overall cost, and the firm should try to develope effective repair techniques to improve the system’s stability for raising the efficiency of operatons or production lines.

5. Conclusions

The study aims to provide Bayesian decision analyses for preventive maintenance related to the hybrid competing failure mode. The proposed model not only can proceed with the prior analysis under a lack of historical failure data but also can consider sampling information to proceed with the posterior analysis. It can help firms to make their best PM alternatives to achieve the advantages of lower cost and operational stability. According to the analysis results of Section 4, the managerial insights and contributions can be summarized as follows: (1) The proposed Bayesian analysis will be a feasible solution in evaluating the related repair and PM costs if the firm cannot proceed with accelerating deterioration experiments to gather enough information to estimate the model’s parameters by traditional statistical methods. (2) A highly intensive PM alternative may not have a beneficial influence on reducing the related costs. It depends on whether the benefit of reducing equipment breakdowns is greater than the increment of PM cost. Therefore, the firm needs to evaluate the trade-off between the two influences before the decision-making. (3) A lower intensive PM alternative may bring serious equipment’ breakdowns during the post-phase. The firm may adopt the strategy of shortening the equipment’s lifetime to prevent this disadvantage. (4) It is suggested that the firm may adopt a medium-to-lower intensive PM alternative with a shorter lifetime if the firm encounters a budget or financial issue in practice. In other words, a firm should not adopt a highly intensive PM alternative with a longer lifetime in consideration of the payback period even if the return of the PM’s investment is high. (5) Misjudging the parameter values of the maintainable and non-maintainable modes will lead to inappropriate decisions. Therefore, the firm should be aware of the potential changes in the related costs and the optimal lifetime to react appropriately. (6) The influence of misjudging the maintainable parameters will be minor because the system’s deterioration can be restored to a younger status after a PM action. It means that the risk of misestimating maintainable parameters is less than that of non-maintainable parameters. (7) Misjudging the standard deviation of the parameters will influence the estimation of the confidence intervals of the overall cost. However, it may not result in more risky decisions since the optimal decision will be less affected by the standard deviations.
In future works, the proposed model could be refined by taking into account the deterioration of equipment or facilities on two-dimensional failure variables. Since most equipment or facilities’ deterioration depends on time and usage, considering only one of these factors could lead to inaccurate estimates. In such cases, a two-dimensional failure model would be suitable for estimating equipment or facilities’ deterioration. Accordingly, the two-factor joint probability distribution (time, usage) will be constructed by using a bivariate-Weibull probability distribution. Besides, the approaches of condition-based PM can be also applied in this scenario. Therefore, decision-makers can utilize the extended approaches to improve their PM strategies.

Author Contributions

Conceptualization, C.-C.F., C.-C.H. and J.-H.L.; Formal Analysis, C.-C.F., C.-C.H. and J.-H.L.; Funding Acquisition, C.-C.F.; Investigation, C.-C.F.; Methodology, C.-C.F., C.-C.H. and J.-H.L.; Resources, C.-C.F., C.-C.H. and J.-H.L.; Supervision, C.-C.F.; Writing—Review and Editing, C.-C.F., C.-C.H. and J.-H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was sponsored by Zhaoqing University’s Science Foundation and Guangdong Basic and Applied Basic Research Foundation, China [grant number 2020A1515010892].

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.

References

  1. Park, D.H.; Jung, G.M.; Yum, J.K. Cost minimization for periodic maintenance policy of a system subject to slow degradation. Reliab. Eng. Syst. Saf. 2000, 68, 105–112. [Google Scholar] [CrossRef]
  2. Yeh, R.H.; Lo, H.C. Optimal preventive-maintenance warranty policy for repairable products. Eur. J. Oper. Res. 2001, 134, 59–69. [Google Scholar] [CrossRef]
  3. Jung, G.M.; Park, D.H. Optimal maintenance policies during the post-warranty period. Reliab. Eng. Syst. Saf. 2003, 82, 173–185. [Google Scholar] [CrossRef]
  4. Seo, J.H.; Bai, D.S. An optimal maintenance policy for a system under periodic overhaul. Math. Comput. Model. 2004, 39, 373–380. [Google Scholar] [CrossRef]
  5. Yeh, R.H.; Chang, W.L. Optimal threshold value of failure-rate for Leased products with preventive maintenance actions. Math. Comput. Model. 2007, 46, 730–737. [Google Scholar] [CrossRef]
  6. Das, A.N.; Sarmah, A.N. Preventive replacement models: An overview and their application in process industries. Eur. J. Ind. Eng. 2010, 4, 280–307. [Google Scholar] [CrossRef]
  7. Yeh, R.H.; Kao, K.C.; Chang, W.L. Preventive-maintenance policy for leased products under various maintenance costs. Expert Syst. Appl. 2011, 38, 3558–3562. [Google Scholar] [CrossRef]
  8. Bouguera, S.; Chelbi, A.; Rezg, N. A decision model for adopting an extended warranty under different maintenance policies. Int. J. Prod. Econ. 2012, 135, 840–849. [Google Scholar] [CrossRef]
  9. Chang, W.L.; Lin, J.H. Optimal maintenance policy and length of extended warranty within the life cycle of products. Comput. Math. Appl. 2012, 63, 144–150. [Google Scholar] [CrossRef] [Green Version]
  10. Beaurepaire, P.; Valdebenito, M.A.; Schuëller, G.I.; Jensen, H.A. Reliability-based optimization of maintenance scheduling of mechanical components under fatigue. Comput. Methods Appl. Mech. Eng. 2012, 221–222, 24–40. [Google Scholar] [CrossRef]
  11. Schutz, J.; Rezg, N. Maintenance strategy for leased equipment. Comput. Ind. Eng. 2013, 66, 593–600. [Google Scholar] [CrossRef]
  12. Kim, B.S.; Ozturkoglu, Y. Scheduling a single machine with multiple preventive maintenance activities and position-based deteriorations using genetic algorithms. J. Adv. Manuf. Technol. 2013, 67, 1127–1137. [Google Scholar] [CrossRef]
  13. Khojandi, A.; Maillart, L.M.; Prokopyev, O.A. Optimal planning of life-depleting maintenance activities. IIE Trans. 2014, 46, 636–652. [Google Scholar] [CrossRef]
  14. Yuan, X.; Lu, Z. Efficient approach for reliability-based optimization based on weighted importance sampling approach. Reliab. Eng. Syst. Saf. 2014, 132, 107–114. [Google Scholar] [CrossRef]
  15. Lu, B.; Zhou, X.; Li, Y. Joint modeling of preventive maintenance and quality improvement for deteriorating single-machine manufacturing systems. Comput. Ind. Eng. 2016, 91, 188–196. [Google Scholar] [CrossRef]
  16. Wang, K.; Djurdjanovic, D. Joint Optimization of Preventive Maintenance, Spare Parts Inventory and Transportation Options for Systems of Geographically Distributed Assets. Machines 2018, 6, 55. [Google Scholar] [CrossRef] [Green Version]
  17. Zhou, Y.; Kou, G.; Xiao, H.; Peng, Y.; Alsaadi, F.E. Sequential imperfect preventive maintenance model with failure intensity reduction with an application to urban buses. Reliab. Eng. Syst. Saf. 2020, 198, 106871. [Google Scholar] [CrossRef]
  18. García, F.J.Á.; Salgado, D.R. Analysis of the Influence of Component Type and Operating Condition on the Selection of Preventive Maintenance Strategy in Multistage Industrial Machines: A Case Study. Machines 2022, 10, 385. [Google Scholar] [CrossRef]
  19. Diatte, K.; O’Halloran, B.; Van Bossuyt, D.L. The Integration of Reliability, Availability, and Maintainability into Model-Based Systems Engineering. Systems 2022, 10, 101. [Google Scholar] [CrossRef]
  20. Paulsen, J.; Cooke, R.; Nyman, R. Comparative evaluation of maintenance performance using subsurvival functions. Reliab. Eng. Syst. Saf. 1997, 58, 157–163. [Google Scholar] [CrossRef]
  21. Cooke, R.; Paulsen, J. Concepts for measuring maintenance performance and methods for analysing competing failure models. Reliab. Eng. Syst. Saf. 1997, 55, 135–141. [Google Scholar] [CrossRef]
  22. Fang, C.C.; Hsu, C.C.; Liu, J.H. The Decision-Making for the Optimization of Finance Lease with Facilities’ Two-Dimensional Deterioration. Systems 2022, 10, 210. [Google Scholar] [CrossRef]
  23. Salinas-Torres, V.H.; Pereira, C.A.B.; Tiwari, R.C. Bayesian nonparametric estimation in a series system or a competing risk model. J. Nonparametr. Stat. 2002, 14, 449–458. [Google Scholar] [CrossRef]
  24. Yousef, M.M.; Hassan, A.S.; Alshanbari, H.M.; El-Bagoury, A.H.; Almetwally, E.M. Bayesian and Non-Bayesian Analysis of Exponentiated Exponential Stress–Strength Model Based on Generalized Progressive Hybrid Censoring Process. Axioms 2022, 11, 455. [Google Scholar] [CrossRef]
  25. Wang, J.; Miao, Y. Optimal preventive maintenance policy of the balanced system under the semi-Markov model. Reliab. Eng. Syst. Saf. 2021, 213, 107690. [Google Scholar] [CrossRef]
  26. Alotaibi, R.; Nassar, M.; Ghosh, I.; Rezk, H. Elshahhat, A. Inferences of a Mixture Bivariate Alpha Power Exponential Model with Engineering Application. Axioms 2022, 11, 459. [Google Scholar] [CrossRef]
  27. Liu, T.; Zhang, L.; Jin, G.; Pan, Z. Reliability Assessment of Heavily Censored Data Based on E-Bayesian Estimation. Mathematics 2022, 10, 4216. [Google Scholar] [CrossRef]
  28. Yousef, M.M.; Hassan, A.S.; Al-Nefaie, A.H.; Almetwally, E.M. Bayesian Estimation Using MCMC Method of System Reliability for Inverted Topp–Leone Distribution Based on Ranked Set Sampling. Mathematics 2022, 10, 3122. [Google Scholar] [CrossRef]
  29. Zequeira, R.I.; Berenguer, C. Periodic imperfect preventive maintenance with two categories of competing failure modes. Reliab. Eng. Syst. Saf. 2006, 91, 460–468. [Google Scholar] [CrossRef]
  30. El-Ferik, S.; Ben-Daya, M. Age-based hybrid model for imperfect preventive maintenance. IIE Trans. 2007, 38, 365–385. [Google Scholar] [CrossRef]
  31. Kahrobaee, S.; Asgarpoor, S. A hybrid analytical-simulation approach for maintenance optimization of deteriorating equipment: Case study of wind turbines. Electr. Power Syst. Res. 2013, 104, 80–86. [Google Scholar] [CrossRef]
  32. Rafiee, K.; Feng, Q.; Coit, D.W. Condition-based maintenance for repairable deteriorating systems subject to a generalized mixed shock model. IEEE Trans. Reliab. 2015, 64, 1164–1174. [Google Scholar] [CrossRef]
  33. Zhou, X.; Wu, C.; Li, Y.; Xi, L. A preventive maintenance model for leased equipment subject to internal degradation and external shock damage. Reliab. Eng. Syst. Saf. 2016, 154, 1–7. [Google Scholar] [CrossRef]
  34. Yang, L.; Zhao, Y.; Peng, R.; Ma, X. Hybrid preventive maintenance of competing failures under random environment. Reliab. Eng. Syst. Saf. 2018, 174, 130–140. [Google Scholar] [CrossRef]
  35. Cao, Y. Modeling the effects of dependence between competing failure processes on the condition-based preventive maintenance policy. Appl. Math. Model. 2021, 99, 400–417. [Google Scholar] [CrossRef]
  36. Liu, J.; Zhuang, X.; Pang, H. Reliability and hybrid maintenance modeling for competing failure systems with multistage periods. Probabilistic Eng. Mech. 2022, 68, 103254. [Google Scholar] [CrossRef]
  37. Basílio, M.P.; Pereira, V.; Costa, H.G.; Santos, M.; Ghosh, A. A Systematic Review of the Applications of Multi-Criteria Decision Aid Methods (1977–2022). Electronics 2022, 11, 1720. [Google Scholar] [CrossRef]
  38. Huang, Y.-S.; Bier, V.M. A Natural Conjugate Prior for the Nonhomogeneous Poisson Process with a Power Law Intensity Function. Commun. Stat.-Simul. Comput. 1998, 27, 525–551. [Google Scholar] [CrossRef]
Figure 1. Timeline of the preventive maintenance model of an equipment or facility.
Figure 1. Timeline of the preventive maintenance model of an equipment or facility.
Axioms 11 00734 g001
Figure 2. Preventive maintenances between perfect recovery and imperfect recovery under a hybrid deterioration.
Figure 2. Preventive maintenances between perfect recovery and imperfect recovery under a hybrid deterioration.
Axioms 11 00734 g002
Figure 3. Heuristic algorithm for obtaining the minimal cost by setting N * and M P q * .
Figure 3. Heuristic algorithm for obtaining the minimal cost by setting N * and M P q * .
Axioms 11 00734 g003
Figure 4. Flowchart for Bayesian decision process.
Figure 4. Flowchart for Bayesian decision process.
Axioms 11 00734 g004
Figure 5. The design of the computerized information system.
Figure 5. The design of the computerized information system.
Axioms 11 00734 g005
Figure 6. The overall costs per unit and year for all PM alternatives estimated by the prior analysis.
Figure 6. The overall costs per unit and year for all PM alternatives estimated by the prior analysis.
Axioms 11 00734 g006
Figure 7. The overall costs per unit and year for all PM alternatives estimated by the prior and posterior analyses.
Figure 7. The overall costs per unit and year for all PM alternatives estimated by the prior and posterior analyses.
Axioms 11 00734 g007
Figure 8. The impacts of μ α o ,   μ α p , μ β o , and μ β p on the overall cost per unit and year.
Figure 8. The impacts of μ α o ,   μ α p , μ β o , and μ β p on the overall cost per unit and year.
Axioms 11 00734 g008
Figure 9. The Impacts of C F q , τ q , C m r , and C p l on the overall cost per unit and year.
Figure 9. The Impacts of C F q , τ q , C m r , and C p l on the overall cost per unit and year.
Axioms 11 00734 g009
Table 1. The notations.
Table 1. The notations.
T L : the lifetime of a equipment or facility.
t : the age of a equipment or facility.
x : the time interval between two PMs.
t k : the effective age of a equipment or facility before the time point of the kth PM.
t k + : the effective age of a equipment or facility after the time point of the kth PM.
α o : the scale factor of the intensity function of non-maintainable failure mode.
β o : the shape factor of the intensity function of non-maintainable failure mode.
α p : the scale factor of the intensity function of maintainable failure mode
β p : the shape factor of the intensity function of maintainable failure mode.
f ( α , β ) : the prior probability distribution of the power-law intensity function.
g ( α , β ) : the posterior probability distribution of the power-law intensity function.
δ : the age reduction factor, where δ [ 0 , 1 ] .
λ o ( t | α o , β o ) : the intensity function of non-maintainable failure mode of the system deterioration.
λ p ( t | α p , β p ) : the intensity function of maintainable failure mode of the system deterioration.
λ h ( t | α o , β o , α p , β p ) : the intensity function of the hybrid mode of the system deterioration.
N : the number of PM action during the whole system lifetime.
N m r ( ) : the expected number of performing minimal repairs of the system.
C m r : the average cost to perform a minimal repair.
C p m k : the cost to perform the kth PM.
C r p : the cost of the overall replacement of a equipment or facility.
Φ ( t r ) : the probability density function of the time for performing a minimal repair.
C p l : the penalty cost if the actual repair time over the time threshold φ
φ : the time threshold for performing a minimal repair.
C F : the base cost for a PM action, which is influenced by the degree of PM.
τ : the increasing rate of PM base cost
Table 2. The detailed information of all candidate PM alternatives.
Table 2. The detailed information of all candidate PM alternatives.
Parameters for the two categories deterioration, which were judged by experts μ α o = 0.6 ,   μ β o = 0.15 ,   σ α o = 1.25 ,   σ β o = 0.3125 ;   μ α p = 0.85 ,   μ β p = 0.2125 ,   σ α p = 1.75 ,   σ β p = 0.4375
Interval between two PM actions x = 0.5   years
PM’s Base cost of the five candidate PM alternatives C F q = { $ 780 , $ 790 , $ 800 , $ 880 , $ 890 }
Age reduction factors of the five candidate PM alternatives δ p m q = { 0.8 , 0.85 , 0.9 , 0.95 , 1.0 }
Periodically increasing rates of PM cost of the five candidate PM alternatives τ q = { 0.19 , 0.195 , 0.2 , 0.235 , 0.24 }
Replacement cost C r p = $ 20 , 000
Expected cost of performinga minimal repair C m r = $250
Penalty cost if the repair time
exceed the time limit φ
C p l = $ 90
Expected value and standard deviation of performing
a minimal repair
E ( t r ) = 5   h ,   σ ( t r ) = 3   h
The limit of tolerable waiting time
for performing a minimal repair
φ = 6.5   h
Table 3. The expected preventive cost, repair cost, penalty cost, replacement cost and overall cost per unit and year for PM Alternatives 2, 3, 4 estimated by the prior analysis.
Table 3. The expected preventive cost, repair cost, penalty cost, replacement cost and overall cost per unit and year for PM Alternatives 2, 3, 4 estimated by the prior analysis.
T L C p m q ( C F q , τ q , x , T L ) T L C m r E P r i [ N m r ( · ) ] T L C p l E [ t r | ω , η ] E P r i [ N m r ( · ) ] T L C r p T L E P r i [ C ( N | M P q ) ]
M P 2 M P 3 M P 4 M P 2 M P 3 M P 4 M P 2 M P 3 M P 4 M P 2 M P 3 M P 4
0.500026126126121921921940,00040,48040,48040,480
186788098329529128724724424020,00021,40921,41521,510
1.512071227138032531730827226625813,33315,13815,14215,280
214161440163035334032829628527510,00012,06412,06612,233
2.5157216001822379363346318304290800010,26910,26710,458
31702173319844043843643393223056667911191069319
3.51816185121294294053803593393195714831883108543
41922196022644524253973793563325000775377417993
4.52021206223924764444123993723464444734073237594
52115216025154994634284183893584000703270127301
5.52207225526345224824434374043713636680267777084
62296234727515445014574564203833333663066016925
6.52383243728655675194724754353953077650264686809
72468252629785895374864944504072857640863716729
7.52553261330906115555005124664192667634363016676
82637270032016335735145314804312500630062546646
8.52719278633116555915285494954422353627662256634
92802287134206776085415675124542222626862146637
9.5 *288329563529699626555586525465210562736212 *6654
102964304036377206435686045394762000628962236681
10.53045312437447426615816225544871905631462436718
113126320738517646785956405684991818634862726763
11.53206329039587866956086595835101739638963086815
123286337340658077126216775975201667643663506873
12.53365345641718297306346956125311600648963976936
133445353842778517476477136265421538654764497005
13.53524362143838727646607316405531481660965067077
143603370344898947816727506555641429667665677154
14.53682378545959167986857686695741379674566317233
153761386747009388156987866835851333681866987316
Table 4. The expected preventive cost, repair cost, penalty cost, replacement cost and overall cost per unit and year for PM Alternative for the prior and posterior analyses.
Table 4. The expected preventive cost, repair cost, penalty cost, replacement cost and overall cost per unit and year for PM Alternative for the prior and posterior analyses.
T L C p m q ( C F q , τ q , x , T L ) T L C m r E [ N m r ( · ) ] T L C p l E [ t r | ω , η ] E [ N m r ( · ) ] T L C r p T L E P r i [ C ( N | M P q ) ] E P o s [ C ( N | M P q ) ]
PriorPosteriorPriorPosterior
0.5026125221921240,00040,48040,464
188029128724424020,00021,41521,407
1.5122731731326626213,33315,14215,135
2144034033428528010,00012,06612,055
2.51600363354304297800010,26710,251
31733384372322312666791069083
3.51851405388339325571483108279
41960425404356339500077417702
4.52062444419372351444473237276
52160463433389363400070126956
5.52255482446404374363667776712
62347501460420385333366016525
6.52437519472435396307764686382
72526537485450406285763716274
7.52613555497466416266763016193
82700573508480426250062546134
8.52786591520495436235362256094
92871608531512445222262146069
9.5 *295662654252545421056212 *6057
10 **3040643553539463200062236056 **
10.53124661563554472190562436064
113207678574568481181862726080
11.53290695584583489173963086103
123373712594597498166763506132
12.53456730604612506160063976166
133538747614626514153864496205
13.53621764623640522148165066248
143703781633655530142965676295
14.53785798642669538137966316345
153867815651683546133366986398
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fang, C.-C.; Hsu, C.-C.; Liu, J.-H. Bayesian Statistical Method Enhance the Decision-Making for Imperfect Preventive Maintenance with a Hybrid Competing Failure Mode. Axioms 2022, 11, 734. https://doi.org/10.3390/axioms11120734

AMA Style

Fang C-C, Hsu C-C, Liu J-H. Bayesian Statistical Method Enhance the Decision-Making for Imperfect Preventive Maintenance with a Hybrid Competing Failure Mode. Axioms. 2022; 11(12):734. https://doi.org/10.3390/axioms11120734

Chicago/Turabian Style

Fang, Chih-Chiang, Chin-Chia Hsu, and Je-Hung Liu. 2022. "Bayesian Statistical Method Enhance the Decision-Making for Imperfect Preventive Maintenance with a Hybrid Competing Failure Mode" Axioms 11, no. 12: 734. https://doi.org/10.3390/axioms11120734

APA Style

Fang, C. -C., Hsu, C. -C., & Liu, J. -H. (2022). Bayesian Statistical Method Enhance the Decision-Making for Imperfect Preventive Maintenance with a Hybrid Competing Failure Mode. Axioms, 11(12), 734. https://doi.org/10.3390/axioms11120734

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