MODELACIN HIDRULICA DE UN CANAL UTILIZANDO HEC-RAS Y

COMPARACIN ENTRE LOS MODELOS INTERPOLADOS Y SIN INTERPOLAR

Alexander Garcs Rojas

Emmanuel Gmez Restrepo

Hidrulica de Canales.
Ingeniera Civil Facultad de Ingeniera
Medelln, 02 de junio de 2017

Introduccin
Siendo los fenmenos hidrulicos unos de los ms impredecibles y variados dentro de la
naturaleza, el tratamiento de estos por medio de mtodos matemticos suele verse limitado
ante la complejidad de los diferentes flujos que se pueden presentar; en la actualidad
disponemos de numerosas tcnicas y herramientas que facilitan esta labor.
A continuacin se presenta un estudio somero de un canal natural mediante la modelacin
digital de este usando el software HEC-RAS (Hydrologic Engineering Center's River
Analysis System) desarrollado por el Cuerpo de Ingenieros de la Armada de Estados
Unidos. Se muestran los principales parmetros hidrulicos en diferentes secciones de un
tramo de 500 metros del canal para los caudales asociados a cuatro periodos de retorno
diferentes, luego se compara analticamente las diferencias de los resultados entre el
modelo para las secciones originales y los resultados para el modelo con secciones
intermedias interpoladas generadas automticamente por el software.

Consideraciones de rugosidad para la modelacin: Para la modelacin del canal se

tuvieron en cuenta los siguientes valores de rugosidades:
Margen izquierda Centro del canal Margen derecha
n Maning 0.013 0.028 0.03
Tabla 1. Rugosidades de Maning del canal.

Margen izquierda: definida por un muro en concreto, para lo cual la rugosidad ha sido
determinada segn valores recomendados en la literatura (Chow, 1994).
 MODELACIN HIDRULICA HEC-RAS
HIDRULICA DE CANALES INGENIERA CIVIL

Margen derecha: definida por una llanura de inundacin ocupada con potreros, la cual
tambin ha sido determinada su rugosidad con base a valores recomendados en la literatura
(Chow, 1994).
Centro del canal: configurado por un lecho en grava con un dimetro promedio de 40 mm.
La rugosidad para esta parte del canal se determina con base a la ecuacin de Strickler
(Duarte, 2009).

Modelo sin interpolar

La siguiente imagen muestra la ubicacin de las secciones analizadas a lo largo del eje, al
igual que la variacin en planta del canal:

Imagen 1. Eje del canal en planta y ubicacin de las secciones analizadas.

La flecha indica el sentido de direccin del flujo en el canal. La seccin 1 (k0+1000)

corresponde a la ltima seccin aguas abajo del canal en el tramo analizado, y la 12
(k0+500) a la primera estacin aguas arriba.
Las siguientes tablas presentan la variacin de algunos parmetros hidrulicos para cada
caudal en su respectivo perodo de retorno:
Estaci
v (m/s) A (m2) T (m) F Y (m) Yc (m) E (m) (Pa) Pm (m) Rh (m)
n
k0+500 2.53 6.10 6.51 0.78 1.45 1.31 1.77 49.00 7.36 0.83
k0+550 1.04 22.62 33.62 0.24 2.07 2.12 6.69 35.84 0.63
k0+600 0.81 30.77 40.98 0.15 3.10 3.12 3.53 43.89 0.70
k0+650 0.42 54.34 34.18 0.08 3.27 3.27 0.95 37.82 1.44

2
 MODELACIN HIDRULICA HEC-RAS
HIDRULICA DE CANALES INGENIERA CIVIL

