Micromechanics Of Fatigue Crack-Growth At Low Stress Intensities In A High-Strength Steel

An attempt to systematically investigate the effects of microstructural parameters in influencing the resistance to fatigue crack growth (FCG) in the near-threshold region under three different temper levels has been made for a high strength low alloy steel to observe in general, widely different trends in the dependence of both the total threshold stress intensity range, DELTA-K(th) and the intrinsic or effective threshold stress intensity range, DELTA-K(eff-th) on the prior austenitic grain size (PAGS). While a low strain hardening microstructure obtained by tempering at high temperatures exhibited strong dependence of DELTA-K(th) on the PAGS by virtue of strong interactions of crack tip slip with the grain boundary, a high strength, high strain hardening microstructure as a result of tempering at low temperature exhibited a weak dependence. The lack of a systematic variation of the near-threshold parameters with respect to grain size in temper embrittled structures appears to be related to the wide variations in the amount of intergranular fracture near threshold. Crack closure, to some extent provides a basis on which the increases in DELTA-K(th) at larger grain sizes can be rationalised. This study, in addition, provides a wide perspective on the relative roles of slip behaviour embrittlement and environment that result in the different trends observed in the grain size dependence of near-threshold fatigue parameters, based on which the inconsistency in the results reported in the literature can be clearly understood. Assessment of fracture modes through extensive fractography revealed that prior austenitic grain boundaries are effective barriers to cyclic crack growth compared to martensitic packet boundaries, especially at low stress intensities. Fracture morphologies comprising of low energy flat transgranular fracture can occur close to threshold depending on the combinations of strain hardening behaviour, yield strength and embrittlement effects. A detailed consideration is given to the discussion of cyclic stress strain behaviour, embrittlement and environmental effects and the implications of these phenomena on the crack growth behaviour near threshold.

Effects of temper level on the dependence of fatigue crack growth threshold and crack closure on the prior austenitic grain size

An attempt has been made to systematically investigate the effects of microstructural parameters, such as the prior austenite grain size (PAGS), in influencing the resistance to fatigue crack growth (FCG) in the near-threshold region under three different temper levels in a quenched and tempered high-strength steel. By austenitizing at various temperatures, the PAGS was varied from about 0.7 to 96 μm. The microstructures with these grain sizes were tempered at 200 °C, 400 °C, and 530 °C and tested for fatigue thresholds and crack closure. It has been found that, in general, three different trends in the dependence of both the total threshold stress intensity range, ΔK th , and the intrinsic threshold stress intensity range, ΔK eff, th , on the PAGS are observable. By considering in detail the factors such as cyclic stress-strain behavior, environmental effects on FCG, and embrittlement during tempering, the present observations could be rationalized. The strong dependence of ΔK th and ΔK eff, th on PAGS in microstructures tempered at 530 °C has been primarily attributed to cyclic softening and thereby the strong interaction of the crack tip deformation field with the grain boundary. On the other hand, a less strong dependence of ΔK th and ΔK eff, th on PAGS is suggested to be caused by the cyclic hardening behavior of lightly tempered microstructures occurring in 200 °C temper. In both microstructures, crack closure influenced near-threshold FCG (NTFCG) to a significant extent, and its magnitude was large at large grain sizes. Microstructures tempered at the intermediate temperatures failed to show a systematic variation of ΔKth and ΔKeff, th with PAGS. The mechanisms of intergranular fracture vary between grain sizes in this temper. A transition from “microstructure-sensitive” to “microstructure-insensitive” crack growth has been found to occur when the zone of cyclic deformation at the crack tip becomes more or less equal to PAGS. Detailed observations on fracture morphology and crack paths corroborate the grain size effects on fatigue thresholds and crack closure.

A general procedure for modified crack closure integral in 3D problems with cracks

