m_3d_beam_assembly.m

来自「结构力学中的有限元例子,包含了7个分类文件夹」· M 代码 · 共 26 行

function [M] = M_3D_beam_assembly (in_data)

% assembly of mass matrix M for 3D-beam elements

for i=1:size(in_data.MASS)        % check if masses exist
    t4 = in_data.MASS(i,1)*6-5;
    if in_data.MASS(i,2)>0
        M(t4,t4)  = in_data.MASS(i,2);
    end;
    if in_data.MASS(i,3)>0
        M(t4+1,t4+1) = in_data.MASS(i,3);
    end;
    if in_data.MASS(i,4)>0
        M(t4+2,t4+2) = in_data.MASS(i,4);   % mass matrix
    end;
    if in_data.MASS(i,5)>0
        M(t4+3,t4+3) = in_data.MASS(i,5);
    end;
    if in_data.MASS(i,6)>0
        M(t4+4,t4+4) = in_data.MASS(i,6);
    end;
    if in_data.MASS(i,7)>0
        M(t4+5,t4+5) = in_data.MASS(i,7);
    end;
end;