*HEADING
 DEEP DRAWING OF CYLINDRICAL CUP WITH MGAX1
*RESTART,WRITE,FREQUENCY=25
*NODE
*NGEN,NSET=BOT
*NSET,NSET=WRKPC
 BOT,
*NODE,NSET=DIE
*NODE,NSET=PUNCH
*NODE,NSET=HOLDER
*NSET,NSET=TOOLS
 PUNCH,DIE,HOLDER
*NSET,NSET=CENTER
*ELEMENT,TYPE=MGAX1,ELSET=BLANK
*ELGEN,ELSET=BLANK
*MEMBRANE SECTION,MATERIAL=STEEL,ELSET=BLANK
 0.00082,
*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
*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
*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
*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
*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
*CONTACT PAIR,INTERACTION=FRIC4
 PCH_PNT,PUNCH
*SURFACE INTERACTION,NAME=FRIC4
*FRICTION
 0.0,
**
*INITIAL CONDITIONS,TYPE=STRESS
BLANK,1000.0,1000.0
*NSET,NSET=NALL,GENERATE
**
*STEP,UNSYMM=YES,NLGEOM
 STEP 1: PUSH THE BLANKHOLDER DOWN BY A PRESCRIBED DISPLACEMENT
*STATIC
 1.,1.
*BOUNDARY
 NALL,5,5
 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
 U,RF
*EL PRINT,ELSET=BLANK,FREQUENCY=500
 S,E
*CONTACT PRINT,FREQUENCY=500
*CONTACT FILE,FREQUENCY=500
*NODE FILE,NSET=TOOLS,FREQUENCY=10
 U,RF
*END STEP
**
**
*STEP,NLGEOM, UNSYMM=YES
 STEP 2: APPLY THE PRESCRIBED FORCE ON THE BLANKHOLDER AND RELEASE THE DISPLACEMENT
*STATIC
 1.,1.
*BOUNDARY,OP=NEW
 CENTER,1,1
 CENTER,5,5
 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.
**
*NODE PRINT,NSET=TOOLS,FREQUENCY=0
*NODE FILE,NSET=TOOLS,FREQUENCY=0
**
*END STEP
**
**
*STEP,INC=500,NLGEOM,UNSYMM=YES
 STEP 3: MOVE THE PUNCH DOWN
*STATIC
 0.01,1.,1.E-6
*BOUNDARY,OP=NEW
 CENTER,1,1
 CENTER,5,5
 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.
*CONTACT CONTROLS, STABILIZE
*NODE PRINT,NSET=TOOLS,FREQUENCY=1
 U,RF
*NODE FILE,NSET=TOOLS,FREQUENCY=1
 U,RF
*OUTPUT,FIELD,FREQ=0
*END STEP
**
**
*STEP,INC=1,NLGEOM, UNSYMM=YES
 STEP 4: FIX ALL NODES AND REMOVE THE IRS ELEMENTS
*STATIC
 1.,1.,1.,1.
*BOUNDARY,FIXED
 WRKPC,1
 WRKPC,2
 WRKPC,5
*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
 STEP 5: REPLACE THE BOUNDARY CONDITIONS BY THE REGULAR SET
*STATIC
 0.1,1.,1.E-6
*BOUNDARY,OP=NEW
 CENTER,1,1
 CENTER,5,5
 181,2,2
**
*END STEP