Search results

Domain integrals, known as volume potentials in 3D elasticity problems, exist in many boundary-type methods, such as the boundary element method (BEM) for inhomogeneous partial differential equations. The purpose of this paper is to develop an accurate and reliable technique to effectively evaluate the volume potentials in 3D elasticity problems.

An adaptive background cell-based domain integration method is proposed for treatment of volume potentials in 3D elasticity problems. The background cells are constructed from the information of the boundary elements based on an oct-tree structure, and the domain integrals are evaluated over the cells rather than volume elements. The cells that contain the boundary elements can be subdivided into smaller sub-cells adaptively according to the sizes and levels of the boundary elements. The fast multipole method (FMM) is further applied in the proposed method to reduce the time complexity of large-scale computation.

The method is a boundary-only discretization method, and it can be applied in the BEM easily. Much computational time is saved by coupling with the FMM. Numerical examples demonstrate the accuracy and efficiency of the proposed method..

Boundary elements are used to create adaptive background cells, and domain integrals are evaluated over the cells rather than volume elements. Large-scale computation is made possible by coupling with the FMM.

Gives introductory remarks about chapter 1 of this group of 31 papers, from ISEF 1999 Proceedings, in the methodologies for field analysis, in the electromagnetic community. Observes that computer package implementation theory contributes to clarification. Discusses the areas covered by some of the papers ‐ such as artificial intelligence using fuzzy logic. Includes applications such as permanent magnets and looks at eddy current problems. States the finite element method is currently the most popular method used for field computation. Closes by pointing out the amalgam of topics.

The boundary element method is a useful method for the analysis of field problems involving unbounded regions. Therefore, the method can be used advantageously in combination with the finite element method. This is sometimes called a combination method and it is suitable as a picture‐frame technique. Although this technique attains good accuracy, the matrix of the discretized equation is not banded, since it is a dense matrix. In this paper, we propose an infinite boundary element which divides the unbounded region radially. By the use of this element, the bandwidth of the discretized system matrix does not increase beyond that of the finite element region and its original matrix structure is maintained. The infinite boundary element can also be applied to homogeneous unbounded field problems, for which the Green's function of the mirror image is difficult to use. To illustrate the validity of the proposed technique, some numerical calculations are demonstrated and the results are compared with those of the usual combination method and the method using the hybrid‐type infinite element.

This paper aims to propose a boundary element analysis of two-dimensional linear elasticity problems by a new expanding element interpolation method.

The expanding element is made up based on a traditional discontinuous element by adding virtual nodes along the perimeter of the element. The internal nodes of the original discontinuous element are referred to as source nodes and its shape function as raw shape function. The shape functions of the expanding element constructed on both source nodes and virtual nodes are referred as fine shape functions. Boundary variables are interpolated by the fine shape functions, while the boundary integral equations are collocated on source nodes.

The expanding element inherits the advantages of both the continuous and discontinuous elements while overcomes their disadvantages. The polynomial order of fine shape functions of the expanding elements increases by two compared with their corresponding raw shape functions, while the expanding elements still keep independence to each other as the original discontinuous elements. This feature makes the expanding elements able to naturally and accurately interpolate both continuous and discontinuous fields.

Numerical examples are presented to verify the proposed method. Results have demonstrated that the accuracy, efficiency and convergence rate of the expanding element method.

When explicit integral analysis is performed on a numerical model with viscoelastic artificial boundary elements, an instability phenomenon is likely to occur in the boundary area, reducing the computational efficiency of the numerical calculation and limiting the use of viscoelastic artificial boundary elements in the explicit dynamic analysis of large-scale engineering sites. The main purpose of this study is to improve the stability condition of viscoelastic artificial boundary elements.

A stability analysis method based on local subsystems was adopted to analyze and improve the stability conditions of three-dimensional (3D) viscoelastic artificial boundary elements. Typical boundary subsystems that can represent the localized characteristics of the overall model were established, and their analytical stability conditions were derived with an analysis based on the spectral radius of the transfer matrix. Then, after analyzing the influence of each physical parameter on the analytical-stability conditions, a method for improving the stability condition of the explicit algorithm by increasing the mass density of the artificial boundary elements was proposed.

Numerical wave propagation simulations in uniform and layered half-space models show that, on the premise of ensuring the accuracy of the viscoelastic artificial boundary, the proposed method can effectively improve the numerical stability and the efficiency of the explicit dynamic calculations for the overall system.

The stability improvement method proposed in this study are significant for improving the applicability of viscoelastic artificial boundary elements in explicit dynamic calculations and the calculation efficiency of wave analysis at large-scale engineering sites.

