*heading
*node
   1, -1., -1., -1.
   2,  1., -1., -1.
   3,  1.,  1., -1.
   4, -1.,  1., -1.
   5, -1., -1.,  1.
   6,  1., -1.,  1.
   7,  1.,  1.,  1.
   8, -1.,  1.,  1.
  11, -1., -1., -1.2
  12,  0., -1., -1.2
  13,  1., -1., -1.2
  14, -1.,  0., -1.2
  15,  0.,  0., -1.2
  16,  1.,  0., -1.2
  17, -1.,  1., -1.2
  18,  0.,  1., -1.2
  19,  1.,  1., -1.2
  21, -1., -1.,  1.2
  22,  0., -1.,  1.2
  23,  1., -1.,  1.2
  24, -1.,  0.,  1.2
  25,  0.,  0.,  1.2
  26,  1.,  0.,  1.2
  27, -1.,  1.,  1.2
  28,  0.,  1.,  1.2
  29,  1.,  1.,  1.2
  31, -1.,-1.2, -1.
  32,  0.,-1.2, -1.
  33,  1.,-1.2, -1.
  34, -1.,-1.2,  0.
  35,  0.,-1.2,  0.
  36,  1.,-1.2,  0.
  37, -1.,-1.2,  1.
  38,  0.,-1.2,  1.
  39,  1.,-1.2,  1.
  41, 1.2, -1., -1.
  42, 1.2,  0., -1.
  43, 1.2,  1., -1.
  44, 1.2, -1.,  0.
  45, 1.2,  0.,  0.
  46, 1.2,  1.,  0.
  47, 1.2, -1.,  1.
  48, 1.2,  0.,  1.
  49, 1.2,  1.,  1.
  51, -1., 1.2, -1.
  52,  0., 1.2, -1.
  53,  1., 1.2, -1.
  54, -1., 1.2,  0.
  55,  0., 1.2,  0.
  56,  1., 1.2,  0.
  57, -1., 1.2,  1.
  58,  0., 1.2,  1.
  59,  1., 1.2,  1.
  61,-1.2, -1., -1.
  62,-1.2,  0., -1.
  63,-1.2,  1., -1.
  64,-1.2, -1.,  0.
  65,-1.2,  0.,  0.
  66,-1.2,  1.,  0.
  67,-1.2, -1.,  1.
  68,-1.2,  0.,  1.
  69,-1.2,  1.,  1.
 101,  0.,  0., -1.2
 102,  0.,  0.,  1.2
 103,  0.,-1.2,  0.
 104, 1.2,  0.,  0.
 105,  0., 1.2,  0.
 106,-1.2,  0.,  0.
 201,  0.,  0., -1.2
 202,  0.,  0.,  1.2
 203,  0.,-1.2,  0.
 204, 1.2,  0.,  0.
 205,  0., 1.2,  0.
 206,-1.2,  0.,  0.
 301,  0.,  0., -1.2
 302,  0.,  0.,  1.2
 303,  0.,-1.2,  0.
 304, 1.2,  0.,  0.
 305,  0., 1.2,  0.
 306,-1.2,  0.,  0.
 401,  0.,  0., -1.2
 402,  0.,  0.,  1.2
 403,  0.,-1.2,  0.
 404, 1.2,  0.,  0.
 405,  0., 1.2,  0.
 406,-1.2,  0.,  0.
*nset, nset=cube_corners, generate
 1, 8, 1
*nset, nset=shell_1_corners
 11,13,17,19
*nset, nset=shell_2_corners
 21,23,27,29
*nset, nset=shell_3_corners
 31,33,37,39
*nset, nset=shell_4_corners
 41,43,47,49
*nset, nset=shell_5_corners
 51,53,57,59
*nset, nset=shell_6_corners
 61,63,67,69
*nset, nset=ref_1
 101,
*nset, nset=ref_2
 102,
*nset, nset=ref_3
 103,
*nset, nset=ref_4
 104,
*nset, nset=ref_5
 105,
*nset, nset=ref_6
 106,
*mpc
 beam, shell_1_corners, ref_1
 beam, shell_2_corners, ref_2
 beam, shell_3_corners, ref_3
 beam, shell_4_corners, ref_4
 beam, shell_5_corners, ref_5
 beam, shell_6_corners, ref_6
*element, type=c3d8, elset=cube
 1, 1,2,3,4,5,6,7,8
*solid section, elset=cube, material=mat1
*element, type=s4r, elset=shell_1
 11, 11,14,15,12
 12, 12,15,16,13
 13, 14,17,18,15
 14, 15,18,19,16
*element, type=s4r, elset=shell_2
 21, 21,22,25,24
 22, 22,23,26,25
 23, 24,25,28,27
 24, 25,26,29,28
*element, type=s4r, elset=shell_3
 31, 31,32,35,34
 32, 32,33,36,35
 33, 34,35,38,37
 34, 35,36,39,38
*element, type=s4r, elset=shell_4
 41, 41,42,45,44
 42, 42,43,46,45
 43, 44,45,48,47
 44, 45,46,49,48
*element, type=s4r, elset=shell_5
 51, 51,54,55,52
 52, 52,55,56,53
 53, 54,57,58,55
 54, 55,58,59,56
*element, type=s4r, elset=shell_6
 61, 61,64,65,62
 62, 62,65,66,63
 63, 64,67,68,65
 64, 65,68,69,66
*elset, elset=shell
 shell_1,shell_2,shell_3,shell_4,shell_5,shell_6
*shell section, elset=shell, material=mat1
 0.1, 5
*element, type=conn3d2, elset=connect_1
 101, 201, 301
*element, type=conn3d2, elset=connect_2
 102, 202, 302
*element, type=conn3d2, elset=connect_3
 103, 203, 303
*element, type=conn3d2, elset=connect_4
 104, 204, 304
*element, type=conn3d2, elset=connect_5
 105, 205, 305
*element, type=conn3d2, elset=connect_6
 106, 206, 306
*orientation, name=oriFast_1
 0.,1.,0., 1.,0.,0.
*orientation, name=oriFast_2
 0.,1.,0., -1.,0.,0.
*orientation, name=oriFast_3
 0.,0.,1., -1.,0.,0.
*orientation, name=oriFast_4
 0.,1.,0., 0.,0.,1.
*orientation, name=oriFast_5
 0.,0.,1., 1.,0.,0.
*orientation, name=oriFast_6
 0.,1.,0., 0.,0.,-1.
*orientation, name=oriConn_1
 1.,0.,0., 0.,1.,0.
*orientation, name=oriConn_2
 1.,0.,0., 0.,1.,0.
*orientation, name=oriConn_3
 1.,0.,0., 0.,1.,0.
*orientation, name=oriConn_4
 1.,0.,0., 0.,1.,0.
*orientation, name=oriConn_5
 1.,0.,0., 0.,1.,0.
*orientation, name=oriConn_6
 1.,0.,0., 0.,1.,0.
*connector section, elset=connect_1
 translator,
 oriConn_1,
*connector section, elset=connect_2
 translator,
 oriConn_2,
*connector section, elset=connect_3
 translator,
 oriConn_3,
*connector section, elset=connect_4
 translator,
 oriConn_4,
*connector section, elset=connect_5
 translator,
 oriConn_5,
*connector section, elset=connect_6
 translator,
 oriConn_6,
*surface, name=surfShell_1
 shell_1, spos
*surface, name=surfShell_2
 shell_2, spos
*surface, name=surfShell_3
 shell_3, spos
*surface, name=surfShell_4
 shell_4, spos
*surface, name=surfShell_5
 shell_5, spos
*surface, name=surfShell_6
 shell_6, spos
*surface, name=surfCube_1
 cube, s1
*surface, name=surfCube_2
 cube, s2
*surface, name=surfCube_3
 cube, s3
*surface, name=surfCube_4
 cube, s4
*surface, name=surfCube_5
 cube, s5
*surface, name=surfCube_6
 cube, s6
*fastener, interaction=int-1, elset=connect_1, property=prop, orient=oriFast_1
 surfShell_1, surfCube_1
*fastener, interaction=int-1, elset=connect_2, property=prop, orient=oriFast_2
 surfShell_2, surfCube_2
*fastener, interaction=int-1, elset=connect_3, property=prop, orient=oriFast_3
 surfShell_3, surfCube_3
*fastener, interaction=int-1, elset=connect_4, property=prop, orient=oriFast_4
 surfShell_4, surfCube_4
*fastener, interaction=int-1, elset=connect_5, property=prop, orient=oriFast_5
 surfShell_5, surfCube_5
*fastener, interaction=int-1, elset=connect_6, property=prop, orient=oriFast_6
 surfShell_6, surfCube_6
*fastener property, name=prop
 1.,
*material, name=mat1
*elastic
 30.e6, 0.3
**
** --------------------------------------------------
**
*step, nlgeom=yes, name=step-1
 static 
*static
 0.05, 1.
*boundary
 cube_corners, 1,3
 ref_1,2,2,1.
 ref_2,2,2,1.
 ref_3,3,3,1.
 ref_4,2,2,1.
 ref_5,3,3,1.
 ref_6,2,2,1.
*node print, freq=0
*output,field
*node output
 u,
*element output, directions=yes
 cu
*end step