Model Reduction for the Material Point Method via Learning the Deformation map and its Spatial-temporal Gradients

Peter Yichen Chen, Columbia University

Maurizio Chiaramonte, Facebook Reality Labs Research

Eitan Grinspun, University of Toronto, Columbia University

Kevin Carlberg, Facebook Reality Labs Research

The reduced-order elasticity simulation contains 3 million particles, runs in real-time, and is an order-of-magnitude faster than its full-order MPM counterpart.

Abstract

This work proposes a model-reduction approach for the material point method on nonlinear manifolds. The technique approximates the kinematics by approximating the deformation map in a manner that restricts deformation trajectories to reside on a low-dimensional manifold expressed from the extrinsic view via a parameterization function. By explicitly approximating the deformation map and its spatial-temporal gradients, the deformation gradient and the velocity can be computed simply by differentiating the associated parameterization function. Unlike classical model reduction techniques that build a subspace for a finite number of degrees of freedom, the proposed method approximates the entire deformation map with infinite degrees of freedom. Therefore, the technique supports resolution changes in the reduced simulation, attaining the challenging task of zero-shot super-resolution by generating material points unseen in the training data. The ability to generate material points also allows for adaptive quadrature rules for stress update. A family of projection methods is devised to generate dynamics, i.e., at every time step, the methods perform three steps: (1) generate quadratures in the full space from the reduced space, (2) compute position and velocity updates in the full space, and (3) perform a least-squares projection of the updated position and velocity onto the low-dimensional manifold and its tangent space. Computational speedup is achieved via hyper-reduction, i.e., only a subset of the original material points are needed for dynamics update. Large-scale numerical examples with millions of material points illustrate the method's ability to gain an order-of-magnitude computational-cost saving -- indeed real-time simulations in some cases -- with negligible errors.