Weibull modulus explained

The Weibull modulus is a dimensionless parameter of the Weibull distribution. It represents the width of a probability density function (PDF) in which a higher modulus is a characteristic of a narrower distribution of values. Use case examples include biological and brittle material failure analysis, where modulus is used to describe the variability of failure strength for materials.

Definition

The Weibull distribution, represented as a cumulative distribution function (CDF), is defined by:

F(x)=1-\exp(-(x-xu
x0

)m)

in which m is the Weibull modulus.[1]

x0

is a parameter found during the fit of data to the Weibull distribution and represents an input value for which ~67% of the data is encompassed. As m increases, the CDF distribution more closely resembles a step function at

x0

, which correlates with a sharper peak in the probability density function (PDF) defined by:
f(x)=(m)(
x0
x-xu
x0

)m-1\exp(-(

x-xu
x0

)m)

Failure analysis often uses this distribution,[2] as a CDF of the probability of failure F of a sample, as a function of applied stress σ, in the form:
F(\sigma)=1-\exp\left[-(\sigma
\sigma0

)m\right]

Failure stress of the sample, σ, is substituted for the

x

property in the above equation. The initial property

xu

is assumed to be 0, an unstressed, equilibrium state of the material.

In the plotted figure of the Weibull CDF, it is worth noting that the plotted functions all intersect at a stress value of 50 MPa, the characteristic strength for the distributions, even though the value of the Weibull moduli vary. It is also worth noting in the plotted figure of the Weibull PDF that a higher Weibull modulus results in a steeper slope within the plot.

The Weibull distribution can also be multi-modal, in which there would be multiple reported

x0

