Latest revision as of 12:27, 1 July 2024

Abstract

In the classical theory of two-body contact, a single shared contact interface is considered between two continuum bodies, and is further discretized as such in the finite element setting. In general, however, the finite element mesh topology of two contacting bodies will be non-conforming at this shared interface, requiring the definition of a preferred or intermediate surface over which integral constraints may be evaluated. The specification of this interface is deemed to be somewhat arbitrary, but in practice the numerical solution of contact problems may exhibit sensitivity to the particular choice of intermediate surface. A further complication concerns the need to establish projective mappings between the discretized finite element surfaces and the chosen intermediate surface, particularly for the sake of evaluating the contact gap function between pairs of points on each of the two bodies. In this work, a new methodology for the enforcement of contact constraints in the context of finite element analyses is proposed. The method entails an alternative representation of contact surface integrals by equivalently integrating over the interstitial – albeit degenerate gap volume between two contacting bodies. An auxiliary indicator field is defined on each body, and is used to represent the degenerate interstitial volume as a non-degenerate hyper-dimensional gap volume. Over this domain, the gradient of the continuously interpolated displacement field with respect to the indicator field yields the oriented displacement gap, which may be used in the formulation of contact inequality constraints. Discretization of the hyper-dimensional gap volume into conforming finite elements is explored, and is observed to offer several advantages over existing contact discretization methods: the proposed method does not require the computation of geometric intersections or projections; it exploits conventional Gaussian quadrature schemes to integrate the hyper-dimensional gap integrals with a sufficient degree of accuracy; and may be naturally and efficiently extended to represent contact between higher-order surfaces. The efficacy of the method is demonstrated on several benchmark problems. Continuing and future work is also discussed, with a focus on intended applications and extensions of the method.

Full Paper