Estaci
v (m/s) A (m2) T (m) F Y (m) Yc (m) E (m) (Pa) Pm (m) Rh (m)
n
k0+700 3.05 5.14 5.78 0.99 1.02 1.02 1.48 72.36 6.95 0.74
k0+750 6.81 2.23 5.82 3.43 0.52 1.02 2.86 483.65 6.19 0.36
k0+800 3.03 5.02 5.52 0.97 1.25 1.25 1.71 72.41 6.48 0.77
k0+850 3.83 3.94 4.86 1.32 1.43 1.61 2.17 124.19 5.80 0.68
k0+900 0.63 22.30 12.50 0.14 2.64 1.14 2.67 2.43 14.20 1.57
k0+907 0.70 22.26 12.49 0.14 2.64 2.67 2.87 14.19 1.57
k0+950 1.47 10.48 7.83 0.39 1.53 1.64 14.66 9.83 1.07
k0+1000 2.37 5.94 5.71 0.66 1.63 1.48 2.01 42.21 7.43 0.80
Tabla 2. Parmetros hidrulicos para el perodo de retorno de 2.33 aos y caudal de 15 m 3/s.
V: Velocidad del agua.
A: rea mojada.
T: Ancho superficial.
F: Nmero de Froude.
Y: Profundidad o altura de la lmina de agua.
E: Energa especfica.
: Esfuerzo cortante del fondo del canal.
Pm: Permetro mojado.
Rh: Radio hidrulico.

Se observa que en la mayora de las estaciones a lo largo del canal se presenta un flujo
subcrtico con base en el nmero de froude, a excepcin de la estacin 5 en la abscisa
k0+850. Se observa tambin que las velocidades son mayores para las secciones de canal
con rea menor, guardando relacin a la ecuacin de continuidad. A continuacin se
presenta una imagen en la cual se observa la estacin con flujo supercrtico:

Imagen 2. Estacin 5 (k0+850) para un perodo de retorno de 2.33 aos y caudal de 15 m 3/s.

La lnea roja corresponde a la lnea de profundidades crticas para dicha seccin, mostrando
as un flujo supercrtico. La lnea verde es la lnea de energa.

3
 MODELACIN HIDRULICA HEC-RAS
HIDRULICA DE CANALES INGENIERA CIVIL

La siguiente tabla presenta los parmetros hidrulicos para el perodo de retorno de 25

aos:
Estaci
v (m/s) A (m2) T (m) F Y (m) Yc (m) E (m) (Pa) Pm (m) Rh (m)
n
k0+500 2.66 7.41 6.87 0.76 1.65 1.48 2.00 51.27 7.88 0.94
k0+550 0.98 30.44 33.80 0.21 2.31 2.34 5.74 36.36 0.84
k0+600 0.76 40.24 40.98 0.14 3.33 3.35 3.04 44.36 0.91
k0+650 0.46 62.59 38.60 0.08 3.49 3.50 1.12 42.49 1.47
k0+700 3.03 6.76 8.89 0.88 1.27 1.27 1.71 65.92 10.39 0.65
k0+750 7.24 2.67 5.89 3.35 0.59 1.16 3.23 516.92 6.36 0.42
k0+800 3.36 5.75 5.64 1.01 1.38 1.40 1.95 85.47 6.76 0.85
k0+850 4.06 4.75 5.07 1.28 1.59 1.79 2.42 131.47 6.22 0.76
k0+900 0.69 28.51 37.42 0.14 2.90 1.28 2.93 2.76 39.39 0.72
k0+907 0.83 28.11 37.42 0.16 2.89 2.93 3.82 39.36 0.71
k0+950 1.56 12.70 12.05 0.38 1.78 1.90 15.80 14.51 0.88
k0+1000 2.58 7.06 8.01 0.67 1.81 1.68 2.27 47.98 9.94 0.71
Tabla 3. Parmetros hidrulicos para el perodo de retorno de 25 aos y caudal de 19 m 3/s.

En este caso se presenta flujo supercrtico en 2 secciones del canal (k0+800 y k0+850,
estaciones 6 y 5 respectivamente). La siguiente imagen muestra la seccin 6:

Imagen 3. Estacin 6 (k0+800) para un perodo de retorno de 25 aos y caudal de 19 m 3/s.

Ac tambin se observa el flujo supercrtico ya que la lmina de agua est por debajo de la
lnea de profundidades crticas.
La siguiente tabla corresponde al perodo de retorno de 50 aos, para un caudal de 25 m3/s:
Estaci
v (m/s) A (m2) T (m) F Y (m) Yc (m) E (m) (Pa) Pm (m) Rh (m)
n
k0+500 3.39 7.67 6.94 0.95 1.69 1.69 2.26 82.28 7.99 0.96
k0+550 1.09 35.69 35.14 0.23 2.46 1.48 2.50 6.92 37.95 0.94
k0+600 0.84 46.55 40.98 0.15 3.48 3.51 3.68 44.67 1.04

