*HEADING
 Basic verification for COUPLING option
 Test a uniform distributed load applied through the
 COUPLING option to a surface with underlying elements CAX8R
*NODE
  1,  5.0, 0.0
  2,  6.0, 0.0
  3,  8.0, 0.0
  4, 11.0, 0.0
  5, 15.0, 0.0
  6, 20.0, 0.0
  7, 26.0, 0.0
  8,  5.0, 0.5
  9,  6.0, 0.5
 10,  8.0, 0.5
 11, 11.0, 0.5
 12, 15.0, 0.5
 13, 20.0, 0.5
 14, 26.0, 0.5
 15,  5.0, 1.0
 16,  6.0, 1.0
 17,  8.0, 1.0
 18, 11.0, 1.0
 19, 15.0, 1.0
 20, 20.0, 1.0
 21, 26.0, 1.0
 22,  5.0, 2.0
 23,  6.0, 2.0
 24,  8.0, 2.0
 25, 11.0, 2.0
 26, 15.0, 2.0
 27, 20.0, 2.0
 28, 26.0, 2.0
 29,  5.0, 3.0
 30,  6.0, 3.0
 31,  8.0, 3.0
 32, 11.0, 3.0
 33, 15.0, 3.0
 34, 20.0, 3.0
 35, 26.0, 3.0
 36,  5.0, 4.5
 37,  6.0, 4.5
 38,  8.0, 4.5
 39, 11.0, 4.5
 40, 15.0, 4.5
 41, 20.0, 4.5
 42, 26.0, 4.5
 43,  5.0, 6.0
 44,  6.0, 6.0
 45,  8.0, 6.0
 46, 11.0, 6.0
 47, 15.0, 6.0
 48, 20.0, 6.0
 49, 26.0, 6.0
 50,  5.0, 8.0
 51,  6.0, 8.0
 52,  8.0, 8.0
 53, 11.0, 8.0
 54, 15.0, 8.0
 55, 20.0, 8.0
 56, 26.0, 8.0
 57,  5.0, 10.0
 58,  6.0, 10.0
 59,  8.0, 10.0
 60, 11.0, 10.0
 61, 15.0, 10.0
 62, 20.0, 10.0
 63, 26.0, 10.0
*ELEMENT, TYPE=CAX4, ELSET=EALL
 1, 1,2,9,8
*ELGEN, ELSET=EALL
 1, 6, 1, 1, 8, 7, 6
*SOLID SECTION, MATERIAL=ELASTIC, ELSET=EALL
*MATERIAL, NAME=ELASTIC
*ELASTIC
 30.0e6,0.3
*ELSET, ELSET=FACE
  1,2,3,4,5,6
*SURFACE, NAME=LSURFACE, TRIM=YES
  FACE,
*NODE, NSET=NREF
 9999, 17.870968, -5.0
*COUPLING, CONSTRAINT NAME=LOAD_DISTRIBUTOR, SURFACE=LSURFACE, REF NODE=NREF
*DISTRIBUTING, WEIGHT=UNIFORM
*BOUNDARY
 57, 2, 2
 58, 2, 2
 59, 2, 2
 60, 2, 2
 61, 2, 2
 62, 2, 2
 63, 2, 2
*STEP, NLGEOM
 Apply uniform load
*STATIC
 ,1.0
*CLOAD
 NREF, 2, 2.e6
*EL PRINT, FREQ=10000
 S
*NODE PRINT, FREQ=10000
 U, RF
*NODE FILE, FREQ=10000
 U, RF
*OUTPUT, FIELD, FREQ=1000
*NODE OUTPUT
 U, RF
*END STEP