A boundary layer framework considering material gradation effects

This paper describes the development and application of a novel modified boundary layer (MBL) model for graded nonhomogeneous materials, e.g. functionally graded materials (FGMs). The proposed model is based on a middle-crack tension, M(T), specimen with traction boundary conditions applied to the top and lateral edges of the model. Finite element analyses are performed using WARP3D, a fracture mechanics research finite element code, which incorporates elements with graded elastic and plastic properties. Elastic crack-tip fields obtained from the proposed MBL model show excellent agreement with those obtained from full models of the cracked component for homogeneous and graded nonhomogeneous materials. The K-T dominance of FGMs is investigated by comparing the actual stress fields with the asymptotic stress fields (the Williams’ solution). The examples investigated in the present study consider a crack parallel to the material gradient. Results of the present study provide insight into the K-T dominance of FGMs and also show the range of applicability of the proposed MBL model. The MBL model is applied to analyze the elastic-plastic crack-tip response of a Ti/TiB FGM SE(T) specimen. The numerical results demonstrate that the proposed MBL model captures the effect of T-stress on plastic zone size and shape, constraint effects, etc. for such configurations.     

Mixed mode asymptotic crack tip fields in orthotropic materials: Derivation and range of dominance

In the first part of this work the asymptotic stress field distribution surrounding a crack in a generally orthotropic solid (i.e. one in which the material and loading axes do not coincide) is derived. The resulting crack tip stress intensity factors are also related to the energy release rate of the cracked solid. In any physical situation, however, the range of dominance of these asymptotic fields will be limited, and will depend on specific geometry and loading parameters. Thus, the second part of this work deals with the range of dominance of the derived stress fields in edge cracked bending loaded fiber reinforced composite plates. The lower and upper limits of the range of dominance of the solution are respectively determined from three-dimensional and two-dimensional full field finite element solutions. The comparison of full field solutions with the asymptotic result provides information on the latter's range of dominance. In all cases the effects of mixed mode loading are considered.     

On higher order parameters in cracked composite plates under far‐field pure shear

This paper presents the analytical solution of the crack tip fields as well as the crack parameters in an infinitely large composite plate with a central crack subjected to pure shear loading. To this end, the complex variable method is employed to formulate an asymptotic solution for the crack tip fields in an anisotropic plane. Using a stress‐based definition of the crack tip modes of loading, only the mode II crack parameters are found to be non‐zero under pure shear load. Special focus is given to the determination of the higher order parameters of the crack tip asymptotic field, particularly the first non‐singular term, ie, the T‐stress. Unlike the isotropic materials, in which the T‐stress is zero under pure shear, it is found that the T‐stress is non‐zero for the case of anisotropic materials, being the only material‐dependent crack tip stress parameter. The veracity of our exact crack tip fields is assessed and verified through a comparison made with respect to the finite element (FE) solution. Finally, we demonstrate the significance of the T‐stress on stresses near the crack tip in composite plates under pure shear loads.     

Analysis of elastoplastic sharp notches

In this paper the elastoplastic solutions with higher-order terms for apex V-notches in power-law hardening materials have been discussed. Two-term expansions of the plane strain and the plane stress solutions have been obtained. It has been shown that the leading-order singularity approaches the value for a crack when the notch angle is not too large. In plane strain cases the elasticity does not enter the second-order solutions when the notch opening angle is too small. For a large notch angle, the two-term expansions of the plane strain near-tip fields are described by a single amplitude parameter. The plane stress solutions generally contain the elasticity terms. The boundary layer formulations based on the small-strain plasticity theory confirm that a dominance zone exists ahead of the notch tip. Finite element results give good agreement to the asymptotic solutions under both plane strain and plane stress conditions. The second-order terms cannot improve the predictions significantly. The near-tip fields are dominated by a single parameter. Finite element calculations under the finite strain J ₂-flow plasticity theory revealed that the finite strains can only affect local characterization of the asymptotic solution. The asymptotic solution has a large dominance zone around the notch tip. For an apex notch bounded to a rigid substrate the leading-order singularity falls with the notch angle significantly more slowly than in the homogeneous material. It vanishes at the notch angle about 135° for all power-hardening exponents. The elasticity effects enter the second-order solutions when the notch angle becomes large enough. The tip fields are characterized by the hydrostatic stress and the shear stress ahead of the notch.     

