%计算普通梁单元刚度矩阵函数
function [Ke] = FrameElementKe1(A,E,I,R,BarLength)
   L=BarLength;
   ke=[E*A/L       0             0          -E*A/L       0                0
       0        12*E*I/L^3      6*E*I/L^2     0       -12*E*I/L^3   6*E*I/L^2;
       0        6*E*I/L^2       4*E*I/L       0        -6*E*I/L^2	2*E*I/L ;
       -E*A/L       0             0          E*A/L       0                0
       0        -12*E*I/L^3     -6*E*I/L^2     0       12*E*I/L^3   -6*E*I/L^2;
       0        6*E*I/L^2       2*E*I/L       0        -6*E*I/L^2	4*E*I/L ;];
   
   Ke=R'*ke*R;
end