有限元方程的求解

skyfing12

各位兄弟:
在对一个空间三维实体进行六面体单元划分后,单元刚度矩阵也已经求出,在下一步如何形成总体刚度矩阵,形成后,如何求解?
各位大侠帮忙!

xinyuxf

%%%%%%%%%装配单元刚度矩阵%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
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年。

图书馆中应该有

skyfing12

强烈感谢院长xinyuxf!
:lol :lol :lol :lol
感谢!!
真正的感谢!!

fgp1

原帖由 xinyuxf 于 2006-12-28 19:12 发表
%%%%%%%%%装配单元刚度矩阵%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function SK = SheetASM(K,EK,lmd)  %K总刚度矩阵，EK单元刚度矩阵，lmd节点编号矩阵
for i = 1:1:4
    for j = 1:1 ...

good

水中鱼儿

顶一下，谢谢楼主哈！以后多分享一些这样得源程序就好了。

linqus

总刚组好后，
需要引入边界、外荷载向量等等

建议先做个小算例（手算），都清楚明白后再编程计算。
推荐两个人的书：
赵更新
刘尔烈
上面有一些编程的实例，可以借鉴。

songyong

可以先编一些，节点位移是零的小程序，在编写节点位移不是零的，会好一些，如果只要实现功能的话，这样的程序还是很好编写的

songyong

从院长的程序中获益，感谢院长