secondye 发表于 2010-10-28 11:44

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

请教了,这方面不是很懂,能不能解释具体一点?
无阻尼时=eig(M,k);
有阻尼时怎么就变成=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));

=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=;
% % outputs are accelerations.
% CC1=[-invM*Khat -invM*C];
%outputs are displacements.
CC1 = ;
Ca=CChat*CC1;

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

VibrationMaster 发表于 2010-10-28 12:50

1. 阻尼比较小的时候,阻尼对振型和固有频率影响比较小,求后者可用无阻尼的代替
2. 通常计算中阻尼设置成比例阻尼,对这种特殊的阻尼形式,有无阻尼两种情形的振型和固有频率精确相等
页: [1]
查看完整版本: 有阻尼时结构的频率和振型与无阻尼时差很远?