Integral formulation for elastodynamic T-stresses

In this paper, a path independent integral formulation is presented for the computation of dynamic T-stresses in a two-dimensional body with a stationary crack. The mutual M-integral expressed through dynamic Ĵ-integrals provides sufficient information for determining T-stresses on the basis of the relationship found between the M-integral and T-stresses. The elastodynamic fields required for the evaluation of the Ĵ-and M-integrals are obtained by the boundary element method. The time-domain approach is used for the solution of the boundary value crack problem and numerical results for two crack problems are presented. In the first a rectangular plate with a central crack is considered and in the second two cracks at a hole in an infinite sheet. A comparison is made with the results obtained by the boundary layer and displacement field methods based on the asymptotic expansions of stresses and displacements at a crack tip vicinity. This revised version was published online in August 2006 with corrections to the Cover Date.     

J-Dominance in Mixed Mode Ductile Fracture Specimens

The asymptotic mixed mode crack tip fields in elastic-plastic solids are scaled by the J-integral and parameterized by a near-tip mixity parameter, M _{_p} . In this paper, the validity and range of dominance of these fields are investigated. To this end, small strain elastic-plastic finite element analyses of mixed mode fracture are first performed using a modified boundary layer formulation. Here, a two term expansion of the elastic crack tip field involving the stress intensity factor |K| the elastic mixity parameter M _{_e} as well as the T-stress is prescribed as remote boundary conditions. The analyses are conducted for different values of M _{_e} and the T-stress. Next, several commonly used mixed mode fracture specimens such as Compact Tension Shear (CTS), Four Point Bend (4PB), and modified Compact Tension specimen are considered. Here, the complete range of loading from contained yielding to large scale yielding is analyzed. Further, different crack to width ratios and strain hardening exponents are considered. The results obtained establish that the mixed mode asymptotic fields dominate over physically relevant length scales in the above geometries, except for predominantly mode I loading and under large scale yielding conditions. This revised version was published online in August 2006 with corrections to the Cover Date.     

Interfacial crack-tip constraints and <Emphasis Type="Italic">J</Emphasis>-integral for bi-materials with plastic hardening mismatch

This paper investigates interfacial crack tip stress fields and the J-integral for bi-materials with plastic hardening mismatch via detailed elastic-plastic finite element analyses. For small scale yielding, the modified boundary layer formulation with the elastic T-stress is employed. For fully plastic yielding, plane strain single-edge- cracked specimens under pure bending are considered. Interfacial crack tip stress fields are explained by modified Prandtl slip-line fields. It is found that, for bi-materials consisting of two elastic-plastic materials, increasing plastic hardening mismatch increases both crack-tip stress constraint in the lower hardening material and the J-contribution there. The implication of asymmetric J-integral in bi-materials is also discussed.     

The dual boundary element method: Effective implementation for crack problems

The present paper is concerned with the effective numerical implementation of the two-dimensional dual boundary element method, for linear elastic crack problems. The dual equations of the method are the displacement and the traction boundary integral equations. When the displacement equation is applied on one of the crack surfaces and the traction equation on the other, general mixed-mode crack problems can be solved with a single-region formulation. Both crack surfaces are discretized with discontinuous quadratic boundary elements; this strategy not only automatically satisfies the necessary conditions for the existence of the finite-part integrals, which occur naturally, but also circumvents the problem of collocation at crack tips, crack kinks and crack-edge corners. Examples of geometries with edge, and embedded crack are analysed with the present method. Highly accurate results are obtained, when the stress intensity factor is evaluated with the J-integral technique. The accuracy and efficiency of the implementation described herein make this formulation ideal for the study of crack growth problems under mixed-mode conditions.     

Stress intensity factor analyses of three-dimensional interface cracks using tetrahedral finite elements

A numerical method is proposed for evaluating the stress intensity factors of a three-dimensional bimaterial interfacial crack using tetrahedral finite elements. This technique is based on the M₁-integral method and employs the moving least-squares approximation. Stress or strain in the M₁-integral equation is automatically approximated from the nodal displacements obtained by the finite element analysis using the moving least-squares method. Therefore, the presented method needs no elemental information from the finite element analysis. In this study, stress intensity factor analyses of some typical three-dimensional interface crack problems using the tetrahedral finite elements are demonstrated.     

