*HEADING
 This input deck verifies the basic modes of the continuum shell (SC8R).  
 All dof are controlled.  Small displacement analysis. The element shape 
 is regular i.e mid-plane is a 2x2 square and has a uniform thickness of 0.1.
 Perturbation steps.
*PREPRINT, ECHO=NO, HISTORY=NO, MODEL=NO
*NODE, NSET=NSC8R
  1, -1.0, -1.0, -0.05
  2,  1.0, -1.0, -0.05
  3,  1.0,  1.0, -0.05
  4, -1.0,  1.0, -0.05
  5, -1.0, -1.0,  0.05
  6,  1.0, -1.0,  0.05
  7,  1.0,  1.0,  0.05
  8, -1.0,  1.0,  0.05
*NODE, NSET=NS4R
  11, -1.0, -1.0, 0.0
  12,  1.0, -1.0, 0.0
  13,  1.0,  1.0, 0.0
  14, -1.0,  1.0, 0.0
*NODE, NSET=NC3D8R
  21, -1.0, -1.0, -0.05
  22,  1.0, -1.0, -0.05
  23,  1.0,  1.0, -0.05
  24, -1.0,  1.0, -0.05
  25, -1.0, -1.0,  0.05
  26,  1.0, -1.0,  0.05
  27,  1.0,  1.0,  0.05
  28, -1.0,  1.0,  0.05
*ELEMENT,TYPE=SC8R, ELSET=ELSET_SC8R
  1, 1, 2, 3, 4, 5, 6, 7 , 8
*ELEMENT,TYPE=S4R, ELSET=ELSET_S4R
 10, 11,12,13,14
*ELEMENT,TYPE=C3D8R, ELSET=ELSET_C3D8R
 20, 21, 22, 23, 24, 25, 26, 27 , 28
*ELSET, ELSET=SHELLS
 ELSET_SC8R,ELSET_S4R
*SHELL SECTION,MATERIAL=A1, ELSET=ELSET_SC8R, POISSON=ELASTIC
 0.1, 
*SHELL SECTION,MATERIAL=A1, ELSET=ELSET_S4R, POISSON=0.0
 0.1, 
*SOLID SECTION,MATERIAL=A1,ELSET=ELSET_C3D8R
*MATERIAL,NAME=A1
*ELASTIC,TYPE=ISOTROPIC
 3000.0,0.0
*DENSITY
 1.0
**
**
** *******************************************************************
**      MEMBRANE HOURGLASS MODE (TOP AND BOTTOM SURFACE IN-SYNC)
** *******************************************************************
**
** Note: Comparison are made with S4R element. Artificial energies 
**       should be identical. Strains/Stresses should be zero.
**       SC8R reaction forces should be 1/2 reported by S4R
**
*STEP, PERT
 Step 1: Membrane Hourglass Mode (Top and Bottom in Sync in y-dir)
*STATIC
*BOUNDARY, OP=NEW
 NSC8R, 1, 1, 0.0
 NSC8R, 3, 3, 0.0
 1, 2, 2, +0.1
 2, 2, 2, -0.1
 3, 2, 2, +0.1
 4, 2, 2, -0.1
 5, 2, 2, +0.1
 6, 2, 2, -0.1
 7, 2, 2, +0.1
 8, 2, 2, -0.1
*BOUNDARY, OP=NEW
 NS4R,1,1,  0.0
 NS4R,3,5,  0.0
 11, 2, 2, +0.1
 12, 2, 2, -0.1
 13, 2, 2, +0.1
 14, 2, 2, -0.1
*BOUNDARY, OP=NEW
 NC3D8R, 1,3, 0.0
**
** ODB OUTPUT
**
*OUTPUT,FIELD, FREQ=1
*ELEMENT OUTPUT, ELSET=SHELLS
 E,
 S,
 SF,
 SE,
 STH
*NODE OUTPUT
 U,
 RF,
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT, ELSET=SHELLS
**
** DAT OUTPUT
**
*EL PRINT, ELSET=SHELLS, FREQ=1
 E,
 S,
 SF,
 SE,
 STH,
*NODE PRINT, FREQ=1
 U,
 RF,
