zbMATH — the first resource for mathematics

A local hypersingular boundary integral equation method using a triangular background mesh. (English) Zbl 1232.74120
Summary: A new meshless local hypersingular boundary integral equation method is presented for the analysis of two-dimensional elastostatic problems. The elastic domain is discretized by placing arbitrarily nodes on its boundary and interior. Given this set of nodes, the corresponding map of background triangles is constructed through a common triangulation algorithm. The local domain of each node consists of the union of triangles that this point lies, thus, creating a polygonal line of its local boundary. The local boundary integral equations of both displacements and stresses of the conventional boundary elements method are taken into account. The interpolation of the unknown fields is performed by taking each face of a triangle of the local domain of a source point as a one-dimensional line element. The essential boundary conditions can be directly implemented easily because the interpolation functions possess the Kronecker delta-function property. After constructing the final linear system of equations, the only unknowns are displacements and stresses of all nodal points. Thus, leading to a banded stiffness matrix as in the Finite Element Method. The effectiveness and efficiency of the proposed method is demonstrated with three elastostatic problems in two-dimensions. Excellent agreement between the numerical results and the exact solutions is found. The numerical examples also show that the accuracy of the proposed method is as good as that of the Boundary Elements Method.

74S15 Boundary element methods applied to problems in solid mechanics
65N38 Boundary element methods for boundary value problems involving PDEs
65R20 Numerical methods for integral equations
45F99 Systems of linear integral equations
Full Text: DOI