*HEADING
 DEEP DRAWING OF CYLINDRICAL CUP WITH CAX4R
*RESTART,WRITE,FREQUENCY=25
*NODE
*NGEN,NSET=BOT
*NGEN,NSET=TOP
*NSET,NSET=WRKPC
BOT,TOP
*NODE,NSET=DIE
*NODE,NSET=PUNCH
*NODE,NSET=HOLDER
*NSET,NSET=TOOLS
PUNCH,DIE,HOLDER
*NSET,NSET=CENTER
*ELEMENT,TYPE=CAX4R,ELSET=BLANK
*ELGEN,ELSET=BLANK
*ELSET,ELSET=ALL
BLANK,
*SOLID SECTION,MATERIAL=STEEL,ORIENTATION=LOCAL,
ELSET=BLANK
*ORIENTATION,NAME=LOCAL
*MATERIAL,NAME=STEEL
*ELASTIC
2.1E11,0.3
*PLASTIC,HARDENING=ISOTROPIC
  0.91294E+08, 0.00000E+00
  0.10129E+09, 0.21052E-03
  0.11129E+09, 0.52686E-03
  0.12129E+09, 0.97685E-03
  0.13129E+09, 0.15923E-02
  0.14129E+09, 0.24090E-02
  0.15129E+09, 0.34674E-02
  0.16129E+09, 0.48120E-02
*ELSET,ELSET=EDIE,GENERATE
*ELSET,ELSET=EHOLDER,GENERATE
*ELSET,ELSET=EPUNCH,GENERATE
*RIGID BODY,ANALYTICAL SURFACE=BSURF,REF NODE=100
*SURFACE,TYPE=SEGMENTS,NAME=BSURF
START,0.05125,-0.060
LINE,0.05125,-0.005
CIRCL,0.05625,0.0,0.05625,-0.005
LINE,0.1,0.0
*RIGID BODY,ANALYTICAL SURFACE=DSURF,REF NODE=300
*SURFACE,TYPE=SEGMENTS,NAME=DSURF
START,0.1,0.00082
LINE,0.05630,0.00082
CIRCL,0.05625,.00087,.05630,.00087
LINE,0.05625,.06
*RIGID BODY,ANALYTICAL SURFACE=FSURF,REF NODE=200
*SURFACE,TYPE=SEGMENTS,FILLET RADIUS=.013,
NAME=FSURF
START,0.05,0.060
LINE,0.05,2.250782E-3
CIRCL,0.0,0.001,0.0,1.001
LINE,-0.001,0.001
*SURFACE,NAME=ASURF
EDIE,S1
*SURFACE,NAME=CSURF
EHOLDER,S3
*SURFACE,NAME=ESURF
EPUNCH,S3
*CONTACT PAIR,INTERACTION=ROUGH1
ASURF,BSURF
*CONTACT PAIR,INTERACTION=ROUGH2
CSURF,DSURF
*CONTACT PAIR,INTERACTION=ROUGH3
ESURF,FSURF
*SURFACE INTERACTION,NAME=ROUGH1
*FRICTION
0.1,
*SURFACE INTERACTION,NAME=ROUGH2
*FRICTION
0.1,
*SURFACE INTERACTION,NAME=ROUGH3
*FRICTION
0.25,
*STEP,NLGEOM, UNSYMM=YES
PUSH THE BLANKHOLDER DOWN BY A PRESCRIBED
DISPLACEMENT
*STATIC
1.,1.
*BOUNDARY
CENTER,1,1
DIE,1,1
DIE,2,2
DIE,6,6
PUNCH,1,1
PUNCH,2,2
PUNCH,6,6
HOLDER,1,1
HOLDER,2,2,-1.75E-8
HOLDER,6,6
*MONITOR,NODE=200,DOF=2
*PRINT,CONTACT=YES
*NODE PRINT,NSET=TOOLS,FREQUENCY=100
COORD,U,RF
*EL PRINT,ELSET=ALL,FREQUENCY=500
S,E
*NODE FILE,NSET=TOOLS,FREQUENCY=10
U,RF
*CONTACT FILE,SLAVE=ASURF,FREQUENCY=10
*CONTACT FILE,SLAVE=CSURF,FREQUENCY=10
*CONTACT FILE,SLAVE=ESURF,FREQUENCY=10
*END STEP
*STEP,NLGEOM
APPLY PRESCRIBED FORCE ON BLANKHOLDER AND
RELEASE DISPLACEMENT
*STATIC
1.,1.
*BOUNDARY,OP=NEW
CENTER,1,1
DIE,1,1
DIE,2,2
DIE,6,6
PUNCH,1,1
PUNCH,2,2
PUNCH,6,6
HOLDER,1,1
HOLDER,6,6
*CLOAD
HOLDER,2,-100000.
*END STEP
*STEP,INC=500,NLGEOM
MOVE THE PUNCH DOWN
*STATIC
.01,1.,1.E-6
*CONTROLS,ANALYSIS=DISCONTINUOUS
*BOUNDARY,OP=NEW
CENTER,1,1
DIE,1,1
DIE,2,2
DIE,6,6
PUNCH,1,1
PUNCH,2,2,-.06
PUNCH,6,6
HOLDER,1,1
HOLDER,6,6
*CLOAD
HOLDER,2,-100000.
*END STEP
*STEP,NLGEOM
FIX ALL NODES AND REMOVE THE CONTACT SURFACES
*STATIC
1.,1.,1.,1.
*BOUNDARY,FIXED
WRKPC,1,2
*MODEL CHANGE,TYPE=CONTACT PAIR,REMOVE
ASURF,BSURF
CSURF,DSURF
ESURF,FSURF
*CLOAD,OP=NEW
HOLDER,2,0.
*END STEP
*STEP,NLGEOM, UNSYMM=NO
REPLACE THE BOUNDARY CONDITIONS BY THE REGULAR SET
*STATIC
.1,1.,1.E-6
*BOUNDARY,OP=NEW
181,2
CENTER,1,1
*END STEP