!   TITLE:  harm_solv3.mac (SOLVE) 

/COM,ANSYS RELEASE  5.4    UP19970828         16:53:05    01/06/2000
/input,start,ans     ,/export/cae/ansys54/docu/,,,,,,,,,,,,,,,,1
!  
/prep7
!  
ET,1,BEAM4
!  
/units,bin
!  
UIMP,1,EX, , ,30e6,
UIMP,1,DENS, , ,.101/386.4,
UIMP,1,NUXY, , ,0.33,
!
n,1,0
n,21,20
fill
E,1,2
EGEN,20,1,1
!
pi=3.1415927
rad=0.5
!
a=pi*rad**2
izz=(1/4)*pi*rad**4
thkz=2*rad
thky=thkz
shrz=10/9
shry=10/9
!
R,1,a,izz,iyy,thkz,thky,0
rmore,0,1/12,shrz,shry
!
*get,prkey_,active,0,prkey
/nopr
!   
/VIEW, 1 ,1,1,1
/ANG, 1
/REP,FAST
/DIST, 1 ,1.371742,1
/REP,FAST
/DIST, 1 ,1.371742,1
/REP,FAST
/DIST, 1 ,1.371742,1
/REP,FAST
/DIST, 1 ,1.371742,1
/REP,FAST
/FOC, 1 ,-0.300000,,,1
/REP,FAST
!  
/title,Harmonic Analysis Trial
eplot
!
/solu
antype,harmic       ! harmonic response element
hropt,reduc     ! full harmonic response
hrout,off       ! print results as amplitudes & phase angles
outpr,basic,1
nsubst,500      ! 50 intervals in range
harfrq,0,10000      ! frequency range 0 to 100 Hz
kbc,1           ! step boundary condition
m,all,uy
!
f,1,fy,1.00
nsel,all
solve
finish
!
/post26
!
csys,0              ! Selected node 1
nsel,s,loc,x,20,,
*get,n1,node,,num,max
!
nsel,s,loc,x,4,,        ! Selected node 2
*get,n2,node,,num,max
!
nsel,s,loc,x,8,,        ! Selected node 3
*get,n3,node,,num,max
!
nsel,s,loc,x,12,,       ! Selected node 4
*get,n4,node,,num,max
!
nsel,s,loc,x,16,,       ! Selected node 5
*get,n5,node,,num,max
!
nsel,all
!
nsol,2,n1,u,y,disp1
nsol,3,n2,u,y,disp2
nsol,4,n3,u,y,disp3
nsol,5,n4,u,y,disp4
nsol,6,n5,u,y,disp5

*if,prkey_,eq,1,then
    /go
*endif

!plcplx,3
/grid,1
/axlab,x,freq
/axlab,y,disp
plvar,2,3,4,5,6
!   
finish