24小时热门版块排行榜

>论坛更新日志 (3959)
>导师招生 (720)
>论文投稿 (397)
>休闲灌水 (218)
>硕博家园 (217)
>文献求助 (162)
>考博 (138)
>新药研发 (98)
>虫友互识 (92)
>学术会议 (86)
>公派出国 (84)
>考研 (71)
>基金申请 (68)
>论文道贺祈福 (52)
>博后之家 (46)
>教师之家 (46)

返回列表

海鸥的飞翔

新虫 (初入文坛)

应助: 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

回复此楼

» 猜你喜欢

售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有13人回复
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有18人回复
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急已经有18人回复
2027博士申请-临床检验诊断学、病原生物学、微生物方向已经有4人回复
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急已经有10人回复
售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有22人回复
越来越担心硕士毕业后没有竞争力了已经有14人回复
论文的挂名怎么评判贡献？基金申请中，排名第一的挂名文章已经有3人回复
申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请已经有3人回复
闲聊已经有6人回复

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

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

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

返回列表