Smooth interpolation of orientations with angular velocity constraints using quaternions