Fracture analysis of magnetoelectroelastic solids by using path independent integrals

A solution scheme based on the fundamental solution for a generalized edge dislocation in an infinite magnetoelectroelastic solid is presented to analyze problems involving single, multiple and slowly growing impermeable cracks. The fundamental solution for a generalized dislocation is obtained by extending the complex potential function formulation used for anisotropic elasticity. The solution for a continuously distributed dislocation is derived by integrating the solution for an edge dislocation. The problem of a system of cracks subjected to remote mechanical, electric and magnetic loading is formulated in terms of set of singular integral equations by applying the principle of superposition and the solution for a continuously distributed dislocation. The singular integral equation system is solved by using a numerical integration technique based on Chebyshev polynomials. The J_i and M-integrals for single crack and multi-cracks problems are derived and their dependence on the coordinate system is investigated. Selected numerical results for the M-integral, total energy release rate and mechanical energy release rate are presented for single, double and multiple crack problems. The case of a slowly growing crack interacting with a stationary crack is also considered. It is found that M-integral presents a reliable and physically acceptable measure for assessment of fracture behaviour and damage of magnetoelectroelastic materials.     

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.     

A non-singular boundary integral formula for determining the T-stress for cracks of arbitrary geometry

A simple, yet accurate 2-D boundary integral equation (BIE) for determining the T-stress for cracks of arbitrarily geometry is introduced in this paper. The formulation is based upon the asymptotic expansion for the stress field in the vicinity of a crack tip. It can be conveniently implemented in the post-processing stage of a boundary element fracture analysis. As demonstrated in this work, the proposed BIE is non-singular, and thus it can directly be collocated at the crack tip under consideration. The technique requires a similar computational effort as that used in calculating the stress components at an interior point of a domain. Consequently, this new approach is very computationally effective and accurate for evaluating the elastic T-stress. Five test examples, involving straight, kink and curved cracks, are studied to validate the proposed technique and to assess its accuracy.     

Boundary element analysis for an interface crack between dissimilar elastoplastic materials

The boundary element method (BEM) is presented for elastoplastic analysis of cracks between two dissimilar materials. The boundary integral equations and integral representation of stress rates are written in such a form that all integrals can be evaluated by the regular Gaussian quadrature rule. An advanced multidomain BEM formulation is suggested for the solution of analysed problems where the substantial reduction of stiffness matrix is observed. The elastoplastic behaviour is modelled through the use of an approximation for the plastic component of the stresses. The boundary and the yielding zone are discretized by elements with quadratic approximations. In numerical examples the path independence of the J- and L-integrals for a straight interface crack and a circular arc-shaped interface crack are investigated, respectively. The influence of the different values of Young's modulus on the J-integral, shape and size of plastic zones is treated too.     

A new boundary element method for mixed boundary value problems involving cracks and holes: Interactions between rigid inclusions and cracks

This paper derives a new boundary integral equation (BIE) formulation for plane elastic bodies containing cracks and holes and subjected to mixed displacement/ traction boundary conditions, and proposes a new boundary element method (BEM) based upon this formulation. The basic unknown in the formulation is a complex boundary function H(t), which is a linear combination of the boundary traction and boundary displacement density. The present BIE formulation can be related directly to Muskhelishvili's formalism. Singular interpolation functions of order r ^–1/2 (where r is the distance measured from the crack tip) are introduced such that singular integrand involved at the element level can be integrated analytically. By applying the BEM, the interaction between a rigid circular inclusion and a crack is investigated in details. Our results for the stress intensity factor are comparable with those given by Erdogan and Gupta (1975) and Gharpuray et al. (1990) for a crack emanating from a stiff inclusion, and with those by Erdogan et al. (1974) for a crack in the neighborhood of a stiff inclusion.     

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.     

A new M-integral parameter for mixed-mode crack growth in orthotropic viscoelastic material