4
 MODELACIN HIDRULICA HEC-RAS
HIDRULICA DE CANALES INGENIERA CIVIL

Estaci
v (m/s) A (m2) T (m) F Y (m) Yc (m) E (m) (Pa) Pm (m) Rh (m)
n
k0+650 0.55 68.89 42.31 0.09 3.65 3.66 1.55 46.37 1.49
k0+700 2.37 13.50 37.66 0.60 1.65 1.65 1.89 36.98 39.67 0.34
k0+750 7.85 3.24 5.98 3.31 0.69 1.34 3.79 571.86 6.57 0.49
k0+800 3.88 6.69 7.22 1.10 1.53 1.93 2.28 109.24 8.46 0.79
k0+850 4.33 5.93 5.37 1.24 1.82 2.03 2.76 139.90 6.81 0.87
k0+900 0.72 41.19 52.24 0.14 3.17 1.45 3.20 2.95 54.92 0.75
k0+907 0.83 40.91 52.18 0.16 3.17 3.20 3.75 54.85 0.75
k0+950 1.57 19.65 38.02 0.36 2.06 2.17 15.08 41.02 0.48
k0+1000 2.73 9.58 11.74 0.66 2.06 2.06 2.55 50.91 13.89 0.69
3
Tabla 4. Parmetros hidrulicos para el perodo de retorno de 50 aos y caudal de 25 m /s.

Para este caudal se presenta flujo supercrtico en las mismas estaciones que para el perodo
de retorno de 25 aos. La siguiente imagen muestra la seccin 9 ubicada en la abscisa
k0+650:

Imagen 4. Estacin 6 (k0+650) para un perodo de retorno de 50 aos y caudal de 25 m 3/s.

Esta seccin presenta la mayor altura de la lmina de agua en la zona de inundacin para
este caudal. Tambin se determina que ocurre un flujo subcrtico (ver Tabla 4). Esta seccin
tambin corresponde a la de mayor energa especfica, la cual fue determinada como la
distancia entre el fondo del canal y la lnea de energa. La lnea de energa se encuentra
ubicada justo en la superficie del agua. Esto ocurre ya que la velocidad es la ms baja en
esta seccin, obtenindose una energa de velocidad muy baja.
En la siguiente tabla se presentan los parmetros hidrulicos para el perodo de 100 aos y
caudal de 33 m3/s:
Estaci
v (m/s) A (m2) T (m) F Y (m) Yc (m) E (m) (Pa) Pm (m) Rh (m)
n
k0+500 3.67 9.43 7.39 0.94 1.93 1.93 2.59 90.94 8.65 1.09
k0+550 1.32 39.24 37.72 0.27 2.56 1.97 2.61 9.96 40.71 0.96
k0+600 1.01 50.53 40.98 0.17 3.58 3.61 5.22 44.87 1.13
k0+650 0.68 73.15 42.84 0.11 3.75 3.76 2.32 47.01 1.56
k0+700 2.13 22.92 60.00 0.51 1.83 1.83 2.00 28.80 62.37 0.37
k0+750 7.92 4.28 6.14 2.93 0.86 1.91 3.99 532.83 6.95 0.62

5
 MODELACIN HIDRULICA HEC-RAS
HIDRULICA DE CANALES INGENIERA CIVIL

Estaci
v (m/s) A (m2) T (m) F Y (m) Yc (m) E (m) (Pa) Pm (m) Rh (m)
n
k0+800 4.27 8.32 14.07 1.13 1.72 2.01 2.59 126.87 15.44 0.54
k0+850 4.44 7.92 7.74 1.13 2.16 2.68 3.14 135.94 9.49 0.83
k0+900 0.92 43.45 52.73 0.18 3.22 1.64 3.26 4.71 55.53 0.78
k0+907 1.05 43.02 52.63 0.19 3.21 3.26 5.99 55.41 0.78
k0+950 2.35 15.95 26.66 0.55 1.95 1.41 2.20 34.44 29.44 0.54
k0+1000 2.30 20.34 50.67 0.51 2.38 2.38 2.58 34.34 53.29 0.38
Tabla 5. Parmetros hidrulicos para el perodo de retorno de 100 aos y caudal de 33 m3/s.

