求质量矩阵集等
哪位大虾有质量矩阵集成以及单元质量矩阵的构成的小程序,有劳传一下谢谢
[ 本帖最后由 mjhzhjg 于 2006-12-29 12:57 编辑 ] 去计算力学版问吧 现在这种单纯小功能的程序很少了,实在需要,你可以用ABQUS或ANSYS计算一下! 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 编辑 ]
页:
[1]