In the most general case, the damped vibration modes of the system are the eigenvalues and vectors of the quadratic eigenproblem ##(m + \lambda d + \lambda^2k)x = 0##. The solutions will be of the form ##x(t) = X e^{(-\sigma + i \omega) t} ## where ##X## is a complex eigenvector, ##-\sigma + i \omega## is the complex eigenvalue, and the physical solution ##x(t)## is the real part of the right hand side of the equation. ##\sigma## represents the amount of damping, and ##\omega## the damped oscillation frequency (which is 0 for an overdamped system)
If you get lucky, the eigenvectors will be real, and the same as the eigenvectors of the undamped system ##(m - \omega_u^2 k)x = 0## (but of course the eigenvalues will be different, and the undamped frequency ##\omega_u## is not equal to the damped frequency ##\omega##).