Simulation of plasticity-induced fatigue crack closure in part-through cracked geometries using finite element analysis

The part-through semi-elliptical surface flaw is commonly encountered in engineering practice. Proper characterization of plasticity-induced crack closure is necessary to predict both flaw growth and flaw shape evolution under cyclic loading. Three-dimensional elastic-plastic finite element analyses are used to model the plasticity-induced closure developed along the surface flaw crack front, and the subsequent crack opening behavior under constant amplitude loading. Resulting crack opening stresses are compared with results from a strip-yield model and with experimentally measured values reported in the literature. It was found that the computed values were larger than those measured.     

Measurement of fatigue crack growth in large compact tension specimens using an AC magnetic field method

Fatigue crack growth rate data are required in order to carry out a numerical analysis of the fatigue performance of complex structural components. These data are obtained by measuring crack growth in standard fracture mechanics specimens. A new method for measuring fatigue crack growth in compact tension specimens has been developed. The technique is based on the measurement of the surface magnetic fields produced when passing a high-frequency alternating current through the specimen. Fatigue crack growth data recorded using this method indicated an accuracy of ±0.02 mm when compared with optical measurements. The technique is suitable for computer-controlled operation and could easily be applied to other standard specimen geometries.     

Utilization of partial crack closure for fatigue crack growth modeling

Load ratio effects are of prime concern when modeling of fatigue crack growth (FCG) rate is required as a prerequisite for a reliable life prediction. The majority of research efforts regarding the load ratio effects are based on Elber's ΔK_eff approach. However, there are intrinsic difficulties encountered with its consistent application to FCG prediction. In this paper two popular crack-growth-life prediction codes FASTRAN and AFGROW are modified utilizing the enhanced partial crack closure model. The proposed utilization aggregates apparent closure mechanisms involved and demonstrates a better correlation and a significant scatter reduction of FCG data taken from literature, especially in the near-threshold region.     

Numerical predictions and experimental measurements of residual stresses in fatigue crack growth specimens

Controlling macro residual stress fields in a material while preserving a desired microstructure is often a challenging proposition. Processing techniques which induce or reduce residual stresses often also alter microstructural characteristics of the material through thermo-mechanical processes. A novel mechanical technique able to generate controlled residual stresses was developed. The method is based on a pin compression approach, and was used to produce well-controlled magnitudes and distributions of residual stresses in rectangular coupons and compact tension specimens typically used in fatigue crack growth testing. Residual stresses created through this method were first computationally modeled with finite element analysis, and then experimentally reproduced with various levels of pin compression. The magnitudes and distributions of residual stresses in experimental specimens were independently assessed with fracture mechanics methods and good correspondence was found between residual stresses produced using the pin compression and processing techniques. Fatigue crack growth data generated from specimens with low residual stresses, high residual stresses resulting from processing, and high residual stresses introduced through the new pin compression technique were compared and validated. The developed method is proposed to facilitate the acquisition and analysis of fatigue crack growth data generated in residual stresses, validate residual stress corrective models, and verify fatigue crack growth simulations and life predictions in the presence of residual stresses.     

Numerical and experimental analysis of residual stresses for fatigue crack growth

In this paper computational and experimental results are presented concerning residual stress effects on fatigue crack growth in a Compact Tension Shear (CTS) specimen under cyclic mode I loading. For a crack of constant length it is found that hardly any compressive residual stresses or crack closure effects are generated along the crack surfaces behind the crack tip through the considered cyclic mode I loading with a load ratio of R=0.1. Only if fatigue crack growth is modelled during the simulation of the cyclic loading process these well-known effects are found. On the other hand it is shown that they have hardly any influence on the residual stresses ahead of the crack tip and thus on further fatigue crack growth. For all cases considered the computational finite element results agree well with the experimental findings obtained through X-ray diffraction techniques.     

Variable amplitude crack growth in notched specimens

Predictions of variable amplitude fatigue crack growth in notched tensile specimens at an elevated temperature and relative high loads were conducted. The predictions were based on constant amplitude fatigue data and mean crack growth rates were calculated. Crack closure levels were calculated by comparison of crack closure free experimental data obtained with the so called constant K_I,max-method and nominal data. Also, numerical simulations of the closure level throughout the specimens were performed. Although experimental data showed rather large a scatter it was concluded that variable amplitude crack growth rates can be predicted with reasonable agreement using constant amplitude fatigue data. Also, crack growth rates cannot be predicted by LEFM in the inelastic notch affected zone, even if absence of crack closure is assumed.     

Crack growth retardation following the application of an overload cycle using a strip-yield model