Este caudal presenta las situaciones hidrulicas menos favorables ya que ocurren las
mayores inundaciones, por lo cual la mayora de los anlisis se concentran ac.
Al igual que para los dos perodos de retorno anteriores, en este caudal se presenta flujo
supercrtico en las mismas estaciones ya mencionadas anteriormente (k0+800 y k0+850).
La velocidad mnima se da para el nmero de Froude ms bajo en la seccin 9 (k0+650)
presentando flujo subcrtico, correspondiente a la profundidad de circulacin mayor con
3,75 metros (ver Tabla 5). En la siguiente imagen se muestra una seccin con flujo
subcrtico:

Imagen 5. Estacin 11 (k0+550) para un perodo de retorno de 100 aos y caudal de 33 m3/s.

La lnea de profundidades crticas es la lnea punteada que se encuentra dentro del agua. En
esta imagen se observa que la lmina de agua se encuentra por arriba de la lnea de
profundidades crticas que arroja el programa, deducindose as la presencia de un flujo
subcrtico. La lnea de energa se observa en color verde por encima de la lmina de agua.
Como es de esperarse, todos los parmetros hidrulicos, a excepcin del esfuerzo cortante,
permetro mojado y radio hidrulico, son mayores para los caudales y perodos de retorno
mayores ya que estos son funciones directas del caudal (ver Tablas 2, 3, 4 y 5).
Anlisis esfuerzo cortante: El esfuerzo cortante es una medida de la fuerza de friccin del
agua en movimiento contra la superficie del lecho por el cual circula, siendo este
directamente proporcional al radio hidrulico y la pendiente hidrulica.

6
MODELACIN HIDRULICA HEC-RAS
HIDRULICA DE CANALES INGENIERA CIVIL

Esfuerzo cortante (Pa)

Izquierd
Estacin a Lecho Derecha
k0+500 23.16 90.94 22.92
k0+550 0.98 9.83 3.59
k0+600 1.37 5.17 1.49
k0+650 0.33 2.31 1.16
k0+700 4.83 28.80 3.21
k0+750 89.00 376.74 94.66
k0+800 7.26 76.81 6.22
k0+850 5.37 116.98 12.91
k0+900 0.81 4.16 0.76
k0+907 0.91 5.20 1.14
k0+950 2.80 22.95 1.92
k0+1000 5.56 34.34 5.83
Tabla 6. Esfuerzo cortante en el canal para un caudal de 33m3/s, relacionado a un periodo de
retorno de 100 aos.

En la Tabla 6 se evidencia que el mayor esfuerzo cortante se presenta en la estacin ubicada

en la abscisa k0+750, a continuacin se muestra la seccin transversal correspondiente a
esta estacin.

Imagen 6. Seccin transversa en k0+750.

Este parmetro es notoriamente mayor en la estacin k0+750 ya que esta presenta la mayor
pendiente como se puede apreciar en la Imagen 13; en la Tabla 6 se muestra que el mayor
esfuerzo cortante se presenta en el lecho central del canal, esto debido a que en esta parte se
presenta el mayor coeficiente de rugosidad y que tiene un mayor radio hidrulico
comparado con las bancas de izquierda y derecha, la banca derecha presenta un esfuerzo
cortante levemente mayor a la izquierda a pesar de sus reas de contacto con el agua ser

7
MODELACIN HIDRULICA HEC-RAS
HIDRULICA DE CANALES INGENIERA CIVIL

similares, esto debido a la diferencia entre coeficientes de rugosidad entre estas dos, el
potrero en la banca derecha ejerce una mayor resistencia al avance que el concreto de la
izquierda.

Modelacin con secciones interpoladas

A continuacin se presenta una vista en planta del canal en la cual se observan las secciones
que el programa interpol:

Imagen 7. Eje del canal interpolado en planta y ubicacin de las secciones analizadas e
interpoladas distanciadas cada 10 metros.

