*HEADING
 TEST OF *PARAMETER SHAPE VARIATION, AXISYMMETRIC, CONTACT, HYBRID ---   
*PARAMETER
par1 = 1.
*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=OBJECT
 NTOP,NBOT,4,10
*NSET, NSET=NCENTER
 1,11,21,31,41
*NSET, NSET=NEND
 5,15,25,35,45
** SELECT THE LAST INCR OF A STATIC ANALYSIS USING CONTACT
** TO DEFINE GEOMETRIC VARIATION
*PARAMETER SHAPE, PARAMETER=par1, FILE=ximpc1, STEP=1, NSET=OBJECT
*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
*STATIC
*NODE FILE, NSET=OBJECT
 COORD,
*END STEP