*HEADING
 Basic element modes for large deformation [SC8R]
**
*PREPRINT, ECHO=NO, MODEL=NO, HISTORY=NO
**
**
*PART, NAME=ELEMENTS
*NODE, NSET=NSC8R
  900000001, -1.0, -1.0, -0.05
  900000002,  1.0, -1.0, -0.05
  900000003,  1.0,  1.0, -0.05
  900000004, -1.0,  1.0, -0.05
  900000005, -1.0, -1.0,  0.05
  900000006,  1.0, -1.0,  0.05
  900000007,  1.0,  1.0,  0.05
  900000008, -1.0,  1.0,  0.05
*ELEMENT,TYPE=SC8R, ELSET=ELSET_SC8R
  100000000, 900000001, 900000002, 900000003, 900000004, 900000005, 900000006, 900000007 , 900000008
*SHELL SECTION,MATERIAL=MAT_1, ELSET=ELSET_SC8R, POISSON=ELASTIC
 0.1, 
*NODE, NSET=NS4R
  900000011, -1.0, -1.0, 0.0
  900000012,  1.0, -1.0, 0.0
  900000013,  1.0,  1.0, 0.0
  900000014, -1.0,  1.0, 0.0
*ELEMENT,TYPE=S4R, ELSET=ELSET_S4R
  200000000, 900000011,900000012,900000013,900000014
*SHELL SECTION,MATERIAL=MAT_1, ELSET=ELSET_S4R, POISSON=0.0
 0.1, 
*NODE, NSET=NC3D8R
  900000021, -1.0, -1.0, -0.05
  900000022,  1.0, -1.0, -0.05
  900000023,  1.0,  1.0, -0.05
  900000024, -1.0,  1.0, -0.05
  900000025, -1.0, -1.0,  0.05
  900000026,  1.0, -1.0,  0.05
  900000027,  1.0,  1.0,  0.05
  900000028, -1.0,  1.0,  0.05
*ELEMENT,TYPE=C3D8R, ELSET=ELSET_C3D8R
  300000000, 900000021, 900000022, 900000023, 900000024, 900000025, 900000026, 900000027, 900000028
*SOLID SECTION,MATERIAL=MAT_1,ELSET=ELSET_C3D8R
*END PART
*MATERIAL,NAME=MAT_1
*ELASTIC,TYPE=ISOTROPIC
 3000.0,0.0
*DENSITY
 1.0
**
**
*ASSEMBLY, NAME=ASSEMBLY
*INSTANCE, NAME=E11, PART=ELEMENTS
 0.0, 0.0, 0.0
*END INSTANCE
*INSTANCE, NAME=E22, PART=ELEMENTS
 5.0, 0.0, 0.0
*END INSTANCE
*INSTANCE, NAME=E12, PART=ELEMENTS
 10.0, 0.0, 0.0
*END INSTANCE
*INSTANCE, NAME=K11, PART=ELEMENTS
 0.0, -5.0, 0.0
*END INSTANCE
*INSTANCE, NAME=K22, PART=ELEMENTS
 5.0, -5.0, 0.0
*END INSTANCE
*INSTANCE, NAME=K12, PART=ELEMENTS
 10.0, -5.0, 0.0
*END INSTANCE
*INSTANCE, NAME=T, PART=ELEMENTS
 0.0, -10.0, 0.0
*END INSTANCE
*INSTANCE, NAME=TGRAD1, PART=ELEMENTS
 5.0, -10.0, 0.0
*END INSTANCE
*INSTANCE, NAME=TGRAD2, PART=ELEMENTS
 10.0, -10.0, 0.0
*END INSTANCE
*INSTANCE, NAME=TVS1, PART=ELEMENTS
 5.0, -15.0, 0.0
*END INSTANCE
*INSTANCE, NAME=TVS2, PART=ELEMENTS
 10.0, -15.0, 0.0
*END INSTANCE
*INSTANCE, NAME=HGBUTTERFLY, PART=ELEMENTS
 0.0, -20.0, 0.0
*END INSTANCE
*INSTANCE, NAME=HGCROP, PART=ELEMENTS
 5.0, -20.0, 0.0
*END INSTANCE
*INSTANCE, NAME=HGMEMBRANE1, PART=ELEMENTS
 0.0, -25.0, 0.0
*END INSTANCE
*INSTANCE, NAME=HGMEMBRANE2, PART=ELEMENTS
 5.0, -25.0, 0.0
*END INSTANCE
*INSTANCE, NAME=HGROT, PART=ELEMENTS
 0.0, -30.0, 0.0
*END INSTANCE
*INSTANCE, NAME=HGWMODE, PART=ELEMENTS
 0.0, -35.0, 0.0
*END INSTANCE
*NSET, NSET=NOUT_MEMBRANE
 E11.NSC8R, E22.NSC8R, E12.NSC8R,
 E11.NS4R , E22.NS4R , E12.NS4R
*ELSET, ELSET=EOUT_MEMBRANE
 E11.ELSET_SC8R, E22.ELSET_SC8R, E12.ELSET_SC8R,
 E11.ELSET_S4R , E22.ELSET_S4R , E12.ELSET_S4R
*NSET, NSET=NOUT_CURVATURE
 K11.NSC8R, K22.NSC8R, K12.NSC8R,
 K11.NS4R , K22.NS4R , K12.NS4R
*ELSET, ELSET=EOUT_CURVATURE
 K11.ELSET_SC8R, K22.ELSET_SC8R, K12.ELSET_SC8R,
 K11.ELSET_S4R , K22.ELSET_S4R , K12.ELSET_S4R
*NSET, NSET=NOUT_THICK
 T.NSC8R, T.NC3D8R,  TGRAD1.NSC8R, TGRAD1.NC3D8R, TGRAD2.NSC8R, TGRAD2.NC3D8R 
*ELSET, ELSET=EOUT_THICK
 T.ELSET_SC8R, T.ELSET_C3D8R, TGRAD1.ELSET_SC8R, TGRAD1.ELSET_C3D8R, TGRAD2.ELSET_SC8R, TGRAD2.ELSET_C3D8R 
*NSET, NSET=NOUT_TVS
 TVS1.NSC8R, TVS2.NSC8R,
 TVS1.NS4R , TVS1.NS4R 
*ELSET, ELSET=EOUT_TVS
 TVS1.ELSET_SC8R, TVS2.ELSET_SC8R,
 TVS1.ELSET_S4R , TVS2.ELSET_S4R 
