MCEN90014 · Materials Engineering
Statistical Analysis: Fracture, Bending & Weibull
This is the Week 3 reliability core of MCEN90014 Materials Engineering at the University of Melbourne, where brittle strength stops being a single number and becomes a distribution. It joins linear-elastic fracture mechanics (K_IC), the bending test used to measure ceramic strength, and the Weibull statistics that describe the scatter — all on the subject's Process→Structure→Property spine. Get it right and the ceramics-and-reliability marks in the final exam become a matter of picking the correct formula and carrying units cleanly.
What this chapter covers
- 01Explain why brittle strength is flaw-controlled and therefore scatters, unlike a metal's yield stress
- 02State the fracture criterion K = Yσ√(πa) with correct units (MPa·√m) and the geometry factor Y
- 03Rearrange for the critical crack a_cr = (1/π)(K_IC/Yσ)² and the critical stress σ_cr = K_IC/(Y√(πa))
- 04Distinguish a centre crack (Y ≈ 1, half-length a) from a surface/edge crack (Y ≈ 1.12, full depth a)
- 05Compute bending stress σ = Mc/I, with M_max = PL/4 in three-point bending
- 06Use I = πd⁴/64 for a round bar (giving σ = 32M/πd³) and I = bh³/12 for a rectangular bar
- 07Read a normal distribution's survival fractions at μ−1σ, μ−2σ, μ−3σ
- 08Fit a Weibull distribution: P_f = 1 − exp[−(σ/σ₀)^m], with modulus m and characteristic strength σ₀
- 09Estimate P_f = n/(N+1), linearise to ln[ln(1/(1−P_f))] = m·lnσ − m·lnσ₀, and read m as the slope
Critical fracture stress of a flawed alumina bar
- +1Choose the right rearrangement. Setting K = K_IC in K = Yσ√(πa) and solving for stress gives σ_cr = K_IC / (Y√(πa)); for a surface flaw Y = 1.12.
- +1Convert the crack depth to metres: a = 20 μm = 20×10⁻⁶ m (keep SI so K_IC in MPa·√m divided by √m returns MPa).
- +1Evaluate the root (π INSIDE): √(πa) = √(π × 20×10⁻⁶) = √(6.28×10⁻⁵) = 7.93×10⁻³ √m.
- +1Form the denominator: Y√(πa) = 1.12 × 7.93×10⁻³ = 8.88×10⁻³ √m.
- +1Divide: σ_cr = 3.5 / 8.88×10⁻³ = 3.94×10² MPa ≈ 394 MPa.
- +1Interpret: below ≈ 394 MPa the 20 μm flaw is stable; above it the crack runs and the part fails without warning. A real design would sit well under this to allow for scatter.
Key terms
- Fracture toughness (K_IC)
- A material constant giving the crack-tip driving force at which fast fracture starts, K_IC = Yσ√(πa) at failure, in MPa·√m. A higher K_IC means the material tolerates a larger flaw or a higher stress before a crack runs.
- Stress-intensity factor (K)
- A measure of how hard the crack tip is being driven by the remote stress, K = Yσ√(πa). The crack propagates when K reaches K_IC; below that it stays put.
- Geometry factor (Y)
- A dimensionless factor for the crack/loading geometry: Y ≈ 1 for a through-thickness centre crack of total length 2a in a wide plate (a = half-length), Y ≈ 1.12 for a surface/edge crack of full depth a.
- Critical crack length (a_cr)
- The flaw size that just triggers fracture at a given stress, a_cr = (1/π)(K_IC/Yσ)². Any flaw longer than a_cr runs; anything shorter is stable — it sets the detection limit for inspection.
- Flexural strength (modulus of rupture)
- The outer-fibre bending stress at fracture, σ = Mc/I, used to strength-test brittle materials in a three-point (M = PL/4) or four-point bend, in MPa. Usually higher than the direct tensile strength of the same ceramic.
- Second moment of area (I)
- The section's resistance to bending, in m⁴: I = πd⁴/64 for a solid round bar (diameter d) and I = bh³/12 for a rectangle (b×h). Note it is /64, not /32 — the 32 belongs to σ = 32M/(πd³).
- Weibull modulus (m)
- The dimensionless shape parameter of the Weibull strength distribution; it is the slope of ln[ln(1/(1−P_f))] against lnσ. A higher m means less scatter and a more reliable material (ceramics m ≈ 5–20; metals m ≈ 50–100).
- Characteristic strength (σ₀)
- The scale parameter in P_f = 1 − exp[−(σ/σ₀)^m]; the stress at which P_f = 1 − 1/e ≈ 0.632 (about 63% have failed). It is a reference marker, not the mean strength.
Statistical Analysis: Fracture, Bending & Weibull FAQ
Why are ceramics tested in bending instead of pulling them in tension?
Ceramics are brittle and flaw-controlled: their strength is set by the largest crack in the loaded volume, and they are far weaker in tension than in compression. Gripping a brittle bar in a tensile machine tends to crack it at the grips and gives wildly scattered results. A three-point bend test loads the bottom face in tension in a repeatable way and, combined with a Weibull treatment of the scatter, yields a usable flexural strength — which is why this chapter pairs the bending formula σ = Mc/I with the statistics.
What is the difference between K_IC and flexural strength?
They answer different questions. Fracture toughness K_IC (MPa·√m) is a deterministic material constant that, together with a known flaw size and the criterion K = Yσ√(πa), tells you whether a specific crack will run. Flexural strength (MPa) is the outer-fibre stress at which a bend specimen actually broke — a whole-specimen result that reflects that specimen's worst flaw, so it scatters and is described statistically with a Weibull modulus. Toughness is about a named crack; flexural strength is about the batch.
Can AI help me with fracture toughness and Weibull analysis in MCEN90014?
Yes — Sia is an AI tutor that explains the method step by step: it can walk you through rearranging K = Yσ√(πa) for the critical crack or stress, check that you kept the √π inside the root and the crack length in metres, or set out the Weibull table (rank → P_f = n/(N+1) → the double log → slope = m) so you understand each line. It is a study aid for learning the reasoning, not a source of ready-made answers, and it will not sit an assessment or guarantee a grade. Always follow the University of Melbourne's academic-integrity and generative-AI rules for MCEN90014 and confirm what is permitted on Canvas.
Studying with AI? Sia — free AI mechanical engineering tutor works through MCEN90014 step by step.
Exam move
Build this chapter around three formulas you can write from a blank page: the fracture criterion K = Yσ√(πa) (and its two rearrangements for a_cr and σ_cr), the bending stress σ = Mc/I with M = PL/4 and I = πd⁴/64, and the Weibull failure law P_f = 1 − exp[−(σ/σ₀)^m] with its linearised form. Because a formula sheet is provided in the final exam, spend practice time on choosing and using the right relation — matching the geometry factor Y to a centre versus surface crack, keeping the √π inside the root, converting crack and beam dimensions to metres — rather than memorising equations. Drill the sign and direction habits until they are automatic: a longer crack raises K, a higher K_IC or a smaller flaw raises the safe stress, and a higher Weibull m means less scatter. For the Weibull fit, always rank the data and use P_f = n/(N+1) (never n/N), then take the double log and read m as the slope, not the intercept. The final exam is 10 questions of 10 marks each (100 marks, all compulsory) worth 50% of the subject with an exam hurdle, so treat every question as equally weighted and confirm the exam duration on the timetable in Canvas.
