*HEADING
 DEEP DRAWING OF CYLINDRICAL CUP WITH SAX1  
*RESTART,WRITE,FREQUENCY=25
*NODE
101, 
181,0.1
*NGEN,NSET=BOT
101,181,2
*NSET,NSET=WRKPC,GENERATE
121,181,2
*NODE,NSET=DIE
100,0.,-0.05
*NODE,NSET=PUNCH
200,0.,.05
*NODE,NSET=HOLDER
300,0.,0.05
*NSET,NSET=TOOLS
PUNCH,DIE,HOLDER
*NSET,NSET=CENTER
101, 
*ELEMENT,TYPE=SAX1,ELSET=BLANK
201,101,121
202,121,123
*ELGEN,ELSET=BLANK
202,30,2,2 
*SHELL SECTION,MATERIAL=STEEL,ELSET=BLANK
.00082,5
*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
  0.17129E+09, 0.64921E-02
  0.18129E+09, 0.85618E-02
  0.19129E+09, 0.11080E-01
  0.20129E+09, 0.14110E-01
  0.21129E+09, 0.17723E-01
  0.22129E+09, 0.21991E-01
  0.23129E+09, 0.26994E-01
  0.24129E+09, 0.32819E-01
  0.25129E+09, 0.39556E-01
  0.26129E+09, 0.47301E-01
  0.27129E+09, 0.56159E-01
  0.28129E+09, 0.66236E-01
  0.29129E+09, 0.77648E-01
  0.30129E+09, 0.90516E-01
  0.31129E+09, 0.10497E+00
  0.32129E+09, 0.12114E+00
  0.33129E+09, 0.13916E+00
  0.34129E+09, 0.15919E+00
  0.35129E+09, 0.18138E+00
  0.36129E+09, 0.20588E+00
  0.37129E+09, 0.23287E+00
  0.38129E+09, 0.26252E+00
  0.39129E+09, 0.29502E+00
  0.40129E+09, 0.33054E+00
  0.41129E+09, 0.36929E+00
  0.42129E+09, 0.41147E+00
  0.43129E+09, 0.45729E+00
  0.44129E+09, 0.50696E+00
  0.45129E+09, 0.56073E+00
  0.46129E+09, 0.61881E+00
  0.47129E+09, 0.68145E+00
  0.48129E+09, 0.74890E+00
  0.49129E+09, 0.82142E+00
  0.50129E+09, 0.89928E+00
  0.51129E+09, 0.98274E+00
  0.52129E+09, 0.10721E+01
*RIGID BODY,ANALYTICAL SURFACE=HOLDER,REF NODE=300
*SURFACE,TYPE=SEGMENTS,NAME=HOLDER,
FILLET RADIUS=0.001
START,0.12,1.75E-8
LINE,0.05625,1.75E-8
LINE,0.05625,.06
*RIGID BODY,ANALYTICAL SURFACE=DIE,REF NODE=100
*SURFACE,TYPE=SEGMENTS,NAME=DIE,
FILLET RADIUS=0.00541
START,0.05125,-0.060
LINE,0.05125,0.0
LINE,0.12,0.0                                        
*RIGID BODY,ANALYTICAL SURFACE=PUNCH,REF NODE=200
*SURFACE,TYPE=SEGMENTS,NAME=PUNCH,
FILLET RADIUS=.01341
START,0.05,0.060
LINE,0.05,2.250782E-3
CIRCL,0.0,.0001,0.0,1.0001
LINE,-0.001,0.0001
**
** Contact with holder
** 
*ELSET,ELSET=HOLD_CON,GENERATE
222,260,2
*SURFACE,NAME=HLD_SURF
HOLD_CON,SPOS
*CONTACT PAIR,INTERACTION=FRIC1
HLD_SURF,HOLDER
*SURFACE INTERACTION,NAME=FRIC1
*SURFACE BEHAVIOR, 
PRESSURE-OVERCLOSURE=EXPONENTIAL         
1.75E-8,4.5E6
*FRICTION
0.1, 
**
** Contact with die
**
*ELSET,ELSET=DIE_CON,GENERATE
212,260,2
*SURFACE,NAME=DIE_SURF
DIE_CON,SNEG
*CONTACT PAIR,INTERACTION=FRIC2
DIE_SURF,DIE
*SURFACE INTERACTION,NAME=FRIC2
*FRICTION
0.1, 
**
** Contact with punch
**
*ELSET,ELSET=PUN_CON2,GENERATE
202,230,2
*SURFACE,NAME=PCH_SURF
PUN_CON2,SPOS
*CONTACT PAIR,INTERACTION=FRIC3
PCH_SURF,PUNCH
*SURFACE INTERACTION,NAME=FRIC3
*FRICTION
0.25,
*SURFACE,TYPE=NODE,NAME=PCH_PNT
101, 
*CONTACT PAIR,INTERACTION=FRIC4
PCH_PNT,PUNCH
*SURFACE INTERACTION,NAME=FRIC4
*FRICTION
0.0, 
**
**
*STEP,UNSYMM=YES,NLGEOM
PUSH THE BLANKHOLDER DOWN BY A PRESCRIBED 
DISPLACEMENT
*STATIC
1.,1.
*BOUNDARY
CENTER,1,1
CENTER,6,6
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=BLANK,FREQUENCY=500
S,E
STH, 
*CONTACT PRINT,FREQUENCY=500
*CONTACT FILE,FREQUENCY=500
*NODE FILE,NSET=TOOLS,FREQUENCY=10
U,RF
*END STEP
*STEP,NLGEOM, UNSYMM=YES
APPLY PRESCRIBED FORCE ON BLANKHOLDER AND 
RELEASE DISPLACEMENT
*STATIC
1.,1.
*BOUNDARY,OP=NEW
CENTER,1,1
CENTER,6,6
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,UNSYMM=YES
MOVE THE PUNCH DOWN
*STATIC
.01,1.,1.E-6
*CONTROLS,ANALYSIS=DISCONTINUOUS
*BOUNDARY,OP=NEW
CENTER,1,1
CENTER,6,6
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,INC=1,NLGEOM, UNSYMM=YES
FIX ALL NODES AND REMOVE THE CONTACT SURFACES
*STATIC
1.,1.,1.,1.
*BOUNDARY,FIXED
WRKPC,1,6
*BOUNDARY,FIXED
CENTER,1,6
*MODEL CHANGE,TYPE=CONTACT PAIR,REMOVE
HLD_SURF,HOLDER
PCH_PNT,PUNCH
PCH_SURF,PUNCH
DIE_SURF,DIE
*CLOAD,OP=NEW
HOLDER,2,0.
*END STEP
*STEP,NLGEOM, UNSYMM=YES
REPLACE THE BOUNDARY CONDITIONS BY THE REGULAR SET
*STATIC
.1,1.,1.E-6
*BOUNDARY,OP=NEW
CENTER,1,1
CENTER,6,6
181,2,2
*END STEP