m_2d_beam_assembly.m

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

function [M] = M_2D_beam_assembly (in_data)

% assembly of mass matrix M for 2D-beam elements
dof_per_node = 3;

for i=1:size(in_data.MASS)        % check if masses exist
    t4 = in_data.MASS(i,1)*dof_per_node-(dof_per_node-1);
    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;
end;