求质量矩阵集等

mengjiehuang

哪位大虾有质量矩阵集成以及单元质量矩阵的构成的小程序，有劳传一下
谢谢

[ 本帖最后由 mjhzhjg 于 2006-12-29 12:57 编辑 ]

AaronSpark

去计算力学版问吧

飞天一笑

现在这种单纯小功能的程序很少了,实在需要,你可以用ABQUS或ANSYS计算一下!

xinyuxf

1.你说的是什么单元的？
2.你需要的是一致质量矩阵还是集中质量矩阵？推荐集中质量矩阵
3.装配总体质量矩阵的程序都差不多的，下面是4节点矩形体单元的总体刚度矩阵装配方法
%%%%%%%%%装配单元质量矩阵%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function SM = RectAMM(M,EM,lmd)   %M整体质量矩阵、EM单元质量矩阵、lmd每个单元的节点坐标矩阵
for i = 1:1:4
    for j = 1:1:4
        for p = 1:1:2
            for q = 1:1:2
                r = (i-1) * 2+ p;
                s = (j-1) * 2 + q;
                R = (lmd(i)-1) *2+ p;
                S = (lmd(j)-1) * 2 + q;
                M(R,S) = M(R,S) + EM(r,s);
            end
        end
    end
end
SM=M;
return;

如果是其他单元，则只需将程序中的4和2替换掉即可（4表示单元有4个节点，2表示每个节点有2个自由度）

[ 本帖最后由 xinyuxf 于 2006-9-16 20:07 编辑 ]