如何求结构的自振频率和振型

liyangcdm

现在已经得到结构的刚度和质量矩阵，如何利用matlab来计算结构的自振频率和振型

[ 本帖最后由 eight 于 2007-11-6 14:50 编辑 ]

VibrationMaster

eig

appleseed05

[V,D]=eig(m\k)

sogooda

无阻尼的情况比较简单

1. [v,w]=eig(m\k);
   
2.   w=sqrt(w);

复制代码

有阻尼稍微郁闷点儿

1. 
   
2.   if norm(d/m*k-k/m*d) < 1e-8*norm(k/m*d)
   
3.     disp(比例阻尼 振型为实数.')
   
4.     [v,w]=eig(m\k);
   
5.     w=sqrt(w);
   
6.    else
   
7.     disp('非比例阻尼, 振型为复数.')
   
8.     a=[0*k eye(length(k));-m\k -m\d];
   
9.     [v,w1]=eig(a);
   
10.     w=abs(w1);
    
11.   end
    
12. 
    

复制代码

liyangcdm

为何会求出的振型和频率有复数形式啊
如何再做进一步处理呢

VibrationMaster

M，K应对称正定。 你随便给个数据未必能保证，就会出现复数形式

[ 本帖最后由 eight 于 2007-11-7 10:54 编辑 ]

liyangcdm

在计算结构的自振频率和振型时,如何在刚度矩阵中考虑位移约束条件