Experiments have shown that the application of an overload cycle can act to retard crack growth and even potentially lead to crack arrest. This paper describes a new method for investigating fatigue crack growth after the application of an overload cycle under plane stress conditions. The developed method is based on the concept of plasticity-induced crack closure and utilises the distributed dislocation technique and a modified strip-yield model. The present results are compared to previous experimental data for several materials. A good agreement is found, with the predictions showing the same trends in the various stages of post-overload crack growth.     

Elastic–plastic finite element analysis of the effect of compressive loading on crack tip parameters and its impact on fatigue crack propagation rate

The effects of applied compressive stress on the crack tip parameters and their implications on fatigue crack growth have been studied. Four center-cracked panel specimens with different crack lengths are analysed using finite element method. The results show that under tension–compression loading the local crack tip parameters are determined by two loading parameters. The two loading parameters are the maximum stress intensity and the maximum compressive stress in the applied stress cycle. Predictions of fatigue crack propagation behaviour based on the maximum stress intensity and maximum compressive stress agree well with experimental observations.     

Singular and non-singular approaches for predicting fatigue crack growth behavior

In this work, three classes of mechanisms that can cause load sequence effects on fatigue crack growth are discussed: mechanisms acting before, at or after the crack tip. After reviewing the crack closure idea, which is based on what happens behind the crack tip, quantitative models are proposed to predict the effects at the crack tip due to crack bifurcation. To predict the behavior ahead of the crack tip, a damage accumulation model is proposed. In this model, fatigue cracking is assumed caused by the sequential failure of volume elements or tiny εN specimens in front of the crack tip, calculated by damage accumulation concepts. The crack is treated as a sharp notch with a small, but not zero radius, avoiding the physically unrealistic singularity at its tip. The crack stress concentration factor and a strain concentration rule are used to calculate the notch root strain and to shift the origin of a modified HRR field, resulting in a non-singular model of the strain distribution ahead of the crack tip. In this way, the damage caused by each load cycle, including the effects of residual stresses, can be calculated at each element ahead of the crack tip using the correct hysteresis loops caused by the loading. The proposed approach is experimentally validated and extended to predict fatigue crack growth under variable amplitude loading, assuming that the width of the volume element broken at each cycle is equal to the region ahead of the crack tip that suffers damage beyond its critical value. The reasonable predictions of the measured fatigue crack growth behavior in steel specimens under service loads corroborate this simple and clear way to correlate da/dN and εN properties.     

History effect in fatigue crack growth under mixed-mode loading conditions

Plastic deformation within the crack tip region introduces internal stresses that modify subsequent behaviour of the crack and are at the origin of history effects in fatigue crack growth. Consequently, fatigue crack growth models should include plasticity-induced history effects. A model was developed and validated for mode I fatigue crack growth under variable amplitude loading conditions. The purpose of this study was to extend this model to mixed-mode loading conditions. Finite element analyses are commonly employed to model crack tip plasticity and were shown to give very satisfactory results. However, if millions of cycles need to be modelled to predict the fatigue behaviour of an industrial component, the finite element method becomes computationally too expensive. By employing a multiscale approach, the local results of FE computations can be brought to the global scale. This approach consists of partitioning the velocity field at the crack tip into plastic and elastic parts. Each part is partitioned into mode I and mode II components, and finally each component is the product of a reference spatial field and an intensity factor. The intensity factor of the mode I and mode II plastic parts of the velocity fields, denoted by dρ_I/dt and dρ_II/dt, allow measuring mixed-mode plasticity in the crack tip region at the global scale. Evolutions of dρ_I/dt and dρ_II/dt, generated using the FE method for various loading histories, enable the identification of an empirical cyclic elastic–plastic constitutive model for the crack tip region at the global scale. Once identified, this empirical model can be employed, with no need of additional FE computations, resulting in faster computations. With the additional hypothesis that the fatigue crack growth rate and direction can be determined from mixed-mode crack tip plasticity (dρ_I/dt and dρ_II/dt), it becomes possible to predict fatigue crack growth under I/II mixed-mode and variable amplitude loading conditions. To compare the predictions of this model with experiments, an asymmetric four point bend test system was setup. It allows applying any mixed-mode loading case from a pure mode I condition to a pure mode II. Initial experimental results showed an increase of the mode I fatigue crack growth rate after the application of a set of mode II overload cycles.     

Predicting fatigue crack growth rate in a welded butt joint: The role of effective R ratio in accounting for residual stress effect

