*HEADING
WANG AND BUDIANSKY'S SPHERICAL PUNCH WITH SAX1, 50 ELEMENTS
*PREPRINT,ECHO=YES
*PART,NAME=BLANK
*NODE
*NGEN
*ELEMENT,TYPE=SAX1,ELSET=BLANK
*ELGEN,ELSET=BLANK
*ELSET,ELSET=CENTER,GEN
*SHELL SECTION,ELSET=BLANK,MATERIAL=SAMP
0.00085,
*SURFACE,TYPE=ELEMENT,NAME=TOP
BLANK,SPOS
*SURFACE,TYPE=ELEMENT,NAME=BOTTOM
BLANK,SNEG
*END PART
*MATERIAL,NAME=SAMP
*DENSITY
7850.,
*ELASTIC
206.8E9,0.3
*PLASTIC
170.0E6,      0.0000000E+00
180.0E6,      1.7205942E-03
190.0E6,      3.8296832E-03
200.0E6,      6.3897874E-03
210.0E6,      9.4694765E-03
220.0E6,      1.3143660E-02
230.0E6,      1.7493792E-02
240.0E6,      2.2608092E-02
**
*PART,NAME=PUNCH
*NODE,NSET=PUNCH
*NGEN,NSET=PUNCH,LINE=C
*ELEMENT,TYPE=RAX2,ELSET=PUNCH
*ELGEN,ELSET=PUNCH
*NSET,NSET=REF_NODE
*RIGID BODY,ELSET=PUNCH,REF NODE=REF_NODE
*END PART
*PART,NAME=DIE
*NODE,NSET=DIE
*NGEN,NSET=DIE,LINE=C
*ELEMENT,TYPE=RAX2,ELSET=DIE
*ELGEN,ELSET=DIE
*NSET,NSET=REF_NODE
*RIGID BODY,ELSET=DIE,REF NODE=REF_NODE
*END PART
*ASSEMBLY,NAME=ASSEMBLY-1
*INSTANCE,NAME=BLANK-1,PART=BLANK
*NSET,NSET=NOUT
*ELSET,ELSET=EOUT
*END INSTANCE
*INSTANCE,NAME=PUNCH-1,PART=PUNCH
***NSET,NSET=PUNCH
*** 1000,
*SURFACE,TYPE=ELEMENT,NAME=PUNCH_BOT
PUNCH,SNEG
*END INSTANCE
*INSTANCE,NAME=DIE-1,PART=DIE
*SURFACE,TYPE=ELEMENT,NAME=DIE_TOP
DIE,SPOS
*END INSTANCE
*END ASSEMBLY
*BOUNDARY
ASSEMBLY-1.BLANK-1.1,1,1
ASSEMBLY-1.BLANK-1.1,6,6
ASSEMBLY-1.BLANK-1.501,1,2
ASSEMBLY-1.PUNCH-1.1000,1,1
ASSEMBLY-1.PUNCH-1.1000,1,6
ASSEMBLY-1.DIE-1.2000,1,2
ASSEMBLY-1.DIE-1.2000,6,6
*AMPLITUDE,NAME=LOAD,TIME=TOTAL
0.,0.,1.24E-3,1.,1.57E-3,1.,1.97E-3,0.,
3.97E-3,-.25
*RESTART,WRITE,NUM=2,TIMEMARKS=NO
*STEP
*DYNAMIC,EXPLICIT
,1.24E-3
*BOUNDARY,TYPE=VELOCITY,AMP=LOAD
ASSEMBLY-1.PUNCH-1.1000,2,2,-30.
*SURFACE INTERACTION,NAME=PUNCH_TOP
*FRICTION
 0.275,
*CONTACT PAIR, INTERACTION=PUNCH_TOP,
CPSET=PUNCH_DIE
ASSEMBLY-1.PUNCH-1.PUNCH_BOT,ASSEMBLY-1.BLANK-1.TOP
ASSEMBLY-1.DIE-1.DIE_TOP,ASSEMBLY-1.BLANK-1.BOTTOM
*******
*NSET, NSET=QA_TEST
ASSEMBLY-1.BLANK-1.NOUT,
*ELSET, ELSET=QA_TEST
ASSEMBLY-1.BLANK-1.EOUT,
*FILE OUTPUT,NUMBER INTERVAL = 1
*EL FILE, ELSET=QA_TEST
PEEQ,
*NODE FILE, NSET=QA_TEST
U,
*ENERGY FILE
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTER=1
*ELEMENT OUTPUT, ELSET=QA_TEST
PEEQ,
*NODE OUTPUT, NSET=QA_TEST
U,
******
*OUTPUT, HIST, FREQ=999
*ENERGY OUTPUT, VAR=PRESELECT
*OUTPUT, FIELD, NUMBER INTER=4
*NODE OUTPUT
U,
*ELEMENT OUTPUT
PEEQ,
***OUTPUT,HISTORY,TIME INTERVAL=0.
***NODE OUTPUT,NSET=ASSEMBLY-1.PUNCH-1.PUNCH
**U2,RF2
***NODE OUTPUT,NSET=ASSEMBLY-1.BLANK-1.NOUT
**U,V
***ELEMENT OUTPUT,ELSET=ASSEMBLY-1.BLANK-1.EOUT
**STH,
**MISES,S,LE
***ENERGY OUTPUT
**ALLKE,ALLSE,ALLWK,ALLPD,ALLIE,ALLVD,ETOTAL,ALLAE,ALLCD,ALLFD,DT
*END STEP
*STEP
*DYNAMIC,EXPLICIT
,.33E-3
*END STEP
*STEP
*DYNAMIC,EXPLICIT
,.40E-3
*END STEP
*STEP
**
** Unloading step
**
*DYNAMIC,EXPLICIT
,2.E-3
*DLOAD,OP=NEW
ASSEMBLY-1.BLANK-1.BLANK,VP,0.35E6
*END STEP