ansys提取剛度矩陣案例（命令流及矩陣文件）

本文格式為Word版，下載可任意編輯——ansys提取剛度矩陣案例（命令流及矩陣文件）交流：516351659

ansys提取剛度矩陣的三種方法

注：本案例借鑒了王新敏老師編著的《ansys工程結構數(shù)值分析》部分內(nèi)容

結構剛度矩陣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提取原始剛度矩陣和節(jié)點載荷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

其意義同上，但數(shù)值有變化。結構剛度矩陣為6行6列，有15個非零元素，6個節(jié)點載荷，

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