Preprint
Review

A Literature Review of for Incorporating Crack Tip Plasticity into Fatigue Crack Growth Models

Altmetrics

Downloads

173

Views

96

Comments

1

A peer-reviewed article of this preprint also exists.

Submitted:

22 November 2023

Posted:

23 November 2023

You are already at the latest version

Alerts
Abstract
This paper conducts an extensive literature review focusing on the utilisation of crack tip plasticity as a crucial parameter in determining and enhancing crack growth models. The review encompasses a comprehensive analysis of various methodologies, predominantly emphasising numerical simulations of crack growth models while also considering analytical approaches. Although experimental investigations are not the focus of this review, their relevance and interplay with numerical and analytical methods are acknowledged. The paper critically examines these methodologies, providing insights into their advantages and limitations. Ultimately, this review aims to offer a holistic understanding of the role of crack tip plasticity in the development of effective crack growth models, highlighting the synergies and gaps between theoretical, experimental, and simulation-based approaches.
Keywords: 
Subject: Engineering  -   Mechanical Engineering

1. Introduction

During the last third of the last century and so far, this century, the publication of reviews on crack growth models has increased considerably. Fatigue failure has been a widely researched topic in academia over the past few years, and has been discussed in several publications that cite key contributions from researchers such as Bathias and Baïlon [1] or Suresh [2]. Crack closure and crack opening models have been discussed, [3], focusing on the volume-based strain energy density approach applied to welded structures with notches [4] and the orientation of the critical plane in multiaxial fatigue [5]. Pippan and Hohenwarter presented a review on the controversial phenomenon of crack closure and its physical implications and consequences [6], a study focusing on the theory of critical distances [7], reviews of the application of the Digital Image Correlation (DIC) technique to fatigue [8,9] and reviews of the use of Finite Element Method (FEM) in the study of crack growth [10,11,12]. In another compilation work [13] on notched components, different types of approaches have been presented to characterise notch fatigue behaviours under complex loading histories. Among them, the consideration of the stress gradient effect is a popular perspective, together with the concepts of fatigue damage zone and critical distance.
Fatigue crack growth has been extensively studied using the curves da/dN - ΔK. ΔK is essentially an elastic parameter, while the fatigue crack growth process is related to non-linear and irreversible plastic phenomena at the crack edge. This explains the limitation of da/dN-ΔK models to predict the effect of stress ratio or load amplitude variation. Among the limitations are a) such curves are completely phenomenological, not derived from physics, and the fitting parameters have no physical justification; b) such curves are only valid in the small-scale yield range; and c) da/dN depends on other parameters, including stress ratio and loading history.
Fatigue damage is reflected in plastic phenomena that are not analysed by elastic models, so it is necessary to introduce different mechanisms to analyse the damage that occurs in each cycle. One of them is the cyclic plastic deformation in which the damage accumulation process is strongly influenced by the evolution of the dislocation substructure [14]. The presence of surface oxides in the base material can act as crack initiation points and can significantly decrease the fatigue life of the material [15]. They can also act as brittle inclusions and behave as cracks when loaded normal to their long axes [16].
Brittle solids often exhibit similar micro mechanisms for the growth of static and cyclic cracks at low temperatures [2]. However, the application of a cyclic load on a static mean stress can result in variable lifetimes and, in some cases, the cyclic frequency can influence the crack growth rate. Furthermore, under fully compressive cyclic loading, crack initiation can occur solely due to cyclic loading, even in regions without stress concentrations. These phenomena can significantly influence the growth and failure of brittle solids. Another failure mechanism is 'Fatigue micro voids', which occur because of the separation of non-coherent secondary particles from the matrix [17]. The void formation process is characterised by the initiation, growth and coalescence of multiple interfacial cracks around the particle.
Subsequently, the phenomenon of crack closure is introduced, which was proposed by [18] and has been used to explain the influence of plasticity on fatigue crack growth. It has also been used to determine numerically, experimentally and analytically the effect of stress ratio [19], short cracks [20,21] and model widths [22]. From the first in-depth studies [23] evolved to the concept of partial crack closure [24]. Afterward [25] proposed a new model for the stress around the crack based on four parameters: KF (SIF opening mode), KS (SIF shear), KR (SIF retardation) and T-strain. KR characterises the effect of crack tip protection arising due to plasticity in both the crack tip and the wake. The use of non-linear parameters was introduced to validate the concept of crack closure and identify the most accurate parameters [26]. Experimentally, a relationship between crack closure and crack growth has been observed with different techniques, in [27] DIC is used for this, and subsequently, a model is established to reproduce the results, obtaining information of the closing stresses. In another study [28] a direct relationship between overloads, crack closure and how it affects crack growth is observed. In the works [29,30] a relationship is established between plastic wake sizes, Crack Opening Displacement (COD) and effective intensity factor at the crack tip. Following that, these parameters have also been related to crack closure with DIC experimental techniques [31]. The DIC technique has been shown to be a very useful tool to determine the elastic or elastoplastic displacement field in the crack environment, to later, through mathematical, analytical or numerical models, calculate the key parameters of the Fatigue Crack Growth Rate (FCGR), and establish the possible driving forces of crack growth. Only in this publication [32] the authors evaluate 5 different methods to locate the crack tip position by applying DIC. The methods are: 2 constrained Newton, the trust-region reflective and the quasi-Newton; the Nelder-Mead Simplex; a constrained Genetic Algorithm (GA) and finally a constrained Pattern Search (PS).
The X-DIC technique [33] has also been shown to be efficient to evaluate the plastic state of the crack tip. Today, there is still a lot of controversy surrounding the crack closure phenomenon, on the one hand, many researchers deny that it has any influence on crack growth. On the other hand, Pippian and Grosinger [34] proposed an original approach to crack closure phenomenon, extending its importance and relevance, at the time of the research publication, the acceptance of crack closure phenomenon and its effect on FCGR was restricted to Low Cycle Fatigue (LCF) processes. Based on images from a scanning electron microscope (Figure 1c) and the interpretation of ΔJ and ΔJeff in the cyclic stress-strain hysteresis curves characteristic of fatigue processes (Figure 1a,b), the authors demonstrate that crack closure also occurs in High Cycle Fatigue (HCF) processes and affects ΔJ and ΔJeff, consequently affecting the FCGR.
The study of the effects of elasto-plastic behaviour on crack growth is approached from many points of view and mechanisms. For example, the effect of the sharpening and hardening that occurs at the crack tip on the crack development itself is addressed [35]. The effect of crack tip sharpening and hardening under maximum and minimum loads has been investigated [36]. This phenomenon has also been experimentally verified [37] with the use of micro fractography. A mathematical model is proposed [38] based on time differential equations describing the thermodynamic process of plastic energy dissipation during the hardening and sharpening of the crack tip growth. Another analytical approach [39] is based on the use of the Smith-Watson-Topper parameter assuming that in the first loading cycles a large and irreversible sharpening of the crack tip occurs.
Another important element to consider is the identification of the crack tip plastic zone. During fatigue crack propagation, two plastic zones are created at the crack tip [40]. One is constant during loading, while the other is cyclic during unloading. Research has shown that the size and shape of the plastic zone significantly influence the crack growth behaviour.
It should also be mentioned that there are numerous papers whose focus is not on the relationship between crack growth and crack tip plasticity, but which provide evidence of the connection. In the work [41]whose main objective is to validate crack growth models (Fatemi-Socie) under biaxial loads and overloads, it was clearly observed that the combined application of axial and torsion almost eliminated crack closure, thereby increasing the crack growth rate. In the following work [42] a finite element model is carried out with Abaqus and the effect of the position of holes drilled in a plate to reduce the growth rate of future cracks is evaluated. The optimum location and size are sought, reducing the rate of crack growth by 50%. To see the effect of the circular holes drilled in the crack tip region, a direct evaluation of the ratchet limit and the plastic deformation range of the crack tip is made. This evaluation is achieved using the new Linear Matching Method (LMM). There is a situation that optimises the slowing down of crack growth because the plastic deformation decays by 50%.
This article systematically reviews recent advances in notch fatigue analysis in relation to stress gradient effects. Despite these advances, a comprehensive analysis of published methods that use crack tip plasticity to determine or improve crack growth models has not yet been carried out. Next, different researches, that propose methods for determining crack growth, will be presented. Works that establish quantifiable relationships between plasticity and crack growth will also be presented. Although they do not present a complete model, they detect and quantify parameters of interest and are the seed of possible future models when all the parameters are found, which, as they are unknown, do not allow the formulation of the complete model. The works have been grouped according to the theories used, as well as the damage parameters. The organization of the review sections is based on the chronological order of the publication year of the works under examination. Additionally, within each section, the presented research is arranged in chronological order to facilitate the comprehension of each model's timeline.

2. Cumulative Damage Theories Based on Stress and/or Strain History

Building upon the foundational theories in the field of crack growth, a variety of models have been proposed, each extending and refining the understanding of the underlying mechanisms. These models, though varied, share a common lineage in their approach to the study of cumulative damage based on the stress and strain history in fatigue crack growth.
Focusing on the primary aspects of cumulative damage theories, one of the earliest works was proposed by Shih [43]. This work presents a relation between the J-integral and the Crack Tip Opening Displacement (CTOD) obtained, based on the description of the stress and strain state proposed by Hutchinson-Rice-Rosengren with singularity at the crack tip. A similar expression is formulated for the CTOD. The author arrives at a relationship between CTOD and J in the equations. The coefficient d, which relates J, depends on the deformation properties of the material and is independent under small-scale plasticisation conditions. But the obtained relations are not valid for large-scale plasticisation due to the singularity. Based on deformation theory, the author arrives at a relation between J and the parameter d, eq. 1, which is a candidate for a more robust model. In fact, after numerical correlations, the model is found to be good at all scales.
d δ / d a = n n + 1 d n σ 0 d J d a
Building on this foundation, Noroozi et al. [39] introduced a new perspective by presenting a FCGR model that hinges on the elastoplastic stress-strain history at the crack tip. The driving force of crack development is derived based on the damage parameter used by Smith-Watson-Topper. The internal effect of stresses induced by the inverse plasticity cycle was also analysed. All these derives to eq. 2, where the final model is presented as a function of the two parameters and where the authors present an extensive network of mathematical developments and balances of forces and stresses. Evidently, the complexity of this model lies in expressing this history of stresses and strains, so the focus is on the maximum and minimum stresses and strains at the crack tip (in eqs. 8-16 in the original publication). Special attention is paid to the contribution of compressive stresses (in eqs. 17-24 in original publication), while elastoplastic stresses and strains in the crack tip environment are described (in eqs. 27-32 in original publication). The calculation of the residual intensity factor Kr is addressed (in eqs. 33-34 in original publication) and all other necessary intensity factors are listed (in eqs. 35-40 in original publication). With all the above, the fatigue crack driving force ΔK and the fatigue crack growth expression da/dN-ΔK are derived analytically (in eqs. 41-60 in original publication). Figure 2 shows a high correlation of the model results with the experimental ones. The authors point out that the fact that being able to relate the Smith-Watson-Topper damage parameter and the stress-strain history of the crack allowed introducing residual stresses into the model has improved its goodness of fit.
d a / d N = f ( ( K m a x , t o t ) ρ ( K t o t ) 1 ρ ) γ
In a similar vein, de Matos and Nowell [44] extended previous models by integrating the plane strain Plastic-Induced Crack Closures (PICC) model for more complex geometries, such as flat plate with a circular hole and two radially symmetric cracks, and residual stress fields. The present model is improved because it can be extended to any geometry and previous residual stress fields. The analytical model starts from expressing the stress state as well as possible residual stresses with boundary elements and plane stress, the objective function is set as the crack growth. A finite element model is implemented and solved for both elastic and elasto-plastic cracks, with and without residual stresses, presenting the results in Figure 3, which correlate the results of the analytical models with experimental.
Extending the exploration of cumulative damage theories to the dynamic interplay between plastic zones and crack growth, Chikh et al. [40] presented a study focusing on this relationship. Since two plastic zones are produced under cyclic loading, one monotonic (rm) and one cyclic (rc ), this work relates this second zone rc to the FCGR, which is produced by the local compression generated at the crack closure of each cycle. Several models can be used to relate the plastic zone to crack propagation. In this study, a simplified model based on the effect of damage accumulated in front of the crack tip da/dN=f(rc) is used. Taking the expression of the stress value around the Linear Elastic Fracture Mechanics (LEFM) crack, it implies that a plastic deformation zone is formed around the crack which has an extension relative to the value of the applied load and the crack length. Furthermore, if the plastic zone is large, a large amount of energy is absorbed during crack propagation, thus, if the zone is small, crack propagation requires less energy, the size of the plastic zone is directly proportional to the hardness of the material. Therefore, the mechanical properties of the material and the stress state govern the size and shape of the plastic zone. It is concluded that the plastic zone depends on the stress field, the applied stress, the specimen thickness and the crack length, and that the crack size is an ideal fatigue crack propagation parameter to determine crack growth, better than any other Elasto-Plastic Fracture Mechanics (EPFM) parameter. The proposed model is based on da/dN=B(rc) 2 where B=3.2E-6(R) e1,3, being Re the yield stress. Experimental tests were carried out on 12NC6 steel yielding showing consistent and conservative results.
In the context of high-temperature conditions, Tong et al. [45,46] bring a novel perspective by incorporating finite element numerical modelling. They also present further results from Finite Element numerical modelling at high temperatures (650 °C). The modelling of the material follows Chaboche's laws [47] and isotropic kinematic hardening (Power Law) and they show the results after the application of 30 loading cycles under different values of load, frequency and load Ratios (R). Among the main conclusions, (a) the ratchetting deformation increases as the R decreases, being higher in negative cases, (b) low frequencies in the loading cycle increase the accumulated deformation, as well as the crack growth rate and (c) maintaining an initial load (Dwell) decreases the fatigue life of the specimen; for example, a Dwell of 100 s at maximum initial load reduces the fatigue life by one third. An interesting factor studied [45] is to use meshing to define different grain sizes and orientations.
Tong et al. [48] further elaborate on the significance of Ratchetting Strain as a damage parameter in their subsequent work. They establish the Ratchetting Strain as a driving force for fatigue crack growth. For this purpose, several models, both independent and time-dependent, are developed using elasto-plastic, visco-plastic and crystal-plastic models [49]. The authors argue that Ratcheting stress always occurs near the crack tip, leading to the progressive accumulation of stress strain normal to the crack growth plane. It seems plausible that this tensile strain, or Ratcheting strain, may be responsible for the material separation leading to crack growth. As mentioned above, the authors introduce several different material models. Subsequently, finite element analyses of the elasto-plastic, visco-plastic and crystal-plastic models are performed, and the results are shown. The correlation of the results obtained with the application of the models with experimental results is presented [49], showing a special agreement when used for the cumulative plastic deformation model. The results obtained for the analytical and numerical models also correlate well. Further experimental results supporting this theory have been published [50].
Complementing this body of work, other research [51] offers a unique perspective by considering crack growth as a process of successive restarts, akin to the initiation stage of a crack. A crack growth model based on accumulated damage was presented in a similar way as a fatigue analysis of stationary notches. Crack growth is considered as the process of successive cracks restarting, as if it was always in its initiation stage. The driving force of crack growth is based on the local stresses and strains around the crack tip, using the Smith-Watson-Topper (SWT) fatigue damage parameter D=σmax -Δε/2. The driving force of crack growth can be expressed as a function of two parameters. On the one hand, the residual cyclic plastic stresses σres (χ) and on the other hand, the residual intensity factor Kr which are expressed in different equations. The study shows the correlation of the experimental results and those obtained with the model, validating the use of the model which is based on modelling the crack as a succession of blunt notches of radius rp, and not as a sharp crack, therefore with realistic stresses and strains around a rounded crack tip.
A different approach is seen in the work [52] where the focus is on identifying decisive parameters in fatigue crack growth based on accumulated plastic deformation. In this work [52] a numerical study based on the accumulated plastic deformation is presented to determine which parameters are decisive in fatigue crack growth. As main results, firstly, it is established the need to stabilise the crack growth with initial cycles where the da/dN is not stable. Secondly, it is compared with experimental results by defining the percentage of accumulated plastic deformation as the driving parameter of crack growth.
Finally, another paper [53] utilises the DIC technique to offer a unique perspective in determining the displacement field around the fatigue crack tip. In detail, DIC technique was used to determine the displacement field around the fatigue crack tip for a constant value of the Stress Intensity Factor (SIF). The data obtained by DIC was used in the evaluation of the T-strain and to see its influence on the FCGR. As an original contribution of the study, the higher order terms of the Williams Expansion (WE) are determined. The displacements of a set of points outside the plastic zone were selected for the application of the Over-Deterministic Method (ODM) to obtain several initial WE term. A numerical model of boundary conditions is also made. The calculated experimental T-stress values show a good agreement with finite element analysis and literature. It was also shown that the level of restraint influences the fatigue crack propagation rate, the higher the T-stress, the lower the growth.
It becomes evident that the explored models, while distinct in their theoretical underpinnings, converge in their emphasis on the stress and/or strain history's role in fatigue crack growth. Initiating with Shih's [43] foundational work that intricately links the J-integral and CTOD, this section traverses through various models, each underscoring different aspects of crack propagation. The intricate elastoplastic stress-strain history analysis by Noroozi et al. [39], the extended application of the PICC model by de Matos and Nowell [44], and Chikh et al.'s [40] investigation into the plastic zone's influence on crack growth rates, collectively illustrate the complexity and multifaceted nature of crack propagation studies. Further, the models by Tong et al. [45,46,48] and the novel applications of the DIC technique [53] highlight innovative approaches to understanding the dynamic interactions at the crack tip.

