Interaction between cracks and rigid-line inclusions by an integral equation approach

 An integral equation approach is presented to investigate the interaction between cracks and rigid-line inclusions embedded in an infinite isotropic elastic matrix subject to remote loading. The relevant fundamental solutions in the integral formulation are presented. Special tip elements are used to simulate the variation of the discontinuous displacements over the crack surfaces, and the axial and shear forces along the rigid-line inclusions. The stress intensity factors at the crack tips and at the ends of the rigid-line inclusions are computed and compared with available solutions. Received: 6 August 2002 / Accepted: 3 February 2003 The work described in this paper was partially supported by a grant from the Research Grant Council of the Hong Kong Special Administration Region, China (Project No.: HKU 7011/01E). The authors would like to thank two reviewers for their constructive comments and suggestions to the paper. The comments of Professor H.P. Hong at the Department of Civil and Environmental Engineering of the University of Western Ontario of Canada are also appreciated.     

Thermal analysis of carbon-nanotube composites using a rigid-line inclusion model by the boundary integral equation method

The boundary integral equation (BIE) method is applied for the thermal analysis of fiber-reinforced composites, particularly the carbon-nanotube (CNT) composites, based on a rigid-line inclusion model. The steady state heat conduction equation is solved using the BIE in a two-dimensional infinite domain containing line inclusions which are assumed to have a much higher thermal conductivity (like CNTs) than that of the host medium. Thus the temperature along the length of a line inclusion can be assumed constant. In this way, each inclusion can be regarded as a rigid line (the opposite of a crack) in the medium. It is shown that, like the crack case, the hypersingular (derivative) BIE can be applied to model these rigid lines. The boundary element method (BEM), accelerated with the fast multipole method, is used to solve the established hypersingular BIE. Numerical examples with up to 10,000 rigid lines (with 1,000,000 equations), are successfully solved by the BEM code on a laptop computer. Effective thermal conductivity of fiber-reinforced composites are evaluated using the computed temperature and heat flux fields. These numerical results are compared with the analytical solution for a single inclusion case and with the experimental one reported in the literature for carbon-nanotube composites for multiple inclusion cases. Good agreements are observed in both situations, which clearly demonstrates the potential of the developed approach in large-scale modeling of fiber-reinforced composites, particularly that of the emerging carbon-nanotube composites.     

Degenerate scale for the analysis of circular thin plate using the boundary integral equation method and boundary element methods

In this paper, the degenerate scale for plate problem is studied. For the continuous model, we use the null-field integral equation, Fourier series and the series expansion in terms of degenerate kernel for fundamental solutions to examine the solvability of BIEM for circular thin plates. Any two of the four boundary integral equations in the plate formulation may be chosen. For the discrete model, the circulant is employed to determine the rank deficiency of the influence matrix. Both approaches, continuous and discrete models, lead to the same result of degenerate scale. We study the nonunique solution analytically for the circular plate and find degenerate scales. The similar properties of solvability condition between the membrane (Laplace) and plate (biharmonic) problems are also examined. The number of degenerate scales for the six boundary integral formulations is also determined. Tel.: 886-2-2462-2192-ext. 6140 or 6177     

Boundary integral equation method for two dissimilar elastic inclusions in an infinite plate

This paper provides a numerical solution for an infinite plate containing two dissimilar elastic inclusions, which is based on complex variable boundary integral equation (CVBIE). The original problem is decomposed into two problems. One is an interior boundary value problem (BVP) for two elastic inclusions, while other is an exterior BVP for the matrix with notches. After performing discretization for the coupled boundary integral equations (BIEs), a system of algebraic equations is formulated. The inverse matrix technique is suggested to solve the relevant algebraic equations, which can avoid using the assembling of some matrices. Several numerical examples are carried out to prove the efficiency of suggested method and the hoop stress along the interface boundary is evaluated.     

The integral equation formulations of an infinite elastic medium containing inclusions, cracks and rigid lines

The integral equation formulations of an infinite homogeneous isotropic medium containing various inclusions, cracks and rigid lines are presented. The present integral equation formulations contain the displacements (no tractions) over the inclusion-matrix interfaces, the discontinuous displacements over crack surfaces and the axial and the shear forces along rigid-line inclusions. Besides, the sub-domain boundary element method is also used in the present research. Numerical results from the present method and the sub-domain boundary element method are compared and discussed.     

Identification of cracks and cavities using the topological sensitivity boundary integral equation

