*HEADING
 DEEP DRAWING OF CYLINDRICAL CUP WITH S4R/S3R 
*RESTART,WRITE,FREQUENCY=25
*PREPRINT,CONTACT=YES
** Define end nodes in x-z plane
*NODE,NSET=ENDS
101, 
181,9.951847E-2,-9.801714E-3,0.
** Generate end nodes rotated by theta=11.25 degrees
*NCOPY, OLD SET=ENDS, CHANGE NUMBER=1000, SHIFT, NEW SET=ENDS2
0.,0.,0.
0.,0.,0.,0.,0.,1.,11.25
**Generate nodes for blank
*NGEN,NSET=BOT
101,181,2
1101,1181,2
*NSET,NSET=PLANE1,GENERATE
103,181,2
*NSET,NSET=PLANE2,GENERATE
1101,1181,2
*NSET,NSET=PLANE
PLANE1,PLANE2
*NSET,NSET=WRKPC
BOT, 
**Define rigid body nodes for the die, punch and holder
*NODE,NSET=DIE
100,0.,0.,-0.05
*NODE,NSET=PUNCH
200,0.,0.,.05
*NODE,NSET=HOLDER
300,0.,0.,0.05
*NSET,NSET=TOOLS
PUNCH,DIE,HOLDER
*NSET,NSET=CENTER
101, 
**Put a triangle in the centre of the wedge
*ELEMENT,TYPE=S3R,ELSET=TRIA
201,101,121,1121
**Create quads for the rest of the blank
*ELEMENT,TYPE=S4R
202,121,123,1123,1121
*ELGEN,ELSET=QUADS
202,30,2,2
*ELSET,ELSET=BLANK
TRIA,QUADS
**Define contact surfaces between the blank and the holder
*ELSET,ELSET=HOLDCON,GENERATE
222,260,2
**Define interface elements between each node of the blank and the die
*ELSET,ELSET=DIECON,GENERATE
212,260,2
**Define interface elements between each node of the blank and the punch
*ELSET,ELSET=PUNCHCON,GENERATE
201,201,1
202,230,2
**Define the thickness and material properties for the blank
*SHELL SECTION,MATERIAL=STEEL,ELSET=BLANK,POISSON=0.5
.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
**Define rigid surfaces (master surfaces) for the die, holder and punch
*rigid body,analytical surface=RDIE,REFNODE=100
*surface,TYPE=REVOLUTION,NAME=RDIE,FILLET RADIUS=0.00541
0.,0.,0.,0.,0.,1.
START,0.05125,-0.06041
LINE,0.05125,-0.00041
LINE,0.12,-0.00041
*rigid body,analytical surface=RHOLDER,REFNODE=300
*surface,TYPE=REVOLUTION,NAME=RHOLDER,FILLET RADIUS=0.001
0.,0.,0.,0.,0.,1.
START,0.12,0.0004100175
LINE,0.05625,0.0004100175
LINE,0.05625,.06041
*rigid body,analytical surface=RPUNCH,REFNODE=200
*surface,TYPE=REVOLUTION,NAME=RPUNCH,FILLET RADIUS=.01341
0.,0.,0.,0.,0.,1.
START,0.05,0.06041
LINE,0.05,2.91782E-3
CIRCL,0.0,.00051,0.0,1.00051
** define slave surfaces
*surface,NAME=DIE
DIECON,S2
*surface,NAME=HOLDER
HOLDCON,S1
*surface,NAME=PUNCH
PUNCHCON,S1
** define contact pairs
*CONTACT PAIR,INTERACTION=DIESURF,TYPE=SURFACE
DIE,RDIE
*CONTACT PAIR,INTERACTION=HOLDSURF,TYPE=SURFACE
HOLDER,RHOLDER
*CONTACT PAIR,INTERACTION=PUNSURF,TYPE=SURFACE
PUNCH,RPUNCH
**Define surface interaction properties. Initally set friction to 0. This is changed
**later in the STEP definitions
*SURFACE INTERACTION,NAME=DIESURF
*FRICTION
0.1, 
*SURFACE INTERACTION,NAME=HOLDSURF
*FRICTION
0.1, 
*SURFACE INTERACTION,NAME=PUNSURF
*FRICTION
0.25, 
**Apply a transform to all nodes at theta=(+/-)5.625
*TRANSFORM,TYPE=C,NSET=PLANE
0.,0.,0.,0.,0.,1.
******************************************************
*STEP, NLGEOM, INC=10, UNSYMM=YES
 STEP 1: PUSH THE BLANKHOLDER DOWN BY A PRESCRIBED DISPLACEMENT