*ENERGY PRINT, ELSET=ELSET_SC8R, FREQ=1
*ENERGY PRINT, ELSET=ELSET_S4R, FREQ=1
**
** FIL OUTPUT
**
*EL FILE, ELSET=SHELLS, FREQ=1
 E,
 S,
 SF,
 SE,
 STH,
*NODE FILE, FREQ=1
 U,
 RF,
*ENERGY FILE, ELSET=ELSET_SC8R, FREQ=1
*ENERGY FILE, ELSET=ELSET_S4R, FREQ=1
*END STEP
**
**
*STEP, PERT
 Step 2: Membrane Hourglass Mode (Top and Bottom in Sync in x-dir)
*STATIC
*BOUNDARY, OP=NEW
 NSC8R, 2, 3, 0.0
 1, 1, 1, -0.1
 2, 1, 1, +0.1
 3, 1, 1, -0.1
 4, 1, 1, +0.1
 5, 1, 1, -0.1
 6, 1, 1, +0.1
 7, 1, 1, -0.1
 8, 1, 1, +0.1
*BOUNDARY, OP=NEW
 NS4R,2,5,  0.0
 11, 1, 1, -0.1
 12, 1, 1, +0.1
 13, 1, 1, -0.1
 14, 1, 1, +0.1
*BOUNDARY, OP=NEW
 NC3D8R, 1,3, 0.0
*END STEP
**
**
** *******************************************************************
**   DIRECTOR HOURGLASS MODES (Similar to bending hg modes in S4R)
** *******************************************************************
**
** Note: For SC8R element this mode is handled via the top and bottom
**       membrane hourglass modes. All strains/stress should be zero.
**       Energies and reaction forces may differ between SC8R and
**       S4R, since these are handled by different hourglass control.
**       This mode was added for comparison purpuses only.
**
*STEP, PERT
 STEP 3: Director Hourglass Mode: y-dir
*STATIC
*BOUNDARY, OP=NEW
 NSC8R, 1, 1, 0.0
 NSC8R, 3, 3, 0.0
 1, 2, 2, +0.055357
 2, 2, 2, -0.055357
 3, 2, 2, +0.055357
 4, 2, 2, -0.055357
 5, 2, 2, -0.055357
 6, 2, 2, +0.055357
 7, 2, 2, -0.055357
 8, 2, 2, +0.055357
*BOUNDARY, OP=NEW
 NS4R,1,3,  0.0
 NS4R,5,6,  0.0
 11, 4, 4, +1.107149
 12, 4, 4, -1.107149
 13, 4, 4, +1.107149
 14, 4, 4, -1.107149
*BOUNDARY, OP=NEW
 NC3D8R, 1,3, 0.0
*END STEP
**
**
*STEP, PERT
 Step 4: Director Hourglass Mode: x-dir
*STATIC
*BOUNDARY, OP=NEW
 NSC8R, 2, 3, 0.0
 1, 1, 1, -0.055357
 2, 1, 1, +0.055357
 3, 1, 1, -0.055357
 4, 1, 1, +0.055357
 5, 1, 1, +0.055357
 6, 1, 1, -0.055357
 7, 1, 1, +0.055357
 8, 1, 1, -0.055357
*BOUNDARY, OP=NEW
 NS4R,1,4,  0.0
 NS4R,6,6,  0.0
 11, 5, 5, +1.107149
 12, 5, 5, -1.107149
 13, 5, 5, +1.107149
 14, 5, 5, -1.107149
*BOUNDARY, OP=NEW
 NC3D8R, 1,3, 0.0
*END STEP
**
**
** *******************************************************************
**                  OUT-OF-SYN W-MODE HOURGLASS
** *******************************************************************
**
**  Note: Strains/Stresses should be zero. Reaction forces and 
**        artificial energy of SC8R are matched to C3D8R to give 
**        equivalent response.
**
*STEP, PERT
 STEP 5: Out-of-Sync W-mode Hourglass Mode
*STATIC
*BOUNDARY, OP=NEW
 NSC8R, 1, 2, 0.0
 1, 3, 3, +0.01
 2, 3, 3, -0.01
 3, 3, 3, +0.01
 4, 3, 3, -0.01
 5, 3, 3, -0.01
 6, 3, 3, +0.01
 7, 3, 3, -0.01
 8, 3, 3, +0.01
