Clase 04 Práctica HEC RAS

Descargar como ppt, pdf o txt
Descargar como ppt, pdf o txt
Está en la página 1de 70

ESPECIALIZACIÓN

EN
INGENIERÍA VIAL
FACULTAD DE INGENIERIA
UNIVERSIDAD NACIONAL DE CUYO

“PRACTICA de Modelación Hidráulica


de Ríos con HEC RAS”
1
Ing. Victor H. Burgos
Temario
1 Introducción. Conceptos. Hipótesis
2 Características generales de HEC RAS
3 Desarrollo de una simulación con HEC RAS
4 Estructuras: Análisis de puentes / Alcant. / Defensas
5 Erosión local en puentes
6 Ejercicios
7 Flujo impermanente
8 HEC Geo RAS
9 Conclusiones
10 Bibliografía
Ing. Victor H. Burgos 2
INA / UTN
Introducción 1

Hipótesis de partida
• HEC-RAS (Hydrologic Engineering Center - River Analysis System) es una
aproximación UNIDIMENSIONAL al flujo en lámina libre.
• Resuelve el FGV a partir de Ecuac. de BALANCE de
ENERGIA (trinomio de Bernoulli) entre dos secciones dadas,
excepto para estructuras (puentes, vertederos, etc) donde
utiliza Ec. de conservac. de CANTIDAD DE
MOVIMIENTO o Ec. Empíricas (erosión)
• 1D: La única componente de velocidad considerada es en
la dirección del flujo. Dir. vert. y transv. se desprecian.
• Sólo p/ pendientes pequeñas < 1:10 (vert/horiz)
Ing. Victor H. Burgos 3
INA / UTN
Introducción
Ecuaciones básicas (Ec.1)
teorema de Bernoulli
conservación de la energía.
Pérdidas de Energía: 1V12  2V22
Ec. Manning he  LS f  C  (Ec.2)
Coef. contraccion/expansion 2g 2g

Ing. Victor H. Burgos 4


INA / UTN
Introducción
Coeficiente de Manning

Ing. Victor H. Burgos 5


INA / UTN
Introducción
Ecuaciones básicas
Conservación (utilizada cuando ocurren cambios de
Cant. Movimiento: regimen, p.e. puentes y confluencias)

P2  P1  Wx  F f  Q· ·V
Donde:
Pi: Presión hidrostática en la sección i
Wx: Peso del agua en la dirección x
Ff: Fuerza de fricción
Q: Caudal
: Densidad del agua
Vx: Cambio de velocidad de la sección 2 a la 1en la dirección x

Ing. Victor H. Burgos 6


INA / UTN
Introducción
Modelos
• Caracterización de la realidad
• Concepto de modelo
• Tipos:
• Solución analítica / aproximación 12

11

• Empirismo 10

Cr i t i ca l C
• Modelación Numérica o Física
9
Upper Reach
8

r.
7

Ing. Victor H. Burgos 2 7


INA / UTN
1
Características Generales de HEC RAS 2

Capacidades:
• Cálculo 1D en régimen permanente e impermanente
• P/ Cauces naturales (ríos, estuarios,etc.) y canales
artificiales (riego, derivación, etc.)
• P/ análisis de efectos de estructuras: compuertas,
alcantarillas, puentes, etc.
Propósitos:
• Dimensionamiento de Obras Hidráulicas
• Simulación hidráulica
• Uso en condiciones nominales
• Uso en eventos extremos
Ing. Victor H. Burgos 8
INA / UTN
Características Generales de HEC RAS

Estructuras
1. Canalización
2. Puente
3. Alcantarilla
4. Estructuras Laterales
5. Defensas

Ing. Victor H. Burgos 9


INA / UTN
Características Generales de HEC RAS

La pérdida de carga (he) es considerada como la suma de


Ecuaciones básicas la pérdida friccional entre dos secciones y la diferencia
(Ec.1) entre cargas de velocidad entre secciones ponderadas por
un coeficiente de pérdida por contracción o expansión.

