UNIVERSIDAD DE ORIENTE

NÚCLEO DE ANZOÁTEGUI
ESCUELA DE INGENIERÍA Y CIENCIAS APLICADAS
DEPARTAMENTO DE INGENIERÍA CIVIL

CURSO DE HIDRÁULICA DE
CANALES ABIERTOS

Realizado Por:

ING. LUIS BELTRÁN GONZÁLEZ G.

Puerto La Cruz; julio de 2023

 ii

Primera Edición agosto de 2006

Impreso por Ameco
Telf. (0281) 2675447- 2672769

Primera Reimpresión Corregida febrero de 2013

Impreso por Techcolor Graphics, C.A.
Telf. (0281) 2675447-2672769

Prohibida la Reproducción parcial o total

de esta obra sin la autorización previa del Autor
LF: 589 2000 6201056
ISBN: 980-234-104-5
 iii

PROLOGO

El origen de este trabajo, se remonta al año de 1982, año en el que recién ingresado
como Instructor a la Universidad de Oriente, preparé unos apuntes personales para dictar la
asignatura Ingeniería Hidráulica en el Departamento de Ingeniería Civil. Dichos apuntes
fueron mejorados con la experiencia docente y con la realización de un Post-Grado en Obras
Hidráulicas desarrollado en 1989 en la Universidad de los Andes. Por otra parte, la
complejidad, lo disperso de la información y el alto costo que presentan textos reconocidos
como Rouse, Chow, French y Aguirre en Venezuela me condujeron a preparar estos apuntes
adaptados al contenido programático de la asignatura Ingeniería Hidráulica en el
Departamento de Ingeniería Civil.

El presente trabajo titulado “Curso de Hidráulica de Canales Abiertos”, no pretende

competir con excelentes textos como los mencionados anteriormente, simplemente es el
resultado de un esfuerzo personal el cual ha dado lugar a un material apropiado para los
estudiantes de Ingeniería Hidráulica de la Carrera de Ingeniería Civil en la Universidad de
Oriente.

También quiero agradecer las recomendaciones y sugerencias permanentemente

brindadas a este material de apoyo por los colegas profesores Hilda Morales y Munir Bou
Ghannan, mil gracias.

Por último, este material fue revisado y mejorado en febrero del 2021.

González Luis Beltrán

 iv

ÍNDICE
Pág.

PROLOGO iii

CAPITULO I. Introducción 1

CAPITULO II. Flujo Uniforme y Permanente en Canales 2

2.1.- Características y resistencia del flujo uniforme en canales abiertos 2

2.2.- El Cálculo del Flujo Uniforme 6

2.3.- Diseño Optimo de Canales Rectangulares 10

2.4.- Diseño óptimo de canales trapezoidales. Otras secciones 11

2.5.- Energía Específica en Canales Rectangulares 14

2.6. Energía Específica en Canales de Cualquiera Geometría. Sección Irregular 17

2.7.- Ocurrencia del Flujo Crítico y Secciones de Control 22

2.8.- Flujo en Canales con Rugosidad Compuesta 26

CAPITULO III. Diseño de Canales 30

3.1.- Tipos de Diseño 30

3.2.- Canales Revestidos 31

3.3.- Canales no revestidos 42

3.4.- Proyecto de un canal 57

CAPITULO IV. Orificios y Vertederos 60

4.1.- Orificios 60

4.2.-Vertederos 78
 v

CAPITULO V. Flujo Gradualmente Variado en Canales 86

5.1 Introducción. 86

5.2 Curvas de Remanso y Tipos de Perfiles. 86

5.3 Ecuación Dinámica del Flujo Gradualmente Variado. 92

5.4 Integración de la Ecuación diferencial del flujo gradualmente variado. 98

5.5. Métodos numéricos para el cálculo del flujo gradualmente variado 101

5.6. Resalto Hidráulico. 106

5.7. Aporte de Caudal hacia un Canal. 133

Bibliografía 136

Apéndice 1. Programa en Basic para el Cálculo de Curvas de Remanso 138

 1

CAPITULO I
INTRODUCCIÓN

El flujo con superficie libre ha sido tradicionalmente conocido como flujo en canales
abiertos, aunque también ocurre en conductos cerrados parcialmente llenos como lo son las
cloacas, en ambos casos sobre la superficie del agua está actuando la presión atmosférica y el
movimiento del fluido lo gobierna la pendiente del fondo del canal. Aunque la mayoría de las
ecuaciones hidráulicas se obtienen de la experimentación no está demás decir que se
fundamentan en leyes físicas tales como: La ecuación de cantidad de movimiento, energía y
conservación de volumen.
En este curso se analizan los casos de flujo permanente en sus modalidades más
usuales: Uniforme, gradualmente variado y rápidamente variado (orificios, vertederos,
Resalto Hidráulico y transiciones).
El Capítulo II presenta los fundamentos hidráulicos básicos en flujo uniforme y
permanente para diseñar canales y el capítulo III explica los métodos para diseñar los mismos
ya sea el canal revestido o no revestido.
La intención del Capítulo IV es presentar un resumen de manera sencilla de la forma
como determinar el caudal que escurre por un orificio o un vertedero que son los tipos de
medidores de flujo más usuales en canales, haciendo más énfasis en los orificios de pared
delgada contracción total y descarga libre y en los vertederos de pared delgada y descarga
libre por ser estos los más usuales.
No siempre el flujo en un canal es uniforme y permanente, tal vez con el estudio
hidrológico se podría estimar un flujo permanente, pero habrán situaciones tales como: La
colocación de un medidor de flujo, producir un resalto hidráulico, desarrollar una transición
entre dos secciones diferentes o simplemente un cambio de pendiente que originarán tramos
en que el flujo deje de ser uniforme, por ese motivo el Capítulo V presenta la manera de
como estudiar el flujo gradualmente variado y en conjunto con el Capítulo IV la manera de
como estudiar el flujo rápidamente variado, para así poder desarrollar lo que se conoce como
el "perfil de continuidad del fluido" que debe ser un requisito indispensable a la hora de
realizar un proyecto hidráulico de esta naturaleza.
 2

CAPITULO II
Flujo Uniforme y Permanente en Canales
2.1.- Características y resistencia del flujo uniforme en canales abiertos

Según el volumen de control señalado en la figura 2.1 en flujo uniforme y permanente,

las únicas fuerzas actuantes y que completan el equilibrio son: la componente del peso en la
dirección del movimiento ALSen, y la de rozamiento entre el fluido y el contorno liso sólido
0PL, donde "0" es el esfuerzo cortante y "P" el perímetro mojado. La ecuación de equilibrio
es:

 0 PL  ALSena  Si R  A / P  Radio Hidráulico   0  RSena (2.1)

energía total

Fig. 2.1 Esquema de definición correspondiente a la ecuación de resistencia

En una tubería a sección llena R = D/4
L V2
como hf  LSena  f y reemplazando Sen de la Ec (2.1) en esta ecuación
4R 2 g
obtenemos:
V2
0   f (2.2)
8
dónde :  = /g
Igualando (2.1) y (2.2)

8g
V  RSena ; Sena  S 0 (2.3)
f
 3

La Ec. (2.3) se conoce como Ecuación de Chezy y se puede expresar de la siguiente manera:

(2.4)
C según Manning y multiplicando la Ec. (2.4) por el área se obtiene la Ecuación de
Manning
(2.5)
donde:
Q = Caudal (m3/seg)
n = Factor de Fricción de Manning (Tabla 2.1)
A = Área mojada de la sección transversal (m2)
R = Radio Hidráulico (A/P) (mts)
S0 = Pendiente del fondo del canal

Tabla 2.1 Valores de la Rugosidad n de Manning

n Superficie

0,01 Muy lisa; vidrio; plástico, cobre

0,011 Concreto muy liso

0,013 Madera suave, metal, concreto

0,017 Canales de tierra, buenas condiciones

0,02 Canales naturales de tierra, libres de vegetación

0,025 Canales naturales con alguna vegetación y piedras esparcidas en el fondo

0,035 Canales naturales con abundante vegetación

0,04 Arroyos de montaña con muchas piedras

Nota: En un canal poco profundo y muy ancho la altura del agua se puede considerar igual al radio hidráulico.

Si utilizamos la ecuación de Chezy (2.3 y 2.4) observamos que aparece el factor de fricción f
de la ecuación de Darcy - Weisbach el cual está en función de la rugosidad relativa (/D) y
del número de Reynolds (Re); si definimos la velocidad de corte como U   0 y

*
sustituyendo  0 de la Ec. 2.2 y V de la Ec. 2.3 tenemos:

0 f
U*  V  gRS 0 (2.6)
 8

U *
 Re*  (2.7)

 4

 
R e  2000  f  64 / Re  (2.8)
Zona de Flujo   )))))0))
 
si Re*  4  Hidráulicamente  2000  Re  10  f  0,316 / Re 
5 1 4
))))
(2.9)
liso  
 R  10 5  1  2 log 
Re f  

 e  2,51   (2.10)
 f   

Zona de
 
si 4  Re*  100  Régimen de   1  2 log 
   2,51 
  12 R  (2.11)
transición 
 f  Re f 


Zona de Régimen 
 1  12 R 

si Re*  100  Completamente   2,03 log   (2.12)

 f   

Rugoso

Donde: VD V (4 R) VD
Re    Re*
v v U *

U* = velocidad de corte
Re* = Nº Reynolds de Corte
Re = Nº Reynolds
R = Radio Hidráulico
 = Rugosidad Absoluta del fluido
 = densidad del fluido
v = Viscosidad Cinemática del fluido

Nota: Estas ecuaciones sirven tanto para tuberías o canales que fluyan por gravedad; si se devuelve el cambio R = D/4
servirán únicamente para tuberías a presión.

Las ecuaciones 2.8, 2.9, 2.10, 2.11 y 2.12 se encuentran representadas en el diagrama de
Moody (Fig. 2.2)
 5

Fig. 2.2 Diagrama de Moody

 6

2.2.- El Cálculo del Flujo Uniforme

La fórmula de Manning se ha convertido en el principal instrumento para el cálculo

práctico de canales.
En la figura 2.3, presentada por Chow (1994) en 1959 y en el supuesto de tenerse como
dato la relación y/b y el talud, o la relación y/d si la forma a adoptar es la circular, el factor
de forma correspondiente permite determinar el caudal de circulación en flujo uniforme ya
que:

AR 2 / 3 Qn AR 2 / 3 Qn
8/3
 y 
b ( S 0 )1 / 2 b 8 / 3 d 8/3
( S 0 )1 / 2 d 8 / 3 (2.13)

Fig. 2.3 Curvas para la determinación de la profundidad normal (Aguirre, 1993)

Ejemplo 2.1 Un canal de concreto, n = 0,013, tiene una anchura de 3 m. en el fondo con
unos taludes de pendiente 4 vertical a 1 horizontal. Su pendiente es de 3,8 m. por Kilómetro.
¿Cuál será la profundidad de flujo cuando el caudal sea de 5 m3/s?
 7

De la Ecuación de Manning

Aplicando el Método de Newton  y = 0,58 m

otra forma es

Qn 5 x0,013 y
1/ 2 8 / 3
 1/ 2 8/3
 m 
 0,0563 fig2.3   0,195  y  0,58m
0.25

(S0 ) b (0,0038) (3) b

Ejemplo 2.2 Se aforó un río ancho y poco profundo y tenía una velocidad media igual a 0,4
m/s y la Prof. y = 1 m. la pendiente del río = 0,000064 en régimen uniforme y permanente
¿Determinar n?

1 1
Q  VA  AR 2 / 3 ( S 0 )1 / 2  0,4  (1) 2 / 3 (0,000064)1 / 2  n  0,02
n n

Ejemplo 2.3 Para una profundidad de 2m., las características de un río son: A=60m2, P=90 m,
V = 0,95 m/s y S0 = 0,003 en régimen uniforme y permanente. Calcular Q para y = 3m, A =
120 m2 y el P = 120 m del mismo río.

2/3
1  60 
0,95 * 60  (60)  (0.003)1 / 2  n  0,04399
n  90 

Secciones Circulares
El cálculo de canales de sección circular, o de tuberías en régimen de superficie libre se
refiere con alguna frecuencia a la condición de sección llena. Si se supone que al variar la
profundidad del flujo uniforme, en un canal de sección circular no varía la rugosidad, entonces
la ecuación de Manning puede dividirse por sí misma para la condición de sección llena,
obteniéndose que:

(2.14)

(2.15)
 8

Graficando las ecuaciones (2.14) y (2.15) se obtiene la figura 2.4

Fig. 2.4 Características de flujo en una sección circular (Aguirre, 1993)

Con el cálculo del caudal y la velocidad para la sección llena, la Fig. 2.4 permite
encontrar las correspondientes cantidades en una sección parcialmente llena. Los niveles a los
cuales se presenta el caudal y la velocidad máxima corresponden a y=0,94d e y=0,81d
respectivamente.
Ejemplo 2.4 Por una alcantarilla de drenaje fluye un caudal de 56 l/s en flujo uniforme. La
alcantarilla tiene una pendiente de 1 por mil, un diámetro de 45 cm. y un valor n= 0,015.
¿Calcular el nivel y la velocidad?
Este ejemplo se puede resolver en forma similar al ejemplo anterior por medio de la Fig.
2.3 y 2.5, pero con el fin de mostrar la utilidad de la Fig. 2.4 haremos uso de ella en la solución:
1 x0,45 2  0,45 
2/3

Qlleno    (0,001)1 / 2  0,079m 3 / s

0,015 4  4 
Q 56 y
  0,719 Figura
 2.4
  0,63  y  28,4cm
Qlleno 78 d

V  1  0,45  2 / 3 
De la misma figura 2.4:  1,1  V  1,1   (0,001)   0, 54 m/s
1/ 2

Vlleno  0,015  4  
 9

Otras Secciones
Cualquier sea la geometría de la sección transversal la ecuación de Manning gobierna
el movimiento del fluido por esta razón la Figura 2.5 tomada de Bolinaga (1979) recoge los
elementos geométricos de las secciones más usuales.

Fig. 2.5 Elementos geométricos de las secciones más usuales.

Ejemplo 2.5 Una alcantarilla circular de concreto con n=0,014, tiene un diámetro de 1,8
m, cuando la profundidad es de 1,1 la alcantarilla conduce un caudal de 2m 3/s en flujo
uniforme y permanente.
Determinar: So, V y Q llena.

y 1,1 Qn
  0,61   1
Figura2.3
8
 0,21 So  7,7328x10 4
D 1,8 So 2 D 3
0.2
Cos    77.16   205,67
0.9
  180
x  205,67 } x  3,589rad  A 
1
8
(3,589  sen 205,67) 1,8 2  1,63m 2
 10

2
V  1,226m / s
1,6289

1 
2/3

Qlleno  1,82  1,8  7,73x10 

4
1
2
 2,96m 3 s
0,014 4  4 

2.3. Diseño Óptimo de Canales Rectangulares

En un canal para un caudal, rugosidad y pendiente dada se puede demostrar a partir de
la fórmula de Manning que cuando la sección transversal es un mínimo, el perímetro mojado
también es un mínimo (Ec. 2.16) y entonces la excavación y el recubrimiento se aproximan a
su mínimo costo para las mismas dimensiones del canal. Entonces:

1 2/3 QnP 2 / 3 CP 2 / 3
Q AR ( S 0 )1 / 2  A  1/ 2 2 / 3
 2/3
 A 5 / 3  CP 2 / 3  A = C'P2/5 (2.16)
n (S 0 ) A A

P = b+2y (2.17)

A = by (2.18)
Donde C'

Sustituyendo b de (2.17) en (2.18) y (2.18) en (2.16)  (P - 2y)y = C'P2/5 (2.19)

Para hallar el valor de y que hace el Pmin se obtiene = 0 de la Ec.(2.19)

0 0
 dP  2 dP
entonces   2  y  ( P  2 y )  C ' P 3 / 5  (2.20)
 dy  5 dy

b=2y (2.21)

A=2y2 (2.22)

R=y/2 (2.23)
 11

2.4 Diseño óptimo de canales trapezoidales. Otras secciones

m x
T = b+2my Tgα = m    x  my;
1 y

por Pitágoras (𝒎𝒚)𝟐 + 𝒚𝟐 = 𝒛𝟐 => 𝒛 = 𝒚√𝟏 + 𝒎𝟐 (2.24)

P = b + 2z = (2.25)

(2.26)
Donde m = máx pendiente para la cual la tierra húmeda no se desmorona

Repitiendo los pasos de canales rectangulares: A = C' P2/5


b de (2.25) en (2.26) (2.26) en (2.16)  P  2 y 1  m y  my  C ' P
2
 2 2/5
(2.27)

P
Haciendo m constante el valor de y que hace P min se obtiene de  0 de la Ec. (2.27)
y
obteniéndose:
0 0
 P 
  2
  2 1  m2  y  P  2 y 1  m2  2my  C' p 3 / 5
 y  5
P
y
 (2.28)

La Ec. (2.28) en (2.25)  (2.29)

La Ec. (2.29) en (2.26)  (2.30)

Dividiendo (2.30) entre (2.28)  (2.31)

P
Haciendo y constante el valor de m que hace Pmin se obtiene de = 0 de la
m
Ec. (2.27) obteniéndose:
 12

0 0

 P 2 y  2 P 2m 3
  (1  m 2 ) 1 / 2 2m  y  y 2  C' P 3 / 5  1 m   Tg 30º
 m 2  5 m 1  m2 3

3
Sustituyendo m = en las Ec. (2.28), (2.29) y (2.30) se obtiene el perímetro mínimo para
3
el talud óptimo, lo que se conoce como máxima eficiencia hidráulica donde se puede inscribir
un semicírculo y dicho trapecio tiene lados iguales (semihexágono) como se puede ver en la
ecuación (2.33).

(2.32)
(2.33)
(2.34)
(2.35)

Ejemplo 2.6 En un canal muy liso; n=0,01; trapezoidal, con pendiente del 1 por mil y cuyos
taludes están en la relación de 2 horizontal a 1 vertical, trabajando en régimen uniforme bajo
condiciones hidráulicas óptimas, se determina el esfuerzo cortante (esfuerzo tractivo)
promedio o= 0,5 x 10-4 Kgf/cm2

a.- ¿Cuáles son las dimensiones b, y del canal?

b.- ¿Cuál es el gasto?
c.- ¿Cuál es la C de Chézy?
d.- ¿Cuál es el factor de fricción f (de Darcy) correspondiente?
e.- ¿Cuál es el diámetro  de los granos de arena (en el diagrama de Moody) que producen
una rugosidad equivalente?
τ 0,5 y
a) 0 = RS0 R 0   0,5m ; R   y  1m
 S 0 1000 x0,001 2

b = 2(1) 1  2 2  2 x2 x1  0,47m

 
2/3
1 1 1
b) Q  AR ( S 0 )  2(1) 2 1  2 2  2 x12   (0,001)1 / 2  4,92m 3 / s
2/3 1/ 2

n 0.01 2
 13

c)

8 g 8 x9,81
d) f    0,0098
C 2 (89,1)

 3,98x10 6 y f = 0,00988 en Moody    0,00001

Q(4 R) 4,92(4 x0,5)
e) Re   6
A v 2,47 x10 D

Por lo tanto  = (4R)0,00001 = (4x0,5)x0,00001 m = 0,02 mm

Verificación

* 
0,5
102
 0,07 m / s } V 
Q 4,92

A 2,47
 1,99m / s

*  gx0,5 x0,001  0,07m / s

f 8 2 8 x0,07 2
*  v  f  2x   0,0098
8 V 1,99 2

 *  0,07 x0,00002
Re *    1,4
v 10 6

Re 
V (4 R) 1,99 x 4 x0,5
10 6

10 6
 3,98 x10 6 } 1f  21g Re2,51f 
  3,98x10 6 f  
21g  1  0
  2,51  f 
 
Newton Raphson : f  0,0093

Otras Secciones
Los elementos geométricos para seis secciones hidráulicas óptimas se muestran en la
Tabla 2.2, pero no siempre estas secciones son prácticas, debido a dificultades en la
construcción y en el uso de material. En general, una sección de canal debe diseñarse para
cumplir con una eficiencia hidráulica óptima, pero debe modificarse para tener en cuenta
aspectos constructivos. Desde un punto de vista práctico, nótese que la sección
hidráulicamente óptima es la sección que da el área mínima para un canal determinado, pero
no necesariamente la mínima excavación.
 14

Tabla 2.2. Secciones Hidráulicas Optimas (Chow, 1994)

Perímetro Radio Ancho Profundidad Factor de

Sección Área
Mojado Hidráulico Superficial Hidráulica Sección
Transversal A
P R T D Z

Trapecio, medio 1 4 3 3 2, 5
hexágono 3 y2 2 3y y 3y y y
2 3 4 2
Rectángulo, 1
2y2 4y y 2y y 2y2,5
medio cuadrado 2
Triángulo, medio 2 1 1 2 2,5
y 2 2y 2y 2y y y
cuadrado 4 2 2

Semicírculo  2 1   2,5
Parábola
y y y 2y y y
2 2 4 4

8 1 2 8
T  2 2y 4/3 2y2 2y 2 2y y 3 y 2,5
3 2y 3 9
Catenaria
1,39586y2 2,9836y 0,46784y 1,917532y 0,72795y 1,19093y2,5
hidrostática

2.5.- Energía Específica en Canales Rectangulares

Se define la energía específica E como la energía del agua sobre el fondo del canal, la
cual se mantiene constante siempre y cuando el flujo sea uniforme y permanente. Si definimos
el caudal por unidad de ancho q= Q/b y recordamos que V= q/y, la energía específica se puede
escribir de la forma:
(2.36)
Manteniendo q constante y dando diferentes valores a y obtenemos la figura 2.6 (a)

Figura 2.6
 15

La figura 2.6 (a) muestra que existen dos profundidades para un valor de energía las
cuales reciben el nombre de profundidades alternas. yc es la profundidad que determina si el
régimen es tranquilo o rápido; el cambio de un flujo tranquilo o rápido o viceversa implica un
flujo no uniforme y permanente. Para pequeños valores de "y" la curva se hace tangente al eje
"E" en el infinito, mientras que para grandes valores de "y" la altura de velocidad es
despreciable y la curva se aproxima a la línea de 45º donde E=y.

El valor de la profundidad que hace Emin, se llama profundidad crítica y se obtiene de

hallar dE/dy = 0 en la Ec. (2.36).

dE 1  q2 2y  q2
 1    1  0q  gy 3 C (2.37)
dy 2g  y 4  gy 3

(2.38)

Eliminando q entre (2.36) y (2.37) (2.39)

Como q = vy y sustituyendo en la Ec. (2.37) tenemos (2.40)

De la Ec. (2.40) (flujo critico F=1) y tomando en cuenta la fig 2.6 (a)
tenemos:
Condiciones:

V y
a) Si  1ó  1  Régimen tranquilo o subcrítico
yc g yC

V  y
b) Si 1ó 1 Régimen rápido o supercrítico
yc g yC

Alternativamente si despejamos q de la Ec. (2.36) tenemos que:

Q
 q  2g( y 2 E  y3 ) (2.41)
b

Manteniendo E constante y graficando y vs. q se obtiene la Fig. 2.6 (b) donde se observa
que la profundidad para el gasto máximo se llama profundidad crítica y se obtiene de hacer
dq/dy=0 en la Ec. (2.41) con E constante.
 16

dq 2
 2 g (2 Ey  3 y 2 )  3 y 2  2 Emin y  y  Emin  yC
dy 3 (2.42)

Lo que nos dice que qmax se obtiene cuando la profundidad normal es igual a yc

Los gráficos de la figura 2.6 se pueden generalizar y hacer adimensionales. Si dividimos

la Ec. (2.36) por yc de la Ec. (2.38) tenemos la Ec. (2.43), graficada se obtiene la Fig. 2.7.

(2.43)

Fig. 2.7. Diagrama Adimensional de Energía Específica

Ejemplo 2.7 En un canal rectangular Q= 7,3 m3/s. b= 9 m, y= 0,9 m. ¿Hallar la energía

específica y el tipo de régimen?

2
1  7,3 
E  0,9     0,94m
2 g  9 x0,9 

1/ 3
q1 / 3  (7,3 / 9) 2 
yC      0,40m  y  Régimen Subcritico
g  9 x 0,9 
 17

Ejemplo 2.8 En un canal rectangular Q= 5,4 m3/s, b= 3,6 m, n= 0,019. Hallar:

a) yc= ? b) Vc= ? c) ¿Qué pendiente producirá la velocidad crítica?

VC  9,81x0,612  2,45m / s

2/3
1  3,6 x0,612 
2,45    ( S C )1 / 2  S C  0,0062
0,019  (3,6  2 x0,612) 

Ejemplo 2.9 En un diseño óptimo de un canal rectangular se observa que el régimen es crítico
en flujo uniforme y permanente. n= 0,02, Vc= 2 m/s. ¿Hallar la pendiente de la línea
piezométrica?

V 2c 22
Vc  gyC  yC    0,40m
g 9,81

2
 V n2 2 / 3 
2
 2 x0,02 x 2 2 / 3 
2/3
1 y 
VC   C   S C   C 2 / 3      0,0137
1/ 2
SC
n 2   yC   0,4
2/3


Sc = 0,0137= 1,37%

2.6. Energía Específica en Canales de Cualquiera Geometría. Sección Irregular

