如何求单元的一致质量矩阵

ZFLY

悬臂梁的一致刚度矩阵和一致质量矩阵怎么算啊？？有什么书上有参考没有？
最好能帮我举个例子，把梁分为四个单元。

谢谢了！！

zsq-w

1  请lz修改标题  标题必须显示主题
2  这个问题一般的有限元入门书籍均有，可以参考。有些振动力学书籍也有说。

ZFLY

2# zsq-w

哪位大侠帮我看看我这个求悬臂梁一致质量矩阵的matlab程序对不对。
悬臂梁分为四个单元，每个节点只考虑水平位移。每个单元的长度为h(i).
谢谢了！
M=zeros(4);
M(4,4)=156*h(4);
for i=1:3
         M(i,i)=156*(h(i)+h(i+1));
         M(i,i+1)=54*h(i+1);
         M(i+1,i)=54*h(i+1);
end
M=density*pi*0.16^2/420*M; %直径为320mm