This paper deals with a new independent path integral which provides the mixed-mode during a creep crack growth process in viscoelastic orthotropic media. The developments are based on an energetic approach using conservative laws. The mixed-mode fracture separation is introduced according to the generalization of the virtual work principle. The fracture algorithm is implemented in a finite element software and coupled with an incremental viscoelastic formulation and an automatic crack growth simulation. This M-integral provides the computation of stress intensity factors and energy release rate for each fracture mode. A numerical validation, in terms of energy release rate and stress intensity factors, is carried out on a CTS specimen under mixed-mode loading for different crack growth speeds.     

Continuum Shape Sensitivity analysis of a mode-I fracture in functionally graded materials

This paper presents a new method for conducting a continuum shape sensitivity analysis of a crack in an isotropic, linear-elastic, functionally graded material. This method involves the material derivative concept from continuum mechanics, domain integral representation of the J-integral and direct differentiation. Unlike virtual crack extension techniques, no mesh perturbation is needed to calculate the sensitivity of stress-intensity factors. Since the governing variational equation is differentiated prior to the process of discretization, the resulting sensitivity equations are independent of approximate numerical techniques, such as the meshless method, finite element method, boundary element method, or others. In addition, since the J-integral is represented by domain integration, only the first-order sensitivity of the displacement field is needed. Several numerical examples are presented to calculate the first-order derivative of the J-integral, using the proposed method. Numerical results obtained using the proposed method are compared with the reference solutions obtained from finite-difference methods for the structural and crack geometries considered in this study.     

Formulation and implementation of a high-order 3-D domain integral method for the extraction of energy release rates

This article presents a three dimensional (3-D) formulation and implementation of a high-order domain integral method for the computation of energy release rate. The method is derived using surface and domain formulations of the J-integral and the weighted residual method. The J-integral along 3-D crack fronts is approximated by high-order Legendre polynomials. The proposed implementation is tailored for the Generalized/eXtended Finite Element Method and can handle discontinuities arbitrarily located within a finite element mesh. The domain integral calculations are based on the same integration elements used for the computation of the stiffness matrix. Discontinuities of the integrands across crack surfaces and across computational element boundaries are fully accounted for. The proposed method is able to deliver smooth approximations and to capture the boundary layer behavior of the J-integral using tetrahedral meshes. Numerical simulations of mode-I and mixed mode benchmark fracture mechanics examples verify expected convergence rates for the computed energy release rates. The results are also in good agreement with other numerical solutions available in the literature.     

J resistance behavior in functionally graded materials using cohesive zone and modified boundary layer models

This paper describes elastic–plastic crack growth resistance simulation in a ceramic/metal functionally graded material (FGM) under mode I loading conditions using cohesive zone and modified boundary layer (MBL) models. For this purpose, we first explore the applicability of two existing, phenomenological cohesive zone models for FGMs. Based on these investigations, we propose a new cohesive zone model. Then, we perform crack growth simulations for TiB/Ti FGM SE(B) and SE(T) specimens using the three cohesive zone models mentioned above. The crack growth resistance of the FGM is characterized by the J-integral. These results show that the two existing cohesive zone models overestimate the actual J value, whereas the model proposed in the present study closely captures the actual fracture and crack growth behaviors of the FGM. Finally, the cohesive zone models are employed in conjunction with the MBL model. The two existing cohesive zone models fail to produce the desired K–T stress field for the MBL model. On the other hand, the proposed cohesive zone model yields the desired K–T stress field for the MBL model, and thus yields J _R curves that match the ones obtained from the SE(B) and SE(T) specimens. These results verify the application of the MBL model to simulate crack growth resistance in FGMs.     

Generalization of <Emphasis Type="Italic">T</Emphasis> and <Emphasis Type="Italic">A</Emphasis> integrals to time-dependent materials: analytical formulations

