If you include an arbitrary damping matrix in your model, the modes will no longer be orthogonal with respect to the eigenproblem.
Simply speaking, the “normal” modal solution doesn’t exist. (In addition to this your FE package would have to solve a polynomial eigenproblem, but that could be...