Importance and use of rigid body mode in boundary element method

Abstract

SUMMARYThe superposition of a rigid body mode on a body should result in a corresponding change in displacement values but should not affect the stresses. However, in the numerical solution by the boundary element method (BEM) large errors may be obtained for displacements and stresses if a rigid body mode is present in the input data. To eliminate the effects of the rigid body mode on the numerical accuracy of the solution, the fundamental solutions for displacements must be correctly interpreted and used. The rigid body mode may be unknowingly present in the boundary condition data. It may be present because the boundary data are not known accurately. Or it may be present if the displacement values at the Support have been computed from a separate analysis. A rigid body mode may arise due to the collocation nature of satisfying the boundary conditions. The point values of the applied load at the collocation point may not satisfy equilibrium. Or the point values of the specified displacements may not satisfy the condition of zero translation and rotation. For bodies under pure traction, we know that the analytical solution can contain an arbitrary amount of rigid body mode. Numerically, however, some unknown value is assigned to this rigid body mode. It might be desirable (for example in limit analysis) to eliminate the rigid body mode from the displacements to obtain deformation of a point with respect to a point on the body. In addition, knowledge and elimination of the rigid body mode is necessary for the implementation of a scheme described by this author in an earlier work. The importance of the earlier work is that it reduced the sensitivity of the BEM to changes and errors in the input data. In this paper the causes, and the effects of the rigid body mode on the BEM, the correct interpretation of the fundamental solution for displacements and an algorithm for determining and accounting for the rigid body mode are discussed. A numerical example validates the ideas in this paper for the indirect version. The algorithm for the direct version is presented without a numerical example in the Appendices.

A necessary and sufficient boundary integral formulation for plane elasticity problems

“…The necessary and sufficient condition for the over-determined problem (l), (22) to have solutions is…”

Coupling of BEM with a large displacement and rotation algorithm

“…Starting from the Fredholm theory of linear operators with zero index 18, 19, different mathematical approaches, referred to as Methods F by Blazquez et al 15, have been proposed by various authors 9, 15, 18, 20–28 with the aim of removing RBMs in elastostatic BEM problems. Methods F can be subdivided into two categories.…”

“…Methods F can be subdivided into two categories. (a) Method F1, also called the augmenting method or the bordering method, has been considered by various researchers9, 15, 18, 20, 21; (b) Method F2, sometimes called the completion method, has been considered by various authors 15, 23–28.…”

“…In engineering practice, additional point supports are directly placed on the unrestrained body when solving TBVPs with the aim of removing non‐uniqueness 31. However, the approach gives inaccurate results for some examples 15, 28, 32.…”

Solution of viscous flow problems by using the boundary element method

“…Studies carried out to improve the methodology and the tools of the BEM include the work of Vable [33] which concerns the rigid body term of the integral equation and determination of the appropriateness of given boundary conditions. Other areas of study include the evaluation of boundary and domain integrals.…”