Strength of Materials/Weakest link determination by use of three parameter Weibull statistics

When for a particular property a material shows a large scatter in measured data (e.g. strength measurement results of brittle materials, or of threads with important variations in thickness as listed in reference [1]), a characteristic value for, or the weakest link of such a property can be obtained by analyzing the measured data on the basis of three parameter Weibull statistics.

Introduction edit

When a material is subject for strength determination, series of strength measurements will show a scatter. When the scatter is low compared with the average value of the series, a minimum strength value   may be defined by:


in which   is the average strength,   is the standard deviation of the measurement series, and   is the factor related to the probability of not failing (e.g. relates to a fail probability of 0.0001)

However, a number of materials, e.g. threads with important variations in thickness, brittle materials as annealed glass, will show a large, or even a very large scatter compared with the average value of the series. And even the scatter of the average values of a number of series of measurements may show such an important scatter compared with the overall average value, that applying expression (1) leads to a nonsensical minimum strength.

In those cases, strength measurements series may be described by two parameter Weibull distributions that indicate that the minimum value approaches to zero, or:


in which   is the cumulated probability of failing,   the measured strength value,   the first or shape parameter, and   the second or scale parameter.

The minimum applicable strength value   can be determined by the following expression:


in which the cumulated probability   of failing has a low value, e.g. 0.0001.

When the scatter of strength measurements is extremely important, or the minimum strength value is certainly larger than zero, straight application of the two parameter Weibull distribution is of no meaning and should be replaced by a more complex method leading to the use of three parameter Weibull distributions, expressed as:


in which   is the third parameter, or the minimum value of the strength. Or in other words, the weakest link.

It should also be taken into account that the parameters shall be considered as distributions   of random variables with a mean value and a variance:

(5) 1st, or shape parameter:  

(6) 2nd, or scale parameter:  

(7) 3rd, or location parameter:  

Quantification method edit

In order to quantify the three parameters, and especially the third one, the weakest link, a relative large number of data should be available, grouped in series. Those series have to be obtained from measurement pieces drawn from lots, where each lot should be considered as more or less homogeneous. The estimation of the parameters then may follow the following iterative procedure:

step i: Use the three parameter Weibull distribution, however declare the third parameter  . Put each series of data into that distribution by means of a best fitted method. This gives for each series the parameters   and  . For a best fitting method, see reference [2], and possibly the Wikipedia contribution "Weibull modulus".

step ii: With those parameters   and  , calculate a specific value and range, e.g. the average   value and the range   with possibly   and  .

step iii: Place all pairs of specific range/specific value, here as example the  /  pairs, into a Cartesian diagram with as abscise the specific range and as ordinate the specific value. The result should be a stretched cloud pointing to the ordinate.

step iv: Calculate the trend, or regression line and calculate the specific value when the specific range equals to zero. The value obtained,  , is also the average third parameter  .

step v: Calculate the best estimation of the residual standard deviation  , needed for the estimation of the lower value of the third parameter, the real weakest link  , with the expressions:

(8)   and (9)  

where   is the square sum of the average   values,   is the number of pairs, and   is the correlation factor.

step vi: Calculate the lower value of the third parameter by means of


where   is the eccentricity of the normal distribution, e.g.   for a probability of 0.05.

step vii: Replace the used third parameter in step i by the   value estimated in step vi. Repeat steps i to vii until the difference between last and fore last value is sufficient small.

Conditions for representation edit

In order to select a most representative number of strength measurement series of the entire, more or less inhomogeneous population, but ensuring within each lot a most representative homogeneous number of measurement pieces, some conditions should be taken into account.

a) The number of series should represent a sufficient time interval of production of the concerned material. A population consisting out of up to five series is certainly too small, a population from six to ten series may lead to doubtful results.

b) The number of measurement pieces per series should be recommended to have at least 10 pieces. The test pieces for one series should be selected randomly from one lot manufactured within a production interval that may be assumed to be as homogeneous as possible.

c) The measurement pieces per series should be of one size: the measurement pieces of different series should differ in size.

Examples of possible minimum numbers of measurement series:

Threads: From each of e.g. four production intervals or lots, spread over a production period of one month, where each production of threads is as homogeneous as possible, a first series of tensile measurements should consist out of measurement pieces of length  , a second series should consist out of measurement pieces of 3 times   and a third series should consist out of measurement pieces of 10 times  . So in total 12 series of tensile measurement pieces.

