Cálculo de Ejes Hidráulicos en Lechos Macro Rugosos - Patricio Pastene

UNIVERSIDAD DE CHILE
FACULTAD DE CIENCIAS FSICAS Y MATEMTICAS

DEPARTAMENTO DE INGENIERA CIVIL
DISEO E IMPLEMENTACIN DE SOFTWARE DE ACOPLE A HEC-

RAS: CLCULO DE EJES HIDRULICOS EN LECHOS MACRO
RUGOSOS
MEMORIA PARA OPTAR AL TTULO DE INGENIERO CIVIL

PATRICIO ANTONIO PASTENES ARAOS
PROFESOR GUA:
ALBERTO DE LA FUENTE STRANGER
MIEMBROS DE LA COMISIN:
ALDO TAMBURRINO TAVANTZIS
MARCELO OLIVARES ALVEAL
SANTIAGO DE CHILE
JUNIO 2012
RESUMEN DE LA MEMORIA
PARA OPTAR AL TTULO DE
INGENIERO CIVIL.
POR: PATRICIO PASTENES A.
FECHA: 20/06/2012.
PROF. GUA: Sr. ALBERTO DE LA FUENTE.
DISEO E IMPLEMENTACIN DE SOFTWARE DE ACOPLE A HEC-RAS:

CLCULO DE EJES HIDRULICOS EN LECHOS MACRO RUGOSOS
En la presente memoria se brinda una herramienta computacional de acople a HEC-RAS,

denominada MR-HR, la cual es capaz de resolver escurrimientos sobre lechos tanto rugosos
como macro rugosos en cauces naturales, mediante la implementacin de diferentes expresiones
no lineales que dan cuenta del efecto de la macro rugosidad.
La primera etapa del estudio, consisti en la recopilacin de antecedentes generalizados de

prdida de carga en cauces naturales. Posteriormente se recopilaron los antecedentes de prdida
de carga asociados con lechos macro rugosos, particularmente, las parametrizaciones del nmero
de Strickler de Keulgean (1938), Limerinos (1970), Parker y Peterson (1980) y Ayala y Oyarce
(1993). Estas expresiones fueron ajustadas para lograr una transicin suave entre lechos macro
rugosos y rugosos, forzando a que el nmero de Strickler en lechos rugosos fuera igual a 0.12.
Luego se desarroll una aplicacin en un caso real en un tramo particular del ro Baker, de
148 km de extensin, comprendidos entre Baker 1 y la desembocadura en caleta Tortel. Esta
aplicacin se dividi en dos partes. En la primera, se simplific el lecho del cauce asumiendo un
slo dimetro caracterstico para todo ste, lo que implic que la macro rugosidad vari
nicamente por el desarrollo del eje hidrulico. En la segunda, se abord el problema lo ms
general posible, admitiendo que la macro rugosidad variar tanto en funcin del desarrollo del eje
hidrulico como en funcin de la variacin longitudinal del tamao de los sedimentos del lecho.
El sentido de dividir la aplicacin en dos partes fue que en la primera etapa, simplificada, se
busc analizar los factores de los que dependan los ejes hidrulicos proporcionados por las
distintas parametrizaciones, y en la segunda etapa, se busc corroborar lo anterior y comparar los
resultados con el eje hidrulico proporcionado por un coeficiente de Manning uniforme para todo
el cauce que obvi los efectos de la macro rugosidad.
Los resultados de la primera etapa sugirieron que el desarrollo de ejes hidrulicos

considerando las parametrizaciones aqu revisadas, es un problema que carece de generalidad en
trminos de que no siempre es posible esperar que una parametrizacin brinde alturas de
escurrimiento mayores que las dems. Esto se debe a que aparte de la parametrizacin
seleccionada, el desarrollo del eje hidrulico est influenciado por la distribucin espacial de los
nmeros de Strickler que a su vez dependen de otros factores como el caudal, la geometra
particular del cauce y las condiciones de borde del problema. En la segunda etapa, se
corroboraron las hiptesis anteriores para un caudal en particular, y con respecto a la
comparacin entre considerar un lecho como macro rugoso y no considerarlo, se obtuvieron
diferencias significativas dada la comparacin desigual realizada.
Finalmente se menciona que los resultados obtenidos son particulares y que por ningn
motivo deben extrapolarse a otros problemas sometidos a distintas condiciones, y se
proporcionan consideraciones que se deberan tener en cuenta al momento de evaluar los
resultados brindados por MR-HR.
i
A MI HERMANA CAROLINA
A QUIEN AMO
ii
AGRADECIMIENTOS
Quiero expresar mis agradecimientos a todas las personas que contribuyeron en mi proceso
de formacin acadmica, en particular:
A mi familia y polola por su incesable apoyo a lo largo de toda mi formacin.
A mi profesor gua, Sr. Alberto de la Fuente S., por su dedicacin, comprensin,

disposicin, contribucin al desarrollo de este trabajo y su aporte en mi formacin
acadmica.
A los miembros que integran esta comisin y al Sr. Yarko Nio, por sus aportes y
sugerencias al desarrollo de este trabajo.
A CEA limitada, por las facilidades brindadas para poner trmino a esta memoria.
iii
INDICE DE CONTENIDOS
CAPTULO 1. INTRODUCCIN ...................................................................................1
1.1 Motivacin.................................................................................................................1
1.2 Objetivos y Alcances .................................................................................................3
1.2.1 Objetivo Principal......................................................................................................3
1.2.2 Objetivos Especficos ................................................................................................3
1.2.3 Alcances ....................................................................................................................4
1.3 Contenido del Informe...............................................................................................4
CAPTULO 2. REVISIN DE ANTECEDENTES .......................................................5
2.1 Prdida de Carga en Lechos Macro Rugosos ............................................................5

2.2 Diferencia entre Lechos Rugosos y Macro Rugosos...............................................11
2.3 Parametrizaciones Consideradas .............................................................................13
2.3.1 Keulegan (1938) ......................................................................................................13
2.3.2 Limerinos (1970) .....................................................................................................15
2.3.3 Parker y Peterson (1980) .........................................................................................18
2.3.4 Ayala y Oyarce (1993) ............................................................................................20
2.4 Resumen de Parametrizaciones Consideradas.........................................................22
CAPTULO 3. IMPLEMENTACIN DE LA HERRAMIENTA: MR-HR .............26
3.1 Ajuste de Parametrizaciones....................................................................................26

3.1.1 Ajuste Lineal ...........................................................................................................27
3.1.2 Ajuste Parablico.....................................................................................................29
3.1.3 Ajuste Ayala y Oyarce.............................................................................................31
3.1.4 Resumen Coeficientes Lineales y Parablicos ........................................................32
3.1.5 Ajustes Inverso de Rugosidad Relativa Mnimo .....................................................35
3.1.6 Resumen de Parametrizaciones Ajustadas ..............................................................36
3.2 Atributos MR-HR ....................................................................................................36
3.3 Relaciones entre MR-HR y HEC-RAS ...................................................................38
3.4 Esquema de Funcionamiento...................................................................................40
3.5 Resultados Proporcionados .....................................................................................41
CAPTULO 4. APLICACIN AL RO BAKER. ........................................................43
4.1 Antecedentes Generales Ro Baker .........................................................................43

4.2 Tramo Particular Estudiado .....................................................................................46
4.3 Caracterizacin Granulomtrica ..............................................................................47
4.4 Clculo de Ejes Hidrulicos Mediante MR-HR ......................................................55
CAPTULO 5. DISCUSIN Y ALCANCES ................................................................70
BIBLIOGRAFA ...................................................................................................................72
ANEXO A. MANUAL DE USO MR-HR .......................................................................... I
A.1 Generalidades ..............................................................................................................I

A.2 Instalacin....................................................................................................................I
A.3 Requisitos .................................................................................................................. II
A.3.1 Computacionales ....................................................................................................... II
A.3.2 Proyecto Importado desde HEC-RAS ......................................................................III
A.4 Componentes .......................................................................................................... VII
A.4.1 Ventana Principal .................................................................................................. VIII
A.4.2 Submen Archivo.....................................................................................................IX
A.4.3 Submen Configuracin......................................................................................... XII
A.4.4 Submen Ejecutar................................................................................................ XVII
A.4.5 Submen Resultados .......................................................................................... XVIII
A.4.6 Submen Ayuda .................................................................................................XXIV
A.5 Ejemplo de Clculo de MR-HR .........................................................................XXVI
A.6 Errores Conocidos .............................................................................................XXXII
ANEXO B. ARCHIVOS DIGITALES ...................................................................XXXIII
B.1 MR-HR ............................................................................................................ XXXIII

B.2 Cdigo Fuente MR-HR ................................................................................... XXXIII
B.3 Ejemplo Aplicacin al Ro Baker.................................................................... XXXIV
B.4 Ejemplo de Clculo Manual de Uso ................................................................ XXXIV
B.5 Resultados Completos Aplicacin al Ro Baker ............................................. XXXIV
ANEXO C. HEC-RAS Y VISUAL BASIC EXPRESS 2010 .................................XXXVI

Clculos de Ejes Hidrulicos en Lechos Macro Rugosos
Patricio Pastenes Araos Captulo 1: Introduccin
Captulo 1. Introduccin
1.1 Motivacin
Este trabajo trata acerca del clculo de ejes hidrulicos en cauces fluviales con lechos
macro rugosos, es decir, lechos con rugosidades relativas, definidas como el cuociente entre el
dimetro caracterstico de sedimento grueso y el radio hidrulico, extremadamente altas
(superiores a 0.1). Este tipo de lechos son muy habituales en Chile, en particular, en
escurrimientos de alta montaa con altas pendientes donde se presentan escurrimientos
torrenciales en los cuales el radio hidrulico se hace comprable con el tamao del sedimento
grueso. Tambin se extiende su aplicabilidad a cauces con lechos rugosos, aquellos con
rugosidades relativas inferiores a 0.1.
La macro rugosidad se caracteriza por la presencia de grandes partculas de sedimento que

generan despegues y estelas en el flujo en torno a stas. Esto provoca un aumento considerable en
el coeficiente de friccin que en este trabajo de titulo es representado por el coeficiente de
Manning. El nmero de Strickler representa un parmetro adimensional que relaciona el
coeficiente de Manning con el tamao de los sedimentos gruesos y el radio hidrulico. Por lo
tanto, para conocer el coeficiente de Manning se necesita conocer el nmero de Strickler y para
conocer este ltimo se necesita conocer previamente el radio hidrulico o en forma equivalente el
desarrollo del eje hidrulico y algn dimetro caracterstico del sedimento grueso. En cauces
naturales, en general, el tamao de los sedimentos de fondo son variables en funcin de su
desarrollo y adems los ejes hidrulicos no son uniformes. Esto produce un problema iterativo
que tanto softwares comerciales como no comerciales no han abordado del todo. Por lo que hoy
en da el calculista debe realizar las iteraciones en forma manual, lo cual se incrementa en
funcin del nmero de secciones transversales a las que se vea enfrentado. Esto implica que, en la
prctica, el calculista no realice ninguna iteracin, proponiendo un valor uniforme para el
coeficiente de Manning y adopta como verdadero el valor del eje hidrulico asociado a ese
Manning.
Distintos investigadores, entre ellos Keulegan (1938), Limerinos (1970), Parker y Peterson
(1980) y Ayala y Oyarce (1993), han analizado este tipo de lechos macro rugoso, y han planteado
1
distintas relaciones no lineales para el clculo del nmero de Strickler, de aqu en adelante
llamadas parametrizaciones del nmero de Strickler, en funcin del tamao de las partculas de
sedimento grueso y el radio hidrulico.
Para el clculo de ejes hidrulicos en lechos de granulometras gruesas existen en general

dos enfoques principales. El primero es denominado mtodo de manual, en donde el ingeniero
recurre a la literatura y selecciona un coeficiente de Manning tabulado por diversos
investigadores o si dispone de granulometra calcula coeficientes de Manning bsicos n0 y luego
se modifican para hacerlos representativos del cauce. El segundo enfoque, el cual se aborda en
este trabajo de ttulo, trata acerca de la clculo del coeficiente de Manning mediante un enfoque
analtico, el cual utiliza tanto informacin granulomtrica como informacin del escurrimiento
del cauce y se resuelve en base a iteraciones sucesivas considerando el efecto de la macro
rugosidad.
Ambos enfoques tienen asociados grados de subjetividad. En el primero se presenta la

decisin de la seleccin de coeficientes de Manning a partir de tabulaciones, lo cual puede variar
en funcin del criterio particular de quien aborda el problema. En el segundo enfoque, la
subjetividad surge a la hora de seleccionar alguna parametrizacin del nmero de Strickler por
sobre las dems. Esta seleccin debe ser acorde con condiciones similares bajo las cuales fueron
deducidas. Sin embargo, el grado de similitud entre cauces es lo que convierte a este enfoque en
uno tambin subjetivo. Para mejorar esto, se consider un parmetro de calibracin de Manning
el cual se adiciona a los coeficientes brindados por las diferentes parametrizaciones y tiene como
objetivo dar cuenta de particularidades del cauce, como obstrucciones, vegetacin, curvas, etc.
Por ejemplo, en el contexto de Cowan (Chow, 1994) el coeficiente de Manning bsico sera el
representado por las diferentes parametrizaciones y la calibracin dara cuenta de las
particularidades del cauce. En este mismo contexto si se cuenta con informacin granulomtrica a
lo largo del cauce, una seleccin adecuada de alguna parametrizacin y una buena calibracin, el
segundo enfoque tiene la ventaja de reflejar la variabilidad del coeficiente de Manning bsico en
estos tipos de lechos.
Entonces, la motivacin de este trabajo de titulo es presentar una herramienta que sea capaz
de extender la aplicacin de la ecuacin de Manning a un rango para el cual no fue deducida,
mediante la implementacin de diferentes parametrizaciones del nmero de Strickler. Adems se
2
desprenden una serie de consecuencias asociadas al clculo adecuado de coeficientes de

Manning, por ejemplo un clculo del eje hidrulico ms preciso que podr ser considerado para el
diseo de obras de ingeniera civil, definir caudales ecolgicos u otros.
1.2 Objetivos y Alcances
1.2.1 Objetivo Principal
El objetivo principal de este trabajo es facilitar, mediante el diseo e implementacin de un

software, el proceso iterativo de clculo de ejes hidrulicos en cauces fluviales con lechos macro
rugosos, en los cuales se considera que el nmero de Strickler es una funcin no lineal que
depende del previo conocimiento del clculo del eje hidrulico y la granulometra del lecho.
1.2.2 Objetivos Especficos
Disear e implementar un software capaz de acoplarse al software HEC-RAS y realizar

los clculos en forma iterativa, sin necesidad que el usuario realice el proceso iterativo de
clculo del nmero de Strickler en forma manual y consecuentemente el clculo del
coeficiente de Manning.
Presentar una herramienta que ample el espectro de toma decisiones en relacin a los
clculos de ejes hidrulicos, debido a que se presentan distintas alternativas, representadas
mediante diferentes parametrizaciones del nmero de Strickler, para el clculo de ejes
hidrulicos.
Realizar una aplicacin real a un cauce de secciones transversales conocidas y de lecho

macro rugoso, en donde se comparen los resultados del eje hidrulico entre las distintas
parametrizaciones del nmero de Strickler propuestas por los diferentes autores.
3
1.2.3 Alcances
Los alcances a los que est ligada la herramienta de acople a HEC-RAS, denominada MR-
HR, son:
No modifica el cdigo fuente de HEC-RAS, sino que lo llama mediante Windows.
No implementa algoritmo de clculo de ejes hidrulicos, sino que los obtiene mediante la
interaccin con HEC-RAS.
Resuelve problemas en rgimen permanente, con flujos turbulentos plenamente

desarrollados, lechos con pared hidrodinmicamente rugosa, granulometras gruesas
(gravas), lechos sin formas de fondo, flujos gradualmente variados y uniformes.
Slo utiliza las siguientes parametrizaciones del nmero de Strickler: Keulegan (1938),
Limerinos (1970), Parker y Peterson (1980) y Ayala y Oyarce (1993).
Adems est limitado por la interaccin con HEC-RAS, lo cual se desarrolla en forma
extensa en la seccin (3.3) la cual trata acerca de las relaciones existentes entre MR-HR y HEC-
RAS. Los requerimientos computacionales se muestran en el anexo (A.3).
1.3 Contenido del Informe
El resto del presente informe consta de las siguientes partes: En el Captulo 2 se realiza una
revisin de antecedentes de prdidas de cargas debido al efecto del sedimento del lecho y, en
particular, se enfoca en las parametrizaciones del nmero de Strickler que son implementadas en
el software MR-HR. En el Captulo 3 se define la implementacin de la herramienta MR-HR
junto con su forma de operacin y ajustes a las parametrizaciones anteriormente sealadas en el
Captulo 2. En el Captulo 4 se realiza una aplicacin a un tramo del ro Baker, comparando los
resultados brindados por las diferentes parametrizaciones seleccionadas y, finalmente, en el
Captulo 5 se desarrollan las conclusiones y alcances del trabajo.
4
Patricio Pastenes Araos Captulo 2: Revisin de Antecedentes
Captulo 2. Revisin de Antecedentes
2.1 Prdida de Carga en Lechos Macro Rugosos
Los antecedentes recopilados son vlidos para granulometras gruesas, en particular, para
lechosos gravosos con rugosidades relativas extremadamente altas, lo cual se denomina lecho
macro rugoso. Este tipo de lechos se presenta, principalmente, en zonas montaosas de gran
pendiente, donde la rugosidad relativa toma una valoracin especial dado que se hace comparable
el radi hidrulico con algn dimetro caracterstico de sedimentos gruesos. El dimetro
caracterstico del sedimento grueso que genera la macro rugosidad es el de la coraza, sin embargo
tambin es posible relacionarlo con algn dimetro caracterstico de la granulometra integral del
lecho.
El lmite que define un lecho macro rugoso, depende de la rugosidad relativa, la cual se
define como:
ds
Rugosidad Relativa = (2.1)
Rh
Donde:
Rh : corresponde al radi hidrulico.
d s : corresponde al dimetro caracterstico del sedimento grueso, tpicamente D90 o D84 para
granulometras integrales del lecho o D50 para granulometras de la coraza, dependiendo del
autor.
Cuando la rugosidad relativa es mayor o igual 0.1 se define un lecho como macro rugoso.
Se pueden diferenciar dos efectos dominantes que afectan a la resistencia hidrulica en este tipo
de lechos: la resistencia friccional clsica que se refiere a la resistencia hidrulica producida por
el sedimento del lecho propiamente tal, y una resistencia de forma que est asociada a distintos
patrones de formas del lecho, como por ejemplo: rizos, dunas, anti dunas, secuencias rpidos-
5
pozas, etc. Las primeras tres formas de fondo anteriormente sealadas se producen por acarreo de
sedimento fino como la arena, y las secuencias rpidos-pozas se asocian a lechos gravosos. Sin
embargo, en este trabajo no se consider el efecto de formas de fondo y slo se centr en el
efecto que produce el sedimento propiamente tal. Bajo este enfoque, el efecto que produce la
macro rugosidad es que se generan despegues y estelas en el flujo en torno a los sedimentos
gruesos, que implica la existencia de una resistencia de forma debido a la fuerza de arrastre que
se genera sobre los elementos de rugosidad del lecho que se agrega a la resistencia friccional
clsica de un lecho hidrodinmicamente rugoso. En estas condiciones se ha determinado
experimentalmente que la resistencia hidrulica del lecho aumenta considerablemente (Nio,
2003).
Con el propsito de facilitar la presentacin de los antecedentes recopilados, a continuacin

se define la nomenclatura utilizada en la descripcin de las variables que intervienen en las
diversas expresiones propuestas en la literatura para predecir la prdida de carga:
Rh : corresponde al radio hidrulico. [m]
ds : corresponde al dimetro caracterstico del sedimento grueso, tpicamente D90

o D84 para granulometras integrales o D50 para granulometras de la coraza,
dependiendo del autor. [m]
S : corresponde a la pendiente del terreno. [m m]
kg
o : corresponde al esfuerzo de corte medio sobre el lecho.
m s
2
U : corresponde a la velocidad media del escurrimiento. [m s ]
u* : corresponde a la velocidad friccional o de corte. [m s ]
n : [
corresponde a el coeficiente de Manning. s m1/ 3 ]
St : corresponde al nmero adimensional de Strickler.
6
Cf : corresponde a un coeficiente o factor de friccin adimensional.
Cchezy : corresponde al coeficiente de Chezy. [ m1 / 2 s ]
y : corresponde a la altura del escurrimiento. [m]
g : corresponde a la aceleracin de gravedad. m s 2 [ ]

: corresponde a la densidad del agua. kg m 3 [ ]
y0 : corresponde a una constante de integracin. [m]
: corresponde a una constante universal caracterstica de la turbulencia o tambin

