Time efficient stiffness model computation for a parallel haptic mechanism via the virtual joint method