%计算弹性地基梁单元刚度矩阵函数
function [Ke] = FrameElementKe2(A,E,I,R,BarLength)
    global m D
   L=BarLength;
   k=m*D*L;
   k1=k;
   k2=k;
   ke=[E*A/L       0             0          -E*A/L       0                0
       0        12*E*I/L^3+4*k1*L/15+k2*L/12     6*E*I/L^2     0       -12*E*I/L^3+k1*L/12+k2*L/15   6*E*I/L^2;
       0        6*E*I/L^2+k1*L^2/30+k2*L^2/60      4*E*I/L       0        -6*E*I/L^2+k1*L^2/60+k2*L^2/60	2*E*I/L ;
       -E*A/L       0             0          E*A/L       0                0
       0        -12*E*I/L^3+k1*L/15+k2*L/12     -6*E*I/L^2     0       12*E*I/L^3+k1*L/12+4*k2*L/15   -6*E*I/L^2;
       0        6*E*I/L^2-k1*L^2/60-k2*L^2/60      2*E*I/L       0        -6*E*I/L^2-k1*L^2/60-k2*L^2/30	4*E*I/L ];
   Ke=R'*ke*R;
end