The purpose of this communication is to present a novel approach to compute the so called Topological Sensitivity (TS) of any variable or functional in elasticity using Boundary Integral Equations (BIEs), and its use as a tool for identification of defects, by itself or in conjunction with zero-order methods, like Genetic Algorithms. The TS of a cost functional provides a measure of the susceptibility of a defect being at a given location. The main contributions are summarized in the following points:     

Boundary value problems of elastodynamics under stationary moving forces using boundary integral equation method

The boundary value problem of elastodynamics is considered in cylindrical domains when acting boundary forces are moving with constant speed parallel to cylinder's axis. The fundamental solutions, boundary integral equations and integral representation of the solutions are constructed using distribution theory for three kinds of speed: subsonic, transonic and supersonic.     

Solution for Eshelby's elastic inclusions in a finite plate using boundary integral equation method

This paper provides a solution for Eshelby's elastic inclusions in a finite plate based on the complex variable boundary integral equation (CVBIE) method. In the formulation, an inclusion with Eshelby's eigenstrains is embedded in an elliptic plate, and the exterior boundary is applied by some static loading. Two BIEs are suggested in the present study. One of BIEs is used for the finite inclusion region, and the other is used for region bounded by interface and the exterior boundary. After the discretization of BIEs, a numerical solution is suggested. In the solution, an inverse matrix technique is suggested which can eliminate one unknown vector in advance. Three numerical examples under different generalized loadings are provided. Interaction between the prescribed eigenstrains and the static loading along the exterior boundary is studied in detail.     

Solution of the second-order one-dimensional hyperbolic telegraph equation by using the dual reciprocity boundary integral equation (DRBIE) method

In this paper, we use a numerical method based on the boundary integral equation (BIE) and an application of the dual reciprocity method (DRM) to solve the second-order one space-dimensional hyperbolic telegraph equation. Also the time stepping scheme is employed to deal with the time derivative. In this study, we have used three different types of radial basis functions (cubic, thin plate spline and linear RBFs), to approximate functions in the dual reciprocity method (DRM). To confirm the accuracy of the new approach and to show the performance of each of the RBFs, several examples are presented. The convergence of the DRBIE method is studied numerically by comparison with the exact solutions of the problems.     

Numerical solutions of hypersingular integral equation for curved cracks in circular regions

In this paper, numerical solutions of a hypersingular integral equation for curved cracks in circular regions are presented. The boundary of the circular regions is assumed to be traction free or fixed. The suggested complex potential is composed of two parts, the principle part and the complementary part. The principle part can model the property of a curved crack in an infinite plate. For the case of the traction free boundary, the complementary part can compensate the traction on the circular boundary caused by the principle part. Physically, the proposed idea is similar to the image method in electrostatics. By using the crack opening displacement (COD) as the unknown function and traction as right hand term in the equation, a hypersingular integral equation for the curved crack problems in the circular regions is obtained. The equation is solved by using the curve length coordinate method. In order to prove that the suggested method can be used to solve more complicated cases of the curved cracks, several numerical examples are given.     

Numerical solution of the variation boundary integral equation for inverse problems

In this paper a procedure to solve the identification inverse problems for two‐dimensional potential fields is presented. The procedure relies on a boundary integral equation (BIE) for the variations of the potential, flux, and geometry. This equation is a linearization of the regular BIE for small changes in the geometry. The aim in the identification inverse problems is to find an unknown part of the boundary of the domain, usually an internal flaw, using experimental measurements as additional information. In this paper this problem is solved without resorting to a minimization of a functional, but by an iterative algorithm which alternately solves the regular BIE and the variation BIE. The variation of the geometry of the flaw is modelled by a virtual strainfield, which allows for greater flexibility in the shape of the assumed flaw. Several numerical examples demonstrate the effectiveness and reliability of the proposed approach. Copyright © 2000 John Wiley & Sons, Ltd.     

Numerical analysis of the inclusion-crack interactions using an integral equation

 A general-purpose integral formulation is proposed for the analysis of the interaction between inclusions and cracks embedded in an elastic isotropic homogeneous infinite medium subjected to a remote loading. This formulation is tailored for the inclusions of arbitrary shapes with the presence of cracks. The discretization is limited to the inclusions (with continuous quadratic triangular and quadrilateral elements) and the cracks (using discontinuous quadratic elements). For the calculation of the stress intensity factors at the crack tips, special crack tip elements are used to model the variation of the displacements near the crack tips. Maximum circumferential stress criterion is adopted to determine the crack propagating direction. Numerical results of benchmark examples are compared with other available methods. Received: 8 January 2002 / Accepted: 24 September 2002     

Modeling of interface cracks in fiber-reinforced composites with the presence of interphases using the boundary element method