*BOUNDARY, OP=NEW
 NS4R, 1, 6,  0.0
*BOUNDARY, OP=NEW
 NC3D8R, 1, 2, 0.0
 21, 3, 3, +0.01
 22, 3, 3, -0.01
 23, 3, 3, +0.01
 24, 3, 3, -0.01
 25, 3, 3, -0.01
 26, 3, 3, +0.01
 27, 3, 3, -0.01
 28, 3, 3, +0.01
**
** ODB OUTPUT
**
*OUTPUT,FIELD, FREQ=1
*ELEMENT OUTPUT, ELSET=SHELLS
 E,
 S,
 SF,
 SE,
 STH
*NODE OUTPUT
 U,
 RF,
**
** DAT OUTPUT
**
*EL PRINT, ELSET=ELSET_SC8R, FREQ=1
 E,
 S,
 SF,
 SE,
 STH,
*EL PRINT, ELSET=ELSET_C3D8R, FREQ=1
 E,
 S,
*NODE PRINT, FREQ=1
 U,
 RF,
*ENERGY PRINT, ELSET=ELSET_SC8R, FREQ=1
*ENERGY PRINT, ELSET=ELSET_C3D8R, FREQ=1
**
** FIL OUTPUT
**
*EL FILE, ELSET=ELSET_SC8R, FREQ=1
 E,
 S,
 SF,
 SE,
 STH,
*EL FILE, ELSET=ELSET_C3D8R, FREQ=1
 E,
 S,
*NODE FILE, FREQ=1
 U,
 RF,
*ENERGY FILE, ELSET=ELSET_SC8R, FREQ=1
*ENERGY FILE, ELSET=ELSET_C3D8R, FREQ=1
*END STEP
**
**
** *******************************************************************
**    Transverse Shear Moment Mode  (Symmetric Shear/Shear Couple)
** *******************************************************************
**
**   Note: Strains/Stresses should be zero. C3D8R included just
**         for comparison. Solutions are not expected to match.
**
*STEP, PERT
 STEP 6:  Transverse shear moment mode : x-dir
*STATIC
*BOUNDARY, OP=NEW
 NSC8R, 1, 2, 0.0
 1, 3, 3, +0.01
 2, 3, 3, -0.01
 3, 3, 3, -0.01
 4, 3, 3, +0.01
 5, 3, 3, -0.01
 6, 3, 3, +0.01
 7, 3, 3, +0.01
 8, 3, 3, -0.01
*BOUNDARY, OP=NEW
 NS4R, 1, 6,  0.0
*BOUNDARY, OP=NEW
 NC3D8R, 1, 2, 0.0
 21, 3, 3, +0.01
 22, 3, 3, -0.01
 23, 3, 3, -0.01
 24, 3, 3, +0.01
 25, 3, 3, -0.01
 26, 3, 3, +0.01
 27, 3, 3, +0.01
 28, 3, 3, -0.01
*END STEP
**
**
*STEP, PERT
 STEP 7:  Transverse shear moment mode : y-dir
*STATIC
*BOUNDARY, OP=NEW
 NSC8R, 1, 2, 0.0
 1, 3, 3, +0.01
 2, 3, 3, +0.01
 3, 3, 3, -0.01
 4, 3, 3, -0.01
 5, 3, 3, -0.01
 6, 3, 3, -0.01
 7, 3, 3, +0.01
 8, 3, 3, +0.01
*BOUNDARY, OP=NEW
 NS4R, 1, 6,  0.0
*BOUNDARY, OP=NEW
 NC3D8R, 1, 2, 0.0
 21, 3, 3, +0.01
 22, 3, 3, +0.01
 23, 3, 3, -0.01
 24, 3, 3, -0.01
 25, 3, 3, -0.01
 26, 3, 3, -0.01
 27, 3, 3, +0.01
 28, 3, 3, +0.01
*END STEP
**
**
** *******************************************************************
**                      Thickness Mode
** *******************************************************************
**
** Note:  SC8R and C3D8R results should match.
**
*STEP, PERT
 STEP 8:thickness mode
