*Heading
** testing general traction loads on continuum solid shell elements
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=Part-1
*End Part
**
** ASSEMBLY
**
*Assembly, name=Assembly
***  
*Instance, name=Part-1-2, part=Part-1
         1.1,           0.,           0.
*Node
      1,          0.5,          0.5,           1.
      2,          0.5,         -0.5,           1.
      3,          0.5,          0.5,           0.
      4,          0.5,         -0.5,           0.
      5,         -0.5,          0.5,           1.
      6,         -0.5,         -0.5,           1.
      7,         -0.5,          0.5,           0.
      8,         -0.5,         -0.5,           0.
*Element, type=CSS8
1, 6, 2, 4, 8, 5, 1, 3, 7
*Node
      9,           0.,           0.,          0.5
*Nset, nset=Part-1-2-RefPt_, internal
9, 
** Region: (Section-1:Picked)
*Elset, elset=_PickedSet3, internal
 1,
** Section: Section-1
*Solid Section, elset=_PickedSet3, material=Material-1
*End Instance
*Nset, nset=_PickedSet8, internal, instance=Part-1-2, generate
 1,  8,  1
*Elset, elset=_PickedSet8, internal, instance=Part-1-2
 1,
*Nset, nset=_PickedSet9, internal, instance=Part-1-2
 9,
*Nset, nset=_PickedSet10, internal, instance=Part-1-2
 9,
*Elset, elset=_SurfPosX_S4, internal, instance=Part-1-2
 1,
*Surface, type=ELEMENT, name=SurfPosX
_SurfPosX_S4, S4
*Elset, elset=_SurfNegX_S6, internal, instance=Part-1-2
 1,
*Surface, type=ELEMENT, name=SurfNegX
_SurfNegX_S6, S6
*Elset, elset=_SurfPosY_S2, internal, instance=Part-1-2
 1,
*Surface, type=ELEMENT, name=SurfPosY
_SurfPosY_S2, S2
*Elset, elset=_SurfNegY_S1, internal, instance=Part-1-2
 1,
*Surface, type=ELEMENT, name=SurfNegY
_SurfNegY_S1, S1
*Elset, elset=_SurfPosZ_S3, internal, instance=Part-1-2
 1,
*Surface, type=ELEMENT, name=SurfPosZ
_SurfPosZ_S3, S3
*Elset, elset=_SurfNegZ_S5, internal, instance=Part-1-2
 1,