3. Damage Theories Based on Crack Growth Concepts. Ctod, Plastic Zone Size, rp, ΔK,

Following the review, a spectrum of damage theories based on crack growth concepts is presented, emphasizing the interplay between plastic zone size, CTOD, and SIF. Each study contributes to a more integrated understanding of fatigue crack growth, weaving together different theoretical approaches with experimental insights.
This was one of the first works [54] to successfully use plastic parameters to model the FCGR. Specifically, in this work, the size of the plastic zone rpc is used as a key parameter. Starting from the equation of rpc based on the Dugdale model, and applying plane strain conditions, eq. 3 is obtained for rpc, and with plane stress conditions, eq. 4 is obtained. That is, a specimen of zero thickness would be ideal for plane strain and one of infinite thickness for plane strain conditions. The FCGR is expressed as da/dN=Ce (ΔKeff) mc. Where Ce and mc are material constants. Figure 4a,b show the correlation of the experimental FCGR with that obtained in the model for different specimen thicknesses, as a function of ΔK and ΔKeff. Figure 4c,d show the FCGR correlations assuming plane strain and plane stress respectively, it is observed that plane strain conditions favour the FCGR, so it increases with specimen thickness. Following the emphasis on the plastic zone size in [54], the study in [22] extends the discussion by measuring crack length and fatigue crack closure to understand the variation of crack growth rate with SIF under different loading conditions.
r p c = π 18 K 2 σ y 2 1 2 + 2 3 2 3 σ a p p σ y
r p c = π 8 K 2 σ y 2 1 2 σ a p p σ y 1 3 4 σ a p p σ y 2
The objectives of the work [22] were to measure crack length and fatigue crack closure to analyse the variation of crack growth rate with Kmax or ΔKeff for different values of R (stress ratio) and B (specimen thickness), to obtain the da/dN-ΔKeff curves, and to propose an empirical model that correlates a new normalised loading parameter U with R, B and ΔK. The results are obtained for constant load amplitude in tension with 3 stress ratios: R=0, 0.2 and 0.4 and 3 specimen thicknesses: B=6, 12 and 24 mm. To measure the crack opening values, gauges are used, and with them ΔKeff. can be calculated. The authors start from a set of equations (in eqs. 1-3 in original publication) to relate in first instance U, Kop and ΔKeff. . In view of the experimental results (shown in Figure 4, Figure 5 and Figure 6 in the original publication) an experimental relation between da/dN and ΔKeff and a readjustment of U is proposed. As a final model, the relation between da/dN and ΔKeff in eq. 5 remains. An important conclusion obtained in the work is that, despite some scatter in obtaining the da/dN - ΔKeff curves, the results of the stress ratios and thicknesses are presented by a simple curve for two parameters of the crack growth rate relationship applied to CK45 steel. A model of U as a function of R, B and ΔK has been proposed based on the experimental data of the crack closure. This model is limited to 0<R<0.4 and 6<B<24 mm.
d a / d N = 9 · 10 9 K e f f 2.95
Building on the understanding of crack closure and stress intensity factors from [22], at the end of the last century [55] Donald and Paris introduced the Adjusted Compliance Method (ACM) and the Crack Wake Influence (CWI), further refining the predictive models for FCGR. In this research, a method for calculating the ΔKeff was established based on the ACM. Subsequently, the CWI is used to understand the intrinsic limitations of the ACM, and by applying adjustments to account for the CWI and its mechanisms, develop an improved model. Starting from the expression of ΔKeff =K -Kmaxop, corrections are applied to arrive at eq. 6 where the parameters exhibit a high degree of similarity (ΔKeff, Kmax and Kop,) but with different nuances.
Δ K e f f = Δ K o p + Δ K m a x Δ K o p · 1 2 π
Δ K e f f = Δ K a p p · 1 V c w i V · I F R
The application of ACR leads to the expression ΔKeff =ACR-ΔKapp, where ACR is a coefficient. The CWI method will be modelled by subtracting Kcwi from K to obtain Keff, Kcwi is related to Vcwi which is the midline shift in a series of equations to obtain the final expression of Keff in eq. 7 Two further modifications are proposed based on the adjusted compliance ratio method and the adjusted compliance ratio/opening load blend method. From this point on, the study focuses on representing the FCGR versus the ΔK, both experimental and those obtained in the models, an example of which is shown in Figure 5. The FCGR will be the da/dN, expressed in mm/Cycle, and the ΔK, the diverse intensity factors proposed in MPa(m)0.5. In total 7 different methods of estimating the effective intensity factor are compared. The ACR method is easy to apply but does not reflect crack closure well. The CWI is mathematically more complex to apply, limiting its applicability. A second version, ACRn2, improves the problems of ACR without adding mathematical difficulty and finally, the empirical approach for AOP (simply an empirical combination of the ACR and OP methods) makes it a simpler model.
The evolution of the partial crack closure model by Donald and Paris [55] is further advanced by [24], who proposes an improvement to integrate the effects of the R-ratio on crack growth, thereby addressing a crucial aspect that had been systematically neglected. The researcher [24] proposes an improvement to the partial crack closure model proposed by Donald and Paris. The improvement consists of a transition function that evolves from the Donald and Paris partial closure model [55,56] to the model in the region close to the Elber growth threshold [57]. Elber proposed the ΔKeff =K-Kmaxop, establishing an effective intensity factor that in turn depended on the intensity factor at crack opening. The model proposed in this paper is developed (in equations 3-10 in the original publication), resulting in eq. 8 for ΔKeffM, which adequately correlates the effects of R-ratio on crack growth at low and intermediate load levels in aluminium alloys. Until the publication of this research the effect of R-ratio had been systematically neglected, it is important to emphasize that the R at the crack tip does not necessarily coincide with that applied to the specimen. It is noteworthy that the key parameter in this model, ΔKeffM, is related to the crack opening profile which obviously depends on the plastic deformation at the tip. Experimental tests show that the proposed model shows improvements to the model in the Elber growth threshold region.
Δ K e f f M = K m a x K o p 1 + 2 π 1 e K m a x K m a x T M 1
Complementing this approach, Zhang et al. [58] introduced a model based on dual parameters, da/dN and da/dS, to define the crack propagation velocity, thus broadening the scope of parameters considered in crack growth modelling. Zhang et al. [58] presented a model based on two parameters, the classical da/dN and the novel da/dS that defines the crack propagation velocity with respect to the applied stress at any time of the cycle. Experimental methods were used to predict the crack propagation behaviour. The mentioned parameters are related and, developing mathematically, eq. 9 is reached where da/dN is expressed as a function of ΔK, α, β and, where α and β can be estimated from the results of two constant stress amplitude fatigue tests. The correlations of the experimental results with those proposed by the models are shown in Figure 6.
d a / d N = C Δ K 2 β + α + 1 1 R 2 α + 1 1 R 2 α + 1
Further contributing to this comprehensive analysis, Lopez-Crespo et al. [59] undertake an analytical-experimental study focusing on the impact of crack closure and plasticity at the crack tip, thereby offering insights into Mode I and II intensity factors. The experimental laboratory work obtained the displacement field at the crack tip and its environment with DIC techniques on two types of Centre-Cracked Plate (CCP) and Wedge Opening Load (WOL) specimens that had been previously initiated with a crack itself, not a notch. The analytical development is based on formulating the displacement field on the specimen surface as two-dimensional, u and v being the horizontal and vertical displacements respectively and expressing them as a complex number z=x+iy. In turn, these displacements can also be obtained from two analytical Muskhelishvili functions φ = φ(z) ψ=ψ(z) (eq. 10). These expressions can be derived by obtaining expressions as a function of u, v and the coefficients (C, D, E, F, α and β). Applying basic tensor algebra, it can be expressed in a matrix form Aχ=b where b is the vector that depends on the displacements u and v. From there it is straightforward to calculate the coefficients (C, D, E, F, α and β) of the matrix A and back in eq. 10 calculate the analytical functions φ = φ(z) and ψ=ψ(z). Ultimately the complex stress intensity factor K=KI -iKII can be calculated with the above data using eq. 11. Figure 7 plots the intensity factor K, the effective intensity factor Keff and the real and imaginary components of the intensity factor KI and KII. The difference in values obtained with the model and those obtained experimentally exists as the model lacks the crack closure effect, while experimentally there were clear indications of its existence and possible effects.
z = ω ζ = R ζ + ζ m ; φ ζ = + a k ζ k ; ψ ζ = + b k ζ k ;
K = 2 2 π lim ζ 1 φ ' ζ ω ζ ω 1 ω ' ζ
On a similar note, Zhang [60] develops a fatigue crack propagation model emphasising the effect of compressive stress, thereby adding another dimension to the understanding of fatigue crack growth under tension-compression loading. The stress near the crack tip, the displacement and the size of the plastic zone are obtained from a kinematically hardened material. The results show that the above three parameters continue to change with the change of applied compressive force. Experimental methods were used to predict the crack propagation behaviour. The study is based on analyses made by Silva [61,62] who found that materials exhibiting strong cyclic hardening and a high Bauschinger effect (such as Al 2024-T351, ck45, 0.4% mild steel) were strongly affected by the applied compressive load while materials not exhibiting cyclic hardening (such as Ti6Al4V and Al7175) were relatively insensitive to the applied compressive load. Silva also concludes that models based on the concept of fatigue crack closure are not suitable, while models based on materials with cyclic-plastic properties are suitable for describing fatigue crack propagation behaviour under tension-compression loading. Based on this, Zhang et al. develop a model where they identify the internal and external parameters that control the fatigue crack growth behaviour and describe the fatigue crack propagation rate under tension-compression loading. The paper concludes that using a kinematically stiffened material subjected to tension-compression loading: the compression part of the applied cyclic load has a significant effect on the crack tip stress, displacement and plastic deformation field, and that the maximum applied stress intensity factor kmax and σmax-comp (maximum applied compressive stress) are the two external loading parameters that determine the above three variables. As an inner parameter, rpc, the reverse plastic zone size is appropriate to correlate the fatigue crack propagation rate, see also Figure 8. Finally, he proposes an equation for the fatigue crack growth rate under tension-compression loading, with special emphasis that it has nothing to do with the concept of crack closure, but with plasticity at the crack tip and based on the plastic damage theory [58].
Additionally, Pokluda [63] presents a discrete model for small plastic deformations, offering a unique perspective on induced crack closure components, such as plasticity and roughness, and their influence on crack growth. The researcher [63] presented a discrete model for small plastic deformations under plane strain conditions that allows direct evaluation of the magnitude of induced crack closure components, such as PICC and Roughness-Induced Crack Closure (RICC). It also determines the influence of the microstructure on the induced roughness. The advantage of the dislocation-based method over those based on continuous mechanisms is its "physical transparency" which allows the evaluation of both induced plasticity and roughness at the same time. Based on measurements taken on an array symmetrical to the crack, experimental measurements are taken, and plane strain is allowed for the interpretation of the measurements. In this case, key is to know the stress state on the crack flanks (position of the array), this stress state is the cause of the PICC and RICC. A model is proposed under certain hypotheses for the calculation of the PICC and RICC components, ΔK, Kmax, Keff, Kcl, rs, Kcl, rl and Kcl, p. The subscripts cl stands for closure, rs range-short, rl range-long, p plasticity. Once these characteristic parameters of RICC and PICC are calculated, they can be related to crack growth as deemed best. Finally, ΔKeff = Kmax -Kcl; Kmax =ΔK/(1-R), and the expression is given in eq. 12.
Δ K e f f = 1 C η R A 2 1 3 η ( R A 1 ) 2 6 + 6 ( R A 1 ) 2 C Δ K 1 R
To talk about numerical models of crack growth and their relation to plastic phenomena in the crack tip environment is to talk about Fernando Antunes. Building on the discrete approach by Pokluda [63], Antunes [64] presents an elasto-plastic numerical model, which provides a different perspective by excluding the contact between the crack flanks, thereby focusing on the relationship between CTOD, plastic phenomena at the crack tip, and crack growth rate. The work [64] presents an elasto-plastic numerical model where the contact between the crack flanks is obviated. In all simulations, the CTOD and CTODP were measured. A relationship between da/dN and CTODp was established for the AA6082-T6 aluminium alloy, which is a quantifiable relationship between plasticity and crack growth rate, on the other hand, the relationship between CTODP and other parameters of interest such as ΔK was quantified. When modelling without contact (visible in Figure 4) both plastic and elastic CTOD are higher, i.e., without contact, without the effect of crack closure, the effect of loading is higher. Without contact, without crack closure, the crack propagation only depends on ΔK. The CTODp is plotted with respect to the deformation, linear relationship, and with respect to the energy, obtaining a quadratic relationship. In this work [26] the main objective is to verify the effectiveness of the crack closure concept, therefore, the crack flank contact is modelled, as well as the following non-linear parameters in the crack tip environment are calculated: cyclic plastic strain range, both elastic and plastic CTOD, size of the plastic inverted zone, rpc in Figure 8, and the plastic energy dissipated per cycle. Direct relationships of these four parameters with the crack growth rate and ΔK are shown throughout the paper, which can be said to form a plasticity-based growth model. The paper also demonstrates the effect of crack closure on the linear parameters, as well as the effect of mesh size on the linear parameters. In another paper [65] the same numerical finite element model is basically used as in the two previous works, which will be now described briefly: the numerical modelling is assumed to be elastoplastic with: 1-Isotropy and 2.- The Von Mises criterion is followed in plasticity together with the Voce isotropic hardening law [66] and the Lemaitre-Chaboche’s kinematic law of hardening [47]. The numerical model is implemented with the DD3IMP software. Only 1/8 of the specimen is modelled by symmetry with the consequent computational saving. The main result of the work is the relationship between the COTDP with da/dN for 7050-T6 aluminium. Specifically, there is a linear relationship, just multiplying by 0.5245. Following the basis of the numerical modelling of the previous publications, in the following work [67] the relationship between plasticity and crack growth is further explored, the relationship between CTODP and da/dN for Aluminium 2050-T8 is also obtained. As a novelty in this work, the effect of the number of load cycles that elapse until another node is released from the crack tip with the consequent crack growth is studied. The higher the number of load cycles between node opening, the higher the CTODP. The da/dN is plotted against the numerical CTODp, and the strong influence of the numerical parameters can be seen. In plane strain, the higher the number of cycles, the lower the CTODp, unlike in plane deformation. In the work [64] the methodology for the calculation of the Plastic CTOD in numerical tests is established in detail, as well as the parameters that will affect this calculation, such as the crack propagation distance, the size of the crack, the distance of the points with respect to the edge, the size of the elements and the number of previous cycles. The mesh size and the software used are those mentioned earlier in Antunes' work. The distance of the first released node after the crack edge is the most sensitive parameter together with the mesh size.
Continuing the theme of elasto-plastic analysis, Kawabata [35] proposes a method for CTOD calculation that incorporates the impact of strain hardening at the crack tip, further enhancing the understanding of plastic effects on crack growth. This method for CTOD calculation [35] is established on the crack tip edge loss due to strain hardening. The researcher introduces a new factor f to account for the plastic factor. This f-factor is a function of the plastic radius and the thickness of the specimen. The result of the correction is verified numerically and experimentally. Starting from the proposed and revised formula [68,69,70,71] for the CTOD calculation, they concluded that the CTOD calculation was overestimated due to the presence of the plastic deformation at the crack tip, which was expected with respect to the edge loss. The CTOD is the sum of an elastic and a plastic term, the elastic term is calculated based on the intensity factor K, and the plastic term based on the Hinge plastic model. For the calculation of the CTOD there are several parameters to adjust: in the plastic part m(YR) is calibrated with eq. 13, the factor f (YR, B) of the plastic part is calibrated in eq. 14, which gives as final result the calculation of the CTOD by eq. 15. The goodness of the model is checked numerically and experimentally, the correlation is shown in Figure 9.
m Y R = 4.9 3.5 Y R
f Y R , B = f B f Y R = 0.8 + 0.2 e 0.019 B 25 1.4 Y R 2 + 2.8 Y R 0.35
δ j w e s = δ e l δ p l = K 2 1 ν 2 m Y R E σ Y S + f ( Y R , B ) r p W a 0 r p W a 0 + a 0 V p
Further expanding on the theme of CTOD, Vasco-Olmo et al. [72] conduct an experimental study to evaluate the capability of CTOD in characterising fatigue crack growth, thereby exploring the practical applications of the theories discussed in previous works. The research [72] conducted an experimental study of CTOD, it is performed to evaluate the ability of this parameter to characterise fatigue crack growth. A methodology is developed to measure and analyse CTOD from experimental data. Displacements are measured by implementing DIC in the fatigue crack growth, in this way the CTOD can be determined. Fatigue tests were performed with R of 0.1 and 0.6 on Compact Tension (CT) specimens made of commercial pure titanium. A sensitivity analysis was performed to explore the effect of the selected position behind the crack tip on the CTOD measurement. Analysis of a full load cycle allowed the elastic and plastic components of the CTOD to be identified. For the plastic, CTOD was found to be directly related to the plastic deformation at the crack. Furthermore, a linear relationship between da/dN and plastic CTOD was observed in both tests (Figure 10). The results show that CTOD can be used as a viable alternative method to ΔK in characterising fatigue crack propagation because the parameter considers the fatigue threshold. This work aims to contribute to a better understanding of the different mechanisms driving fatigue crack growth and the direction of fatigue crack growth, a controversy associated with plasticity-induced fatigue crack closure.
Lastly, Medhi et al. [73] integrate DIC with analytical solutions to monitor structural health, focusing on the stress intensity factors under cyclic loading, thus connecting the theoretical models with practical structural health monitoring applications. This methodology for structural health monitoring is based on a combination of DIC and an analytical elastic solution. Experimentally (DIC), the full displacement field around the crack tip in a CT specimen made of 2024-T351 aluminium alloy under cyclic loading at different load levels is obtained. Then, an over deterministic multipoint method is used to calculate the SIF in Modes I and II (KI and KII in eq. 5 in the original publication), to apply the method it is necessary to adjust the experimental displacement field obtained with the series of William [74]. The parameter was determined to facilitate the difference between the SIF result of the (nominal) model and the experimental one. It was observed that, with increasing amplitude in cyclic loading, the difference between nominal and experimental K increased. This may be because of crack tip plasticity, which was not considered in the nominal evaluations. To account for plasticity, Irwin's approximation was used in the analytical model. The plastic term is added to the crack length, acorr =a+ry, where ry is the plastic radius given by the Irwin approximation, r y = 1 / 2 π · K I σ y s 2 . The results showed a better agreement in the evolution of K under cyclic conditions, increasing the load up to the sharp fracture of the specimen.
This section intricately navigates through diverse damage theories centred on crack growth, weaving together the significance of plastic zone size, CTOD, and stress intensity factors. Commencing with the pioneering use of plastic parameters by [54], the discussion evolves, with [22] exploring crack length and closure. This narrative progresses with [55] refining FCGR models through the ACM and CWI methods. The evolution continues with [24] enhancing crack closure models by integrating R-ratio effects, and Zhang et al. [58] expanding the parameters in crack growth modelling. Complementary studies by Lopez-Crespo et al. [59], Zhang [60], Pokluda [63] and Antunes [64] delve into the nuances of crack closure, plasticity, and microstructural influences, while Kawabata [35] and Vasco-Olmo et al. [72] focus on CTOD methodologies. Lastly, Medhi et al. [73] bridge these theoretical models with practical structural health monitoring applications. Collectively, these studies form a cohesive narrative, significantly enriching the understanding of fatigue crack propagation and the importance of CTOD and crack tip plasticity that are demonstrated through experimental work [64,72,73].