A simple and efficient method is presented in this paper for predicting fatigue crack growth rate in welded butt joints. Three well-known empirical crack growth laws are employed using the material constants that were obtained from the base material coupon tests. Based on the superposition rule of the linear elastic fracture mechanics, welding residual stress effect is accounted for by replacing the nominal stress ratio (R) in the empirical laws by the effective stress intensity factor ratio (R_eff). The key part of the analysis process is to calculate the stress intensity factor due to the initial residual stress field and also the stress relaxation and redistribution due to crack growth. The finite element method in conjunction with the modified virtual crack closure technique was used for this analysis. Fatigue crack growth rates were then calculated by the empirical laws and comparisons were made among these predictions as well as against published experimental tests, which were conducted under either constant amplitude load or constant stress intensity factor range. Test samples were M(T) geometry made of aluminium alloy 2024-T351 with a longitudinal weld by the variable polarity plasma arc welding process. Good agreement was achieved.     

Evaluation of fatigue crack growth behavior of GLARE3 fiber/metal laminates using a compliance method

The objectives of this study were to investigate the effectiveness of a compliance method for analyzing the fatigue crack growth of GLARE3 fiber/metal laminates. The materials tested were GLARE3-5/4 (2.6 mm thick) and GLARE3-3/2 (1.4 mm thick). Centrally notched specimens with two kinds of notch length and two kinds of fiber orientation were fatigue tested under constant amplitude loading. The expression of the experimental stress intensity factor, K_exp, for the 2024-T3 aluminum-alloy layers of a GLARE3 is formulated and K_exp were obtained from the relationship between crack length and specimen compliance. The test results clarified the following: (1) da/dN–ΔK_exp relationships roughly show the linear relationship independent of the maximum stress level, specimen thickness, notch length, and fiber orientations, (2) the da/dN–ΔK_exp relationships approximately agree with the linear part and its extension of Paris–Erdogan’s law obtained for the da/dN–ΔK relationship of the 2024-T3 aluminum-alloy, (3) the compliance method is effective for analyzing fatigue crack growth in GLARE3 laminates.     

A three-dimensional (3D) numerical study of fatigue crack growth using remeshing techniques

Numerical analyses based on the finite element (FE) method and remeshing techniques have been employed in order to develop a damage tolerance approach to be used for the design of aeroengines shaft components. Preliminary experimental tests have permitted the calculation of fatigue crack growth parameters for the high strength alloy steel adopted in this research. Then, a robust numerical study have been carried out to understand the influence of various factors (such as: crack shape, crack closure) on non-planar crack evolution in solid and hollow shafts under mixed-mode loading. The FE analyses have displayed a satisfactory agreement compared to experimental data on compact specimens (CT) and solid shafts.     

Evolution of residual stresses with fatigue loading and subsequent crack growth in a welded aluminium alloy middle tension specimen

Neutron diffraction has been used to measure the evolution of the residual stresses in a VPPA welded Al-2024 alloy middle tension (M(T)) specimen with fatigue loading and subsequent crack growth. The measurements were carried out on the diffractometer ENGIN-X, a time-of-flight instrument, at the ISIS Pulsed Neutron Source. Fatigue crack growth was performed in situ and strain measurements averaged through the thickness of the specimen were made along two orthogonal directions as the crack grew, allowing the stresses in the specimen to be calculated assuming plane stress. 2D finite element simulation of the evolution of the initial residual stress field with crack growth, using an elastic model produced predictions that were in reasonable agreement with the experimental results. The results further indicate that some re-distribution of the residual stress field occurred due to the crack tip plasticity associated with the fatigue loading.     

Role of elasto-plastic analysis under cyclic loading in fatigue crack growth studies

Linear Elastic Fracture Mechanics (LEFM) has been widely used in the past for fatigue crack growth studies, but this is acceptable only in situations which are within small scale yielding (SSY). In many practical structural components, conditions of SSY could be violated and one has to look for fracture criteria based on elasto-plastic analysis. Crack closure phenomenon, one of the most striking discoveries based on inelastic deformations during crack growth, has significant effect on fatigue crack growth rate. Numerical simulation of this phenomenon is computationally intensive and involved but has been successfully implemented. Stress intensity factors and strain energy release rates lose their meaning,J-integral (or its incremental) values are applicable only in specific situations, whereas alternate path independent integrals have been proposed in the literature for use with elasto-plastic fracture mechanics (EPFM) based criteria. This paper presents certain salient features of two independent finite element (numerical) studies of relevance to fatigue crack growth, where elasto-plastic analysis becomes significant. These problems can only be handled in the current day computational environment, and would have been only a dream just a few years ago. The work presented in this paper is supported by sponsored research projects of the Aeronautics R & D Board, Government of India and their support is acknowledged.     

Fatigue crack closure determination by means of finite element analysis

