|
%%%%%%%%%装配单元刚度矩阵%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function SK = SheetASM(K,EK,lmd) %K总刚度矩阵,EK单元刚度矩阵,lmd节点编号矩阵
for i = 1:1:4
for j = 1:1:4
for p = 1:1:3
for q = 1:1:3
m = (i-1) * 3 + p;
n = (j-1) * 3 + q;
M = (lmd(i)-1) * 3 + p;
N = (lmd(j)-1) * 3 + q;
K(M,N) = K(M,N) + EK(m,n);
end
end
end
end
SK=K;
return;
skyfing12 ,你好。上次你PM给我时,正在帮助整理论坛,没有及时回,后来竟给忘了,对不住。以前在论坛有组装刚度矩阵的程序,应该能搜到。现在直接贴给你,免得搜了。
至于如何求解,具体问题具体对待,动态和静态都不一样。建议你找些书看看。推荐一本书:
MATLAB有限元分析与应用/(德) P. I. Kattan著 ,清华大学出版社 2004年。
图书馆中应该有 |
评分
-
1
查看全部评分
-
|