*NSET, NSET=NOUT_HGMODES
 HGBUTTERFLY.NSC8R, HGCROP.NSC8R, HGMEMBRANE1.NSC8R,, HGMEMBRANE2.NSC8R,, HGROT.NSC8R,
 HGBUTTERFLY.NS4R, HGCROP.NS4R, HGMEMBRANE1.NS4R,, HGMEMBRANE2.NS4R,, HGROT.NS4R,
*ELSET, ELSET=EOUT_HGMODES
 HGBUTTERFLY.ELSET_SC8R, HGCROP.ELSET_SC8R, HGMEMBRANE1.ELSET_SC8R,, HGMEMBRANE2.ELSET_SC8R,, HGROT.ELSET_SC8R
 HGBUTTERFLY.ELSET_S4R, HGCROP.ELSET_S4R, HGMEMBRANE1.ELSET_S4R,, HGMEMBRANE2.ELSET_S4R,, HGROT.ELSET_S4R
*NSET, NSET=NOUT_HGWMODE
 HGWMODE.NSC8R, HGWMODE.NC3D8R
*END ASSEMBLY
*AMPLITUDE, NAME=K22_uy, DEFINITION=PERIODIC
 1, 0.78539816, 0., 0.
 0., 1.0
*AMPLITUDE, NAME=K22_uz, DEFINITION=PERIODIC
 1, 0.78539816, 0., 1.0
 -1.0, 0.
*AMPLITUDE, NAME=K11_ux, DEFINITION=PERIODIC
 1, 0.78539816, 0., 0.
 0., 1.0
*AMPLITUDE, NAME=K11_uz, DEFINITION=PERIODIC
 1, 0.78539816, 0., 1.0
 -1.0, 0.
*AMPLITUDE, NAME=K12_ux, DEFINITION=PERIODIC
 1, 0.78539816, 0., 0.
 0., 1.0
*AMPLITUDE, NAME=K12_uy, DEFINITION=PERIODIC
 1, 0.78539816, 0., 0.
 0., 1.0
*AMPLITUDE, NAME=K12_uz, DEFINITION=PERIODIC
 1, 0.78539816, 0., 1.0
 -1.0, 0.
*AMPLITUDE, NAME=TVS1_ux, DEFINITION=PERIODIC
 1, 0.78539816, 0., 0.
 0., 1.0
*AMPLITUDE, NAME=TVS1_uz, DEFINITION=PERIODIC
 1, 0.78539816, 0., 1.0
 -1.0, 0.
*AMPLITUDE, NAME=TVS2_uy, DEFINITION=PERIODIC
 1, 0.78539816, 0., 0.
 0., 1.0
*AMPLITUDE, NAME=TVS2_uz, DEFINITION=PERIODIC
 1, 0.78539816, 0., 1.0
 -1.0, 0.
*AMPLITUDE, NAME=CROP_ux, DEFINITION=PERIODIC
 1, 0.78539816, 0., 0.
 0., 1.0
*AMPLITUDE, NAME=CROP_uy, DEFINITION=PERIODIC
 1, 0.78539816, 0., 0.
 0., 1.0
*AMPLITUDE, NAME=CROP_uz, DEFINITION=PERIODIC
 1, 0.78539816, 0., 1.0
 -1.0, 0.
*AMPLITUDE, NAME=HGROT_uy, DEFINITION=PERIODIC
 1, 0.78539816, 0., 0.
 0., 1.0
*AMPLITUDE, NAME=HGROT_uz, DEFINITION=PERIODIC
 1, 0.78539816, 0., 1.0
 -1.0, 0.
**
**
**
*STEP, NLGEOM=YES, INC=20
*STATIC
 0.05, 1.0, ,0.05
**
** Membrane loading e11
**
** SC8R and S4R results are expected to match.
**
*BOUNDARY, OP=NEW
 ASSEMBLY.E11.NSC8R, 2, 3, 0.0
 ASSEMBLY.E11.900000001, 1, 1, -0.8
 ASSEMBLY.E11.900000002, 1, 1, +0.8
 ASSEMBLY.E11.900000003, 1, 1, +0.8
 ASSEMBLY.E11.900000004, 1, 1, -0.8
 ASSEMBLY.E11.900000005, 1, 1, -0.8
 ASSEMBLY.E11.900000006, 1, 1, +0.8
 ASSEMBLY.E11.900000007, 1, 1, +0.8
 ASSEMBLY.E11.900000008, 1, 1, -0.8
*BOUNDARY, OP=NEW
 ASSEMBLY.E11.NS4R,2,6,  0.0
 ASSEMBLY.E11.900000011, 1, 1, -0.8
 ASSEMBLY.E11.900000012, 1, 1, +0.8
 ASSEMBLY.E11.900000013, 1, 1, +0.8
 ASSEMBLY.E11.900000014, 1, 1, -0.8
*BOUNDARY, OP=NEW
 ASSEMBLY.E11.NC3D8R, 1, 3, 0.0
**
** Membrane loading e22
**
** SC8R and S4R results are expected to match.
**
*BOUNDARY, OP=NEW
 ASSEMBLY.E22.NSC8R, 1, 1, 0.0
 ASSEMBLY.E22.NSC8R, 3, 3, 0.0
 ASSEMBLY.E22.900000001, 2, 2, -0.8
 ASSEMBLY.E22.900000002, 2, 2, -0.8
 ASSEMBLY.E22.900000003, 2, 2, +0.8
 ASSEMBLY.E22.900000004, 2, 2, +0.8
 ASSEMBLY.E22.900000005, 2, 2, -0.8
 ASSEMBLY.E22.900000006, 2, 2, -0.8
 ASSEMBLY.E22.900000007, 2, 2, +0.8
 ASSEMBLY.E22.900000008, 2, 2, +0.8
*BOUNDARY, OP=NEW
 ASSEMBLY.E22.NS4R,1,1,  0.0
 ASSEMBLY.E22.NS4R,3,6,  0.0
 ASSEMBLY.E22.900000011, 2, 2, -0.8
 ASSEMBLY.E22.900000012, 2, 2, -0.8
 ASSEMBLY.E22.900000013, 2, 2, +0.8
 ASSEMBLY.E22.900000014, 2, 2, +0.8
*BOUNDARY, OP=NEW
 ASSEMBLY.E22.NC3D8R, 1, 3, 0.0
