Computational model for short-fiber composites with eigenstrain formulation of boundary integral equations

A computational model is proposed for short-fiber reinforced materials with the eigenstrain formulation of the boundary integral equations （BIE） and solved with the newly developed boundary point method （BPM）. The model is closely derived from the concept of the equivalent inclusion Of Eshelby tensors. Eigenstrains are iteratively determined for each short-fiber embedded in the matrix with various properties via the Eshelby tensors, which can be readily obtained beforehand either through analytical or numerical means. As unknown variables appear only on the boundary of the solution domain, the solution scale of the inhomogeneity problem with the model is greatly reduced. This feature is considered significant because such a traditionally time-consuming problem with inhomogeneity can be solved most cost-effectively compared with existing numerical models of the FEM or the BEM. The numerical examples are presented to compute the overall elastic properties for various short-fiber reinforced composites over a representative volume element （RVE）, showing the validity and the effectiveness of the proposed computational modal and the solution procedure.     

Computational model for short-fiber composites with eigenstrain formulation of boundary integral equations

A computational model is proposed for short-fiber reinforced materials with the eigenstrain formulation of the boundary integral equations(BIE)and solved with the newly developed boundary point method(BPM).The model is closely derived from the concept of the equivalent inclusion of Eshelby tensors.Eigenstrains are iteratively determined for each short.fiber embedded in the matrix with various properties via the Eshelby tensors,which can be readily obtained beforehand either through analytical or numerical means.As unknown variables appear only on the boundary of the solution domain,the solution scale of the inhomogeneity problem with the model is greatly reduced.This feature is considered significant because such a traditionally time-consuming problem with inhomogeneity can be solved most cost-effectively compared with existing numerical models of the FEM or the BEM.The numerical examples are presented to compute the overall elastic properties for various short-fiber reinforced composites over a representative volume element(RVE),showing the validity and the effectiveness of the proposed computational modal and the solution procedure.     

Boundary element method for thermoelastic analysis of three-dimensional transversely isotropic solids

Thermal effects are well known to manifest themselves as additional volume integral terms in the direct formulation of the boundary integral equation (BIE) for linear elastic solids when using the boundary element method (BEM). This domain integral has been successfully transformed in an exact manner to surface ones only in isotropy and in 2D anisotropy, thereby restoring the BEM as a truly boundary solution technique. The difficulties with extending it to 3D general anisotropic solids lie in the mathematical complexity of the Green’s function and its derivatives for such materials. These quantities are required items in the BEM formulation. In this paper, the exact, analytical transformation of the volume integral associated with thermal effects to surface ones is achieved for a transversely isotropic material using a similar approach which the authors have previously employed for the same task in BEM for 2D general anisotropy. A numerical scheme, however, needs to be employed to evaluate some of the new terms introduced in the surface integrals that arise from this process here. The mathematical soundness of the formulation is demonstrated by a few examples; the numerical results obtained are checked by alternative means, including those obtained from the commercial FEM code, ANSYS.     

Boundary element analysis for elastic and elastoplastic problems of 2D orthotropic media with stress concentration

Both the orthotropy and the stress concentration are common issues in modern structural engineering. This paper introduces the boundary element method (BEM) into the elastic and elastoplastic analyses for 2D orthotropic media with stress concentration. The discretized boundary element formulations are established, and the stress formulae as well as the fundamental solutions are derived in matrix notations. The numerical procedures are proposed to analyze both elastic and elastoplastic problems of 2D orthotropic media with stress concentration. To obtain more precise stress values with fewer elements, the quadratic isoparametric element formulation is adopted in the boundary discretization and numerical procedures. Numerical examples show that there are significant stress concentrations and different elastoplastic behaviors in some orthotropic media, and some of the computational results are compared with other solutions. Good agreements are also observed, which demonstrates the efficiency and reliability of the present BEM in the stress concentration analysis for orthotropic media. The project supported by the Basic Research Foundation of Tsinghua University, the National Foundation for Excellent Doctoral Thesis (200025) and the National Natural Science Foundation of China (19902007). The English text was polished by Keren Wang.     

Multiple reciprocity boundary element analysis of two-dimensional anisotropic thermoelasticity involving an internal arbitrary non-uniform volume heat source

