The Weibull Distribution

The Weibull distribution is one of the most widely used lifetime distributions in reliability engineering. It is a versatile distribution that can take on the characteristics of other types of distributions, based on the value of the shape parameter, [math]\displaystyle{ {\beta} \,\! }[/math]. This chapter provides a brief background on the Weibull distribution, presents and derives most of the applicable equations and presents examples calculated both manually and by using ReliaSoft's Weibull++ software.

Weibull Probability Density Function

The 3-Parameter Weibull

The 3-parameter Weibull pdf is given by:

[math]\displaystyle{ f(t)={ \frac{\beta }{\eta }}\left( {\frac{t-\gamma }{\eta }}\right) ^{\beta -1}e^{-\left( {\frac{t-\gamma }{\eta }}\right) ^{\beta }} \,\! }[/math]

where:

[math]\displaystyle{ f(t)\geq 0,\text{ }t\geq \gamma \,\! }[/math]

[math]\displaystyle{ \beta\gt 0\ \,\! }[/math]

[math]\displaystyle{ \eta \gt 0 \,\! }[/math]

[math]\displaystyle{ -\infty \lt \gamma \lt +\infty \,\! }[/math]

and:

[math]\displaystyle{ \eta= \,\! }[/math] scale parameter, or characteristic life

[math]\displaystyle{ \beta= \,\! }[/math] shape parameter (or slope)

[math]\displaystyle{ \gamma= \,\! }[/math] location parameter (or failure free life)

The 2-Parameter Weibull

The 2-parameter Weibull pdf is obtained by setting [math]\displaystyle{ \gamma=0 \,\! }[/math], and is given by:

[math]\displaystyle{ f(t)={ \frac{\beta }{\eta }}\left( {\frac{t}{\eta }}\right) ^{\beta -1}e^{-\left( { \frac{t}{\eta }}\right) ^{\beta }} \,\! }[/math]

The 1-Parameter Weibull

The 1-parameter Weibull pdf is obtained by again setting [math]\displaystyle{ \gamma=0 \,\! }[/math] and assuming [math]\displaystyle{ \beta=C=Constant \,\! }[/math] assumed value or:

[math]\displaystyle{ f(t)={ \frac{C}{\eta }}\left( {\frac{t}{\eta }}\right) ^{C-1}e^{-\left( {\frac{t}{ \eta }}\right) ^{C}} \,\! }[/math]

where the only unknown parameter is the scale parameter, [math]\displaystyle{ \eta\,\! }[/math].

Note that in the formulation of the 1-parameter Weibull, we assume that the shape parameter [math]\displaystyle{ \beta \,\! }[/math] is known a priori from past experience with identical or similar products. The advantage of doing this is that data sets with few or no failures can be analyzed.

Weibull Distribution Functions

The Mean or MTTF

The mean, [math]\displaystyle{ \overline{T} \,\! }[/math], (also called MTTF) of the Weibull pdf is given by:

[math]\displaystyle{ \overline{T}=\gamma +\eta \cdot \Gamma \left( {\frac{1}{\beta }}+1\right) \,\! }[/math]

where

[math]\displaystyle{ \Gamma \left( {\frac{1}{\beta }}+1\right) \,\! }[/math]

is the gamma function evaluated at the value of:

[math]\displaystyle{ \left( { \frac{1}{\beta }}+1\right) \,\! }[/math]

The gamma function is defined as:

[math]\displaystyle{ \Gamma (n)=\int_{0}^{\infty }e^{-x}x^{n-1}dx \,\! }[/math]

For the 2-parameter case, this can be reduced to:

[math]\displaystyle{ \overline{T}=\eta \cdot \Gamma \left( {\frac{1}{\beta }}+1\right) \,\! }[/math]

Note that some practitioners erroneously assume that [math]\displaystyle{ \eta \,\! }[/math] is equal to the MTTF, [math]\displaystyle{ \overline{T}\,\! }[/math]. This is only true for the case of: [math]\displaystyle{ \beta=1 \,\! }[/math] or:

