Invariant subspace reduction for linear dynamic analysis of finite-dimensional viscoelastic structures