Search

Maxwell's methodology is combined with the boundary element method (BEM) for evaluating the two-dimensional effective elastic properties of composite, porous, and microcracked isotropic materials with periodic or random structure. The approach is based on the idea that the effective properties of the material can be deduced from the effects that a cluster of fibers, pores, or cracks embedded in an infinite matrix has on the far-fields. The fibers, pores, or cracks can have arbitrary shapes, sizes, and elastic properties, provided that the overall behavior is isotropic, and their effects are assumed to be the same as those of an equivalent circular inhomogeneity. The key aspect of the approach is to precisely account for the interactions between all the constituents in the cluster that represent the material in question. This is done by using the complex-variables version of the BEM to solve the problem of a finite cluster of fibers, pores or cracks embedded an infinite isotropic, linearly elastic matrix. The effective properties of the material are evaluated by comparing the far-field solutions for the cluster with that of the equivalent inhomogeneity. It is shown that the model adequately captures the influence of the micro-structure of the material on its overall properties.

This paper presents a numerical method for solving the two-dimensional problem of a polygonal linear viscoelastic domain containing an arbitrary number of non-overlapping circular holes of arbitrary sizes. The solution of the problem is based on the use of the correspondence principle. The governing equation for the problem in the Laplace domain is a complex hypersingular boundary integral equation written in terms of the unknown transformed displacements on the boundaries of the holes and the exterior boundaries of the finite body. No specific physical model is involved in the governing equation, which means that the method is capable of handling a variety of viscoelastic models. A truncated complex Fourier series with coefficients dependent on the transform parameter is used to approximate the unknown transformed displacements on the boundaries of the holes. A truncated complex series of Chebyshev polynomials with coefficients dependent on the transform parameter is used to approximate the unknown transformed displacements on the straight boundaries of the finite body. A system of linear algebraic equations is formed using the overspecification method. The viscoelastic stresses and displacements are calculated through the viscoelastic analogs of the Kolosov–Muskhelishvili potentials, and an analytical inverse Laplace transform is used to provide the time domain solution. Using the concept of representative volume, the effective viscoelastic properties of an equivalent homogeneous material are then found directly from the corresponding constitutive equations for the average field values.

A fast and accurate algorithm is presented to increase the computational efficiency of a Galerkin boundary integral method for solving two-dimensional elastostatics problems involving numerous straight cracks and circular inhomogeneities. The efficiency is improved by computing the combined influences of groups, or blocks, of elements—with each element being an inclusion, a hole, or a crack—using asymptotic expansions, multiple shifts, and Taylor series expansions. The coefficients in the asymptotic and Taylor series expansions are computed analytically. Implementation of this algorithm involves a single- or multi-level grid, a clustering technique, and a tree data structure. An iterative procedure is adopted to solve the coefficients in the series expansions of boundary unknowns block by block. The elastic fields in each block are calculated by superposition of the direct influences from the nearby elements and the grouped far-field influences from all the other elements. This fast multipole algorithm is considerably more efficient for large-scale practical problems than the conventional approach.

A direct boundary integral method in the time domain is presented to solve the problem of an infinite, isotropic Boltzmann viscoelastic plane containing a large number of randomly distributed, non-overlapping circular holes and perfectly bonded elastic inclusions. The holes and inclusions are of arbitrary size and the elastic properties of all of the inclusions can, in general, be different. The method is based on a direct boundary integral approach for the problem of an infinite elastic plane containing multiple circular holes and elastic inclusions described by Crouch and Mogilevskaya [1], and a time marching strategy for viscoelastic analysis described by Mesquita and Coda [2–8]. Benchmark problems and numerical examples are included to demonstrate the accuracy and efficiency of the method.

This paper describes a numerical procedure for solving two-dimensional elastostatics problems with multiple circular holes and elastic inclusions in a finite domain with a circular boundary. The inclusions may have arbitrary elastic properties, different from those of the matrix, and the holes may be traction free or loaded with uniform normal pressure. The loading can be applied on all or part of the finite external boundary. Complex potentials are expressed in the form of integrals of the tractions and displacements on the boundaries. The unknown boundary tractions and displacements are approximated by truncated complex Fourier series. A linear algebraic system is obtained by using Taylor series expansion without boundary discretization. The matrix of the linear system has diagonal submatrices on its diagonal, which allows the system to be effectively solved by using a block Gauss-Seidel iterative algorithm.