An acoustic Riemann solver for large strain computational contact dynamics

Abstract

AbstractThis article presents a vertex‐centered finite volume algorithm for the explicit dynamic analysis of large strain contact problems. The methodology exploits the use of a system of first order conservation equations written in terms of the linear momentum and a triplet of geometric deformation measures (comprising the deformation gradient tensor, its co‐factor, and its Jacobian) together with their associated jump conditions. The latter can be used to derive several dynamic contact models ensuring the preservation of hyperbolic characteristic structure across solution discontinuities at the contact interface, a clear advantage over the standard quasi‐static contact models where the influence of inertial effects at the contact interface is completely neglected. Taking advantage of the conservative nature of the formalism, both kinetic (traction) and kinematic (velocity) contact interface conditions are explicitly enforced at the fluxes through the use of appropriate jump conditions. Specifically, the kinetic condition is enforced in the usual linear momentum equation, whereas the kinematic condition can now be easily enforced in the geometric conservation equations without requiring a computationally demanding iterative algorithm. Additionally, a total variation diminishing shock capturing technique can be suitably incorporated in order to improve dramatically the performance of the algorithm at the vicinity of shocks. Moreover, and to guarantee stability from the spatial discretization standpoint, global entropy production is demonstrated through the satisfaction of semi‐discrete version of the classical Coleman–Noll procedure expressed in terms of the time rate of the so‐called Hamiltonian energy of the system. Finally, a series of numerical examples is examined in order to assess the performance and applicability of the algorithm suitably implemented in OpenFOAM. The knowledge of the potential contact loci between contact interfaces is assumed to be known a priori.