Near-tip strain evolution under cyclic loading: In situ experimental observation and numerical modelling

The concept of ratchetting strain as a crack driving force in controlling crack growth has previously been explored at Portsmouth using numerical approaches for nickel-based superalloys. In this paper, we report the first quantitative experimental evidence of near-tip strain ratchetting with cycles, as captured in situ by digital image correlation (DIC) technique on a compact tension specimen of stainless steel 316L, using both Stereo and SEM systems. The evolution of the near-tip strains with loading cycles was monitored whilst the crack tip was kept stationary. The strains normal to the crack plane were examined over selected distances from 6 to 57 μm to the crack tip for a number of cycles. The results show that strain ratchetting occurs with loading cycles, and is particularly evident close to the crack tip and under higher loads. 3D finite element models have also been developed to simulate the experiments and the results from the simulation are compared with those from the DIC measurements.This is the first time that near-tip strain ratchetting has been captured in situ at the peak loads during cyclic loading.     

Plastic strain-controlled short crack growth and fatigue life

Constant plastic strain-controlled and constant stress-controlled tests were performed on smooth and lightly notched specimens machined from a massive forging of 42CrMo4 steel. Comparison of the fatigue life curves plotted as function of the plastic strain amplitude and stress amplitude shows a decisive role of plastic strain amplitude. Crack initiation and the kinetics of short crack growth were studied in constant plastic strain amplitude loading and the relation between the crack growth coefficient and plastic strain amplitude was established. This is equivalent to the Coffin–Manson law and shows that the Coffin–Manson law can be interpreted in terms of short crack growth.     

Simulation of Fatigue Crack Propagation in Ductile Metals by Blunting and Re-sharpening

Laird and Smith [(1962). Philosophical Magazine 8, 847–857] proposed a plastic sliding-off mechanism for the stage II fatigue crack growth via striation formation. In their view, the fatigue crack extension results solely from the changing character of deformation at the crack tip during loading and unloading. In particular, the crack tip blunts during the loading stage and folds into a double notch during the unloading stage, resulting in striation formation. In order to verify Laird’s plastic blunting mechanism for ductile polycrystals as well as for ductile fcc single crystals, FE calculations were performed for a rectangular plate with an initially sharp crack under plane strain conditions. The plate was subjected to a fully reversed tension-to-pressure cyclic load perpendicular to the crack plane (Mode 1). In the single crystal case the crack propagation simulations were carried out for cracks with crack plane (001) for two different crack growth orientations [110] and [100]. No initial radius for the crack tip was assumed. The actual shape of the crack tip followed from an initially sharp crack by repeated remeshing. To model the constitutive behavior typical for polycrystalline ductile metals, J₂ hypo-elasto-plasticity model with Armstrong–Frederick kinematic hardening was used. To model the constitutive behavior typical for ductile fcc single crystals, a geometrically nonlinear version of Cailletaud’s model based on the multiplicative elasto-plastic decomposition of the deformation gradient was implemented into the FE program ABAQUS. For simplicity, only octahedral slip systems were considered. Using repeated remeshing for severely distorted elements at the advancing crack tip, deformation patterns in the sense of Laird’s mechanism for fatigue crack propagation with striation formation were obtained in the case of the polycrystal simulation as well as in the case of the single crystal simulation for [110] crack growth direction. The simulation for [100] crack growth direction with the same stress level as for [110] direction also yielded crack extension by progressive large deformations but without striation formation. The dependence of the fatigue striation formation on the crack growth direction as predicted by the simulation of crack propagation in single crystals is verified by the experimental results of Neumann [(1974). Acta Metallurgica 22, 1155–1165] on pure copper single crystals.     

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.     

In situ observation of electric-field-induced domain switching and crack propagation in poled PMNT62/38 single crystals

It is important to in situ observe the crack propagation and the corresponding domain switching in ferroelectrics subjected to electric loading since domain switching has been widely assumed to play a critical role in the electric-field-induced crack growth. In this investigation, we in situ observed the crack propagation and the domain switching in PMNT62/38 single crystals poled along the [001] orientation. An experimental setup was designed and constructed to investigate the crack propagation and the domain switching in thin plate specimens with pre-crack subjected to electric field by using polarized light microscope (PLM). The pre-crack began to propagate forward accompanied by the appearance of domain switching zones near the crack tip and the disappearance of switched zones behind the crack tip at the unipolar electric field of E = 0.8E_C. The results indicate that the structure mismatch of the adjacent switched zones with different polarizations stimulated by the intensive electric field near the crack tip results in the electric-field-induced crack growth.     

