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

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

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

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

結(jié)構(gòu)剛度矩陣K

交流：516351659

(1)單元?jiǎng)傟嚨奶崛≡斠妏age356/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文件可得到單元?jiǎng)偠染仃?/p>

(2)用HBMAT提取原始剛度矩陣和節(jié)點(diǎn)載荷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格式記錄的文件（僅含非零項(xiàng)的上三角陣，剛度矩陣為對稱矩陣）。矩陣有12行12列，33個(gè)非零元素

（3）用HBMAT提取結(jié)構(gòu)剛度矩陣

與（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ù)值有變化。結(jié)構(gòu)剛度矩陣為6行6列，有15個(gè)非零元素，6個(gè)節(jié)點(diǎn)載荷，

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)用超單元提取結(jié)構(gòu)剛度矩陣，與（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