*Heading
*Node
*Element, type=CPE4, elset=block1
*Element, type=CPE4, elset=block2
*Elset, Elset = blocks
 block1, block2
*Solid Section, elset=blocks, material=Mat1
 1.,
*Element,Type=B21,Elset=beam11
*Element,Type=B21,Elset=beam12
*Element,Type=B21,Elset=beam21
*Element,Type=B21,Elset=beam22
*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
*Nset, nset=drivenodes
*Material, name=Mat1
*Elastic
 1000., 0.3
*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
*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-shear
*Static
 .5, 1., 1e-05, 1.
*Boundary
 refnodes,1,1,0.
 refnodes,2,2,0.3
*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