Prediction of crack growth in a nickel-based superalloy under fatigue-oxidation conditions

Prediction of oxidation-assisted crack growth has been carried out for a nickel-based superalloy at elevated temperature based on finite element analyses of oxygen diffusion, coupled with viscoplastic deformation, near a fatigue crack tip. The material constitutive behaviour, implemented in the finite element code ABAQUS via a user-defined material subroutine (UMAT), was described by a unified viscoplastic model with non-linear kinematic and isotropic hardening rules. Diffusion of oxygen was assumed to be controlled by two parameters, the oxygen diffusivity and deformation-assisted oxygen mobility. Low frequencies and superimposed hold periods at peak loads significantly enhanced oxygen concentration near the crack tip. Evaluations of near-tip deformation and oxygen concentration were performed, which led to the construction of a failure envelop for crack growth based on the consideration of both oxygen concentration and accumulated inelastic strain near the crack tip. The failure envelop was then utilised to predict crack growth rates in a compact tension (CT) specimen under fatigue-oxidation conditions for selected loading ranges, frequencies and dwell periods. The predictions from the fatigue-oxidation failure envelop compared well with the experimental results for triangular and dwell loading waveforms, with marked improvements achieved over those predicted from the viscoplastic model alone. The fatigue-oxidation predictions also agree well with the experimental results for slow-fast loading waveforms, but not for fast-slow waveforms where the effect of oxidation is much reduced.     

Biaxial cyclic deformation of an epoxy resin: Experiments and constitutive modeling

Biaxial (proportional and non-proportional) cyclic tests were conducted on thin-walled tubular specimens to investigate deformation behavior of an epoxy resin, Epon 826/Epi-Cure Curing Agent 9551. The focus was placed on the biaxial stress-strain response and their dependency on the load control mode, stress or strain range and loading path. Experimental results indicated that under strain-controlled equi-biaxial (proportional) cyclic loading, mean stress relaxation occurred in both axial and hoop directions, whereas under stress-controlled equi-biaxial cyclic loading, ratcheting strains accumulated in both principal directions. Under strain- or stress-controlled non-proportional cyclic loading, anisotropy in stress-strain responses was induced in both axial and hoop directions, and the axial and hoop hysteresis loops rotated in opposite directions. This was particularly evident at high stress or strain levels. The experimental results were further used to evaluate the predictive capabilities of a nonlinear viscoelastic constitutive model. Qualitative and quantitative comparison with the test data indicated a good agreement in predicting the complex stress-strain response under biaxial cyclic loading with various loading paths, applied stress or strain ranges and loading control modes.     

A crystal plasticity study of cyclic constitutive behaviour, crack-tip deformation and crack-growth path for a polycrystalline nickel-based superalloy

Crystal plasticity has been applied to model the cyclic constitutive behaviour of a polycrystalline nickel-based superalloy at elevated temperature using finite element analyses. A representative volume element, consisting of randomly oriented grains, was considered for the finite element analyses under periodic boundary constraints. Strain-controlled cyclic test data at 650 °C were used to determine the model parameters from a fitting process, where three loading rates were considered. Model simulations are in good agreement with the experimental results for stress–strain loops, cyclic hardening behaviour and stress relaxation behaviour. Stress and strain distributions within the representative volume element are of heterogeneous nature due to the orientation mismatch between neighbouring grains. Stress concentrations tend to occur within “hard” grains while strain concentrations tend to locate within “soft” grains, depending on the orientation of grains with respect to the loading direction. The model was further applied to study the near-tip deformation of a transgranular crack in a compact tension specimen using a submodelling technique. Grain microstructure is shown to have an influence on the von Mises stress distribution near the crack tip, and the gain texture heterogeneity disturbs the well-known butterfly shape obtained from the viscoplasticity analysis at continuum level. The stress–strain response near the crack tip, as well as the accumulated shear deformation along slip system, is influenced by the orientation of the grain at the crack tip, which might dictate the subsequent crack growth through grains. Individual slip systems near the crack tip tend to have different amounts of accumulated shear deformation, which was utilised as a criterion to predict the crack growth path.     

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.     

