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.     

A super‐element for crack analysis in the time domain

A super‐element for the dynamic analysis of two‐dimensional crack problems is developed based on the scaled boundary finite‐element method. The boundary of the super‐element containing a crack tip is discretized with line elements. The governing partial differential equations formulated in the scaled boundary co‐ordinates are transformed to ordinary differential equations in the frequency domain by applying the Galerkin's weighted residual technique. The displacements in the radial direction from the crack tip to a point on the boundary are solved analytically without any a priori assumption. The scaled boundary finite‐element formulation leads to symmetric static stiffness and mass matrices. The super‐element can be coupled seamlessly with standard finite elements. The transient response is evaluated directly in the time domain using a standard time‐integration scheme. The stress field, including the singularity around the crack tip, is expressed semi‐analytically. The stress intensity factors are evaluated without directly addressing singular functions, as the limit in their definitions is performed analytically. Copyright © 2004 John Wiley & Sons, Ltd.     

Time domain boundary element method for dynamic stress intensity factor computations

This paper presents a procedure for transient dynamic stress intensity factor computations using traction singular quarter-point boundary elements in combination with the direct time domain formulation of the Boundary Element Method. The stress intensity factors are computed directly from the traction nodal values at the crack tip. Several examples of finite cracks in finite domains under mode-I and mixed mode dynamic loading conditions are presented. The computed stress intensity factors are represented versus time and compared with those obtained by other authors using different methods. The agreement is very good. The results are reliable and little mesh dependent. These facts allow for the analysis of dynamic crack problems with simple boundary discretizations. The versatile procedure presented can be easily applied to problems with complex geometry which include one or several cracks.     

Efficient solution of multiple cracks in great number using eigen COD boundary integral equations with iteration procedure

A newly developed computational approach is proposed in the paper for the analysis of multiple crack problems based on the eigen crack opening displacement (COD) boundary integral equations. The eigen COD particularly refers to a crack in an infinite domain under fictitious traction acting on the crack surface. With the concept of eigen COD, the multiple cracks in great number can be solved by using the conventional displacement discontinuity boundary integral equations in an iterative fashion with a small size of system matrix to determine all the unknown CODs step by step. To deal with the interactions among cracks for multiple crack problems, all cracks in the problem are divided into two groups, namely the adjacent group and the far-field group, according to the distance to the current crack in consideration. The adjacent group contains cracks with relatively small distances but strong effects to the current crack, while the others, the cracks of far-field group are composed of those with relatively large distances. Correspondingly, the eigen COD of the current crack is computed in two parts. The first part is computed by using the fictitious tractions of adjacent cracks via the local Eshelby matrix derived from the traction boundary integral equations in discretized form, while the second part is computed by using those of far-field cracks so that the high computational efficiency can be achieved in the proposed approach. The numerical results of the proposed approach are compared not only with those using the dual boundary integral equations (D-BIE) and the BIE with numerical Green's functions (NGF) but also with those of the analytical solutions in literature. The effectiveness and the efficiency of the proposed approach is verified. Numerical examples are provided for the stress intensity factors of cracks, up to several thousands in number, in both the finite and infinite plates.     

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.     

Automatic crack tip detection and stress intensity factors estimation of curved cracks from digital images

With the development of full‐field measurement techniques, it has been possible to analyze crack propagation experimentally with an increasing level of robustness. However, the analysis of curved cracks is made difficult and almost unexplored because the possible analysis domain size decreases with crack curvature, leading to an increasing uncertainty level. This paper proposes a digital image correlation technique, augmented by an elastic regularization, combining finite element kinematics on an adapted mesh and a truncated Williams' expansion. Through the analysis of two examples, the proposed technique is shown to be able to address the experimental problems of crack tip detection and stress intensity factors estimation along a curved crack path. Copyright © 2015 John Wiley & Sons, Ltd.     

