Boundary element analysis of three-dimensional crack problems in two joined transversely isotropic solids

The present paper presents a boundary element analysis of penny-shaped crack problems in two joined transversely isotropic solids. The boundary element analysis is carried out by incorporating the fundamental singular solution for a concentrated point load in a transversely isotropic bi-material solid of infinite space into the conventional displacement boundary integral equations. The conventional multi-region method is used to analyze the crack problems. The traction-singular elements are employed to capture the singularity around the crack front. The values of the stress intensity factors are obtained by using crack opening displacements. The numerical scheme results are verified with the closed-form solutions available in the literature for a penny-shaped crack parallel to the plane of the isotropy of a homogeneous and transversely isotropic solid of infinite extent. The new problem of a penny-shaped crack perpendicular to the interface of a transversely isotropic bi-material solid is then examined in detail. The crack surfaces are subject to the three normal tractions and the uniform shear traction. The associated stress intensity factor values are obtained and analyzed. The present results can be used for the prediction of the stability of composite structures and the hydraulic fracturing in deep rock strata and reservoir engineering.     

Stress intensity factors of square crack inclined to interface of transversely isotropic bi-material

This paper analyzes a square crack in a transversely isotropic bi-material solid by using dual boundary element method. The square crack is inclined to the interface of the bi-material. The fundamental solution for the bi-material solid occupying an infinite region is incorporated into the dual boundary integral equations. The square crack can have an arbitrary angle with respect to the plane of isotropy of the bi-material occupying either finite or infinite regions. The stress intensity factor (SIF) values of the modes I, II, and III associated with the square crack are calculated from the crack opening displacements. Numerical results show that the properties of the anisotropic bi-material have evident influences on the values of the three SIFs. The values of the three SIFs are further examined by taking into account the effect of the external boundary of the internally cracked bi-material.     

Application of finite-part integrals to three-dimensional fracture problems for piezoelectric media Part II: Numerical analysis

In this paper, elliptical cracks and rectangular cracks embedded in a three-dimensional infinite transversely isotropic piezoelectric solid are analyzed under combined mechanical tension and electric fields. The hypersingular integral equation method is used to solve the mentioned problems. The unknown function in the hypersingular integral equations is approximated with a product of the fundamental density function and polynomials. The hypersingular integrals can be numerically evaluated by using a method of Taylor series expansion. Therefore, the hypersingular integral equations for the crack problems can be solved immediately. Finally, numerical examples of the stress and electric displacement intensity factors as well as the energy release rates for these crack configurations are presented. The numerical results demonstrate the present approach to be very efficient.     

Boundary element analysis of three‐dimensional cracks in anisotropic solids

This paper presents a boundary element analysis of linear elastic fracture mechanics in three‐dimensional cracks of anisotropic solids. The method is a single‐domain based, thus it can model the solids with multiple interacting cracks or damage. In addition, the method can apply the fracture analysis in both bounded and unbounded anisotropic media and the stress intensity factors (SIFs) can be deduced directly from the boundary element solutions. The present boundary element formulation is based on a pair of boundary integral equations, namely, the displacement and traction boundary integral equations. While the former is collocated exclusively on the uncracked boundary, the latter is discretized only on one side of the crack surface. The displacement and/or traction are used as unknown variables on the uncracked boundary and the relative crack opening displacement (COD) (i.e. displacement discontinuity, or dislocation) is treated as a unknown quantity on the crack surface. This formulation possesses the advantages of both the traditional displacement boundary element method (BEM) and the displacement discontinuity (or dislocation) method, and thus eliminates the deficiency associated with the BEMs in modelling fracture behaviour of the solids. Special crack‐front elements are introduced to capture the crack‐tip behaviour. Numerical examples of stress intensity factors (SIFs) calculation are given for transversely isotropic orthotropic and anisotropic solids. For a penny‐shaped or a square‐shaped crack located in the plane of isotropy, the SIFs obtained with the present formulation are in very good agreement with existing closed‐form solutions and numerical results. For the crack not aligned with the plane of isotropy or in an anisotropic solid under remote pure tension, mixed mode fracture behavior occurs due to the material anisotropy and SIFs strongly depend on material anisotropy. Copyright © 2000 John Wiley & Sons, Ltd.     