Las lneas en amarillo corresponden a las secciones interpoladas por el programa, y las de
color caf a las introducidas por nosotros.
Estacin v (m/s) A (m2) T (m) F Y (m) Yc (m) E (m) (Pa) Pm (m) Rh (m)
k0+500 3.67 9.43 7.39 0.94 1.93 1.93 2.59 90.94 8.65 1.09
k0+550 1.31 39.48 37.82 0.27 2.56 1.97 2.62 9.83 40.82 0.97
k0+600 1.00 50.74 40.98 0.17 3.58 3.62 5.17 44.87 1.13
k0+650 0.67 73.34 42.87 0.11 3.75 3.77 2.31 47.03 1.56
k0+700 2.13 22.92 60.00 0.51 1.83 1.83 2.00 28.80 62.37 0.37
k0+750 6.82 4.98 6.25 2.35 0.97 1.91 3.29 376.74 7.20 0.69
k0+800 3.38 12.21 35.04 0.85 1.87 2.01 2.36 76.81 36.64 0.33

8
MODELACIN HIDRULICA HEC-RAS
HIDRULICA DE CANALES INGENIERA CIVIL

Estacin v (m/s) A (m2) T (m) F Y (m) Yc (m) E (m) (Pa) Pm (m) Rh (m)
k0+850 4.16 8.78 10.79 1.03 2.25 2.68 3.11 116.98 12.60 0.70
k0+900 0.87 46.72 53.43 0.17 3.28 3.31 4.16 56.39 0.83
k0+907 0.98 46.39 53.36 0.18 3.27 3.31 5.20 56.30 0.82
k0+950 1.94 21.45 40.80 0.44 2.11 2.27 22.95 43.86 0.49
k0+1000 2.30 20.34 50.67 0.51 2.38 2.38 2.58 34.34 93.47 0.22
Tabla 7. Parmetros hidrulicos modelo interpolado para el perodo de retorno de 100 aos y
caudal de 33 m3/s.

Comparando las Tablas 5 y 7 se observan algunas variaciones no muy significativas pero

que si pueden ser apreciables para clculos con base en estos resultados. Se observa
tambin que en el modelo interpolado el programa arroja menos valores de profundidades
crticas.
En general, los valores que arroja el programa con el modelo interpolado son un poco
menores a los obtenidos anteriormente sin interpolar secciones.
A continuacin se presenta una seccin interpolada y la siguiente ms prxima que fue
introducida:

Imagen 8. Estacin 4.4 (k0+880) interpolada.

9
MODELACIN HIDRULICA HEC-RAS
HIDRULICA DE CANALES INGENIERA CIVIL

Imagen 9. Estacin 4 (k0+900).

Se observan similitudes geomtricas en ambas secciones, las interpoladas buscan formas

parecidas a las reales y el programa hace una aproximacin gradual a la forma de estas.

Altura mnima de diques o muros para prevenir inundaciones

Como es de esperarse, las mayores inundaciones ocurren para los caudales mayores, en
nuestro caso el correspondiente al perodo de retorno de 100 aos (33 m 3/s). La siguiente
tabla presenta las alturas mnimas necesarias para evitar inundaciones a cada lado del canal:
Tr (aos) Q (m3/s) H mn. izquierdo (m) H mn. derecho (m)
2.33 15 0.26 1.25
25 19 0.5 1.5
50 25 0.62 1.66
100 33 0.75 2.12
Tabla 8. Alturas mnimas de muros o diques a cada borde del canal.

Tr: perodo de retorno.

Q: Caudal.
H mn.: Altura mnima de la estructura (dique o muro) para evitar inundacin en la llanura.

Las alturas mnimas se presentan para cada borde del canal. La seccin ms crtica para el
lado izquierdo del canal corresponde a la estacin 10 (k0+600), y en el lado derecho en la
estacin 9 (k0+650). A continuacin se presentan imgenes para estas secciones
correspondientes al perodo de retorno de 100 aos:

10
MODELACIN HIDRULICA HEC-RAS
HIDRULICA DE CANALES INGENIERA CIVIL

Imagen 10. Estacin 9 (k0+650). Inundacin ms crtica al lado derecho del canal para un
perodo de retorno de 100 aos.

Imagen 11. Estacin 10 (k0+650). Inundacin ms crtica al lado izquierdo del canal para un
perodo de retorno de 100 aos.

11
MODELACIN HIDRULICA HEC-RAS
HIDRULICA DE CANALES INGENIERA CIVIL

Comparacin perfiles de flujo y velocidad

