for performance reasons instead of handling contacts globally, such as in continuum methods or solving exact Coulomb friction. Although penalty forces can cause instabilities, using a semi-implicit integration method combined with tapering spring constants during contact release produces stable results.