*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, SHIFT
 0.0, 0.0, 1.6666
 
*NCOPY, CHANGE NUMBER=200, OLD SET=A, SHIFT
 0.0, 0.0, 3.0

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

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

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

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

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

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

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

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

*ELEMENT, ELSET=ELEM, TYPE=S4
 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
*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=MID
 506,
**
*COUPLING, CONSTRAINT NAME=LOAD, REF NODE=9999, SURFACE=TOP
*DISTRIBUTING
**
*STEP
*STATIC
*BOUNDARY
 EDGE_1, 1, 6 
 EDGE_2, 1, 6
 EDGE_3, 1, 6
 EDGE_4, 1, 6
 9999, 1,1
 9999, 3,6
*CLOAD
 9999, 2, -100
*NODE FILE, FREQ=1000
 U,
 RF,
*NODE PRINT, FREQ=1000, TOTALS=YES, NSET=MID
 U,
*OUTPUT, FIELD, FREQ=1000
*NODE OUTPUT
 U,
 RF,
*END STEP