V V
2 2
he  LS f  C 2 2  1 1
2g 2g

(Ec.2)

(Ec.3)
Donde:
L = largo representativo del tramo
Sf = pendiente de fricción representativa entre dos secciones
C = coeficiente de pérdida por expansión o contracción

L se calcula como la suma de las distancia entre las subsecciones (planicies


laterales y canal principal) ponderadas por sus respectivos caudales
promedios y dividido por el caudal promedio total

Ing. Victor H. Burgos 10


INA / UTN
Características Generales de HEC RAS

Ecuaciones básicas
Dado que el modelo es unidimensional es necesario utilizar solo
una altura de velocidad, este coeficiente compatibiliza las diferencias
entre las alturas de velocidad de las distintas subsecciones

Se define en función de las áreas de flujo y la capacidad de


conducción

 At  
2   K lob  3

 K ch  3

 K rob  3

2 
 A
 lob  2
 A  2
 A  
 ch rob

 Kt  3

Ing. Victor H. Burgos 11


INA / UTN
Características Generales de HEC RAS

Ecuaciones básicas
Se define la capacidad de conducción de una subsección como el
caudal dividido por la raíz de la pérdida de carga, y se introduce
en la ecuación de Manning.
Q
K (m3 / s) 1
Sf K AR 2 / 3 (m 3 / s )
n

donde:
K = capacidad de transporte de una subdivisión
n = coeficiente de Manning de una subdivisión
A = área de flujo de una subdivisión
R = radio hidráulico de una subdivisión

Ing. Victor H. Burgos 12


INA / UTN
Características Generales de HEC RAS

Subsecciones
Por defecto HEC-RAS considera subsecciones en donde existe
variación del coeficiente de rugosidad de Manning, sin embargo, por
compatibilidad con HEC-2 puede modificarse la definición de
subsecciones en cada nodo del perfil (Estilo HEC-2)
La conducción total de una sección será la suma de la conducción de
las subsecciones

Ing. Victor H. Burgos 13


INA / UTN
Características Generales de HEC RAS

Proceso computacional (reg.permanente)

Ing. Victor H. Burgos 14


INA / UTN
Características Generales de HEC RAS

Módulos y Componentes
• Módulo 1: Cálculo de perfil de lámina de agua en REGIMEN
PERMANENTE (p/ simular tramos simples o dendríticos o red de
canales). Pérdidas de energía por Manning y por Coef. Contracc. y
Expansión
• Módulo 2: Simulación en REGIMEN IMPERMANENTE
• Módulo 3: Estimación de transporte de sedimentos en lecho móvil.
• SubMódulo: Erosión local / Diseño Hidráulico
• Módulo 4: Análisis de Calidad
• Módulo SIG: RAS MAPPER
• 3 Componentes: Geométrico / Datos de flujo / Plan (Escenario)

Ing. Victor H. Burgos 15


INA / UTN
Características Generales de HEC RAS

Los datos de flujo permanente necesarios para calcular el


perfil de la superficie del agua. Son básicamente:
• Régimen del flujo
• Condiciones de borde
• Caudales del sistema

El modelo realiza sus cálculos a partir de una sección


transversal con condiciones iniciales conocidas o supuestas
y siguen aguas arriba para flujo subcrítico o aguas abajo
para flujo supercrítico.

Ing. Victor H. Burgos 16


INA / UTN
Características Generales de HEC RAS

Para régimen subcrítico el modelo restringe como


resultados correctos valores mayores a la altura crítica.

Para régimen supercrítico el modelo restringe como


resultados correctos valores menores a la altura crítica.

En casos donde el régimen de flujo cambia de subcrítico a


supercrítico, o viceversa, el usuario debe ejecutar el
programa en modo de régimen de flujo mixto

Ing. Victor H. Burgos 17


INA / UTN
Características Generales de HEC RAS

Condiciones de Borde
1. Nivel conocido de agua
2. Tirante crítico
3. Altura normal (pendiente de energía)
4. Curva de descarga (H-Q)

