欧阳中华 发表于 2010-8-30 14:03
.
将特征向量中各个节点位移项找出来,加上边界点上的值就可以描出各阶的振型了. . . .
请问系统刚度矩阵为51阶的,节点矩阵为17阶的,我想画振型图,是不是在特征向量的某一阶中取出y节点的位移项(即2:2:51),因为我的单元刚度矩阵是(u1,v1,theta1,u2,v2,theta2),,
不要检查边界条件,消去相应的自由度吗,因为如果不消的话,我得到的特征向量就只有u的值
%求特征值和特征向量(圆频率和振型)
[V,D]=eig(kk2,mm2);
[lambda,ki]=sort(diag(D)); % 排序
omega=sqrt(lambda);
omega1=sqrt(lambda)/(2*pi);
V=V(:,ki);
%--------------------------------------------------------------------------
%--------------------------------------------------------------------------
%--------------------------------------------------------------------------
%画振型图
jk=3; %选择显示第3阶模态
Vi=V(:,jk); %振型向量的第jk列
x(1)=0
for i=1:16
x(i+1)=x(i)+gcoord(i,1);
end
i=1;
for j=2:3:45
y1(i)=Vi(j,1);
i=i+1;
end
y=[0,y1,0]
plot(x,y);我画出的振型图是错的 |