Modelling and Simulation in Engineering

Volume 2018, Article ID 1241856, 12 pages

https://doi.org/10.1155/2018/1241856

## Reliability-Centered Maintenance: Analyzing Failure in Harvest Sugarcane Machine Using Some Generalizations of the Weibull Distribution

Institute of Mathematical and Computer Sciences, São Paulo University, São Carlos, SP, Brazil

Correspondence should be addressed to Pedro L. Ramos; rb.psu@somarlordep

Received 7 December 2017; Revised 22 January 2018; Accepted 14 February 2018; Published 1 April 2018

Academic Editor: Farouk Yalaoui

Copyright © 2018 Pedro L. Ramos et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We considered five generalizations of the standard Weibull distribution to describe the lifetime of two important components of sugarcane harvesting machines. The harvesters considered in the analysis harvest an average of 20 tons of sugarcane per hour and their malfunction may lead to major losses; therefore, an effective maintenance approach is of main interest for cost savings. For the considered distributions, mathematical background is presented. Maximum likelihood is used for parameter estimation. Further, different discrimination procedures were used to obtain the best fit for each component. At the end, we propose a maintenance scheduling for the components of the harvesters using predictive analysis.

#### 1. Introduction

The arrival of the sugarcane culture in Brazil has had a significant impact on the national economy, which led the country to become the largest producer in the world [1]. Its subproducts are used in the food and chemical industries, as well as in electricity generation and fuel production. Mechanized harvesting is one of the most important stages in the sugar and ethanol mills, since it provides the raw material with quality, time, and competitive costs for later processing. Among the used machines in the mechanized harvest, the harvesters stand out for having a large number of corrective stops, given the functionality in such extreme environmental conditions. In addition, its operation is in a regime of 24 hours on the workdays, having impact on fatigue and wear of their parts. During operation, the harvester processes an average of 20 tons of sugarcane per hour and its malfunction may lead to major losses; therefore, an effective maintenance approach is of keen interest [2].

Reliability-centered maintenance consists of determining the most effective maintenance approach [3, 4]. This process was firstly developed in the aviation industry for deciding what maintenance work is needed to keep aircraft airborne, driven by the need to improve reliability, while reducing the cost of maintenance [5]. Reliability analysis can be used to estimate time-related parameters to the next machine stop [6], providing information to manage and control the preventive maintenance of harvesters which could result in increased production and has potential for cost savings.

In reliability, common procedures are usually based on the assumption that the data follows a Weibull distribution. Introduced by Weibull [7], this distribution has convenient mathematical properties and its physiological failure process arises in many areas (see Manton and Yashin [8]). Additionally, McCool [9] provided an extensive discussion about its use in reliability. However, this distribution cannot be used to describe data with nonmonotone hazard function (bathtub, upside-down bathtub, to list a few). To overcome this problem, many generalizations of the standard Weibull distribution have been proposed. Murthy et al. [10] presented the application of some generalized Weibull distributions for modeling complex failure data sets. Pham and Lai [11] discussed recent generations of Weibull-related lifetime. Further, Lai [12] reviewed more than 25 generalizations of the Weibull distribution, and Tahir and Cordeiro [13] cited more than 30 compounded Weibull models.

In this paper, we consider five important generalized Weibull distributions with three parameters to describe the lifetime of two important components of the sugarcane harvesting machines. Our main goal here is to correctly predict the next failure of the components, not to present an extensive review of the generalizations of Weibull distribution. In the cited papers, the authors only described mathematical properties of the distributions and conducted the fit for different data. However, reliability is more about correctly predicting the future than describing the past; in this sense, no predictive analysis was presented considering such generalizations.

The distributions considered are the gamma-Weibull distribution [14], generalized Weibull (GW) distribution [15], exponentiated Weibull (EW) distribution [16], Marshall-Olkin Weibull (MOW) distribution [17], and the extended Poisson-Weibull (EPW) distribution [18]. While the first three distributions are the most common three parameter generalization of the Weibull, the MOW and the EPW arise in the competitive and complementary risk scenario (see Louzada [19], for a detailed discussion). In these cases, the latent variables follow, respectively, a geometric and a zero-truncated Poisson distribution and each of components in risk came from a Weibull baseline distribution.

