Dynamic fracture mechanics analysis of failure mode transitions along weakened interfaces in elastic solids

Dynamic fracture mechanics theory was employed to analyze the crack deflection behavior of dynamic mode-I cracks propagating towards inclined weak planes/interfaces in otherwise homogenous elastic solids. When the incident mode-I crack reached the weak interface, it kinked out of its original plane and continued to propagate along the weak interface. The dynamic stress intensity factors and the non-singular T-stresses of the incident cracks were fitted, and then dynamic fracture mechanics concepts were used to obtain the stress intensity factors of the kinked cracks as functions of kinking angles and crack tip speeds. The T-stress of the incident crack has a small positive value but the crack path was quite stable. In order to validate fracture mechanics predictions, the theoretical photoelasticity fringe patterns of the kinked cracks were compared with the recorded experimental fringes. Moreover, the mode mixity of the kinked crack was found to depend on the kinking angle and the crack tip speed. A weak interface will lead to a high mode-II component and a fast crack tip speed of the kinked mixed-mode crack.     

Mode I cracks subjected to large T-stresses

There are several criteria for predicting brittle fracture in mode I and mixed mode loading. In this paper, the modified maximum tangential stress criterion originally proposed for mixed mode loading, is employed to study theoretically brittle fracture for mode I cracks. In particular, the effect of the non-singular term of stress, often known as the T-stress, on the angle of initiation of fracture and the onset of crack growth is explored. The T-stress component of the tangential stress vanishes along the crack line. Therefore, it is often postulated for linear elastic materials that the effect of T-stress on mode I brittle fracture can be ignored. However, it is shown here that the maximum tangential stress is no longer along the line of initial crack when the T-stress exceeds a critical value. Thus, a deviation in the angle of initiation of fracture can be expected for specimens having a large T-stress. It is shown that the deviation angle increases for larger values of T-stress. Theoretical results show that the apparent fracture toughness decreases significantly when a deviation in angle occurs. Earlier experimental results are used to corroborate the findings. The effect of large T-stresses is also explored for a crack specimen undergoing moderate scale yielding. The elastic-plastic investigation is conducted using finite element analysis. The finite element results reveal a similar deviation in the angle of maximum tangential stress for small to moderate scale yielding.     

Variation of elastic T-stresses along slightly wavy 3D crack fronts

The non-singular terms in the series expansion of the elastic crack-tip stress field, commonly referred to as the elastic T-stresses, play an important role in fracture mechanics in areas such as the stability of a crack path and the two-parameter characterization of elastic-plastic crack-tip deformation. In this paper, a first order perturbation analysis is performed to study some basic properties of the T-stress variation along a slightly wavy 3D crack front. The analysis employs important properties of Bueckner-Rice 3D weight function fields. Using the Boussinesq-Papkovitch potential representation for the mode I weight function field, it is shown that, for coplanar cracks in an infinite isotropic and homogeneous linear elastic body, the mean T-stress along an arbitrary crack front is independent of the shape and size of the crack. Further, a universal relation is discovered between the mean T-stress and the stress field at the same crack front location under the same loading but in the absence of a crack. First-order-accurate solutions are given for the T-stress variation along a slightly wavy crack front with nearly circular or straight confifurations. Specifically, cosine wave functions are adopted to describe smooth polygonal and slightly undulating planar crack shapes. The results indicate that T ₁₁, the 2D T-stress component acting normal to the crack front, increases with the curvature of the crack front as it bows out but T ₃₃, acting parallel to the crack front, decreases with the crack front curvature.     

Fractal effects of crack propagation on dynamic stress intensity factors and crack velocities

