*HEADING
 TEST OF *IMPERFECTION, AXISYMMETRIC, CONTACT, HYBRID (INITIAL GEOM)       
 MULTIPLE NSET DEFINITIONS ARE PRINTED TO THE RESULTS FILE
*RESTART, WRITE, FREQUENCY=1
*NODE, NSET=NCORNER
 1,  0., 2.0
 5,  4., 2.0
 41, 0., 0.0
 45, 4., 0.0
*NODE, NSET=N99
 99, 0.,4.5
*NGEN, NSET=NTOP
 1,5
*NGEN, NSET=NBOT
 41,45
*NFILL, NSET=NALL
 NTOP,NBOT,4,10
*NSET, NSET=NCENTER
 1,11,21,31,41
*NSET, NSET=NEND
 5,15,25,35,45
*ELEMENT, TYPE=CAX4H, ELSET=MAT1
 1,11,12,2,1
*ELGEN, ELSET=MAT1
 1,4,1,1,4,10,5
*ELSET, ELSET=ETOP
 1,4,1
*rigid body,analytical surface=BSURF,REF NODE=99
*surface,TYPE=SEGMENTS, NAME=BSURF
 START,2.5,4.5
 CIRCL,0.0,2.0,0,4.5
 CIRCL,-2.5,4.5,0,4.5
*surface, NAME=ASURF
 ETOP,S3
*CONTACT PAIR, INTERACTION=ROUGH
 ASURF,BSURF
*SURFACE INTERACTION, NAME=ROUGH
*FRICTION
 0.8, 
*SOLID SECTION, ELSET=MAT1, MATERIAL=MAT1, ORIENTATION=RECT
*ORIENTATION, NAME=RECT, SYSTEM=RECTANGULAR
 1.,0.,0., 0.,1.,0.
 3,0.
*MATERIAL, NAME=MAT1
*DENSITY
 1.0, 
*ELASTIC
 200,0.15
*BOUNDARY
 NBOT,1,2     
 NCENTER,1
 N99,1
 N99,6
*AMPLITUDE, NAME=RAMP
 0.,0.,1.,1.
*STEP, INC=20, NLGEOM
*STATIC
 0.1,1.0
*BOUNDARY, AMP=RAMP
 N99,2,2,-0.5
*NODE PRINT, GLOBAL=YES
 U,
 COORD,
**
** PRINT MULTIPLE NSET DEFINITIONS
**
*NODE FILE, NSET=NTOP, GLOBAL=YES
 U,COORD
*OUTPUT,FIELD
*NODE OUTPUT,NSET=NTOP
U,COORD
*NODE FILE, NSET=NALL, GLOBAL=YES
 U, 
*OUTPUT,FIELD
*NODE OUTPUT,NSET=NALL
U,
*NODE FILE, NSET=NCORNER, GLOBAL=YES
 U, 
*OUTPUT,FIELD
*NODE OUTPUT,NSET=NCORNER
U,
*END STEP