Families of asymptotic tensile crack tip stress fields in elastic-perfectly plastic single crystals

In this work, two families of asymptotic near-tip stress fields are constructed in an elastic-ideally plastic FCC single crystal under mode I plane strain conditions. A crack is taken to lie on the (010) plane and its front is aligned along the [(1) over bar 01] direction. Finite element analysis is first used to systematically examine the stress distributions corresponding to different constraint levels. The general framework developed by Rice (Mech Mater 6:317-335, 1987) and Drugan (J Mech Phys Solids 49:2155-2176, 2001) is then adopted to generate low triaxiality solutions by introducing an elastic sector near the crack tip. The two families of stress fields are parameterized by the normalized opening stress (tau(A)(22)/tau(o)) prevailing in the plastic sector in front of the tip and by the coordinates of a point where elastic unloading commences in stress space. It is found that the angular stress variations obtained from the analytical solutions show good agreement with finite element analysis.

Mixed mode (I and II) crack tip fields in bulk metallic glasses

Stationary crack tip fields in bulk metallic glasses under mixed mode (I and II) loading are studied through detailed finite element simulations assuming plane strain, small scale yielding conditions. The influence of internal friction or pressure sensitivity on the plastic zones. notch deformation, stress and plastic strain fields is examined for different mode mixities. Under mixed mode loading, the notch deforms into a shape such that one part of its surface sharpens while the other part blunts. Increase in mode If component of loading dramatically enhances the normalized plastic zone size, lowers the stresses but significantly elevates the plastic strain levels near the notch tip. Higher internal friction reduces the peak tangential stress but increases the plastic strain and stretching near the blunted part of the notch. The simulated shear bands are straight and extend over a long distance ahead of the notch tip under mode II dominant loading. The possible variations of fracture toughness with mode mixity corresponding to failure by brittle micro-cracking and ductile shear banding are predicted employing two simple fracture criteria. The salient results from finite element simulations are validated by comparison with those from mixed mode (I and II) fracture experiments on a Zr-based bulk metallic glass.

Stationary crack tip fields in elastic-plastic solids: an overview of recent numerical simulations

In this paper, an overview of some recent numerical simulations of stationary crack tip fields in elastic-plastic solids is presented. First, asymptotic analyses carried out within the framework of 2D plane strain or plane stress conditions in both pressure insensitive and pressure sensitive plastic solids are reviewed. This is followed by discussion of salient results obtained from recent computational studies. These pertain to 3D characteristics of elastic-plastic near-front fields under mixed mode loading, mechanics of fracture and simulation of near-tip shear banding process of amorphous alloys and influence of crack tip constraint on the structure of near-tip fields in ductile single crystals. These results serve to illustrate several important features associated with stress and strain distributions near the crack tip and provide the foundation for understanding the operative failure mechanisms. The paper concludes by highlighting some of the future prospects for this field of study.

A numerical study of crack tip constraint in ductile single crystals

In this work, the effect of crack tip constraint on near-tip stress and deformation fields in a ductile FCC single crystal is studied under mode I, plane strain conditions. To this end, modified boundary layer simulations within crystal plasticity framework are performed, neglecting elastic anisotropy. The first and second terms of the isotropic elastic crack tip field, which are governed by the stress intensity factor K and T-stress, are prescribed as remote boundary conditions and solutions pertaining to different levels of T-stress are generated. It is found that the near-tip deformation field, especially, the development of kink or slip shear bands, is sensitive to the constraint level. The stress distribution and the size and shape of the plastic zone near the crack tip are also strongly influenced by the level of T-stress, with progressive loss of crack tip constraint occurring as T-stress becomes more negative. A family of near-tip fields is obtained which are characterized by two terms (such as K and T or J and a constraint parameter Q) as in isotropic plastic solids.

Crack Tip Fields in a Single Edge Notched Aluminum Single Crystal Specimen