**
** Membrane loading e12
**
** SC8R and S4R results are expected to match.
**
*BOUNDARY, OP=NEW
 ASSEMBLY.E12.NSC8R, 3, 3, 0.0
 ASSEMBLY.E12.900000001, 2, 2, -0.8
 ASSEMBLY.E12.900000002, 2, 2, +0.8
 ASSEMBLY.E12.900000003, 2, 2, +0.8
 ASSEMBLY.E12.900000004, 2, 2, -0.8
 ASSEMBLY.E12.900000001, 1, 1, -0.8
 ASSEMBLY.E12.900000002, 1, 1, -0.8
 ASSEMBLY.E12.900000003, 1, 1, +0.8
 ASSEMBLY.E12.900000004, 1, 1, +0.8
 ASSEMBLY.E12.900000005, 2, 2, -0.8
 ASSEMBLY.E12.900000006, 2, 2, +0.8
 ASSEMBLY.E12.900000007, 2, 2, +0.8
 ASSEMBLY.E12.900000008, 2, 2, -0.8
 ASSEMBLY.E12.900000005, 1, 1, -0.8
 ASSEMBLY.E12.900000006, 1, 1, -0.8
 ASSEMBLY.E12.900000007, 1, 1, +0.8
 ASSEMBLY.E12.900000008, 1, 1, +0.8
*BOUNDARY, OP=NEW
 ASSEMBLY.E12.NS4R,3,6,  0.0
 ASSEMBLY.E12.900000011, 2, 2, -0.8
 ASSEMBLY.E12.900000012, 2, 2, +0.8
 ASSEMBLY.E12.900000013, 2, 2, +0.8
 ASSEMBLY.E12.900000014, 2, 2, -0.8
 ASSEMBLY.E12.900000011, 1, 1, -0.8
 ASSEMBLY.E12.900000012, 1, 1, -0.8
 ASSEMBLY.E12.900000013, 1, 1, +0.8
 ASSEMBLY.E12.900000014, 1, 1, +0.8
*BOUNDARY, OP=NEW
 ASSEMBLY.E12.NC3D8R, 1, 3, 0.0
**
** Curvature loading K11
**
** SC8R and S4R results are expected to match.
** Note: May want to decrease time steps to converge to
** proper solution.
**
*BOUNDARY, OP=NEW
 ASSEMBLY.K11.NSC8R, 2, 2, 0.0
*BOUNDARY, OP=NEW, AMP=k11_ux
 ASSEMBLY.K11.900000001, 1, 1, -0.05
 ASSEMBLY.K11.900000002, 1, 1, +0.05
 ASSEMBLY.K11.900000003, 1, 1, +0.05
 ASSEMBLY.K11.900000004, 1, 1, -0.05
 ASSEMBLY.K11.900000005, 1, 1, +0.05
 ASSEMBLY.K11.900000006, 1, 1, -0.05
 ASSEMBLY.K11.900000007, 1, 1, -0.05
 ASSEMBLY.K11.900000008, 1, 1, +0.05
*BOUNDARY, OP=NEW, AMP=k11_uz
 ASSEMBLY.K11.900000001, 3, 3, +0.05
 ASSEMBLY.K11.900000002, 3, 3, +0.05
 ASSEMBLY.K11.900000003, 3, 3, +0.05
 ASSEMBLY.K11.900000004, 3, 3, +0.05
 ASSEMBLY.K11.900000005, 3, 3, -0.05
 ASSEMBLY.K11.900000006, 3, 3, -0.05
 ASSEMBLY.K11.900000007, 3, 3, -0.05
 ASSEMBLY.K11.900000008, 3, 3, -0.05
*BOUNDARY, OP=NEW
 ASSEMBLY.K11.NS4R,1,4,  0.0
 ASSEMBLY.K11.NS4R,6,6,  0.0
 ASSEMBLY.K11.900000011, 5, 5, +0.78539816
 ASSEMBLY.K11.900000012, 5, 5, -0.78539816
 ASSEMBLY.K11.900000013, 5, 5, -0.78539816
 ASSEMBLY.K11.900000014, 5, 5, +0.78539816
*BOUNDARY, OP=NEW
 ASSEMBLY.K11.NC3D8R, 1, 3, 0.0
**
** Curvature loading K22
**
** SC8R and S4R results are expected to match.
** Note: May want to decrease time steps to converge to
** proper solution.
**
*BOUNDARY, OP=NEW
 ASSEMBLY.K22.NSC8R, 1, 1, 0.0
*BOUNDARY, OP=NEW, AMP=k22_uy
 ASSEMBLY.K22.900000001, 2, 2, +0.05
 ASSEMBLY.K22.900000002, 2, 2, +0.05
 ASSEMBLY.K22.900000003, 2, 2, -0.05
 ASSEMBLY.K22.900000004, 2, 2, -0.05
 ASSEMBLY.K22.900000005, 2, 2, -0.05
 ASSEMBLY.K22.900000006, 2, 2, -0.05
 ASSEMBLY.K22.900000007, 2, 2, +0.05
 ASSEMBLY.K22.900000008, 2, 2, +0.05
*BOUNDARY, OP=NEW, AMP=k22_uz
 ASSEMBLY.K22.900000001, 3, 3, +0.05
 ASSEMBLY.K22.900000002, 3, 3, +0.05
 ASSEMBLY.K22.900000003, 3, 3, +0.05
 ASSEMBLY.K22.900000004, 3, 3, +0.05
 ASSEMBLY.K22.900000005, 3, 3, -0.05
 ASSEMBLY.K22.900000006, 3, 3, -0.05
 ASSEMBLY.K22.900000007, 3, 3, -0.05
 ASSEMBLY.K22.900000008, 3, 3, -0.05
*BOUNDARY, OP=NEW
 ASSEMBLY.K22.NS4R,1,3,  0.0
 ASSEMBLY.K22.NS4R,5,6,  0.0
 ASSEMBLY.K22.900000011, 4, 4, +0.78539816
 ASSEMBLY.K22.900000012, 4, 4, +0.78539816
 ASSEMBLY.K22.900000013, 4, 4, -0.78539816
 ASSEMBLY.K22.900000014, 4, 4, -0.78539816
*BOUNDARY, OP=NEW
 ASSEMBLY.K22.NC3D8R, 1, 3, 0.0
