Open access peer-reviewed chapter

Spare Parts Forecasting Based on Reliability

By Nataša Kontrec and Stefan Panić

Submitted: November 27th 2016Reviewed: May 5th 2017Published: December 20th 2017

DOI: 10.5772/intechopen.69608

Downloaded: 1074


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.


  • reliability
  • spare parts
  • forecasting methods
  • Rayleigh model
  • Weibull model
  • failure rate
  • inventory

1. Introduction

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 [1].

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 [2].

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 [5]. This method was most widely used in industry and is still a part of many software packages for spare parts forecasting. Rao [6] and then Schultz [7] studied Croston’s method and proposed certain alterations but with no effect on execution results. Willemain et al. [8], and later Johnston and Boylan [9] 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 [16]. 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 [17], 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 [2225]. An overview of the abovementioned models in spare parts forecasting with some new approaches proposed has been given in Ref. [26].

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 (Tut) expressed in hours, μ, as a stochastic process modeled with Rayleigh’s distribution [27]. Probability density function (PDF) of Rayleigh’s model is stated in Eq. (1)


where σis the parameter of Rayleigh’s distribution determined by relation E(μ2) = 2σ2. E(μ) is the mathematical expectation of Rayleigh’s random variable μ. Based on that, Tut can be presented as follows:


With replacement μ2/(2σ2) = x, Eq. (2) is transformed into


Γ(a) is Gamma function [28] and Γ(3/2)=π/2, so the average life span can be presented as Tut=σπ2.

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 μ2π/(4Tut2)=uin 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. [29]. This approach stresses the stochastic nature of variable Tut. By observing the expected number of variations of this random variable in time interval (Tut + dTut) for a given slope T˙utin its specific environment dTut for Nunits of stated component, we can then determine the number of components most likely to fail as follows:


By implementing and substituting T˙ut2/(2σ˙2)=u, the previous equation is reduced to


where T˙utis the Gauss random variable with variance V(T˙ut)=σ˙.

Now, given that μis a Rayleigh’s random variable with mathematical expectation E(μ) = Tut and variance V(μ) = 2 Tut2/π, then the average number of components n, which will be subject to failure in time Tut, can be determined as


Moreover, the number of spare parts required in inventory can be determined by observing the total time when random variable μis below Tut:


This model can be a part of software for spare parts forecasting as in Ref. [30].

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 [31]. 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 βand characteristic life η, has the following form:


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. [32] 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


and implement



x= lnt, a= βand b= −βlnη, then Eq. (19) can be expressed in the following linear form:


If t1,t2,,tnare observed values of random variable t, the procedure for parameter estimation is as follows [33]:

  1. Data should be arranged in an ascending order: t(1)t(2)t(n).

  2. Value of empirical function of distribution is then calculated: F^(t(i)),i=1,2,,n.

  3. Calculating yi=ln[ln[1F^(t(i))]],i=1,2,,n.

  4. Calculating xi=lnt(i),i=1,2,,n.

  5. Drawing points Ti=(xi,yi)in the coordinate system and mathematically determine direction y=β^x+b^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 a^and b^are the least-squares estimates of aand b, and can be determined with the following formulae:




  • Based on the equations, it is now simple to determine parameters β i η.

  • In cases where data sample is greater than 15, parameters are determined by maximum likelihood estimator method [33].

    By observing t1, t2, t3,…, tn data sample of size npopulation with given PDF equation (13), then the joint density of likelihood function can be determined as a product of the densities of each point


    Taking natural logarithm of both sides gives


    By differentiating partially previous equation, we obtain




    From the previous equation, we can estimate ηas


    By replacing Eq. (29) into Eq. (27), we get


    Eq. (30) does not have a closed-form solution, so we can estimate shape parameter βby using Newton-Raphson’s method or any other numeric procedure. After we determined β^, by replacing it into Eq. (29) we can calculate η^.

    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. [34]. 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 σ=η2, while μ=tβ2, wherein μis Rayleigh’s random and tis Weibull’s random variable.


    |J| represents Jacobian transformation expressed with Eq. (32)


    By substituting Eq. (32) in Eq. (31), we obtain Eq. (33)


    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 (μ, μ+) for a given slope μ˙in specific environment dμ, then the number of spare parts that are most likely to fail can be determined as follows:


    Based on the previous equations, the number of spare parts that are subject to failure in time period of [0, t] is


    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, t]. We will use the approach presented in paper [29] but in this case the quantity will be calculated as a quotient of CDF function of Weibull’s random variable and size n. More accurately, as parameter ηmarks the time in which 63.2% of units will fail and is approximately equal to MTTF [35], we estimate the time interval in which tis below ηas


    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 [36], we can determine the underage costs, that is, cost per unit of unsatisfied demand as


    where Φ-1 presents inverse distribution function (complementary error function), cu are underage costs, and co are the overage costs, which in our case is spare part price.

    4. Case study

    Depending on data availability, one of the two recommended approaches will be implemented. If the total unit time (Tut) is available for a specific part, we opt to use model for spare parts forecasting based on Rayleigh’s distribution. When historic failure/censored data are available, we use model based on Weibull’s distribution.

    Table 1 provides data on specific aircraft model’s Windshield, taken from paper [37]. 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.


    Table 1.

    Windshield failure time.

    Weibull probability plot for data from previous table is presented in Figure 1.

    Figure 1.

    Weibull probability plot forTable 1data.

    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 β= 2.28, and characteristic life parameter η= 3450.54.

    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.

    Figure 2.

    Windshield reliability.

    Figure 3.

    Failure rate function of windshield.

    Finally, we can determine the quantity of spare parts required to be kept in inventory in a given time interval based on Eq. (36). A graphical representation is presented in Figure 4.

    Figure 4.

    Quantity of spare parts in function of time.

    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.

    5. Conclusions

    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.

    © 2017 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

    How to cite and reference

    Link to this chapter Copy to clipboard

    Cite this chapter Copy to clipboard

    Nataša Kontrec and Stefan Panić (December 20th 2017). Spare Parts Forecasting Based on Reliability, System Reliability, Constantin Volosencu, IntechOpen, DOI: 10.5772/intechopen.69608. Available from:

    chapter statistics

    1074total chapter downloads

    1Crossref citations

    More statistics for editors and authors

    Login to your personal dashboard for more detailed statistics on your publications.

    Access personal reporting

    Related Content

    This Book

    Next chapter

    Model Development for Reliability Cannibalization

    By Bernard Tonderayi Mangara

    Related Book

    First chapter

    Microassembly Using Water Drop

    By Taksehi Mizuno

    We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

    More About Us