谐响应算例

 

/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