ansys提取刚度矩阵案例（命令流及矩阵文件）

本文格式为Word版，下载可任意编辑——ansys提取刚度矩阵案例（命令流及矩阵文件）交流：516351659

ansys提取刚度矩阵的三种方法

注：本案例借鉴了王新敏老师编著的《ansys工程结构数值分析》部分内容

结构刚度矩阵K

交流：516351659

(1)单元刚阵的提取详见page356/PREP7ET,1,BEAM3MP,EX,1,2E5

R,1,1E-2,32E-5,0.5N,1N,2,0,4N,3,4,4N,4,4,0EN,1,2,3EN,2,1,2EN,3,4,3F,2,FX,5

SFBEAM,1,1,PRES,10,,,,2,-1SFBEAM,2,1,PRES,3D,1,ALLD,4,ALL/SOLU

/OUTPUT,ELEMSTIFF,TXT/DEBUG,-1,,,1SOLVE/OUTPUTFINISH

用任一文本编辑器开启ELEMSTIFF.TXT文件可得到单元刚度矩阵

(2)用HBMAT提取原始刚度矩阵和节点载荷FINISH/CLEAR

/FILNAME，HBFILE/PREP7ET,1,BEAM3MP,EX,1,2E5

R,1,1E-2,32E-5,0.5N,1N,2,0,4N,3,4,4N,4,4,0EN,1,2,3EN,2,1,2EN,3,4,3F,2,FX,5

SFBEAM,1,1,PRES,10,,,,2,-1

ELEMSTIFF.TXT交流：516351659

SFBEAM,2,1,PRES,3/SOLUWRFULL,1SOLVEFINISH/AUX2

FILE,HBFILE,FULL

HBMAT,HBFILE,TXT,ASCII,STIFF,YESFINISH

用文本编辑器开启HBFILE.TXT可看到用Harwell-Boeing格式记录的文件（仅含非零项的上三角阵，刚度矩阵为对称矩阵）。矩阵有12行12列，33个非零元素

（3）用HBMAT提取结构刚度矩阵

与（2）一致，但施加约束条件即可。所生成的HBMAT.TXT的前5行为/PREP7ET,1,BEAM3MP,EX,1,2E5

R,1,1E-2,32E-5,0.5N,1N,2,0,4N,3,4,4N,4,4,0EN,1,2,3EN,2,1,2EN,3,4,3F,2,FX,5

SFBEAM,1,1,PRES,10,,,,2,-1SFBEAM,2,1,PRES,3D,1,ALLD,4,ALL/SOLUWRFULL,1SOLVEFINISH/AUX2

FILE,文件名缺省,FULL

HBMAT,文件名缺省,TXT,ASCII,STIFF,YESFINISH

其意义同上，但数值有变化。结构刚度矩阵为6行6列，有15个非零元素，6个节点载荷，

HBFILE.TXTfile.TXT

交流：516351659

（5）用超单元提取原始刚度矩阵。/PREP7ET,1,BEAM3MP,EX,1,2E5

R,1,1E-2,32E-5,0.5N,1N,2,0,4N,3,4,4N,4,4,0EN,1,2,3EN,2,1,2EN,3,4,3F,2,FX,5

SFBEAM,1,1,PRES,10,,,,2,-1SFBEAM,2,1,PRES,3/SOLUANTYPE,7

SEOPT,SUBMAT,1M,ALL,ALLSOLVE

SELIST,SUBMAT,3

该原始刚度矩阵为12行12列，文件按行ROW1~ROW12显示

(6)用超单元提取结构刚度矩阵，与（5）一致但要施加约束条件/PREP7ET,1,BEAM3MP,EX,1,2E5

R,1,1E-2,32E-5,0.5N,1N,2,0,4N,3,4,4N,4,4,0EN,1,2,3EN,2,1,2EN,3,4,3F,2,FX,5

SFBEAM,1,1,PRES,10,,,,2,-1SFBEAM,2,1,PRES,3D,1,ALLD,4,ALL/SOLUANTYPE,7

SELIST(原始刚度).txt交流：516351659

SEOPT,SUBMAT,1M,ALL,ALLSOLVE

SELIST,SUBMAT,3