*Surface, type=ELEMENT, name=SurfNegZ
_SurfNegZ_S5, S5
*Surface, type=NODE, name=_PickedSet8_CNS_, internal
_PickedSet8, 1.
** Constraint: Constraint-2
*Coupling, constraint name=Constraint-2, ref node=_PickedSet9, surface=_PickedSet8_CNS_
*Kinematic
*End Assembly
** 
** MATERIALS
** 
*Material, name=Material-1
*Elastic
1000., 0.3
** 
** BOUNDARY CONDITIONS
** 
** Name: BC-1 Type: Displacement/Rotation
*Boundary
_PickedSet10, 1, 1
_PickedSet10, 2, 2
_PickedSet10, 3, 3
_PickedSet10, 4, 4
_PickedSet10, 5, 5
_PickedSet10, 6, 6
** 
** ----------------------------------------------------------------
** 
** STEP: Step-1
** 
*Step, name=Step-1, perturbation
*Static
*load case, name=trvec_px_100_ds
*Dsload
SurfPosX,trvec, 1, 1, 0, 0
*end load case
*load case, name=trvec_px_100_d
*Dload
_SurfPosX_S4,trvec4, 1, 1, 0, 0
*end load case
**
*load case, name=trvec_px_010_ds
*Dsload
SurfPosX,trvec, 1, 0, 1, 0
*end load case
*load case, name=trvec_px_010_d
*Dload
_SurfPosX_s4,trvec4, 1, 0, 1, 0
*end load case
**
*load case, name=trvec_px_001_ds
*Dsload
SurfPosX,trvec, 1, 0, 0, 1
*end load case
*load case, name=trvec_px_001_d
*Dload
_SurfPosX_s4,trvec4, 1, 0, 0, 1
*end load case
**
*load case, name=trvec_px_011_ds
*Dsload
SurfPosX,trvec, 1.41421356,  0., 1., 1.
*end load case
*load case, name=trvec_px_011_d
*Dload
_SurfPosX_s4,trvec4, 1.41421356,  0., 1., 1.
*end load case
**
*load case, name=trvec_px_111_ds
*Dsload
SurfPosX,trvec, 1.73205081,  1., 1., 1.
*end load case
*load case, name=trvec_px_111_d
*Dload
_SurfPosX_s4,trvec4, 1.73205081,  1., 1., 1.
*end load case
** 
**
*load case, name=trshr_px_100_ds
*Dsload
SurfPosX,trshr, 1, 1, 0, 0
*end load case
*load case, name=trshr_px_100_d
*Dload
_SurfPosX_s4,trshr4, 1, 1, 0, 0
*end load case
**
*load case, name=trshr_px_010_ds
*Dsload
SurfPosX,trshr, 1, 0, 1, 0
*end load case
*load case, name=trshr_px_010_d
*Dload
_SurfPosX_s4,trshr4, 1, 0, 1, 0
*end load case
**
*load case, name=trshr_px_001_ds
*Dsload
SurfPosX,trshr, 1, 0, 0, 1
*end load case
*load case, name=trshr_px_001_d
*Dload
_SurfPosX_s4,trshr4, 1, 0, 0, 1
*end load case
**
*load case, name=trshr_px_011_ds
*Dsload
SurfPosX,trshr, 1.41421356,  0., 1., 1.
*end load case
*load case, name=trshr_px_011_d
*Dload
_SurfPosX_s4,trshr4, 1.41421356,  0., 1., 1.
*end load case
**
*load case, name=trshr_px_111_ds
*Dsload
SurfPosX,trshr, 1.41421356, 1., 1.,1.
*end load case
*load case, name=trshr_px_111_d
*Dload
_SurfPosX_s4,trshr4, 1.41421356, 1., 1.,1.
*end load case
**
**
** FIELD OUTPUT: F-Output-1
** 
*Output, field, variable=PRESELECT
** 
*Node print, nset=_PickedSet10
RF,
*OUTPUT,FIELD
*ELEMENT OUTPUT
 TRNOR,TRSHR
*End Step
** 
** STEP: Step-2
** 
*Step, name=Step-2, perturbation
*Static
*load case, name=trvec_nx_100_ds
*Dsload
SurfNegX,trvec, 1, 1, 0, 0
*end load case
*load case, name=trvec_nx_100_d
*Dload
_SurfNegX_S6,trvec6, 1, 1, 0, 0
*end load case
**
*load case, name=trvec_nx_010_ds
*Dsload
SurfNegX,trvec, 1, 0, 1, 0
*end load case
*load case, name=trvec_nx_010_d
*Dload
_SurfNegX_s6,trvec6, 1, 0, 1, 0
*end load case
**
*load case, name=trvec_nx_001_ds
*Dsload
SurfNegX,trvec, 1, 0, 0, 1
*end load case
*load case, name=trvec_nx_001_d
*Dload
_SurfNegX_s6,trvec6, 1, 0, 0, 1
*end load case
**
*load case, name=trvec_nx_011_ds
*Dsload
SurfNegX,trvec, 1.41421356,  0., 1., 1.
*end load case
*load case, name=trvec_nx_011_d
*Dload
_SurfNegX_s6,trvec6, 1.41421356,  0., 1., 1.
*end load case
**
*load case, name=trvec_nx_111_ds
*Dsload
SurfNegX,trvec, 1.73205081,  1., 1., 1.
*end load case
*load case, name=trvec_nx_111_d
*Dload
_SurfNegX_s6,trvec6, 1.73205081,  1., 1., 1.
*end load case
** 
**
*load case, name=trshr_nx_100_ds
*Dsload
SurfNegX,trshr, 1, 1, 0, 0
*end load case
*load case, name=trshr_nx_100_d
*Dload
_SurfNegX_s6,trshr6, 1, 1, 0, 0
*end load case
**
*load case, name=trshr_nx_010_ds
*Dsload
SurfNegX,trshr, 1, 0, 1, 0
*end load case
*load case, name=trshr_nx_010_d
*Dload
_SurfNegX_s6,trshr6, 1, 0, 1, 0
*end load case
**
*load case, name=trshr_nx_001_ds
*Dsload
SurfNegX,trshr, 1, 0, 0, 1
*end load case
*load case, name=trshr_nx_001_d
*Dload
_SurfNegX_s6,trshr6, 1, 0, 0, 1
*end load case
**
*load case, name=trshr_nx_011_ds
*Dsload
SurfNegX,trshr, 1.41421356,  0., 1., 1.
*end load case
*load case, name=trshr_nx_011_d
*Dload
_SurfNegX_s6,trshr6, 1.41421356,  0., 1., 1.
*end load case
**
*load case, name=trshr_nx_111_ds
*Dsload
SurfNegX,trshr, 1.41421356, 1., 1.,1.
*end load case
*load case, name=trshr_nx_111_d
*Dload
_SurfNegX_s6,trshr6, 1.41421356, 1., 1.,1.
*end load case
*OUTPUT,FIELD
*ELEMENT OUTPUT
 TRNOR,TRSHR