**
** Curvature loading K12
** Note: May want to decrease time steps to converge to
** proper solution.
**
*BOUNDARY, OP=NEW, AMP=k12_ux
 ASSEMBLY.K12.900000001, 1, 1, +0.0353553
 ASSEMBLY.K12.900000002, 1, 1, +0.0353553
 ASSEMBLY.K12.900000003, 1, 1, -0.0353553
 ASSEMBLY.K12.900000004, 1, 1, -0.0353553
 ASSEMBLY.K12.900000005, 1, 1, -0.0353553
 ASSEMBLY.K12.900000006, 1, 1, -0.0353553
 ASSEMBLY.K12.900000007, 1, 1, +0.0353553
 ASSEMBLY.K12.900000008, 1, 1, +0.0353553
*BOUNDARY, OP=NEW, AMP=k12_uy
 ASSEMBLY.K12.900000001, 2, 2, +0.0353553
 ASSEMBLY.K12.900000002, 2, 2, -0.0353553
 ASSEMBLY.K12.900000003, 2, 2, -0.0353553
 ASSEMBLY.K12.900000004, 2, 2, +0.0353553
 ASSEMBLY.K12.900000005, 2, 2, -0.0353553
 ASSEMBLY.K12.900000006, 2, 2, +0.0353553
 ASSEMBLY.K12.900000007, 2, 2, +0.0353553
 ASSEMBLY.K12.900000008, 2, 2, -0.0353553
*BOUNDARY, OP=NEW, AMP=k12_uz
 ASSEMBLY.K12.900000001, 3, 3, +0.05
 ASSEMBLY.K12.900000002, 3, 3, +0.05
 ASSEMBLY.K12.900000003, 3, 3, +0.05
 ASSEMBLY.K12.900000004, 3, 3, +0.05
 ASSEMBLY.K12.900000005, 3, 3, -0.05
 ASSEMBLY.K12.900000006, 3, 3, -0.05
 ASSEMBLY.K12.900000007, 3, 3, -0.05
 ASSEMBLY.K12.900000008, 3, 3, -0.05
*BOUNDARY, OP=NEW
 ASSEMBLY.K12.NS4R,1,3,  0.0
 ASSEMBLY.K12.NS4R,6,6,  0.0
 ASSEMBLY.K12.900000011, 4, 4, +0.55536036
 ASSEMBLY.K12.900000012, 4, 4, -0.55536036
 ASSEMBLY.K12.900000013, 4, 4, -0.55536036
 ASSEMBLY.K12.900000014, 4, 4, +0.55536036
 ASSEMBLY.K12.900000011, 5, 5, -0.55536036
 ASSEMBLY.K12.900000012, 5, 5, -0.55536036
 ASSEMBLY.K12.900000013, 5, 5, +0.55536036
 ASSEMBLY.K12.900000014, 5, 5, +0.55536036
*BOUNDARY, OP=NEW
 ASSEMBLY.K12.NC3D8R, 1, 3, 0.0
**
**
** Thickness direction loading 
** Solution is expected to match C3D8R.
**
*BOUNDARY, OP=NEW
 ASSEMBLY.T.NSC8R, 1, 2, 0.0
 ASSEMBLY.T.900000001, 3, 3, +0.04
 ASSEMBLY.T.900000002, 3, 3, +0.04
 ASSEMBLY.T.900000003, 3, 3, +0.04
 ASSEMBLY.T.900000004, 3, 3, +0.04
 ASSEMBLY.T.900000005, 3, 3, -0.04
 ASSEMBLY.T.900000006, 3, 3, -0.04
 ASSEMBLY.T.900000007, 3, 3, -0.04
 ASSEMBLY.T.900000008, 3, 3, -0.04
*BOUNDARY, OP=NEW
 ASSEMBLY.T.NS4R, 1, 6, 0.0
*BOUNDARY, OP=NEW
 ASSEMBLY.T.NC3D8R, 1, 2, 0.0
 ASSEMBLY.T.900000021, 3, 3, +0.04
 ASSEMBLY.T.900000022, 3, 3, +0.04
 ASSEMBLY.T.900000023, 3, 3, +0.04
 ASSEMBLY.T.900000024, 3, 3, +0.04
 ASSEMBLY.T.900000025, 3, 3, -0.04
 ASSEMBLY.T.900000026, 3, 3, -0.04
 ASSEMBLY.T.900000027, 3, 3, -0.04
 ASSEMBLY.T.900000028, 3, 3, -0.04
**
**
** Thickness gradient in 1-dir
** Strains and stresses should be zero. Solution
** is not expected match C3D8R. 
**
*BOUNDARY, OP=NEW
 ASSEMBLY.TGRAD1.NSC8R, 1, 2, 0.0
 ASSEMBLY.TGRAD1.900000001, 3, 3, +0.010517
 ASSEMBLY.TGRAD1.900000002, 3, 3, -0.009516
 ASSEMBLY.TGRAD1.900000003, 3, 3, -0.009516
 ASSEMBLY.TGRAD1.900000004, 3, 3, +0.010517
 ASSEMBLY.TGRAD1.900000005, 3, 3, -0.009516
 ASSEMBLY.TGRAD1.900000006, 3, 3, +0.010517
 ASSEMBLY.TGRAD1.900000007, 3, 3, +0.010517
 ASSEMBLY.TGRAD1.900000008, 3, 3, -0.009516
*BOUNDARY, OP=NEW
 ASSEMBLY.TGRAD1.NS4R, 1, 6, 0.0
*BOUNDARY, OP=NEW
 ASSEMBLY.TGRAD1.NC3D8R, 1, 2, 0.0
 ASSEMBLY.TGRAD1.900000021, 3, 3, +0.010517
 ASSEMBLY.TGRAD1.900000022, 3, 3, -0.009516
 ASSEMBLY.TGRAD1.900000023, 3, 3, -0.009516
 ASSEMBLY.TGRAD1.900000024, 3, 3, +0.010517
 ASSEMBLY.TGRAD1.900000025, 3, 3, -0.009516
 ASSEMBLY.TGRAD1.900000026, 3, 3, +0.010517
 ASSEMBLY.TGRAD1.900000027, 3, 3, +0.010517
 ASSEMBLY.TGRAD1.900000028, 3, 3, -0.009516