*STATIC
*BOUNDARY, OP=NEW
 NSC8R, 1, 2, 0.0
 1, 3, 3, +0.04
 2, 3, 3, +0.04
 3, 3, 3, +0.04
 4, 3, 3, +0.04
 5, 3, 3, -0.04
 6, 3, 3, -0.04
 7, 3, 3, -0.04
 8, 3, 3, -0.04
*BOUNDARY, OP=NEW
 NS4R, 1, 6,  0.0
*BOUNDARY, OP=NEW
 NC3D8R, 1, 2, 0.0
 21, 3, 3, +0.04
 22, 3, 3, +0.04
 23, 3, 3, +0.04
 24, 3, 3, +0.04
 25, 3, 3, -0.04
 26, 3, 3, -0.04
 27, 3, 3, -0.04
 28, 3, 3, -0.04
*END STEP
**
**
**********************************************************************
**                     Transverse Shear 
**********************************************************************
**
** Note: S4R and SC8R results should match In plane disp. causing tvs 
**       shear. This does not cause in-plane elongation.
**
*STEP, PERT
 Step 9: Transverse Shear in 2-dir
*STATIC
*BOUNDARY, OP=NEW
 NSC8R, 1, 1, 0.0
 NSC8R, 3, 3, 0.0
 1, 2, 2, 0.0
 2, 2, 2, 0.0
 3, 2, 2, 0.0
 4, 2, 2, 0.0
 5, 2, 2, 0.13734
 6, 2, 2, 0.13734
 7, 2, 2, 0.13734
 8, 2, 2, 0.13734
*BOUNDARY, OP=NEW
 NS4R, 1, 3, 0.0
 NS4R, 5, 6, 0.0
 11, 4, 4, -1.3734
 12, 4, 4, -1.3734
 13, 4, 4, -1.3734
 14, 4, 4, -1.3734
*BOUNDARY, OP=NEW
 NC3D8R, 1,3, 0.0
** ODB OUTPUT
**
*OUTPUT,FIELD, FREQ=1
*ELEMENT OUTPUT, ELSET=SHELLS
 E,
 S,
 SF,
 SE,
 STH
*NODE OUTPUT
 U,
 RF,
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT, ELSET=SHELLS
**
** DAT OUTPUT
**
*EL PRINT, ELSET=SHELLS, FREQ=1
 E,
 S,
 SF,
 SE,
 STH,
*NODE PRINT, FREQ=1
 U,
 RF,
*ENERGY PRINT, ELSET=ELSET_SC8R, FREQ=1
*ENERGY PRINT, ELSET=ELSET_S4R, FREQ=1
**
** FIL OUTPUT
**
*EL FILE, ELSET=SHELLS, FREQ=1
 E,
 S,
 SF,
 SE,
 STH,
*NODE FILE, FREQ=1
 U,
 RF,
*ENERGY FILE, ELSET=ELSET_SC8R, FREQ=1
*ENERGY FILE, ELSET=ELSET_S4R, FREQ=1
*END STEP
**
**
*STEP, PERT
 Step 10: Transverse Shear in x-dir
*STATIC
*BOUNDARY, OP=NEW
 NSC8R, 2, 3, 0.0
 1, 1, 1, 0.0
 2, 1, 1, 0.0
 3, 1, 1, 0.0
 4, 1, 1, 0.0
 5, 1, 1, -0.13734
 6, 1, 1, -0.13734
 7, 1, 1, -0.13734
 8, 1, 1, -0.13734
*BOUNDARY, OP=NEW
 NS4R,1,4,  0.0
 NS4R,6,6,  0.0
** Linear geometry
 11, 5, 5, -1.3734
 12, 5, 5, -1.3734
 13, 5, 5, -1.3734
 14, 5, 5, -1.3734
*BOUNDARY, OP=NEW
 NC3D8R, 1,3, 0.0
*END STEP
**
**
**********************************************************************
**                    Butterfly Hourglass Mode
**********************************************************************
**
** Note:  S4R and SC8R results should match.
**
*STEP, PERT
 Step 11: Butterfly Hourglass Mode