Crack extension paths are often irregular, producing rough fracture surfaces which have a fractal geometry. In this paper, crack tip motion along a fractal crack trace is analysed. A fractal kinking model of the crack extension path is established to describe irregular crack growth. A formula is derived to describe the effects of fractal crack propagation on the dynamic stress intensity factor and on crack velocity. The ratio of the dynamic stress intensity factor to the applied stress intensity factor K(L(D, t), V)/K(L(t), 0), is a function of apparent crack velocity V_o, microstructure parameter d/a (grain size/crack increment step length), fractal dimension D, and fractal kinking angle of crack extension path . For fractal crack propagation, the apparent (or measured) crack velocity V_o, cannot approach the Rayleigh wave speed Cr. Why V_o is significantly lower than Cr in dynamic fracture experiments can be explained by the effects of fractal crack propagation. The dynamic stress intensity factor and apparent crack velocity are strongly affected by the microstructure parameter (d/a), fractal dimension D, and fractal kinking angle of crack extension path . This is in good agreement with experimental findings.     

A theoretical note on mode-I crack branching and kinking

An energy-based fracture mode has been derived for the mode-I crack branching and kinking. The classic J_i-integral has been further explored by a new partial integral path and the analytical solution of the energy release rate for crack branching and kinking from a mode-I crack tip has been established. The crack branching/kinking angle has also been analytically derived. It shows that the Griffith’s theorem and conservation law can be applied to both mode-I crack extension and mode-I crack branching and kinking. The branching mechanism for quasi-static mode-I crack has been theoretically investigated. The branching toughness and the K-based criterion for crack branching have been defined. The crack branching phenomena predicted by the present model are in well agreement with the experimental observations reported in the literatures.     

On the full set of elastic T-stress terms of internal circular cracks under mixed-mode loading conditions

Analytical expressions for all non-singular stress terms of the asymptotic crack tip field, the so-called T-stresses of internal mixed mode circular (penny shaped) cracks in infinite homogeneous and isotropic elastic solids are addressed. To solve the problem the mixed mode crack problem is divided into sub-problems using the superposition method, and each sub-problem is then solved for the asymptotic stress field. Considering the expansion of the local stress field at the crack front, the elastic T-stress terms are derived for each sub-problem. The results are superimposed to give the analytical expressions of the so far missing elastic T-stresses of internal mixed mode penny shaped cracks.The effect of the T-stresses on the size and shape of the plastic zone at the crack tip is discussed, and analytical results are compared to the ones from finite element analyses, both for the T-stress components and the size of the plastic zone. For an accurate prediction of the plastic zone all the singular terms and the constant terms (T-stresses) in the stress expansion formulae should be considered. It is observed that negative T-stresses increase the size of the plastic zone, while positive ones reduce it.     

T-stress near the tips of a cruciform crack with unequal arms

The T-stress near the tips of a crack of cross shape embedded in an isotropic elastic solid is analyzed. The integral transform technique is employed to convert the associated boundary value problem to a system of singular integral equations. According to the stress difference method, T-stresses can be expressed as a sum of an integral involving crack opening displacement (COD) and applied loading at infinity. Obtained results indicate that, in addition to applied loading, T-stresses at the horizontal (vertical) crack tips depend on the COD of the vertical (horizontal) crack surface. COD plays a leading role in determining T-stresses, in particular for a cruciform crack of two crack-arm lengths of the same order. Moreover, T-stresses for a single-crack limiting case can be recovered from the present results as the length of one arm approaches zero. For a biaxial tension of the same magnitude, T-stresses are present for a cruciform crack, but absent for a single crack. Finally, for several cases of interest, T-stresses for a cruciform crack are evaluated and compared with those for a single crack, and the influence of the ratio of two crack-arm lengths b/a and the COD on the T-stress of a cruciform crack is examined.     

Maximum Average Tangential Stress Criterion for Prediction of the Crack Path

The concept of the average stress has been employed to propose the maximum average tangential stress (MATS) criterion for predicting the direction of fracture angle. This criterion states that a crack grows when the maximum average tangential stress in the fracture process zone ahead of the crack tip reaches its critical value and the crack growth direction coincides with the direction of the maximum average tangential stress along a constant radius around the crack tip. The tangential stress is described by the singular and non-singular (T-stress) terms in the Williams series solution. The predicted directions of fracture angle are consistent with the experimental data for the mixed mode I/II crack growth behavior of Guiting limestone.     