Low-cycle fatigue crack growth prediction by strain intensity factor

The validity of the strain intensity factor for representing driving force of fatigue crack growth under the large scale yielding condition was shown in this study. First, a crack growth test technique using plate specimens was developed in order to apply a fully-reversed cyclic load and to measure the global strain for calculating the strain intensity factor. Fatigue crack growth tests using Type 316 stainless steel revealed that the growth rates correlated well with the strain and effective strain intensity factor ranges for various stress or strain ranges, specimen geometries and loading modes (stress- and strain-controlled). Based on the test results, validity and physical meaning of the strain intensity factor were discussed. It was concluded that the strain intensity factor was an effective parameter for predicting crack growth under the small and large scale yielding conditions.     

Numerical simulations of dynamic crack growth along an interface

Dynamic crack growth is analyzed numerically for a plane strain bimaterial block with an initial central crack. The material on each side of the bond line is characterized by an isotropic hyperelastic constitutive relation. A cohesive surface constitutive relation is also specified that relates the tractions and displacement jumps across the bond line and that allows for the creation of new free surface. The resistance to crack initiation and the crack speed history are predicted without invoking any ad hoc failure criterion. Full finite strain transient analyses are carried out, with two types of loading considered; tensile loading on one side of the specimen and crack face loading. The crack speed history and the evolution of the crack tip stress state are investigated for parameters characterizing a PMMA/Al bimaterial. Additionally, the separate effects of elastic modulus mismatch and elastic wave speed mismatch on interface crack growth are explored for various PMMA-artificial material combinations. The mode mixity of the near tip fields is found to increase with increasing crack speed and in some cases large scale contact occurs in the vicinity of the crack tip. Crack speeds that exceed the smaller of the two Rayleigh wave speeds are also found.     

Modeling of uniaxial and biaxial ratcheting behavior of 1026 Carbon steel using the simplified Viscoplasticity Theory Based on Overstress (VBO)

Summary.  Hassan and Kyriakides [1], Hassan et al. [2], and Corona et al. [3] performed uniaxial and biaxial ratcheting experiments on heat-treated 1026 Carbon steel. The loading histories performed with uniaxial and tubular specimens were selected to simulate those encountered in nuclear reactor vessels. The stress-strain diagram of 1026 Carbon steel was used to determine the material constants in a simplified version of VBO. Small rate dependence was allowed. The model represents some nonlinear rate dependence, kinematic hardening and cyclically neutral behavior. The set of material constants determined only from uniaxial tests was used throughout the paper. Numerical experiments included are: (i) uniaxial stress-controlled cycling with various mean stresses, (ii) strain-controlled axial tests with tubular specimens under constant and variable internal pressure, and (iii) examination and variation of certain material constants of the VBO model that can influence ratcheting. Very good agreement with the experiments is found for the uniaxial case. However the ratchet strain accumulation during biaxial cycling is over-predicted by VBO and other constitutive models. Received November 5, 2001; revised June 12, 2002 Published online: January 16, 2003     

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.     

Atomistic modeling of crack growth in magnesium single crystal

The analysis of crack growth in magnesium single crystal was performed using molecular dynamics simulation with Embedded Atom Method (EAM) potentials. The twinning process at the crack tip was analyzed. Four specimens with increasing sizes were used to investigate the influences of material length scale on crack growth of magnesium single crystals. Furthermore, the effects of temperature, and the loading strain rate were also verified. The specimens were subjected to uniaxial tension strain up to the total strain level of 0.2 with a constant strain rate. In the simulation of each specimen, the averaged stress strain curve was monitored. The simulation results show that the specimen size, loading strain rate, and temperature strongly influence the peak stress point at which the twin nucleated and subsequently the crack grew. The initial slope of the averaged stress strain curve is independent of the loading strain rate and temperature. Moreover, high temperatures induce increased atomic mobility, and thereby atom reorganization, which, in turn, releases the stress at the crack tip.     

Cyclic loadings and crystallization of natural rubber: An explanation of fatigue crack propagation reinforcement under a positive loading ratio