**
**
** Thickness gradient in 2-dir
** Strains and stresses should be zero. Solution
** is not expected match C3D8R. 
**
*BOUNDARY, OP=NEW
 ASSEMBLY.TGRAD2.NSC8R, 1, 2, 0.0
 ASSEMBLY.TGRAD2.900000001, 3, 3, +0.010517
 ASSEMBLY.TGRAD2.900000002, 3, 3, +0.010517
 ASSEMBLY.TGRAD2.900000003, 3, 3, -0.009516
 ASSEMBLY.TGRAD2.900000004, 3, 3, -0.009516
 ASSEMBLY.TGRAD2.900000005, 3, 3, -0.009516
 ASSEMBLY.TGRAD2.900000006, 3, 3, -0.009516
 ASSEMBLY.TGRAD2.900000007, 3, 3, +0.010517
 ASSEMBLY.TGRAD2.900000008, 3, 3, +0.010517
*BOUNDARY, OP=NEW
 ASSEMBLY.TGRAD2.NS4R, 1, 6, 0.0
*BOUNDARY, OP=NEW
 ASSEMBLY.TGRAD2.NC3D8R, 1, 2, 0.0
 ASSEMBLY.TGRAD2.900000021, 3, 3, +0.010517
 ASSEMBLY.TGRAD2.900000022, 3, 3, +0.010517
 ASSEMBLY.TGRAD2.900000023, 3, 3, -0.009516
 ASSEMBLY.TGRAD2.900000024, 3, 3, -0.009516
 ASSEMBLY.TGRAD2.900000025, 3, 3, -0.009516
 ASSEMBLY.TGRAD2.900000026, 3, 3, -0.009516
 ASSEMBLY.TGRAD2.900000027, 3, 3, +0.010517
 ASSEMBLY.TGRAD2.900000028, 3, 3, +0.010517
**
** Transverse shear in 1-direction
** Results for SC8R are S4R are supposed to match.
**
*BOUNDARY, OP=NEW
 ASSEMBLY.TVS1.NSC8R, 2, 2, 0.0
*BOUNDARY, OP=NEW, AMP=tvs1_ux
 ASSEMBLY.TVS1.900000001, 1, 1, +0.05
 ASSEMBLY.TVS1.900000002, 1, 1, +0.05
 ASSEMBLY.TVS1.900000003, 1, 1, +0.05
 ASSEMBLY.TVS1.900000004, 1, 1, +0.05
 ASSEMBLY.TVS1.900000005, 1, 1, -0.05
 ASSEMBLY.TVS1.900000006, 1, 1, -0.05
 ASSEMBLY.TVS1.900000007, 1, 1, -0.05
 ASSEMBLY.TVS1.900000008, 1, 1, -0.05
*BOUNDARY, OP=NEW, AMP=tvs1_uz
 ASSEMBLY.TVS1.900000001, 3, 3, +0.05
 ASSEMBLY.TVS1.900000002, 3, 3, +0.05
 ASSEMBLY.TVS1.900000003, 3, 3, +0.05
 ASSEMBLY.TVS1.900000004, 3, 3, +0.05
 ASSEMBLY.TVS1.900000005, 3, 3, -0.05
 ASSEMBLY.TVS1.900000006, 3, 3, -0.05
 ASSEMBLY.TVS1.900000007, 3, 3, -0.05
 ASSEMBLY.TVS1.900000008, 3, 3, -0.05
*BOUNDARY, OP=NEW
 ASSEMBLY.TVS1.NS4R,1,4,  0.0
 ASSEMBLY.TVS1.NS4R,6,6,  0.0
 ASSEMBLY.TVS1.900000011, 5, 5, -0.78539816
 ASSEMBLY.TVS1.900000012, 5, 5, -0.78539816
 ASSEMBLY.TVS1.900000013, 5, 5, -0.78539816
 ASSEMBLY.TVS1.900000014, 5, 5, -0.78539816
*BOUNDARY, OP=NEW
 ASSEMBLY.TVS1.NC3D8R, 1, 3, 0.0
**
** Transverse shear in 2-direction
** Results for SC8R are S4R are supposed to match.
**
*BOUNDARY, OP=NEW
 ASSEMBLY.TVS2.NSC8R, 1, 1, 0.0
*BOUNDARY, OP=NEW, AMP=tvs2_uy
 ASSEMBLY.TVS2.900000001, 2, 2, -0.05
 ASSEMBLY.TVS2.900000002, 2, 2, -0.05
 ASSEMBLY.TVS2.900000003, 2, 2, -0.05
 ASSEMBLY.TVS2.900000004, 2, 2, -0.05
 ASSEMBLY.TVS2.900000005, 2, 2, +0.05
 ASSEMBLY.TVS2.900000006, 2, 2, +0.05
 ASSEMBLY.TVS2.900000007, 2, 2, +0.05
 ASSEMBLY.TVS2.900000008, 2, 2, +0.05
*BOUNDARY, OP=NEW, AMP=tvs2_uz
 ASSEMBLY.TVS2.900000001, 3, 3, +0.05
 ASSEMBLY.TVS2.900000002, 3, 3, +0.05
 ASSEMBLY.TVS2.900000003, 3, 3, +0.05
 ASSEMBLY.TVS2.900000004, 3, 3, +0.05
 ASSEMBLY.TVS2.900000005, 3, 3, -0.05
 ASSEMBLY.TVS2.900000006, 3, 3, -0.05
 ASSEMBLY.TVS2.900000007, 3, 3, -0.05
 ASSEMBLY.TVS2.900000008, 3, 3, -0.05
*BOUNDARY, OP=NEW
 ASSEMBLY.TVS2.NS4R,1,3,  0.0
 ASSEMBLY.TVS2.NS4R,5,6,  0.0
 ASSEMBLY.TVS2.900000011, 4, 4, -0.78539816
 ASSEMBLY.TVS2.900000012, 4, 4, -0.78539816
 ASSEMBLY.TVS2.900000013, 4, 4, -0.78539816
 ASSEMBLY.TVS2.900000014, 4, 4, -0.78539816
*BOUNDARY, OP=NEW
 ASSEMBLY.TVS2.NC3D8R, 1, 3, 0.0
**
** Butterfly Hourglass Mode
** Section strains and stresses should be zero. Reaction
** forces on SC8R should be 1/2 of S4R at corresponding 
** edges.
**
*BOUNDARY, OP=NEW
 ASSEMBLY.HGBUTTERFLY.NSC8R, 1, 1, 0.0
 ASSEMBLY.HGBUTTERFLY.NSC8R, 2, 2, 0.0
 ASSEMBLY.HGBUTTERFLY.900000001, 3, 3, +0.1
 ASSEMBLY.HGBUTTERFLY.900000002, 3, 3, -0.1
 ASSEMBLY.HGBUTTERFLY.900000003, 3, 3, +0.1
 ASSEMBLY.HGBUTTERFLY.900000004, 3, 3, -0.1
 ASSEMBLY.HGBUTTERFLY.900000005, 3, 3, +0.1
 ASSEMBLY.HGBUTTERFLY.900000006, 3, 3, -0.1
 ASSEMBLY.HGBUTTERFLY.900000007, 3, 3, +0.1
 ASSEMBLY.HGBUTTERFLY.900000008, 3, 3, -0.1