Experimentelle Charakterisierung und zweidimensionale Simulation der mikrostrukturbestimmten Ermüdungsrissausbreitung in einer β‐Titanlegierung

In this project the initiation and propagation of short fatigue cracks in the metastable β‐titanium alloy TIMETAL®LCB is investigated. By means of an interferometric strain/displacement gauge system (ISDG) to measure the crack opening displacement (COD) and the electron back scattered diffraction technique (EBSD) to determine the orientation of individual grains the microstructural influence on short crack initiation and growth can be characterized. Finite element calculations show a high influence of the elastic anisotropy on the initiation sites of cracks. Crack propagation takes place transgranulary along slip planes as well as intergranulary along grain boundaries. The crack growth rate depends strongly on the active mechanism at the crack tip which in turn is influenced by crack length, the applied stress and the orientation of the grains involved. The value of the steady state crack closure stress changes from a positive value at low applied stresses (roughness induced) to a negative one at higher applied stresses (due to plastic deformations at the crack tip). The crack growth simulation is realised by a two‐dimensional boundary element technique, which contains the ideas of Navarro und de los Rios. The model includes the sequence of the applied stress amplitude as well as the experimental measured roughness induced crack closure.     

Analyzing interaction between coplanar square cracks using an efficient boundary element‐free method

This paper examines the interaction between coplanar square cracks by combining the moving least‐squares (MLS) approximation and the derived boundary integral equation (BIE). A new traction BIE involving only the Cauchy singular kernels is derived by applying integration by parts to the traditional boundary integral formulation. The new traction BIE can be directly applied to a crack surface and no displacement BIE is necessary because all crack boundary conditions (both upper and lower ones) are incorporated. A boundary element‐free method is then developed by combining the derived BIE and MLS approximation, in which the crack opening displacement is first expressed as the product of weight functions and the characteristic terms, and the unknown weight is approximated with the MLS approximation. The efficiency of the developed method is tested for isotropic and transversely isotropic media. The interaction between two and three coplanar square cracks in isotropic elastic body is numerically studied and the case of any number of coplanar square cracks is deduced and discussed. Copyright © 2012 John Wiley & Sons, Ltd.     

Transient analysis of the DSIFs and dynamic T-stress for particulate composite materials—Numerical vs. experimental results

The fracture behavior of particulate composite materials when subjected to dynamic loading has been a great concern for many industrial applications as these materials are particularly susceptible to impact loading conditions. As a result, many numerical and experimental techniques have been developed to deal with this class of problems. In this work, the fracture behavior of particulate composites under impact loading conditions is numerically studied via the two most important fracture parameters: dynamic stress intensity factors (DSIFs) and dynamic T-stress (DTS), and the results are compared with the experimental data obtained in Refs. [1,2]. Here, micromechanics models (self-consistent, Mori–Tanaka, …) or experimental techniques need to be employed first to determine the effective material properties of particulate composites. Then, the symmetric-Galerkin boundary element method for elastodynamics in the Fourier-space frequency domain is used in conjunction with displacement correlation technique to evaluate the DSIFs and stress correlation technique to determine the DTS. To obtain transient responses of the fracture parameters, fast Fourier transform (FFT) and inverse FFT are subsequently used to convert the DSIFs and DTS from the frequency domain to the time domain. Test examples involving free–free beams made of particulate composites are considered in this study. The numerical results are found to agree very well with the experimental ones.     

Near‐crack contour behaviour and extraction of log‐singular stress terms of the self‐regular traction boundary integral equation

This work contains an analytical study of the asymptotic near‐crack contour behaviour of stresses obtained from the self‐regular traction‐boundary integral equation (BIE), both in two and in three dimensions, and for various crack displacement modes. The flat crack case is chosen for detailed analysis of the singular stress for points approaching the crack contour. By imposing a condition of bounded stresses on the crack surface, the work shows that the boundary stresses on the crack are in fact zero for an unloaded crack, and the interior stresses reproduce the known inverse square root behaviour when the distance from the interior point to the crack contour approaches zero. The correct order of the stress singularity is obtained after the integrals for the self‐regular traction‐BIE formulation are evaluated analytically for the assumed displacement discontinuity model. Based on the analytic results, a new near‐crack contour self‐regular traction‐BIE is proposed for collocation points near the crack contour. In this new formulation, the asymptotic log‐singular stresses are identified and extracted from the BIE. Log‐singular stress terms are revealed for the free integrals written as contour integrals and for the self‐regularized integral with the integration region divided into sub‐regions. These terms are shown to cancel each other exactly when combined and can therefore be eliminated from the final BIE formulation. This work separates mathematical and physical singularities in a unique manner. Mathematical singularities are identified, and the singular information is all contained in the region near the crack contour. Copyright © 2003 John Wiley & Sons, Ltd.     

Boundary element method for <Emphasis Type="Italic">J</Emphasis>-integral and stress intensity factor computations in three-dimensional interface cracks

A general numerical tool for the analysis of three–dimensional bimaterial interface cracks is presented in this paper. The proposed tool is based on a multidomain formulation of the Boundary Element Method (BEM), with the crack located at the interface of the domain. Mixed mode stress intensity factors are computed along the three-dimensional crack fronts using the Energy Domain Integral (EDI) methodology and decoupled via the Interaction Integral. The capability of the procedure is demonstrated by solving a number of examples. The last of these examples consists in a thick centre cracked panel for which the behaviour of the J-integral and the mixed-mode stress intensity factors along the crack front is studied as a function of the material mismatch.     

A dual BIE approach for large‐scale modelling of 3‐D electrostatic problems with the fast multipole boundary element method

A dual boundary integral equation (BIE) formulation is presented for the analysis of general 3‐D electrostatic problems, especially those involving thin structures. This dual BIE formulation uses a linear combination of the conventional BIE and hypersingular BIE on the entire boundary of a problem domain. Similar to crack problems in elasticity, the conventional BIE degenerates when the field outside a thin body is investigated, such as the electrostatic field around a thin conducting plate. The dual BIE formulation, however, does not degenerate in such cases. Most importantly, the dual BIE is found to have better conditioning for the equations using the boundary element method (BEM) compared with the conventional BIE, even for domains with regular shapes. Thus the dual BIE is well suited for implementation with the fast multipole BEM. The fast multipole BEM for the dual BIE formulation is developed based on an adaptive fast multiple approach for the conventional BIE. Several examples are studied with the fast multipole BEM code, including finite and infinite domain problems, bulky and thin plate structures, and simplified comb‐drive models having more than 440 thin beams with the total number of equations above 1.45 million and solved on a PC. The numerical results clearly demonstrate that the dual BIE is very effective in solving general 3‐D electrostatic problems, as well as special cases involving thin perfect conducting structures, and that the adaptive fast multipole BEM with the dual BIE formulation is very efficient and promising in solving large‐scale electrostatic problems. Copyright © 2007 John Wiley & Sons, Ltd.     

Stress intensity factors for cracks emanating from a triangular or square hole in an infinite plate by boundary elements

This paper concerns stress intensity factors of cracks emanating from a triangular or square hole in an infinite plate subjected to internal pressure calculated by means of a boundary element method, which consists of constant displacement discontinuity element presented by Crouch and Starfied and crack tip displacement discontinuity elements proposed by the author. In the boundary element implementation the left or the right crack tip displacement discontinuity element is placed locally at corresponding left or right crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. Numerical examples are included to show that the method is very efficient and accurate for calculating stress intensity factors of plane elasticity crack problems. Specifically, the numerical results of stress intensity factors of cracks emanating from a triangular or square hole in an infinite plate subjected to internal pressure are given.     

Comparative study of time and frequency domain BEM approaches in frictional contact problem for antiplane crack under harmonic loading

