Randomized parallel simulation of constrained multibody systems for VR/haptic applications