Crack growth equation

A crack growth equation is used for calculating the size of a fatigue crack growing from cyclic loads. The growth of fatigue cracks can result in catastrophic failure, particularly in the case of aircraft. A crack growth equation can be used to ensure safety, both in the design phase and during operation, by predicting the size of cracks. In critical structure, loads can be recorded and used to predict the size of cracks to ensure maintenance or retirement occurs prior to any of the cracks failing.

Fatigue life can be subdivided into an initiation or nucleation period and a crack growth period.[1] The initiation period consists of crack initiation and microcrack growth which leads to the next phase of macrocrack growth. Crack growth equations are used to predict the crack size starting from a given initial flaw and are typically based on experimental data obtained from constant amplitude fatigue tests.

One of the earliest crack growth equations based on the stress intensity factor range of a load cycle (${\displaystyle \Delta K}$) is the Paris-Erdogan equation[2]

${\displaystyle {da \over dN}=C(\Delta K)^{m}}$

where ${\displaystyle a}$ is the crack length and ${\displaystyle {\rm {d}}a/{\rm {d}}N}$ is the fatigue crack growth for a single load cycle ${\displaystyle N}$. A variety of crack growth equations similar to the Paris-Erdogan equation have been developed to include factors that affect the crack growth rate such as stress ratio, overloads and load history effects.

The stress intensity range can be calculated from the maximum and minimum stress intensity for a cycle

${\displaystyle \Delta K=K_{\text{max}}-K_{\text{min}}}$

A geometry factor ${\displaystyle \beta }$ is used to relate the far field stress ${\displaystyle \sigma }$ to the crack tip stress intensity using

${\displaystyle K=\beta \sigma {\sqrt {\pi a}}}$.

There are standard references containing the geometry factors for many different configurations.[3][4][5]

History of crack propagation laws

Many crack propagation laws have been proposed over the years to improve prediction accuracy and incorporate a variety of effects. The works of Head,[6] Frost and Dugdale,[7] McEvily and Illg,[8] and Liu[9] on fatigue crack-growth behaviour laid the foundation in this topic. The general form of these crack propagation equations may be expressed as

${\displaystyle {da \over dN}=f(\Delta \sigma ,a,C_{i}),}$

where, the crack length is denoted by ${\displaystyle a}$, the number of cycles of load applied is given by ${\displaystyle N}$, the stress range by ${\displaystyle \Delta \sigma }$, and the material parameters by ${\displaystyle C_{i}}$. For symmetrical configurations, the length of the crack from the line of symmetry is defined as ${\displaystyle a}$ and is half of the total crack length ${\displaystyle 2a}$.

Crack growth equations of the form ${\displaystyle da/dN}$ are not a true differential equation as they do not model the process of crack growth in a continuous manner throughout the loading cycle. As such, separate cycle counting or identification algorithms such as the commonly used rainflow-counting algorithm, are required to identify the maximum and minimum values in a cycle. Although developed for the stress/strain-life methods rainflow counting has also been shown to work for crack growth.[10] There have been a small number of true derivative fatigue crack growth equations that have also been developed.[11][12]

Crack growth rate in different regimes

Figure 1 shows a typical plot of the rate of crack growth as a function of the alternating stress intensity or crack tip driving force ${\displaystyle \Delta K}$ plotted on log scales. The crack growth rate behaviour with respect to the alternating stress intensity can be explained in different regimes (see, figure 1) as follows

Regime A: At low growth rates, variations in microstructure, mean stress (or load ratio), and environment have significant effects on the crack propagation rates. It is observed at low load ratios that the growth rate is most sensitive to microstructure and in low strength materials it is most sensitive to load ratio.[13]

To predict the crack growth rate at the near threshold region, the following relation[14] has been used

${\displaystyle {da \over dN}=A\left(\Delta K-\Delta K_{\text{th}}\right)^{p}.}$

Regime B: At mid-range of growth rates, variations in microstructure, mean stress (or load ratio), thickness, and environment have no significant effects on the crack propagation rates.

Paris-Erdogan Equation

To predict the crack growth rate in this intermediate regime, the Paris' law is used[2]

${\displaystyle {da \over dN}=C\left(\Delta K\right)^{m}.}$

Regime C: At high growth rates, crack propagation is highly sensitive to the variations in microstructure, mean stress (or load ratio), and thickness. Environmental effects have relatively very less influence.

Forman Equation

In 1967, Forman proposed the following relation to account for the increased growth rates due to stress ratio and when approaching the fracture toughness ${\displaystyle K_{\text{c}}}$[15]

