A numerical study of fatigue crack closure induced by plasticity

Crack closure delays the intrinsic mechanisms responsible for crack growth, therefore, it must be considered in fatigue crack growth modelling. The objective of this work is to develop a numerical procedure to predict crack closure induced by plasticity. First the crack closure was experimentally measured on M(T) 6082‐T6 aluminium alloy specimens of 3 mm thickness. A pin microgauge was used with the compliance technique. Then different parameters of the numerical procedure were analysed, namely the finite element mesh and the crack propagation scheme. The size of crack‐tip elements has an important influence and it is recommended to be of the same order of cyclic plastic zone. Crack‐opening levels only 10% lower than experimental results were obtained considering kinematic hardening and two load cycles in each increment.     

A numerical analysis of the mechanisms behind plasticity induced crack closure: Application to variable amplitude loadings

The effect of loading parameters on fatigue crack growth has been explained using the concept of crack closure. Plasticity induced crack closure (PICC) is linked to the crack tip plastic deformation, which becomes residual with crack propagation. The objective here is to identify the main mechanisms behind PICC, and for that different loading cases were considered namely overloads and load blocks. An analytical model was used to isolate the effect of residual plastic deformation on PICC, however significant differences were obtained relatively to finite element results. A second mechanism, which is crack tip blunting, was used to explain the transient behaviour observed after overloads and load blocks. For overloads and low–high load sequences there is a sudden increase of crack tip blunting with load increase which explains the sudden decrease of crack opening level. For high–low load sequences there is a sudden decrease of crack tip blunting which enhances the effect of residual plastic wake. Finally, the partial closure concept was tested looking to non-linear crack tip parameters but no evidences of Donald’s effect were found for the material studied.     

Influence of through-thickness crack shape on plasticity induced crack closure

In most of the previous three‐dimensional (3D) numerical studies of plasticity induced crack closure (PICC), ideal shapes have been assumed for the cracks. The aim of present paper is to study the effect of crack shape on PICC. With this objective a 3D numerical model was developed to predict PICC in middle‐tension (MT) specimens with different thicknesses and crack shapes. The radial size of crack tip elements and the stabilization of closure level were studied to ensure the quality of numerical predictions. Simultaneously, an independent numerical model was developed to predict crack shape evolution, stable crack shapes and corresponding K distributions. Crack closure was found to produce a significant tunnelling effect, with maximum values of ΔK and K_max at the surface. The curved crack presented significant plastic deformation near the free surface which has a high impact on the computation time, compared to the straight crack. The modification of ΔK and K_max with crack shape produced a variation of 38% in opening values at the interior positions, but relatively small variations at the surface. Considering the great influence of crack shape on PICC, it is fundamental to model realistic crack shapes.     

A parameter for quantitative analysis of plasticity induced crack closure

Numerical models have been successfully developed to predict plasticity induced crack closure (PICC). However, despite the large research effort a full understanding of the links between physical parameters, residual plastic wake and PICC has not been achieved yet. The plastic extension of material behind crack tip, Δy_p, obtained by the integration of vertical plastic deformation perpendicularly to crack flank, is proposed here to quantify the residual plastic field. The values of Δy_p and PICC were obtained numerically in a M(T) specimen using the finite element method. An excellent correlation was found between PICC and Δy_p which indicates that this parameter controls the phenomenon, and can be used to quantify the effect of physical parameters. An empirical model was developed to predict PICC assuming that the residual plastic field is a set of vertical plastic wedges, that the linear superposition principle applies and that the influence of a particular wedge exponentially decreases with distance to crack tip. The model was applied successfully to predict PICC for different residual plastic fields which provided an additional validation of Δy_p as the parameter controlling PICC.     

Plastic mismatch effect on plasticity induced crack closure: Fatigue crack propagation perpendicularly across a plastically mismatched interface

In the present work, comprehensive investigation of both theoretical analysis and numerical simulation was carried out to investigate the plastic mismatch effect on plasticity induced crack closure (PICC) behavior and effective fatigue crack tip driving force. During the process of crack tip approaching interface, crack tip load and crack tip load ratio will change, resulting in the change of PICC degree. When the crack propagates towards higher strength side, K_op/K_max increases; when the crack propagates towards lower strength side, K_op/K_max decreases firstly and then increases. The two mechanisms of “interface plastic mismatch effect on nominal fatigue crack tip driving force” and “interface plastic mismatch effect on PICC degree” were compared. The second mechanism must be considered when building crack tip driving force model for describing fatigue crack crossing plastically mismatched interface, because it is more physically factual and maybe more important than the first mechanism.     