*BOUNDARY, OP=NEW
 ASSEMBLY.HGBUTTERFLY.NS4R,1,2,  0.0
 ASSEMBLY.HGBUTTERFLY.NS4R,4,6,  0.0
 ASSEMBLY.HGBUTTERFLY.900000011, 3, 3,  0.1
 ASSEMBLY.HGBUTTERFLY.900000012, 3, 3, -0.1
 ASSEMBLY.HGBUTTERFLY.900000013, 3, 3,  0.1
 ASSEMBLY.HGBUTTERFLY.900000014, 3, 3, -0.1
*BOUNDARY, OP=NEW
 ASSEMBLY.HGBUTTERFLY.NC3D8R,1,3,  0.0
**
** Crop circle hourglass mode
** Strains and stress should be zero. Reaction forces
** should correspond to S4R.
**
*BOUNDARY, OP=NEW, AMP=CROP_ux
 ASSEMBLY.HGCROP.900000001, 1, 1, -0.0353553
 ASSEMBLY.HGCROP.900000002, 1, 1, -0.0353553
 ASSEMBLY.HGCROP.900000003, 1, 1, +0.0353553
 ASSEMBLY.HGCROP.900000004, 1, 1, +0.0353553
 ASSEMBLY.HGCROP.900000005, 1, 1, +0.0353553
 ASSEMBLY.HGCROP.900000006, 1, 1, +0.0353553
 ASSEMBLY.HGCROP.900000007, 1, 1, -0.0353553
 ASSEMBLY.HGCROP.900000008, 1, 1, -0.0353553
*BOUNDARY, OP=NEW, AMP=CROP_uy
 ASSEMBLY.HGCROP.900000001, 2, 2, +0.0353553
 ASSEMBLY.HGCROP.900000002, 2, 2, -0.0353553
 ASSEMBLY.HGCROP.900000003, 2, 2, -0.0353553
 ASSEMBLY.HGCROP.900000004, 2, 2, +0.0353553
 ASSEMBLY.HGCROP.900000005, 2, 2, -0.0353553
 ASSEMBLY.HGCROP.900000006, 2, 2, +0.0353553
 ASSEMBLY.HGCROP.900000007, 2, 2, +0.0353553
 ASSEMBLY.HGCROP.900000008, 2, 2, -0.0353553
*BOUNDARY, OP=NEW, AMP=CROP_uz
 ASSEMBLY.HGCROP.900000001, 3, 3, +0.05
 ASSEMBLY.HGCROP.900000002, 3, 3, +0.05
 ASSEMBLY.HGCROP.900000003, 3, 3, +0.05
 ASSEMBLY.HGCROP.900000004, 3, 3, +0.05
 ASSEMBLY.HGCROP.900000005, 3, 3, -0.05
 ASSEMBLY.HGCROP.900000006, 3, 3, -0.05
 ASSEMBLY.HGCROP.900000007, 3, 3, -0.05
 ASSEMBLY.HGCROP.900000008, 3, 3, -0.05
*BOUNDARY, OP=NEW
 ASSEMBLY.HGCROP.NS4R,1,3,  0.0
 ASSEMBLY.HGCROP.NS4R,6,6,  0.0
 ASSEMBLY.HGCROP.900000011, 4, 4, +0.55536036
 ASSEMBLY.HGCROP.900000012, 4, 4, -0.55536036
 ASSEMBLY.HGCROP.900000013, 4, 4, -0.55536036
 ASSEMBLY.HGCROP.900000014, 4, 4, +0.55536036
 ASSEMBLY.HGCROP.900000011, 5, 5, +0.55536036
 ASSEMBLY.HGCROP.900000012, 5, 5, +0.55536036
 ASSEMBLY.HGCROP.900000013, 5, 5, -0.55536036
 ASSEMBLY.HGCROP.900000014, 5, 5, -0.55536036
*BOUNDARY, OP=NEW
 ASSEMBLY.HGCROP.NC3D8R, 1, 3, 0.0
**
** Membrane Hourglass mode top and bottom faces
** Section strains and stresses should be zero. Reaction
** forces on SC8R should be 1/2 of S4R at corresponding 
** edges.
**
*BOUNDARY, OP=NEW
 ASSEMBLY.HGMEMBRANE1.NSC8R, 2, 2, 0.0
 ASSEMBLY.HGMEMBRANE1.NSC8R, 3, 3, 0.0
 ASSEMBLY.HGMEMBRANE1.900000001, 1, 1, +0.5
 ASSEMBLY.HGMEMBRANE1.900000002, 1, 1, -0.5
 ASSEMBLY.HGMEMBRANE1.900000003, 1, 1, +0.5
 ASSEMBLY.HGMEMBRANE1.900000004, 1, 1, -0.5
 ASSEMBLY.HGMEMBRANE1.900000005, 1, 1, +0.5
 ASSEMBLY.HGMEMBRANE1.900000006, 1, 1, -0.5
 ASSEMBLY.HGMEMBRANE1.900000007, 1, 1, +0.5
 ASSEMBLY.HGMEMBRANE1.900000008, 1, 1, -0.5
*BOUNDARY, OP=NEW
 ASSEMBLY.HGMEMBRANE1.NS4R,2, 5, 0.0
 ASSEMBLY.HGMEMBRANE1.900000011, 1, 1, +0.5
 ASSEMBLY.HGMEMBRANE1.900000012, 1, 1, -0.5
 ASSEMBLY.HGMEMBRANE1.900000013, 1, 1, +0.5
 ASSEMBLY.HGMEMBRANE1.900000014, 1, 1, -0.5
*BOUNDARY, OP=NEW
 ASSEMBLY.HGMEMBRANE1.NC3D8R,1,3,  0.0