*STATIC
*BOUNDARY, OP=NEW
 NSC8R, 1, 1, 0.0
 NSC8R, 2, 2, 0.0
 1, 3, 3, +0.1
 2, 3, 3, -0.1
 3, 3, 3, +0.1
 4, 3, 3, -0.1
 5, 3, 3, +0.1
 6, 3, 3, -0.1
 7, 3, 3, +0.1
 8, 3, 3, -0.1
*BOUNDARY, OP=NEW
 NS4R,1,2,  0.0
 NS4R,4,6,  0.0
 11, 3, 3,  0.1
 12, 3, 3, -0.1
 13, 3, 3,  0.1
 14, 3, 3, -0.1
*BOUNDARY, OP=NEW
 NC3D8R, 1,3, 0.0
*END STEP
**
**
**********************************************************************
**                        Crop Circle Mode
**********************************************************************
**
** Note: The results of S4R and SC8R should match. 
**
*STEP, PERT
 Step 12: Crop Circle Hourglass Mode
*STATIC
*BOUNDARY, OP=NEW
 NSC8R, 3, 3, 0.0
 1,1,1,-3.5355E-03 
 2,1,1,-3.5355E-03 
 3,1,1, 3.5355E-03 
 4,1,1, 3.5355E-03 
 5,1,1, 3.5355E-03 
 6,1,1, 3.5355E-03 
 7,1,1,-3.5355E-03 
 8,1,1,-3.5355E-03 
 1,2,2, 3.5355E-03 
 2,2,2,-3.5355E-03 
 3,2,2,-3.5355E-03 
 4,2,2, 3.5355E-03 
 5,2,2,-3.5355E-03 
 6,2,2, 3.5355E-03 
 7,2,2, 3.5355E-03 
 8,2,2,-3.5355E-03 
*BOUNDARY, OP=NEW
 NS4R,1,3,  0.0
 NS4R,6,6,  0.0
 11, 4,4,  0.07045
 12, 4,4, -0.07045
 13, 4,4, -0.07045
 14, 4,4,  0.07045
 11, 5,5, +0.07045
 12, 5,5, +0.07045
 13, 5,5, -0.07045
 14, 5,5, -0.07045
*BOUNDARY, OP=NEW
 NC3D8R, 1,3, 0.0
*END STEP
**
**
** *******************************************************************
**                        REGULAR MEMBRANE MODES
** *******************************************************************
**
** Note: Results should match that of S4R.
**
*STEP, PERT
 STEP 13: In-plane 11-strain.
*STATIC
*BOUNDARY, OP=NEW
 NSC8R, 2, 3, 0.0
 1, 1, 1, -0.8
 2, 1, 1, +0.8
 3, 1, 1, +0.8
 4, 1, 1, -0.8
 5, 1, 1, -0.8
 6, 1, 1, +0.8
 7, 1, 1, +0.8
 8, 1, 1, -0.8
*BOUNDARY, OP=NEW
 NS4R,2,6,  0.0
 11, 1, 1, -0.8
 12, 1, 1, +0.8
 13, 1, 1, +0.8
 14, 1, 1, -0.8
*BOUNDARY, OP=NEW
 NC3D8R, 1,3, 0.0
*END STEP
**
**
*STEP, PERT
 STEP 14: In-plane 22-strain.
*STATIC
*BOUNDARY, OP=NEW
 NSC8R, 1, 1, 0.0
 NSC8R, 3, 3, 0.0
 1, 2, 2, -0.8
 2, 2, 2, -0.8
 3, 2, 2, +0.8
 4, 2, 2, +0.8
 5, 2, 2, -0.8
 6, 2, 2, -0.8
 7, 2, 2, +0.8
 8, 2, 2, +0.8
*BOUNDARY, OP=NEW
 NS4R,1,1,  0.0
 NS4R,3,6,  0.0
 11, 2, 2, -0.8
 12, 2, 2, -0.8
 13, 2, 2, +0.8
 14, 2, 2, +0.8
*BOUNDARY, OP=NEW
 NC3D8R, 1,3, 0.0
*END STEP
**
**
*STEP, PERT
 STEP 15: In-plane 12-strain.
