clc
clear;
zr=textread('Mode_Re.txt');
zi=textread('Mode_Im.txt');
r=abs(zr+ i*zi);
r=r/max(r);
e=angle(zr+ i*zi);
x=;
y=[(10:10:90); (10:10:90); (10:10:90)];
z=zeros(3,9);
frame_num=100;
for i=0:frame_num
z1=sin((i*2*pi/frame_num)+e).*r*0.7;
for i=1:3
for j=1:9
z(i,j)=z1(j+(i-1)*9);
end
end
surf(x,y,z)
axis()
pause(0.1)
end
dw04116 发表于 2012-3-14 09:30 static/image/common/back.gif
我计算梁的振型时,是把每个节点的横向位移提取出来当做纵坐标,横坐标按梁的节点位移坐标,然后就绘制出来 ...
你用的是一个节点两个自由度的梁单元吧,是求出特征向量以后,取出特征向量V某一列然后提取(v,theta)中V的所有数据,再画图吗 {:{03}:}挺好的帖子啊
页:
1
[2]