*HEADING
WANG AND BUDIANSKY'S SPHERICAL PUNCH WITH SAX1
ELEMENTS
PUNCH AND DIE ARE ANALYTICAL RIGID SEGMENTS
*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, NAME=TOP
BLANK,SPOS
*SURFACE, 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=REF_NODE
1000,0.,.051225
*RIGID BODY, REF NODE=REF_NODE,
ANALYTICAL SURFACE=PUNCH_BOT
*SURFACE, NAME=PUNCH_BOT, TYPE=SEGMENTS
START, .0508,.051225
CIRCL, 0.,0.000425, 0.,.051225
*END PART
*PART,NAME=DIE
*NODE,NSET=REF_NODE
2000, .05918,-.006775
*RIGID BODY, REF NODE=REF_NODE,
ANALYTICAL SURFACE=DIE_TOP
*SURFACE, NAME=DIE_TOP, TYPE=SEGMENTS
START, .05283,-.030425
LINE, .05283,-.006775
CIRCL, .05918,-0.000425, .05918,-.006775
LINE, .05930,-0.000425
*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,
*END INSTANCE
*INSTANCE,NAME=DIE-1,PART=DIE
*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
*OUTPUT, FIELD, NUMBER INTERVAL=4, TIMEMARKS=NO
*CONTACT OUTPUT, CPSET=PUNCH_DIE, 
VARIABLE=PRESELECT
*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