%计算单元刚度矩阵函数
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])