p/ flujo subcrítico: condic.borde aguas abajo


p/ flujo supercrítico: condic.borde aguas arriba
p/flujo mixto: ambas

Ing. Victor H. Burgos 18


INA / UTN
Características Generales de HEC RAS

Módulos y Componentes MODULOS

COMPONENTES:
Escenario

Geometría

Tipo de
flujo

Íconos de acceso
Ing. Victor H. Burgos 19
INA / UTN
Características Generales de HEC RAS

Módulos y Componentes Ver secciones


Íconos de acceso: Ver resultados
Ver Perfil long.
Guardar proyecto tabulares
Ver Curvas
Abrir proyecto Geometría Reg. de flujo Ver 3D
Mod.Sedimentos Ver hidrogramas
Mod.Calidad /limn Ver
database

Ver mensajes
error / avisos

Unidades
Correr Diseño Hidr
Correr Reg. Correr A.Calidad
Correr Reg. Impermanente Archivos
Permanente Correr A.Sedimentos
Ing. Victor H. Burgos 20
INA / UTN
Características Generales de HEC RAS

Instalación y configuración inicial

Ing. Victor H. Burgos 21


INA / UTN
Características Generales de HEC RAS

Unidades

Ing. Victor H. Burgos 22


INA / UTN
Desarrollo de una simulación con HEC RAS 3

Ing. Victor H. Burgos 23


INA / UTN
Desarrollo de una simulación con HEC RAS

Ing. Victor H. Burgos 24


INA / UTN
Desarrollo de una simulación con HEC RAS

Esquema de río
Morfología de la red
Conectividad entre
cauces
Definición de puntos de
intersección
Nomenclatura (Cauces,
tramos e intersecciones)
Secciones transversales
Topografía del cauce
Distancia entre
secciones
Singularidades
Coeficientes de pérdida
de energía Ing. Victor H. Burgos 25
INA / UTN
Desarrollo de una simulación con HEC RAS

Ing. Victor H. Burgos 26


INA / UTN
Desarrollo de una simulación con HEC RAS

Ing. Victor H. Burgos 27


INA / UTN
Desarrollo de una simulación con HEC RAS

Ing. Victor H. Burgos 28


INA / UTN
Desarrollo de una simulación con HEC RAS

Ing. Victor H. Burgos 29


INA / UTN
Desarrollo de una simulación con HEC RAS

Ing. Victor H. Burgos 30


INA / UTN
Desarrollo de una simulación con HEC RAS

Ing. Victor H. Burgos 31


INA / UTN
Desarrollo de una simulación con HEC RAS

Ing. Victor H. Burgos 32


INA / UTN
Desarrollo de una simulación con HEC RAS

Ing. Victor H. Burgos 33


INA / UTN
Desarrollo de una simulación con HEC RAS

Ing. Victor H. Burgos 34


INA / UTN
Desarrollo de una simulación con HEC RAS

Ing. Victor H. Burgos 35


INA / UTN
Desarrollo de una simulación con HEC RAS

Ing. Victor H. Burgos 36


INA / UTN
Desarrollo de una simulación con HEC RAS

Ing. Victor H. Burgos 37


INA / UTN
Desarrollo de una simulación con HEC RAS

Legend

WS 22FEB1999 0500
Ground
Bank Sta
Ground
Levee
Ineff

Ing. Victor H. Burgos 38


INA / UTN
Desarrollo de una simulación con HEC RAS

Flujo perpendicular a la sección


Cálculos SOLO en secciones.
Geometría para definir áreas de flujo y perímetro mojado

Ing. Victor H. Burgos 39


INA / UTN
Desarrollo de una simulación con HEC RAS

Ing. Victor H. Burgos 40


INA / UTN
Estructuras: Análisis de puentes 4

Ing. Victor H. Burgos 41


INA / UTN
Estructuras: Análisis de puentes

- Sección 1: A una distancia L aguas arriba del puente.


- Sección 2: Sección aguas arriba del puente.
- Sección 3: Sección aguas abajo del puente.
- Sección 4: A una distancia 4 L aguas abajo del puente.

