!   TITLE:  aaav_harmonic_solv.mac (AAAV)
!
/prep7
!
/title,AAAV Harmonic Analysis / Aluminum
aclear,all
ddele,all
numcmp,all
modmsh,detach
emid,remove,all
ndele,all
numcmp,all
eplot
finish
!
/solu
!
maxf=500            ! MAXIMUM FREQUENCY IN RANGE
minf=100            ! MINIMUM FREQUENCY IN RANGE
fstep=40            ! NUMBER OF FREQUENCIES PER STEP
numinc=(maxf-minf)/fstep    ! NUMBER OF INCREMENTS IN FREQUENCY RANGE (INTEGER)
g=386.4             ! ACCELERATION OF GRAVITY (in/sec^2)
numg=1              ! NUMBER OF g's ACCELERATION
!
antype,harmic       ! harmonic response analysis
hropt,full      ! full harmonic response
hrout,off       ! print results as amplitudes & phase angles
outpr,basic,1
nsubst,numinc       ! 100 intervals in range
harfrq,minf,maxf    ! frequency range 0 to 100 Hz
kbc,1
dmprat,0.05     ! step boundary condition
!
csys,gcs
nsel,s,loc,z,-8.600,-8.600,,
csys,0
d,all,uy,1.00
nsel,all
solve
finish
!
/post26
!
csys,0
!
nsel,s,loc,z,-.10,.10,          ! Housing top surface
nsel,r,loc,y,-.10,0,
nsel,r,loc,x,2.57-.10,2.57+.10,,
*get,n1,node,,num,max
!
nsel,s,node,,44958          ! Back of mirror
*get,n2,node,,num,max
!
nsel,s,node,,33786          ! Bottom along gimbal axis
*get,n3,node,,num,max
!
nsel,s,loc,z,3.280-.20,3.280+.20,   ! Right end of elevation axis
nsel,r,loc,y,-4.592-.20,-4.592+.20,
nsel,r,loc,x,1.781-.20,1.781+.20,
*get,n4,node,,num,max
!
nsel,s,node,,40561          ! Left end of elevation axis
*get,n5,node,,num,max
!
nsel,s,node,,42311          ! Center of mirror along axis
*get,n6,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
nsol,7,n6,u,y,disp6
!
/grid,1
/axlab,x,freq (Hz)
/axlab,y,displacement (y)
prnsol,u,y
plvar,2,3,4,5,6,7
!
finish