/TITLE, Your name and date here, cylind1.a
  ! Axial load on cylinder
  ! PLANE42 is 2D, 4-node structural solid
  ! 2in dia on top, 2in dia on bottom, 12in high
  ! axial node numbering
  ! exact answer is

/PREP7               ! preprocessor section
ET,1,PLANE42,,,1     ! axisymmetric option
for = 20e3           ! force at upper end
len = 12             ! length
dia = 2              ! diameter
pi = 4*atan(1)
area = pi/4*dia**2
press = for/area     ! force conv to pressure
elax = 8             ! axial number of elements
elra = 4             ! radial number of elements
nax = elax+1         ! axial number of nodes
nrad = elra+1        ! radial number of nodes

MP,EX,1,29E6         ! Young-s modulus, steel
MP,GXY,1,11.2E6      ! shear modulus

N,1,0                ! bottom center
N,nax,0,len          ! upper-left node
FILL,1,nax           ! axial nodes
NGEN,nrad,10,1,nax,1,(dia/2/elra),0 ! axial nodes
/PNUM,node,1          ! plot the nodes
! NPLOT
APLOT,all           ! plot all areas
/PLOPTS,INFO,OFF    ! Removes entire legend column
/PLOPTS,FRAME,OFF   ! Removes frame

TYPE,1               ! element type (ET) 1
E,1,11,12,2          ! lower elem

EPLOT