[math]\displaystyle{ \begin{align} \overline{T} &= \eta \cdot \Gamma \left( {\frac{1}{1}}+1\right) \\ &= \eta \cdot \Gamma \left( {\frac{1}{1}}+1\right) \\ &= \eta \cdot \Gamma \left( {2}\right) \\ &= \eta \cdot 1\\ &= \eta \end{align} \,\! }[/math]

The Median

The median, [math]\displaystyle{ \breve{T}\,\! }[/math], of the Weibull distribution is given by:

[math]\displaystyle{ \breve{T}=\gamma +\eta \left( \ln 2\right) ^{\frac{1}{\beta }} \,\! }[/math]

The Mode

The mode, [math]\displaystyle{ \tilde{T} \,\! }[/math], is given by:

[math]\displaystyle{ \tilde{T}=\gamma +\eta \left( 1-\frac{1}{\beta }\right) ^{\frac{1}{\beta }} \,\! }[/math]

The Standard Deviation

The standard deviation, [math]\displaystyle{ \sigma _{T}\,\! }[/math], is given by:

[math]\displaystyle{ \sigma _{T}=\eta \cdot \sqrt{\Gamma \left( {\frac{2}{\beta }}+1\right) -\Gamma \left( {\frac{1}{ \beta }}+1\right) ^{2}} \,\! }[/math]

The Weibull Reliability Function

The equation for the 3-parameter Weibull cumulative density function, cdf, is given by:

[math]\displaystyle{ F(t)=1-e^{-\left( \frac{t-\gamma }{\eta }\right) ^{\beta }} \,\! }[/math]

This is also referred to as unreliability and designated as [math]\displaystyle{ Q(t) \,\! }[/math] by some authors.

Recalling that the reliability function of a distribution is simply one minus the cdf, the reliability function for the 3-parameter Weibull distribution is then given by:

[math]\displaystyle{ R(t)=e^{-\left( { \frac{t-\gamma }{\eta }}\right) ^{\beta }} \,\! }[/math]

The Weibull Conditional Reliability Function

The 3-parameter Weibull conditional reliability function is given by:

[math]\displaystyle{ R(t|T)={ \frac{R(T+t)}{R(T)}}={\frac{e^{-\left( {\frac{T+t-\gamma }{\eta }}\right) ^{\beta }}}{e^{-\left( {\frac{T-\gamma }{\eta }}\right) ^{\beta }}}} \,\! }[/math]

or:

[math]\displaystyle{ R(t|T)=e^{-\left[ \left( {\frac{T+t-\gamma }{\eta }}\right) ^{\beta }-\left( {\frac{T-\gamma }{\eta }}\right) ^{\beta }\right] } \,\! }[/math]

These give the reliability for a new mission of [math]\displaystyle{ t \,\! }[/math] duration, having already accumulated [math]\displaystyle{ T \,\! }[/math] time of operation up to the start of this new mission, and the units are checked out to assure that they will start the next mission successfully. It is called conditional because you can calculate the reliability of a new mission based on the fact that the unit or units already accumulated hours of operation successfully.

The Weibull Reliable Life

The reliable life, [math]\displaystyle{ T_{R}\,\! }[/math], of a unit for a specified reliability, [math]\displaystyle{ R\,\! }[/math], starting the mission at age zero, is given by:

[math]\displaystyle{ T_{R}=\gamma +\eta \cdot \left\{ -\ln ( R ) \right\} ^{ \frac{1}{\beta }} \,\! }[/math]

This is the life for which the unit/item will be functioning successfully with a reliability of [math]\displaystyle{ R\,\! }[/math]. If [math]\displaystyle{ R = 0.50\,\! }[/math], then [math]\displaystyle{ T_{R}=\breve{T} \,\! }[/math], the median life, or the life by which half of the units will survive.

The Weibull Failure Rate Function

The Weibull failure rate function, [math]\displaystyle{ \lambda(t) \,\! }[/math], is given by:

[math]\displaystyle{ \lambda \left( t\right) = \frac{f\left( t\right) }{R\left( t\right) }=\frac{\beta }{\eta }\left( \frac{ t-\gamma }{\eta }\right) ^{\beta -1} \,\! }[/math]