fini
/cle
/batch,list
/PREP7
!*
et,1,beam188

mp,ex,1,2.07e11
mp,nuxy,1,0.28
mp,dens,1,7860
!* Parameters
thick1 = 0.03810
thick  = thick1
hthick = thick/2.0
radius = 0.09525
inradius = radius - hthick
outradius = radius + hthick
length = 0.254
!*
pi=acos(-1)
!di=inradius*2
!da=outradius*2
!d2=da*da-di*di
!d4=da**4-di**4
!r,1,pi*d2/4,pi*d4/64,pi*d4/64,da,da
sectype,1,beam,ctube,,20
secdata,inradius,outradius,20
secnum,1

k,1
k,2,0,0,length
l,1,2
mat,1
esize,0.01
real,1
lmesh,1

/esha,1 $ eplot

finish
!*
/SOLU
antype,modal
modopt,qrdamp,20
omega,0,0,314.16,1
!
! Apply Coriolis effect and use default
!   rotating reference frame
coriolis,on,,,off
!
solve
finish
/post1
set,list
!set,1,1
!pldi,1