声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2535|回复: 7

[计算力学] 有限元方程的求解

[复制链接]
发表于 2006-12-28 17:12 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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

使用道具 举报

发表于 2006-12-28 19:12 | 显示全部楼层
%%%%%%%%%装配单元刚度矩阵%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
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

查看全部评分

 楼主| 发表于 2006-12-28 20:53 | 显示全部楼层

感谢

强烈感谢院长xinyuxf!
:lol :lol :lol :lol
感谢!!
真正的感谢!!
发表于 2006-12-29 09:40 | 显示全部楼层
原帖由 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
发表于 2007-1-5 10:50 | 显示全部楼层
顶一下,谢谢楼主哈!以后多分享一些这样得源程序就好了。
发表于 2007-1-8 10:10 | 显示全部楼层
总刚组好后,
需要引入边界、外荷载向量等等

建议先做个小算例(手算),都清楚明白后再编程计算。
推荐两个人的书:
赵更新
刘尔烈
上面有一些编程的实例,可以借鉴。
发表于 2007-1-8 18:37 | 显示全部楼层
可以先编一些,节点位移是零的小程序,在编写节点位移不是零的,会好一些,如果只要实现功能的话,这样的程序还是很好编写的
发表于 2007-1-8 18:40 | 显示全部楼层
从院长的程序中获益,感谢院长
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 09:34 , Processed in 0.082382 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表