For each distribution, the mathematical background is reviewed and the parameters estimators are presented using the maximum likelihood estimators. Further, different discrimination procedures are used to obtain the best fit for each component. At the end, we propose a maintenance scheduling for the components of the harvesters using predictive analysis.

The remainder of this paper is organized as follows. Section 2 presents the literature review related to the survival models adopted. Section 3 exposes the data collection and empirical analysis, as well as carrying out the predictive analysis based on the parametric models. Finally, in Section 4, we present some final remarks related to the contribution of this study.

#### 2. Theoretical Background

In this section, we present the statistical background on the adopted distributions and its parameter estimation procedures. The following distributions are considered: gamma-Weibull, generalized Weibull, exponentiated Weibull, Marshall-Olkin-Weibull, and Marshall-Olkin-Weibull. Their choice is based on their flexibility to accommodate lifetime dataset with hazard functions with different shapes, for instance, constant, increasing, decreasing, bathtub, and upside-down bathtub.

##### 2.1. The Gamma-Weibull Distribution

Introduced by Stacy [14], the gamma-Weibull distribution with three parameters is a flexible model for reliability data due to its ability to accommodate various forms of the hazard function. This distribution is also known as generalized gamma (GG) distribution as it generalizes the two-parameter gamma distribution; hereafter, we will refer to this model as GG distribution to avoid confusion with the GW distribution. A random variable has GG distribution if its probability density function (PDF) is given bywhere , , and . The mean and variance of GG are given bySome relevant distributions are special cases such as the Weibull distribution (when ), the distribution gamma (), log-normal (case limit when ), and the generalized normal distribution (). For example, the generalized normal distribution is also a distribution that includes several distributions known as half-normal (), Rayleigh (), Maxwell-Boltzmann () and chi (). The cumulative distribution function (CDF) is given bywhere is the lower incomplete gamma function. The survival function iswhere is the upper incomplete gamma function.

The hazard function of the GG distribution iswhere the hazard function has constant, increasing, decreasing, bathtub, and upside-down bathtub hazard rate.

For parameter estimation, let be a random sample of size , where . Then, the likelihood function related to the PDF (1) is given by

The log-likelihood is given by

Setting the partial derivatives , , and equal to , we obtain the following maximum likelihood estimators:where . The solution provides the maximum likelihood estimates (MLEs). See, for instance, Ramos et al. [20, 21] and Achcar et al. [22] for a detailed discussion.

Under mild conditions, the estimators become unbiased for large samples and asymptotically efficient. Moreover, such estimators have asymptotically normal joint distribution given bywhere is the Fisher information matrix: that is,and is the trigamma function.

##### 2.2. The Generalized Weibull Distribution

Introduced by Mudholkar et al. [15], the generalized Weibull distribution has PDF given bywhere , , and . The CDF and the survival function are, respectively, given by

The hazard function of the GW distribution is

This model is very flexible to describe lifetime data, since it has the hazard function with constant, increasing, decreasing, bathtub, and upside-down bathtub hazard rate. The quantile function of the GW distribution has closed form and is given by

For parameter estimation, let be a random sample of size , where . Then, the likelihood function related to the PDF (11) is given by

The log-likelihood is given by

Setting the partial derivatives equal to , we obtain the maximum likelihood estimators. Here, we follow Mudholkar et al. [15] which considers the direct maximization of (16). Under mild conditions, the obtained estimators are consistent and efficient with an asymptotically normal joint distribution given bywhere is the Fisher information matrix associated with the vector of parameters and is the Fisher information elements in and given by

Since the Fisher information matrix does not have closed-form expression for some terms, an alternative is to consider the observed information matrix, where the terms are given by

Hereafter, we considered the same approach to obtain the confidence intervals for the parameters from other distributions.

##### 2.3. The Exponentiated Weibull Distribution

Introduced by Mudholkar et al. [16], the exponentiated Weibull distribution with PDF is given bywhere , and .

The exponentiated Weibull distribution includes the Weibull distribution () and the exponentiated exponential distribution (). The survival function is given by

The hazard function of the GG distribution is

The shapes of the hazard function are analogous to the GG and GW distribution. Additionally, the quantile function of the EW distribution has closed form and is given by