梯度复合材料裂纹扩展路径和起裂载荷的有限元分析

为了模拟功能梯度材料(FGM)在工程应用中可能会出现的断裂问题并计算相应的开裂载荷,通过编写用户自定义UEL子程序将梯度扩展单元嵌入到ABAQUS软件中模拟功能梯度材料的物理场,并编写交互能量积分后处理子程序计算裂纹尖端的混合模式应力强度因子(SIF),采用最大周向应力准则编写子程序计算裂纹的偏转角,并模拟了裂纹扩展路径,计算了裂纹的起裂载荷。讨论了材料梯度参数对裂纹扩展路径以及起裂载荷的影响规律。通过与均匀材料的对比,验证了功能梯度材料断裂性能的优越性。研究表明:外载平行于梯度方向时,垂直梯度方向的初始裂纹朝着等效弹性模量小的方向扩展,且偏转角在梯度指数线性时出现峰值,并随着组分弹性模量比的增加而变大;当外载和初始裂纹均平行于梯度方向时,材料等效弹性模量和断裂韧性的增加或者梯度指数的减小都导致起裂载荷变大。     

XFEM for direct evaluation of mixed mode SIFs in homogeneous and bi‐materials

The extended finite element method (XFEM) is improved to directly evaluate mixed mode stress intensity factors (SIFs) without extra post‐processing, for homogeneous materials as well as for bimaterials. This is achieved by enriching the finite element (FE) approximation of the nodes surrounding the crack tip with not only the first term but also the higher order terms of the crack tip asymptotic field using a partition of unity method (PUM). The crack faces behind the tip(s) are modelled independently of the mesh by displacement jump functions. The additional coefficients corresponding to the enrichments at the nodes of the elements surrounding the crack tip are forced to be equal by a penalty function method, thus ensuring that the displacement approximations reduce to the actual asymptotic fields adjacent to the crack tip. The numerical results so obtained are in excellent agreement with analytical and numerical results available in the literature. Copyright © 2004 John Wiley & Sons, Ltd.     

Analytic and numerical studies on mode I and mode II fracture in elastic-plastic materials with strain gradient effects

This paper presents a study on fracture of materials at microscale (∼1 μm) by the strain gradient theory (Fleck and Hutchinson, 1993; Fleck et al., 1994). For remotely imposed classical K fields, the full-field solutions are obtained analytically or numerically for elastic and elastic-plastic materials with strain gradient effects. The analytical elastic full-field solution shows that stresses ahead of a crack tip are significantly higher than their counterparts in the classical K fields. The sizes of dominance zones for mode I and mode II near-tip asymptotic fields are 0.3l and 0.5l,while strain gradient effects are observed within land 2l to the crack tip, respectively, where l is the intrinsic material length in strain gradient theory and is on the order of microns in strain gradient plasticity (Fleck et al., 1994; Nix and Gao, 1998; Stolken and Evans, 1997). The Dugdale–Barenblatt type plasticity model is obtained to provide an estimation of plastic zone size for mode II fracture in materials with strain grain effects. The finite element method is used to investigate the small-scale-yielding solution for an elastic-power law hardening solid. It is found that the size of the dominance zone for the near-tip asymptotic field is the intrinsic material lengthl. For mode II fracture under the small-scale-yielding condition, transition from the remote classical K _IIfield to the near-tip asymptotic field in strain gradient plasticity goes through the HRR field only when K _IIis relatively large such that the plastic zone size is much larger than the intrinsic material length l. For mode I fracture under small-scale-yielding condition, however, transition from the remote classical K _I field to the near-tip asymptotic field in strain gradient plasticity does not go through the HRR field, but via a plastic zone. This revised version was published online in July 2006 with corrections to the Cover Date.     

An Interaction Integral Method for Computing Fracture Parameters in Functionally Graded Magnetoelectroelastic Composites

