NivaarExam PrepOfficial exam papers ↗

22-Mec-B9 Advanced Engineering Structures · May 2016

Question 5 of 8: Paris-law inspection interval for an edge crack

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

Notes on this paper

Paper format. National Exams, May 2016 — 07-Mec-B9 Advanced Engineering Structures. Three hours, open book, any non-communicating calculator permitted. Eight questions of equal total value (20 marks each); the rubric states that any five constitute a complete paper. All eight are solved here, because the complete set is the study resource.

Reference texts.

Axes and sign convention used throughout. A right-handed set is used for every thin-walled question: the span axis runs along the member, $Y$ is vertical (upward) in the cross-section and $Z$ is horizontal (to the right). Section constants are written $I_{ZZ}=\int Z^{2}\,dA$, $I_{YY}=\int Y^{2}\,dA$ and $I_{YZ}=\int ZY\,dA$ about the centroid, and the bending stress is carried as the linear field $\sigma = aZ + bY$, whose coefficients follow from the pair

$$\begin{aligned} a\,I_{ZZ} + b\,I_{YZ} &= M_{Y}, \\ a\,I_{YZ} + b\,I_{YY} &= M_{Z} \end{aligned}$$

with $M_{Y}\equiv\int\sigma Z\,dA$ and $M_{Z}\equiv\int\sigma Y\,dA$. This one pair replaces the Megson fraction and handles symmetric and unsymmetrical sections alike. Shear flows are counted positive in the direction of the walk stated with each figure; a negative answer simply means the flow runs the other way round the circuit. Moments and torques are anticlockwise-positive in the $Z$–$Y$ plane.

Question 5: Paris-law inspection interval for an edge crack (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 semi-infinite plate carrying a single edge crack under constant-amplitude cyclic tension, with a Paris exponent of exactly four.

Given data — Question 5
QuantitySymbolValue
Initial edge crack length$a_{0}$0.25 mm
Constant-amplitude stress range$\Delta\sigma$130 N/mm2
Fracture toughness$K_{IC}$2 950 N/mm3/2
Paris coefficient$C$$34\times10^{-15}$ mm/cycle
Paris exponent$m$4
Edge-crack geometry factor$Y$1.12

Find. The number of cycles taken for the crack to grow from 0.25 mm to half its critical length, which is the longest permissible interval between inspections.

Approach. Get the critical crack length from the toughness, halve it to set the detection target, then integrate the Paris law between the two lengths — which for $m = 4$ has a closed form.

  1. Write the stress intensity factor for an edge crack. For a single edge crack in a semi-infinite plate the standard result is $$K = Y\sigma\sqrt{\pi a} = 1.12\,\sigma\sqrt{\pi a}$$ The 1.12 factor is the free-surface correction; a centre crack of the same length would give $Y = 1$ and a longer life.
  2. Find the critical crack length. Fast fracture occurs when the peak stress intensity reaches the toughness, $1.12\Delta\sigma\sqrt{\pi a_{c}} = K_{IC}$, so $$a_{c} = \frac{1}{\pi}\left(\frac{K_{IC}}{1.12\Delta\sigma}\right)^{2} = \frac{1}{\pi}\left(\frac{2950}{1.12\times130}\right)^{2} = \frac{1}{\pi}\left(20.261\right)^{2} = \boxed{130.67\ \text{mm}}$$ so the detection target is $$a_{f} = \frac{a_{c}}{2} = \boxed{65.33\ \text{mm}}$$
  3. Set up the Paris integral. Substituting the stress intensity range into $da/dN = C\left(\Delta K\right)^{4}$ gives $$\begin{aligned} \frac{da}{dN} &= C\left(1.12\Delta\sigma\right)^{4}\pi^{2}a^{2} = C'a^{2}, \\ C' &= C\left(1.12\Delta\sigma\right)^{4}\pi^{2} \end{aligned}$$ Because the exponent is exactly four, the crack length appears only as $a^{2}$ and the integral is elementary.
  4. Evaluate the lumped constant. $$C' = 34\times10^{-15}\times\left(145.6\right)^{4}\times\pi^{2} = 34\times10^{-15}\times4.4941\times10^{8}\times9.8696 = \boxed{1.5081\times10^{-4}}$$ in units of reciprocal millimetre-cycles.
  5. Integrate between the two crack lengths. Separating variables, $$N = \int_{a_{0}}^{a_{f}}\frac{da}{C'a^{2}} = \frac{1}{C'}\left(\frac{1}{a_{0}}-\frac{1}{a_{f}}\right)$$ $$N = \frac{1}{1.5081\times10^{-4}}\left(\frac{1}{0.25}-\frac{1}{65.33}\right) = 6631.2\times\left(4.0000-0.0153\right) = \boxed{26\,422\ \text{cycles}}$$
  6. Sanity-check the number against the full life. Repeating the integral to the full critical length gives 26 473 cycles, so the last half of the crack’s growth — from 65 mm to 131 mm — takes only 51 cycles, two tenths of one per cent of the life. The stress intensity range rises from 129 N/mm3/2 at $a_{0}$ to 2086 N/mm3/2 at $a_{f}$, and since the growth rate goes as the fourth power of that, the final stage is effectively instantaneous.
  7. State the maintenance interval. The inspection interval must be no longer than $$\boxed{N_{\text{inspect}} = 26\,400\ \text{cycles}}$$ rounded down from 26 422, and in practice a factor of two or more would be applied to allow for scatter in $C$, for the detection threshold of the chosen inspection method, and for the fact that a crack shorter than 0.25 mm may already exist at the start of service.
0 27 55 82 110 137 0k 7k 13k 20k 26k a critical = 130.7 mm half critical = 65.3 mm N = 26 422 cycles N crack length a (mm) the axis stops at the maintenance interval; growth is flat then runs away
Figure 5.1 — crack length against cycles from the Paris law. The curve is almost flat for the first twenty thousand cycles and then turns sharply upward, which is why the inspection interval must be set well before the target length is reached.
Question 5 — results
QuantityResult
Critical crack length130.67 mm
Detection target (half critical)65.33 mm
Lumped Paris constant $C'$$1.5081\times10^{-4}$ mm−1 per cycle
$\Delta K$ at the initial crack129.0 N/mm3/2
$\Delta K$ at the detection target2 086 N/mm3/2
Cycles from 0.25 mm to 65.33 mm26 422 cycles
Cycles to full critical length26 473 cycles
Recommended maintenance interval26 400 cycles (before any additional scatter factor)