Effect of material plasticity on fatigue crack propagation under complex stress state

The maximum energy release rate criterion, i.e., G _max criterion, is commonly used for crack propagation analysis. However, this fracture criterion is based on the elastic macroscopic strength of materials. In the present investigation, a modification has been made to G _max criterion to implement the consideration of the plastic strain energy. This criterion is extended to study the fatigue crack growth characteristics of mixed mode cracks in steel pipes. To predict crack propagation due to fatigue loads, a new elasto-plastic energy model is presented. This new model includes the effects of material properties like strain hardening exponent n, yield strength σ_y and fracture toughness and stress intensity factor ranges. The results obtained are compared with those obtained using the commonly employed crack growth law and the experimental data.     

Numerical evaluation of the quarter-point crack tip element

This paper attempts to answer two commonly raised questions during the preparation of a finite element mesh, for the linear elastic fracture analysis of cracked structure: how to set up the finite element mesh around the crack tip, and what level of accuracy is to be expected from such a modelling. Two test problems, with known analytical expressions for their stress intensity factors, are analysed by the finite element method using the isoparametric quadratic singular element. The modified parameters were the order of integration, aspect ratio, number of elements surrounding the crack tip, use of transition elements, the singular element length over the total crack length, the symmetry of the mesh around the crack tip. Based on these analyses, a data base is created and various plots produced. The results are interpreted, the accuracy evaluated and recommendations drawn. Contrary to previous reports, it is found that the computed stress intensity factor (SIF) remains within engineering accuracy (10 per cent) throughout a large range of l/a (singular element length over crack length) for problems with a uniform non-singular stress distribution ahead of the crack tip (i.e. double edge notch), and l/a should be less than 0·1 for problems with a non-singular stress gradient (i.e three-point bend). Also, it is found that the best results are achieved by using at least four singular elements around the crack tip, with their internal angles around 45 degrees, and a reduced (2 × 2) numerical integration.     

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₂ 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.     

SINGULAR p-VERSION FINITE ELEMENTS FOR STRESS INTENSITY FACTOR COMPUTATIONS

The finite element analysis of linear elastic fracture mechanics problems is complicated by the presence of the singular and finite non-singular stress distributions in the crack tip region. The availability of a constant stress term in addition to the singular term in the standard h-version singular finite elements is insufficient to model the finite nonsingular stress zone. A p-version singular finite element capable of modelling the higher-order non-singular stress terms in addition to the singular term and the constant term is presented. The formulation for the displacement substitution technique for computing the stress intensity factors using singular p-version triangular finite elements is developed. Unlike the standard h-version formulation, the stress intensity factors computed using the p-version displacement substitution technique do not depend on the specific arrangement and length of the quarter point elements, and require simple mesh designs as well as fewer number of degrees of freedom. Numerical studies comparing the convergence of the stress intensity factors computed by the p-version method against other available alternatives such as the h-version method and the contour integral method are presented to demonstrate the effectiveness of the present developments. © 1997 by John Wiley & Sons, Ltd.     

T-stress effects on steady crack growth in a thin, ductile plate under small-scale yielding conditions: Three-dimensional modeling

The non-singular T-stress provides a first-order estimate of geometry and loading mode, e.g. tension vs. bending, effects on elastic–plastic, crack-front fields under mode I conditions. The T-stress has a pronounced effect on measured crack growth resistance curves for ductile metals – trends most computational models confirm using a two-dimensional setting. This work examines T-stress effects on three-dimensional (3D), elastic–plastic fields surrounding a steadily advancing crack for a moderately hardening material in the framework of a 3D, small-scale yielding boundary-layer model. A flat, straight crack front advances at a constant quasi-static rate under near invariant local and global mode I loading. The boundary-layer model has thickness B that defines the only geometric length-scale. The material flow properties and (local) toughness combine to limit the in-plane plastic-zone size during steady growth to at most a few multiples of the thickness (conditions obtainable, for example, in large, thin aluminum components). The computational model requires no crack growth criterion; rather, the crack front extends steadily at constant values of the plane-stress displacements imposed on the remote boundary for the specified far-field stress intensity factor and T-stress. The specific numerical results presented demonstrate similarity scaling of the 3D near-front stresses in terms of two non-dimensional loading parameters. The analyses reveal a strong effect of T-stress on key stress and strain quantities for low loading levels and less effect for higher loading levels, where much of the plastic zone experiences plane-stress conditions. To understand the combined effects of T-stress on stresses and plastic strain levels, normalized values from a simple void-growth model, computed over the crack plane for low loading, clearly reveal the tendency for crack-front tunneling, shear-lip formation near the outside surfaces, and a minimum steady-state fracture toughness for T = 0 loading.     