Los conceptos desarrollados de energía específica en canales rectangulares se pueden
generalizar para cualquiera geometría, así la ecuación (2.36) se puede escribir como:

(2.44)

Manteniendo Q constante y dando diferentes valores a y obtenemos la Fig. 2.8 (a).

 18

Fig. 2.8

En una sección cualquiera el valor del área que hace la Emin, se llama área crítica,
entonces haciendo dE/dy=0 en la Ec. (2.44) y tomando en cuenta la Fig. 2.9 tenemos:

Fig. 2.9 Sección transversal de cualquiera Geometría

Como dA= tcdy y sustituyendo en la ecuación anterior tenemos: (2.45)

 19

dividiendo la Ec. (2.45) entre AC2 

(2.46)

De la Ec. (2.46) (flujo critico F=1) y tomando en cuenta la figura 2.8

(a) tenemos las siguientes condiciones:

V y
a) Si 1 ó  1  Régimen tranquilo o subcrítico
gAC / Tc yc

V y
b) Si 1 ó  1  Régimen rápido o supercrítico
gAC / Tc yc

Q  VC AC  AC gAc / Tc
(2.47)

VC2 A
E min  y c   yc  c (2.48)
2g 2TC

Esta ecuación según la figura 2.8 (b) muestra que el mínimo de energía se presenta
cuando la carga de velocidad es la mitad del tirante promedio A/T.
Para el cálculo de profundidades críticas en canales de cualquiera geometría, es
frecuente el empleo del factor de sección FS definido como el producto del área y la raíz
cuadrada de la profundidad hidráulica media

A
FS  A y m  A (2.49)
T

Si esta expresión se sustituye en la Ec. (2.45) se obtiene para un flujo critico (F=1) que

(2.50)

Las relaciones (2.49) y (2.50) son de gran utilidad en el cálculo y análisis del flujo crítico
en canales abiertos. Para un caudal dado, las ecuaciones anteriores dan el factor de sección
crítico FSc y por consiguiente la profundidad crítica yc. En la figura 2.10 tomada de
 20

Chow (1994), se presentan con el fin de facilitar el cálculo, curvas adimensionales que dan la
relación entre la profundidad y el factor de sección FS para canales rectangulares,
trapezoidales y circulares.

Valores de FS/d2.5 para canales trapezoidales

Valores de FS/d2.5 para secciones circulares

Fig. 2.10. Curvas para la Determinación de la Profundidad Crítica (Aguirre, 1993)

Ejemplo 2.10 Calcular la profundidad y velocidad crítica en un canal trapezoidal que

tiene una base b= 6m y taludes 2H: 1V, para un gasto Q= 12000 l/s.
 21

Q 2 AC3 Q 2T 2 12 2 x(6  2 x2 xyC ) Aplicando

   1  0  1  0  yC  0,684m
g TC gAC3
9,81(6 yC  2 yC 2
Newton  Raphson

AC  6 x,684  2 x0,684 2  5,039m 2 TC  6  4 x0,684  8,736m

𝐴 5,039
𝑉𝐶 = √𝑔 𝑇𝐶 = √9,81𝑥 = 2,378𝑚/𝑠
𝐶 8,73

Otra forma:
Q 12 FS 3,83 m2 y
FS    3,83  2,5  2,5  0,0434  C  0,115
g 9,81 b 6 Fig .2.10 b

yc  0,115x6  0,69m

12
VC   2,38m / s
5,039

Sección Irregular
A) Si en la sección transversal no conocemos la profundidad crítica.

1) Se calculan las áreas y anchos superficiales para diferentes prof. en la sección

Y A T Fs  A A / T
Y1 A1 T1 Fs1
.
.
.
.
Ymax A Total T Fs
 22

2) Se gráfica Y Vs. FS

3) Como Q es un dato se calcula el FS= Q 2 / g y en la figura anterior se halla el valor

de yc y luego se pueden utilizar las ecuaciones (2.45), (2.46), (2.47) y (2.48)
Nota: Si la sección es compuesta el procedimiento sería el mismo.

2.7.- Ocurrencia del Flujo Crítico y Secciones de Control

El control de flujo en un canal abierto se define de muchas maneras. Tal como se utiliza
aquí, el término significa el establecimiento de una condición definitiva de flujo en un canal
o, más específicamente, una relación definitiva entre el nivel y el caudal del flujo. Cuando el
control de flujo se alcanza en una cierta sección de canal, esta sección es una sección de
control. Más adelante se demostrará que la sección de control regula el flujo de tal modo que
restringe la transmisión de efectos de cambios en la condición del flujo, ya sea en una
dirección hacia aguas arriba o en una dirección hacia aguas abajo, según el estado de flujo en
el canal. Como la sección de control mantiene una relación nivel-caudal definitiva, siempre
es un lugar adecuado para una estación de aforo y para el desarrollo de una curva de
calibración de caudales, la cual es una curva que representa la relación profundidad-caudal en
la estación de aforo.

En el estado crítico de flujo puede establecerse una relación definitiva nivel-caudal, y

representarla mediante la ecuación (2.45). Esta ecuación muestra que la relación nivel-caudal
es teóricamente independiente de la rugosidad del canal y de otras circunstancias no
controladas. Por consiguiente, una sección de flujo crítico es una sección de control.
 23

La localización de la sección de control en un canal prismático por lo general está

gobernada por el estado de flujo, el cual a su vez se determina mediante la pendiente del canal.
Tómese como ejemplo un canal prismático recto largo en el cual se crea un embalse mediante
una prensa a través del canal y el agua fluye por encima de la prense a través de un vertedero
de rebose (Figura 2.11). Se muestran tres condiciones de flujo en el canal, que representan los
flujos subcrítico, crítico y supercrítico, respectivamente. Las pendientes del canal en los tres
casos son, suave o subcrítica, crítica y empinada o supercrítica, de manera equivalente

Figura 2.11. Condiciones de Flujo en un Canal Pismático Largo

Si el canal tiene una pendiente crítica (esquema intermedio de la figura 2.11), entonces
en principio el flujo es uniforme y crítico a través del canal. Sin embargo, debido a la presencia
de la presa, el flujo a través del embalse será subcrítico y la superficie del embalse será casi
horizontal. En el extremo de aguas abajo se desarrollará una curva de caída, que se extiende
aguas arriba desde una sección cercana a la cresta del vertedero y se vuelve asintótica al nivel
del embalse.
 24

Si el canal tiene una pendiente subcrítica (esquema superior de la figura 2.11) en

principio el flujo es subcrítico. Debido a la presencia de la presa, el nivel del embalse se
extenderá aguas arriba a lo largo de una gran distancia formando una curva conocida como
curva de remanso. La profundidad de agua adicional se requiere para obtener la altura
suficiente con el fin de alcanzar el incremento de velocidad necesario para pasar el agua por
encima del vertedero. Este efecto de represar el agua por detrás de la presa se conoce como
efecto de remanso. En el extremo de aguas abajo la curva de remanso se conecta con una curva
de caída suave, que lleva el agua por encima del vertedero.

Si el canal tiene una pendiente supercrítica (esquema inferior de la figura 2.11), en

principio el flujo es supercrítico. Debido a la presencia de la presa, el efecto de remanso
originado en el embalse no se extiende muy lejos aguas arriba. Por el contrario, el flujo en el
canal de aguas arriba continuará en la dirección de aguas abajo en estado supercrítico, hasta
que el perfil de la superficie de flujo se encuentre por debajo del nivel del embalse; luego
aumentará abruptamente hasta la elevación del embalse a través de un resalto hidráulico. El
efecto de remanso no se extenderá aguas arriba a través del resalto hidráulico. El flujo aguas
arriba del resalto está completamente gobernado por las condiciones de aguas arriba.

Los ejemplos anteriores explican la importancia del hecho de que en pendientes

subcríticas el efecto de cambio en la elevación de la superficie del agua del lado de aguas
abajo se transmite aguas arriba a través de una curva de remanso, en tanto que en pendientes
supercríticas no puede transmitirse lejos aguas arriba. La condición de flujo en un canal
subcrítico se afecta por las condiciones de aguas abajo; pero, en un canal supercrítico o en el
lugar donde el agua entra al canal, la condición de flujo depende por completo de las
condiciones aguas arriba. En consecuencia, se dice que el control de flujo se localiza en el
extremo de aguas abajo para canales con pendientes subcríticas, y en extremo de aguas arriba
para canales con pendientes supercríticas.

Cuando el canal tiene una pendiente subcrítica, una sección de control en el extremo de
aguas abajo puede ser una sección crítica, como la creada en la parte superior de un vertedero
de rebose. En una pendiente supercrítica, la sección de control en el extremo de aguas arriba
también puede ser una sección crítica, tal como se muestra en la figura. Una compuerta
deslizante, o un orificio o cualquier otra estructura de control también pueden utilizarse para
crear una sección de control. Nótese que, si la pendiente del canal es crítica, subcrítica o
supercrítica, dependerá no sólo de la medida de la pendiente real, sino también del caudal o
en la profundidad de flujo.
 25

Se puede obtener la profundidad crítica en puntos de un canal abierto, en los que el

fondo del canal se eleva por la construcción de una protuberancia baja, o en donde se estrecha
el canal moviendo hacia adentro las paredes laterales. Tales contracciones (que por lo general
tienen una elevación en el fondo del canal) se conocen en general como canalones de Venturi,
y en la medición de agua para irrigación se usan ampliamente diseños específicos de estos
canalones. Se puede hacer un análisis preliminar de tales contracciones (que presentan flujo
acelerado) suponiendo un flujo unidimensional y despreciando las pérdidas de carga, pero los
diseños finales requieren información más refinada, procedente, ya sea de otros diseños o de
experimentos de laboratorio. Las ventajas de tales canalones consisten en su habilidad para
dejar pasar agua cargada de sedimentos sin que éstos se depositen, y en el pequeño cambio
neto de nivel de agua requerido entre los canales de entrada y de salida.

Ejemplo 2.11 En un canal rectangular largo de 3 m de ancho que tiene una n de Manning de
0,015 y está dispuesto con pendiente de 0,001, ocurre un flujo uniforme a profundidad de 1,5
m. Calcúlense (a) la altura mínima de la protuberancia que se puede construir en el piso de
este canal para producir la profundidad crítica, y (b) la anchura máxima de contracción que
puede producir profundidad crítica.

Solución:

0,89
E min  3 *  1,33m
2

Una altura de protuberancia menor que x hará descender la superficie del agua sobre la
protuberancia, pero no podrá producir profundidad crítica; una altura de protuberancia mayor
que ésta, producirá profundidad crítica, siendo entonces la protuberancia un vertedor de cresta
amplia. Sin embargo, la última condición elevará la línea de energía en la vecindad de la
protuberancia y (siendo este un flujo subcrítico) aumentará la profundidad corriente arriba de
este punto. Por lo anterior, la altura mínima de la protuberancia para producir la
 26

profundidad crítica, será aquella que lo haga sin elevar la línea de energía. Del diagrama, se
ve que ésta es 1,65 - 1,33 o 0,32 m.

(b) Siendo desconocida la anchura de la contracción, b, pero estando especificada la

profundidad crítica,

De la parte (a) se puede ver que las anchuras de contracción más pequeñas que la
requerida para apenas producir la profundidad crítica, lo harán así elevando la línea de energía
y la profundidad del agua corriente arriba de la contracción. De acuerdo con lo anterior, la
condición deseada será aquella que produzca la profundidad crítica en la contracción, sin
elevar la línea de energía. Como la contracción es de sección transversal rectangular.
1,65
yC  2 x  1,1m
3

Igualando esta expresión con la precedente y resolviendo para b, b = 2,17 m.

2.8.- Flujo en Canales con Rugosidad Compuesta

En algunos canales la rugosidad puede ser diferente en distintos tramos del perímetro
de su sección transversal. Cuando se trata de aplicar la Ecuación de Manning en tales canales
es necesario determinar un valor equivalente de n para toda la sección y que tenga el mismo
efecto de los coeficientes parciales n1, n2, ...nN (Ver figura 2.12). Para encontrar el valor de n,
hay que dividir convenientemente el área del agua imaginaria en N partes de tal manera que
los perímetros P1, P2, ...Pn, correspondan a los coeficientes n1, n2, ...nN. Asumiendo una
velocidad uniforme en toda la sección igual a la media e igual a la velocidad en cada área
imaginaria:
 27

Fig. 2.12 Canales con varios revestimientos (Bolinaga. 1985)

Como:
A   A1  
n
A
3/ 2
P
Pn 
i i
A
3/ 2

n P
 P n 
3/ 2 i
3/ 2
i

n
P n 1/ 5
i i  P2 n12/ 5  ...Pn n1n/ 5 
2/3 (2.51)
P2/ 3
Si la forma de la sección transversal del canal es tal que no permita suponer que la
velocidad sea única, como por ejemplo en el caso de canales aluvionales en los cuales se
producen inundaciones en épocas de alta precipitación, entonces la rugosidad equivalente se
puede hallar haciendo la sumatoria de los gastos correspondientes a los canales parciales.
Consideremos, en la fig. 2.13, el canal compuesto por las secciones 1 y 2 con parámetro de
subíndice 1 y 2.

Fig. 2.13. Canal de Rugosidad Compuesta y Velocidad No Uniforme

De la fórmula de Manning para el gasto total Q = Q1+Q2 se obtiene que:

   
   
1A 5 / 3 i 2
 1( Ai ) 5 / 3  En General A 5/3
 1  A 5/3
 1 
     n  2 / 3 
2/3      
np 2 / 3 i 1  ni ( Pi )  p  i  N  A5 / 3   p p  i  N  A5 / 3  
5 / 3 1

 
 2 / 3   5 / 3 1 
 
 i 1  np  i   i 1  np p  i 
Simplificando
 28

PR 5 / 3
n  iN (2.52)
 PR 5 / 3 
 
i 1 

n i
Ejemplo 2.12.- Un canal trapecial de 4 m de ancho en la base y lados 1:1, tienen varios
revestimientos como se muestra en la figura. Para una altura de agua de 6 m, calcule la n
equivalente por los dos criterios antes comentados y el gasto aproximado. (S0 = 0,0012).

Parte Ai (m2) Pi (m) Ri (m) n Q0 (m3/s)

1 4,5 4,24 1,06 0,024 6,75
2 13,5 4,24 3,18 0,038 26,61
3 24 4 6 0,015 183
4 13,5 4,24 3,18 0,038 26,61
5 4,5 4,24 1,06 0,024 6,75
Totales 60 20,96 249,72

2/3
 2 x4,24 x0,0241,5  2 x4,24 x0,0381,5  4 x0.0151,5 
n     0,029
 20,96 
  60 
5/ 3

 20,96 x  
  20,96  
n 5/3 5/3 5/3   0,017
 2 x 4,24 x1,06  2 x 4,24 x3,18  4 x6 
 0,024 0,038 0,015 
 
2
1  60  3 1 3
Qaprox  (60)  (0,0012) 2  144,49 m s
0,029  20,96 
Puede notarse la diferencia entre los dos criterios, el autor recomienda según la forma
geométrica de la sección, el uso del primer criterio (velocidad uniforme en toda la sección).
Aplicando la ecuación de Manning con el segundo criterio n=0,017 el gasto es 246,49 m3/s
valor muy parecido al arrojado por la tabla.
 29

Ejemplo 2.13.- El canal indicado en la figura tiene una pendiente de fondo de 0,0012.
Calcúlese el gasto suponiendo flujo uniforme y altura normal 2,5 m.

n1,5 = no se pueden hallar por el criterio del gasto por ser una sola sección entonces:

2/3
 2 x0,041,5  1,5x0,0131,5 
n1 y 5     0,03
 3 ,5 

2/3
 2 x0,0151,5  1,41x0,041,5 
n2 y 4     0,027
 3,41 

Parte Ai (m2) Pi (m) Ri (m) n Q0 (m3/s)

1 3 3,5 0,86 0,03 3,12
2 7 3,41 2,05 0,027 14,50
3 13 4,82 2,7 0,015 58,17
4 7 3,41 2,05 0,027 14,50
5 3 3,50 0,86 0,03 3,12
Totales 33 18,64 93,41

La sumatoria de los gastos resulta en un gasto total de 93,4 m3/seg

Puede notarse la combinación de los dos criterios debido a que en la sección 1 y 5 el n1,5
no se puede hallar por el criterio del gasto por ser una sola sección, sin embargo, el gasto total
es el calculado con el criterio del gasto.
Si el canal reviste mucha importancia se recomienda investigar su rugosidad mediante
el uso de modelos físicos.
 30

CAPITULO III
DISEÑO DE CANALES

3.1.- Tipos de Diseño

El problema del diseño de un canal generalmente se presenta teniendo como datos el
gasto que debe transportar, la pendiente disponible de acuerdo con la topografía del terreno y
la rugosidad de sus paredes. Con estos datos es posible determinar, a partir de la Ec. (2.13),
un único valor para el módulo de sección AR2/3. Sin embargo, como se observa en la Fig. 2.3,
el mismo factor se puede satisfacer con distintas formas de la sección, unas más eficientes que
otras, lo que implica más de una solución.

Una de las soluciones consistiría en elegir la forma y dimensiones adecuadas que debe
tener la sección, de modo que se pueda adaptar a la topografía del terreno donde se va a
excavar el canal, y que sea lo más económica posible. Sin embargo, de acuerdo con el material
en que se excave el canal, y no existiendo revestimiento, habrá tramos en que la velocidad del
agua, erosione los taludes y la plantilla modificando la sección escogida. Por ello conviene
diferenciar entre los canales revestidos y canales no revestidos. Los primeros comprenden a
los canales que se revisten con un material resistente a la acción erosiva del agua (concreto,
mampostería, madera, plástico, etc.), o bien que se excaven en un material de iguales
características (cimentación firme, roca sana, etc.). Los segundos comprenden a los canales
excavados en un material que resiste a la acción erosiva mientras la velocidad o el esfuerzo
tangencial de fricción ejercido por el agua sobre los granos no rebasen a una magnitud,
prefijada de acuerdo con las características del material.

Lógicamente, esta diferenciación cambia el criterio de diseño. En el proyecto de un

canal revestido se calculan las dimensiones óptimas de la sección que proporcionen máxima
eficiencia hidráulica, mínimo costo o ambas. En cambio, en el diseño no revestido rigen los
criterios de velocidad permisible o de esfuerzo tangencial crítico los cuales dependen del tipo
de material en que se excava la sección del canal y que determinan también la rugosidad, la
velocidad mínima permisible para evitar el depósito (si el agua transporta sedimentos), taludes
de la sección, pendiente longitudinal, el borde libre y la sección óptima. En ambos casos, la
tarea del proyectista será minimizar el costo del canal.

Se describirán a continuación los métodos de cálculo indicados, tanto para los canales
revestidos como para los no revestidos.
 31

3.2.- Canales Revestidos

El revestimiento de un canal tiene por objeto prevenir la erosión, evitar las
infiltraciones y disminuir la rugosidad de las paredes.

El volumen de excavación y la superficie de revestimiento son los factores más

importantes en el costo del canal. El primero depende del área de la sección y la segunda del
perímetro mojado. La optimización de estos dos factores reducirá el costo al "mínimo".

La sección de máxima eficiencia hidráulica será la de mínimo perímetro mojado para

un área dada ya que en ella se tendrá la mínima resistencia al escurrimiento, así como el
mínimo costo de revestimiento (o en su defecto, la mínima superficie de infiltración), aunque
no necesariamente la mínima excavación.

Para un gasto dado, la sección hidráulica "óptima" sería aquella para la cual el área es
mínima; esto implica que la velocidad sea máxima. Según las fórmulas de Chezy y Manning,
esto significaría que el radio Hidráulico R=A/P fuera el máximo. Para ello será necesario
minimizar también el perímetro mojado. Los factores que se consideran en el diseño son:

3.2.1.- Taludes de la Sección Transversal

La sección trapecial es la más usada en canales. Normalmente el talud de la sección
no se elige arbitrariamente ya que está supeditado al que pueda resistir el material de
excavación y en la tabla 3.1 (Chow, 1994) se dan algunas recomendaciones al respecto.

Tabla 3.1 Pendientes laterales apropiadas para canales construidos

en diferentes clases de materiales

Material Pendiente Lateral

Roca Aproximadamente vertical
Estiércol y suelos de turba ¼:1
Arcilla rígida o tierra con recubrimiento de concreto ½: 1 a 1:1
Tierra con recubrimiento de piedras o tierra en canales grandes 1:1
Arcilla firme o tierra en canales pequeños 1½:1
Tierra arenosa suelta 2:1
Marga arenosa o arcilla porosa 3:1
 32

3.2.2.- Coeficiente de Rugosidad

Según las Normas e Instructivos para el proyecto de alcantarillados (INOS, 1975) los
valores del coeficiente de rugosidad "n" según el material en tramos rectos se pueden
obtener en tabla 3.2

Tabla 3.2 Valores de n en colectores y canales

Material "n"

a) Colectores Cerrados Prefabricados

PVC y Fiberglass 0,011
Asbesto - Cemento, Acero, HF, HF Dúctil 0,012
Arcilla Vitrificada, Concreto Liso 0,013
Concreto Rugoso 0,015

b) Colectores Cerrados Vaciados en Sitio

Concreto 0,014

c) Canales
Revestidos de Asfalto o Concreto 0,015
Excavados en Tierra 0,022 – 0,03
Lechos Pedregosos y Taludes con Grama 0,035

Si en el canal existen irregularidades superficiales, variación en la forma y tamaño de

la sección, obstrucciones vegetación y meandros se podría utilizar el método de Cowan y los
valores obtenidos en la tabla 3.2 representan los valores de n0 en la fórmula de Cowan:

n = (n0 + n1+ n2 + n3 + n4) m5 (3.1)

Donde:
n0 = valor básico del coeficiente establecido según el material para un canal recto, liso, de
sección uniforme.
n1 = corrección adicional para tomar en cuenta los efectos de irregularidad superficial.
n2 = corrección adicional para tomar en cuenta las variaciones en forma y tamaño de la
sección.
n3 = corrección adicional para tomar en cuenta los efectos de las obstrucciones producidas por
arrastrar, raíces, troncos flotantes.
 33

n4 = corrección adicional para tomar en cuenta los efectos de la vegetación en el cauce.

m5 = factor de corrección para tomar en cuenta el efecto de los meandros.
Los valores de n0, n1, n2, n3, n4, y m5 aparecen en la tabla 3.3

Tabla 3.3 Método de Cowan

Condiciones del Canal Valores
Tierra 0,02
Roca Excavada 0,025
Material n0
Grava Fina 0,024
Grava Gruesa 0,028
Ninguno, como una superficie lisa 0,000
Grado de Escaso, como en los canales bien dragados o aquellos con 0,005
Irregularidad los lados algo socavados n1
Moderado, como en los canales con taludes socavados 0,01
Superficial
Notable, como en los canales muy socavados, o aquellos 0,020
excavados en roca, sin ningún acabado
Gradual 0,000
Variación en forma y
Ocasional n2 0,005
tamaño de la Sección
Frecuente 0,010-0,015
Despreciable 0,000
Escaso, como cuando los arrastres obstruyen algo el paso
0,010-0,015
Efecto de del agua.
n3
Obstrucciones Apreciable, como cuando se ven raíces y árboles que
0,020-0,030
obstruyen el paso.
Notables, como cuando hay grandes troncos atascados, 0,040-0,060
carameras grandes, etc.
Baja: hierbas flexibles, la profundidad del agua es tres veces
0,005-0,010
la altura de la vegetación predominante
Media: hierbas, arbustos; la profundidad del agua es como 0,010-0,025
Vegetación dos veces la altura de la vegetación predominante. n4
Alta: hierbas, pequeños árboles; la profundidad del agua es 0,025-0,050
comparable a la altura de la vegetación predominante.
Muy Alta: la profundidad del agua es como la mitad de la
0,050-0,100
altura de la vegetación predominante.
Escasos: Lm/Ls entre 1,0 y 1,2 1,00
Efectos de
Apreciable: Lm/Ls entre 1,2 y 1,5 m5 1,15
los meandros (*)
Notable: Lm/Ls mayor de 1,5 1,3
(*) Lm es la longitud del tramo medida a lo largo del curso del agua .Ls es la longitud del tramo medida en línea recta.
Nota: Para el cálculo de Lm=Dc y Ls=Cl en Curvas Circulares Simétricas ( Ver Curvatura en Canales (3.2.7))
 34

Si el canal en su totalidad no presenta irregularidad superficial, variación en la forma ni

tamaño de la sección, obstrucciones ni meandros, en la tabla 3.4 se obtiene el valor total de n.

Tabla 3.4 Coeficiente "n" de Manning

