提取整体刚度矩阵、质量矩阵及阻尼矩阵的简单方法

glise · 发表于 2006-3-11 19:39

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

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

x

本帖最后由 Rainyboy 于 2011-4-11 11:35 编辑

其原理很简单，即使用ansys的超单元即可解决问题。定义超单元，然后列出超单元的刚度矩阵即可。
面是一个小例题，自可明白。
/prep7
k,1
k,2,3000
l,1,2
et,1,beam3
mp,ex,1,2e5
mp,prxy,1,0.3
r,1,5000,2e7,200
lesize,all,,,10
lmesh,all
finish
!----以上正常建立模型，不必施加约束和荷载
/solu
antype,7 !substructuring分析类型
seopt,matname,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等)
nsel,all !选择所有节点
m,all,all !定义所有节点自由度为主自由度
solve !求解
selist,matname,3 !列出整体刚度矩阵

shiwencai · 发表于 2006-5-11 19:30

本帖最后由 wdhd 于 2016-3-4 09:58 编辑

　　

我试过这个算例了。。。。。
但是selist列出来的数据有些乱阿，看不懂啊。。。。[em04]
请指教。。。。

另外对不单元节点较多时，则质量刚度矩阵岂不是很大吗?这样我们为了计算，怎么处理呢?

[em06]

sysh320 · 发表于 2006-5-11 19:35

本帖最后由 wdhd 于 2016-3-4 09:59 编辑

　　

楼主能不能再举几个例子?

zhangload · 发表于 2007-1-11 21:21

原帖由 glise 于 2006-3-11 19:39 发表
其原理很简单，即使用ansys的超单元即可解决问题。定义超单元，然后列出超单元的刚度矩阵即可。 面是一个小例题，自可明白。 /prep7 k,1 k,2,3000 l,1,2 &l ...

得出的数据不知道使什么意思？

再求解时没设置固定边界。。所以对所求的刚度正确性表示质疑。。。。（如果设置了他会报错）

另外，，对于网格稍大些的情况就直接推出了，比较费内存。。。

zhenxinghe · 发表于 2008-3-27 22:45

请问搂主：为什么上面的例子里面无法输出整体质量矩阵啊，
我把seopt,matname,1 改为seopt,matname,2 也不行

风自南来 · 发表于 2008-3-29 14:51

这个例子中没有定义材料密度
就没有质量矩阵的吧

账号		自动登录	找回密码
密码			我要加入

[综合讨论] 提取整体刚度矩阵、质量矩阵及阻尼矩阵的简单方法

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

本帖被以下淘专辑推荐:

怎么读懂selist的数据阿。。。。

回复 4楼的帖子

[综合讨论] 提取整体刚度矩阵、质量矩阵及阻尼矩阵的简单方法

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

本帖被以下淘专辑推荐:

怎么读懂selist的数据阿。。。。

回复 4楼 的帖子

回复 4楼的帖子