A Penalty Function Method for Modelling Frictional Contact in MPM

In the present paper, the boundary of a continuum object is discretised separately from the MPM discretisation and traced accurately during the solution advancement. Therefore, contact forces among different boundaries are evaluated using a penalty function method commonly implemented in Lagrangian analyses. These forces are then mapped to the computational mesh, where the momentum equation is solved, as an external force. The suggested approach is validated with a benchmark problem, where the closed-form solution is available, and compared to the classical MPM contact algorithm. Furthermore, an application of two bodies collision with large deformation is demonstrated.