Canales Abiertos (Machado, 1979)
I. Canales recubiertos
I.1 Concreto, según los diferentes acabados
I.1.1 Acabado liso 0,013
I.1.2 Acabado con cepillo 0,015
I.1.3 Acabado con cepillo con algo de grava en el fondo 0,017
I.1.4 Acabado rústico 0,017
I.1.5 "Gunite", sección bien acabada 0,019
1.1.6 "Gunite", sección ondulada 0,022
1.2 Fondo de grava y lados según se especifican:
1.2.1 De concreto 0,020
1.2.2 Zampeado 0,023
1.2.3 Gaviones 0,033
1.3 Asfalto
1.3.1 Liso 0,013
1.3.2 Rugoso 0,016
1.4 De hierba, sin podar
1.4.1 Grama bermuda o similar 0,07
1.4.2 Hierbas 0,30 m altura 0,11
1.4.3 Hierbas muy altas 0,14
II. Canales excavados
II.1 En tierra, sección uniforme
II.1.1 Limpios, sin vegetación 0,022
II.1.2 Hierba corta, algunas malezas 0,027
II.2 En tierra, sección aproximadamente uniforme
II.2.1 Sin vegetación 0,025
II.2.2 Hierba corta, algunas malezas 0,030
II.2.3 Plantas acuáticas, enea, etc. 0,035
II.2.4 Fondos de cantos rodados y taludes sin vegetación 0,040
II.3 Dragados
II.3.1 Sin vegetación 0,028
II.3.2 Vegetación ligera en los taludes 0,05
II.4 Sin mantenimiento
II.4.1 Gran cantidad de malezas 0,080
II.4.2 Gran cantidad de malezas, arbustos, y vegetación en las márgenes 0,100
II.5 En roca
II.5.1 Sección uniforme y lisa 0,035
II.5.2 Sección irregular y escabrosa 0,040
III. Canales naturales
III.1 En terrenos con poca pendiente (menor de 1%)
III.1.1 Limpios, rectos, sin pozos, uniforme 0,030
III.1.2 Sección uniforme, rectos, sin pozos, con algo de piedra y malezas 0,035
III.1.3 Tramos con densa vegetación, árboles, arbustos 0,100
III.2 Torrentes, terrenos montañosos, sin vegetación, taludes escarpados
III.2.1 Fondo de grava, escasos cantos rodados 0,040
III.2.2 Fondo con abundantes cantos rodados 0,050
IV. Conductos
IV.1 Tubos de concreto
IV.1.1 Prefabricados 0,012
IV.1.2 Vaciados en sitio 0,014
IV.2 Tubos de metal corrugado galvanizado
IV.2.1 Corriente 0,024
IV.2.2 Con 25% de su periferia pavimentada 0,021
IV.2.3 Con 50% de su periferia pavimentada 0,018
IV.3 Hierro fundido 0,013
IV.4 Acero 0,010
 35

El método de Cowan puede ser empleado tanto para canales revestidos como no
revestidos con la claridad de que los valores de n0 (tipo de material) son diferentes motivado
al revestimiento.

3.2.3.- Velocidades máximas y mínimas

Las normas del INOS (INOS, 1975) recomiendan una velocidad mínima para
acueductos y cloacas de 0,6 m/s, en cambio para drenaje recomiendan 0,75 m/s. El segundo
punto consiste en comprobar que la velocidad circulante sea menor que la velocidad máxima.
A continuación mostramos según el tipo de revestimiento sus velocidades máximas
recomendadas (tabla 3.5):

Tabla 3.5 Velocidades máximas recomendadas en canales revestidos

Material Vmax (m/s)
Concreto Rcc28 = 210 Kgf/cm2 5

Concreto Rcc28 = 280 Kgf/cm2 6

Concreto Rcc28 = 350 Kgf/cm2 7,50

Concreto Rcc28 = 420 Kgf/cm2 9,5

Arcilla Vitrificada 6

Asbesto - Cemento 4,5

PVC (Cloruro de Polivinilo) 4,5

Hierro Fundido Sin límite

Para canales no revestidos la velocidad mínima también se puede considerar 0,75 m/s
y las máximas recomendadas en encuentran en la tabla 3.6 (MOP, 1967) (Aguirre, 1980).
 36

Tabla 3.6 Velocidades máximas recomendadas en canales no revestidos

Tipo de Suelo Vmax (m/s)

Arena fina - no coloidal 0,75
Greda arenosa - no coloidal 0,75
Greda limosa - no coloidal 0,90
Greda firme 1
Grava fina 1,2
Arcilla dura - muy coloidal 1,4
Limos aluvionales - coloidales 1,4
Limos aluvionales - no coloidales 0,9
Materiales Gradados - no coloidales
Greda a Grava 1,4
Limo a Grava 1,6
Esquisto Arcilloso 1,8
Grava 1,8
Grava gruesa 2
Grava a Cantos Rodados 2,3

Cuando el fluido transporte material en suspensión tanto las velocidades máximas y

mínimas en canales revestidos y no revestidos deberían ser estudiadas con mayor profundidad
(Aguirre, 1980 y Chow, 1994).

3.2.4.- Juntas de Dilatación

Las juntas de dilatación en losas de concreto armado en canales se colocan con listones
de madera separados una distancia entre veinte y cien del espesor de la losa (20ed100e),
una vez que fragua el concreto se retira el listón y el espacio vacío se rellena con asfalto o con
el mismo material, no es normativo.

3.2.5.- Confluencia de Corriente en Canales

En la confluencia tanto de canales revestidos como no revestidos el ángulo interior
debe ser  20º, si en lo posible la corriente menor debe descargar a un nivel superior al N.A.N.
de la corriente mayor. En caso de no ser posible el desnivel se deben buscar ángulos < a 20º,
realizar protecciones debidas a los taludes y un estudio de la situación.
 37

3.2.6.- Borde Libre

Según el U.S. Bureau of Reclamation y tomado de Bolinaga (1979) el borde libre tanto
para canales revestidos como no revestidos se obtiene de la figura 3.1 y representa la altura
desde el N.A.N de diseño hasta el tope de bancada. El topo de recubrimiento de roca o concreto
representa la altura desde el N.A.N. de diseño hasta donde mínimo debe de llegar en altura el
material de revestimiento y el tope de recubrimiento de tierra hasta el equivalente hasta donde
debe estar bien compactada la superficie del canal no revestido.

Gasto(M3/s)
Figura 3.1 Borde Libre Recomendable

3.2.7.- Curvatura en Canales

El alineamiento horizontal en canales tanto revestidos como no revestidos se traza
siguiendo los mismos principios del alineamiento de un eje de una vialidad; de hecho, gran
cantidad de canales de drenajes mantienen un alineamiento paralelo al de la vialidad tributaria.
Los elementos geométricos de las curvas (generalmente circulares) se calculan siguiendo las
expresiones conocidas:
 0 
TE  TS  R0Tang   (3.2)
 2 
  0  
E  R0   Sec   1
 (3.3)
  2  
 0 
CL  LS  2R0 Sen  (3.4)
 2 
R0
DC  LM  R0  rad
(3.5)
180
 38

Cuando las aguas siguen una trayectoria curva, la fuerza centrífuga origina una
elevación en el contorno exterior y una depresión en el interior generando un flujo secundario
y una pérdida local de energía. La sobreelevación se puede estimar como:

CV 2 B
y   Borde libre (3.6)
gR0

Donde:
C= Según Bolinaga (1979) tabla 3.7
V= Velocidad media de la sección (m/s)
B= Ancho promedio de la sección transversal del agua (m)
R0 = Radio del tramo de Curva Circular Simétrica

Tabla 3.7 - Coeficientes de Sobre-Elevación en Curvas

Sección
Tipo de Régimen Tipo de Curva C
Transversal
Subcrítico Rectangular Circular 1,0
Subcrítico Trapecial Circular 1,0
Supercrítico Rectangular Circular 2,0
Supercrítico Trapecial Circular 2,0
Supercrítico Rectangular Transiciones Espirales 1,0
Supercrítico Trapecial Transiciones Espirales 2,0
Supercrítico Rectangular Espirales e inclinación de fondo 1,0

Experimentalmente existe una relación Q vs Radio de Curvatura que pueden ser

útiles en el diseño práctico:

Q (m3/s) 20 15 5 1 0,5
Rmin (m) 100 80 20 10 5

3.2.8.- Espesores de revestimiento

Según Torres (1981) los espesores en canales revestidos se obtienen en la fig. 3.2 y
en canales no revestidos se obtienen en la figura 3.3. En los canales revestidos en concreto
con el área de acero mínimo es suficiente motivado a que generalmente las fuerzas actuantes
son pequeñas.
 39

Fig. 3.2 Espesores de Revestimientos de Superficie Dura para Canales

t t1 t2 Notas:
metros m m Si el material de revestimiento
0,60 ó menor 0,3 0,9 requiere una capa protectora de grava
o roca para prevenir socavación o
t 0,60 a 1,20 0,45 1,2 erosión, la excavación se debe
t 1,20 a 1,80 0,6 1,8 ampliar para acomodar el espesor
designado del revestimiento mas el
mayor de 1,80 0,6 2,4 espesor de la capa de grava o roca

Fig 3.3 Sección Típica para Canales con Revestimiento de Tierra

 40

3.2.9.- Procedimiento para el diseño de canales revestidos

a) Se obtienen Q (estudio hidrológico), n (método de Cowan), S0 (topografía).
b) Según la topografía se decide la forma geométrica de la sección transversal. Si la sección
es trapecial o triangular “m” se obtiene por la tabla 3.1 ó figura 3.6. (θ<Φ).
c) Se sustituyen la ecuación óptima en la ecuación de Manning y encontramos "y" y “b”
d) Se chequea que Vmin  Vdiseño Vmax
e) Se obtiene el Borde Libre (BL), el tope de recubrimiento de roca o concreto (TR r ó c) y
llevar el tope de recubrimiento de tierra (TRt) hasta TRr ó c (ver fig 3.1).
f) Se chequean las curvas (y  Borde Libre) si las hay.
g) Se calculan los espesores de revestimiento en fig 3.2 (revestidos) o fig. 3.3 (revestidos en
tierra)
h) Se harían ajustes si cumple la velocidad y y en curvas, dependiente de condiciones
operativas durante la construcción.

Ejemplo 3.1: ¿Calcular las dimensiones de la sección del gasto máximo (sección óptima) en
un canal que debido a la topografía y a la clase de suelo se le fijó un A=9m2 y un talud m=1,5?

𝐴
𝐴 = 2𝑦 2 √1 + 𝑚2 − 𝑚𝑦 2 = 𝑦 2 (2√1 + 𝑚2 − 𝑚) ⇒ 𝑦 = √ = 2,07 𝑚
(2√1 + 𝑚2 − 𝑚)

b  2 y 1  m 2  2my  2 x 2,07 x 1  1,5 2  2 x1,5 x 2,07  1,25m

Ejemplo 3.2 Un canal de fuerza de sección rectangular debe conducir un gasto Q=3 m3/s con
una velocidad V = 1,2 m/s. ¿Calcular las dimensiones de la sección óptima y la pendiente
necesaria si se reviste de concreto (n = 0,017)?

Q 3 2,5
A   2,5m 2  2 y 2  y   1,12m b  2 y  2,24m
V 1,2 2

2
 
 
 1,2 x0,017 
2
1  Vn 
Q AR 2 / 3 ( S 0 )1 / 2  S0   2 / 3    2/3 
 0,0009  0,9 0/00
n R    1,12  
 2 
  
 41

Ejemplo 3.3 Diseñar un canal trapecial interurbano revestido en concreto RCC28 = 280
Kgf/cm2, el material del suelo es una arcilla rígida la cual será excavada. Hay una curva
circular simétrica con R=100 m, ∆° = 40º33'30", S0 = 3 %o, Q=20 m3/s, el efecto de las
obstrucciones es apreciable y sin mantenimiento (vegetación baja).

a) n=(n0+ n1 + n2+ n3+ n4) ms= (0,015 + 0,025 + 0.0075) x 1= 0,0475

  40,5583 
DC  Lm  (100)( 40,5583)  70,78m ; CL  2 x100 xsen   69,32m
180  2 
Lm
  1,02  m5  1
Ls
b) (Según la tabla 3.1 coloco máxima eficiencia hidráulica).

 
2/3
1  y 2 3y
c) 20  3y 2   (0,003)1 / 2  y  2,82m ; b   3,26m
0,0475 2 3

d) V  20
3 ( 2 ,8 2 ) 2
 1 , 4 5 m /s 0,75 m/s < 1,45 m/s < 6 m/s.

e) BL = 1 m. TRc= TRt= 0.4 m.

CV 2 B 1x1,45 2 x 4,88
f) y    0,0105 m  B.L.
gR0 9,81x100

20 2 (3,26 y c  0,577 y c2 ) 3
3
Q 2 Ac
    y c  1,43m  y n  R. subcritico  C  1
g Tc 9,81 3,26  2 x0,577 y c

T  b  2 m y = 3 , 2 6 + 2 x 0 ,5 7 7 x 2 ,8 2 = 6 , 5 1 m
3,26  6,51
B  4,88 m.
2
g) emin=10,5 cm. en concreto reforzado.

0.4 m 1.00 m

m
e =10,5 Cm.
3 3

3 3

m
 42

Ejemplo 3.4 En un diseño óptimo de un canal trapecial y revestido se tiene:

b = 2,5 m ; n = 0,05 ; S0 = 1 % ¿Seleccionar el material más económico de revestimiento?

b  2y 1  m2  2my =>y  b
2 1m 2 2m
 2, 74 m.

2/3
1  2,74 
V   (0,01)1 / 2  2,46m / s
0,05  2 

Según la tabla 3.5 tomo asbesto-cemento ó PVC Vmax=4,5 m/s

3.3 Canales no revestidos

La fórmula de flujo uniforme utilizada en el diseño de canales revestidos es
insuficiente en el caso de los no revestidos debido a que el diseño es esencialmente un
problema de estabilidad de la sección. Si el canal transporta sedimentos o está excavado en
material erosionable, es necesario que no ocurra depósito ni erosión, esto es, el canal debe
estar en equilibrio con respecto al transporte de sedimentos, de manera que la cantidad total
de los mismos que transporta sea a lo largo del canal, o bien, impedir dicho transporte.

Existen fundamentalmente dos tipos de problemas en el diseño de canales erosionables

de acuerdo con las condiciones que deben cumplir y para su estabilidad.

a) Canales transportando agua limpia o material fino en suspensión.

b) Canales transportando material sólido de arrastre sobre el fondo.

Aquí sólo se tratará del diseño de canales erosionables que presentan las características
del primer tipo. Los del segundo tipo son objeto de estudio detallado en la hidráulica fluvial.

En el caso de canales transportando agua limpia o con material fino en suspensión las
condiciones del canal exigen que no se deposite dicho material y que la capacidad erosiva del
flujo sea tal que no erosione el lecho y paredes del canal.

Pueden mencionarse dos métodos que sirven de guía para el diseño de canales en
estas condiciones:
a) Método de la velocidad máxima permisible.
b) Método de la "fuerza tractiva".
 43

3.3.1. Método de la Velocidad Máxima Permisible (M.V.M.P.)

La velocidad máxima permisible o la velocidad no erosionante es la mayor velocidad
que no causará sería erosión en los contornos del canal. Esta velocidad es muy imprecisa y
variable y sólo se puede estimar con experiencia previa y buen juicio. En general los canales
viejos soportarán mayores velocidades que los nuevos, ya que en aquellos el fondo está mejor
estabilizado debido a la deposición de sustancias coloidales. Un canal profundo podrá
conducir el agua a mayores velocidades, sin erosión que otro menos profundo, ya que las
velocidades del fondo son mayores en este que en aquel cuando la velocidad media es la
misma.

Maza (1967) recopiló y actualizó experiencias de otros investigadores y presentó las

tablas 3.8 y 3.9 de velocidades máximas no erosivas para suelos no cohesivos y cohesivos
respectivamente.

Tabla 3.8
Velocidades Medias de la Corriente del Agua, que son Admisibles (No Erosivas) para
Suelos no Cohesivos, en m/s
Suelos no Diámetro
Medio Profundidad Media de la Corriente en m.
Cohesivos de las
Partículas 0,40 1,00 2,00 3,00 5,00 Más de 10
en mm.
Polvo y limo 0,005-0,05 0,15-0,2 0,2-0,30 0,25-0,40 0,30-0,45 0,40-0,55 0,45-0,65
Arena fina 0,05-0,25 0,20-0,35 0,3-0,45 0,40-0,55 0,45-0,60 0,55-0,70 0,65-0,80
Arena media 0,25-1,0 0,35-0,50 0,45-0,60 0,55-0,70 0,60-0,75 0,70-0,85 0,80-0,45
Arena gruesa 1,0-2,5 0,50-0,65 0,60-0,75 0,70-0,80 0,75-0,90 0,80-1,00 0,95-1,20
Gravilla fina 2,5-5,0 0,65-0,80 0,75-0,85 0,80-1,00 0,90-1,10 2,00-1,20 1,20-1,50
Gravilla media 5,0-10 0,80-0,90 0,85-1,05 1,00-1,15 1,10-1,30 1,20-1,45 1,50-1,75
Gravilla gruesa 10-15 0,90-1,10 1,05-1,20 1,15-1,35 1,30-1,50 1,45-1,65 1,75-2,00
Grava fina 15-25 1,10-1,25 1,20-1,45 1,35-1,65 1,50-1,85 1,65-2,00 2,00-2,30
Grava media 25-40 1,25-1,50 1,45-1,85 1,65-2,10 1,85-2,30 2,00-2,45 2,30-2,70
Grava gruesa 40-75 1,50-2,00 1,85-2,40 2,10-2,75 2,30-3,10 2,45-3,30 2,70-3,60
Guijarro fino 75-100 2,00-2,45 2,40-2,80 2,75-3,20 3,10-3,50 3,30-3,80 5,60-4,20
Guijarro medio 100-150 2,45-3,00 2,80-3,35 3,20-3,75 3,50-4,10 3,80-4,40 4,20-4,50
Guijarro grueso 150-200 3,00-3,50 3,35-3,80 3,75-4,30 4,10-4,65 4,40-5,00 4,50-5,40
Canto rodado fino 200-300 3,50-3,85 3,80-4,35 4,30-4,70 4,65-4,90 5,00-5,50 5,40-5,90
Canto rodado medio 300-400 4,35-4,75 4,70-4,95 4,90-5,30 5,50-5,60 5,90-6,00
Canto rodado grueso 400-500 4,95-5,35 5,30-5,50 5,60-6,00 6,00-6,20
o más
 44

Tabla 3.9 Velocidades medias permisibles (no erosivas) para suelos cohesivos,
en m/seg
Denominación de Porcentaje del Suelos poco Suelos medianamente Suelos compactos, Suelos muy compactos,
los suelos contenido de compactos, peso compactados, peso peso volumétrico del peso volumétrico del
partículas volumétrico del volumétrico del material seco, 1,66 a material seco, 2,04 a
material seco hasta material seco, 1,20 a 2,04 ton/m3 2,14 ton/m3
1,20 ton/m3 1,66 ton/m3
Tirantes medios, en m
<0,005 0,005-0,05
0,4 1,0 2,0 3,0 0,4 1,0 2,0 3,0 0,4 1,0 2,0 3,0 0,4 1,0 2,0 3,0
(mm) (mm)
Arcillas y tierras 30-50 70-50
0,35 0,4 0,45 0,5 0,7 0,85 0,95 1,1 1 1,2 1,4 1,5 1,4 1,7 1,9 2,1
muy arcillosas 20-30 80-70
Tierras ligeramente
10-20 90-80 0,35 0,4 0,45 0,5 0,65 0,8 0,9 1,0 0,95 1,2 1,4 1,5 1,4 1,7 1,9 2,1
arcillosas
Suelos de aluvión
0,6 0,7 0,8 0,85 0,8 1,0 1,2 1,3 1,1 1,3 1,5 1,7
y arcillas margosas
Tierras arenosas 5-10 20-40 Según la fig. 3.4 de acuerdo con el tamaño de las fracciones arenosas

3.3.1.1.- Procedimiento para diseñar por el M.V.M.P.

a) Se obtienen Q (estudio hidrológico), n (método de Cowan; donde n0 se podría obtener
por la fórmula n0= (d75)1/6/26 en caso de tener la granulometría) y S0 (topografía). d75
se expresa en metros (diámetro al cual se acumula el 75% del material en peso desde
los diámetros menores a los mayores).
b) Según la topografía se decide la forma geométrica de la sección transversal. Si la
sección es trapecial o triangular “m” se obtiene de la tabla 3.1 ó fig 3.6. (θ<Φ).
c) Con la granulometría (d50) y tipo de suelo se asume una profundidad y una velocidad
máxima en la tabla 3.8 ó fig. 3.4 y tabla 3.9 dependiendo del tipo de material; en su
defecto se obtiene Vmax en la tabla 3.6.

Empleando las expresiones de A y P en función del ancho b de la base y la profundidad

y (figura 2.5), resolver para "b" y "y". Si “y” es diferente al yasumido , con el valor de “y”
en la tabla 3.8 ó fig. 3.4 y tabla 3.9 se obtiene una nueva Vmax y se va al paso c) hasta
que aproximadamente yultimo sea igual a ypenúltimo. Si la Vmax se obtuvo de la tabla 3.6 (no
hay granulometría) simplemente la primera solución real del sistema de ecuaciones
(Área y Perímetro) representa la solución “b” y “y”.
d) Continuar igual a canales revestidos
 45

Ejemplo 3.5 Construir la tabla 3.8 en forma gráfica para crear la figura 3.4

Figura 3.4 Velocidad permisible en materiales no cohesivos

Ejemplo 3.6 ¿Determinar el tipo de material que resista la velocidad V = 1,2 m/s y y=1,12 m?

Para y = 1,12 m y material cohesivo, de la tabla 3.9 se admitiría, por ejemplo, un suelo
arcilloso de 1,66 a 2,04 tf/m3 de peso volumétrico seco (los 2 primeros). Para material no
cohesivo, de la figura 3.4 esta debería ser grava gruesa con d50 = 15 mm.

Ejemplo 3.7 (El mismo 3.3 pero no revestido). Diseñar un canal trapecial interurbano no
revestido, el material del suelo es una arcilla rígida la cual será excavada. Hay una curva
circular simétrica con R = 100 m y º = 40º33'30", S0 = 3 %o, Q=20 m3/s, el efecto de las
obstrucciones es apreciable y sin mantenimiento (vegetación baja). n0= 0,026
 46

a) n = (n0+n1+n2+n3+n4)m5=(0,026+0,025+0,0075) x 1= 0,0585

b) (tabla 3.1)  1
2  m  1  OK

c) Vmax = 1,4 m/s (tabla 3.6)

3/ 2
V n 
3/ 2
 1,4 x0,0585 
R permisible   max1 / 2 
   
1/ 2 
 1,828m
 S 0    (0,003) 

20 A 14,28
A  14,28m 2 Pmin    7,812m
1,4 R permisible 1,828

14,28  by  0,577 y 2
 1,732 y 2  7,812 y  14,28  0  Ambas Raices Imaginaria s
7,812  b  2 y 1  0,577 2

Bajamos la velocidad máxima en un 25%  Vmax = 1,4*0,75 = 1,05 m/s

3/ 2
 1,05x0,0585 
R permisible   
1/ 2 
 1,19m
 (0,003) 

A=20/1,05=19,05m2 Pmin= 19,05/1,19=16,01m

19,05  by  0,577 y 2  y  7,84m  b  2,09m

 1,732 y 2  16,01y  19,05  0   1
16,01  b  2 y 1  0,577 2  y 2  1,40m  b  12,80m
T=b+2my= 14,42m; B=(T+b)/2= 13,61m

d) 0,75 m/s < V = 1,05 m/s = Vmax

e) En la fig 3.1  BL= 1 m; TRt = 0,26 m

CV 2 B 1x1,05 2 x13,61
f) y    0,015m  BL  1m .
gR 9,81x100
Q 2 Ac 3 20 2 (12,80 y c  0,577 y c2 ) 3
    y c  0,62m y n  Rsubcritico  C  1 (Tabla3.7)
g Tc 9,81 (12,80  2 x0,577 y c )
g) De la figura 3.3 tenemos t1=0,6m y t2=1,8m
 47

Ejemplo 3.8 Diseñar un canal trapecial estable por el método de la velocidad máxima
permisible que conduzca un Q = 4 m3/s. El suelo en que se realizaría la excavación tiene un
d50 = d75 = 5 mm y un S0 = 0,001. Las partículas del material del fondo son muy angulares
(tierra arenosa suelta).

a) n = (d75)1/6/26 = 0,0051/6/26 = 0,016

b) De la tabla 3.1 (m = 2)

c) De la fig 3.4 con d50 y asumiendo y = 1 m, tenemos Vmax = 0,85 m/s; entonces

Rp (m)
Fig 3.4 V A (m2) P (m) y (m) b (m)
y (m) MAX n 3 / 2
Vmax ( 1/ 2
) Q/VMAX A/Rp (m  2 1  m ) y  py  A  0
2 2
P  2 y 1  m2
(so )
1 0,85 0,282 4,71 16,70 0,295 15,38
0,295 0,75 0,23 5,33 23,17 0,24 22,10
0,24 0,74 0,229 5,40 23,58 0,234 22,53

d) La velocidad de Diseño  0,74 m/s no cumple con la velocidad mínima (Vmin= 0,75m/s)

e) De la figura 3.1 BL= 0,70m; y TRt = 0,15 m

3.3.2.- Método de la fuerza tractiva o esfuerzo tangencial crítico

El esfuerzo cortante promedio que produce el flujo de agua sobre el perímetro mojado
de un canal está dado por le Ec. (2.1) 0 = RS0; donde:  (Peso específico del fluido [Kgf/m3]),
R (radio hidráulico [m]) y S0 es la pendiente del fondo del canal.

Por otra parte Shields en 1936 determinó que el esfuerzo cortante que produce el inicio
del movimiento en el material de fondo para material suelto, o no cohesivo está dado por:

     c = 0,056(s-)d75 (3.7)

Donde: c= Esfuerzo crítico (Kgf/m2)

s= peso específico del suelo saturado (Kgf/m3)
 = peso específico del agua (Kgf/m3)
 48