In the direct boundary element method (BEM) formulation of anisotropic thermoelasticity, thermal loads manifest themselves as additional volume integral terms in the boundary integral equation (BIE). Conventionally, this requires internal cell discretisation throughout the whole domain. In this paper, the multiple reciprocity method in BEM analysis is employed to treat the general 2D thermoelasticity problem when the thermal loading is due to an internal non-uniform volume heat source. By successively performing the “volume-to-surface” integral transformation, the general formulation of the associated BIE for the problem is derived. The successful implementation of such a scheme is illustrated by three numerical examples.     

Wideband fast multipole boundary element method: Application to acoustic scattering from aerodynamic bodies

A wideband adaptive multi‐level fast multipole method (MLFMM) is used to accelerate the matrix–vector products arising from a boundary element method (BEM) formulation which solves the Burton–Miller boundary integral equation (BIE). The wideband MLFMM presented here applies a plane wave expansion formulation with fast interpolation and filtering for calculations in the high‐frequency regime and a partial wave expansion formulation with rotation‐coaxial translation in the low‐frequency regime. The iterative solvers GMRES, Bi‐CGSTAB and CGS are tested and compared and a block diagonal preconditioner is used to improve the condition number of the BEM matrices and to accelerate the convergence of the iterative solvers. Details on the implementation of the formulations are described, including the treatment of singular integrals. Results for acoustic scattering from a wing plus engine nacelle configuration for a prescribed source in a subsonic uniform flow are presented for a broad range of frequencies in order to assess the implemented capability. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical comparison of two boundary element methods for plane harmonic functions

A direct boundary element method (BEM) has been studied in the paper based on a set of sufficient and necessary boundary integral equations (BIE) for the plane harmonic functions. The new sufficient and necessary BEM leads to accurate results while the conventional insufficient BEM will lead to inaccurate results when the conventional BIE has multiple solutions. Theoretical and numerical analyses show that it is beneficial to use the sufficient and necessary BEM, to avoid hidden dangers due to non-unique solution of the conventional BIE.     

Continuum Dyson's equation and defect Green's function in a heterogeneous anisotropic solid

A continuum Dyson's equation and a defect Green's function (GF) in a heterogeneous, anisotropic and linearly elastic solid under homogeneous boundary conditions have been introduced. The continuum Dyson's equation relates the point-force Green's responses of two systems of identical geometry and boundary conditions but of different media. Given the GF of either system (i.e., a reference), the GF of the other (i.e., a defect system with “defect” change of materials property relative to the reference) can be obtained by solving the Dyson's equation. The defect GF is applied to solve the eigenstrain problem of a heterogeneous solid. In particular, the problem of slightly inhomogeneous inclusions is examined in detail. Based on the Dyson's equation, approximate schemes are proposed to efficiently evaluate the elastic field. Numerical results are reported for inhomogeneous inclusions in a semi-infinite substrate with a traction-free surface to demonstrate the validity of the present formulation.     

Elastic field of an isotropic matrix with a nanoscale elliptical inhomogeneity

In traditional continuum mechanics, the effect of surface energy is ignored as it is small compared to the bulk energy. For nanoscale materials and structures, however, the surface effects become significant due to the high surface/volume ratio. In this paper, two-dimensional elastic field of a nanoscale elliptical inhomogeneity embedded in an infinite matrix under arbitrary remote loading and a uniform eigenstrain in the inhomogeneity is investigated. The Gurtin–Murdoch surface/interface elasticity model is applied to take into account the surface/interface stress effects. By using the complex variable technique of Muskhelishvili, the analytic potential functions are obtained in the form of an infinite series. Selected numerical results are presented to study the size-dependency of the elastic field and the effects of surface elastic moduli and residual surface stress. It is found that the elastic field of an elliptic inhomogeneity under uniform eigenstrain is no longer uniform when the interfacial stress effects are taken into account.     

弹性力学问题中一个新的边界积分方程——自然边界积分方程

位移导数边界积分方程一直存在着超奇异积分计算的障碍,该文提出以符号算子δye和εye作用于位移导数边界积分方程,施用一系列变换将边界位移、面力和位移导数转成为新的边界张量,从而得到一个新的边界积分方程－－自然边界积分方程,自然边界积分方奇异性为强奇性,文中给出了相应的Cauchy主值积分算式,自然边界积分方程与位移边界积分方程联合可直接获取边界应力,几个算例表明了自然边界积分方程的正确性。     