**
** Membrane Hourglass mode top and bottom faces
** Section strains and stresses should be zero. Reaction
** forces on SC8R should be 1/2 of S4R at corresponding 
** edges.
*BOUNDARY, OP=NEW
 ASSEMBLY.HGMEMBRANE2.NSC8R, 1, 1, 0.0
 ASSEMBLY.HGMEMBRANE2.NSC8R, 3, 3, 0.0
 ASSEMBLY.HGMEMBRANE2.900000001, 2, 2, +0.5
 ASSEMBLY.HGMEMBRANE2.900000002, 2, 2, -0.5
 ASSEMBLY.HGMEMBRANE2.900000003, 2, 2, +0.5
 ASSEMBLY.HGMEMBRANE2.900000004, 2, 2, -0.5
 ASSEMBLY.HGMEMBRANE2.900000005, 2, 2, +0.5
 ASSEMBLY.HGMEMBRANE2.900000006, 2, 2, -0.5
 ASSEMBLY.HGMEMBRANE2.900000007, 2, 2, +0.5
 ASSEMBLY.HGMEMBRANE2.900000008, 2, 2, -0.5
*BOUNDARY, OP=NEW
 ASSEMBLY.HGMEMBRANE2.NS4R,1,1,  0.0
 ASSEMBLY.HGMEMBRANE2.NS4R,3,5,  0.0
 ASSEMBLY.HGMEMBRANE2.900000011, 2, 2, +0.5
 ASSEMBLY.HGMEMBRANE2.900000012, 2, 2, -0.5
 ASSEMBLY.HGMEMBRANE2.900000013, 2, 2, +0.5
 ASSEMBLY.HGMEMBRANE2.900000014, 2, 2, -0.5
*BOUNDARY, OP=NEW
 ASSEMBLY.HGMEMBRANE2.NC3D8R,1,3,  0.0
**
** Rotational hourglass mode
** Section strains and stresses should be zero.
** Results for S4R are given for comparison purpose. Reaction
** forces between S4R and SC8R do not have to match since
** this mode is controlled by to different hourglass modes.
**
*BOUNDARY, OP=NEW
 ASSEMBLY.HGROT.NSC8R, 1, 1, 0.0
*BOUNDARY, OP=NEW, AMP=hgrot_uy
 ASSEMBLY.HGROT.900000001, 2, 2, -0.05
 ASSEMBLY.HGROT.900000002, 2, 2, +0.05
 ASSEMBLY.HGROT.900000003, 2, 2, -0.05
 ASSEMBLY.HGROT.900000004, 2, 2, +0.05
 ASSEMBLY.HGROT.900000005, 2, 2, +0.05
 ASSEMBLY.HGROT.900000006, 2, 2, -0.05
 ASSEMBLY.HGROT.900000007, 2, 2, +0.05
 ASSEMBLY.HGROT.900000008, 2, 2, -0.05
*BOUNDARY, OP=NEW, AMP=hgrot_uz
 ASSEMBLY.HGROT.900000001, 3, 3, +0.05
 ASSEMBLY.HGROT.900000002, 3, 3, +0.05
 ASSEMBLY.HGROT.900000003, 3, 3, +0.05
 ASSEMBLY.HGROT.900000004, 3, 3, +0.05
 ASSEMBLY.HGROT.900000005, 3, 3, -0.05
 ASSEMBLY.HGROT.900000006, 3, 3, -0.05
 ASSEMBLY.HGROT.900000007, 3, 3, -0.05
 ASSEMBLY.HGROT.900000008, 3, 3, -0.05
*BOUNDARY, OP=NEW
 ASSEMBLY.HGROT.NS4R,1,3,  0.0
 ASSEMBLY.HGROT.NS4R,5,6,  0.0
 ASSEMBLY.HGROT.900000011, 4, 4, +0.78539816
 ASSEMBLY.HGROT.900000012, 4, 4, -0.78539816
 ASSEMBLY.HGROT.900000013, 4, 4, +0.78539816
 ASSEMBLY.HGROT.900000014, 4, 4, -0.78539816
*BOUNDARY, OP=NEW
 ASSEMBLY.HGROT.NC3D8R,1,3,  0.0
**
** W-mode Hourglass mode
** Section strains and stresses should be zero.
** Reaction forces should match C3D8R.
*BOUNDARY, OP=NEW
 ASSEMBLY.HGWMODE.NSC8R, 1, 2, 0.0
 ASSEMBLY.HGWMODE.900000001, 3, 3, +0.010517
 ASSEMBLY.HGWMODE.900000002, 3, 3, -0.009516
 ASSEMBLY.HGWMODE.900000003, 3, 3, +0.010517
 ASSEMBLY.HGWMODE.900000004, 3, 3, -0.009516
 ASSEMBLY.HGWMODE.900000005, 3, 3, -0.009516
 ASSEMBLY.HGWMODE.900000006, 3, 3, +0.010517
 ASSEMBLY.HGWMODE.900000007, 3, 3, -0.009516
 ASSEMBLY.HGWMODE.900000008, 3, 3, +0.010517
*BOUNDARY, OP=NEW
 ASSEMBLY.HGWMODE.NS4R, 1, 6, 0.0
*BOUNDARY, OP=NEW
 ASSEMBLY.HGWMODE.NC3D8R, 1, 2, 0.0
 ASSEMBLY.HGWMODE.900000021, 3, 3, +0.01
 ASSEMBLY.HGWMODE.900000022, 3, 3, -0.01
 ASSEMBLY.HGWMODE.900000023, 3, 3, +0.01
 ASSEMBLY.HGWMODE.900000024, 3, 3, -0.01
 ASSEMBLY.HGWMODE.900000025, 3, 3, -0.01
 ASSEMBLY.HGWMODE.900000026, 3, 3, +0.01
 ASSEMBLY.HGWMODE.900000027, 3, 3, -0.01
 ASSEMBLY.HGWMODE.900000028, 3, 3, +0.01
** OUTPUT FOR MEMBRANE 
**
*OUTPUT,FIELD, FREQ=1000
*ELEMENT OUTPUT, ELSET=ASSEMBLY.EOUT_MEMBRANE
 LE,
 S,
 SF,
 SE,
 STH,
*NODE OUTPUT, NSET=NOUT_MEMBRANE
 U,
 RF,
*EL FILE, FREQ=1000, ELSET=ASSEMBLY.EOUT_MEMBRANE
 LE,
 S,
 SF,
 SE,
 STH,
*NODE FILE, FREQ=1000, NSET=NOUT_MEMBRANE
 U,
 RF,
*EL PRINT, FREQ=1000, ELSET=ASSEMBLY.EOUT_MEMBRANE
 LE,
 S,
 SF,
 SE,
 STH,
