Model reduction of contact dynamics simulation using a modified Lyapunov balancing method