The th moment of the EW distribution is given bywhere , . The proof of this equality is presented by Choudhury [23].

For parameter estimation, let be a random sample of size , where . Then, the likelihood function related to the PDF (20) is given by

The log-likelihood is given by

Setting the partial derivatives , , and equal to , we obtain the following maximum likelihood estimators:

##### 2.4. The Marshall-Olkin-Weibull Distribution

Marshall and Olkin [17] presented a new procedure for introducing an additional parameter into a family of distribution. In this case, the authors applied such procedure in the Weibull distribution. The obtained PDF of the MOW distribution is given bywhere , , and . The MOW distribution arises naturally in competing risks scenarios. Let , where is a random variable with geometrical distribution and are assumed to be independent and identically distributed according to a Weibull distribution; then the has a PDF given by (28). Cordeiro and Lemonte [24] derived many properties and the parameter estimators for the MOW distribution; the following results were obtained from the cited work. The survival function is given by

The hazard function of the MOW distribution iswhere its behavior is constant, increasing, decreasing, bathtub, and unimodal. Moreover, the quantile function of the MOW distribution has closed form and is given by

For parameter estimation, let be a random sample of size , where . Then, the likelihood function related to the PDF (28) is given by

The log-likelihood is given by

Setting the partial derivatives , , and equal to , we obtain the following maximum likelihood estimators:for more details, see Cordeiro and Lemonte [24].

##### 2.5. The Extended Poisson-Weibull Distribution

Ramos et al. [18] introduced the extended Poisson-Weibull (EPW) distribution as a generalization of Weibull-Poisson distribution (see Hemmati et al. [25]), where its PDF is given bywhere , , and . Analogously to the MOW distribution, the EWP model arises naturally in competing risks scenarios. Let , where is a random variable with a zero-truncated Poisson distribution and are assumed to be independent and identically distributed according to a Weibull distribution, then the has a PDF given by (35). The survival function is given by

The hazard function of the GG distribution is

For the EWP distribution, the hazard function has different shapes such as constant, increasing, decreasing, bathtub, and upside-down bathtub. Furthermore, the quantile function of the EPW distribution has closed form and is given by

For parameter estimation, let be a random sample of size , where . Then, the likelihood function related to the PDF (20) is given by

The log-likelihood is given by

Setting the partial derivatives , , and equal to , we obtain the following maximum likelihood estimators:

##### 2.6. Goodness of Fit

Firstly, in order to verify the behavior of the empirical data, the Total Time on Test plot (TTT-plot) was considered (Barlow and Campo [26]). The TTT-plot is obtained through the plot of , where, , and is the ordered data. For data with concave (convex) curve, the hazard function has increasing (decreasing) shape. If the behavior starts convex and then becomes concave (concave and then convex), the hazard function has bathtub (inverse bathtub) shape.

The goodness of fit is checked considering the Kolmogorov-Smirnov (KS) test. This procedure is based on the KS statistic , where is the supremum of the set of distances, is the empirical distribution function, and is CDF. A hypothesis test is conducted at the level of significance to test whether or not the data comes from . In this case, the null hypothesis is rejected if the returned value is smaller than .

The following discrimination criterion methods were adopted: Akaike information criteria (AIC) and the corrected AIC (AICc) computed, respectively, by and , where is the number of parameters to be fitted and is MLEs of . For a set of candidate models for , the best one provides the minimum values.

#### 3. Data Collection and Empirical Analysis

The dataset came from two sources: a manual stop system, which brings the history of revisions and corrective stops of two sugarcane harvesters; and data from the onboard computers of the harvesters, which provide information on the operation of the machine. The data were collected from January 2015 to August 2017, a period corresponding to 2.5 harvests (crops), that is, a period of thirty months of activity.

##### 3.1. Empirical Analysis

Firstly, considering all the stops and their reasons, records of the performance of the predictive maintenance are required to be observed. In total, 1347 stops were observed, of which 186 were preventive and 1161 corrective stops. Thus, it is possible to observe the superior amount of unplanned stops, thus questioning the effectiveness of preventive maintenance. Table 1 shows the failure among the harvests, considering both machines analysis.