*HEADING 
 DEEP DRAWING OF A SQUARE BOX 
 (1/4 SYMMETRY MODEL)
**
**  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.
*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
*ELEMENT,TYPE=MASS,ELSET=EMASS
 20000,299991
*MASS,ELSET=EMASS
 6.396E-1,
**
**  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
*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
**
** Define surfaces
**
*SURFACE, NAME=PUNCH
PUNCH,SNEG
*SURFACE, NAME=HOLDER
HOLDER,SPOS
*SURFACE, NAME=DIE
DIE,SPOS
*SURFACE, NAME=BLANK
BLANK,
**
**  Apply symmetry boundary conditions
**
*BOUNDARY
D1,2,2
P1,2,2
H1,2,2
D2,1,1
P2,1,1
H2,1,1
B1,YSYMM
B4,XSYMM
199991, 1,2
199991, 4,6
299991,1,2
299991,4,6
599991,1,3
599991,4,6
**
**  Simulate the deep drawing operation
**
*RESTART,WRITE,NUMBER INTERVAL=2, TIMEMARKS=NO
**
*STEP
*DYNAMIC,EXPLICIT
,.0036
*AMPLITUDE,DEFINITION=SMOOTH STEP,NAME=SMOOTH
0.,0., 0.0036, 1.
*BOUNDARY,AMPLITUDE=SMOOTH
199991,3,3,-0.036
*CLOAD
 299991, 3, -2.287E4
**
*SURFACE INTERACTION,NAME=PUNCH_TOP
*FRICTION
 0.25,
*CONTACT PAIR,INTERACTION=PUNCH_TOP, 
MECHANICAL CONSTRAINT=PENALTY
PUNCH,BLANK
*CONTACT PAIR, MECHANICAL CONSTRAINT=PENALTY
HOLDER,BLANK
*SURFACE INTERACTION,NAME=DIE_BOTT
*FRICTION
 0.125,
*CONTACT PAIR,INTERACTION=DIE_BOTT,
MECHANICAL CONSTRAINT=PENALTY
DIE,BLANK
**
** Output requests for postprocessing
**
*NSET,NSET=NPUNCH
 199991,
*NSET,NSET=HOLDER
 299991,
*ELSET,ELSET=ELHIST
 3541,
*OUTPUT,FIELD,VAR=PRESELECT
*ELEMENT OUTPUT
STH,PEEQ
*OUTPUT,HISTORY,TIME INTERVAL=1.8E-4
*NODE OUTPUT,NSET=NPUNCH
U3,RF3
*NODE OUTPUT,NSET=HOLDER
U3,
*ELEMENT OUTPUT,ELSET=ELHIST
STH,
PEEQ,
**
** Output requests for qa testing
**
*ELSET,ELSET=QA_TESTEL,GEN
3701,3735,1
*NSET,NSET=QA_TESTN,ELSET=QA_TESTEL
*NSET,NSET=QA_TESTRB
199991,299991,599991
*ELSET,ELSET=QA_TESTELH
3647,
*OUTPUT,FIELD,NUM=1
*ELEMENT OUTPUT,ELSET=QA_TESTEL
 STH, PEEQ
*NODE OUTPUT,NSET=QA_TESTN
 U, 
*OUTPUT,FIELD,NUM=4,TIMEMARKS=YES
*NODE OUTPUT,NSET=QA_TESTRB
U,RF
*ELEMENT OUTPUT,ELSET=QA_TESTELH
STH,
PEEQ,
*OUTPUT,HISTORY,TIME INTERVAL=3.E-4
*ENERGY OUTPUT,VAR=PRESELECT
*END STEP