plotmode.m

有限元平面问题的求解相关程序。 包括质量矩阵

%---------------mode------------
V=zeros(2*nodnum,2*freenod);
V(2*(line+1)+1:2*nodnum,1:2*freenod)=40*v(:,:)/Tmax;  %将模态矩阵数值归一化
Ve=zeros(2*subnodnum,2*subfreenod);
Ve(2*(subline+1)+1:2*subnodnum,1:2*subfreenod)=10*vG(:,:)/tmax;    %将模态矩阵数值归一化
number=10;       % input('please input a number:');
modeaixs=zeros(nodnum,2*number);   %各节点坐标矩阵
for j=1:number
    for i=1:nodnum
        modeaixs(i,2*j-1)=V(2*i-1,j)+aixs(i,1);
        modeaixs(i,2*j)=V(2*i,j)+aixs(i,2);
    end
end
modeaixsmx=zeros(line+1,ar+1);
modeaixsmy=zeros(line+1,ar+1);
for numb=1:number
    figure;
    t=1;
    for j=1:ar+1
        for i=1:line+1
            modeaixsmx(i,j)=modeaixs(t,2*numb-1);
            modeaixsmy(i,j)=modeaixs(t,2*numb);
            t=t+1;
        end
    end
    mesh(modeaixsmx,modeaixsmy,aixsmz);
    colormap(white);
    title(['No. ',num2str(numb),' mode']);
    view(2);
    axis off;
end
submodeaixs=zeros(subnodnum,2*number);   %各节点坐标矩阵
for j=1:number
    for i=1:subnodnum
        modeaixs(i,2*j-1)=Ve(2*i-1,j)+subaixs(i,1);
        modeaixs(i,2*j)=Ve(2*i,j)+subaixs(i,2);
    end
end
modeaixsmx=zeros(subline+1,subar+1);
modeaixsmy=zeros(subline+1,subar+1);
for numb=1:number
    figure;
    t=1;
    for j=1:subar+1
        for i=1:subline+1
            modeaixsmx(i,j)=modeaixs(se1(t,1),2*numb-1);
            modeaixsmy(i,j)=modeaixs(se1(t,1),2*numb);
            t=t+1;
        end
    end
    mesh(modeaixsmx,modeaixsmy,subaixsmzg);
    colormap(white);
    title(['No. ',num2str(numb),' mode']);
    view(2);
    axis off;
end
%------------------------------------------------