*HEADING
 WANG AND BUDIANSKY'S SPHERICAL PUNCH WITH M3D9R  
*PREPRINT,ECHO=YES,model=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
*NCOPY,OLD SET=METND,NEW SET=METND,CHANGE NUMBER=1000,SHIFT,MULTIPLE=2
 0., 0., 0.,
 0.,-1.,0.,0.,1.,0.,5.
*NSET,NSET=METND2,GENERATE
11,501,10
1001,1501,10
2001,2501,10
*NSET,NSET=LOWER,GENERATE
 11,501,10
*NSET,NSET=UPPER,GENERATE
 2011,2511,10
*ELEMENT,TYPE=M3D9R,ELSET=BLANK
 1,21,2021,1,1,1021,2011,1,11,1011
 2, 21,41,2041,2021,31,1041,2031,1021,1031
*ELGEN,ELSET=BLANK
 2, 24,20,2
*ELSET,ELSET=MET_SURF,GENERATE
2,48,2
*MEMBRANE SECTION,ELSET=BLANK,MATERIAL=SAMP,POISSON=0.5,ORIENTATION=OR1, MEMBRANE THICKNESS=DISTTHICK
 1.,
*distribution, name=distThick, table=distTableThick, location=element 
  ,0.1
 blank, .00085
*ORIENTATION,NAME=OR1,SYSTEM=CYLINDRICAL
0.0,-1.0,0.0, 0.0,1.0,0.0
3,0.0
*SURFACE, NAME=TOP
MET_SURF,SPOS
*SURFACE, NAME=BOTTOM
MET_SURF,SNEG
*END PART
*distribution table, name=distTableThick
 length
*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=REVOLUTION
0.,0.,0.,0.,0.,0.
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=REVOLUTION
0.,0.,0.,0.,0.,0.
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
*TRANSFORM,NSET=METND2,TYPE=C
0.,0.,0.,0.,1.,0.
*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
*INITIAL CONDITIONS,TYPE=STRESS
ASSEMBLY-1.BLANK-1.BLANK,8.5E6,8.5E6
*BOUNDARY
ASSEMBLY-1.BLANK-1.1,1,1
ASSEMBLY-1.BLANK-1.1,3,3
**ASSEMBLY-1.BLANK-1.1,4,6
ASSEMBLY-1.BLANK-1.501,1,3
ASSEMBLY-1.BLANK-1.1501,1,3
ASSEMBLY-1.BLANK-1.2501,1,3
ASSEMBLY-1.BLANK-1.LOWER,YSYMM
ASSEMBLY-1.BLANK-1.UPPER,YSYMM
ASSEMBLY-1.PUNCH-1.1000,1,1
ASSEMBLY-1.PUNCH-1.1000,3,3
ASSEMBLY-1.PUNCH-1.1000,4,6
ASSEMBLY-1.DIE-1.2000,1,2
ASSEMBLY-1.DIE-1.2000,3,3
ASSEMBLY-1.DIE-1.2000,4,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
*NSET,NSET=QA_TEST_NSET1,GENERATE
ASSEMBLY-1.BLANK-1.21,ASSEMBLY-1.BLANK-1.111,1
*ELSET,ELSET=QA_TEST_ELSET1
ASSEMBLY-1.EALL
*NSET,NSET=QA_TEST_NSET2
ASSEMBLY-1.NALL
*EL PRINT,FREQUENCY=0
*NODE FILE,FREQUENCY=1000
U,RF
COORD, 
*CONTACT FILE,SECONDARY=ASSEMBLY-1.BLANK-1.TOP,FREQUENCY=1000
*CONTACT FILE,SECONDARY=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=1000
*CONTACT OUTPUT,VARIABLE=PRESELECT,SECONDARY=ASSEMBLY-1.BLANK-1.TOP
*CONTACT OUTPUT,VARIABLE=PRESELECT,SECONDARY=ASSEMBLY-1.BLANK-1.BOTTOM
*NODE OUTPUT,VARIABLE=PRESELECT
*ELEMENT OUTPUT
S,LE,NE
*OUTPUT,HISTORY,FREQUENCY=1000
*ENERGY OUTPUT,VARIABLE=PRESELECT
*OUTPUT,FIELD,FREQUENCY=1000
*CONTACT OUTPUT,VARIABLE=PRESELECT,SECONDARY=ASSEMBLY-1.BLANK-1.TOP,NSET=QA_TEST_NSET1
*CONTACT OUTPUT,VARIABLE=PRESELECT,SECONDARY=ASSEMBLY-1.BLANK-1.BOTTOM,NSET=QA_TEST_NSET1
*NODE OUTPUT,VARIABLE=PRESELECT,NSET=QA_TEST_NSET2
*ELEMENT OUTPUT,ELSET=QA_TEST_ELSET1
S,LE,NE
*OUTPUT,HISTORY,FREQUENCY=1000
*ENERGY OUTPUT,VARIABLE=PRESELECT,ELSET=QA_TEST_ELSET1
*END STEP
*STEP,INC=2000,NLGEOM,UNSYMM=YES
*STATIC
0.05,100.,1E-5
*EL FILE,ELSET=ASSEMBLY-1.BLANK-1.BLANK,FREQUENCY=1000 
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,3
**ASSEMBLY-1.BLANK-1.METND,4,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 
S, 
*END STEP
**
** STEP 5 IS NOT DONE IN THIS ANALYSIS BECAUSE AT THE END OF THE
** SPRINGBACK STEP THE SYSTEM BECOMES SINGULAR DUE TO THE MEMBRANES
** NOT HAVING ANY BENDING STIFFNESS.
**