DISEÑO Y CONSTRUCCIÓN DE CANALES 1

Comportamiento Hidraulico de Canales en

HEC-RAS, Matlab, Python y Excel
Juan Camilo Mojica Arias
[email protected]

Palabras clave: Flujo uniforme, Energía específica, Número B. Objetivo Especifico:

de Froude, Profundidad crítica, Tirante hidráulico, Flujo sub-
Realizar 2 conclusiones sobre el comportamiento hidraulico
crítico, Flujo supercrítico.
de cada sistema.

I. R ESUMEN :
III. M ETODOLOGIA
El comportamiento hidráulico de canales se refiere al análi-
sis del flujo de agua en canales abiertos, como ríos o canales En primer lugar, se plantea el ejercicio con el objetivo de
artificiales. Se estudian diferentes tipos de flujo, como el flujo verificar las condiciones del canal, así como las variables
uniforme, donde la velocidad y el caudal son constantes, y el conocidas y las solicitadas. El ejercicio se desarrolla
flujo no uniforme, en el que las características del flujo varían inicialmente de forma manual utilizando el editor de
a lo largo del canal. Los factores claves incluyen la pendiente ecuaciones de Lyx, con las ecuaciones que se enumeran a
del canal, la rugosidad del fondo (a menudo descrita con el continuación:
coeficiente de Manning), y el perfil de flujo, que determina si
el flujo es subcrítico (lento) o supercrítico (rápido), lo cual se
evalúa con el número de Froude. Ec. N°1 Energia Especifica (m): Ee(m) = Y c + V2
:
2g

Ec. N°2 Ar. Hidr. (m2): A = 12 ZaY c2 + bY c + 12 ZbY c2 :

Ec. N°3 Espejo (m): T = ZaY c + b + ZbY c :

√ √
Ec. N°4 P. Moj (m): P = Y c Za2 + 1 + b + Y c Zb2 + 1 :

0.5(ZaY c2 +2bY c+ZbY c2 )