*STATIC
1.,1.
*BOUNDARY
CENTER,1,2
CENTER,4,6
WRKPC,2,2
WRKPC,4,4
WRKPC,6,6
DIE,1,6
PUNCH,1,6
HOLDER,1,2
HOLDER,3,3,-1.77E-8
HOLDER,4,6
*MONITOR,NODE=200,DOF=3
*PRINT,CONTACT=YES
*NODE PRINT,NSET=TOOLS,FREQUENCY=100
COORD,U
RF, 
*EL PRINT,ELSET=BLANK,FREQUENCY=500
S,E
STH, 
*CONTACT PRINT,SLAVE=HOLDER,MASTER=RHOLDER,FREQUENCY=500
*CONTACT PRINT,SLAVE=DIE,MASTER=RDIE,FREQUENCY=500
*CONTACT PRINT,SLAVE=PUNCH,MASTER=RPUNCH,FREQUENCY=500
*OUTPUT, FIELD, FREQ=10
*NODE OUTPUT,NSET=TOOLS
 U,RF
*END STEP
******************************************************
*STEP,INC=10,NLGEOM, UNSYMM=YES
 STEP 2: APPLY THE PRESCRIBED FORCE ON THE BLANKHOLDER AND RELEASE THE DISPLACEMENT
*STATIC
1.,1.
*BOUNDARY,OP=NEW
CENTER,1,2
CENTER,4,6
WRKPC,2,2
WRKPC,4,4
WRKPC,6,6
DIE,1,6
PUNCH,1,6
HOLDER,1,2
HOLDER,4,6
*CLOAD
HOLDER,3,-3125.
*END STEP
******************************************************
*STEP,INC=500,NLGEOM,UNSYMM=YES
 STEP 3: MOVE THE PUNCH DOWN
*STATIC
.01,1.,1.E-6
*controls,parameter=time
,,,,,,20
*contact controls, auto
*BOUNDARY,OP=NEW
CENTER,1,2
CENTER,4,6
WRKPC,2,2
WRKPC,4,4
WRKPC,6,6
DIE,1,6
PUNCH,1,2
PUNCH,3,3,-.06
PUNCH,4,6
HOLDER,1,2
HOLDER,4,6
*CLOAD
HOLDER,3,-3125.
*END STEP
******************************************************
*STEP,INC=1,NLGEOM, UNSYMM=YES
 STEP 4: FIX ALL NODES AND REMOVE THE CONTACT SURFACES
*STATIC
1.,1.,1.,1.
*BOUNDARY,FIXED
WRKPC,1,6
*MODEL CHANGE,TYPE=CONTACT PAIR,REMOVE
DIE,RDIE
HOLDER,RHOLDER
PUNCH,RPUNCH
*CLOAD,OP=NEW
HOLDER,3,0.
*END STEP
******************************************************
*STEP,INC=50,NLGEOM, UNSYMM=YES
 STEP 5: REPLACE THE BOUNDARY CONDITIONS BY THE REGULAR SET
*STATIC
.1,1.,1.E-6
*BOUNDARY,OP=NEW
CENTER,1,2
CENTER,4,6
WRKPC,2,2
WRKPC,4,4
WRKPC,6,6
181,3,3
*END STEP