This paper presents a new algorithm to identify the occurrence of impact (contact) among structures. Structures are modeled here by finite elements and the resulting formulation is employed to solve two-dimensional frame impact problems by the Lagrange Multiplier technique. The proposed algorithm is based on potential theory, largely employed in various fields of physics. The potential theory is used together with integral equations making possible to compare relative positions among nodes of different structures involved in a collision process. The present technique may be used to model the contact prevision among two or three-dimensional structures or between structures and rigid obstacles. In order to show the applicability of the proposed technique, numerical answers for 2D problems are compared with both analytical solutions and numerical results given by other authors.
IntroductionThe dynamic contact among structures is a subject of great interest in various engineering fields, such as mechanical, civil, aeronautic etc. Over the last two decades there has been a great advance in numerical modeling of impact, contact and crash problems, mainly using FEM. The first difficulty that emerges when leading with numerical modeling of impact is the proper identification, at a specific instant, of the contact occurrence. The main objective of this study is to present an alternative algorithm for this purpose.Two dimensional straight frame finite elements have been used to model structures in order to demonstrate the applicability of the technique. However, the algorithm can also be used with membrane elements, mixed problems and 3D simple applications. Virtual boundary elements (VBE) are used to integrate potential fields in order to identify the contact/impact occurrence. These VBE are placed over the boundary of the target structure. For frame or shell structures they are geometrically coincident with finite elements. Rigid and fixed bodies are also discretized here by straight virtual boundary elements and do not need finite element discretization.In the specialized literature there are algorithms that identify the contact by local procedures, as for example [1][2][3][4][5]. In this case the contact identification is performed for any point of the projectile structure considering circular (or other geometry) (2D) or spherical (3D) regions defined in the neighborhood of the target structure nodes, called here influence regions. This kind of prediction may present some problems when the geometry of the studied structures is very complicated. In this situation the referred methodologies can present superposition between influence regions leading to difficulties in the choice of the contact element. The proposed technique does not depend upon the local geometry of the colliding bodies, overcoming this kind of difficulty for contact identification.Other sophisticated and important algorithms can be found in literature; for example, references [6-10] are related to contact identification and momentum conservation and re...