*HEADING 
 DEEP DRAWING OF A SQUARE BOX (1/4 SYMMETRY MODEL) AND THE SPRING BACK, BOTH 
 DONE IN /STANDARD.  FACETED RIGID SURFACE IS USED FOR THE RIGID BODIES. 
**
*PREPRINT,ECHO=NO
**  Generate the mesh for the PUNCH
**
**
*NODE
 101,.0 ,.0,0.00041
 107,.04,.0,0.00041
 113,.05,.0 ,.01
 119,.05,.0 ,.07
 199,.04,.0 ,.01
 701,.0 ,.04,0.00041
 707,.04,.04,0.00041
 713,.05,.04,.01
 719,.05,.04,.07
 799,.04,.04,.01
1307,.04,.04,.00041
1313,.04,.05,.01
1319,.04,.05,.07
1399,.04,.04,.01
1907,.0 ,.04,.00041
1913,.0 ,.05,.01
1919,.0 ,.05,.07
1999,.0 ,.04,.01
199991,0.,0.,0.00041
*NGEN,NSET=P1
101,107,1
113,119,1
*NGEN,NSET=P1,LINE=C
107,113,1,199
*NGEN,NSET=P2
701,707,1
713,719,1
*NGEN,NSET=P2,LINE=C
707,713,1,799
*NSET,NSET=P2C,GEN
707,719,1
*NGEN,NSET=P3
1313,1319,1
*NGEN,NSET=P3,LINE=C
1307,1313,1,1399
*NGEN,NSET=P4
1913,1919,1
*NGEN,NSET=P4,LINE=C
1907,1913,1,1999
*NFILL
P1,P2,6,100
P3,P4,6,100
*NCOPY,CHANGE NUMBER=100,OLD SET=P2C,SHIFT,MULTIPLE=6
0.,0.,0.
.04,.04,0., .04,.04,1.,15.
*NSET,NSET=P4,GEN
101,701,100
*ELEMENT,TYPE=R3D4,ELSET=PUNCH
  1, 101,102,202,201