In linear elastic fracture mechanics (LEFM), Irwin's crack closure integral (CCI) is one of the signficant concepts for the estimation of strain energy release rates (SERR) G, in individual as well as mixed-mode configurations. For effective utilization of this concept in conjunction with the finite element method (FEM), Rybicki and Kanninen [Engng Fracture Mech. 9, 931 938 (1977)] have proposed simple and direct estimations of the CCI in terms of nodal forces and displacements in the elements forming the crack tip from a single finite element analysis instead of the conventional two configuration analyses. These modified CCI (MCCI) expressions are basically element dependent. A systematic derivation of these expressions using element stress and displacement distributions is required. In the present work, a general procedure is given for the derivation of MCCI expressions in 3D problems with cracks. Further, a concept of sub-area integration is proposed which facilitates evaluation of SERR at a large number of points along the crack front without refining the finite element mesh. Numerical data are presented for two standard problems, a thick centre-cracked tension specimen and a semi-elliptical surface crack in a thick slab. Estimates for the stress intensity factor based on MCCI expressions corresponding to eight-noded brick elements are obtained and compared with available results in the literature.

Modified Crack Closure Integral (MCCI) For 3-D Problems Using 20-Noded Brick Elements

The Modified Crack Closure Integral (MCCI) technique based on Irwin's crack closure integral concept is very effective for estimation of strain energy release rates G in individual as well as mixed-mode configurations in linear elastic fracture mechanics problems. In a finite element approach, MCCI can be evaluated in the post-processing stage in terms of nodal forces and displacements near the crack tip. The MCCI expressions are however, element dependent and require a systematic derivation using stress and displacement distributions in the crack tip elements. Earlier a general procedure was proposed by the present authors for the derivation of MCCI expressions for 3-dimensional (3-d) crack problems modelled with 8-noded brick elements. A concept of sub-area integration was proposed to estimate strain energy release rates at a large number of points along the crack front. In the present paper a similar procedure is adopted for the derivation of MCCI expressions for 3-d cracks modelled with 20-noded brick elements. Numerical results are presented for centre crack tension and edge crack shear specimens in thick slabs, showing a comparison between present results and those available in the literature.

Modified crack closure integral (MCCI) for 3-d problems using 20 noded brick elememnts

The Modified Crack Closure Integral (MCCI) technique based on Irwin's crack closure integral concept is very effective for estimation of strain energy release rates G in individual as well as mixed-mode configurations in linear elastic fracture mechanics problems. In a finite element approach, MCCI can be evaluated in the post-processing stage in terms of nodal forces and displacements near the crack tip. The MCCI expressions are however, element dependent and require a systematic derivation using stress and displacement distributions in the crack tip elements. Earlier a general procedure was proposed by the present authors for the derivation of MCCI expressions for 3-dimensional (3-d) crack problems modelled with 8-noded brick elements. A concept of sub-area integration was proposed to estimate strain energy release rates at a large number of points along the crack front. In the present paper a similar procedure is adopted for the derivation of MCCI expressions for 3-d cracks modelled with 20-noded brick elements. Numerical results are presented for centre crack tension and edge crack shear specimens in thick slabs, showing a comparison between present results and those available in the literature.

Asymptotic and finite element analyses of mode III dynamic crack growth at a ductile-brittle interface

In this work, dynamic crack growth along a ductile-brittle interface under anti-plane strain conditions is studied. The ductile solid is taken to obey the J(2) flow theory of plasticity with linear isotropic strain hardening, while the substrate is assumed to exhibit linear elastic behavior. Firstly, the asymptotic near-tip stress and velocity fields are derived. These fields are assumed to be variable-separable with a power singularity in the radial coordinate centered at the crack tip. The effects of crack speed, strain hardening of the ductile phase and mismatch in elastic moduli of the two phases on the singularity exponent and the angular functions are studied. Secondly, full-field finite element analyses of the problem under small-scale yielding conditions are performed. The validity of the asymptotic fields and their range of dominance are determined by comparing them with the results of the full-field finite element analyses. Finally, theoretical predictions are made of the variations of the dynamic fracture toughness with crack velocity. The influence of the bi-material parameters on the above variation is investigated.

Re-examination of crack opening model of interface fracture

