*HEADING
: DEFORMATION OF PLATE BY A RIGID PUNCH 
MODEL 1 CONTAINS USER SUBROUTINE DEFINED RIGID SURFACE
MODEL 2 CONTAINS AN ANALYTICAL RIGID SURFACE
CALLS USER SUBROUTINE A
** 
**     MODEL DEFINITION 1
**
*NODE
1,0,0,
11,10,0
101,0,1
111,10,1
*NSET,NSET=N601
 601,
*NGEN, NSET=BODY1
1,11
101,111
*ELEMENT, TYPE=CAX4, ELSET=BODY1
1,1,2,102,101
*ELGEN, ELSET=BODY1
1,10,1,1
*SOLID SECTION, ELSET=BODY1, MATERIAL=ELASTIC1
1., 
*MATERIAL, NAME=ELASTIC1
*ELASTIC
3.E5, 0.3
*NODE, NSET=CYLSURF1
99991, 0.0, 10.0
*surface, NAME=SECONDARY1
BODY1,S3
*rigid body,analytical surface=MAIN1, REFNOD=99991 
*surface,NAME=MAIN1, TYPE=USER
*CONTACT PAIR, INTERACTION=SURF1
SECONDARY1,MAIN1
*SURFACE INTERACTION, NAME=SURF1
*BOUNDARY
99991,1,1
99991,3,6
11,1,3
**
**     MODEL DEFINITION 2
**
*NODE
501,0,10,
511,10,10
601,0,11
611,10,11
*NGEN, NSET=BODY2
501,511
601,611
*ELEMENT, TYPE=CAX4, ELSET=BODY2
51,501,502,602,601
*ELGEN, ELSET=BODY2
51,10,1,1
*SOLID SECTION, ELSET=BODY2, MATERIAL=ELASTIC2
1., 
*MATERIAL, NAME=ELASTIC2
*ELASTIC
3.E5, 0.3
*NODE, NSET=CYLSURF2
99992, 0.0, 20.0
*surface, NAME=SECONDARY2
BODY2,S3
*rigid body,analytical surface=MAIN2, REFNOD=99992
*surface,NAME=MAIN2, TYPE=SEGMENTS
START,6.8301,17.8318
LINE, 4.3301,13.5
CIRCL,0,11,0,16.
*CONTACT PAIR, INTERACTION=SURF2,ADJUST=N601
SECONDARY2,MAIN2
*SURFACE INTERACTION, NAME=SURF2
*BOUNDARY
99992,1,1
99992,3,6
511,1,3
**
** STEP DEFINITION
**
*STEP, INC=40, NLGEOM
*STATIC
0.025,1.
*CLOAD
99991,2,-150000
99992,2,-150000
*NODE PRINT
U, 
*NODE FILE
U, 
*ELPRINT
S, 
*ELFILE
S, 
*END STEP