*HEADING
 SHEET METAL FORMING WITH ANNEALING
 WITH AN ANALYTICAL RIGID SURFACE
 DIE3 SLIDES DOWN FROM ABOVE TO AVOID INITIAL OVERCLOSURE
*NODE,NSET=BLANK
1,0.,0.0005
41,.150,0.0005
*NGEN,NSET=BLANK
1,41,1
*ELEMENT,TYPE=SAX1,ELSET=BLANK
1, 1,2
*ELGEN,ELSET=BLANK
1, 40,1,1
*SHELL SECTION,ELSET=BLANK,MATERIAL=STEEL, SECTION INTEGR=GAUSS
.001,5
*MATERIAL,NAME=STEEL
*DENSITY
7800., 
*ELASTIC
2.1E11,0.3
*PLASTIC
  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
**
*NODE
199991, 0., -0.05
299991, .100, -.060
399991,.0,0.0044
399992,.0,0.0044
**
*NSET,NSET=OUTN,GENERATE
1,41,1
*NSET,NSET=OUTN
199991,299991,399991
**
*BOUNDARY
1,XSYMM
199991,1,6
299991,1,6
399991,1,2
399991,6,6
399992,1,2
399992,6,6
*AMPLITUDE,NAME=R1,TIME=STEP TIME
0.,0.,  .8E-3,1.5E6, 1.7E-3,1.5E6, 3.E-3,1.E7,
3.5E-3,1.E7
*AMPLITUDE,NAME=R2,TIME=STEP TIME
0.,1.E7,   1.E-3,0.
*AMPLITUDE,NAME=R3A,TIME=STEP TIME,DEFINITION=SMOOTH STEP 
0.,0.,   1.E-3, 1.0
*AMPLITUDE,NAME=R4,TIME=STEP TIME
0.,0.,   1.E-3,.6E7
*AMPLITUDE,NAME=R5,TIME=STEP TIME
0.,.6E7,   1.E-3,0.
*RESTART,WRITE,NUM=1
**
** First downward pressure loading
**
*SURFACE,TYPE=ELEMENT,NAME=TOP
BLANK,SPOS
*SURFACE,TYPE=ELEMENT,NAME=TOP1
BLANK,SPOS
*SURFACE, NAME=DIE1, TYPE=SEGMENTS
START,  0.,-0.05
LINE, .090,-0.05
CIRCL, .100,-.040, .090,-.040
*SURFACE, NAME=DIE2, TYPE=SEGMENTS
START, .100,-.060
LINE, .100,-.010
CIRCL, .110, 0., .110,-.010
LINE, .160, 0.
*SURFACE,TYPE=ELEMENT,NAME=BOTTOM
BLANK,SNEG
*SURFACE, NAME=DIE3, TYPE=SEGMENTS
START,  .0988, .0405
LINE, .0988,0.001
CIRCL, .0888,-.009, .0888,0.001
LINE, .050,-.009
CIRCL,  .0, 0.0044, .0,-0.0956
*SURFACE, NAME=DIE4, TYPE=SEGMENTS
START,  .099, .0005
LINE, .099,-.039
CIRCL, .089,-.049, .089,-.039
LINE, .050,-.049
CIRCL,  .0, -.0356,.0,-.1356
*RIGID BODY, REF NODE=399992, ANALYTICAL SURFACE =DIE4
*RIGID BODY, REF NODE=399991, ANALYTICAL SURFACE =DIE3
*RIGID BODY, REF NODE=199991, ANALYTICAL SURFACE =DIE1
*RIGID BODY, REF NODE=299991, ANALYTICAL SURFACE =DIE2
*ELSET, ELSET=QA_TEST, GEN
5, 40, 5
*ELSET, ELSET=QA_TEST
1,
*NSET, NSET=QA_TEST, GEN
1, 41, 5
*STEP
*DYNAMIC,EXPLICIT
,3.5E-3
*DLOAD,AMP=R1
BLANK,P,-1.0
*SURFACE INTERACTION,NAME=FRICT
*FRICTION
 0.1,
*CONTACT PAIR,INTERACTION=FRICT
DIE2,BOTTOM
DIE1,BOTTOM
**
**************
*FILE OUTPUT,NUM=1, TIMEMARKS=YES
*EL FILE, ELSET=QA_TEST
STH,
PEEQ
*NODE FILE, NSET=QA_TEST
U,
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTERVAL=1
*ELEMENT OUTPUT, ELSET=QA_TEST
STH, PEEQ
*NODE OUTPUT, NSET=QA_TEST
U,
*OUTPUT, FIELD, NUMBER INTERVAL=10
*ELEMENT OUTPUT
PEEQ,S
*NODE OUTPUT
U,V,A
*OUTPUT, HISTORY
*ENERGY OUTPUT, VAR=PRESELECT
*END STEP
**
** First springback 
**
*STEP
*DYNAMIC,EXPLICIT
,1.E-3
*DLOAD,AMP=R2,OP=NEW
BLANK,P,-1.0
*END STEP
*******************************
**
** anneal 
**
*STEP
*ANNEAL
*END STEP
** 
** gradually slide die3 into position
**
*STEP
*DYNAMIC,EXPLICIT
,1.0E-3
*BOUNDARY, AMP=R3A  
399991,2,2,-0.04
***DLOAD,OP=NEW
**BLANK, VP, 100.0
*CONTACT PAIR,INTERACTION=FRICT
DIE3,TOP
*CONTACT PAIR,OP=DELETE
DIE1,BOTTOM
*END STEP
**
** Second upward pressure loading
**
*STEP
*DYNAMIC,EXPLICIT
,1.0E-3
*ELSET,ELSET=LOAD,GEN
1,26,1
*DLOAD,AMP=R4,OP=NEW
LOAD,P,1.0
*END STEP
**
** Second springback 
**
*STEP
*DYNAMIC,EXPLICIT
,1.E-3
*DLOAD,AMP=R5,OP=NEW
LOAD,P,1.0
*END STEP