*HEADING
 2D Draw Bending problem for Numisheet'93 - Explicit analysis
 Punch velocity 30m/s
**
**   BLANK
*NODE
    1, 0.000, 0.000
  351, 0.175, 0.000
 1001, 0.000, 0.005
 1351, 0.175, 0.005
*NGEN,NSET=BLANK
    1, 351, 1
 1001,1351, 1
*NSET,NSET=CENTER
 1, 1001
*ELEMENT,TYPE=S4R
 1, 1, 3, 1003, 1001
*ELGEN,ELSET=BLANK
 1, 175, 2, 2
*SHELL SECTION,ELSET=BLANK,MATERIAL=STEEL
 0.00078, 5
*MATERIAL,NAME=STEEL
*ELASTIC
 206.E9, 0.3
*PLASTIC
 0.16700E+09, 0.00000E+00
 0.19410E+09, 0.10000E-01
 0.21913E+09, 0.20000E-01
 0.23803E+09, 0.30000E-01
 0.25348E+09, 0.40000E-01
 0.26668E+09, 0.50000E-01
 0.27826E+09, 0.60000E-01
 0.28864E+09, 0.70000E-01
 0.29807E+09, 0.80000E-01
 0.30674E+09, 0.90000E-01
 0.31477E+09, 0.10000E+00
 0.32226E+09, 0.11000E+00
 0.32930E+09, 0.12000E+00
 0.33594E+09, 0.13000E+00
 0.34224E+09, 0.14000E+00
 0.34822E+09, 0.15000E+00
 0.35393E+09, 0.16000E+00
 0.35940E+09, 0.17000E+00
 0.36464E+09, 0.18000E+00
 0.36968E+09, 0.19000E+00
 0.37454E+09, 0.20000E+00
 0.37922E+09, 0.21000E+00
 0.38375E+09, 0.22000E+00
 0.38814E+09, 0.23000E+00
 0.39239E+09, 0.24000E+00
 0.39652E+09, 0.25000E+00
 0.40053E+09, 0.26000E+00
 0.40443E+09, 0.27000E+00
 0.40822E+09, 0.28000E+00
 0.41193E+09, 0.29000E+00
 0.41554E+09, 0.30000E+00
 0.41906E+09, 0.31000E+00
 0.42251E+09, 0.32000E+00
 0.42588E+09, 0.33000E+00
 0.42917E+09, 0.34000E+00
 0.43240E+09, 0.35000E+00
 0.43556E+09, 0.36000E+00
 0.43866E+09, 0.37000E+00
 0.44169E+09, 0.38000E+00
 0.44468E+09, 0.39000E+00
 0.44760E+09, 0.40000E+00
 0.45047E+09, 0.41000E+00
 0.45330E+09, 0.42000E+00
 0.45607E+09, 0.43000E+00
 0.45880E+09, 0.44000E+00
 0.46148E+09, 0.45000E+00
 0.46413E+09, 0.46000E+00
 0.46672E+09, 0.47000E+00
 0.46928E+09, 0.48000E+00
 0.47181E+09, 0.49000E+00
 0.47429E+09, 0.50000E+00
*POTENTIAL
 1., 1.0402, 1.24897, 1.07895, 1., 1.
*DENSITY
 7800., 
**
**   SURFACES
*SURFACE,NAME=TOP
 BLANK, SPOS
*SURFACE,NAME=BOT
 BLANK, SNEG
**
**   DIE
*NODE
 2000, 0.026, 0.0025, -0.080
*SURFACE,TYPE=CYLINDER,NAME=DIE
 0.0, 0.0, 0.0, 1.0, 0.0, 0.0
 0.0, 1.0, 0.0
 START, 0.026, -0.08000
 LINE,  0.026, -0.00539
 CIRCL, 0.031, -0.00039, 0.031, -0.00539 
 LINE,  0.180, -0.00039
*RIGID BODY,REF NODE=2000,ANALYTICAL SURFACE=DIE
**
**   PUNCH
*NODE
 3000, 0.000, 0.0025, 0.001
*SURFACE,TYPE=CYLINDER,NAME=PUNCH
 0.0, 0.0, 0.0, 1.0, 0.0, 0.0
 0.0, 1.0, 0.0
 START, 0.025, 0.081
 LINE,  0.025, 0.006
 CIRCL, 0.020, 0.001, 0.020, 0.006
 LINE,  0.000, 0.001
*RIGID BODY,REF NODE=3000,ANALYTICAL SURFACE=PUNCH
**
**   HOLDER
*NODE
 4000, 0.0585, 0.0025, 0.010
*SURFACE,TYPE=CYLINDER,NAME=HOLDER
 0.0, 0.0, 0.0, 1.0, 0.0, 0.0
 0.0, 1.0, 0.0
 START, 0.086, 0.01000
 LINE,  0.086, 0.00539
 CIRCL, 0.081, 0.00039, 0.081, 0.00539
 LINE,  0.031, 0.00039
 LINE,  0.031, 0.01000
*RIGID BODY,REF NODE=4000,ANALYTICAL SURFACE=HOLDER
**
*ELEMENT,TYPE=MASS,ELSET=EMASS
 1000, 4000
*MASS,ELSET=EMASS
 2.5e-2, 
**
**   BOUNDARY CONDITIONS
*BOUNDARY
 CENTER, 1, 1
 CENTER, 5, 5
 BLANK,  2, 2
 BLANK,  4, 4
 BLANK,  6, 6
 2000,   1, 6
 3000,   1, 6
 4000,   1, 2
 4000,   4, 6
**
*AMPLITUDE,NAME=APUNCH, DEFINITION=SMOOTH STEP
0.,0., 0.00233333333, 1., 0.004666666667, 0.
*AMPLITUDE,NAME=SMOOTH, DEFINITION=SMOOTH STEP
 0.0, 0.0, 0.001, 1.0
*ELSET,ELSET=QA_TEST
 57, 151
*NSET,NSET=QA_TEST
 1, 20, 57, 2000, 3000, 4000
**
**   APPLY BLANKHOLDER FORCE
**
*STEP
*DYNAMIC,EXPLICIT
 ,0.001
*CLOAD,AMPLITUDE=SMOOTH
 4000, 3, -175.0
**
*CONTACT PAIR,INTERACTION=CONTACT
 PUNCH, TOP
 HOLDER, TOP
 DIE, BOT
*SURFACE INTERACTION,NAME=CONTACT
*FRICTION
 0.144, 
*RESTART,WRITE,NUMBER INTERVAL=1
**
*FILE OUTPUT,NUMBER INTERVAL=1
*EL FILE,ELSET=QA_TEST
 PEEQ,
 STH,
*NODE FILE, NSET=QA_TEST
 U,
*OUTPUT,FIELD,NUMBER INTERVAL=1
*ELEMENT OUTPUT,ELSET=QA_TEST
 PEEQ, STH
*NODE OUTPUT,NSET=QA_TEST
 U,
*OUTPUT,HISTORY,FREQ=1000,VARIABLE=PRESELECT
*OUTPUT,FIELD,NUMBER INTERVAL=2
*ELEMENT OUTPUT
 PEEQ, STH
*NODE OUTPUT
 U,
*END STEP
**
**   MOVE PUNCH DOWN
**
*STEP
*DYNAMIC,EXPLICIT
 ,0.004666666667
*BOUNDARY,AMPLITUDE=APUNCH,TYPE=VELOCITY
 3000,3,3,-30.0
*RESTART,WRITE,NUMBER INTERVAL=1
*END STEP