A technique combining finite elements and boundary elements is promising for unbounded field problems. A hypothetical boundary is assumed in the unbounded domain, and the usual finite element method is applied to the inner region, while the boundary element method is applied to the outer infinite region. On the coupling boundary, therefore, both potential and flux must be compatible. In the finite element method, the flux is defined as the derivative of the potential for which a trial function is defined. In the boundary element method, on the other hand, the same polynomial function is chosen for the potential and the flux. Thus, the compatibility cannot be satisfied unless a special device is considered. In the present paper, several compatibility conditions are discussed concerning the total flux or energy flow continuity across the coupling boundary. Some numerical examples of Poisson and Helmholtz problems are demonstrated.

The spectral boundary element method for solving two‐dimensional transient heat conduction problems is developed. This method is combined with the fast Fourier transform (FFT) to convert the solution between the time and frequency domains. The fundamental solutions in the frequency domain, required for the present method, are discussed. The resulting line integrations in the frequency domain are discretized using constant boundary elements and used in a Fortran boundary element program. Three examples are used to illustrate the accuracy and effectiveness of the method in both the frequency and time domains. First, the frequency domain solution procedure is verified using the steady‐state example of a semi‐infinite half space with a heat flux applied to a patch of the surface. This spectral boundary element method is then applied to the problem of a circular hole in an infinite solid subjected to a time‐varying heat flux, and solutions in both the frequency and time domains are presented. Finally, the method is used to solve the circular hole problem with a convection boundary condition. The accurary of these results leads to the conclusion that the spectral boundary element method in conjunction with the FFT is a viable option for transient problems. In addition, this spectral approach naturally produces frequence domain information which is itself of interest.

The purpose of this paper is to evaluate the efficiency of the fast multipole boundary element method (FMBEM) in the analysis of stress and effective properties of 3D linear elastic structures with cavities. In particular, a comparison between the FMBEM and the finite element method (FEM) is performed in terms of accuracy, model size and computation time.

The developed FMBEM uses eight-node Serendipity boundary elements with numerical integration based on the adaptive subdivision of elements. Multipole and local expansions and translations involve solid harmonics. The proposed model is used to analyse a solid body with two interacting spherical cavities, and to predict the homogenized response of a porous material under linear displacement boundary condition. The FEM results are generated in commercial codes Ansys and MSC Patran/Nastran, and the results are compared in terms of accuracy, model size and execution time. Analytical solutions available in the literature are also considered.

FMBEM and FEM approximate the geometry with similar accuracy and provide similar results. However, FMBEM requires a model size that is smaller by an order of magnitude in terms of the number of degrees of freedom. The problems under consideration can be solved by using FMBEM within the time comparable to the FEM with an iterative solver.

The present results are limited to linear elasticity.

This work is a step towards a comprehensive efficiency evaluation of the FMBEM applied to selected problems of micromechanics, by comparison with the commercial FEM codes.

The purpose of this paper is to develop an easy-to-implement and accurate fast boundary element method (BEM) for solving large-scale elastodynamic problems in frequency and time domains.

A newly developed kernel-independent fast multipole method (KIFMM) is applied to accelerating the evaluation of displacements, strains and stresses in frequency domain elastodynamic BEM analysis, in which the far-field interactions are evaluated efficiently utilizing equivalent densities and check potentials. Although there are six boundary integrals with unique kernel functions, by using the elastic theory, the authors managed to accelerate these six boundary integrals by KIFMM with the same kind of equivalent densities and check potentials. The boundary integral equations are discretized by Nyström method with curved quadratic elements. The method is further used to conduct the time-domain analysis by using the frequency-domain approach.

Numerical results show that by the fast BEM, high accuracy can be achieved and the computational complexity is brought down to linear. The performance of the present method is further demonstrated by large-scale simulations with more than two millions of unknowns in the frequency domain and one million of unknowns in the time domain. Besides, the method is applied to the topological derivatives for solving elastodynamic inverse problems.

An efficient KIFMM is implemented in the acceleration of the elastodynamic BEM. Combining with the Nyström discretization based on quadratic elements and the frequency-domain approach, an accurate and highly efficient fast BEM is achieved for large-scale elastodynamic frequency domain analysis and time-domain analysis.

The penalty boundary method (PBM) is a new method for performing finite element analysis using a regular overlapping mesh that does not have to coincide with the geometric boundaries. The PBM uses CAD solid geometry directly to generate element matrix equations and apply boundary conditions, removing the need for a separate representation of the geometry. The preliminary results show that the PBM can significantly reduce the time and manual intervention required to prepare finite element models and perform analyses. This paper presents the PBM approach for representing the problem domain on an overlapping mesh that results in a more traditional method for applying natural boundary conditions.