文档介绍：1 / 5
anｓyｓ提取刚度矩阵的三种方法
注:本案例借鉴了王新敏老师编著的《aｎsys工程结构数值分析》部分内容
结构刚度矩阵K
2 / 5
（1)单元刚阵的提取 --－-－详见pａge35６
/ＰＲEP7
ET，1，BＥAM3
MP，EX,１,2E5
R,1，1E-２，32E-5,
N,1
N，2，０,４
N,3，４,４
Ｎ,4,4,０
EN,1,2,3
EN,2,1，２
EN，3,４，３
Ｆ,2,FX，５
ＳＦＢＥAＭ,1,1,ＰＲEＳ,10,,,,２,-1
SFBEAＭ，2，1，ＰREＳ,3
D，1，ＡLL
Ｄ，4，ALＬ
/SＯLU
/OUTPUT,ELEMＳTＩFＦ,TXT
/DＥＢUＧ，-1，,，1
SOＬVE
/OUTPＵT
FＩＮIＳＨ
（２)用ＨＢＭＡT提取原始刚度矩阵和节点载荷
ＦINISH
/ＣLEＡR
/FIＬＮＡＭE，ＨＢＦIＬE
／PRＥP７
ＥＴ,１,BEＡＭ3
MP，EＸ，１,2Ｅ5
R,１,１Ｅ-2，3２E－5,0．５
N，1
Ｎ,2，0，4
N,3，4,4
N，4,4，0
EN，１,２，３
ＥN，2,１,2
ＥN,3,4,３
F,2,FX,5
SFBEAM,1，1,ＰＲＥS，1０,，,,2,-1
3 / 5
SFBＥＡＭ,2，1,PＲES，3
/ＳOLU
WＲFＵLL,1
SOＬVE
FINISH
／ＡUX2
ＨＢＭAＴ,HB
FIＮISＨ
用文本编辑器打开HB可看到用Haｒｗｅll-Boｅｉng格式记录的文件（仅含非零项的上三角阵,刚度矩阵为对称矩阵)。矩阵有12行12列，3３个非零元素
(3）用HBMAT提取结构刚度矩阵
与(2）相同,但施加约束条件即可。
/ＰＲＥP7
ＥT，１,BEＡM３
MＰ,EＸ,1，2Ｅ5
R，１,１Ｅ－2,３2E－５,0．５
N,1
Ｎ,2,0，4
N，3,4,4
Ｎ,4,4，０
EN,1，２,3
EN，2,１,2
EN,3,4，3
Ｆ,2,FX,５
ＳＦＢEＡM,1，１,PＲEＳ,10,,，,2，-１
ＳFＢＥAM，2，１，PＲＥS,３
D，1,ＡLL
Ｄ,４,ALL
/SＯLU
WRＦULL，1
SOLＶE
FＩNISH
／AＵX2
FＩＬE,文件名缺省,FULL
HＢＭAＴ,文件名缺省,ＴXT,ASＣII,ＳＴIＦF,YES
ＦINIＳH
其意义同上，但数值有变化。结构刚度矩阵为6行６列,有１５个非零元素,6个节点载荷,
4 / 5
（5）用超单元提取原始刚度矩阵。
/PＲＥP7
ET,1，BEＡＭ3
MP，EX,１,2E5
R,1,1Ｅ-２,3２E-5,０．5
N,1
N,2,０，4
N,3,4，４