*End Step
** 
** STEP: Step-3
** 
*Step, name=Step-3, perturbation
*Static
*load case, name=trvec_py_100_ds
*Dsload
SurfPosY,trvec, 1, 1, 0, 0
*end load case
*load case, name=trvec_py_100_d
*Dload
_SurfPosY_S2,trvec2, 1, 1, 0, 0
*end load case
**
*load case, name=trvec_py_010_ds
*Dsload
SurfPosY,trvec, 1, 0, 1, 0
*end load case
*load case, name=trvec_py_010_d
*Dload
_SurfPosY_s2,trvec2, 1, 0, 1, 0
*end load case
**
*load case, name=trvec_py_001_ds
*Dsload
SurfPosY,trvec, 1, 0, 0, 1
*end load case
*load case, name=trvec_py_001_d
*Dload
_SurfPosY_s2,trvec2, 1, 0, 0, 1
*end load case
**
*load case, name=trvec_py_101_ds
*Dsload
SurfPosY,trvec, 1.41421356,  1., 0., 1.
*end load case
*load case, name=trvec_py_101_d
*Dload
_SurfPosY_s2,trvec2, 1.41421356,  1., 0., 1.
*end load case
**
*load case, name=trvec_py_111_ds
*Dsload
SurfPosY,trvec, 1.73205081,  1., 1., 1.
*end load case
*load case, name=trvec_py_111_d
*Dload
_SurfPosY_s2,trvec2, 1.73205081,  1., 1., 1.
*end load case
** 
**
*load case, name=trshr_py_100_ds
*Dsload
SurfPosY,trshr, 1, 1, 0, 0
*end load case
*load case, name=trshr_py_100_d
*Dload
_SurfPosY_s2,trshr2, 1, 1, 0, 0
*end load case
**
*load case, name=trshr_py_010_ds
*Dsload
SurfPosY,trshr, 1, 0, 1, 0
*end load case
*load case, name=trshr_py_010_d
*Dload
_SurfPosY_s2,trshr2, 1, 0, 1, 0
*end load case
**
*load case, name=trshr_py_001_ds
*Dsload
SurfPosY,trshr, 1, 0, 0, 1
*end load case
*load case, name=trshr_py_001_d
*Dload
_SurfPosY_s2,trshr2, 1, 0, 0, 1
*end load case
**
*load case, name=trshr_py_101_ds
*Dsload
SurfPosY,trshr, 1.41421356,  1., 0., 1.
*end load case
*load case, name=trshr_py_101_d
*Dload
_SurfPosY_s2,trshr2, 1.41421356,  1., 0., 1.
*end load case
**
*load case, name=trshr_py_111_ds
*Dsload
SurfPosY,trshr, 1.41421356, 1., 1.,1.
*end load case
*load case, name=trshr_py_111_d
*Dload
_SurfPosY_s2,trshr2, 1.41421356, 1., 1.,1.
*end load case
**
*OUTPUT,FIELD
*ELEMENT OUTPUT
 TRNOR,TRSHR