A numerical study of plasticity induced crack closure under plane strain conditions

The level of plasticity induced crack closure (PICC) is greatly affected by stress state. Under plane strain conditions, however, the level and even the existence of PICC still are controversial. The objective here is to study the influence of the main numerical parameters on plane strain PICC, namely the total crack propagation, the number of load cycles between crack increments, the finite element mesh and the parameter used to quantify PICC. The PICC predictions were included in a parallel numerical study of crack propagation, in order to quantify the impact of plane strain values on fatigue life. The results indicate that literature may be overestimating plane strain PICC due to incorrect numerical parameters. The number of load cycles usually considered is unrealistically small, and its increase was found to vanish crack closure, particularly for kinematic hardening. This effect was linked to the ratcheting effect observed at the crack tip. The total crack increment, Δa, must be large enough to obtain stabilized PICC values, but this may imply a huge numerical effort particularly for 3D models. The size of crack tip plastic zone may be overestimated in literature, which means that the meshes used may be too large. Additionally, the crack propagation study showed that the plane strain PICC has usually a dominant effect on fatigue life, and plane stress PICC is only relevant for relatively thin geometries.     

A numerical investigation on the effect of an overload on fatigue crack opening and closure behaviour

The effect of a superimposed single tensile overload or a compressive overload in a block of constant-amplitude cycles on the crack opening and closing stresses is investigated using an elastic–plastic finite element analysis. The results obtained are in basic agreement with the experimental observations. Following an applied tensile overload cycle, the crack opening and closing stresses increase instantaneously, while the imposition of a compressive overload cycle results in a small decrease of the crack opening and closing stresses. A detailed discussion of the residual stress and strain patterns caused by the plastic deformation during the overload cycle and the corresponding crack profiles is included. The main governing mechanism resulting in the change of crack growth due to an overload is pointed out.     

Fatigue crack closure: a review of the physical phenomena

Plasticity‐induced, roughness‐induced and oxide‐induced crack closures are reviewed. Special attention is devoted to the physical origin, the consequences for the experimental determination and the prediction of the effective crack driving force for fatigue crack propagation. Plasticity‐induced crack closure under plane stress and plane strain conditions require, in principle, a different explanation; however, both types are predictable. This is even the case in the transition region from the plane strain to the plane stress state and all types of loading conditions including constant and variable amplitude loading, the short crack case or the transition from small‐scale to large‐scale yielding. In contrast, the prediction of roughness‐induced and oxide‐induced closures is not as straightforward.     

Plasticity induced crack closure in Middle‐Crack Tension specimen: numerical versus experimental

Numerical studies played a major role in the current understanding of plasticity induced crack closure (PICC), however it is well known that the numerical models simplify the reality by considering discrete crack propagations, relatively high fatigue crack growth rates, sharp cracks and that propagation occurs at a well‐defined load. The objective of this paper is to find if despite these limitations, the numerical predictions of PICC are close to the experimental results. Experimental work was developed to obtain crack opening values in 1‐mm‐thick middle‐tension (M(T)) specimens of AA6016‐T4 aluminium alloy. The opening level was determined from remote compliance data captured by a pin microgauge placed at the centre of the specimen. A numerical model was developed and optimized replicating the experimental test in terms of sample geometry, loading parameters, material behaviour and crack opening measurement technique. The adoption, in the numerical analysis, of a remote measurement technique for determining PICC was found to be robust and adequate for mesh refinement studies. A good agreement was found between experimental and numerical results.     

Numerical simulation of plasticity induced crack closure: Identification and discussion of parameters

