*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 *DENSITY 0.2836 *ELASTIC 30.0E6, 0.3 *BOUNDARY NALL,1,6 *RESTART, WRITE *STEP Normal distributed drag force FDD. *DYNAMIC, EXPLICIT ,.2 *DLOAD, OP=NEW EALL, FDD, 1.0, 2.0, 1.2, 1.0 *OUTPUT, FIELD *NODE OUTPUT RF *OUTPUT, HISTORY *NODE OUTPUT, NSET=NALL RF *END STEP