Two different boundary element methods (BEM) for crack analysis in two dimensional (2-D) antiplane, homogeneous, isotropic and linear elastic solids by considering frictional contact of the crack edges are presented. Hypersingular boundary integral equations (BIE) in time-domain (TD) and frequency domain (FD), with corresponding elastodynamic fundamental solutions are applied for this purpose. For evaluation of the hypersingular integrals involved in BIEs a special regularization process that converts the hypersingular integrals to regular integrals is applied. Simple regular formulas for their calculation are presented. For the problems solution while considering frictional contact of the crack edges a special iterative algorithm of Udzava's type is elaborated and used. Numerical results for crack opening, frictional contact forces and dynamic stress intensity factors (SIFs) are presented and discussed for a finite III-mode crack in an infinite domain subjected to a harmonic crack-face loading and considering crack edges frictional contact interaction using the TD and FD approaches.     

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.     

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.     

2D SH modelling of ultrasonic testing for cracks near a non-planar surface

A model of 2D SH ultrasonic nondestructive testing for interior strip-like cracks near a non-planar back surface in a thick-walled elastic solid is presented. The model employs a Green's function to reformulate the 2D antiplane wave scattering problem as two coupled boundary integral equations (BIE): a displacement BIE for the back surface displacement and a hypersingular traction BIE for the crack opening displacement (COD). The integral equations are solved by performing a boundary element discretization of the back surface and expanding the COD in a series of Chebyshev functions which incorporate the correct behaviour at the crack edges. The transmitting ultrasonic probe is modelled by prescribing the traction underneath it, enabling the consequent calculation of the incident field. An electromechanical reciprocity relation is used to model the action of the receiving probe. A few numerical examples which illustrate the influence of the non-planar back surface are given.     

A boundary integral equation method for an inverse problem related to crack detection

This paper discusses an application of a boundary integral equation method (BIEM) to an inverse problem of determining the shape and the location of cracks by boundary measurements. Suppose that a given body contains an interior crack, the shape and the location of which are unknown. On the exterior boundary of this body one carries out measurements which are interpreted mathematically as prescribing Dirichlet data and measuring the corresponding Neumann data, or vice versa, for a field governed by Laplace's equation. The inverse problem considered here attempts to determine the geometry of the crack from these experimental data. We propose to solve this problem by minimizing the error of a certain boundary integral equation (BIE). The process of this minimization, however, is shown to require solutions of certain are proposed. Several 2D and 3D numerical examples are given in order to test the performance of the present method.     

A time‐dependent formulation of dual boundary element method for 3D dynamic crack problems

In this paper, the dual boundary element method in time domain is developed for three‐dimensional dynamic crack problems. The boundary integral equations for displacement and traction in time domain are presented. By using the displacement equation and traction equation on crack surfaces, the discontinuity displacement on the crack can be determined. The integral equations are solved numerically by a time‐stepping technique with quadratic boundary elements. The dynamic stress intensity factors are calculated from the crack opening displacement. Several examples are presented to demonstrate the accuracy of this method. Copyright © 1999 John Wiley & Sons, Ltd     

An effective method of stress intensity factor calculation for cracks emanating from a triangular or square hole under biaxial loads

This paper is concerned with stress intensity factors for cracks emanating from a triangular or square hole under biaxial loads by means of a new boundary element method. The boundary element method consists of the constant displacement discontinuity element presented by Crouch and Starfied and the crack‐tip displacement discontinuity elements proposed by the author. In the boundary element implementation, the left or the right crack‐tip displacement discontinuity element is placed locally at the corresponding left or right crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. The method is called a Hybrid Displacement Discontinuity Method (HDDM). Numerical examples are included to show that the method is very efficient and accurate for calculating stress intensity factors for plane elastic crack problems. In addition, the present numerical results can reveal the effect of the biaxial loads on stress intensity factors.