Simulation of Granular Media with the Material Point Method

We propose an extension to the Material Point Method for the simulation of granular media. We model the dynamics with an elastoplastic, continuum assumption. The behavior of the granular media is captured by the Drucker-Prager yield criterion that naturally represents the frictional relationship between shear and normal stresses. We develop a stress projection algorithm that is well suited for both explicit and implicit time integration, and uses a non-associative flow rule to ensure volume preservation. Our approach is able to recreate the dynamics of granular media undergoing large deformations, topological changes, and collisions.