DISEO HIDRAULICO DE UNA RAPIDA

Ubicacin: Prog. 0+670 Km

1) Datos del canal aguas arriba

3
m
Caudal : Q 2.25
s
Ancho base del canal: b 1.20 m

Talud : z 0.5

Pendiente : S 0.001

1
Rugosidad : 3
n 0.014s m

2

3
( b z y) y
( b z y ) y S
b 2 1 z2 y
y rootQ
y 0.01 m 1.5m
n

yn y 0.963 m ;Tirante normal

2
A ( b z y) y 1.619 m ;Area

Q m
V 1.39 ;Velocidad
A s

T b 2 z y 2.163 m ;Espejo de agua

V
Fr 0.513 ;Nmero de Froude
A
g
T

if ( Fr 1 "flujo supercrtico" "flujo subcrtico" ) "flujo subcrtico"

2
V
En y 1.061 m
2g

2. Datos del canal aguas abajo

Ancho base del canal : b5 1.80 m

Talud : z5 0

Pendiente : S5 0.001

1
3
Rugosidad : n5 0.014s m

2

3
( b5 z5 y5) y5
( b5 z5 y5) y5 S5
b5 2 1 z52 y5
y5 rootQ
y50.01 m 1.5m
n5

y5 0.932 m ;Tirante normal

 2
A ( b5 z5 y5) y5 1.677 m ;Area

Q m
v5 1.342 ;Velocidad
A s

T b5 2 z5 y5 1.8 m ;Espejo de agua

v5
F 0.444 ;Nmero de Froude
A
g
T

if ( F 1 "flujo supercrtico" "flujo subcrtico" ) "flujo subcrtico"

2
v5
E y5 1.023 m
2g

3. Datos de elevaciones del perfil longitudinal de la Rpida

PROGI 670m ;Progresiva al inicio del canal de la rpida - seccin de control

C 1 2835m ;Elevacin al inicio del canal de la rpida - seccin de control

C 12 0 m ;Elevacin al inicio del 2do tramo del canal de la rpida

C 23 0 m ;Elevacin al inicio del 3er tramo del canal de la rpida

C 2 2822.00m ;Elevacin al final del canal de la rpida - inicio de la trayectoria

C 5 2821.60m ;Elevacin a la salida de la poza - canal aguas abajo

4) Canal de la rpida
0.4
s 0.4
br 0.765 Q 1.058 m ;ancho del canal de la rpida, segn la frmula de Dadenkov
.2
m
Ancho asumido: br 1.00 m

Talud: zr 0

Pendiente canal de la rpida: Sr 0.34

1
Rugosidad : 3
nr 0.014s m

5) Seccin de Control

Tirante crtico (yc)

Q2 [ ( br zr yc) yc]
3
yc root yc 0 m 1 m
g br 2 zr yc

yc 0.802 m ;Tirante crtico

2
Ac ( br zr yc) yc 0.802 m ;Area

Q m
Vc 2.805 ;Velocidad
Ac s

Tirante al inicio del canal de la rpida (yo)

yc if ( Fr 1 ynyc)

yc 0.802 m

Altura de solera (h)

Figura 1. Seccin de control

;Ecuacin de la conservacin de energa en la seccin de control

Vc
2
h if Fr 1 0 En yc
2g

h 0.142 m
6) Tirantes (Y) de escurrimiento en el canal de la rpida

Figura 2. Secciones consideradas, mtodo de los incrementos finitos

;Ecuacin de la conservacin de la energa

en para el tramo: i, i+1

;hf prdida de energa por friccin en el tramo

;Sf pendiente de friccin

C1 C2 ;L, longitud del canal de la rpida

L 38.235 m
Sr

N 20 ;N, nmero de tramos a subdividir la longitud total del canal de la rpida

L
L 1.912 m
N

x L 1 Sr2 2.019 m
M 0 1 N

DX M PROGI M L
Y Y0 yc

for i 0 N

Ai br zr Yi Yi
Q
Vi
Ai

2
Pi br 2 Yi 1 zr

Ai
Ri
Pi
2
nr Vi 2
Sf i
4

Ri 3
2 Q
2
nr
Sf [ ( br zr Yi1) Yi1]
i 4
3
2 ( br zr Yi1) Yi1
Q
Vi 2 br 2 Yi1 1 zr2
Yi 1 rootL Sr Yi Yi1
[ ( br zr Yi1) Yi1] x Yi1 0.001m yc
2g 2g 2
Yi 1 Yi 1

AM br zr YM YM ;Area mojada

Q
VM ;Velocidad
AM

TM br 2 zr YM ;Espejo de agua

VM
FrM
A
;Nmero de Froude
M
g
TM

Aireacin del flujo

