Gilles Daviet
We present a family of mixed Material Point Methods well suited to CFL-rate simulation of stiff elastoviscoplastic materials, up to the incompressible limit. Our work builds on the mixed discretization from Daviet and Bertails-Descoubes [2016a] and extends it to handle finite-strain viscoelasticity and more general flow rules, allowing the simulation of a much wider range of materials. Our implicit integration scheme yields a well-posed, symmetric optimization problem with compact stencils, together with an efficient GPU solver. We demonstrate our method on a variety of examples ranging from granular materials and snow to elastic solids, including two-way coupling with rigid-body solvers.