*STATIC
*BOUNDARY, OP=NEW
 NSC8R, 3, 3, 0.0
 1, 2, 2, -0.8
 2, 2, 2, +0.8
 3, 2, 2, +0.8
 4, 2, 2, -0.8
 1, 1, 1, -0.8
 2, 1, 1, -0.8
 3, 1, 1, +0.8
 4, 1, 1, +0.8
 5, 2, 2, -0.8
 6, 2, 2, +0.8
 7, 2, 2, +0.8
 8, 2, 2, -0.8
 5, 1, 1, -0.8
 6, 1, 1, -0.8
 7, 1, 1, +0.8
 8, 1, 1, +0.8
*BOUNDARY, OP=NEW
 NS4R,3,6,  0.0
 11, 2, 2, -0.8
 12, 2, 2, +0.8
 13, 2, 2, +0.8
 14, 2, 2, -0.8
 11, 1, 1, -0.8
 12, 1, 1, -0.8
 13, 1, 1, +0.8
 14, 1, 1, +0.8
*BOUNDARY, OP=NEW
 NC3D8R, 1,3, 0.0
*END STEP
**
**
** *******************************************************************
**                   REGULAR BENDING MODES
** *******************************************************************
**
**  Note: Results should match that of S4R.
**
*STEP, PERT
 STEP 16: In-plane 11-curvature.
*STATIC
*BOUNDARY, OP=NEW
 NSC8R, 2, 3, 0.0
 1, 1, 1, -0.0754189
 2, 1, 1, +0.0754189
 3, 1, 1, +0.0754189
 4, 1, 1, -0.0754189
 5, 1, 1, +0.0754189
 6, 1, 1, -0.0754189
 7, 1, 1, -0.0754189
 8, 1, 1, +0.0754189
*BOUNDARY, OP=NEW
 NS4R,1,4,  0.0
 NS4R,6,6,  0.0
 11, 5, 5, +1.5083775
 12, 5, 5, -1.5083775
 13, 5, 5, -1.5083775
 14, 5, 5, +1.5083775
*BOUNDARY, OP=NEW
 NC3D8R, 1,3, 0.0
*END STEP
**
**
*STEP, PERT
 STEP 17:  22-curvature.
*STATIC
*BOUNDARY, OP=NEW
 NSC8R, 1, 1, 0.0
 NSC8R, 3, 3, 0.0
 1, 2, 2, -0.0754189
 2, 2, 2, -0.0754189
 3, 2, 2, +0.0754189
 4, 2, 2, +0.0754189
 5, 2, 2, +0.0754189
 6, 2, 2, +0.0754189
 7, 2, 2, -0.0754189
 8, 2, 2, -0.0754189
*BOUNDARY, OP=NEW
 NS4R,1,3,  0.0
 NS4R,5,6,  0.0
 11, 4, 4, -1.5083775
 12, 4, 4, -1.5083775
 13, 4, 4, +1.5083775
 14, 4, 4, +1.5083775
*BOUNDARY, OP=NEW
 NC3D8R, 1,3, 0.0
*END STEP
*STEP, PERT
 STEP 18: 12-curvature.
*STATIC
*BOUNDARY, OP=NEW
 NSC8R, 3, 3, 0.0
 1, 1, 1, +0.0535
 1, 2, 2, +0.0535
 2, 1, 1, +0.0535
 2, 2, 2, -0.0535
 3, 1, 1, -0.0535
 3, 2, 2, -0.0535
 4, 1, 1, -0.0535
 4, 2, 2, +0.0535
 5, 1, 1, -0.0535
 5, 2, 2, -0.0535
 6, 1, 1, -0.0535
 6, 2, 2, +0.0535
 7, 1, 1, +0.0535
 7, 2, 2, +0.0535
 8, 1, 1, +0.0535
 8, 2, 2, -0.0535
*BOUNDARY, OP=NEW
 NS4R,1,3,  0.0
 NS4R,6,6,  0.0
 11, 4, 4, +1.07
 12, 4, 4, -1.07
 13, 4, 4, -1.07
 14, 4, 4, +1.07
 11, 5, 5, -1.07
 12, 5, 5, -1.07
 13, 5, 5, +1.07
 14, 5, 5, +1.07
*BOUNDARY, OP=NEW
 NC3D8R, 1,3, 0.0
*END STEP