Perfiles de flujo: Las siguientes imgenes presentan los perfiles de flujo para el modelo
interpolado y sin interpolar, correspondientes ambos al perodo de retorno de 100 aos y
caudal de 33 m3/s.

Imagen 12. Perfil longitudinal interpolado del canal para un perodo de retorno de 100 aos y
caudal de 33 m3/s.

12
MODELACIN HIDRULICA HEC-RAS
HIDRULICA DE CANALES INGENIERA CIVIL

Imagen 13. Perfil longitudinal sin interpolar del canal para un perodo de retorno de 100 aos y
caudal de 33 m3/s.

La lnea roja corresponde a la lnea de profundidades crticas, y la verde a la lnea de

energa. El flujo circula de derecha a izquierda en el perfil. Para ambos perfiles se observan
las partes del canal en donde se presenta flujo subcrtico y supercrtico y como es esta
variacin a lo largo de la distancia. El programa no calcula la lnea de profundidades
crticas para todas las secciones, solo para algunas. Se observan cambios ms suaves a lo
largo del perfil interpolado, las variaciones en la lmina de agua son ms graduales y el
segundo resalto hidrulico aguas abajo ms notorio. En el inicio de los perfiles tambin se
observa un resalto hidrulico ms notorio para el interpolado. Las lneas de profundidades
crticas y de energa tambin presentan cambios ms suaves para el perfil interpolado.

Perfiles de velocidad a lo largo del canal: Las siguientes imgenes muestran los perfiles
de velocidad para los modelos interpolados y sin interpolar.

Imagen 14. Perfil de velocidad interpolado para el perodo de 100 aos y caudal de 33 m3/s.

13
MODELACIN HIDRULICA HEC-RAS
HIDRULICA DE CANALES INGENIERA CIVIL

Imagen 15. Perfil de velocidad sin interpolar para el perodo de 100 aos y caudal de 33 m3/s.

La lnea azul corresponde a la velocidad del canal en su centro, la negra a la parte izquierda
y la verde punteada a la parte derecha. En general se observa que la velocidad es menor en
el modelo interpolado.
Para ambos modelos se observa que la velocidad es mayor en el centro del canal, adems
de que la velocidad de la parte derecha del tramo es menor a la de la izquierda. Esto ocurre
ya que el esfuerzo cortante es mayor en la parte derecha del canal, haciendo esto que la
velocidad en dicha zona sea la menor.
En el modelo interpolado se presenta una variacin de velocidad con cambios ms suaves
en comparacin con el modelo sin interpolar. Tambin se observa en el interpolado que en
la distancia del eje de 200 metros, correspondiente a la abscisa k0+800, la velocidad
disminuye para luego subir notoriamente. Esto se da ya que en esta parte del canal ocurre
un resalto hidrulico el cual se observa mejor en el modelo interpolado. Tambin se observa
que en el inicio del canal la velocidad del interpolado sube y luego cae suavemente, en
cambio en el modelo sin interpolar simplemente desciende desde su inicio.

Conclusiones
Como era de esperarse a mayor caudal mayor es la mancha de inundacin,
consideracin importante a la hora de realizar cualquier obra hidrulica pero no
siendo este el nico factor influyente; no es prudente disear para un perodo de
retorno de 100 aos una obra cuya vida til sea menor a este, entonces se pueden
evitar sobrecostos diseando utilizando como criterio los caudales correspondientes
a perodos de retorno menores.

14
MODELACIN HIDRULICA HEC-RAS
HIDRULICA DE CANALES INGENIERA CIVIL

El sector de mayores esfuerzos cortantes se presenta entre la abscisa k0+700 y la

abscisa k0+800, en esta parte del canal se presentan dos cadas y la topografa forma
secciones encaonadas, el agua gana velocidad y energa y el rea de flujo es menor
a la de las dems secciones en general.
La modelacin con secciones interpoladas determina cambios ms suaves en los
perfiles de flujo y de velocidad. Se aprecian mejor los resaltos hidrulicos presentes
en el canal, y en general la velocidad para modelos interpolados es menor que para
modelos sin interpolar.

Bibliografa

Chow, V. T. (1994). Hidrulica de Canales Abiertos. Bogot: McGraw - Hill.

Duarte, C. (2009). Introduccin a la Hidrulica de Canales. Bogot: Universidad Nacional de
Colombia Sede Bogot.

15