A contour integral method is developed for the computation of stress intensity, electric and magnetic intensity factors for cracks in continuously nonhomogeneous magnetoelectroelastic solids under a transient dynamic load. It is shown that the asymptotic fields in the crack-tip vicinity in a continuously nonhomogeneos medium are the same as in a homogeneous one. A meshless method based on the local Petrov-Galerkin approach is applied for the computation of the physical fields occurring in the contour integral expressions of intensity factors. A unit step function is used as the test functions in the local weak-form. This leads to local integral equations (LIEs) involving only contour-integrals on the surfaces of subdomains. The moving least-squares (MLS) method is adopted for approximating the physical quantities in the LIEs. The accuracy of the present method for computing the stress intensity factors (SIF), electrical displacement intensity factors (EDIF) and magnetic induction intensity factors (MIIF) are discussed by comparison with numerical solutions for homogeneous materials.     

Crack tip fields in ductile crystals

Results on the asymptotic analysis of crack tip fields in elastic-plastic single crystals are presented and some preliminary results of finite element solutions for cracked solids of this type are summarized. In the cases studied, involving plane strain tensile and anti-plane shear cracks in ideally plastic f c c and b c c crystals, analyzed within conventional small displacement gradient assumptions, the asymptotic analyses reveal striking discontinuous fields at the crack tip.For the stationary crack the stress state is found to be locally uniform in each of a family of angular sectors at the crack tip, but to jump discontinuously at sector boundaries, which are also the surfaces of shear discontinuities in the displacement field. For the quasi-statically growing crack the stress state is fully continuous from one near-tip angular sector to the next, but now some of the sectors involve elastic unloading from, and reloading to, a yielded state, and shear discontinuities of the velocity field develop at sector boundaries. In an anti-plane case studied, inclusion of inertial terms for (dynamically) growing cracks restores a discontinuous stress field at the tip which moves through the material as an elastic-plastic shock wave. For high symmetry crack orientations relative to the crystal, the discontinuity surfaces are sometimes coincident with the active crystal slip planes, but as often lie perpendicular to the family of active slip planes so that the discontinuities correspond to a kinking mode of shear.The finite element studies so far attempted, simulating the ideally plastic material model in a small displacement gradient type program, appear to be consistent with the asymptotic analyses. Small scale yielding solutions confirm the expected discontinuities, within limits of mesh resolution, of displacement for a stationary crack and of velocity for quasi-static growth. Further, the discontinuities apparently extend well into the near-tip plastic zone. A finite element formulation suitable for arbitrary deformation has been used to solve for the plane strain tension of a Taylor-hardening crystal panel containing, a center crack with an initially rounded tip. This shows effects due to lattice rotation, which distinguishes the regular versus kinking shear modes of crack tip relaxation. and holds promise for exploring the mechanics of crack opening at the tip.     

K‐dominance of static crack tip in functionally gradient materials

K‐dominance of static crack tip in functionally gradient materials (FGMs) with a crack oriented along the direction of the elastic gradient is studied through coherent gradient sensing (CGS), digital speckle correlation method (DSCM) and finite element method (FEM). In the direction of crack propagation, the shear modulus has a linear variation with constant mass density and Poisson's ratio. First, the CGS and DSCM governing equations related to the measurements and the elastic solutions at mode I crack in FGMs are obtained in terms of the stress intensity factor, material constants and graded index. Secondly, two kinds of FGMs specimens and one homogenous specimen are prepared to observe the influences of the property variation on the K‐dominance. Then, CGS and DSCM experiments using three‐point‐bending of FGMs and homogenous beams are performed. Thirdly, based on the results of the experiments, the stress intensity factors of three kinds of specimens are calculated by CGS and DSCM. Meanwhile, the stress intensity factors are obtained by FEM. Finally, comparing the results from CGS, DSCM and FEM, the K‐dominance of mode‐I static crack tip in FGMs is discussed in detail.     

A finite element investigation of quasi-static and dynamic asymptotic crack-tip fields in hardening elastic-plastic solids under plane stress

