24小时热门版块排行榜

返回列表

海鸥的飞翔

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 17.4
散金: 30
帖子: 47
在线: 27.5小时
虫号: 3875028
注册: 2015-05-17
专业: 机械动力学

[求助] 这样组装平面桁架的质量矩阵对不对呢？

%单元刚度矩阵
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

回复此楼

» 猜你喜欢

聘U V热熔胶研究人员已经有10人回复
职称评审没过，求安慰已经有7人回复
求助文献已经有3人回复
垃圾破二本职称评审标准已经有10人回复
投稿返修后收到这样的回复，还有希望吗已经有8人回复
三无产品还有机会吗已经有6人回复
博士申请都是内定的吗？已经有14人回复
谈谈两天一夜的“延安行” 已经有13人回复
氨基封端PDMS和HDI反应快速固化已经有11人回复
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛已经有3人回复

1楼2015-09-14 16:05:44

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主海鸥的飞翔的主题更新

返回列表