d75 = Corresponde al diámetro donde se acumula el 75% del material en peso cuando
este se va acumulando desde los diámetros pequeños a los mayores (m) [Rango d 75
 1 mm]
(s-) = peso específico del suelo sumergido.

La fig. 3.5 presenta los esfuerzos cortantes críticos en suelos cohesivos y la fig. 3.6
ayuda a la tabla 3.1 a la hora de obtener la inclinación del talud

Si conocemos la granulometría n0 se puede obtener por la expresión:

(3.8)

Esta información permite calcular, sin considerar los esfuerzos en los taludes, un
canal de fondo estable siguiendo los siguientes pasos:

a) Igual al método de la velocidad máxima permisible.

b) Igual al método de la velocidad máxima permisible.
 49

0,056 s   d 75
c) Hallar Rhmax 
S 0

; ;

Empleando las expresiones de A y P en función de "b" y "y" se resuelve el sistema. Si las

raíces son Imaginarias bajamos el Rhmax en 25% y repetimos el proceso hasta que tenga
solución el sistema de ecuaciones.

d) Continuar igual a canales revestidos.

Ejemplo 3.9 Un canal debe transportar 2800 l/s a través de un terreno compuesto de granzón
grueso con d75 = 5 cm. La excavación para la cual tiene una pendiente S0 = 0,005 las orillas
están cubiertas con hierba y protegidas de la erosión. ¿Calcular las dimensiones mínimas para
que el fondo sea estable utilizando el método de la fuerza tractiva? s = 2650 Kgf/m3

a)

1
b) d75 = 5cm en la fig 3.6 (Poco angular) Φ= 37º Tg37º =  m  1,32  1,35  θ<Φ
m

0,056(2650  1000)0,05
c) Rhmax   0,924m
1000 x0,005

m
 50

 
0,9589  by  1,35 y 2   2,01y  1,37 y  0,9589  0
2

 Raices imaginarias
1,037  b  2 y 1  1,352  b  1,037  3,36 y 
 
Rhmax = 0,75*0,924 = 0,70 m

 
1,16  by  1,35 y 2   2,01y  1,66 y  1,16  0
2

 Raices imaginarias
1,66  b  2 y 1  1,35 2  b  1,66  3,36 y 
 
Rhmax = 0,75*0,70 = 0,525 m

V
1
0,5252 3 0,005 12  2m / s
0,023

1,4
A
2,8
2
 1,4 m2 P  2,66m
0,525

 
1,4  by  1,35 y 2   2,01y  2,66 y  1,4  0
2

 Raices imaginarias
2,66  b  2 y 1  1,35 2  b  2,66  3,36 y 
 

Rhmax = 0,75*0,525 = 0,39375 m

 
1,69  by  1,35 y 2   2,01y  4,29 y  1,69  0
2
 y1  1,61m  b1  1,12m
 
4,29  b  2 y 1  1,352  b  4,29  3,36 y  y 2  0,52m  b 2  2,54m
 
d) V = 1,65 m/s > 0,75 m/s
e) De la fig 3.1 TRt = 0,15 m y BL = 0,62 m
 51

3.3.2.1.- El Problema de los taludes

Si el ángulo de reposo de un determinado material es  el criterio de equilibrio para
ese material en un talud sobre el que actúa agua en reposo (sin movimiento), se establece,
con referencia a la fig 3.7 como sigue

Fig. 3.7 Condición de equilibrio de un talud

Si no hay movimiento PSen = NTg  PSen = PCosTg   = , donde P es el

peso del agua por unidad de área, sobre una partícula.

Si hay movimiento, sobre las partículas que descansan en los taludes de un canal
trapecial actúan dos fuerzas por unidad de área: La fuerza tangencial de arrastre 0 y la
componente del peso en la dirección de la pendiente máxima del talud PSen, las cuales hacen
que la partícula tienda a moverse. La resultante de estas fuerzas está equilibrada por la fuerza
de fricción ejercida sobre ella, que es igual al producto de la componente normal al talud
correspondiente al peso de la partícula (PCos)multiplicada por el coeficiente de fricción
interno: Tg. En el caso límite, cuando la partícula está a punto de rodar, se establece el
siguiente equilibrio:
Tg 2
PCosTg  P 2 Sen   02 , Despejando  0 (talud )  PCosTg 1  (3.9)
Tg 2

En el caso de partículas descansando en el fondo del canal  = 0 la ecuación será:

c = PTg, llamando k a la relación entre el esfuerzo tangencial crítico en los taludes 0 y el

esfuerzo tangencial de arrastre en el fondo c, se tiene:
 52

0 Tg 2
K   Cos 1  2 (3.10)
c Tg 

Expresión esta que permite hallar el esfuerzo cortante crítico en el talud (0), dados el ángulo
del mismo y el ángulo de reposo del material que lo compone.

El U.S. Bureau of Reclamation ha realizado extensos estudios para determinar el

ángulo de reposo de los distintos materiales y la distribución de esfuerzos cortantes en los
taludes, ya que no es satisfactorio usar un esfuerzo cortante único igual al promedio en la
sección. Como producto de estas investigaciones en relación con el ángulo de reposo en
material no cohesivo se presenta la fig 3.6.

En relación con la distribución de esfuerzos, es evidente que estos deben variar con el
ángulo  y la relación b/y. Pero para los propósitos prácticos de diseño se puede considerar
que el esfuerzo cortante en el fondo es 0 = RS0 y que el máximo en los taludes está dado por
0 = 0,75RS0. El procedimiento de diseño se puede resumir en los siguientes pasos:

a) Igual al método de la velocidad máxima permisible.

b) Igual al método de la velocidad máxima permisible.

c) Determinar la relación de esfuerzos, 0/c ;

0 0,75RS 0 Tg 2
  Cos 1  2 (3.11)
 c 0,056( s   )d 75 Tg 

Calcular Rh (máximo) ó S (máximo) de la Ecuación 3.11

Calcular “b” y “y”, empleando la Ecuación de Manning, mediante un tanteo hasta que el
Rcal ≤ Rh (máximo).

d) Continuar igual a canales revestidos.

 53

Ejemplo 3.10.- Resolver el ejemplo 3.8 por el método de la fuerza tractiva para que los taludes
sean estables. Datos: Q = 4 m3/s; d50 = d75 = 5 mm; S0 = 0,001; partículas de fondo muy
angulares

a) n = n0

1
b) De la fig 3.6 con D75 = 5 mm,  32º, m   1,6  2  Tg = 1/2 
tg32º
 = Arctg 0,5 = 26,56º   < 

0,056(2650  1000) x0,005 Tg 2 (26,56)

c) Rh max  Cos26,56 1   0,33m
0,75x1000 x0,001 Tg 2 (32)

con la Ecuación de Manning

θ

Dando valores de “b” obtengo “y” y chequeo que Rhmax 0,33 m.

b (m) y (m) A=(by+2y2) Rcal =A/P Rcrítico

10 0,38 4,08 11,69 0,35> 0,33
14 0,31 4,53 15,38 0,29< 0,33

Respuesta b = 14 m, y = 0,31 m ,

4
d) V   0,88m / s  0,75m / s
14 x0,31  2(0,31) 2

e) De la figura 3.1 tenemos BL = 0,7 m, TRt = 0,15 m, la dimensión del revestimiento en

tierra se obtiene de la figura 3.3.

Ejemplo 3.11 Resolver el ejemplo 3.8 por el método de la fuerza tractiva para que el fondo
sea estable.
Del problema anterior tomamos: Q = 4 m3/s; d50 = d75 = 5 mm; S0 = 0,001; n = 0,016; m =2

0,056( s   )d 75 0,056(2650  1000)0,005

entonces Rh max   0,462m
S 0 1000 x0,001
 54

3,38  by  2 y 2  2,47 y 2  7,33 y  3,38  0 y  0,57m

 
7,33  b  2 y 1  2 2  b  7,33  4,47 y  b  4,78m

Las dimensiones finales de los problemas 3.8, 3.10 y 3.11 son las siguientes respectivamente:

y = 0,234 m b = 22,53 m Rp = 0,229 m (M.V.M.P.)

y = 0,31 m b = 14 m Rhmax = 0,29 m (M.F.T. para que el talud sea estable)
y = 0,57 m b = 4,78 m Rhmax = 0,462 m (M.F.T. para que el fondo sea estable)

3.3.3.- Método propuesto por el autor.

El autor propone para diseñar canales no revestidos seguir el procedimiento de canales

revestidos (diseño óptimo) respetando el ángulo de reposo del material ( <  y
posteriormente chequear los siguientes pasos para evitar la erosión:

1) Para la velocidad máxima permisible  Vmin ≤ V ≤ Vmax

2) Para el método de la fuerza tractiva con fondo estable  0 ≤ c

donde, 0 = RhS0
c = 0,056( s   )d 75

3) Para el método de la fuerza tractiva con taludes estable  0 ≤ c

donde, 0 = 0,75RhS0
Tg 2
c = 0,056( s   )d 75 Cos 1  2
Tg 

Ejemplo 3.12(Ej.3.7propuesto por el autor) n=0,0585 ;S0 = 3 %o; Q=20

m3/s
a) Q=20 m3/s; n = 0,0585; S0 = 3 %o

b) (tabla 3.1)
 55

2
 y 3
1
1 2 3. y
c) 20  3. y 2   0,003 2  y  3,05m  b   3,52m
0,0585 2 3
20 m m m m
d) V   1,24  0,75  1,24  1,4 (tabla3.6)
3.3,05 s s s s

e) BL = 1m; TRt = 0,26m

20 2 (3,52 y c  0,577 y c ) 3
2
Q 2 Ac 3
f)     y c  1,37m y n  3,05m  R.Subcrítico  c  1
g Tc 9,81 3,52  (2  0,577 y c )

Como T=b+2my=7,04m y b=3,52m entonces B= (b+T) / 2=5,28m

cV 2 B 1  124
, 2  5, 28
y    0,008 m  BL  1m
gR 9,81 100
g) De la Figura 3.3 t1=0,6m t2=2,4m

1) Para la velocidad máxima permisible, ver chequeo en el método propuesto por el autor

2) Para el método de la fuerza tractiva con fondo estable:

A 3.3,05 2
Rh    0,50 m, Asumo d50 = d75 = 50mm (muy angular)
P 2 3.3,05

0 = *0,5*0,003 = 1,5 Kgf/m2

0 ≤ c
c =  = 4,62 Kgf/m 2

3) Para el método de la fuerza tractiva con taludes estables:

fig .3.6
Con d75 = 50mm  =40°  =90º - 30º = 60º   >  No cumple.
muy..angular Ver nota.

0 =  = 1,13 Kgf/m2 No cumple.

0 > c Ver nota.
Tg 2 60º
c = 4,62* Cos 60 º 1  = Raiz imaginaria
Tg 2 40º
Nota: Se podría recalcular el método con un “m” mayor de tal forma que  < . En estos casos
también se podría recubrir con bentonita los taludes para aumentar su estabilidad.
 56

Ejemplo 3.13. Diseñar un canal no revestido Q=5 m3/s; d50 = d75 = 50mm (redondo);
S0 = 0,001; recto uniforme.
1
6
0,05
n0   0,0233
26

fig.3.6 1  1 
Con d75= 50mm =37º  m   1,327  1,35    arctg   36,53º( < 
redondo tg  1,35 

 
2
 y
1
1 3
5 2 y 2 1  1,35 2  1,35 y 2   0,0012  y=1,49m; b=0,98m; A= 4,46m2, P=5,99m
0,023 2
V=1,12 m/sfig. 3.4 0,75m/s ≤ 1,12m/s ≤ 2,10m/s

Chequeo del método de la fuerza tractiva para que el fondo sea estable:

0 = *(4,46/5,99)*0,001 = 0,74 Kgf/m2

0 ≤ c
c =  = 4,62 Kgf/m 2

Chequeo del método de la fuerza tractiva para que los taludes sean estables:

0 = 0,75**(4,46/5,99)*0,001 = 0,56 Kgf/m2

0 ≤ c
Tg 2 36,53º
c = 4,62* Cos36,53º 1  = 0,68 Kgf/m2
Tg 2 37º
Ejemplo 3.14. Determine la máxima pendiente permisible en un canal trapecial. Aplicar la
fuerza tractiva considerando la ecuación de estabilidad de los taludes. El canal tiene las
siguientes características B = 20 m, y = 2 m, d = 3,5 mm,  36o, m= 2. Encuentre además el
caudal transportado.

A  by  my 2  20 x2  2(2) 2  48m 2

P  b  2 y 1  m 2  20  2 x2 x 1  2 2  28,94m
 57

  arctg (1 / 2)  26,56º    36º

0 0,75RS 0 tg 2 0,056 x1650 x0,0035 tg 2 26,56

  Cos 1  2  S 0  Cos(26,56) 1 
 c 0,056( s   )d 75 tg  0,75 x1000 x1,658 tg 2 36

48
S 0  1,68 x10 4  Q  (1,658) 2 / 3 (0,000168)1 / 2  58m 3 / s
0,0149

3.4. Proyecto de un canal

El diseño de una sección dada simplemente es una parte del proyecto de un canal, un
buen estudio tanto cualitativo como cuantitativo de la situación daría un mejor proyecto, en la
realidad hay diversos factores que influyen en la elaboración de un proyecto por lo cual se
debe realizar un anteproyecto.

Mencionaremos los pasos más importantes (desde el punto de vista técnico) en la

elaboración de un proyecto:

1.- Recopilación de información.

2.- Recorrido e inspección del terreno para detectar cualquier anomalía.
3.- Trazado del eje tentativo en planta buscando la ruta lógica, económica y de máximo
aprovechamiento.
4.- Trazar las curvas horizontales y perfil longitudinal del terreno.
5.- Trazar la rasante del canal.

La rasante y la profundidad depende de la condición del terreno, lo que si se debe

procurar es compensar los volúmenes de corte y relleno por motivos de economía teniendo en
cuenta que es preferible volúmenes sobrantes que tenerlos que comprar; si a los lados del canal
se van a hacer vías paralelas a dicho canal, el material sobrante de esta rasante sirve para
relleno de dicha vía (Ver Fig. 3.8). Tomar progresivas cada 100 mts.
 58

Progresiva
C. Terreno
C. Rasante
Ord. Corte
Ord. Relleno
Alineamiento
Fig. 3.8. Rasante en el Perfil Longitudinal

6.- Dibujar las secciones transversales (papel milimetrado)

Esc. H=V=1:200
7.- Ubicar la sección típica diseñada (óptima)

8.- Calcular las áreas y volúmenes

9.- Determinar el movimiento de tierra y diagrama de masas.

Se determina para obtener un equilibrio y volumen de material necesario o sobrante
según el resultado del diagrama de masas. La distancia de libre acarreo es aquella en que el
material puede ser transportado sin ningún costo adicional. La distancia sobre acarreo se paga,
pero habría que hacer un estudio para ver si es más económico que comprar y transportar de
otro material de préstamo.
 59

3.4.1. Propiedades del diagrama de masas

1.- Cualquier línea horizontal paralela a la línea base me permite identificar un volumen.
2.- Una rama ascendente en la curva de masas significa un corte en el perfil longitudinal.
3.- Una rama descendente en la curva de masas indica un relleno en el perfil longitudinal.
4.- Cuando las secciones transversales son a media ladera una rama ascendente del
diagrama de masa indica un dominio del corte sobre el terraplén, en caso contrario
indica un mayor dominio del terraplén sobre el corte.
5.- En general un punto máximo o mínimo en el diagrama de masas indica en el perfil
longitudinal un paso de corte a relleno o de relleno a corte respectivamente. (Fig. 3.9).

Perfil longitudinal

ota

Diagrama de Masa

Fig. 3.9. Perfil Longitudinal y Diagrama de Masas respectivamente

 60

CAPITULO IV
ORIFICIOS Y VERTEDEROS

4.1. Orificios
En el presente capítulo se resume de una manera sencilla la forma de determinar el
caudal que sale por un orificio. El capítulo abarca todos los tipos de orificios, pero se extiende
en forma más detallada a orificios de pared delgada, en vista de que a partir de ellos se pueden
determinar los gastos de los orificios de pared gruesa, por supuesto teniendo en cuenta las
diferencias existentes entre ambos; y además, teniendo en cuenta que suele ser este tipo de
orificio el más característico por su costo y funcionabilidad tan sencilla.

El tema de orificio, es muy amplio, y se necesitaría todo un texto para lograr abarcarlo
completo. Sin embargo, este capítulo hace un resumen, que a nivel estudiantil será muy útil
ya que servirá de guía para entender mejor la solución de los problemas que suelen
presentarse.

En general, trataremos de ser lo más objetivo posible en la elaboración de este capítulo

y realizarlo de una manera tal, que el estudiante de hidráulica lo utilice como guía a la hora
de calcular el gasto que escurre por un orificio.

4.1.1 Definiciones

4.1.1.1. Orificio: En lo que a hidráulica se refiere, es una abertura de forma regular practicada
en paredes, o en el fondo de un recipiente y por el cual escurre agua en forma continua y en
contacto con toda su sección. Esto implica el hecho de que el contorno del orificio debe estar
totalmente sumergido.

4.1.1.2. Vena Contracta: Es la parte de menor sección de la vena líquida que escurre a través
de un orificio. Su importancia radica en que en esa sección de vena líquida, las líneas de
corriente son paralelas y la sección es mínima. En ese punto se termina de producir la
contracción y por ende la pérdida de energía por la salida del orificio.
 61

4.1.1.3. Tipos de Orificios

• Orificio de Pared Delgada (e/h <0,67): Cuando la vena líquida fluye a través de una arista
viva. Eso puede suceder en dos casos:
a) Cuando el espesor de la pared es pequeño en comparación con la carga de agua
(Fig. 4.1 b).
b) Cuando las paredes del orificio son biseladas, este es el caso en el cual se puede
asegurar que el orificio es de pared delgada (Fig. 4.1 a).

• Orificio de Pared Gruesa (e/h  0,67): Cuando la vena líquida que escurre a través del
orificio se pone en contacto con las paredes gruesas del mismo (Fig. 4.1 b).

Fig. 4.1 Orificio de pared delgada y gruesa

• Orificio con Contracción Parcial: Si el orificio está en contacto con el fondo o una pared
lateral se suprime la curvatura de las líneas de corriente que se encuentran de ese lado (Fig.
4.2 a).
• Orificio con Contracción Incompleta: (D/2<X<3D): Es cuando la contracción de las
líneas de corriente puede hacerse menos pronunciada mediante adecuadas disposiciones (Fig.
4.2 b).

Fig. 4.2. Orificios con contracción parcial e incompleta

 62

• Orificio Totalmente Sumergido: Es cuando el nivel aguas abajo está por encima del borde
superior del orificio (Fig. 4.3a).

• Orificio Parcialmente Sumergido: Cuado el nivel aguas abajo está comprendido entre el
borde superior y el borde inferior del orifiio (Fig. 4.3 b).

Fig. 4.3 Orificio totalmente y parcialmente sumergido

Estos dos últimos tipos de orificios presentan la característica de que la vena líquida
rige la "Ley Hidrostática" con la altura piezométrica correspondiente a la masa líquida que la
rodea. Mientras que en los otros casos cuando la vena desagua a la atmósfera, la presión que
la rodea es prácticamente constante. Tomando en cuenta que las dimensiones de la vena son
casi siempre reducidas, se puede aceptar que la presión que la rodea rige en toda la vena
contracta. En orificios de pared gruesa o en disposiciones equivalentes, la vena se contrae en
la entrada y se ensancha después hasta llenar todo el tubo, ocasionándose así una pérdida de
carga.

Los orificios de pared delgada quedan caracterizados por la contracción; los orificios
de pared gruesa se caracterizan por la pérdida de carga en la entrada, pérdida por fricción y
contracción final, circunstancia que pueden faltar o aminorarse según las disposiciones. En la
tabla (4.1) se presenta un resumen de los tipos de orificios según su espesor, tipo de
contracción y tipo de descarga.
 63

Tabla 4.1 Tipos de Orificios

Característica Tipos Límite de Aplicación
Según su Espesor Pared Delgada e/h < 0,67 ó aristas vivas
Pared Gruesa e/h  0,67
Según la Contracción Parcial En contacto con el fondo o
pared lateral
Incompleta D/2 < X < 3D
Total X  3D
Según la Descarga Totalmente Sumergido
Parcialmente Sumergido
Descarga Libre

4.1.2. Orificios de pared delgada. Contracción total. Descarga libre

Cuando la vena líquida sale del orificio, se forma en una parte de ella un achatamiento,
es decir se presenta una sección de la vena líquida muy contraída, a esta sección se le denomina
"Vena Contracta". Esta vena se caracteriza por lo siguiente: Los filetes exteriores de la vena
líquida que sale de un orificio de pared delgada escurren aguas arriba de la línea de contacto
aplicados a la pared del depósito, después lo hacen al aire libre; estas trayectorias pasan
rápidamente de la dirección tangencial a la pared, a una dirección prácticamente normal a ella,
es decir, en ese cambio brusco de dirección, los filetes presentan una curvatura muy fuerte,
pero no infinita, o sea, que presentan un radio de curvatura finito. Este razonamiento se hace
debido a que las fuerzas actúan sobre las moléculas de estos filetes, no pueden producir una
discontinuidad en su dirección o velocidad. La tensión superficial tiene en la vena contacta
una parte muy secundaria; esto se comprueba suprimiendo perfectamente la contracción en
parte del perímetro de un orificio que desagua al aire con tabiques interiores que impidan
únicamente la curvatura de las trayectorias de las partículas afluentes. O sea, que en la vena
contracta se debe aceptar el cumplimiento de la ley hidrostática, porque es allí donde los filetes
o líneas de corrientes de la vena líquida son totalmente paralelos, produciéndose también la
sección mínima de la vena líquida. Esto se debe cumplir para todos los casos, ya sean en
contracciones parcialmente suprimidas, incompletas, sumergidas o parcialmente sumergidas.

Para hallar la velocidad en la vena contracta, para un orificio que descarga a la

atmósfera, aplicamos Bernouilli entre la sección 1 y 2 de la figura (4.4.); entonces:
 64

h V2
2g =>V  2gh (4.1)

Figura 4.4

En la ecuación (4.1) de Toricelli h es la distancia desde el centro del orificio a la

superficie del depósito. Esta velocidad sería la velocidad teórica o ideal, en ella se consideran
a las pérdidas despreciables. La gravedad influye sobre el fluido al salir por el orificio de la
siguiente manera: a medida que el flujo va saliendo del depósito, la altura h, va variando
haciendo seguir a la vena líquida una variación en su forma. En la figura (4.5) se puede
observar como varía la vena líquida a medida que va saliendo del orificio, notándose la
variación de los filetes a lo largo de la vena.

Figura 4.5
 65

Q= Qr= Cd Qi=Cd (Vi Ao)= Cc Cv (Vi Ao)= (Cv Vi) (Cc Ao)= Vr Ar, entonces:

(4.2)

Q= Qr= Caudal que sale del orificio

Qi= Caudal ideal que sale del orificio
Cd=Cc Cv= Coeficiente de descarga (figura 4.6)
Cc= Coeficiente de contracción (figura 4.6)
Cv= Coeficiente de velocidad (figura 4.6)
Vi= Velocidad ideal
Ao= Ai= Area del orificio
Ar=Ac=Ao Cc= Area real en la vena contracta
Vr= Vi Cv= Velocidad real en la vena contraída.

El coeficiente de contracción es la razón entre la sección de la vena contraída (Ac) y

la sección del orificio (Ao), es decir, Cc= Ac/Ao.

Si la contracción es total este coeficiente vale aproximadamente 0,64 y si es incompleta

aproximadamente 1.

La manera más completa de obtener los coeficientes de contracción, velocidad y

descarga es mediante la figura (4.6); utilizable para la contracción total. Si la forma del orificio
es rectangular en contracción total se podría utilizar su menor dimensión, en lugar de el
diámetro en el número de Reynolds obteniendo resultados bastante ajustados.

Figura 4.6. Variación de los coeficientes de velocidad, contracción y gasto, con el número
de Reynolds en un orificio circular (Sotelo, 1982).
 66

El coeficiente de velocidad expresa la influencia de la fricción y de la viscosidad,

generalmente su valor varía entre 0,96 y 0,99.

La pérdida de energía en el orificio es proporcional a la carga de velocidad media en

la sección contraída. El coeficiente de pérdida K no tiene dimensiones y es función sólo del
coeficiente de velocidad.

La ecuación (4.3) nos permite calcular la pérdida de energía en un orificio de pared

delgada en contracción total y descarga libre.

(4.3)

1
K=( − 1)
Cv 2

En la práctica teniendo un par de coordenadas de la lámina de agua que escurre por un

orificio, midiendo el caudal y teniendo la geometría de un recipiente elevado como muestra
la figura (4.7), podemos obtener el coeficiente de velocidad, velocidad real en la vena
contraída y comparar los coeficientes y caudal obtenidos en el laboratorio con los obtenidos
utilizando la fórmula (4.2) y la figura (4.6).

Fig.4.7
 67

Aplicando las ecuaciones de movimiento de un cuerpo en el espacio en la vena

contraída:
gt 2 gt 2 x x x 2g
X  Vox t  ; y  Voy t   Vox  Vr   (4.4)
2 2 t 2y 2y
g

como Vr=Vi Cv  Cv= (4.5)

