有阻尼时结构的频率和振型与无阻尼时差很远？

secondye

请教了，这方面不是很懂，能不能解释具体一点？
无阻尼时[U,S]=eig(M,k);
有阻尼时怎么就变成[U,S]=eig(inv(M)*K)了？这样解出来的频率跟上面差很多哦，好像不对啊

在这个帖子上http://forum.vibunion.com/thread-43483-1-1.html
20052731 兄说，转换到状态空间，系统矩阵的零点就是特征值。
本人最近在读一程序，算法似乎跟20052731 兄说的一样，但还不是很理解。
程序是这样的：对一桁架模型，先用静力凝聚法缩减刚度矩阵，下面的看不懂了，好像是求系统的A，B，C，D矩阵，这部分贴出来:

Mhat=M_In(Inu,Inu); % 经缩减后的质量矩阵
Khat=Kuu-Kus*AA;    % 经缩减后的刚度矩阵
invM=Mhatdiag(ones(length(Mhat),1));

[UV2,w2]=eig(invM*Khat);  
omega=sqrt(diag(w2));
ff=omega/2/pi;
zeta=0.005;              
Cq=diag(2*zeta*omega);        %%% 是否阻尼矩阵对角线值？
C=Mhat*UV2*Cq*inv(UV2);    %%% 是否构造阻尼矩阵？
Ndof=length(Inu);
%return
Aa=[zeros(Ndof) eye(Ndof);
   -invM*Khat     -invM*C  ];
Ba=[zeros(Ndof,length(Phat(1,:))); invM*Phat];
% % outputs are accelerations.
% CC1=[-invM*Khat -invM*C];
%outputs are displacements.
CC1 = [eye(Ndof) zeros(Ndof)];
Ca=CChat*CC1;

本人是新手啊，还请各位多指教

VibrationMaster

1. 阻尼比较小的时候，阻尼对振型和固有频率影响比较小，求后者可用无阻尼的代替
2. 通常计算中阻尼设置成比例阻尼，对这种特殊的阻尼形式，有无阻尼两种情形的振型和固有频率精确相等