109, 708,709,809,808
241, 707, 708, 808,707
242, 707, 808, 908,707
243, 707, 908,1008,707
244, 707,1008,1108,707
245, 707,1108,1208,707
246, 707,1208,1308,707
247, 707,1308,1408,706
248, 706,1408,1508,705
249, 705,1508,1608,704
250, 704,1608,1708,703
251, 703,1708,1808,702
252, 702,1808,1908,701
*ELGEN,ELSET=PUNCH
1, 18,1,1, 6,100,18
109, 11,1,1, 12,100,11
*RIGID BODY,ELSET=PUNCH,REF NODE=199991
**
**  Generate the mesh for the HOLDER
**
*NODE
20101,.06,.0,.01
20102,.06,.0,.00041
20108,.13,.0,.00041
20109,.13,.0,.01
20701,.06,.05,.01
20702,.06,.05,.00041
20708,.13,.05,.00041
20709,.13,.05,.01
21301,.05,.06,.01
21302,.05,.06,.00041
21308,.05,.13,.00041
21309,.05,.13,.01
21901,.0,.06,.01
21902,.0,.06,.00041
21908,.0,.13,.00041
21909,.0,.13,.01
299991,0.06,0.,0.01
*NGEN,NSET=H1
20101,20102,1
20102,20108,1
20108,20109,1
*NGEN,NSET=H2
20701,20702,1
20702,20708,1
20708,20709,1
*NGEN,NSET=H3
21301,21302,1
21302,21308,1
21308,21309,1
*NGEN,NSET=H4
21901,21902,1
21902,21908,1
21908,21909,1
*NFILL
H1,H2,6,100
H3,H4,6,100
*NCOPY,CHANGE NUMBER=100,OLD SET=H2,SHIFT,MULTIPLE=6
0.,0.,0.
.05,.05,0., .05,.05,1.,15.
*NSET, NSET=HLDREDG1, GEN 
20108, 21908, 100 
*NSET, NSET=HLDREDG2, GEN 
20102, 21902, 100 
*ELSET, ELSET=HLDREDG1, GEN
20007, 20143, 8 
*ELSET, ELSET=HLDREDG2, GEN
20002, 20138, 8
*NORMAL, TYPE=CONTACT SURFACE 
HOLDER, HLDREDG1, 0., 0., -1.
HOLDER, HLDREDG2, 0., 0., -1.
*ELEMENT,TYPE=R3D4,ELSET=HOLDER
20001, 20101,20201,20202,20102
*ELGEN,ELSET=HOLDER
20001, 8,1,1, 18,100,8
*RIGID BODY,ELSET=HOLDER,REF NODE=299991
**  Generate the mesh for the DIE
**
**
*NODE
599991, .05125, 0., -0.06
*NODE,NSET=DIE
50101, .05125,0.,-.06
50107, .05125,0.,-.01
50113, .06125,0.,-0.00041
50119, .13000,0.,-0.00041
50199, .06125,0.,-.01
**
50701, .05125,.04125,-.06
50707, .05125,.04125,-.01
50713, .06125,.04125,-0.00041
50719, .13000,.04125,-0.00041
50799, .06125,.04125,-.01
**
51301, .04125,.05125,-.06
51307, .04125,.05125,-.01
51313, .04125,.06125,-0.00041
51319, .04125,.13000,-0.00041
51399, .04125,.06125,-.01
**
51901, 0.,.05125,-.06
51907, 0.,.05125,-.01
51913, 0.,.06125,-0.00041
51919, 0.,.13000,-0.00041
51999, 0.,.06125,-.01
*NGEN,NSET=D1
50101,50107,1
50113,50119,1
*NGEN,NSET=D1,LINE=C
50107,50113,1,50199
*NGEN,NSET=D2
50701,50707,1
50713,50719,1
*NGEN,NSET=D2,LINE=C
50707,50713,1,50799
*NGEN,NSET=D3
51301,51307,1
51313,51319,1
*NGEN,NSET=D3,LINE=C
51307,51313,1,51399
*NGEN,NSET=D4
51901,51907,1
51913,51919,1
*NGEN,NSET=D4,LINE=C
51907,51913,1,51999
*NFILL
D1,D2,6,100
D3,D4,6,100
*NCOPY,CHANGE NUMBER=100,OLD SET=D2,SHIFT,MULTIPLE=6
0.,0.,0.
.04125,.04125,0., .04125,.04125,1., 15.
*ELEMENT,TYPE=R3D4
50001, 50101,50102,50202,50201
*ELGEN,ELSET=DIE
50001, 18,1,1, 18, 100,18
*RIGID BODY,ELSET=DIE,REF NODE=599991 
**
**
**  Generate the mesh for the BLANK
**
*NODE
90001,.0,.0,0.
90036,.1,.0,0.
91261,.0,.1,0.
91296,.1,.1,0.
*NGEN,NSET=B1
90001,90036,1
*NGEN,NSET=B2
91261,91296,1
*NFILL
B1,B2,35,36
*NSET,NSET=B4,GEN
90001,91261,36
*ELEMENT,TYPE=S4R,ELSET=BLANK
3001, 90001,90002,90038,90037
*ELGEN,ELSET=BLANK
3001, 35,1,1, 35,36,35
*SHELL SECTION,MATERIAL=STEEL,ELSET=BLANK
.00082,5
**
**  Define material properties for STEEL
**
*MATERIAL,NAME=STEEL
*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
**
**  Apply symmetry boundary conditions
**
*BOUNDARY
B1,YSYMM
B4,XSYMM
199991,1,6
299991,1,2
299991,4,6
599991,1,6
**
***************** contact definition
**
*surface, NAME=PUNCH
PUNCH,SNEG
*surface, NAME=HOLDER
HOLDER,SPOS
*surface, NAME=DIE
DIE,SPOS
*surface, NAME=TOP
BLANK,SPOS
*surface, NAME=BOTTOM
BLANK,SNEG
*CONTACT PAIR,TYPE=SURFACE TO SURFACE,INTERACTION=PUNCH_TOP
TOP,PUNCH
*CONTACT PAIR,TYPE=SURFACE TO SURFACE,INTERACTION=HLDR_TOP
TOP, HOLDER
*CONTACT PAIR,TYPE=SURFACE TO SURFACE,INTERACTION=DIE_BOTT 
BOTTOM, DIE
*SURFACE INTERACTION,NAME=HLDR_TOP
*surface behavior, penalty
,,0.01
*SURFACE INTERACTION,NAME=DIE_BOTT
*surface behavior, penalty
,,0.01
*FRICTION
 0.125,
