*HEADING
 WANG AND BUDIANSKY'S SPHERICAL PUNCH WITH M3D4R   
*PREPRINT,ECHO=YES
*RESTART,WRITE,FREQUENCY=40
*PART,NAME=BLANK
   *NODE,NSET=MID
    1,0.0,0.0,0.
   *NODE,NSET=REFD
    401,50.59,0.,0.
   *NFILL,BIAS=1.0,NSET=METND
    MID,REFD,40,10
   *NODE,NSET=END
    501,59.18,0.,0.
   *NFILL,BIAS=1.0,NSET=METND
    REFD,END,10,10
   *NCOPY,OLD SET=METND,NEW SET=METND,CHANGE NUMBER=1000,SHIFT,MULTIPLE=2
    0., 0., 0.,
    0., 0., -1., 0., 0., 1., 5.
   *NSET,NSET=NODWR,GENERATE
    1,501,10
   *NSET,NSET=NODWR1,GENERATE
    1001,1501,10
   *NSET,NSET=MIDDLE,GENERATE
    1011,1501,10
   *NSET,NSET=NODWR2,GENERATE
    2001,2501,10
   *NSET,NSET=LOWER,GENERATE
    11,501,10
   *NSET,NSET=UPPER,GENERATE
    2011,2511,10
   *ELEMENT,TYPE=M3D3,ELSET=TRIS
    1, 11,1011,1
    1001, 1011,2011,1
   *ELEMENT,TYPE=M3D4R,ELSET=QUADS
    2, 11,21,1021,1011
   *ELGEN,ELSET=QUADS
    2, 49,10,1,2,1000,1000
   *ELSET,ELSET=METAL
    TRIS,QUADS
   *ELSET,ELSET=METAL1,GENERATE
    1,40,1
    1001,1040,1
   *ELSET,ELSET=METAL2,GENERATE
    41,50,1
    1041,1050,1
   *MEMBRANE SECTION,ELSET=METAL,MATERIAL=SAMP,ORIENTATION=OR1
    0.85, 
   *ORIENTATION,NAME=OR1,SYSTEM=CYLINDRICAL
    0.0,0.0,-1.0, 0.0,0.0, 1.0
    3,0.0
*END PART
*MATERIAL,NAME=SAMP
*ELASTIC
 206.8,0.3
*PLASTIC
 0.1700000,      0.0000000E+00
 0.1800000 ,     1.7205942E-03
 0.1900000  ,    3.8296832E-03
 0.2000000 ,     6.3897874E-03
 0.2100000,      9.4694765E-03
 0.2200000,      1.3143660E-02
 0.2300000,      1.7493792E-02
 0.2400000,      2.2608092E-02
 0.2500000,      2.8581845E-02
 0.2600000,      3.5517555E-02
 0.2700000,      4.3525275E-02
 0.2800000,      5.2722659E-02
 0.2900000,      6.3235357E-02
 0.3000000,      7.5197279E-02
 0.3100000,      8.8750519E-02
 0.3200000,      0.1040458
 0.3300000,      0.1212430
 0.3400000,      0.1405106
 0.3500000,      0.1620263
 0.3600000,      0.1859779
 0.3700000,      0.2125620
 0.3800000,      0.2419857
 0.3900000,      0.2744660
 0.4000000,      0.3102303
 0.4100000,      0.3495160
 0.4200000,      0.3925720
 0.4300000,      0.4396578
 0.4400000,      0.4910434
 0.4500000,      0.5470111
 0.4600000,      0.6078544
 0.4700000,      0.6738777
 0.4800000,      0.7453985
 0.4900000,      0.8227461
 0.5000000,      0.9062610
 0.5100000 ,     0.9962980
*PART,NAME=PUNCH
   *SURFACE,TYPE=REVOLUTION,NAME=PUNCH

    START,0.0,0.0
    CIRCL,50.8,-50.80,0.0,-50.80
   *NODE,NSET=PUNCH
    1000,0.,0.,0.
   *RIGID BODY,ANALYTICAL SURFACE=PUNCH,REFNODE=1000
*END PART
*PART,NAME=DIE
   *SURFACE,TYPE=REVOLUTION,NAME=DIE
   
    START,61.00,0.05
    LINE,59.18,0.05
    CIRCL,52.83,6.4,59.18,6.4
    LINE,52.83,8.
   *NODE,NSET=DIE
    2000,59.18,0.05,0.0
   *RIGID BODY,ANALYTICAL SURFACE=DIE,REFNODE=2000
*END PART
*ASSEMBLY,NAME=FORM
   *INSTANCE,NAME=BLANK-1,PART=BLANK
      *TRANSFORM,NSET=UPPER,TYPE=C
        0.,0.,0., 0., 0., 1.
      *TRANSFORM,NSET=LOWER,TYPE=C
        0.,0.,0., 0., 0., 1.
      *TRANSFORM,NSET=MIDDLE,TYPE=C
        0.,0.,0., 0., 0., 1.
      *SURFACE, NAME=METAL1
       METAL1,SNEG
      *SURFACE, NAME=METAL2
       METAL2,SPOS
   *END INSTANCE
   *INSTANCE,NAME=PUNCH-1,PART=PUNCH

    0.,0.,0.,1.,0.,0.,90.
   *END INSTANCE
   *INSTANCE,NAME=DIE-1,PART=DIE
   
    0.,0.,0.,1.,0.,0.,90.
   *END INSTANCE
   *NSET,NSET=PUNK
    PUNCH-1.PUNCH,DIE-1.DIE
   *ELSET,ELSET=AFEW
    BLANK-1.19,BLANK-1.20,BLANK-1.21,BLANK-1.22,BLANK-1.23,BLANK-1.24