A frequency-domain BEM for 3D non-synchronous crack interaction analysis in elastic solids

A frequency-domain boundary element method (BEM) is presented for non-synchronous crack interaction analysis in three-dimensional (3D), infinite, isotropic and linear elastic solids with multiple coplanar cracks. The cracks are subjected to non-synchronous time-harmonic crack-surface loading with contrast frequencies. Hypersingular frequency-domain traction boundary integral equations (BIEs) are applied to solve the boundary value problem. A collocation method is adopted for solving the BIEs numerically. The local square-root behavior of the crack-opening-displacements at the crack-front is taken into account in the present method. For two coplanar penny-shaped cracks of equal radius subjected to non-synchronous time-harmonic crack-surface loading, numerical results for the dynamic stress intensity factors are presented and discussed.     

Thermoelectroelastic analysis of cracks in piezoelectric half-plane by BEM

The Green's function and the boundary element method for analysing fracture behaviour of cracks in piezoelectric half-plane are presented in this paper. By combining Stroh formalism and the concept of perturbation, a general thermoelectroelastic solution for half-plane solid subjected to point heat source and/or temperature discontinuity has been derived. Using the proposed solution and the potential variational principle, a boundary element model (BEM) for 2-D half-plane solid with multiple cracks has been developed and used to calculate the stress intensity factors of the multiple crack problem. The method is available for multiple crack problems in both finite and infinite solids. Numerical results for a two-crack system are presented and compared with those from finite element method (FEM).     

BEM for crack-hole problems in thermopiezoelectric materials

The boundary element formulation for analysing interaction between a hole and multiple cracks in piezoelectric materials is presented. Using Green's function for hole problems and variational principle, a boundary element model (BEM) for a 2-D thermopiezoelectric solid with cracks and holes has been developed and used to calculate stress intensity factors of the crack-hole problem. In BEM, the boundary condition on the hole circumference is satisfied a priori by Green's function, and is not involved in the boundary element equations. The method is applicable to multiple-crack problems in both finite and infinite solids. Numerical results for stress and electric displacement intensity factors at a particular crack tip in a crack-hole system of piezoelectric materials are presented to illustrate the application of the proposed formulation.     

Computational modeling of strong and weak discontinuities using the extended finite element method

In this article, the extended finite element method is employed to solve problems, including weak and strong discontinuities. To this end, a level set framework is used to represent the discontinuities location, and the Heaviside and Branch function are included in the standard finite element method. The case of two arbitrary curved cracks is solved numerically and stress intensity factor (SIF) values at the crack tips are calculated based on the evaluation of the crack tip opening displacement. Afterwards, J-integral methodology is adopted to evaluate the SIFs for isotropic and anisotropic bi-material interface crack problems. Numerical results are verified with those presented in the literature.     

Analysis of cracked transversely isotropic and inhomogeneous solids by a special BIE formulation

In this paper, a special boundary integral equation (BIE) formation is proposed to analyze the fracture problem in transversely isotropic and inhomogeneous solids. In this formulation, the single-domain boundary element method (BEM) is utilized to discretize the cracked matrix and the displacement BEM to the surface of the embedded inhomogeneity. The two regions are then connected through the continuity conditions along their joint interface. The conventional and three special nine-node quadrilateral elements are utilized to discretize the inhomogeneity–matrix interface and the crack surface. From the crack-opening displacements on the crack surface, the mixed-mode stress intensity factors (SIFs) are calculated, using the well-known asymptotic expression in terms of the Barnett–Lothe tensor. In the numerical analysis, the distance between the inhomogeneity and the crack as well as the orientation of the isotropic plane of the transversely isotropic media is varied to show their influences on the mixed-mode SIFs along the crack fronts.     

Generalized Kelvin solution based boundary element method for crack problems in multilayered solids