*SURFACE INTERACTION,NAME=PUNCH_TOP
*surface behavior, penalty
,,0.01
*FRICTION
 0.25,
**
**  Simulate the deep drawing operation
**
**  grip the blank between the holder and the die in the next two steps.
**
*STEP, NLGEOM, CONVERT SDI=YES
*STATIC 
0.0000012,.0000012
*model change,remove,type=contact pair
 TOP,PUNCH
*RESTART,WRITE,FREQ=100
*PRINT, CONTACT=YES 
*BOUNDARY
299991,3,3,-0.000011
*NSET,NSET=NOUT
 B1, 199991, 299991, 599991 
*NSET,NSET=NPUNCH
 199991,
*NSET,NSET=HOLDER
 299991,
*EL PRINT, FREQ=0
*NODE PRINT, FREQ=0 
*ELSET,ELSET=ELHIST
 3541,
*OUTPUT,FIELD,FREQ=20
*ELEMENT OUTPUT,VARIABLE=PRESELECT
1,3,5
STH,
*NODE OUTPUT,VARIABLE=PRESELECT
*CONTACT OUTPUT,VARIABLE=PRESELECT
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=NPUNCH
U3,RF3
*ELEMENT OUTPUT,ELSET=ELHIST
STH,
PEEQ,
*energy output, var=all
*END STEP
**
*STEP
*STATIC 
,.0000012
*BOUNDARY,OP=NEW 
B1,YSYMM
B4,XSYMM
199991,1,6
299991,1,2
299991,4,6
599991,1,6
*CLOAD
 299991, 3, -2.287E4
**
*END STEP
**
**  move the punch down to deform the blank 
**
*STEP, NLGEOM, UNSYMM=YES, INC=10000
*STATIC
0.01, 1.0, 1.0e-10,
*MODEL CHANGE, ADD, TYPE=CONTACT PAIR
 TOP,PUNCH
*BOUNDARY
199991,3,3,-0.009
*contact controls,stabilize=0.01,tangent=0.0
*END STEP
**
**  move the punch down further to deform the blank 
**
*STEP, NLGEOM, UNSYMM=YES, INC=10000
*STATIC,stabilize,allsdtol
0.01, 1.0, 1.0e-10,
*BOUNDARY
199991,3,3,-0.036
**
*END STEP
**
**  spring back analysis 
** 
*STEP, NLGEOM
*STATIC 
0.1, 1.0
*MODEL CHANGE, REMOVE, TYPE=CONTACT PAIR
TOP,PUNCH
*CONTACT CONTROLS, STABILIZE=0.1
*END STEP
**
*STEP, NLGEOM
*STATIC 
0.1, 1.0
*BOUNDARY, FIXED
90001, 3, 3
299991, 3, 3
*MODEL CHANGE, REMOVE, TYPE=CONTACT PAIR
TOP, HOLDER
BOTTOM, DIE
*CLOAD, OP=NEW
*CONTACT CONTROLS, STABILIZE=0.1
*END STEP