Ing. Victor H. Burgos 42


INA / UTN
Estructuras: Análisis de puentes

En algunos puentes los estribos son muy cortos longitudinalmente,


por lo que se puede optar por situar las secciones 1 y 4 a distancias
de 0,5 B y 2 B respectivamente

Ing. Victor H. Burgos 43


INA / UTN
Estructuras: Alcantarillas

Ing. Victor H. Burgos 44


INA / UTN
Estructuras: Alcantarillas

Ing. Victor H. Burgos 45


INA / UTN
Estructuras Laterales

Ing. Victor H. Burgos 46


INA / UTN
Estructuras: Defensas (Levees)

Ing. Victor H. Burgos 47


INA / UTN
Erosión local en puentes 5

Una vez corrido (y calibrado) el modelo hidráulico se pueden


estimar con el módulo Hydraulic Design Function/ Scour at
Bridges las erosiones locales y por contracción.

La erosión por contracción puede ser analizada con la ec. de


Laursen para agua clara (1963) o para lecho móvil (1960)

Se debe ingresar K1 (coef. temperatura) y d50

La erosión vendría dada por tanto por: e=y2-y1

Ing. Victor H. Burgos 48


INA / UTN
Erosión local en puentes

La erosión local por pilas puede ser analizada con la ec. FHWA,
Richarson, 1975

Se debe ingresar K1 (coef. forma de pila) y d50


Donde:
e= erosión máxima (m)
B= anchura de la pila (m)
k1= constante de forma de la pila (1.0 para la pila circular, 1.1 para la
pila rectangular)
k2= constante de ángulo de ataque, que puede omitirse si se usa la
anchura B* de la pila proyectada perpendicularmente a la corriente en
lugar de B.
y1= calado aguas arriba
Fr1= número de Froude aguas arriba
Ing. Victor H. Burgos 49
INA / UTN
EJERCICIOS 6
 Ejercicio 1: Estimar el tirante y la velocidad de flujo en el badén de
RP86 (Tupungato) en función de los resultados obtenidos por
modelación hidrológica de la cuenca del Aº Anchayuyo para una
recurrencia de 50 años. 2) Dibujar el tramo
1) cargar imagen y vectores

Aº Anchayuyo Baden

3) Cargar secciones (con tabla anexa)


Villa Bastías
(TUP)
Ing. Victor H. Burgos 50
INA / UTN
EJERCICIOS
 Ejercicio 1
Agregar secciones (con Option/Add new cross section) indicando la
progresiva (seg/tabla anexa)
Cargar las progresivas y cotas con tabla anexa (copiando y pegando)
Cargar los coef de rugosidad de Manning (seg/tabla n)
Cargar las distancias hacia aguas abajo (seg/ tabla anexa)

Ing. Victor H. Burgos 51


INA / UTN
EJERCICIOS
 Ejercicio 1
El perfil longitudinal debería quedar así:

Con botón derecho se pueden medir pendientes (para condiciones de borde)

Ing. Victor H. Burgos 52


INA / UTN
EJERCICIOS
 Ejercicio 1
 Cargar las condiciones de flujo (permanente) según caudales
resultantes, indicando las condiciones de borde, y correr el escenario
(plan)
 Utilizar distribución de velocidades
 Analizar resultados

Ing. Victor H. Burgos 53


INA / UTN
EJERCICIOS
 Ejercicio 1

Análisis de
secciones

Ing. Victor H. Burgos 54


INA / UTN
EJERCICIOS
 Ejercicio 1

Análisis de
Tablas

Ing. Victor H. Burgos 55


INA / UTN
EJERCICIOS
 Ejercicio 1

Análisis de
Tablas

Ing. Victor H. Burgos 56


INA / UTN
EJERCICIOS
Ejercicio 2:
Estimar la erosión local por pilas y contracción del puente Tierras
Blancas sobre RN40, utilizando los resultados de la modelación
hidrológica de la cuenca interviniente para una recurrencia de 50 años.

