*HEADING
 SHEET METAL FORMING WITH ANNEALING
 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,NSET=DIE1
1001,  0.,-0.05
1002,.090,-0.05
1011,.100,-.040
199991, 0., -0.05
*NGEN,NSET=DIE1,LINE=C
1002,1011,1,,.090,-.040,0.
*ELEMENT,TYPE=RAX2,ELSET=DIE1
1001, 1001,1002
*ELGEN,ELSET=DIE1
1001, 10,1,1
*NODE,NSET=DIE2
2001,.100,-.060
2002,.100,-.010
2011,.110, 0.
2012,.160, 0.
299991, .100, -.060
*NGEN,NSET=DIE2,LINE=C
2002,2011,1,,.110,-.010,0.,0.,0.,-1.
*ELEMENT,TYPE=RAX2,ELSET=DIE2
2001, 2001,2002
*ELGEN,ELSET=DIE2
2001, 11,1,1
*NODE,NSET=DIE3
**raised by 0.05 from original 
**shift outer surface inward by 
**then lower by 0.005 (half the radius of curvature)
**further lower by 0.005
3001,.0  ,0.0044
3014,.050,-0.009
3015,.0888,-0.009
3024,.0988,0.001
3025,.0988,0.0405
399991,.0,0.0044
*NGEN,NSET=DIE3,LINE=C
3001,3014,1,, 0.0,-0.0956,0.,0.,0.,-1. 
3015,3024,1,,.0888, 0.001, 0.
*ELEMENT,TYPE=RAX2,ELSET=DIE3
3001, 3001,3002
*ELGEN,ELSET=DIE3
3001, 24,1,1
*BOUNDARY
1,XSYMM
199991,1,6
299991,1,6
399991,1,2
399991,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=DIE1
DIE1,SPOS
*SURFACE,TYPE=ELEMENT,NAME=DIE2
DIE2,SPOS
*SURFACE,TYPE=ELEMENT,NAME=DIE3
DIE3,SNEG
*SURFACE,TYPE=ELEMENT,NAME=BOTTOM
BLANK,SNEG
*RIGID BODY,ELSET=DIE3,REF NODE=399991
*RIGID BODY,ELSET=DIE1,REF NODE=199991
*RIGID BODY,ELSET=DIE2,REF NODE=299991
*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
*BULK VISCOSITY
0.02, 0.4
*DLOAD,AMP=R1
BLANK,P,-1.0
*SURFACE INTERACTION,NAME=FRICT
*FRICTION
 0.1,
*CONTACT PAIR,INTERACTION=FRICT,MECHANICAL CONSTRAINT=PENALTY
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
PEEQ,STH
*NODE OUTPUT, NSET=QA_TEST
U,
*OUTPUT,HISTORY,VAR=PRESELECT,TIME INTERVAL=1E-3
*OUTPUT, HISTORY
*ENERGY OUTPUT, VAR=PRESELECT
**************
***OUTPUT,HISTORY,TIME INTERVAL=3.5E-6
*NSET,NSET=NOUT
 1,
***NODE OUTPUT,NSET=NOUT
**U,V
*ELSET,ELSET=EOUT
 24,25,26
***ELEMENT OUTPUT,ELSET=EOUT
**STH,
**PEEQ,MISES,S,LE,PE
***ENERGY OUTPUT
**ALLKE,ALLSE,ALLWK,ALLPD,ALLAE,ALLCD,ALLFD,ALLIE,ALLVD,ETOTAL,DT
*END STEP
**
** First springback 
**
*STEP
*DYNAMIC,EXPLICIT
,1.E-3
*BULK VISCOSITY
0.02, 0.4
*DLOAD,AMP=R2,OP=NEW
BLANK,P,-1.0
*END STEP
**
** anneal 
**
*STEP
*ANNEAL
*END STEP
**
** gradually slide die 3 into position 
**
*STEP
*DYNAMIC,EXPLICIT
,1.0E-3
*BULK VISCOSITY
0.02, 0.4
*BOUNDARY, AMP=R3A  
399991,2,2,-0.04
***DLOAD,OP=NEW
**BLANK, VP, 100.0
*CONTACT PAIR,INTERACTION=FRICT,MECHANICAL CONSTRAINT=PENALTY
DIE3,TOP
*CONTACT PAIR,OP=DELETE
DIE1,BOTTOM
*END STEP
**
** Second upward pressure loading
**
*STEP
*DYNAMIC,EXPLICIT
,1.0E-3
*BULK VISCOSITY
0.02, 0.4
*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
*BULK VISCOSITY
0.02, 0.4
*DLOAD,AMP=R5,OP=NEW
LOAD,P,1.0
*END STEP