| 查看: 3349 | 回复: 0 | ||
[求助]
这样组装平面桁架的质量矩阵对不对呢？
|
|
%单元刚度矩阵 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;]; %单元质量矩阵 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); %组装质量矩阵 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) %该函数进行单元刚度矩阵的组装 %输入单元刚度矩阵 k，单元的节点编号 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;  qq.png |
回复此楼» 猜你喜欢
专硕085601求调剂
已经有6人回复
-
282求调剂 不挑专业 求收留
已经有6人回复
-
材料专硕 调剂
已经有7人回复
-
285求调剂
已经有6人回复
-
280求调剂
已经有17人回复
-
工科 267求调剂
已经有6人回复
-
调剂
已经有5人回复
-
材料调剂
已经有13人回复
-
315求调剂
已经有3人回复
-
270求调剂
已经有4人回复