4. Energy-Based Theories of Damage

Expanding on the concept of energy accumulation and/or dissipation as a fundamental factor in fatigue crack growth, this section delves into various models that utilise energy-based theories to elucidate the damage process in ductile solids. Each model, while grounded in the same theoretical domain, offers unique perspectives and methodologies, contributing to a broader understanding of crack propagation dynamics.
In one of the earliest researches, Klingbeil [75] proposed a new theory of fatigue crack growth in ductile solids based on the total plastic energy dissipation per cycle ahead of the crack. The fundamental hypothesis of the theory proposes a unified criterion for crack growth under monotonic conditions and fatigue loading, so that the fatigue crack growth rate is given explicitly in terms of the total plastic dissipation per cycle and the monotonic plane strain or plane stress, and the fracture toughness of the material. Here, the total plastic dissipation per cycle is obtained by 2-D elastic-plastic finite element analysis of a stationary crack with constant amplitude loading, Mode I CT specimens. Initially, the aim is to relate exclusively the fatigue crack growth to the total plastic energy dissipated, emulating the concept of the J-integral as a crack growth parameter. For the calculation of the dissipated energy, a stationary numerical model is proposed in which the possibility of crack closure is not given. The FEM was implemented with ABAQUS, under small elastoplastic deformations and Von-Mises criterion, reduced integration to abolish shear stiffness and minimum element size less than 0.05 µm. The dissipated energy is directly calculated by ABAQUS as shown in Figure 11b, and Figure 11a shows how the parameters affect the energy. It was obtained as a result that the dissipated energy decreased with hardening. As a result, eq. 16 fits a crack growth model as a function of dissipated energy.
d a d N * = d a d N K I c σ y 2 d W d N *   a n d   K * = K / K I c
Transitioning from the specific focus on plastic dissipation in [75], Pommier and Risbert's study [38] introduces a more holistic approach by considering a supplementary state variable for the crack, enhancing the understanding of the energy dynamics in fatigue crack growth. Pommier and Risbert [38] presented a paper where the main objective was to propose a set of equations time derived for fatigue crack growth in order to avoid any deconstruction of the cycle. The model is based on the thermodynamics of dissipative processes. Its main originality lies in the introduction of a supplementary state variable for the crack, which allows the crack state to be described continuously throughout any complex loading sequence. The crack state is fully characterised on a global scale by its length a, its plastic blunting ρ and its elastic opening ∆CTOD. The plastic blunting is calculated by integrating the difference of the elastoplastic and elastic displacement at the crack tip along d, and dividing the result by d. The model starts from an energy balance of the crack where D is the dissipated energy which is formulated in eq. 17. Together with the Clausius Duhem inequality they are combined in the expression of da/dN in eq. 18. The model finally consists of two laws: a crack propagation law, which is a relation between d ρ / d t and da/dt and observes the inequality derived from the Clausius Duhem inequality, and an elasto-plastic constitutive behaviour of the cracked structure, which gives d ρ / d t versus load and is derived from the energy balance equation. By solving the set of equations with the appropriate boundary conditions, the evolution of the different crack parameters with time is obtained. The model was implemented and tested. It successfully reproduces the main features of fatigue crack growth as reported in the literature, such as the Paris law, the stress ratio effect and the overload retardation effect.
D = ϕ a e f f d a d t + ϕ ρ e f f * d ρ d t + ϕ a χ d a d t + ϕ ρ χ * d ρ d t
d a d N = t ( N ) t ( N + 1 ) d a d t d t = α 2 ρ m a x ρ m i n α 2 C T O D 2
Further advancing the energy-based perspective, the subsequent study [76] shifts focus to the cumulative change in the cyclic strain energy, offering a complementary view on how energy parameters drive crack growth under various loading conditions. The aim of the work [76] was to demonstrate that fatigue crack growth under elasto-plastic conditions can be perfectly correlated using the concept of physically finding a driving force parameter. In this work the selected parameter is the cumulative change in the cyclic strain energy of the net section, thereby finding a physical approach independent of geometry and loading conditions, establishing a parameter as the driving force of the crack, which is calculated analytically from the elastoplastic behaviour of the material and the relative sizes of the crack sections in their corresponding planes. The author demonstrates and calculates the stresses (σ), strains (ε), elastic and plastic energies (U), and idem with the energies of the net cross-sections (UL, in eqs. 1-14 in the original publication). The applied loads will be cyclic in both, stress and strain control. The author presents an approach based on the treatment of elasto-plastic deformations in the net section, from the material resistance point of view. The physical driving parameters of the crack are calculated as the difference of the energy of the net section and the energy before the crack appeared. ΔC, eq. 19, for stress-controlled fatigue and ΔC , eq. 20, for strain-controlled fatigue. Both parameters are expressed in terms of material hardening and deformation parameters (k, n), crack size (a) and specimen width (W). It achieves excellent correlations between the crack growth ratios as a function of the strain energy parameters in the net section as seen in Figure 12. Furthermore, the correlations have been obtained without any consideration of the concept of crack closure. As a conclusion, it was also added that this approach can be used in many inelastic deformation situations, such as creep crack growth and crack growth in viscoelastic materials.
Δ C p σ s L = n σ n + 1 n σ n + 1 n o K 1 n n + 1 W W a 1 n 1
Δ C p e L = K n + 1 p n + 1 o n + 1 a W a
Finally, [77] propels the energy-based discussion forward by proposing a comprehensive energy equation for fatigue crack growth. This approach seeks to quantify the various energy contributions per cycle, thus providing a more detailed understanding of the energetic interactions during crack propagation. It is proposed [77] an energy equation for fatigue crack growth. It equates the total external work per cycle (dW/dN) to the sum of the plastic dissipation (dUpl /dN), the increase in stored strain energy (dUe /dN) and the energy dissipated in the formation of new crack surface (dUa /dN). Experimental measurements of fatigue crack growth were performed to obtain the relationship between fatigue crack growth rate (da/dN) and energy variables. The result shows that dUpl /dN and dUe /dN are not directly related to da/dN. The dUa /dN, whose value cannot be obtained experimentally with sufficient regression, may be the variable directly related to da/dN within the test range.
A comprehensive examination of energy-based theories in the context of fatigue crack growth is presented, with each model contributing uniquely to the collective understanding of this complex phenomenon. The works of [38,75,76,77] collectively underscore the pivotal role of energy dissipation and strain energy parameters in delineating crack propagation dynamics. From [75]'s exploration of total plastic energy dissipation to [38]'s innovative introduction of a supplementary state variable for crack characterisation, and onto [76]'s focus on the cumulative change in cyclic strain energy and [77]'s holistic energy equation, these studies reflect a concerted effort to quantify and understand the energetic interactions at play in fatigue crack growth. This cohesive body of work not only reinforces the significance of energy considerations in crack growth studies but also illuminates the relation of various energy forms in influencing crack propagation under diverse loading conditions.

