Weibull modulus

Last updated

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.

Contents

Definition

CDF of Weibull distribution for the example of predicting failure in materials, s0= 50 MPa Weibull CDF.png
CDF of Weibull distribution for the example of predicting failure in materials, σ0= 50 MPa

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

in which m is the Weibull modulus. [1] 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 , which correlates with a sharper peak in the probability density function (PDF) defined by:

PDF of Weibull distribution for the example of predicting failure in materials, s0= 50 MPa Weibull PDF.png
PDF of Weibull distribution for the example of predicting failure in materials, σ0= 50 MPa
PDF of a bimodal Weibull distribution with Weibull Moduli of 4 and 10 and characteristic strengths of 40 and 120 MPa Bimodal Weibull PDF.png
PDF of a bimodal Weibull distribution with Weibull Moduli of 4 and 10 and characteristic strengths of 40 and 120 MPa

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:

Failure stress of the sample, σ, is substituted for the property in the above equation. The initial property 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 values and multiple reported moduli, m. The CDF for a bimodal Weibull distribution has the following form, [3] when applied to materials failure analysis:

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:

CDF of a bimodal Weibull distribution with Weibull Moduli of 4 and 10 and characteristic strengths of 40 and 120 MPa Bimodal Webull CDF.png
CDF of a bimodal Weibull distribution with Weibull Moduli of 4 and 10 and characteristic strengths of 40 and 120 MPa

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:

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

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:

The linearization of the Weibull CDFs shown above. Linearization of Weibull CDF.png
The linearization of the Weibull CDFs shown above.

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

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 .

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 [4] uses the following formula:

where is the specimen number as ranked and is the total number of specimens in the sample. From there 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] [9] They have also been applied to other fields as well such as meteorology where wind speeds are often described using Weibull statistics. [10] [11] [12]

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

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
MaterialWeibull Modulus
Silicon Nitride 20 [9]
Silicon Carbide 18 [9]
Aluminum Nitride 10 [9]
Aluminum Oxide (98%)10 [9]
Aluminum Oxide (99%)10 [9]
Zirconia Toughened Alumina (ZTA)13 [9]
Boron Carbide 12 [9]
Titanium Diboride 11 [9]
Zirconia 15 [9]

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.

Related Research Articles

<span class="mw-page-title-main">Normal distribution</span> Probability distribution

In probability theory and statistics, a normal distribution or Gaussian distribution is a type of continuous probability distribution for a real-valued random variable. The general form of its probability density function is

<span class="mw-page-title-main">Pareto distribution</span> Probability distribution

The Pareto distribution, named after the Italian civil engineer, economist, and sociologist Vilfredo Pareto, is a power-law probability distribution that is used in description of social, quality control, scientific, geophysical, actuarial, and many other types of observable phenomena; the principle originally applied to describing the distribution of wealth in a society, fitting the trend that a large portion of wealth is held by a small fraction of the population. The Pareto principle or "80-20 rule" stating that 80% of outcomes are due to 20% of causes was named in honour of Pareto, but the concepts are distinct, and only Pareto distributions with shape value of log45 ≈ 1.16 precisely reflect it. Empirical observation has shown that this 80-20 distribution fits a wide range of cases, including natural phenomena and human activities.

<span class="mw-page-title-main">Log-normal distribution</span> Probability distribution

In probability theory, a log-normal (or lognormal) distribution is a continuous probability distribution of a random variable whose logarithm is normally distributed. Thus, if the random variable X is log-normally distributed, then Y = ln(X) has a normal distribution. Equivalently, if Y has a normal distribution, then the exponential function of Y, X = exp(Y), has a log-normal distribution. A random variable which is log-normally distributed takes only positive real values. It is a convenient and useful model for measurements in exact and engineering sciences, as well as medicine, economics and other topics (e.g., energies, concentrations, lengths, prices of financial instruments, and other metrics).

<span class="mw-page-title-main">Weibull distribution</span> Continuous probability distribution

In probability theory and statistics, the Weibull distribution is a continuous probability distribution. It models a broad range of random variables, largely in the nature of a time to failure or time between events. Examples are maximum one-day rainfalls and the time a user spends on a web page.

In engineering, deformation refers to the change in size or shape of an object. Displacements are the absolute change in position of a point on the object. Deflection is the relative change in external displacements on an object. Strain is the relative internal change in shape of an infinitesimal cube of material and can be expressed as a non-dimensional change in length or angle of distortion of the cube. Strains are related to the forces acting on the cube, which are known as stress, by a stress-strain curve. The relationship between stress and strain is generally linear and reversible up until the yield point and the deformation is elastic. Elasticity in materials occurs when applied stress does not surpass the energy required to break molecular bonds, allowing the material to deform reversibly and return to its original shape once the stress is removed. The linear relationship for a material is known as Young's modulus. Above the yield point, some degree of permanent distortion remains after unloading and is termed plastic deformation. The determination of the stress and strain throughout a solid object is given by the field of strength of materials and for a structure by structural analysis.

