FEM-Course-Matlab/12+13框架结构动力分析（模态分析+时程分析）/FrameElementMe.m

%计算单元刚度矩阵函数
function [Me] = FrameElementMe(A,rho,R,BarLength)
   L=BarLength;
   me=rho*A*L/2*[1     0     0     0     0     0
                 0     1     0     0     0     0
                 0     0     0     0     0     0
                 0     0     0     1     0     0
                 0     0     0     0     1     0
                 0     0     0     0     0     0];
   
   Me=R'*me*R;
end
% diag([1 1 0 1 1 0])