*HEADING
 FORMING OF A SHEET AROUND A CYLINDER - RSURFU 
*RESTART, WRITE, F=10
**
** SHEET DEFINITION
**
*NODE
1,    0.0,  5.0,  15.0
11,  10.0,  5.0,  15.0
101,  0.0,  5.0,  20.0
111, 10.0,  5.0,  20.0
*NGEN, NSET=EDGE1
1,101,20
*NGEN, NSET=EDGE2
11,111,20
*NSET, NSET=BOTEDG, GENERATE
1,11
*NFILL, NSET=ALL
EDGE1,EDGE2,10,1
*ELEMENT, TYPE=S4R, ELSET=SHELL
1,1,21,22,2
*ELGEN, ELSET=SHELL
1,10,1,1,5,20,100
*SHELL SECTION, ELSET=SHELL, MATERIAL=ELASTIC
0.25, 
*MATERIAL, NAME=ELASTIC
*ELASTIC
3.0E6, 0.3
**
** RIGID surface
**
*NODE, NSET=CYLSURF
10000, 0.0, 0.0, 20.0
*NSET,NSET=CYL,GENERATE
2,11
22,31
42,51
62,71
82,91
102,111
*surface,type=node, NAME=SLAVE
CYL,
*CONTACT PAIR, INTERACTION=SURF, HCRIT=1.0E10
SLAVE,MASTER
*SURFACE INTERACTION, NAME=SURF
*rigid body,analytical surface=MASTER,REFNOD=10000 
*surface,NAME=MASTER, TYPE=USER
**
** BOUNDARY CONDITIONS
**
*BOUNDARY
EDGE1,1,6
CYLSURF,1,6
**
** STEP DEFINITION
**
*STEP, NLGEOM
*STATIC
0.1,1.0
*DLOAD
SHELL,P,-25.0
*PRINT, CONTACT=YES
*NODE PRINT,FREQ=999
*EL PRINT,FREQ=999
*NODE FILE, NSET=BOTEDG, FREQ=999
U, 
*ENDSTEP
*STEP, INC=40, NLGEOM
*STATIC
0.01,1.0
*DLOAD,op=mod
SHELL,P,-700.0
*PRINT, CONTACT=YES
*NODE PRINT,FREQ=999
*EL PRINT,FREQ=999
*NODE FILE, NSET=BOTEDG, FREQ=999
U, 
*ENDSTEP