What is π doing in the Stress Intensity expression?

Today it is π-day (3/14), a perfect day to wonder what π is doing in the stress intensity factor (SIF) equation. Let us have a look to this SIF equation first:

K=\beta \cdot S\cdot \sqrt{\pi \cdot a}

In this equation, K is the stress intensity factor, β is a geometry factor, S is the gross stress and a is the (half) crack length. This SIF describes the severity of the stress field around a crack tip. With increasing SIF, the fatigue crack growth rate will increase as well:

\frac{da}{dn}=C\cdot \left (\Delta K \right )^{n}

It seems that the SIF equation could do quite well without the π. To understand why the π was introduced we have to go back in history.

Westergaard

Westergaard [1] introduced an Airy stress function of complex numbers to describe the stress field in an infinite plate with a crack. The equations are:

\sigma _{xx}=Re\left (Z \right )-y\cdot Im\left ( Z' \right )

\sigma _{yy}=Re\left (Z \right )+y\cdot Im\left ( Z' \right )

\tau _{xy}=-y\cdot Re\left (Z' \right )

Z and Z’ are:

Z\left ( z \right )=\frac{S}{\sqrt{1-\left ( \frac{a}{z} \right )^{2}}}

Z'\left ( z \right )=\frac{-S\cdot a^{2}}{z^{3} \left [(1-\left ( \frac{a}{z} \right )^{2} \right ]^{\frac{3}{2}}}

a is the half crack length and z equals x+iy. For y=0, so in the plane of the crack:

\sigma _{yy}=\frac{S}{\sqrt{1-\left ( \frac{a}{x} \right )^{2}}}

At a=x, there is a singularity; σyy is infinite.

Irwin

Later Irwin [2] simplified Westergaard’s solution. He proposed the following expression for z:

z=a+r\cdot e^{^{i\theta }}

By doing so, Irwin found expressions for the stress field in the close vicinity of the crack tip using the polar coordinates r and θ. After applying a lot of algebra, the following equation for y=0 is found:

\sigma _{yy}=\frac{S\sqrt{a}}{\sqrt{2r}}\cdot cos\frac{\theta }{2}\left (1+sin\frac{\theta }{2}sin\frac{3\theta }{2} \right )

Still without a π. In the years after Irwin published his paper, the term √(π/π) was introduced in the equation:

\sigma _{yy}=\frac{S\sqrt{\pi a}}{\sqrt{2\pi r}}\cdot cos\frac{\theta }{2}\left (1+sin\frac{\theta }{2}sin\frac{3\theta }{2} \right )

With \inline K=S\cdot \sqrt{\pi \cdot a}  the stress field can be expressed as:

\sigma _{xx}=\frac{K}{\sqrt{2\pi r}}\cdot f\left ( \theta \right )\;\;\;\;\; \sigma _{yy}=\frac{K}{\sqrt{2\pi r}}\cdot g\left ( \theta \right )\;\;\;\;\; \tau _{xy}=\frac{K}{\sqrt{2\pi r}}\cdot h\left ( \theta \right )

This is only an approximation of the exact solution of Westergaard. Close to the crack tip the approximation is quite good. The advantage of Irwin’s approximation is that it uses the stress intensity factor K. This SIF is the key parameter in linear elastic fracture mechanics.

But we still haven’t seen why we need π in the equation for K.

Griffith

Early last century, Griffith [3] worked on an energy-based criterion to predict the failure behaviour of a plate with a crack. in short, the criterion implies that failure (unstable crack growth) occurs when the release of strain energy is larger than the required energy to break the atomic bonds in a material. Griffith critical energy release rate is written as:

G_{c}=\frac{S_{c}^{2}\pi a}{E}

Decades later Irwin introduced the SIF, leading to:

G=\frac{K^{2}}{E}

And here we have the point where the π was introduced.

References

  1. Westergaard, H.M., Bearing Pressures and Cracks, Journal of Applied Mechanics, Vol. 6, pp. A49-53, 1939.
  2. Irwin, G.R., Analysis of Stresses and Strains Near the End of a Crack Traversing a Plate, Journal of Applied Mechanics, Vol. 24, pp. 361-364, 1957.
  3. Griffith, A.A., The Phenomena of Rupture and Flow in Solids, Philosophical Transactions, Series A, Vol. 221, pp. 163-198, 1920.

Do you want to know more about fatigue?

Subscribe to our newsletter and we will keep you informed about fatigue and Fatec Engineering.