/PREP7
ANTYPE,STATIC
ET,1,solid45
MP,EX,1,2.1E11 ! STEEL
MP,NUXY,1,0
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,7800
BLC4, , ,0.01,0.01,0.1
vsweep,1
/SOLUTION
DA,1,ALL,
ANTYPE,2
EQSLV,SPAR,1e-5,
MXPAND,4, , ,1
MODOPT,LANB,4,0,0, ,OFF
solve
finish
ffre1=700
ffre2=1000
! harmonic analysis
/SOLU
ANTYPE,3 ! HARMONIC RESPONSE ANALYSIS
HROPT,MSUP ! MODE SUPERPOSITION
LUMPM,1
HROPT,MSUP,3,1,
HROUT,OFF,ON,0
ACEL,100,0,0,
KBC,1 ! STEP CHANGE FORCE
HARFRQ,ffre1,ffre2 ! OBTAIN FREQUENCY EVERY EIGHT HERTZ
NSUBST,20
DMPRAT,0.04
OUTRES,NSOL,NONE,
OUTRES,NSOL,ALL,
OUTRES,STRS,ALL,
SOLVE
FINISH
/solu
EXPASS,ON
NUMEXP,ALL,
HREXP,ALL,
SOLVE
/post1
hrcplx,1,14,360,1 !将830.9HZ的实部虚部结果转换成幅值相位结果
plnsol,s,eqv,2
/post26
ANSOL,2,172,S,EQV,SEQV_2
NSOL,3,172,U,X, UX_3
PLCPLX,0 !应力幅值
XVAR,1
PLVAR,2,
PLCPLX,2 !实频图
PLVAR,3,
PLCPLX,3 !虚频图
PLVAR,3,
fini
/clear
摘自:http://blog.sina.com.cn/s/blog_4d81e5a001000btk.html