Ec. N°5 Rad. Hidrau. (m): R = √ √
Y c( Za2 +1+b+Y c Zb2 +1
:

Q
Ec. N°6 Velocidad V = 0.5(ZaY c2 +2bY c+ZbY c2 ) :

Ec. N°7 Flujo Critico gA3 T −1 − Q2 = 0 :

Ec. N°8 Caudal Q = V A

Fig. 1. Canal La Albina en el barrio Gustavo Restrepo, kr 14 Bis B con calle Ec. N°9 NF N F = pV
gDh
31 A Sur, en Bogotá. q
2
Ec. N°10 Yc (Canal rectangular): Y c 3 bQ
2 ∗g

El análisis incluye la determinación de profundidades críti-

cas, curvas de energia especifica, y pérdidas de energía, que Posteriormente, se implementa en HEC-RAS, MATLAB,
son fundamentales para diseñar canales que sean eficientes PYTHON y EXCEL. Finalmente, se extraen dos conclusiones
y seguros. Además, los cambios en las condiciones de flujo, sobre el comportamiento hidráulico del sistema.
como transiciones de ancho o pendiente, deben gestionarse
para evitar erosión o ineficiencias en el transporte de agua.
IV. P ROCEDIMIENTO
II. O BJETIVOS A. Ejercicio 1:
A. Objetivo General: En un canal trapezoidal de ancho b= 0.3 m y talud z=
Diseñar, en Hec Ras, Matlab, Python y Excel el compor- 1, determinar el caudal que transporta para una energía
tamiento hidráulico para cada sistema propuesto por el Ing. específica mínima igual a 0.48 mkg/kg.
Edgar Ladino para el Taller #1 de Diseño y Construccion de
Canales. Solucion:
DISEÑO Y CONSTRUCCIÓN DE CANALES 2

Variables conocidas: b (Base), z (Talud), Ee (Energia las profundidades alternas. Graficar curva EE
Especifica)

Variables solicitadas: Q(Caudal) Solucion:

Partimos de la Ec. N°1
V2
Ee(m) = Y c + 2g

De la Ec. N°8, despejamos Velocidad y reemplazamos:

Q2
Ee(m) = Y c + A2 2g

Reemplazamos el Caudal por la Ec. N°7 Variables conocidas: Q (Caudal), Ee (Energia Especifica)
3 −1
gA T
Ee(m) = Y c + A2 2g

Simplificamos la Ecuacion:
A
Ee(m) = Y c + 2T

Reemplazamos A:(Ec. N°2) y T:(Ec. N°3)

1
2 ZaY c2 +bY c+ 12 ZbY c2
Ee(m) = Y c + 2(ZaY c+b+ZbY c) Variables solicitadas: Profundidades Alternas, Curva EE
Reemplazamos b, Ee, Za y Zb Partimos de la Ec. N°1
1
2 (1)Yc2 +(0.3m)Y c+ 12 (1)Y c2
0.48m = Y c + 2((1)Y c+0.3m+(1)Y c)

Simplificamos la Ecuacion: De la Ec. N°8, despejamos VelocidadVy2 reemplazamos:

2
Ee(m) = Y c + 2g
Y c +(0.3m)Y c
0.48m = Y c + 2(Y c+0.3m)

Utilizando el Solver de la Calculadora: Reemplazamos el Caudal por la Ec. N°7

2
Ee(m) = Y c + AQ2 2g
Y c =: 0.3627m
De la Ec. N°2, calculamos el A: Simplificamos la Ecuacion: gA3 T −1
Ee(m) = Y c + A2 2g
1 2 1 2
A= 2 (1)(0.3627m) + (0.3m)(0.3627m) + 2 (1)(0.3627m)

Resolviendo para A: Reemplazamos A y T: A

Ee(m) = Y c + 2T
2
A = 0.2404m
De la Ec. N°3, calculamos el T: Cancelamos b, y obtenemos Yc by
Ee(m) = Y c + 2b
T = (1)(0.3627m) + 0.3m + (1)(0.3627)
Resolviendo para T: Simplificamos la Ecuacion:
2
3 Ee(m) = Y c
T = 1.0256m
De la Ec. N°7, despejamos Q: Utilizando el Solver de2 la Calculadora:
p 3 (1.829m) = Y c
gA3 T −1 = Q :
Reemplazamos en valores g, A y T: De la Ec. N°7, despejamos Q:
√ Y c =: 1.219m
9.81 ∗ 0.24043 ∗ 1.0256−1 = Q :
Resolvemos para Q: Reemplazamos A, T, Q 3y g:
3
gA T −1 − Q2 = 0 :
Q= 0.3646 ms :
Conclusiones: Resolvemos para b: 3 −1 3

Se obtuvo que para un canal trapezoidal de ancho de 0.3 m y 9.81(by) b − (0.879 ms )2 = 0

talud de 1, el caudal que transporta para una energía específica
mínima igual a 0.48 mkg/kg es de 0.3646 m r 3
Resolvemos b: (0.879 ms )
3 =b
B. Ejercicio 2: (9.81 ms )∗(1.219m)3

Un canal rectangular transporta un caudal de 0.85+0.029

m3/s-m, si la energía específica es de 1.80+0.029m, calcular Con ayuda de excel, graficamos la Curva EE (Anexo)
b = 0.2085m
DISEÑO Y CONSTRUCCIÓN DE CANALES 3

Asi, se obtuvieron las siguientes profundidades alternas, y se

evidencia que el flujo critico es de 0.768:

Fig. 2. Curva de Energia Especifica

Fig. 4. Profundidades Alternas

C. Ejercicio 3:
Un canal rectangular transporta un flujo con una D. Ejercicio 4:
profundidad de 1.20 m, con un Número de Froude de
Un canal rectangular transporta un flujo con un Número
0.5-0.0051 B a) Cuál es la profundidad crítica para este flujo,
de Froude de 0.6 para el cual corresponde una profundidad
b) Cual es la profundidad alterna correspondiente.
crítica de 0.50+0.015 m. Calcular el caudal específico y la
profundidad de flujo.
Variables conocidas: Dh (Profundidad Hidraulica), NF
(Numero de Froude) Variables conocidas: Yc (Profundidad Critica), NF (Numero
de Froude)
Variables solicitadas: Profundidad Critica, Profundidad
Alterna Variables solicitadas: Caudal Especifico, Profundidad de
De la ecuacion N°9, obtenemos: Flujo

N F = pV Reemplazamos el Caudal por la Ec. N°7

gDh
gA3 T −1
Ahora, despejamos la velocidad Ee(m) = Y c + A2 2g
√
V = N F gDh De la Ec. N°7, despejamos el caudal:
Obtenemos V= gA3 T −1 = Q2 :
V = 1.698 m
s Reemplazamos el Caudal:
Ahora, suponiendo una base del canal de 1 metro, y con gA3 T −1 = V 2 ∗ A2 :
ayuda de excel, graficamos la curva de Ee:
Simplificamos la Ecuacion:
gA
T =V2
Reemplazamos A y T, teniendo en cuenta que es un canal
rectangular:
g∗b∗y
b =V2 :
Simplificamos y obtenemos:
√
V = gy
Simplificamos la Ecuacion:
V = (9.81 sm2 )(0.515m) = 2.2476 m
p
s

Ahora, suponiendo una base de canal de 1, obtenemos un

caudal especifico:
Q=V ∗A
Fig. 3. Curva de Energia Especifica Reemplazamos V y A, y obtenemos el Caudal Especifico:
DISEÑO Y CONSTRUCCIÓN DE CANALES 4

3
Q = (2.2476 m m
s )(0.515m)(1m) = 1.1575 s Igualando E1 y E2:
V 22
De la Ecuacion de Froude, despejamos la profundidad 0.8385m = 0.75789m + 2g
hidraulica:
Resolviendo para V2:
V2
Dh = N F ∗g V 22
0.8385m − 0.75789m = 2g
Obtenemos V=
2
Obtenemos:
(2.2476 m
s )
Dh = (0.515m)∗(9.81 m2 )
s V 2 = 1.257 m
s
Asi, obtenemos la profundidad de flujo: El área de la sección final es:
Dh = 0.858m A2 = Q
= 1.211
= 0.9638m2
V2 1.257

E. Ejercicio 5: La nueva profundidad se obtiene dividiendo el área por el

Un canal rectangular se transporta un caudal de 1.20 ancho:
m3/s m con una profundidad de 0.80 m El fondo del canal A2 0.9638
h2 = b = 1.87 = 0.5154m
experimenta un suave cambio de profundidad que produce un
decremento de la elevación absoluta de la superficie del agua Finalmente, la variación en el fondo del canal es:
de 0.0477 m Calcular la variación experimentada en el fondo ∆hf = h1 − h2 = 0.80559 − 0.5154 = 0.2902m
del canal:
Por lo tanto, la variación experimentada en el fondo del
canal es de 0.2902 m (aproximadamente 29 cm).
Variables conocidas: Q (Caudal), h1 (Profundidad),
Decremento en la elevación

Variables solicitadas: Variacion Fondo del Canal:

De la ecuacion N°10, despejamos la base:

q
c3 ∗g
b = YQ 2

Reemplazamos:
r 3
(0.80559)3 ∗(9.81 ms )
b= 3 :
(1.211 ms )2

Asi, obtenemos:
b = 1.87m
El área de la sección transversal con el nuevo ancho es:
A1 = b ∗ h1 = 1.87 ∗ 0.80559 = 1.5065m2
La velocidad inicial v1 es:
Q 1.496
v1 = A1 = 1.20 ≈ 0.8037m/s
La energía específica en la sección inicial es:
V 12
E1 = h1 + 2g

Sustituyendo los valores:

(0.8037)2
E1 = 0.80559 + 2∗9.81

Obtenemos:
E1 = 0.8385m
Con una disminución de la elevación de 0.0477 m, la nueva
profundidad de la superficie es
h2 = 0.80559m − 0.0477m = 0.75789m
La energía en esta nueva sección es:
V 22
E2 = h2 + 2g