欧阳中华 发表于 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);我画出的振型图是错的 |
发表于 2010-8-25 17:01
提升卡
变色卡
显身卡
发表于 2010-8-30 14:03
楼主
