Windshield failure time.
Stochastic models for spare parts forecasting have not been widely researched in scientific literature from the aspect of their reliability. In this chapter, the authors present models which analyze standard reliability parameters of technical systems’ parts/components. By analyzing system reliability and failure rate, we estimate the required number of spare parts in the moment of expected failure or when reliability falls below the predefined level. Two different approaches based on data availability are presented herewith.
- spare parts
- forecasting methods
- Rayleigh model
- Weibull model
- failure rate
Technical systems, for example, aircrafts or weapon systems, are typically highly complex and composed of a large number of components and parts. Maintaining those systems adheres to strict rules and procedures for specific part or component. In order for such a system to be successfully maintained, the efficient management of spare part inventory is required, that is, a specific part should be provided in the right place at the right time .
Unpredictability of future events, that is, equipment and parts failure, has major impact on this problem. One way to reduce the level of unpredictability is to maintain a sufficient number of spare parts in inventory which results in the increase of warehouse costs and capital trapped in spare parts; another way is the assessment of spare parts inventory by using one of the available models for spare parts forecasting .
Forecasting is an essential skill on predicting future events and represents a foundation for every valid assessment. Although forecasts are often deemed expensive and time consuming, a vast number of researchers have been involved in finding novel methods with improved and more accurate assessment in recent years.
Concerning the demand for spare parts, exponential smoothing model and Croston’s model are still among the most utilized due to their simplicity. The exponential smoothing model described in Refs. [3, 4] is based on pounder-moving average and can be implemented straightforwardly. In 1972, Croston presented a model based on exponential smoothing but far superior to it . This method was most widely used in industry and is still a part of many software packages for spare parts forecasting. Rao  and then Schultz  studied Croston’s method and proposed certain alterations but with no effect on execution results. Willemain et al. , and later Johnston and Boylan  proved that for majority of cases Croston’s method gives better results than exponential smoothing method. Syntetos and Boyland [10, 11] made a considerable contribution to this scientific field by proposing a modified method based on criticism of Croston’s model, calling it biased with reference to spare parts demand per time unit.
Another commonly used method is Bootstrap method—a computer statistical technique which, based on available data sample, creates a large number of new samples of same range as the original sample, by random sampling with replacement from original data set. From the aspect of inventory management, this model was examined in detail in Refs. [12, 13]. In addition, Refs. [14, 15] examine spare parts forecasting using Poisson’s model. A conclusion was reached that traditional statistical methods based on analysis of time series can incorrectly assess functional form related to dependent and independent variables.
The nature of the spare parts demand is stochastic, and previously stated models do not always provide most accurate assessments . For that reason, the number of models dealing with assessing the required number of spare parts, based on parameters such as spare part reliability, maintenance possibilities, life span, maintenance costs, and so on, has greatly increased during the past decade.
Due to the increase of system complexity , presenting reliability as a quantitative measure is commonly considered. By analyzing system’s reliability and failure rate, we can estimate the required number of spare parts in the moment of expected failure or when reliability falls below the predefined level. There are numerous papers on the topic of determining the required number of spare parts, particularly as a part of logistic support [18, 19]. Refs. [20, 21] mostly deal with repairable spare parts or inventory managements with the aim of achieving previously set system reliability. On the other side, quantitative theories based on the theory of reliability were used to estimate the failure rate in order to precisely determine demand for a specific spare part [22–25]. An overview of the abovementioned models in spare parts forecasting with some new approaches proposed has been given in Ref. .
In this chapter, we present two models for spare parts forecasting based on the analysis of reliability parameter. Each model can be used depending on data availability.
2. Spare parts forecasting using Rayleigh model
Spare parts manufacturers provide only basic information on the part they produce. In this case, we observe the average life span of a part/component (
Based on the aforementioned, the PDF of Rayleigh’s model can now be presented as follows:
Cumulative density function (CDF) can similarly be determined as
If we substitute in the previous equation, it can be reduced to
The function of part’s reliability can further be determined with Eq. (7)
Finally, based on the previous equations, we can define the failure function as a probability that the examined part will cease to perform its function in a specific time interval
The essence of this model’s implementation is in an approach of determining the quantity of spare parts in inventory and costs that occur due to negative level of inventory at the end of usability period of examined part (underage costs) originally presented in Ref. . This approach stresses the stochastic nature of variable
By implementing and substituting , the previous equation is reduced to
where is the Gauss random variable with variance .
Now, given that
Moreover, the number of spare parts required in inventory can be determined by observing the total time when random variable
This model can be a part of software for spare parts forecasting as in Ref. .
3. Spare parts forecasting using Weibull’s model
For spare parts forecasting in cases where the total unit time of a product is unavailable but data on previous failure rates are available, Weibull’s model can be used . This model is most often used when reliability of a technical system is being determined. The PDF of Weibull’s model with two parameters, shape parameter
The function of reliability can be expressed by the following formula:
while CDF is
The failure rate function of Weibull’s model is
In relation to Weibull’s model, a problem can arise while estimating the distribution parameters. There are multiple ways for parameter estimation but based on Ref.  if the available data sample is less than 15, then linear regression for parameter estimation is used, and in the opposite case maximum likelihood estimator (MLE) gives best results.
If linear regression is used, we start from CDF function of Weibull’s model expressed as follows:
Taking the logarithm of previous formula, we get
If we take the logarithm of the previous expression once again
If are observed values of random variable
Data should be arranged in an ascending order:
Value of empirical function of distribution is then calculated: .
Drawing points in the coordinate system and mathematically determine direction which best approximates obtained data set. The most frequently used is the least-square method or linear regression which minimizes the vertical deviation of observed points from the given direction, in the following way:
where and are the least-squares estimates of
Based on the equations, it is now simple to determine parameters
In cases where data sample is greater than 15, parameters are determined by maximum likelihood estimator method .
Taking natural logarithm of both sides gives
By differentiating partially previous equation, we obtain
From the previous equation, we can estimate
Eq. (30) does not have a closed-form solution, so we can estimate shape parameter
Once the parameters of Weibull’s model are determined, we can estimate the required number of spare parts that should be in inventory in the given time interval. In order to achieve that, we use the approach presented in Ref. . The PDF of Weibull’s distributed failure time is given by Eq. (13), while the PDF of Rayleigh’s distribute failure time is given by Eq. (1). Based on those two equations, we can conclude that , while , wherein
If we stress the stochastic nature of a specific part’s failure rate by observing the expected number of variations of Rayleigh’s random variable in interval (
Based on the previous equations, the number of spare parts that are subject to failure in time period of [0,
As in the previous section, it is now necessary to determine the number of spare parts required to be kept in inventory in time interval [0,
When we determine the quantity of spare parts that need to be kept in inventory, then in case when we know the price per unit of product based on Newsvendor model , we can determine the underage costs, that is, cost per unit of unsatisfied demand as
4. Case study
Depending on data availability, one of the two recommended approaches will be implemented. If the total unit time (
Table 1 provides data on specific aircraft model’s Windshield, taken from paper . Table 1 consists of 88 records on the part’s failure time and 65 records on censored time expressed in flight hours. Censored time (or service time) means that the Windshields have not failed at the time of observation.
Weibull probability plot for data from previous table is presented in Figure 1.
As noted in Figure 1, records from Table 1 follow Weibull’s distribution, and as their number is greater than 15, based on the previous research, we use maximum likelihood estimation method for the estimation of parameters in Weibull’s distribution, described extensively in the previous section.
Implementation of this method results in the value of shape parameter
Now, when values of these parameters are known we can determine the reliability of part windshield based on Eq. (15), while its failure function can be determined based on Eq. (16). A graphical representation of these functions is presented in Figures 2 and 3.
Based on Figure 4, it can be concluded that it is necessary to have one spare part in inventory after 2000 flight hours. The results of this analysis, that is, data on spare part reliability, failure time, and required quantity in inventory can be of great value to decision makers, on questions related to what to keep in inventory and in what quantity, and when to plan maintenance activities in order to prevent the occurrence of failure.
The aforementioned models provide possibility for taking into consideration underage costs during decision-making process by using Eq. (37), for those cases when price of the part in issue is available. These costs are difficult to determine objectively. Some consequences of lack of spare parts, for example, damage to company’s reputation due to delays, are difficult to express quantitatively.
This section elaborates on determining spare parts required to be kept in inventory from the aspect of reliability analysis of that part. Depending on data availability, either Rayleigh’s or Weibull’s method can be used. In case of Weibull’s method, two approaches for the assessment of parameters are given, depending on data sample size. Information obtained with this analysis can have a major role in the process of supply management. Based on them, it is possible to reduce costs that occur due to delays, unplanned cancellations, and so on. These models can serve as a solid foundation for the creation of software for spare parts forecasting. Although the emphasis was placed on planning maintenance activities and avoiding delays, all the aforementioned leads to limiting the consequences of suboptimal supply management, that is, minimize the spare parts overstocking. In the case when we are dealing with reparable systems, researches should focus on the determination or minimization of the repair rate of such system.