On the Solution of an Elliptical Inhomogeneity in Plane Elasticity by the Equivalent Inclusion Method

Stress analysis of an elliptical inhomogeneity in an infinite isotropic elastic plane is a classical elasticity problem, which is usually solved by means of the complex variable formulation. In this work, we demonstrate that an alternative method of solution for such a problem, via the equivalent inclusion method, may be more convenient and straightforward without recourse to complex potentials or curvilinear coordinates. The explicit analytical solution can be derived through simple algebraic manipulation, although the longitudinal eigenstrain component should be handled with care in the case of plane strain. Since the exterior Eshelby tensor for an elliptical inclusion is available in closed-form, the present study provides a full field stress solution expressed in Cartesian coordinates. Furthermore, the in-plane stress components are represented in terms of Dundurs’ parameters. The solution methodology and the convenient formulae of the stress concentration may be of practical use to the engineers in developing benchmarks for design evaluation.     

Time-harmonic BEM for 2-D piezoelectricity applied to eigenvalue problems

We will derive the fundamental generalized displacement solution, using the Radon transform, and present the direct formulation of the time-harmonic boundary element method (BEM) for the two-dimensional general piezoelectric solids. The fundamental solution consists of the static singular and the dynamics regular parts; the former, evaluated analytically, is the fundamental solution for the static problem and the latter is given by a line integral along the unit circle. The static BEM is a component of the time-harmonic BEM, which is formulated following the physical interpretation of Somigliana’s identity in terms of the fundamental generalized line force and dislocation solutions obtained through the Stroh–Lekhnitskii (SL) formalism. The time-harmonic BEM is obtained by adding the boundary integrals for the dynamic regular part which, from the original double integral representation over the boundary element and the unit circle, are reduced to simple line integrals along the unit circle.The BEM will be applied to the determination of the eigen frequencies of piezoelectric resonators. The eigenvalue problem deals with full non-symmetric complex-valued matrices whose components depend non-linearly on the frequency. A comparative study will be made of non-linear eigenvalue solvers: QZ algorithm and the implicitly restarted Arnoldi method (IRAM). The FEM results whose accuracy is well established serve as the basis of the comparison. It is found that the IRAM is faster and has more control over the solution procedure than the QZ algorithm. The use of the time-harmonic fundamental solution provides a clean boundary only formulation of the BEM and, when applied to the eigenvalue problems with IRAM, provides eigen frequencies accurate enough to be used for industrial applications. It supersedes the dual reciprocity BEM and challenges to replace the FEM designed for the eigenvalue problems for piezoelectricity.     

A boundary element application for mixed modeloading idealized sawtooth fracture surface

In this paper, a 2-D elastic-plastic BEM formulation predicting the reduced mode IIand the enhanced mode I stress intensity factors are presented. The dilatant boundary conditions (DBC) are assumed to be idealized uniform sawtooth crack surfaces and an effective Coulombsliding law. Three types of crack face boundary conditions, i.e. (1) BEM sawtooth model-elasticcenter crack tip; (2) BEM sawtooth model-plastic center crack tip; and (3) BEM sawtoothmodel-edge crack with asperity wear are presented. The model is developed to attempt todescribe experimentally observed non-monotonic, non-linear dependence of shear crack behavioron applied shear stress, superimposed tensile stress, and crack length. The asperity sliding isgoverned by Coulombs law of friction applied on the inclined asperity surface which hascoefficient of friction μ. The traction and displacement Greens functions which derive fromNaviers equations are obtained as well as the governing boundary integral equations for an infiniteelastic medium. Accuracy test is performed by comparison stress intensity factors of the BEMmodel with analytical solutions of the elastic center crack tip. The numerical results show thepotential application of the BEM model to investigate the effect of mixed mode loading problemswith various boundary conditions and physical interactions.     

Development of plasticity and damage in vibrating structural elements performing guided rigid-body motions

