Physics:Paris' law

Short description: Formula in materials science

Typical plot of crack growth rate with respect to the stress intensity range where the Paris–Erdogan equation fits the central, linear region of Regime B.

Paris' law (also known as the Paris–Erdogan equation) is a crack growth equation that gives the rate of growth of a fatigue crack. The stress intensity factor [math]\displaystyle{ K }[/math] characterises the load around a crack tip and the rate of crack growth is experimentally shown to be a function of the range of stress intensity [math]\displaystyle{ \Delta K }[/math] seen in a loading cycle. The Paris equation is

[math]\displaystyle{ \begin{align} {da \over dN} &= C\left(\Delta K\right)^{m}, \end{align} }[/math]

where [math]\displaystyle{ a }[/math] is the crack length and [math]\displaystyle{ {\rm d}a/{\rm d}N }[/math] is the fatigue crack growth for a load cycle [math]\displaystyle{ N }[/math]. The material coefficients [math]\displaystyle{ C }[/math] and [math]\displaystyle{ m }[/math] are obtained experimentally and also depend on environment, frequency, temperature and stress ratio. The stress intensity factor range has been found to correlate the rate of crack growth from a variety of different conditions and is the difference between the maximum and minimum stress intensity factors in a load cycle and is defined as

[math]\displaystyle{ \Delta K= K_{\text{max}}-K_{\text{min}}. }[/math]

Being a power law relationship between the crack growth rate during cyclic loading and the range of the stress intensity factor, the Paris–Erdogan equation can be visualized as a straight line on a log-log plot, where the x-axis is denoted by the range of the stress intensity factor and the y-axis is denoted by the crack growth rate.

The ability of ΔK to correlate crack growth rate data depends to a large extent on the fact that alternating stresses causing crack growth are small compared to the yield strength. Therefore crack tip plastic zones are small compared to crack length even in very ductile materials like stainless steels.

The equation gives the growth for a single cycle. Single cycles can be readily counted for constant-amplitude loading. Additional cycle identification techniques such as rainflow-counting algorithm need to be used to extract the equivalent constant-amplitude cycles from a variable-amplitude loading sequence.

History

In a 1961 paper, P. C. Paris introduced the idea that the rate of crack growth may depend on the stress intensity factor. Then in their 1963 paper, Paris and Erdogan indirectly suggested the equation with the aside remark "The authors are hesitant but cannot resist the temptation to draw the straight line slope 1/4 through the data" after reviewing data on a log-log plot of crack growth versus stress intensity range. The Paris equation was then presented with the fixed exponent of 4.

Domain of applicability

Stress ratio

Higher mean stress is known to increase the rate of crack growth and is known as the mean stress effect. The mean stress of a cycle is expressed in terms of the stress ratio [math]\displaystyle{ R }[/math] which is defined as

[math]\displaystyle{ R = {K_{\text{min}} \over K_{\text{max}}}, }[/math]

or ratio of minimum to maximum stress intensity factors. In the linear elastic fracture regime, [math]\displaystyle{ R }[/math] is also equivalent to the load ratio

[math]\displaystyle{ R \equiv {P_{\text{min}} \over P_{\text{max}}}. }[/math]

The Paris–Erdogan equation does not explicitly include the effect of stress ratio, although equation coefficients can be chosen for a specific stress ratio. Other crack growth equations such as the Forman equation do explicitly include the effect of stress ratio, as does the Elber equation by modelling the effect of crack closure.

Intermediate stress intensity range

The Paris–Erdogan equation holds over the mid-range of growth rate regime, but does not apply for very low values of [math]\displaystyle{ \Delta K }[/math]approaching the threshold value [math]\displaystyle{ \Delta K_{\text{th}} }[/math], or for very high values approaching the material's fracture toughness, [math]\displaystyle{ K_{\text{Ic}} }[/math]. The alternating stress intensity at the critical limit is given by [math]\displaystyle{ \begin{align} \Delta K_{\text{cr}} &= (1-R)K_{\text{Ic}} \end{align} }[/math].

The slope of the crack growth rate curve on log-log scale denotes the value of the exponent [math]\displaystyle{ m }[/math] and is typically found to lie between [math]\displaystyle{ 2 }[/math] and [math]\displaystyle{ 4 }[/math], although for materials with low static fracture toughness such as high-strength steels, the value of [math]\displaystyle{ m }[/math] can be as high as [math]\displaystyle{ 10 }[/math].

Long cracks

Because the size of the plastic zone [math]\displaystyle{ (r_{\text{p}} \approx K_{I}^2/\sigma_{y}^2) }[/math] is small in comparison to the crack length, [math]\displaystyle{ a }[/math] (here, [math]\displaystyle{ \sigma_{y} }[/math] is yield stress), the approximation of small-scale yielding applies, enabling the use of linear elastic fracture mechanics and the stress intensity factor. Thus, the Paris–Erdogan equation is only valid in the linear elastic fracture regime, under tensile loading and for long cracks.