*HEADING
 DEEP DRAWING OF A SQUARE BOX (1/4 SYMMETRY MODEL)
**
**  MESH FOR THE PUNCH
**
*PREPRINT, CONTACT=YES, MODEL=YES, HISTORY=YES
*RESTART,WRITE,NUMBER INTERVAL=2,TIMEMARKS=NO
***PREPRINT,ECHO=NO,MODEL=NO,HISTORY=NO
**
*NODE
*NGEN
*NGEN,LINE=C
*NSET,NSET=P1,GENERATE
*NSET,NSET=P2,GENERATE
*NFILL
 P1,P2,12,37
**
*NSET,NSET=P2C,GENERATE
*NCOPY,CHANGE NUMBER=25,OLD SET=P2C,SHIFT,MULTIPLE=12
*NSET,NSET=P3,GENERATE
*NSET,NSET=P4,GENERATE
*NFILL
 P3,P4,12,25
**
*NSET,NSET=PXZ
 P1,
*NSET,NSET=PYZ,GENERATE
**
**
*ELEMENT,TYPE=R3D4,ELSET=PUNCH
*ELGEN,ELSET=PUNCH
*RIGID BODY,ELSET=PUNCH,REF NODE=99901
**
**  MESH FOR THE HOLDER
**
*NODE
**
*NGEN,NSET=H1
*NGEN,NSET=H2
*NFILL
 H1,H2,12,17
*NCOPY,CHANGE NUMBER=17,OLD SET=H2,SHIFT,MULTIPLE=12
*NSET,NSET=H3,GENERATE
*NGEN,NSET=H4
*NFILL
 H3,H4,12,17
**
*ELEMENT,TYPE=R3D4
*ELGEN,ELSET=HOLDER
*RIGID BODY,ELSET=HOLDER,REF NODE=99902
**
*ELEMENT,TYPE=MASS,ELSET=EMASS
*MASS,ELSET=EMASS
 0.6396,
**
**  MESH FOR THE DIE
**
*NODE
*NGEN
*NGEN,LINE=C
*NSET,NSET=D1,GENERATE
*NSET,NSET=D2,GENERATE
*NFILL
 D1,D2,12,37
*NCOPY,CHANGE NUMBER=37,OLD SET=D2,SHIFT,MULTIPLE=12
*NSET,NSET=D3,GENERATE
*NSET,NSET=D4,GENERATE
*NFILL
 D3,D4,12,37
**
*ELEMENT,TYPE=R3D4
*ELGEN,ELSET=DIE
*RIGID BODY,ELSET=DIE,REF NODE=99903
**
**  MESH FOR THE BLANK (70 X 70)
**
*NODE
*NGEN,NSET=B1
*NGEN,NSET=B2
*NFILL
 B1,B2,70,71
*NSET,NSET=B4,GENERATE
**
*ELEMENT,TYPE=S4R
*ELGEN,ELSET=BLANK
*SHELL SECTION,MATERIAL=STEEL,ELSET=BLANK,SECTION INTEGRA=GAUSS,
CONTROLS=HGC
 .00082,5
*SECTION CONTROLS,NAME=HGC,HOURGLASS=ENHANCED
**
**  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
**
** Define surfaces
**
*surface, NAME=PUNCH
 PUNCH,SNEG
*surface, NAME=DIE
 DIE,SPOS
*surface, NAME=BLANK
 BLANK,
**
**  Apply symmetry boundary conditions
**
*BOUNDARY
 D1,YSYMM
 PXZ,YSYMM
 H1,YSYMM
 D4,XSYMM
 PYZ,XSYMM
 H4,XSYMM
**
 B1,YSYMM
 B4,XSYMM
**
 99901,1,2
 99901,4,6
 99902,1,2
 99902,4,6
 99903,1,6
**
*ELSET,ELSET=ELHIST
*ELSET, ELSET=QA_TEST
ELHIST,
*NSET, NSET=QA_TEST
**
**  Simulate the deep drawing operation
**
*************************************************************
**General contact definition
*CONTACT
*CONTACT INCLUSIONS, ALL ELEMENT BASED
*CONTACT PROPERTY ASSIGNMENT
 PUNCH, BLANK, PUNCH_T
 DIE,   BLANK, DIE_B
*SURFACE PROPERTY ASSIGNMENT, PROPERTY=THICKNESS
 BLANK,THINNING
*CONTACT FORMULATION, TYPE=PURE MASTER-SLAVE
 PUNCH, BLANK, MASTER
 DIE,   BLANK, MASTER
*SURFACE INTERACTION, NAME = PUNCH_T
*FRICTION
 0.25,
*SURFACE INTERACTION, NAME = DIE_B
*FRICTION
 0.125,
*************************************************************
*STEP
*DYNAMIC,EXPLICIT
 ,.0036
*AMPLITUDE,DEFINITION=SMOOTH STEP,NAME=SMOOTH
0.,0., 0.0036, 1.
*BOUNDARY, AMPLITUDE=SMOOTH
 99901,3,3,-0.036
*CLOAD
 99902,3,-2.287E4
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTER=1
*ELEMENT OUTPUT, ELSET=QA_TEST
PEEQ, STH
*NODE OUTPUT, NSET=QA_TEST
U,
*OUTPUT, HIST, FREQ=999
*ENERGY OUTPUT, VAR=PRESELECT
**
** output
**
*NSET,NSET=NPUNCH
*NSET,NSET=HOLDER
*OUTPUT,FIELD,VAR=PRESELECT, NUMBER INTERVAL=2
*ELEMENT OUTPUT
STH,PEEQ
*OUTPUT,HISTORY,TIME INTERVAL=1.8E-4
*ELEMENT OUTPUT,ELSET=ELHIST
STH,
PEEQ,
*NODE OUTPUT,NSET=NPUNCH
U3,RF3
*NODE OUTPUT,NSET=HOLDER
U3,
*ENERGY OUTPUT
ALLKE,ALLIE,ALLSE,ALLAE,ALLVD,ALLCD,ALLPD,ALLFD,ALLWK,ETOTAL
*END STEP