*END ASSEMBLY
**
**Begin QA_TEST set definitions
**
*ELSET,ELSET=QA_TEST_ELSET1,GENERATE
FORM.BLANK-1.1,FORM.BLANK-1.16,1
FORM.BLANK-1.26,FORM.BLANK-1.50,1
*NSET,NSET=QA_TEST_NSET1
 FORM.BLANK-1.1,FORM.BLANK-1.231
**
**End QA_TEST set definitions
**
*CONTACT PAIR, INTERACTION=SMOOTH
 FORM.BLANK-1.METAL1,FORM.PUNCH-1.PUNCH
*CONTACT PAIR, INTERACTION=SMOOTH
 FORM.BLANK-1.METAL2,FORM.DIE-1.DIE
*SURFACE INTERACTION,NAME=SMOOTH
*FRICTION
 0.275, 
*INITIAL CONDITIONS,TYPE=STRESS
 FORM.BLANK-1.METAL,0.0085,0.0085
*BOUNDARY
 FORM.BLANK-1.1,1,2
 FORM.BLANK-1.1,4,6
 FORM.BLANK-1.501,1,3
 FORM.BLANK-1.1501,1,3
 FORM.BLANK-1.2501,1,3
 FORM.BLANK-1.LOWER,YSYMM
 FORM.BLANK-1.UPPER,YSYMM
 FORM.BLANK-1.MIDDLE,YSYMM
 FORM.PUNCH-1.1000,1,2
 FORM.PUNCH-1.1000,4,6
 FORM.DIE-1.2000,1,2
 FORM.DIE-1.2000,4,6
*STEP,NLGEOM, UNSYMM=YES
*STATIC
1.,1.
*CONTROLS, PARAMETER=FIELD, FIELD=DISPLACEMENT
0.001,
**
** AUTOMATIC TOLERANCES ARE ENTERED TO OVERCOME DIFFERENCES
** IN ROUND OFF CONTACT FORCES IN DIFFERENT PLATFORMS
**
**CONTACT CONTROLS, AUTOMATIC TOLERANCES
*BOUNDARY
FORM.PUNCH-1.1000,3,3
FORM.DIE-1.2000,3,3,-0.05
*PRINT,RESIDUAL=NO,FREQUENCY=10,CONTACT=YES
*EL PRINT,FREQUENCY=100
*NODE FILE,NSET=FORM.PUNK,FREQUENCY=1000
U,
COORD,
*OUTPUT,HISTORY,FREQUENCY=0 
*OUTPUT,FIELD,FREQUENCY=1000
*ELEMENT OUTPUT,ELSET=FORM.BLANK-1.METAL 
S,LE,PEEQ
*NODE OUTPUT
U,
*END STEP
*STEP,NLGEOM,UNSYMM=YES
*STATIC
1.0,100.
*BOUNDARY
FORM.PUNCH-1.1000,3,3,18.6
FORM.DIE-1.2000,3,3,-0.05
*MONITOR,NODE=FORM.PUNCH-1.1000,DOF=3
*NODE PRINT,FREQUENCY=100
U,
RF,
*OUTPUT,HISTORY,FREQUENCY=0
*OUTPUT,FIELD,FREQUENCY=1000
*ELEMENT OUTPUT,ELSET=FORM.BLANK-1.METAL 
S,LE,PEEQ
*NODE OUTPUT
U,
*END STEP
*STEP,NLGEOM,INC=200, UNSYMM=YES
*STATIC
1.0,100.,
*BOUNDARY
FORM.PUNCH-1.1000,3,3,28.5
FORM.DIE-1.2000,3,3,-0.05
*EL FILE,ELSET=FORM.AFEW,FREQUENCY=1000
 S,
 E,
 PE,
*OUTPUT,HISTORY,FREQUENCY=0
*OUTPUT,FIELD,FREQUENCY=1000
*ELEMENT OUTPUT,ELSET=FORM.AFEW
 S,
 E,
 PE, 
*ELEMENT OUTPUT,ELSET=FORM.BLANK-1.METAL 
S,LE,PEEQ
*NODE OUTPUT
U,
*END STEP
*STEP,INC=2000,NLGEOM, UNSYMM=YES
*STATIC
1.0,100.
*BOUNDARY
FORM.PUNCH-1.1000,3,3,34.5
FORM.DIE-1.2000,3,3,-0.05
*OUTPUT,HISTORY,FREQUENCY=0
*OUTPUT,FIELD,FREQUENCY=1000
*ELEMENT OUTPUT,ELSET=FORM.BLANK-1.METAL 
S,LE,PEEQ
*NODE OUTPUT
U,
**
*END STEP
*STEP,INC=2000,NLGEOM, UNSYMM=YES
*STATIC
100.,100.
*BOUNDARY,FIXED
FORM.BLANK-1.METND,1,6
*MODEL CHANGE,REMOVE,TYPE=CONTACT PAIR
FORM.BLANK-1.METAL1,FORM.PUNCH-1.PUNCH
FORM.BLANK-1.METAL2,FORM.DIE-1.DIE
*EL FILE,ELSET=FORM.BLANK-1.METAL,FREQUENCY=1000 
S,
*OUTPUT,HISTORY,FREQUENCY=0 
*OUTPUT,FIELD,FREQUENCY=1000
*ELEMENT OUTPUT,ELSET=FORM.BLANK-1.METAL 
S,LE,PEEQ
*NODE OUTPUT
U,
** Begin QA_TEST output request
** Too sensitive across platforms. 
**ELEMENT OUTPUT,ELSET=QA_TEST_ELSET1
**LE,
*NODE OUTPUT,NSET=QA_TEST_NSET1
U,
**
** End QA_TEST output request
**
*END STEP
**
** STEP 6 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.
**