<span class="mw-page-title-main">Stress–strain curve</span> Curve representing a materials response to applied forces

In engineering and materials science, a stress–strain curve for a material gives the relationship between stress and strain. It is obtained by gradually applying load to a test coupon and measuring the deformation, from which the stress and strain can be determined. These curves reveal many of the properties of a material, such as the Young's modulus, the yield strength and the ultimate tensile strength.

In mathematics, a Gaussian function, often simply referred to as a Gaussian, is a function of the base form

<span class="mw-page-title-main">Fracture</span> Split of materials or structures under stress

Fracture is the appearance of a crack or complete separation of an object or material into two or more pieces under the action of stress. The fracture of a solid usually occurs due to the development of certain displacement discontinuity surfaces within the solid. If a displacement develops perpendicular to the surface, it is called a normal tensile crack or simply a crack; if a displacement develops tangentially, it is called a shear crack, slip band, or dislocation.

<span class="mw-page-title-main">Compressive strength</span> Capacity of a material or structure to withstand loads tending to reduce size

In mechanics, compressive strength is the capacity of a material or structure to withstand loads tending to reduce size. In other words, compressive strength resists compression, whereas tensile strength resists tension. In the study of strength of materials, tensile strength, compressive strength, and shear strength can be analyzed independently.

Thermal shock is a phenomenon characterized by a rapid change in temperature that results in a transient mechanical load on an object. The load is caused by the differential expansion of different parts of the object due to the temperature change. This differential expansion can be understood in terms of strain, rather than stress. When the strain exceeds the tensile strength of the material, it can cause cracks to form and eventually lead to structural failure.

In probability theory and statistics, the generalized extreme value (GEV) distribution is a family of continuous probability distributions developed within extreme value theory to combine the Gumbel, Fréchet and Weibull families also known as type I, II and III extreme value distributions. By the extreme value theorem the GEV distribution is the only possible limit distribution of properly normalized maxima of a sequence of independent and identically distributed random variables. Note that a limit distribution needs to exist, which requires regularity conditions on the tail of the distribution. Despite this, the GEV distribution is often used as an approximation to model the maxima of long (finite) sequences of random variables.

Estimation theory is a branch of statistics that deals with estimating the values of parameters based on measured empirical data that has a random component. The parameters describe an underlying physical setting in such a way that their value affects the distribution of the measured data. An estimator attempts to approximate the unknown parameters using the measurements. In estimation theory, two approaches are generally considered:

The Anderson–Darling test is a statistical test of whether a given sample of data is drawn from a given probability distribution. In its basic form, the test assumes that there are no parameters to be estimated in the distribution being tested, in which case the test and its set of critical values is distribution-free. However, the test is most often used in contexts where a family of distributions is being tested, in which case the parameters of that family need to be estimated and account must be taken of this in adjusting either the test-statistic or its critical values. When applied to testing whether a normal distribution adequately describes a set of data, it is one of the most powerful statistical tools for detecting most departures from normality. K-sample Anderson–Darling tests are available for testing whether several collections of observations can be modelled as coming from a single population, where the distribution function does not have to be specified.

<span class="mw-page-title-main">Flexural strength</span> Material property

Flexural strength, also known as modulus of rupture, or bend strength, or transverse rupture strength is a material property, defined as the stress in a material just before it yields in a flexure test. The transverse bending test is most frequently employed, in which a specimen having either a circular or rectangular cross-section is bent until fracture or yielding using a three-point flexural test technique. The flexural strength represents the highest stress experienced within the material at its moment of yield. It is measured in terms of stress, here given the symbol .

Expected shortfall (ES) is a risk measure—a concept used in the field of financial risk measurement to evaluate the market risk or credit risk of a portfolio. The "expected shortfall at q% level" is the expected return on the portfolio in the worst of cases. ES is an alternative to value at risk that is more sensitive to the shape of the tail of the loss distribution.

<span class="mw-page-title-main">Nakagami distribution</span> Statistical distribution

The Nakagami distribution or the Nakagami-m distribution is a probability distribution related to the gamma distribution. It is used to model physical phenomena, such as those found in medical ultrasound imaging, communications engineering, meteorology, hydrology, multimedia, and seismology.

In financial mathematics, tail value at risk (TVaR), also known as tail conditional expectation (TCE) or conditional tail expectation (CTE), is a risk measure associated with the more general value at risk. It quantifies the expected value of the loss given that an event outside a given probability level has occurred.

