This work develops robust contact algorithms capable of dealing with multibody nonsmooth contact
geometries for which neither normals nor gap functions can be defined. Such situations arise
in the early stage of fragmentation when a number of angular fragments undergo complex collision
sequences before eventually scattering. Such situations precludes the application of most contact
algorithms proposed to date