Ejemplo 4.1: En la figura (4.7) D= 0,1 m y desaloja un gasto Q= 29,5 l/s con una carga h= 2
m. Si tenemos las coordenada (x,y)=(3m, 1,15m). Calcular los coeficientes y comprobarlos
con la figura (4.6).

; Cv=

0,0295 Cd
Cd   0,6 ; Cd   0,607
 Cv
(0,1) 2 2 gx2
4

2hg xD 2 gx2 x0,1

R  6
 6,26 x10 5 de la figura 4.6
v 10

Cd= 0,6 ; Cc= 0,605 ; Cv= 0,992

Ejemplo 4.2. El recipiente está descargando agua como muestra la figura se pregunta A) Si el
flujo es ideal ¿hallar "h"; el diámetro (D) del flujo a la elevación de 30 cm. después de salir
del orificio y los valores de "x" y "V" al final de la caída del agua? B) Si C d=0,6; Cv= 0,98;
Cc= 0,612. ¿Hallar "h", "x", "v" y "hf" ?.

Recordatorio
at 2 at 2
x  V0 x t  y  V0 y t 
2 2

V f  V0 t  at (Acelerado)

(Retardado)
 68

Respuesta A) V o  2g4  8, 85m/s ; Vox= 7,67 m/s ; Voy= 4,429 m/s

V y  V oy  9, 81t  t  0, 4515 seg.

9,81
h=4,29x0,4515- 2 (0, 4515 )
2
 0, 99m

x1=7,62 (0,4515)=3,46m

2x9,81x0,3=(4,429)2-Vy 2  Vy  3, 7m/s  V  7, 67 2  3, 7 2  8, 51m/s

Q  8, 85x 4 (0, 05 ) 2  0, 0174  8, 51x 4 D 2 => D  0, 051m

9,81t 2
20, 99  2
 t  2, 0686 seg.  x  7, 67(2, 0686 )  15, 86m

V y  9, 81(2, 0686 )  20, 29m/s  V  20, 29 2  7, 67 2  2 1, 69m/s

Respuesta B) Vo=8,85 x 0,98=8,673 m/s ; Vox= 7,51 m/s ; Voy=4,3365 m/s

 69

9,81
V y  V oy  gt  t  0, 442 seg. h  4, 3365(0, 442 )   0, 9583m
 2
2 (0, 442 )

9,81t 2
20, 95  2
 t  2, 067 seg.  x  7, 51(2, 067 )  15, 52m

V y  9, 81x2, 067  20, 27m/s  V  7, 51 2  20, 27 2  21, 62m / s

8,673 2
4 2g  hf  h f  0, 166m

Ejemplo 4.3. En el ejemplo anterior, calcular la altura máxima que alcanza un chorro que sale
de un plano inclinado en función de H y . Despreciar las pérdidas

gt 2 gt 2
x  V ox t  2  2gh sen  t 1 ; y= 2gh cos  t- 2 2 t de 1 en 2
;
2
y   4H Sx en 2   x cot g dy
 0   2H Sxen2  cot g  x  2 H Sen Cos
dx 

Sustituyendo x en y  y  H Cos 2 

Otra forma será aplicando Bernouilli Hy V 

2g V  2g(H  y )

como la velocidad en el eje x es la misma 2gH Sen = 2g(H  y )  y  HCos 2 

4.1.3. Orificio de pared delgada. Contracción incompleta. Descarga libre

Cuando las paredes o el fondo del recipiente se encuentran a distancias (medida al

centro del orificio) inferiores a 3D (D es el diámetro de los orificios) o bien, a 3a (a, dimensión
mínima en orificios rectangulares), se dice que la contracción es incompleta, entonces:

(4.6)

donde Cd es el coeficiente de gasto del orificio; Cdo el coeficiente de gasto del mismo orificio
con contracción total; Ao el área del orificio; At el área de la pared del recipiente en contacto
con el agua (Fig.4.8).
 70

Fig. 4.8 Contracción Incompleta en un Orificio

Ejemplo 4.4. El orificio de pared delgada de la fig 4.9, es cuadrado (a = 0,18 m) y trabaja con
una carga h = 0,5 m. Sobre la superficie libre del líquido actúa una presión de po= 1,45
Kgf/cm2. Determinar el gasto que descarga del orificio.υ= 0,0175 cm2/s

Fig. 4.9 del Problema 4.4

x < 3a (contracción incompleta) (ambos lados) 0,15 y 0,19 < 3 x 0,18 =0,54  Contracción
Incompleta
18x 2x981x1500
h  50 cm+
1450 gr/cm2
1 gr/cm 3
=1500 cm  R = 0,0175 cm 2 /s  1, 76 x 10 6 => C do  0, 60

0,18x0,18 2
Cd  0, 60 1  0, 641 0,3x 0,69  0, 609

Q  0, 609(0, 18x0, 18 ) 2 x g x 15  0, 338 m 3 /s

 71

4.1.4. Orificio de pared delgada. Contracción total. Descarga sumergida y parcialmente

sumergida
Cuando el orificio descarga a otro tanque cuyo nivel está por arriba del canto inferior
del orificio, se dice que la descarga es ahogada. El ahogamiento puede ser total (Fig. 4.10 a)
o parcial (Fig. 4.10 b).

Fig. 4.10 Orificio con Descarga Ahogada (a) y Parcialmente Sumergida (b)

En el caso de descarga ahogada total se puede derivar una ecuación análoga a la general
(4.2) con la única diferencia que la altura h es entonces H (diferencia de niveles entre los dos
recipientes); el gasto es entonces:

Q= Cd A 2gH (4.7)

Cuando el ahogamiento es parcial, como en la figura 4.10 b, el gasto total descargado

por el orificio se puede expresar como la suma Q1 y Q2, donde Q1 es el gasto correspondiente
a la porción del orificio con descarga ahogada, es decir:

Q 1  Cd 1 A 1 2gH (4.8)

y Q2 es el gasto de la porción del orificio con descarga libre, es decir:

Q2  Cd 2 A 2 2gHm (4.9)

No hay investigaciones confiables acerca de los coeficientes de gasto Cd1 y Cd2; al

respecto, Schlag (1957) propone que Cd1= 0,70 y Cd2= 0,675, en el caso de que el orificio
tenga un umbral en el fondo, como en la figura 4.10 b.
 72

4.1.5.- Orificio de pared delgada. Contracción parcial

Una compuerta consiste en una placa móvil, plana o curva, que al levantarse permite
graduar la altura del orificio que se va descubriendo, y a la vez controlar la descarga producida.
El orificio generalmente se hace entre el piso de un canal y el borde inferior de la compuerta,
por lo que su ancho coincide con el del canal. (Fig.4.11 y 4.12) (Sotelo, 1982).

Fig. 4.11 Compuerta plana inclinada Fig. 4.12 Compuerta radial

4.1.5.1. Compuerta plana

Si se establece la ecuación de Bernouilli en una compuerta plana entre la sección 1 y
la vena contraída (sección 2) se obtiene la ecuación 4.10, la ecuación de continuidad entre
estas secciones (EC. 4.11) sustituida en la ecuación 4.10, tenemos:

(4.10)

(4.11)

C ca
2
C ca 2 V 2 V22 V22 y 1 C ca 1 y1
y1   Cc a   
y1
y1 2g 2g => 2g 1
C ca 2 1
C ca
y1 1
C ca
y1
y1

Por tanto, la velocidad media real en la sección contraída es la ecuación 4.12 en que
Cv es el coeficiente de velocidad. El gasto o caudal se obtiene al multiplicar la velocidad media
real por el área contraída (ecuación 4.13).

(4.12)
 73

(4.13)

(4.14)

(4.15)

Experimentalmente para compuertas planas verticales y descarga libre o ahogada

(Límite Superior Cv=1) (4.16)

Las figuras 4.13 y 4.14 sirven para calcular los Cd en compuertas planas

Figura 4.13. Coeficientes de gasto para compuertas planas inclinadas con descarga libre.

Figura 4.14. Coeficiente de gasto de una compuerta plana vertical según Cofré y Buchheister
(Sotelo, 1982)
 74

Ejemplo 4.5. En la compuerta rectangular (Fig. 4.13), 90º; y1= 2,4 m; b=3m. Calcular: a)
La abertura que debe tener para descargar un gasto de 7m3/s; b) Con esta misma abertura,
calcular el gasto que descarga cuando el tirante, aguas abajo, y3= 1,8 m.
Q
asumiendo Cd  0, 56 => a  Cdb 2gy1
 7
0,56x3x 2gx2,4
 0, 607

y1 2,4 Fig. 4.14 Ec. 4.16 Ec. 4.15

a  0,607  3, 95  Fig    Cd  0, 57 =>a  0, 60m      Cv  0, 985      Cc  0, 62
. 4.13

y2= 0,62x0,60 = 0,372 m.

4.1.5.2. Compuerta radial

Para compuertas radiales Knapp (Sotelo, 1982) encontró una ecuación semejante a la
(4.16) para calcular el coeficiente de velocidad donde Cv tiene nuevamente como límite
superior Cv=1. Dicha ecuación modificada es:

C v  0, 960  (0, 001615o  0, 0475 ) y 1

a
(4.17)
En la fig. 4.15 se presentan los valores de Cd en compuertas radiales con descarga libre
que son aplicables a la Ec. 4.13.

Figura 4.15. Coeficientes de gasto en compuertas radiales con descarga libre, según
Gentilini (Sotelo, 1982)
 75

El coeficiente de contracción se podría obtener por la Ec. 4.15 o por la Ec. 4.18.

C c  1  0 , 0 0 8 3 3 0  0 , 0 0 0 0 4 4 5 ( 0 )
2
(4.18)
En la Fig. 4.16 se presentan los valores del coeficiente de gasto, los cuales se han
relacionado con y1/ r, a/r, y3/r para cada valor de h/r, tanto en descarga libre como en
sumergida. Los Cc y Cv se podría obtener igual a la descarga libre.

Figura 4.16. Coeficiente de Gasto de una Compuerta Radial

En la obtención del coeficiente de contracción hay que tomar en cuenta que

Cos  hr  a
r (4.19)
y1  y1/r
a a/r
(4.20)

Ejemplo 4.6. Una compuerta radial de 4,5 m de radio y altura del perno h= 4m debe descargar
un gasto por unidad de ancho q=2,6 m3/sm, con un tirante, aguas arriba, y1= 4,5m y otro, aguas
abajo, y3=3,45m. a) Calcular la abertura de la compuerta para las condiciones de descarga
ahogada; b) Calcular el gasto, por unidad de ancho de la compuerta, con la misma abertura si
la descarga es libre, así como los coeficientes de contracción y gasto correspondiente.

y1 y3
h
r  4
4,5  0, 89; r  1; r  0, 767; asumo a
r  0, 2   Fig.4.16C   C d  0, 41 =>a  0, 675m.
Extrapolacio´n
 76

Un valor más exacto a  0,675  0, 15   F


ig.4.16C
    Cd  0, 37 =>a  0, 75m.
r 4,5 Extrapolacio´n
q  C d a 2gy 1  C d x0, 75 2g4, 5

q  0, 68x0, 75 2gx4, 5  4, 79 m3 /sm.

cos   hr  ar  0, 89  0, 167  0, 723          44 o

y de la Ec. 4.15 Cc=0,0416+

el Cc también se podría obtener por la ecuación 4.18.

4.1.6. Orificios de pared gruesa (e/h  0,67).

Cuando la pared en el contorno de un orificio no tiene aristas afiladas, el orificio es

de pared gruesa o tubo corto (Fig. 4.17).

Fig.4.17. Descarga a través de un tubo corto

 77

En este tipo de orificio se observa que el chorro, una vez que ha pasado la sección
contraída, tiene todavía espacio dentro del tubo para expandirse y llenar la totalidad de la
sección. Entre la sección contraída y la final ocurre un rápido descenso de la velocidad
acompañado de turbulencia y fuerte pérdida de energía. Análogamente a los orificios de pared
delgada el caudal se puede calcular por la ecuación 4.2 y los Cd están representados en la
Tabla 4.2.
Tabla 4.2
Coeficientes de gasto para diferentes tipos de orificios de pared gruesa
 78

Ejemplo 4.7. a) Calcular el gasto en un orificio circular de pared delgada practicado en una
pared vertical de un recipiente con D= 0,1m, desaloja un gasto Q=29,5 l/s y h=3 m; b) Calcular
el nuevo gasto si el orificio de pared delgada se reemplaza por un tubo corto de 40 cm. de
largo; c) Si se sustituye por un tubo divergente con redondeamiento de sus aristas cuyo ángulo
es = 6º.

Solución a)
0,1x 2g3 Fig. 4.6  x 0,1 2
R= 10 6  7, 67x10 5      Q  Cd Ao 2gy 1  0, 60x 4 2x9, 81x3  0, 036m 3 /s
C 0,6
d

Solución b)
 4       C d  0, 805 =>Q  0, 805x x 4
Tabla 4.2 0,1 2
e
D  40
10 2x9, 81x3  0, 0485m3 /s

Solución c)
Ta bla 4 .2 x0,1 2
e
d  4 ;   6 o
       C d  0 , 8 6 =>Q  0 , 86 x 4 2x9, 8 1x3  0, 0518m 3 /s

4.2.Vertederos
Cuando la descarga del líquido se efectúa por encima de un muro o una placa y a
superficie libre, la estructura hidráulica en la que ocurre se llama vertedor; éste puede
 79

presentar diferentes formas según las finalidades a que se destine. Así, cuando la descarga se
efectúa sobre una placa con perfil de cualquier forma, pero con arista aguda, el vertedor se
llama de pared delgada; por el contrario, cuando el contacto entre la pared y la lámina vertiente
es más bien toda una superficie, el vertedor es de pared gruesa.

4.2.1. Vertedor de pared delgada (e/h < 0,67)

Fig.4.18 Vertedor de Pared Delgada en Forma General

h= Desnivel entre la superficie inalterada del agua y la cresta.

Vo= Velocidad media de la sección inalterada del agua.

Aplicando Bernouilli entre la sección 0 y 1 de la Fig. 4.18, despreciando las pérdidas

de energía y la velocidad en la sección inalterada, tenemos la ecuación del gasto para
cualquiera geometría (Ec.4.21).

V 2o V2
h 2g y 2g  V  2g(h  y )

dQ  VdA  2g(h  y ) (xdy2 )       Q  2 2g   o x (h  y) dy

Integrando h

Q  2 2g   o x (h  y ) dy
h
(4.21)

Donde  considera el efecto de contracción de la lámina vertiente.

 80

4.2.2. Vertedor de Pared Delgada Rectangular

Si en la Ec.4.21 x= b/2, donde b es la longitud de la cresta e integramos la Ec. 4.21
tenemos:

h
Q  2 2g  o b2 (h  y ) dy  23 b 2g (h  y )
h 1/2 3/2
o

Q 2
3 2g  bh 3/2 (4.22)

La Ec. 4.22 es la ecuación general para calcular el gasto en un vertedor rectangular

cuya carga de velocidad de llegada es despreciable, motivado a esto diferentes científicos han
investigado y obtenido fórmulas experimentales para evaluar  , en los cuales se incluya dicho
efecto para permitir que tales ecuaciones conserven su validez (Tabla 4.3) (Sotelo, 1982).

Cuando el vertedor rectangular sin contracciones laterales tiene una inclinación  con
respecto a la horizontal (Fig.4.19), el coeficiente de gasto  de la tabla 4.3 debe multiplicarse
por un coeficiente C que depende del ángulo de inclinación  (Ec. 4.23).

Fig.4.19 Vertedor Inclinado

o
C   1 , 19 5 1  0 , 39 0 2 1 8 0 o (4.23)
Q  C 23 2g  bh 3/2 (4.24)

Las relaciones de la tabla 4.3 tienen ciertas limitaciones en común para su mayor
precisión:
 e  2 mm.
 La carga sobre el vertedero se debe medir en un punto, aguas arriba, por lo menos
cuatro veces la carga máxima aguas arriba.
 Contracciones equidistantes a las paredes laterales del canal.
 81

Tabla 4.3. Fórmulas experimentales para determinar el coeficiente de gasto  aplicable a la

Ec. 4.22 para vertedores rectangulares con contracciones laterales o sin ellas. En el caso de
vertedores sin contracciones laterales haga b= B, en las fórmulas.

Autor Fórmula Límites de Observaciones

Aplicación

  0, 6075  0, 045( x

Bb 0,10m h 0,60m
)
Hegly (Ref.30) 0,0041 El primer límite de aplicación es el más
(1921) B h 0,50m b 2,00m importante. Para h/b>0,13 tiene mayor precisión
0,20m w 1,13m que la fórmula SIAS

 0, 578  0, 037( B )  1000h1.6 x

Sociedad de b 3.6153 b/B
2 ( )2 0,025h0,80m Para vertedores sin contracciones laterales los
Ingenieros y b0,3B límites son:
Arquitectos Suizos w0,30m
(1924) (Ref. 9)

x 1  0 .5 ( ) ( )
0,025mh0,80m
(Fórmula SIAS) b 4 h 2 h/w1 en el caso 0,30mw
B h w de contracciones h/w1
laterales
Para h/b 0,13; es más precisa que la de Hegly

  0 .6 1 6 ( 1 
Hamilton-Smith b 0,075mh0,60m Si B(h+w)< 10 bh, se deberá reemplazar en la Ec.
10B ) 0,30mb 4.22 el valor de h por h' donde
0,30mw

hw/2

b(B-2h)

h/b0,5

es la velocidad de llegada

Francis (Ref.31) V2o 3/2 V2o 3/2 0,18mh0,50m

 0.623[1  0.1n b ] 1  2gh 
h
2gh 2,40mb3,00m
0,60mw1,50m velocidad de llegada.
b3h n = 2 en vertedores con contracciones laterales
n = 0 en vertedores sin contracciones laterales
0,01mh0,80m
Vale sólo para vertedores sin contracciones
h0.0011
 [0.6035 0.0813( w )][1  h ]
Rehbock (1929) 0.0011 3/2 b0,30m
laterales. Es muy precisa y de las más utilizadas
(Ref. 9) w0,06m
por su sencillez.
h/w1

Ejemplo 4.8. Calcular el gasto en un vertedor rectangular de pared delgada en un canal del
mismo ancho de la cresta b=2,5m, que trabaja con una carga h=0,42 m, cuya cresta se
encuentra a W=1m del piso del canal.

Solución: Utilizando la fórmula de Hegly y Rehbock de la Tabla 4.3

0,42 2
  0, 6075 
0,0041
(Hegly) 0,42 x 1  0, 55 0,421  0, 647

Q  2g x0, 647x2, 5x0, 42 3/2  1, 299m3 /s

2
3

0,420,0011 0,0011 3/2

  0, 6035  0, 0813 1  0, 640
(Rehbock) 1 0,42

 Q  1 , 2 8 6 m /s
3
 82

4.2.3. Vertederos de Pared Delgada con Descarga Sumergida

