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

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.

No abstract

A first order hyperbolic framework for large strain computational solid dynamics. Part II: Total Lagrangian compressible, nearly incompressible and truly incompressible elasticity Antonio J. Gil, Chun Hean Lee, Javier Bonet, Rogelio Ortigosa PII: S0045-7825(15) 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 first order hyperbolic framework for large strain computational solid dynamics. AbstractIn Part I of this series, Bonet et al.[1] introduced a new computational framework for the analysis of large strain isothermal fast solid dynamics, where a mixed set of Total Lagrangian conservation laws was presented in terms of the linear momentum and an extended set of strain measures, namely the deformation gradient, its co-factor and its Jacobian. The main aim of this paper is to expand this formulation to the case of nearly incompressible and truly incompressible materials. The paper is further enhanced with three key novelties. First, the use of polyconvex nearly incompressible strain energy functionals enables the definition of generalised convex entropy functions and associated entropy fluxes. Two variants of the same formulation can then be obtained, namely, conservation-based and entropy-based, depending on the unknowns of the system. Crucially, the study of the eigenvalue structure of the system is carried out in order to demonstrate its hyperbolicity and, thus, obtain the correct time step bounds for explicit time integrators. Second, 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. Third, an adapted fractional step method, built upon the work presented in Gil et al. range of applications towards the incompressibility limit. Finally, a series of numerical examples are presented in order to assess the applicability and robustness of the proposed formulation. The overall scheme shows excellent behaviour in compressible, nearly incompressible and truly incompressible scenarios, yielding equal order of convergence for velocities and stresses.

A vertex centred Jameson-Schmidt-Turkel (JST) finite volume algorithm was recently introduced by the authors (Aguirre et al., 2014 [1]) in the context of fast solid isothermal dynamics. The spatial discretisation scheme was constructed upon a Lagrangian two-field mixed (linear momentum and the deformation gradient) formulation presented as a system of conservation laws [2][3][4]. In this paper, the formulation is further enhanced by introducing a novel upwind vertex centred finite volume algorithm with three key novelties. First, a conservation law for the volume map is incorporated into the existing two-field system to extend the range of applications towards the incompressibility limit (Gil et al., 2014 [5]). Second, the use of a linearised Riemann solver and reconstruction limiters is derived for the stabilisation of the scheme together with an efficient edge-based implementation. Third, the treatment of thermo-mechanical processes through a Mie-Grüneisen equation of state is incorporated in the proposed formulation. For completeness, the study of the eigenvalue structure of the resulting system of conservation laws is carried out to demonstrate hyperbolicity and obtain the correct time step bounds for non-isothermal processes. A series of numerical examples are presented in order to assess the robustness of the proposed methodology. The overall scheme shows excellent behaviour in shock and bending dominated nearly incompressible scenarios without spurious pressure oscillations, yielding second order of convergence for both velocities and stresses.

This paper presents a new Smooth Particle Hydrodynamics (SPH) computational framework for large strain explicit solid dynamics. A mixed-based set of Total Lagrangian conservation laws [1,2] is presented in terms of the linear momentum and an extended set of geometric strain measures, comprised of the deformation gradient, its co-factor and the Jacobian. Taking advantage of this representation, the main aim of this paper is the adaptation of the very efficient Jameson-Schmidt-Turkel (JST) algorithm [3], extensively used in computational fluid dynamics, to a SPH based discretisation of the mixed-based set of conservation laws, with three key distinct novelties. First, a conservative JST-based SPH computational framework is presented with emphasis in nearly incompressible materials. Second, the suppression of numerical instabilities associated with the non-physical zero-energy modes is addressed through a well-established stabilisation procedure. Third, the use of a discrete angular momentum projection algorithm is presented in conjunction with a monolithic Total Variation Diminishing Runge-Kutta time integrator in order to guarantee the global conservation of angular momentum. For completeness, exact enforcement of essential boundary conditions is incorporated through the use of a Lagrange multiplier projection technique. A series of challenging numerical examples (e.g. in the near incompressibility regime) are examined in order to assess the robustness and accuracy of the proposed algorithm. The obtained results are benchmarked against a wide spectrum of alternative numerical strategies.

Citation for this version held on GALA:Haider, Jibran, Lee, Chun Hean, Gil, Antonio J. and Bonet, Javier (2016) A first order hyperbolic framework for large strain computational solid dynamics: An upwind cell centred Total Lagrangian scheme. London: Greenwich Academic Literature Archive. , a first order hyperbolic system of conservation laws was introduced in terms of the linear momentum and the deformation gradient tensor of the system, leading to excellent behaviour in two dimensional bending dominated nearly incompressible scenarios. The main aim of this paper is the extension of this algorithm into three dimensions, its tailor-made implementation into OpenFOAM and the enhancement of the formulation with three key novelties. First, the introduction of two different strategies in order to ensure the satisfaction of the underlying involutions of the system, that is, that the deformation gradient tensor must be curl-free throughout the deformation process. Second, the use of a discrete angular momentum projection algorithm and a monolithic Total Variation Diminishing Runge-Kutta time integrator combined in order to guarantee the conservation of angular momentum. Third, and for comparison purposes, an adapted Total Lagrangian version of the Hyperelastic-GLACE nodal scheme of Kluth and Després [2] is presented. A series of challenging numerical examples are examined in order to assess the robustness and accuracy of the proposed algorithm, benchmarking it against an ample spectrum of alternative numerical strategies developed by the authors in recent publications.