Para canales que tienen velocidad de flujo alto, la superficie de agua puede ser ligeramente superior
que lo esperado debido al ingreso de aire, que debe tenerse en cuenta al calcular la altura de los
muros del canal de la rpida.

El HEC-RAS considera sta situacin a travs de las siguientes dos ecuaciones:

;para nmeros de Froude menor o igual a 8,2

;para nmeros de Froude mayor a 8,2

Donde:

Ya = tirante de agua con el ingreso de aire

Y = tirante de agua sin el ingreso de aire
F = Nmero de Froude

0.061 Fr
M
0.1051 Fr
M
Ya M if FrM 8.2 0.906 YM e 0.620 YM e
 DX M YM Ya M VM
670.00 m 0.802 m 0.773 m 2.805 m
671.91 0.433 0.458 5.193 s
673.82 0.359 0.399 6.262
675.74 0.319 0.368 7.058
677.65 0.292 0.350 7.693
679.56 0.274 0.337 8.214
681.47 0.260 0.328 8.648
683.38 0.250 0.321 9.014
685.29 0.241 0.316 9.323
687.21 0.235 0.313 9.586
689.12 0.229 0.310 9.810
691.03 0.225 0.307 10.000
692.94 0.221 0.305 10.163
694.85 0.218 0.304 10.302
696.76 0.216 0.303 10.422
698.68 0.214 0.302 10.523
700.59 0.212 0.301 10.611
702.50 0.211 0.300 10.685
704.41 0.209 0.300 10.749
706.32 0.208 0.299 10.804
708.24 0.207 0.299 10.851

Perfil de la superficie de agua

1
Tirante
Tirantes(m)

0.8
Y Tirante aireado
M
0.6
Ya
M
0.4

0.2
670 680 690 700
DX
M

Distancia-Progresiva(m)

Velocidad del flujo de agua en el canal de la rpida

12
Velocidad(m/s)

10

8
V
M
6

4 Velocidad
2
670 680 690 700
DX
M

Distancia-Progresiva(m)
7) Estabilidad del flujo en el canal de la rpida

El flujo uniforme se vuelve inestable cuando la velocidad de flujo es muy alta o la pendiente del canal es
muy empinada. Cuando esto ocurre, la inestabilidad de la superficie libre se caracteriza por la formacion
de una serie de ondas de remolino.

En el diseo de una rpida debe evitarse la formacion de ondas flujo pulstil en el canal de la rpida. La
amplitud de las olas puede superar la altura de los muros del canal, provocando desbordes y asimismo
generar un resalto hidrulico inestable en la poza de disipacin.

Para verificar la posible existencia de rgimen inestable y flujo pulstil en el canal de la rpida, se aplican
los criterios del Nmero de Vedernikov y de Montuori.

El Nmero de Vedernikov:

La ecuacin bsica que caracteriza la condicin de flujo de acuerdo a la estabilidad de la superficie libre
es la que expresa el Nmero de Vedernikov ( Nv ) (Chow, 1959).

Tabla 1. Clasificacion del rgimen de flujo en canales abiertos

Para un canal de seccin trapezoidal

y Fr es el nmero de Froude

El Nmero de Montouri (M2 )

Para que no se produzcan trenes de onda en el canal de la rpida, el valor mnimo del Nmero de
Montouri (M2 ) debe ser 0.2.

Donde :

g = Aceleracin de la gravedad,
L = Longitud del tramo en consideracin,
s = Pendiente media del gradiente de energa, s=tan()
V = Velocidad y
= Angulo de inclinacin del gradiente de energa.
Verificacin por flujo inestable y pulstil en el canal de la rpida
Durante la operacin del canal el caudal puede variar en un rango de caudales, por lo tanto es necesario
verificar el diseo de la rpida para asegurar su operacin en dicho rango.

Calculados los Numeros de Vedernikov y de Mountori, se ingresa a la Figura 3 y se verifica si alguno de

los valores cae dentro de la zona de flujo pulsante, se deben verificar puntos intermedios para determinar
el punto en el que las ondas comienzan a formarse.

Figura 3. Criterio para flujo estable (sin ondas)

M 1 N

DX M M L
;DX es la longitud del tramo en consideracin
2
Vc
E1 C 1 yc ;Energa al inicio del canal de la rpida
2g

VM 2 ;Energa a lo largo del canal de la rpida

EM C1 Sr DXM YM
2g

E1 EM
M atan ;, ngulo de inclinacin del gradiente de energa
DXM

VM 2 ;M2, Nmero de Mountori

M2 M
g E1 EM cos M

Y YM
2
1 2 zr M 2 zr 1 zr
2
2 br
Nv M
br
;Nv, Nmero de Vedernikov
3 FrM
YM 2 YM
1 2 1 zr 1 2 zr
br br