An advanced boundary element method (BEM) with thin-body capabilities was developed recently for the study of interphases in fiber-reinforced composite materials (Y.J. Liu, N. Xu and J.F. Luo, Modeling of interphases in fiber-reinforced composites under transverse loading using the boundary element method, ASME J. Appl. Mech. 67 (2000) 41–49). In this BEM approach, the interphases are modeled as thin elastic layers based on the elasticity theory, as opposed to spring-like models in the previous BEM and some FEM work. In the present paper, this advanced BEM is extended to study the interface cracks at the interphases in the fiber-reinforced composites. These interface cracks are curved cracks between the fiber and matrix, with the presence of the interphases. Stress intensity factors (SIFs) for these interface cracks are evaluated based on the developed models. The BEM approach is validated first using available analytical and other numerical results for curved cracks in a single material and straight interface cracks between two materials. Then, the interface cracks at the interphases of fiber-reinforced composites are studied and the effects of the interphases (such as the thickness and materials) on the SIFs are investigated. As a special case, results of the SIFs for sub-interface cracks are also presented. It is shown that the developed BEM is very accurate and efficient for the interface crack analyses, and that the properties of the interphases have significant influences on the SIFs for interface cracks in fiber-reinforced composites.     

Wave propagation in the presence of empty cracks in an elastic medium

This paper proposes the use of a traction boundary element method (TBEM) to evaluate 3D wave propagation in unbounded elastic media containing cracks whose geometry does not change along one direction. The proposed formulation is developed in the frequency domain and handles the thin-body difficulty presented by the classical boundary element method (BEM). The empty crack may have any geometry and orientation and may even exhibit null thickness. Implementing this model yields hypersingular integrals, which are evaluated here analytically, thereby surmounting one of the drawbacks of this formulation. The TBEM formulation enables the crack to be modelled as a single line, allowing the computation of displacement jumps in the opposing sides of the crack. Furthermore, if this formulation is combined with the classical BEM formulation the displacements in the opposing sides of the crack can be computed by modelling the crack as a closed empty thin body.     

On solvability of a boundary integral equation of the first kind for Dirichlet boundary value problems in plane elasticity

The solution of a Dirichlet boundary value problem of plane isotropic elasticity by the boundary integral equation (BIE) of the first kind obtained from the Somigliana identity is considered. The logarithmic function appearing in the integral kernel leads to the possibility of this operator being non-invertible, the solution of the BIE either being non-unique or not existing. Such a situation occurs if the size of the boundary coincides with the so-called critical (or degenerate) scale for a certain form of the fundamental solution used. Techniques for the evaluation of these critical scales and for the removal of the non-uniqueness appearing in the problems with critical scales solved by the BIE of the first kind are proposed and analysed, and some recommendations for BEM code programmers based on the analysis presented are given.     

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.     

Numerical expansion-iterative method for analysis of integral equation models arising in one- and two-dimensional electromagnetic scattering

Most integral equations of the first kind are ill-posed, and obtaining their numerical solution needs often to solve a linear system of algebraic equations of large condition number. So, solving this system may be difficult or impossible. Since many problems in one- and two-dimensional scattering from perfectly conducting bodies can be modeled by Fredholm integral equations of the first kind, this paper presents an effective numerical expansion-iterative method for solving them. This method is based on vector forms of block-pulse functions. By using this approach, solving the first kind integral equation reduces to solve a recurrence relation. The approximate solution is most easily produced iteratively via the recurrence relation. Therefore, computing the numerical solution does not need to directly solve any linear system of algebraic equations and to use any matrix inversion. Also, the method practically transforms solving of the first kind Fredholm integral equation which is inherently ill-posed into solving second kind Fredholm integral equation. Another advantage is low cost of setting up the equations without applying any projection method such as collocation, Galerkin, etc. To show convergence and stability of the method, some computable error bounds are obtained. Test problems are provided to illustrate its accuracy and computational efficiency, and some practical one- and two-dimensional scatterers are analyzed by it.     

Transient heat conduction analysis in functionally graded materials by the meshless local boundary integral equation method

Advanced computational method for transient heat conduction analysis in continuously nonhomogeneous functionally graded materials (FGM) is proposed. The method is based on the local boundary integral equations with moving least square approximation of the temperature and heat flux. The initial-boundary value problem is solved by the Laplace transform technique. Both Papoulis and Stehfest algorithms are applied for the numerical Laplace inversion to obtain the time-dependent solutions. Numerical results are presented for a finite strip and a hollow cylinder with an exponential spatial variation of material parameters.