*HEADING
 2D Draw Bending problem for Numisheet'93
 Explicit analysis
*prePRINT, MODEL=YES
**
** The blank, made of mild steel, initially
** measures 350mm by 35mm and is 0.78mm
** thick. The total blank holder force is
** 2.45kN, and a mass of 5kg is attached to
** the blank holder.
**
** The problem is symmetric about a plane
** through the center of the punch, thus only
** half of the problem is modeled. The blank
** is modeled with a single row of 175
** first-order shell elements (S4R) with an
** element size of 1mm by 5mm. Therefore,
** the blank holder force used in this input
** file is 2.45kN * 5mm / 35mm = 175N, and
** the mass attached to the blank holder is
** 5.0kg / 2 = 2.5kg. The mass attached to
** the blank holder is reduced to 0.025kg in
** order to prevent a large mass mismatch
** between the contact surfaces.
**
**   BLANK
*NODE
 900000001, 0.000, 0.000
 900000351, 0.175, 0.000
 900001001, 0.000, 0.005
 900001351, 0.175, 0.005
*NGEN,NSET=BLANK
 900000001, 900000351, 1
 900001001, 900001351, 1
*NSET,NSET=CENTER
 900000001, 900001001
*ELEMENT,TYPE=S4R
 800000001, 900000001, 900000003, 900001003, 900001001
*ELGEN,ELSET=BLANK
 800000001, 175, 2, 2
*SHELL SECTION,ELSET=BLANK,MATERIAL=STEEL,CONTROLS=HGC
 0.00078, 5
*SECTION CONTROLS,NAME=HGC,HOURGLASS=ENHANCED
*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
 900002000, 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=900002000,ANALYTICAL SURFACE=DIE
**
**   PUNCH
*NODE
 900003000, 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=900003000,ANALYTICAL SURFACE=PUNCH
**
**   HOLDER
*NODE
 900004000, 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=900004000,ANALYTICAL SURFACE=HOLDER
**
*ELEMENT,TYPE=MASS,ELSET=EMASS
 1000, 900004000
*MASS,ELSET=EMASS
 2.5e-2, 
**
**   BOUNDARY CONDITIONS
*BOUNDARY
 CENTER, 1, 1
 CENTER, 5, 5
 BLANK,  2, 2
 BLANK,  4, 4
 BLANK,  6, 6
 900002000,   1, 6
 900003000,   1, 6
 900004000,   1, 2
 900004000,   4, 6
**
*AMPLITUDE,NAME=APUNCH, DEFINITION=SMOOTH STEP
 0.0, 0.0, 0.00466666667, 1.0, 0.00933333333, 0.0
*AMPLITUDE,NAME=SMOOTH, DEFINITION=SMOOTH STEP
 0.0, 0.0, 0.001, 1.0
*ELSET,ELSET=QA_TEST
 800000057, 800000151
*NSET,NSET=QA_TEST
 900000001, 900000020, 900000057, 
 900002000, 900003000, 900004000
*SURFACE INTERACTION,NAME=CONTACT
*FRICTION
 0.144, 
**
**   APPLY BLANKHOLDER FORCE
**
*STEP
*DYNAMIC,EXPLICIT
 ,0.001
*CLOAD,AMPLITUDE=SMOOTH
 900004000, 3, -175.0
**
*CONTACT
*CONTACT INCLUSIONS
,
*CONTACT PROPERTY ASSIGNMENT
,,CONTACT
*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.00933333333
*BOUNDARY,AMPLITUDE=APUNCH,TYPE=VELOCITY
 900003000,3,3,-15.0
*RESTART,WRITE,NUMBER INTERVAL=1
*END STEP