模型&结果
finish
/clear
/PREP7
g=9.8 ! gravity m/sec^2
pi=3.141592654 ! pi
/title,DIRECT DISPLACEMENT METHOD (UNIT LOAD) 386.4 IN/SEC^2 BASE INPUT
!!!!Pre-processing!!!!
RECTNG,20,,.25,,
ET,1,PLANE42
MP,EX,1,30e6,
MP,NUXY,1,.3,
MP,DENS,1,.000729,
LESIZE,2, , ,4, , , , ,1
LESIZE,4, , ,4, , , , ,1
LESIZE,1, , ,80, , , , ,1
LESIZE,3, , ,80, , , , ,1
MSHAPE,0,2D
MSHKEY,1
AMESH,all
nsel,s,loc,x,-10,.0001
d,all,all,0
allsel,all
FINISH
!!!!Modal Analysis!!!!
/SOLU
ANTYPE,2
MODOPT,LANB,6
EQSLV,SPAR !Sparse direct equation solver
MXPAND,6, , ,0
MODOPT,LANB,6,0,0, ,OFF,
SOLVE
FINISH
!!!!Postprocess Modal Analysis!!!!
/POST1
*GET,freq1,MODE,1,FREQ
*GET,freq2,MODE,2,FREQ
*GET,freq3,MODE,3,FREQ
finish
!!!!Solution Process Full Harmonic Analysis!!!!
/solu
ANTYPE,3
HROPT,FULL
HROUT,ON !
!*
EQSLV,FRONT,0, !直接求解
DDELE,1,ALL
nsel,s,loc,x,-10,.0001
d,all,ux,0
d,all,uy,1
allsel,all
HARFRQ,freq2*.5,freq3*1.2,
NSUBST,50,
KBC,1
!*
DMPRAT,.0003, !常熟阻尼比
!*
outres,all,all
solve
!!!!Postprocess Harmonic Analysis!!!!
finish
/POST26
NSOL,2,83,U,Y, ! raw results
!PROD, IR, IA, IB, IC, Name, --, --, FACTA, FACTB, FACTC
!IR = (FACTA x IA) x (FACTB x IB) x (FACTC x IC)
!3=2*pi*1x2*pi*1
PROD,3,1,1, , , , ,2*pi,2*pi,1, !Omega^2 圆频率
!
PROD,4,2,3 , , , , ,1,1, ! accelerations: raw results*accel
PROD,5,2,3 , , , , ,1/g,1 ! g's: (raw results*accel)/g