Implicit dynamics in the material-point method

Abstract A time-implicit discretization is derived and validated for the material-point method (MPM). The resulting non-linear, discrete equations are solved using Newton's method combined with either the conjugate gradient method or the generalized minimum residual method. These Newton–Krylov solvers are implemented in a matrix-free fashion for numerical efficiency. A description of the algorithms and evaluation of their performance is presented. On all test problems, if the time step is chosen appropriately, the implicit solution technique is more efficient than an explicit method without loss of desired features in the solutions. In a dramatic example, time steps 10,000 times the explicit step size are possible for the large deformation compression of a cylindrical billet at 1.2% the computational cost.

[1]  C. Kelley.  Iterative Methods for Linear and Nonlinear Equations , 1987 .

[2]  D. Sulsky,et al.  THE MATERIAL POINT METHOD FOR SIMULATION OF THIN MEMBRANES , 1999 .

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

[4]  Howard L. Schreyer,et al.  Modeling delamination as a strong discontinuity with the material point method , 2002 .

[5]  Richard Barrett,et al.  Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods , 1994, Other Titles in Applied Mathematics.

[6]  Jose Luis Silveira,et al.  An algorithm for shakedown analysis with nonlinear yield functions , 2002 .

[7]  Giovanni Lapenta,et al.  Control of the number of particles in fluid and MHD particle in cell methods , 1995 .

[8]  Howard L. Schreyer,et al.  Fluid–membrane interaction based on the material point method , 2000 .

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

[10]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[11]  Howard L. Schreyer,et al.  Axisymmetric form of the material point method with applications to upsetting and Taylor impact problems , 1996 .

[12]  Lee A. Segel,et al.  Mathematics applied to deterministic problems in the natural sciences , 1974, Classics in applied mathematics.

[13]  G. Golub,et al.  Gmres: a Generalized Minimum Residual Algorithm for Solving , 2022 .

[14]  Ivo Babuška,et al.  Approximation properties of the h-p version of the finite element method , 1996 .

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

[16]  Dana A. Knoll,et al.  Enhanced Nonlinear Iterative Techniques Applied to a Nonequilibrium Plasma Flow , 1998, SIAM J. Sci. Comput..

[17]  J. Brackbill,et al.  Flip: A low-dissipation, particle-in-cell method for fluid flow , 1988 .

[18]  J. Brackbill,et al.  An implicit particle-in-cell method for granular materials , 2002 .

[19]  O. C. Zienkiewicz,et al.  A new cloud-based hp finite element method , 1998 .

[20]  James E. Guilkey,et al.  An implicit time integration strategy for use with the material point method , 2001 .

[21]  D. Sulsky.  Erratum: Application of a particle-in-cell method to solid mechanics , 1995 .

[22]  J. Brackbill,et al.  The material-point method for granular materials , 2000 .

[23]  Dana A. Knoll,et al.  A Multilevel Iterative Field Solver for Implicit, Kinetic, Plasma Simulation , 1999 .