A Penalty Function Method for Modelling Frictional Contact in MPM