Ing. Victor H. Burgos 57


INA / UTN
EJERCICIOS
Ejercicio 2: Datos:
22.85 22.85 22.85
CR:897.05
9.17

D=1.20m
8.39 CFV:895.39
887.0

874.4 12.6 874.4 D=1.50m

13.25

Ing. Victor H. Burgos 58


INA / UTN
EJERCICIOS
Ejercicio 2:
Luego de cargar la imagen satelital (con botón de
mapita) y trazar el eje del arroyo (con botón
River Reach), ir a botón cross section y luego
Option / Add a new cross section, indicando la
Progresiva (Station) de la planilla.

Cargar las secciones transversales copiando y


pegando en la ventana de edición de secciones
(Cross Section Data)

Colocar las distancias hacia aguas abajo

Agregar los coef. Manning (seg/tabla n)

Ing. Victor H. Burgos 59


INA / UTN
EJERCICIOS
Ejercicio 2: Agregando el puente
Desde Geometría, con botón Brdg/Culv

Agregar Puente:
Option/Add Bridge
(Progresiva entre las dos
secciones más próximas
al puente)

Botones de geometría
del puente. Indicar
ancho tablero, cotas
rasante, fondo viga,
pilas, estribos, etc.

Ing. Victor H. Burgos 60


INA / UTN
EJERCICIOS
Ejercicio 2:
Debería quedar así el perfil longitudinal

En componente de flujo,
elegimos flujo perm,
e indicamos condiciones de
borde y caudal a simular.
(Desde perfil long. Podemos
medir pendientes con botón
derecho) Ing. Victor H. Burgos 61
INA / UTN
EJERCICIOS
Ejercicio 2:
Corremos la simulación y analizamos resultados.

Ing. Victor H. Burgos 62


INA / UTN
EJERCICIOS
Ejercicio 2:
Con módulo de análisis hidráulico, estimamos erosión local.
Se anexa estudio de suelos (obtener d50, d95)
Agregar coef. de forma de pilas, temperatura y coef. acorazamiento

Ing. Victor H. Burgos 63


INA / UTN
EJERCICIOS
Ejercicio 2:
Calculamos con compute y analizamos reporte

Hydraulic Design Data


Contraction Scour
Left Channel Right
Input Data
Average Depth (m): 0.72 6.39 0.14
Approach Velocity (m/s):0.56 3.48 0.21
Br Average Depth (m): 4.53
Ing. Victor H. Burgos 64
BR Opening Flow (m3/s): 1647.00
INA / UTN
BR Top WD (m): 54.68
EJERCICIOS
Ejercicio 2: Results
Scour Depth Ys (m): 4.67
Hydraulic Design Data Froude #: 0.70
Contraction Scour Equation: CSU equation
Left Channel Right
Input Data Abutment Scour
Average Depth (m): 0.72 6.39 0.14 Left Right
Approach Velocity (m/s): 0.56 3.48 Input 0.21
Data
Br Average Depth (m): 4.53 Station at Toe (m): 166.46 227.05
BR Opening Flow (m3/s): 1647.00 Toe Sta at appr (m): 156.46 220.41
BR Top WD (m): 54.68 Abutment Length (m): 48.62 198.45
Grain Size D50 (mm): 3.20 3.20 3.20 Depth at Toe (m):0.30 0.30
Approach Flow (m3/s): 17.13 1624.27 5.60 K1 Shape Coef: 1.00 - Vertical abutment
Approach Top WD (m): 42.97 72.93 195.11 Degree of Skew (degrees): 90.00 90.00
K1 Coefficient: 0.590 0.640 0.590 K2 Skew Coef: 1.00 1.00
Results Projected Length L' (m): 48.62 198.45
Scour Depth Ys (m): 3.25 Avg Depth Obstructed Ya (m): 1.10 0.19
Critical Velocity (m/s): Flow Obstructed Qe (m3/s): 70.32 21.82
Equation: Live Area Obstructed Ae (m2): 53.71 38.11
Results
Pier Scour Scour Depth Ys (m): 1.72 1.98
All piers have the same scour depth Froude #: 0.49 0.74
Input Data Equation: HIRE HIRE
Pier Shape: Round nose
Pier Width (m): 1.50 Combined Scour Depths
Grain Size D50 (mm): 3.20000
Depth Upstream (m): 6.30 Pier Scour + Contraction Scour (m):
Velocity Upstream (m/s): 5.49 Channel: 7.92
K1 Nose Shape: 1.00
Pier Angle: 0.00 Left abutment scour + contraction scour (m): 4.97
Pier Length (m): 13.25 Right abutment scour + contraction scour (m): 5.22
K2 Angle Coef: 1.00
K3 Bed Cond Coef: 1.10
Grain Size D90 (mm): 15.00000
K4 Armouring Coef: 1.00
Ing. Victor H. Burgos 65
INA / UTN
Flujo impermanente 7