Verificacin M if Nv M 9.187 M2M 1.12 "Flujo estable" "Flujo inestable y pulstil"
ProgresivaM DX M PROGI

ProgresivaM M2 M Nv M 0
671.912 m 61.273 0.900 0 "Flujo estable"
673.824 27.624 1.294 1 "Flujo estable"
675.735 17.257 1.625 2 "Flujo estable"
677.647 12.252 1.911 3 "Flujo estable"
679.559 9.332 2.159 4 "Flujo estable"
681.471 7.438 2.374 5 "Flujo estable"
683.382 6.120 2.562 6 "Flujo estable"
685.294 5.156 2.725 7 "Flujo estable"
687.206 4.426 2.867 8 "Flujo estable"
689.118 3.856 2.989 Verificacin M 9 "Flujo estable"
691.029 3.401 3.095 10 "Flujo estable"
692.941 3.031 3.187 11 "Flujo estable"
694.853 2.725 3.266 12 "Flujo estable"
696.765 2.469 3.335 13 "Flujo estable"
698.676 2.252 3.394 14 "Flujo estable"
700.588 2.066 3.444 15 "Flujo estable"
702.500 1.905 3.488 16 "Flujo estable"
704.412 1.765 3.526 17 "Flujo estable"
706.324 1.643 3.558 18 "Flujo estable"
708.235 1.535 3.586 19 "Flujo estable"

8. Trayectoria de la Rpida

Curva parablica

Figura 4. Trayectoria parablica.

= ngulo de inclinacion del canal de la rpida en el origen de la trayectoria (tg=So )

Vo = Velocidad al incio de la trayectoria,
yo = Tirante al incio de la trayectoria,
Vmax= 1.5Vo , para asegurar una presin positiva del flujo sobre el piso y
Z = Talud de la rampa de ingreso a la poza
So Sr 0.34
m
Vo VN 10.851
s
m
Vmax 1.5 Vo 16.276
s

yo YN 0.207 m

Z 2

2
g x
f ( x) C2 x So
2
1 So ;Ecuacin de la trayectoria parablica de la rpida
2
2 ( Vmax)

1 So
Z
X 2.25 ( Vo) 2 20.34 m f ( X) 2806.542 m ;Coordenada del P.T. =(X,f(X))

g 1 So2

j 0 10

X
xj j
10

Coordenadas de la trayectoria:
xj PROGI L f xj Perfil de la trayectoria de la rpida
708.235 m 2822.00 m 2825

710.269 2821.22
712.303 2820.28 2820
Elevacin(m)

714.337 2819.16
716.371 2817.87
f x j 2815

718.405 2816.41
2810
720.439 2814.78
722.473 2812.97
2805
724.507 2811.00 700 710 720 730

726.541 2808.86 x PROGI L

j
728.575 2806.54 Distancia-x(m)
Curva circular
Tambin se puede adoptar una curva vertical circular para unir la pendiente ltima de la rpida con el
plano inclinado del principio del cuenco amortiguador.

El radio de la curva circular de transicin est definido por la ecuacin siguiente (Meshgin and Moore,
1970).
Donde:

r = radio de la curva de transicin, m

yo = tirante al inicio de la curva , m
Fr = Nmero de Froude al inicio de la curva
yo
r r 7.901 m
1.5

FrN
2
Figura 5. Trayectoria circular
e 1

9. Clculo de los tirantes conjugados y3, y4 del resalto hidralico en la poza de disipacin

Figura 6. Seccin longitudinal de la poza

Estableciendo balance de energas (E) al inicio de la trayectoria y al pie de la poza, puntos 2 y 3:

;hf prdida de energa

( Ecuacin 1 )

Estableciendo balance de energas entre los puntos 4 y 5, y considerando un %K de ahogamiento para

asegurar un resalto sumergido dentro de la poza, se tiene:

( Ecuacin 2 )

Aplicando la ecuacin de momentum entre los puntos 3 y 4, lugar donde se desarrolla el resalto
hidrulico, se tiene:
 ( Ecuacin 3 )

Por otro lado el ancho de la poza B, debe ser menor que el ancho que alcanza un chorro natural en
pendiente a una distancia igual a la longitud total de la poza, es decir:

( Ecuacin 4 )

( Ecuacin 5 )

Ordenando y simplificando las Ecuaciones 1, 2, 3 y 5 asimismo igualando las cotas C3 = C4, se obtiene
un sistema de 3 Ecuaciones No lineales con tres variables (y3, y4, B):

K 30% ;% de ahogamiento del resalto hidrulico en la poza

Fo FrN
m
v2 VN 10.851
s
y2 YN 0.207 m

