ansys提取刚度矩阵案例(命令流及矩阵文件)

交流QQ：516351659ansys提取刚度矩阵的三种方法注：本案例借鉴了王新敏老师编著的《ansys工程结构数值分析》部分内容结构刚度矩阵K(1)单元刚阵的提取详见page356/PREP7ET,1,BEAM3MP,EX,1,2E5R,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,5SFBEAM,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,2E5R,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,5SFBEAM,1,1,PRES,10,,,,2,-1SFBEAM,2,1,PRES,3/SOLUWRFULL,1SOLVEFINISH/AUX2FILE,HBFILE,FULLHBMAT,HBFILE,TXT,ASCII,STIFF,YESFINISH用文本编辑器打开HBFILE.TXT可看到用Harwell-Boeing格式记录的文件（仅含非零项的上三角阵，刚度矩阵为对称矩阵）。矩阵有12行12列，33个非零元素（3）用HBMAT提取结构刚度矩阵与（2）相同，但施加约束条件即可。所生成的HBMAT.TXT的前5行为/PREP7ET,1,BEAM3MP,EX,1,2E5R,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,5SFBEAM,1,1,PRES,10,,,,2,-1SFBEAM,2,1,PRES,3D,1,ALLD,4,ALL/SOLUWRFULL,1SOLVEFINISH/AUX2FILE,文件名缺省,FULLHBMAT,文件名缺省,TXT,ASCII,STIFF,YESFINISH其意义同上，但数值有变化。结构刚度矩阵为6行6列，有15个非零元素，6个节点载荷，（5）用超单元提取原始刚度矩阵。/PREP7ET,1,BEAM3MP,EX,1,2E5R,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,5SFBEAM,1,1,PRES,10,,,,2,-1SFBEAM,2,1,PRES,3/SOLUANTYPE,7SEOPT,SUBMAT,1M,ALL,ALLSOLVESELIST,SUBMAT,3该原始刚度矩阵为12行12列，文件按行ROW1~ROW12显示(6)用超单元提取结构刚度矩阵，与（5）相同但要施加约束条件/PREP7ET,1,BEAM3MP,EX,1,2E5R,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,5SFBEAM,1,1,PRES,10,,,,2,-1SFBEAM,2,1,PRES,3D,1,ALLD,4,ALL/SOLUANTYPE,7SEOPT,SUB