%單元剛度矩陣
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+08;
-0.2832000E+08 -0.2124000E+08 0.2832000E+08 0.2124000E+08;];
k4=[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;];
%單元質(zhì)量矩陣
m1=[0.9466667E-01 0.0000000E+00 0.4733333E-01 0.0000000E+00;
0.0000000E+00 0.9466667E-01 0.0000000E+00 0.4733333E-01;
0.4733333E-01 0.0000000E+00 0.9466667E-01 0.0000000E+00;
0.0000000E+00 0.4733333E-01 0.0000000E+00 0.9466667E-01;];
m2=[0.7100000E-01 0.0000000E+00 0.3550000E-01 0.0000000E+00;
0.0000000E+00 0.7100000E-01 0.0000000E+00 0.3550000E-01;
0.3550000E-01 0.0000000E+00 0.7100000E-01 0.0000000E+00;
0.0000000E+00 0.3550000E-01 0.0000000E+00 0.7100000E-01;];
m3=[0.1183333E+00 0.0000000E+00 0.5916667E-01 0.0000000E+00;
0.0000000E+00 0.1183333E+00 0.0000000E+00 0.5916667E-01;
0.5916667E-01 0.0000000E+00 0.1183333E+00 0.0000000E+00;
0.0000000E+00 0.5916667E-01 0.0000000E+00 0.1183333E+00;];
m4=[0.9466667E-01 0.0000000E+00 0.4733333E-01 0.0000000E+00;
0.0000000E+00 0.9466667E-01 0.0000000E+00 0.4733333E-01;
0.4733333E-01 0.0000000E+00 0.9466667E-01 0.0000000E+00;
0.0000000E+00 0.4733333E-01 0.0000000E+00 0.9466667E-01;];
%組裝剛度矩陣
KK=zeros(8,8);
KK=Bar2D2Node_Assembly(KK,k1,1,2);
KK=Bar2D2Node_Assembly(KK,k2,3,2);
KK=Bar2D2Node_Assembly(KK,k3,1,3);
KK=Bar2D2Node_Assembly(KK,k4,3,4);
%組裝質(zhì)量矩陣
MM=zeros(8,8);
MM=Bar2D2Node_Assembly(MM,m1,1,2);
MM=Bar2D2Node_Assembly(MM,m2,3,2);
MM=Bar2D2Node_Assembly(MM,m3,1,3);
MM=Bar2D2Node_Assembly(MM,m4,3,4);
function z = Bar2D2Node_Assembly(KK,k,i,j)
%該函數(shù)進(jìn)行單元剛度矩陣的組裝
%輸入單元剛度矩陣 k,單元的節(jié)點編號 i、 j
%輸出整體剛度矩陣 KK
%--------------------------------------------------------
DOF(1)=2*i-1;
DOF(2)=2*i;
DOF(3)=2*j-1;
DOF(4)=2*j;
for n1=1:4
for n2=1:4
KK(DOF(n1),DOF(n2))= KK(DOF(n1),DOF(n2))+k(n1,n2);
end
end
z=KK;
![這樣組裝平面桁架的質(zhì)量矩陣對不對呢?]()
qq.png |