${\displaystyle {da \over dN}={\frac {C(\Delta K)^{n}}{(1-R)K_{\text{c}}-\Delta K}}}$

McEvily-Groeger Equation

McEvily and Groeger[16] proposed the following power-law relationship which considers the effects of both high and low values of ${\displaystyle \Delta K}$

${\displaystyle {da \over dN}=A(\Delta K-\Delta K_{\text{th}})^{2}{\Big [}1+{\frac {\Delta K}{K_{\text{IC}}-K_{\text{max}}}}{\Big ]}}$.

NASGRO Equation

The NASGRO equation is used in the crack growth programs AFGROW, FASTRAN and NASGRO software.[17] It is a general equations that covers the lower growth rates near the threshold and increased growth rates near the fracture toughness, along with the mean stress effect

${\displaystyle {\frac {da}{dN}}=C\left[\left({\frac {1-f}{1-R}}\right)\Delta K\right]^{n}{\left(1-{\frac {\Delta K_{\text{th}}}{\Delta K}}\right)^{p} \over \left(1-{\frac {K_{\max }}{K_{\text{crit}}}}\right)^{q}}}$

where ${\displaystyle C}$, ${\displaystyle f}$, ${\displaystyle n}$, ${\displaystyle p}$, ${\displaystyle q}$, ${\displaystyle \Delta K_{\text{th}}}$ and ${\displaystyle K_{\text{crit}}}$ are the equation coefficients.

McClintock Equation

In 1967, McClintock developed an equation for the upper limit of crack growth based on the cyclic crack tip opening displacement ${\displaystyle \Delta {\text{CTOD}}}$[18]

