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.

A two‐grid iterative method for 3D linear elasticity problems, discretized using quadratic tetrahedral elements is proposed. The conjugate‐gradient method is used as smoother. As compared to the conjugate‐gradient alone, it is shown, via numerical examples, that the method is much more efficient on the basis of computing time and memory allocation. The convergence property of the method is sensitive to the regularity of the problem.

This paper gives a bibliographical review of the finite element and boundary element parallel processing techniques from the theoretical and application points of view. Topics include: theory – domain decomposition/partitioning, load balancing, parallel solvers/algorithms, parallel mesh generation, adaptive methods, and visualization/graphics; applications – structural mechanics problems, dynamic problems, material/geometrical non‐linear problems, contact problems, fracture mechanics, field problems, coupled problems, sensitivity and optimization, and other problems; hardware and software environments – hardware environments, programming techniques, and software development and presentations. The bibliography at the end of this paper contains 850 references to papers, conference proceedings and theses/dissertations dealing with presented subjects that were published between 1996 and 2002.

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.

Based on the recent advances of hybrid stress finite elements, a series of alternative stress assumptions for these elements are investigated. Several new element models are proposed by using different concepts for the stress interpolation. Under a unified formulation presented in this paper for Hellinger—Reissner principle based hybrid stress element models, the element series 5β‐family for plane stress and 18β‐family for three‐dimensional problems are discussed. The extra incompatible displacements sometimes also added are not introduced in this unified formulation. A number of popular benchmark elastic problems are examined for both two element families. In each family, the element model presented in this paper using normalized transformed higher order stress trials usually gives better predictions than the others.

The purpose of this paper is to present special nine‐node quadrilateral elements to discretize the un‐cracked boundary and the inclined surface crack in a transversely isotropic cuboid under a uniform vertical traction along its top and bottom surfaces by a three‐dimensional (3D) boundary element method (BEM) formulation. The mixed‐mode stress intensity factors (SIFs), K_I, K_II and K_III, are calculated.

A 3D dual‐BEM or single‐domain BEM is employed to solve the fracture problems in a linear anisotropic elastic cuboid. The transversely isotropic plane has an arbitrary orientation, and the crack surface is along an inclined plane. The mixed 3D SIFs are evaluated by using the asymptotical relation between the SIFs and the relative crack opening displacements.

Numerical results show clearly the influence of the material and crack orientations on the mixed‐mode SIFs. For comparison, the mode‐I SIF when a horizontal rectangular crack is embedded entirely within the cuboid is calculated also. It is observed that the SIF values along the crack front are larger when the crack is closer to the surface of the cuboid than those when the crack is far away from the surface.

The FORTRAN program developed is limited to regular surface cracks which can be discretized by the quadrilateral shape function; it is not very efficient and suitable for irregular crack shapes.

The evaluation of the 3D mixed‐mode SIFs in the transversely isotropic material may have direct practical applications. The SIFs have been used in engineering design to obtain the safety factor of the elastic structures.

This is the first time that the special nine‐node quadrilateral shape function has been applied to the boundary containing the crack mouth. The numerical method developed can be applied to the SIF calculation in a finite transversely isotropic cuboid within an inclined surface crack. The computational approach and the results of SIFs are of great value for the modeling and design of anisotropic elastic structures.

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.

3D steady heat conduction analysis considering heat source is conducted on the fundamental of the fast multipole method (FMM)-accelerated line integration boundary element method (LIBEM).

Due to considering the heat source, domain integral is generated in the traditional heat conduction boundary integral equation (BIE), which will counteract the well-known merit of the BEM, namely, boundary-only discretization. To avoid volume discretization, the enhanced BEM, the LIBEM with dimension reduction property is introduced to transfer the domain integral into line integrals. Besides, owing to the unsatisfactory performance of the LIBEM when it comes to large-scale structures requiring massive computation, the FMM-accelerated LIBEM (FM-LIBEM) is proposed to improve the computation efficiency further.

Assuming N and M are the numbers of nodes and integral lines, respectively, the FM-LIBEM can reduce the time complexity from O(NM) to about O(N+ M), and a full discussion and verification of the advantage are done based on numerical examples under heat conduction.

(1) The LIBEM is applied to 3D heat conduction analysis with heat source. (2) The domain integrals can be transformed into boundary integrals with straight line integrals by the LIM. (3) A FM-LIBEM is proposed and can reduce the time complexity from O(NM) to O(N+ M). (4) The FM-LIBEM with high computational efficiency is exerted to solve 3D heat conduction analysis with heat source in massive computation successfully.

At present numerical simulation of seismicity, used in mining and hydraulic fracturing practice, is quite time expensive what hampers its combined employing with observed seismicity in real time. The purpose of this paper is to suggest a mean for drastic speeding up numerical modeling seismic and aseismic events.

The authors propose the means to radically decrease the time expense for the bottleneck stage of simulation: calculations of stresses, induced by a large group of already activated flaws (sources of events), at locations of flaws of another large group, which may be activated by the stresses. This is achieved by building a hierarchical tree and properly accounting for the sizes of activated flaws, excluding check of their influence on flaws, which are beyond strictly defined near-regions of strong interaction.

Comparative simulations of seismicity by conventional and improved methods demonstrate high efficiency of the means developed. When applied to practical mining and hydrofracturing problems, it requires some two orders less time to obtain practically the same output results as those of conventional methods.

The proposed improvement provides a means for simulation of seismicity in real time of mining steps and hydrofracture propagation. It can be also used in other applications involving seismic and aseismic events and acoustic emission.

This paper aims to develop a new isogeometric boundary element formulation based on non-uniform rational basis splines (NURBS) curves for solving Reissner’s shear-deformable plates.

The generalized displacements and tractions along the problem boundary are approximated as NURBS curves having the same rational B-spline basis functions used to describe the geometrical boundary of the problem. The source points positions are determined over the problem boundary by the well-known Greville abscissae definition. The singular integrals are accurately evaluated using the singularity subtraction technique.

Numerical examples are solved to demonstrate the validity and the accuracy of the developed formulation.

This formulation is considered to preserve the exact geometry of the problem and to reduce or cancel mesh generation time by using NURBS curves employed in computer aided designs as a tool for isogeometric analysis. The present formulation extends such curves to be implemented as a stress analysis tool.