This paper presents a new boundary element method (BEM) for linear elastic fracture mechanics in three-dimensional multilayered solids. The BEM is based on a generalized Kelvin solution. The generalized Kelvin solution is the fundamental singular solution for a multilayered elastic solid subject to point concentrated body-forces. For solving three-dimensional elastic crack problems in a finite region, a multi-region method is also employed in the present BEM. For crack problems in an infinite space, a large finite body is used to approximate the infinite body. In addition, eight-node traction-singular boundary elements are used in representing the displacements and tractions in the vicinity of a crack front. The incorporation of the generalized Kelvin solution into the boundary integral formulation has the advantages in elimination of the element discretization at the interfaces of different elastic layers. Three numerical examples are presented to illustrate the proposed method for the calculation of stress intensity factors for cracks in layered solids. The results obtained using the proposed method are well compared with the existing results available in the relevant literature.     

An Altgernating Method Based on the VNA Solution for Analysis of Damaged Solid Containing Arbitrarily Distributed Elliptical Microcracks

This paper presents the development of an alternating method for the interaction analysis of arbitrary distributed numerous elliptical microcracks. The complete analytical solutions (VNA solutions) for a single elliptical crack in an infinite solid, subject to arbitrary crack-face tractions, are implemented in the present alternating method, together with the coordinate transformations for stress tensors. First, the present method is verified by solving the problems of two interacting cracks for which accurate numerical solutions have been obtained previously. Next, the present method demonstrates obtaining efficient and accurate solutions for the problems of many interacting elliptical cracks, which cannot be solved in a practical sense by the ordinary numerical methods such as the finite element method. Furthermore, damaged solids containing periodically distributed elliptical microcracks are analyzed by the present alternating method. The effective elastic moduli are evaluated for varying microcrack density. Detailed structures of the interactions in the damaged solids are visualized and clarified. This revised version was published online in August 2006 with corrections to the Cover Date.     

An alternating method based on the VNA solution for analysis of damaged solid containing arbitrarily distributed elliptical microcracks

This paper presents the development of an alternating method for the interaction analysis of arbitrary distributed numerous elliptical microcracks. The complete analytical solutions (VNA solutions) for a single elliptical crack in an infinite solid, subject to arbitrary crack-face tractions, are implemented in the present alternating method, together with the coordinate transformations for stress tensors. First, the present method is verified by solving the problems of two interacting cracks for which accurate numerical solutions have been obtained previously. Next, the present method demonstrates obtaining efficient and accurate solutions for the problems of many interacting elliptical cracks, which cannot be solved in a practical sense by the ordinary numerical methods such as the finite element method. Furthermore, damaged solids containing periodically distributed elliptical microcracks are analyzed by the present alternating method. The effective elastic moduli are evaluated for varying microcrack density. Detailed structures of the interactions in the damaged solids are visualized and clarified. This revised version was published online in July 2006 with corrections to the Cover Date.     

Boundary integral equation method for conductive cracks in two and three-dimensional transversely isotropic piezoelectric media

Using the fundamental solutions and the Somigliana identity of piezoelectric medium, the boundary integral equations are obtained for a conductive planar crack of arbitrary shape in three-dimensional transversely isotropic piezoelectric medium. The singular behaviors near the crack edge are studied by boundary integral equation approach, and the intensity factors are derived in terms of the displacement discontinuity and the electric displacement boundary value sum near the crack edge on crack faces. The boundary integral equations for two dimensional crack problems are deduced as a special case of infinite strip planar crack. Based on the analogy of the obtained boundary integral equations and those for cracks in conventional isotropic elastic material and for contact problem of half-space under the action of a rigid punch, an analysis method is proposed. As an example, the solution to conductive Griffith crack is derived.     

Analytical solution for embedded elliptical cracks,and finite element alternating method for elliptical surface cracks,subjected to arbitrary loadings

The complete solution for an embedded elliptical crack in an infinite solid and subjected to arbitrary tractions on the crack surface is rederived from Vijayakumar and Atluri's general solution procedure. The general procedure for evaluating the necessary elliptic integrals in the generalized solution for elliptical crack is also derived in this paper. The generalized solution is employed in the Schwartz alternating technique in conjunction with the finite element method. This finite element-alternating method gives an inexpensive way to evaluate accurate stress intensity factors for embedded or elliptical cracks in engineering structural components.     