We report a combined experimental and computational study of a low constraint aluminum single crystal fracture geometry and investigate the near-tip stress and strain fields. To this end, a single edge notched tensile (SENT) specimen is considered. A notch, with a radius of 50 µm, is taken to lie in the (010) plane and its front is aligned along the [101] direction. Experiments are conducted by subjecting the specimen to tensile loading using a special fixture inside a scanning electron microscope chamber. Both SEM micrographs and electron back-scattered diffraction (EBSD) maps are obtained from the near-tip region. The experiments are complemented by performing 3D and 2D plane strain finite element simulations within a continuum crystal plasticity framework assuming an isotropic hardening response characterized by the Pierce–Asaro–Needleman model. The simulations show a distinct slip band forming at about 55 deg with respect to the notch line corresponding to slip on (11-bar 1)[011] system, which corroborates well with experimental data. Furthermore, two kink bands occur at about 45 deg and 90 deg with respect to the notch line within which large rotations in the crystal orientation take place. These predictions are in good agreement with the EBSD observations. Finally, the near-tip angular variations of the 3D stress and plastic strain fields in the low constraint SENT fracture geometry are examined in detail.

Numerical simulations of crack tip fields in polycrystalline plastic solids

In this paper, modes I and II crack tip fields in polycrystalline plastic solids are studied under plane strain, small scale yielding conditions. Two different initial textures of an Al–Mg alloy, viz., continuous cast AA5754 sheets in the recrystallized and cold rolled conditions, are considered. The former is nearly-isotropic, while the latter displays distinct anisotropy. Finite element simulations are performed by employing crystal plasticity constitutive equations along with a Taylor-type homogenization as well as by using the Hill quadratic yield theory. It is found that significant texture evolution occurs close to the notch tip which profoundly influences the stress and plastic strain distributions. Also, the cold rolling texture gives rise to higher magnitude of plastic strain near the tip.

Numerical simulations of crack tip fields in polycrystalline plastic solids

In this paper, modes I and II crack tip fields in polycrystalline plastic solids are studied under plane strain, small scale yielding conditions. Two different initial textures of an Al-Mg alloy, viz.,continuous cast AA5754 sheets in the recrystallized and cold rolled conditions, are considered. The former is nearly-isotropic, while the latter displays distinct anisotropy. Finite element simulations are performed by employing crystal plasticity constitutive equations along with a Taylor-type homogenization as well as by using the Hill quadratic yield theory. It is found that significant texture evolution occurs close to the notch tip which profoundly influences the stress and plastic strain distributions. Also, the cold rolling texture gives rise to higher magnitude of plastic strain near the tip. (C) 2010 Elsevier Ltd. All rights reserved.

Influence of crack tip constraint on void growth in ductile FCC single crystals

Effect of constraint (stress triaxiality) on void growth near a notch tip in a FCC single crystal is investigated. Finite element simulations within the modified boundary layer framework are conducted using crystal plasticity constitutive equations and neglecting elastic anisotropy. Displacement boundary conditions based on model, elastic, two term K-T field are applied on the outer boundary of a large circular domain. A pre-nucleated void is considered ahead of a stationary notch tip. The interaction between the notch tip and the void is studied under different constraints (T-stress levels) and crystal orientations. It is found that negative T-stress retards the mechanisms of ductile fracture. However, the extent of retardation depends on the crystal orientation. Further, it is found that there exists a particular orientation which delays the ductile fracture processes and hence can potentially improve ductility. This optimal orientation depends on the constraint level. (C) 2010 Published by Elsevier B.V.

Effect of lattice orientation on crack tip constraint in ductile single crystals

In this work, the effect of lattice orientation on the fields prevailing near a notch tip is investigated pertaining to various constraint levels in FCC single crystals. A modified boundary layer formulation is employed and numerical solutions under mode I, plane strain conditions are generated by assuming an elastic-perfectly plastic FCC single crystal. The analysis is carried out corresponding to different lattice orientations with respect to the notch line. It is found that the near-tip deformation field, especially the development of kink or slip shear bands is sensitive to the constraint level. The stress distribution and the size and shape of the plastic zone near the notch tip are also strongly influenced by the level of T-stress. The present results clearly establish that ductile single crystal fracture geometries would progressively lose crack tip constraint as the T-stress becomes more negative irrespective of lattice orientation. Also, the near-tip field for a range of constraint levels can be characterized by two-parameters such as K-T or J-Q as in isotropic plastic solids.

Quasi-static crack tip fields in rate-sensitive FCC single crystals