Plasticity-induced crack closure is an observed phenomenon during fatigue crack growth. However, accurate determination of fatigue crack closure has been a complex task for years. It has been approached by means of experimental and numerical methods. The finite element method (FEM) has been the principal numerical tool employed. In this paper the results of a broad study of fatigue crack closure in plane stress and plane strain by means of FEM are presented. The effect of three principal factors has been analysed in depth, the maximum load, the crack length and the stress ratio. It has been found that the results are independent of maximum load and the crack length, and there exists a direct influence of the stress ratio. This relation has been numerically correlated and compared with experimental results. Differences have also been established between opening and closure points and between the different criteria employed to compute crack closure.     

Investigation of fatigue crack closure using multiscale image correlation experiments

Two full-field macroscale methods are introduced for estimating fatigue crack opening levels based on digital image correlation (DIC) displacement measurements near the crack tip. Crack opening levels from these two full-field methods are compared to results from a third (microscale) method that directly measures opening of the crack flanks immediately behind the crack tip using two-point DIC displacement gages. Of the two full-field methods, the first one measures effective stress intensity factors through the displacement field (over a wide region behind and ahead of the crack tip). This method reveals crack opening levels comparable to the limiting values (crack opening levels far from the crack tip) from the third method (microscale). The second full-field method involves a compliance offset measurement based on displacements obtained near the crack tip. This method delivers results comparable to crack tip opening levels from the microscale two-point method. The results of these experiments point to a normalized crack tip opening level of 0.35 for R ∼ 0 loading in grade 2 titanium. This opening level was found at low and intermediate ΔK levels. It is shown that the second full-field macroscale method indicates crack opening levels comparable to surface crack tip opening levels (corresponding to unzipping of the entire crack). This indicates that effective stress intensity factors determined from full-field displacements could be used to predict crack opening levels.     

Predicting corner crack fatigue propagation from cold worked holes

We predict the fatigue propagation of corner cracks from cold worked holes using three dimensional finite element models. The models account for the through thickness variation in residual stress left after cold working. The predictions are compared to experimental results in aluminum 2024-T351 and 7075-T651. The models show the evolution of P-shaped crack fronts similar to those observed in experiments. Predictions based on the initial residual stress field left after cold working were nonconservative, predicting either slower than experimental crack growth or crack growth that arrests. Predictions based on an estimate of the stable relaxed residual stress field near the hole were conservative, and predicted 5-10 times greater life than the current Department of Defense reduced initial flaw size approach.     

Plastic dissipation in mixed-mode fatigue crack growth along plastically mismatched interfaces

A new theory of fatigue crack growth in ductile solids has recently been proposed based on the total plastic energy dissipation per cycle ahead of the crack. This and previous energy based approaches in the literature suggest that the total plastic dissipation per cycle can be closely correlated with fatigue crack growth rates under mode I loading. In a recent paper, the authors have extended the dissipated energy approach to the case of fatigue crack growth in a homogeneous material under sustained mixed-mode loading conditions. The goal of the current study is to further extend the approach to mixed-mode fatigue delamination of ductile interfaces in layered materials. Attention is restricted to material combinations with identical elastic properties, but with mismatches in plastic properties (both yield strength and hardening modulus) across the interface. Such systems can occur in brazing, soldering, welding, and a variety of layered manufacturing applications, where high-temperature material deposition can result in a mismatch in mechanical properties between the deposited material and the substrate. In this study, the total plastic dissipation per cycle is obtained through plane strain elastic–plastic finite element analysis of a stationary crack in a general layered specimen geometry under constant amplitude, mixed-mode loading. Numerical results for a dimensionless plastic dissipation per cycle are presented over the full range of relevant material combinations and mixed-mode loading conditions. Results suggest that while applied mode-mix ratio is the dominant parameter, mismatches in yield strength and hardening modulus can have a significant effect on the total plastic dissipation per cycle, which is dominated by the weaker/softer material.     

A procedure for the simulation of fatigue crack growth in adhesively bonded joints based on the cohesive zone model and different mixed-mode propagation criteria

This work deals with the simulation of the fatigue crack growth (FCG) in bonded joints. In particular a cohesive damage model is implemented in the commercial software Abaqus, in order to take into account for the damage produced by fatigue loading. The crack growth rate is evaluated with different Paris-like power laws expressed in terms of strain energy release rate. The crack growth rate is then translated into a variation of the damage distribution over the cohesive zone setting an equivalence between the increment of crack length and the increment of damage. The model takes also into account mixed mode I/II conditions. In this work the validity of the model is tested by comparison with theoretical trends for conditions of pure mode I, pure mode II and mixed mode loading. In the case of mixed mode conditions, different models are implemented for the crack growth rate computation. The results of the model are in very good agreement with the expected trends, therefore the model is adequate to simulate the fatigue crack growth behaviour of bonded joint.