*heading
 Test of a small spotweld inside a 4-node shell element
 Rotary inertia added at nodes to test the anisotropic rotations
     
          -->|   |<-- delx
    21 o-----o---+---------o
       |     |22           |23
    |  |     |             |
    -  o-----o-------------o
       |11   |12           |13 
    -  +     |   X         |   
    ^  |     |             |       
    |  |     |             |       X - spotweld at (0,0,0), always inside
  dely |     |             |           element 2
       |     |             |
       o-----o---+---------o
       1     2             3
*preprint, echo=yes, model=yes
*parameter
 rota = 1.391e-10
 rotaoff = rota*0.5
 step_time=1.e-3
 time_interval = step_time*0.01
 square_side = 20.0
 corner      = square_side/2.0
** Skew the middle node in the x- and y-direction
 offcenter   = 0.0
 delx = 0.0
 dely = 0.0
 x1  = -corner
 y1  = -corner
 x3  =  corner
 y3  = -corner
 x23 =  corner
 y23 =  corner
 x21 = -corner
 y21 =  corner
 x2  = -delx
 y2  = -corner
 x11 = -corner
 y11 =  dely
 x12 = -delx
 y12 =  dely
 x13 =  corner
 y13 =  dely
 x22 = -delx
 y22 =  corner
 z = -3.0
**xsw  = 0.5*(x12+x13)-offcenter
**ysw  = 0.5*(y2+y12)+offcenter
xsw  = x12
ysw = y12
*node,nset=nall-top
  1,  <x1>,  <y1>, 0.0
  2,  <x2>,  <y2>, 0.0
  3,  <x3>,  <y3>, 0.0
 11, <x11>, <y11>, 0.0
 12, <x12>, <y12>, 0.0
 13, <x13>, <y13>, 0.0
 21, <x21>, <y21>, 0.0
 22, <x22>, <y22>, 0.0
 23, <x23>, <y23>, 0.0
 1001, <xsw>, <ysw>, 0.0
*node,nset=nall-bot
 101,  <x1>,  <y1>, <z>
 102,  <x2>,  <y2>, <z>
 103,  <x3>,  <y3>, <z>
 111, <x11>, <y11>, <z>
 112, <x12>, <y12>, <z>
 113, <x13>, <y13>, <z>
 121, <x21>, <y21>, <z>
 122, <x22>, <y22>, <z>
 123, <x23>, <y23>, <z>
 1002, <xsw>, <ysw>, <z>
*element,type=rotaryI, elset=top
10001,1
10002,2
10003,3
10004,11
10005,12
10006,13
10007,21
10008,22
**10009,23
*rotary inertia, elset=top
<rota>,<rota>,<rota>,<rotaoff>,<rotaoff>,<rotaoff>
*element, type=s4r, elset=shell
 1, 1, 2, 12, 11
 2, 2, 3, 13, 12
 3, 11, 12, 22, 21
 4, 12, 13, 23, 22
*shell general section, elset=shell, material=mat
 3.0,
*element, type=s4r, elset=base
 101, 101, 102, 112, 111
 102, 102, 103, 113, 112
 103, 111, 112, 122, 121
 104, 112, 113, 123, 122
*shell general section, elset=base, material=mat
 3.0,
*material, name=mat
*elastic
 1.e6, 0.0
*density
  0.7820E-08,
*surfaceinteraction, name=INTPROP-1
*************************************************
*nset, nset=n1001
 1001,
*surface, name=top, type=element
 shell, sneg
*surface, name=bot, type=element
 base, spos
*element, type=conn3d2, elset=Fastener_elements-2
155268,      1001,      1002
*CONNECTOR SECTION, ELSET=Fastener_elements-2
BEAM,
*COUPLING, CONSTRAINT NAME=up  , REF NODE=1001, SURFACE=top
*DISTRIBUTING,coupling=structural
*COUPLING, CONSTRAINT NAME=down, REF NODE=1002, SURFACE=bot
*DISTRIBUTING,coupling=structural
*nset,nset=left_edge
 1,11,21
*nset,nset=center
 2,12,22
*boundary
101, 1,6
111, 1,6
121, 1,6
102, 1,6
112, 1,6
122, 1,6
103, 1,6
113, 1,6
123, 1,6
*************************************************
*step,nlgeom=yes
*DYNAMIC, EXPLICIT, DIRECT
**, Element by element
2.e-7 ,<step_time>
*amplitude,name=smooth,definition=smooth
0,0,<step_time>,1
*contact, op=NEW
*contactinclusions, allelementbased
*contactpropertyassignment
, , INTPROP-1
*cload,amplitude=smooth
  3, 3, 1.25e5
 13, 3, 2.5e5
 23, 3, 1.25e5
*OUTPUT,FIELD,NUMBER INTERVAL=20
*node output
 u, rf
*contact output
Cforce,
*nset,nset=QA_TEST
3,13,23
*OUTPUT, HISTORY,time interval = <time_interval>
*node output,nset=QA_TEST
 u,
*ENERGY OUTPUT
ALLKE,ALLAE,ALLIE,ALLWK,ETOTAL
*end step