conocida como la constante de Von Karman, adimensional.
Como se mencion anteriormente, la condicin para que un lecho sea macro rugoso es que
la rugosidad relativa sea mayor a 0.1. En forma equivalente, en este trabajo se opt por el inverso
de la rugosidad relativa, por lo tanto se dice que un escurrimiento se desarrolla sobre un lecho
macro rugoso si:
Rh
10 (2.2)
ds
Para poder caracterizar escurrimientos en lechos macro rugosos, es necesaria la

introduccin de algn coeficiente de friccin que permita dar cuenta de los efectos que produce el
lecho en el escurrimiento. Desde el punto de vista de la ingeniera Civil, el coeficiente de
Manning es el ms usado a pesar que es un coeficiente dimensional. De esta forma, el presente
trabajo de ttulo aborda el clculo de ejes hidrulicos mediante el coeficiente de Manning.
De acuerdo a la ecuacin de Strickler, el coeficiente de Manning depende del tamao del

sedimento del lecho de la siguiente forma:
1/ 6
ds
n= (2.3)
26.4
7
En base a la relacin (2.3), se define un parmetro adimensional que relaciona el

coeficiente de Manning y el dimetro caracterstico del sedimento grueso. Este parmetro se
denomina el nmero de Strickler, S t , y se define como:
g n
St = 1/ 6
(2.4)
ds
Reemplazando la ecuacin (2.3) en la ecuacin (2.4), resulta que el nmero de Strickler es

un nmero igual a 0.12, independiente del tipo de lecho. Este valor fue determinado por Strickler
(1923) y es vlido para lechos rugosos y de granulometras gruesas considerando ds=d90.
Diversos investigadores han estudiado escurrimientos en lechos macro rugosos (Keulegan,

1938; Limerinos, 1970; Parker y Peterson, 1980; Ayala y Oyarce, 1993) y han concluido que, a
diferencia de Strickler de (2.3), en lechos macro rugosos, el nmero de Strickler est fuertemente
influenciado por una variable adicional al dimetro caracterstico de los sedimentos gruesos, la
cual es el radio hidrulico. Esto provoca una variabilidad en el nmero de Strickler.
Entonces, para entender las relaciones que se presentan en estos escurrimientos sobre
lechos macro rugosos, se recurre a una ley de resistencia hidrulica general, la cual plantea la
siguiente relacin:
o = C f U 2 (2.5)
Adems, se tiene que a partir del esfuerzo de corte sobre el fondo se define una velocidad
friccional o de corte, como:
o
u* = (2.6)

Por lo tanto, manipulando las ecuaciones (2.5) y (2.6) se obtiene que:
8
U 1
*
= (2.7)
u Cf
Luego, considerando equilibrio de fuerzas en canales abiertos, para flujos permanentes y

uniformes se tiene que:
o = g Rh S (2.8)
Por lo tanto, reemplazando la ecuacin (2.8) en (2.6) se obtiene la siguiente relacin para la
velocidad friccional o de corte:
u * = g Rh S (2.9)
Adems, considerando las mismas condiciones anteriores, es decir, flujo permanente y

uniforme, se tiene la ecuacin de Manning:
S 2/3
U= Rh (2.10)
n
Reemplazando los resultados de las ecuaciones (2.10) y (2.9) en la ecuacin (2.7) se

obtiene:
2
g n
C f = 1/ 6 (2.11)
R
h
Finalmente, es posible sustituir la ecuacin (2.4) en la ecuacin (2.11), obtenindose:
1/ 6
1 U 1 R
= * = h (2.12)
Cf u St d s
Notar que en la deduccin de la ecuacin (2.12), el coeficiente de friccin adimensional es

una propiedad que solo depende del inverso de la rugosidad relativa. Esto es vlido slo para
flujos turbulentos plenamente desarrollados y, adems, con pared hidrodinmicamente rugosa.
9
Adems, es interesante notar de esta ecuacin que si es posible conocer la distribucin de

velocidades, entonces es posible conocer el nmero de Strickler y con ayuda de la ecuacin (2.4)
es posible conocer el coeficiente de Manning. Precisamente es la distribucin de velocidades, la
cual vara segn los autores que estudiaron este fenmeno: Keulegan (1938); Limerinos (1970);
Parker y Peterson (1980); Ayala y Oyarce (1993). En general, la mayora de los autores sealados
anteriormente concuerdan en la forma del cuociente de velocidades, basado el la ley de
distribucin de velocidades logartmica de Prandtl (1926) von Karman (1930) (ambos citados en
Keulegan, 1938). Sin embargo, difieren entre si en los valores de las constantes asociadas a esta
ley. Esto se debe a que los valores de las constantes fueron determinadas de manera experimental
sujetas a diferentes condiciones de terreno y/o laboratorio. En algunos casos, los estudios se
limitaron a correlacionar datos suponiendo la forma del cuociente de velocidades y se
seleccionaron los coeficientes que mejor se ajustaban a los datos, en relacin a los coeficientes de
Manning calculados con la ayuda tanto de mediciones de terreno y/o laboratorio como de la
ecuacin (2.10).
Por otra parte y en estricto rigor, la aplicabilidad de la ecuacin de Manning est sujeta a
flujos permanentes y uniformes, es decir, a problemas donde la altura de escurrimiento no vara
ni espacialmente, ni temporalmente. Esto implica: cauces prismticos y sin variacin de ninguna
propiedad tanto espacialmente como temporalmente. Sin embargo, ms adelante en la seccin
(3.3) se extender su validez a problemas con escurrimientos gradualmente variados y a cauces
no prismticos.
Con respecto a la forma cmo se calcula el coeficiente de Manning en este trabajo, como se
mencion anteriormente, slo es vlido para flujos turbulentos plenamente desarrollados y con
pared hidrodinmicamente rugosa. Si no se cumple est ltima condicin, el coeficiente de Chezy
y, consecuentemente, el coeficiente de Manning dependern del nmero de Reynolds, dado que
las fuerzas viscosas no son despreciables en la cercana de la pared.
10
2.2 Diferencia entre Lechos Rugosos y Macro Rugosos
Vale la pena aclarar, antes de comenzar con la diferencia entre lechos, que el nmero de
Strickler tambin suele ser interpretado como la constante de Strickler. La diferencia entre ambos
es que el primero es representado por la ecuacin:
g n
St = 1/ 6
ds
Mientras que el segundo es representado por el cuociente:
n
1/ 6
= Constante de Strickler
ds
Notar que ambas definiciones son equivalentes, al nmero de Strickler se le agrega la raz
cuadrada de la gravedad para hacerlo adimensional.
La principal diferencia existente entre lechos rugosos y macro rugosos es la variabilidad de

la constante de Strickler en lechos macro rugosos. En el primer tipo de lecho, con inversos de
rugosidades relativas mayores que 10, es posible asumir que la constante de Strickler se comporta
como tal, y su valor es igual a 0.0342 (Strickler, 1923). En el segundo tipo de lecho, la constante
de Strickler es variable y depende del inverso de la rugosidad relativa. Datos recolectados de
algunos ros del noroeste de Estados Unidos y el ro Mississippi dan cuenta de la diferencia entre
ambos lechos (Chow, 1994), lo cual se muestra en la figura (2.1). Adems se aprecia que el eje de
la abscisas representa el cuociente entre el radio hidrulico y la altura de rugosidad (adimensional
considerando tanto la altura de rugosidad como el radio hidrulico en pies), y en el eje de las
ordenadas el cuociente entre el coeficiente de Manning y la altura de rugosidad elevada a un
sexto. Si se supone que la altura de rugosidad k puede ser representada por algn dimetro
caracterstico ( k = d s ), entonces se recupera la ecuacin (2.2), es decir, el inverso de la rugosidad
relativa.
11
Figura 2.1: Representacin Grfica Constante de Strickler

Fuente: Chow (1994). Hidrulica de Canales Abiertos. Pgina 202.
Es interesante notar que la constante de Strickler para inversos de rugosidades relativas

menores o iguales que 10 aumenta considerablemente, lo cual es cuantificado por las expresiones
de Keulegan (1938), Limerinos (1970), Parker y Peterson (1980) y Ayala y Oyarce (1993). Por
otra parte, para inversos de rugosidades relativas mayores que 10, se tiene que Strickler se vuelve
constante igual al valor 0.0342. Este valor, propuesto por Strickler (1923), surge de
observaciones reales realizadas en ros de Suiza, y adems de considerar la altura de rugosidad
igual al dimetro medio de los sedimentos de fondo, es decir, k = D50 . Estudios posteriores
realizados en Waterways Experiment Station (WES) (Maynard, 1991 citado en USACE, 2010b),
mostraron que el valor de la constante de Strickler, considerando k = D90 y para consideraciones
de clculo de caudales en cauces de lechos gravosos, debe ser 0.038, lo cual se muestra a
continuacin:
Reordenando la ecuacin (2.3), resulta:
n 1
1/ 6
= (2.13)
ds 26.4
con d s en metros. En forma equivalente resulta:
n
1/ 6
= 0.038 (2.14)
ds
12
En la ecuacin (2.14) se recupera el valor de la constante de Strickler 0.038. Luego

reemplazando (2.14) en (2.4) resulta:
S t = 0.12 (2.15)
2.3 Parametrizaciones Consideradas
2.3.1 Keulegan (1938)
Keulegan (1938) propuso la siguiente expresin para el cuociente entre las velocidades
media y friccional:
U R
*
= 2.5 ln12 h (2.16)
u ds
con d s = D90 .
La ecuacin (2.16), puede ser expresada en trminos del nmero de Strickler, mediante la
ayuda de la ecuacin (2.12) y resulta:
1/ 6
Rh

St = ds (2.17)
R
2.5 ln12 h
ds
Rh
con d s = D90 . El valor de S t versus se muestra en la figura (2.2).
ds
La ecuacin (2.16) est basada en la ley de distribucin de velocidades logartmica de

Prandtl (1926) von Karman (1930) (ambos citados en Keulegan, 1938). Est ley es la siguiente:
13
U 1 y
= ln (2.18)
u * y0
Consideraciones tericas indican que la constante debera ser independiente de la

naturaleza de la superficie de la pared. Esto ha sido verificado en casos de flujos en tuberas
circulares por von Karman (Karman, 1930 citado en Keulegan, 1938). Con respecto a la
constante de integracin yo , Keulegan la denomin "caracterstica de la superficie", sta vara en
funcin de la velocidad friccional y de acuerdo a la rugosidad de la pared (Keulegan, 1938).
Para determinar la constante yo , Keulegan se apoy en estudios experimentales

proporcionados por Nikuradse (1933) (citado en Keulegan, 1938), quien realiz
experimentaciones en tuberas circulares con arena. Finalmente, Keulegan corrobor la ecuacin
(2.16) mediante experimentaciones realizadas por Bazin (1865) (citado en Keulegan, 1938) en
canales abiertos, entre los aos 1855 al 1860. La experimentacin mediante la cual se valid la
ecuacin (2.16), comprende los siguientes rangos: pendientes entre [0.0015 - 0.0089], inversos de
rugosidad relativas entre [6.3 - 1,030], radios hidrulicos entre [0.048 - 0.396] m, velocidades
medias entre [0.39 - 2.61] m/s, temperaturas entre [7.5 - 19] C, secciones transversales
prismticas rectangulares y de formas variadas. De esta forma, los requisitos que debe cumplir,
simultneamente, la ecuacin (2.16) para ser vlida son: lechos con granulometras finas1
(arenosas), flujos turbulentos plenamente desarrollados, lechos con pared hidrodinmicamente
rugosa y ausencia de formas de fondo.
1
En estricto rigor, Keulegan basado en experimentacin de Nikuradse con arenas, determin en su publicacin el
concepto de rugosidad de arena equivalente. Con esto es posible evaluar problemas con diferentes tipos de
granulometras, aparte de arenas. Sin embargo, en este trabajo al fijar d s = 1 D90 solo tiene sentido para
granulometras arenosas.
14
Figura 2.2: Representacin Grfica Ecuacin (2.17)

Fuente: Elaboracin Propia
2.3.2 Limerinos (1970)
Limerinos (1970) propuso la siguiente expresin para el cuociente entre las velocidades
media y friccional:
U R
*
= 2.5 ln 3.8 h (2.19)
u ds
con d s = D84 .
15
1/ 6
Rh

St = ds (2.20)
R
2.5 ln 3.8 h
ds
Rh
ds
La ecuacin (2.19) est basada en un estudio experimental realizado en terreno. La

expresin original propuesta por Limerinos fue:
n 0.0926
1/ 6
= (2.21)
Rh R
1.16 + 2.0 log h
ds
La expresin (2.21) est en unidades inglesas, sin embargo es posible transformarla al

sistema internacional y expresarla en trminos del cuociente de velocidades mediante las
siguientes expresiones:
1/ 6
1.486 Rh
Cchezy = (2.22)
n
g
Cf = 2
(2.23)
Cchezy
Finalmente utilizando la ecuacin (2.7) es posible deducir la expresin (2.19), la cual est
en trminos del cuociente de velocidades.
El estudio experimental realizado por Limerinos abarc once cauces naturales en

California, Estados Unidos, donde se midieron: caudales, alturas de escurrimiento, velocidades,
etc. La deduccin de la expresin (2.21) o su equivalente (2.19) se bas en que Limerinos saba la
forma que debera tener la relacin buscada, debido a la revisin bibliogrfica realizada
previamente. Por lo tanto, se limit a correlacionar los datos y escogi los coeficientes que mejor
16
se ajustaban ellos. Es interesante destacar que estos coeficientes, si son vistos del punto de vista
terico segn la ecuacin (2.18), corresponden a la constante de von Karman y a la constante de
integracin que es caracterstica de la superficie (Keulegan, 1938). Bajo este punto de vista, solo
se aprecia una diferencia con respecto a la ecuacin brindada por Keulegan (2.16) en la constante
de integracin, lo cual se debe a que Keulegan (1938) determin el valor de esta constante para
arenas y Limerinos para gravas.
La experimentacin de Limerinos mediante la cual se valid la ecuacin (2.19), comprende

los siguientes rangos: rango de inversos de rugosidades relativas entre [1.06 - 68.05], caudales
entre [5.97 - 254.85] m3/s y dimetros caractersticos entre [0.02 - 0.75] m. Los requisitos que
debe cumplir, simultneamente, la ecuacin (2.19) para ser vlida son: lechos con granulometras
gruesas (gravas), flujos turbulentos plenamente desarrollados, lechos con pared
hidrodinmicamente rugosa, ausencia de formas de fondo, cauces rectos, con poca o nula
vegetacin, cauces relativamente anchos de seccin trapezoidal y cauces sin irregularidades tanto
de fondo como en las riberas. Limerinos en su estudio escogi cauces con todas las caractersticas
descritas anteriormente, debido a que buscaba una relacin que cuantificara el efecto aislado de la
resistencia que produce el sedimento de fondo. Su idea se basaba en que si lograba calcular en
forma aislada el efecto del sedimento de fondo, luego podra corregir dicho valor en presencia de
cauces: sinuosos, con vegetacin, irregularidades, etc., mediante una expresin propuesta en la
literatura, como por ejemplo, la de Cowan (Chow, 1994).
17

2.3.3 Parker y Peterson (1980)
Parker y Peterson (1980) propusieron la siguiente expresin para el cuociente entre las
velocidades media y friccional, debida solamente al efecto en la resistencia producida por el
sedimento de fondo (sin incluir formas de fondo):
U R
*
= 2.46 ln 5.5 h (2.24)
u ds
con d s = D90 .
18
1/ 6
Rh

St = ds (2.25)
R
2.46 ln 5.5 h
ds
Rh
ds
El estudio de Parker y Peterson se centra en predecir la prdida de carga producida por una
forma de fondo particular, la cual obedece a secuencias de rpidos-pozas. Un problema recurrente
que se enfrenta en estos estudios es que diversos factores intervienen en la prdida de carga,
aparte de los efectos de las formas de fondo y del sedimento del lecho. Para solucionar esto,
Parker y Peterson corroboraron que en los cauces naturales estudiados en Alberta, Canad, los
efectos dominantes que contribuyen a la resistencia son la forma de fondo rpido-poza y el
sedimento del lecho. Entonces proponen linealizar ambos efectos y basarse en estudios previos,
para una expresin que de cuenta de la prdida de carga producida por el sedimento del lecho.
Con esto, y por diferencia, son capaces de cuantificar el efecto producido por las formas de fondo
rpidos-pozas.
Los estudios previos en los que se basaron Parker y Peterson para la deduccin de la
ecuacin (2.24), fueron los proporcionados por: Keulegan (1938), Limerinos (1970) y Kamphuis
(1974) (citado en Parker y Peterson, 1980).
Por otra parte, para la prediccin de la expresin que de cuenta de la prdida de carga por la
forma de fondo particular, Parker y Peterson se basaron en mediciones de terreno brindados por
Kellerhals, Neill y Bray (1972) (citado en Parker and Peterson, 1980). Sin embargo, esta
informacin es representativa de cauces con secuencias de rpidos-pozas, las cuales no son
consideradas en este trabajo de ttulo. La expresin (2.24) est sujeta a condiciones similares las
propuestas por Limerinos (1970).
19
Los requisitos que debe cumplir, simultneamente, la ecuacin (2.24) para ser vlida son:
lechos con granulometras gruesas (gravas), flujos turbulentos plenamente desarrollados, lechos
con pared hidrodinmicamente rugosa y ausencia de formas de fondo.

2.3.4 Ayala y Oyarce (1993)
Ayala y Oyarce en 1993 propusieron la siguiente expresin para el cuociente entre las
velocidades media y friccional:
0.57
U R
*
= 3.3 h (2.26)
u ds
con d s = D90 .
20
1/ 6
Rh

d
S t = s 0.57 (2.27)
R
3.3 h
ds
Rh
ds
El estudio experimental realizado por Ayala y Oyarce consisti en la produccin de un

modelo fsico a escala 1:30 de un prototipo representativo de cauces cordilleranos Andinos. El
prototipo utilizado fue un tramo del ro Mapocho de 300 metros de longitud, ubicado en un sector
denominado los Almendros, Chile, con una pendiente media de 0.0154. Luego de construido el
modelo se realizaron experiencias sin arrastre de sedimentos, sin arrastre de sedimentos con lecho
previamente acorazado y con arrastre generalizado. Con los resultados de estas experiencias
propusieron la expresin (2.26), basada en correlacin de datos, la cual da cuenta de las
experiencias anteriormente sealadas. Es interesante notar que a diferencia de las ecuaciones
propuestas anteriormente por diferentes autores, la ecuacin (2.26) es del tipo monomial y no
logartmica. Aguirre y Fuentes (1987) (citado en Oyarce, 1993) postulan que para lechos macro
rugosos, la ecuacin (2.18) que representa la ley de distribucin de velocidades logartmica de
Prandtl (1926) von Karman (1930) (ambos citados en Keulegan, 1938) deja de ser vlida a una
determinada profundidad debido a la presencia de grandes sedimentos en el fondo. Estos
generaran un perfil de velocidades constante en sus cercanas debido a las estelas de flujo que se
produce en torno a estos.
La experimentacin de Ayala y Oyarce mediante la cual se calibr la ecuacin (2.26),

comprende los siguientes rangos correspondientes al prototipo: inversos de rugosidades relativas
entre [2.05 - 7.23], caudales entre [11.1 - 204] m3/s y dimetros caractersticos entre [0.04 - 0.47]
m. Los requisitos que debe cumplir, simultneamente, la ecuacin (2.26) para ser vlida son:
lechos con granulometras gruesas (gravas), flujos turbulentos plenamente desarrollados, lechos
21
con pared hidrodinmicamente rugosa, ausencia de formas de fondo, rgimen de escurrimiento

torrencial y fuertes pendientes tpicas de cauces cordilleranos andinos.

2.4 Resumen de Parametrizaciones Consideradas
En general, las expresiones relacionadas con la prdida de carga tienen un comportamiento

logartmico comn, el cual est basado en la ley de distribucin de velocidades logartmicas,
salvo la expresin propuesta por Ayala y Oyarce (1993), la cual es del tipo monomial. Adems
son, en general, vlidas para granulometras gruesas y extendidas, salvo Keulegan (1938) que es
vlida para arenas.
Todas las expresiones comparten los siguientes requisitos para ser vlidas: flujos
turbulentos plenamente desarrollados, lechos con pared hidrodinmicamente rugosa y ausencia de
formas de fondo, esto ltimo debido a que las ecuaciones solo consideran el efecto en la prdida
de carga producida por el sedimento del lecho. En estricto rigor, el efecto del sedimento es
atribuible a la coraza de los lechos, sin embargo, una medida indirecta del sedimento acorazado
22
es la medicin de un dimetro caracterstico del sedimento del lecho original, lo cual es