In this work, the effects of loading rate, material rate sensitivity and constraint level on quasi-static crack tip fields in a FCC single crystal are studied. Finite element simulations are performed within a mode I, plane strain modified boundary layer framework by prescribing the two term (K-T) elastic crack tip field as remote boundary conditions. The material is assumed to obey a rate-dependent crystal plasticity theory. The orientation of the single crystal is chosen so that the crack surface coincides with the crystallographic (010) plane and the crack front lies along 101] direction. Solutions corresponding to different stress intensity rates K., T-stress values and strain rate exponents m are obtained. The results show that the stress levels ahead of the crack tip increase with K. which is accompanied by gradual shrinking of the plastic zone size. However, the nature of the shear band patterns around the crack tip is not affected by the loading rate. Further, it is found that while positive T-stress enhances the opening and hydrostatic stress levels ahead of crack tip, they are considerably reduced with imposition of negative T-stress. Also, negative T-stress promotes formation of shear bands in the forward sector ahead of the crack tip and suppresses them behind the tip.

Universal crack closure integral for SIF estimation

Numerical analysis of cracked structures often involves numerical estimation of stress intensity factors (SIFs) at a crack tip/front. A newly developed formulation called universal crack closure integral (UCCI) for the evaluation of potential energy release rates (PERRs) and the corresponding SIFs is presented in this paper. Unlike the existing element dedicated forms of crack closure integrals (MCCI, VCCI) with application limited to finite element analysis, this new numerical SIF/PERR estimation technique is independent of the basic stress analysis procedure, making it universally applicable. The second merit of this procedure is that it avoids the generally error-producing zones close to the crack tip/front singularity. The UCCI procedure, based on Irwin's original CCI, is formulated and explored using a simple 2D problem of a straight crack in an infinite sheet. It is then applied to some three-dimensional crack geometries with the stresses and displacements obtained from a boundary element program.

Biaxial-load effects on isopachics for a crack at the interface of dissimilar media

The plane problem of two dissimilar materials, bonded together and containing a crack along their common interface, which were subjected to a biaxial load at infinity, is examined by giving a closed-form expression for the first stress invariant of the normal stresses, which is equally valid everywhere, near to, and far from, the crack-tip region. This exact expression for the first-stress invariant is compared by constructing the respective isopachic-fringe patterns, to the approximate expression with non-singular terms, due to the biaxiality factor, for the same quantity. Significant differences between respective isopachic-patterns were found and their dependence on the elastic properties of both materials and the applied loads was demonstrated. The relative errors between the computedK I - andK II -components by using the approximate expression for the first stress-invariant and the accurate one, derived from closed-form solution along either isopachic-fringes or along circles and radii from the crack-tip have been given, indicating in some cases large discrepancies between exact and approximate solutions.

New a posteriori error estimator and adaptive mesh refinement strategy for 2-D crack problems

A posteriori error estimation and adaptive refinement technique for fracture analysis of 2-D/3-D crack problems is the state-of-the-art. The objective of the present paper is to propose a new a posteriori error estimator based on strain energy release rate (SERR) or stress intensity factor (SIF) at the crack tip region and to use this along with the stress based error estimator available in the literature for the region away from the crack tip. The proposed a posteriori error estimator is called the K-S error estimator. Further, an adaptive mesh refinement (h-) strategy which can be used with K-S error estimator has been proposed for fracture analysis of 2-D crack problems. The performance of the proposed a posteriori error estimator and the h-adaptive refinement strategy have been demonstrated by employing the 4-noded, 8-noded and 9-noded plane stress finite elements. The proposed error estimator together with the h-adaptive refinement strategy will facilitate automation of fracture analysis process to provide reliable solutions.

Observation of crack closure using a crack mouth opening displacement gauge

A crack mouth opening displacement (CMOD) gauge has been designed to estimate the crack length as well as the crack opening stress level in an automated fatigue crack propagation test programme. The CMOD gauge accurately predicts crack tip opening levels in K-controlled tests with constant K, K increasing or K decreasing. In all three K-controlled tests with a single 100% overload cycle, the CMOD gauge does not measure the crack tip opening stress level over a large range of crack lengths after the application of the overload. The CMOD gauge measures the stress level at which the overload plastic zone site opens. Caution should thus be exercised in using the crack opening stress level, estimated by the CMOD gauge, to explain fatigue crack propagation under arbitrary load sequences from the viewpoint of crack closure phenomena.