%单元刚度矩阵
k1=[0.7375000E+08 0.0000000E+00 -0.7375000E+08 0.0000000E+00;
0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00;
-0.7375000E+08 0.0000000E+00 0.7375000E+08 0.0000000E+00;
0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00;];
k2=[0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00;
0.0000000E+00 0.9833333E+08 0.0000000E+00 -0.9833333E+08;
0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00;
0.0000000E+00 -0.9833333E+08 0.0000000E+00 0.9833333E+08;];
k3=[0.3776000E+08 0.2832000E+08 -0.3776000E+08 -0.2832000E+08;
0.2832000E+08 0.2124000E+08 -0.2832000E+08 -0.2124000E+08;
-0.3776000E+08 -0.2832000E+08 0.3776000E+08 0.2832000E&