expresado por los distintos autores como D90 o D84 , en particular, Ayala y Oyarce (1993)
determinaron que D50 coraza = D90 . Se repite que tanto D90 como D84 deben estar asociados a la
granulometra integral del lecho, mientras que D50 coraza debe estar asociado a la granulometra de
la coraza.
Si bien, los diferentes autores sugieren amplios rangos de validez para sus
parametrizaciones del nmero de Strickler, en trminos del inverso de la rugosidad relativa, no
existe un consenso en cuanto al valor del nmero de Strickler en el lmite entre lechos rugosos y
macro rugosos, lo que se aprecia en la siguiente tabla.
Tabla 2.1: Valores Correspondientes a Lmites entre Lechos Rugosos y Macro Rugosos
R
Autor S t h = 10
ds
Keulegan (1930) 0.123
Limerinos (1970) 0.161
Parker y Peterson (1980) 0.149
Ayala y Oyarce (1993) 0.120
En la tabla (2.1) se puede apreciar que tanto Limerinos (1970) como Parker y Peterson
(1980) brindan valores superiores a 0.12 para la transicin entre lechos rugosos y macro rugosos,
mientras que Keulegan (1930) converge a dicho valor. Entonces en funcin de los antecedentes
recopilados en este trabajo de ttulo, se adopt como valor de transicin entre lechos 0.12 y se
mantiene dicho valor constante para lechos rugosos. En este contexto, las figuras: (2.2), (2.3),
(2.4) y (2.5) que representan el nmero de Strickler en funcin del inverso de la rugosidad
relativa, fueron modificadas para lograr una transicin suave entre lechos macro rugosos y
rugosos propuesta en este trabajo de ttulo. Esto se describe en la seccin (3.1) mediante
modificaciones lineales y parablicas. El sentido de modificar dichas parametrizaciones y
extender su aplicabilidad fuera de su rango de validez recae en la transicin que deben cumplir
segn la recopilacin realizada en este trabajado de ttulo.
23
En la figura (2.6) se resumen grficamente las parametrizaciones del nmero de Strickler

revisadas en este captulo. En la tabla (2.2) se muestra un resumen con estas parametrizaciones
indicando su autor, expresin y rango de validez.
Figura 2.6: Resumen Parametrizaciones del Nmero de Strickler Originales Segn Autor.
Tabla 2.2: Resumen Parametrizaciones del Nmero de Strickler.
Rango de Validez
Autor Expresin Rh Supuesto
( )
ds
Granulometras
Strickler (1928) 0.12 >10
Gruesas (gravas)
1/ 6
Rh

Keulegan (1930) St = ds [6.3 - 1,030]
Granulometras
Finas (arenas)
R
2.5 ln12 h
ds
24
Tabla 2.2 (continuacin): Resumen Parametrizaciones del Nmero de Strickler.
Rango de Validez
Autor Expresin Rh Supuesto
( )
ds
1/ 6
Rh

Limerinos (1970) St = ds [1.06 - 68.05]
Granulometras
Gruesas (gravas)
R
2.5 ln 3.8 h
ds
1/ 6
Rh

Parker y Peterson
St = ds [1.06 - 68.05]
Granulometras
(1980) Gruesas (gravas)
R
2.46 ln 5.5 h
ds
1/ 6
Rh

d
Ayala y Oyarce (1993) S t = s 0.57 [2.05 - 7.23]
Granulometras
Gruesas (gravas)
Rh
3.3
ds
25
Patricio Pastenes Araos Captulo 3: Implementacin MR-HR
Captulo 3. Implementacin de la Herramienta: MR-HR
Este captulo surge de la necesidad de comprender los fundamentos de MR-HR como su

interaccin con HEC-RAS. Para esto se explican las modificaciones que se realizaron a las
parametrizaciones consideradas en la seccin (2.3). Luego, se mencionan tanto los atributos como
los alcances, desde el punto de vista hidrulico, a los que est sujeto MR-HR. Para finalmente,
explicar el esquema de operacin de este ltimo programa y como este interacta con HEC-RAS,
junto a los resultados proporcionados por MR-HR.
3.1 Ajuste de Parametrizaciones
Con el propsito de facilitar la presentacin de las deducciones realizadas en este captulo,

a continuacin se muestra la nomenclatura utilizada:
f (x) : corresponde a la parametrizacin del nmero de Strickler de Keulegan (1938), Limerinos

(1970), Parker y Peterson (1980) y Ayala y Oyarce (1993).
g (x) : corresponde a la parametrizacin lineal del nmero de Strickler aqu propuesta.
h(x) : corresponde a la parametrizacin parablica del nmero de Strickler aqu propuesta.
x : corresponde al inverso de la rugosidad relativa.
Zi : corresponde a la cota de fondo del lecho de la seccin i-sima, medida con respecto a un
datum de referencia. [m]
Yi : corresponde a la altura de escurrimiento de la seccin i-sima, medida con respecto a la

cota de fondo del lecho. [m]
Vi : corresponde a la velocidad media de la seccin i-sima. [m s ]
ai : corresponde al coeficiente de Coriolis asociado a la seccin i-sima.
Se propone realizar ajustes a las parametrizaciones del nmero de Strickler, expuestas en la

seccin (2.3), para lograr una transicin suave entre los escurrimientos desarrollados en lechos
26
macro rugosos y rugosos, lo cual se justific en la seccin (2.4). Esto se realiz haciendo tanto un
ajuste lineal como otro parablico, apoyado en la generalizacin de la forma de las ecuaciones
de: Keulegan (1938), Limerinos (1970) y Parker y Peterson (1980), dado que siguen una ley
logartmica2, es decir, la forma genrica de estas ecuaciones es:
1 x1 / 6
f ( x ) = (3.1)
k1 ln (k 2 x )
3.1.1 Ajuste Lineal
Se proyect una lnea recta, tal que en un punto x * la recta compartiera tanto el mismo
valor de la parametrizacin, como la misma derivada, sujeto a que en el fin de la zona macro
rugosa x = Rh ds = 10 se cumple que la recta debe tener el valor del Strickler constante e igual a
0.12. Grficamente, se tiene lo siguiente:
Figura 3.1: Representacin Grfica Ajuste Lineal

2
Notar que Ayala y Oyarce obedece a una ley de Potencia, por lo que no se incluye en este anlisis.
27
El planteamiento formal, es el que se muestra a continuacin:
Sea g (x) la recta que representa la parametrizacin lineal del nmero de Strickler:
df
g ( x) = x+b (3.2)
dx x = x*
Con las siguientes condiciones de borde:
df
g (10) = 0.12 10 + b = 0.12 (3.3)
dx x = x*
df
g ( x*) = f ( x*) x * +b = f ( x*) (3.4)
dx x = x*
Restando la ecuacin (3.3) a la ecuacin (3.4), resulta:
df
(10 x *) = 0.12 f ( x *) (3.5)
dx x = x*
Por otra parte, derivando la ecuacin (3.1), se obtiene:
df 1 1 1 1
= 5 / 6 5 / 6 (3.6)
dx 6 k1 x ln (k 2 x ) k1 x ln (k 2 x )2
Finalmente, reemplazando la ecuacin (3.6) y (3.1) en la ecuacin (3.5), resulta:
1 1 1 1 1/ 6
(10 x *) = 0.12 1 x *
6 k x *5 / 6 ln (k x *) k x *5 / 6 ln (k x *)2 k ln (k x *)
1 2 1 2 1 2
(3.7)
Resolviendo la ecuacin (3.7), se encuentra el punto x * , donde nace la recta que comparte
tanto la misma derivada como el mismo valor de la parametrizacin. Luego de conocido el punto
28
x*, basta con reemplazar en la ecuacin (3.3) o en (3.4) para conocer el coeficiente de corte de la
lnea recta con el eje de las ordenadas, es decir, el coeficiente b .
3.1.2 Ajuste Parablico
Con la finalidad de realizar una transicin suave en las cercanas del inverso de la
rugosidad relativa (Rh ds ) igual a 10, se propuso una transicin parablica3. Est transicin
ocurre en el intervalo [8,12] en el dominio del inverso de la rugosidad relativa. La parbola debe
ser tal que, en su extremo inferior ( x = 8 ), debe tener el mismo valor de la recta proyectada
anteriormente. Adems en su extremo superior ( x = 12 ) debe tener un valor igual a 0.12 y su
derivada debe ser nula. Grficamente, la idea es la siguiente:
Figura 3.2: Representacin Grfica Ajuste Parablico