The asymptotic structures of crack-tip stress and deformation fields are investigated numerically for quasi-static and dynamic crack growth in isotropic linear hardening elastic-plastic solids under mode I, plane stress, and small-scale yielding conditions. An Eulerian type finite element scheme is employed. The materials are assumed to obey the von Mises yield criterion and the associated flow rule. The ratio of the crack-tip plastic zone size to that of the element nearest to the crack tip is of the order of 1.6 × 10⁴. The results of this study strongly suggest the existence of crack-tip stress and strain singularities of the type r ^s (s < 0) at r=0, where r is the distance to the crack tip, which confirms the asymptotic solutions of variable-separable type by Amazigo and Hutchinson [1] and Ponte Castañeda [2] for quasi-static crack growth, and by Achenbach, Kanninen and Popelar [3] for dynamic crack propagation. Both the values of the parameter s and the angular stress and velocity field variations from the present full-field finite element analysis agree very well with those from the analytical solutions. It is found that the dominance zone of the r ^s-singularity is quite large compared to the size of the crack-tip active plastic zone. The effects of hardening and inertia on the crack-tip fields as well as on the shape and size of the crack-tip active plastic zone are also studied in detail. It is discovered that as the level of hardening decreases and the crack propagation speed increases, a secondary yield zone emerges along the crack flank, and kinks in stress and velocity angular variations begin to develop. This dynamic phenomenon observable only for rapid crack growth and for low hardening materials may explain the numerical difficulties, in obtaining solutions for such cases, encountered by Achenbach et al. who, in their asymptotic analysis, neglected the existence of reverse yielding zones along the crack surfaces.     

Analytical forms of higher-order asymptotic elastic-plastic crack-tip fields in a linear hardening material under antiplane shear

An analytical study of the higher-order asymptotic solutions of the stress and strain fields near the traction-free crack tip under antiplane shear in a linear hardening material is investigated. The results show that every term of the asymptotic fields is controlled by both elasticity and plasticity and all the higher-order asymptotic fields are governed by linear nonhomogeneous equations. The first four term solutions are presented analytically and the first four terms are described by two independent parameters J and K ₂. The amplitude of the second order term solution is only dependent on the material properties, but independent of loading and geometry. This paper focuses on the case with traction-free crack surface boundary conditions. The effects of different crack surface boundary conditions, such as clamped and mixed surfaces, on the crack-tip fields are also presented. Comparison of multi-term solution with leading term solution, and finite element solution in an infinite strip with semi-infinite crack under constant displacements along the edges is provided.     

A domain-independent interaction integral for linear elastic fracture analysis of micropolar materials

This paper establishes a domain-independent interaction integral (DII-integral) for linear elastic fracture mechanics of micropolar elastic solids. The DII-integral has three amazing features that make it effective for solving the fracture parameters of complex micropolar materials. The first one is that the DII-integral can decouple the stress intensity factors (SIFs) and couple stress intensity factors (CSIFs) both of which are the key fracture parameters charactering the crack-tip asymptotic singular fields. In details, the DII-integral is derived from the J-integral by superimposing an actual field and an auxiliary field. By assigning the fracture parameters in the auxiliary field with different values, the SIFs and CSIFs of different crack opening modes can be obtained separately through the DII-integral. The second important feature is that the DII-integral is domain-independent for material nonhomogeneity and discontinuity. Thanks to this feature, the DII-integral becomes extremely effective for the micropolar materials with arbitrary nonhomogeneous properties or complex interfaces. The third feature is that the DII-integral does not contain any derivatives of material properties, which feature facilitate the practical implementation of the DII-integral on complex micropolar materials. Finally, the DII-integral combined with the extended finite element method (XFEM) is employed to solve four representative crack problems and the results show good validity of the DII-integral for complex micropolar materials.     

Three-dimensional analysis of interface cracks in rubber materials

Three-dimensional and plane stress finite element analyses were carried out to investigate the stress fields and fracture parameters of interface cracks in rubber materials. The tearing energy computations for any arbitrarily shaped 3D crack front in dissimilar materials were obtained using the virtual crack extension method. The finite element results obtained are validated against existing alternate solutions and experimental data for cracks in homogeneous as well as bimaterial cases. The effects of different rubber material models, tearing energy distributions, crack extension angles, and three-dimensional regions at the crack tip for interface cracks are presented and discussed. It is shown that nonlinear materials have larger three-dimensional effects near the interface crack front, and that this effect increases as the material mismatch increases.     

