ansys质量矩阵刚度矩阵提取.docx

文档介绍：精品文档，仅供学****与交流，如有侵权请联系网站删除
【精品文档】第 1 页
ansys质量矩阵刚度矩阵提取
看了这么久了都没人回，查了一些质料终于找到答案了，，
下面提供三种方法：方便与其他程序进行接口编程1.
Which matrix you would like? element stiffness matrix or full stiffness matrix?
element stiffness is within . full stiffness matrix is within
A simple way to dump the matrix is as follow:
/aux2
fileaux2,file,emat
form,long
dump,all
2.
可以使用/DEBUG命令来得到。详细步骤参见下面的宏文件
finish
/clear
PI=
w1=3
w2=10
w3=6
w4=
r=.8
t=
/PREP7
ET,1,SHELL63
R,1,t
精品文档，仅供学****与交流，如有侵权请联系网站删除
【精品文档】第 2 页
ET,2,MASS21
R,2,500,500,500,2000,2000,2000,
UIMP,1,EX, , ,2e11
UIMP,1,NUXY, , ,,
UIMP,1,DAMP, , ,,
UIMP,1,DENS, , ,7800,
BLC4,0,0,w2,w1
ESIZE,,0,
AMESH,all
NSEL,S,LOC,X,
D,all, , , , , ,ALL, , , , ,
allsel,all
SFA,all,1,PRES,12
FINISH
/OUTPUT,cp,out,, !
/debug,-1,,,1 ! 指定输出单元矩阵
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 将输出信息送到output windows中
! ，可以看到按单元写出的质量、刚度等矩阵
3.
其原理很简单，即使用ansys的超单元即可解决问题。定义超单元，然后列出超单元的刚度矩阵即可。
面是一个小例题，自可明白。
精品文档，仅供学****与交流，如有侵权请联系网站删除
【精品文档】第 3 页
/prep7
k,1
k,2,3000
l,1,2
et,1,beam3
mp,ex,1,2e5
mp,prxy,1,
mp,dens,1,2e3
r,1,5000,2e7,200
lesize,all,,,10
lmesh,all
finish
!----以上正常建立模型，不必施加约束和荷载
/solu
antype,7 !substructuring分析类型
seopt,matname,2 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等)
nsel,all !选择所有节点
m,all,all !定义所有节点自由度为主自由度
solve !求解
selist,matname,3 !列出整体刚度矩阵
ANSYS单元和整体刚度矩阵的提取
一、单元刚度矩阵的提取
/DEBUG命令
详细说明：
finish
精品文档，仅供学****与交流，如有侵权请联系网站删除
【精品文档】第 4 页
/clear
PI=
w1=3
w2=10
w3=6
w4=
r=.8
t=
/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, , ,,
UIMP,1,DAMP, , ,,
UIMP,1,DENS, , ,7800,
BLC4,0,0,w2,w1
ESIZE,,0,
AMESH,all
NSEL,S,LOC,X,
D,all, , , , , ,ALL, , , , ,
allsel,all
SFA,all,1,PRES,12
FINISH
精品文档，仅供学****与交流，如有侵权请联系网站删除
【精品文档】第 5