5. Hybrid Damage Theories or Parameter Definition

This section explores a range of hybrid damage theories and parameter definitions, emphasising the intricate interaction between microstructural components, energy dissipation, and crack growth dynamics in various materials. Each study, while distinct in its approach, contributes to a deeper understanding of fatigue damage, highlighting the importance of multi-dimensional analysis in fatigue crack growth models.
The following work [15] focuses on the high cycle fatigue (HCF) life in cast Al-Mg-Si alloys, which is particularly sensitive to the combination of microstructural components, inclusions and stress concentrators. Inclusions can range from large-scale shrinkage porosities to trapped oxides introduced during casting. When controlling for shrinkage porosity, the relevant microstructure at nucleation sites are often the larger Si particles within the eutectic regions. In this paper, an HCF model is introduced that recognises multiple scales of inclusions for crack formation. This ambitious work presents a model that addresses the role of constrained microplasticity around detached particles, shrinkage pores in the formation, and the microstructural growth of small fatigue cracks. It is based on cyclic crack tip displacement rather than linear elastic fracture toughness factor. The model is applied to an aluminium alloy A356-T6. A relationship between da/dN, Δσ, Su, A and ΔCTOD (Su is the maximum force) is presented, not seeming complicated on the surface, the authors try to make a collection of all possibilities and steps (eq. 1 in original publication). For example, the number of cycles is expressed as NT=NINC+NMSC+NPSC+NLC. Where INC is incubation, MSC propagation of a microstructurally small crack, PSC during the transition from MSC status to that of a dominant long crack and LC long crack. These cycle numbers are related to the porosity size using various coefficients, for the MSC/PSC regimes the cycle numbers are also expressed as a function of the most characteristic parameters of the process according to the authors. The whole process detailed so far is repeated for five types of inclusion models: (a) Distributed microporosity and Si particles; no significant pores or oxides (b) high levels of microporosity--shrinkage or gas pores with maximum diameter Dp < 3DCS but greater than the maximum Si particle diameter (c) large pores near the free surface (Dp > 3DCS) (d) large pores near the free surface (Dp > 3DCS) (e) large oxide films. The model results are correlated with the experimental A356-T6 plate results (in Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11 of the original publication). The work ends with some interesting graphs (Figure 13 and Figure 14 in the original publication) where the size of the inclusions is related to the crack growth mechanisms, as well as to the fatigue life of the materials. Transitioning from the microstructural emphasis in [15], the subsequent study by Khelil et al. [78] introduces an energetic perspective, illustrating how different approaches can converge to enhance the understanding of fatigue crack growth.
Khelil et al. [78] presented an energetic approach to Fatigue Crack Growth (FCG). This approach is based on the numerical determination of the plastic zone by introducing a new form of plastic radius. Experimental results on two aluminium alloys of the types 2024-T351 and 7075-T7351 were used to validate the developed numerical model. A good agreement between the two types of results has been found. Starting from eq. 21 where ΔWp is formulated as the plastic energy of the hysteresis cycle characteristic of cyclic loads and U as the Specific Energy. U, eq. 22, is expressed as a function of da/dN and subdivided into the three phases of the Kinetics Fatigue Failure Diagrams (KFFD). On the other hand, taking as valid the value of the angular strain γ=ε(3)0.5 ,near the crack tip one can express RN (ϴ) as the singular dominant term approximating the elastoplastic boundary. This allows further development of the ΔWp which in eq. 21 is left as an integral of stress and strain terms along rp (plastic radius), Spz (plasticised surface) and Q, eq. 23, (total energy dissipated).
W ( P Z ) = 2 1 N ' 1 N ' ' 1 + N ' ' σ 0 ϵ 0 K I σ 0 4 0 π 2 f N θ 2 d θ
U = M 2 B · d a / d N
Q = W ( P Z )   ·   B
These expressions are a function of parameters that must be calculated iteratively as indicated in the flow diagram (Figure 7 in the original publication). Figure 13 (a and c) shows the correlations of dissipated energy and ΔK and Figure 13 (b and d) does the same for da/dN and Q with acceptable correlations. Finally (in eqs. 35-37 in the original publication), the relationship of da/dN and ΔK for the 3 phases of the KFDD is established.
Building upon the energetic approach of [78], the next study [79] broadens the scope by integrating experimental studies and modelling techniques, thus providing a comprehensive view of the fatigue crack growth rates in aerospace alloys. The following work [79] apart from the crack growth model due to dissipated plastic energy, also presents an overview of experimental studies and modelling of fatigue crack growth rates in the aerospace titanium alloy Ti-6Al-4V. Firstly, numerous experimental tests were carried out on CT specimens of the alloy, from which the following was obtained: the FCGR under constant load. The influence of various parameters such as crack length and the effect of overload (retardation) were also evaluated experimentally. The crack opening was measured using DIC with a digital camera equipped with a Questar wide-range telescope, the set-up can record at 0.1 µm pixel quality. The Residual SIF (KR) was determined using X-ray diffraction techniques with a synchrotron, which also served to map all the deformation in the vicinity of the crack tip. The modelling techniques used for FCGR prediction are based on three concepts: (i) Crack advancement is controlled by the damage processes occurring within the fracture zone located ahead of the crack tip. This zone is embedded within the plastic zone of the crack tip, which in turn is surrounded by the zone of dominance of the elastic mechanics stress field solution (Kfield). (ii) Prediction of crack advance should be possible by knowing the deformation of a small volume around the crack tip, including residual stress and damage accumulation. (iii) Strictly speaking, a distinction must be made between plasticity and damage: although both are dissipative mechanisms, ultimately it is the damage component that determines failure. However, in many metals, the two parameters are related, i.e., it can be assumed that the damage at each point in the material is a function of the plastic deformation. Therefore, FCGR must be correlated with the plastic deformation processes at the crack tip. To have a good growth model it is essential to choose well the parameters, in this study it is first resorted to the crack tip blunting ρ , and the crack length a, obtaining the first relationship between both. The FCGRs is related to the residual intensity factor, the ratio of cyclic loads R and a series of adjustable constants. The adjustment and calculation of the parameters and constants leads to postulate the equation where the FCGRs is related to the dissipated plastic energy and a coefficient β, this equation implies the relation that a part of the dissipated plastic energy is converted into damage causing propagation. Finally, the equations ρ p (crack tip blunting) and ΔWp (plastic dissipate energy) are chosen as the constituents of the model whose results are compared with the experimental ones shown in Figure 14, offering a more than satisfactory correlation with better results of the blunting parameter.
Further exploring the realm of fatigue crack growth, the research by [80] delves into the complexities of plane strain cracks, challenging some prevailing notions and suggesting new directions for future modelling efforts. The following paper does not present a model per se, but rather evaluates some parameters as possible future keys in models and discards others. The paper [80] present the finite deformation elastoplastic analysis of plane strain cracks subjected to cyclic mode I loads with constant amplitude at various load ranges and ratios, as well as with overload and underload cycles. The Laird-Smith mechanism of crack growth by cyclic blunting and re-sharpening, which transfers material from the crack tip to the crack flanks, is visualised. In the present model, crack closure has never been detected. Furthermore, the assumed origin of PICC is ruled out. However, the simulated Fatigue Crack Growth (FCG) by blunting and regrinding reproduces the key experimental trends related to the effects of ΔK and single peaks of overload/underload. The calculated curves are kinked despite the absence of PICC, raising doubts about their reliability as a means of crack closure detection and assessment. Thus, the modelling performed manifests ambiguities with respect to PICC as the universal intrinsic factor capable of uniquely controlling FCG. On the other hand, in the absence of PICC, the calculated stress-strain responses near the tip, being the driving forces of fatigue damage and bond-breaking crack advance, manifest affinities with experimental FCG trends without the intervention of PICC. This implies both independent parameters of cyclic crack-tip loading under small-scale creep, such as the pair Kmax and Kmin, or Kmax and ΔK, or other equivalent, as well as the indispensable variables that drive FCG directly without PICC mediation. One of the major problems of looking at the influence of PICC on crack growth is that under plane strain conditions PICC is not observable. Therefore, numerical modelling appears to be a good tool. The finite element models use a perfect elastoplastic model with Von-Mises criterion associated with perfect plastic flow. No contact between the flanks was introduced, nor was meshing used. To maintain this criterion, the number of cycles was lowered, and different mesh models with very small element sizes were used. The most remarkable thing according to the authors is that the effect of PICC on FCG seems to be observable, even though PICC has not been modelled, so it is questionable whether PICC affects FCG or not.
Ultimately, the study in [81] evaluates the application of strip yields models, further enriching the discourse on fatigue crack growth under random loading, thereby showcasing the diversity and depth of hybrid damage theories. In this work [81] the application of strip yields models implemented in NASGRO (software) to estimate fatigue crack growth under random loading is evaluated. The two different strip yields model options (constant constraint-loss (CCL) and variable constraint-loss (VCL)) are evaluated and compared. The CCL is very similar to Newman's FASTRAN fatigue crack closure model. In this model, the tensile restraint factor α is constant along the plastic zone elements, but its value depends on the stress state, ranging from plane strain to plane stress. This constraint-loss is based on the observation that cracks that initially start with a flat face will eventually grow in a tilted-face mode. Newman proposed that the transition occurs when the size of the cyclic plastic zone (calculated from ΔKeff) reaches a percentage of the thickness. In the second, VCL, the tensile restraint factor α varies along the plastic zone according to a parabolic expression. The constraint decays from its value at the crack tip (αtip) to a flat stress value of 1.15 at the end in front of the plastic zone. The loss of restraint is also evaluated, but unlike the previous model, the value of αtip in both strain and plane strain is calculated from the ratio of the size of the plastic zone to the thickness of the specimen. The FCGR is formulated (eq. 24), where C, n, p and q are parameters that can be adjusted thanks to the NASMAT module. Three combinations of these parameters were considered, and the results compared with data from the NASGRO materials database and the literature. The ability of the models to estimate fatigue life and variability is analysed by comparing the simulated results with experimental fatigue crack growth data under different stationary Gaussian random loading processes on 2024-T351 aluminium CT specimens. The analysis showed that the two models with their three configurations provide good fatigue life predictions, with very similar results. The variability of the results due to the randomness of the loading was also analysed, and in this case, the CCL model provides a better estimate than the VCL. Finally, it should be noted that the best correlation with the experimental results was achieved with one of the combinations of constants proposed by the authors. This combination implemented with the NASMAT module improved any combination of the NASGRO database.
d a / d N = C ( 1 f 1 R K ) n 1 K t h K p 1 K m a x K c q
The models presented offer a variety of hybrid damage theories and parameter definitions, each contributing a unique lens through which the complex phenomenon of fatigue crack growth can be understood. From the detailed exploration of microstructural influences on high cycle fatigue in Al-Mg-Si alloys [15], to the energetic approaches of Khelil et al. [78] and the comprehensive experimental and modelling techniques applied to aerospace titanium alloys [79], these studies collectively underscore the multifaceted nature of fatigue damage. The critical evaluation of prevailing notions in [80], alongside the implementation of strip yield models in [81], further exemplifies the dynamic connection between theoretical modelling and practical applications. This diverse array of research not only enhances the knowledge in fatigue crack growth across different materials and conditions but also emphasises the necessity of integrating various theoretical and experimental approaches to develop more robust and comprehensive fatigue damage models.

6. Theory of Critical Distances (TCD)