Natural rubber is known to have excellent fatigue properties. Fatigue crack propagation studies show that, under uniaxial tension loading, fatigue crack growth resistance increases with the loading ratio, even if the peak stress increases. Studies dealing with crack initiation confirm this trend. If strain induced crystallization is believed to play a major role in this reinforcement process, it is not clear yet by which mechanism this reinforcement takes place. Using SEM investigation, it is shown here that the reinforcement process is associated with strong crack branching in the crack tip region. From experimental results it is shown that under particular reinforcing loading condition a cyclic strain hardening process can be observed on the natural rubber which is able to overcome classically observed softening effects. A cumulative strain induced crystallization process is proposed to explain the stress ratio effect on fatigue crack initiation and propagation properties of natural rubber.     

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.     

Assessment of small fatigue crack growth driving forces in single crystals with and without slip bands

Strain localization under low amplitude cyclic loading is a manifestation of plastic irreversible deformation associated with early crack growth. However, traditional constitutive models cannot usually reproduce strain localization in smooth single crystals, which can affect crack growth predictions for crystallographic fatigue cracks. This work analyzes the influence of bands of localized plastic shear strain on the cyclic crack tip displacement and on a fatigue indicator parameter by making special provision of a crack along the interface of a deformation band. Furthermore, the quality of local and volume-averaged fatigue indicator parameters are assessed using finite element models of a Cu single crystal cycled to induce plastic deformation under multiple loading conditions.     

Crack observation methods, their application and simulation of curved fatigue crack growth

The observation of cracks with curved crack path has to be done optically. It is shown that a modified commercial flat bed scanner and a combination of a high-resolution scanner camera system with 10,000 × 10,000 pixels with a telecentric lens are appropriate for high-resolution (up to 8 μm/pixel) optical recording of multiple crack ends on sample areas up to 210 mm × 290 mm. The high-resolution photographs are suitable for determination of crack lengths. It is also possible to observe crack paths or geometrical crack tip parameters and strain fields by image correlation. The method is used to determine static crack resistance and cyclic crack growth curves on center crack tension and biaxial cruciform samples.Furthermore, the paper presents an improved finite element technique for the simulation of curved fatigue crack growth in a multiple arbitrarily pre-cracked isotropic sheet under biaxial plane stress loading applying a predictor-corrector procedure in combination with the modified virtual crack closure integral (MVCCI) method including the consideration of the plastic limit loads. For this, the program PCCS-2D was extended to analyse the crack growth and the plastic limit load for each crack propagation step in a fully automatic simulation. The proposed solution algorithm provides a powerful tool for flaw assessment with the failure assessment diagram procedure in combination with a numerical crack path simulation. Finally, the simulation is verified experimentally.     

Dynamic fracture of a beam or plate under tensile loading

The dynamic fracture response of a long beam of brittle elastic material under tensile loading is studied. If the magnitude of the applied loading is increased to a critical value, a crack is assumed to propagate across the beam cross section. In a parallel analysis to [t] the crack length and applied loading at the fracture face are determined as functions of time measured from fracture initiation. The results of the analysis are shown in graphs of crack length, crack tip speed and fracturing section tensile loading vs time. As found in [1], the crack tip accelerates very quickly to a speed near the characteristic terminal speed for the material, travels at this speed through most of the beam thickness, and then decelerates rapidly in the final stage of the process. Finally, by appropriate change of the elastic modulus, the results may be applied to plane strain fracture of a plate under pure tensile loading.     

THE SIGNIFICANCE OF CRACK TIP DEFORMATION FOR SHORT AND LONG FATIGUE CRACKS

Abstract— The behaviour of physical short mode I cracks under constant amplitude cyclic loading was investigated both numerically and experimentally. A dynamic two-dimensional elastic-plastic finite element technique was utilised to simulate cyclic crack tip plastic deformation. Different idealisations were investigated. Both stationary and artificially advanced long and short cracks were analysed. A parameter which characterises the plastically deformed crack tip zone, the strain field generated within that zone and the opening and closure of the crack tip were considered. The growth of physically short mode I cracks under constant amplitude fully reversed fatigue loading was investigated experimentally using conventional cast steel EN-9 specimens. Based on a numerical analysis, a crack tip deformation parameter was devised to correlate fatigue crack propagation rates.