Summary  A numerical algorithm for studying the development of plastic and damaged zones in a vibrating structural element with a large, guided rigid-body motion is presented. Beam-type elements vibrating in the small-strain regime are considered. A machine element performing rotatory motions, similar to an element of a slider-crank mechanism, is treated as a benchmark problem. Microstructural changes in the deforming material are described by the mesolevel variables of plastic strain and damage, which are consistently included into a macroscopic analysis of the overall beam motion. The method is based on an eigenstrain formulation, considering plastic strain and damage to contribute to an eigenstrain loading of a linear elastic background structure. Rigid-body coordinates are incorporated into this beam-type structural formulation, and an implicit numerical scheme is presented for iterative computation of the eigenstrains from the mesolevel constitutive behavior. Owing to the eigenstrain formulation, any of the existing constitutive models with internal variables could in principle be implemented. Linear elastic/perfectly plastic behavior is exemplarily treated in a numerical study, where plastic strain is connected to the Kachanov damage parameter by a simple damage law. Inelastic effects like plastic shakedown and damage-induced low-cycle rupture are shown to occur in the examplary problems. Received 1 September 1999; accepted for publication 9 March 2000     

A NOVEL BOUNDARY INTEGRAL EQUATION METHOD FOR LINEAR ELASTICITY--NATURAL BOUNDARY INTEGRAL EQUATION

The boundary integral equation (BIE) of displacement derivatives is put at a disadvantage for the difficulty involved in the evaluation of the hypersingular integrals. In this paper, the operators δij and εij are used to act on the derivative BIE. The boundary displacements, tractions and displacement derivatives are transformed into a set of new boundary tensors as boundary variables. A new BIE formulation termed natural boundary integral equation (NBIE) is obtained. The NBIE is applied to solving two-dimensional elasticity problems. In the NBIE only the strongly singular integrals are contained. The Cauchy principal value integrals occurring in the NBIE are evaluated. A combination of the NBIE and displacement BIE can be used to directly calculate the boundary stresses. The numerical results of several examples demonstrate the accuracy of the NBIE.     

A boundary element method for calculating the SH waves scattering from arbitrarily shaped holes in an anisotropic medium

The scattering problem of elastic wave by arbitrarily shaped cavities in an infinite anisotropic medium is investigated by the boundary integral equation (BIE) method. The formulations of BIE are derived with the help of generalized Green's formula. The discretization of BIE is based upon constant elements. After confirmation of the accuracy of the present method, some numerical examples are given for various cavities in a full space, in which an isotropic body with a circular cylinder hole is used for comparison and good agreement is observed. It has been proved that the method developed in this paper is effective.     

Degenerate scale for multiply connected Laplace problems

The degenerate scale in the boundary integral equation (BIE) or boundary element method (BEM) solution of multiply connected problem is studied in this paper. For the mathematical analysis, we use the null-field integral equation, degenerate kernels and Fourier series to examine the solvability of BIE for multiply connected problem in the discrete system. Two treatments, the method of adding a rigid body term and CHEEF concept (Combined Helmholtz Exterior integral Equation Formulation), are applied to remedy the non-unique solution due to the critical scale. The efficiency and accuracy of the two regularizations are also addressed. For simplicity without loss of generality, the eccentric case is considered to demonstrate the occurring mechanism of degenerate scale.     

各向异性体内含任意孔洞对反平面波散射的边界元方法

本文借助于广义格林公式导出了用位移表示的各向异性介质中SH波入射时的边界积分方程.根据本文作者在文献[8]给出的基本解,求解了各向异性介质中孔洞对SH波的散射问题.边界积分方程的离散基于常数元模式.文中给出了一个圆柱、一个椭圆柱和两个椭圆柱形式的孔洞周围的位移场和应力场的数值结果.最后,对入射波频率较高时的情形作了说明.     

Analysis of FRP composites under frictional contact conditions

This work studies numerically the tribological behavior of fiber-reinforced plastics (FRP) under different frictional contact conditions, using a boundary element methodology. The formulation uses the Boundary Element Method (BEM) with an explicit approach for fundamental solutions evaluation, for computing the elastic influence coefficients. To enforce the contact constraints on the potential contact zone: Signorini’s contact conditions and an orthotropic law of friction, contact operators over the augmented Lagrangian are considered in the formulation. The methodology and the proposed algorithm are applied to study two types of glass FRP and two types of carbon FRP, with the same fiber volume fraction, under frictional contact. In these studies, it can be observed how the fiber orientation and sliding orientation affect the normal and tangential contact compliance, as well as the contact traction distribution. Furthermore, the formulation considers a micromechanics model for FRP that allows also to study the influence of fiber volume fraction on normal and tangential contact compliances.