3
La seleccin de los puntos del dominio de la parbola x [8,12] fueron realizados en forma arbitraria para lograr una
transicin suave y no sobre estimar en demasa los valores del nmero de Strickler en lechos rugosos, es decir, para
inversos de rugosidades relativas mayores que 10.
29
El planteamiento formal, es el que se muestra a continuacin:
Sea h(x) la parbola:
h( x ) = a x 2 + b x + c (3.8)
Las condiciones de borde a las que est sujeta la parbola, son las siguientes:
dh
=0 (3.9)
dx x = x ''
h(x ' ) = g (x ' ) a (x ' ) + b (x ' ) + c = g (x ' )

2
(3.10)
h(x '' ) = g (10) a (x '' ) + b (x '' ) + c = g (x '' )

2
(3.11)
Donde:
x' = 8 x '' = 12 g (10 ) = 0.12
Derivando la ecuacin (3.8) resulta:
dh
= 2a x +b (3.12)
dx
Luego, reemplazando la ecuacin (3.12) en la ecuacin (3.9), resulta:
b = 2 a x ' ' (3.13)
Restando la ecuacin (3.10) a la ecuacin (3.11), resulta:
( 2 2
)
g (x ' ) g (10) = a (x ' ) (x '' ) + b (x ' x '' ) (3.14)
30
Reemplazando la ecuacin (3.13) en la ecuacin (3.14), resulta:
a=
( )
g x ' g (10 )
(3.15)
((x ) (x ) ) 2 x (x x )
' 2 '' 2 '' ' ''
Donde:
x' = 8 g (x ' ) = Es un valor conocido de la parametrizacin lineal anterior.
x '' = 12 g (10 ) = 0.12
Con estas ecuaciones, el problema ya est resuelto, ya que una vez conocido a se
reemplaza en la ecuacin (3.13) y se conoce b . Finalmente para conocer c , basta con reemplazar
a y b en alguna condicin de borde, es decir, en las ecuaciones (3.10) o (3.11).
3.1.3 Ajuste Ayala y Oyarce
Dado que Ayala y Oyarce no obedecen a una ley logartmica, sino ms bien a una ley de
potencia, slo se ajust en trminos parablicos, utilizando las mismas ecuaciones anteriores. La
nica modificacin es en la ecuacin (3.10), que modificada resulta:
h(x ' ) = f (x ' ) a (x ' ) + b (x ' ) + c = f (x ' )

2
(3.16)
Ntese que se cambi la funcin g (x ' ) por f (x ' ) , dado que no es necesario un ajuste lineal
y, por ende, slo se suaviza mediante una parbola en las cercanas de Rh ds = 10 . Para esto se
utilizan las mismas ecuaciones anteriores, salvo que se reemplaza la condicin de borde, definida
por la ecuacin (3.10) por la ecuacin (3.16).
31
3.1.4 Resumen Coeficientes Lineales y Parablicos
En la siguiente tabla, se resumen los valores de los ajustes lineales a las parametrizaciones
sealadas en la seccin (3.1.1):
Tabla 3.1: Resumen de parmetros de ajustes lineales

Parametrizacin k1 k2 x* df b
dx x = x*
Keulegan 2.50 12.00 5.800 -0.001504010 0.134967918

Limerinos 2.50 3.80 3.691 -0.010814549 0.228100000
Parker y Peterson 2.46 5.50 3.675 -0.007585932 0.195900000
Ayala y Oyarce - - - - -
Fuente: Elaboracin Propia.
En la siguiente tabla, se resumen los valores de los ajustes parablicos a las

parametrizaciones sealadas en la seccin (3.1.2):
Tabla 3.2: Resumen de parmetros de ajustes parablicos

Parametrizacin a b c
Keulegan 0.000183490 -0.004403752 0.146422513
Limerinos 0.001348975 -0.032375411 0.314252469
Parker y Peterson 0.000950784 -0.022818811 0.256912867
Ayala y Oyarce 0.000686906 -0.016485734 0.218914406
Las tablas (3.1) y (3.2) representan todos esos decimales debido a que se deja constancia de
la forma de operacin de MR-HR, con cuatro decimales es suficiente ya que se consigue una
variacin en el cuarto decimal del coeficiente de Manning, lo cual es imperceptible en la
variacin de alturas de escurrimiento en metros. A continuacin se muestran grficamente los
ajustes tanto lineales como parablicos realizados a todas las parametrizaciones descritas en las
tablas anteriores:
32
Figura 3.3: Ajustes Parametrizacin Keulegan

Figura 3.4: Ajustes Parametrizacin Limerinos

33
Figura 3.5: Ajustes Parametrizacin Parker y Peterson

Figura 3.6: Ajustes Parametrizacin Ayala y Oyarce

34
3.1.5 Ajustes Inverso de Rugosidad Relativa Mnimo
Es importante destacar que las parametrizaciones de Keulegan (1938), Limerinos (1970) y

Parker y Peterson (1980) obedecen a una ley logartmica, la cual, eventualmente, podra
convertirse en nmero negativo si el argumento del logaritmo natural es inferior a 1 (ntese que
para Rh ds << 1 , el problema carece de sentido desde el punto de vista ingenieril). Por lo tanto,
se realiz un ltimo ajuste a las parametrizaciones y para unificar criterio, tambin se realiz para
Ayala y Oyarce (1993), este ajuste consiste en que para todos los inversos de rugosidades
relativas pertenecientes al intervalo (0,1], cada parametrizacin tiene un nico valor constante a
lo largo de todo el intervalo, el cual se obtiene de evaluar la parametrizacin particular en el
inverso de rugosidad relativa igual a 1. En la tabla (3.3) se muestran los resultados obtenidos. El
sentido de forzar el funcionamiento de las parametrizaciones fuera de su rango de validez de
inversos de rugosidades relativas menores a uno, es netamente numrico, es decir, impedir que
los logaritmos naturales se conviertan en nmeros negativos por un uso inadecuado de las
parametrizaciones. Por lo tanto, en estos casos los resultados entregados por el software deben ser
interpretados con criterio.
Tabla 3.3: Resumen de Nmeros de Strickler vlidos para todo el intervalo de inversos de
rugosidades relativas (0,1]
R
S t h = 1
Parametrizacin ds
Keulegan 0.160971842
Limerinos 0.299625229
Parker y Peterson 0.238454037
Ayala y Oyarce 0.303030303
35
3.1.6 Resumen de Parametrizaciones Ajustadas
A continuacin, se muestra grficamente el resumen de parametrizaciones ajustadas, segn

los criterios anteriormente descritos. Es importante destacar que bajo estas parametrizaciones, la
herramienta MR-HR realizar sus clculos:
Figura 3.7: Resumen Parametrizaciones Ajustadas

Las parametrizaciones del nmero de Strickler que se observan en la figura (3.7) se

obtuvieron a partir de las tablas (3.1), (3.2), (3.3) y las ecuaciones (2.17), (2.20), (2.25), (2.27).
3.2 Atributos MR-HR
En el diseo e implementacin de MR-HR se pens en un software que fuese capaz de

calcular los coeficientes de Manning en escurrimientos sobre lechos macro rugosos, mediante el
clculo iterativo de nmeros de Strickler y ejes hidrulicos. Para esto, MR-HR interacta con
HEC-RAS, dado que HEC-RAS tiene una vasta utilizacin en el mundo ingenieril, con un
36
respaldo de aos de desarrollo y validacin por parte de ingenieros de todo el mundo, incluido
Chile. Por lo tanto, los atributos considerados en el diseo e implementacin de MR-HR fueron:
Clculo iterativo: En base a iteraciones sucesivas del nmero de Strickler, MR-HR logra
la convergencia de los coeficientes de Manning y, consecuentemente, el desarrollo de ejes
hidrulicos.
Parametrizaciones del Nmero de Strickler: como se vio en la seccin (3.1.6), MR-HR

incorpora diferentes parametrizaciones del nmero de Strickler. Esto se realiz con la
finalidad de brindar un abanico de posibilidades en cuanto al clculo de ejes hidrulicos.
Perfiles de Flujo: Si bien, MR-HR se basa en iteraciones sucesivas del nmero de

Strickler, tambin es posible realizar iteraciones en trminos del caudal. Esto ltimo se
implement con la finalidad de conseguir ejes hidrulicos para diferentes perfiles de flujo.
Resultados Grficos: Con la finalidad de facilitar la interpretacin de los resultados

brindados por MR-HR, se incorporaron diferentes opciones grficas, las cuales pueden ser
consultadas en detalle en el Anexo (A.4.5).
Acoplamiento HEC-RAS: El acoplamiento con HEC-RAS se debe a que fue necesario un

software que resolviera ejes hidrulicos para coeficientes de Manning dados.
Las opciones de seleccin consideras en la implementacin de MR-HR, fueron las

siguientes:
Seleccin de Parametrizacin: Es posible la seleccin de la parametrizacin del nmero

de Strickler deseada, entre las opciones presentadas en la seccin (3.1.6).
Dimetro Caracterstico: Es posible ingresar cualquier magnitud de dimetro

caracterstico deseado, asociado a cada perfil transversal. En caso de no poseer la
granulometra completa a lo largo del cauce, el programa interpolar de forma lineal.
37
Nmero de Iteraciones: Es posible ingresar el nmero de iteraciones que se desee,

referentes al clculo de nmeros de Strickler y ejes hidrulicos, para lograr la
convergencia de los coeficientes de Manning.
Rango de Caudales: Es posible ingresar cualquier rango de caudales deseado, solamente

aguas arriba, para lograr diferentes perfiles de flujo y analizar los diferentes ejes
hidrulicos proporcionados por MR-HR. Adems, si se desea, se puede hacer funcionar
MR-HR para solo un caudal aguas arriba.
3.3 Relaciones entre MR-HR y HEC-RAS
Las limitaciones existentes en HEC-RAS son ntegramente traspasadas a MR-HR, dado que
ste se apoya en los clculos brindados por HEC-RAS.
Hasta ahora slo se ha visto que las ecuaciones expresadas en la seccin (2.1) son vlidas
para flujos permanentes y uniformes, los cuales no son tpicamente encontrados en los problemas
de la ingeniera hidrulica. Por el contrario, los problemas comnmente encontrados son flujos
permanentes y gradualmente variados, es decir, la altura de escurrimiento vara paulatinamente a
lo largo de los cauces. HEC-RAS aborda esta clase de problemas mediante la resolucin de la
ecuacin de la energa, la cual es la siguiente:
2 2
a V a V
Z 2 + Y2 + 2 2 = Z1 + Y1 + 1 1 + he (3.17)
2 g 2 g
Donde he representa las prdidas de energa tanto por efectos friccinales como expansiones o
contracciones del flujo.
2 2
a V a V
he = L S f + C 2 2 1 1 (3.18)
2 g 2 g
38
Donde S f representa el gradiente o pendiente de friccin, el cual debe ser representativo entre
las dos secciones transversales consideradas.
Para el clculo de S f se utiliza la ecuacin de Manning (2.10). Adems debe ser
representativo del tramo, es decir, debe considerar las pendientes de ambas secciones
transversales que conforman el tramo estudiado. Si bien la ecuacin de Manning es vlida para
problemas con flujos uniformes y permanentes donde la pendiente de fondo es idntica tanto a la
pendiente de energa como a la pendiente de la superficie libre del agua, se considera que es
localmente aplicable en escurrimientos gradualmente variados al utilizar S f en lugar de S en la
ecuacin (2.10). El algoritmo de resolucin que utiliza HEC-RAS es el mtodo del paso estndar
en vez del mtodo del paso directo. El primer mtodo obtiene una altura de escurrimiento dadas
la altura de escurrimiento de la seccin anterior y la distancia entre secciones. El segundo obtiene
la distancia entre dos secciones en base a alturas predeterminadas. Esta implementacin permite
enfrentar problemas gradualmente variados en cauces no prismticos donde varan la geometra
de las secciones transversales, las pendientes de fondo, los coeficientes de Manning en el eje
longitudinal, etc.
Es importante notar que en la ecuacin (3.17) est expresada en forma explicita la ley de
distribucin hidrosttica de presiones, representada por la carga piezomtrica. Este supuesto es
vlido para pendientes menores o iguales que el 10 %.
La ampliacin de flujos uniformes a gradualmente variados en canales no prismticos tiene

un impacto considerable, debido a que el primer tipo de flujo es inmensamente restrictivo, dado
que no pueden variar las secciones transversales y tampoco ninguna propiedad del flujo en
trminos espaciales. En cambio, en el segundo tipo de flujo se admite que vare la altura de
escurrimiento y por ende, se pueden abordar problemas ms generales donde cambien las
velocidades, geometras de las secciones transversales, pendientes de fondo, etc. Estos tipos de
problemas son recurrentes en el quehacer cotidiano del ingeniero hidrulico.
Otro aspecto interesante de HEC-RAS es que en nodos (puntos donde confluyen o divergen
2 o ms ros) utiliza igualdad de energa. Adems, y en caso que se requiera, decide cual de los
39
dos tramos controla mediante la aplicacin de la funcin de momentum. El que tenga mayor
momentum o fuerza especfica ser la condicin de borde que definir la cota de la superficie
libre en el nodo. Este mtodo es el que HEC-RAS utiliza por defecto y con el cual se trabaja en
MR-HR.
3.4 Esquema de Funcionamiento
El funcionamiento de MR-HR est ntimamente ligado con HEC-RAS 4.1.0 debido a que
los clculos de ejes hidrulicos se desarrollan a travs del ltimo programa. Sin embargo, MR-
HR realiza una serie de clculos en base a los resultados de HEC-RAS, y adems realiza una
estructura lgica mediante la cual decide si terminar o no la ejecucin.
Los pasos que realiza la herramienta son los siguientes:
1. Introduce coeficientes de Manning a todas las secciones transversales, inicialmente todos igual
a 0.030 (valor propuesto en Chow (1994) para cauces naturales de granulometras gravosas) y
adems prepara el primer perfil de flujo, es decir, el caudal aguas arriba para el cual se realizarn
las iteraciones.
2. Calcula mediante HEC-RAS el respectivo eje hidrulico, y consecuentemente los radios

hidrulicos relativos a todas las secciones transversales.
3. Gracias a los resultados del punto (2), recalcula los coeficientes de Manning, de cada seccin
transversal en particular, a travs de la parametrizacin del nmero de Strickler seleccionada.
4. Evala si quedan iteraciones pendientes por realizar, previamente definidas por el usuario. Esto
ltimo es equivalente a preguntarse Convergieron los coeficientes de Manning, es decir, se
realizaron todas las iteraciones ingresadas, previamente, por el usuario relativas a los coeficientes
de Manning?
4.1 No, entonces introduce los Coeficientes de Manning recalculados y vuelve a iterar desde el
punto (2).
40
4.2 S, entonces evala Quedan perfiles restantes, es decir, otros caudales para los cuales el
usuario quiera calcular nuevamente los ejes hidrulicos?.
4.2.1. No, entonces fin de la ejecucin.
4.2.2. S, entonces vuelve al punto (1) con el nuevo perfil y realiza todos los pasos.
MR-HR HEC RAS
INTRODUCE COEFICIENTES
DE MANNING INICIAL Y
CAUDAL ESPECIFICO
CONVERGIERON LOS RECALCULA COEFICIENTES CALCULA EJE

COEFICIENTES DE DE MANNING A TRAVES HIDRAULICO Y RADIOS
MANNING? DEL NDE STRICKLER HIDRAULICOS
NO INTRODUCE
COEFICIENTES DE
MANNING RECALCULADOS
SI
RESULTADOS FINALES
SE INCREMENTA EL
SI
CAUDAL?
NO
FIN DE EJECUCION
Figura 3.8: Esquema de Funcionamiento MR-HR.

Finalmente, es posible apreciar todos los pasos anteriormente descritos, resumidos en la

figura (3.8).
3.5 Resultados Proporcionados
El principal aporte de MR-HR es que extiende la aplicacin de la ecuacin de Manning a

un rango para el cual no fue deducida, pero que se utiliza dada la ausencia de un mtodo de
clculo. Este aporte lo refleja mediante el clculo de coeficientes de Manning, a travs de las
41
diferentes parametrizaciones del nmero de Strickler. En la figura (3.9) se aprecia la ventana

principal del programa.
Figura 3.9: Ventana Principal MR-HR.

En el captulo Anexo (A.4.5) se muestra la forma en que se presentan los resultados MR-
HR y adems, las representaciones grficas que la herramienta proporciona.
42
Patricio Pastenes Araos Captulo 4: Aplicacin al Ro Baker
Captulo 4. Aplicacin al ro Baker.
En el presente captulo se desarroll una aplicacin, a travs de la herramienta MR-HR, a

un tramo del Ro Baker, ubicado en la XI regin de Aysn del General Carlos Ibez del Campo.
En una primera etapa se presentan antecedentes generales del ro Baker y su cuenca,

basados en su totalidad en el Estudio de Impacto Ambiental del Proyecto Hidroelctrico Aysn
(Centrales Hidroelctricas de Aysn S.A., 2008). Luego se describe un tramo del ro Baker, que
fue considerado para la ejecucin de MR-HR junto con su caracterizacin granulomtrica. Esto
ltimo se materializ basado en un estudio realizado para HidroAysn por el Grupo de
Investigacin en Procesos de Transporte y Flujos Ambientales de la Universidad de Chile
(Universidad de Chile, 2010). Finalmente, se tomaron diversos caudales de inters para la
realizacin de pruebas con MR-HR.
Es importante destacar que el presente captulo no tiene como objetivo definir que
parametrizacin es mejor para el ro Baker, ni tampoco ahondar en los aspectos particulares del
escurrimiento en el ro, sino que el objetivo principal es la implementacin de la herramienta
computacional aqu diseada a un caso real. Es por esto que la aplicacin, aqu desarrollada, est
orientada a ver las diferencias en los ejes hidrulicos calculados usando las diferentes
parametrizaciones del nmero de Strickler.
4.1 Antecedentes Generales Ro Baker
El Ro Baker drena una superficie de 26.487 km2, de los cuales 17.159 km2 estn en
territorio nacional. Su caudal medio est entre 1.200 y 1.100 m3/s, siendo el ro ms caudaloso
del pas. Tiene su origen en el lago Bertrand, ubicado en la cota 203 m.s.n.m., que se encuentra
aguas abajo del lago General Carrera, cuerpo lacustre de 1.047,5 km2, el ms extenso del pas y
de carcter internacional, llamado Buenos Aires en el lado argentino. Los principales ros
aportantes al lago General Carrera son: ro Ibez, ro Avellanos, ro Murta, ro Delta o Leones y
ro Jeinemeni, todos de orden 4 a excepcin del ro Ibez que es de orden de jerarqua 5.
43
En su recorrido de 175 km hasta desembocar en las inmediaciones de Tortel, el ro Baker

recibe numerosos e importantes afluentes tales como: Chacabuco, Cochrane, Del Salto y De Los
adis por el oriente; Nef, De La Colonia y Ventisquero por el occidente.
Figura 4.1: Red de Drenaje Subcuenca del Ro Baker.

Fuente: Centrales Hidroelctricas de Aysn S.A. (2008)
Las precipitaciones en la regin de Aysn aumentan en forma considerable de este a oeste.

En Chile Chico, en el lmite con Argentina, la precipitacin media anual es del orden de 240 mm,
en Puerto Gaudal de 790 mm y en Puerto Bertrand, en el desage del lago, de 1330 mm.
El rgimen del ro Baker es nivo-pluvial con flujos mximos en pocas de deshielo. En su

nacimiento, el ro Baker es muy uniforme debido al efecto regulador del lago General Carrera y
presenta aguas cristalinas debido al efecto decantador del mismo lago. Los primeros afluentes son
los ros Nef y Chacabuco, que generan un aumento en la turbiedad sobretodo por los aportes del
Ro de la Colonia que, por provenir directamente del derretimiento de nieve y hielos de los
44
glaciares de Campo de Hielo Norte, presentan un gran arrastre de slidos suspendidos y de

arrastre de fondo.
En general se distinguen 3 tipos de zonas morfolgicas en el ro Baker, las cuales son:
Zonas tipo cascada: se presentan situaciones aguas arriba de la confluencia del el ro Nef
con el ro baker, aguas arriba de la junta entre los ros Baker y Chacabuco y aguas arriba
de la futura central Baker 2, la cual estar emplazada en la angostura El Saltn a unos 2
km aguas arriba de la confluencia de los ros Baker y El Saltn. Este tipo de forma se
caracteriza por alta pendiente de terreno, sustrato grueso tipo rocas que obstruyen el
escurrimiento formando zonas de aguas blancas.
Regiones con canal recto con ausencia de rpidos y remansos: se ubica principalmente
aguas abajo de la confluencia entre los ros Baker y Chacabuco y algunos sitios
intermedios. Esta regin presenta una variabilidad importante en las condiciones del
escurrimiento en el eje longitudinal, dando origen a sectores de altas velocidades (rpidos
ubicados aguas abajo de secciones con escurrimientos crtico) y zonas de bajas
velocidades ubicadas aguas arriba de la seccin crtica (remansos).
Tramos trenzados: esta clase de escurrimiento se presenta en amplias regiones del ro

Baker y sus afluentes, y se generan por la disminucin de la pendiente de terreno, hecho
que favorece el desarrollo lateral, permitiendo as la existencia de ms de un brazo del ro.
Adems, dados los procesos de erosin de las riberas de los cauces, la variabilidad
temporal de la forma del cauce es amplia, resultando que el trenzado del ro vara de
forma en el transcurso de un ao. Desde el punto de vista hidrulico, el rgimen en estas
zonas es subcrtico (ro), caracterizado por bajas velocidades y altas alturas de
escurrimiento. Este tipo de escurrimiento se verifica en el tramo del ro Baker
comprendido desde 8 Km. aguas arriba de la confluencia del ro Del Salto hasta dos
kilmetros aguas abajo de la confluencia con el ro De Los adis.
45
Figura 4.2: Ubicacin de las Centrales Baker 1 y 2 en el Eje Longitudinal y la Distribucin

Longitudinal de los Tipos Morfolgicos Presentes en el Ro Baker.
Fuente: CEA (2006) citado en Centrales Hidroelctricas de Aysn S.A. (2008)
4.2 Tramo Particular Estudiado
La aplicacin se desarroll en el ro Baker aguas abajo del futuro emplazamiento de Baker

1 (central hidroelctrica que construir HidroAysn), hasta su desembocadura en caleta Tortel.
Este tramo tiene una longitud de 148 km y est caracterizado por 740 perfiles transversales, los
cuales se detallan en el Anexo (B.3). En la figura (4.3) se puede apreciar el tramo particular
estudiado.
La central hidroelctrica Baker 1 estar emplazada en la angostura Chacabuco, a unos 1000

m aguas arriba de la confluencia entre los ros Baker y Chacabuco. Baker 2 estar emplazada en
la angostura El Saltn, unos 2 km aguas arriba de la confluencia de los ros baker y El Saltn.
46
Figura 4.3 : Tramo Estudiado en la Presente Aplicacin.

Fuente: Centrales Hidroelctricas de Aysn S.A. (2008)
Notar que la figura (4.3) es referencial dado que ninguna de las dos centrales ha sido
construida y las superficies de inundacin, marcadas en negro, son las superficies estimadas que
tendrn ambas centrales. Por lo tanto, para efectos de la aplicacin de MR-HR no existen zonas
de inundacin, y slo se considera el cauce natural desde Baker 1 a la desembocadura en caleta
Tortel.
4.3 Caracterizacin Granulomtrica
El trabajo realizado por Universidad de Chile (2010) consisti en mejorar la caracterizacin

de la granulometra de los ros Baker y Pascua. Para esto, se realizaron 5 visitas a terreno
comprendidas entre los aos 2006 al 2010. Los sitios muestreados se resumen en la siguiente
tabla:
47
Tabla 4.1: Ubicacin de Lugares de Muestreos Granulomtricos, Ro Baker.
Coordenadas UTM
Tipo de
Sitio Ubicacin (Datum WGS84)
muestreo
N E
Ro Baker en islote bajo
SLB1 Integral 4.788.440 665.288
sector desage L. Bertrand
Ro Baker en tramo
SLB2 Integral 4.785.155 667.109
Bertrand - Nef
Ro Baker antes junta Ro
SLB3 Integral 4.779.055 667.877
Nef
Ro Baker bajo Angostura
LB-S1 Integral 4.777.528 673.999
Chacabuco
Ro Baker en sector Balsa
LB-S2 Integral 4.769.355 679.328
Baker
Ro Baker en isla sector
RB2A Integral 4.762.407 665.091
Valle Grande
Ro Baker bajo junta Ro de
RB2B Superficial 4.749.013 659.955
la Colonia
Ro Baker antes junta Ro de
RB2C Superficial 4.743.142 654.808
Los adis
LB-S4 Ro Baker bajo junta Vargas Integral 4.713.594 634.305
Fuente: Universidad de Chile (2010)
El muestreo superficial corresponde a la granulometra de la coraza del lecho, mientras que

el muestro integral representa la granulometra integral del lecho sin incluir su coraza. Slo deben
considerarse de la tabla (4.1) los sitios: LB-S1, LB-S2, RB2A, RB2B, RB2C, LB-S4, debido a
que estas ubicaciones se encuentran aguas abajo de la futura central hidroelctrica Baker 1
(angostura Chacabuco). A continuacin se muestran los esquemas de ubicacin de dichos sitios.
48
Figura 4.4: Esquema de ubicacin de calicatas Ro Baker (inicio) y Ro Nef.

Figura 4.5: Esquema de ubicacin de calicatas Ro Baker (Valle Grande) y Ro de la Colonia.

49
Figura 4.6: Esquema de ubicacin de calicatas Ro Baker (bajo adis hasta Ro Vargas) y Ro de
los adis.
Las campaas de muestreo asociadas a los distintos sitios, se muestran a continuacin:
Tabla 4.2: Campaas de terreno de muestreo granulomtrico.
Fechas Sitios muestreados
21-25 Nov. 2006 LB-S2, LB-S3, LB-S4, LB-S5

19-21 Ene. 2007 LB-S1, CB1, CB2, CB3
RB2A, RB2B, RB2C, RNAD, RC4,
12-14 Nov. 2009
RC5G, RC6, RN1, RN2, RN3
16-17 Dic. 2009 SLB1, SLB2, SLB2
RP1, RP2, RP3, RP3a, RP4, RP5,
13-15 Feb. 2010
RP6, RP7, RP8
50
Finalmente, los resultados que proporcion el trabajo granulomtrico fueron:
Tabla 4.3: Tamaos caractersticos del sustrato, Ro Baker (mm).

SLB1 SLB2 SLB3 LB-S1 LB-S2 RB2A LB-S4
D90 131.2 20.1 62.3 72.3 45.0 22.2 18.4
D84 116.4 12.7 51.0 57.1 36.2 17.5 13.0
D50 22.3 1.0 18.1 9.8 15.4 5.1 3.9
D16 2.1 0.2 0.6 0.4 1.7 0.4 0.7
En la tabla (4.3) se excluyen los resultados de los sitios RB2B y RB2C debido a que en la
confeccin de esa tabla se buscaba caracterizar el sustrato y esos puntos fueron muestreados
superficialmente. Sin embargo, tambin es posible agregar informacin granulomtrica de la
coraza, en particular el dimetro caracterstico medio de la coraza ya que es vlido ligar ste con
el dimetro caracterstico 90% y 84% de la granulometra integral del lecho. En este anlisis, en
particular, no se incluyeron dichos sitios de muestreo, pero se podran haber incluido sin
problema. A continuacin se adjunta la curva que representa la variacin longitudinal del
sedimento del sustrato en el Ro Baker.
Figura 4.7: Variacin longitudinal de los tamaos caractersticos de los sedimentos del lecho del
Ro Baker (Sustrato)
51
De la figura (4.7) se deben extraer los datos representativos del cauce natural particular
estudiado, es decir, dado que el tramo particular en estudio nace aguas abajo de la angostura
Chacabuco, slo se consideran los muestreos desde LB-S1 en adelante, o en forma equivalente,
desde los 20 km en adelante. Adems es posible observar un cambio de tendencia en la
distribucin de los sedimentos gruesos D90 , debido a que en la estaciones LB-S1, LB-S2 y RB2A
se observa una tendencia decreciente y luego, entre las estaciones RB2A y LB-S4, se observa una
tendencia casi constante.
Para efectos del anlisis en MR-HR, se necesitan los dimetros caractersticos D90 y D84 .
En un primer anlisis se consider slo un dimetro caracterstico para todo el cauce estudiado.
De esta forma, la macro rugosidad vari en funcin del desarrollo del eje hidrulico. Luego, se
compararon los resultados en la forma ms general posible, es decir, se consider la variacin de
los dimetros caractersticos del lecho a lo largo del cauce junto con el efecto del desarrollo del
eje hidrulico en la macro rugosidad versus proponer un slo coeficiente de Manning para todo el
cauce, basado en el dimetro caracterstico representativo de todo el cauce calculado en la etapa
anterior.
En la primera etapa, se plante la siguiente metodologa para definir el dimetro

caracterstico representativo de todo el cauce:
1. Se realiz un ajuste lineal entre los puntos de muestreo LB-S1 y RB2A, es decir, desde el
km. 20 al km. 60 segn la figura (4.7). Con esto se obtuvo una representacin lineal de la
distribucin de los sedimentos entre dichos puntos.
2. Se realiz otro ajuste lineal entre los puntos de muestreo RB2A y LB-S4, es decir, desde
el km. 60 al km. 140 segn la figura (4.7). Con esto se obtuvo otra representacin lineal
de la distribucin de los sedimentos entre dichos puntos.
52
3. Con esto se procedi a calcular el valor representativo del sedimento, mediante el valor
medio ponderado (VMP):
60 140
y(x ) dx + y(x ) dx
20 60
DS (VMP ) = 140
(4.1)
dx
20
Donde:
y(x ) : representa la parametrizacin lineal de la distribucin de sedimentos, entre las

estaciones LB-S1 y RB2A. Es funcin de x dado que vara a lo largo del eje
longitudinal. (mm)
y ( x ) : representa la parametrizacin lineal de la distribucin de sedimentos, entre las

estaciones RB2A y LB-S4. Es funcin de x dado que vara a lo largo del eje
longitudinal. (mm)
x : representa el eje longitudinal, es decir, la distancia. (km)
A continuacin se muestra el resumen de clculos realizados para la obtencin de la

caracterizacin granulomtrica.
Tabla 4.4: Coeficientes de Correlacin D90 (del sustrato), Ro Baker.
Rectas: c1 x + c2
Recta c1 c2 R2
y (x) -1.2525 96.6 0.9974
y (x) -0.0487 25.1231 1.0
53
Tabla 4.5: Coeficientes de Correlacin D84 (del sustrato), Ro Baker.
Rectas: c1 x + c2
Recta c1 c2 R2
y ( x) -0.99 76.5333 0.9990
y ( x) -0.0577 20.9615 1.0
Tabla 4.6: Resumen de Valores Medios Ponderados (VMP), Ro Baker.
Dimetro Caracterstico VMP [mm]

D90 29.0
D84 22.4
Finalmente, se muestra en forma grfica los resultados de las tablas (4.4), (4.5) y (4.6).
Figura 4.8: Parametrizacin lineal de la variacin longitudinal de los tamaos caractersticos de

los sedimentos del lecho del Ro Baker (Sustrato D90)
54
Figura 4.9: Parametrizacin lineal de la variacin longitudinal de los tamaos caractersticos de

los sedimentos del lecho del Ro Baker (Sustrato D84)
4.4 Clculo de Ejes Hidrulicos Mediante MR-HR
Para el clculo de ejes hidrulicos mediante MR-HR en una primera etapa, se defini un
rango de caudales entre [1 - 50] m3/s, con un paso de 1 m3/s. Adems, se seleccionaron los
dimetros caractersticos:
Valor Medio Ponderado (VMP)
Valor Mnimo
Valor Mximo
La seleccin de estos dimetros se debe a que es posible obtener un eje hidrulico medio, a
travs del valor medio ponderado (VMP), y adems obtener un eje hidrulico mximo y otro
mnimo, representados por los dimetros caractersticos mximo y mnimo respectivamente.
55
La siguiente tabla resume los valores de los dimetros caractersticos anteriormente

sealados:
Tabla 4.7: Resumen de dimetros caractersticos considerados
D90 D84
Valor Considerado
[mm] [mm]
Valor Mnimo 18.4 13.0
Valor Medio Ponderado (VMP) 29.0 22.4
Valor Mximo 72.3 57.1
Para definir una mtrica de cul es el dimetro caracterstico que genera mayor macro
rugosidad, se emple el concepto de nmero de Strickler promedio de todo el cauce, para un
determinado caudal:
740
St(
i =1
i ,q )
St promedio = (4.2)
740
Donde:
St (i ,q ) : representa el nmero de Strickler, para una parametrizacin particular, de la seccin
transversal i-sima y para el caudal q-simo.
La ecuacin (4.2) es interesante desde el punto de vista de poder seleccionar el dimetro

caracterstico que genere un mayor grado de macro rugosidad, representado por el nmero de
Strickler promedio de todo el cauce para un determinado caudal, para poder centralizar el anlisis
en ese dimetro. Adems esta manera de cuantificar la macro rugosidad permite ver la evolucin
de sta en funcin del caudal. Con esto ltimo, tambin se podr centrar la atencin en los
caudales que generen mayor macro rugosidad.
La importancia de focalizar el anlisis en las macro rugosidades ms elevadas recae en que

si se generan muchos sectores con lechos rugosos, pierde un poco de sentido el uso de MR-HR,
56
debido a que no se necesitar realizar iteraciones en esas zonas pues el valor del nmero de
Strickler es constante.
Entonces se ejecut este problema, adjunto en Anexo (B.3), en MR-HR para todas las
parametrizaciones del nmero de Strickler, considerando los valores de dimetros caractersticos
indicados en la tabla (4.7) y con un rango de caudales [1 - 50] m3/s, con un paso de 1 m3/s. Los
resultados completos se adjuntan en el Anexo (B.5). En este trabajo se presentan los resultados
resumidos mediante figuras.
Figura 4.10: Variacin del Strickler promedio vs. caudal, para diferentes dimetros.
57
De la figura (4.10) se aprecia el correcto funcionamiento de MR-HR ya que para todas las
parametrizaciones, el dimetro caracterstico mximo es que el produce mayor macro rugosidad.
Esto se debe a que el inverso de la rugosidad relativa es menor para dimetros ms grandes, lo
que implica nmeros de Strickler mayores. Adems se aprecia que para caudales menores el
nmero de Strickler promedio es mayor, esto se debe a que a menor caudal disminuyen las alturas
de escurrimiento y, consecuentemente, los radios hidrulicos, lo que implica que disminuya el
inverso de la rugosidad relativa.
Por otra parte, se aprecia para dimetros caractersticos mnimos y medios que pasado el
umbral de 30 m3/s el Strickler promedio vale 0.12. Esto indica que bajo dichas condiciones los
problemas son de lecho rugoso, por lo que no se consideraron en este anlisis.
Por lo tanto, la aplicacin se centra en el rango de caudales [1-50] m3/s y utilizando el

dimetro caracterstico mximo. En particular, se escogi el caudal 5 m3/s y se graficaron los ejes
hidrulicos brindados por las diferentes parametrizaciones, lo cual se muestra en la figura (4.11).
Vale la pena destacar de esta figura, que los ejes hidrulicos calculados mediante las distintas
parametrizaciones aparentan ser iguales entre si, pero no lo son. Esta distorsin ocurre porque
para este caso en particular, los ejes hidrulicos difieren en un orden de magnitud de centmetros
y la escala de la figura (4.11) es demasiado grande para apreciar estas diferencias4. Entonces, no
habra motivo para pensar que pudiesen presentarse diferencias apreciables en esta escala para
caudales mayores, dado que la figura (4.10) muestra que en general el lecho tiende a tornarse
rugoso, es decir, deja de ser macro rugoso para valores de nmeros de Strickler iguales a 0.12, a
medida que se incrementa el caudal. Para solucionar el tema de escalas, se recurre a realizar el
anlisis a secciones transversales particulares, donde se maximicen las diferencias de alturas de
escurrimiento y a una escala apreciable visualmente. Adems es de especial inters seleccionar
un par de secciones transversales que representen tanto rpidos como pozas, lo cual es
determinado por el nmero de Froude. Con esta comparacin se podr apreciar las variaciones
entre las diferentes parametrizaciones.
4
Para consultar en detalle los resultados computados por MR-HR, vase el Anexo (B.5), el cual representa los
resultados completos de la presente aplicacin.
58
Figura 4.11: Resultados MR-HR: Ejes hidrulicos considerando diferentes parametrizaciones,

para caudal 5 m3/s.
Las secciones transversales seleccionadas para el anlisis fueron:
Tabla 4.8: Secciones Transversales Consideradas.
Tipo Seccin N Froude Distancia desde Baker 1

Transversal [km]
Rpido 75831.2 >1 75.8
Poza 127168.5 <1 21.0
59
Figura 4.12: Resultados MR-HR: Cotas superficie libre del agua vs. Caudal (Seccin Subcrtica)
De la figura (4.12) es interesante notar que en un comienzo, para caudales menores a 25

3
m /s, la parametrizacin de Ayala y Oyarce brinda las mayores cotas de superficie libre, luego
para caudales superiores a 25 m3/s, es la parametrizacin de Parker y Peterson la que brinda las
mayores cotas. Tambin, es interesente notar que para caudales sobre los 40 m3/s, Limerinos
brinda cotas de superficie libre similares a las de Keulegan. Segn la figura (3.7), la cual
representa el resumen de parametrizaciones ajustadas en funcin del inverso de rugosidad
relativa, Limerinos brinda valores del nmero de Strickler, por lo general, mayores que la
parametrizacin de Keulegan. En este contexto, debera brindar mayores valores de alturas de
escurrimiento, salvo en las cercanas del lecho rugoso en donde todas las parametrizaciones
convergen a 0.12. Entonces surgen las interrogantes si Keulegan y Limerinos convergen al lecho
rugoso por qu Ayala y Oyarce y Parker y Peterson no?, deberan brindar todas las
parametrizaciones la misma altura de escurrimiento si es que estn convergiendo al lecho
rugoso?, por qu se alternan las mayores cotas de superficie libre en funcin del caudal?. Estas
60
preguntas sugieren analizar como evoluciona espacialmente el nmero de Strickler en funcin de

la longitud del cauce para diferentes caudales, lo cual se muestra en las figuras (4.14) y (4.15).
Figura 4.13: Resultados MR-HR: Cotas superficie libre del agua vs. Caudal (Supercrtico)
Para la condicin de flujo supercrtico, en la figura (4.13) se aprecia que para caudales
pequeos, comprendidos en el rango 1 a 10 m3/s, se observan diferencias considerables entre las
diferentes parametrizaciones. Luego, en el rango de 20 a 30 m3/s, las cotas se tornan casi
idnticas. Finalmente, en el rango de caudales mayores a 35 m3/s, se observa que tienden a
diferir, levemente, las parametrizaciones entre s. Esto nuevamente sugiere el anlisis de la
evolucin espacial del nmero de Strickler considerando distintos caudales, dado que las
parametrizaciones son distintas, pero segn la figura brindan valores similares aun cuando los
lechos son macro rugosos en los rango de caudales bajos de 1 a 30 m3/s.
61
Para una mejor interpretacin de la evolucin espacial del nmero de Strickler en funcin
de diferentes caudales, se seleccion un tramo particular del cauce en estudio. Este tramo
corresponde al comprendido entre los kilmetros [15.5 - 18.5], medidos con respecto al origen en
Baker 1. Se adjunta en el Anexo (B.5) dos videos que dan cuenta de est evolucin completa.
Ntese que la evolucin espacial del nmero de Strickler conduce a la compresin de los
comportamientos de las figuras (4.12) y (4.13) debido a que implcitamente se analiza como
evolucionan los coeficientes de Manning y por tanto las alturas de escurrimiento.
Figura 4.14: Resultados MR-HR: Evolucin Espacial de Nmeros de Strickler (Caudal 10 m3/s)
Fuente: Elaboracin propia
De la figura (4.14) se aprecia que no existe una nica parametrizacin que entregue los
mayores valores del nmero de Strickler, sino que ms bien se observa una alternancia entre
stas. Entre las distancias 15.5 km y 16.5 km se tiene que los valores del nmero de Strickler, en
orden decreciente, son: Ayala y Oyarce, Limerinos, Parker y Peterson y Keulegan. En el
kilmetro 17 se observa el siguiente ordenamiento decreciente: Parker y Peterson, Ayala y
Oyarce, Limerinos y Keulegan. Hacia el final del tramo, entre los kilmetros 18 y 18.5, los
nmeros de Strickler en orden decreciente son Limerinos, Ayala y Oyarce, Parker y Peterson y
62
Keulegan, aun cuando Limerinos, Ayala y Oyarce, Parker y Peterson pueden considerarse
iguales.
De la figura (4.15) se aprecian cambios considerables con respecto a la figura (4.14).

Primero que todo, las amplitudes del nmero de Strickler han disminuido, lo cual es correcto
debido a que a mayores caudales el escurrimiento tiende a ser rugoso, como se apreci en la
figura (4.10). Entre las distancias 15.5 km y 16 km, se tiene un tramo donde la parametrizacin
que brinda el mayor nmero de Strickler es Parker y Peterson. Lo mismo ocurre en el tramo
comprendido entre los kilmetros 18 y 18.5. En este ltimo tramo tambin se observa que la
parametrizacin de Limerinos converge ms rpido que la de Keulegan al escurrimiento rugoso.
Los comportamientos descritos en base a las figuras (4.14) y (4.15), alternancia en los
valores del nmero de Strickler a lo largo del cauce proporcionados por las distintas
parametrizaciones y velocidad de convergencia a escurrimientos sobre lechos rugosos, no se
explican totalmente en trminos de la figura (3.7) que representa el resumen de parametrizaciones
63
ajustadas en funcin del inverso de la rugosidad relativa. Esto se debe a que la evolucin espacial
del nmero de Strickler en funcin del caudal, est ntimamente ligado con la parametrizacin
seleccionada, la forma del cauce particular estudiado, las condiciones de borde y el caudal. Esto
ocurre debido a que el nmero de Strickler es variable a lo largo de todo el cauce y entonces,
dependiendo de la parametrizacin seleccionada, los nmeros de Strickler pueden provocar
secciones de control hidrulico, cambiando los tipos de escurrimientos y por ende las alturas de
escurrimiento.
Para ejemplificar an mejor los conceptos anteriormente descritos, se presenta la figura

(4.16). Esta figura muestra la seccin transversal nmero 117610.3, la cual es rugosa para todo el
rango de caudales, es decir, tiene un Strickler constante e igual a 0.12. Esto implica a que en este
punto del cauce en particular, los coeficientes de Manning brindados por las distintas
parametrizaciones sean idnticos. Sin embargo las cotas de superficie libre del agua no lo son.
Figura 4.16: Resultados MR-HR: Cotas superficie libre del agua vs. Caudal
64
Las cotas de superficie libre del agua no son idnticas a pesar de ser un escurrimiento
rugoso en sta seccin, esto se debe a la influencia que ejerce en este punto la distribucin
espacial de los nmeros de Strickler y, consecuentemente, los coeficientes de Manning a lo largo
del cauce, segn la parametrizacin seleccionada.
En la segunda etapa se plante el anlisis de las diferencias brindadas por las diferentes
parametrizaciones en su forma ms general, es decir, permitiendo que la macro rugosidad
dependiera tanto del desarrollo del eje hidrulico como de la variacin longitudinal de los
dimetros caractersticos a lo largo del cauce versus proponer un coeficiente de Manning
uniforme para todo el cauce, basado en el clculo previo del dimetro caracterstico asociado al
valor medio ponderado.
La interpolacin lineal de los dimetros caractersticos se basaron en los datos

proporcionados por la tabla (4.3). Esta interpolacin se resume en la siguiente figura5:
Figura 4.17: Interpolacin lineal realizada a los sedimentos del sustrato del lecho
5
Para consultar en detalle las interpolaciones realizadas por MR-HR, vase el Anexo (B.5).
65
Luego se definieron coeficientes de Manning uniformes basados en los resultados

obtenidos en la tabla (4.7). Para esto se consideraron los dimetros caractersticos asociados al
valor medio ponderado, con la finalidad de que stos fueran representativos del cauce en estudio.
Adems se asumi que el coeficiente de Manning dependa slo del dimetro caracterstico de los
sedimentos y no del desarrollo del eje hidrulico, es decir, se asumi que el lecho era rugoso con
un nmero de Strickler uniforme a lo largo del cauce igual a 0.12. Entonces, con la ayuda de la
ecuacin (2.4) se calcularon dos coeficientes de Manning para todo el cauce, los cuales se
muestran en la siguiente tabla:
Tabla 4.9: Coeficientes de Manning Uniformes Basados en Valores Medios Ponderados de

Dimetros Caractersticos.
Coeficientes de Manning
(sm-1/3)
n (D90) n (D84)
0.0212 0.0204
Finalmente, para cuantificar el efecto de considerar un lecho rugoso con un nico dimetro
caracterstico versus uno macro rugoso, se trabaj con las diferencias en las alturas de
escurrimiento. Estas diferencias son con respecto al caso de Manning uniforme correspondientes
a D90 de la tabla anterior. Se escogi ese coeficiente para minimizar las diferencias esperadas e
intentar ser objetivo en la comparacin. En la figura (4.18) se muestran las diferencias
anteriormente sealadas para un caudal de 5 m3/s. Bajo stas condiciones el problema presenta,
en su mayora, rgimen de escurrimiento subcrtico, es decir, nmeros de Froude menores a la
unidad, salvo en las secciones transversales donde ocurren rpidos que corresponden a las
diferencias nulas entre las diferentes parametrizaciones y la proposicin del coeficiente de
Manning uniforme a lo largo de todo el cauce. Estas secciones se caracterizan por tener gran
pendiente y un nmero de Froude muy cercano a la unidad, produciendo crisis y desconectando
hidrulicamente los escurrimientos entre los diferentes tramos del ro. Para fijar idea, la figura
(4.11) muestra las zonas donde se producen las crisis, las cuales corresponden a las mismas
secciones mostradas en la figura (4.18). En los extremos es posible apreciar que, para este caso en
particular, las condiciones de borde no ejercen demasiada influencia en el cauce dado que
cercano a ellas se producen crisis.
66
Figura 4.18: Diferencia entre Ejes Hidrulicos con Respecto a Manning Uniforme (Ds90), Q=5
m3/s
Con respecto a las diferencias, como se dijo anteriormente, no existe un predominio

exclusivo de alguna parametrizacin por sobre las dems, sino que se alternan. En los primeros
70 kilmetros se observa un aumento paulatino de las diferencias debido a que en este tramo la
interpolacin lineal de los sedimentos del lecho es creciente a medida que se acerca al origen del
cauce, lo que implica que se agudiza la macro rugosidad. Por el contrario, posterior a los 70
kilmetros los dimetros caractersticos son decrecientes y junto con ello tambin la macro
rugosidad. Las diferencias se tornan negativas dado que, en sus tramos correspondientes, el lecho
se comporta tanto macro rugoso como rugoso. Si bien al comportarse como macro rugoso el
nmero de Strickler es superior a 0.12, el dimetro caracterstico es inferior al medio ponderado
por lo que se produce una competencia entre el nmero de Strickler y la raz sexta de los tamaos
de sedimentos en la determinacin de los coeficientes de Manning, lo cual se puede apreciar en la
ecuacin (2.4). Esta competencia conduce a que los coeficientes de Manning sean inferiores a los
67
proporcionados por el valor medio ponderado y consecuentemente las alturas de escurrimiento.

En la zonas donde el lecho se comporta como rugoso es evidente que los coeficientes de Manning
sean inferiores a los proporcionados por el valor medio ponderado Ds90, ya que el tamao de los
sedimentos es inferior. La siguiente figura da cuenta de cuando el lecho se comporta como macro
rugoso y cuando no.
De la figura (4.19) se observa que aguas abajo, en la condicin de borde km. 148, se genera
macro rugosidad. La condicin de borde es de altura normal con la respectiva pendiente del
terreno de esa zona. Inmediatamente aguas arriba de la condicin de borde no se observa macro
rugosidad dado el cambio de geometra desde una muy llana a una desnivelada. Para el caudal 5
m3/s el escurrimiento en esa zona es subcrtico por lo que la altura de escurrimiento es
68
determinada por la macro rugosidad, lo que se puede interpretar en la figura (4.18) como las
diferencias de alturas de escurrimiento brindadas por las diferentes parametrizaciones.
Los resultados anteriormente brindados no deben interpretarse como absolutos, debido a

que son referentes al cauce particular estudiado y a las condiciones a las que est sometido. Por
ejemplo, en un cauce con escurrimiento subcrtico en donde no se generen crisis, el efecto de la
macro rugosidad podra remontarse aguas arriba del cauce, llegando a su inicio. Esto podra
agudizar aun ms las diferencias en las alturas de escurrimiento en funcin de cuan macro rugoso
se comporte el lecho. Por el contrario, las diferencias podran ser menores si el lecho se comporta
en la transicin entre la rugosidad y la macro rugosidad. Sin embargo, las diferencias entre alturas
de escurrimientos aplicado en este anlisis son slo para fijar idea, debido a que no son el tema de
fondo sino la variacin de los coeficiente de Manning, ya que por ejemplo para una geometra de
un cauce llano, en trminos de sus secciones transversales, el efecto de la macro rugosidad no se
vera reflejado en forma tan marcada como en el caso anterior. Por otra parte, la seleccin del
coeficiente de Manning uniforme para todo el cauce, se realiz en pro de minimizar las
diferencias en base a la experiencia del autor de este trabajo. Sin embargo, no es justo comparar
un slo coeficiente de Manning uniforme a lo largo de todo el cauce, basado en un dimetro
caracterstico "representativo", con un enfoque que considera toda la variabilidad de los
sedimentos del lecho en el eje longitudinal. Ms bien la comparacin debera incluir coeficientes
de Manning uniformes, pero variables por tramos, lo cual tendera a minimizar las diferencias
anteriormente mostradas. Por lo tanto, el ejemplo aqu propuesto tiene como objetivo fijar idea de
lo que ocurre en trminos generales y por ningn motivo deben extrapolarse sus resultados a otro
tipo de cauces sometidos a diferentes condiciones.
69
Patricio Pastenes Araos Captulo 5: Discusin y Alcances
Captulo 5. Discusin y Alcances
En lo que respecta al diseo e implementacin de MR-HR, es importante reconocer que las

principales limitaciones, vistas desde el punto de vista hidrulico, estn condicionadas por la
interaccin con HEC-RAS. Sin embargo y segn los antecedentes revisados en este trabajo, se
concluye que MR-HR es capaz de enfrentar problemas frecuentes de la ingeniera hidrulica.
Esto se refleja en que se pueden abordar: cauces no prismticos, problemas en rgimen
turbulento, problemas con pared hidrodinmicamente rugosa, problemas con pendientes
inferiores al 10%, etc. Estas son condiciones que se dan en general en los cauces naturales y son
muy recurrentes en el que hacer cotidiano del ingeniero hidrulico.
Con respecto a la aplicacin de MR-HR al ro Baker, vale la pena mencionar que se

utilizaron caudales tan bajos, de 1 a 50 m3/s, dado que en este rango el lecho se comportaba como
macro rugoso, y para valores superiores del caudal la aplicacin de MR-HR perda un poco de
sentido al comportarse el lecho como rugoso, debido a que la altura de escurrimiento en el cauce
provocaba inversos de rugosidades relativas superiores a 10.
En cuanto a los resultados brindados por MR-HR, considerando slo un dimetro

caracterstico para todo el cauce, fue posible concluir que las alturas de escurrimiento estn
fuertemente influenciadas por la distribucin espacial de los nmeros de Strickler que a su vez
depende de la parametrizacin seleccionada, el caudal, las condiciones de borde del problema y
la geometra. Esto produce una carencia de generalidad en trminos de los resultados esperados
en funcin de las diversas parametrizaciones, es decir, no es posible asumir que siempre la
seleccin de una determinada parametrizacin brindar mayores o menores alturas de
escurrimiento en una determinada posicin del cauce en comparacin con las dems. La
influencia que ejerce la distribucin espacial de los nmeros de Strickler en un punto particular
del cauce, podra producir secciones de control hidrulico que alteren considerablemente el
desarrollo del eje hidrulico entre una parametrizacin y otra. Sin embargo, en general es posible
esperar que el nmero de Strickler promediado a lo largo de todo el cauce se comporte en forma
acorde a las parametrizaciones revisadas en este trabajo, considerando su dependencia nica y
exclusivamente con el inverso de la rugosidad, es decir, se esperan nmeros de Strickler
promedio mayores segn la seleccin de la parametrizacin, pero esto no implica que se pueda
70
Patricio Pastenes Araos Captulo 5: Discusin y Alcances
extender al desarrollo del eje hidrulico. Por lo tanto, y dada la falta de generalizacin, en caso de
no contar con una seleccin adecuada de alguna parametrizacin, se recomienda evaluar los
resultados de todas las parametrizaciones en la aplicacin de MR-HR y adoptar como vlido el
eje hidrulico en funcin del tipo de problema que se enfrenta y su criterio.
En el contexto de comparar los resultados proporcionados por MR-HR, considerando toda

la informacin granulomtrica del lecho, versus la seleccin de un coeficiente de Manning
uniforme a lo largo de todo el cauce, se apreci que en particular para el cauce estudiado las
condiciones de borde del problema no juegan un rol predominante en el desarrollo de todo el eje
hidrulico, dadas las desconexiones hidrulicas que ocurren a lo largo del cauce. Sin embargo,
esto ocurre particularmente en el cauce aqu estudiado y no debe ser siempre interpretado de la
misma manera. En lo que respecta a la influencia que ejerce la distribucin espacial de los
nmeros de Strickler en una determinada posicin, se observ nuevamente una carencia de
generalidad en cuanto a las alturas de escurrimiento brindadas por las diferentes
parametrizaciones en trminos que no existe una parametrizacin que siempre brinde mayores
alturas de escurrimiento por sobre las dems. Con respecto a las diferencias de alturas de
escurrimiento obtenidas, stas son slo para fijar idea, ya que estn sujetas a las condiciones
particulares del cauce, y no representan el problema de fondo el cual es reconocer que la ecuacin
de Manning no es vlida para lechos macro rugosos, por lo que MR-HR si resuelve el problema
de fondo que es la determinacin correcta de los coeficientes de Manning con la ayuda de HEC-
RAS, y consecuentemente los respectivos ejes hidrulicos.
Por lo tanto, se recomienda no generalizar los resultados obtenidos en la aplicacin

particular de MR-HR realizada en este trabajo y en caso que se disponga de alturas de
escurrimiento, compararlas con los resultados brindados tanto por MR-HR como los resultados
proporcionados a partir de tabulaciones de coeficientes de Manning para decidir cual representa
mejor el cauce en estudio.
71
Patricio Pastenes Araos Bibliografa
Bibliografa
1. AYALA, L. y OYARCE, O. (1993). Prdida de carga de cauces cordilleranos sujetos a

acorazamiento y arrastre En: XI CONGRESO Chileno de Ingeniera Hidrulica. Concepcin,
Universidad de Concepcin. Vol. 2, pp. 1-15.
2. BRADLEY, J. C. AND MILLSPAUGHASDF A. (2009). Advanced Programming Using

Visual Basic 2008. 4th Edition USA, McGraw-Hill. 672p.
3. BRADLEY, J. C. AND MILLSPAUGHASDF A. (2010). Programming in Visual Basic

2010. USA, McGraw-Hill. 704p.
4. CENTRALES HIDROELECTRICAS DE AYSEN S.A. (2008). Estudio de Impacto

Ambiental: Proyecto Hidroelctrico Aysn. 1973p
5. CHOW V. T. (1994). Hidrulica de Canales Abiertos. Santaf de Bogot, McGraw-Hill.

667p.
6. KEULEGAN, G. (1938). Laws of turbulent flow in open channels. National Bureau of

Standards, Research Paper 1151 (U.S.). 21(6):707-741.
7. LIMERINOS, J. (1970). Determination of the Manning Coefficient from Measured Bed

Roughness in Natural Channels. Geological Survey Water-Supply Paper 1898-B.
8. MICROSOFT (2011a). Developing Applications with Visual Basic. [en lnea]

<http://msdn.microsoft.com/en-us/library/k6a6etxs>. [consulta : 25 septiembre 2011].
9. MICROSOFT (2011b). Visual Basic Programming Guide. [en lnea] <

http://msdn.microsoft.com/en-us/library/y4wf33f0>. [consulta : 30 septiembre 2011].
72
10. NIO, Y. (2003). Resistencia Hidrulica en Canales Aluviales. En: Apuntes del Curso CI61F
Transporte Hidrulico de Slidos, Departamento de Ingeniera Civil, Universidad de Chile. pp.
24-28.
11. OYARCE, O. (1993). Diseo, Construccin y Operacin del Modelo Fsico de un Prototipo
Representativo de Cauces Cordilleranos Andinos. Memoria Para Optar al Ttulo de Ingeniero
Civil, Universidad de Chile, Facultad de Ciencias Fsicas y Matemticas. 159p.
12. PARKER, G. AND PETERSON, A. (1980). Bar Resistance Of Gravel-Bed Streams. Journal
Of The Hydraulics Division, ASCE. 106(HY10):1559-1575.
13. STRICKLER, A. (1923). Beitrge zur Frage der Geschwindigkeitsformel und der
Rauhigkeitszahlen fr Strme, Kanle und geschlossene Leitungen (Some contributions to the
problem of velocity formula and roughness factors for rivers, canals, and closed conduits).
Mitteilungen des eidgenssischen Amtes fr Wasserwirtschaft, N16, Bern Switzerland.
14. UNIVERSIDAD DE CHILE (2010). Antecedentes Complementarios a la Caracterizacin

Granulomtrica de los Ros Baker y Pascua. Grupo de Investigacin en Procesos de Transporte y
Flujos Ambientales, Departamento de Ingeniera Civil, Universidad de Chile.
15. U.S. ARMY CORPS OF ENGINEERS (USACE) (1994). Methods for Predicting n Values
for the Manning Equation. En: Engineering and Design: Hydraulic Design of Flood Control
Channels, Engineer Manual No. 1110-2-1601. Washington, USA. pp. 5.1-5.14.
16. U.S. ARMY CORPS OF ENGINEERS (USACE) (2010a). HEC-RAS River Analysis System:
Applications Guide, version 4.1. California, USA, Institute for Water Resources, Hydrologic
Engineering Center (HEC). 351p.
17. U.S. ARMY CORPS OF ENGINEERS (USACE) (2010b). HEC-RAS River Analysis System:
Hydraulic Reference Manual, version 4.1. California, USA, Institute for Water Resources,
Hydrologic Engineering Center (HEC). 417p
73
18. U.S. ARMY CORPS OF ENGINEERS (USACE) (2010c). HEC-RAS River Analysis System:
User's Manual, version 4.1. California, USA, Institute for Water Resources, Hydrologic
Engineering Center (HEC). 790p
74
Patricio Pastenes Araos ANEXO A: Manual de Uso MR-HR
Anexo A. Manual de uso MR-HR
A.1 Generalidades
MR-HR es una herramienta de acople a HEC-RAS, que calcula los coeficientes de

Manning mediante diferentes parametrizaciones del nmero de Strickler. Para esto realiza las
iteraciones que el usuario desee hasta lograr la convergencia de los coeficientes. Si bien, MR-HR
est principalmente enfocado en problemas de lechos macro rugosos, tambin resuelve problemas
de lechos rugosos mediante transiciones suaves realizadas a las parametrizaciones de Keulegan
(1938), Limerinos (1970), Parker y Peterson (1980) y Ayala y Oyarce (1993).
A.2 Instalacin
MR-HR consiste en tres archivos que deben estar localizados en una misma ruta sin
importar esta. En la siguiente figura se muestra un ejemplo:
Figura A-1: Ejemplo de Instalacin MR-HR.

I
Los tres archivos que se aprecian en la figura (A-1) son: Iterop.RAS41.dll, MR-HR.exe y
Manual de Uso.pdf. El primero representa el enlace entre MR-HR y HEC-RAS 4.1.0. El segundo,
que tiene por icono una gota de agua, representa el archivo ejecutable que abre el programa y el
tercero este Manual de Uso.
A.3 Requisitos
A.3.1 Computacionales
Los requerimientos computacionales que necesita MR-HR para su correcto funcionamiento

son los siguientes:
Tener instalado el sistema operativo Microsoft Windows6 XP, Microsoft Windows

Vista o Microsoft Windows 7, ya sea en sus versiones de 32-bits como la de 64-bits.
Tener instalado Microsoft .NET Framework 7, versin 2.0 o superior (requisito que
se cumple por defecto al tener instalado Microsoft Windows XP o superior).
Tener instalado HEC-RAS 4.1.0,8 y crear los proyectos en la misma versin.
Tener instalado algn lector de archivos PDF para la visualizacin del Manual de
Uso, como por ejemplo Adobe Reader9.
6
La configuracin regional de Windows debe tener como separador de decimales el punto.
7
Programa de distribucin gratuita. Ms informacin: http://www.microsoft.com/
8
Programa de distribucin gratuita. Ms informacin: http://www.hec.usace.army.mil/software/hec-ras/
9
Programa de distribucin gratuita. Ms informacin: http://get.adobe.com/es/reader/
II
A.3.2 Proyecto Importado desde HEC-RAS
Para el correcto funcionamiento de MR-HR, los proyectos creados en HEC-RAS 4.1.0

deben contener slo un tramo y un slo ro. Este requisito no implica prdida de generalidad
desde el punto de vista del cauce principal, dado que HEC-RAS aplica en nodos (puntos donde
confluyen o divergen dos o ms ros) igualdad de energa y en caso que se requiera la funcin de
momentum. Por lo tanto, teniendo en cuenta estos principios es posible modelar slo el cauce
principal e ingresar las magnitudes de caudales que ingresan o egresan del cauce principal sin
afectar el desarrollo del eje hidrulico.
Figura A-2: Proyecto creado en HEC-RAS no compatible con MR-HR

La figura (A-2), no es compatible con MR-HR dado que contiene tanto un afluente como
un efluente en la seccin del tramo principal 3030.*. Para resolver esto se debe modelar slo un
cauce principal y hacer un balance de caudal en la seccin donde converge el afluente y donde
diverge el efluente. Esto se denomina nodo fantasma. Entonces en este nodo fantasma (seccin
III
3030.*) se debe ingresar el balance total de caudal, es decir, el caudal proveniente de aguas arriba
ms el caudal aportado por el afluente menos el caudal del efluente.
A continuacin se muestra un proyecto correctamente creado en HEC-RAS, para su

posterior utilizacin en MR-HR:
Figura A-3: Proyecto creado en HEC-RAS compatible con MR-HR

Es posible observar de la figura (A-3) que se eliminaron los afluentes y efluentes. Para fijar
idea, supongamos que desde aguas arriba (seccin 5010) provienen 5 m3/s, el afluente aporta 2
m3/s y el efluente extrae 3 m3/s. Entonces el balance de caudal en el nodo fantasma es de 4 m3/s.
Esto se refleja en la siguiente figura:
IV
Figura A-4: Proyecto creado en HEC-RAS compatible con MR-HR.

La figuras (A-3) y (A-4) producen un eje hidrulico idntico al de la figura (A-2), desde el
punto de vista del cauce principal, mediante el ingreso de balances de caudal en los nodos
fantasmas. Recordar que para esta aplicacin se debe tener en mente los principios de
conservacin de la Energa y la funcin de Momentum, que son los principios que HEC-RAS
utiliza por defecto.
Otro requisito del proyecto creado en HEC-RAS, es que solo puede existir un perfil de
flujo, si se quieren evaluar ms de un perfil de flujo, MR-HR permite la evaluacin de mltiples
perfiles de flujo va configuracin adentro de MR-HR. A continuacin se muestra una
configuracin de perfiles de flujo no compatible con MR-HR, representada por la figura (A-5), y
una configuracin de perfil de flujo que es compatible con MR-HR, representada por la figura
(A-6).
V
Figura A-5: Proyecto creado en HEC-RAS no compatible con MR-HR.

Figura A-6: Proyecto creado en HEC-RAS compatible con MR-HR.

En el punto (A.4.3) que representa el Submen de Configuracin de MR-HR se explicar

como configurar proyectos, adentro de MR-HR, para uno o ms perfiles de flujo.
VI
Finalmente, para el correcto funcionamiento de MR-HR los proyectos creados en HEC-

RAS deben tener asociados un PLAN, GEOMETRA y FLUJO.
Figura A-7: Proyecto Creado en HEC-RAS compatible con MR-HR.

De la figura (A-7) es posible observar que el proyecto creado en HEC-RAS, contiene asociado:
PLAN : Plan Permanente

GEOMETRA : Geometra
FLUJO : flujo
Una vez creada tanto la geometra, el flujo, como el plan, estos se deben guardar en el
proyecto HEC-RAS y recin entonces, estarn asociados al proyecto.
A.4 Componentes
El diseo de MR-HR fue pensado en ser lo ms intuitivo posible, es por esto que sus
componentes se van habilitando uno tras otro, es decir, una vez que se cumpla un requisito de la
configuracin se habilitara el prximo, hasta la eventual ejecucin del programa y la habilitacin
de los resultados brindados por ste.
VII
A.4.1 Ventana Principal
A continuacin se muestra la ventana principal de MR-HR y se enumeran los principales

componentes que este programa contiene.
Figura A-8: Ventana Principal de MR-HR.

El significado de la numeracin descrita en la figura (A-8) se explica a continuacin:
1. Barra de Men Contextual: mediante esta barra se importa el proyecto creado en HEC-
RAS, se configura la parametrizacin seleccionada, el o los dimetros caractersticos, el
nmero de iteraciones que se deseen realizar, se define el o los perfiles de flujo para los
cuales se va a iterar, se ejecuta el proyecto una vez configurado, se muestran los
resultados y se recurre a ayuda en caso que se requiera. Cada elemento (Archivo /
Configuracin / Ejecutar / Resultados / Ayuda) despliega un submen donde se accede a
configurar y/o mostrar lo anteriormente sealado.
VIII
2. Proyecto Cargado: Este cuadro de texto brinda la informacin del proyecto cargado
(importado) de HEC-RAS. Proporciona el nombre del proyecto, plan, geometra y el
archivo de flujo asociado. Una vez cargado (importado) el proyecto este cuadro de texto
se actualiza automticamente.
3. Parametrizacin: Este cuadro de texto brinda la informacin de la configuracin

referente a la parametrizacin del nmero de Strickler seleccionada, el o los dimetros
caractersticos ingresados, el nmero de iteraciones seleccionado y la calibracin del
coeficiente de Manning. La calibracin del coeficiente de Manning se refiere a un
desplazamiento ya sea positivo o negativo de todos los coeficientes de Manning en una
determinada magnitud ingresada por el usuario. El objetivo de este parmetro es calibrar
los coeficientes de Manning en caso de mucha irregularidad en el cauce, vegetacin,
obstrucciones, etc., adems si se dispone de alturas de escurrimiento medidas en terreno y
las parametrizaciones seleccionadas no brindan resultados acordes a los de terreno,
entonces la calibracin del coeficiente de Manning sumar o restar una magnitud
predefinida por el usuario a todos los coeficientes de Manning para intentar ajustar el eje
hidrulico a los datos obtenidos en terreno. Mediante se configure el proyecto, el cuadro
de texto se ir actualizando constantemente, desplegando la informacin ingresada por el
usuario.
4. Configuracin Caudal Variable: Este cuadro de texto brinda la informacin de la

configuracin referente a el o los perfiles de flujo para los cuales se va a iterar en funcin
del coeficiente de Manning. La seccin de cambio de caudal, siempre es la aguas arriba,
donde comienza el cauce. Una vez configurado el caudal, este cuadro de texto se
actualizara automticamente.
A.4.2 Submen Archivo
El submen archivo es el primer submen, de izquierda a derecha, que aparece en la barra

contextual (ver figura A-8: Ventana Principal de MR-HR). Al hacer click en el submen archivo
se despliegan las opciones:
IX
1. Cargar Proyecto : Es donde se importa el proyecto creado en HEC-RAS.

2. Salir : Se abandona MR-HR.
Ntese que el submen Archivo y Ayuda son los nicos habilitados, dado que no tiene
sentido configurar un proyecto si antes no ha sido cargado o importado. Lo mismo ocurre con los
submens Ejecutar y Resultados que tambin carecen de sentido si el proyecto no ha sido cargado
o importado.
Figura A-9: Submen Archivo

Al hacer click en Cargar Proyecto, MR-HR lanzar el siguiente mensaje: "Se recomienda
guardar y cerrar todas las instancias de HEC-RAS actualmente abiertas, ya que se cerraran
automticamente". Esto se debe a que a partir de este punto MR-HR controlar por completo a
HEC-RAS y es posible que si no ha guardado algn trabajo este se pierda.
Luego debe navegar hasta seleccionar el proyecto creado en HEC-RAS. Ntese que solo es
posible seleccionar archivos *.prj que es la extensin de los proyectos creados en HEC-RAS. La
figura (A-10) muestra dicha situacin.
X
Figura A-10: Importacin proyecto creado en HEC-RAS en MR-HR.

Una vez importado su proyecto, MR-HR desplegar el siguiente mensaje: "MR-HR

realizar cambios en su proyecto creado en HEC-RAS. Se recomienda crear copias de respaldo
de sus archivos". Esto se debe a que MR-HR cambiar los perfiles de flujo y coeficientes de
Manning a medida que vayan ocurriendo las iteraciones, en el proyecto original importado. La
importancia que tiene esto es que una vez finalizados los clculos, el usuario puede recurrir al
proyecto modificado por MR-HR, el cual contendr la informacin tanto del flujo como de los
coeficientes de Manning de la ltima iteracin. Entonces podr ejecutarlo normalmente en HEC-
RAS y podr consultar las variables que le parezcan de inters, no consideradas en MR-HR, tales
como: el nmero de Froude, pendientes energticas, etc.
Finalmente, luego que el proyecto haya sido cargado o importado se actualizar la ventana
principal de MR-HR y se habilitar el submen Configuracin de la barra de men Contextual.
Esto se puede apreciar en la figura (A-11).
XI
Figura A-11: Actualizacin de MR-HR una vez cargado (importado) el proyecto de HEC-RAS
Ntese de la figura (A-11) que se actualiz el campo Proyecto Cargado, que contiene la
informacin del proyecto importado, adems se habilit el submen Configuracin. Sin embargo,
los submens Ejecutar y Resultados an aparecen inactivos, dado que no se ha configurado el
proyecto en MR-HR, por lo que carece de sentido la Ejecucin y tambin la visualizacin de
Resultados.
A.4.3 Submen Configuracin
Una vez cargado o importado el proyecto, se habilita el submen Configuracin. Este

contiene cuatro opciones las cuales son:
4.3.1 Parametrizacin
4.3.2 Dimetros Caractersticos
4.3.3 Criterio Convergencia
4.3.4 Caudal Variable
XII
Figura A-12: Submen Configuracin

Es posible apreciar de la figura (A-12) que solo la opcin Parametrizacin est disponible.
Luego de completada esta opcin se irn habilitando sucesivamente: Dimetros Caractersticos,
Criterio de Convergencia y Caudal Variable.
4.3.1 Parametrizacin
Figura A-13: Submen Configuracin / Parametrizacin

XIII
En la figura (A-13) se aprecia que es posible seleccionar la parametrizacin del nmero de

Strickler deseada, es decir: Keulegan, Limerinos, Parker y Peterson, Ayala y Oyarce. Existe la
opcin de ingresar una magnitud para calibrar los coeficientes de Manning, denominada
Calibracin Manning. Esta magnitud puede ser tanto positiva como negativa y tiene como
objetivo sumar o restar a todos los coeficientes de Manning una cantidad arbitraria definida por el
usuario para ajustar el eje hidrulico. Esto tiene sentido si se tienen mediciones de terreno de
alturas de escurrimiento y si ninguna parametrizacin se ajusta bien a las mediciones de terreno
y/o el cauce presenta muchas singularidades, es decir, mucha vegetacin, obstrucciones, curvas,
etc. El valor por defecto es cero y se debe ser cuidadoso en la magnitud de Calibracin ingresada,
pues un valor negativo muy pequeo podra tornar los coeficientes de Manning negativos, lo cual
carece de sentido fsico.
En particular en la figura (A-13) se seleccion la parametrizacin de Parker y Peterson con

Calibracin nula. Una vez configurada la Parametrizacin se habilita la opcin Dimetros
Caractersticos.
4.3.2 Dimetros Caractersticos
Figura A-14: Configuracin / Dimetros Caractersticos

XIV
En la figura (A-14) se aprecia que se reconoce previamente la parametrizacin

seleccionada. Se resumen todas las secciones transversales asociadas al proyecto junto con su
distancia aguas abajo a la prxima seccin transversal medida con respecto al centro del cauce.
Adems, se tiene la opcin de ingresar el dimetro caracterstico de cada seccin transversal. Si
no se conoce todo el detalle en cada seccin transversal, se debe interpolar linealmente. En
particular, se ingres slo un dimetro caracterstico en la seccin transversal 2010.* igual a 0.09
m, al interpolar ste dimetro se convierte en uniforme para todo el resto de las secciones
transversales. Una vez configurados los Dimetros Caractersticos se habilita la opcin Criterio
de Convergencia.
4.3.3 Criterio de Convergencia
Figura A-15: Submen Configuracin / Criterio de Convergencia

En la figura (A-15) se define el nmero de iteraciones que se realizarn para lograr la

convergencia de los coeficientes de Manning, el valor por defecto es 5. Sin embargo, es posible
aumentar este nmero cuanto se desee, pero trae consigo consecuencias en cuanto al tiempo de
clculo. Una vez configurado el Criterio de Convergencia se habilita la opcin Caudal Variable.
4.3.3 Caudal Variable
En la figura (A-16) se aprecia que se reconoce automticamente y no se pueden cambiar el

ro, el tramo y la seccin de aguas arriba donde nace el cauce. En esta seccin (aguas arriba) se
XV
realizarn las variaciones de caudal. Abajo se definen el o los perfiles de flujo para los cuales se
iterar en funcin del coeficiente de Manning.
Para definir un nico perfil de flujo el Caudal Mnimo debe ser igual al Caudal Mximo y
el Incremento de Caudal debe ser cero. Para definir ms de un perfil de flujo, el Caudal Mnimo
debe ser diferente del Caudal Mximo y el incremento de Caudal debe ser distinto de cero.
Figura A-16: Submen Configuracin / Caudal Variable

En particular, en la figura (A-16) se observa que se defini: Caudal Mnimo 5 m3/s, Caudal
Mximo 15 m3/s e incremento de Caudal 0.5 m3/s. Esto quiere decir que se realizaran 21 perfiles
de flujo comenzando desde el caudal 5 m3/s, en la seccin de aguas arriba donde nace el cauce,
incrementndose en 0.5 m3/s hasta llegar a los 15 m3/s.
Una vez completada toda la configuracin, la ventana principal de MR-HR se actualizara

con los datos ingresados y se habilitara el submen Ejecutar. Esto se aprecia en la figura (A-17).
XVI
Figura A-17: Configuracin Finalizada

A.4.4 Submen Ejecutar
Una vez configurado por completo el proyecto en MR-HR se habilita el submen Ejecutar.
Al hacer click en este submen MR-HR preguntar donde desea guardar los resultados. Luego
comienza el proceso de clculo.
En la figura (A-18) es posible apreciar la ventana Clculos, la cual indica el progreso en

porcentaje de los clculos totales que se deben realizar. Una vez finalizada la ejecucin, es decir,
progreso igual al 100%, entonces se habilita el submen Resultados. Ntese que anteriormente ya
se decidi en que lugar quedaran guardados los resultados. Estos resultados estn en formato de
texto plano (*.txt) y contienen toda la memoria de clculo realizada por MR-HR en conjunto con
HEC-RAS.
XVII
Figura A-18: MR-HR Ejecutndose

A.4.5 Submen Resultados
Una vez finalizados los clculos se habilita el submen Resultados. En este se despliega
otro submen que contiene tanto resultados numricos como grficos. En la siguiente figura se
aprecia esto.
XVIII
Figura A-19: Submen Resultados

En la figura (A-19) se puede apreciar dos opciones pertenecientes al submen Resultados,

estas son:
1. Ver Resultados : Esta opcin muestra los resultados numricos, tanto de la ltima
iteracin como de toda la memoria de clculo realizada.
2. Grficos : Esta opcin despliega grficos de alturas de escurrimiento,

velocidad promedio en el canal y coeficientes de Manning vs. Caudales. Cada grfico es
referente a una seccin transversal particular seleccionada por el usuario.
1. Ver Resultados
La opcin Ver Resultados, contiene los resultados numricos realizados por MR-HR. Se
puede escoger entre desplegar los resultados de la ltima iteracin para cada caudal o desplegar
XIX
los resultados completos donde se incluye toda la memoria de clculo. Esto se muestra en la
figura (A-20).
Figura A-20: Submen Resultados / Ver Resultados

Independientemente de la opcin seleccionada, es decir, Resultados Completos o Ultima

Iteracin, las variables que se despliegan son las mismas. Estas son:
Caudal [m3/s]
N Iteracin
Seccin Transversal
Manning [s*m-1/3]
Altura de Escurrimiento [m]
Velocidad Promedio [m/s]
Radio Hidrulico [m]
Diferencia Manning con respecto a la iteracin anterior
N de Strickler
Inverso Rugosidad Relativa (Rh/ds)
Dimetro Caracterstico
En las figuras (A-21) y (A-22) se muestran las opciones de Resultados Completos y Ultima
Iteracin, respectivamente:
XX
Figura A-21: Submen Resultados / Ver Resultados / Resultados Completos

Figura A-22: Submen Resultados / Ver Resultados / Ultima Iteracin

XXI
De la figura (A-21) se puede apreciar que muestra todas las iteraciones para cada uno de los
caudales definidos en los perfiles. Esto se puede apreciar para el caudal 5 m3/s y N de iteracin,
dado que parte desde la primera iteracin hasta la quinta (recordar que como criterio de
convergencia se seleccionaron 5 iteraciones por caudal). La ventaja de esto es que se puede
apreciar como fueron evolucionando o convergiendo los coeficientes de Manning.
De la figura (A-22) solo se aprecia la ltima iteracin para cada caudal, es decir, la
iteracin nmero 5. Si bien no se aprecia la evolucin de los coeficientes de Manning, la columna
que representa la diferencia con respecto al Manning de la iteracin anterior, brinda la
informacin necesaria para saber si convergieron o no los Manning. En caso que la diferencia con
respecto a los Manning anteriores sea demasiado grande, se sugiere aumentar el nmero de
iteraciones. Adems existen opciones de copiar la informacin total o parcialmente al
portapapeles o exportarla en formato de texto plano (*.txt).
Es posible que en algunos casos la informacin no pueda ser desplegada en pantalla dada la
gran extensin que esta podra representar. Esto podra provocar que la memoria disponible del
computador no de abasto para desplegar toda la informacin. Es por esto que, anteriormente,
cuando se realizan los clculos, MR-HR pregunta donde quiere guardar los resultados en formato
de texto plano (*.txt). Si bien no es posible desplegar toda la informacin en pantalla el archivo
de texto puede ser manipulado con softwares como Matlab y Excel para generar grficas o lo que
quiera el usuario.
XXII
2. Grficos
Figura A-23: Submen Resultados / Grficos

Figura A-24: Grfica H vs. Q seccin 3010.*

XXIII
El grfico de la figura (A-24) proviene del submen Resultados / Grficos / Altura vs.
Caudal, figura (A-23). Para generar el grfico se debe seleccionar la seccin transversal que se
desea graficar, en este caso 3030.*, luego se hace click en Graficar y se despliega la grfica
mostrada en la figura (A-24). El procedimiento es anlogo para las grficas de Velocidad vs.
Caudal y Manning vs. Caudal.
Adems de desplegar los grficos en pantalla, MR-HR cuanta con otras opciones, las cuales
son:
1. Copiar : Copia el grfico al portapapeles para luego pegarlo en otro programa.

2. Exportar : Exporta la imagen en formato de mapa de bits (*.bmp).
Es interesante notar, que si en la configuracin de Caudal Variable perteneciente al

submen Configurar, se seleccion solo 1 perfil de flujo, es decir, solo un caudal, estos grficos
estarn deshabilitados dado que son en funcin del caudal.
A.4.6 Submen Ayuda
El submen Ayuda consta con las opciones de: Manual de Uso y Acerca de. Esto se
muestra en la siguiente figura.
XXIV
Figura A-25: Submen Ayuda

Figura A-26: Submen Ayuda / Acerca de..

XXV
Como se puede apreciar de la figura (A-26), la opcin Acerca de.. muestra informacin
referente al desarrollo de MR-HR, los derechos de autor y la informacin de contacto en caso que
se requiera.
La opcin Manual de Uso, hace referencia a este Manual ejecutndolo en un lector de

archivos PDF.
A.5 Ejemplo de Clculo de MR-HR
En esta seccin se muestra como MR-HR realiza los clculos, lo cual puede ser de especial
inters en caso de querer verificar sus resultados. Para esto se utiliz un ejemplo sencillo adjunto
en el Anexo (B.4). Este ejemplo, llamado PruebaPequeo.prj, consta con 6 secciones
transversales, las cuales se muestra en la siguiente Figura:
XXVI
RS = 5010 RS = 4010.*
.0261 .0261 .0261 .0262 .0262 .0262
110 Legend 109 Legend
Ground Ground
109 108
Bank Sta Bank Sta
Elevation (m)
Elevation (m)
108 107
107 106
106 105
105 104
-10 -5 0 5 10 -10 -5 0 5 10
Station (m) Station (m)
RS = 3010.* RS = 2010.*
.0262 .0262 .0262 .0263 .0263 .0263
108 107 Legend
Legend
Ground Ground
107 106
Bank Sta Bank Sta
Elevation (m)
Elevation (m)
106 105
105 104
104 103
103 102
-10 -5 0 5 10 -10 -5 0 5 10
RS = 1010.* RS = 10
.0257 .0257 .0257 .0284 . .0284
106 Legend 105 0 Legend
2
8
Ground 4 Ground
105 104
Bank Sta Bank Sta
Elevation (m)
Elevation (m)
104 103
103 102
102 101
101 100
-10 -5 0 5 10 -10 -5 0 5 10
Figura A-27: Secciones Transversales Utilizadas en Ejemplo

XXVII
De la figura (A-27) se puede apreciar las 6 secciones transversales, de las cuales 5 son
trapezoidales. Esto se realiz por simplicidad, sin embargo, la forma de las secciones
transversales no es un impedimento para la ejecucin de MR-HR. El orden de las secciones
transversales desde aguas arriba hacia aguas abajo son: 5010, 4010.*, 3010.*, 2010.*, 1010.*, 10.
En la figura (A-28) se muestra un esquema del orden de las secciones transversales.
Figura A-28: Esquema de Orden de las Secciones Transversales Utilizadas en Ejemplo

La condicin de borde impuesta aguas abajo es la de altura crtica y el rgimen del flujo es
Subcrtico. Para ms informacin consultar el Anexo (B.4).
El resumen de la configuracin utilizada en MR-HR, se muestra en la figura (A-29).
XXVIII
Figura A-29: Resumen Configuracin del Ejemplo en MR-HR.

Ntese que la configuracin seleccionada fue Parker y Peterson, por ende la

parametrizacin del nmero de Strickler corresponde a la siguiente ecuacin:
0.23845403 7 Rh d s 1

Rh
1/ 6

ds
1 < Rh d s 3.675
2.46 ln 5.5 Rh
d s
R
S t h =
ds
Par ker y Peterson
- 0.00758593 2 R + 0.19590000 0 3.675 < Rh d s 8
h

2
0.00095078 4 R h - 0.02281881 1 R h + 0.25691286 7 8 < Rh d s 12

0.12 Rh d s > 12
(A-1)
XXIX
Los resultados brindados por MR-HR se muestran en la siguientes tablas:
Tabla A-1: Resultados 1ra Iteracin

Diferencia
Caudal Altura Velocidad Radio
N Seccin Manning (n) Manning
(Q) Escurrimiento Promedio Hidrulico N Strickler Rh/ds ds [m]
Iteracin Transversal [s*m^(-1/3)] Iteracin
[m3/s] (h) [m] (v) [m/s] (Rh) [m]
Anterior
5 1 5010 0.03 105.5614 0.8432418 0.5116991 - - 2.558 0.2

5 1 4010.* 0.03 104.4257 1.126402 0.3961863 - - 1.980 0.2
5 1 3010.* 0.03 103.5601 0.8455431 0.5104873 - - 2.552 0.2
5 1 2010.* 0.03 102.639 0.7356053 0.5756877 - - 2.878 0.2
5 1 1010.* 0.03 101.6055 0.778744 0.5481948 - - 2.740 0.2
5 1 10 0.03 101.0981 2.072777 0.4799535 - - 2.399 0.2
En la tabla (A-1) se puede apreciar que en la primera iteracin todos los coeficientes de
Manning son constantes e iguales a 0.030 sm-1/3 (este valor es tpico de cauces naturales de
granulometras gravosas y siempre el mismo para comenzar las iteraciones). En base a esto MR-
HR ejecuta HEC-RAS y obtiene el eje hidrulico asociado a estos coeficientes de Manning junto
con sus respectivos radios hidrulicos, para luego calcular el inverso de la rugosidad relativa
resultante, Rh/ds. Luego, con los resultados de los inversos de la rugosidad relativa de la primera
iteracin es posible calcular los coeficientes de Strickler mediante la ayuda de la ecuacin (A-1).
Ntese que los nmeros de Strickler diferirn entre secciones transversales dado que los inversos
de rugosidades relativas asociados a cada seccin son distintos. Por lo tanto, una vez calculados
los nmeros de Strickler, es posible calcular los coeficientes de Manning asociados a cada
seccin transversal mediante la ecuacin (2.4), la cual es:
g n
St = 1/ 6
ds
Despejando "n" (el coeficiente de Manning) de la ecuacin (2.4) es posible calcular los
nuevos coeficientes de Manning asociados a la segunda iteracin de MR-HR. Con estos nuevos
coeficientes, MR-HR los ingresa a HEC-RAS para el clculo del eje hidrulico y sus respectivos
radios hidrulicos, para volver a calcular los inversos de la rugosidad relativa asociados a cada
seccin transversal. Esto se muestra en la tabla (A-2).
XXX
Tabla A-2: Resultados 2da Iteracin

Diferencia
[m3/s] (h) [m] (v) [m/s] (Rh) [m]
Anterior
5 2 5010 0.0439203 105.6918 0.6759859 0.6186205 1.39E-02 0.17979324 3.093 0.2

5 2 4010.* 0.04659544 104.5605 0.8446297 0.5109695 1.66E-02 0.19074426 2.554 0.2
5 2 3010.* 0.04394235 103.6999 0.6677873 0.6250367 1.39E-02 0.17988352 3.125 0.2
5 2 2010.* 0.04288051 102.8012 0.5778531 0.7054375 1.29E-02 0.17553673 3.527 0.2
5 2 1010.* 0.04329935 101.7112 0.6563797 0.634192 1.33E-02 0.17725129 3.170 0.2
5 2 10 0.04453265 101.0981 2.072777 0.4799535 1.45E-02 0.18229994 2.399 0.2
Con los nuevos inversos de rugosidades relativas se procede en forma anloga hasta
completar las iteraciones totales definidas anteriormente, en este caso son 5.
Tabla A-3: Resultados 3ra, 4ta y 5ta Iteracin
Diferencia
[m3/s] (h) [m] (v) [m/s] (Rh) [m]
Anterior
5 3 5010 0.04229619 105.6752 0.6936757 0.6052205 1.62E-03 0.17314473 3.026 0.2

5 3 4010.* 0.04393357 104.5423 0.8744072 0.4957663 2.66E-03 0.17984754 2.478 0.2
5 3 3010.* 0.04221528 103.6836 0.6847824 0.6118829 1.73E-03 0.17281353 3.059 0.2
5 3 2010.* 0.04131748 102.7844 0.591127 0.692274 1.56E-03 0.16913825 3.461 0.2
5 3 1010.* 0.04210247 101.705 0.6625214 0.6292295 1.20E-03 0.1723517 3.146 0.2
5 3 10 0.04453265 101.0981 2.072777 0.4799535 0.00E+00 0.18229994 2.399 0.2
5 4 5010 0.04247028 105.6773 0.6913764 0.6069289 1.74E-04 0.17385741 3.034 0.2
5 4 4010.* 0.04421835 104.544 0.8714683 0.4972261 2.85E-04 0.18101333 2.486 0.2
5 4 3010.* 0.04238284 103.6848 0.6835253 0.6128367 1.68E-04 0.17349945 3.064 0.2
5 4 2010.* 0.04145128 102.7852 0.5904808 0.6929033 1.34E-04 0.169686 3.464 0.2
5 4 1010.* 0.04216324 101.705 0.6625214 0.6292295 6.08E-05 0.17260048 3.146 0.2
5 4 10 0.04453265 101.0981 2.072777 0.4799535 0.00E+00 0.18229994 2.399 0.2
5 5 5010 0.04244769 105.677 0.6916391 0.6067333 2.26E-05 0.17376493 3.033 0.2
5 5 4010.* 0.04419029 104.5435 0.8723464 0.496789 2.81E-05 0.18089846 2.483 0.2
5 5 3010.* 0.04237047 103.6848 0.6835253 0.6128367 1.24E-05 0.17344879 3.064 0.2
5 5 2010.* 0.04144478 102.7852 0.5904808 0.6929033 6.50E-06 0.16965938 3.464 0.2
5 5 1010.* 0.04216324 101.705 0.6625214 0.6292295 0.00E+00 0.17260048 3.146 0.2
5 5 10 0.04453265 101.0981 2.072777 0.4799535 0.00E+00 0.18229994 2.399 0.2
XXXI
De la tabla (A-3) se puede apreciar que para la ltima iteracin, la nmero 5, las diferencias
entre los coeficientes de Manning de las iteraciones 4 y 5 difieren entre el 5to y 6to decimal. Si se
busca mayor precisin o no se est satisfecho con esas diferencias, es posible aumentar el nmero
de iteraciones para lograr minimizar dichas diferencias.
A.6 Errores Conocidos
El nico error registrado de MR-HR a ocurrido a nivel de ejecucin del cdigo fuente a
travs del compilador de Visual Basic y no desde el archivo ejecutable (*.exe). Este error se
muestra en la figura (A-30) a continuacin:
Figura A-30: Error Asociado a Lectura de Archivos de Salida de HEC-RAS

El error registrado en la figura (A-30) tiene relacin con la lectura del archivo de salida
generado por HEC-RAS. Esto, probablemente, ocurre por conflictos de lectura de archivos de
Windows. Sin embargo, el programa MR-HR no arroja ningn error, puesto que este error se
genera fuera del cdigo fuente de MR-HR, es decir, dentro de HEC-RAS. Se vuelve a reiterar,
que este error ha sido registrado a nivel de ejecucin del cdigo fuente a travs del compilador de
Visual Basic y no desde al archivo ejecutable (*.exe). Para su solucin, basta con re compilar el
cdigo fuente y se evita el error.
XXXII
Patricio Pastenes Araos ANEXO B: Archivos Digitales
Anexo B. Archivos Digitales
Dada la gran extensin de resultados generados por MR-HR, tanto como la escritura del
cdigo fuente, no es factible agregarlo en papel. Es por esto que se adjunta, en el final de esta
memoria10, un DVD-ROM con dicha informacin. La estructura del DVD-ROM se muestra en la
figura (B-1).
Figura B-1: Contenido DVD-ROM con archivos digitales

Se aprecia de la figura (B-1) que dentro del DVD-ROM existen 5 carpetas con los prefijos
B1, B2, B3, B4 y B5. El contenido de estas carpetas se detalla a continuacin:
B.1 MR-HR
Esta carpeta contiene en su interior el programa MR-HR.
B.2 Cdigo Fuente MR-HR
Esta carpeta contiene en su interior el cdigo fuente de MR-HR. Este fue programado en
Visual Basic Express 2010. Sin embargo el cdigo se puede leer con cualquier editor de texto.
10
Los archivos Digitales estn tambin disponibles en el enlace: http://www.ingcivil.uchile.cl/rhma/?page_id=1051
XXXIII
B.3 Ejemplo Aplicacin al Ro Baker
Esta carpeta contiene en su interior el proyecto HEC-RAS con el cual se realiz la

aplicacin al ro Baker. Aqu se puede consultar variables que no considera MR-HR, tanto como
la geometra de cada una de las secciones transversales (740 en total).
B.4 Ejemplo de Clculo Manual de Uso
Esta carpeta contiene en su interior el proyecto HEC-RAS con el cual se realiz la seccin
(A.5), es decir, el ejemplo de clculo de MR-HR.
B.5 Resultados Completos Aplicacin al Ro Baker
Los resultados completos generados por MR-HR, son archivos de texto plano (*.txt) y
tienen la siguiente estructura:
[Parametrizacin] [Ro] [Dimetro Caracterstico Utilizado].[txt]
Por ejemplo, el archivo que contiene los resultados completos del ro Baker, utilizando el
dimetro caracterstico mximo y la parametrizacin de Limerinos, tiene por nombre:
Limerinos Baker MAX.txt
Los archivos de video que dan cuenta de la evolucin espacial del nmero de Strickler en
funcin del caudal son:
Variacin Strickler Tramo Completo.avi: Muestra la variacin espacial del coeficiente de

Strickler en funcin del caudal, desde 1 m3/s a 50 m3/s, con un paso de 1 m3/s y para todo
el tramo del ro Baker estudiado.
XXXIV
Variacin Strickler Tramo Particular.avi: Muestra la variacin espacial del coeficiente de

Strickler en funcin del caudal, desde 1 m3/s a 50 m3/s, con un paso de 1 m3/s y para un
tramo acotado del ro Baker estudiado, comprendido entre los kilmetros 15.5 a 18.5.
XXXV
Patricio Pastenes Araos ANEXO C: HEC-RAS y Visual Basic Express 2010
Anexo C. HEC-RAS y Visual Basic Express 2010
En esta parte se detallan las indexaciones de las diferentes variables proporcionadas por
HEC-RAS v 4.1.0, controlado mediante Visual Basic Express 2010.
La funcin mediante la cual se leen los resultados proporcionados por HEC-RAS a travs
de Visual Basic Express 2010 es la siguiente:
hecras.Output_ReachOutput(riv,rch,prof,nVar,nRS,RS,ChannelDist,Value)11
Donde:
riv := Como Integer asociado a cada ro de un proyecto en particular. Este es

asignado por HEC-RAS a cada ro del proyecto.
rch := Como Integer asociado a cada tramo de cada ro de un proyecto en

particular. Este es asignado por HEC-RAS a cada tramo del proyecto.
prof := Como Integer asociado al perfil de flujo definido previamente en HEC-RAS.
nVar := Como Integer relacionado con qu variable deseo obtener para un ro y

tramo en particular. Su valor vara entre 1 y 266 dependiendo de la variable buscada.
nRS := Como Integer nmero de estaciones de ros o nmero de perfiles

transversales.
RS := Arreglo proporcionado por la funcin con todas las secciones

transversales del ro y tramo predefinidas anteriormente.
ChannelDist := Arreglo proporcionado por la funcin con las distancias entre secciones
transversales.
Value := Arreglo de valores de la variable definida previamente en nVar,

proporcionado por la funcin, asociado a todas las secciones transversales del ro y tramo
11
Para ver la aplicacin de esta funcin, vase Anexo (B.2) el que contiene el cdigo fuente de MR-HR.
XXXVI
predefinidas anteriormente en conjunto con el perfil de flujo para el cual se desea

ejecutar HEC-RAS.
En la siguiente tabla se indica en detalle toda la indexacin de nVar en conjunto con su

descripcin.
Tabla C-1: Indexacin, Significado y Descripcin de Variable nVar.
nVar Significado Descripcin

1 Profile Profile number.
2 W.S. Elev Calculated water surface from energy equation.
3 E.G. Elev Energy gradeline for given WSEL.
4 Max Chl Dpth Maximum main channel depth.
5 Min Ch El Minimum channel elevation.
6 Q Left Flow in left overbank.
7 Q Channel Flow in main channel.
8 Q Right Flow in right overbank.
9 Q Total Total flow in cross section.
10 Flow Area Total area of cross section active flow.
11 Flow Area L Area of left overbank active flow.
12 Flow Area Ch Area of main channel active flow.
13 Flow Area R Area of right overbank active flow.
14 W.P. Total Wetted perimeter of total cross section.
15 W.P. Left Wetted perimeter of left overbank.
16 W.P. Channel Wetted perimeter of main channel.
17 W.P. Right Wetted perimeter of right overbank.
18 Conv. Total Conveyance of total cross section.
19 Conv. Left Conveyance of left overbank.
20 Conv. Chnl Conveyance of main channel.
21 Conv. Right Conveyance of right overbank.
22 Vel Head Velocity head.
23 Vel Total Average velocity of flow in total cross section.
24 Vel Left Average velocity of flow in left overbank.
25 Vel Chnl Average velocity of flow in main channel.
26 Vel Right Average velocity of flow in right overbank.
27 Alpha Alpha - energy weighting coefficient.
28 Beta Beta - momentum weighting coefficient.
29 Top Wdth Act Top width of the wetted cross section, not including ineffective flow.
30 E.G. Slope Slope of the energy grade line at a cross section.
31 Volume Cumulative volume of water from the downstream end of the reach (including ineffective areas).
32 Area Flow area of the entire cross section including ineffective flow.
33 Area Left Flow area of the left overbank including ineffective flow.
34 Area Channel Flow area of the main channel including ineffective flow.
35 Area Right Flow area of the right overbank including ineffective flow.
36 Sta W.S. Lft Left station where water intersects the ground.
37 Sta W.S. Rgt Right station where water intersects the ground.
38 Left Sta Eff Furthest left station where there is effective flow.
XXXVII
Tabla C-1: Indexacin, Significado y Descripcin de Variable nVar. (continuacin 1)
39 Rght Sta Eff Furthest right station that still has effective flow.
40 Length Wtd. Weighted length based on flow distribution, in left bank, channel, and right bank.
41 Length Left Downstream reach length of the left overbank.
42 Length Chnl Downstream reach length of the main channel to next XS (unless BR is d/s, then this is the distance to the deck/roadway).
43 Length Rght Downstream reach length of the right overbank.
44 Mann Wtd Left Conveyance weighted Manning's n for the left overbank.
45 Mann Wtd Chnl Conveyance weighted Manning's n for the main channel.
46 Mann Wtd Rght Conveyance weighted Manning's n for the right overbank.
47 Mann Comp Mannings n value for main channel based on composite roughness equation.
48 Froude # Chl Froude number for the main channel.
49 Froude # XS Froude number for the entire cross section.
50 Trvl Tme Avg Cumulative travel time based on the average velocity of the entire cross section per reach.
51 Trvl Tme Chl Cumulative travel time based on the average velocity of the main channel per reach.
52 Conv. Ratio Ratio of the conveyance of the current cross section to the conveyance of the downstream cross section.
53 Specif Force The specific force for this cross section at the computed water surface elevation.
54 Spc Force PR Specific force prime. For mixed flow, the specific force at this cross section for the flow regime that does not control.
55 W.S. Prime Water surface prime. For mixed flow, the water surface of the flow regime that does not control.
56 Crit W.S. Critical water surface elevation. Water surface corresponding to the minimum energy on the energy versus depth curve.
57 Crit E.G. Critical energy elevation. Minumum energy on the energy versus depth curve.
58 Crit Depth Critical depth. Corresponds to critical water surface.
59 Frctn Loss Friction loss between two cross sections.
60 C & E Loss Contraction or expansion loss between two cross sections.
61 Headloss Total energy loss between two cross sections.
62 Top Width Top width of the wetted cross section.
63 Top W Left Top width of the left overbank. Does not include ìslands', but it does include ineffective flow.
64 Top W Chnl Top width of the main channel. Does not include ìslands', but it does include ineffective flow.
65 Top W Right Top width of the right overbank. Does not include ìslands', but it does include ineffective flow.
66 Num Trials Current number (or final number) of trials attempted before the energy equation is balanced.
Standard step method used to determine WSEL (1 = successful convergence, 2 = minimum error,
67 Std Stp Case 3 = resorted to critical depth).
68 Frctn Slope Representative friction slope between two cross sections.
69 Frctn Slp Md Friction slope averaging method used.
70 Min Error The minimum error, between the calculated and assumed water surfaces when balancing the energy equation.
71 Delta WS Change in water surface through culvert(s) and Bridge(s).
72 Delta EG Change in energy grade line through culvert(s) and Bridge(s).
73 Q Culv Group Flow through all barrels in a culvert.
74 Q Barrel Flow through one barrel in a culvert.
Upstream water surface elevation upstream of bridge, culvert or weir (specific to that opening,
75 W.S. US. not necessarily the energy weighted average).
76 Clv EG No Wr Energy grade elevation at the culvert that was calculated without the weir.
77 E.G. US. Upstream energy grade elevation at bridge or culvert (specific to that opening, not necessarily the weighted average).
78 E.G. IC Upstream energy gradeline based on inlet control.
79 E.G. OC Upstream energy gradeline based on outlet control.
80 Culv Nml Depth Normal depth for this culvert (and flow).
81 Culv Vel DS Velocity in culvert at defined downstream.
82 Culv Vel US Velocity in culvert at defined upstream.
83 Culv Frctn Ls Friction loss through the culvert.
84 Culv Entr Loss Entrance loss (energy loss due only to entrance).
XXXVIII
85 Culv Exit Loss Exit loss (energy loss due to exit).

86 Culv Full Len The length that the culvert flows full.
87 Culv Crt Depth Critical depth inside the culvert.
88 Culv Inv El Up Culvert invert elevation upstream.
89 Culv Inv El Dn Culvert invert elevation downstream.
90 Culv EG Inlet Energy gradeline inside the culvert at the inlet.
91 Culv EG Outlet Energy gradeline inside the culvert at the outlet.
92 Culv WS Inlet Water surface elevation inside the culvert at the inlet.
93 Culv WS Outlet Water surface elevation inside the culvert at the outlet.
94 Q Weir Flow over the weir.
95 Weir Flow Area Area of the flow going over the weir.
96 Weir Sta Lft Station where flow starts on the left side.
97 Weir Sta Rgt Station where flow ends on the right side.
98 Weir Max Depth The maximum depth over the weir.
99 Weir Avg Depth The average depth over the weir.
100 Weir Submerg The ratio of the downstream depth above the weir to the upstream depth above the weir.
101 Min El Weir Flow Elevation where weir flow begins.
102 Wr Top Wdth Top width of water over the weir.
103 Energy/Wr WS Water surface elevation upstream of bridge for low flow energy method and weir flow.
104 Yarnell WS Water surface elevation upstream of bridge for Yarnell method.
105 WSPRO WS Water surface elevation upstream of bridge for the WSPRO method.
106 Prs/Wr WS Water surface elevation upstream of bridge for pressure and/or weir method.
107 Energy WS Water surface elevation upstream of bridge for energy only method.
108 Momen. WS Water surface elevation upstream of bridge for momentum method.
109 Prs O WS Water surface elevation upstream of bridge for pressure only method.
110 Energy/Wr EG Energy grade elevation upstream of bridge for energy method .
111 Yarnell EG Energy grade elevation upstream of bridge for Yarnell method.
112 WSPRO EG Energy grade elevation upstream of bridge for the WSPRO method.
113 Prs/Wr EG Energy grade elevation upstream of bridge for pressure and/or weir method.
114 Energy EG Energy grade elevation upstream of bridge for energy only method.
115 Momen. EG Energy grade elevation upstream of bridge for momentum method.
116 Prs O EG Energy grade elevation upstream of bridge for pressure only method.
117 Br Sel Method Selected bridge method.
118 Min El Prs Elevation at the bridge when pressure flow begins.
119 Crit Num Number of critical depths found.
120 Crit W.S. 1 Water surface elevation of first critical depth.
121 Crit W.S. 2 Water surface elevation of second critical depth.
122 Crit W.S. 3 Water surface elevation of third critical depth.
123 Crit Enrgy 1 Energy associated with first critical depth.
124 Crit Enrgy 2 Energy associated with second critical depth.
125 Crit Enrgy 3 Energy associated with third critical depth.
126 Hydr Depth Hydraulic depth for cross section.
127 Hydr Depth L Hydraulic depth in left over bank.
128 Hydr Depth C Hydraulic depth in channel.
129 Hydr Depth R Hydraulic depth for right over bank.
130 Deck Width Width of Deck.
131 # Barrels Number of barrels in a culvert.
XXXIX
132 Q Bridge Flow through a bridge opening.

Cumulative volume of water in the left overbank from the downstream end of the reach
133 Vol Left (including ineffective areas).
Cumulative volume of water in the channel from the downstream end of the reach
134 Vol Chan (including ineffective areas).
Cumulative volume of water in the right overbank from the downstream end of the reach
135 Vol Right (including ineffective areas).
136 Min El Minimum overall section elevation.
137 Enc Val 1 Target for encroachment analysis.
138 Enc Val 2 Second target for encroachment analysis.
139 Enc Sta L Left station of encroachment.
140 Enc Sta R Right station of encroachment.
141 Dist Center L Distance from center of channel to left encroachment.
142 Dist Center R Distance from center of channel to right encroachment.
143 K Perc L Conveyance reduction from left encroachment.
144 K Perc R Conveyance reduction from right encroachment.
145 Q Perc L Percent of flow in left overbank.
146 Q Perc Chan Percent of flow in main channel.
147 Q Perc R Percent of flow in right overbank.
148 Prof Delta WS Difference in WS between current profile and WS for first profile.
149 Prof Delta EG Difference in EG between current profile and EG for first profile.
150 Shear Total Shear stress in total section.
151 Shear LOB Shear stress in left overbank.
152 Shear Chan Shear stress in main channel.
153 Shear ROB Shear stress in right overbank.
154 Power Total Total stream power.
155 Power LOB Total stream power in left overbank.
156 Power Chan Total stream power in main channel.
157 Power ROB Total stream power in right overbank.
158 Ch Sta L Left station of channel.
159 Ch Sta R Right station of channel.
160 Base WS Water surface for first profile (used in comparison of encroachments).
161 Center Station Center station of main channel.
162 XS Delta WS Change in water surface between current section and next one downstream.
163 XS Delta EG Change in energy gradeline between current section and next one downstream.
164 SA Total Cumulative surface area for entire cross section (including ineffective areas) from the downstream end of the reach.
165 SA Left Cumulative surface area for left overbank (including ineffective areas) from the downstream end of the reach.
166 SA Chan Cumulative surface area for main channel (including ineffective areas) from the downstream end of the reach.
167 SA Right Cumulative surface area for right overbank (including ineffective areas) from the downstream end of the reach.
168 Enc Method Encroachment method.
169 Q Gate Group Flow through all gate openings in a gate group.
170 Gate Open Ht Height of gate opening.
171 Gate #Open The number of gates opened in the current group.
172 Gate Area The flow area in an opened gate.
173 Gate Submerg The ratio of the downstream depth above the gate to the upstream depth above the gate.
174 Gate Invert Gate spillway invert elevation.
175 Q Gates Total flow through all of the gate groups of an inline/lateral structure.
176 BR Open Area Total area of the entire bridge opening.
177 Coef of Q WSPRO bridge method coefficient of discharge.
XL
178 Cum Ch Len Cumulative Channel Length from the downstream end of the reach.
179 Enc WD Encroachment Width.
180 Obs WS Observed Water Surface.
181 WS Air Entr. Water surface elevation accounting for air entrainment.
182 BR Open Vel Average velocity inside the bridge opening (Maximum of BU and BD).
183 Ice Thick LOB Ice thickness in the left overbank.
184 Ice Thick Chan Ice thickness in the main channel.
185 Ice Thick ROB Ice thickness in the right overbank.
186 Ice Vol Total Cummulative volume of ice in an ice jam.
187 Ice Vol. LOB Cummulative volume of ice in the left overbank for an ice jam.
188 Ice Vol. Chan Cummulative volume of ice in the main channel for an ice jam.
189 Ice Vol. ROB Cummulative volume of ice in the right overbank for an ice jam.
190 Ice Top LOB The top elevation of ice in the left overbank.
191 Ice Top Chan The top elevation of ice in the main channel.
192 Ice Top ROB The top elevation of ice in the right overbank.
193 Ice Btm LOB The bottom elevation of ice in the left overbank.
194 Ice Btm Chan The bottom elevation of ice in the main channel.
195 Ice Btm ROB The bottom elevation of ice in the right overbank.
196 Invert Slope The slope from the invert of this cross section to the next cross section downstream.
197 LOB Elev The ground elevation at the left bank of the main channel.
198 ROB Elev The ground elevation at the right bank of the main channel.
199 L. Freeboard The freeboard in the main channel at the left bank.
200 R. Freeboard The freeboard in the main channel at the right bank.
201 Levee El Left The elevation of the left levee.
202 Levee El Right The elevation of the right levee.
203 Ineff El Left The elevation of the left ineffective area.
204 Ineff El Right The elevation of the right ineffective area.
205 L. Levee Frbrd The freeboard before the left levee is over-topped.
206 R. Levee Frbrd The freeboard before the right levee is over-topped.
207 Mann Wtd Total Mannings n value for the total main cross section.
208 Hydr Radius Hydraulic radius for cross section.
209 Hydr Radius L Hydraulic radius in left over bank.
210 Hydr Radius C Hydraulic radius in channel.
211 Hydr Radius R Hydraulic radius for right over bank.
212 Hydr Rad 2/3 Hydraulic radius for cross section to the 2/3 power.
213 W.S. DS Water surface downstream.
214 E.G. DS Energy elevation downstream.
215 Min Weir El Minimum weir elevation.
216 Perc Q Leaving Percentage of flow leaving through a lateral structure.
217 Q US Flow in cross section upstream of a lateral structure.
218 Q DS Flow in cross section downstream of lateral structure.
219 Weir Sta US Upstream station for weir flow starts.
220 Weir Sta DS Downstream station where weir flow ends.
221 Q Leaving Total Total flow leaving in a lateral structure including all gates, culverts and lateral rating curves.
222 SA Min El Minimum elevation of a storage area.
223 SA Area Surface area of a storage area.
224 SA Volume Storage volume of a storage area.
XLI
225 Top W Act Left Top width of the wetted left bank, not including ineffective flow.
226 Top W Act Chan Top width of the wetted channel, not including ineffective flow.
227 Top W Act Right Top width of the wetted right bank, not including ineffective flow.
228 Culv Depth Blocked Depth of fill in a culvert.
229 Culv Inlet Mann n The composite n value at the culvert inlet.
230 Culv Outlet Mann n The composite n value at the culvert outlet.
231 Ice WS Err Convergence error in water surface for dynamic ice jam.
232 Ice Err Convergence error in ice thickness for dynamic ice jam.
233 Piping Flow Flow from piping weir failure.
234 Breach CL Center line of weir breach.
235 Breach WD Bottom width of weir breach.
236 Breach Bottom El Bottom Elevation of weir breach.
237 Breach Top El Top Elevation of weir breach.
238 Breach SSL Left side slope of weir breach.
239 Breach SSR Right side slope of weir breach.
240 Q Pump Group Pump group flow.
241 Q Lat RC Lateral rating curve flow.
242 Q Culv Total flow in all culvert groups.
243 Culv Length Length of the culvert barrel.
244 Q Pump Station Total flow in all pump groups in a pump station.
245 WS Inlet WS at the inlet of a pump station.
246 WS Outlet WS at the outlet of a pump station.
247 Pumping Head Pumping head for the pump station.
248 Inflow Net inflow into a storage area.
249 Outflow Net outflow into a storage area.
250 Net Flux Net inflow - outflow for a storage area.
251 Enc Offset L Minimum setback from the left overbank station.
252 Enc Offset R Minimum setback from the right overbank station.
253 Min Ch Pilot Minimum channel elevation (including pilot channels).
254 Diff Difference between the previous two columns.
255 Min Ch El Sta Station of the minimum channel elevation.
256 Culv Area DS Cross sectional flow area in culvert at defined downstream.
257 Culv Area US Cross sectional flow area in culvert at defined upstream.
258 Gate Weir Coef Coefficient used in weir flow over the gate.
259 Weir Coef Coefficient used in weir flow.
260 Q Breach Flow through a breach.
261 Breach Avg Velocity Average flow velocity through a breach.
262 Breach Flow Area Flow area through a breach.
263 Left Station Left station of the cross section.
264 Right Station Right station of the cross section.
265 Levee Sta Left Left levee station.
266 Levee Sta Right Right levee station.
Fuente: Elaboracin Propia en Base a HEC-RAS v 4.1.0
XLII

Cálculo de Ejes Hidráulicos en Lechos Macro Rugosos - Patricio Pastene

Cargado por

Copyright:

Formatos disponibles

Cálculo de Ejes Hidráulicos en Lechos Macro Rugosos - Patricio Pastene

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Cálculo de Ejes Hidráulicos en Lechos Macro Rugosos - Patricio Pastene

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD DE CHILE

FACULTAD DE CIENCIAS FSICAS Y MATEMTICAS

DISEO E IMPLEMENTACIN DE SOFTWARE DE ACOPLE A HEC-

MEMORIA PARA OPTAR AL TTULO DE INGENIERO CIVIL

DISEO E IMPLEMENTACIN DE SOFTWARE DE ACOPLE A HEC-RAS:

En la presente memoria se brinda una herramienta computacional de acople a HEC-RAS,

La primera etapa del estudio, consisti en la recopilacin de antecedentes generalizados de

Los resultados de la primera etapa sugirieron que el desarrollo de ejes hidrulicos

A mi familia y polola por su incesable apoyo a lo largo de toda mi formacin.

A mi profesor gua, Sr. Alberto de la Fuente S., por su dedicacin, comprensin,

CAPTULO 1. INTRODUCCIN ...................................................................................1

CAPTULO 2. REVISIN DE ANTECEDENTES .......................................................5

2.1 Prdida de Carga en Lechos Macro Rugosos ............................................................5

CAPTULO 3. IMPLEMENTACIN DE LA HERRAMIENTA: MR-HR .............26

3.1 Ajuste de Parametrizaciones....................................................................................26

CAPTULO 4. APLICACIN AL RO BAKER. ........................................................43

4.1 Antecedentes Generales Ro Baker .........................................................................43

CAPTULO 5. DISCUSIN Y ALCANCES ................................................................70

A.1 Generalidades ..............................................................................................................I

ANEXO B. ARCHIVOS DIGITALES ...................................................................XXXIII

B.1 MR-HR ............................................................................................................ XXXIII

ANEXO C. HEC-RAS Y VISUAL BASIC EXPRESS 2010 .................................XXXVI

La macro rugosidad se caracteriza por la presencia de grandes partculas de sedimento que

Para el clculo de ejes hidrulicos en lechos de granulometras gruesas existen en general

Ambos enfoques tienen asociados grados de subjetividad. En el primero se presenta la

desprenden una serie de consecuencias asociadas al clculo adecuado de coeficientes de

1.2 Objetivos y Alcances

1.2.1 Objetivo Principal

El objetivo principal de este trabajo es facilitar, mediante el diseo e implementacin de un

1.2.2 Objetivos Especficos

Disear e implementar un software capaz de acoplarse al software HEC-RAS y realizar

Realizar una aplicacin real a un cauce de secciones transversales conocidas y de lecho

No modifica el cdigo fuente de HEC-RAS, sino que lo llama mediante Windows.

Resuelve problemas en rgimen permanente, con flujos turbulentos plenamente

1.3 Contenido del Informe

Captulo 2. Revisin de Antecedentes

2.1 Prdida de Carga en Lechos Macro Rugosos

Rh : corresponde al radi hidrulico.

Con el propsito de facilitar la presentacin de los antecedentes recopilados, a continuacin

Rh : corresponde al radio hidrulico. [m]

ds : corresponde al dimetro caracterstico del sedimento grueso, tpicamente D90

S : corresponde a la pendiente del terreno. [m m]

U : corresponde a la velocidad media del escurrimiento. [m s ]

u* : corresponde a la velocidad friccional o de corte. [m s ]

Cf : corresponde a un coeficiente o factor de friccin adimensional.

Cchezy : corresponde al coeficiente de Chezy. [ m1 / 2 s ]

y : corresponde a la altura del escurrimiento. [m]

g : corresponde a la aceleracin de gravedad. m s 2 [ ]

: corresponde a una constante universal caracterstica de la turbulencia o tambin

Para poder caracterizar escurrimientos en lechos macro rugosos, es necesaria la

De acuerdo a la ecuacin de Strickler, el coeficiente de Manning depende del tamao del

En base a la relacin (2.3), se define un parmetro adimensional que relaciona el

Reemplazando la ecuacin (2.3) en la ecuacin (2.4), resulta que el nmero de Strickler es

Diversos investigadores han estudiado escurrimientos en lechos macro rugosos (Keulegan,

Por lo tanto, manipulando las ecuaciones (2.5) y (2.6) se obtiene que:

Luego, considerando equilibrio de fuerzas en canales abiertos, para flujos permanentes y

Adems, considerando las mismas condiciones anteriores, es decir, flujo permanente y

Reemplazando los resultados de las ecuaciones (2.10) y (2.9) en la ecuacin (2.7) se

Finalmente, es posible sustituir la ecuacin (2.4) en la ecuacin (2.11), obtenindose: