*HEADING
 Basic verification for *COUPLING
 Test tributary area of surface with underlying element s4
**
** Uniformly Loaded plate clamped at all 4 edges
** Solution according to Theory of Plates and Shells, 
**                       Timoshenko and Woinowsky-Krieger, 
**
**
**        0.00126 (a/b) (t)^4
** wmax = -------------------  = 0.00459
**                 D
**
** where 
**
**               Eh^3
**        D = ---------  = 2747.25;  (a/b) = 1.0;  t=0.1
**            12(1-v^2)
**
**
*NODE, NSET=A
  1,  0.0000,  0.0,  0.0
  2,  1.6666,  0.0,  0.0
  3,  3.0000,  0.0,  0.0
  4,  4.0000,  0.0,  0.0
  5,  4.6666,  0.0,  0.0
  6,  5.0000,  0.0,  0.0
  7,  5.3333,  0.0,  0.0
  8,  6.0000,  0.0,  0.0
  9,  7.0000,  0.0,  0.0
 10,  8.3333,  0.0,  0.0
 11, 10.0000,  0.0,  0.0
*NCOPY, CHANGE NUMBER=100, OLD SET=A, NEW SET=A1,SHIFT
 0.0, 0.0, 1.6666
 
*NCOPY, CHANGE NUMBER=200, OLD SET=A, NEW SET=A2,SHIFT
 0.0, 0.0, 3.0

*NCOPY, CHANGE NUMBER=300, OLD SET=A, NEW SET=A3,SHIFT
 0.0, 0.0, 4.0

*NCOPY, CHANGE NUMBER=400, OLD SET=A, NEW SET=A4,SHIFT
 0.0, 0.0, 4.6666

*NCOPY, CHANGE NUMBER=500, OLD SET=A, NEW SET=A5,SHIFT
 0.0, 0.0, 5.0

*NCOPY, CHANGE NUMBER=600, OLD SET=A, NEW SET=A6,SHIFT
 0.0, 0.0, 5.3333

*NCOPY, CHANGE NUMBER=700, OLD SET=A, NEW SET=A7,SHIFT
 0.0, 0.0, 6.0

*NCOPY, CHANGE NUMBER=800, OLD SET=A, NEW SET=A8,SHIFT
 0.0, 0.0, 7.0

*NCOPY, CHANGE NUMBER=900, OLD SET=A, NEW SET=A9,SHIFT
 0.0, 0.0, 8.66666

*NCOPY, CHANGE NUMBER=1000,OLD SET=A, NEW SET=A10,SHIFT
 0.0, 0.0, 10.0

*ELEMENT, ELSET=ELEM, TYPE=S4R
 1, 1, 2, 102, 101
*ELGEN, ELSET=ELEM
 1, 10, 1, 1, 10, 100, 11
*SHELL SECTION, ELSET=ELEM, MATERIAL=ELASTIC
 0.1,
*MATERIAL, NAME=ELASTIC
*ELASTIC
 30e6, 0.3
*DENSITY
.1,
*ELEMENT, TYPE=MASS, ELSET=MASS
10000,9999
*MASS, ELSET=MASS
1.E-6,
*ELEMENT, TYPE=ROTARYI, ELSET=ROTARY
10001,9999
*ROTARY INERTIA, ELSET=ROTARY
1.E-6,1.E-6,1.E-6
*NSET, NSET=EDGE_1, GEN
 1, 11, 1
*NSET, NSET=EDGE_2, GEN
 111, 1011, 100
*NSET, NSET=EDGE_3, GEN
 1001, 1010, 1
*NSET, NSET=EDGE_4, GEN
 101, 901, 100
*NSET, NSET=MID
 1,
*SURFACE, NAME=TOP
 ELEM, SPOS
*NODE
 9999, 5.0, 5.0, 5.0
*NSET, NSET=QA_TEST_MID
 506,
**
*COUPLING, CONSTRAINT NAME=LOAD, REF NODE=9999, SURFACE=TOP
*DISTRIBUTING
*AMPLITUDE,NAME=AMP,DEFINITION=SMOOTH STEP
0.,0.,.1,1,
**
*NSET, NSET=QA_TEST_ALLNODES
A,A1,A2,A3,A4,A5,A6,A7,A8,A9,A10
*STEP
*DYNAMIC, EXPLICIT
,1.E-1
**,1.0
*BOUNDARY
 EDGE_1, 1, 6 
 EDGE_2, 1, 6
 EDGE_3, 1, 6
 EDGE_4, 1, 6
 9999, 1,1
 9999, 3,6
*CLOAD,AMPLITUDE=AMP
 9999, 2, -1.E+2
*OUTPUT, FIELD, NUMBER=1
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
 U, 
**OUTPUT, HISTORY,FREQ=1
**NODE OUTPUT,NSET=QA_TEST_MID
**U,
*END STEP