*Heading
*Node
      1,  0., 0.
      2,  1., 0.
      3,  1., 1.
      4,  0., 1.
     11,  2., 0.
     12,  3., 0.
     13,  3., 1.
     14,  2., 1.
    101, 0.5, 1.
    102, 0.5, 2.
    103,  0., 2.
    104,  1., 2.
    111, 2.5, 1.
    112, 2.5, 2.
    113,  2., 2.
    114,  3., 2.
   1000,  0., 0.
*Element, type=CPE4, elset=block1
  1,    1,    2,    3,    4
*Element, type=CPE4, elset=block2
  2,   11,   12,   13,   14
*Elset, Elset = blocks
 block1, block2
*Solid Section, elset=blocks, material=Mat1
 1.,
*Element,Type=B21,Elset=beam11
 101, 101, 102
*Element,Type=B21,Elset=beam12
 102, 102, 103
 103, 102, 104
*Element,Type=B21,Elset=beam21
 111, 111, 112
*Element,Type=B21,Elset=beam22
 112, 112, 113
 113, 112, 114
*Elset, Elset=beam1
 beam11,beam21
*Elset, Elset=beam2
 beam12,beam22
*Coupling,constraint name=beam2block1,surface=top1,ref node=101
*Distributing			      		   
 1,2
 6,6				      		   
*Coupling,constraint name=beam2block2,surface=top2,ref node=111
*Distributing			      		   
 1,2
*Beam Section, Material=Mat2,Section=Circ,Elset=beam1
 1.0, 
 0., 0., -1.
*Beam Section, Material=Mat2,Section=Circ,Elset=beam2
 1.0, 
 0., 0., -1.
*Surface, type=ELEMENT, name=top1
 block1, S3
*Surface, type=ELEMENT, name=top2
 block2, S3
*Surface, type=ELEMENT, name=bot
 blocks, S1
*Surface, type=segments, name=plate
 start, 0., 0.
 line,  3., 0.
*Rigid body, analytical surface=plate, ref node=1000
*Tie, name=blocksToPlate
 bot, plate
*Nset, nset=refnodes
 101, 111
*Nset, nset=drivenodes
 102,112
*Material, name=Mat1
***Elastic
** 1000., 0.3
*Hyperelastic, neo hooke
 1000., .1
*Material, name=Mat2
*Elastic
 30.e6, 0.3
** ----------------------------------------------------------------
** 
** STEP: Step-1
** 
*Step, name=Step-1, perturbation
 perturbation analyses
*Static
*Boundary
 1000,1,2 
 1000,6,6
*Output, field, freq=1
*Node Output
 U,
*El Print, freq=9999
*Node Print, freq=9999
**
** LOAD CASES
**
*Load Case,name=shear_x
*Boundary
 refnodes,1,1,0.3
 refnodes,6,6
*End Load Case
*Load Case,name=pull_y
*Boundary
 refnodes,2,2,0.3
 refnodes,6,6
*End Load Case
*Load Case,name=twist_z
*Boundary
 refnodes,6,6,0.3
*End Load Case
*Load Case,name=all
*Boundary
 refnodes,1,2,0.3
 refnodes,6,6,0.3
*End Load Case
*End Step
** ----------------------------------------------------------------
** 
** STEP: Step-2
** 
*Step, name=Step-2, nlgeom
 x-shear 
*Static
 .5, 1., 1e-05, 1.
*Boundary
 1000, 1, 2
 1000, 6, 6
 refnodes,1,1,0.3
 refnodes,6,6
*Output, field, freq=1
*Node Output
 U,
*El Print, freq=0
*Node Print, freq=0
*End Step
** ----------------------------------------------------------------
** 
** STEP: Step-3
** 
*Step, name=Step-3, nlgeom
 y-pull
*Static
 .5, 1., 1e-05, 1.
*Boundary
 refnodes,1,1,0.
 refnodes,2,2,0.3
 refnodes,6,6
*End Step
** ----------------------------------------------------------------
** 
** STEP: Step-4
** 
*Step, name=Step-4, nlgeom
 z-rotate
*Static
 .5, 1., 1e-05, 1.
*Boundary
 refnodes,6,6,0.3
*End Step
** ----------------------------------------------------------------
** 
** STEP: Step-5
** 
*Step, name=Step-5, nlgeom
 z-rotate back
*Static
 .5, 1., 1e-05, 1.
*Boundary
 refnodes,1,2,0.
 refnodes,6,6,0.
*End Step
** ----------------------------------------------------------------
** 
** STEP: Step-6
** 
*Step, name=Step-6, nlgeom
 all
*Static
 .5, 1., 1e-05, 1.
*Boundary
 refnodes,1,2,0.3
 refnodes,6,6,0.3
*End Step