*NODE PRINT, FREQ=1000, NSET=NOUT_MEMBRANE
 U,
 RF,
**
** OUTPUT FOR CURVATURE 
**
*OUTPUT,FIELD, FREQ=1000
*ELEMENT OUTPUT, ELSET=ASSEMBLY.EOUT_CURVATURE
 LE,
 S,
 SF,
 SE,
 STH,
*NODE OUTPUT, NSET=NOUT_CURVATURE
 U,
 RF,
*EL FILE, FREQ=1000, ELSET=ASSEMBLY.EOUT_CURVATURE
 LE,
 S,
 SF,
 SE,
 STH,
*NODE FILE, FREQ=1000, NSET=NOUT_CURVATURE
 U,
 RF,
*EL PRINT, FREQ=1000, ELSET=ASSEMBLY.EOUT_CURVATURE
 LE,
 S,
 SF,
 SE,
 STH,
*NODE PRINT, FREQ=1000, NSET=NOUT_CURVATURE
 U,
 RF,
**
** OUTPUT FOR SHEAR 
**
*OUTPUT,FIELD, FREQ=1000
*ELEMENT OUTPUT, ELSET=ASSEMBLY.EOUT_TVS
 LE,
 S,
 SF,
 SE,
 STH,
*ELEMENT OUTPUT, ELSET=ASSEMBLY.EOUT_TVS
 3,
 CTSHR,
*NODE OUTPUT, NSET=NOUT_TVS
 U,
 RF,
*EL FILE, FREQ=1000, ELSET=ASSEMBLY.EOUT_TVS
 LE,
 S,
 SF,
 SE,
 STH,
*NODE FILE, FREQ=1000, NSET=NOUT_TVS
 U,
 RF,
*EL PRINT, FREQ=1000, ELSET=ASSEMBLY.EOUT_TVS
 LE,
 S,
 SF,
 SE,
 STH,
*EL PRINT, FREQ=1000, ELSET=ASSEMBLY.EOUT_TVS
 3,
 CTSHR,
 TSHR,
*NODE PRINT, FREQ=1000, NSET=NOUT_TVS
 U,
 RF,
**
** HOURGLASS MODES
**
*OUTPUT,FIELD, FREQ=1000
*ELEMENT OUTPUT, ELSET=ASSEMBLY.EOUT_HGMODES
 LE,
 S,
 SF,
 SE,
 STH,
*NODE OUTPUT, NSET=NOUT_HGMODES
 U,
 RF,
*EL FILE, FREQ=1000, ELSET=ASSEMBLY.EOUT_HGMODES
 LE,
 S,
 SF,
 SE,
 STH,
*NODE FILE, FREQ=1000, NSET=NOUT_HGMODES
 U,
 RF,
*EL PRINT, FREQ=1000, ELSET=ASSEMBLY.EOUT_HGMODES
 LE,
 S,
 SF,
 SE,
 STH,
*NODE PRINT, FREQ=1000, NSET=NOUT_HGMODES
 U,
 RF,
**
** OUTPUT FOR THICKNESS 
**
*OUTPUT,FIELD, FREQ=1000
*ELEMENT OUTPUT, ELSET=ASSEMBLY.T.ELSET_SC8R
 LE,
 S,
 SF,
 SE,
 STH,
*ELEMENT OUTPUT, ELSET=ASSEMBLY.TGRAD1.ELSET_SC8R
 LE,
 S,
 SF,
 SE,
 STH,
*ELEMENT OUTPUT, ELSET=ASSEMBLY.TGRAD2.ELSET_SC8R
 LE,
 S,
 SF,
 SE,
 STH,
*ELEMENT OUTPUT, ELSET=ASSEMBLY.T.ELSET_C3D8R
 LE,
 S,
*NODE OUTPUT, NSET=NOUT_THICK
 U,
 RF,
*EL FILE, FREQ=1000, ELSET=ASSEMBLY.T.ELSET_SC8R
 LE,
 S,
 SF,
 SE,
 STH,
*EL FILE, FREQ=1000, ELSET=ASSEMBLY.TGRAD1.ELSET_SC8R
 LE,
 S,
 SF,
 SE,
 STH,
*EL FILE, FREQ=1000, ELSET=ASSEMBLY.TGRAD2.ELSET_SC8R
 LE,
 S,
 SF,
 SE,
 STH,
*EL FILE, ELSET=ASSEMBLY.T.ELSET_C3D8R
 LE,
 S,
*NODE FILE, FREQ=1000, NSET=NOUT_THICK
 U,
 RF,
*EL PRINT, FREQ=1000, ELSET=ASSEMBLY.T.ELSET_SC8R
 LE,
 S,
 SF,
 SE,
 STH,
*EL PRINT, FREQ=1000, ELSET=ASSEMBLY.TGRAD1.ELSET_SC8R
 LE,
 S,
 SF,
 SE,
 STH,
*EL PRINT, FREQ=1000, ELSET=ASSEMBLY.TGRAD2.ELSET_SC8R
 LE,
 S,
 SF,
 SE,
 STH,
*EL PRINT, ELSET=ASSEMBLY.T.ELSET_C3D8R, FREQ=1000
 LE,
 S,
*NODE PRINT, FREQ=1000, NSET=NOUT_THICK
 U,
 RF,
**
** OUTPUT FOR W-Hourglass Mode 
**
*OUTPUT,FIELD, FREQ=1000
*ELEMENT OUTPUT, ELSET=ASSEMBLY.HGWMODE.ELSET_SC8R
 LE,
 S,
 SF,
 SE,
 STH,
*NODE OUTPUT, NSET=NOUT_HGWMODE
 U,
 RF,
*EL FILE, FREQ=1000, ELSET=ASSEMBLY.HGWMODE.ELSET_SC8R
 LE,
 S,
 SF,
 SE,
 STH,
*EL FILE, ELSET=ASSEMBLY.HGWMODE.ELSET_C3D8R
 LE,
 S,
*NODE FILE, FREQ=1000, NSET=NOUT_HGWMODE
 U,
 RF,
*EL PRINT, FREQ=1000, ELSET=ASSEMBLY.HGWMODE.ELSET_SC8R
 LE,
 S,
 SF,
 SE,
 STH,
*EL PRINT, ELSET=ASSEMBLY.HGWMODE.ELSET_C3D8R, FREQ=1000
 LE,
 S,
*NODE PRINT, FREQ=1000, NSET=NOUT_HGWMODE
 U,
 RF,
*END STEP