Mustaine,thanks!
what I want is the eigenvalue and vector of tangent stiffness matrix, so it should be solved like (K – λ I) ψ=0, where K is the tangent stiffness matrix, λ- eigenvalue, I- identity matrix, ψ- eigenvector...
But if I do frequency analysis for each configuration, it will...