Nicht-glatte Mehrkörpersysteme werden entscheidend von den Zuständen der einseitigen Kontakte beeinflusst. Die Unstetigkeiten in den Zustandsgrößen oder deren Ableitungen sowie die Strukturvarianz führen zu Problemen bei der numerischen Integration der Bewegungsgleichungen.