This paper presents a new computational formulation for large strain polyconvex elasticity. The formulation is based on using the deformation gradient (the fibre map), its adjoint (the area map) and its determinant (the volume map) as independent kinematic variables of a convex strain energy function. Compatibility relationships between these variables and the deformed geometry are enforced by means of a multi-field variational principle with additional constraints. This process allows the use of different approximation spaces for each variable. The paper also introduces conjugate stresses to these kinematic variables which can be used to define a generalised convex complementary energy function and a corresponding complementary energy principle of the Hellinger-Reissner type, where the new conjugate stresses are primary variables together with the deformed geometry. Both compressible and incompressible or nearly incompressible elastic models are considered. A key element to the developments presented in the paper is the definition of a new tensor cross product which facilitates the algebra associated with the adjoint of the deformation gradient. Two particular choices of interpolation spaces are considered to illustrate the formulation. The first is based on quadratic displacements, constant pressure discontinuous across element edges and linear discontinuous element stresses conjugate to the deformation gradient and its adjoint. The second choice involves linear displacement and continuous linear pressure stabilised using a Petrov-Galerkin technique.

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 [1], 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 use of a new tensor cross product 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-1 Corresponding author: j.bonet@swansea.ac.uk 2 Corresponding author: a.j.gil@swansea.ac.uk 3 Corresponding author: c.h.lee@swansea.ac.uk Preprint submitted to Computer Methods in Applied Mechanics and EngineeringSeptember 16, 2014Rivlin 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.

No abstract

Please cite this article in press as: M. Aguirre et al., A vertex centred Finite Volume Jameson-Schmidt-Turkel (JST) algorithm for a mixed conservation formulation in solid dynamics, Journal of Computational Physics (2013), http://dx.doi.org/10. 1016/j.jcp.2013.12.012 This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain. A mixed formulation based upon the use of the linear momentum, the deformation gradient tensor and the total energy as conservation variables is discretised in space using linear triangles and tetrahedra in two-dimensional and three-dimensional computations, respectively. The scheme is implemented using central differences for the evaluation of the interface fluxes in conjunction with the Jameson-Schmidt-Turkel (JST) artificial dissipation term. The discretisation in time is performed by using a Total Variational Diminishing (TVD) two-stage Runge-Kutta time integrator. The JST algorithm is adapted in order to ensure the preservation of linear and angular momenta. The framework results in a low order computationally efficient solver for solid dynamics, which proves to be very competitive in nearly incompressible scenarios and bending dominated applications.

A mixed second order stabilised Petrov-Galerkin finite element framework was recently introduced by the authors (C.H.Lee, A.J.Gil and J.Bonet. "Development of a stabilised Petrov-Galerkin formulation for conservation laws in Lagrangian fast solid dynamics", CMAME, 268:40-64, 2014). The new mixed formulation, written as a system of conservation laws for the linear momentum and the deformation gradient, performs extremely well in bending dominated scenarios (even when linear tetrahedral elements are used) yielding equal order of convergence for displacements and stresses. In this paper, this formulation is further enhanced for nearly and truly incompressible deformations with three key novelties. First, a new conservation law for the Jacobian of the deformation is added into the system providing extra flexibility to the scheme. Second, a variationally consistent PetrovGalerkin stabilisation methodology is derived. Third, an adapted fractional step method is presented for both incompressible and nearly incompressible materials in the context of nonlinear elastodynamics. For completeness and ease of understanding, these three improvements are presented both in small and large strain regimes, studying the eigenstructure of the resulting systems. A series of numerical examples are presented in order to demonstrate the robustness of the enhanced methodology with respect to the work previously published by the authors.