Partially stiffened elastic half-plane with an edge crack

A technique, using the Brazilian disk specimen, for measuring the fracture toughness of unidirectional fiber-reinforced composites, over the entire range of crack-tip mode mixities, was developed. The fracture toughness of a graphite/epoxy fiber-reinforced composite was measured, under both mode-I and mode-II loading conditions. We found that for certain material orientations the mode-II fracture toughness is substantially higher than the mode-I toughness. The complete dependence of the fracture toughness on the crack-tip mixity was determined for particular material orientations and the phenomenological fracture toughness curves were constructed. Using the Brazilian disk specimen, together with a hydraulic testing machine, the fracture toughness of the composite under moderate loading rates was measured. We observed that the mode-I fracture toughness was not sensitive to the loading rate at the crack tip, K, while the mode-II ‘dynamic’ fracture toughness increased approximately 50 percent over the quasi-static fracture toughness. A qualitative explanation of the dependency of fracture toughness on crack-tip loading rate is discussed. Finally, a mechanical fracture criterion, at the microscopic level, which governs the crack initiation under mixed-mode loading conditions is presented; these theoretical predictions closely follow the trend of experimental measurements. This revised version was published online in August 2006 with corrections to the Cover Date.     

Mechanical properties of the solid Li-ion conducting electrolyte: Li0.33La0.57TiO3

Li_0.33La_0.57TiO₃ (LLTO) is a potential Li-ion conducting membrane for use in aqueous Li-air batteries. To be in this configuration its mechanical properties must be determined. Dense LLTO was prepared using a solid-state (SS) or sol–gel (SG) procedure and was hot-pressed to yield a high relative density material (>95 %). Young’s modulus, hardness, and fracture toughness of the LLTO-SS and sol–gel LLTO-SG materials was determined and compared to other solid Li-ion conducting electrolytes. The Young’s modulus for LLTO-SG and LLTO-SS was 186 ± 4 and 200 ± 3 GPa, respectively. The Vickers hardness of LLTO-SG and LLTO-SS was 9.7 ± 0.7 and 9.2 ± 0.2 GPa, respectively. The fracture toughness, K _IC, of both LLTO-SG and LLTO-SS was ~1 MPa m^1/2; the fracture toughness of LLTO-SG was slightly higher than that of LLTO-SS. Both LLTO-SG and LLTO-SS have a Young’s modulus and hardness greater than the other possible solid Li-ion conducting membranes; Li₇La₃Zr₂O₁₂ and Li_1+x+y Al _x Ti_2−x Si _y P_3−y O₁₂. Based on modulus and hardness hot-pressed LLTO exhibits sufficient mechanical integrity to be used as a solid Li-ion conducting membrane in aqueous Li-air batteries but, its fracture toughness needs to be improved without degrading its ionic conductivity.     

Frictional contact induced crack initiation in incompressible substrate

Surface crack initiation in an incompressible substrate induced through frictional fretting contact is analyzed using an energy-based fracture mechanics model. A closed-form energy release rate for surface crack initiation at the contact boundary has been derived with the crack growth angle determined by the mixed mode singular stress field at the contact edge. The driving forces in the form of J_i-integral, the critical energy release rate and the critical load for crack initiation from the crack free surface have been formulated. The relations between the friction coefficient and crack initiation angle, critical load have been specified.     

