*HEADING
 WANG AND BUDIANSKY'S SPHERICAL PUNCH WITH M3D4R
*PREPRINT,ECHO=YES
*RESTART,WRITE,FREQUENCY=40
*PART,NAME=BLANK
   *NODE,NSET=MID
   *NODE,NSET=REFD
   *NFILL,BIAS=1.0,NSET=METND
    MID,REFD,40,10
   *NODE,NSET=END
   *NFILL,BIAS=1.0,NSET=METND
    REFD,END,10,10
   *NCOPY,OLD SET=METND,NEW SET=METND,CHANGE NUMBER=1000,SHIFT,MULTIPLE=2
   *NSET,NSET=NODWR,GENERATE
   *NSET,NSET=NODWR1,GENERATE
   *NSET,NSET=MIDDLE,GENERATE
   *NSET,NSET=NODWR2,GENERATE
   *NSET,NSET=LOWER,GENERATE
   *NSET,NSET=UPPER,GENERATE
   *ELEMENT,TYPE=M3D3,ELSET=TRIS
   *ELEMENT,TYPE=M3D4R,ELSET=QUADS
   *ELGEN,ELSET=QUADS
   *ELSET,ELSET=METAL
    TRIS,QUADS
   *ELSET,ELSET=METAL1,GENERATE
   *ELSET,ELSET=METAL2,GENERATE
   *MEMBRANE SECTION,ELSET=METAL,MATERIAL=SAMP,ORIENTATION=OR1
    0.85,
   *ORIENTATION,NAME=OR1,SYSTEM=CYLINDRICAL
*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
*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
   *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
   *RIGID BODY,ANALYTICAL SURFACE=DIE,REFNODE=2000
*END PART
*ASSEMBLY,NAME=FORM
   *INSTANCE,NAME=BLANK-1,PART=BLANK
      *TRANSFORM,NSET=UPPER,TYPE=C
      *TRANSFORM,NSET=LOWER,TYPE=C
      *TRANSFORM,NSET=MIDDLE,TYPE=C
      *SURFACE, NAME=METAL1
       METAL1,SNEG
      *SURFACE, NAME=METAL2
       METAL2,SPOS
   *END INSTANCE
   *INSTANCE,NAME=PUNCH-1,PART=PUNCH
   *END INSTANCE
   *INSTANCE,NAME=DIE-1,PART=DIE
   *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
**
*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.
**