*HEADING
WANG AND BUDIANSKY'S SPHERICAL PUNCH WITH SAX1, 50 ELEMENTS
*PREPRINT,ECHO=YES
*PART,NAME=BLANK
*NODE
1,0.0,0.0
401,.05059,0.
501,.05918,0.
*NGEN
  1,401,10
401,501,10
*ELEMENT,TYPE=SAX1,ELSET=BLANK
1,1,11
*ELGEN,ELSET=BLANK
1,50,10
*ELSET,ELSET=CENTER,GEN
1,10,1
*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
250.0E6,      2.8581845E-02
260.0E6,      3.5517555E-02
270.0E6,      4.3525275E-02
280.0E6,      5.2722659E-02
290.0E6,      6.3235357E-02
300.0E6,      7.5197279E-02
310.0E6,      8.8750519E-02
320.0E6,      0.1040458
330.0E6,      0.1212430
340.0E6,      0.1405106
350.0E6,      0.1620263
360.0E6,      0.1859779
370.0E6,      0.2125620
380.0E6,      0.2419857
390.0E6,      0.2744660
400.0E6,      0.3102303
410.0E6,      0.3495160
420.0E6,      0.3925720
430.0E6,      0.4396578
440.0E6,      0.4910434
450.0E6,      0.5470111
460.0E6,      0.6078544
470.0E6,      0.6738777
480.0E6,      0.7453985
490.0E6,      0.8227461
500.0E6,      0.9062610
510.0E6,      0.9962980
**
*PART,NAME=PUNCH
*NODE,NSET=PUNCH
1000,0.,.051225
1001,0.,0.000425
1031,.0508,.051225
*NGEN,NSET=PUNCH,LINE=C
1001,1031,1,1000
*ELEMENT,TYPE=RAX2,ELSET=PUNCH
1001, 1001,1002
*ELGEN,ELSET=PUNCH
1001, 30,1,1
*NSET,NSET=REF_NODE
1000,
*RIGID BODY,ELSET=PUNCH,REF NODE=REF_NODE
*END PART
*PART,NAME=DIE
*NODE,NSET=DIE
2000, .05918,-.006775
2001, .05283,-.030425
2002, .05283,-.006775
2022, .05918,-0.000425
2023, .05930,-0.000425
*NGEN,NSET=DIE,LINE=C
2002,2022,1,2000,,,,0,0,-1.
*ELEMENT,TYPE=RAX2,ELSET=DIE
2001,2001,2002
*ELGEN,ELSET=DIE
2001,22,1,1
*NSET,NSET=REF_NODE
2000,
*RIGID BODY,ELSET=DIE,REF NODE=REF_NODE
*END PART
*ASSEMBLY,NAME=ASSEMBLY-1
*INSTANCE,NAME=BLANK-1,PART=BLANK
*NSET,NSET=NOUT
 1,
*ELSET,ELSET=EOUT
 22,23,24,25
*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,

PE,
*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