*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
101,-8.,-8,0
111,8,-8,0
1101,-8,8,0
1111,8,8,0
*NGEN, NSET=FRONT1
101,111
*NGEN, NSET=BEHIND1
1101,1111
*NFILL,NSET=NALL1
FRONT1, BEHIND1, 10, 100
*NSET, NSET=SIDES1,GEN
101,1101,100
111,1111,100
*NSET, NSET=BOUND1
FRONT1,BEHIND1,SIDES1
*ELEMENT, TYPE=S4R, ELSET=EALL1
1,101,102,202,201
*ELGEN, ELSET=EALL1
1,10,1,1,10,100,10
*SHELL SECTION, ELSET=EALL1, MATERIAL=ELASTIC1
1., 
*MATERIAL, NAME=ELASTIC1
*ELASTIC
3.0E5, 0.3
*NODE, NSET=CYLSURF1
99991, 0.0, 10.0
*surface, NAME=SLAVE1
EALL1,SPOS
*rigid body,analytical surface=MASTER1, REFNOD=99991 
*surface,NAME=MASTER1, TYPE=USER
*CONTACT PAIR, INTERACTION=SURF1
SLAVE1,MASTER1
*SURFACE INTERACTION, NAME=SURF1
*BOUNDARY
99991,1,2
99991,4,6
BOUND1,1,3
**
**
**   MODEL DEFINITION 2
*NODE
2101,12,-8,0
2111,28,-8,0
3101,12,8,0
3111,28,8,0
*NGEN, NSET=FRONT2
2101,2111
*NGEN, NSET=BEHIND2
3101,3111
*NFILL,NSET=NALL2
FRONT2, BEHIND2, 10, 100
*NSET, NSET=SIDES2,GEN
2101,3101,100
2111,3111,100
*NSET, NSET=BOUND2
FRONT2,BEHIND2,SIDES2
*ELEMENT, TYPE=S4R, ELSET=EALL2
201,2101,2102,2202,2201
*ELGEN, ELSET=EALL2
201,10,1,1,10,100,10
*SHELL SECTION, ELSET=EALL2, MATERIAL=ELASTIC2
1., 
*MATERIAL, NAME=ELASTIC2
*ELASTIC
3.0E5, 0.3
*NODE, NSET=CYLSURF2
99992, 12.0, 10.0
*surface, NAME=SLAVE2
EALL2,SPOS
*rigid body,analytical surface=MASTER2, REFNOD=99992
*surface,NAME=MASTER2, TYPE=REVOLUTION
20,0,0,20,0,1
START,6.8301,6.8301
LINE,4.3301,2.5
CIRCL,0,0,0,5
*CONTACT PAIR, INTERACTION=SURF2
SLAVE2,MASTER2
*SURFACE INTERACTION, NAME=SURF2
*BOUNDARY
99992,1,2
99992,4,6
BOUND2,1,3
**
** STEP DEFINITION
**
*STEP, INC=40, NLGEOM
*STATIC
0.25,1.
*CLOAD
99991,3,-30000
99992,3,-30000
*NODE PRINT
U, 
*NODE FILE,NSET=NALL1,FREQ=9999
U, 
*ELPRINT
S, 
*ELFILE,FREQ=9999
S, 
*END STEP