*Heading
** testing general tractions on 3d rigid 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-1, part=Part-1
*Node
      1,         -0.5,          0.5,           0.
      2,         -0.5,         -0.5,           0.
      3,          0.5,          0.5,           0.
      4,          0.5,         -0.5,           0.
*Element, type=R3D4
1, 1, 2, 4, 3
*Node
      5,           0.,           0.,           0.
*Nset, nset=Part-1-1-RefPt_, internal
5, 
*End Instance
**  
*Instance, name=Part-1-2, part=Part-1
          2.,           0.,           0.
*Node
      1,         -0.5,          0.5,           0.
      2,         -0.5,         -0.5,           0.
      3,          0.5,         -0.5,           0.
      4,          0.5,          0.5,           0.
*Element, type=R3D3
1, 2, 3, 1
2, 1, 3, 4
*Node
      5,           0.,           0.,           0.
*Nset, nset=Part-1-2-RefPt_, internal
5, 
*End Instance
*Nset, nset=_PickedSet6, internal, instance=Part-1-1
 5,
*Nset, nset=_PickedSet7, internal, instance=Part-1-1, generate
 1,  4,  1
*Elset, elset=_PickedSet7, internal, instance=Part-1-1
 1,
*Nset, nset=_PickedSet8, internal, instance=Part-1-2
 5,
*Nset, nset=_PickedSet9, internal, instance=Part-1-2, generate
 1,  4,  1
*Elset, elset=_PickedSet9, internal, instance=Part-1-2
 1, 2
*Nset, nset=_PickedSet10, internal, instance=Part-1-1
 5,
*Nset, nset=_PickedSet10, internal, instance=Part-1-2
 5,
*Elset, elset=_SurfPos_SPOS, internal, instance=Part-1-1
 1,
*Elset, elset=_SurfPos_SPOS, internal, instance=Part-1-2
 1, 2
*Surface, type=ELEMENT, name=SurfPos
_SurfPos_SPOS, SPOS
*Elset, elset=_SurfNeg_SNEG, internal, instance=Part-1-1
 1,
*Elset, elset=_SurfNeg_SNEG, internal, instance=Part-1-2
 1, 2
*Surface, type=ELEMENT, name=SurfNeg
_SurfNeg_SNEG, SNEG
** Constraint: Constraint-1
*Rigid Body, ref node=_PickedSet6, elset=_PickedSet7
** Constraint: Constraint-2
*Rigid Body, ref node=_PickedSet8, elset=_PickedSet9
*End Assembly
** ----------------------------------------------------------------
** 
** 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
** 
** 
** OUTPUT REQUESTS
** FIELD OUTPUT: F-Output-1
** 
*Output, field, variable=PRESELECT
** 
*Node Print, freq=999999, nset=_PickedSet10
RF,
**
**
*Load case, name=trvec100_pos_ds
*dsload
surfPos, trvec, 1, 1, 0, 0
*End load case
*Load case, name=trvec100_pos_d
*dload
_SurfPos_SPOS, trvec, 1, 1, 0, 0
*End load case
**
*Load case, name=trvec010_pos_ds
*dsload
surfPos, trvec, 1, 0, 1, 0
*End load case
*Load case, name=trvec010_pos_d
*dload
_SurfPos_SPOS, trvec, 1, 0, 1, 0
*End load case
**
*Load case, name=trvec001_pos_ds
*dsload
surfPos, trvec, 1, 0, 0, 1
*End load case
*Load case, name=trvec001_pos_d
*dload
_SurfPos_SPOS, trvec, 1, 0, 0, 1
*End load case
**
**
*Load case, name=trshr100_pos_ds
*dsload
surfPos, trshr, 1, 1, 0, 0
*End load case
*Load case, name=trshr100_pos_d
*dload
_SurfPos_SPOS, trshr, 1, 1, 0, 0
*End load case
**
*Load case, name=trshr010_pos_ds
*dsload
surfPos, trshr, 1, 0, 1, 0
*End load case
*Load case, name=trshr010_pos_d
*dload
_SurfPos_SPOS, trshr, 1, 0, 1, 0
*End load case
**
*Load case, name=trshr001_pos_ds
*dsload
surfPos, trshr, 1, 0, 0, 1
*End load case
*Load case, name=trshr001_pos_d
*dload
_SurfPos_SPOS, trshr, 1, 0, 0, 1
*End load case
**
*OUTPUT,FIELD
*ELEMENT OUTPUT
 TRNOR,TRSHR
*End Step
**
**
**
*Step, name=Step-2, pert
*Static
**
**
*Load case, name=trvec100_neg_ds
*dsload
surfNeg, trvec, 1, 1, 0, 0
*End load case
*Load case, name=trvec100_neg_d
*dload
_SurfNeg_SNEG, trvec, 1, 1, 0, 0
*End load case
**
*Load case, name=trvec010_neg_ds
*dsload
surfNeg, trvec, 1, 0, 1, 0
*End load case
*Load case, name=trvec010_neg_d
*dload
_SurfNeg_SNEG, trvec, 1, 0, 1, 0
*End load case
**
*Load case, name=trvec001_neg_ds
*dsload
surfNeg, trvec, 1, 0, 0, 1
*End load case
*Load case, name=trvec001_neg_d
*dload
_SurfNeg_SNEG, trvec, 1, 0, 0, 1
*End load case
**
**
*Load case, name=trshr100_neg_ds
*dsload
surfNeg, trshr, 1, 1, 0, 0
*End load case
*Load case, name=trshr100_neg_d
*dload
_SurfNeg_SNEG, trshr, 1, 1, 0, 0
*End load case
**
*Load case, name=trshr010_neg_ds
*dsload
surfNeg, trshr, 1, 0, 1, 0
*End load case
*Load case, name=trshr010_neg_d
*dload
_SurfNeg_SNEG, trshr, 1, 0, 1, 0
*End load case
**
*Load case, name=trshr001_neg_ds
*dsload
surfNeg, trshr, 1, 0, 0, 1
*End load case
*Load case, name=trshr001_neg_d
*dload
_SurfNeg_SNEG, trshr, 1, 0, 0, 1
*End load case
**
*OUTPUT,FIELD
*ELEMENT OUTPUT
 TRNOR,TRSHR
*End Step