This section delves into various applications and adaptations of the Theory of Critical Distances (TCD), a key concept in understanding fatigue crack growth. While each model presented here is based on TCD, they offer unique adaptations or applications of this theory, showcasing the versatility and broad applicability of TCD in different contexts.
In the following works [82,83] a numerical model based on the theory of critical distances (TCD) is proposed, this theory is not new [7,84] and basically consists of ensuring that crack propagation will occur when a certain physical quantity, at a given location in front of the crack front, reaches a specific value. In this work, this theory is partially modified to model the crack growth rate, the point method (PM) is adopted, and the physical quantity will be the Plastic Energy dissipation in front of the crack tip, at the critical point, similar to Figure 8. The numerical model is developed with Abaqus, under plane strain conditions due to the thin thickness of the specimens, a total of 7667 4-node elements with reduced integration were used and the minimum mesh size in the crack is less than 10 µm. Crack closure is modelled with the contact of the fronts, the linear penalty method coupled with an augmented Lagrangian algorithm was implemented to improve the accuracy of the overpressure during closure, also the frictionless sliding of the crack plane along the x-axis of contact was allowed. One of the critical points of this work, evidently, is the calculation of the plastic energy dissipated in each cycle for its parametric use as damage, and, therefore, causing the fatigue crack advance. For this purpose, an equation is proposed that basically integrates the stress over the plastic strain differential. To model the crack propagation, the mesh size must be considered in the order of 2.5-7.5 µm with the actual growth rate being around 0.1 µm /cycle, it quickly follows that many cycles must accumulate in the numerical model before the next node is released. During this accumulation of cycles, there is a stationary regime, which becomes transient when the energy accumulates at the critical distance and the crack propagates releasing the corresponding node. The modelling of crack propagation requires a specific post-processing program. Two are the key input parameters, the LCD, the Critical Distance Length, and the EC, the Critical Plastic Energy. The determination of these parameters is of vital importance. For this purpose, the da/dN vs ΔK curves are determined experimentally. Figure 15 plots the numerical results obtained for the da/dN vs ΔK curves for different values of LCD and EC. Finally, the parameters are fitted with the least-square method at LCD =17.5µm and EC =0.55 J/mm3. Good correlations between the model and the experimental results are observed.
Building on this foundation, subsequent studies further explore the intricacies of TCD, examining how factors like mesh size and microstructure size, along with variations in loading conditions, influence the outcomes predicted by the TCD-based models. An important step is to study how mesh size and microstructure size affect, the latter was done because the adjusted critical distance is of the order of the microstructure size, and several authors have previously focused on this detail [85,86]. Subsequently, in another paper [83] the same authors tested the goodness of the method against varying conditions of overload, underload, acceleration and deceleration of the loading cycles, as well as the moment (crack length) at which this is done. The model responds satisfactorily with respect to the experimental results obtained. The aim of this work [87] is to predict the fatigue life of 6201-T81 aluminium alloy wires containing geometric discontinuities using the TCD. The equivalent stress was evaluated by the PM and the Volume Method (VM), based on the maximum principal stress (σI). Regarding the latter method, the relationship between the critical distance and the number of cycles to failure was calibrated using two different methodologies. When life estimations were performed based on the first calibration, the results were considered unsatisfactory. However, when the second calibration method was used, almost all predictions fell within scatter bands with respect to the experimental data used to calibrate the model.
LCF prediction is critical to ensure the structural integrity of an engine. It is developed [88] a concept of strain energy gradient and a general procedure is established that combines the energy concept with the critical distance theory for life prediction of components with notches in LCF. There is a relationship of the strain energy distributed within the effective damage zone and fatigue crack growth. For specimens with GH4169 and TC4 notches, as well as a case study of a high-pressure turbine disc, the proposed procedure provides better life correlations than the Fatemi-Socie and Smith-Watson-Topper models.
As a continuous advance in TCD, this article [89] aims to compare the predictive capabilities of TCD based on a damage parameter, specifically four different ones were analysed. It also relates local strains and stresses to evaluate the fatigue life of notched components subjected to bending-torsion. The fatigue damage parameters tested are defined using well-known stress-based, strain-based, SWT-based, and energy-based relationships (Table 2). Multiaxial cyclic plasticity in the notch zone has been modelled with 3D Finite Elements, using three local stress-strain approaches: Neuber's rule, equivalent strain Energy Density rule (ESED) and the modified rule (ESEDM). Neuber's rule always led to more conservative results, and the ESEDM rule gave slightly better fatigue life predictions compared to the original ESED rule. As for the fatigue damage parameters, the energy-based models were more accurate.
In essence, the studies in this section demonstrate a comprehensive and multifaceted approach to applying the TCD, reflecting its adaptability and relevance across a spectrum of scenarios and materials. The diverse methodologies, ranging from the nuanced analysis of plastic energy dissipation to the integration of strain energy gradients and various damage parameters, collectively enrich the understanding of TCD and its practical applications in predicting fatigue life and crack growth.

7. Conclusions

In concluding this comprehensive review, it has been thoroughly explored the role of crack tip plasticity in shaping fatigue crack growth models. This paper seamlessly integrates numerical simulations and analytical methods, meticulously evaluating various approaches to emphasize their merits and constraints. The vital synergy of experimental data with theoretical and simulation-based models is illuminated, providing a detailed perspective on the progression of fatigue crack growth analysis.
In the realm of methodological advancements, it was observed a substantial shift towards integrating more sophisticated computational models, particularly Finite Element Modelling, with experimental approaches. This hybrid methodological approach has enhanced the precision of fatigue life predictions, allowing for more accurate simulations of real-world material behaviour under various stress conditions. The application of Finite Element Modelling has been particularly noteworthy in the context of complex loading conditions and the study of microscopic fracture mechanisms.
In these advances, it is of huge interest to pay attention to the work of de Matos and Nowell [44], who proposed modified models for diverse geometry and the possibility of incorporating residual stress fields. The importance lies in the ability to model fatigue problems with variables that in most cases cannot be neglected in real engineering applications. Models that are able to reproduce better real-signs are highly valuable. Due to the different approach taken by Mikheevskiy et al. [51], this research is worth to mention. New approaches and ideas might be the key ideas that are needed in the field to give real solution to a problem that, till today, is not yet fully understood.
For simplification purposes, summary tables (Table 3, Table 4 and Table 5) have been drawn up for the most used parameters in the literature. The information shown in the tables consists of the bibliographical reference with its first author and date of publication, as well as a brief description of the contribution with respect to the parameter in question. It also shows the methodology used in the research, or the methodologies in case there are several of them. The authors of this review subjectively suggest how positive (+), very positive (++) or extraordinarily positive (+++) the contribution of the research has been to advance in the knowledge and relationship that the parameter has with crack growth. This suggestion has an objective basis based on the number of citations of the paper, and another part of subjective evaluation in which the authors of this review have scored according to their personal experience. Table 3 condenses the information on the CTOD and the CTODP, this parameter has been questioned to replace even the CTOD, the advantages are varied, from the unity of the parameter, m, in the case of the CTOD, and above all the elastoplastic nature of the CTOD as opposed to the elastic approach of ΔK. At present, it is a well-established hypothesis that the origin and first stage of the crack has a plastic and microstructural basis, while the second growth phase is of an elastic nature and much more influential on a macro level. Any numerical or analytical approach needs experimental evidence of it. In this direction, works of Vaso-Olmo et al. [72] and Medhi et al. [73] give insights on the concept of CTOD and plasticity in the crack tip. Both works are able to explain the fatigue crack growth behaviour with the introduction of the crack tip plasticity in different ways. It is important to highlight latest Antunes et al. studies [64,65,67] for its exceptional contribution to advancing the understanding of the implications of crack tip plasticity in fatigue crack growth predictions based on numerical models with validation with experimental results. In summary, this research successfully established a robust method for the numerical calculation of CTOD and CTODP, emphasizing their dependency on key parameters. The study initially focuses on the 7050-T6 aluminium alloy, where cyclic plastic deformation is meticulously determined through experimental tests and subsequently modelled analytically. The development of a 3D numerical model to predict CTODP marks a significant advancement in this area. Experimental investigations provide a foundational understanding of the relationship between CTOD and FCGR. This is further complemented by numerical predictions of CTODP, which are tailored for varying crack lengths and da/dN values, thereby enhancing the precision and applicability of the model in predicting fatigue crack growth.
The study conducted by Pommier and Risbert [38] deserves recognition for its novel approach as it represents a significant stride in fatigue crack growth modelling by introducing plastic behaviour as an energy-based parameter evolving over time. Their novel approach, rooted in the thermodynamics of dissipative processes, uniquely incorporates a supplementary state variable for the crack. This allows for a continuous and detailed description of the crack state through complex loading sequences, a concept widely acknowledged and supported within the research community. The integration of plastic blunting, an established and well-regarded concept, into this framework further strengthens the model's applicability. By correlating the thermodynamic process of plastic energy dissipation with the hardening and sharpening of the crack tip growth, the model adeptly captures key aspects of fatigue crack progression, such as the Paris law and overload retardation effects. This work not only advances the understanding of fatigue crack dynamics but also sets a new precedent for incorporating time-dependent plastic behaviour into predictive modelling.
A significant portion of this review delved into the theoretical foundations underpinning the study of fatigue and fracture in materials. The TCD, a pivotal concept in understanding crack propagation and fatigue life, received considerable attention. Particularly Branco [89], in the context of assessing the fatigue life of notched components subjected to bending-torsion loading, a comparative analysis was conducted to evaluate various local approaches. It was observed that the implementation of ESEDM techniques in 3D finite element models resulted in a higher degree of precision in the prediction of fatigue life.
In Table 4, geometrical parameters that are used to measure some property or definition of the plastic zone are collected, suggestion column is also added here. Without obtaining a valid model for all situations, there are many investigations that manage to obtain very reliable correlations between the geometric parameters that define the plastic zone and crack growth.
Incorporating hybrid damage theories in this review is essential to explore alternative research avenues beyond the established ones. It is important to mention McDowell et al. [15] contribution about the sensitivity to the combination of microstructural components, inclusions and stress concentrators. Meanwhile Khelil et al. [78] introduces a energy approach through a numerical model, using of a new form of plastic radius that has been validate with experimental results. Korsunsky et al. [79] presents an energetic approach through the integration of modelling techniques and experimental studies in aerospace applications.
Table 5 presents a summary of the main lines of research that have in common some aspect of plastic energy as a driving parameter in crack growth. It is wise to mention the differences between the terms in Table 5. They can be divided into groups, terms related to Plastic Strain Energy, Elastic Strain Energy, Strain energy gradient, Total strain energy and Specific Energy. The terms related to Plastic Strain Energy are Plastic CTOD [64], plastic energy in the point used in PM [82,83], Kujawsky and Ellyin plastic energy [85,86], ΔWp [75,78,79] , Upl [77] and Q [78]. All of them consider a threshold of the plastic energy that triggers the crack growth when surpassed. Nevertheless, there are some differences between them; while [75,78,79] refer to the same term ΔWp related to load cycle, other terms like [85,86], Upl [77] and Q [78] are related to crack or fracture cycle release. Quan and Alderliesten [77] also refer to the variation in elastic strain energy in one entire crack propagation cycle, which is included inside the Elastic Strain Energy group. The Strain Energy Gradient is mentioned in [85,86,88] with the difference that Zhu et al. reflect stress and strain effects. Regarding the group where Total Strain Energy is mentioned, terms like ΔWt [89], UL, ΔCpσ and ΔCpε [76], and Ua [77] are included. While ΔWt [89], UL, and ΔCpσ [76] are stress-related parameters, ΔCpε [76] is strain-related, and Ua [77] is surface related. The last term is linked to Specific Energy, U [78], energy dissipated per unit volume during fatigue crack growth.
The existing problems are based on the absence of a physical basis for most of the proposed models. This makes models based on experimental results representative only if the operational conditions are the same. Something similar, but with a different nature, occurs in Finite Element Models, where boundary conditions, mesh size, analysable parameters, etc. are imposed. FEMs must be validated with experimental results and serve to extrapolate new test conditions; numerical simulations are much more efficient and economical than experimental tests. But one must not lose sight of the main objective of fatigue research, which is mathematical modelling based on physical phenomena that are common to fatigue cracks. As previously mentioned, in the problem of fracture and fatigue, the microstructural and macrostructural levels, as well as the plastic and/or elastic domain, intersect, representing a handicap for its physical and mathematical modelling.
The difficulty in this case lies in how to measure or estimate experimentally such specific aspects of the energy in such a small space as the plastic zone at the vertex of the crack. Also, analytically it presents an enormous challenge, because repeating concepts, applying mechanics of continuous media or particle systems, to date there has not been published (except for error or omission) an energetic study with Lagrangian or Hamiltonian mechanics, which were developed for an energetic approach to physics [90].

8. Future Directions

The field continues to advance through the interplay of theoretical development, methodological innovation, and practical application. Looking into the future, the ongoing challenges and emerging research directions promise to further deepen our understanding of material behaviour. The synergy between experimental research and computational modelling will be crucial in this research field, enabling researchers to predict, understand, and ultimately control material behaviour under the most demanding conditions.

Advancements in Multiscale Approaches

The development of universally applicable models for predicting fatigue life across diverse materials and loading conditions remains a primary challenge in material science. The various theories in the field of fracture and fatigue mechanics have generated a disparity that distances the link between the theory of continuous media (macroscopic approach) and the formulation of particle systems (microscopic approach), hindering the search for a point of convergence. In this sense, the need to explore a rapprochement between the two perspectives is evident, given the potential benefit that could be derived from the synergy between these two apparently disparate approaches.
It is interesting to note how, in a similar context, the expert Pradl achieved an important breakthrough in the field of fluid mechanics by bringing apparently divergent theories closer together.

Challenges in Experimental Measurement and Analytical Modelling

One of the major challenges in the field is the accurate measurement and estimation of energy within the plastic zone at the crack tip. The small scale of this area makes it difficult to assess, despite its critical role in understanding crack propagation. Future experimental efforts should not only focus on advanced microscopy and nano-indentation methods but also incorporate the use of high-energy X-rays, such as those available in synchrotron facilities. These high-energy X-rays can provide invaluable information about the bulk properties of materials, offering insights that are not attainable through conventional methods. This approach will significantly enhance our understanding of energy dynamics at the crack tip and facilitate the validation of computational models.
Incorporating energy-based concepts like Lagrangian or Hamiltonian mechanics into the study of fatigue and fracture in materials presents a promising yet largely unexplored avenue. These mechanics, known for their robust energy analysis framework, could transform our understanding of material behaviour under stress. Future research should focus on integrating these mechanics into material science models to deepen our theoretical understanding of energy dynamics within materials.

Energetic modelling procedure

A change in the numerical modelling of fracture and fatigue mechanics is proposed, moving away from conventional models based on Newton's theory and focusing on energy models. This new approach seeks not only to quantify the energy invested in the transformation of matter but also to integrate the energy dissipated during the test, considering phenomena such as friction and molecular friction.
This energy approach finds its roots in Griffith's theory, whose fundamental contribution was the introduction of the concept of the plastic radius, originally developed to explain the failure of brittle materials. It thus became a valuable theoretical framework for developing the relationship between theoretical and practical results.
The future directions in researching fatigue life prediction and fracture mechanics underscore the imperative to reduce the gap between macroscopic and microscopic approaches in material science. The challenge lies in developing universally applicable models that transcend diverse materials and loading conditions. Inspired by the successful integration of apparently divergent theories in fluid mechanics by expert Pradl, a similar rapprochement is sought between the theory of continuous media and particle systems in the fracture and fatigue domain. Overcoming challenges in experimental measurement and analytical modelling, particularly in accurately assessing energy within the plastic zone at the crack tip, necessitates advanced techniques such as high-energy X-rays available in synchrotron facilities. Future efforts should also explore the integration of energy-based concepts like Lagrangian or Hamiltonian mechanics, with their robust energy analysis framework, into material science models to revolutionize our understanding of material behaviour under stress. Additionally, a proposed shift in numerical modelling towards energy models, departing from conventional Newtonian-based models, holds promise in quantifying both the invested and dissipated energy during tests, with roots in Griffith's theory providing a foundational framework for this innovative approach. This multifaceted strategy aims to advance our theoretical understanding and practical applications in the field of fatigue and fracture mechanics.

Acknowledgements

Authors would like to acknowledge the following financial support: “FEDER Programa Operativo” by “Junta de Andalucía” (Spain) – grant reference UMA18-FEDERJA-250 and “IIPPITD” by “Universidad de Málaga” (Spain).

Nomenclature

