*HEADING
 SHEET METAL FORMING WITH ANNEALING
 WITH AN ANALYTICAL RIGID SURFACE
 DIE3 SLIDES DOWN FROM ABOVE TO AVOID INITIAL OVERCLOSURE
*NODE,NSET=BLANK
*NGEN,NSET=BLANK
*ELEMENT,TYPE=SAX1,ELSET=BLANK
*ELGEN,ELSET=BLANK
*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
**
*NODE
*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, 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
*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
*ELSET, ELSET=QA_TEST
*NSET, NSET=QA_TEST, GEN
*STEP
*DYNAMIC,EXPLICIT
,3.5E-3
*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,VAR=PRESELECT
*NODE OUTPUT, NSET=QA_TEST
U,
*OUTPUT, HISTORY, TIME INTERVAL=0.5E-3
*ENERGY OUTPUT
ALLAE,ALLPD,ALLCD,ALLFD
**************
***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
*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
*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
*ELSET,ELSET=LOAD,GEN
*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