Cuando es sumergida la descarga (existe h' abajo) de los vertederos de pared delgada,
de cualquier forma (rectangular, triangular, trapecial, parabólico y proporcionales) la ecuación
4.25 proporciona un método simple para evaluar el efecto de la sumersión.

(4.25)

En dicha ecuación, Q es el gasto del vertedor con descarga sumergida; Q1 el gasto del
mismo vertedor con igual carga en el supuesto de descarga libre, S es la relación de sumersión
(h'abajo/ h arriba); y n el exponente de la carga h sobre la cresta en la ecuación correspondiente
a Q1 (3/2 para vertedor rectangular, 5/2 para triangular, etc.). La carga de sumersión deberá
medirse desde la superficie, aguas abajo (fuera de la zona de disturbios) hasta la cresta.

4.2.4. Vertederos de Pared Gruesa Rectangular (e/h  0,67)

Para calcular el caudal en un vertedor rectangular de pared gruesa y aristas vivas (Fig.
4.20) se utiliza la fórmula 4.22 afectada por un factor 1. (Ec. 4.26).

Fig. 4.20 Vertedor rectangular de pared gruesa

Q=1 3 2g  bh  1 C bh
2 3/2 3/2
(4.26)

Donde C es el coeficiente adecuado de un vertedor de pared delgada sin contracciones

laterales y descarga libre que se evalúa de la Tabla 4.3. El coeficiente 1 depende de la relación
e/h y se obtiene de la Fig. 4.21, si la descarga es ahogada la fórmula 4.27 sirve para obtener
dicho gasto y 2 se obtiene de la fig.4.22.
 83

Q  12 23 2g  bh 3/2  1 2 C bh 3/2 (4.27)

Figura 4.21. Coeficiente 1 de la Ec. 4.26 para vertedores de pared gruesa con descarga libre,
según F. Dominguez (1959).

Figura 4.22. Coeficiente 2 para vertedores de pared gruesa con descarga ahogada, según F.
Dominguez (1959).
 84

4.2.5. Vertedores con Cresta Redondeada

Cuando la cresta del vertedor se redondea, el coeficiente de gasto C aumenta

considerablemente respecto del calculado para uno de pared gruesa. Esto se explica por una
baja en la contracción del chorro sobre el vertedor, pues actúa sobre las partículas una
aceleración centrífuga debido a la curvatura de las líneas de corriente.

Figura 4.23 Vertedor de cresta Figura 4.24 Vertedor de cresta

circular y talud inclinado, circular y talud vertical, aguas
aguas abajo abajo

En el caso del vertedor de la Fig. 4.23 el coeficiente  de la Ec.4.22 será:

  0, 312  0, 09 w  0, 3  0, 01(5  h/r )

h 2
(4.28)

que vale para:

h / r  (6 0 
2 0r
w 3 r )
h/w 1 w>r>0,2m

En el caso de vertedores como el de la Fig. 4.24, Kramer encontró

experimentalmente:

  1, 02 
1,015
h/r2,08  0, 04(h/r  0, 19 ) 2  0, 0223 r/w (4.29)

que vale para:

h/r 4,2 ; h/r 0,4+0,32 w/r+ 0,06 r  12, 5

 85

Este tipo de vertedor es importante en el caso de compuertas cilíndricas con vertido

superior.

Figura 4.25.Vertedor de cresta elíptica y talud inclinado, aguas abajo

Para vertedores como el de la Fig. 4.25, Kramer propone el empleo de la ecuación

anterior (4.29), considerando que el radio r en ella se calcule con la expresión:

(4.30)

donde a y b son las longitudes de los ejes de la elipse. Esta fórmula ha sido verificada para
valores de a=b/2, a=b, a=2b y a=6b.
 86

CAPITULO V
Flujo Gradualmente Variado en Canales

5.1 Introducción.
Como ya se ha señalado, el flujo uniforme no es representativo de las condiciones
reales de un canal, ya que existen ciertas fuerzas oponentes al movimiento que originan que
el flujo de agua no sólo no sea paralelo al fondo, sino que además tendrá diferentes formas y
posiciones con el tiempo cuando el flujo no es permanente.

Generalmente lo que se hace es suponer un flujo permanente gradualmente variado

diseñado para el gasto máximo.

Un flujo gradualmente variado es aquel en el cual el tirante o profundidad varía

gradualmente a lo largo del canal, permitiendo suponer paralelismo en las líneas de corriente.
La suposición fundamental en que se basa el análisis es suponer que las pérdidas de energía
en un tramo corto del canal son las mismas que en un flujo uniforme tomando para el cálculo
los promedios de los radios hidráulicos y las velocidades en los extremos del tramo.

Cualquier proceso de cálculo de flujo gradualmente variado requiere conocer

previamente las profundidades normales (yn ó yo ) y crítica (yc ). La primera la podemos
encontrar por la conocida ecuación de Manning (Ec. 2.5) y la segunda por la (Ec. 2.45) que
define al flujo crítico que simplificada para una sección rectangular tenemos yc = (q2 / g)1/3.

5.2 Curvas de Remanso y Tipos de Perfiles.

En un canal pueden existir secciones donde la profundidad del agua está forzada a
tener un valor diferente al tirante normal. En estos casos el perfil de la superficie líquida
variará desde estos valores a la altura normal. Esta forma adoptada por la superficie libre del
agua donde ocurre la transición entre los dos valores se denominan curva de Remanso
correspondiente a un flujo gradualmente variado.

Se ha comprobado experimentalmente, que de acuerdo a las relaciones entre yn y yc la

variación de y con respecto a x (Fig. 5.1) sólo pueden ocurrir de acuerdo a ciertas formas o
perfiles de flujo en relación a los siguientes tipos de pendientes: Pronunciadas (S), Suaves
(M), Críticas (C), Adversas (A) y Horizontales. Ver figura 5.2
 87

Figura 5.1.

Figura 5.2. Clasificación de los Perfiles de Flujo Gradualmente Variado

 88

Pendiente Pronunciada (S): Corresponde al caso en que la altura normal es menor que la altura
crítica, o sea que si el nivel fuera uniforme el régimen sería supercrítico. Las líneas de fondo,
altura normal y altura crítica definen tres zonas donde pueden formarse las curvas de remanso,
que se denominan según el tipo de pendiente (S) y con sub-índice según la zona donde se
desarrollen dependiendo de la altura real del flujo.

Pendiente Suave (M): Correspondería en flujo uniforme a un régimen subcrítico o lento,

debido a que la línea de profundidad normal estará por encima de la línea crítica. Al igual que
el caso anterior, habrá 3 zonas donde podrán desarrollarse perfiles de flujo.

Pendiente Crítica (C): Estos perfiles representan las condiciones de transición entre los
perfiles M y S. En este tipo de pendiente yn = yc y sólo existirán dos perfiles posibles.

Pendiente Horizontal (H): ( S0 = 0 y yn =Estos son los casos limitantes de perfiles M

cuando el fondo del canal se vuelve horizontal. Los perfiles H2 y H3 corresponden a los
perfiles M2 y M3 , pero un perfil H1 no puede establecerse en realidad, debido a que yn es
infinito.

Pendiente Adversa (A): (S0 < 0 y yn = no existe). Se caracteriza porque la pendiente del canal
sube en la dirección del flujo. Tampoco puede originarse un flujo uniforme y sólo la
profundidad crítica servirá de referencia para establecer la posición de una curva de remanso.
Giles (1994) presenta la tabla 5.1 donde se presenta un resumen de los tipos de variantes de
flujo. 

Con el propósito de facilitar el análisis de continuidad en los perfiles superficiales, se

presentan en la fig. 5.3 los casos más frecuentes en la variación de los perfiles superficiales
cuando existe un quiebre en la pendiente del terreno. (Aguirre, 1980).
 89

Tabla 5.1. Tipos de Variantes de Flujo

dy Prof. En el
Pendiente del Relación de
sent. De la Simbolo Tipo de flujo Forma del perfil
canal profundidad dl corriente

y>yN>yc + Aumenta M1 Subcrítico

Suave
yN > y >yc - Disminuye M2 Subcrítico
0<S<Sc

yN > yc >
+ Aumenta M3 Supercrítico
y

y > yc - Disminuye H2 Subcrítico

Horizontal
S=0
yN = 
yc>y + Aumenta H3 Supercrítico

Subcrítico
y>yc=yN + Aumenta C1
Crítica
SN= Sc Uniforme
yc=y=yN Constante C2
yN= Yc crítico
yc=yN>y + Aumenta C3 Supercrítico

yc>y>yN + Aumenta S1 Subcrítico

Pronunciada
S>Sc>0 yc>y>yn - Disminuye S2 Supercrítico

yc>yN>y + Aumenta S3 Supercrítico

Adversa y>yc - Disminuye A2 Subcrítico

S< 0
yn =  yc > y + Aumenta A3 Supercrítico
 90

Fig 5.3 Perfiles de flujo gradualmente variado en canales prismáticos de gran longitud, con
quiebre en la pendiente de fondo.
 91

Fig 5.3 Perfiles de flujo gradualmente variado en canales prismáticos de gran longitud, con
quiebre en la pendiente de fondo (Continuación).(Chow.1994 y Bolinaga.1985)
 92

5.3 Ecuación Dinámica del Flujo Gradualmente Variado.

Se define el flujo gradualmente variado como un flujo permanente y no uniforme que

sufre pequeñas variaciones de profundidad y velocidad con la posición, por efecto del
desequilibrio entre las fuerzas de rozamiento y las fuerzas de gravedad.

En el flujo permanente y uniforme la pendiente de la línea de energía S, coincide con

la pendiente S0 del fondo del canal, mientras que en flujo gradualmente variado, las
variaciones de profundidad y velocidad establecen una pendiente S diferente.

Fig. 5.4 Esquema de definición de flujo gradualmente variado.

Si, en la fig 5.4 designamos por z la elevación del fondo con respecto a un datum
arbitrario, la energía total H, de la ecuación de Bernoulli estaría dada por:

H  z  y cos   A22g
Q2
(5.1)

donde  es el ángulo de inclinación del fondo,  es el coeficiente de energía, y es la

profundidad tomada perpendicularmente al fondo, Q es el gasto y A la sección transversal del
canal en estudio. La derivada de la ec. 5.1 con respecto a la dirección x correspondiente a la
dirección del fondo del canal, se obtiene:
 93

 cos  dx  gA3
dy Q2
dH
dx  dz
dx
dA
dx
(5.2)

en donde dH/dx = - S; dz/dx = - S0 ; y , dA/dx = T dy/dx. Así, podemos escribir:

𝑑𝑦 𝑆0 −𝑆
= 𝑄2 𝑇
(5.3)
𝑑𝑥 𝑐𝑜𝑠∅−𝛼
𝑔𝐴3

𝑑𝑦 𝑑𝐸/𝑑𝑥
y también, = 𝑄2 𝑇
(5.4)
𝑑𝑥 𝑐𝑜𝑠∅−𝛼
𝑔𝐴3

debido a que S0 - S = d(H - z )/dx = dE/dx.

Las ecs. 5.3 y 5.4 conocidas como las ecuaciones de resistencia del flujo gradualmente
variado, son ecuaciones diferenciales en y = f (x), que en general no son explícitamente
solubles, pero existen numerosos métodos numéricos para obtener soluciones. El problema
principal con estas ecuaciones, es que requieren modificaciones para mostrar la variación
cualitativa del perfil superficial. Es deseable desarrollar un proceso sistemático que nos
permita definir el tipo de variación que sufrirá un perfil superficial, si en una determinada
sección se conocen sus características hidráulicas. Es frecuente recurrir a los canales
rectangulares de gran anchura, con distribución uniforme de velocidades tal que  = 1 y con
inclinación tal que cos 1 para el análisis cualitativo de los perfiles superficiales. Entonces
la pendiente S de la línea de energía puede obtenerse de la ecuación de Manning, así

q2 n2
S= y 10/3 (5.5)

donde q es el gasto por unidad de ancho. Así la ec. 5.3 para las condiciones señaladas, se
puede escribir:

q 2n 2
dy S0
y10/3
dx = q2
1
gy3 (5.6)

Introducimos la profundidad crítica yc en q2/g = yc3 ; así :

1- q2/(gy3) = 1-(yc /y)3. Por lo tanto la ec. 5.6 se puede entonces transformar en:
 94

q 2n2
dy S0 - y10/3

dx  yc3 (5.7)
1- y3

En la cual S0 es también la pendiente de la línea de energía del flujo uniforme de

profundidad normal yn . Para condiciones de flujo uniforme, de la ecuación de Manning se
obtiene:

q2 n 2
S0 = y 1n0/3 (5.8)

En esta ecuación se puede despejar q2n2 , y reemplazarlo en la ec. 5.7 para obtener:

yn 10/3
dy 1 y
dx = S0 yc 3
1 y (5.9)

Esta ecuación permite estudiar el cambio longitudinal de la profundidad en función de

tres parámetros fundamentales: la pendiente del fondo, la relación entre la profundidad normal
del flujo uniforme y la profundidad existente, y la relación entre la profundidad crítica y la
profundidad existente.

El interés primario en el análisis cualitativo del perfil superficial reposa en el signo de

𝑑𝑦
, ya que si él es positivo la profundidad crece hacia aguas abajo y si es negativo la
𝑑𝑥
profundidad decrece en la dirección de aguas abajo. Se puede mostrar que la variación de los
perfiles superficiales corresponde a lo ilustrado en la Fig. 5.5 en la cual se presentan las curvas
posibles para pendiente de flujo normal subcrítico dando lugar a los perfiles M; para pendiente
de flujo normal supercrítico originando los perfiles S; para pendiente crítica los perfiles C,
pendiente adversa con perfiles A y para pendiente horizontal con S0 = 0 los perfiles H. (Ver
tabla 5.1 y Fig. 5.6).
 95

Fig. 5.5 Perfiles Superficiales en flujo gradualmente variado.

El análisis cualitativo de la variación de profundidades, en un canal dado, se realiza de

la siguiente forma: Se calculan los parámetros yn e yc para el gasto y las características del
canal dados y se grafican en un diagrama de perfiles en que las escalas de las longitudes en la
dirección del movimiento están reducidas con respecto a las de profundidades. Luego se
calculan las profundidades en las secciones cuyas características de flujo están
predeterminadas, es decir, en las secciones de control, tomas, vertederos, compuertas, caídas
libres y quiebres en la pendiente del canal. Se toma en cuenta que los perfiles sobre las líneas
de profundidad crítica están determinados por los controles de aguas abajo; y que los perfiles
bajo las líneas de profundidad crítica están determinados por los controles de aguas arriba.

Luego se dibujan, sobre el diagrama del canal, las porciones apropiadas de los perfiles
de la Fig. 5.5. Si la profundidad tiene que pasar de un valor supercrítico a otro subcrítico, se
genera un resalto hidráulico situado en posición compatible con la relación de profundidad
propia del resalto. Después de haber resuelto el problema cualitativamente se procede a la
integración de la ecuación diferencial del flujo gradualmente variado para determinar y = f(x).
 96

Figura 5.6. Ejemplos de Perfiles de Flujo (Chow, 1994)

5.3.1 Secciones de Control:

Cuando un conducto prismático está formado por diferentes secciones o cambios de
pendientes se presentarán distintos tipo de perfiles de flujo, que unidos forman lo que se
conoce como perfil de continuidad de un canal. Para conocer el perfil de continuidad es
necesario determinar los distintos puntos o secciones de control y la relación entre yn y yc .
 97

Una sección de control es aquella donde se presuma que existirá una altura de agua
conocida. Entre las secciones de control más comunes están:

a) Controles de Yn : Todo canal que tienda al infinito, (canal largo) tenderá a flujo uniforme.
Asimismo, todo aquel que provenga del infinito, vendrá de flujo uniforme, basado en esto:

a.1) Extremo de aguas arriba de un tramo largo de una pendiente M, donde se conoce que
"y" estará muy cerca de yn (Fig 5.7b).
a.2) Extremos de aguas abajo de un tramo largo de una pendiente S, donde la altura del
agua tendrá el valor de la altura normal. En la fig 5.7b se puede observar que las
curvas de remanso tienden a coincidir con la línea de altura normal aguas abajo,
mientras para una pendiente M, lo son aguas arriba. En ambos casos, las curvas
tienden a ser perpendiculares a la línea de altura crítica (comprobar con la ec. 5.9).

b) Controles de yc: Todo sitio donde el flujo pase con su energía específica mínima (fig 2.8),
poseerá profundidad crítica. Los lugares donde esto puede ocurrir son: los escalones
(ejemplo 2.11), los estrechamientos, las caídas y los cambios de pendiente M a S (fig 5.7b).

c) Controles artificiales: Son mecanismos u obstáculos donde se obliga al líquido a tener una
altura determinada que se fija a voluntad. Los más usuales son los vertederos, los orificios,
las compuertas de fondo, etc. (fig. 5.7c).

d.) Controles de niveles: Existen sitios donde los niveles del líquido son conocidos para un
gasto determinado, por ejemplo, una estación de aforos, o bien un cuerpo estacionario de
agua como un embalse, un lago o el mar (fig 5.7a).

Figura 5.7 Secciones de Control en Canales

 98

Las perturbaciones producidas por una sección de control sólo podrán viajar aguas
arriba si la velocidad del flujo es menor que la celeridad, en caso contrario influirá aguas
abajo. En otras palabras, cuando el flujo sea supercrítico (y < yc ó V > celeridad Velocidad
crítica) las perturbaciones no influirán aguas arriba y el cálculo se iniciará desde aguas arriba
hacia aguas abajo. Si en cambio el flujo es lento (y>yc ó V<Vc ), cualquier control situado
aguas abajo influirá en la altura del agua, aguas arriba; por lo que el cálculo se iniciará desde
aguas abajo.

Con relación a lo anteriormente señalado, cabe observar que en un flujo gradualmente

variado, el tipo de régimen (subcrítico o supercrítico) no dependerá del tipo de pendiente
(suave o pronunciada) ni de la relación de yn a yc , sino de la relación de la altura real del flujo
a la altura crítica. Así, un perfil M3 que se desarrolla en una pendiente suave (yn >yc )
corresponde a un régimen supercrítico ya que la altura real del agua (y) será menor a la altura
crítica; por su parte, un perfil S1 corresponderá a un flujo subcrítico a pesar de originarse sobre
una pendiente fuerte (yc > yn ).

Como conclusión, en la fig 5.2, se pueden apreciar los puntos de partida para el cálculo
de cada uno de los perfiles según el flujo sea subcrítico o supercrítico.

5.4 Integración de la Ecuación Diferencial del Flujo Gradualmente Variado.

La ecuación general del flujo gradualmente variado, ec. 5.3, se puede integrar para
obtener la distancia x1 - 2 necesaria para pasar de una profundidad y1 a otra profundidad y2 .

En efecto, si reescribimos la ec. 5.3 para cosy = 1 en la forma

1 - Q2 T / (gA3 )
dx= S0 - n 2 Q 2 / (A2 R4/3 ) dy
(5.10)

y si se integra entre y1 , e y2 se tiene que

1Q2 T/ gA3 
x 12   y21
y

S 0 n2 Q2 /(A2R4/3 )
(5.11)

Cuando el numerador de la función integral es cero, el flujo es crítico, indicando que

no hay cambio en x al variar y, despreciando la curvatura de las líneas de corriente y la
distribución de presiones no hidrostática en la sección. Cuando esto ocurre no se puede decir
que el flujo sea gradualmente variado y las fórmulas anteriores no dan resultados que
 99

describan con precisión la variación del flujo cuando éste está próximo a su valor crítico.
Cuando el denominador de la función integral es cero, el flujo es uniforme y la profundidad
alcanza su valor normal.
Para un canal de sección transversal no variable, con valores de n y s0 constantes, la
función integral depende únicamente de y. Si esta función se denomina F(y) tal que

1Q2 T/ (gA3 )
F (y) = S0 n2 Q2 /(A2 R4/3 )

entonces se tiene que

X 12   y1
y2
F (y) d y (5.12)

La representación gráfica de F (y) en un esquema como el de la fig. 5.8 nos permite

hallar el valor del área bajo la curva F (y) entre y1 e y2 la cual es igual a la distancia x1-2.
Evidentemente, se puede hallar una aproximación numérica de la integral dividiendo el área
en prismas de base y, y altura F (y) que permitan obtener una sumatoria que será tanto más
precisa cuanto menor sean los incrementos y.

Fig. 5.8 Integración gráfica de la ecuación del flujo gradualmente variado.

Ejemplo 5.1. Bajo la compuerta que se indica en la figura, fluye un caudal de 7m3/seg/m hacia
un canal horizontal muy ancho, con una rugosidad n= 0,017. Calcular la distancia a la que
debe ubicarse la caída, para que la profundidad antes del resalto sea de 1,10 m.

La profundidad crítica yc está dada por

yc  3 q 2 /g  3
7 2 /9, 81  1,71 m .
 100

Aguas abajo de la compuerta se produce un perfil H3, luego un resalto hidráulico y un perfil
H2 hasta la profundidad crítica en la caída.

La profundidad secuente del resalto hidráulico es, según la ec. 5.17.

y1
y2  1  1  8(F1 ) 2  1  1  8(7 2 /9, 81 x 1, 1 3 )
1,10
2 2

y2 = 2,51 m

Para la determinación del perfil será necesario desarrollar la ecuación diferencial

5.10:

1 - Q2 T / (gA3 )
dx= S 0 - n 2 Q2 /(A2 R4/3 ) dy

Integrando entre y1 e y2 , se tiene:

1Q2 T/gA3 
x12   y21
y
S 0 n 2 Q2 /(A2 R4/3 ) dy
 101

Para canales horizontales de gran anchura S0 = 0 y R = y, quedando:

y 2 1q 2 gy 3
x 12    y 1
y 1/3 y 10/3
/
n 2q 2/y 10/3 dy   yy12 gn 2  n2q2 dy

resolviendo la integral y sustituyendo límites, se tiene:

Aplicando la ecuación anterior entre los puntos y1 = 0,8 m e y2 = 1,10 m se determina

la distancia horizontal existente entre la compuerta y el resalto hidráulico:

x12   0, 84/3 )   0, 8 13/3 )

3 4/3 3 13/3
(
4x0,0172 x9,81 1, 10
(
13x0,0172 x72 1, 10

= 85,49 m

La caída constituye una sección de control, con profundidad crítica, de forma que
volviendo a aplicar la ecuación anterior y3= 2,51 m e y4 = 1,71 m se establece la distancia
horizontal entre el resalto y la caída:

x34   2, 514/3 )  13x0,03172 x72 (1, 7113/3  2, 5113/3 )

3 4/3
4x0,0172 x9,81 (1, 71

= 350,97 m

En definitiva, la longitud total del canal entre la compuerta y la caída debe ser
(despreciando la longitud del resalto hidráulico):

xtotal = 85,49+350,97 = 436,46 m

5.5. Métodos Numéricos para el Cálculo del Flujo Gradualmente Variado

Se puede realizar una aproximación numérica para el cálculo de los perfiles
superficiales mediante una tabulación sistemática. Una de las formas más difundidas
corresponde al llamado método de las aproximaciones sucesivas o método directo por etapas
el cual supone que las pérdidas de energía en un tramo corto del canal son las mismas
 102

que en un flujo uniforme tomando para el cálculo los promedios de los radios hidráulicos y
las velocidades en los extremos del tramo. Aplicando Bernouilli la sección 1 y 2 de la figura
5.9 tenemos:

Fig 5.9 Esquema de definición.

_
Como Z1 = LS0 y hf = LS ; sustituyendo

V2 V 21
__ y 2  2g2  y 1 
V21 V22 2g E 2 E1
LS 0  y 1  2g  y2  2g  L S ; L  __
S 0 S
 __
S0  S

E2 E1
L  __
S0 S (5.13)

El método de cálculo se puede resumir de la siguiente manera:

1.- Se establece el tipo de perfil según los controles y según la variación de yn y yc.
2.- Comenzando en una sección de y conocida se divide el tramo en sub-tramos con
incrementos o disminuciones de y según el tipo de régimen en que se encuentre el perfil; el
cálculo se iniciará hacia aguas arriba para regimenes subcríticos o hacia aguas abajo para
regimenes supercríticos.
3.- Se estima la velocidad del agua para la altura establecida (V=Q/A).
4.- Se obtiene E = y+V2/2g para cada valor de y.
5.- Usando la Ec.de Manning se obtiene el valor de la pendiente S para cada valor de la altura

del agua .

6.- Se obtiene el valor de para un intervalo de longitud L dado por

 103

7.- Se obtiene S0 - y E2-E1.

E2 E1
8.- Luego L 
__
S0 S , la suma de los valores de L dará la longitud total entre los límites de
variación de y.

Ejemplo 5.2 En un canal rectangular muy largo ocurre un quiebre de pendiente, aguas arriba
la pendiente es suave y aguas abajo rápida. Si Q = 31 m3/s, b=2,5m; n = 0,03; S0= 0,077 (aguas
abajo) determinar y estudiar la curva de remanso originada aguas abajo del quiebre de
pendiente.

(2,5y ) 2/3
31  1
0,03 (2, 5y ) 2,52y (0, 077 ) 1/2 y n  1, 675 m

q2 1/3 (31/2,5 ) 2 1/3

yc  g  9,81  2, 5 m  y c  y n  r. sup ercr´itico

La curva que se originará será del tipo S y se extenderá desde el control donde y = yc
hasta aproximarse a la altura normal, luego el perfil será del tipo S2. El flujo es supercrítico y
las perturbaciones no influirán aguas arriba. Por lo cual el cálculo se iniciará desde aguas
arriba hacia aguas abajo.
 104

Puede verse que la distancia más próxima a la real es la de 45,69 m ya que se discretiza
en el cálculo del perfil lo cual implica una pendiente en el intervalo (S) más próxima a la real
y por ende una distancia L más precisa.

Ejemplo 5.3 Un río de fondo ancho y poco profundo, casi rectangular, tiene una pendiente
uniforme de 0,0004 y presenta una profundidad normal de 3 m. ¿Determinar la curva de
remanso producida por una represa que origina una profundidad de 7,5 m. n = 0,022?

Vby q2 1/3 5,67 2 1/3

 Vy  1, 89 x 3  5, 67 m 3 /Sm  y c 
Q
q b  VA
b  b g  9,81

= 1,48 m

Como yn> yc la pendiente será del tipo M, y debido a la altura del agua en la represa
será del tipo M1 donde y > yn > yc, entonces la represa será una sección de control para el
régimen aguas arriba y el cálculo se iniciará de aguas abajo (y = 7,5 m) hacia aguas arriba (y
= 3 m). Dicho cálculo del perfil se podría hacer ya que q = q/b=vy y R=y tanto con la ecuación
5.11 como con la ecuación 5.13 (ver ejemplo 5.1)

Ejemplo 5.4. La profundidad aguas arriba es 0,5 m y la velocidad V = 10 m/s. El canal es

muy ancho y poco profundo (rectangular) y tiene una contrapendiente S0=0,005 ¿Calcular la
profundidad en el otro extremo? n=0,015
 105

q  V 1 y 1  1 0 x 0, 5  5 mSm  E 1  0 , 5  10 2
 5, 6 m ;
3
2g

2 105/y 2 2
__ __
0,015
S 
nV 2
__ 2/3
R 0,5y 2 2/3
2 ; sustituyendo en la Ec. 5.13

52
y 2 5,6
2gy 22 Aplicando
45  105/y 2 2                   y 2  0, 64 m
M. Newton
0,015 2
0,005
0, 5y 2 2/3
2

5 2 1/ 3
yc  9, 81  1, 36 m  y 2  No hay resalto hidra´ uli co

Ejemplo 5.5. Un canal trapezoidal de n= 0,015; Q= 15 m3/seg; b= 3m; m= 1 y So=0,001, este

en su final experimenta un cambio brusco de pendiente (Tipo S). ¿Conocer y dibujar el perfil
de la superficie libre aguas arriba?

(3yy 2 ) (3yy 2 ) 2/3 Aplicando el

15  0,015 (32y 2 ) (0, 001 ) 2  M. Newton  y n  1, 58m

Q2 A3c (3y cy 2c ) 3 Aplicando el

 91,581 
2
g  Tc (32y c )        yc  1, 19m  y n Tipo M
M. Newton
 106

Lo que implica aguas arriba un perfil M2 y el cálculo de perfil se comenzará desde la

sección de control yc = 1,19 m hasta yn= 1,58 m.

y A P R ¯R V V2/2g ¯V E S E2-E1 ¯S S0- ¯ S L L

(m) (m2) (m) (m) (m) (m/s) (m) (m/s) (m) (m) (m) (m)
1,19 4,98 6,36 0,78 0,81 3 0,46 2,84 1,65 0,0028 0,0157 0,0024 0,0014 11,12 11,12
1,30 5,59 6,67 0,83 2,68 0,36 1,66 0,002
1,40 6,16 6,95 0,88 0,86 2,43 0,30 2,55 1,70 0,0015 0,0352 0,0018 0,00079 44,10 55,22
1,50 6,75 7,24 0,93 0,90 2,22 0,25 2,32 1,75 0,0012 0,0494 0,0013 0,00038 127,99 183,22
1,58 7,23 7,46 0,96 0,95 2,07 0,21 2,14 1,79 0,001 0,0473 0,0011 0,00011 428,17 611,4

5.6. Resalto Hidráulico.

El flujo gradualmente variado podría convertirse en flujo fuertemente variado si el

tirante varía en forma abrupta en una corta longitud del canal. En la figura 5.10 (a), se observa
un cambio brusco de una pendiente suave a una pronunciada. La línea crítica será paralela al
fondo aguas arriba y la altura del agua teóricamente será mayor a la crítica y la energía mayor
a la mínima, el tirante disminuye a medida que el flujo avanza y la energía disminuye por
efecto de fricción hasta llegar a ser igual a la energía mínima; esto ocurre cerca del cambio de
pendiente, de allí en adelante el tirante continúa disminuyendo con el incremento de la
pendiente, el cual abastece de mayor energía al flujo, a la que disipa por fricción. De esta
manera, se concluye que la transición de un régimen subcrítico a uno supercrítico es gradual,
acompañada de poca turbulencia y pérdida de energía debido sólo a la fricción en el
movimiento.

Fig. 5.10

Estudiaremos ahora el proceso inverso de transición de régimen supercrítico a

subcrítico representado en la figura 5.10 (b). Para explicar el proceso de transición,
podemos recurrir a un análisis semejante al del caso anterior. El flujo inicialmente en el
 107

régimen supercrítico se frena por efecto de la fricción y disminución de pendiente hasta

alcanzar la condición crítica. Toda vez, que aguas abajo existe un régimen subcrítico. Esto
implicaría que la energía específica del flujo, debe ser mayor que la mínima.

Para encontrar la forma de transición, se puede apelar a la evidencia experimental, la

cual demuestra con toda claridad que por el contrario al caso anterior, la transferencia del
régimen supercrítico a subcrítico, es violenta y se acompaña de mucha turbulencia y gran
pérdida de energía. Después de un crecimiento irregular y brusco de la superficie del agua
hasta alcanzar un tirante aproximadamente igual a la altura normal de un tramo relativamente
corto, el frente turbulento se normaliza de manera inmediata y continúa libremente un régimen
subcrítico.

El fenómeno antes descrito se conoce como Onda Estacionaria o Resalto Hidráulico.

Además de su gran mérito como disipador natural de energía , el salto hidráulico tiene muchos
otros usos prácticos entre los cuales se puede mencionar:
a) Prevención de la socavación aguas abajo de las estructuras hidráulicas donde sea necesario
disipar la energía. Ejemplos: Caídas, canales con grandes pendientes, salidas de alcantarillas,
etc.
b) El mezclado eficiente de fluidos o de sustancias químicas usadas en la purificación del
agua (Plantas de Tratamiento).
c) El aireamiento del agua destinada al abastecimiento (Plantas de Tratamiento).
d) La recuperación de cargas aguas abajo de un aforador.

