Analysis of temperature distribution near the crack tip under constant amplitude loading

An analytical/numerical method has been developed to find the temperature rise near the crack tip under fatigue loading. The cyclic plastic zone ahead of the crack tip is assumed to be the shape of the source of heat generation and some fraction of plastic work done in cyclic plastic zone as heat generation. Plastic work during fatigue load was found by obtaining stress and strain distribution within the plastic zone by Hutchinson, Rice and Rosengren (HRR) crack tip singularity fields applied to small scale yielding on the cyclic stress strain curve. A two‐dimensional conduction heat transfer equation, in moving co‐ordinates, was used to obtain temperature distribution around the crack tip. Temperature rise was found to be a function of frequency of loading, applied stress intensity factor and thermal properties of the material. A power–law relation was found between the rise in temperature at a fixed point near the crack tip and range of stress intensity factor.     

Cyclic crack opening displacement behavior during high-amplitude block loading

Load-crack opening displacement hysteresis behavior was monitored for fifty cycles of highamplitude loading which followed fatigue pre-cracking at low stress intensity factor levels. The material studied was quenched and tempered (400°C) AISI 4140 steel which showed pronounced cyclic softening. Despite this softening behavior, cycle-to-cycle decreases in load-COD hysteresis were observed during the initial cycles of high-amplitude loading. Steady state load-COD hysteresis behavior was attained by fifty loading cycles in each case and the fifty-cycle hysteresis loop widths agreed well with those for continuously cycled (non-pre-cracked) samples for equivalent loading conditions. The cycles during which the load-COD hysteresis decreased most dramatically represented fatigue crack growth distances equal to approximately 30% of the calculated plane strain monotonic plastic zone size. Greater percentage reductions in load-COD hysteresis were observed for lower stress intensity factor ranges. The observed behavior was in general agreement with that predicted by finite element fatigue crack closure models in the literature. In addition, the level of prior loading was found to have a pronounced effect on subsequently measured fracture toughness values for this material.     

Fatigue crack growth model for constant amplitude loading

A fatigue crack growth model under constant amplitude loading has been developed considering energy balance during growth of the crack. The plastic energy dissipated during growth of a crack within cyclic plastic zone and area below cyclic stress–strain curve was used in the energy balance. The near crack tip elastic–plastic stress and strain were calculated on the basis of Hutchinson, Rice and Rosengren (HRR) formulations. Fatigue crack growth rate in linear and near threshold region of da/dN versus ΔK curve can be determined on the basis of the proposed model in terms of low cycle fatigue (LCF) properties determined on smooth specimen. The predictions of the model have been compared with the experimental and theoretical results available in the literature using mechanical and fatigue properties. The model compares well in the threshold and intermediate region of the da/dN versus ΔK curve for wide range of material tested.     

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.     

Mechanisms and modeling of low cycle fatigue crack propagation in a pressure vessel steel Q345

The fatigue process near crack is governed by highly concentrated strain and stress in the crack tip region. Based on the theory of elastic–plastic fracture mechanics, we explore the cyclic J-integral as breakthrough point, an analytical model is presented in this paper to determine the CTOD for cracked component subjected to cyclic axial in-plane loading. A simple fracture mechanism based model for fatigue crack growth assumes a linear correlation between the cyclic crack tip opening displacement (ΔCTOD) and the crack growth rate (da/dN). In order to validate the model and to calibrate the model parameters, the low cycle fatigue crack propagation experiment was carried out for CT specimen made of Q345 steel. The effects of stress ratio and crack closure on fatigue crack growth were investigated by elastic–plastic finite element stress–strain analysis of a cracked component. A good comparison has been found between predictions and experimental results, which shows that the crack opening displacement is able to characterize the crack tip state at large scale yielding constant amplitude fatigue crack growth.     

Thermal effect of plastic dissipation at the crack tip on the stress intensity factor under cyclic loading

Plastic dissipation at the crack tip under cyclic loading is responsible for the creation of an heterogeneous temperature field around the crack tip. A thermomechanical model is proposed in this paper for the theoretical problem of an infinite plate with a semi-infinite through crack under mode I cyclic loading both in plane stress or in plane strain condition. It is assumed that the heat source is located in the reverse cyclic plastic zone. The proposed analytical solution of the thermo-mechanical problem shows that the crack tip is under compression due to thermal stresses coming from the heterogeneous stress field around the crack tip. The effect of this stress field on the stress intensity factor (its maximum and its range) is calculated analytically for the infinite plate and by finite element analysis. The heat flux within the reverse cyclic plastic zone is the key parameter to quantify the effect of dissipation at the crack tip on the stress intensity factor.     

Plastic strain distribution at the tip of a fatigue crack. application to fatigue crack closure in the threshold regime

