*HEADING
 WANG AND BUDIANSKY'S SPHERICAL PUNCH WITH SC8R  
*PREPRINT,ECHO=YES
*PART,NAME=BLANK
*NODE
     1,.00000,-0.000425
   401,.05059,-0.000425
   501,.05918,-0.000425
 50001,.00000,+0.000425
 50401,.05059,+0.000425
 50501,.05918,+0.000425
*NGEN,NSET=METND
       1,401,10
     401,501,10
 50001,50401,10
 50401,50501,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
  50011,50501,10
  51001,51501,10
  52001,52501,10
*NSET,NSET=LOWER,GENERATE
 11,501,10
 50011,50501,10
*NSET,NSET=UPPER,GENERATE
 2011,2511,10
 52011,52511,10
*ELEMENT,TYPE=SC6R,ELSET=BLANK
    1,  11,1011,1,50011,51011,50001
 1001,1011,2011,1,51011,52011,50001
*ELEMENT,TYPE=SC8R,ELSET=BLANK
 2, 11,21,1021,1011, 50011,50021,51021,51011
*ELGEN,ELSET=BLANK
 2, 49,10,1,2,1000,1000
*SHELL SECTION,ELSET=BLANK,MATERIAL=SAMP,ORIENTATION=OR1,POISSON=ELASTIC
 0.00085,9
*ORIENTATION,NAME=OR1,SYSTEM=CYLINDRICAL
 0.0,-1.0,0.0, 0.0,1.0,0.0
 3,0.0
*SURFACE, NAME=TOP
 BLANK,S2
*SURFACE, NAME=BOTTOM
 BLANK,S1
*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=REVOLUTION

 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

 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.0,0.
*NSET,NSET=PUNCH
 1000,
*END INSTANCE
*INSTANCE,NAME=DIE-1,PART=DIE
 0.,0.0,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
*EQUATION
 2
 BLANK-1.501,  1, 1.0, BLANK-1.50501,  1, 1.0
 2
 BLANK-1.1501, 1, 1.0, BLANK-1.51501,  1, 1.0
 2
 BLANK-1.2501, 1, 1.0, BLANK-1.52501, 1, 1.0
*END ASSEMBLY
*BOUNDARY
 ASSEMBLY-1.BLANK-1.1,1,1
 ASSEMBLY-1.BLANK-1.1,3,3
 ASSEMBLY-1.BLANK-1.501,2,3
 ASSEMBLY-1.BLANK-1.1501,2,3
 ASSEMBLY-1.BLANK-1.2501,2,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=24
*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, 
*MONITOR,NODE=ASSEMBLY-1.PUNCH-1.1000,DOF=2
*OUTPUT,FIELD,FREQUENCY=1000
*CONTACT OUTPUT,SLAVE=ASSEMBLY-1.BLANK-1.TOP,NSET=QA_TEST_NSET1
 CSTRESS,
*CONTACT OUTPUT,SLAVE=ASSEMBLY-1.BLANK-1.BOTTOM,NSET=QA_TEST_NSET1
 CSTRESS,
*NODE OUTPUT,NSET=QA_TEST_NSET2
 U,
*ELEMENT OUTPUT,ELSET=QA_TEST_ELSET1
 1,5,9
 S,LE,NE
*END STEP
*STEP,INC=2000,NLGEOM,UNSYMM=YES
*STATIC
 0.05,100.,1E-5
*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
*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,3,3
 ASSEMBLY-1.BLANK-1.501,1,3
 ASSEMBLY-1.BLANK-1.1501,2,3
 ASSEMBLY-1.BLANK-1.2501,2,3
 ASSEMBLY-1.BLANK-1.LOWER,YSYMM
 ASSEMBLY-1.BLANK-1.UPPER,YSYMM
 ASSEMBLY-1.PUNCH-1.1000,1,2
 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
*MONITOR,NODE=ASSEMBLY-1.BLANK-1.1,DOF=2
*END STEP