On some unconditionally stable, higher order methods for the numerical solution of the structural dynamics equations