A first order hyperbolic framework for large strain computational solid dynamics. Part II: Total Lagrangian compressible, nearly incompressible and truly incompressible elasticity

Abstract In Part I of this series, Bonet et al. (2015) 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. (2014), is presented to extend the 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.

[1]  Guglielmo Scovazzi,et al.  Stabilized shock hydrodynamics: II. Design and physical interpretation of the SUPG operator for Lagrangian computations☆ , 2007 .

[2]  A. J. Gil,et al.  A first‐order hyperbolic framework for large strain computational solid dynamics: An upwind cell centred Total Lagrangian scheme , 2017 .

[3]  O. C. Zienkiewicz,et al.  Triangles and tetrahedra in explicit dynamic codes for solids , 1998 .

[4]  Antonio J. Gil,et al.  A two-step Taylor-Galerkin formulation for fast dynamics , 2014 .

[5]  A. Huerta,et al.  Finite Element Methods for Flow Problems , 2003 .

[6]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[7]  P. Neff,et al.  Invariant formulation of hyperelastic transverse isotropy based on polyconvex free energy functions , 2003 .

[8]  R. D. Wood,et al.  Nonlinear Continuum Mechanics for Finite Element Analysis , 1997 .

[9]  Michael A. Puso,et al.  A highly efficient enhanced assumed strain physically stabilized hexahedral element , 2000 .

[10]  J. C. Simo,et al.  Variational and projection methods for the volume constraint in finite deformation elasto-plasticity , 1985 .

[11]  Rogelio Ortigosa,et al.  On a tensor cross product based formulation of large strain solid mechanics , 2016 .

[12]  Guglielmo Scovazzi,et al.  Lagrangian shock hydrodynamics on tetrahedral meshes: A stable and accurate variational multiscale approach , 2012, J. Comput. Phys..

[13]  Antonio J. Gil,et al.  A coupled hp-finite element scheme for the solution of two-dimensional electrostrictive materials , 2012 .

[14]  Antonio J. Gil,et al.  Development of a cell centred upwind finite volume algorithm for a new conservation law formulation in structural dynamics , 2013 .

[15]  J. Peraire,et al.  A variationally consistent mesh adaptation method for triangular elements in explicit Lagrangian dynamics , 2010 .

[16]  Mikhail Shashkov,et al.  A multi-scale Q1/P0 approach to langrangian shock hydrodynamics. , 2006 .

[17]  Guglielmo Scovazzi,et al.  A generalized view on Galilean invariance in stabilized compressible flow computations , 2010 .

[18]  T. Hughes,et al.  A new finite element formulation for computational fluid dynamics: V. Circumventing the Babuscka-Brezzi condition: A stable Petrov-Galerkin formulation of , 1986 .

[19]  Peter Wriggers,et al.  A new Mixed Finite Element based on Different Approximations of the Minors of Deformation Tensors , 2011 .

[20]  Clark R. Dohrmann,et al.  A uniform nodal strain tetrahedron with isochoric stabilization , 2009 .

[21]  T. Hughes.  Multiscale phenomena: Green's functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origins of stabilized methods , 1995 .

[22]  Kenji Amaya,et al.  A locking‐free selective smoothed finite element method using tetrahedral and triangular elements with adaptive mesh rezoning for large deformation problems , 2014 .

[23]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[24]  Guglielmo Scovazzi,et al.  A discourse on Galilean invariance, SUPG stabilization, and the variational multiscale framework , 2007 .

[25]  T. Hughes.  Generalization of selective integration procedures to anisotropic and nonlinear media , 1980 .

[26]  Mikhail Shashkov,et al.  Multi-Scale Lagrangian Shock Hydrodynamics on Q1/P0 Finite Elements: Theoretical Framework and Two-dimensional Computations. , 2008 .

[27]  Antonio J. Gil,et al.  An upwind vertex centred Finite Volume solver for Lagrangian solid dynamics , 2015, J. Comput. Phys..

[28]  William J. Rider,et al.  A conservative nodal variational multiscale method for Lagrangian shock hydrodynamics , 2010 .

[29]  Guglielmo Scovazzi,et al.  Galilean invariance and stabilized methods for compressible flows , 2007 .