Thermoelastic analysis for an opening crack in an orthotropic material

The thermoelastic analysis of an opening crack embedded in an orthotropic material is made under applied uniform heat flux and mechanical loadings. To simulate the case of an opening crack filled with a medium, a thermal-medium crack model is proposed. The thermally permeable and impermeable cracks are the limiting ones of the proposed thermal-medium one. The crack-tip thermoelastic fields induced by a crack in an orthotropic material are determined in closed forms. The elastic T-stress can be also obtained explicitly. The effects of applied mechanical loadings and the thermal conductivity of crack interior on the heat flux at the crack surfaces and the mode-II stress intensity factor are investigated through numerical computations. The obtained results reveal that an increase of the thermal conductivity of crack interior decreases the mode-II stress intensity factor. And when an applied mechanical loading is increasing, the mode-II stress intensity factor is rising.     

On the full set of elastic T-stress terms of internal elliptical cracks under mixed-mode loading condition

Analytical expressions for the elastic constant stress terms of the asymptotic field, the so called T-stresses, for internal mixed-mode elliptical cracks in infinite homogeneous and isotropic elastic solids are addressed. To solve the problem the mixed-mode crack problem is divided into sub-problems using the superposition method, and each sub-problem is then solved for the asymptotic stress field. Considering the expansion of the local stress field at the crack front, the elastic T-stress terms are derived for each sub-problem. The results are superimposed to give the analytical expressions of the so far missing elastic T-stresses for mixed-mode elliptical cracks.The effect of the T-stresses on the size and shape of the plastic zone at the crack tip is discussed, and analytical results are compared to the ones from finite element analyses, both for the T-stress components and the size of the plastic zone. For an accurate prediction of the plastic zone all singular and constant terms (T-stresses) in the stress expansion formulae should be considered. It is observed that negative T-stresses increase the size of the plastic zone, while positive ones reduce it.     

Effect of particles on interlaminar crack growth in cross-plied carbon-fibre/epoxy laminates

This paper discusses the effect of particulate additions on the mode-I and mode-II interlaminar fracture toughness of a cross-plied, carbon-fibre-reinforced, epoxy-resin laminate. Particles of graphite, silicon carbide and polyethylene were mixed with the epoxy resin prior to laminating with woven carbon-fibre cloth. Tests have been performed on double cantilever beam (DCB) and end-notched-flexure (ENF) specimens to obtain the critical-strain energy-release rates,G _IC andG _IIC, for the laminates with and without particulate additions. The dependences of the values ofG _IC andG _IIC on the crack length are also considered. The results indicate that the interlaminar-fracture-toughness (mode-I and mode-II) values of the CFRP laminate increase with increases in the particle content up to about 3%, and thereafter they decrease with further increases in the particle content. This was thought to be due to an increase in multiple-crack formation in the middle region of the cracked-plate samples. Furthermore, mode-I tests indicate that the propagation values ofG _IC are dependent on the crack length.     

On Influence of Non-Singular Stress States on Brittle Fracture

In the case of sufficiently brittle material the use of stress intensity factor as a fracture parameter alone is well justified within the Linear-Elastic Fracture Mechanics. This is because the singular stress field associated with the stress intensity factor is dominant near the crack tip. However, there are numerous experimental evidences that the critical stress intensity factor to cause fracture initiation (or fracture toughness) can be affected by the specimen geometry as well as loading conditions. To address this issue a number of twoparameter criteria have been proposed in the past, which often utilise non-singular terms of the classical asymptotic expansion of the stress field near the crack tip. Therefore, there is a problem of the selection of an appropriate parameter in addition to the stress intensity factor, which could account for various effects induced by the specimen geometry and loading on initiation of brittle fracture. This short paper demonstrates that brittle fracture conditions can be successfully predicted with various two-parameter criteria, and there are no clear advantages in the use of T-stress as the additional parameter in fracture criterion, in comparison with the next non-singular term, A ₁, of the asymptotic expansion.