*heading
 Test of a small spotweld inside a 4-node shell element and compare with VUGENS results
     
          -->|   |<-- 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
 step_time=1.e-3
 square_side = 20.0
 corner      = square_side/2.0
 xsw         = 0.0
 ysw         = 0.0
 field_time = step_time/20
 hist_time  = step_time/200
** Skew the middle node in the x- and y-direction
** delx        = 1.0
** dely        = 2.0
** delx = 3.33333333333
** dely = 2.
** delx = 3.333333333333
** dely = 3.333333333333
 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
 xshift = square_side + corner
 xxsw = xsw + xshift     
 xx1  = x1 + xshift
 xx3  = x3 + xshift
 xx23 = x23 + xshift
 xx21 = x21 + xshift
 xx2  = x2 + xshift
 xx11 = x11 + xshift
 xx12 = x12 + xshift
 xx13 = x13 + xshift
 xx22 = x22 + xshift
*node,nset=nall-top-user
  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-user
 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>
*node,nset=nall-top
 201,  <xx1>,  <y1>, 0.0
 202,  <xx2>,  <y2>, 0.0
 203,  <xx3>,  <y3>, 0.0
 211, <xx11>, <y11>, 0.0
 212, <xx12>, <y12>, 0.0
 213, <xx13>, <y13>, 0.0
 221, <xx21>, <y21>, 0.0
 222, <xx22>, <y22>, 0.0
 223, <xx23>, <y23>, 0.0
 1002,<xxsw>, <ysw>, 0.0
*node,nset=nall-bot
 301,  <xx1>,  <y1>, <z>
 302,  <xx2>,  <y2>, <z>
 303,  <xx3>,  <y3>, <z>
 311, <xx11>, <y11>, <z>
 312, <xx12>, <y12>, <z>
 313, <xx13>, <y13>, <z>
 321, <xx21>, <y21>, <z>
 322, <xx22>, <y22>, <z>
 323, <xx23>, <y23>, <z>
*element, type=s4r, elset=shell-user
 1, 1, 2, 12, 11
 2, 2, 3, 13, 12
 3, 11, 12, 22, 21
 4, 12, 13, 23, 22
*element, type=s4r, elset=shell
 201, 201, 202, 212, 211
 202, 202, 203, 213, 212
 203, 211, 212, 222, 221
 204, 212, 213, 223, 222
*shell general section, elset=shell-user, density = 0.7820E-08, POISSON=0.0, USER, PROPERTIES=21,VARIABLES=7,DELETE=7
3.
3.E+06, 0., 3.E+06,0.,0.,1.5E+06,0.,0.,
0.,2250000,0.,0.,0.,0.,2250000,0.,
0.,0.,0.,0.,1125000.
*TRANSVERSE SHEAR STIFFNESS
1249653., 1249653.,  0.0000 
*shell general section, elset=shell, density = 0.7820E-08, POISSON=0.0, PROPERTIES=21,DELETE=7
3.E+06, 0., 3.E+06,0.,0.,1.5E+06,0.,0.,
0.,2250000,0.,0.,0.,0.,2250000,0.,
0.,0.,0.,0.,1125000.
*TRANSVERSE SHEAR STIFFNESS
1249653., 1249653.,  0.0000 
*element, type=s4r, elset=base-user
 101, 101, 102, 112, 111
 102, 102, 103, 113, 112
 103, 111, 112, 122, 121
 104, 112, 113, 123, 122
*element, type=s4r, elset=base
 301, 301, 302, 312, 311
 302, 302, 303, 313, 312
 303, 311, 312, 322, 321
 304, 312, 313, 323, 322
*shell general section, elset=base-user, density = 0.7820E-08, POISSON=0.0, USER, PROPERTIES=21,VARIABLES=7,DELETE=7
3.
3.E+06, 0., 3.E+06,0.,0.,1.5E+06,0.,0.,
0.,2250000,0.,0.,0.,0.,2250000,0.,
0.,0.,0.,0.,1125000.
*TRANSVERSE SHEAR STIFFNESS
1249653., 1249653.,  0.0000 
*shell general section, elset=base, density = 0.7820E-08, POISSON=0.0, PROPERTIES=21,DELETE=7
3.E+06, 0., 3.E+06,0.,0.,1.5E+06,0.,0.,
0.,2250000,0.,0.,0.,0.,2250000,0.,
0.,0.,0.,0.,1125000.
*TRANSVERSE SHEAR STIFFNESS
1249653., 1249653.,  0.0000 
*surfaceinteraction, name=INTPROP-1
*************************************************
*nset, nset=n1001
 1001,
*nset, nset=n1002
 1002,
*surface, name=top-user, type=element
 shell-user, sneg
*surface, name=bot-user, type=element
 base-user, spos
*surface, name=top, type=element
 shell, sneg
*surface, name=bot, type=element
 base, spos
*fastener, interaction name=spot, reference node set=n1001,property=spotprop,coupling=structural
 top-user, bot-user
*fastener, interaction name=spot2, reference node set=n1002,property=spotprop,coupling=structural
 top, bot
*fastener property, name=spotprop
 0.01,
 1,6
*nset,nset=left_edge
 1,11,21
*nset,nset=center
 2,12,22
*boundary
 nall-bot, 1,6
 nall-bot-user, 1,6
*nset, nset=qa_test
nall-top, nall-top-user 
*************************************************
*step,nlgeom=yes
*DYNAMIC, EXPLICIT, Element by element
 ,<step_time>
*amplitude,name=smooth,definition=smooth
0,0,<step_time>,1
*contact, op=NEW
*contactinclusions, allexterior
*contactpropertyassignment
, , INTPROP-1
*cload,amplitude=smooth
  3, 3, 1.25e5
 13, 3, 2.5e5
 23, 3, 1.25e5
 203, 3, 1.25e5
 213, 3, 2.5e5
 223, 3, 1.25e5
*output,field, ,TIME INTERVAl=<field_time>
*node output
 u, rf
*contact output
Cforce,
*OUTPUT, HISTORY,,TIME INTERVAl=<hist_time>
*node output, nset=qa_test
u,rf
*OUTPUT, HISTORY,TIME INTERVAl=<hist_time>
*ENERGY OUTPUT, VARIABLE=PRESELECT
*end step