values and multiple reported moduli, m. The CDF for a bimodal Weibull distribution has the following form,[3] when applied to materials failure analysis:
F(\sigma)=1-\phi\exp[(-(\sigma
\sigma01
m1
)

]-(1-\phi)\exp[-(

\sigma
\sigma02
m2
)

]

This represents a material which fails by two different modes. In this equation m1 is the modulus for the first mode, and m2 is the modulus for the second mode. Φ is the fraction of the sample set which fail by the first mode. The corresponding PDF is defined by:

f(\sigma)=\phi(m1)(
\sigma01
\sigma
\sigma01
m1-1
)\exp[-(
\sigma
\sigma01
m1
)]+(1-\phi)(
m2)(
\sigma02
\sigma
\sigma02
m2-1
)\exp[-(
\sigma
\sigma02
m2
)

]

Examples of a bimodal Weibull PDF and CDF are plotted in the figures of this article with values of the characteristic strength being 40 and 120 MPa, the Weibull moduli being 4 and 10, and the value of Φ is 0.5, corresponding to 50% of the specimens failing by each failure mode.

Linearization of the CDF

The compliment of the cumulative Weibull distribution function can be expressed as:

P=1-F

Where P corresponds to the probability of survival of a specimen for a given stress value. Thus, it follows that:

P(\sigma)=1-\left[1-\exp\left[-(

\sigma
\sigma0
m\right]\right]=\exp\left[-(\sigma
\sigma0
)

)m\right]

where m is the Weibull modulus. If the probability is plotted vs the stress, we find that the graph is sigmoidal, as shown in the figure above. Taking advantage of the fact that the exponential is the base of the natural logarithm, the above equation can be rearranged to:

ln\left[ln\left(

1
1-F

\right)\right]=ln\left[\left(

\sigma
\sigma0

\right)m\right]

Which using the properties of logarithms can also be expressed as:

ln\left[ln\left(

1
1-F

\right)\right]=mln(\sigma)-mln(\sigma0)

When the left side of this equation is plotted as a function of the natural logarithm of stress, a linear plot can be created which has a slope of the Weibull modulus, m, and an x-intercept of

ln(\sigma0)

.

Looking at the plotted linearization of the CDFs from above it can be seen that all of the lines intersect the x-axis at the same point because all of the functions have the same value of the characteristic strength. The slopes vary because of the differing values of the Weibull moduli.

Measurement

Standards organizations have created multiple standards for measuring and reporting values of Weibull parameters, along with other statistical analyses of strength data:

When applying a Weibull distribution to a set of data the data points must first be put in ranked order. For the use case of failure analysis specimens' failure strengths are ranked in ascending order, i.e. from lowest to greatest strength. A probability of failure is then assigned to each failure strength measured, ASTM C1239-13 uses the following formula:

F(\sigma)=i-0.5
N

where

i

is the specimen number as ranked and

N

is the total number of specimens in the sample. From there

F

can be plotted against failure strength to obtain a Weibull CDF. The Weibull parameters, modulus and characteristic strength, can be obtained from fitting or using the linearization method detailed above.

Example uses from published work

Weibull statistics are often used for ceramics and other brittle materials.[8] They have also been applied to other fields as well such as meteorology where wind speeds are often described using Weibull statistics.[9] [10] [11]

Ceramics and brittle materials

For ceramics and other brittle materials, the maximum stress that a sample can be measured to withstand before failure may vary from specimen to specimen, even under identical testing conditions. This is related to the distribution of physical flaws present in the surface or body of the brittle specimen, since brittle failure processes originate at these weak points. Much work has been done to describe brittle failure with the field of linear elastic fracture mechanics and specifically with the development of the ideas of the stress intensity factor and Griffith Criterion. When flaws are consistent and evenly distributed, samples will behave more uniformly than when flaws are clustered inconsistently. This must be taken into account when describing the strength of the material, so strength is best represented as a distribution of values rather than as one specific value.

Consider strength measurements made on many small samples of a brittle ceramic material. If the measurements show little variation from sample to sample, the calculated Weibull modulus will be high, and a single strength value would serve as a good description of the sample-to-sample performance. It may be concluded that its physical flaws, whether inherent to the material itself or resulting from the manufacturing process, are distributed uniformly throughout the material. If the measurements show high variation, the calculated Weibull modulus will be low; this reveals that flaws are clustered inconsistently, and the measured strength will be generally weak and variable. Products made from components of low Weibull modulus will exhibit low reliability and their strengths will be broadly distributed. With careful manufacturing processes Weibull moduli of up to 98 have been seen for glass fibers tested in tension.[12]

A table is provided with the Weibull moduli for several common materials. However, it is important to note that the Weibull modulus is a fitting parameter from strength data, and therefore the reported value may vary from source to source. It also is specific to the sample preparation and testing method, and subject to change if the analysis or manufacturing process changes.

Table of Weibull Moduli for Common Materials!Material!Weibull Modulus
Silicon Nitride20[13]
Silicon Carbide18
Aluminum Nitride10
Aluminum Oxide (98%)10
Aluminum Oxide (99%)10
Zirconia Toughened Alumina (ZTA)13
Boron Carbide12
Titanium Diboride11
Zirconia15

Organic materials

Studies examining organic brittle materials highlight the consistency and variability of the Weibull modulus within naturally occurring ceramics such as human dentin and abalone nacre. Research on human dentin[14] samples indicates that the Weibull modulus remains stable across different depths or locations within the tooth, with an average value of approximately 4.5 and a range between 3 and 6. Variations in the modulus suggest differences in flaw populations between individual teeth, thought to be caused by random defects introduced during specimen preparation. Speculation exists regarding a potential decrease in the Weibull modulus with age due to changes in flaw distribution and stress sensitivity. Failure in dentin typically initiates at these flaws, which can be intrinsic or extrinsic in origin, arising from factors such as cavity preparation, wear, damage, or cyclic loading.

Studies on the abalone shell illustrate its unique structural adaptations, sacrificing tensile strength perpendicular to its structure to enhance strength parallel to the tile arrangement. The Weibull modulus of abalone nacre samples[15] is determined to be 1.8, indicating a moderate degree of variability in strength among specimens.

Quasi-brittle materials

The Weibull modulus of quasi-brittle materials correlates with the decline in the slope of the energy barrier spectrum, as established in fracture mechanics models. This relationship allows for the determination of both the fracture energy barrier spectrum decline slope and the Weibull modulus, while keeping factors like crack interaction and defect-induced degradation in consideration. Temperature dependence and variations due to crack interactions or stress field interactions are observed in the Weibull modulus of quasi-brittle materials. Damage accumulation leads to a rapid decrease in the Weibull modulus, resulting in a right-shifted distribution with a smaller Weibull modulus as damage increases.[16]

Quality analysis

Weibull analysis is also used in quality control and "life analysis"[17] for products. A higher Weibull modulus allows for companies to more confidently predict the life of their product for use in determining warranty periods.

Other methods of characterization for brittle materials

A further method to determine the strength of brittle materials has been described by the Wikibook contribution Weakest link determination by use of three parameter Weibull statistics.

Notes and References

  1. Weibull . Waloddi . 1951 . A Statistical Distribution Function of Wide Applicability . Journal of Applied Mechanics. 18 . 3 . 293 . 10.1115/1.4010337 . 1951JAM....18..293W .
  2. Book: Chiang . Yet-ming . Physical ceramics . Birnie . Dunbar P. . Kingery . W. D. . 1997 . J. Wiley . 978-0-471-59873-2 . MIT series in materials science & engineering . New York.
  3. Loidl . Dieter . Paris . Oskar . Rennhofer . H. . Müller . Martin . Peterlik . Herwig . November 2007 . Skin-core structure and bimodal Weibull distribution of the strength of carbon fibers . Carbon . 45 . 14 . 2801–2805 . 10.1016/j.carbon.2007.09.011 . 2007Carbo..45.2801L . 0008-6223.
  4. 2013 . "ASTM C1239-13: Standard Practice for Reporting Uniaxial Strength Data and Estimating Weibull Distribution Parameters for Advanced Ceramics" . ASTM . ASTM Compass.
  5. 2021 . Standard Practice for Reporting Uniaxial Strength Data and Estimating Weibull Distribution Parameters for Advanced Graphites . ASTM . ASTM Compass.
  6. March 2019 . ISO 20501:2019 Fine ceramics (advanced ceramics, advanced technical ceramics) - Weibull statistics for strength data . International Organization for Standardization . Iso.org.
  7. 2007 . Advanced Technical Ceramics - Mechanical Properties of Monolithic Ceramics at Room Temperature - Part 5: Statistical Analysis . American National Standards Institute.
  8. Quinn . JB . Quinn . GD . February 2010 . A practical and systematic review of Weibull statistics for reporting strengths of dental materials . Dental Materials . 26 . 2 . 135–147 . 10.1016/j.dental.2009.09.006 . 19945745 . 3086645 .
  9. Conradsen . K. . Nielsen . L. B. . Prahm . L. P. . 1984 . Review of Weibull Statistics for Estimation of Wind Speed Distributions . Journal of Climate and Applied Meteorology . 23 . 8 . 1173–1183 . 10.1175/1520-0450(1984)023<1173:ROWSFE>2.0.CO;2 . 26181388 . 1984JApMe..23.1173C . 0733-3021. free .
  10. Pavia . Edgar G. . O'Brien . James J. . 1986-10-01 . Weibull Statistics of Wind Speed over the Ocean . Journal of Applied Meteorology and Climatology . EN . 25 . 10 . 1324–1332 . 10.1175/1520-0450(1986)025<1324:WSOWSO>2.0.CO;2 . 1520-0450 . free. 1986JApMe..25.1324P .
  11. Carta . J.A. . Ramirez . P. . 2007 . Analysis of two-component mixture Weibull statistics for estimation of wind speed distributions . Renewable Energy . 32 . 3 . 518–531 . 10.1016/j.renene.2006.05.005 . Elsevier Science Direct.
  12. Matthewson . M. John . Kurkjian . Charles R. . Gulati . Suresh T. . November 1986 . Strength Measurement of Optical Fibers by Bending . Journal of the American Ceramic Society . en . 69 . 11 . 815–821 . 10.1111/j.1151-2916.1986.tb07366.x . 0002-7820.
  13. Book: Meyers . Marc . Mechanical Behavior of Materials . Chawla . Krishan . Cambridge University Press . 1998 . 978-0-521-86675-0 . 2nd . Cambridge, United Kingdom.
  14. Quinn . Janet B. . Quinn . George D. . February 2010 . A practical and systematic review of Weibull statistics for reporting strengths of dental materials . Dental Materials . 26 . 2 . 135–147 . 10.1016/j.dental.2009.09.006 . 0109-5641 . 3086645 . 19945745.
  15. Lin . Albert Yu-Min . Meyers . Marc André . December 2009 . Interfacial shear strength in abalone nacre . Journal of the Mechanical Behavior of Biomedical Materials . 2 . 6 . 607–612 . 10.1016/j.jmbbm.2009.04.003 . 19716105 . 1751-6161.
  16. Qi . Chengzhi . Lu . Chunsheng . Chanyshev . A.I. . Li . Xiaozhao . Qu . Xiaolei . June 2023 . Preliminary study on the determination of the Weibull modulus of strength distribution in quasi-brittle materials . Geohazard Mechanics . 1 . 2 . 103–109 . 10.1016/j.ghm.2023.05.002 . 2949-7418. free . 2023GeohM...1..103Q .
  17. Web site: 2018-02-20 . Weibull Analysis Quality-One . 2024-04-27 . en-US.