The complex singularity associated with a crack at the interface between two dissimilar, isotropic and homogeneous materials leads to mathematical artefacts, such as stress oscillations and crack face interpenetrations in the vicinity of the crack tip. To avoid these unrealistic features, Sinclair (Sinclair GB. On the stress singularity at an interface crack. International Journal of Fracture 1980;16(2):111-9) assumed a finite crack opening angle (COA) such that the singularity lambda became real equal to 1/2. This paper extends the COA model by considering real singularities not necessarily equal to 1/2. When COA is 0 degrees: the interface crack singularity is complex with a real part equal to 1/2. On increasing COA, the imaginary part of the singularity decreases and becomes zero at a threshold value of COA; at this point, the singularity is a real, repeated value. A further increase in COA results in a pair of real singularities. Different crack opening configurations and material combinations are studied, and results presented for threshold COAs and associated values of singularity. Stress analyses for these three regimes: (a) complex, (b) real pair and (c) real repeated singularities, are reported. It is seen that additional complexities are present in the last case. Typical results for stress fields are also included for comparing with standard fields. (C) 1999 Elsevier Science Ltd. All rights reserved.

The effect of T-stress on plane strain dynamic crack growth in elastic-plastic materials

In this paper, the effects of T -stress on steady, dynamic crack growth in an elastic-plastic material are examined using a modified boundary layer formulation. The analyses are carried out under mode I, plane strain conditions by employing a special finite element procedure based on moving crack tip coordinates. The material is assumed to obey the J (2) flow theory of plasticity with isotropic power law hardening. The results show that the crack opening profile as well as the opening stress at a finite distance from the tip are strongly affected by the magnitude and sign of the T -stress at any given crack speed. Further, it is found that the fracture toughness predicted by the analyses enhances significantly with negative T -stress for both ductile and cleavage mode of crack growth.

Evaluation of Fracture Parameters by Double-G, Double-K Models and Crack Extension Resistance for High Strength and Ultra High Strength Concrete Beams

This paper presents the advanced analytical methodologies such as Double- G and Double - K models for fracture analysis of concrete specimens made up of high strength concrete (HSC, HSC1) and ultra high strength concrete. Brief details about characterization and experimentation of HSC, HSC1 and UHSC have been provided. Double-G model is based on energy concept and couples the Griffith's brittle fracture theory with the bridging softening property of concrete. The double-K fracture model is based on stress intensity factor approach. Various fracture parameters such as cohesive fracture toughness (4), unstable fracture toughness (K-Ic(c)), unstable fracture toughness (K-Ic(un)) and initiation fracture toughness (K-Ic(ini)) have been evaluated based on linear elastic fracture mechanics and nonlinear fracture mechanics principles. Double-G and double-K method uses the secant compliance at the peak point of measured P-CMOD curves for determining the effective crack length. Bi-linear tension softening model has been employed to account for cohesive stresses ahead of the crack tip. From the studies, it is observed that the fracture parameters obtained by using double - G and double - K models are in good agreement with each other. Crack extension resistance has been estimated by using the fracture parameters obtained through double - K model. It is observed that the values of the crack extension resistance at the critical unstable point are almost equal to the values of the unstable fracture toughness K-Ic(un) of the materials. The computed fracture parameters will be useful for crack growth study, remaining life and residual strength evaluation of concrete structural components.

Fatigue Crack Growth Study and Remaining Life Assessment of High Strength and Ultra High Strength Concrete Beams

This paper presents the details of crack growth study and remaining life assessment of concrete specimens made up of high strength concrete (HSC, HSC1) and ultra high strength concrete (UHSC). Flexural fatigue tests have been conducted on HSC, HSC1 and UHSC beams under constant amplitude loading with a stress ratio of 0.2. It is observed from the studies that (i) the failure patterns of HSC1 and UHSC beams indicate their ductility as the member was intact till the crack propagated up to 90% of the beam depth and (ii) the remaining life decreases with increase of notch depth (iii) the failure of the specimen is influenced by the frequency of loading. A ``Net K'' model has been proposed by using non-linear fracture mechanics principles for crack growth analysis and remaining life prediction. SIF (K) has been computed by using the principle of superposition. SIP due to the cohesive forces applied on the effective crack face inside the process zone has been obtained through Green's function approach by applying bi-linear tension softening relationship to consider the cohesive the stresses acting ahead of the crack tip. Remaining life values have been have been predicted and compared with the corresponding experimental values and observed that they are in good agreement with each other.

Finite element simulations of notch tip fields in magnesium single crystals

