Development of an implicit contact technique for the material point method

Abstract

An implicit contact algorithm for the material point method (MPM) has been developed to simulate contact. This allows recently developed implicit MPM codes to simulate large-scale deformations and interaction with external bodies. The performance of the method has been investigated and compared to an existing explicit method using benchmark and geotechnical examples. In particular, the proposed formulation has been shown to conserve energy in a similar way to the explicit formulation and reach similar results. The method typically converges to the analytical solution when an adequate time step and mesh size are used, with the time step generally being around ten times larger than the explicit method, although during the contact phase this does not always result in faster computation due to the iterative solution procedure.