*End Step
** 
** STEP: Step-4
** 
*Step, name=Step-4, perturbation
*Static
*load case, name=trvec_ny_100_ds
*Dsload
SurfNegY,trvec, 1, 1, 0, 0
*end load case
*load case, name=trvec_ny_100_d
*Dload
_SurfNegY_S1,trvec1, 1, 1, 0, 0
*end load case
**
*load case, name=trvec_ny_010_ds
*Dsload
SurfNegY,trvec, 1, 0, 1, 0
*end load case
*load case, name=trvec_ny_010_d
*Dload
_SurfNegY_s1,trvec1, 1, 0, 1, 0
*end load case
**
*load case, name=trvec_ny_001_ds
*Dsload
SurfNegY,trvec, 1, 0, 0, 1
*end load case
*load case, name=trvec_ny_001_d
*Dload
_SurfNegY_s1,trvec1, 1, 0, 0, 1
*end load case
**
*load case, name=trvec_ny_101_ds
*Dsload
SurfNegY,trvec, 1.41421356,  1., 0., 1.
*end load case
*load case, name=trvec_ny_101_d
*Dload
_SurfNegY_s1,trvec1, 1.41421356,  1., 0., 1.
*end load case
**
*load case, name=trvec_ny_111_ds
*Dsload
SurfNegY,trvec, 1.73205081,  1., 1., 1.
*end load case
*load case, name=trvec_ny_111_d
*Dload
_SurfNegY_s1,trvec1, 1.73205081,  1., 1., 1.
*end load case
** 
**
*load case, name=trshr_ny_100_ds
*Dsload
SurfNegY,trshr, 1, 1, 0, 0
*end load case
*load case, name=trshr_ny_100_d
*Dload
_SurfNegY_s1,trshr1, 1, 1, 0, 0
*end load case
**
*load case, name=trshr_ny_010_ds
*Dsload
SurfNegY,trshr, 1, 0, 1, 0
*end load case
*load case, name=trshr_ny_010_d
*Dload
_SurfNegY_s1,trshr1, 1, 0, 1, 0
*end load case
**
*load case, name=trshr_ny_001_ds
*Dsload
SurfNegY,trshr, 1, 0, 0, 1
*end load case
*load case, name=trshr_ny_001_d
*Dload
_SurfNegY_s1,trshr1, 1, 0, 0, 1
*end load case
**
*load case, name=trshr_ny_101_ds
*Dsload
SurfNegY,trshr, 1.41421356,  1., 0., 1.
*end load case
*load case, name=trshr_ny_101_d
*Dload
_SurfNegY_s1,trshr1, 1.41421356,  1., 0., 1.
*end load case
**
*load case, name=trshr_ny_111_ds
*Dsload
SurfNegY,trshr, 1.41421356, 1., 1.,1.
*end load case
*load case, name=trshr_ny_111_d
*Dload
_SurfNegY_s1,trshr1, 1.41421356, 1., 1.,1.
*end load case
**
*OUTPUT,FIELD
*ELEMENT OUTPUT
 TRNOR,TRSHR
*End Step
** 
** STEP: Step-5
** 
*Step, name=Step-5, perturbation
*Static
*load case, name=trvec_pz_100_ds
*Dsload
SurfPosZ,trvec, 1, 1, 0, 0
*end load case
*load case, name=trvec_pz_100_d
*Dload
_SurfPosZ_S3,trvec3, 1, 1, 0, 0
*end load case
**
*load case, name=trvec_pz_010_ds
*Dsload
SurfPosZ,trvec, 1, 0, 1, 0
*end load case
*load case, name=trvec_pz_010_d
*Dload
_SurfPosZ_s3,trvec3, 1, 0, 1, 0
*end load case
**
*load case, name=trvec_pz_001_ds
*Dsload
SurfPosZ,trvec, 1, 0, 0, 1
*end load case
*load case, name=trvec_pz_001_d
*Dload
_SurfPosZ_s3,trvec3, 1, 0, 0, 1
*end load case
**
*load case, name=trvec_pz_110_ds
*Dsload
SurfPosZ,trvec, 1.41421356,  1., 1., 0.
*end load case
*load case, name=trvec_pz_110_d
*Dload
_SurfPosZ_s3,trvec3, 1.41421356,  1., 1., 0.
*end load case
**
*load case, name=trvec_pz_111_ds
*Dsload
SurfPosZ,trvec, 1.73205081,  1., 1., 1.
*end load case
*load case, name=trvec_pz_111_d
*Dload
_SurfPosZ_s3,trvec3, 1.73205081,  1., 1., 1.
*end load case
** 
**
*load case, name=trshr_pz_100_ds
*Dsload
SurfPosZ,trshr, 1, 1, 0, 0
*end load case
*load case, name=trshr_pz_100_d
*Dload
_SurfPosZ_s3,trshr3, 1, 1, 0, 0
*end load case
**
*load case, name=trshr_pz_010_ds
*Dsload
SurfPosZ,trshr, 1, 0, 1, 0
*end load case
*load case, name=trshr_pz_010_d
*Dload
_SurfPosZ_s3,trshr3, 1, 0, 1, 0
*end load case
**
*load case, name=trshr_pz_001_ds
*Dsload
SurfPosZ,trshr, 1, 0, 0, 1
*end load case
*load case, name=trshr_pz_001_d
*Dload
_SurfPosZ_s3,trshr3, 1, 0, 0, 1
*end load case
**
*load case, name=trshr_pz_110_ds
*Dsload
SurfPosZ,trshr, 1.41421356,  1., 1., 0.
*end load case
*load case, name=trshr_pz_110_d
*Dload
_SurfPosZ_s3,trshr3, 1.41421356,  1., 1., 0.
*end load case
**
*load case, name=trshr_pz_111_ds
*Dsload
SurfPosZ,trshr, 1.41421356, 1., 1.,1.
*end load case
*load case, name=trshr_pz_111_d
*Dload
_SurfPosZ_s3,trshr3, 1.41421356, 1., 1.,1.
*end load case
**
*OUTPUT,FIELD
*ELEMENT OUTPUT
 TRNOR,TRSHR