Z 0 ;Talud de la poza para seccin rectangular

Z1 2 ;Talud de la rampa de ingreso a la poza

Valores iniciales para iniciar el clculo:

y3 0.3m ;Tirante conjugado menor

y4 1 m ;Tirante conjugado mayor

B 0.7m ;Ancho de la poza

Dado

Q
2
Q
2

( B Z y4) y4 y3 1.1 ( B Z y3) y3 C y5
2 2
v5 v2
( 1 K) y4 5 C2 y2 =0
2g 2g 2g 2g

Q Q
Q Q
( B Z y4) y4 B

Z y4 2 ( B Z y3) y3 B Z y3 2
y4 y3 =0
g 2 3 g 2 3

Q
2

v5
2
( B Z y4) y4
C C5 y5 ( 1 K) y4
2
2
2 g 2g
B br 2 1 Z1 =0
3 Fo
 0.137
Y Find( y3y4B) Y 1.927 m

1.378

y3 Y0 0.137 m ;Tirante conjugado menor

y4 Y1 1.927 m ;Tirante conjugado mayor

B Y2 1.378 m ;Ancho de la base de la poza, calculado

2.5
m Q
Bp 18.78 2.279 m ;Bp, ancho de la base de la poza, calculado por la frmula
0.5 3
s m emprica
10.11 Q
s

B 1.50 m ;ancho de poza asumido

Reclculo de los tirantes conjugados, con el ancho de poza asumido B:

Dado

Q
2
Q
2

( B Z y4) y4 y3 1.1 ( B Z y3) y3 C y5
2 2
v5 v2
( 1 K) y4 5 C2 y2 =0
2g 2g 2g 2g

Q Q
Q Q
( B Z y4) y4

B Z y4 2 ( B Z y3) y3 B Z y3 2
y4 y3 =0
g 2 3 g 2 3

0.126 m
Y Find( y3y4) Y
1.843
y3 Y0 0.126 m ;Tirante conjugado menor recalculado

y4 Y1 1.843 m ;Tirante conjugado mayor recalculado

Q m
v3 11.869
( B Z y3) y3 s

Q m
v4 0.814
( B Z y4) y4 s

Cota de fondo de la poza (C3 )

v5
2 v4
2
C 3 C 5 y5 ( 1 K) y4 2820.183 m
2 g 2g
Profundidad de la poza(h)

h C 5 C 3 1.417 m

Nmero de Froude (F3)

v3
F3 10.661
g y3
 10. Longitud de la poza (Lp)

Lp 5 ( y4 y3) 8.585 m ;Long. resalto segn Siechin

y3
0.1 y3 160 tanh
F3 12 ""
if 8.346 m ;Long. resalto segn Hager (1990) 2<F1<16
B 20
11. Longitud de transicin de salida (LTs)

B b5
LTS 0.677 m

2 tan12.5


180

12. Verificacin del funcionamiento de la poza

N 1 2 10

QN
( 11 N ) Q
p C2 C3 yo
v2
2
r
1.1 Q2 p
2
2 p
3
27 r
10 2g 2 AA BB
2 g B 3 27

BB
2
TH acos AA TH 4 p
Y3 2 cos

AA
3 3 3 3

3

QN QN

Q QN
N ( B Z y4) y4 B Z y4 B Z Y3N Y3N B Z Y3N 2
Y3N y4
2
Y4 N root y4
g 2 3 g 2 3

TW4N C3 Y4N

2

3
( b5 z5 y5) y5
( b5 z5 y5) y5 S5
b5 2 1 z52 y5
Y5 N rootQN
y5
n5
TW5N C5 Y5N

Verificacin N if TW4N TW5N "Ok" "----"
N QN Y3 N Y4 N TW4N TW5N 0
3
1 2.25 m 0.126 m 1.843 m 2822.027 m 2822.532 m 0 "Ok"
2 2.025 s 0.114 1.753 2821.936 2822.461 1 "Ok"
3 1.8 0.101 1.656 2821.839 2822.389 2 "Ok"
4 1.575 0.088 1.553 2821.736 2822.315 3 "Ok"
5 1.35 0.076 1.441 2821.624 2822.240 Verificacin N 4 "Ok"
6 1.125 0.063 1.319 2821.502 2822.162 5 "Ok"
7 0.9 0.05 1.183 2821.366 2822.080 6 "Ok"
8 0.675 0.038 1.028 2821.211 2821.994 7 "Ok"
9 0.45 0.025 0.842 2821.026 2821.899 8 "Ok"
10 0.225 0.013 0.598 2820.782 2821.790 9 "Ok"
 2823

2822.5

TW4 2822
N

TW5
N
2821.5

2821

2820.5
0 1 2 3
Q
N

Resultado del diseo: