*Heading
*Node
      1, -1.,   0., -1.
      2,  1.,   0., -1.
      3,  1.,   2., -1.
      4, -1.,   2., -1.
      5, -1.,   0.,  1.
      6,  1.,   0.,  1.
      7,  1.,   2.,  1.
      8, -1.,   2.,  1.
     11, -1.,   4., -1.
     12,  1.,   4., -1.
     13,  1.,   6., -1.
     14, -1.,   6., -1.
     15, -1.,   4.,  1.
     16,  1.,   4.,  1.
     17,  1.,   6.,  1.
     18, -1.,   6.,  1.
     21, -1.,   8., -1.
     22,  1.,   8., -1.
     23,  1.,  10., -1.
     24, -1.,  10., -1.
     25, -1.,   8.,  1.
     26,  1.,   8.,  1.
     27,  1.,  10.,  1.
     28, -1.,  10.,  1.
     31, -1.,  12., -1.
     32,  1.,  12., -1.
     33,  1.,  14., -1.
     34, -1.,  14., -1.
     35, -1.,  12.,  1.
     36,  1.,  12.,  1.
     37,  1.,  14.,  1.
     38, -1.,  14.,  1.
     41,  3.,   0., -1.
     42,  5.,   0., -1.
     43,  5.,   2., -1.
     44,  3.,   2., -1.
     45,  3.,   0.,  1.
     46,  5.,   0.,  1.
     47,  5.,   2.,  1.
     48,  3.,   2.,  1.
     51,  3.,   4., -1.
     52,  5.,   4., -1.
     53,  5.,   6., -1.
     54,  3.,   6., -1.
     55,  3.,   4.,  1.
     56,  5.,   4.,  1.
     57,  5.,   6.,  1.
     58,  3.,   6.,  1.
     61,  3.,   8., -1.
     62,  5.,   8., -1.
     63,  5.,  10., -1.
     64,  3.,  10., -1.
     65,  3.,   8.,  1.
     66,  5.,   8.,  1.
     67,  5.,  10.,  1.
     68,  3.,  10.,  1.
     71,  3.,  12., -1.
     72,  5.,  12., -1.
     73,  5.,  14., -1.
     74,  3.,  14., -1.
     75,  3.,  12.,  1.
     76,  5.,  12.,  1.
     77,  5.,  14.,  1.
     78,  3.,  14.,  1.
    101,  0.,   1.,  1.
    102,  0.,   1.,  2.
    103,  0.,   0.,  2.
    104,  0.,   2.,  2.
    111,  0.,   5.,  1.
    112,  0.,   5.,  2.
    113,  0.,   4.,  2.
    114,  0.,   6.,  2.
    121,  0.,   9.,  1.
    122,  0.,   9.,  2.
    123,  0.,   8.,  2.
    124,  0.,  10.,  2.
    131,  0.,  13.,  1.
    132,  0.,  13.,  2.
    133,  0.,  12.,  2.
    134,  0.,  14.,  2.
    141,  4.,   1.,  1.
    142,  4.,   1.,  2.
    143,  4.,   0.,  2.
    144,  4.,   2.,  2.
    151,  4.,   5.,  1.
    152,  4.,   5.,  2.
    153,  4.,   4.,  2.
    154,  4.,   6.,  2.
    161,  4.,   9.,  1.
    162,  4.,   9.,  2.
    163,  4.,   8.,  2.
    164,  4.,  10.,  2.
    171,  4.,  13.,  1.
    172,  4.,  13.,  2.
    173,  4.,  12.,  2.
    174,  4.,  14.,  2.
   1000,  0.,   0., -1.
*Element, type=C3D8, elset=block1
  1,    1,    2,    3,    4,    5,    6,    7,    8
*Element, type=C3D8, elset=block2
  2,   11,   12,   13,   14,   15,   16,   17,   18
*Element, type=C3D8, elset=block3
  3,   21,   22,   23,   24,   25,   26,   27,   28
*Element, type=C3D8, elset=block4
  4,   31,   32,   33,   34,   35,   36,   37,   38
*Element, type=C3D8, elset=block5
  5,   41,   42,   43,   44,   45,   46,   47,   48
*Element, type=C3D8, elset=block6
  6,   51,   52,   53,   54,   55,   56,   57,   58   
*Element, type=C3D8, elset=block7
  7,   61,   62,   63,   64,   65,   66,   67,   68  
*Element, type=C3D8, elset=block8
  8,   71,   72,   73,   74,   75,   76,   77,   78  
*Elset, Elset = blocks
 block1, block2, block3, block4,
 block5, block6, block7, block8
*Solid Section, elset=blocks, material=Mat1
 1.,
*Element,Type=B31,Elset=beam11
 101, 101, 102
*Element,Type=B31,Elset=beam12
 102, 102, 103
 103, 102, 104
*Element,Type=B31,Elset=beam21
 111, 111, 112
*Element,Type=B31,Elset=beam22
 112, 112, 113
 113, 112, 114
*Element,Type=B31,Elset=beam31
 121, 121, 122
*Element,Type=B31,Elset=beam32
 122, 122, 123
 123, 122, 124
*Element,Type=B31,Elset=beam41
 131, 131, 132
*Element,Type=B31,Elset=beam42
 132, 132, 133
 133, 132, 134
*Element,Type=B31,Elset=beam51
 141, 141, 142
*Element,Type=B31,Elset=beam52
 142, 142, 143
 143, 142, 144
*Element,Type=B31,Elset=beam61
 151, 151, 152
*Element,Type=B31,Elset=beam62
 152, 152, 153
 153, 152, 154
*Element,Type=B31,Elset=beam71
 161, 161, 162
*Element,Type=B31,Elset=beam72
 162, 162, 163
 163, 162, 164
*Element,Type=B31,Elset=beam81
 171, 171, 172
*Element,Type=B31,Elset=beam82
 172, 172, 173
 173, 172, 174
*Elset, Elset=beam1
 beam11,beam21,beam31,beam41,
 beam51,beam61,beam71,beam81
*Elset, Elset=beam2
 beam12,beam22,beam32,beam42,
 beam52,beam62,beam72,beam82
*Element, type=spring1, elset=spring4
 104, 101
 114, 111
 124, 121
 134, 131
 144, 141
 154, 151
 164, 161
 174, 171
*Spring, elset=spring4
 4,
 1., 
*Element, type=spring1, elset=spring5
 105, 101
 115, 111
 125, 121
 135, 131
 145, 141
 155, 151
 165, 161
 175, 171
*Spring, elset=spring5
 5,
 1., 
*Element, type=spring1, elset=spring6
 106, 101
 116, 111
 126, 121
 136, 131
 146, 141
 156, 151
 166, 161
 176, 171
*Spring, elset=spring6
 6,
 1., 
*Orientation,name=yzx,system=rectangular
 0.,1.,0., 0.,0.,1.
*Coupling,constraint name=beam2block1,surface=top1,ref node=101
*Distributing			      		   
 1,6				      		   
*Coupling,constraint name=beam2block2,surface=top2,ref node=111
*Distributing			      		   
 1,5				      		   
*Coupling,constraint name=beam2block3,surface=top3,ref node=121
*Distributing			      		   
 1,4				      		   
*Coupling,constraint name=beam2block4,surface=top4,ref node=131
*Distributing			      		   
 1,3				      		   
*Coupling,constraint name=beam2block5,surface=top5,ref node=141,orientation=yzx
*Distributing			      		   
 1,6				      		   
*Coupling,constraint name=beam2block6,surface=top6,ref node=151,orientation=yzx
*Distributing			      		   
 1,4
 6,6				      		   
*Coupling,constraint name=beam2block7,surface=top7,ref node=161,orientation=yzx	   
*Distributing			      		   
 1,3
 6,6				      		   
*Coupling,constraint name=beam2block8,surface=top8,ref node=171,orientation=yzx
*Distributing
 1,3
*Beam Section, Material=Mat2,Section=Circ,Elset=beam1
 1.0, 
 1.,0.,0., 0.,1.,0.
*Beam Section, Material=Mat2,Section=Circ,Elset=beam2
 1.0, 
  0.,0.,1., 0.,1.,0.
*Surface, type=ELEMENT, name=top1
 block1, S2
*Surface, type=ELEMENT, name=top2
 block2, S2
*Surface, type=ELEMENT, name=top3
 block3, S2
*Surface, type=ELEMENT, name=top4
 block4, S2
*Surface, type=ELEMENT, name=top5
 block5, S2
*Surface, type=ELEMENT, name=top6
 block6, S2
*Surface, type=ELEMENT, name=top7
 block7, S2
*Surface, type=ELEMENT, name=top8
 block8, S2
*Surface, type=ELEMENT, name=bot
 blocks, S1
*Surface, type=cylinder, name=plate
 0., 0., -1., 1., 0., -1.
 0., 1., -1.
 start, -1., 0.
 line,   5., 0.
*Rigid body, analytical surface=plate, ref node=1000
*Tie, name=blocksToPlate
 bot, plate
*Nset, nset=refnodes, generate
 101,171,10
*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
*Output, field, freq=1
*Node Output
 U,
*El Print, freq=0
*Node Print, freq=0
*Boundary
 1000,1,6
**
** LOAD CASES
**
*Load Case,name=shear_x
*Boundary
 refnodes,1,1,0.3
 refnodes,4,6
*End Load Case
*Load Case,name=shear_y
*Boundary
 refnodes,2,2,0.3
 refnodes,4,6
*End Load Case
*Load Case,name=pull_z
*Boundary
 refnodes,3,3,0.3
 refnodes,4,6
*End Load Case
*Load Case,name=twist_x
*Boundary
 refnodes,4,4,0.3
*End Load Case
*Load Case,name=twist_y
*Boundary
 refnodes,5,5,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,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, 6
 refnodes,1,1,0.3
 refnodes,4,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-pull
*Static
 .5, 1., 1e-05, 1.
*Boundary
 refnodes,2,2,0.
 refnodes,3,3,0.3
*End Step
** ----------------------------------------------------------------
** 
** STEP: Step-5
** 
*Step, name=Step-5, nlgeom
 x-rotate
*Static
 .5, 1., 1e-05, 1.
*Boundary
 refnodes,3,3,0.
*Boundary
 refnodes,4,4,0.3
 refnodes,5,6,0.
*End Step
** ----------------------------------------------------------------
** 
** STEP: Step-6
** 
*Step, name=Step-6, nlgeom
 x-rotate back
*Static
 .5, 1., 1e-05, 1.
*Boundary
 refnodes,1,6,0.
*End Step
** ----------------------------------------------------------------
** 
** STEP: Step-7
** 
*Step, name=Step-7, nlgeom
 y-rotate
*Static
 .5, 1., 1e-05, 1.
*Boundary
 refnodes,4,4,0.
 refnodes,5,5,0.3
 refnodes,6,6,0.
*End Step
** ----------------------------------------------------------------
** 
** STEP: Step-8
** 
*Step, name=Step-8, nlgeom
 y-rotate back
*Static
 .5, 1., 1e-05, 1.
*Boundary
 refnodes,1,6,0.
*End Step
** ----------------------------------------------------------------
** 
** STEP: Step-9
** 
*Step, name=Step-9, nlgeom
 z-rotate
*Static
 .5, 1., 1e-05, 1.
*Boundary
 refnodes,4,5,0.
 refnodes,6,6,0.3
*End Step
** ----------------------------------------------------------------
** 
** STEP: Step-10
** 
*Step, name=Step-10, nlgeom
 z-rotate back
*Static
 .5, 1., 1e-05, 1.
*Boundary
 refnodes,1,6,0.
*End Step
** ----------------------------------------------------------------
** 
** STEP: Step-11
** 
*Step, name=Step-11, nlgeom
 all
*Static
 .125, 1., 1e-05, .125
*Boundary
 refnodes,1,6,0.3
*End Step