Crack kinking under high pressure in an elastic-plastic material

Directional crack growth criteria in compressed elastic–plastic materials are considered. The conditions at the crack tip are evaluated for a straight stationary crack. Remote load is a combined hydrostatic stress and pure shear, applied via a boundary layer assuming small scale yielding. Strains and deformations are assumed to be small. Different candidates for crack path criteria are examined. Maximum non-negative hoop stress to judge the risk of mode I and maximum shear stress for mode II extension of the crack are examined in some detail. Crack surfaces in contact are assumed to develop Coulumb friction from the very beginning. Hence, a condition of slip occurs throughout the crack faces. The plane in which the crack extends is calculated using a finite element method. Slip-line solutions are derived for comparison with the numerically computed asymptotic field. An excellent agreement between numerical and analytical solutions is found. The agreement is good in the region from the crack tip to around halfway to the elastic–plastic boundary. The relation between friction stress and yield stress is varied. The crack is found to extend in a direction straight ahead in shear mode for sufficiently high compressive pressure. At a limit pressure a kink is formed at a finite angle to the crack plane. For lower pressures the crack extends via a kink forming an angle to the parent crack plane that increases with decreasing pressure.     

Estimation of fracture parameters and stress field for edge cracks in finite elastically graded plates using boundary collocation

Summary The fracture parameters, stress intensity factor and T-stress are obtained for edge cracks aligned along the gradient in finite size elastically graded plates using the technique of boundary collocation. A scheme for extending the recently derived crack tip stress field for elastically graded materials is proposed. Using this extended stress field, the fracture parameters are evaluated for edge cracks subjected to far field tension and three point bending. The results for far field tension agreed well with published theoretical results over a good range of elastic gradients. The maximum shear stress calculated over the entire domain of the cracked plate using boundary collocation agrees very well with that obtained from finite element analysis. The efficacy of the extended stress field in capturing the effects of the elastic gradient on the stresses and fracture parameters is thus established in this study.     

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

A hybrid finite element analysis of interface cracks

A versatile hybrid finite element scheme consisting of special crack-tip elements and crack face contact elements is developed to analyse a partially closed interface crack between two dissimilar anisotropic elastic materials. The crack-tip element incorporates higher-order asymptotic solutions for an interfacial crack tip. These solutions are obtained from complex variable methods in Stroh formalism. For a closed interfacial crack tip, a generalized contact model in which the crack-tip oscillation is eliminated is adopted in the calculation. The hybrid finite element modelling allows the stress singularity at an open and closed crack tip to be accurately treated. The accuracy and convergence of the developed scheme are tested with respect to the known interface crack solutions. Utilizing this numerical scheme, the stress intensity factors and contact zone are calculated for a finite interface crack between a laminated composite material.     

On the evaluation of the J-integral by a plastic crack tip element

Two crack tip elements are formulated for a stationary, mode I plastic crack in planar structures using hybrid assumed stress approach, based on the secant modulus and the Newton-Raphson schemes, respectively. The stress distribution in the crack tip element is assumed to be the HRR field superimposed by the regular polynomial terms. The formulated (hybrid) crack tip elements are compatible with the isoparametric element so that they can be used conveniently along with the conventional displacement-based finite elements. The intensity of the HRR stress field, the J-integral, is determined directly from the finite element equations together with the nodal displacements. The dominance of the HRR stress field at the crack tip is pertinent to the present approach, which depends on geometry and loading conditions. Since the J-integral is globally path-independent for nonlinear elastic materials (deformation plasticity model), in order to assess the accuracy and efficiency of the methodology as compared to the contour integration approach, numerical studies of common plane-stress cracked configurations are performed for these materials. The results indicate that for a sufficiently small crack tip element size, J from the present approach correlates well, within 6 percent difference, with that from the contour integration for a wide range of material hardening coefficients if the HRR zone exists at the crack tip. These highly accurate results for J from the crack tip stresses could not be achieved without using (newly) modified variational principles and a refined numerical technique. It should be emphasized that the present methodology also can be applied to cracks in J ₂ flow materials under HRR dominance. In such case, the J integral may not be globally path independent, and hence it now must be determined from the stress and strain fields near the crack tip.