Annealed flat glass: From each of e.g. four production intervals or lots, spread over a production period of some months, where each production of annealed flat glass is as homogeneous as possible (e.g. series could be taken from one very large sheet of glass), a first series should consist out of measurement pieces of surface  , a second series should consist out of measurement pieces of surface  , and a third series should consist out of measurement pieces of surface  . For various surfaces, refer to the series of international standards CEN ISO 1288 (reference [3]). So in total 12 series

Computer generated example edit

Figure 1

An example with computer generated data is given below.

Measurement pieces:

- 4 lots, where each lot contains 3348 strength numbers;

- each lot deliver 3 series of measurement pieces, where the pieces contain 121, 25 or 9 strength numbers, in respectively the first, second and third series;

- each series consists out of 12 measurement pieces;

- total 4 x 3 x 12 = 144 measurement pieces;

- strength numbers   are randomly distributed following the expression


where   is a random variable related to lots with values between 0 and 1, and   and   are random variables related to any stress number, also with values between 0 and 1;

Third parameter estimation:

The steps i to vii have been performed two times. The result is graphically given in figure 1. The stability of the method has been investigated by repeating the whole process 20 times. The averaged minimum third parameter value   was 16.2 by a standard deviation  =2.2.

Background edit

Symbols edit

  = occurrence, e.g. breaking stress

  = subscript that indicates mean value

  = a dimension, e.g. length, surface, ..

  = probability that an occurrence occurs

  = intensity of irregularities that influences, most negatively, an occurrence

  = range, or value of difference between two occurrences

  = variance of within a series of values

  = first parameter of Weibull distributions

  = second or scale parameter of Weibull distributions

  = third parameter of Weibull distributions

Weibull distribution with parameters subject of random variation edit

The following sections give the information for the development from expressions in the “Introduction” to those in step v and step vi in the “Quantification method”.

From expression (4) the stress   belonging to a certain occurrence probability F is


From expression (i) and expressions (5) to (7) it can be concluded that

at first the variance of the stress with a determined probability to occur, is

(ii)   where + means "combined with".

When the probability   is selected in such a way that   , then expression (i) will be simplified to

(iii)   which can be detailed as

(iv.a) mean value   and

(iv.b) variance  

and secondly the value of the stress   with a determined probability   to occur varies linear with the value of the second parameter   . When the value of   can be reduced to zero, the remaining expression will be

(v)   which can be detailed as

(vi.a) mean value   and

(vi.b) variance  

Properties of the second or scale parameter   edit

Failing of an object under stress of inhomogeneous objects such as threads, or brittle materials such as annealed glass, depends on the size   of the objects as well as the intensity   of those irregularities in or on the object from where failing starts. For instance for threads the size is characterized by the length, the intensity of the irregularities by a.o. variations in thickness. In the case of annealed glass the size is characterized by the surface of the object, the intensity of the irregularities by the presence and depth of flaws on that surface.

This dependency can be expressed as


so that the expression (6) can be more detailed by


This expression demonstrates that when the size   approaches the infinitive and, or the irregularities are such that   approaches zero, the value of the second parameter   approaches zero.

Reduction of second parameter   to zero edit

For a certain measurement series, the occurrence   occurring with a probability  , and the occurrence   occurring with a probability  , are

(ix)   and


When stating that the range   for that particular measurement series is


then the second parameter   can be expressed as


When substituting   in expression (ix) by expression (xii), the result is that


Further, when taking into account the variability of   by using expression (viii), it can be concluded that the occurrence   with a probability of   is


This means that when the size   becomes an infinitive large, or when the intensity of the irregularities are so important that   approaches zero, the occurrence equals to the value of the third parameter, or


to be detailed as

(xvi.a) mean value:  

(xvi.b) and variance:  

Conclusion edit

On the basis of (xiii), the measured values of pairs   and   can be placed in a Cartesian diagram, which should give a scatter of points where the scatter is determined by the distributions of  ,   and  . When the specimens per measurement series do not differ in dimension and, or when the series are drawn from lots with no significant difference in the intensity of irregularities, the scatter will be an undetermined cloud of points. No conclusion can be drawn. When on the other hand these differences have been respected, the scatter allows to draw a regression line of   as function of the range  , and there where the range   equals to zero, to apply expression (xv), and thus the quantification method as exposed above.

References edit

[1] Strength of Materials/Unsorted topics

[2] Typical glass type bending strength - Jacob Zwart - International Glass review, Flat Glass, Issue 3, 1999.

[3] European and International Standard series EN ISO 1288: Glazing in building - Determination of the bending strength of glass