5.6.1. Características.
Un resalto hidráulico tiene varias características: Sus profundidades alternas ó
conjugadas, su pérdida de energía, su longitud, su altura y su eficiencia:

5.6.1.1. Profundidades alternas o conjugadas.

La profundidad aguas arriba en régimen supercrítico, tiene su conjugada en el régimen
subcrítico aguas abajo. Estas alturas son función del número de Froude y de la geometría del
conducto donde se produce el resalto cuando la pendiente de fondo no influye sensiblemente
en los cálculos.

Para el caso general de cualquier sección prismática, la ecuación que expresa la

relación de las alturas de un resalto puede obtenerse a partir de las ecuaciones de continuidad
y de cantidad de movimiento.
 108

De la Segunda Ley de Newton, podemos obtener la ecuación de cantidad de

movimiento, aplicada a un fluido de un volumen de control:


V V
F m t
 
t
 QV  V (V d A) (5.14)

Figura 5.11. Esquema de definición para determinar la relación de Profundidades en resalto

Hidráulico con So= 0.

que aplicada a la figura 5.11 y recordando = /g y el fluido es permanente:

F 1  F 2  Q
Q2
A2 
Q1
A1
 A 1 y 1  A 2 y 2  
g
Q2
A2 
Q2
A1 , entonces :

(5.15)

La ecuación 5.15 representa la ecuación general que relaciona las alturas conjugadas
en un resalto hidráulico para una sección cualquiera.

La tabla 5.2 tomada de Bolinaga (1979), se puede utilizar para obtener ¯y en canales
trapezoidales.
 109

Tabla 5.2 Distancia del Centro de Gravedad a la superficie en canales trapeciales

Z=m y¯/y
y/b Z= 1/8 Z= ¼ Z= ½ Z= ¾ Z=1 Z=1 ½ Z=2 Z=2 ½ Z=3 Z=4
0,05 0,499 0,498 0,496 0,494 0,492 0,488 0,485 0,481 0,478 0,472
0,10 0,498 0,496 0,492 0,488 0,485 0,478 0,472 0,467 0,462 0,462
0,15 0,497 0,494 0,488 0,483 0,478 0,469 0,462 0,455 0,448 0,438
0,20 0,496 0,492 0,485 0,478 0,472 0,462 0,452 0,444 0,438 0,426
0,25 0,495 0,49 0,481 0,474 0,467 0,455 0,444 0,436 0,429 0,417
0,30 0,494 0,488 0,478 0,469 0,462 0,44 0,438 0,429 0,421 0,409
0,35 0,493 0,487 0,475 0,465 0,457 0,443 0,431 0,422 0,415 0,403
0,40 0,492 0,485 0,472 0,462 0,452 0,438 0,426 0,417 0,409 0,397
0,45 0,491 0,483 0,469 0,458 0,448 0,433 0,421 0,412 0,404 0,393
0,50 0,49 0,481 0,467 0,455 0,444 0,429 0,417 0,407 0,40 0,389
0,60 0,488 0,478 0,462 0,44 0,438 0,421 0,409 0,40 0,398 0,382
0,70 0,487 0,475 0,457 0,443 0,431 0,415 0,403 0,394 0,387 0,377
0,80 0,485 0,472 0,452 0,43 0,426 0,409 0,397 0,389 0,382 0,373
0,90 0,483 0,469 0,448 0,435 0,421 0,404 0,393 0,385 0,378 0,37
1,00 0,481 0,467 0,444 0,429 0,417 0,4 0,389 0,381 0,375 0,367
1,10 0,48 0,464 0,441 0,425 0,413 0,396 0,385 0,378 0,372 0,364
1,20 0,478 0,462 0,438 0,421 0,409 0,393 0,382 0,375 0,37 0,362
1,30 0,477 0,459 0,434 0,418 0,406 0,39 0,38 0,373 0,367 0,36
1,40 0,475 0,457 0,431 0,415 0,403 0,387 0,377 0,37 0,365 0,359
1,50 0,474 0,455 0,429 0,412 0,4 0,385 0,375 0,368 0,364 0,357
1,60 0,472 0,452 0,426 0,409 0,397 0,382 0,373 0,367 0,362 0,356
1,70 0,471 0,45 0,423 0,407 0,395 0,38 0,371 0,365 0,361 0,355
1,80 0,469 0,448 0,421 0,404 0,393 0,378 0,37 0,364 0,359 0,364
1,90 0,468 0,446 0,419 0,402 0,391 0,377 0,368 0,362 0,358 0,353
2,00 0,467 0,444 0,417 0,4 0,389 0,375 0,367 0,361 0,357 0,352

Ejemplo 5.6. Un canal trapecial de b= 3 m, m= 1, conduce un gasto de 8 m3/seg a una

profundidad de 0,3 m. Se desea calcular la altura conjugada.

__
__
De la tabla 5.2  by  0,3
3  0, 1 y m=1 
y
y  0, 485  y  0, 145m

8,08
A  3 x 0 , 3  0, 3 2  0, 99m 2 F V
  4, 92  1 (R.Ra´ pido)
g AT 9,81x
0,99
3,6
 110

T  3  2 x 1 x 0, 3  3, 6 m

Q2
g 
A3c
Tc
y c  0, 8 m

__ Q2
0, 145 x 0, 99   y 2 A2 
Aplicando la Ec. (5.15) 82
9,81x 0,99 gA2

por Tanteo: y2 = 1,69 m

¯y2 = 0,745 m
A2 = 7,93 m2

8,08 2 2 1
E 1  0, 3  2x 9,81  3, 62 m ; E 2  1, 69  8
7,93 2x 9,81  1, 74 m

H1-2 = 3,62 - 1,74 = 1,88 m

Ejemplo 5.7. En un canal abierto de forma de V con lados de pendiente a 45º, ocurre un flujo
supercrítico de 2,8 m3/s, a una profundidad de 0,9 m. Calcúlese la profundidad apenas
corriente abajo de un salto hidráulico en este flujo.

Q2 Q2
 0,39  y 32
y1 y2 3 2,8 2 2,8 2
y 21 3  y 21 g
 y 22 3  y 22 g 0,9 2 x9,81  3  y 22 x9,81  y 2  1, 32m
,

Sección Rectangular

Si el canal es rectangular, sustituyendo en la Es. 5.15 y Q permanente

 b2y 1 
2 V 21 b 2 y 21 by 22 V 21 b 2 y 21
y1 Q2 y2 Q2
by 1 2  by 1 g
 by 2 2  by 2 g by 1 g
 2  by 2 g

multiplicando ambos miembros por dos (2) y recordando que F  V/ gy , o sea F2 =V2/gy,
la ecuación se transforma en:
 111

y3
y 21  2F 21 y 21  y 22  2F 21 y12  y 21 y 2  2F 21 y 21 y 2  y 32  2F 21 y 31 , dividiendo por y31
3 2
y2 y y  V2  q  1
 2F12 2   2   2F12 ; comoF12  1    y sustituyendo en la Ec. anterior
y1 y1  y1  gy1  y1  gy1

(5.16)

Como , factorizando tenemos

y2 1 1 2 4x 2xF21 1 18F21

Luego y1  2x 1  2

De donde
y2 1

y1 2

1  8F1  1
2
 (5.17)

De la ecuación 5.15 tenemos:

, que simplificada es

que multiplicada por y1/y22 y repitiéndolos pasos subsecuentes tenemos:

(5.18)

5.6.1.2. Pérdida de Energía.

La pérdida de energía en el salto es igual a la diferencia entre las energías específicas
antes y después del salto. Se puede demostrar analíticamente que dicha pérdida es igual a:

(y 2 y 1 )
E  E 1  E 2  4y 1 y 2 (5.19)

La relación /E1 se conoce como pérdida de energía relativa

 112

5.6.1.3. Longitud del salto

La longitud del salto, ha recibido gran atención de los investigadores pero hasta ahora
no se ha desarrollado un procedimiento satisfactorio para su cálculo. Sin duda, esto se debe al
hecho de que el problema no ha sido analizado teóricamente, así como a las complicaciones
prácticas derivadas de la inestabilidad general del fenómeno y la dificultad en definir las
secciones de inicio y fin del salto.

Se acepta comúnmente que la longitud LR del salto se defina como la distancia medida
entre la sección de inicio y la sección inmediatamente aguas abajo en que termina la zona
turbulenta (Fig 5.11).

Según el U.S. Bureau of Reclamation (1955), la longitud del salto hidráulico en canales
rectangulares horizontal, varía de acuerdo con la tabla 5.3, o bien con la curva So= 0 en la fig.
5.21.

Tabla 5.3. Longitud del resalto en canales rectangulares.

F 1  V 1 / gy 1 1,7 2 2,5 3 3,5 4 5 6 8 10
L R /y 2 4 4,35 4,85 5,28 5,55 5,8 6 6,1 6,12 6,12

La longitud de un resalto en un canal trapecial, es algo mayor debido a la asimetría

que se produce por efecto de la distribución no uniforme de la velocidad, la LR se puede
calcular como LR= A (y2- y1); donde A depende del talud del canal. En la tabla 5.4 se pueden
obtener distintos valores de A.

Tabla 5.4. Coeficientes A para LR en canales trapezoidales

Talud 0 0,5 0,75 1 1,25 1,5
A 5 7,9 9,2 10,6 12,6 15

Otros autores recomiendan estimar la longitud del resalto como LR=6(y2-y1); en

ambos casos y1 y y2 representan las profundidades alternas antes y después del salto
respectivamente.

5.6.1.4. Altura del Salto.

La diferencia entre las profundidades antes y después del resalto, es la altura del salto,
ó hj= y2 - y1. Al expresar cada término como la relación con respecto a la energía específica
inicial,
 113

(5.20)

donde hJ /E1, es la altura relativa, y1/E1, es la profundidad inicial relativa y y2/E1 es la

profundidad secuente relativa.

5.6.1.5. Eficiencia.
La relación entre la energía específica después y antes del resalto, se define como la
eficiencia del resalto, entonces E2/E1, es la eficiencia del mismo.

Como la pérdida relativa, la eficiencia, la altura relativa y las profundidades inicial

relativa y secuente relativa de un resalto hidráulico en un canal rectangular horizontal, son
funciones de F1, pueden graficarse contra F1 , dando como resultado un conjunto de curvas
características (figura 5.12). Con referencia a estas curvas nótense los siguientes aspectos:

Figura 5.12. Curvas Características de Resaltes Hidráulicos en Canales Rectangulares

Horizontales (Chow, 1994)

1.- La máxima altura relativa hj/E1, es 0.507, la cual ocurre cuando F1=2.77.
2.- La máxima profundidad relativa y2/E1, es 0.8, la cual ocurre en y1/E1=0.4 y F1= 1.73. Los
experimentos han demostrado que la transición desde un resalto ondulante a uno directo
ocurre aproximadamente en el punto F1= 1.73.
3.- Cuando F1= 1, el flujo es crítico y y1= y2=2/3 E1.
4.- Cuando F1 aumenta, los cambios en todas las relaciones características se vuelven
graduales.
 114

5.6.2. Localización del Resalto.

El resalto hidráulico ocurre en un flujo supercrítico, cuando su profundidad cambia
abruptamente a su profundidad secuente. En teoría, el resalto ocurrirá en un canal horizontal
rectangular, si las profundidades inicial y secuente y el número de Froude de aproximación
satisfacen la ecuación (5.17). Esta condición teórica por lo general se utiliza para localizar la
posición de un resalto. Para un cálculo más exacto de la posición del resalto, sin embargo,
debe considerarse la longitud del mismo. A continuación se ilustrará la localización de un
resalto hidráulico en tres casos comunes (Figura 5.13):

Figura 5.13. Localización de un Resalto Hidráulico

El caso A, muestra el resalto por debajo de una compuerta deslizante en un canal de

pendiente suave. Los perfiles AB Y CD pueden identificarse con facilidad como del tipo M3
 115

y M2, respectivamente. Los métodos de cálculo de estos perfiles se estudiaron en este capítulo.
La curva A'B, es una gráfica de la profundidad secuente con respecto a AB. Mediante la
posición F', puede estimarse la longitud del resalto. Mediante ensayo y error, puede
encontrarse una intersección horizontal entre las curvas A'B y CD, igual a la longitud del
resalto. Por ejemplo, la distancia horizontal EF es igual a dicha longitud, correspondiente a la
profundidad y2 en F. Es claro que el resalto se formará entre G y F, debido a que la profundidad
en F es secuente a la profundidad en G y a que la distancia EF mide la longitud del resalto.
Nótese que, si la longitud del resalto no se hubiera considerado en el análisis, se consideraría
que el resalto se forma con el punto F' aguas arriba, dando como resultado un error
representado por F'F. En el caso A, puede verse que al incrementar la profundidad del agua
hacia aguas abajo o al subir la curva CD, el resalto puede moverse hacia aguas arriba. La
profundidad de aguas abajo puede subirse hasta una altura para la cual el resalto
eventualmente se ahogue al frente de la compuerta deslizante. Al bajar la profundidad de aguas
abajo o al disminuir CD, el resalto se moverá hacia aguas abajo.

El análisis anterior, se aplica también a la localización de un resalto formado a la

salida de un vertedero de cresta delgada o de uno de rebose.

Cuando existe un resalto hidráulico por debajo de una compuerta deslizante, el flujo
que sale de la compuerta formará un chorro que posee una vena contracta. La distancia Le
desde la vena contracta hasta la apertura de la compuerta deslizante a menudo es corta. Con
respecto a esta distancia, existe una regla general adoptada, la cual establece que la vena
contracta se localiza aproximadamente a una distancia h desde la apertura de la compuerta
deslizante.

El caso B, muestra el resalto en un canal que tiene un quiebre en su pendiente de fondo

que cambia de empinada a suave. Por simplicidad, se supone que el flujo es uniforme a lo
largo del canal con excepción del tramo entre el resalto y el quiebre.

El resalto puede ocurrir en el canal empinado o en el suave, según si la profundidad

aguas abajo y2, es menor o mayor que la profundidad y1', secuente a la profundidad de aguas
arriba y1. Si la profundidad y2 , es mayor que y1' , el resalto ocurrirá en región empinada.
Luego la superficie curva OC es del tipo S1. La línea A'P, indica la profundidad secuente a
la línea AR. Ahora, determine una intersección horizontal IJ entre A'P y CO, el cual es igual
a la longitud del resalto. Es claro que un resalto HJ empezará en la sección que contiene a I.
Si la profundidad y2 se baja aproximadamente a algo menor que y1', el resalto empezará a
 116

moverse dentro del canal suave. En este caso el resalto se localiza como se describió en el
caso A.

El caso C muestra el resalto por detrás de una barrera de rebose. En teoría se formará
un resalto si la profundidad en la barrera es mayor que la profundidad secuente y1',
correspondiente a la profundidad supercrítica de aproximación y1. La localización del resalto
es igual a la del caso B, si éste ocurre en la región empinada. Al aumentar la altura de la
barrera, el resalto se moverá hacia aguas arriba. Al disminuirla se moverá hacia aguas abajo.
Cuando la profundidad de la barrera es menor que la secuente y1', la barrera será sobrepasada
por una oleada estacionaria en la forma de un ascenso superficial ondular solitario, a la que
no seguirán ondulaciones adicionales.

Ejemplo 5.8. En un canal rectangular de 6 m de ancho, transporta 11m3/s y descarga en una

solera protectora de 6 m de ancho, de pendiente nula, la velocidad inicial del resalto es Vi= 6
m/s. ¿Hallar HJ, E y yc?

q Q
y1  V1  bV 1  11
6x6  0, 306 m

y1 8.V21 0,306
y2  1 1  1  1  1, 353 m
8x6 2
2 gy 1 2 9,81x 0,306

Hj= 1,353-0,306= 1,047 m

62 2 1
E 1  0, 306  2xg = 2,14m E 2  1, 353 
11
6x1,353 2g  1, 44 m E  0, 69m

Ejemplo 5.9. Al pie de la caída de un canal rectangular, la velocidad del agua es de 6,4 m/s,
el espesor del chorro es de 0,42 m. ¿Calcular la profundidad conjugada y la longitud del resalto
utilizando la tabla 5.3 o fig. 5.15?

y2 
0,42
2
1
8x 6,4 2
 1 1, 67 m; F 1  3, 15 L R  5, 45 x 1, 67  9, 1m
9,
0 81x 0,42
 117

5.6.3. Tipos de Resalto y Cuencos Amortiguadores.

Los resaltos hidráulicos en fondos horizontales se clasifican en varias clases. De
acuerdo con los estudios del U.S. Bureau of Reclamation (1955), éstos pueden clasificarse
convenientemente según el número de Froude F1 del flujo entrante (figura 5.14) Como sigue:

Para F1 = 1 a 1,7, la superficie del agua muestra ondulaciones, no son necesarios cuencos
amortiguadores y se presenta el resalto ondulante.

Para F1 = 1,7 a 2,5, se desarrolla una serie de remolino sobre la superficie del resalto, pero la
superficie de aguas hacia aguas abajo permanece uniforme. La velocidad a través de la sección
es razonablemente uniforme y la pérdida de energía es baja. No se necesita cuenco
amortiguador y se presenta el resalto débil.

Figura 5.14. Diferentes Tipos de Resaltos Hidráulicos.

Para F1 = 2,5 a 4,5, existe un chorro oscilante o de transición que entra desde el fondo
del resalto hasta la superficie y se devuelve sin ninguna periodicidad. Cada oscilación produce
una onda grande con período irregular, muy común en canales, que puede viajar a lo largo de
varias millas causando daños ilimitados a bancas de tierra y a enrocados de protección.
Cuando el número de Froude está en este intervalo la figura 5.15, representa el cuenco más
eficiente que debe colocarse, el resalto que se presenta en este intervalo se llama resalto
oscilante.
 118

Figura 5.15. Características del Cuenco Amortiguador para Números de Froude entre 2,5 y 4,5

Para F1 = 4,5 a 9, la extremidad de aguas abajo el remolino superficial y el punto sobre

el cual el chorro de alta velocidad tiende a dejar el flujo, ocurren prácticamente en la misma
sección vertical. La acción y la posición de este resalto son menos sensibles a la variación el
la profundidad de aguas abajo. El resalto se encuentra bien balanceado y su comportamiento
es el mejor. La disipación de energía varía de 45% a 70% de E1. Se presenta entonces el resalto
estable.
 119

Para F1 =9 y mayores, el chorro de alta velocidad choca con paquetes de agua intermitentes
que corren hacia abajo a lo largo de la cara frontal del resalto, generando ondas hacia aguas
abajo, y puede prevalecer una superficie rugosa. La acción del resalto es brusca pero efectiva,
debido a que la disipación de energía puede alcanzar un 85% de E1 . Se produce entonces el
resalto fuerte.

En las figuras 5.16 y 5.17, se presentan las características del cuenco amortiguador
para números de Froude superior a 4,5.

Figura 5.16. Características del Cuenco Amortiguador para Números de Froude

Superior a 4,5 y Velocidad de llegada inferior a 15 m/s.
 120

Figura 5.17. Características del cuenco amortiguador para números de

Froude mayores de 4,5.

Nótese que los rangos del número de Froude dados arriba para los diferentes tipos de
resalto no están claramente demarcados sino que se traslapan en cierto modo según las
condiciones locales.
La tabla 5.5 (Azevedo; N y Acosta; G. 1976) nos da los calados o profundidades al final del
resalto para canales rectangulares, triangulares y trapezoidales y representa a la Ecuación 5.15;
en donde D1 = y1 , H1 = E1 y Z=m.
 121

Tabla 5.5 Calado más alto de igual energía y resalto hidráulico para rectangulares,
triangulares y trapeciales.
D1 = profundidad antes del resalto; z = pendiente de los costados del canal, horizontal a vertical;
H1 = carga de energía antes del resalto; b= ancho del fondo del canal. Para cada valor de D 1/H1
hay dos hileras de cifras; las cifras de la hilera superior multiplicadas por H1 dan los calados más
altos de igual energía; las cifras de la hilera inferior multiplicadas por H1 dan las profundidades
después del resalto. La primera columna es para canales rectangulares; la última, para canales
triangulares; el resto, para los trapeciales.
D1/H1 zH1/b
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1 10 
0,05 0,99 0,99
0,41 0,41 0,40 0,39 0,39 0,38 0,38 0,37 0,37 0,36 0,36 0,31 0,24

0,1 0,99 0,99 0,99 0,99 0,99

0,56 0,55 0,55 0,54 0,53 0,53 0,52 0,52 0,51 0,51 0,5 0,42 0,37

0,15 0,98 0,98 0,99 0,99 0,99 0,99 0,99 0,99 0,99 0,99 0,99
0,65 0,64 0,64 0,63 0,62 0,62 0,61 0,61 0,61 0,60 0,60 0,52 0,48

0,2 0,97 0,97 0,98 0,98 0,98 0,98 0,99 0,99 0,99 0,98 0,99
0,71 0,7 0,7 0,69 0,69 0,69 0,68 0,68 0,68 0,67 0,67 0,60 0,56

0,25 0,95 0,95 0,96 0,97 0,97 0,97 0,98 0,98 0,98 0,98 0,99
0,75 0,75 0,74 0,74 0,74 0,73 0,73 0,73 0,73 0,73 0,72 0,66 0,63

0,3 0,92 0,93 0,94 0,95 0,95 0,96 0,97 0,97 0,97 0,97 0,97 0,99
0,78 0,78 0,78 0,78 0,77 0,77 0,77 0,77 0,77 0,77 0,77 0,72 0,69

0,35 0,90 0,91 0,92 0,93 0,94 0,94 0,95 0,95 0,95 0,95 0,96 0,98 0,99
0,79 0,80 0,80 0,80 0,80 0,80 0,80 0,80 0,80 0,80 0,80 0,76 0,75

0,4 0,87 0,89 0,90 0,91 0,92 0,92 0,93 0,93 0,94 0,94 0,95 0,98 0,99
0,80 0,80 0,81 0,81 0,81 0,82 0,82 0,82 0,82 0,82 0,82 0,80 0,79

0,45 0,84 0,86 0,87 0,88 0,89 0,90 0,91 0,91 0,92 0,92 0,93 0,97 0,98
0,80 0,80 0,81 0,82 0,82 0,83 0,83 0,83 0,83 0,83 0,83 0,84 0,83

0,5 0,81 0,83 0,85 0,86 0,87 0,87 0,88 0,89 0,90 0,90 0,91 0,96 0,97
0,78 0,79 0,80 0,80 0,81 0,82 0,83 0,83 0,83 0,84 0,84 0,86 0,85

0,55 0,77 0,79 0,81 0,82 0,83 0,85 0,86 0,87 0,87 0,88 0,88 0,94 0,95
0,76 0,77 0,78 0,79 0,80 0,81 0,82 0,83 0,83 0,83 0,84 0,87 0,87

0,6 0,73 0,75 0,77 0,79 0,80 0,82 0,82 0,83 0,84 0,85 0,85 0,92 0,93
0,73 0,74 0,76 0,77 0,79 0,80 0,81 0,81 0,82 0,82 0,83 0,87 0,88

0,65 0,69 0,71 0,73 0,74 0,76 0,78 0,79 0,80 0,81 0,81 0,82 0,89 0,91
0,69 0,71 0,73 0,74 0,76 0,77 0,78 0,79 0,80 0,80 0,81 0,87 0,88

0,7 ------ ------ ------ ------ 0,72 0,73 0,74 0,75 0,76 0,77 0,78 0,86 0,88
------ ------ ------ ------ 0,72 0,73 0,74 0,75 0,76 0,77 0,78 0,85 0,87

0,75 ------ ------ ------ ------ ------ ------ ------ ------ ------ ------ ------ 0,83 0,85
------ ------ ------ ------ ------ ------ ------ ------ ------ ------ ------ 0,83 0,84

D0 0,67 0,68 0,69 0,70 0,71 0,72 0,72 0,73 0,73 0,74 0,74 0,79 0,80

Ejemplo 5.10. Un canal rectangular de 15 m de ancho se inicia al pie de un cimacio que tiene
una altura de 4,27 m, dicho cimacio tiene la misma longitud de cresta que el ancho del canal
 122

y con una carga h = 2,43 m sobre la misma, deberá descargar un gasto Q = 112,5 m 3/s, n =
0,025. El régimen aguas abajo en flujo uniforme debe ser subcrítico. Determinar la pendiente
necesaria en el canal para que el salto hidráulico principie justo al pie de la caída, así como la

longitud L de la zona que debe protegerse utilizando la tabla 5.3 ó fig. 5.15. hf0-1 =

112,5 7,5 V 20 7,5 2 1/3

q  7, 5 Sxmm V 0   1, 119 ms ,  0, 064m; yc 
3
15 6,7 2g 9,81

= 1,79 m

7,5 V21 10,373

V1  0,723  10, 373 ms ; 2g  5, 49m; E 1  6, 213m; F1  9,81x0,723
 3, 879

0,723 7,5 V 22
y2  2 1  8 x 3, 897 2  1  3, 639m; V 2  3,639  2, 061 ms ; 2g  0, 217m

H = E1-E2 = 6,213-(3,639+0,217) = 2,357 m

A2 = 15x3,638 = 54,58 m2; P2 = 22,278m; Rh = 2,45m

F 1  3, 897  y2R  5, 75 L R  20, 92  21m (tabla 5.3 o figura 5.15)

L

 

Ejemplo 5.11 En el problema anterior determinar la posición del salto hidráulico si la