HEC-RAS permite la simulación del flujo en régimen crítico y


subcrítico en la misma simulación.

Para ello utiliza el método LPI (Local Partial Inertia), que


añade un término a las ecuaciones para reducir los términos
de inercia cuando el número de Froude se acerca a la
unidad y por tanto habrá un cambio de régimen y una posible
discontinuidad.

Este método hace posible que el modelo no se vuelva


inestable ante una discontinuidad del flujo a costa de una
pérdida de precisión de resultados donde ésta se produzca.

Ing. Victor H. Burgos 66


INA / UTN
HEC Geo RAS 8
El modelo tiene capacidad para interactuar con modelos digitales de
terreno (DEM o TIN) mediante entorno GIS.

Para trabajar con HEC-RAS existe una aplicación, HEC-GeoRAS,


que trabaja en entorno GIS (ArcInfo, ArcView o ArcGIS). A partir de
un modelo digital del terreno puede crearse un archivo para importar
a HEC-RAS la información geométrica necesaria para llevar a cabo
una simulación (cauce, secciones, etc.).

De la misma forma después de haber realizado una simulación con


HEC-RAS se pueden pasar los resultados al entorno GIS para
realizar mapas de inundación.

Ing. Victor H. Burgos 67


INA / UTN
CONCLUSIONES 9
 HEC-RAS es un modelo unidimensional, sin embargo incorpora la
posibilidad de realizar simulaciones cuasibidimensionales
permitiendo la representación de zonas donde se considera que no
hay flujo o con zonas de almacenamiento de agua.
 Los resultados en tramos donde el flujo sea marcadamente
bidimensional (grandes llanuras de inundación) no serán todo lo
rigurosos que sería deseable.
 Para el estudio de los niveles y velocidades de agua en ríos, la
aproximación que sin duda se ha utilizado más es la de flujo
unidimensional y régimen permanente gradualmente variado.
• ……….

Ing. Victor H. Burgos 68


INA / UTN
BIBLIOGRAFIA 10

• E. Bladé, M. Sánchez–Juny, H.P. Sánchez, D. Niñerola, M. Gómez (2009) “Modelación


numérica en ríos en régimen permanente y variable” Ediciones UPC, 2009
• Brunner, G. W. (1995). HEC-RAS River Analysis System. Hydraulic Reference Manual. Version
1.0. HYDROLOGIC ENGINEERING CENTER DAVIS CA.
• Brunner, G. W. (2002). HEC-RAS River Analysis System: User's Manual. US Army Corps of
Engineers, Institute for Water Resources, Hydrologic Engineering Center.
• Richardson, E. V., & Davis, S. R. (2001). HEC-18: Evaluating Scour at Bridges. FHWA
Hydraulic Engineering Circulars (HEC).

Ing. Victor H. Burgos 69


INA / UTN
BIBLIOGRAFIA

MUCHAS GRACIAS
POR SU ATENCION

Datos de contacto:

Victor Burgos
vhburgos@yahoo.com.ar
vburgos@ina.gob.ar

Ing. Victor H. Burgos 70


INA / UTN

También podría gustarte