[30]  Peter Wriggers,et al.  Approximation of incompressible large deformation elastic problems: some unresolved issues , 2013 .

[31]  D. Benson.  Computational methods in Lagrangian and Eulerian hydrocodes , 1992 .

[32]  B. Carnes,et al.  A simple, stable, and accurate linear tetrahedral finite element for transient, nearly, and fully incompressible solid dynamics: a dynamic variational multiscale approach , 2016 .

[33]  Oubay Hassan,et al.  An averaged nodal deformation gradient linear tetrahedral element for large strain explicit dynamic applications , 2001 .

[34]  T. Hughes,et al.  B¯ and F¯ projection methods for nearly incompressible linear and non-linear elasticity and plasticity using higher-order NURBS elements , 2008 .

[35]  Jintai Chung,et al.  A Time Integration Algorithm for Structural Dynamics With Improved Numerical Dissipation: The Generalized-α Method , 1993 .

[36]  J. Ball.  Convexity conditions and existence theorems in nonlinear elasticity , 1976 .

[37]  Rogelio Ortigosa,et al.  A first order hyperbolic framework for large strain computational solid dynamics. Part I: Total Lagrangian isothermal elasticity , 2015 .

[38]  D. Owen,et al.  Design of simple low order finite elements for large strain analysis of nearly incompressible solids , 1996 .

[39]  Clark R. Dohrmann,et al.  Uniform Strain Elements for Three-Node Triangular and Four-Node Tetrahedral Meshes , 1999 .

[40]  T. Hughes,et al.  A new finite element formulation for computational fluid dynamics: I. Symmetric forms of the compressible Euler and Navier—Stokes equations and the second law of thermodynamics , 1986 .

[41]  Oubay Hassan,et al.  A discrete geometric conservation law (DGCL) for a cell vertex finite‐volume algorithm on moving domains , 2010 .

[42]  Antonio J. Gil,et al.  Development of a stabilised Petrov–Galerkin formulation for conservation laws in Lagrangian fast solid dynamics , 2014 .

[43]  Reint Boer,et al.  Vektor- und Tensorrechnung für Ingenieure , 1982 .

[44]  Antonio J. Gil,et al.  A stabilised Petrov-Galerkin formulation for linear tetrahedral elements in compressible, nearly incompressible and truly incompressible fast dynamics , 2014 .

[45]  P. Thomas,et al.  Geometric Conservation Law and Its Application to Flow Computations on Moving Grids , 1979 .

[46]  Michael A. Puso,et al.  A stabilized nodally integrated tetrahedral , 2006 .

[47]  Bruno Després,et al.  Discretization of hyperelasticity on unstructured mesh with a cell-centered Lagrangian scheme , 2010, J. Comput. Phys..

[48]  Antonio J. Gil,et al.  A vertex centred Finite Volume Jameson-Schmidt-Turkel (JST) algorithm for a mixed conservation formulation in solid dynamics , 2014, J. Comput. Phys..

[49]  Thomas J. R. Hughes,et al.  Stabilized shock hydrodynamics: I. A Lagrangian method , 2007 .

[50]  T. Belytschko,et al.  A uniform strain hexahedron and quadrilateral with orthogonal hourglass control , 1981 .

[51]  M. Chial,et al.  in simple , 2003 .

[52]  Rogelio Ortigosa,et al.  A computational framework for polyconvex large strain elasticity for geometrically exact beam theory , 2015, Computational Mechanics.

[53]  J. Bonet,et al.  A simple average nodal pressure tetrahedral element for incompressible and nearly incompressible dynamic explicit applications , 1998 .

[54]  A. Chorin.  Numerical Solution of the Navier-Stokes Equations* , 1989 .

[55]  Ngoc Cuong Nguyen,et al.  Hybridizable discontinuous Galerkin methods for partial differential equations in continuum mechanics , 2012, J. Comput. Phys..

[56]  Antonio J. Gil,et al.  An hp-fem framework for the simulation of electrostrictive and magnetostrictive materials , 2014 .

[57]  T. Hughes,et al.  A new finite element formulation for computational fluid dynamics: II. Beyond SUPG , 1986 .

[58]  J. Bonet,et al.  Stability and comparison of different linear tetrahedral formulations for nearly incompressible explicit dynamic applications , 2001 .