mode.m

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

number=11;       % input('please input a number:');
modeaixs=zeros(nodnum,2);   %各节点坐标矩阵
for j=number:number
    V=zeros(2*nodnum,1);
    V(2*(line+1)+1:2*nodnum,1)=4*v(:,number)/max(abs(v(:,number)));      %将模态矩阵数值归一化
    for i=1:nodnum
        modeaixs(i,2*j-1)=V(2*i-1,1)+aixs(i,1);
        modeaixs(i,2*j)=V(2*i,1)+aixs(i,2);
    end
end
modeaixsmx=zeros(line+1,ar+1);
modeaixsmy=zeros(line+1,ar+1);
figure;
for numb=number:number
    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=number:number
    Ve=zeros(2*subnodnum,1);
    Ve(2*(subline+1)+1:2*subnodnum,1)=-4*vG(:,number)/max(abs(vG(:,number)));      %将模态矩阵数值归一化
    for i=1:subnodnum
        modeaixs(i,2*j-1)=Ve(2*i-1,1)+subaixs(i,1);
        modeaixs(i,2*j)=Ve(2*i,1)+subaixs(i,2);
    end
end
modeaixsmx=zeros(subline+1,subar+1);
modeaixsmy=zeros(subline+1,subar+1);
figure;
for numb=number:number
    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
%------------------------------------------------