The fatigue crack growth rate behaviour of a Co-33 wt pct Ni alloy was investigated at room temperature down to the threshold regime using CT specimens for two load ratios 0.1 and 0.3. Cyclic equivalent plastic strain distributions along an axis normal to the crack plane were experimentally determined over the whole range of crack growth rates using two techniques microhardness and a quantitative metallographic technique applied to twins, both calibrated on low cycle fatigue specimens. These experimental values were compared with theoretical curves as obtained from the monotonic plane strain finite element analysis of Tracey and adapted to cyclic loading according to the procedure proposed by Rice. A good agreement was found in stage II crack growth in the vicinity of the crack tip but a discrepancy was observed in the low crack growth rate regime, indicative of crack closure. It was possible to determine the effective amplitude of the stress intensity factor which accounts for this discrepancy and an intrinsic crack growth law was obtained which obeys Paris equation and which applies in the whole crack growth rate range independently of the load ratio.     

Model of fatigue crack propagation by damage accumulation at the crack tip

Modelling fatigue crack propagation by damage accumulation at the crack tip originally proposed by McClintock is reworked out using metallurgical considerations. On a physical basis it is shown that the validity of such models is actually confined to the low crack growth rate range corresponding to stage I growth such as along crystallographic planes. At higher crack growth rates there is a transition to plastic stretching mechanisms which could be described by crack opening displacement models. An evaluation of two models has been carried out on a ³³Co Ni alloy where extensive information was available: both are based on the strain singularity as computed by Tracey from a finite element analysis of plane strain small scale yielding and as adapted to cyclic loading under Rice's hypothesis and taking the grain as the critical element below which size continuum mechanics do no longer apply. It is pointed out that the use of this strain singularity which is not able to account for crack closure, renders models unable to predict experimental crack growth rates curves but only intrinsic curves relating the growth rate to the effective stress intensity amplitude as defined by Elber. The first model which assumes fatal cracking of a grain with an average uniform cyclic equivalent strain, underestimates the crack growth rates and in addition it yields results very sensitive to the shape and size of the grain. The second model which assumes progressive cracking in a grain with a cyclic equivalent strain gradient is able to predict the intrinsic crack growth rate curve in the ³³Co Ni alloy and to yield predictions consistent with other experimental data.     

Biaxial fatigue for proportional and non‐proportional loading paths

In order to study the use of a local approach to predict crack‐initiation life on notches in mechanical components under multiaxial fatigue conditions, the study of the local cyclic elasto‐plastic behaviour and the selection of an appropriate multiaxial fatigue model are essential steps in fatigue‐life prediction. The evolution of stress–strain fields from the initial state to the stabilized state depends on the material type, loading amplitude and loading paths. A series of biaxial tension–compression tests with static or cyclic torsion were carried out on a biaxial servo‐hydraulic testing machine. Specimens were made of an alloy steel 42CrMo4 quenched and tempered. The shear stress relaxations of the cyclic tension–compression with a steady torsion angle were observed for various loading levels. Finite element analyses were used to simulate the cyclic behaviour and good agreement was found. Based on the local stabilized cyclic elastic–plastic stress–strain responses, the strain‐based multiaxial fatigue damage parameters were applied and correlated with the experimentally obtained lives. As a comparison, a stress‐invariant‐based approach with the minimum circumscribed ellipse (MCE) approach for evaluating the effective shear stress amplitude was also applied for fatigue life prediction. The comparison showed that both the equivalent strain range and the stress‐invariant parameter with non‐proportional factors correlated well with the experimental results obtained in this study.     

Digital image correlation measurement of near‐tip fatigue crack displacement fields: constant amplitude loading and load history effects

Recent work by de Matos and colleagues employed digital image correlation to measure near tip displacement fields for fatigue cracks in 6082 T6 aluminium alloy. The main focus of this work was to directly measure fatigue crack closure, but the measurements can also be used to examine conditions at and ahead of the crack tip. In this paper, the results are re‐analysed and compared to two crack‐tip deformation models. The first assumes simple elastic deformation (according the Westergaard solution). This allows the history of crack‐tip stress intensity to be examined. Reasonable agreement with the elastic model is obtained, although there is a residual stress intensity caused by the plastic wake, which gives rise to crack closure. The second model examined is a simple elastic–plastic assumption, proposed by Pommier and colleagues. This can be applied to constant amplitude loading, although the results obtained here are very similar to the elastic case. A slightly more complex load case (a single overload in an otherwise constant amplitude variation of load) gives a much more complicated crack‐tip history. Here, the importance of crack‐tip plastic displacement, represented by the second term in Pommier's model becomes much clearer. Load history effects are captured by the residual value of this term and its associated displacement fields as well as by stress intensity factor. The implications for further modelling and experimental work are discussed.     

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.     

