!       TITLE:  harmonic_solv.mac (SOLVE)
/solu
!
nsel,s,loc,z,0,0,,
ddele,all
!
antype,harmic           ! initiate harmonic analysis
hropt,msup,10,1         ! mode-superposition option
hrout,off           ! results as amplitude & phase angle
outpr,basic,1
nsubst,1000         ! 1000 intervals in frequency range
harfrq,0,2500           ! frequency range from 0 to 2500 Hz
kbc,1               ! step boundary condition
d,all,uy,2.00
!nsel,s,loc,z,20,,
!nsel,r,loc,y,0.5,,
!*get,nnum1,node,,num,max
!f,all,fy,20.0          ! apply load amplitude
nsel,all
solve
!
finish
!
/post26
!
csys,0              ! Selected node 1
nsel,s,loc,z,20,,
nsel,r,loc,y,0.5,,
*get,nnum1,node,,num,max
!
nsel,s,loc,z,4,,        ! Selected node 2
nsel,r,loc,y,0.5,,
*get,nnum2,node,,num,max
!
nsel,s,loc,z,8,,        ! Selected node 3
nsel,r,loc,y,0.5,,
*get,nnum3,node,,num,max
!
nsel,s,loc,z,12,,       ! Selected node 4
nsel,r,loc,y,0.5,,
*get,nnum4,node,,num,max
!
nsel,s,loc,z,16,,       ! Selected node 5
nsel,r,loc,y,0.5,,
*get,nnum5,node,,num,max
!
file,harmonic,rfrq
nsol,2,nnum1,u,y,disp1
nsol,3,nnum2,u,y,disp2
nsol,4,nnum3,u,y,disp3
nsol,5,nnum4,u,y,disp4
nsol,6,nnum5,u,y,disp5
plcplx,3
/grid,1
/axlab,x,freq
/axlab,y,disp
plvar,2,3,4,5,6
finish
!