Numerical studies play a major role in the understanding and prediction of plasticity induced crack closure (PICC). However, the available numerical models can be considered simplifications of reality as they consider discrete crack propagations, relatively high fatigue crack growth rates (FCGR), sharp cracks, and propagation occurring at well-defined loads. Besides, there are a great number of numerical and physical parameters affecting the predictions of PICC. The aim of this paper is to discuss the numerical study of PICC. The numerical parameters affecting the accuracy of the numerical simulations, and the dependent parameters used to characterise the plastic wake and the closure level, are identified. The influence of the radial size of crack front elements and crack propagation is analysed. An extrapolation model is proposed, with excellent results. An intrinsic uncertainty is associated with the number of load cycles between crack increments and the definition of crack closure level. Finally, the effect of the stress ratio (R) on crack closure level is analysed.     

Some experimental observations on crack closure and crack-tip plasticity

This study presents a methodology for evaluating crack closure and the effect of crack-tip plasticity on stress intensity. Full-field displacement maps obtained by digital image correlation are used to obtain the mixed-mode, crack-driving force. The methodology allows the quantification of the effect of a range of contact phenomena: effects arising from interlocking, plastic deformation of crack face asperities and wedging generated as a consequence of sliding displacements of fatigue cracks have been identified. By evaluating the effective crack-tip stress intensity factor, crack opening levels can be quantified for both mode I and mode II. Moreover, the approach can take into account plasticity effects local to the crack in determining the stress intensity factor. All the information can be extracted in a non-contacting fashion with equipment that can be easily incorporated into industrial environments.     

Analytical and numerical modelling of plasticity-induced crack closure in cold-expanded holes

The aim of this research was the development of an analytical model for plasticity-induced fatigue crack closure for cold expanded holes. This paper extends Nowell's plane stress model of plasticity-induced crack closure for a plate with a circular hole and two radial symmetric cracks. The possibility of existence of an initial residual stress field is also taken into account. This model has potential to be applied to other cracked geometries and arbitrary residual stress fields, although the paper is focused on the study of cold-expanded holes. Hole cold-expansion is widely used in aircraft industry, for improving the fatigue performance of rivet holes by delaying fatigue crack propagation. 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. The analytical model developed, was compared with finite element analyses of plasticity-induced crack closure with and without residual stresses. Finally, the model was used to predict fatigue lives for some experiments recently reported in the literature for fatigue crack propagation from cold-expanded holes. Predicted fatigue lives correlate well with experimental data.     

Simulating small crack growth behaviour using crystal plasticity theory and finite element analysis

Predictions of small crack growth under cyclic loading in aluminium alloy 7075 are performed using finite element analysis (FEA), and results are compared with published experimental data. A double‐slip crystal plasticity model is implemented within the analyses to enable the anisotropic nature of individual grains to be approximated. Small edge‐cracks in a single grain with a starting length of 6 μm are incrementally grown following a node‐release scheme. Crack‐tip opening displacements (CTOD) and crack opening stresses are calculated during the simulated crack growth, and da/dN against ΔK diagrams are computed. Interactions between the crack tip and a grain boundary are also considered. The computations are shown to accurately capture the magnitude and the variability normally observed in small crack fatigue data.     

Empirical model for plasticity‐induced crack closure based on Kmax and ΔK

The mean stress has a significant effect on crack propagation life and must be included in prediction models. However, there is no consensus in the fatigue community regarding the dominant mechanism explaining the mean stress effect. The concept of crack closure has been widely used and several empirical models can be found in literature. The stress ratio, R, is usually the main parameter of these models, but present numerical results showed a significant influence of K_max. A new empirical model is therefore proposed here, dependent on K_max and ΔK, with four empirical constants. The model also includes the effect of material's yield stress, and two additional parameters were defined to account for stress state and crack closure parameter. A comparison was made with Kujawski's and Glinka's parameters, for a wide range of loading conditions. ΔK_eff lies between Kujawski's and Glinka's parameters, and some agreement is evident, although the concepts are quite different. The crack opening model was applied to literature results and was able to collapse da/dN–ΔK curves for different stress ratios to a single master curve.     

Numerical modelling of plane strain plasticity induced crack closure effects for bimaterial interfacial cracks

The effects of plane strain plasticity induced crack closure on fatigue cracks located at the interface of dissimilar steel materials are presented using finite element modelling. Based on the study, it has been observed that bimaterial cracks produced unsymmetrical residual plastic strains and crack profiles in the crack wakes. It is seen that Young’s modulus and yield stress mismatch have profound effects on the development of unsymmetrical residual plastic strain and crack profiles, whereas the effect of Poisson’s ratio is insignificant. However, it has been found that for the material properties considered, low value of crack closure levels have been identified.     

Finite-element analysis of fatigue crack closure under plane strain conditions: stabilization behaviour and mesh size effect

An elastic–plastic finite‐element analysis of fatigue crack closure is performed for plane strain conditions. The stabilization behaviour of crack opening level and the effect of mesh size on the crack opening stress are investigated. It has been well reported that the crack opening level under plane stress conditions becomes stable after the crack advances beyond the initial monotonic plastic zone. In order to obtain a stabilized crack opening level for plane strain conditions, the crack must be advanced through approximately four times the initial monotonic plastic zone. The crack opening load tends to increase with the decrease of mesh size. The mesh size nearly equal to the theoretical plane strain cyclic plastic zone size may provide reasonable numerical results comparable with experimental crack opening data. The crack opening behaviour is influenced by the crack growth increment and discontinuous opening behaviour is observed.     

Estimation of the effects of plasticity and resulting crack closure during small fatigue crack growth

On the basis of a plastic-strip model and the method of singular integral equations, a closed-form analytical solution of the problem of an elastic-plastic plate containing a rectilinear fatigue crack is considered. The solution is used for the prediction of fatigue growth of `mechanically-small' crack by accounting for reverse plastic yielding and plasticity-induced crack closure in the material. The main effects of these factors on the crack-growth rate are analyzed, and the predicted results are compared with experimental data on small fatigue-crack growth in a aluminum-lithium alloy 2091-T351 and Fe-3% Si alloy.     

An integrated methodology for separating closure and residual stress effects from fatigue crack growth rate data

To properly interpret the results of standard fatigue crack growth tests it is often necessary to incorporate corrective techniques to the ΔK applied data. This is especially true in the near‐threshold regime where long crack data need to be closure corrected to predict small crack behaviour. It is also an issue in the presence of residual stress. A methodology to separate the influence of sample size, geometry, crack length and residual stress from the standard crack growth test data to obtain a true material response is presented. Stress ratio and residual stress contributions from known combinations of assumed crack size, applied stress and residual stress are also addressed and incorporated in the fatigue crack growth behaviour.     

Effect of crack closure on non-linear crack tip parameters

Crack closure concept has been widely used to explain different issues of fatigue crack propagation. However, some authors have questioned the relevance of crack closure and have proposed alternative concepts. The main objective here is to check the effectiveness of crack closure concept by linking the contact of crack flanks with non-linear crack tip parameters. Accordingly, 3D-FE numerical models with and without contact were developed for a wide range of loading scenarios and the crack tip parameters usually linked to fatigue crack growth, namely range of cyclic plastic strain, crack tip opening displacement, size of reversed plastic zone and total plastic dissipation per cycle were investigated. It was demonstrated that: (i) LEFM concepts are applicable to the problem under study; (ii) the crack closure phenomenon has a great influence on crack tip parameters decreasing their values; (iii) the ΔK_eff concept is able to explain the variations of crack tip parameters produced by the contact of crack flanks; and (iv) the analysis of remote compliance is the best numerical parameter to quantify the crack opening level. Therefore the crack closure concept seems to be valid. Additionally, the curves of crack tip parameters against stress intensity factor range obtained without contact may be seen as master curves.     

Numerical simulation of plasticity induced fatigue crack opening and closure for autofrettaged intersecting holes

The autofrettage of intersecting holes leads to extremely high compressive residual stress fields. These stresses in combination with the plastic deformations decelerate fatigue cracks initiated at the hole intersection notch. Simulations of plasticity induced crack closure of such cracks are presented based on the strip yield and a finite element model. The strip yield model has been extended to allow for an input of residual stresses coming from elsewhere, e.g. from a finite element calculation or measurements. The calculations are applied for constant as well as variable amplitude loading. The numerical expense of the finite element based modelling for variable amplitude loading is still too high if millions of cycles have to be considered. Therefore, a new approximation method is proposed introducing compensatory load sequences. Simulation results are compared to experimentally determined results showing good agreement. However, the accuracy of crack initiation life estimates has turned out to provide a high potential for further improvement.