*HEADING
 2D Draw Bending problem for Numisheet'93 - Explicit analysis
 Punch velocity 5m/s 
** Analysis using the general contact capability
**
*PREPRINT, CONTACT=YES, MODEL=YES, HISTORY=YES
**
**   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., 
**
**   DIE
*NODE
 2000, 0.026, 0.0025, -0.080
*NODE, NSET=DIE
 2001, 0.026, -0.001, -0.080
 2002, 0.026, -0.001, -0.00539
 2008, 0.031, -0.001, -0.00039
 2009, 0.180, -0.001, -0.00039
*NGEN, NSET=DIE, LINE=C
 2002, 2008, 1, , 0.031, -0.001, -0.00539
*NCOPY, CHANGE=100, OLD=DIE, NEW=DIE, SHIFT
 0.0, 0.007

*ELEMENT, TYPE=R3D4
 2001, 2001, 2002, 2102, 2101
*ELGEN, ELSET=DIE
 2001, 8, 1, 1
*RIGID BODY, ELSET=DIE, REF NODE=2000
**
**   PUNCH
*NODE
 3000, 0.000, 0.0025, 0.001
*NODE, NSET=PUNCH
 3001, 0.000, -0.001, 0.001
 3002, 0.020, -0.001, 0.001
 3008, 0.025, -0.001, 0.006
 3009, 0.025, -0.001, 0.081
*NGEN, NSET=PUNCH, LINE=C
 3002, 3008, 1, , 0.02, -0.001, 0.006
*NCOPY, CHANGE=100, OLD=PUNCH, NEW=PUNCH, SHIFT
0.0, 0.007

*ELEMENT, TYPE=R3D4
 3001, 3001, 3002, 3102, 3101
*ELGEN, ELSET=PUNCH
 3001, 8, 1, 1
*RIGID BODY, ELSET=PUNCH, REF NODE=3000
**
**   HOLDER
*NODE
 4000, 0.0585, 0.0025, 0.010
*NODE, NSET=HOLDER
 4001, 0.031, -0.001, 0.010
 4002, 0.031, -0.001, 0.00039
 4003, 0.081, -0.001, 0.00039
 4009, 0.086, -0.001, 0.00539
*NGEN, NSET=HOLDER, LINE=C
 4003, 4009, 1, , 0.081, -0.001, 0.00539
*NCOPY, OLD=HOLDER, NEW=HOLDER, CHANGE=100, SHIFT
 0, 0.007

*ELEMENT, TYPE=R3D4
 4001, 4001, 4002, 4102, 4101
*ELGEN, ELSET=HOLDER
 4001, 8, 1, 1
*RIGID BODY, ELSET=HOLDER, REF NODE=4000
**
*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.014, 1., 0.028, 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
**
** GENERAL CONTACT DEFINITION
**
*CONTACT
*CONTACT INCLUSIONS, ALL EXTERIOR
*CONTACT PROPERTY ASSIGNMENT
 , , CONTACT
*SURFACE INTERACTION,NAME=CONTACT
*FRICTION
 0.144, 
**
**   APPLY BLANKHOLDER FORCE
**
*STEP
*DYNAMIC,EXPLICIT
 ,0.001
*CLOAD,AMPLITUDE=SMOOTH
 4000, 3, -175.0
**
*RESTART,WRITE,NUMBER INTERVAL=1
**Output request for QA
*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,var=preselect,time interval=1e-3
**
*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.028
*BOUNDARY,AMPLITUDE=APUNCH,TYPE=VELOCITY
 3000,3,3,-5.0
*RESTART,WRITE,NUMBER INTERVAL=1
*END STEP