Recent experiments using three point bend specimens of Mg single crystals have revealed that tensile twins of {10 (1) over bar2}-type form profusely near a notch tip and enhance the fracture toughness through large plastic dissipation. In this work, 3D finite element simulations of these experiments are carried out using a crystal plasticity framework which includes slip and twinning to gain insights on the mechanics of fracture. The predicted load-displacement curves, slip and tensile twinning activities from finite element analysis corroborate well with the experimental observations. The numerical results are used to explore the 3D nature of the crack tip stress, plastic slip and twin volume fraction distributions near the notch root. The occurrence of tensile twinning is rationalized from the variation of normal stress ahead of the notch tip. Further, deflection of the crack path at twin-twin intersections observed in the experiments is examined from an energy standpoint by modeling discrete twins close to the notch root.

Wallner lines, crack velocity and mechanisms of crack nucleation and growth in a brittle bulk metallic glass

Mode I fracture experiments were conducted on brittle bulk metallic glass (BMG) samples and the fracture surface features were analyzed in detail to understand the underlying physical processes. Wollner lines, which result from the interaction between the propagating crack front and shear waves emanating from a secondary source, were observed on the fracture surface and geometric analysis of them indicates that the maximum crack velocity is similar to 800 m s(-1), which corresponds to similar to 0.32 times the shear wave speed. Fractography reveals that the sharp crack nucleation at the notch tip occurs at the mid-section of the specimens with the observation of flat and half-penny-shaped cracks. On this basis, we conclude that the crack initiation in brittle BMGs is stress-controlled and occurs through hydrostatic stress-assisted cavity nucleation ahead of the notch tip. High magnification scanning electron and atomic force microscopies of the dynamic crack growth regions reveal highly organized, nanoscale periodic patterns with a spacing of similar to 79 nm. Juxtaposition of the crack velocity with this spacing suggests that the crack takes similar to 10(-10) s for peak-to-peak propagation. This, and the estimated adiabatic temperature rise ahead of the propagating crack tip that suggests local softening, is utilized to critically discuss possible causes for the nanocorrugation formation. Taylor's fluid meniscus instability is unequivocally ruled out. Then, two other possible mechanisms, viz. (a) crack tip blunting and resharpening through nanovoid nucleation and growth ahead of the crack tip and eventual coalescence, and (b) dynamic oscillation of the crack in a thin slab of softened zone ahead of the crack-tip, are critically discussed. (C) 2014 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

An automatic algorithm for mixed mode crack growth rate based on drop potential method

A new automatic algorithm for the assessment of mixed mode crack growth rate characteristics is presented based on the concept of an equivalent crack. The residual ligament size approach is introduced to implementation this algorithm for identifying the crack tip position on a curved path with respect to the drop potential signal. The automatic algorithm accounting for the curvilinear crack trajectory and employing an electrical potential difference was calibrated with respect to the optical measurements for the growing crack under cyclic mixed mode loading conditions. The effectiveness of the proposed algorithm is confirmed by fatigue tests performed on ST3 steel compact tension-shear specimens in the full range of mode mixities from pure mode Ito pure mode II. (C) 2015 Elsevier Ltd. All rights reserved.

Photoelastic determination of stress intensity factors in patched cracked plates

n this paper, the influence of patch parameters on stress intensity factors in edge cracked plates is studied by employing transmission photoelasticity. Edge cracked plates made of photo-elastic material are patched on one side only by E glass-epoxy and carbon-epoxy unidirectional composites. The patch is located on the crack in such a way that the crack tip is not covered. Magnified isochromatic fringes are obtained by using a projection microscope of magnification 50, converted into a polariscope. Irwin's method is used to compute stress intensity factors from photoelastic data. The reduction in stress intensity factors is presented in graphical form as a function of patch parameters, namely stiffness, location and length. An empirical equation connecting reduction in stress intensity factor and these patch parameters is presented.

On an analytical-numerical procedure for the analysis of cylindrical shells with arbitrarily oriented cracks

An analytical-numerical procedure for obtaining stress intensity factor solutions for an arbitrarily oriented crack in a long, thin circular cylindrical shell is presented. The method of analysis involves obtaining a series solution to the governing shell equation in terms of Mathieu and modified Mathieu functions by the method of separation of variables and satisfying the crack surface boundary conditions numerically using collocation. The solution is then transformed from elliptic coordinates to polar coordinates with crack tip as the origin through a Taylor series expansion and membrane and bending stress intensity factors are computed. Numerical results are presented and discussed for the pressure loading case.