Application of the distributed dislocation technique for calculating cyclic crack tip plasticity effects

This paper describes a method for modelling cyclic crack tip plasticity effects based on the distributed dislocation technique (DDT). A strip‐yield model is utilised to allow for the determination of the crack opening displacement, size of the plastic zones and in the case of a fatigue crack, the wake of plasticity. The DDT can be easily implemented for a wide range of cracked geometries with reliable control over the accuracy and convergence. Thickness effects can also be incorporated through a recently obtained solution for an edge dislocation in an infinite plate of finite thickness. Results for finite length cracks that have had limited growth, such that there is no plastic wake, are presented for a range of applied loads and R‐ratios. Further results are provided for a steady‐state fatigue crack in a plate of finite thickness. The present results are compared with analytical solutions and they show an excellent agreement.     

Acoustic emission behaviour during stage II fatigue crack growth in an AISI type 316 austenitic stainless steel

Acoustic emission (AE) behaviour during fatigue crack growth (FCG) in a ductile AISI type 316 austenitic stainless steel is reported. The two substages in the stage II Paris regime of FCG could be distinguished by a change in the rate of acoustic activity with increase in crack growth rate. The transition point in the cumulative ringdown count plot coincides with that in the da/dn plot. The AE activity increases with increase in ΔK during stage IIa and decreases during stage IIb. The major source of AE during stage IIa is found to be the plastic deformation within the cyclic plastic zone (CPZ) as compared to the phenomena such as monotonic plastic zone (MPZ) expansion, ductile crack growth, crack closure, etc. The increase in AE activity with increase in ΔK during stage IIa is attributed to the increase in the size of the CPZ which is generated and developed only under plane strain conditions. The decrease in AE activity during stage IIb is attributed to the decrease in the size of the CPZ under plane stress condition. The high acoustic activity during the substage IIa is attributed to irreversible cyclic plasticity with extensive multiplication and rearrangement of dislocations taking place within the CPZ. The AE activity is found to strongly depend on the optimum combination of the volume of the CPZ, average plastic strain range and the number of cycles before each crack extension. Based on this, an empirical relationship between the cumulative RDC and ΔK has been proposed and is found to agree well with experimentally observed values.     

Local crack plasticity and its influences on the global elastic stress field

This paper presents the background and development of a novel ‘plastic inclusion’ approach for dealing with the local plasticity which occurs at the tip of a growing fatigue crack. Localised plasticity arises from crack growth mechanisms and essentially blunts the crack, creates a reversed cyclic plastic zone, and induces shear along the crack flanks, along with the possible generation of wake contact stresses which act on the applied elastic stress field at the boundary of the elastic–plastic enclave surrounding the crack. The paper outlines the development of a meso-scale model of the elastic stress field around a growing crack that explicitly incorporates these interaction effects. The outcome is a modified crack tip stress intensity factor that includes some aspects of the magnitude of plastic wake-induced crack tip shielding and which the authors propose has the potential to help resolve some long-standing controversies associated with plasticity-induced closure. A full-field approach is developed for stress using photoelasticity and also for displacement using digital image correlation.     

Effects of grain size on cyclic plasticity and fatigue crack initiation in nickel

Fatigue experiments were conducted on polycrystalline nickel of two grain sizes, 24 and 290 μm, to evaluate the effects of grain size on cyclic plasticity and fatigue crack initiation. Specimens were cycled at room temperature at plastic strain amplitudes ranging from 2.5×10⁻⁵ to 2.5×10⁻³. Analyses of the cyclic stress–strain response and evolution of hysteresis loop shape indicate that the back stress component of the cyclic stress is significantly affected by grain size and plastic strain amplitude, whereas these parameters have little effect on friction stress. A nonlinear kinematic hardening framework was used to study the evolution of back stress parameters with cumulative plastic strain. These are related to substructural evolution features. In particular, long range back stress components are related to persistent slip bands. The difference in cyclic plasticity behavior between the two grain sizes is related to the effect of grain size on persistent slip band (PSB) morphology, and the effect this has on long range back stress. Fine grain specimens had a much longer fatigue life, especially at low plastic strain amplitude, as a result of the influence of grain size on fatigue crack initiation characteristics. At low plastic strain amplitude (2.5×10⁻⁴), coarse grain specimens initiated cracks where PSBs impinged on grain boundaries. Fine grain specimens formed cracks along PSBs. At high plastic strain amplitude (2.5×10⁻³), both grain sizes initiated cracks at grain boundaries.     

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.     

Finite deformation cohesive polygonal finite elements for modeling pervasive fracture