This paper deals with the generalization of T-integral to crack growth process in viscoelastic materials. In order to implement this expression in a finite element software, a modelling form of this integral, called Aθ, is developed. The analytical formulation is based on conservative law, independent path integral, and a combination of real, virtual displacement fields, and real, virtual thermal fields introducing, in the same time, a bilinear form of free energy density F. According to the generalization of Noether’s method, the application of Gauss Ostrogradski’s theorem combined with curvilinear cracked contour, T _v is obtained. By introducing a volume domain around crack tip, the modelling expression Aθ is also defined.. Finally, the viscoelastic generalization through a thermodynamic approach, called A _v , is introduced by using a discretisation of the creep tensor according to a generalized Kelvin Voigt representation.     

Application of a conservation integral to an interface crack interacting with singularities

A mutual integral, which has the conservation property is applied to the problem of an interfacial crack. The stress intensity factors K ₁, K ₂, K ₃ and T-stress for the problem in an infinite medium are easily obtained by using the mutual integral without solving the boundary value problem. In doing so the auxiliary solutions are required and they are taken from the known asymptotic solutions. This method is amenable to numerical evaluation of the stress intensity factors and T-stress if the interfacial crack in a finite medium is considered.     

Stress intensity factor analysis of an interface crack between dissimilar anisotropic materials under thermal stress using the finite element analysis

New numerical methods were presented for stress intensity factor analyses of two-dimensional interfacial crack between dissimilar anisotropic materials subjected to thermal stress. The virtual crack extension method and the thermal M-integral method for a crack along the interface between two different materials were applied to the thermoelastic interfacial crack in anisotropic bimaterials. The moving least-squares approximation was used to calculate the value of the thermal M-integral. The thermal M-integral in conjunction with the moving least-squares approximation can calculate the stress intensity factors from only nodal displacements obtained by the finite element analysis. The stress intensity factors analyses of double edge cracks in jointed dissimilar isotropic semi-infinite plates subjected to thermal load were demonstrated. Excellent agreement was achieved between the numerical results obtained by the present methods and the exact solution. In addition, the stress intensity factors of double edge cracks in jointed dissimilar anisotropic semi-infinite plates subjected to thermal loads were analyzed. Their results appear reasonable.     

Applications of the element-free Galerkin method for singular stress analysis under strain gradient plasticity theories

It is known that the plasticity models affect characterization of the crack tip fields. To predict failure one has to understand the crack tip stress field and control the crack. In the present work the element-free Galerkin methods for gradient plasticity theories have been developed and implemented into the commercial finite element code ABAQUS and used to analyze crack tip fields. Based on the modified boundary layer formulation it is confirmed that the stress singularity in the gradient plasticity theories is significantly higher than the known HRR solution and seems numerically to equal to 0.78, independently of the strain-hardening exponent. The strain singularity is much lower than the known HRR one. The crack field in gradient plasticity under small-scale yielding condition consists of three zones: The elastic K-field, the plastic HRR-field dominated by the J-integral and the hyper-singular stress field. Even under gradient plasticity there exists an HRR-zone described by the known J-integral, whereas the hyper-singular zone cannot be characterized by J. The hyper-singular zone is very small (r ? J/σ₀) and contained by the HRR zone in the infinitesimal deformation framework. The finite strains under the gradient plasticity will not eliminate the stress singularity as r → 0, in contrast to the known finite strain results under the Mises plasticity. Numerically no significant changes in characterization of the stress field were found in comparison with the infinitesimal deformation theory. Since the hyper-singular stress field is much smaller than the HRR zone and in the same size as the fracture process zone, one may still use the known J concept to control the crack in the gradient plasticities. In this sense the gradient plasticity will not change characterization of the crack.     

Fracture mechanics analysis of anisotropic plates by the boundary element method

This paper presents an analysis of mixed mode fracture mechanics problems arising in anisotropic composite laminates. The boundary element method (BEM) and the J _k integral are presented as accurate techniques to compute the stress intensity factors K _I and K _II of two dimensional anisotropic bodies. Using function of a complex variable a decoupling procedure is derived to obtain the stress intensity factors. The procedure is based on the computation of the J ₁-integral and of the ratio of relative displacements at the crack faces, near the crack tip. Applications are presented for unidirectional and symmetric laminates of glass, boron and graphite-epoxy materials. Numerical examples of problems of pure mode I and mixed mode deformations are given, in order to demonstrate the accuracy of the method.