A first order hyperbolic framework for large strain computational solid dynamics. Part I: Total Lagrangian isothermal elasticity

This paper introduces a new computational framework for the analysis of large strain fast solid dynamics. The paper builds upon previous work published by the authors (Gil etal., 2014) [48], where a first order system of hyperbolic equations is introduced for the simulation of isothermal elastic materials in terms of the linear momentum, the deformation gradient and its Jacobian as unknown variables. In this work, the formulation is further enhanced with four key novelties. First, the use of a new geometric conservation law for the co-factor of the deformation leads to an enhanced mixed formulation, advantageous in those scenarios where the co-factor plays a dominant role. Second, the use of polyconvex strain energy functionals enables the definition of generalised convex entropy functions and associated entropy fluxes for solid dynamics problems. Moreover, the introduction of suitable conjugate entropy variables enables the derivation of a symmetric system of hyperbolic equations, dual of that expressed in terms of conservation variables. Third, the new use of a tensor cross product [61] greatly facilitates the algebraic manipulations of expressions involving the co-factor of the deformation. Fourth, the development of a stabilised Petrov-Galerkin framework is presented for both systems of hyperbolic equations, that is, when expressed in terms of either conservation or entropy variables. As an example, a polyconvex Mooney-Rivlin material is used and, for completeness, the eigen-structure of the resulting system of equations is studied to guarantee the existence of real wave speeds. Finally, a series of numerical examples is presented in order to assess the robustness and accuracy of the new mixed methodology, benchmarking it against an ample spectrum of alternative numerical strategies, including implicit multi-field Fraeijs de Veubeke-Hu-Washizu variational type approaches and explicit cell and vertex centred Finite Volume schemes.