声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1283|回复: 4

[综合讨论] 质量与刚度矩阵的重组

[复制链接]
发表于 2009-10-13 09:48 | 显示全部楼层 |阅读模式

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

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

x
如果现在已经提取刚度矩阵为:
1 2 3  4 ...  10
                        1
                        2
                        3
                        4
                        ...
                        10
整体刚度矩阵按照节点1-10组装
现在想改为6 ,7 ,8,4,1,2,3,5,9,10
请问怎么操作呢,最好有例子说明
还有 在ansys中提取的刚度矩阵默认是按照节点顺序给出的么,请教了,谢谢
回复
分享到:

使用道具 举报

发表于 2010-3-15 16:34 | 显示全部楼层
呃。。。我也想知道:@(
发表于 2010-12-27 14:33 | 显示全部楼层
问一下,本人想提取整体刚度矩阵,楼主可有命令流提供?谢谢!
发表于 2010-12-30 04:01 | 显示全部楼层

来自《振动论坛FAQ(1.0版)之Ansys---pdf版》

2、提取整体刚度矩阵、质量矩阵及阻尼矩阵的三种方法
1.
Which matrix you would like? element stiffness matrix or full stiffness matrix?
element stiffness is within file.emat. full stiffness matrix is within file.full
A simple way to dump the matrix is as follow:
-------------------
/aux2
fileaux2,file,emat
form,long
dump,all
-------------------

2. 可以使用/DEBUG 命令来得到。详细步骤参见下面的宏文件
finish
/clear
PI=3.1415926
w1=3
w2=10
w3=6
w4=1.2
r=.8
t=0.08
/PREP7
!*
ET,1,SHELL63
R,1,t
ET,2,MASS21
R,2,500,500,500,2000,2000,2000,
!*
UIMP,1,EX, , ,2e11
UIMP,1,NUXY, , ,0.3,
UIMP,1,DAMP, , ,0.2,
UIMP,1,DENS, , ,7800,
BLC4,0,0,w2,w1
ESIZE,1.5,0,
AMESH,all
NSEL,S,LOC,X,0.0
D,all, , , , , ,ALL, , , , ,
allsel,all
SFA,all,1,PRES,12
FINISH
/OUTPUT,cp,out,, ! 将输出信息送到cp.out 文件
/debug,-1,,,1 ! 指定输出单元矩阵
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 将输出信息送到output windows 中
! 这时用编辑器打开cp.out 文件,可以看到按单元写出的质量、刚度等矩阵


3.其原理很简单,即使用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 !列出整体刚度矩阵


发表于 2010-12-30 04:16 | 显示全部楼层
zhuweihong2004 发表于 2009-10-13 09:48
如果现在已经提取刚度矩阵为:
1 2 3  4 ...  10
                        1

关于刚度矩阵和节点之间的关系
我记得用超单元法得到的整体刚度矩阵和节点顺序是对应的
其它方法你可以尝试做几个简单的模型试一下,比如2个节点的、3个节点的模型
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-11 18:00 , Processed in 0.061323 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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