Here,we describe the robust and efficient application of the conventional 3D BEM in solving elasticity problems. We have focused on the precise computation of weakly singular integrals. The conformal Duffy-distance tr...Here,we describe the robust and efficient application of the conventional 3D BEM in solving elasticity problems. We have focused on the precise computation of weakly singular integrals. The conformal Duffy-distance transformation was employed to alleviate near singularities caused from two aspects:(1) the large aspect ratio of elements,i.e.,element shape distortions;and(2)the closeness of element boundaries to field points,i.e.,ill-shaped patches. Then,the rigid body motion method was employed to evaluate strongly singular integrals. Numerical solutions of 3D elastostatic problems demonstrated the high accuracy of the proposed method with coarse meshes and high convergence rates with mesh refinement. Compared with the Duffy transformation and original polar coordinate transformations,the proposed method is insensitive to element shapes.展开更多
基金This work was supported by the National Natural Science Foundation of China(Grant Nos.51879245,41731284&11672360)the Fundamental Research Funds for the Central Universities,China University of Geosciences(Wuhan)(Grant Nos.CUGCJ1821&CUG170645)。
文摘Here,we describe the robust and efficient application of the conventional 3D BEM in solving elasticity problems. We have focused on the precise computation of weakly singular integrals. The conformal Duffy-distance transformation was employed to alleviate near singularities caused from two aspects:(1) the large aspect ratio of elements,i.e.,element shape distortions;and(2)the closeness of element boundaries to field points,i.e.,ill-shaped patches. Then,the rigid body motion method was employed to evaluate strongly singular integrals. Numerical solutions of 3D elastostatic problems demonstrated the high accuracy of the proposed method with coarse meshes and high convergence rates with mesh refinement. Compared with the Duffy transformation and original polar coordinate transformations,the proposed method is insensitive to element shapes.
