ansys提取剛度矩陣案例(命令流及矩陣文件)

交流QQ：516351659ansys提取剛度矩陣的三種方法注：本案例借鑒了王新敏老師編著的《ansys工程結(jié)構(gòu)數(shù)值分析》部分內(nèi)容結(jié)構(gòu)剛度矩陣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提取原始剛度矩陣和節(jié)點(diǎn)載荷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格式記錄的文件（僅含非零項(xiàng)的上三角陣，剛度矩陣為對稱矩陣）。矩陣有12行12列，33個非零元素（3）用HBMAT提取結(jié)構(gòu)剛度矩陣與（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其意義同上，但數(shù)值有變化。結(jié)構(gòu)剛度矩陣為6行6列，有15個非零元素，6個節(jié)點(diǎn)載荷，（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)用超單元提取結(jié)構(gòu)剛度矩陣，與（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