*HEADING
: AQUA Spatial variation of steady state current velocity.
**
** AQUA steady state currents are tested with respect to spatial
** variation. A pipe element fully constrained and submerged,
** is subjected a steady state current which varies with respect to 
** the x-coordinates and y-coordinates. For this simple case the total 
** drag force can be analytically determined for each member and 
** compared to the member reaction forces.
**
** EXPECTED RESULTS:
**
**      R = 1/4*rho*Cd*D*vel^2*L
**
** where    R = Member reaction at either node
**        rho = fluid density = 1.987
**         Cd = Drag Coefficient = 1.2
**          D = Effective Diameter = 2.0
**        vel = Steady current velocity = fn(x)
**          L = Member Length = 10
**         vf = 2.896
**
** Expected Results:
**
** Nodes      |    Velocity   |   Ri
** -----------------------------------
**   1,11     |      1 vf     |  100
**   2,12     |      2 vf     |  400
**   3,13     |      3 vf     |  900
**   4,14     |      3 vf     |  900
**   5,15     |    1.5 vf     |  225
**   6,16     |      0 vf     |    0
** -----------------------------------
**     21     |    1.5 vf     |  225
**     22     |      3 vf     |  900
**     23     |    4.5 vf     | 2025
**     24     |    4.5 vf     | 2025
**     25     |   2.25 vf     |  506.25
**     26     |      0 vf     |    0
** -----------------------------------
**  31,41     |      2 vf     |  400
**  32,42     |      4 vf     | 1600
**  33,43     |      6 vf     | 3600
**  34,44     |      6 vf     | 3600
**  35,45     |      3 vf     |  900
**  36,46     |      0 vf     |    0
**
**
** Note: Add 100 for the other member nodes
**
*AQUA
0.0,40.0,32.2,1.987
**
**    vx       vy     vz    elev      x1    x2
**-----------------------------------------------------
5.792358548,  0.0,   0.0,    0.0,  200.0, 100.0
5.792358548,  0.0,   0.0,   40.0,  200.0, 100.0
2.896179274,  0.0,   0.0,    0.0,  800.0, 100.0
2.896179274,  0.0,   0.0,   40.0,  800.0, 100.0
17.37707564,  0.0,   0.0,    0.0,  200.0, 300.0
17.37707564,  0.0,   0.0,   40.0,  200.0, 300.0
8.688537822,  0.0,   0.0,    0.0,  800.0, 300.0
8.688537822,  0.0,   0.0,   40.0,  800.0, 300.0
17.37707564,  0.0,   0.0,    0.0,  200.0, 600.0
17.37707564,  0.0,   0.0,   40.0,  200.0, 600.0
8.688537822,  0.0,   0.0,    0.0,  800.0, 600.0
8.688537822,  0.0,   0.0,   40.0,  800.0, 600.0
0.0,          0.0,   0.0,    0.0,  200.0, 900.0
0.0,          0.0,   0.0,   40.0,  200.0, 900.0
0.0,          0.0,   0.0,    0.0,  800.0, 900.0
0.0,          0.0,   0.0,   40.0,  800.0, 900.0
**
**
*NODE, NSET=NALL
1,   900.0,  50.0,  20.0
2,   900.0, 200.0,  20.0
3,   900.0, 300.0,  20.0
4,   900.0, 450.0,  20.0
5,   900.0, 750.0,  20.0
6,   900.0, 990.0,  20.0
** ----------------------
11,  800.0,  50.0,  20.0
12,  800.0, 200.0,  20.0
13,  800.0, 300.0,  20.0
14,  800.0, 450.0,  20.0
15,  800.0, 750.0,  20.0
16,  800.0, 990.0,  20.0
** ----------------------
21,  500.0,  50.0,  20.0
22,  500.0, 200.0,  20.0
23,  500.0, 300.0,  20.0
24,  500.0, 450.0,  20.0
25,  500.0, 750.0,  20.0
26,  500.0, 990.0,  20.0
** ----------------------
31,  200.0,  50.0,  20.0
32,  200.0, 200.0,  20.0
33,  200.0, 300.0,  20.0
34,  200.0, 450.0,  20.0
35,  200.0, 750.0,  20.0
36,  200.0, 990.0,  20.0
** ----------------------
41,  100.0,  50.0,  20.0
42,  100.0, 200.0,  20.0
43,  100.0, 300.0,  20.0
44,  100.0, 450.0,  20.0
45,  100.0, 750.0,  20.0
46,  100.0, 990.0,  20.0
** ----------------------
101, 900.0,  50.0,  30.0
102, 900.0, 200.0,  30.0
103, 900.0, 300.0,  30.0
104, 900.0, 450.0,  30.0
105, 900.0, 750.0,  30.0
106, 900.0, 990.0,  30.0
** ---------------------
111, 800.0,  50.0,  30.0
112, 800.0, 200.0,  30.0
113, 800.0, 300.0,  30.0
114, 800.0, 450.0,  30.0
115, 800.0, 750.0,  30.0
116, 800.0, 990.0,  30.0
** ----------------------
121, 500.0,  50.0,  30.0
122, 500.0, 200.0,  30.0
123, 500.0, 300.0,  30.0
124, 500.0, 450.0,  30.0
125, 500.0, 750.0,  30.0
126, 500.0, 990.0,  30.0
** ----------------------
131, 200.0,  50.0,  30.0
132, 200.0, 200.0,  30.0
133, 200.0, 300.0,  30.0
134, 200.0, 450.0,  30.0
135, 200.0, 750.0,  30.0
136, 200.0, 990.0,  30.0
** ----------------------
141, 100.0,  50.0,  30.0
142, 100.0, 200.0,  30.0
143, 100.0, 300.0,  30.0
144, 100.0, 450.0,  30.0
145, 100.0, 750.0,  30.0
146, 100.0, 990.0,  30.0
** ----------------------
1000, 0, 990.0, 30.0
** ---------------------
*ELEMENT, TYPE=PIPE31, ELSET=PIPES
1, 1, 101, 1000
2, 2, 102, 1000
3, 3, 103, 1000
4, 4, 104, 1000
5, 5, 105, 1000
6, 6, 106, 1000
** ---------------
11, 11, 111, 1000
12, 12, 112, 1000
13, 13, 113, 1000
14, 14, 114, 1000
15, 15, 115, 1000
16, 16, 116, 1000
** ---------------
21, 21, 121, 1000
22, 22, 122, 1000
23, 23, 123, 1000
24, 24, 124, 1000
25, 25, 125, 1000
26, 26, 126, 1000
** ---------------
31, 31, 131, 1000
32, 32, 132, 1000
33, 33, 133, 1000
34, 34, 134, 1000
35, 35, 135, 1000
36, 36, 136, 1000
** ---------------
41, 41, 141, 1000
42, 42, 142, 1000
43, 43, 143, 1000
44, 44, 144, 1000
45, 45, 145, 1000
46, 46, 146, 1000
** ---------------
*BEAM SECTION, SECTION=PIPE, MATERIAL=STEEL, ELSET=PIPES
1.0, 0.05
*ELSET, ELSET=EALL
PIPES, 
*MATERIAL, NAME=STEEL
*ELASTIC
30.0E9, 0.3
*BOUNDARY
NALL,1,6
*RESTART, WRITE
*STEP
Normal distributed drag force FDD.
*STATIC
.1,.2
*DLOAD, OP=NEW
EALL, FDD, 1.0, 2.0, 1.2, 1.0
*NODE FILE, FREQ=100
RF, 
*EL FILE, FREQ=100
LOADS, 
*NODE PRINT
RF, 
*EL PRINT
LOADS,  
*END STEP