Determination of stress intensity factors for bimaterial interface stationary rigid line inclusions by boundary element method

The stress intensity factors for a rigid line inclusion lying along a bimaterial interface are calculated by the boundary element method with the multiregion and the discontinuous traction singular elements. The relationships between the stress intensity factors and the inclusion surface stresses are derived. The numerically computed stress intensity factors for the bimaterial interface rigid line inclusion in the infinite body are proved to be in good agreement within 3% when compared with the previous exact solutions. In the finite bimaterial models, the stress intensity factors for the center and edge rigid line inclusions at the interface are computed with the variation of the rigid line inclusion length and the shear modulus ratio under the uniaxial and biaxial loading conditions.     

An iterative boundary element method for Cauchy inverse problems

 We are interested in this paper in recovering an harmonic function from the knowledge of Cauchy data on some part of the boundary. A new inversion method is introduced. It reduces the Cauchy problem resolution to the determination of the resolution of a sequence of well-posed problems. The sequence of these solutions is proved to converge to the Cauchy problem solution. The algorithm is implemented in the framework of boundary elements. Displayed numerical results highlight its accuracy, as well as its robustness to noisy data. Received 6 November 2000     

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.     

A fast convergent iterative boundary element method on PVM cluster

This paper reports a fast convergent boundary element method on a Parallel Virtual Machine (PVM) (Geist et al., PVM: Parallel Virtual Machine, A Users' Guide and Tutorial for Networked Parallel Computing. MIT Press, Cambridge, 1994) cluster using the SIMD computing model (Single Instructions Multiple Data). The method uses the strategy of subdividing the domain into a number of smaller subdomains in order to reduce the size of the system matrix and to achieve overall speedup. Unlike traditional subregioning methods, where equations from all subregions are assembled into a single linear algebraic system, the present scheme is iterative and each subdomain is handled by a separate PVM node in parallel. The iterative nature of the overall solution procedure arises due to the introduction of the artificial boundaries. However, the system equations for each subdomain is now smaller and solved by direct Gaussian elimination within each iteration. Furthermore, the boundary conditions at the artificial interfaces are estimated from the result of the previous iteration by a reapplication of the boundary integral equation for internal points. This method provides a consistent mechanism for the specification of boundary conditions on artificial interfaces, both initially and during the iterative process. The method is fast convergent in comparison with other methods in the literature. The achievements of this method are therefore: (a) simplicity and consistency of methodology and implementation; (b) more flexible choice of type of boundary conditions at the artificial interfaces; (c) fast convergence; and (d) the potential to solve large problems on very affordable PVM clusters. The present parallel method is suitable where (a) one has a distributed computing environment; (b) the problem is big enough to benefit from the speedup achieved by coarse-grained parallelisation; and (c) the subregioning is such that communication overhead is only a small percentage of total computation time.     

Acoustic problems analysis of 3D solid with small holes by fast multipole boundary face method

In this paper, an adaptive fast multipole boundary face method is introduced to implement acoustic problems analysis of 3D solids with open-end small tubular shaped holes. The fast multipole boundary face method is referred as FMBFM. These holes are modeled by proposed tube elements. The hole is open-end and intersected with the outer surface of the body. The discretization of the surface with circular inclusions is achieved by applying several special triangular elements or quadrilateral elements. In the FMBFM, the boundary integration and field variables approximation are both performed in the parametric space of each boundary face exactly the same as the B-rep data structure in standard solid modeling packages. Numerical examples for acoustic radiation in this paper demonstrated the accuracy, efficiency and validity of this method.     

A complex variable boundary element method for solving interface crack problems

The problem of finite bimaterial plates with an edge crack along the interface is studied. A complex variable boundary element method is presented and applied to determine the stress intensity factor for finite bimaterial plates. Using the pseudo-orthogonal characteristic of the eigenfunction expansion forms and the well-known Bueckner work conjugate integral and taking the different complex potentials as auxiliary fields, the interfacial stress intensity factors associated with the physical stress-displacement fields are evaluated. The effects of material properties and crack geometry on stress intensity factors are investigated. The numerical examples for three typical specimens with six different combinations of the bimaterial are given. This revised version was published online in July 2006 with corrections to the Cover Date.     

Radial integration boundary element method for acoustic eigenvalue problems

In this paper, the radial integration boundary element method is developed to solve acoustic eigenvalue problems for the sake of eliminating the frequency dependency of the coefficient matrices in traditional boundary element method. The radial integration method is presented to transform domain integrals to boundary integrals. In this case, the unknown acoustic variable contained in domain integrals is approximated with the use of compactly supported radial basis functions and the combination of radial basis functions and global functions. As a domain integrals transformation method, the radial integration method is based on pure mathematical treatments and eliminates the dependence on particular solutions of the dual reciprocity method and the particular integral method. Eventually, the acoustic eigenvalue analysis procedure based on the radial integration method resorts to a generalized eigenvalue problem rather than an enhanced determinant search method or a standard eigenvalue analysis with matrices of large size, just like the multiple reciprocity method. Several numerical examples are presented to demonstrate the validity and accuracy of the proposed approach.     

A semi-analytic boundary element method for parabolic problems

A new semi-analytic solution method is proposed for solving linear parabolic problems using the boundary element method. This method constructs a solution as an eigenfunction expansion using separation of variables. The eigenfunctions are determined using the dual reciprocity boundary element method. This separation of variables-dual reciprocity method (SOV-DRM) allows a solution to be determined without requiring either time-stepping or domain discretisation. The accuracy and computational efficiency of the SOV-DRM is found to improve as time increases. These properties make the SOV-DRM an attractive technique for solving parabolic problems.     

A new boundary element method for resonance of an acoustic enclosure

Abstract

This paper presents a new boundary element formulation in which the eigenvalue appears outside the integral operator, which distinguishes it from the Helmholtz integral equation. Thus, the formation of global matrices need only be assembled once. Since the kernel of the operator used in the new formulation is real‐valued, all calculations can be carried out in a much simpler way in the real domain. The complex acoustic pressure amplitude is considered herein to deivate by a certain amount from a harmonic function. It is an important contribution that an exact relation between the deviator and the complex acoustic pressure amplitude is constructed locally and thus no more approximations are introduced except conventional boundary discretizations. Several examples are given to illustrate the feasibility of an accurate, effective prediction of resonance.     

Evolutionary optimization in thermoelastic problems using the boundary element method

 The paper is devoted to application of evolutionary algorithms and the boundary element method to shape optimization of structures for various thermomechanical criteria, inverse problems of finding an optimal distribution of temperature on the boundary and identification of unknown boundary. Design variables are specified by Bezier curves. Several numerical examples of evolutionary computation are presented. Received 6 November 2000     

A general boundary element method approach to the solution of three-dimensional frictionless contact problems

A direct boundary element method is presented for three-dimensional stress analysis of frictionless contact problems. The isoparametric formulation of the boundary element method is implemented for the general case of contact in the absence of friction, which is limited to linear elastic homogeneous and isotropic materials. An iterative procedure is employed to determine the correct size of the contact zone by finding a boundary solution compatible with the contact condition. The applicability of the procedure is tested by application to three problems of advancing and conforming contact. The computed results are compared with numerical and analytical solutions where possible.     

An alternative collocation boundary element method for static and dynamic problems

A collocation boundary element formulation is presented which is based on a mixed approximation formulation similar to the Galerkin boundary element method presented by Steinbach (SIAM J Numer Anal 38:401–413, 2000) for the solution of Laplace’s equation. The method is also applicable to vector problems such as elasticity. Moreover, dynamic problems of acoustics and elastodynamics are included. The resulting system matrices have an ordered structure and small condition numbers in comparison to the standard collocation approach. Moreover, the employment of Robin boundary conditions is easily included in this formulation. Details on the numerical integration of the occurring regular and singular integrals and on the solution of the arising systems of equations are given. Numerical experiments have been carried out for different reference problems. In these experiments, the presented approach is compared to the common nodal collocation method with respect to accuracy, condition numbers, and stability in the dynamic case.     

Application of the boundary element method to three-dimensional potential problems in heterogeneous media

The present work discusses a solution procedure for heterogeneous media three-dimensional potential problems, involving nonlinear boundary conditions. The problem is represented mathematically by the Laplace equation and the adopted numerical technique is the boundary element method (BEM), here using velocity correcting fields to simulate the conductivity variation of the domain. The integral equation is discretized using surface elements for the boundary integrals and cells, for the domain integrals. The adopted strategy subdivides the discretized equations in two systems: the principal one involves the calculation of the potential in all boundary nodes and the secondary which determines the correcting field of the directional derivatives of the potential in all points. Comparisons with other numerical and analytical solutions are presented for some examples.     

Reliability analysis of plane elasticity problems by stochastic spline fictitious boundary element method

In this paper a stochastic spline fictitious boundary element method (SFBEM) is proposed for reliability analysis of plane elasticity problems in conjunction with the advanced first-order-second-moment (AFOSM) method. The AFOSM method has been demonstrated to be a reliable and practical approach to the structural reliability analysis, yielding results of reasonable accuracy for the engineering applications. And as a modified method for the conventional indirect boundary element method, SFBEM can provide accurate numerical solutions at high efficiency in deterministic analyses. For the purpose of structural reliability analysis, SFBEM is introduced during the iteration process of the AFOSM method, to obtain the required values of structural responses and their derivations with respect to the random variables considered. The use of SFBEM in the formulation of the AFOSM method makes it unnecessary to construct an explicit expression to the implicit limit state function of the problem, leading to a higher efficiency and better accuracy. The present approach is validated by comparing calculated solutions with those of Monte Carlo simulation for a number of example problems and a good agreement of the results is achieved.     

Stochastic spline fictitious boundary element method in elastostatic problems with random fields

Mathematical formulation and computational implementation of the stochastic spline fictitious boundary element method (SFBEM) are presented for the analysis of plane elasticity problems with material parameters modeled with random fields. Two sets of governing differential equations with respect to the means and deviations of structural responses are derived by including the first order terms of deviations. These equations, being in similar forms to those of deterministic elastostatic problems, can be solved using deterministic fundamental solutions. The calculation is conducted with SFBEM, a modified indirect boundary element method (IBEM), resulting in the means and covariances of responses. The proposed method is validated by comparing the solutions obtained with Monte Carlo simulation for a number of example problems and a good agreement of results is observed.     

Fracture mechanics analysis of cracking in plain and reinforced concrete using the boundary element method

Boundary element formulations for modelling the nonlinear behaviour of concrete are reviewed. The analysis based on the dual boundary element method (BEM) to represent the cracks in concrete is presented. The fictitious crack is adopted to represent the fracture process zone in concrete. The influence of reinforcements on the concrete is considered as a distribution of forces over the region of attachment. The yielding of reinforcement is considered when the total force at any section of the reinforcement is greater than the yielding force and is assumed to be broken when the strain reaches the maximum strain. In using the BEM to simulate cracks, the crack path need not be known in advance since it can be calculated during the iteration process and as such remeshing becomes obsolete. The numerical results obtained are compared to the FEM analysis.     

A hybrid boundary element method for three-dimensional fracture analysis

At first, a hybrid boundary element method used for three-dimensional linear elastic fracture analysis is established by introducing the relative displacement fundamental function into the first and the second kind of boundary integral equations. Then the numerical approaches are presented in detail. Finally, several numerical examples are given out to check the proposed method. The numerical results show that the hybrid boundary element method has a very high accuracy for analysis of a three-dimensional stress intensity factor.     

Dual boundary element method for three-dimensional thermoelastic crack problems

This paper describes the formulation and numerical implementation of the three-dimensional dual boundary element method (DBEM) for the thermoelastic analysis of mixed-mode crack problems in linear elastic fracture mechanics. The DBEM incorporates two pairs of independent boundary integral equations; namely the temperature and displacement, and the flux and traction equations. In this technique, one pair is applied on one of the crack faces and the other pair on the opposite one. On non-crack boundaries, the temperature and displacement equations are applied. This revised version was published online in July 2006 with corrections to the Cover Date.     

Automated simulation of fatigue crack propagation for two-dimensional linear elastic fracture mechanics problems by boundary element method

In this paper, automated simulation of multiple crack fatigue propagation for two-dimensional (2D) linear elastic fracture mechanics (LEFM) problems is developed by using boundary element method (BEM). The boundary element method is the displacement discontinuity method with crack-tip elements proposed by the author. Because of an intrinsic feature of the boundary element method, a general growth problem of multiple cracks can be solved in a single-region formulation. In the numerical simulation, for each increment of crack extension, remeshing of existing boundaries is not necessary. Local discretization on the incremental crack extension is performed easily. Further the new adding elements and the existing elements on the existing boundaries are employed to construct easily the total structural mesh representation. Here, the mixed-mode stress intensity factors are calculated by using the formulas based on the displacement fields around crack tip. The maximum circumferential stress theory is used to predict crack stability and direction of propagation at each step. The well-known Paris’ equation is extended to multiple crack case under mixed-mode loadings. Also, the user does not need to provide a desired crack length increment at the beginning of each simulation. The numerical examples are included to illustrate the validation of the numerical approach for fatigue growth simulation of multiple cracks for 2D LEFM problems.     

A boundary element method for the solution of a class of heat conduction problems for inhomogeneous media

A boundary element method is derived for solving the two-dimensional heat equation for an inhomogeneous body subject to suitably prescribed temperature and/or heat flux on the boundary of the solution domain. Numerical results for a specific test problem is given.