A Discontinuous Galerkin Material Point Method for the solution of impact problems in solid dynamics

Abstract An extension of the Material Point Method [1] based on the Discontinuous Galerkin approximation (DG) [2] is presented here. A solid domain is represented by a collection of particles that can move and carry the fields of the problem inside an arbitrary computational grid in order to provide a Lagrangian description of the deformation without mesh tangling issues. The background mesh is then used as a support for the Discontinuous Galerkin approximation that leads to a weak form of conservation laws involving numerical fluxes defined at element faces. Those terms allow the introduction of the characteristic structure of hyperbolic problems within the numerical method by using an approximate Riemann solver [3] . The Discontinuous Galerkin Material Point Method, which can be viewed as a Discontinuous Galerkin Finite Element Method (DGFEM) with modified quadrature rule, aims at meeting advantages of both mesh-free and DG methods. The method is derived within the finite deformation framework for multidimensional problems by using a total Lagrangian formulation. A particular attention is paid to one specific discretization leading to a stability condition that allows to set the CFL number at one. The approach is illustrated and compared to existing or developed analytical solutions on one-dimensional problems and compared to the finite element method on two-dimensional simulations.

[1]  S. Osher.  Riemann Solvers, the Entropy Condition, and Difference , 1984 .

[2]  Robert B. Haber,et al.  A space-time discontinuous Galerkin method for linearized elastodynamics with element-wise momentum balance , 2006 .

[3]  John A. Trangenstein,et al.  Numerical Solution of Hyperbolic Partial Differential Equations , 2009 .

[4]  A. Tekkaya,et al.  Electromagnetic forming—A review , 2011 .

[5]  D. Sulsky,et al.  A particle method for history-dependent materials , 1993 .

[6]  Gerard R. Richter,et al.  An explicit finite element method for the wave equation , 1994 .

[7]  David H. Sharp,et al.  A conservative Eulerian formulation of the equations for elastic flow , 1988 .

[8]  Matthias Möller,et al.  A High Order Material Point Method , 2017 .

[9]  G. Chavent,et al.  A finite-element method for the 1-D water flooding problem with gravity , 1982 .

[10]  Deborah Sulsky,et al.  Mass matrix formulation of the FLIP particle-in-cell method , 1992 .

[11]  A. Harten.  High Resolution Schemes for Hyperbolic Conservation Laws , 2017 .

[12]  W. H. Reed,et al.  Triangular mesh methods for the neutron transport equation , 1973 .

[13]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[14]  Gilles Kluth.  Analyse mathématique et numérique de systèmes hyperélastiques et introduction de la plasticité , 2008 .

[15]  Zhen Chen,et al.  The Material Point Method: A Continuum-Based Particle Method for Extreme Loading Cases , 2016 .

[16]  P. Colella.  Multidimensional upwind methods for hyperbolic conservation laws , 1990 .

[17]  Thomas Heuzé,et al.  Lax-Wendroff and TVD finite volume methods for unidimensional thermomechanical numerical simulations of impacts on elastic-plastic solids , 2017, J. Comput. Phys..

[18]  Wing Kam Liu,et al.  Nonlinear Finite Elements for Continua and Structures , 2000 .

[19]  Bernardo Cockburn,et al.  Discontinuous Galerkin Methods for Convection-Dominated Problems , 1999 .

[20]  J. Brackbill,et al.  FLIP: A method for adaptively zoned, particle-in-cell calculations of fluid flows in two dimensions , 1986 .

[21]  Jie Liu,et al.  ParaView visualization of Abaqus output on the mechanical deformation of complex microstructures , 2017, Comput. Geosci..

[22]  Phillip Colella,et al.  A higher-order Godunov method for modeling finite deformation in elastic-plastic solids , 1991 .

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

[24]  Bruno Després,et al.  Perfect plasticity and hyperelastic models for isotropic materials , 2008 .

[25]  Nicolas Favrie,et al.  Dynamics of shock waves in elastic-plastic solids , 2011 .

[26]  Li-Li Wang,et al.  Foundations of Stress Waves , 2007 .

[27]  Anna Giacomini,et al.  Caveats on the implementation of the generalized material point method , 2008 .

[28]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[29]  Jan S. Hesthaven,et al.  High-order nodal discontinuous Galerkin particle-in-cell method on unstructured grids , 2006, J. Comput. Phys..

[30]  A. Ern,et al.  Mathematical Aspects of Discontinuous Galerkin Methods , 2011 .

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

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

[33]  Z. Więckowski.  The material point method in large strain engineering problems , 2004 .

[34]  Martin Berzins,et al.  Examination and Analysis of Implementation Choices within the Material Point Method (MPM) , 2008 .

[35]  J. Hesthaven,et al.  Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications , 2007 .

[36]  Deborah Sulsky,et al.  An unconditionally stable, energy–momentum consistent implementation of the material-point method , 2006 .

[37]  Utkarsh Ayachit,et al.  The ParaView Guide: A Parallel Visualization Application , 2015 .

[38]  Yulong Xing,et al.  Recent developments in discontinuous Galerkin finite element methods for partial differential equations : 2012 John H Barrett memorial lectures , 2014 .