pendiente del canal vale S0 = 0,0009 (En el cuenco n = 0,015)
 123

(0 , 000 9 ) 1/2  y n  3 , 4 8m
15y n
1 12 , 5  1
0,025 ( 15y n ) 152y n

 0, 369 y 1 
2,155 3,48
F2   1  8 x 0, 369 2  1  0, 775m
V2
gy 2 9,81x3,48 2

V21 7,5 2 1 (7,5/0,775)

E1  y1  2g  0, 775  0,775 2g  5, 553 m ; F1  9,81/ 0,775
 3, 51

Tabla 5.3
 L R  5, 55 x 3, 48  19, 3 m

H= 5,553 - 3,717 = 1,836 m

y A P R ¯R V V2/2g ¯V E S E2-E1 ¯S S0- ¯ S L L

(m) (m2) (m) (m) (m) (m/s) (m) (m/s) (m) (m) (m) (m)
0,723 10,84 16,44 0,659 10,37 5,48 6,20 0,0421
0,78 11,625 16,55 0,70 0,68 9,67 4,77 10,02 5,54 0,0337 0,659 0,03775 0,0368 17,89 17,89

Ejemplo 5.12. Un canal rectangular de b = 2m, consiste en tres tramos de diferentes pendientes
que sirven de descarga a una compuerta de abertura 0,90 m. La finalización del canal es
abrupta. Q = 10 m3/s; n=0,016; S01 =S03 =0,004; S02 = 0,01. ¿Dibujar los tipos de perfiles que
se presentan?
 124

(2y ) 2y 2/3
En el primer y último tramo  10  0,016 22y (0, 004 ) 1/2  y n1  1, 725  y c

(2y ) 2y 2/3
En el segundo tramo  10  0,016 22y (0, 01 ) 1/2  y n2  1, 20  y c

Vamos a comprobar si el primer cambio de pendiente es un control:

y A P R ¯R V V2/2g ¯V E S E2-E1 ¯S S0- ¯ S L L

(m) (m2) (m) (m) (m) (m/s) (m) (m/s) (m) (m) (m) (m)
0,9 1,8 3,8 0,47 0,486 5,55 1,57 5,27 2,473 0,0213 0,1988 0,0186 0,0146 13,6 13,6
1 2 4 0,5 5 1,27 2,27 0,0161
1,1 2,2 4,2 0,52 0,511 4,54 1,05 4,77 2,15 0,0125 0,1211 0,0142 0,0102 11,83 25,43
1,2 2,4 4,4 0,54 0,53 4,16 0,88 4,35 2,08 0,0099 0,068 0,01119 0,00719 9,47 34,91
1,3 2,6 4,6 0,56 0,55 3,84 0,75 4 2,05 0,081 0,030 0,009 0,005 6,17 41,09
1,37 2,74 4,74 0,57 3,64 0,67 2,04 0,007
0,571 3,74 0,005 0,0075 0,0035 1,42 42,5
M3 Si hay control en el primer cambio de pendiente
M2
1,37 2,74 4,74 0,57 0,58 3,64 0,67 3,61 2,04 0,007 0,0012 0,0068 0,0028 0,42 0,42
1,40 2,8 4,8 0,58 3,57 0,65 2,05 0,0066
1,50 3 5 0,6 0,59 3,33 0,56 3,45 2,06 0,0056 0,0162 0,0061 0,0021 7,56 7,98
1,60 3,2 5,2 0,61 0,607 3,12 0,49 3,22 2,09 0,0047 0,0314 0,0051 0,0011 26,49 34,47
1,70 3,4 5,4 0,629 0,62 2,94 0,44 3,03 2,14 0,0041 0,0431 0,0044 0,0004 100,19 134,66
1,725 3,45 5,45 0,633 2,89 0,42 2,15 0,0039
0,63 2,91 0,0123 0,0040 0,000029 412,19 546,86
S2
1,37 2,74 4,74 0,57 0,571 3,649 0,678 3,747 2,04 0,007 0,00507 0,00757 0,00242 2,09 2,09
1,3 2,6 4,6 0,56 3,84 0,753 2,05 0,0081
1,2 2,4 4,4 0,54 0,55 4,16 0,88 4 2,08 0,0099 0,03 0,009 0,00099 30,96 33
M3
1,2 2,4 4,4 0,545 0,561 4,16 0,88 3,9 2,08 0,0099 0,0359 0,0084 0,0044 8,11 8,11
1,37 2,74 4,74 0,57 3,64 0,67 2,04 0,007
 125

Como la conjugada de y1= 0,9 m, es y2= 1,97 m > yn lo que implica que hay resalto
hidráulico y termina en y2= yn= 1,725, entonces:

y1 
y2
2
 1  8F  1  1,725
2
2

2 

 1  8 xV22
gy 2
 1,725  1  8 x 2,89 2
 1 
 2 

gx1,725

 1  1,06m

  

Tabla 5.3
 3, 75  L R  6, 46 m (No colocar cuenco)
V1 4,7169 LR
F1  gy 1  9,8 1x1,06
 1, 46           y2

Ejemplo 5.13. Un canal rectangular muy largo de 3 m de ancho, conduce 1,8m3/s, experimenta
un cambio de rugosidad de n=0,0115 a n= 0,03. So= 0,01. Dibujar el perfil del agua en las
proximidades del cambio de rugosidad. ¿Qué tipo de curvas de remanso se presentan y a qué
distancia del cambio de rugosidad ocurre el resalto?

3y (3y ) 2 /3
1, 8  0,03 (32y )
(0, 01 ) 1/2      y n 2  0, 39 2 m

Se busca la conjugada para yn1= 0,21 m, que es y2= 0,49 m, lo que implica un resalto
aguas abajo del cambio de rugosidad, buscaremos a qué distancia del cambio de rugosidad
termina el m3 y comienza el resalto hidráulico, para eso buscamos la conjugada de yn2=0,392
m que es y1=0,28 m.
 126

y A P R ¯
R V V2/2g V¯ E S E2-E1 S
¯ L L
(m) (m2) (m) (m) (m) (m/s) (m) (m/s) (m) (m) So - S¯ (m) (m)
0,21 0,63 3,42 0,18 2,85 0,41 0,62 0,070
0,21 2,50 0,112 0,045 0,035 3,20 3,20
0,28 0,84 3,56 0,24 2,14 0,23 0,51 0,028

2,14
F1    1, 29  2, 5 (No colocar cuenco amortiguador )
V1
gy 1 9,81x 0,28

Ejemplo 5.14 En un canal rectangular horizontal, existe una compuerta de fondo de pared
delgada vertical a= 0,5m, que descarga Q= 10 m3/s, n= 0,015, b= 2,12m. Calcular y1 para que
salgan los 10 m3/s.

Calculando yc=1,3138 m, la conjugada de y = 0,5m, es y = 2,77 m, vemos que si se corta el

perfil H2 con el perfil conjugado de H3 por lo cual la descarga es libre.

¯ ¯ _
y A P R V V2/2g E E2-E1 _ So - S L L
R V S
(m) (m2) (m) (m) (m) (m) (m) (m) S (m) (m)
(m) (m/s)
0,5 1,06 3,12 0,34 9,43 4,53 5,03 0,084
0,46 6,51 3,065 0,0266 0,0266 115,13 115,13
1,3138 2,78 4,74 0,58 3,59 0,65 1,97 0,005
0,67
2,64 0,946 0,0026 0,0026 356,94 356,94
2,77 5,87 7,66 0,76 1,70 0,147 2,91 0,0009

Si asumo

 10       Cd=0,6  Q  0, 6 x 2, 12 x 0, 5 x 2 x g x 5
y1 Fig. 4.13
a  5
0,5

Para mayor valor de y1 es Cd= 0,6 entonces:

 127

Q 2 1 2 1
y1    12, 61 m.
10
Cd ba 2g 0,6x2,12x0,5 2g

Ejemplo 5.15. Bajo una compuerta sale un q = 6,1 m3/sm, el canal donde ocurre la descarga
es horizontal y n = 0,015. El canal se extiende 600 m aguas abajo, en la vena contraída y =
0,6 m de profundidad. Calcular y dibujar el perfil superficial resultante, si se produce un
resalto, determinar su ubicación. Considerar el canal muy ancho y poco profundo.

Aguas abajo de la compuerta (y<yc), se producirá un perfil H3 y mediante la ecuación

5.11, con So=0 y R= y se puede obtener la longitud donde se alcanzará la altura crítica.
Resumiendo el cálculo, para y1= 0,6m tenemos:

Y2 (m) 0,70 0,80 0.90 1 1,10 1,20 1,30 1,4 1,50

X12 (m) 35,4 72,2 108 143 172,7 202,2 225,7 244,3 254,3

Como se observa, la profundidad crítica se alcanzaría antes de llegar a la caída (Aprox.

260 m). Por consiguiente, ésta se transforma en sección de control, con profundidad crítica y
desde allá hacia aguas arriba se extenderá un perfil H2 hasta una sección tal que tenga una
profundidad secuente a la del perfil H3, produciéndose un resalto. El perfil H2,
correspondiente, está dado a continuación para y2=1,56 m.

y1 (m) 1,7 1,8 1,9 2 2,10 2,20 2,30 2,40

y2 (m) 10,2 30,3 67 122,4 197 283 409,2 546,8
 128

La curva de profundidades secuentes a las del perfil H3, se obtiene aplicando la Ec.
5.17 del Resalto Hidráulico. Esta ecuación produce los siguientes resultados:

y1 (m) 0,6 0,7 0,80 0,90 1 1,10 1,20

y2 (m) 3,26 2,97 2,69 2,48 2,30 2,14 1,98

La representación gráfica de las curvas H2, H3 y de profundidades secuentes del resalto,

nos indica que el resalto se produce en la sección de intersección de las curvas H3 y de
profundidades secuentes del resalto. Esta sección está ubicada a 136 m de la vena contraida.

Ejemplo 5.16. En el canal rectangular mostrado b= 2 m, Q = 10 m3/s, n= 0,013, So= 0,01.

Hallar yn, yc, distancia donde comienza el resalto (usar tabla 5.3 ó fig. 5.15) y dibujar la
superficie libre.
 129

(2y n ) 2y n 2/3
10  0,013 22y n (0, 01 ) 1/2 y n  1, 025 m

y1 8q 2 1,025
y2  2 1 gy 31
1  2 1 8x25
9,81x1,25 3  1  1, 77 m

Y A P R ¯R V V2/2g ¯V E S E2-E1 ¯S S0- ¯ S L L

(m) (m2) (m) (m) (m) (m/s) (m) (m) (m) (m) (m) (m)
2,5 5 7 0,714 2 0,203 2,70 0,001
0,70 2,08 0,163 0,00117 0,00882 18,46 18,46
2,3 4,6 6,6 0,696 2,17 0,24 2,54 0,0012
2,1 4,2 6,2 0,677 0,687 2,38 0,288 2,27 2,38 0,0016 0,1519 0,00144 0,00855 17,76 36,22
1,9 3,8 5,8 0,65 0,66 2,63 0,35 2,506 2,25 0,002 0,1359 0,00182 0,00817 16,63 52,85
1,77 3,54 5,54 0,638 0,64 2,82 0,406 2,72 2,17 0,0024 0,08 0,00224 0,00775 9,83 62,69

V1 q L
F1     1, 53  y R2  3, 75 LR  6, 64 m
5
gy 1 y 1 gy 1 1,025 9,81x1,025
 

Distancia= LR+x= 6,64 + 62,69 = 69,33 m.

Ejemplo 5.17. En los canales con cambios de pendiente mostrados, las longitudes son muy
largas. Dibujar y nombrar a mano alzada los tipos de perfiles que se pueden desarrollar.
 130

5.6.4. El Control de los Resaltos.

La formación y situación de los resaltos puede ser controlada por medio de
sobreelevaciones o caídas en el fondo del canal, por la acción de rebordes de cresta aguda o
de cresta ancha, y por la localización de dientes en un pozo de amortiguación (punto 5.6.3. de
este capítulo).

El efecto de una sobreelevación h en el fondo del canal fue estudiado, analítica y

experimentalmente, por J.W. Foster y R.A. Skrinde, en 1949. Los resultados experimentales
indican que los valores de y3, según se indica en la Fig. 5.18 (Chow, 1994), son menores que
los obtenidos analíticamente por aplicación de la ecuación de cantidad de movimiento entre
las secciones 1 y 3 similar al hecho para un vertedero de cresta ancha (ejemplo 3.2 Chow,
1994), suponiendo que el empuje sobre el escalón se produce por efecto de la distribución
hidrostática de presiones que produce la profundidad y2. La razón de esta discrepancia puede
explicarse por el hecho de que el chorro de aproximación no se uniformiza completamente
después del resalto y produce un impacto sobre la grada, el cual genera una fuerza mayor que
la supuesta. En todos los experimentos el resalto se localizó en una longitud x = 5(h+y3). Los
resultados del trabajo de Forster y Skrinde se presentan en la Fig. 5.18 con sus respectivos
rangos de aplicación (yc< y3 < y2). El resalto siempre se producirá en el lugar que indica el
dibujo de la Fig. 5.18 siempre y cuando yc< y3<y2.

Fig. 5.18. Relaciones experimentales entre F, y3/y1 y h/y1 para

una elevación abrupta.

Ejemplo 5.18: En un canal de flujo supercrítico, la prof. y1= 1,5 m y F1= 6, se desea producir
un resalto hidráulico con una sobrelevación h = 3 m. Cuál es la prof. y2 después del resalto, y3
sobre el escalón y la longitud del resalto.
 131

y1 1,5
y2  2 1  8F 21  1  2 ( 1  8 x 36  1 )  12 m;
Fig. 5.18 y3
F1  6 y h
y1  3
1,5 2   y1  5, 2  y 3  7, 8 m  L R  5(3  7, 8 )  54 m

En una sección de un canal recto sin sobreelevación si y2>y3 el resalto se rueda aguas
abajo y una sobrelevación lo controlaría en un sitio fijo como indica la Fig. 5.18, en caso
contrario y2<y3 el resalto se rueda aguas arriba y una caída lo mantendría en un sitio fijo como
muestra la Figura 5.19.

Cuando el fondo del canal tiene una caída Z0, el cálculo de la relación y2/y1 del resalto
producido varía dependiendo de que el vértice principal entre o no hacia el canal de aguas
arriba. H. Rouse, B.V. Bhoota y E.Y. Hsu, en 1949 (Chow (1994) y Aguirre (1980)),
investigaron el efecto estabilizador de las caídas y encontraron que entre los límites señalados
por el caso a, y el caso b, en la figura 5.19 se presenta una región de transición caracterizada
por un resalto ondular, el cual es altamente inestable pudiendo transformarse en un resalto del
caso a ó del caso b.

Fig. 5.19. Características del Salto Hidráulico Producido por una Caída
 132

Ejemplo 5.19. Determinar la altura del escalón Z0 para que el resalto termine justo en la
caída (caso a) y para que el resalto comience justo en la caída (caso b) para que y2= 6 m, si en
el canal rectangular de aproximación las profundidades y1= 1 m y F1= 3.

y1
y2  1  8F 21  1  2 ( 1  8 x 9  1 )  3, 77 m
1
2

𝑦2 ∆𝑍0
con = 6 𝑦 𝐹1 = 3 ⇒ = 2 ⇒ ∆𝑍0 = 2𝑚
𝑦1 𝑦

El caso b no se desarrolla.

El control de resaltos por vertederos puede ser observado en el texto de Chow (1994)
y Aguirre (1990).

5.6.5. Resalto en Canales Inclinados:

Cuando el fondo del canal no es horizontal las figuras 5.20 y 5.21 nos permiten hacer
la corrección para hallar las profundidades conjugadas (y ó d) y la longitud del resalto
hidráulico (distancia horizontal desde y1 hasta y2). Dichas figuras son tomadas de Chow
(1994).

Figura 5.20. Relaciones experimentales entre F1 y y2/y1 o d2/d1 para resaltos en canales
inclinados.
 133

Figura 5.21. Longitud en términos de la profundidad secuente de resaltos en canales

inclinados.

5.7. Aporte de Caudal hacia un Canal.

Hasta ahora hemos considerado que el caudal que fluye por un canal es conocido
cuando se estudia el perfil longitudinal. En muchos casos, el caudal no se conoce previamente.
Uno de los problemas más frecuentes se presenta cuando desde un depósito se agua en relativo
reposo, de suficiente extensión, parte un canal de dimensiones y pendiente dadas. Si la
pendiente del canal es supercrítica, en la sección de entrada se producirá un flujo crítico, el
cual se prolongará hacia aguas abajo mediante un perfil S2. Evidentemente, para los efectos
del cálculo, será necesario hacer la hipótesis de flujo supercrítico hacia aguas abajo de la
entrada, antes de calcular el caudal suponiendo que el flujo es crítico en ella. Si la entrada está
debidamente construida y se supone que las pérdidas en ella se pueden despreciar, se procede
al cálculo de la profundidad crítica la cual será 2/3E0, para entradas de sección rectangular,
donde E0 es la energía específica en el depósito de agua por encima del nivel de la cresta. La
profundidad crítica permite el cálculo del caudal y éste a su vez el de la profundidad normal.
Si la hipótesis de flujo supercrítico es acertada, la profundidad normal debe corresponder a la
de un flujo supercrítico; en otro caso, será necesario proceder en forma diferente. En los casos
prácticos, los canales de pendiente pronunciada o supercrítica son cortos tal como en los
aliviaderos de embalses o en descargas de aguas de lluvia. Si el canal tiene mucha pendiente,
se desarrolla un flujo ultra rápido, el cual no se puede considerar como flujo permanente.
 134

Cuando la pendiente del canal produce un flujo uniforme de profundidad crítica, ésta se
prolonga hasta la misma entrada. Este tipo de Flujo es inestable y debe evitarse.

El caso de flujo subcrítico es mucho más común. La forma más simple de este
problema se presenta con la descarga desde un lago o embalse de nivel conocido hacia un
canal de pendiente uniforme y suficientemente largo como para que los efectos de aguas abajo
no perturben el flujo de profundidad normal. Dadas estas circunstancias, no pueden producirse
perfiles M1 o M2 que hagan variar la profundidad desde el valor de E0 dado hasta la
profundidad normal. El flujo adopta, por tanto, la profundidad normal apenas entra en el canal.
El caudal Q y la profundidad normal yn se hallan resolviendo simultáneamente las ecuaciones
de Manning y de energía específica. Si la entrada produce pérdidas que deban ser
consideradas, éstas se incorporarán, en términos de la energía cinética o de la profundidad
normal, a la ecuación de energía específica.

Cuando los controles de aguas abajo no permiten que se produzca la profundidad

normal en las cercanías de la entrada hacia el canal, el problema se puede resolver suponiendo
un determinado gasto. Partiendo del control de aguas abajo, se remonta el canal calculando su
perfil hasta llegar al embalse con ciertos valores de V e y para los cuales se calcula la energía
específica. En la fig. 5.22 se ilustran algunos de los casos descritos en el texto precedente.

Fig. 5.22. Flujo desde un lago o embalse hacia un canal

 135

Ejemplo 5.20. Un canal rectangular de 3,70 m de ancho toma agua de un embalse. La toma es
suave y redondeada. El nivel del agua sobre la cresta de entrada es de E 0 = 1,85 m. El canal
de concreto con n = 0,013; es recto y muy largo. La pendiente es S = 0,001. Calcular el caudal
y el tipo de perfil superficial en la entrada del canal si se supone que las pérdidas son
despreciables.

Supongamos que la pendiente es supercrítica. Como el canal es rectangular, la

profundidad crítica será yc= E0 (2/3)= 1,23 m.

V c  y c g  3, 48 m/s Q  3, 48 x (3, 7 x 1, 23 )  15, 84 m 3 /s

2/3
3,7 yn 3,7 yn
15  0,013 3,72 yn (0, 001 ) 1/2 y n  1, 78 m.

Como yn= 1,78 > yc= 1,23 m., la pendiente será subcrítica, para el gasto Q = 15,84
m3/s y la hipótesis supuesta es incorrecta. Entonces, la solución corresponderá a la del caso b,
en la figura 5.22. Resolveremos así, simultáneamente la ecuación de Manning y la de energía
específica. Lo haremos por tanteos, en forma tabulada hasta obtener una yn que origine un
valor de E = E0. En la tabla que sigue se encuentra que yn= 1,60 m.

En consecuencia, el nivel de embalse, sobre la cresta bajará hasta yn=1,60 m. y esta

profundidad se mantendrá en el canal, correspondiente al caso b de la figura 5.22

yn (m) Q (m3/s) A (m2) V (m/s) V2/2g (m) E(m)

(supuesta)
1,7 14,11 6,29 2,24 0,255 1,95
1,65 13,55 6,1 2,22 0,251 1,9
1,6 13 5,92 2,19 0,245 1,85
 136

BIBLIOGRAFÍA

Aguirre, Pe. (1980). Hidráulica de Canales. Universidad de los Andes. Mérida, Venezuela.

Aguirre, Pe. Flores, I. y Macagno, E. (1993). Mecánica de los Fluidos. Universidad de Los
Andes. Mérida, Venezuela.

Azevedo, N. Acosta, G. (1976). Manual de Hidráulica. Editorial Harla, México, D.F.

Bolinaga, J.J. (1979). Drenaje Urbano. M.A.R.N.R., I.N.O.S. Caracas, Venezuela.

Bolinaga, J.J. (1985). Mecánica Elemental de los Fluidos. Fundación Polar. Universidad
Católica Andrés Bello. Caracas, Venezuela.

Chow, W.T. (1994). Hidráulica de Canales Abiertos. McGraw-Hill. Santafé de Bogotá,

Colombia.

Domínguez., F.J. (1934). Hidráulica. Editorial Universitaria, S.A.; Santiago de Chile.

Foster, J.W. y Skrinde, R.A. (1949). "Control of Hydraulic Jump by Sills", Proceedings,
A.S.C.E.

French, R.H. (1988). Hidráulica de Canales Abiertos. McGraw-Hill. México, D.F.

Giles, Ranald. Evett, J. y Cheng, L. (1994). Mecánica de los Fluidos e Hidráulica. Serie
McGraw-Hill. Tercera Edición. Madrid. España.

Henderson, F.M. (1996). "Open Channel Flow", MacMillan Series in Civi Enginnering New
York.

INOS. (1975). Normas e Instructivos para el Proyecto de Alcantarillado. Caracas,

Venezuela.

Machado P, J.L. (1979). Fundamentos del Ciclo Hidrológico. Dpto. de Meterología e

Hidrología. U.C.V. Caracas, Venezuela.
 137

Maza, J.A. (1967). "Erosión del Cauce de un Río en el Cruce de un Puente", Ingeniería
Hidráulica de México, Vol XXI, Nº 1 y 2, México.

MOP. (19679. Manual de Drenajes. División de Vialidad. Caracas, Venezuela.

Rouse, H. (1950). "Hidráulica", Traducción, Editorial Dossat, Madrid, España.

Rouse, H. (1961)., "Fluid Mechanics for Hydraulics Engineers", Dover Publications, New
York.

Schlag, A. (1957). Hidrauliqué Genérale. Dunod. París.

Sotelo A., G. (1982). Hidráulica General. Volumen 1. Editorial Limusa, México.

Torres H., F. (1981). Obras Hidráulicas. Editorial Limusa. México.

U.S. Bureau of Reclamation. (1955)., "Research Sudies on Stilling basins, energy

dissipators, and associated appurtenances", Hydraulic Laboratory Report. Nº Hyd-
399.
 138

APENDICE 1
Programa en Basic para El Cálculo de Curvas de Remanso
5 Cls: Print “Curvas De Remanso”
10 L = 0
11 Rp = 0
12 Vp = 0
13 DE = 0
14 S = 0
15 W = 0
16 Ds = 0
17 Dx = 0
18 L = 0
20 INPUT “ Q = “; Q
30 INPUT “ B = “; B
40 INPUT “ M = “; M
50 INPUT “ N = “; N
60 INPUT “So = “; So
70 INPUT “Y = “; Y
80 A = (B.Y)+M.(Y)^2
90 P = B+(2*SQR(Y^2+(M*Y)^2))
100 R = A/P
105 X = (Q*N/(R^(2/3)*A))^2
110 V =Q/A
115 Z = V^2/(2*9.81)
120 Eo = ABS(Y+(V^2/(2*9.81)))
130 IF W = 0 THEN 205
140 RP = (RP+R)/2
150 VP = (VP+V)/2
160 DE = ABS (Eo-DE)
170 S = ABS ((VP*N/RP^(2/3))^2)
180 DS = ABS (S-So)
190 DX = DE/DS
200 L = L+DX
205 ’BEEP 1
210 PRINT “AREA = ”; A
220 PRINT “PERIMETRO = ”; P
230 PRINT “RADIO H = ”; R
240 PRINT “RADIO H.P = ”; RP
250 PRINT “VELOCIDAD = ”; V
255 PRINT “V^2/(2G) = ”; Z
260 PRINT “Veloc. Prom = ”; VP
270 PRINT “ENERGIA = ”; Eo
275 PRINT “PEND. (S1) = ”; X
280 PRINT “(E2 – E1) = ”; DE
290 PRINT “PROM. (S) = ”; S
300 PRINT “ (So-S) = ”; DS
310 PRINT “DELTA X = ”; DX
320 PRINT “DELTA X ACUMULADO = ”; L
330 W = W+1
340 RP = R
350 VP = V
360 DE = Eo
370 GOTO 70