Multiparametric computation of eigenvalues for linear viscoelastic structures