A0 Parameter dependent of the stress intensity factor RICC Rugosity Induced Crack Closure
a0 Initial crack length Rint Interior radius for the area of interest
acorr Corrected crack length rm Monotonic plastic zone size
ACR Adjusted Compliance Ratio RN Dominant singular term approximation to the elastic-plastic boundary
B Number of additions per column Rout Outer radius for the area of interest
C Matrix function of the polar coordinates in system of equations rp Irwing plastic zone size for plane strain conditions
CCP Centre-Crack Plate specimen rpc Plastic zone size
COD Crack Open Displacement Rpr Reversed plastic zone size
CT Compact Tension specimen ry Plastic radius from Irwing
CTOD Crack Tip Open Displacement S Shift applied to the crack tip coordinates
CTODBS CTOD in BS7448-1:1991 SPZ Plastified surface
CTODJWES CTOD proposed by Japan Welding Engineering Society SWT Smith-Watson-Topper parameter
CTODp Plastic CTOD T-Stress Crack Tip Stress
CWI Crack Wake Influence U Normalized load ratio parameter / fatigue crack energy
d Matrix including εyy information in system of equations Ua Surface energy dissipated through new crack surface formation
D Dissipated energy UL Net section energy
da/dN Fatigue crack growth rate UPl Plastic energy dissipation
dW/dN Plastic work per cycle UTS Ultimate Tensile Strength
E Young's modulus VCWI E displacement at centreline due to concentrated force P on crack surface
Ec Critical plastic energy Vp Plastic component of crack mouth opening displacement
F Detection function for the symmetry axis W Specimen width
f Correction factors for plastic component of CTODJWES WOL Wedge Opening Loading specimen
IFR Influence Ratio X Vector of unknowns in system of equations
K Stress Intensity Factor (SIF) x, y Crack growing coordinate and crack opening coordinate
Kci,p Stress Intensity Factor when first contact between crack flanks occurs, plastic xct X coordinate of the crack-tip
Kcl,rl Stress Intensity Factor when first contact between crack flanks occurs, range-long xmax X coordinate of the maximum εyy value
Kcl,rs Stress Intensity Factor when first contact between crack flanks occurs, range-short xmin X coordinate of the minimum εyy value
KCWI Stress Intensity Factor due to concentrated force YR Yield-to-tensile Ratio, rys/ruts
Kexp Experimental estimation of the stress intensity factor δ SIF Relative error, % SIF
KF SIF opening mode ΔC Change in net-section strain energy for crack extension under plasticity in strain-controlled testing
KFFD Kinetics Fatigue Failure Diagrams ΔC Change in net-section strain energy for crack extension under plasticity in stress-controlled testing
Kfield SIF field Δij Kronecker delta
KI Opening mode SIF ΔJ J-integral
KIc Mode I fracture toughness ΔJeff Effective j-integral
KII SIF mode II ΔK SIF range
Kmax Maximum SIF ΔKeff Effective SIF range
Knom Nominally applied SIF δp Crack tip blunting
Kop SIF open ΔWe Elastic nominal strain energy density
KR SIF retardation ΔWe,eff Effective elastic nominal strain energy density
Kr Residual SIF ΔWpl Plastic nominal strain energy density
KS SIF Shear ΔWpl,eff Effective plastic nominal strain energy density
Kth Threshold Stress Intensity Factor ranges εxx , εyy Elastic strain in x and y directions
LCD Length Critical Distance η Statistical parameter
P Pearson's coefficient ν Poisson's ratio
PICC Plasticity-Induced Crack Closure ρ Plastic blunting
Q Total energy dissipated σapp Applied stress
R Stress ratio σmaxcomp Maximum stress at compression
r, θ Polar coordinates from the crack-tip σxx , σyy Stress in x and y directions
RA Area roughness of crack flanks σy Yield stress
rc Cyclic plastic zone size