*End Step
** 
** STEP: Step-6
** 
*Step, name=Step-6, perturbation
*Static
*load case, name=trvec_nz_100_ds
*Dsload
SurfNegZ,trvec, 1, 1, 0, 0
*end load case
*load case, name=trvec_nz_100_d
*Dload
_SurfNegZ_S5,trvec5, 1, 1, 0, 0
*end load case
**
*load case, name=trvec_nz_010_ds
*Dsload
SurfNegZ,trvec, 1, 0, 1, 0
*end load case
*load case, name=trvec_nz_010_d
*Dload
_SurfNegZ_s5,trvec5, 1, 0, 1, 0
*end load case
**
*load case, name=trvec_nz_001_ds
*Dsload
SurfNegZ,trvec, 1, 0, 0, 1
*end load case
*load case, name=trvec_nz_001_d
*Dload
_SurfNegZ_s5,trvec5, 1, 0, 0, 1
*end load case
**
*load case, name=trvec_nz_110_ds
*Dsload
SurfNegZ,trvec, 1.41421356,  1., 1., 0.
*end load case
*load case, name=trvec_nz_110_d
*Dload
_SurfNegZ_s5,trvec5, 1.41421356,  1., 1., 0.
*end load case
**
*load case, name=trvec_nz_111_ds
*Dsload
SurfNegZ,trvec, 1.73205081,  1., 1., 1.
*end load case
*load case, name=trvec_nz_111_d
*Dload
_SurfNegZ_s5,trvec5, 1.73205081,  1., 1., 1.
*end load case
** 
**
*load case, name=trshr_nz_100_ds
*Dsload
SurfNegZ,trshr, 1, 1, 0, 0
*end load case
*load case, name=trshr_nz_100_d
*Dload
_SurfNegZ_s5,trshr5, 1, 1, 0, 0
*end load case
**
*load case, name=trshr_nz_010_ds
*Dsload
SurfNegZ,trshr, 1, 0, 1, 0
*end load case
*load case, name=trshr_nz_010_d
*Dload
_SurfNegZ_s5,trshr5, 1, 0, 1, 0
*end load case
**
*load case, name=trshr_nz_001_ds
*Dsload
SurfNegZ,trshr, 1, 0, 0, 1
*end load case
*load case, name=trshr_nz_001_d
*Dload
_SurfNegZ_s5,trshr5, 1, 0, 0, 1
*end load case
**
*load case, name=trshr_nz_110_ds
*Dsload
SurfNegZ,trshr, 1.41421356,  1., 1., 0.
*end load case
*load case, name=trshr_nz_110_d
*Dload
_SurfNegZ_s5,trshr5, 1.41421356,  1., 1., 0.
*end load case
**
*load case, name=trshr_nz_111_ds
*Dsload
SurfNegZ,trshr, 1.41421356, 1., 1.,1.
*end load case
*load case, name=trshr_nz_111_d
*Dload
_SurfNegZ_s5,trshr5, 1.41421356, 1., 1.,1.
*end load case
**
*OUTPUT,FIELD
*ELEMENT OUTPUT
 TRNOR,TRSHR
*End Step