<span class="mw-page-title-main">Log-Cauchy distribution</span>

In probability theory, a log-Cauchy distribution is a probability distribution of a random variable whose logarithm is distributed in accordance with a Cauchy distribution. If X is a random variable with a Cauchy distribution, then Y = exp(X) has a log-Cauchy distribution; likewise, if Y has a log-Cauchy distribution, then X = log(Y) has a Cauchy distribution.

<span class="mw-page-title-main">Size effect on structural strength</span> Deviation with the scale in the theories of elastic or plastic structures

According to the classical theories of elastic or plastic structures made from a material with non-random strength (ft), the nominal strength (σN) of a structure is independent of the structure size (D) when geometrically similar structures are considered. Any deviation from this property is called the size effect. For example, conventional strength of materials predicts that a large beam and a tiny beam will fail at the same stress if they are made of the same material. In the real world, because of size effects, a larger beam will fail at a lower stress than a smaller beam.

The four-point flexural test provides values for the modulus of elasticity in bending , flexural stress , flexural strain and the flexural stress-strain response of the material. This test is very similar to the three-point bending flexural test. The major difference being that with the addition of a fourth bearing the portion of the beam between the two loading points is put under maximum stress, as opposed to only the material right under the central bearing in the case of three-point bending.

References

  1. Weibull, Waloddi (1951). "A Statistical Distribution Function of Wide Applicability" (PDF). Journal of Applied Mechanics. 18 (3): 293. Bibcode:1951JAM....18..293W. doi:10.1115/1.4010337.
  2. Chiang, Yet-ming; Birnie, Dunbar P.; Kingery, W. D. (1997). Physical ceramics. MIT series in materials science & engineering. New York: J. Wiley. ISBN   978-0-471-59873-2.
  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. Bibcode:2007Carbo..45.2801L. doi:10.1016/j.carbon.2007.09.011. ISSN   0008-6223.
  4. 1 2 ""ASTM C1239-13: Standard Practice for Reporting Uniaxial Strength Data and Estimating Weibull Distribution Parameters for Advanced Ceramics"". ASTM. 2013 via ASTM Compass.
  5. "Standard Practice for Reporting Uniaxial Strength Data and Estimating Weibull Distribution Parameters for Advanced Graphites". ASTM. 2021 via ASTM Compass.
  6. "ISO 20501:2019 Fine ceramics (advanced ceramics, advanced technical ceramics) - Weibull statistics for strength data". International Organization for Standardization. March 2019 via Iso.org.
  7. "Advanced Technical Ceramics - Mechanical Properties of Monolithic Ceramics at Room Temperature - Part 5: Statistical Analysis". American National Standards Institute. 2007.
  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. doi:10.1016/j.dental.2009.09.006. PMC   3086645 . PMID   19945745.
  9. 1 2 3 4 5 6 7 8 9 10 Meyers, Marc; Chawla, Krishan (1998). Mechanical Behavior of Materials (2nd ed.). Cambridge, United Kingdom: Cambridge University Press. ISBN   978-0-521-86675-0.
  10. 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. Bibcode:1984JApMe..23.1173C. doi:10.1175/1520-0450(1984)023<1173:ROWSFE>2.0.CO;2. ISSN   0733-3021. JSTOR   26181388.
  11. Pavia, Edgar G.; O'Brien, James J. (1986-10-01). "Weibull Statistics of Wind Speed over the Ocean". Journal of Applied Meteorology and Climatology. 25 (10): 1324–1332. Bibcode:1986JApMe..25.1324P. doi: 10.1175/1520-0450(1986)025<1324:WSOWSO>2.0.CO;2 . ISSN   1520-0450.
  12. 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. doi:10.1016/j.renene.2006.05.005 via Elsevier Science Direct.
  13. Matthewson, M. John; Kurkjian, Charles R.; Gulati, Suresh T. (November 1986). "Strength Measurement of Optical Fibers by Bending". Journal of the American Ceramic Society. 69 (11): 815–821. doi:10.1111/j.1151-2916.1986.tb07366.x. ISSN   0002-7820.
  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. doi:10.1016/j.dental.2009.09.006. ISSN   0109-5641. PMC   3086645 . PMID   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. doi:10.1016/j.jmbbm.2009.04.003. ISSN   1751-6161. PMID   19716105.
  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. Bibcode:2023GeohM...1..103Q. doi: 10.1016/j.ghm.2023.05.002 . ISSN   2949-7418.
  17. "Weibull Analysis | Quality-One". 2018-02-20. Retrieved 2024-04-27.