Dynamic crack problems in three-dimensional transversely isotropic solids

In this paper, a general boundary element approach for three-dimensional dynamic crack problems in transversely isotropic bodies is presented for the first time. Quarter-point and singular quarter-point elements are implemented in a quadratic isoparametric element context. The procedure is based on the subdomain technique, the displacement integral representation for elastodynamic problems and the expressions of the time-harmonic point load fundamental solution for transversely isotropic media. Numerical results corresponding to cracks under the effects of impinging waves are presented. The accuracy of the present approach for the analysis of dynamic fracture mechanics problems in transversely isotropic solids is shown by comparison of the obtained results with existing solutions.     

The gauss-hermite numerical integration method for the solution of the plane elastic problem of semi-infinite periodic cracks

The problem of parallel semi-infinite periodic cracks subjected to a transversely directed load in an infinite isotropic and elastic medium under conditions of plane stress or plane strain can be reduced to the solution of a Cauchy-type singular integral equation along one of the cracks. This equation can be transformed into a system of linear equations by means of an approximation of the integrals through the Gauss-Hermite procedure and application of the equation to distinct points along the faces of the crack. Stress intensity factors thus determined for the crack tips under constant load along the cracks are in satisfactory agreement with corresponding values derived previously.     

Boundary element analysis of thermal stress intensity factors for interface griffith and cusp cracks

The thermal stress intensity factors for interface cracks of Griffith and symmetric lip cusp types under vertical uniform heat flow in a finite body are calculated by the boundary element method. The boundary conditions on the crack surfaces are insulated or fixed to constant temperature. The relationship between the stress intensity factors and the displacements on the nodal point of a crack-tip element is derived. The numerical values of the thermal stress intensity factors for an interface Griffith crack in an infinite body are compared with the previous solutions. The thermal stress intensity factors for a symmetric lip cusp interface crack in a finite body are calculated with respect to various effective crack lengths, configuration parameters, material property ratios and the thermal boundary conditions on the crack surfaces. Under the same outer boundary conditions, there are no appreciable differences in the distribution of thermal stress intensity factors with respect to each material property. However, the effect of crack surface thermal boundary conditions on the thermal stress intensity factors is considerable.     

On the problem of two coplanar cracks inside an infinite isotropic elastic solid

A method is proposed for the approximate evaluation of normal displacements and normal stresses on the plane of two coplanar cracks located inside an infinite isotropic elastic solid and subjected to normal internal pressure. The formulation results in a single integral equation for the unknown normal stresses on the plane of the cracks. Numerical results are given for the stress intensity factor K_I of two coplanar circular cracks and two coplanar elliptical cracks opened up under a uniform internal pressure.     

The effect of material combinations and relative crack size to the stress intensity factors at the crack tip of a bi-material bonded strip

Several types of singular stress fields may appear at the corner where an interface between two bonded materials intersects a traction-free edge depending on the material combinations. Since the failure of the multi-layer systems often originates at the free-edge corner, the analysis of the edge interface crack is the most fundamental to simulate crack extension. In this study, the stress intensity factors for an edge interfacial crack in a bi-material bonded strip subjected to longitudinal tensile stress are evaluated for various combinations of materials using the finite element method. Then, the stress intensity factors are calculated systematically with varying the relative crack sizes from shallow to very deep cracks. Finally, the variations of stress intensity factors of a bi-material bonded strip are discussed with varying the relative crack size and material combinations. This investigation may contribute to a better understanding of the initiation and propagation of the interfacial cracks.     

Computation of shear mode weight functions for circular and elliptical cracks

Three-dimensional shear mode fundamental fields in finite bodies with mixed boundary conditions are analyzed by a special finite element method for circular and elliptical cracks. A procedure for determining the Fourier coefficients of the stress intensity factor for circular cracks is presented. A special series is proposed to represent the computed crack face weight functions for elliptical cracks.