*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,NSET=METND
  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,9
*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
0.,-0.000425,0.
*NSET,NSET=PUNCH
 1000,
*END INSTANCE
*INSTANCE,NAME=DIE-1,PART=DIE
0.,0.000425,0.
*END INSTANCE
*NSET,NSET=PUNK
PUNCH-1.1000,DIE-1.2000
*NSET,NSET=NALL
PUNK,BLANK-1.METND
*ELSET,ELSET=EALL
BLANK-1.BLANK
*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.BLANK-1.501,6,6
ASSEMBLY-1.PUNCH-1.1000,1,2
ASSEMBLY-1.PUNCH-1.1000,6,6
ASSEMBLY-1.DIE-1.2000,1,2
ASSEMBLY-1.DIE-1.2000,6,6
*SURFACE INTERACTION,NAME=PUNCH_TOP
*FRICTION
 0.275,
*CONTACT PAIR, INTERACTION=PUNCH_TOP
ASSEMBLY-1.BLANK-1.TOP,ASSEMBLY-1.PUNCH-1.PUNCH_BOT
ASSEMBLY-1.BLANK-1.BOTTOM,ASSEMBLY-1.DIE-1.DIE_TOP
*RESTART,WRITE,FREQUENCY=250
*STEP,INC=2000,NLGEOM,UNSYMM=YES
*STATIC
0.05,100.,1.E-5
*BOUNDARY
ASSEMBLY-1.PUNCH-1.1000,2,2,-0.0186
*EL PRINT,FREQUENCY=0
*NODE FILE,FREQUENCY=1000
U,RF
COORD, 
*CONTACT FILE,SLAVE=ASSEMBLY-1.BLANK-1.TOP,FREQUENCY=1000
*CONTACT FILE,SLAVE=ASSEMBLY-1.BLANK-1.BOTTOM,FREQUENCY=1000
*MONITOR,NODE=ASSEMBLY-1.PUNCH-1.1000,DOF=2
*NODE PRINT,NSET=ASSEMBLY-1.PUNK,FREQUENCY=100
U,RF
COORD, 
*OUTPUT,FIELD,FREQUENCY=4
*CONTACT OUTPUT,VARIABLE=PRESELECT
*NODE OUTPUT,VARIABLE=PRESELECT
*ELEMENT OUTPUT
1,5,9
S,LE,NE
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT
*END STEP
*STEP,INC=2000,NLGEOM,UNSYMM=YES
*STATIC
0.05,100.,1E-5
*EL FILE,ELSET=ASSEMBLY-1.BLANK-1.BLANK,FREQUENCY=1000
5, 
S,E
*BOUNDARY
ASSEMBLY-1.PUNCH-1.1000,2,2,-0.0285
*END STEP
*STEP,INC=2000,NLGEOM,UNSYMM=YES
*STATIC
0.05,100.,1.E-5
*BOUNDARY
ASSEMBLY-1.PUNCH-1.1000,2,2,-0.0345
*END STEP
*STEP,INC=2000,NLGEOM,UNSYMM=YES
*STATIC
100.,100.
*BOUNDARY,FIXED
ASSEMBLY-1.BLANK-1.METND,1,2
ASSEMBLY-1.BLANK-1.METND,6
*MODEL CHANGE,TYPE=CONTACT PAIR,REMOVE
ASSEMBLY-1.BLANK-1.TOP,ASSEMBLY-1.PUNCH-1.PUNCH_BOT
ASSEMBLY-1.BLANK-1.BOTTOM,ASSEMBLY-1.DIE-1.DIE_TOP
*EL FILE,ELSET=ASSEMBLY-1.BLANK-1.BLANK,FREQUENCY=1000
5, 
S, 
*END STEP
*STEP,INC=2000,NLGEOM, UNSYMM=YES
*STATIC
1.,100.
*BOUNDARY,OP=NEW
ASSEMBLY-1.BLANK-1.1,1,1
ASSEMBLY-1.BLANK-1.1,6,6
ASSEMBLY-1.BLANK-1.501,1,2
**ASSEMBLY-1.BLANK-1.501,6,6
ASSEMBLY-1.PUNCH-1.1000,1,2
ASSEMBLY-1.PUNCH-1.1000,6,6
ASSEMBLY-1.DIE-1.2000,1,2
ASSEMBLY-1.DIE-1.2000,6,6
*MONITOR,NODE=ASSEMBLY-1.BLANK-1.1,DOF=2
*EL FILE,ELSET=ASSEMBLY-1.BLANK-1.BLANK,FREQUENCY=1000
5, 
S, 
*END STEP