A boundary weak singularity elimination method for multilayer structures

In the numerical discretization for the models of multilayer thin structures, the elements with narrow strip shape may be appeared. When using the previous methods to integrate these elements, the potentially near singularity will arise in the circumferential direction. In this paper, a weak singularity elimination method for multilayer structures of 3D boundary element method is presented to accurately calculate weakly singular integrals. Its numerical implementation is divided into the following two parts: (1) (alpha, beta) coordinate transformation is applied to deal with the weakly singular integral arising in boundary integral equation and separate out the integral in the circumferential direction (beta direction); (2) Based on the integral form of beta, a corresponding sinh transformation method is developed to eliminate the near singularity in beta direction. Finally, the proposed transformation formulation is extended to mull-domain potential problem and applied to the boundary analysis for multilayer structures. Several numerical examples are given to verify the computational accuracy of the proposed method. And the numerical results are promising and encouraging.

Automated summary

This paper proposes a method for the numerical discretization of multilayer thin structures to eliminate weak singularities, accurately calculate integrals through coordinate transformation and sinh transformation method. The proposed method shows promising numerical accuracy in several examples.