${\displaystyle {da \over dN}\propto \Delta {\text{CTOD}}\approx \beta {(\Delta K)^{2} \over {2\sigma _{0}E'}}}$

where ${\displaystyle \sigma _{0}}$ is the flow stress, ${\displaystyle E'}$ is the Young's modulus and ${\displaystyle \beta }$ is a constant typically in the range 0.1-0.5.

Stress ratio effect

Cycles with higher stress ratio ${\displaystyle R=K_{\text{min}}}/{K_{\text{max}}\equiv P_{\text{min}}/{P_{\text{max}}}}$ have an increased rate of crack growth.[19] This effect is often explained using the crack closure concept which describes the observation that the crack faces can remain in contact with each other at loads above zero. This reduces the effective stress intensity factor range and the fatigue crack growth rate.[20]

Walker Equation

To account for the stress ratio effect, Walker[21] suggested a modified form of the Paris' law

${\displaystyle {da \over dN}=C_{0}{\Big (}{\overline {\Delta K}}{\Big )}^{m}=C_{0}{\bigg (}{\frac {(\Delta K)}{(1-R)^{1-\gamma }}}{\bigg )}^{m}=C_{0}{\big (}K_{\text{max}}(1-R)^{\gamma }{\big )}^{m}}$

where, ${\displaystyle \gamma }$ is a material parameter which represents the influence of stress ratio on the fatigue crack growth rate. Typically, ${\displaystyle \gamma }$ takes a value around ${\displaystyle 0.5}$, but can vary between ${\displaystyle 0.3-1.0}$. In general, it is assumed that compressive portion of the loading cycle ${\displaystyle {\big (}R<0{\big )}}$ has no effect on the crack growth by considering ${\displaystyle \gamma =0,}$ which gives ${\displaystyle {\overline {\Delta K}}=K_{\text{max}}.}$ This can be physically explained by considering that the crack closes at zero load and does not behave like a crack under compressive loads. In very ductile materials like Man-Ten steel[22] compressive loading does contribute to the crack growth according to ${\displaystyle \gamma =0.22}$.

Elber Equation

Elber modified the Paris-Erdogan equation to allow for crack closure with the introduction of the opening stress intensity level ${\displaystyle K_{\text{op}}}$ at which contact occurs. Below this level there is no movement at the crack tip and hence no growth. This effect has been used to explain the stress ratio effect and the increased rate of growth observed with short cracks. Elber's equation is [23]

${\displaystyle \Delta K_{\text{eff}}=K_{\text{max}}-K_{\text{op}}}$
${\displaystyle {da \over dN}=C(\Delta K_{\text{eff}})^{m}}$

Fatigue crack propagation in ductile and brittle materials

The general form of the fatigue-crack growth rate in ductile and brittle materials is given by[18]

${\displaystyle {da \over dN}\propto (K_{\text{max}})^{n}(\Delta K)^{p},}$

where, ${\displaystyle n}$ and ${\displaystyle p}$ are material parameters. Based on different crack-advance and crack-tip shielding mechanisms in metals, ceramics, and intermetallics, it is observed that the fatigue crack growth rate in metals is significantly dependent on ${\displaystyle \Delta K}$ term, in ceramics on ${\displaystyle K_{\text{max}}}$, and intermetallics have almost similar dependence on ${\displaystyle \Delta K}$ and ${\displaystyle K_{\text{max}}}$ terms.

Prediction of fatigue life

Computer programs

There are many computer programs that implement crack growth equations such as Nasgro,[24] AFGROW and Fastran. In addition, there are also programs that implement a probabilistic approach to crack growth that calculate the probability of failure throughout the life of a component.[25][26]

Crack growth programs grow a crack from an initial flaw size until it exceeds the fracture toughness of a material and fails. Because the fracture toughness depends on the boundary conditions, the fracture toughness may change from plane strain conditions for a semi-circular surface crack to plane stress conditions for a through crack. The fracture toughness for plane stress conditions is typically twice as large as that for plane strain. However, because of the rapid rate of growth of a crack near the end of its life, variations in fracture toughness do not significantly alter the life of a component.

Crack growth programs typically provide a choice of:

• cycle counting methods to extract cycle extremes
• geometry factors that select for the shape of the crack and the applied loading
• crack growth equation
• acceleration/retardation models
• material properties such as yield strength and fracture toughness

Analytical solution

The stress intensity factor is given by

${\displaystyle K=\sigma {\sqrt {\pi a}}Y{\Big (}{\frac {a}{W}}{\Big )},}$

where ${\displaystyle \sigma }$ is the applied uniform tensile stress acting on the specimen in the direction perpendicular to the crack plane, ${\displaystyle a}$ and ${\displaystyle W}$represent the crack size and width of the specimen respectively, and ${\displaystyle Y{\Big (}{\frac {a}{W}}{\Big )}}$ is a dimensionless parameter that depends on the geometry of the specimen. The alternating stress intensity becomes

{\displaystyle {\begin{aligned}\Delta K&={\begin{cases}(\sigma _{\text{max}}-\sigma _{\text{min}})Y{\sqrt {\pi a}}=\Delta \sigma Y{\sqrt {\pi a}},\qquad R\geq 0\\\sigma _{\text{max}}Y{\sqrt {\pi a}},\qquad R<0\end{cases}},\end{aligned}}}

where ${\displaystyle \Delta \sigma }$ is the range of the cyclic stress amplitude.

By assuming the initial crack size to be ${\displaystyle a_{0}}$, the critical crack size ${\displaystyle a_{c}}$ before the specimen fails can be computed using ${\displaystyle {\big (}K=K_{\text{max}}=K_{\text{IC}}{\big )}}$ as

{\displaystyle {\begin{aligned}K_{\text{IC}}&=\sigma _{\text{max}}Y{\Big (}{\frac {a_{c}}{W}}{\Big )}{\sqrt {\pi a_{c}}},\\\Rightarrow a_{c}&={\frac {1}{\pi }}{\bigg (}{\frac {K_{\text{IC}}}{\sigma _{\text{max}}Y{\big (}{\frac {a_{c}}{W}}{\big )}}}{\bigg )}^{2}.\end{aligned}}}

The above equation in ${\displaystyle a_{c}}$is implicit in nature and can be solved numerically if necessary.

Case I

For ${\displaystyle R\geq 0.7,}$ crack closure has negligible effect on the crack growth rate[27] and the Paris' law can be used to compute the fatigue life of a specimen before it reaches the critical crack size ${\displaystyle a_{c}}$ as

{\displaystyle {\begin{aligned}{da \over dN}&=C(\Delta K)^{m}=C{\bigg (}\Delta \sigma Y{\Big (}{\frac {a}{W}}{\Big )}{\sqrt {\pi a}}{\bigg )}^{m},\\\Rightarrow N_{f}&={\frac {1}{({\sqrt {\pi }}\Delta \sigma )^{m}}}\int _{a_{0}}^{a_{c}}{\frac {da}{(C{\sqrt {a}}Y{\Big (}{\frac {a}{W}}{\Big )})^{m}}}.\end{aligned}}}
Crack growth model with constant value of ${\displaystyle Y}$and R = 0

For the Griffith-Irwin crack growth model or center crack of length ${\displaystyle 2a}$ in an infinite sheet as shown in the figure 2, we have ${\displaystyle Y=1}$ and is independent of the crack length. Also, ${\displaystyle C}$ can be considered to be independent of the crack length. By assuming ${\displaystyle Y={\text{constant}},}$ the above integral simplifies to

${\displaystyle N_{f}={\frac {1}{C({\sqrt {\pi }}Y\Delta \sigma )^{m}}}\int _{a_{0}}^{a_{c}}{\frac {da}{({\sqrt {a}})^{m}}},}$

by integrating the above expression for ${\displaystyle m\neq 2}$ and ${\displaystyle m=2}$ cases, the total number of load cycles ${\displaystyle N_{f}}$ are given by

{\displaystyle {\begin{aligned}N_{f}&={\frac {2}{(m-2)C({\sqrt {\pi }}\Delta \sigma Y)^{m}}}{\Bigg [}{\frac {1}{(a_{0})^{\frac {m-2}{2}}}}-{\frac {1}{(a_{c})^{\frac {m-2}{2}}}}{\Bigg ]},\qquad m\neq 2,\\N_{f}&={\frac {1}{\pi C(\Delta \sigma Y)^{2}}}\ln {\frac {a_{c}}{a_{0}}},\qquad m=2.\end{aligned}}}

Now, for ${\displaystyle m>2}$ and critical crack size to be very large in comparison to the initial crack size ${\displaystyle {\big (}a_{c}>>a_{0}{\big )}}$ will give

${\displaystyle N_{f}={\frac {2}{(m-2)C({\sqrt {\pi }}\Delta \sigma Y)^{m}}}(a_{0})^{\frac {2-m}{2}}.}$

The above analytical expressions for the total number of load cycles to fracture ${\displaystyle {\big (}N_{f}{\big )}}$ are obtained by assuming ${\displaystyle Y={\text{constant}}}$. For the cases, where ${\displaystyle Y}$ is dependent on the crack size such as the Single Edge Notch Tension (SENT), Center Cracked Tension (CCT) geometries, numerical integration can be used to compute ${\displaystyle N_{f}}$.

Case II

For ${\displaystyle R<0.7,}$ crack closure phenomenon has an effect on the crack growth rate and we can invoke Walker equation to compute the fatigue life of a specimen before it reaches the critical crack size ${\displaystyle a_{c}}$ as

{\displaystyle {\begin{aligned}{da \over dN}&=C{\bigg (}{\frac {\Delta K}{(1-R)^{1-\gamma }}}{\bigg )}^{m}={\frac {C}{(1-R)^{m(1-\gamma )}}}{\bigg (}\Delta \sigma Y{\Big (}{\frac {a}{W}}{\Big )}{\sqrt {\pi a}}{\bigg )}^{m},\\\Rightarrow N_{f}&={\frac {(1-R)^{m(1-\gamma )}}{({\sqrt {\pi }}\Delta \sigma )^{m}}}\int _{a_{0}}^{a_{c}}{\frac {da}{(C{\sqrt {a}}Y{\Big (}{\frac {a}{W}}{\Big )})^{m}}}.\end{aligned}}}

Numerical calculation

This scheme is useful when ${\displaystyle Y}$ is dependent on the crack size ${\displaystyle a}$. The initial crack size is considered to be ${\displaystyle a_{0}}$. The stress intensity factor at the current crack size ${\displaystyle a}$ is computed using the maximum applied stress as

{\displaystyle {\begin{aligned}K_{\text{max}}&=Y\sigma _{\text{max}}{\sqrt {\pi a}}.\end{aligned}}}
If ${\displaystyle K_{\text{max}}}$ is less than the fracture toughness ${\displaystyle K_{\text{IC}}}$, the crack has not reached its critical size ${\displaystyle a_{c}}$ and the simulation is continued with the current crack size to calculate the alternating stress intensity as

${\displaystyle \Delta K=Y{\Big (}{\frac {a}{W}}{\Big )}\Delta \sigma {\sqrt {\pi a}}.}$

Now, by substituting the stress intensity factor in Paris' law, the increment in the crack size ${\displaystyle \Delta a}$ is computed as

${\displaystyle \Delta a=C(\Delta K)^{m}\Delta N,}$

where ${\displaystyle \Delta N}$ is cycle step size. The new crack size becomes

${\displaystyle a_{i+1}=a_{i}+\Delta a,}$

where index ${\displaystyle i}$ refers to the current iteration step. The new crack size is used to calculate the stress intensity at maximum applied stress for the next iteration. This iterative process is continued until

${\displaystyle K_{\text{max}}\geq K_{\text{IC}}.}$

Once this failure criterion is met, the simulation is stopped.

The schematic representation of the fatigue life prediction process is shown in figure 3.

Example

The stress intensity factor in a SENT specimen (see, figure 4) under fatigue crack growth is given by[5]

{\displaystyle {\begin{aligned}K_{I}&=\sigma {\sqrt {\pi a}}\,Y{\Big (}{\frac {a}{W}}{\Big )}=\sigma {\sqrt {\pi a}}{\Bigg [}0.265{\bigg [}1-{\frac {a}{W}}{\bigg ]}^{4}+{\frac {0.857+0.265{\frac {a}{W}}}{{\big [}1-{\frac {a}{W}}{\big ]}^{\frac {3}{2}}}}{\Bigg ]},\\\Delta K_{I}&=K_{\text{max}}-K_{\text{min}}=\Delta \sigma {\sqrt {\pi a}}\,Y{\Big (}{\frac {a}{W}}{\Big )}.\end{aligned}}}

The following parameters are considered for the calculation

{\displaystyle {\begin{aligned}a_{0}&=5{\text{ mm}},\qquad W=100{\text{ mm}},\qquad h=200{\text{ mm}},\qquad K_{\text{IC}}=30{\text{ MPa}}{\sqrt {\text{m}}},\qquad R={\frac {K_{\text{min}}}{K_{\text{max}}}}=0.7,\\\Delta \sigma &=20{\text{ MPa}}\qquad C=4.6774\times 10^{-11}{\frac {\text{m}}{\text{cycle}}}{\frac {1}{({\text{MPa}}{\sqrt {\text{m}}})^{m}}},\qquad m=3.874.\end{aligned}}}

The critical crack length, ${\displaystyle a=a_{c}}$, can be computed when ${\displaystyle K_{\text{max}}=K_{\text{IC}}}$ as

${\displaystyle a_{c}={\frac {1}{\pi }}{\Bigg (}{\frac {0.45}{Y{\Big (}{\frac {a_{c}}{W}}{\Big )}}}{\Bigg )}^{2}.}$

By solving the above equation, the critical crack length is obtained as ${\displaystyle a_{c}=26.7{\text{ mm}}}$.

Now, invoking the Paris' law gives

${\displaystyle N_{f}={\frac {1}{C(\Delta \sigma )^{m}({\sqrt {\pi }})^{m}}}\int _{a_{0}}^{a_{c}}{\frac {da}{a^{\frac {m}{2}}{\Bigg [}0.265{\bigg [}1-{\frac {a}{W}}{\bigg ]}^{4}+{\frac {0.857+0.265{\frac {a}{W}}}{{\big [}1-{\frac {a}{W}}{\big ]}^{\frac {3}{2}}}}{\Bigg ]}^{m}}}}$

By numerical integration of the above expression, the total number of load cycles to failure is obtained as ${\displaystyle N_{f}=1.2085\times 10^{6}{\text{ cycles}}}$.

References

1. Schijve, J. (January 1979). "Four lectures on fatigue crack growth". Engineering Fracture Mechanics. 11 (1): 169–181. doi:10.1016/0013-7944(79)90039-0. ISSN 0013-7944.
2. Paris, P. C.; Erdogan, F. (1963). "A critical analysis of crack propagation laws". Journal of Basic Engineering. 18 (4): 528–534. doi:10.1115/1.3656900..
3. Murakami, Y.; Aoki, S. (1987). Stress Intensity Factors Handbook. Pergamon, Oxford.
4. Rooke, D. P.; Cartwright, D. J. (1976). Compendium of Stress Intensity Factors. Her Majesty’s Stationery Office, London.
5. Tada, Hiroshi; Paris, Paul C.; Irwin, George R. (1 January 2000). The Stress Analysis of Cracks Handbook, Third Edition. Three Park Avenue New York, NY 10016-5990: ASME. doi:10.1115/1.801535. ISBN 0791801535.
6. Head, A.K. (September 1953). "XCVIII. The growth of fatigue cracks". The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science. 44 (356): 925–938. doi:10.1080/14786440908521062. ISSN 1941-5982.
7. Frost, N.E.; Dugdale, D.S. (January 1958). "The propagation of fatigue cracks in sheet specimens". Journal of the Mechanics and Physics of Solids. 6 (2): 92–110. doi:10.1016/0022-5096(58)90018-8. ISSN 0022-5096.
8. McEvily, Arthur J.; Illg, Walter (1960), "A Method for Predicting the Rate of Fatigue-Crack Propagation", Symposium on Fatigue of Aircraft Structures, ASTM International, pp. 112–112–8, doi:10.1520/stp45927s, ISBN 9780803165793
9. Liu, H. W. (1961). "Crack Propagation in Thin Metal Sheet Under Repeated Loading". Journal of Basic Engineering. 83 (1): 23–31. doi:10.1115/1.3658886. ISSN 0021-9223.
10. Sunder, R.; Seetharam, S. A.; Bhaskaran, T. A. (1984). "Cycle counting for fatigue crack growth analysis". International Journal of Fatigue. 6 (3): 147–156. doi:10.1016/0142-1123(84)90032-X.
11. Pommier, S.; Risbet, M. (2005). "Time derivative equations for mode I fatigue crack growth in metals". International Journal of Fatigue. 27 (10–12): 1297–1306. doi:10.1016/j.ijfatigue.2005.06.034.
12. Lu, Zizi; Liu, Yongming (2010). "Small time scale fatigue crack growth analysis". International Journal of Fatigue. 32 (8): 1306–1321. doi:10.1016/j.ijfatigue.2010.01.010.
13. Ritchie, R. O. (1977). "Near-Threshold Fatigue Crack Propagation in Ultra-High Strength Steel: Influence of Load Ratio and Cyclic Strength". Journal of Engineering Materials and Technology. 99 (3): 195–204. doi:10.1115/1.3443519. ISSN 0094-4289.
14. Allen, R. J.; Booth, G. S.; Jutla, T. (March 1988). "A review of fatigue crack growth characterisation by Linear Elastic Fracture Mechanics (LEFM). Part II - Advisory documents and applications within National Standards". Fatigue & Fracture of Engineering Materials and Structures. 11 (2): 71–108. doi:10.1111/j.1460-2695.1988.tb01162.x. ISSN 8756-758X.
15. Forman, R. G.; Kearney, V. E.; Engle, R. M. (1967). "Numerical Analysis of Crack Propagation in Cyclic-Loaded Structures". Journal of Basic Engineering. 89 (3): 459–463. doi:10.1115/1.3609637. ISSN 0021-9223.
16. McEvily, A.J.; Groeger, J. (1978), "On the threshold for fatigue crack growth", Advances in Research on the Strength and Fracture of Materials, Elsevier, pp. 1293–1298, doi:10.1016/b978-0-08-022140-3.50087-2, ISBN 9780080221403
17. Forman, R. G.; Shivakumar, V.; Cardinal, J. W.; Williams, L. C.; McKeighan, P.C. (2005). "Fatigue Crack Growth Database for Damage Tolerance Analysis" (PDF). FAA. Retrieved 6 July 2019.
18. Ritchie, R.O. (1 November 1999). "Mechanisms of fatigue-crack propagation in ductile and brittle solids". International Journal of Fracture. 100 (1): 55–83. doi:10.1023/A:1018655917051. ISSN 1573-2673.
19. Maddox, S. J. (1975). "The effect of mean stress on fatigue crack propagation—A literature review". International Journal of Fracture. 1 (3).
20. Elber, W. (1971), "The Significance of Fatigue Crack Closure", Damage Tolerance in Aircraft Structures, ASTM International, pp. 230–242, doi:10.1520/stp26680s, ISBN 9780803100312
21. Walker, K (1970), "The Effect of Stress Ratio During Crack Propagation and Fatigue for 2024-T3 and 7075-T6 Aluminum", Effects of Environment and Complex Load History on Fatigue Life, ASTM International, pp. 1–14, doi:10.1520/stp32032s, ISBN 9780803100329
22. Dowling, Norman E. (2012). Mechanical behavior of materials : engineering methods for deformation, fracture, and fatigue. Pearson. ISBN 978-0131395060. OCLC 1055566537.
23. Suresh, S. (2004). Fatigue of Materials. Cambridge University Press. ISBN 978-0-521-57046-6.
24. "NASGRO® Fracture Mechanics & Fatigue Crack Growth Software". Retrieved 14 July 2019.
25. "DARWIN Fracture mechanics and reliability assessment software". Retrieved 14 July 2019.
26. Zehnder, Alan T. (2012). Fracture Mechanics. Lecture Notes in Applied and Computational Mechanics. 62. Dordrecht: Springer Netherlands. doi:10.1007/978-94-007-2595-9. ISBN 9789400725942.
27. "Fatigue Crack Growth". Retrieved 6 July 2019.