References

  1. Bathias, C.; Baïlon, J.-P. La Fatigue Des Matériaux et Des Structures. Paris Montréal by Maloine Presses de l’Université de Montréal, 2224; ISBN 222400608X. [Google Scholar]
  2. Suresh, S. Suresh, S. (Subra) Fatigue of Materials. Cambridge University Press, 1998; ISBN 9780521570466. [Google Scholar]
  3. Stoychev, S.; Kujawski, D. Methods for Crack Opening Load and Crack Tip Shielding Determination: A Review. Fatigue Fract Eng Mater Struct 2003, 26, 1053–1067. [Google Scholar] [CrossRef]
  4. Berto, F.; Lazzarin, P. A Review of the Volume-Based Strain Energy Density Approach Applied to V-Notches and Welded Structures. Theoretical and Applied Fracture Mechanics 2009, 52, 183–194. [Google Scholar] [CrossRef]
  5. Karolczuk, A.; Macha, E. A Review of Critical Plane Orientations in Multiaxial Fatigue Failure Criteria of Metallic Materials. Int J Fract 2005, 134, 267–304. [Google Scholar] [CrossRef]
  6. Pippan, R.; Hohenwarter, A. Fatigue Crack Closure: A Review of the Physical Phenomena. Fatigue Fract Eng Mater Struct 2017, 40, 471–495. [Google Scholar] [CrossRef] [PubMed]
  7. Susmel, L. The Theory of Critical Distances: A Review of Its Applications in Fatigue. Eng Fract Mech 2008, 75, 1706–1724. [Google Scholar] [CrossRef]
  8. Hild, F.; Roux, S. Digital Image Correlation: From Displacement Measurement to Identification of Elastic Properties – a Review. Strain 2006, 42, 69–80. [Google Scholar] [CrossRef]
  9. Pan, B.; Qian, K.; Xie, H.; Asundi, A. Two-Dimensional Digital Image Correlation for in-Plane Displacement and Strain Measurement: A Review. Meas Sci Technol 2009, 20, 62001. [Google Scholar] [CrossRef]
  10. Branco, R.; Antunes, F. V; Costa, J.D. A Review on 3D-FE Adaptive Remeshing Techniques for Crack Growth Modelling. Eng Fract Mech 2015, 41, 170–195. [Google Scholar] [CrossRef]
  11. Rege, K.; Lemu, H.G. A Review of Fatigue Crack Propagation Modelling Techniques Using FEM and XFEM. IOP Conf Ser Mater Sci Eng 2017, 276, 012027. [Google Scholar] [CrossRef]
  12. Vikram, N.; Kumar, R. Review on Fatigue-Crack Growth and Finite Element Method. Int J Sci Eng Res 2013, 4, 833–843. [Google Scholar]
  13. Zhu, S.-P.; Ye, W.-L.; Correia, J.A.F.O.; Jesus, A.M.P.; Wang, Q. Stress Gradient Effect in Metal Fatigue: Review and Solutions. Theoretical and Applied Fracture Mechanics 2022, 121, 103513. [Google Scholar] [CrossRef]
  14. Paul, S.K.; Sivaprasad, S.; Dhar, S.; Tarafder, S. Key Issues in Cyclic Plastic Deformation: Experimentation. Mechanics of Materials 2011, 43, 705–720. [Google Scholar] [CrossRef]
  15. McDowell, D.L.; Gall, K.; Horstemeyer, M.F.; Fan, J. Microstructure-Based Fatigue Modeling of Cast A356-T6 Alloy. Eng Fract Mech 2003, 70, 49–80. [Google Scholar] [CrossRef]
  16. Campbell, J. Invisible Macrodefects in Castings. Journal de Physique IV Proceedings 1993, 03, 861–872. [Google Scholar] [CrossRef]
  17. Sunder, R.; Porter, W.J.; Ashbaugh, N.E. Fatigue Voids and Their Significance. Fatigue Fract Eng Mater Struct 2002, 25, 1015–1024. [Google Scholar] [CrossRef]
  18. Elber, W. Fatigue Crack Closure under Cyclic Tension. Eng Fract Mech 1970, 2, 37–45. [Google Scholar]
  19. Blom, A.F.; Holm, D.K. An Experimental and Numerical Study of Crack Closure. Eng Fract Mech 1985, 22, 997–1011. [Google Scholar] [CrossRef]
  20. Donald, K.; Paris, P.C. An Evaluation of DKeff Estimation Procedure on 6061-T6 and 2024-T3 Aluminum Alloys. Journal of Fatigue 1999, 21, S47–S57. [Google Scholar] [CrossRef]
  21. Rao, K.T.V.; Yu, W.; Ritchie, R.O. On the behaviour of Small Fatigue Cracks in Commercial Aluminium Lithium Alloys. Eng Fract Mech 1988, 31, 623–635. [Google Scholar] [CrossRef]
  22. Costa, J.D.M.; Ferreira, J.A.M. Efect of Stress Ratio and Specimen Thickness on Fatigue Crack Growth of CK45 Steel. Theoretical and applied fracture mechanics 1998, 30, 65–73. [Google Scholar] [CrossRef]
  23. Budianske, B.; Hutchinson, W.J. Analysis of Closure in Fatigue Crack Growth. journal of applied mechacins 1978, 45, 267–277. [Google Scholar] [CrossRef]
  24. Kujawski, D. Enhanced Model of Partial Crack Closure for Correlation of R-Ratio Effects in Aluminum Alloys. Int J Fatigue 2001, 23, 95–102. [Google Scholar] [CrossRef]
  25. Christopher, C.J.; James, M.N.; Patterson, E.A.; Tee, K.F. Towards a New Model of Crack Tip Stress Fields. Int J Fract 2007, 148, 361–371. [Google Scholar] [CrossRef]
  26. Antunes, F.V.; Sousa, T.; Branco, R.; Correia, L. Effect of Crack Closure on Non-Linear Crack Tip Parameters. Int J Fatigue 2015, 71, 53–63. [Google Scholar] [CrossRef]
  27. Nowell, D.; de Matos, P.F.P. Application of Digital Image Correlation to the Investigation of Crack Closure Following Overloads. Procedia Eng 2010, 2, 1035–1043. [Google Scholar] [CrossRef]
  28. Yusof, F.; Lopez-Crespo, P.; Withers, P.J. Effect of Overload on Crack Closure in Thick and Thin Specimens via Digital Image Correlation. Int J Fatigue 2013, 56, 17–24. [Google Scholar] [CrossRef]
  29. Lopez-Crespo, P.; Shterenlikht, A.; Patterson, E.A.; Yates, J.R.; Withers, P.J. The Stress Intensity of Mixed Mode Cracks Determined by Digital Image Correlation. J Strain Anal Eng Des 2008, 43, 769–780. [Google Scholar] [CrossRef]
  30. Shterenlikht, A.; Díaz Garrido, F.A.; Lopez-Crespo, P.; Withers, P.J.; Patterson, E.A. Mixed Mode (KI+KII) Stress Intensity Factor Measurement by Electronic Speckle Pattern Interferometry and Image Correlation. Applied Mechanics and Materials 2004, 1–2, 107–112. [Google Scholar] [CrossRef]
  31. Lopez-Crespo, P.; Burguete, R.L.; Patterson, E.A.; Shterenlikht, A.; Withers, P.J.; Yates, J.R. Study of a Crack at a Fastener Hole by Digital Image Correlation. Exp Mech 2009, 49, 551–559. [Google Scholar] [CrossRef]
  32. Zanganeh, M.; Lopez-Crespo, P.; Tai, Y.H.; Yates, J.R. Locating the Crack Tip Using Displacement Field Data: A Comparative Study. Strain 2013, 49, 102–115. [Google Scholar] [CrossRef]
  33. Chen, J.; Zhan, N.; Zhang, X.; Wang, J. Improved Extended Digital Image Correlation for Crack Tip Deformation Measurement. Opt Lasers Eng 2015, 65, 103–109. [Google Scholar] [CrossRef]
  34. Pippan, R.; Grosinger, W. Fatigue Crack Closure: From LCF to Small Scale Yielding. Int J Fatigue 2013, 46, 41–48. [Google Scholar] [CrossRef]
  35. Kawabata, T.; Tagawa, T.; Sakimoto, T.; Kayamori, Y.; Ohata, M.; Yamashita, Y.; Tamura, E.; Yoshinari, H.; Aihara, S.; Minami, F.; et al. Proposal for a New CTOD Calculation Formula. Eng Fract Mech 2016, 159, 16–34. [Google Scholar] [CrossRef]
  36. Tvergaard, V. On Fatigue Crack Growth in Ductile Materials by Crack–Tip Blunting. J Mech Phys Solids 2004, 52, 2149–2166. [Google Scholar] [CrossRef]
  37. Pelloux, R.M.N. Crack Extension by Alternating Shear. Eng Fract Mech 1970, 1, 697–700. [Google Scholar] [CrossRef]
  38. Pommier, S.; Risbet, M. Time Derivative Equations for Mode I Fatigue Crack Growth in Metals. Int J Fatigue 2005, 27, 1297–1306. [Google Scholar] [CrossRef]
  39. Noroozi, A.; Glinka, G.; Lambert, S. A Two Parameter Driving Force for Fatigue Crack Growth Analysis. Int J Fatigue 2005, 27, 1277–1296. [Google Scholar] [CrossRef]
  40. Ould Chikh, B.; Imad, A.; Benguediab, M. Influence of the Cyclic Plastic Zone Size on the Propagation of the Fatigue Crack in Case of 12NC6 Steel. Comput Mater Sci 2008, 43, 1010–1017. [Google Scholar] [CrossRef]
  41. Cruces, A.S.; Mokhtarishirazabad, M.; Moreno, B.; Zanganeh, M.; Lopez-Crespo, P. Study of the Biaxial Fatigue Behaviour and Overloads on S355 Low Carbon Steel. Int J Fatigue 2020, 134, 105466. [Google Scholar] [CrossRef]
  42. Chen, H.; Chen, W.; Li, T.; Ure, J. Effect of Circular Holes on the Ratchet Limit and Crack Tip Plastic Strain Range in a Centre Cracked Plate. Eng Fract Mech 2011, 78, 2310–2324. [Google Scholar] [CrossRef]
  43. Shih, C.F. Relationships between the J-Integral and the Crack Opening Displacement for Stationary and Extending Cracks. J Mech Phys Solids 1981, 29, 305–326. [Google Scholar] [CrossRef]
  44. De-Matos, P.F.P.; Nowell, D. Analytical and Numerical Modelling of Plasticity-Induced Crack Closure in Cold-Expanded Holes. Fatigue Fract Eng Mater Struct 2008, 31, 488–503. [Google Scholar] [CrossRef]
  45. Cornet, C.; Zhao, L.G.; Tong, J. Ratchetting Strain as a Damage Parameter in Controlling Crack Growth at Elevated Temperature. Eng Fract Mech 2009, 76, 2538–2553. [Google Scholar] [CrossRef]
  46. Zhao, L.; Tong, J. A Viscoplastic Study of Crack-Tip Deformation and Crack Growth in a Nickel-Based Superalloy at Elevated Temperature. J Mech Phys Solids 2008, 56, 3363–3378. [Google Scholar] [CrossRef]
  47. Chaboche, J.L. A Review of Some Plasticity and Viscoplasticity Constitutive Theories. Int J Plast 2008, 24, 1642–1693. [Google Scholar] [CrossRef]
  48. Tong, J.; Zhao, L.G.; Lin, B. Ratchetting Strain as a Driving Force for Fatigue Crack Growth. Int J Fatigue 2013, 46, 49–57. [Google Scholar] [CrossRef]
  49. Lin, B.; Zhao, L.G.; Tong, J. A Crystal Plasticity Study of Cyclic Constitutive Behaviour, Crack-Tip Deformation and Crack-Growth Path for a Polycrystalline Nickel-Based Superalloy. Eng Fract Mech 2011, 78, 2174–2192. [Google Scholar] [CrossRef]
  50. Lu, Y.W.; Lupton, C.; Zhu, M.L.; Tong, J. In Situ Experimental Study of Near-Tip Strain Evolution of Fatigue Cracks. Exp Mech 2015, 55, 1175–1185. [Google Scholar] [CrossRef]
  51. Mikheevskiy, S.; Bogdanov, S.; Glinka, G. Analysis of Fatigue Crack Growth under Spectrum Loading – The UniGrow Fatigue Crack Growth Model. Theoretical and Applied Fracture Mechanics 2015, 79, 25–33. [Google Scholar] [CrossRef]
  52. Borges, M.F.; Neto, D.M.; Antunes, F. V. Numerical Simulation of Fatigue Crack Growth Based on Accumulated Plastic Strain. Theoretical and Applied Fracture Mechanics 2020, 108, 102676. [Google Scholar] [CrossRef]
  53. Miarka, P.; Cruces, A.S.; Seitl, S.; Malíková, L.; Lopez-Crespo, P. Influence of the Constraint Effect on the Fatigue Crack Growth Rate in S355 J2 Steel Using Digital Image Correlation. Fatigue Fract Eng Mater Struct 2020, 43, 1703–1718. [Google Scholar] [CrossRef]
  54. Park, H.-B.; Kim, K.-M.; Lee, B.-W. Plastic Zone Size in Fatigue Cracking. International Journal of Pressure Vessels and Piping 1996, 68, 279–285. [Google Scholar] [CrossRef]
  55. Donald, K.; Paris, P. An Evaluation of ΔKeff Estimation Procedures on 6061-T6 and 2024-T3 Aluminum Alloys*1. Int J Fatigue 1999, 21, S47–S57. [Google Scholar] [CrossRef]
  56. Paris, P.; Tada, H.; Donald, K. Service Load Fatigue Damage - a Historical Perspective. Int J Fatigue 1999, 21, S35–S46. [Google Scholar] [CrossRef]
  57. Elber, W. The Significance of Fatigue Crack Closure. Damage tolerance in aircraft structures. ASTM STP 486 1971, 230. [Google Scholar]
  58. Zhang, J.; He, X.; Du, S. Analyses of the Fatigue Crack Propagation Process and Stress Ratio Effects Using the Two Parameter Method. Int J Fatigue 2005, 27, 1314–1318. [Google Scholar] [CrossRef]
  59. Lopez-Crespo, P.; Shterenlikht, A.; Yates, J.R.; Patterson, E.A.; Withers, P.J. Some Experimental Observations on Crack Closure and Crack-Tip Plasticity. Fatigue Fract Eng Mater Struct 2009, 32, 418–429. [Google Scholar] [CrossRef]
  60. Zhang, J.; He, X.D.; Sha, Y.; Du, S.Y. The Compressive Stress Effect on Fatigue Crack Growth under Tension–Compression Loading. Int J Fatigue 2010, 32, 361–367. [Google Scholar] [CrossRef]
  61. Silva, F.S. The Importance of Compressive Stresses on Fatigue Crack Propagation Rate. Int J Fatigue 2005, 27, 1441–1452. [Google Scholar] [CrossRef]
  62. Silva, F.S. Crack Closure Inadequacy at Negative Stress Ratios. Int J Fatigue 2004, 26, 241–252. [Google Scholar] [CrossRef]
  63. Pokluda, J. Dislocation-Based Model of Plasticity and Roughness-Induced Crack Closure. Int J Fatigue 2013, 46, 35–40. [Google Scholar] [CrossRef]
  64. Antunes, F.V.; Díaz, F.A.; Vasco-Olmo, J.M.; Prates, P. Numerical Determination of Plastic CTOD. Fatigue Fract Eng Mater Struct 2018, 41, 2197–2207. [Google Scholar] [CrossRef]
  65. Antunes, F. V.; Branco, R.; Prates, P.A.; Borrego, L. Fatigue Crack Growth Modelling Based on CTOD for the 7050-T6 Alloy. Fatigue Fract Eng Mater Struct 2017, 40, 1309–1320. [Google Scholar] [CrossRef]
  66. Voce, F. The Relationship between Stress and Strain for Homogeneous Deformation. Journal of the Institue of Metals 1948, 74, 537–562. [Google Scholar]
  67. Antunes, F.V.; Serrano, S.; Branco, R.; Prates, P. Fatigue Crack Growth in the 2050-T8 Aluminium Alloy. Int J Fatigue 2018, 115, 79–88. [Google Scholar] [CrossRef]
  68. Tagawa, T.; Kawabata, T.; Sakimoto, T.; Kayamori, Y.; Ohata, M.; Yamashita, Y.; Tamura, E.; Yoshinari, H.; Aihara, S.; Minami, F.; et al. Experimental Measurements of Deformed Crack Tips in Different Yield-to-Tensile Ratio Steels. Eng Fract Mech 2014, 128, 157–170. [Google Scholar] [CrossRef]
  69. Tagawa, T.; Kayamori, Y.; Ohata, M.; Handa, T.; Kawabata, T.; Yamashita, Y.; Tsutsumi, K.; Yoshinari, H.; Aihara, S.; Hagihara, Y. Comparison of CTOD Standards: BS 7448-Part 1 and Revised ASTM E1290. Eng Fract Mech 2010, 77, 327–336. [Google Scholar] [CrossRef]
  70. Kayamori, Y.; Inoue, T.; Tagawa, T. Transformation of BS7448-CTOD to ASTM E1290-CTOD. J Press Vessel Technol 2010, 132, 041401. [Google Scholar] [CrossRef]
  71. Kayamori, Y.; Inque, T.; Tagawa, T. Changes in ISO 15653-Based CTOD for Specimens of A0/W=0.45*. Journal of Solid Mechanics and Materials Engineering 2012, 6, 645–653. [Google Scholar] [CrossRef]
  72. Vasco-Olmo, J.M.; Diaz, F.A.; Antunes, F.V.; James, M.N. Experimental Evaluation of CTOD in Constant Amplitude Fatigue Crack Growth from Crack Tip Displacement Fields. Frattura ed Integrità Strutturale 2017, 11, 157–165. [Google Scholar] [CrossRef]
  73. Mokhtarishirazabad, M.; Lopez-Crespo, P.; Zanganeh, M. Stress Intensity Factor Monitoring under Cyclic Loading by Digital Image Correlation. Fatigue Fract Eng Mater Struct 2018, 41, 2162–2171. [Google Scholar] [CrossRef]
  74. Yates, J.R.; Zanganeh, M.; Tai, Y.H. Quantifying Crack Tip Displacement Fields with DIC. Eng Fract Mech 2010, 77, 2063–2076. [Google Scholar] [CrossRef]
  75. Klingbeil, N. A Total Dissipated Energy Theory of Fatigue Crack Growth in Ductile Solids. Int J Fatigue 2003, 25, 117–128. [Google Scholar] [CrossRef]
  76. Ravi Chandran, K.S. Mechanics of Fatigue Crack Growth under Large-Scale Plasticity: A Direct Physical Approach for Single-Valued Correlation of Fatigue Crack Growth Data. Int J Fatigue 2018, 117, 299–313. [Google Scholar] [CrossRef]
  77. Quan, H.; Alderliesten, R.C. The Energy Dissipation during Fatigue Crack Growth in Metallic Materials. Eng Fract Mech 2022, 269, 108567. [Google Scholar] [CrossRef]
  78. Khelil, F.; Aour, B.; Belhouari, M.; Benseddiq, N. Modeling of Fatigue Crack Propagation in Aluminum Alloys Using an Energy Based Approach. Engineering, Technology & Applied Science Research 2013, 3, 488–496. [Google Scholar] [CrossRef]
  79. Korsunsky, A.M.; Song, X.; Belnoue, J.; Jun, T.; Hofmann, F.; De Matos, P.F.P.; Nowell, D.; Dini, D.; Aparicio-Blanco, O.; Walsh, M.J. Crack Tip Deformation Fields and Fatigue Crack Growth Rates in Ti–6Al–4V☆. Int J Fatigue 2009, 31, 1771–1779. [Google Scholar] [CrossRef]
  80. Toribio, J.; Kharin, V. Simulations of Fatigue Crack Growth by Blunting–Re-Sharpening: Plasticity Induced Crack Closure vs. Alternative Controlling Variables. Int J Fatigue 2013, 50, 72–82. [Google Scholar] [CrossRef]
  81. Moreno, B.; Martin, A.; Lopez-Crespo, P.; Zapatero, J.; Domínguez, J. Estimations of Fatigue Life and Variability under Random Loading in Aluminum Al-2024T351 Using Strip Yield Models from NASGRO. Int J Fatigue 2016, 91, 414–422. [Google Scholar] [CrossRef]
  82. Zheng, X.; Cui, H.; Su, X.; Engler-Pinto, C.C.; Wen, W. Numerical Modeling of Fatigue Crack Propagation Based on the Theory of Critical Distances. Eng Fract Mech 2013, 114, 151–165. [Google Scholar] [CrossRef]
  83. Zheng, X.; Cui, H.; Engler-Pinto, C.C.; Su, X.; Wen, W. Numerical Modeling of Fatigue Crack Propagation Based on the Theory of Critical Distances: Effects of Overloads and Underloads. Eng Fract Mech 2014, 128, 91–102. [Google Scholar] [CrossRef]
  84. Taylor, D. The Theory of Critical Distances. Eng Fract Mech 2008, 75, 1696–1705. [Google Scholar] [CrossRef]
  85. Kujawski, D.; Ellyin, F. A Fatigue Crack Propagation Model. Eng Fract Mech 1984, 20, 695–704. [Google Scholar] [CrossRef]
  86. Chalant, G.; Remy, L. Model of Fatigue Crack Propagation by Damage Accumulation at the Crack Tip. Eng Fract Mech 1983, 18, 939–952. [Google Scholar] [CrossRef]
  87. Adriano, V.S.R.; Martínez, J.M.G.; Ferreira, J.L.A.; Araújo, J.A.; da Silva, C.R.M. The Influence of the Fatigue Process Zone Size on Fatigue Life Estimations Performed on Aluminum Wires Containing Geometric Discontinuities Using the Theory of Critical Distances. Theoretical and Applied Fracture Mechanics 2018, 97, 265–278. [Google Scholar] [CrossRef]
  88. Zhu, S.P.; Liu, Y.; Liu, Q.; Yu, Z.Y. Strain Energy Gradient-Based LCF Life Prediction of Turbine Discs Using Critical Distance Concept. Int J Fatigue 2018, 113, 33–42. [Google Scholar] [CrossRef]
  89. Branco, R.; Costa, J.D.; Borrego, L.P.; Berto, F.; Razavi, S.M.J.; Macek, W. Comparison of Different One-Parameter Damage Laws and Local Stress-Strain Approaches in Multiaxial Fatigue Life Assessment of Notched Components. Int J Fatigue 2021, 151, 106405. [Google Scholar] [CrossRef]
  90. Landau, L.D.; Lifshitz, E.M.; Berestetskii, V.B.; Pitaevskii, L.P. Física Teórica. Mecánica: 1; Editorial Reverté, Ed 2005; ISBN 978-84-291-4081-1.
