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.     

Fatigue life and crack closure in specimens subjected to variable amplitude loads under plane strain conditions

The propagation of fatigue cracks in specimens subjected to variable amplitude loading under plane strain conditions was investigated experimentally and numerically, to find the influence of the variable amplitude loading on the stabilised crack closure level. Experiments on four-point-bend specimens with a Gurney block load scheme, showed that the crack closure level depends on the crack length but not on the stress range of the fluctuations. Numerical simulations performed in the fatigue crack growth program FASTRAN-II showed good agreement with the experimental results. In addition, statistical uncertainty analyses performed on the fatigue life show that, for technical applications, the uncertainties in initial crack length and load levels have a greater influence on the uncertainty in fatigue life, than the fluctuation level of the load.     

The evolution of the stress–strain fields near a fatigue crack tip and plasticity-induced crack closure revisited

The evolution of the stress–strain fields near a stationary crack tip under cyclic loading at selected R‐ratios has been studied in a detailed elastic–plastic finite element analysis. The material behaviour was described by a full constitutive model of cyclic plasticity with both kinematic and isotropic hardening variables. Whilst the stress/strain range remains mostly constant during the cyclic loading and scales with the external load range, progressive accumulation of tensile strain occurs, particularly at high R‐ratios. These results may be of significance for the characterization of crack growth, particularly near the fatigue threshold. Elastic–plastic finite element simulations of advancing fatigue cracks were carried out under plane‐stress, plane‐strain and generalized plane‐strain conditions in a compact tension specimen. Physical contact of the crack flanks was observed in plane stress but not in the plane‐strain and generalized plane‐strain conditions. The lack of crack closure in plane strain was found to be independent of the material studied. Significant crack closure was observed under plane‐stress conditions, where a displacement method was used to obtain the actual stress intensity variation during a loading cycle in the presence of crack closure. The results reveal no direct correlation between the attenuation in the stress intensity factor range estimated by the conventional compliance method and that determined by the displacement method. This finding seems to cast some doubts on the validity of the current practice in crack‐closure measurement, and indeed on the role of plasticity‐induced crack closure in the reduction of the applied stress intensity factor range.     

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.     

On the accurate assessment of crack opening and closing stresses in plasticity-induced fatigue crack closure problems

The accurate calculation of the opening and closing stresses is an important issue in fatigue crack closure problems, since the effective driving force for crack growth is dependent on accurate calculation of the opening stresses. Often numerical methods such as finite element analysis are used to model plasticity-induced fatigue crack closure problems. There are many difficulties associated with this modelling work, since the results may depend on a wide range of parameters such as mesh refinement, node release scheme and modelling of the contact between the crack faces etc. Even after a great deal of modelling work some arbitrariness is evident in the technique used for assessing the opening and closing stresses. A number of techniques have been proposed in the literature and the current work will assess and compare these approaches. The node displacement method, the change in stresses at the crack tip, and the weight function technique will each be applied to a finite element model of a plane stress crack for a range of stress levels. In addition, an analytical model for plasticity-induced crack closure under plane stress conditions will be used to discuss the accuracy of these techniques. The investigation shows that all these techniques are equivalent provided that the displacement and stress at the crack tip are assessed accurately. However, it will be shown that use of the tensile tip stress method, proposed by some authors for assessing the closing stress, is erroneous.     

Experimental and numerical investigation of thickness effects in plasticity-induced fatigue crack closure

The understanding of fatigue crack closure has been proved to be a challenging and controversial topic among the fatigue community over the last three decades. The effect of the specimen (or component) thickness has been shown to have a significant effect on closure behaviour and this seems to be related to the relative size of the plastic zone. Real cracks are inherently three-dimensional; plane stress-like behaviour is found close to the region where the crack front intersects the free surface, whereas most of the crack front will experience something close to plane strain. The aim of the present work is to investigate the influence of specimen thickness on closure behaviour (both close to and remote from the surface) and on fatigue crack propagation. The paper will present results from a simple experimental program, which consists of fatigue testing CT specimens with different thicknesses. Fatigue crack propagation is measured optically. Crack closure is assessed using traditional compliance techniques (clip gauge and back face strain gauge) and Digital Image Correlation methods. Experimental results are compared with two and three-dimensional simulations of plasticity-induced fatigue crack closure. The implications of thickness effects for predicting the propagation of three-dimensional fatigue cracks are discussed.     

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.     

Computations of fatigue crack growth with strain gradient plasticity and an irreversible cohesive zone model

Computations of fatigue crack growth with a first-order strain gradient plasticity (SGP) model and an irreversible cohesive zone model are reported. SGP plays a significant role in the model predictions and leads to increased fatigue crack growth rates relative to predictions with classical plasticity. Increased magnitudes of tractions and material separation at the crack tip together with reduced crack closure appear as the cause for accelerated crack growth in SGP. Under plane strain conditions SGP appears as an essential feature of the development of the crack closure zone. Size effects are explored relative to changes in internal material length scale as well as to structural length scales.     

Effect of overload on crack closure in thick and thin specimens via digital image correlation

The displacement field of compact tension (CT) specimens have been mapped by digital image correlation (DIC) local to growing fatigue cracks to study overload effects for plane stress and plane strain. We have extracted crack opening displacement (ΔCOD) and stress intensity (K) determined by a Muskhelishvili fit to the crack tip displacement field to infer the closure load. In both cases a classical knee was observed upon unloading consistent with closure which disappeared during the accelerated growth following OL, before increasing during retardation. In both cases following OL the crack growth rate is perturbed for a distance similar to the plastic zone.     

Characteristics of fatigue crack propagation behaviour as identified by hysteresis loop at the crack tip

Fatigue crack propagation tests have been carried out under various load conditions. Hysteresis loops denoting the relationship between load and strain at the crack tip are obtained by using local compliance measurement. Crack growth acceleration, delayed retardation and non‐propagation phenomena are investigated by considering the variation of hysteresis loop expansion and hysteresis loop tail. Based on the physical meaning of hysteresis loops, two types of crack closure are ascertained and the effect of crack closure on fatigue crack propagation is studied. Results show that change of the effective amplitude of the stress intensity factor at the crack tip is the reason that crack propagation rates vary.     

Methods for crack opening load and crack tip shielding determination: a review

In this paper a review of the literature on crack closure/opening load and crack tip shielding effects determination methods is presented. Commonly used ‘subjective’ (visual) and ‘non‐subjective’ approaches have been included. Procedures associated with the determination of an effective crack driving force for both Elber type and that of partial (or incremental) crack closure models have been covered. Comparison among different methods of analyses based on compliance and fatigue crack growth rate measurements is discussed together with their implications and difficulties in fatigue crack growth correlations.     

Derivation of crack closure and crack growth rate data from effective-strain fatigue life data for fracture mechanics fatigue life predictions

This study deals with the behavior of short cracks growing out of notches. Three types of load histories are used: (a) a fully-reversed constant amplitude history; (b) a periodic compressive overload history consisting of repeated load blocks containing one fully-reversed constant amplitude yield–stress magnitude cycle (the overload) followed by a group of smaller constant amplitude cycles having the same maximum stress as the overload cycle; (c) and a service strain history. Procedures are presented for deriving crack closure data and crack growth rate vs effective stress intensity factor range data from data obtained by subjecting a small number of smooth laboratory specimens to simple periodic compressive overload tests to obtain closure-free strain-life data. These procedures are illustrated in an example in which fatigue life predictions are made for a service strain history applied to notched plate specimens. The fatigue life predictions based on the measured and the derived crack closure and crack growth rate data are in good agreement with the experimentally determined fatigue lives.     

Plane strain crack closure and cyclic hardening

This study is focused on the understanding of the mechanical effects of cyclic hardening on crack tip plasticity and on plasticity-induced crack closure. Various finite element analyses were conducted using abaqus. Cyclic hardening is found to affect both crack closure and the shape of the plastic zone at the crack tip. Crack growth modelling in plane strain conditions in a cyclically hardening material is discussed. An empirical formula is provided which allows the calculation of the crack tip plastic zone size under plane strain conditions in a cyclically hardening material. The effects of overloads are also examined.     

Fatigue crack growth and crack closure in an AlMgSi alloy

This study reports an experimental investigation of fatigue crack propagation in AlMgSi1-T6 aluminium alloy using both constant and variable load amplitudes. Crack closure was monitored in all tests by the compliance technique using a pin microgauge. For the constant amplitude tests four different stress ratios were analysed. The crack closure parameter U was calculated and related with Δ K and the stress ratio, R . The threshold of the stress intensity factor range, Δ K _th , was also obtained. Fatigue crack propagation tests with single tensile peak overloads have been performed at constant load amplitude conditions. The observed transient post overload behaviour is discussed in terms of the overload ratio, Δ K baseline level and R . The crack closure parameter U trends are compared with the crack growth transients. Experimental support is given for the hypothesis that crack closure is the main factor determining the transient crack growth behaviour following overloads on AlMgSi1-T6 alloy for plane stress conditions.     

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.     

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 visualization of fatigue crack closure in plane stress and plane strain

Elastic-plastic finite element simulations of growing fatigue cracks in both plane stress and plane strain are used as an aid to visualization and analysis of the crack closure phenomenon. Residual stress and strain fields near the crack tip are depicted by both color fringe plots and x-y graphs. Development of the residual plastic stretch in the wake of a growing plane stress fatigue crack is shown to be associated with the transfer of material from the thickness direction to the axial direction. Finite element analyses indicate that crack closure does occur under pure plane strain conditions. The development of the residual plastic stretch in plane strain is shown to be associated with the transfer of material from the in-plane transverse direction to the axial direction. This in-plane contraction also leads to the generation of complex residual stress fields. The total length of closed crack at minimum load in plane strain is shown to be a small fraction of the total crack length, especially for positive stress ratios. This suggests that experimental measurement of plane strain closure would be extremely difficult, and may explain why some investigators have concluded that closure does not occur in plane strain.     

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.     

An investigation of crack closure and the propagation of semi-elliptical fatigue cracks in Q1N (HY80) pressure vessel steel

The results of an experimental investigation of the effect of crack closure on the propagation of semi-elliptical fatigue cracks are presented. Load-shedding fatigue threshold tests were carried out at stress ratios of 0.2, 0.35, 0.5 and 0.7. Crack closure was measured at the surface and depth positions using backface strain gauges, near-tip gauges, and a clip gauge. Differences between the surface and depth growth behaviour are explained by considerations of the effects of the transition from plane stress conditions at the surface to plane strain conditions at the depth. The effects of stress ratios are attributed largely to differences in the crack opening displacement, which result in asperities coming into contact to induce roughness-induced crack closure.