In the present study, mode I crack subjected to cyclic loading has been investigated for plastically compressible hardening and hardening–softening–hardening solids using the crack tip blunting model where we assume that the crack tip blunts during the maximum load and re-sharpening of the crack tip takes place under minimum load. Plane strain and small scale yielding conditions have been assumed for analysis. The influence of cyclic stress intensity factor range (\(\Delta \hbox {K})\), load ratio (R), number of cycles (N), plastic compressibility (\({\upalpha })\) and material softening on near tip deformation, stress–strain fields were studied. The present numerical calculations show that the crack tip opening displacement (CTOD), convergence of the cyclic trajectories of CTOD to stable self-similar loops, plastic crack growth, plastic zone shape and size, contours of accumulated plastic strain and hydrostatic stress distribution near the crack tip depend significantly on \(\Delta \hbox {K}\), R, N, \({\upalpha }\) and material softening. For both hardening and hardening–softening–hardening materials, yielding occurs during both loading and unloading phases, and resharpening of the crack tip during the unloading phase of the loading cycle is very significant. The similarities are revealed between computed near tip stress–strain variables and the experimental trends of the fatigue crack growth rate. There was no crack closure during unloading for any of the load cycles considered in the present study.     

A new model for describing a stable cyclic stress–strain relationship under non-proportional loading based on activation state of slip systems

This paper proposes a new simple model for cyclic incremental plasticity based on activation states of slip systems describing stable cyclic stress–strain relationships under non‐proportional loading. In the model, the magnitude and the direction of incremental plastic strain are estimated by (1+αf_NP) and Q , respectively. Here, α is the constant related to the dependence of material on additional hardening and f_NP the intensity factor expressing the severity of non‐proportional loading. Q is the second‐order tensor describing the activation states of slip systems in polycrystalline metals and is given by the calculation using a virtual specimen. The model was examined by application to the prediction of the stable cyclic stress–strain relationship in extensive non‐proportional low cycle fatigue tests for type 304 stainless steel and 6061 aluminium alloy. The simulated results showed that the model gave a satisfactory prediction of the stable cyclic stress–strain relationship under complex non‐proportional multiaxial loadings for the two materials.     

A study of crack growth retardation due to artificially induced crack surface contact

The contact of the cracked surfaces during a part of a loading cycle generally results in a reduced crack growth rate. A critical experiment was designed to evaluate the influence of the crack surface contact on crack growth. A round compact specimen made of 1070 steel with a round hole at the wake of the fatigue crack was designed. Two mating wedges were inserted into the hole of the specimen while the external load was kept at its maximum in a loading cycle. In this way, the wedges and the hole in the specimen were in firm contact during the entire loading cycle in the subsequent loading. Experiments showed that the addition of the wedges resulted in a reduction of crack growth rate in the subsequent constant amplitude loading. However, crack growth did not arrest. With the increase in the subsequent loading cycles, crack growth rate increased. The traditional crack closure concept cannot explain the experimental phenomenon because the effective stress intensity factor range was zero after the insertion of the wedges. The detailed stress–strain responses of the material near the crack tip were analyzed by using the finite element method with the implementation of a robust cyclic plasticity theory. A multiaxial fatigue criterion was used to determine the fatigue damage based upon the detailed stresses and strains. The crack growth was simulated and the predicted results were in good agreement with the experimental observations. It was confirmed that the stresses and strains near the crack tip governed cracking behavior. Crack surface contact reduced the crack tip cyclic plasticity and the result was the observed retardation in crack growth.     

Fatigue crack closure study based on whole-field displacements

This study is concerned with crack tip strain field fluctuations at loads below the point of crack closure in fatigue cycling. Moiré interferometry was used to investigate crack tip fields in compact tension specimens, cracked under constant stress intensity range and fixed R-ratio conditions. An elastic-plastic finite element model of simulated closure was developed to provide a theoretical cross-reference for the moiré studies. The ‘stretched zone’, which is believed to be the most significant source of closure effects, was simulated by inserting a constant thickness strip of elements into the crack before unloading from the maximum load point. Analysis of the crack tip fields in the experimental and theoretical cases was made in terms of crack face opening profiles, compliance changes and elastic stress intensity parameters. The latter were inferred through stress and displacement measurements made along circular and radial paths relative to the crack tip. Closure on the stretched zone was found to generate non-proportional loading in the crack tip field, so that the resulting stress changes were not well characterized by the asymptotic elastic equations. It is concluded firstly, that significant strain fluctuations occur below the point of closure load and that these should not be ignored in crack propagation studies. Secondly, the effective stress intensity range in fatigue cycling is not simply related to the open-crack stress intensity range and the need therefore remains for R-ratio and geometry effects to be treated as variables in crack propagation data collection programmes.