Figure 1. Schematic illustration of the procedure to determine ΔJ (a) and the ΔJeff (b). From [34] (c) Stress-strain curve of a low cycle fatigue experiment at a plastic strain amplitude of about 0.3%. The SEM micrographs show the closing of the crack. From [34] .
Figure 1. Schematic illustration of the procedure to determine ΔJ (a) and the ΔJeff (b). From [34] (c) Stress-strain curve of a low cycle fatigue experiment at a plastic strain amplitude of about 0.3%. The SEM micrographs show the closing of the crack. From [34] .
Preprints 91137 g001
Figure 2. Fatigue crack growth data for Al 2024 T351 aluminium alloy as a function of the two-parameter driving force, ΔK. From [39].
Figure 2. Fatigue crack growth data for Al 2024 T351 aluminium alloy as a function of the two-parameter driving force, ΔK. From [39].
Preprints 91137 g002
Figure 3. Opening stresses for an infinite plate with a circular hole and two radial symmetric cracks propagating under constant amplitude loading with and without residual stresses due to cold expansion (ur/r = 4%). (a) and (b) are for σ0max=0.4 and 0.6, respectively, with R=0. From [44].
Figure 3. Opening stresses for an infinite plate with a circular hole and two radial symmetric cracks propagating under constant amplitude loading with and without residual stresses due to cold expansion (ur/r = 4%). (a) and (b) are for σ0max=0.4 and 0.6, respectively, with R=0. From [44].
Preprints 91137 g003
Figure 4. (a, b) Crack growth per cycle (da/dN) versus plastic CTOD range (ΔCTODp) for both tests. (c) Comparison of FCGR specimens with relation to cyclic plastic zone Plane Strain conditions. (d) Plane Stress and Strain. From [54].
Figure 4. (a, b) Crack growth per cycle (da/dN) versus plastic CTOD range (ΔCTODp) for both tests. (c) Comparison of FCGR specimens with relation to cyclic plastic zone Plane Strain conditions. (d) Plane Stress and Strain. From [54].
Preprints 91137 g004
Figure 5. FCGR response of alloys (a) 6061-T6 and (b) 2024-T3 showing R=0.7 data and constant Kmax data (Kmax =22 MPa m ) as a function of Kapp. From [55] .
Figure 5. FCGR response of alloys (a) 6061-T6 and (b) 2024-T3 showing R=0.7 data and constant Kmax data (Kmax =22 MPa m ) as a function of Kapp. From [55] .
Preprints 91137 g005
Figure 6. The comparison between the predicted Kmax vs ΔK curve and the experimental obtained data for (a) 2324 aluminium alloy and (b) 6013 aluminium alloy. From [58].
Figure 6. The comparison between the predicted Kmax vs ΔK curve and the experimental obtained data for (a) 2324 aluminium alloy and (b) 6013 aluminium alloy. From [58].
Preprints 91137 g006
Figure 7. (a) Opening (KI ) and (b) shearing (KII ) mode stress intensity factors measured during the loading portion of a fatigue cycle (R = 0) with θ = 30º on a 7010 aluminium alloy centre-cracked plate. The nominal values were calculated neglecting closure effects. The experimental values were determined relative to an image captured at zero load. From [59].
Figure 7. (a) Opening (KI ) and (b) shearing (KII ) mode stress intensity factors measured during the loading portion of a fatigue cycle (R = 0) with θ = 30º on a 7010 aluminium alloy centre-cracked plate. The nominal values were calculated neglecting closure effects. The experimental values were determined relative to an image captured at zero load. From [59].
Preprints 91137 g007
Figure 8. Interesting crack tip region/point, assuming that damage is a local issue Theory of Critical Distance (TCD). Region of reversed plastic zone rpc: if fatigue damage accumulates. Region of forward plastic zone rpm.
Figure 8. Interesting crack tip region/point, assuming that damage is a local issue Theory of Critical Distance (TCD). Region of reversed plastic zone rpc: if fatigue damage accumulates. Region of forward plastic zone rpm.
Preprints 91137 g008
Figure 9. Correlation between CTOD calculated from P-Vg curve and actual CTOD at mid-thickness experimentally measured by silicone rubber casting for specimen with thickness B = 30 mm. (a) CTODJWES. (b) CTODBS. From [35]. .
Figure 9. Correlation between CTOD calculated from P-Vg curve and actual CTOD at mid-thickness experimentally measured by silicone rubber casting for specimen with thickness B = 30 mm. (a) CTODJWES. (b) CTODBS. From [35]. .
Preprints 91137 g009
Figure 10. Crack growth per cycle (da/dN) versus CTOD range (ΔCTODp) for both tests).
Figure 10. Crack growth per cycle (da/dN) versus CTOD range (ΔCTODp) for both tests).
Preprints 91137 g010
Figure 11. Evolution of fatigue crack propagation rate estimated by numerical simulation. (a) Master plot for R=0 and (b) Plastic dissipation per cycle. From [75].
Figure 11. Evolution of fatigue crack propagation rate estimated by numerical simulation. (a) Master plot for R=0 and (b) Plastic dissipation per cycle. From [75].
Preprints 91137 g011
Figure 12. Fatigue crack growth data under stress control for cracks growing from starter holes in round tension specimens of 0.45% C steel, plotted in terms of (a) ΔK and (b) the change in net-section strain energy, ΔC. From [76].
Figure 12. Fatigue crack growth data under stress control for cracks growing from starter holes in round tension specimens of 0.45% C steel, plotted in terms of (a) ΔK and (b) the change in net-section strain energy, ΔC. From [76].
Preprints 91137 g012
Figure 13. Comparison of measured and estimated dissipated energy per cycle for (a) 2024-T351 and (b) 7075-T7351. Comparison of measured Evolution of da/dN with Q. (c)for 2024-T351 (d) for 7075-T7351. From [78].
Figure 13. Comparison of measured and estimated dissipated energy per cycle for (a) 2024-T351 and (b) 7075-T7351. Comparison of measured Evolution of da/dN with Q. (c)for 2024-T351 (d) for 7075-T7351. From [78].
Preprints 91137 g013
Figure 14. (a)The dependence of crack blunting parameter increment (per cycle) on the applied stress intensity factor, DK, that is described well by a power law relationship. (b) The dependence of crack tip plastic energy dissipation (per cycle) on the applied stress intensity factor, ΔK, that is described well by a power law relationship. From [79].
Figure 14. (a)The dependence of crack blunting parameter increment (per cycle) on the applied stress intensity factor, DK, that is described well by a power law relationship. (b) The dependence of crack tip plastic energy dissipation (per cycle) on the applied stress intensity factor, ΔK, that is described well by a power law relationship. From [79].
Preprints 91137 g014
Figure 15. Evolution of fatigue crack propagation rate estimated by numerical simulation. From [82,83].
Figure 15. Evolution of fatigue crack propagation rate estimated by numerical simulation. From [82,83].
Preprints 91137 g015
Table 2. Damage parameters.
Table 2. Damage parameters.
Damage parameter Equation
Stress σa = (σ'fm) (2Nf) b
Strain εa = (σ'f - σm) (1/E) (2N) fb +ε'f (2Nf) c
Energy ΔWt = κt (2Nf) αt +ΔW0t
SWT εa σmax = (σ'f )2 (1/E) (2N) f2b + ε'f (2Nf) c+b
Table 3. CTOD and/or CTODP as a critical parameter used as driving force in the crack grown modelling. CTOD is Crack Tip Open Displacement, where P denotes Plastic.
Table 3. CTOD and/or CTODP as a critical parameter used as driving force in the crack grown modelling. CTOD is Crack Tip Open Displacement, where P denotes Plastic.
Ref Authors and Date Description/Main Contribution Methods Sug.
[26] Antunes et al. 2014 Analysis of remote compliance is the best numerical parameter to quantify the crack opening level Numerical ++
[26] Antunes et al. 2015 Establish an analytical relation between CTOD and da/dN. This relation was tested numerically Numerical ++
[35] Kawabata et al. 2016 A new CTOD method is investigated considering the variation of crack tip blunting (strain hardening). The calculation formula is based on three-dimensional elasto-plastic FEM Numerical-Experimental ++
[43] Shih 1986 Establish the relation between the J-integral and the crack opening displacement by exploiting the dominance of the Hutchinson--Rice-Rosengren singularity in the crack-tip region. Numerical-Experimental ++
[63] Pokluda 2011 A discrete dislocation model of contact effects in small-scale yielding is presented. The model enables to directly assess the magnitude of both plasticity and roughness-induced components of crack closure. Analytical ++
[64] Antunes et al. 2018 To establish a method of numerical calculation of the CTOD and CTODP and the dependence on certain parameters Numerical +++
[65] Antunes et al. 2017 The 7050-T6 aluminium alloy cyclic plastic deformation was determined experimentally and modelled analytically. A 3D numerical model was developed to predict the CTODP Numerical-Experimental-Analytical +++
[67] Antunes et al. 2018 First, experimental tests were conducted to obtain the relation between CTOD and FCGR. Then, numerical predictions of CTODP were obtained for different crack length and da/dN. Numerical-Experimental +++
[68] Tagawa et al. 2014 Numerical and experimental methods to determine a method to calculate CTOD and CTODP Numerical-Experimental +
[69] Tagawa et al. 2009 The CTOD testing methodologies effects on CTOD values were investigated according to tests conducted by the Japan Welding Engineering Society (WES) Numerical-Experimental ++
[70] Kayamori et al. 2010 Experimental investigations and analytical developments into crack tip opening displacement (CTOD) conducted to stablish the relationship between BS7448-CTOD and ASTM E1290-CTOD. Numerical-Experimental +
[71] Kayamori et al. 2012 Two new CTOD calculations were proposed, for deep-notched specimens’ displacement-conversion CTOD, and for shallow-notched specimens, a J-conversion CTOD was proposed. Numerical-Experimental ++
[72] Vasco-Olmo et al. 2017 A methodology is developed to measure and analyse the CTOD and CTODP from experimental data. Experimental +++
[74] Yates et al. 2010 The paper gives an overview of some DIC applications for crack tip characterisation such as CTOD and CTODP measures as well as data obtained. Analytical-Experimental ++
Table 4. Plastic Size as a critical parameter used as driving force in the crack grown modelisation.
Table 4. Plastic Size as a critical parameter used as driving force in the crack grown modelisation.
Ref SpecifiPara. Authors and date Description/Main Contribution Methods Sug.
[26] Rpr Antunes et al. 2015 Reverse Plastic zone Size. the crack closure phenomenon has a great influence on crack tip parameters decreasing their values; Numerical ++
[35] CWI Kawabata et al. 2016 Crack Wake Influence A new factor f is introduced to correct the plastic term. In this factor, the blunted crack tip shape is considered to depend on the strain hardening exponent, and f is given as a function of the yield-to-tensile ratio (YR) of the material and the specimen thickness. Numerical-Experimental +
[38] ρ Pommier and Risbet 2005 In the equations, special attention is paid to the elastic energy stored inside the crack tip plastic zone, sync, in practice, residual stresses at the crack tip are known to considerably influence fatigue crack growth Analytical +++
[39] Noroozi et al. 2005 The results demonstrate the crack closure influence for the LCF behaviour. The change of crack closure from LCF to high cycle fatigue and their consequences for lifetime prediction Analytical-experimental ++
[40] Rpr Ould Chikh et al. 2008 Plastic zone Size The cyclic plastic strain can be the principal parameter for the fatigue crack growth under a cyclic loading. Generally, FCGR is a plastic zone size rc function, and it increases as the plastic zone size increases. Analytical +++
[44] De Matos et al. 2008 This paper shows that the residual stress field due to cold expansion has a strong influence on the closure behaviour and therefore on fatigue crack propagation. Numerical-Analytical ++
[52] εPA Borges et al. 2020 Fatigue crack growth (FCG) is simulated here by node release, which is made when the accumulated plastic strain reaches a critical value. Numerical +++
[54] Rpc Park et al. 1996 Plastic Zone Size Experimental tests showed that Plastic zone Size was an important parameter in crack propagation Experimental +++
[55] ry Donald and Paris 1999 The ACR and CWI methods measure the change in displacement at minimum load due to closure That quantity is less subject to variability than is the measurement of the opening load. Analytical +++
[58] da/dS Zhang et al. 2010 da/dS, defines the fatigue crack propagation rate with the change of the applied stress at any moment of a stress cycle the relationship between this new parameter and the conventional da/dN is given. Numerical ++
[60] Rpc Zhang et al. 2010 Plastic zone Size The results have shown that the near crack tip the reverse plastic zone size continues to change with the change of the applied compressive stress Numerical-Experimental ++
[63] PICC Pokluda 2011 Plasticity induced in the crack closure There is a good qualitative agreement between the plasticity-induced shielding terms employed in the dislocation-based model and the continuum-based multi-parameter model. Analytical +
[64] Rpr Antunes et al. 2018 Reverse Plastic zone Size the increase of crack plastic deformation also produces an increase of crack closure phenomenon, which cancels the increase of plastic deformation Numerical ++
[79] ρ Korsunsky et al. 2009 Plastic blunting, crack tip blunting Two approaches were considered in the present study: the approach based on the consideration of crack tip blunting due to Pommier and Risbet [38], and the presently proposed approach based on the analysis of local energy dissipation in the immediate vicinity of the crack tip. Analytical-Experimental ++
Table 5. Energy as a critical parameter used as driving force in the crack grown modelisation.
Table 5. Energy as a critical parameter used as driving force in the crack grown modelisation.
Ref Specific Parameter Authors and Date Description Definition Methodology
[64] CTODp Antunes et al. 2018 Plastic Energy Dissipated per Cycle Plastic portion of the crack tip opening displacement Plastic CTOD is obtained by subtracting the elastic CTOD from the total. Numerical
[82,83] PM Zheng et al. 2013, 2014 Critical plastic energy, at the point close to crack tip The fatigue damage experienced by a point located within a specific distance from the crack tip can accurately represent the average damage condition at the crack tip area. Numerical
[85] EC Kujawsky and Ellyin 1984 Critical plastic energy Amount of plastic strain energy that a material can dissipate before experiencing fatigue failure Numerical
[86] EC Chalant and Remy 1983 Critical plastic energy The strain gradient inside the grain at the crack tip Numerical
[89] ΔWt Branco et al. 2021 Critical plastic energy Accounts for the mean stress effect, measure of the energy dissipated per cycle and is capable of unifying both the low-cycle and high-cycle fatigue regimes. Numerical-Analytical
[88] χW Zhu et al. 2018 Critical plastic energy Reflects the distribution of both stress and strain gradients within the actual structure. Analytical
[79] ΔWp Konsunsky et al. 2009 Equivalent deformation energy Amount of energy dissipated due to plastic deformation at the crack tip during each loading cycle. Analytical-Experimental
[75] ΔWp Klingbei 2003 Strain energy gradient Change in total plastic dissipation per unit width during a specific cycle. Numerical-Analytical
[76] UL Ravi Chandran 2018 Total dissipated plastic energy Total net-section strain energy in the crack plane of the ligament. Combination of elastic and plastic strain energies due to the increased stress. Analytical-Experimental
[76] ΔCpσ Ravi Chandran 2018 Total dissipated plastic energy Change in net-section strain energy parameter in stress-controlled fatigue.Eq. 19 Analytical-Experimental
[76] ΔCpε Ravi Chandran 2018 Cumulative change in cyclic strain energy of the net section Change in net-section strain energy density in strain-controlled fatigue.Eq. 20 Analytical-Experimental
[77] Upl Quan and Alderliesten 2022 Plastic energy difference in the net-section Energy consumed in the process of crack growth through plastic deformation of the material surrounding the crack. Numerical-Experimental
[77] Ue Quan and Alderliesten 2022 Elastic energy difference in the net-section Variation in elastic strain energy stored throughout one full cycle of crack propagation. Numerical-Experimental
[77] Ua Quan and Alderliesten 2022 Surface energy difference in the net-section Surface energy differential dissipated through new crack surface formation Numerical-Experimental
[78] ΔWp Kheli et all. 2013 Stored deformation dissipation Cyclic plastic strain energy, corresponding to one loading cycle.Eq. 21 Numerical-Experimental
[78] U Kheli et all. 2013 Dissipation new crack surface formation Specific energy, energy dissipated per unit volume during fatigue crack growth.Eq. 22 Numerical-Experimental
[78] Q Kheli et all. 2013 Plastic energy of the hysteresis cycle characteristic of cyclic loads Total dissipated energy in the specimen during fatigue crack growthEq. 23 Numerical-Experimental
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated