NivaarExam PrepOfficial exam papers ↗

22-Mec-B9 Advanced Engineering Structures · December 2019

Question 5 of 7: Fatigue-crack growth in an edge-cracked plate

Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)

Notes on this paper

Paper format. 16-Mec-B9 Advanced Engineering Structures, National Examinations, December 2019. Three hours, open book, any non-communicating calculator permitted. Seven problems of equal value; any five constitute a complete paper. All seven are solved here, because the set is a study resource rather than an examination script.

Reference texts. T. H. G. Megson, Aircraft Structures for Engineering Students, 6th ed. (thin-walled open and closed sections, structural idealisation, shear centre, multi-cell torsion); A. P. Boresi and R. J. Schmidt, Advanced Mechanics of Materials, 6th ed. (curved beams, yield criteria); A. C. Ugural and S. K. Fenster, Advanced Mechanics of Materials and Applied Elasticity, 5th ed.; R. C. Hibbeler, Mechanics of Materials, 10th ed. (statically indeterminate axial members, thermal stress); G. E. Dieter, Mechanical Metallurgy, 3rd ed. (fatigue-crack growth, Paris law).

Question 5: Fatigue-crack growth in an edge-cracked plate (20 marks)

Question text not reproduced: the examination questions are © Engineers and Geoscientists BC. Open the official past paper (linked at the top of this page) to read the question, then follow the worked solution below.

Given. A single edge crack in a plate wide enough to be treated as semi-infinite, under constant-amplitude repeated loading, with a Paris growth law of exponent four.

Fatigue data
QuantityValue
Initial crack length a00.4 mm
Repeated stress range Δσ150 N/mm2
Fracture toughness KIC (as printed)180 N/mm3
Paris coefficient C, exponent m28.5 × 10−15, 4
Edge-crack geometry factor Y1.12

Find. (a) the crack length at which fast fracture occurs, and (b) the number of constant-amplitude cycles taken to grow the crack from 0.4 mm to that length.

aΔσ = 150 N/mm²semi-infinite edge crackaₜ = 36.54 mma₀cycles N (failure at 11 035)crack length agrowth is slow until the last few hundred cycles
Figure 5.1 — Edge-cracked plate and the resulting crack-length history. With an exponent of four the crack spends most of its life almost stationary and then runs away over the last few hundred cycles.

Check: the printed unit for fracture toughness is dimensionally impossible, and at the printed magnitude the question has no answer. Stress intensity has units of stress × √length, so N/mm3 cannot be a toughness; the only self-consistent unit here is N/mm3/2, which is also what the Paris coefficient of 28.5 × 10−15 implies. Taken at face value, KIC = 180 N/mm3/2 gives a critical crack length of 0.365 mm, which is shorter than the 0.4 mm crack the plate starts with: the stress intensity at the initial crack is already 188.3 N/mm3/2, so the plate would fracture on the first application of load and both parts would be vacuous. This solution therefore adopts KIC = 1800 N/mm3/2, and tabulates the two other candidate readings in step 5. The reassuring outcome is that the cycle count barely moves between the physically admissible readings — only the critical crack length does.

Approach. Set the stress-intensity range at the crack tip equal to the toughness to obtain the critical crack length, then integrate the Paris law in closed form between the initial and critical lengths.

  1. Write the stress-intensity range. For a single edge crack in a semi-infinite plate the geometry factor is 1.12, so $$\Delta K=Y\,\Delta\sigma\sqrt{\pi a}=1.12(150)\sqrt{\pi a}$$ At the starting crack this is $\Delta K=188.3\ \text{N/mm}^{3/2}$, a number worth computing early because it is what exposes the units problem discussed above.
  2. Part (a) — find the critical crack length. Fast fracture occurs when the peak stress intensity reaches the toughness. The loading is "repeated", meaning it cycles between zero and the maximum, so the peak intensity equals the range and $$a_c=\frac{1}{\pi}\left(\frac{K_{IC}}{Y\,\Delta\sigma}\right)^{2}=\frac{1}{\pi}\left(\frac{1800}{1.12\times150}\right)^{2}$$ $$\boxed{a_c = 36.54\ \text{mm}}$$
  3. Reduce the Paris law to a single power of a. Substituting the stress-intensity expression into $da/dN=C(\Delta K)^m$ with $m=4$ makes the bracket a clean square: $$\frac{da}{dN}=C\,(Y\Delta\sigma)^4\pi^2a^2=C'a^2,\qquad C'=C(Y\Delta\sigma)^4\pi^2=2.2407{\times}10^{-4}$$ The exponent of four is what makes the integral elementary; any other exponent would leave a fractional power of $a$.
  4. Part (b) — integrate to failure. Separating variables between the initial and critical lengths, $$N=\int_{a_0}^{a_c}\frac{da}{C'a^{2}}=\frac{1}{C'}\left(\frac{1}{a_0}-\frac{1}{a_c}\right)=\frac{1}{2.2407{\times}10^{-4}}\left(\frac{1}{0.4}-\frac{1}{36.54}\right)$$ $$\boxed{N = 11\,035\ \text{cycles}}$$ The bracket is $2.500-0.027$, so the critical length contributes barely one per cent of the answer: the life is set almost entirely by how small the starting crack is.
  5. Show how little the reading of KIC matters. Because the second term in the bracket is so small, the life is nearly independent of the toughness. Taking the alternative admissible reading of 180 MPa·√m — that is 5692 N/mm3/2 — moves the critical length by a factor of ten but changes the life by only one per cent, from 11 035 to 11 145 cycles. Only the literal 180 N/mm3/2 is inadmissible, for the reason set out above.
  6. Sense-check the growth rates. At the initial crack the rate is $3.59{\times}10^{-5}$ mm/cycle; at the critical length it is 0.299 mm/cycle, eight thousand times faster. Practically, the crack takes about ten thousand cycles to reach a few millimetres and then covers the remaining thirty millimetres in a few hundred. This is why inspection intervals for fatigue-critical structure are set from the detectable crack size rather than from the critical size.
Question 5 — results
QuantityAdopted KIC = 1800 N/mm3/2If KIC = 180 MPa·√mLiteral 180 N/mm3/2
(a) Crack length at failure36.54 mm365.4 mm0.365 mm (below a0)
(b) Cycles to failure11 03511 1450 — fractures on loading
Lumped Paris coefficient C′2.2407 × 10−4 mm−1/cycle
ΔK at a0 = 0.4 mm188.3 N/mm3/2
Growth rate at a0 / at ac3.59 × 10−5 / 0.299 mm per cycle