Modelo Tridimensional Del Flujo en Pilas de Puente 01 PDF

UNIVERSIDAD FRANCISCO DE PAULA SANTANDER
BIBLIOTECA EDUARDO COTE LAMUS

RESUMEN TESIS DE GRADO
AUTORES: IVAN HERNANDO RAMIREZ MENDOZA

SAILY ANTONIA TORRES ROSADO
FACULTAD: INGENIERA
PLAN DE ESTUDIOS: INGENIERA CIVIL
DIRECTOR: MEIMER PEARANDA CARRILLO
TITULO DE LA TESIS: MODELO TRIDIMENSIONAL DEL FLUJO DE AGUA EN PILAS

DE PUENTE POR SIMULACIN NUMRICA
.
RESUMEN:
El propsito del siguiente trabajo es modelar el flujo de agua en los tipos

geomtricos comunes de pila, como lo son; rectangular, circular y en punta de
diamante empleando ANSYS 10.0, efectuando el anlisis de la actuacin hidrulica
del cauce en el entorno de las pilas y estimando la magnitud de la socavacin local
en cada geometra estudiada, extendiendo el trabajo al proponer una seccin con
forma hidrodinmica que perturbe en menor grado los parmetros hidrulicos del
flujo en el cauce
CARACTERSTICAS
PAGINAS:
141
PLANOS:
ILUSTRACIONES:
CD ROM: 1
MODELO TRIDIMENSIONAL DEL FLUJO DE AGUA EN PILAS DE PUENTES

POR SIMULACIN NUMRICA
IVAN HERNANDO RAMIREZ MENDOZA


FACULTAD DE INGENIERA
PLAN DE ESTUDIOS DE INGENIERIA CIVL
SAN JOS DE CCUTA
2008
MODELO TRIDIMENSIONAL DEL FLUJO DE AGUA EN PILAS DE PUENTES

POR SIMULACIN NUMRICA
IVAN HERNANDO RAMIREZ MENDOZA

Trabajo de grado presentado como requisito para optar al ttulo de

Ingeniero civil
Director
MEIMER PEARANDA CARRILLO
Ingeniero Mecnico

FACULTAD DE INGENIERA
PLAN DE ESTUDIOS DE INGENIERIA CIVL
SAN JOS DE CCUTA
2008
A mi madre Nhora Esmeralda Mendoza Sanchez por sus oraciones y

pensamientos desde la distancia, dndome el calor de madre que aunque no era
fsico posee la misma intensidad, a mi madre doy gracias por ser el eje
fundamental de mi vida.
A mi padre Pedro Ivan Ramirez Jaimes por ensearme afrontar los problemas con
tranquilidad, igualmente expreso un voluminoso agradecimiento a cada integrante
de mi familia y amistades, que colaboraron con mi noble causa en el transcurso de
mi carrera.
Ivan
A Mariluz Rosado Mendoza porque nunca me ha fallado como madre y padre a la

vez; ya que ese papel lo ha sabido cumplir muy bien, brindndome todo el amor,
cario, compresin del mundo y motivacin de seguir adelante.
Saily
AGRADECIMIENTOS
Los autores del trabajo expresan sus agradecimientos a:
Al semillero de investigacin ACUA, en el cual surgi la idea de centrar nuestra

investigacin, proporcionando buenas bases para comenzar la aplicacin de
software en el estudio hidrulico de las pila de puente.
Al M.Sc. Gustavo Adolfo Carrillo por enfocarnos en el tema especifico de nuestro

trabajo de grado, brindndonos informacin y asesoras adecuada en la
realizacin de la investigacin, inculcndonos que la ingeniera civil es un medio
de continuo aprendizaje.
Al Ingeniero Mecnico Meimer Pearanda Carrillo por sus sugerencias brindadas

para el desarrollo de este trabajo.
Al PhD Carlos Acevedo por asignarnos el espacio en la sala del grupo

investigacin GIDIMA donde se realizo nuestro proyecto de grado .
CONTENIDO
Pg.
INTRODUCCIN
25
1. ANTECEDENTES
26
2. FUNDAMENTOS HIDRAULICOS
30
2.1 FLUJO TRIDIMENSIONAL
30
2.2 PROPIEDADES DEL FLUJO
31
2.2.1 Velocidad
31
2.3 PROPIEDADES DEL FLUIDO
31
2.3.1 Densidad
31
2.3.2 Viscosidad dinmica
31
2.4 TIPOS DE FLUJO
31
2.4.1 Flujo libre
31
2.4.2 Flujo variado permanente
32
2.5 RGIMEN DEL FLUJO
32
2.5.1 Flujo con rgimen turbulento
32
2.5.2 Flujo incompresible
33
2.6 FUNDAMENTOS DE VISUALIZACIN DEL FLUJO
33
2.6.1 Punto de estancamiento
33
2.6.2 Equilibrio
33
2.6.3 Lneas de corriente
34
2.6.4 Graficas de contornos
34
2.7 VORTICIDAD Y ROTACIONALIDAD
35
2.8 FUERZA DE ARRASTRE (DRAG)
36
3. PROCEDIMIENTO GENERAL PARA SOLUCIN DE

PROBLEMAS CFD TRIDIMENSIONAL EN ANSYS
41
3.1 GENERACIN DE LA MALLA ESTRUCTURADA
42
3.2 CONDICIONES DE FRONTERA
42
3.3 CONDICIONES DE FRONTERA EN LA SUPERFICIE SLIDA (PILA)
42
3.4 CONDICIONES DE FRONTERA DE FLUJO DE

ENTRADA O FLUJO DE SALIDA
43
4. PUENTE
44
4.1 PILAS
45
4.2 SOCAVACIN LOCAL EN PILAS
47
5. PLANTEAMIENTO DEL PROYECTO
50
5.1 DOMINIO DEL PROBLEMA
56
5.2 DISCRETIZACIN
57
6. METODOLOGA PARA LA APLICACIN DEL SOFTWARE

(ANSYS) EN LA MODELACIN
58
6.1 PILA RECTANGULAR MACIZA
58
6.1.1 Construccin de la geometra de la pila
61
6.1.2 Discretizacin
67
6.1.3 Condiciones de frontera del flujo de entrada
69
6.1.4 Condicin de frontera del flujo de salida
71
6.1.5 Condicin de frontera en la superficie solida (pila)
71
6.1.6 Solucin
73
6.1.7 Postproceso
75
6.2 PILA CIRCULAR APORTICADA
81
83
6.2.2 Discretizacin
91
93
6.2.4 Condicin de frontera del flujo de salida
94
94
6.2.6 Solucin
95
6.2.7 Postproceso
96
6.3 PILA EN PUNTA DE DIAMANTE
101
103
6.3.2 Discretizacin
109
111
6.3.4 Condicin de frontera de flujo de salida
112
112
6.3.6 Solucin
114
6.3.7 Postproceso
115
7. ANALISIS DE RESULTADOS
119
8. GEOMETRIA PROPUESTA PARA OPTIMIZACIN HIDRULICA DE

PILARES DE PUENTE
126
9. CONCLUSIONES
133
10. RECOMENDACIONES
134
BIBLIOGRAFA
135
ANEXOS
138
LISTA DE FIGURAS
Pg.
Figura 1. Sistema de vrtices en socavacin local en pilas
27
Figura 2. Nmero de Froud en el entorno del puente
27
Figura 3. Modificacin de la geometra del puente
28
Figura 4. Geometra B probada por Roca
29
Figura 5. Campo vectorial para flujo tridimensional
30
Figura 6. Flujo a superficie libre
32
Figura 7. Flujo turbulento
33
Figura 8. Perfil de velocidad en pila cuadrada proveniente

de la simulacin numrica
34
Figura 9. Lneas de corriente a travs de un cuerpo hidrodinmico
35
Figura 10. Diferencia entre el flujo rotacional y el irrotacional
36
Figura 11. Pila circular en lmina libre que ilustre el punto

de estancamiento, separacin y estela turbulenta creada detrs de la pila
37
Figura 12. rea proyectada de la pila cuadrada necesaria

para el clculo del arrastre de presin
40
Figura 13. Dominio computacional
41
Figura 14. Fronteras empleadas en el modelo
43
Figura 15. Seccin longitudinal puente
45
Figura 16. Caracterstica del flujo alrededor de una pila
48
Figura 17. Hidrodinmica y morfologa de un foso de erosin local en planta
49
Figura 18. Geometra de pilas con igual base
50
Figura 19. Seccin transversal del canal natural para anlisis del proyecto
54
Figura 20. Geometra del domino del problema
56
Figura 21. Geometra de FLUID142
57
Figura 22. Herramienta que ofrece ANSYS para el anlisis

en el campo que se necesita
59
Figura 23. Tipo de elemento
60
Figura 24. Propiedades del material, densidad
60
Figura 25. Propiedades del material, viscosidad
61
Figura 26. Tipo de anlisis pila rectangular
62
Figura 27. Vista de los ejes de coordenadas y orientacin

del sistema de coordenadas globales
63
Figura 28. Orientacin de las ventanas de salida, con el

men interactivo de ANSYS
63
Figura 29. Volumen para el rea rectangular
66
Figura 30. Forma del volumen de control por extruccin pila rectangular
67
Figura 31. Equilibrio discretizado en la pila rectangular
69
Figura 32. Giros en cualquier direccin
70
Figura 33. Aplicacin de carga en el rea especificada en la pila rectangular
72
Figura 34. Propiedades del fluido pila rectangular
73
Figura 35. Propiedades del flujo pila rectangular
74
Figura 36. Solucin nodal de la velocidad del fluido en la pila rectangular
75
Figura 37. Wp settings pila rectangular
76
Figura 38. Movimiento de working plane en el lecho del

rio en la pila rectangular
77
Figura 39. Cambio a working plane para mostrar los resultados
77
Figura 40. Campo de velocidad en el lecho del rio en la pila rectangular
78
Figura 41. Vector del gradiente de velocidad en la pila rectangular
78
Figura 42. Puntos de inters donde se pretende simular

el movimiento de las lneas de corriente
79
Figura 43. Plot flow trace pila rectangular
80
Figura 44. Trayectoria de las lneas de corriente en la pila rectangular
80
Figura 45. Perfil de presiones sobre el lecho del rio en la pila rectangular
81
Figura 46. Generacin de puntos para crear rea circular aporticada
85
Figura 47. Generacin de lneas a partir de los puntos
86
Figura 48. Lneas que componen la seccin circular
87
Figura 49. Pila circular aporticada
89
Figura 50. Pila circular aporticada interactuando con el equilibrio
89
Figura 51. Forma del volumen de control por extruccin
90
Figura 52. Equilibrio discretizado en la pila circular aporticada
91
Figura 53. Aplicacin de carga en el rea especificada en la pila circular

aporticada
95
Figura 54. Solucin nodal de la velocidad del fluido en la

pila circular aporticada
96
Figura 55. Campo de velocidad en el lecho del rio en la pila circular aporticada 97
Figura 56. Vector del gradiente de velocidad en la pila circular aporticada
98

99
Figura 58. Plot flow trace pila circular
99
Figura 59. Trayectoria de las lneas de corriente en la pila circular aporticada 100
Figura 60. Perfil de presiones sobre el lecho del rio en

la pila circular aporticada
101
Figura 61. Generacin de puntos para crear el rea de punta diamante
105
106
Figura 63. Pila en punta de diamante
107
Figura 64. Equilibrio interactuando con la pila
108
109
Figura 66. Equilibrio discretizado en la pila punta diamante
111
Figura 67. Aplicacin de carga en el rea especificada

en la pila punta diamante
113
Figura 68. Solucin nodal de la velocidad del fluido en la pila punta diamante 114
Figura 69. Campo de velocidad en el lecho del rio en la pila punta diamante
116
Figura 70. Vector del gradiente de velocidad en la pila punta diamante
116

117
Figura 72. Trayectoria de las lneas de corriente en la pila punta diamante
118
Figura 73. Perfil de presiones sobre el lecho del rio en la pila punta diamante 118
Figura 74. Gradientes de velocidad prxima a la frontera
solida para las formas geomtricas analizadas paralelas al flujo
121
Figura 75. Perfil de presin para las formas geomtricas

analizadas sobre el sobre el lecho del rio
123
Figura 76. Fenmenos hidrulicos
124
Figura 77. Propuesta de las nuevas geometras para optimizacin hidrulica
127
Figura 78. Isocontorno de presin en geometras propuestas
130
Figura 79. Esquema general P2
131
Figura 80. Campo de velocidad y lneas de corriente en P2
131
LISTA DE TABLAS
Pg.
Tabla 1. Resumen de los datos de campo (45 datos) empleados

en el proyecto
51
Tabla 2. Distancia a nivel de banca a banca del Perfil del cauce
53
Tabla 3. Calculo de elemento geomtrico del canal natural
55
Tabla 4. Calculo de perfil logartmico
56
Tabla 5. Coordenadas para dar forma al rea rectangular
64
Tabla 6. Generacin de lneas a partir de los puntos P.R
65
Tabla 7. Generacin de reas entre lneas P.R
65
Tabla 8. Tasa de divisin de lneas que componen el equilibrio P.R
68
Tabla 9. Coordenadas para dar forma al rea de pila circular aporticada
84
Tabla 10. Generacin de lneas a partir de los puntos P.C
86
Tabla 11. Generacin de reas entre lneas P.C
88
Tabla 12. Tasa de divisin de lneas que componen el equilibrio P. C
92
Tabla 13. Coordenadas para dar forma al rea de pila punta diamante
104
Tabla 14. Generacin de lneas a partir de los puntos P.P. D
105
Tabla 15. Generacin de reas entre lneas P.P. D
106
Tabla 16. Tasa de divisin de lneas que componen el equilibrio P.P. D
110
Tabla 17. Clculo de arrastre de forma
125
Tabla 18. Componentes de velocidad aguas arriba en lecho de rio para

nuevas geometras pro las puestas
128
Tabla 19. Estimacin de la socavacin local
132
LISTA DE GRAFICAS
Pg.
Grafica 1. Gradiente de velocidad aguas arriba normal al flujo
119
Grafica 2. Gradiente de velocidad aguas abajo normal al flujo
120
Grafica 3. Variaciones de las velocidades aguas arriba de las pilas
129
LISTA DE FOTOS
Pg.
Foto 1. Alma llena y voladizos
46
Foto 2. Alma llena
46
Foto 3. Pila aporticada
47
LISTA DE ANEXOS
Pg.
Anexo A. Lista de comandos para modelar la seccin p2 utilizando

seccin p2 utilizando ANSYS 10
139
Anexo B. Sitios de inters para el aprendizaje de ANSYS 10
142
GLOSARIO
CALLE DE VRTICE: ordenamiento de vrtices dentro de la estela.

CARENAR: optimizar aerodinmicamente o hidrodinmicamente un cuerpo para
disminuir la friccin con el fluido en el cual se encuentra.
DISCRETIZACIN: es el proceso de transformar el dominio geomtrico continuo
en una malla de elementos finitos.
EQUILIBRIO: regin fija en el espacio donde puede existir flujo de fluido a travs
de sus fronteras.
FLUJO INCOMPRESIBLE: el flujo es incompresible si los cambios de densidad
de un punto a otro son despreciables.
FLUJO PERMANENTE: los parmetros hidrulicos del flujo permanecen
constantes en el tiempo o sea que la velocidad de las partculas que ocupan un
punto dado es la misma para cada instante.
FLUJO TRIDIMENSIONAL: el gradiente del flujo existen en las tres direcciones,
este es el caso ms general del flujo y es el que presenta la mayora de
conducciones o cauces.
FUERZAS DE ARRASTRE: fuerza sobre un cuerpo ocasionada por el fluido que
opone resistencia en la direccin del movimiento del cuerpo, se descompone en
arrastre de forma y arrastre de friccin.
GRADOS DE LIBERTAD: numero de variables principales del problema,
determina el sistema completamente.
MTODO DE ELEMENTO FINITO: se aproxima el dominio geomtrico continuo
por una unin de subdominios ms simples, y a cada uno de estos se le
denomina elemento finito.
NODO: punto donde se va a encontrar la solucin del problema para despus
extenderla al resto del dominio interpolando dicha solucin. Pueden estar en los
vrtices o dentro del elemento.
PUNTO DE ESTANCAMIENTO: se define como un punto en el campo de flujo en
donde la velocidad es idnticamente cero.
SECCIN DE CONTROL: Paso de rgimen lento a rpido.

VELOCIDAD: vector que representa la direccin, sentido y magnitud de la rapidez
de movimiento del fluido.
INTRODUCCIN
En el momento que se construye una pila de una obra de arte vial, como es el
caso de un puente, se produce perturbaciones en el funcionamiento hidrulico del
cauce en el entorno donde se encuentra la pila, debido a que el elemento
estructural se localiza inmerso en el flujo, generando un gran obstculo al cuerpo
de agua, esto provoca que el flujo aguas arriba de dicho puente este desacoplado
del flujo aguas abajo, creando sobre elevaciones al frente de la pila.
En la actualidad las soluciones a los problemas de optimizacin como el descrito

anteriormente son modelos matemticos complejos debido a la variacin
constante del campo de velocidad. Dichos modelos se han podido resolver unidos
al desarrollo de la modelacin numrica en reas de conocimiento prximos que
permiten obtener soluciones aproximadas de flujos, que sean susceptibles de ser
representados por un sistema de ecuaciones diferenciales.
La existencia actualmente de potentes modelos computacionales basados en el

mtodo de elementos finitos (M.E.F) permite producir comportamientos de flujos
de agua en cualquier tipo de contorno. Una de las principales ventajas de utilizar
la aproximacin por elementos finitos, debido a su gran capacidad para ajustar
geometras complejas y permitir un refinamiento local en los lugares del modelo en
que se hace necesario. Sin embargo, la mejora en la solucin es tambin muy
importante, ya que en general, en los problemas prcticos se precisa de una
discretizacin en tres dimensiones, con un nmero de grados de libertad mucho
mayor que en problemas tpicos de estructuras (son usuales en dichos casos de
105 a 106 grados de libertad).
El trabajo presentado en el siguiente proyecto desarrollar esquemas numricos

tridimensionales, que representar el comportamiento del flujo en cada pila,
efectuando el anlisis hidrulico del flujo alrededor de las diferentes geometras de
pilas (rectangular, circular y punta de diamante), y a su vez comparar el impacto
que crean las secciones de pilas en el cauce, conjuntamente se proporcionara
un perfil hidrodinmico que perturbe en menor grado el flujo a travs de, software
comercial llamado ANSYS.
25
1. ANTECEDENTES
En el presente capitulo se hace un repaso de las tcnicas utilizadas por otros

investigadores, para resolver problemas similares al objetivo del trabajo de grado
(modelacin del flujo de agua en pilares de puentes con superficie libre), y se
comentan tambin algunos estudios ms tericos que se aplicarn en la
simulacin numrica.
El profesor Irving H. Shames posee el ttulo de profesor de facultad de ingeniera

y ciencia aplicada de la Universidad del estado de Nueva York en Buffalo, escribi
en su libro de mecnica de fluidos que el arrastre de presin en un cuerpo inmerso
en el agua puede reducirse agregndole material en el lado aguas abajo del
cuerpo en un esfuerzo para retrasar la separacin tanto como sea posible. Esto es
dar una forma aerodinmica o hidrodinmica, este paso disminuir el arrastre de
presin.
El seor R. H. Lorefice en el ao 2006 publico en su articulo de Erosin local de

pilas a nivel mesomecnico por simulacin numrica, que el desprendimiento del
material de fondo (suelo o sedimento) se debe principalmente a la accin de
fuerzas tractivas, las cuales bajo determinadas condiciones del flujo, generan un
estado combinado de tensiones que eventualmente supera a la resistencia
cohesiva/friccional del suelo. Desde el punto de vista de la fsica que gobierna el
problema, el desprendimiento del material del fondo ocurre debido a que las
fuerzas tractivas originadas por la dinmica del flujo vencen la resistencia
intrnseca del material en estado saturado.
Juan Pedro Martn Vide en el ao 2007 manifiesta en su libro Ingeniera de ros

que fsicamente el fenmeno de erosin local consiste que alrededor de la pila se
dan velocidades locales mayores que los de la corriente media, acompaadas de
un sistema de vrtices frontales, laterales y estelas detrs de la pila (Vase la
Figura 1), sealando al sistema de vrtices mencionado anteriormente como el
factor principal la socavacin local.
26
Figura 1. Sistema de vrtices en socavacin local en pilas
Fuente: MARTN VIDE, Juan Pedro. Ingeniera de ros. Barcelona: UPC, 2006.
El Ph.D Ernest Blad Castellet en el ao 2005 estudio la influencia de las pilas

de puentes antiguos en el funcionamiento hidrulico del rio Llierca, dicho estudio
muestra que el puente antiguo representa una seccin de control (paso de
rgimen lento a rpido) basado en el numero de Froude, en las zonas azules el
rgimen de flujo es subcrtico y en las zonas negra el rgimen de flujo es
supercrtico (Vase la Figura 2).
Figura 2. Nmero de Froud en el entorno del puente
Fuente: CASTELLET, Ernest Blad. Modelacin del flujo en lmina libre sobre
cauces naturales. Barcelona: Universidad Politcnica de Catalunya, 2005.
27
Adems del perfil de velocidad en el rio aguas arriba y aguas abajo de la pila,
analiz los resultados entregados por la simulacin numrica y propuso carenar
(adecuar cualquier elemento con forma hidrodinmica o aerodinmica) la pila
aguas arriba con un semicrculo y aguas abajo con un perfil Joukovsky. (Vase la
Figura 3).
Figura 3. Modificacin de la geometra del puente
Fuente: CASTELLET, Ernest Blad. Modelacin del flujo en lmina libre sobre
cauces naturales. Barcelona: Universidad Politcnica de Catalunya, 2005.
Roca prob tres nuevas geometras para el elemento protector contra la

socavacin local denominadas A, B y C. La forma de estas geometras fue
escogida de acuerdo a la forma del cuenco de socavacin que se produce en la
estructura sin proteccin, (Vase la figura 4), se muestra la geometra con el
porcentaje de reduccin de socavacin local de 89[%], el mayor presentado en su
investigacin.
28
Figura 4. Geometra B probada por Roca
Fuente: ROCA, Alex. Socavacin en Pilas de Puente. Bogot: Universidad de los

Andes, 1997.
29
2. FUNDAMENTOS HIDRAULICOS
En estudios hidrulicos, se requiere conocer con detalle el campo de velocidades

de un flujo a superficie libre. Ejemplos de estos estudios son: el efecto en las
caractersticas del flujo al colocar pilas de puente o diques en un rio, La
estimacin de la erosin y el depsito de sedimentos en cauces.
En la actualidad es comn el uso de modelos matemticos en la solucin de

problemas de flujo a superficies libre. Dependiendo del problema que se va a
resolver, se puedan usar modelos unidimensionales, bidimensionales y
tridimensionales.
2.1 FLUJO TRIDIMENSIONAL
Es el flujo ms general en el que las componentes de velocidad u, v, w en

direcciones mutuamente perpendiculares son funciones de las coordenadas del
espacio y del tiempo x, y, z y t (Vase la Figura 5). Los mtodos de anlisis son
generalmente complejos matemticamente. Es frecuente, usar un modelo
tridimensional cuando se necesita conocer la variacin del campo de velocidades
en las tres componentes de velocidad.
Figura 5. Campo vectorial para flujo tridimensional
30
2.2 PROPIEDADES DEL FLUJO
2.2.1 Velocidad. Es un vector que representa la direccin, sentido y magnitud de

la rapidez de movimiento del fluido. Componente principal en el estudio del
comportamiento del flujo en el modelo computacional.
2.3 PROPIEDADES DEL FLUIDO
Las siguientes son algunas de las propiedades de los fluidos. Los valores de stas
pueden depender de otras variables como: temperatura, presin, etc.
2.3.1 Densidad. La densidad () de un fluido es su masa por unidad de volumen.

Propiedad del material (agua) que se debe especificar como dato de entrada para
el modelo.
2.3.2 Viscosidad dinmica. En la prctica se observa que algunos fluidos se

mueven con mayor facilidad que otros. Esto se debe a fuerzas de rozamiento
internas en el fluido. Propiedad del material (agua) que se debe especificar como
dato de entrada para el modelo.
2.4 TIPOS DE FLUJO
El flujo es el movimiento de un fluido y se clasifican segn el movimiento, el

espacio y el tiempo.
2.4.1 Flujo libre. El movimiento del fluido se realiza por conductos abiertos o
cerrados parcialmente llenos, de forma que existe una superficie libre que est en
contacto con la atmsfera(Vase la Figura 6). El movimiento se realiza gracias a la
fuerza de la gravedad.
31
Figura 6. Flujo a superficie libre

Superficie libre del agua
Fondo del canal
Nivel de referencia
El flujo libre tiene lugar en la naturaleza en ros, arroyos y en forma artificial en

canales de conduccin de fluidos, acueductos, alcantarillados, drenajes, etc.
2.4.2 Flujo variado permanente. Los parmetros hidrulicos del flujo varan en el
espacio pero no en el tiempo. Este tipo de flujo puede subdividirse en
gradualmente variado o rpidamente variado.
ANSYS dentro de sus aplicaciones permite elegir el tipo de rgimen para el

modelo, en este caso ser steady state flow analysis, donde el modelo del
proyecto acepta la variacin de la distribucin de velocidad y presin en el
equilibrio mas no en el tiempo.
2.5 RGIMEN DEL FLUJO
2.5.1 Flujo con rgimen turbulento. Se presenta si las fuerzas viscosas son
dbiles con relacin a las fuerzas inerciales. Las partculas del fluido con rgimen
laminar se mueven ordenadamente siguiendo trayectorias definidas, pero al
aumentar la velocidad las partculas del fluido chocan entre s y se desvan
siguiendo trayectorias irregulares que no son suaves ni fijas (Vase la Figura 7), y
que constituyen el flujo turbulento para un Re> 1000[-].
Se hace la anotacin que el rgimen caracterstico de flujo en ros es turbulento,

adems cuando el fluido interacta con obstculos slidos como pilas de puente,
se presentan turbulencias en el entorno de la pila. Por esta razn se especifica el
tipo de anlisis turbulento para el modelo.
32
Figura 7. Flujo turbulento
Fuente: POTTER, Merle. Mecnica de fluidos. Bogot: Cengage Learning, 2002.
2.5.2 Flujo incompresible. El flujo es incompresible si los cambios de densidades

de un punto a otro son despreciables. Los lquidos y gases a bajas temperaturas
se consideran incompresibles. Para el modelo se plantea la incomprensibilidad,
sin variacin de densidad en el equilibrio.
2.6 FUNDAMENTOS DE VISUALIZACIN DEL FLUJO
La visualizacin del flujo permite observar el movimiento del fluido para el modelo
que se desarrolla en el paquete computacional, debido a que la mente humana
esta diseada para procesar con rapidez una cantidad increble de informacin
visual; como se dice: una imagen vale ms que mil palabras. De hecho,
precisamente lo primero que un ingeniero hace cuando utiliza CFD, despus de
obtener una solucin numrica , es simular alguna forma de visualizacin del flujo,
de modo que pueda ver la imagen completa , en vez de solo una lista de nmeros
y datos cuantitativo.
2.6.1 Punto de Estancamiento. Se define como un punto en el campo de flujo en

donde la velocidad es idnticamente cero.
2.6.2 Equilibrio. Se considera, un volumen aislado de fluido en cierto instante de

tiempo. Se entiende por equilibrio de control una regin fija en el espacio donde
puede existir flujo de fluido a travs de sus fronteras. Se estableci la geometra
del domino del problema como un equilibrio.
33
Figura 8. Perfil de velocidad en pila cuadrada proveniente de la simulacin

numrica
2.6.3 Lneas de corriente. Son lneas cuya tangente es en todos los puntos
paralela al vector velocidad en un instante. Las lneas de corriente permiten
simular el movimiento de las partculas del fluido en determinado instante.
(Vase la Figura 9).
2.6.4 Graficas de contornos. Muestra las curvas de valor constante de una

propiedad escalar (o magnitud de una propiedad de un vector) en un determinado
instante. Las graficas de contorno exteriorizan la variacin de la presin, la
magnitud de la velocidad en el modelo computacional. Una grafica de contornos
puede relevar con rapidez las regiones de valores altos (o bajos) de la propiedad
del flujo que se est estudiando.
Una grafica de contorno puede consistir, sencillamente, de curvas que indiquen

varios niveles de la propiedad; esta se conoce como grafica de lneas de
contornos. De modo opcional, los contornos se pueden rellenar con colores o
sombras de gris; esto se conoce como grafica de contornos rellenos. En la CDF a
menudo se presenta las lneas de contorno en llamativos, indicando por lo comn
con rojo el valor ms alto del escalar y con azul el ms bajo. El ojo humano sano
puede identificar con facilidad una regin roja o azul y de este modo, localizar las
regiones de valor alto o bajo de la propiedad del flujo. Debido alas bellas imgenes
producidas por la CDF, a la dinmica de fluidos a todo color.
34
Figura 9. Lneas de corriente a travs de un cuerpo hidrodinmico
Fuente: CENGEL, Yunus. CIMBALA, John. Mecnica de Fluidos. Bogot: Mc

Graw Hill Interamericana, 2006.
2.7 VORTICIDAD Y ROTACIONALIDAD
La vorticidad es una medida de la rotacin de una partcula de fluido. Si la

vorticidad en un punto en un campo de flujo es diferente de cero, la partcula de
fluido que llegue a ocupar ese punto en el espacio est girando; se dice que el
flujo en esa regin es rotacional. De modo semejante, si la vorticidad en una
regin del flujo es cero o (o despreciable pequea) las partculas de fluido all no
estn girando; se dice que el flujo en esa regin es irrotacional. Desde el punto de
vista fsico, las partculas de fluido que estn en una regin rotacional de un flujo
giran extremo sobre extremo a medida que avanzan en ese flujo.
Por ejemplo, partculas de fluido dentro de la capa limite viscosa cercana a una
pared solida son rotacionales (y, por lo tanto, tienen vorticidad diferente de cero),
en tanto que las partculas de fluido que estn afuera de la capa limite son
irrotacionales (y su vorticidad es cero).Estos dos casos se ilustran en la figura 10.
La rotacin de los elementos de fluidos se asocia con las estelas, las capas
limites.
Por consiguiente, si el flujo se origina en una regin irrotacional, contina siendo

irrotacional hasta que algn proceso no uniforme lo altera. Las lneas de corriente
en las condiciones de frontera de entrada, nunca se cortan porque el flujo es
35
permanente, como resultante el flujo que se origina es irrotacional y se mantiene

constante hasta que se presenta el rea frontal de la pila en su trayectoria,
haciendo rotar las partculas de fluido generando los vrtices laterales alrededor
de la pila.
Figura 10. Diferencia entre el flujo rotacional y el irrotacional

2.8 FUERZA DE ARRASTRE (DRAG)
Cuando un flujo se desplaza sobre un cuerpo solid, ejerce fuerzas de presin

normales a la superficie y fuerzas de corte paralelas a la superficie a lo largo de
la superficie exterior del cuerpo. Usualmente existe inters en las resultantes de
las fuerzas de presin y de corte que actan sobre el cuerpo en lugar de los
detalles de las distribuciones de estas fuerzas en la superficie del cuerpo, la
componente de la fuerza de presin y de corte resultante que acta en la
direccin del flujo se llama fuerza de arrastre (o slo de arrastre).En ocasiones
un fluido se mueve sobre un cuerpo estacionario (como el viento que sopla sobre
un edificio o el empuje dinmico del agua sobre las pilas de puente) y otras veces
un cuerpo se mueve a travs de un fluido quieto (como un automvil que se
desplaza por el aire). Estos procesos al parecer diferentes son mutuamente
equivalentes; lo que importa es el movimiento relativo entre el fluido y el cuerpo.
Estos movimientos se analizan de manera conveniente cuando se fija el sistema
de coordenadas sobre el cuerpo; se les conoce como flujos sobre cuerpos o Flujo
externo.
36
Los cuerpos sometidos a flujo de fluidos se clasifican como currentilneos o

romos, esto depende de su forma global. Se dice que un cuerpo es currentilneos
Si se hace un esfuerzo conciso para alinear su forma con las lneas de corriente
anticipadas en el flujo. De otro modo un cuerpo (como una pila de puente) tiende
a bloquear el flujo y se dice que es un romo.
La velocidad del fluido que se aproxima a un cuerpo se llama velocidad de flujo

libre y se denota con V. Para que se visualice la influencia de la presin dinmica
sobre el arrastre (Vase la figura 11), donde se muestra la esfera en una corriente
de fluido. Las lneas de corriente ilustran la trayectoria del fluido conforme ste se
aproxima y fluye alrededor de la esfera. En el punto s de la superficie de la esfera,
la corriente
de fluido est en reposo estancada. El termino punto de
estancamiento se emplea para denotar dicho punto la relacin entre la presin ps
y la presin de la corriente no perturbada del punto 1, se encuentra por medio de
la ecuacin de Bernoulli a lo largo de la lnea de corriente.
p
v
1 s
2g
p1
(1)
Figura 11. Pila circular en lmina libre que ilustre el punto de estancamiento,
separacin y estela turbulenta creada detrs de la pila
P1
Punto de separacin
Estela Turbulenta
Punto de estancamiento PS en la superficie de no deslizamiento
Al resolver para p s , obtenemos
37
p s p1
* v1 2
2g
Debido a que =
p s p1
* v1 2
2
(2)
, tenemos
g
(3)
La presin de estancamiento es mayor que la presin esttica en la corriente libre

por la magnitud de la presin dinmica. La energa cintica de la corriente libre en
movimiento se transforma en un a clase de energa potencial en forma de presin.
Es de esperar que el incremento de presin en el punto de estancamiento
produzca una fuera sobre el cuerpo, opuesta a su movimiento, es decir, una fuerza
de arrastre. Sin embargo, la magnitud de la fuerza depende no solo de la presin
de estancamiento, sino tambin de la presin en lado trasero del cuerpo. Debido a
la dificultad para predecir la variacin real de la presin en el lado posterior, es
comn que se utilice el coeficiente de arrastre.
El arrastre total sobre un cuerpo se debe a dos componentes, el arrastre de

presin (tambin llamado arrastre de forma) se debe a disturbios de la corriente
de flujo
conforme pasa el cuerpo, lo que crea una estela turbulenta. Las
caractersticas de los disturbios dependen de la forma del cuerpo, y a veces del
nmero de Reynolds de flujo y de la rugosidad de la superficie. El arrastre de
friccin se debe a las fuerzas cortantes en la capa delgada de fluido que se
encuentra cerca de la superficie del cuerpo, la cual recibe el nombre de la capa
lmite. Se describen en detalle estos dos tipos de arrastre.
El arrastre de presin es proporcional es proporcional al rea frontal y a la

diferencia de presiones que actan sobre la parte frontal y posterior del cuerpo
sumergido. En consecuencia el arrastre de presin usualmente domina para
cuerpos romos.
La ecuacin para calcular la fuerza de arrastre en funcin del coeficiente de

arrastre se encuentra definida por:
38
v2
FD CD
A
2
(4)
Los trminos en esta ecuacin son los siguientes:
CD es el coeficiente de arrastre. Se trata de un nmero adimensional que

depende de la forma del cuerpo, rugosidad y su orientacin con respecto a la
corriente del fluido, los respectivos coeficientes de arrastre para las formas
geomtricas de pilas (Circular, rectangular y punta de diamante) se extrajeron
de la tabla 11- 1 del libro: CENGEL, Yunus. CIMBALA John. Mecnica de
Fluidos. Bogot: Mc Graw Hill Interamericana, 2006.
es la densidad del fluido. Debido a que la densidad de los lquidos es mayor

que los gases, el orden general de magnitud de las fuerzas de arrastre sobre
objetos que se mueve en el agua o la presin ejercida por el agua a un cuerpo
estacionario es mucho ms grande que para los objetos que se mueven en el
aire.
es la velocidad superficial de la corriente libre del fluido en relacin con el

cuerpo. En general, no importa si el que se mueve es el cuerpo o el fluido. Sin
embargo, la ubicacin de otras superficies cerca del cuerpo de inters afecta
el arrastre.
A es el rea mxima proyectada en un plano normal al flujo, es el rea que

una persona vera si observase al cuerpo desde la direccin del fluido que se
aproxima. Piense cual seria la forma bidimensional mas grande si mirara desde
la corriente aguas arriba la pila (Vase la figura 12); donde el plano yz de color
azul representa la corriente, Ap. define el rea mxima proyectada Esa seria el
rea que se empleara para calcular el arrastre sobre la pila, y h constituye el
tirante del cauce aguas arriba de la pila, el arrastre calculado con la proyeccin
en el plano normal al flujo del rea se denomina arrastre de forma o arrastre de
presin.
39
Figura 12. rea proyectada de la pila cuadrada necesaria para el clculo del
arrastre de presin
Ap.
X
z
40
3. PROCEDIMIENTO GENERAL PARA SOLUCIN DE PROBLEMAS CFD

TRIDIMENSIONAL EN ANSYS
Para resolver de manera numrica los problemas ingenieriles, se efectan los

siguientes pasos.
Se elige un dominio computacional y se genera una malla (conocida tambin

como red de nodos) con elemento FLUID142; el dominio se divide en muchos
elementos pequeos llamados celdas. Para dominios tridimensionales (3-D) las
celdas son volmenes (Vase la Figura 13). Puede considerase a cada celda
como un pequeo volumen de control en el que se resuelven las versiones
separadas de las ecuaciones de conservacin y transporte. La calidad de una
solucin de DFC depende mucho de la calidad de la malla.
Las condiciones de frontera se especifican en cada cara del dominio (flujos 3-D).
Se especifica el tipo de fluido (agua), juntos con las propiedades fsico-qumicas
del fluido (temperatura, densidad, viscosidad, etc.). Se procede a la solucin, las
variables de campo de flujo como la velocidad y la presin se trazan y se analizan
de manera grafica.
Figura 13. Dominio computacional

41
3.1 GENERACIN DE LA MALLA ESTRUCTURADA
Se debe generar una malla que define las celdas en las que se calculan las
variables de flujo (velocidad y presin, entre otras) en todo el dominio
computacional. Para el modelo se empleara Una malla estructurada consta de
celdas volumtricas con seis caras (3-D). Aunque la forma rectangular de las
celdas podra estar distorsionada, cada celda se numera de acuerdo con los
ndices i, j, k que no necesariamente corresponden a las coordenadas x, y, z. En
la capas limites, donde las variables de flujo cambian con rapidez en la direccin
normal a la pared y se requieren mallas de altas resolucin en la cercana de sta,
por se debe elegir la malla estructurada que permiten una resolucin mucho ms
fina en la interaccin de agua- pila.
3.2 CONDICIONES DE FRONTERA
Las condiciones de frontera apropiadas se necesitan para obtener una solucin de

CFD exacta en el modelo. Se establece las condiciones de frontera, la cuales son;
velocidad y presin, el tipo de flujo que se modela, se determina mediante las
condiciones de frontera impuestas. Cada modelo presenta condiciones de
fronteras diferentes, en las descripciones dadas, se emplean las palabras cara o
plano, lo que significa flujo tridimensional. (Vase la figura 14), se muestran las
condiciones de frontera que se deben especificar para el modelo, en la frontera de
entrada se define la velocidad, en la frontera de salida se especifica la presin
manomtrica, la frontera solida en el modelo es la pila, para las cuales la velocidad
normal y tangencial son nulas.
3.3 CONDICIONES DE FRONTERA EN LA SUPERFICIE SLIDA (PILA)
La condicin de frontera ms simple es la que se formula para una superficie

solida (pila). Puesto que el fluido no pueda pasar por una superficie slida, la
componente normal de la velocidad relativa con respecto a la superficie slida se
iguala a cero a lo largo de aquella cara en la cual se establece la condicin de
frontera. Adems debido a la condicin de no deslizamiento, se fija tambin en
cero la componente de la velocidad tangencial a una superficie slida en reposo.
(Vase la Figura 14).
42
Las superficies slidas mviles y superficies con esfuerzos de corte especificados,

como en los casos donde se pretende simular el flujo en la interfaz pila-suelo-agua
donde se desea dejar que el fluido se deslice a lo largo de la pared (a esto se le
conoce como pared invscida).
3.4 CONDICIONES DE FRONTERA DE FLUJO DE ENTRADA O FLUJO DE

SALIDA
Existen varias opciones de definir las condiciones de frontera en aquellas partes

de frontera a travs de las cuales entra el fluido al dominio computacional (flujo de
entrada) o sale del dominio (flujo de salida). Se clasifican en general como
condiciones que especifican la velocidad (entrada o salida de presin).
Figura 14. Fronteras empleadas en el modelo

Frontera Solida (Pila)
Salida
Entrada
Dominio Computacional
En una entrada de velocidad se especifica la velocidad del flujo entrante a lo

largo de la cara de admisin.
En una salida de presin, el fluido sale del dominio computacional. Se especifica

la presin esttica a lo largo de la cara de salida; en muchos casos esta es la
presin atmosfrica (presin manomtrica cero). Las propiedades del flujo, como
la temperatura, y las propiedades de turbulencia se especifican tambin en las
entradas y salidas de presin.
43
La presin no se especifica en la entrada de velocidad, ya que esto dara lugar a

una sobre especificacin matemtica, porque la presin y las velocidades acoplan
en la ecuacin de movimiento. Adems, la presin en una entrada de velocidad se
ajusta por si misma para coincidir con el resto del campo de flujo. De modo similar,
la velocidad no se especifica en una entrada o salida de presin, ya que esto
ocasionara tambin sobre especificacin matemtica. Adems, si se establece
una condicin que especifica presin, la velocidad se ajusta por si misma para
coincidir con el resto del campo de flujo.
Una situacin comn en una aplicacin simple de la DFC es especificar velocidad

en una o mas entradas de flujo al dominio computacional, presin en una o mas
salidas de flujos, y establecer las condiciones de frontera en las paredes que
definen la configuracin geomtrica del resto del dominio computacional.
44
4. PUENTE
Es una obra de arte que permite dar continuidad a una va salvando un obstculo,
el cual puede ser ro, una depresin del terreno u otra va en cuyo caso tenemos
los pasos elevados. El puente est conformado por la superestructura o zona por
donde pasa el trfico y por la infraestructura o elementos que trasmiten la carga al
suelo. La geometra del puente se muestra en una seccin longitudinal (perfil),
donde se indica: la superestructura y la infraestructura (pila y estribo).
Figura 15. Seccin longitudinal puente (perfil)
Superestructura
Pila
Lecho del rio
4.1 PILAS
Son apoyos intermedios que son soportan la superestructura, en trminos

generales no tienen empuje de tierra pero cuando son de gran altura influye
mucho el viento y el agua. Pueden ser de varios tipos segn la seccin transversal
en las Fotografa (1, 2 y 3), se puede observar las formas geomtricas que se
modelaron en el proyecto.
45
Alma llena
Alma llena y voladizos
Aporticada
Foto 1. Alma llena y voladizos
Foto 2. Alma llena
46
Foto 3. Pila aporticada
4.2 SOCAVACIN LOCAL EN PILAS
Cuando se coloca una pila de puente en la corriente de un ro se produce un

cambio en las condiciones hidrulicas de sta, y, por lo tanto, en su capacidad
para producir arrastre slido. Si la capacidad de arrastre supera localmente el
aporte del gasto slido del ro, ocurrir en la pila una socavacin local.
El fenmeno de socavacin o erosin local en pilares es uno de los problemas con

los cuales hay que tomar contacto durante el diseo de obras de lechos de ros, es
el pronstico de erosiones localizadas aguas arriba y aguas abajo de la obra
hidrulica. La profundidad de la erosin localizada, est estrechamente vinculada
a las particularidades hidrulicas del flujo, es decir a la distribucin de los caudales
especficos. Adems de los procesos naturales de modificacin permanente de los
lechos de ros (transporte y arrastre de sedimentos).
La erosin causada por el flujo alrededor de obstculos, como pilares de puente,

se llama erosin local. Fsicamente el fenmeno consiste que alrededor del pilar
se dan velocidades localmente mayores que las medias de la corriente,
acompaadas de un sistema de vrtices frontales, laterales y de estela detrs del
pilar. Aguas arriba de la pila, en la base y los alrededores se presentan vrtices
de eje horizontal y frontal (horseshoe vortex), Aguas abajo de la pila se presentan
vrtices de eje vertical; estela (wake vortex) vistos usualmente como remolinos
aguas debajo de la pila. Por las caractersticas de las lneas de corriente, los
vrtices de eje horizontal son los responsables de la mayor parte de la socavacin
producida. La geometra de una pila influencia significativamente la profundidad de
47
socavacin puesto que esta refleja la intensidad de los vrtices de eje horizontal
en la base de la pila. (Vase la figura 16) representacin grfica del proceso.
Figura 16. Caracterstica del flujo alrededor de una pila
Fuente: U.S DEPARTMENT OF TRANSPORTATION. Evaluating scour at bridges.

Washington: National Highway Institute, 2002.
Este sistema de vrtices es el principal responsable de la socavacin. Los granos

del lecho son aspirados por los vrtices y el fondo parece hervir por el movimiento
de los granos. El foso que se forma rodea al pilar, con la mayor profundidad y
extensin situada en la cara frontal. (Vase la Figura 17) Como en otros
fenmenos de erosin, hay una dependencia mutua entre flujo y el foso de
socavacin, de manera que a largo plazo, si las condiciones hidrulicas son
permanentes, se alcanza un equilibrio en la forma y tamao de los fosos.
Para estimacin de socavacin local en la pila rectangular, punta de diamante,

circular y P2, se utilizara la ecuacin propuesta por Chandra Nalluri, Francisco
Plata , Juan G. Saldarriaga, Revisin de ecuaciones que predicen la socavacin
local alrededor de pilas de puente, articulo 7 revista 9, Universidad de los Andes,
Colombia. Presentan la mejor ecuacin de datos mezclado (Mezcla de datos de
campo y laboratorio) para el clculo de socavacin local, producto del anlisis de
48
18 mtodos obtenidos tanto de campo como de laboratorio que se reportan en la

literatura.
ds
y
0.1176
D
b
0.538
Fr
0.126
b

D
1.25
(5)
ds
D
y
b
Fr
=
=
=
=
Profundidad de socavacin medida desde el

nivel medio del lecho
Dimetro medio de los sedimentos del lecho
Profundidad del flujo
Base o dimetro de la pila
Numero de Froude del flujo
Figura 17. Hidrodinmica y morfologa de un foso de erosin local en planta
Fuente: MARTN VIDE, Juan Pedro. Ingeniera de ros. Barcelona: UPC, 2006.
49
5. PLANTEAMIENTO DEL PROYECTO
Para trabajar con dimensiones reales de pilas, en la simulacin, se midieron las

dimensiones de la pila n 2 en el puente Francisco de Paula Andrade, las
dimensiones en planta fueron las siguientes:
Base = 0.80 [m]
Altura = 3.0 [m]
Las dimensiones extradas se tomaron como base, para la generacin de cada

una de las reas base para la pila (Rectangular, circular y punta de diamante) en
la simulacin. Se debe aclarar que para llevar a cabo la simulacin y comparacin
del comportamiento del flujo en el entorno de las pilas; las pilas poseen
dimensiones equivalentes, se entiende por dimensiones: Base, altura vistas en
planta, el tirante (profundidad del flujo) de agua siempre ser constante. A modo
de ejemplo la base de la pila rectangular es el mismo dimetro de la pila circular y
punta de diamante (Vase la Figura 18).
Figura 18. Geometra de pilas igual base
50
Para el presente anlisis se utilizaron datos de campo (45 datos) medidos

principalmente en puentes ubicados en la India recopilados por Kafi y Alam, se
extrajeron estos datos de india porque el modelo pretende simular el flujo de agua
en diferentes pila en condiciones reales, sin puntualizar ningn caso de puente
construido, los parmetros fueron los siguientes; velocidad mxima, dimetro
medio de los sedimentos del lecho y profundidad media del flujo, (Vase tabla 1).
Tabla 1. Resumen de los datos de campo (45 datos) empleados en el

proyecto
ESTADISTICA
m
seg
[mm]
[cm]
Valor Mnimo
1.19
0.1
206
Valor Mximo
7.03
1836
Promedio
2.08
0.59
674.86
Desviacin Estndar
1.40
1.08
386.32
Fuente: KAFI, M. and ALAM. Modification of Local Scour Equations.Journal of the

Institution of Engineers, Vol 76, No 5, pp 25-28. India: 1995.
La pendiente longitudinal del terreno empleada en el modelo, fue consultada del
informe tcnico presentado por la facultad de ingeniera de la Universidad
Francisco de Paula Santander sobre la creciente histrica del rio pamplonita el 11
de noviembre de 1984, Estos datos se escogieron para trabajar el modelo con
datos reales, algunos de los valores presentados en el informe son:
Pendiente longitudinal del rio: 0.7 [%]
Nivel del agua alcanzado en la avenida: 3.5 [m]
Velocidad de la corriente media: 3.69 m

seg
51
Para el clculo de la velocidad en la corriente superficial se empleara la ley

universal de la distribucin de velocidades, de Prandtl Von Karman. Esta ley ha
sido verificada mediante varios experimentos. Los resultados indican una
similaridad notable entre la distribucin de velocidad observada y la calculada, por
consiguiente, ofrece una justificacin razonable para el uso de esta ley logartmica
para el clculo de la velocidad superficial de la corriente que se utilizara en la
condicin de frontera de entrada para la simulacin.
La ecuacin de distribucin de velocidades en flujos turbulentos sobre superficies

rugosas (lecho de rio) es:
ds
y
0.1176
D
b
0.538
Fr
0.126
b

D
1.25
(6)
Vf = Velocidad de corte
y = Profundidad del flujo
K = Dimetro medio de los sedimentos
La velocidad de corte mencionada anteriormente se calcula mediante la siguiente

ecuacin:
Vf g s R
(7)
R = Radio hidrulico
s = Pendiente longitudinal del cauce
g = Aceleracin gravitacional
Para calcular el radio hidrulico de la seccin a estudiar se necesita el perfil

transversal del rio, en ese caso se desarrollo el perfil transversal propuesto por
A partir de los datos a continuacin genere la seccin transversal de un canal
natural.
52
Tabla 2. Distancia a nivel de banca a banca del perfil del cauce
Distancia desde un punto de referencia

cerca de la banca izquierda, en metros
BANCA
IZQUIERDA
-7
7,2
-5
5,6
-4
4,6
-2
4
0
1,9
1
1,5
2
1
3
0,8
5
0,7
BANCA
DERECHA
7
0,6
9
0,2
11
0
13
0,4
15
1,2
17
2,6
19
3,2
20
4,1
22
6,8
Fuente: CHOW , Ven Te. Hidrulica de canales abiertos. pp 34- 35, numeral [2-5],
Bogot: McGraw Hill Interamericana, 2004.
En la figura 19 se observa el perfil del cauce, la altura de la lamina de agua, que

es 1.7 veces la alcanzada por la avenida del rio pamplonita el 11 de noviembre de
1984, para hallar el rea y el permetro mojado; se grafico el perfil, luego por el
mtodo grafico se calculo el rea del rio. En la tabla 2, se presenta los clculos del
permetro mojado; dividiendo todo su permetro en tringulos rectngulos, pues
conociendo las cotas y las distancias aplicando Pitgoras se calcula la hipotenusa
de cada triangulo, que a posteriori sumndolas sern el permetro del cauce.
53
Figura 19. Seccin transversal del canal natural para anlisis del proyecto
8
SECCION TRANVERSAL DE CANAL NATURAL PARA ANALISIS DEL PROYECTO
NIVEL [m]
TIRANTE DE AGUA EN AVENIDA DE RIO

Ti = 6 [m]
3
SLOPE = 0,70%
0
10
DISTACIA [m]
54
11
12
13
14
15
16
17
18
Tabla 3. Calculo de elemento geomtrico del canal natural
DISTANCIA
-5
-4
2
-2 m
0
1
2
3
5
7
9
11
13
15
17
19
20
NIVEL
AH
AL
PERIMETRO
[m]
[m]
[m]
[m]
5,6
4,6
4
1,9
1,5
1
0,8
0,7
0,6
0,2
0
0,4
1,2
2,6
3,2
4,1
1
0,6
2,1
0,4
0,5
0,2
0,1
0,1
0,4
0,2
-0,4
-0,8
-1,4
-0,6
-0,9
-1
-2
2
1
1
1
2
2
2
2
2
2
2
2
1
1,41
2,09
2,90
1,08
1,12
1,02
2,00
2,00
2,04
2,01
2,04
2,15
2,44
2,09
1,35
PERIMETRO CANAL NATURAL

AREA CANAL NATURAL
RADIO HIDRAULICO CANAL NATURAL
27,74
156,72
5,65
[m]
[m]
El siguiente proceso, es hallar las velocidades que componen la distribucin

logartmica, para tomar el mayor valor de la distribucin, que a posteriori se
introducir al modelo computacional como condiciones de frontera de entrada
(perfil uniforme), la distribucin de velocidad desarrollada para el cauce (rio),
entrego rango de valores que fluctan entre [5 8,9 m/seg], tabla 3 se verifico que
la velocidad mxima calculada para el proyecto, se encuentra en el rango de
valores de datos de campo en puentes ubicados en la India recopilados por Kafi y
Alam (1995).
55
Tabla 4. Calculo de perfil logartmico
y
[ft]
1,64
3,28
4,92
6,56
8,2
9,84
11,48
13,12
14,76
16,4
18,04
19,68
Vf
ft
seg
2,044
2,044
2,044
2,044
2,044
2,044
2,044
2,044
2,044
2,044
2,044
2,044
[ft ]
ft
seg
1,935
1,935
1,935
1,935
1,935
1,935
1,935
1,935
1,935
1,935
1,935
1,935
16,52
20,05
22,12
23,59
24,73
25,66
26,45
27,13
27,73
28,27
28,76
29,20
m
seg

5,03
6,11
6,74
7,19
7,54
7,82
8,06
8,27
8,45
8,62
8,76
8,90
5.1 DOMINIO DEL PROBLEMA
La figura 20 esquematiza la geometra del domino del problema, la frontera

ABCDEFGH particulariza el domino del flujo de agua adecuado para la correcta
modelacin del problema.
Figura 20. Geometra del domino del problema
56
Las caras ACEG, ABEF y BDFH se localizan dos metros de la cara de la pila, se
suponen que a esa distancia el flujo no es perturbado, la altura del dominio es de
6 [m] y la cara CDHG se encuentra distanciada a seis metros de la parte posterior
de la pila donde el flujo se ha estabilizado.
Condicin de contorno
Sobre la cara ABEF = velocidad del Flujo 9 [m/seg] (se especifica un perfil
uniforme de velocidad, velocidad horizontal, equivalente a la velocidad del flujo
de agua, debido a que la coordenada horizontal es predominante en el
movimiento del agua en cauces)
Sobre la cara ACEG y la cara BDFH = velocidad del Flujo 9 [m/seg]
Sobre la frontera solida (pila) V= 0 [m/seg]; (se declara velocidad normal y

tangencial nula).
Sobre la cara CDGH: EFGH: ABCD = 0 (se define presin cero)
5.2 DISCRETIZACIN
Para el enmallado del equilibrio se utiliz como tipo de elemento base a

FLUID142 (Vase la Figura 22) que posee las siguientes caractersticas (las que
utilizan este modelo).
Figura 21. Geometra de FLUID142
Nodos: I, J, K, L, M, N, O, P
Grados de libertad: VX, VY, VZ, ENKE, ENDS
Propiedades del material: Densidad, viscosidad

dinmica.
Cargas: Fuerza (presin, velocidad)

57
6. METODOLOGA PARA LA APLICACIN DEL SOFTWARE (ANSYS) EN LA

MODELACIN
6.1 PILA RECTANGULAR MACIZA
Empleando la siguiente lista de comandos se proporcionar nombre al proyecto,

se indicara el titulo del problema y se establece al sistema internacional (SI) como
patrn de medidas en el proyecto.
/FILNAME, Pila rectangular
Define el nombre del proyecto como

pila Rectangular en la base de datos
Modelacin tridimensional del Indica el titulo del problema

flujo de agua en la pila
/TITLE,
/UNIT,SI
Se
establece el sistema internacional

de unidades
A continuacin se introducen parmetros escalares

que se utilizaran
posteriormente en la simulacin tridimensional del flujo de agua en pilas como lo
son; seccin de la pila y tirante aguas arriba.
Figura en planta de la pila especificando h y b, tirante del agua y control volumen
Ba= 0.8
Se define como base de la seccin de

la pila 1. [m]
He= 3
Altura de la seccin de la pila 3.0 [m]
Ti= 6
Profundidad
aguas arriba
58
del flujo en lamina libre
Seleccionamos del men principal de ANSYS (ANSYS main men) Preferences >
FLOTRAN CFD > OK; Se utiliza para establecer preferencias generales por
dinmica de fluidos computacionales (Computational fluid dynamic), herramienta
que ofrece (ANSYS products) para el anlisis en el campo del flujo de fluidos en
dos dimensiones 2D y tres dimensiones 3D.
Figura 22. Herramienta que ofrece ANSYS para el anlisis en el campo que
se necesita
Se ingresa al preprocesador, se especifica el tipo de elemento para enmallar el

volumen de control posteriormente y se puntualiza las propiedades del material
como son; densidad, viscosidad dinmica en funcin de la temperatura estndar
de 20 C.
59
/ PREP7
Se entra en el preprocesador para la

construccin del modelo en elementos
finitos
ET,1, FLUID 142
Tipo de elemento que se emplea

FLUID 142
Figura 23. Tipo de elemento
MP,DENS,1,998.92
Propiedad del material, densidad, del

elemento fluid 142, para 20 C es
Kg m
de 998,2
3
Figura 24. Propiedades del material, densidad
60
MP,VISC,1,1.10E-03
Propiedad del material, viscosidad, del

N seg
de 1.005E-03
2
Figura 25. Propiedades del material, viscosidad
Se define el tipo de anlisis como steady-state, adiabtico, turbulento e

incomprensible, en el men principal de ANSYS (ANSYS main men)
Preprocessor > FLOTRAN Set Up > Solution Options, cambiando el rgimen
de flujo laminar a turbulento como se observa en el cuadro de lista sombreado
con azul cuadro de dialogo.
6.1.1 Construccin de la geometra de la pila. Previamente al desarrollo de la

geometra de la pila rectangular, se cambiara la vista de los ejes de coordenadas,
la orientacin del sistema de coordenadas globales, para visualizar en forma ms
sencilla el step by step en la generacin de la pila, Se utilizarn los siguientes
comandos:
61
Figura 26. Tipo de anlisis pila rectangular
/VIEW,1,1,1,1
Define la direccin para la ventana de

salida (output) vase la figura 27
/VUP,1,Y
Orientacin del sistema de coordenadas

para el proyecto vase la figura 27
/REPLOT
Dibuja automticamente los cambios

hechos en la ventana de salida
62
Figura 27. Vista de los ejes de coordenadas y orientacin del sistema de

coordenadas globales
Los comandos anteriormente explicados cambian la orientacin de la ventana de

salida, pero ahora se explicara la misma accin trabajando con el men
interactivo de ANSYS. En la ventana de grafica damos click en Isometric view, as
realizaremos la misma tarea de forma simple (Vase la Figura 28).
Figura 28. Orientacin de las ventanas de salida, con el men interactivo de

ANSYS
63
La ubicacin de cada keypoint que dar forma al rea de la pila que

posteriormente cambiara a volumen por extruccin se proporciona en la tabla 5,
aclarando que todas las dimensiones se encuentran en metros.
Genera el keypoint, enumera, se
especifica las coordenadas x, y, z
respectivamente
K, nm., x, y, z
Tabla 5. Coordenadas para dar forma al rea rectangular
keypoint
1
2
3
4
5
6
7
8
x
0
0
2+He+5
2+He+5
2
2
2+He
2+He
y
0
0
-0.07
-0.07
-0.035
-0.035
-0.035
-0.035
z
0
2+Ba+2
0
2+Ba+2
2
2+Ba
2
2+Ba
Si escribi mal las coordenadas de cualquier keypoint puede eliminar dicho

keypoint con el comando KDELE, NP; Donde NP significa nmero del punto a
borrar.
KSEL,ALL
Selecciona todos los keypoints
LSEL,ALL
Selecciona todas las lneas
KPLOT,ALL
Dibuja todos los keypoints
LPLOT,ALL
Dibuja todas las lneas
/Pnum,KP,1
Numera todos los keypoints
/Pnum,LINE,1
Numera todas las lneas

64
Numera todas las areas
/Pnum,AREA,1
Generacin de las lneas a partir de la unin de dos keypoints
Genera lneas entre dos puntos diferentes
L, K1, K2
Tabla 6. Generacin de lneas a partir de los puntos P. R
LINEA
1
2
3
4
5
6
7
8
Ki
1
1
3
4
5
5
7
6
Kj
2
3
4
4
6
7
8
8
Generacin de reas definiendo lneas previamente

Genera reas entre lneas definidas, no es
necesario
numerar
las
lneas;
el
preprocesador de ansys automticamente
las enumera
AL, L1, L2, L3, L2,
Tabla 7. Generacin de reas entre lneas P.R
AREA
2
1
L1
1
5
L2
2
6
L3
3
7
L4
4
8
Construida el rea base para el proyecto, substraemos el rea que en teora va ha

ser la pila rectangular, la suposicin que se hace es la siguiente; se crea volumen
para el rea rectangular y el rea de la seccin de control, utilizando el comando
de extruccin, despus substraemos el volumen de la pila rectangular para que
65
solo se encuentre el rea de agua, (Vase figura 29) donde ocurren los
fenmenos fsicos objeto de estudio del proyecto.
VOFFST,1,7,
Genera un volumen normal al rea

creada, por extruccin en este caso el rea
1, con la altura que se especifica en este
caso 7. este paso se lleva a cabo con las
dems reas
Figura 29. Volumen para el rea rectangular
VOFFST,2,ti,
Genera
VSBV, NV1, NV2, SEPO, KEEP1,

KEEP2
Substraer volumen desde volumen.
un volumen normal al rea creada,

por extruccin en este caso el rea 2, con la
altura que se especifica en este caso Ti. este
paso se lleva a cabo con las dems reas
66
VSBV, 2,1,,,
Substraer
volumen desde volumen. En este

caso substraer del volumen 2, el volumen 1;
para dar forma al volumen de control
Figura 30. Forma del volumen de control por extruccin pila rectangular
Se guarda la el proceso realizado hasta aqu, por medio del siguiente

procedimiento en el men de utilidades (Utility Men) File > Save as Jobname.db
de esta forma guardamos en cualquiera momento en la base de datos de ANSYS
6.1.2 Discretizacin. En la fase del enmallando se transforma el elemento fsico

de estudio con incgnitas infinitas, a uno que tiene un sistema con un numero de
incgnitas finitas relacionada entre si por elementos volumtricos de volmenes
finitos, capaces transferir informacin entre si, por medio de la interaccin de lo
nodos comunes entre los elementos.
Para empezar con esta fase se muestran y se numeran las lneas del proyecto,
posteriormente seleccionamos las lneas que dan forma a la pila rectangular
dividindolas en elementos con longitudes iguales a una tasa reducida; necesaria
67
por refinamiento en el sitio de inters, donde en teora, ocurren los fenmeno

hidrulicos a estudiar, lo mismo hacemos con las lneas del volumen de control
pero con una tasa de divisin menor. Este proceso se lleva a cabo con el siguiente
comando.
Divide lneas especificando la tasa
LESIZE,NL1,SIZE,,,,,,,
/PNUM,
LINE, 1
LPLOT
La tasa de divisin para cada lnea se sirve en la siguiente tabla
Tabla 8. Tasa de divisin de lneas que componen el equilibrio P. R
LINE
SIZE
LINE
SIZE
1
2
3
4
7
17
18
19
20
21
22
23
0.3
0.3
0.3
0.3
0.08
0.3
0.3
0.3
0.3
0.3
0.3
0.3
24
25
26
27
28
29
30
31
32
33
34
35
0.3
0.08
0.08
0.08
0.08
0.08
0.08
0.08
0.08
0.08
0.08
0.08
VPLOT
VSWEEP, 3
Se guarda la el proceso realizado hasta aqu, por medio del siguiente

procedimiento en el men de utilidades (Utility Men) File > Save as Jobname.db
68
Figura 31. Equilibrio discretizado en la pila rectangular
6.1.3 Condiciones de frontera del flujo de entrada. Ahora establecemos las

condiciones de frontera, la cuales son; velocidad, pero primero debemos ver el
volumen por el cara ZY, para ello utilizamos el PlotCtrls > Pam Zomm Rotate. Esta
aplicacin nos sirve para girar en cualquier direccin (x, y, z) hasta llegar al perfil
deseado en este caso ZY, tambin a su vez para acercar y alejar el volumen por
medio de los botones con crculos; el circulo pequeo aleja y el circulo grande
aplica un zoom.
Despus aplicamos la velocidad de entrada, valor mximo de la distribucin

logartmica calculada, aplicando la carga por rea; en el fondo del lecho de rio la
velocidad es cero.
69
Figura 32. Giros en cualquier direccin
Se establece un bucle que aplica la velocidad del flujo en cada rea, gracias al
sistema de numeracin, las areas que compone el dominio del problema, poseen
un nmero de identificacin. El comando ejecutado es el siguiente:
*DO, Par, IVAL, FVAL, INC
Donde , *DO es el inicio de DO Loop,

Par es el contador , IVAL y FVAL son
los valores inicial y final de nodos del
bucle, INC es el incremento entre IVAL
y FVAL
*ENDDO
Finaliza el DO Loop y Empieza el

bucle ( hacerlo hasta FVAL)
70
*DO, T, 9,10, 1
T es el contador ,
DA, T, VX, 9, , , , , , , ,
Aplica
la
carga
en
el
rea
especificada, en este caso ser el
contador I, la velocidad en la dimensin
X, con un valor de 9 [m/seg]
*ENDDO
DA, 12, VX, 9, , , , , , , ,
El comando DA,
6.1.4 Condicin de frontera del flujo de salida. En la salida de presin, el fluido

sale del dominio computacional. Se especifica la presin esttica en el lado ZY de
salida igual, a la presin atmosfrica (presin manomtrica cero). Aplicamos los
siguientes comandos.
/PNUM, AREA,1
Activa el sistema de numeracin para

las reas
APLOT
Dibuja reas
DA, 11,PRES,0
Aplica carga en el rea especificada

en este caso 11, carga aplicar presin,
con un valor de 0
/PNUM, AREA,0
Desactiva el sistema de numeracin

para reas
6.1.5 Condicin de frontera en la superficie solida (pila). Debido a que, la pila

es una superficie solida, el fluido no puede infiltrarse en el alma de concreto, la
componente normal de la velocidad con respecto a la superficie solida a lo largo
de aquella cara en la cual se establece la condicin de frontera, se fija tambin
cero la componente de velocidad tangencial, consecuente a la condicin de no
deslizamiento.
71
*DO, Q,13,16, 1
Q es el contador ,
DA, Q, VX, 0, , , , , , , ,
Aplica
la
carga
en
el
rea
*ENDDO
Fin DO
*DO, P,13,16, 1
p es el contador ,
DA, P, VY, 0, , , , , , , ,
Aplica
la
carga
en
el
rea
Y, con un valor de 0 [m/seg]
*ENDDO
Fin DO
FINISH
Figura 33. Aplicacin de carga en el rea especificada en la pila rectangular
72
6.1.6 Solucin. Se ingresa el procesador solucin, para establecer la aceleracin

gravitacional en el problema debido a que, sin gravedad, el volumen de agua no
se desplazara.
/SOLU
Ingreso al procesador de solucin
ACEL,,-9.81,
Especifica
gravedad
la
aceleracin
de
la
Luego se procede a fijar las propiedades del fluido en el men interactivo ANSYS
Solution > Flotran Set up > Fluid Properties > OK.
En el label Density properity type CONSTANT escribimos la densidad de el agua

a 20C, la cual es 998.2 [Kg/m^3] y en el label Viscosity properity type
CONSTANT escribimos la viscosidad dinamica de el agua a 20C, la cual es
1.005E -03 [Kg/m^3].
Figura 34. Propiedades del fluido pila rectangular
73
Figura 35. Propiedades del flujo pila rectangular
Aumentamos el nmero de iteraciones y establecemos modificaciones en los

factores de relajacin inercial. Estas modificaciones se llevan a cabo, porque; En
el momento de correr el software, surga un mensaje que declaraba que la matriz
de coeficiente posea una diagonal negativa, creando un espurio en la solucin,
entregando resultados no fsicos y alterando la solucin del problema, entonces se
empleo esta tcnica de relajacin para la resolucin de flujos con alto grado de
complejidad. Luego se procede a correr el programa, para resolver el modelo.
FLDATA2, ITER, EXEC, 50
Numero de iteraciones, 50 en total
FLDATA34, MIR, MOME, 0.3
Modificacin de factor para momento,

el rango es de [0 -1], el utilizado en el
modelo es de 0.3
74
SOLVE
stars solucin
6.1.7 Postproceso. Se ingresa al postprocesador. Se lee el ltimo resultado;

General postproc > Read Result > Last set, luego se especifica al software que
entregue la simulacin del contorno de velocidad a travs del men interactivo;
General postproc > Plot Results > Contour plot > Nodal Solu > DOF solution >
Fluid velocity > OK (Vase figura 36).
El area de interes del proyecto, se encuentra cerca del lecho de rio, podemos
puntualizar en esa zona, para ver los efectos del flujo de agua en la pila
rectangular, moviendo el Working Plane a lo largo de la longitud axial en el lecho
del rio empleando el siguiente procedimiento.
En el men de utilidades de ANSYS dar click en WorkPlane > WP Settings a

continuacin seleccionar (Cartesian, Grid Only).
Figura 36. Solucin nodal de la velocidad del fluido en la pila rectangular
Cambiar los valores existentes que por defecto aporta ANSYS wp settings por
los siguientes.
75
Figura 37. Wp settings pila rectangular
Snap incr
: 0.1
Spacing
Minimum
:1
Maximum
: 10
Tolerance
: 0.01
Ir al ANSYS Utiliy Menu > WorkPlane y seleccionar Display Working

Plane.
Ahora que el working plane es seleccionado, ir a ANSYS Main
Menu>WorkPlane >Offset WP by Increments y ajustar el working plane Tl
que se encuentre situado a lo largo del lecho. Cuando est terminando el
movimiento, el plano que debers ver es el siguiente.
76
Figura 38. Movimiento de working plane en el lecho del rio en la pila

rectangular
Una vez el plano es una cuadricula lineal, vamos ANSYS Utiliy Menu
>PlotCntrls>Style>Hidden Line Options. (Vease la Figura 39), En la ventana
cambiar el Type of Plot a Q-Slice Z-buffer, y Cutting Plane is a Working Plane y
damos click en ok. ANSYS ahora nos mostrar los resultados del anlisis con el
working plane, como el dibujo es cortado.
Figura 39. Cambio a working plane para mostrar los resultados
77
El resultado final es el siguiente
Figura 40. Campo de velocidad en el lecho del rio en la pila rectangular
Posteriormente vamos Main Menu > General Postproc > Plot Results>Vector
Plot > Predefined, y en la venta que se presenta a continuacin damos click en
OK para mostrar el vector del gradiente de velocidad.
Figura 41. Vector del gradiente de velocidad en la pila rectangular
78
Para mostrar las lneas de corriente que son lneas tangentes al vector del
gradiente de velocidad, que nos permiten observar el movimiento del flujo en el
entorno de la pila, definimos los punto inciales del trazo del flujo, Main Menu >
General Postproc > Plot Results> Flow trace > Defi trace pt.
En la figura 41, Seleccionamos los puntos de inters donde se pretenda simular

el movimiento de las lneas de corriente. Despus utilizando la siguiente cadena
Main Menu > General Postproc > Plot Results> Flow Trace > Plot Flow Trace.
En la ventana de la figura 43 se elige DOF solution , VSUM. Se cambia el valor
de max numbers of loops por el numero de punto que desee simular, en este
paso se utiliz 50.
Figura 42. Puntos de inters donde se pretende simular el movimiento de las

lneas de corriente
79
Figura 43. Plot flow trace pila rectangular
En el momento que se da click en OK, ANSYS simula el movimiento de las lneas

de corriente. Se puede llevar a cabo, animacin virtual de la simulacin entregada
anteriormente.
Figura 44. Trayectoria de las lneas de corriente en la pila rectangular
80
Ir al ANSYS Utiliy Menu > PlotCtrls > Animate > Particle Flow y seleccionar
de la ventana DOF solution, VSUM. Luego dar click en Ok y se tendr la
animacin virtual del flujo de partculas de la figura 44.
Perfil de presiones sobre el lecho del rio, en el entorno de la pila se puede hallar
mediante la siguiente cadena en el men interactivo de ANSYS. General
postproc > Plot Results > Contour plot > Nodal Solu > DOF solution > Press>
OK.
Figura 45. Perfil de presiones sobre el lecho del rio en la pila rectangular
6.2 PILA CIRCULAR APORTICADA

81
/FILNAME, Pila circular

pila circular en la base de datos
/TITLE, Modelacin tridimensional del Indica el titulo del problema

flujo de agua en la pila circular
/UNITS,SI
Se establece el sistema internacional

de unidades
A continuacin se introducen parmetros escalares

que se utilizaran
D0= 0.80
Se define como dimetro de la pila

0.80 [m]
lu= 2.20
luz entre ejes en las pilas circulares

que componen la pila aporticada
Ti= 6
Profundidad del flujo en lamina libre

aguas arriba
ro=0.40
Se define como radio de la pila

0.40 [m]

de 20 C.
82
/ PREP7

finitos
ET,1, FLUID 142

FLUID 142 cuadro de dialogo
MP,DENS,1,1000
Propiedad del material, densidad, del

Kg m
de 998,2
3
MP,VISC,1,1.10E-03
Propiedad del material, viscosidad, del

elemento fluid 142, para 20 C es de
N seg
1.005E-03
2

incomprensible, en el menu principal de ANSYS (ANSYS main men)
de flujo laminar a turbulento.

geometra de la pila circular, se cambiara la vista de los ejes de coordenadas, la
orientacin del sistema de coordenadas globales, para visualizar en forma ms
sencilla el step by step en la generacin de la pila, Se utilizarn los siguientes
comandos:
/VIEW,1,1,1,1
!*Define la direccin para la ventana de

salida (output)
/VUP,1,Y
!*Orientacin
del
sistema
coordenadas para el proyecto
/REPLOT
!*Dibuja automticamente los cambios

83
de
hechos en la ventana de salida

posteriormente cambiara a volumen por extruccin se proporciona en la tabla 9,
aclarando que todas las dimensiones se encuentran en metros.
Genera el keypoint, enumera, se

especifica las coordenadas x, y, z
respectivamente
K, nm., x, y, z
Tabla 9. Coordenadas para dar forma al rea de pila circular aporticada
KEYPOINT
2+Do+2
2+2*ro+lu+5
-0,07
2+2*ro+lu+5
-0,07
2+Do+2
-0.0196
-0.0196
2+Do
2+Do
-0.0196
2+Do
-0.0196
2+Do
2+ ro
-0.0196
2+ ro
10
2+ro+lu
-0.035
2+ ro
11
2+Do+1.40
-0.035
12
2+Do+1.40
-0.035
2+Do
13
2+Do+1.40+Do
-0.035
14
2+Do+1.40+Do
-0.035
2+Do
84
KSEL,ALL
Selecciona todos los keypoints
LSEL,ALL
Selecciona todas las lneas
KPLOT,ALL
Dibuja todos los keypoints
LPLOT,ALL
Dibuja todas las lneas
/Pnum,KP,1
Numera todos los keypoints
/Pnum,LINE,1
Numera todas las lneas
/Pnum,AREA,1
Numera todas las reas
Figura 46. Generacin de puntos para crear rea circular aporticada
85
Generacin de las lneas a partir de la unin de dos keypoints
Genera lneas
diferentes
L, K1, K2
entre
Tabla 10. Generacin de lneas a partir de los puntos P. C
LINEA
1
2
3
4
5
6
7
8
9
10
11
12
Ki
1
1
3
2
5
5
7
6
11
11
13
12
Kj
2
3
4
4
6
7
8
8
12
13
14
14
86
dos
puntos
(Vase en la figura 47), una serie de puntos keypoint que forman la pila, pero de
geometra cuadrada; la realidad del caso es que se generan keypoints en forma
cuadrada es el primer paso, luego se crea las lneas 5, 6, 7, 8, 9, 10, 11, 12 y
posteriormente con el siguiente cadena de comandos creamos las lneas que
componen la seccin circular para la pila.
LFILLT, NL1, NL2, RAD, PCENT
genera una cinta o curva en la

intercepcin
de
dos
lneas
especificando las lneas que se
interceptan, radio de la curva.
LFILLT, 5, 6, r0, 9
genera la curva en la intercepcin de la

lneas 5 y 6, con radio de 0.4 [m]
introducido como parmetro escalar
anteriormente y eligiendo el keypoint 9
como centro, se realiza el mismo
procedimiento para las dems lneas.
LFILLT, 5, 8, r0, 9
Figura 48. Lneas que componen la seccin circular
87
LFILLT, 6, 7, r0, 9
LFILLT, 7, 8, r0, 9
LFILLT, 9, 10, r0, 10
LFILLT, 12, 9, r0, 10
LFILLT, 10, 11, r0, 10
LFILLT, 12, 11, r0, 10
Genera reas entre lneas definidas
AL, L1, L2, L3, L4,
Tabla 11. Generacin de reas entre lneas P. C
AREA
1
2
3
L1
7
5
1
L2
8
6
2
L3
9
13
3
L4
10
14
4

ser la pila circular a porticada (Vase la figura 49), la suposicin que hacemos es
la siguiente; se crea volumen para las reas circulares y el rea de la seccin de
control, utilizando el comando de extruccin, despus substraemos el volumen de
las pilas circulares para que solo se encuentre el volumen de agua, figura 50
donde ocurren los fenmenos fsicos objeto de estudio del proyecto.
VOFFST,1,7,
Genera un volumen normal al rea creada,

paso se lleva a cabo con las dems reas.
88
VOFFST,2,7,
Figura 49. Pila circular aporticada
VOFFST,3,Ti,
Figura 50. Pila circular aporticada interactuando con el equilibrio
89

KEEP2
VSBV, 3, 1,,,
Substraer volumen desde volumen. En este

caso substraer del volumen 3, el volumen 1;
para dar forma al volumen de control.
VSBV, 4, 2,,,
Se guarda el proceso realizado hasta aqu, por medio del siguiente procedimiento
en el men de utilidades (Utility Men) File > Save as Jobname.db de esta forma
guardamos en cualquiera momento en la base de datos de ANSYS.
90

posteriormente seleccionamos las lneas que dan forma a la pila circular
comando.
/PNUM, LINE, 1
LPLOT
La tasa de divisin para cada lnea se sirve en la siguiente tabla 10
VPLOT
VSWEEP, 3
91
Tabla 12. Tasa de divisin de lneas que componen el equilibrio P.C
LINE
1
2
3
4
29
30
31
32
33
34
35
36
37
38
39
40
41
42
SIZE
0.3
0.3
0.3
0.3
0.3
0.3
0.3
0.3
0.3
0.3
0.3
0.3
0.08
0.08
0.08
0.08
0.08
0.08
LINE
7
8
9
10
11
12
15
16
17
18
19
20
43
44
45
46
47
48
SIZE
0.08
0.08
0.08
0.08
0.08
0.08
0.08
0.08
0.08
0.08
0.08
0.08
0.08
0.08
0.08
0.08
0.08
0.08
Figura 52. Equilibrio discretizado en la pila circular aporticada
92
en el men de utilidades (Utility Men) File > Save as Jobname.db

condiciones de frontera, la cuales son; velocidad, aplicamos la velocidad de
entrada, valor mximo 9[m/s] de la distribucin logartmica calculada aplicando la
velocidad por rea; en el fondo del lecho de ro la velocidad es cero.
sistema de numeracin, las reas que compone el dominio del problema, poseen
un nmero de identificacin. El comando ejecutado es el siguiente.
Donde , *DO es el inicio de DO Loop,

Par es el contador , IVAL y FVAL son
y FVAL
*ENDDO

bucle ( hacerlo hasta FVAL)
*DO, T, 9,10, 1
T es el contador ,
DA, T, VX, 9, , , , , , , ,
Aplica la carga en el rea

*ENDDO
DA, 12, VX, 9, , , , , , , ,
93

/PNUM, AREA,1

las reas
APLOT
Dibuja reas
DA, 17,PRES,0

con un valor de 0.
/PNUM, AREA,0

para reas
6.2.5 Condicin de frontera en la superficie solida (pila). Debido a que, la pila

deslizamiento.
*DO, Q,19,22, 1
Q es el contador ,
DA, Q, VX, 0, , , , , , , ,
Aplica
la
carga
en
el
rea
contador, la velocidad en la dimensin
*ENDDO
Fin DO
*DO, P,19,22, 1
P es el contador ,
94
DA, P, VY, 0, , , , , , , ,
Aplica
la
carga
en
el
rea
Y, con un valor de 0 [m/seg]
*ENDDO
Fin Do
FINISH
Figura 53. Aplicacin de carga en el rea especificada en la pila circular

aporticada

se desplazara.
95
/SOLU
ACEL,,-9.81,
Especifica
gravedad
la
aceleracin
de
la

CONSTANT escribimos la viscosidad dinmica de el agua a 20C, la cual es
1.005E -03 [Kg/m^3], posteriormente damos click en OK.


el rango es de [0 -1], el utilizado en el
modelo es de 0.3
SOLVE
Stars solucin

General postproc > Read Result > Last set, luego se le especifica al software que
entregue la simulacin del contorno de velocidad a travs del men interactivo;
Fluid velocity > OK (Vase la figura 54).
96
puntualizar en esa zona como se realizo en seccion anterior, para ver los efectos
del flujo del agua en la pila punta diamante.
Figura 54. Solucin nodal de la velocidad del fluido en la pila circular

aporticada
Figura 55. Campo de velocidad en el lecho del rio en la pila circular

aporticada
97
Figura 56. Vector del gradiente de velocidad en la pila circular aporticada
Para mostrar las lneas de corriente, que son lneas tangentes al vector del
gradiente de velocidad, permitiendo observar el movimiento del flujo en el entorno
de la pila, definimos los punto inciales del trazo del flujo, Main Menu > General
Postproc > Plot Results> Flow trace > Defi trace pt.
En la figura 56, seleccionamos los puntos de inters donde se pretenda simular el

movimiento de las lneas de corriente. Despus utilizando la siguiente cadena
Main Menu > General Postproc > Plot Results> Flow trace > Plot flow trace.
En la ventana del figura 58, se elige DOF solution, VSUM, y se da click en OK.
98

lneas de corriente
Figura 58. Plot flow trace pila circular
99
ANSYS simula el movimiento de las lneas de corriente. Se puede llevar a cabo,

animacin virtual de la simulacin entregada anteriormente.
animacin virtual del flujo de partculas (Vase la figura 59).
Figura 59. Trayectoria de las lneas de corriente en la pila circular aporticada
Perfil de presiones sobre el lecho del rio y la pila se puede hallar mediante la
siguiente cadena en el men interactivo de ANSYS. General postproc > Plot
Results > Contour plot > Nodal Solu > DOF solution > Press > OK.
100
Figura 60. Perfil de presiones sobre el lecho del rio en la pila circular
aporticada
6.3 PILA EN PUNTA DE DIAMANTE

/FILNAME, Pila en punta de diamante

pila circular en la base de datos.
/TITLE, Modelacin tridimensional del Indica el titulo del problema.

flujo de agua en la pila en punta de
diamante
/UNITS,SI
Se establece el sistema internacional

101
de unidades
A continuacin se introducen parmetros escalares que se utilizaran

Ba= 0.8
Se define como base de la seccin de

la pila 1. [m]
He= 1.62
Altura de la seccin de la pila 3.0 [m]
Ti= 6
Profundidad del flujo en lamina libre

aguas arriba.
Hd= 0.69
Altura del triangulo que genera la forma

de punta de diamante.

de 20 C.
/ PREP7

finitos
ET,1, FLUID 142

FLUID 142 cuadro de dialogo
102
MP,DENS,1,998.2
MP,VISC,1,1.005E-03
Propiedad del material, densidad,

elemento fluid 142, para 20 C
Kg m
de 998,2
3
m
Propiedad del material, viscosidad,
elemento fluid 142, para 20 C
N seg
de 1.005E-03
2
del
es
del
es

incomprensible, en el men principal de ANSYS (ANSYS main men)
de flujo laminar a turbulento.

geometra de la pila en punta de diamante, se cambiara la vista de los ejes de
coordenadas, la orientacin del sistema de coordenadas globales, para visualizar
en forma ms sencilla el step by step en la generacin de la pila, Se utilizarn los
/VIEW,1,1,1,1
Define la direccin para la ventana de

salida (output).
/VUP,1,Y
Orientacin
del
sistema
coordenadas para el proyecto.
/REPLOT
Dibuja automticamente los cambios

hechos en la ventana de salida.
de

posteriormente cambiara a volumen por extruccin se proporciona en la siguiente
tabla 11, aclarando que todas las dimensiones se encuentran en metros.
K, nm., x, y, z
Genera el keypoint, se especifica las

coordenadas x, y, z respectivamente.
103
Tabla 13. rea Coordenadas para dar forma al rea de pila punta diamante
KEYPOINT
1
2
3
4
5
6
7
8
9
10
X
0
0
2+He+2*hd+5
2+He+2*hd+5
2
2+hd
2+hd
2+hd+He
2+hd+He
2+2*hd+He
Y
0
0
-0.07
-0.07
-0.035
-0.035
-0.035
-0.035
-0.035
-0.035
Z
0
2+Ba+2
0
2+Ba+2
2+Ba/2
2
2+Ba
2
2+Ba
2+Ba/2
KSEL,ALL
!*Selecciona todos los keypoints
LSEL,ALL
!*Selecciona todas las lneas
KPLOT,ALL
!*Dibuja todos los keypoints
LPLOT,ALL
!*Dibuja todas las lneas
/Pnum,KP,1
!*Numera todos los keypoints
/Pnum,LINE,1
!*Numera todas las lneas
/Pnum,AREA,1
!*Numera todas las reas
104
Figura 61. Generacin de puntos para crear el rea de punta Diamante
Generacin de las lneas a partir de la unin de dos keypoints.
Genera lneas
diferentes
L, K1, K2
entre
Tabla 14. Generacin de lneas a partir de los puntos P. P. D
LINEA
1
2
3
4
5
6
7
8
9
10
Ki
1
1
3
2
5
5
7
6
8
9
105
Kj
2
3
4
4
6
7
9
8
10
10
dos
puntos
Genera reas entre lneas definidas
AL, L1, L2, L3, L4,
Tabla 15. Generacin de reas entre lneas P. P. D
AREA
L1
L2
L3
L4
L5
L6
10

ser la pila en punta de diamante maciza, (Vase la figura 63), la suposicin que
hacemos es la siguiente; se crea volumen para el rea en punta de diamante y el
rea de la seccin de control, utilizando el comando de extruccin, despus
substraemos el volumen de la pila en punta diamante para que solo se encuentre
el rea de agua, (Vase la figura 64).
106
Donde ocurren los fenmenos fsicos objeto de estudio del proyecto.

VOFFST,1,7,
Genera un volumen normal al rea creada,

paso se lleva a cabo con las dems areas.
Figura 63. Pila en punta de diamante
VOFFST,2,Ti,
107
Figura 64. Equilibrio interactuando con la pila

KEEP2
VSBV, 2, 1,,,
Substraer volumen desde volumen. En

este caso substraer del volumen 2, el
volumen 1; para dar forma al volumen
de control.
VSBV, 4, 2,,,
en el men de utilidades (Utility Men) File > Save as Jobname.db de esta forma
guardamos en cualquiera momento en la base de datos de ANSYS.
108

posteriormente seleccionamos las lneas que dan forma a la pila rectangular
comando.
109
/PNUM, LINE, 1
LPLOT
La tasa de divisin para cada lnea se sirve en la siguiente (Vase la tabla 16).
Tabla 16. Tasa de divisin de lneas que componen el equilibrio P. P. D
LINE
SIZE
LINE
SIZE
0.3
34
0.08
0.3
35
0.08
0.3
36
0.08
0.3
37
0.08
23
0.3
38
0.08
24
0.3
39
0.08
25
0.3
40
0.08
26
0.3
41
0.08
27
0.3
42
0.08
28
0.3
43
0.08
29
0.3
44
0.08
30
0.3
45
0.08
31
0.08
46
0.08
32
0.08
47
0.08
33
0.08
48
0.08
VPLOT
VSWEEP, 3
110
Figura 66. Equilibrio discretizado en la pila punta diamante
en el men de utilidades (Utility Men) File > Save as Jobname.db.

condiciones de frontera, la cuales son; velocidad, aplicamos la velocidad de
entrada, valor mximo 9[m/s] de la distribucin logartmica calculada aplicando la
velocidad por rea; en el fondo del lecho de rio la velocidad es cero.
sistema de numeracin, las areas que compone el dominio del problema, poseen
un nmero de identificacin. El comando ejecutado es el siguiente.
Donde, *DO es el inicio de DO Loop,

Par es el contador, IVAL y FVAL son
y FVAL.
*ENDDO

bucle (hacerlo hasta FVAL).
111
*DO, T, 9,10, 1
T es el contador ,
DA, T, VX, 9, , , , , , , ,
Aplica
la
carga
en
el
rea
X, con un valor de 9 [m/seg].
*ENDDO
DA, 12, VX, 9, , , , , , , ,
El comando DA,

/PNUM, AREA,1

las reas.
PLOT
Dibuja reas
DA, 11,PRES,0

con un valor de 0.
/PNUM, AREA,0

para reas.
6.3.5 Condicin de frontera en la superficie solida (pila). Debido a que la pila

deslizamiento.
112
*DO, Q,13,16, 1
Q es el contador ,
DA, Q, VX, 0, , , , , , , ,
Aplica
la
carga
en
el
rea
X, con un valor de 0 [m/seg].
*ENDDO
Fin DO
*DO, P,13,16, 1
p es el contador ,
DA, P, VY, 0, , , , , , , ,
Aplica
la
carga
en
el
rea
Y, con un valor de 0 [m/seg].
*ENDDO
Fin DO
FINISH
Figura 67. Aplicacin de carga en el rea especificada en la pila

diamante
113
punta

se desplazara.
/SOLU
ACEL,,-9.81,
Especifica
gravedad
la
aceleracin
de
la

CONSTANT escribimos la viscosidad dinmica de el agua a 20C, la cual es
1.005E -03 [Kg/m^3], posteriormente damos click en OK.


el rango es de
[0 -1], el utilizado
en el modelo es de 0.3
SOLVE
stars solucin
114

General postproc > Read Result > Last set, luego se le especifica al software que
entrgue la simulacin del contorno de velocidad a travs del men interactivo;
Fluid velocity > OK. (Vase la figura 68).
puntualizar en esa zona realizando el procedimiento aplicado para la pila
rectangular, para ver los efectos del flujo del agua en la pila punta diamante.
Figura 68. Solucin nodal de la velocidad del fluido en la pila punta diamante
115
Figura 69. Campo de velocidad en el lecho del rio en la pila punta diamante
Figura 70. Vector del gradiente de velocidad en la pila punta diamante
116
Para mostrar las lneas de corriente que son lneas tangentes al vector del
gradiente de velocidad, que nos permiten observar el movimiento del flujo en el
entorno de la pila, definimos los punto inciales del trazo del flujo, Main Menu >
General Postproc > Plot Results> Flow trace > Defi trace pt.
En la figura 70 Seleccionamos los puntos de inters donde se pretenda simular el

movimiento de las lneas de corriente. Despus utilizando la siguiente cadena
Main Menu > General Postproc > Plot Results> Flow trace > Plot flow trace.
En la ventana de la se elige DOF solution, VSUM, y se da click en OK.

lneas de corriente
ANSYS simula el movimiento de las lneas de corriente. Se puede llevar a cabo,

animacin virtual de la simulacin entregada anteriormente.
animacin virtual del flujo de partculas de la figura 71.
117
Figura 72. Trayectoria de las lneas de corriente en la pila punta diamante
Perfil de presiones sobre el lecho del rio y la pila se puede hallar mediante la
siguiente cadena en el men interactivo de ANSYS. General postproc > Plot
Results > Contour plot > Nodal Solu > DOF solution > Press > OK.
Figura 73. Perfil de presiones sobre el lecho del rio en la pila punta diamante
118
7. ANALISIS DE RESULTADOS
Se presenta a continuacin la variacin de las componentes de velocidad y la

fuerza msica (presin), dentro del volumen aislado de fluido en cierto instante de
tiempo, mediante la interfaz grafica. Debido al continuo cambio de los
desplazamientos, es natural que concentremos nuestra atencin en las
velocidades y presin en el entorno de la pila, que se obtuvieron en el men
interactivo ANSYS con operaciones de trayectoria, donde se especifica la
direccin donde se requiere obtener los resultados.
Grafica 1. Gradiente de velocidad aguas arriba normal al flujo
14.000
VELOCIDAD AGUAS ARRIBA [m/seg]
13.000
12.000
11.000
10.000
9.000
8.000
7.000
6.000
5.000
4.000
3.000
2.000
1.000
0.000
0.00
1.00
2.00
3.00
4.00
5.00
DISTANCIA NORMAL AL FLUJO [m]

Punta de Diamante
Circular
Rectangular
En la grafica 1, la velocidad del flujo en el contorno de la pila rectangular y circular,

aguas arriba, se redujo en 96[%] y 86[%] con respecto a la corriente media del
flujo, debido a la presencia de estancamiento del flujo, generando velocidades
bajas por consiguiente presiones dinmicas altas Bernoulli, explicacin lgica a la
119
formacin de vrtices de herradura frontales aguas arriba de la pila, que forman el

principal cuenco de socavacin(Vase la Figura 45, 60 y 73).
Se puede observar la mxima aceleracin del flujo en la pila rectangular; con un

valor de 12.4 [m/seg], consecuencia de la separacin en la trayectoria del flujo
incapaz de seguir la forma geomtrica de la pila. Aguas arriba la pila en punta de
diamante, es la seccin que altero en menor magnitud el flujo del agua
manteniendo un rango de velocidades de [6 9 m/seg].
Grafica 2. Gradiente de velocidad aguas abajo normal al flujo
14.000
13.000
VELOCIDAD AGUAS ABAJO [m/seg]
12.000
11.000
10.000
9.000
8.000
7.000
6.000
5.000
4.000
3.000
2.000
1.000
0.000
0.00
1.00
2.00
3.00
4.00
5.00
DISTANCIA NORMAL AL FLUJO [m]

Punta de Diamante
Rectangular
Circular
(Vase la grafica 2), el comportamiento del flujo, donde se crea una zona de
aguas muertas (velocidades bajas) aguas abajo de las pilas simuladas en el
proyecto. (Vase la Figuras 40, 55, 69). Las pilas con llevan a que se presenten,
cambios en el rgimen de flujo, actuando como una seccin de control (rgimen
rpido a lento).
120
Figura 74. Gradientes de velocidad prxima a la frontera solida para las

formas geomtricas analizadas paralelas al flujo
Figura 74 A Pila rectangular
Figura 74 B Pila circular aporticada
Figura 74 C Pila punta de diamante
121
Se concluye que las mayores velocidades del flujo del agua se presentan en las
zonas ms prximas de lados laterales de las pilas (franja de color rojo), dichas
velocidades son localmente mayores que la velocidad de la corriente media,
siendo el mximo en la pila rectangular, a una distancia de 2 [m], el mnimo en la
geometra circular a la distancia de 2.40 [m] (Vase la Figura 74).
La pila en punta de diamante presenta la mayor velocidad a la distancia de 3 [m],

favorable porque se reduce la vorticidad lateral, logrando adems desplazar los
vrtices en longitud aguas abajo de la pila figura 74C, en las zonas alejadas el
flujo no es alterado. La pila circular aporticada disminuye la capacidad de
transporte del equilibrio (volumen aislado de fluido) al obstaculizar el fluido en la
primera pila, creando vrtices laterales y estelas en medio del sistema aporticado,
que socavan el lecho de rio, al existir mayor rea normal al flujo (Vase la
Figura 74).
(Vase la Figura 75). Como todas las pilas crean, crculos de presiones dinmicas
positivas, de color rojo situados en la cara frontal de cada pila, que extraen
material de lecho del rio, distinguiendo la mayor presin ejercida en magnitud y
extensin es la generada por la pila rectangular.
(Vase la Figura 76). Se verifica como despus del punto de separacin, hay un
incremento continuo de la velocidad de la corriente principal, tambin se crea una
regin gruesa irregular con estelas y corrientes secundarias llamada regin
separada, inicialmente en la zona separada se genera un par de vrtices, estos
vrtices actan como rodillos hidrodinmicos por los cuales fluye la corriente
principal, los vrtices en la estela son componentes del sistema de vrtices que
erosionan el lecho de rio.
El punto se estancamiento es el punto en el campo de flujo donde la velocidad es

cero, las lneas de corriente en las condiciones de frontera de entrada, nunca se
cortan porque el flujo es permanente, por consiguiente el flujo que se origina es
irrotacional y se mantiene constante hasta que se presenta el rea frontal de la
pila en su trayectoria, haciendo rotar las partculas de fluido generando los vrtices
laterales alrededor de la pila. Como muestra la (Vase la Figura 76).
122
Figura 75. Perfil de presin para las formas geomtricas analizadas sobre el
lecho del rio
Figura 75 A Pila rectangular
Figura 75 B Pila circular aporticada
Figura 75 C Pila punta de diamante
123
Figura 76. Fenmenos hidrulicos
Punto de
estancamiento
Corriente secundaria
Estela
Punto de separacin
Cuando un fluido se mueve sobre un cuerpo solido, ejerce fuerzas de presin y

fuerzas de corte paralela a lo largo de la superficie exterior del cuerpo, existe
inters en la resultante de presin exclusivamente, porque en este proyecto se
despreciara los efectos viscosos del agua por consiguiente se asume que las
fuerzas de corte son cero. La resultante que acta en la direccin de flujo sobre el
rea frontal de la pila se llama fuerzas de arrastre de presin o forma.
El calculo de la fuerza de arrastre se llevo a cabo, con la ecuacin 4. La ecuacin

es f (CD, , v, A), cada variable de la ecuacin es descritas en detalle, en la
seccin 2.8. El coeficiente de arrastre depende de la forma del cuerpo y
orientacin del cuerpo respecto al flujo, en la tabla 11-1 del libro YUNUS A.
CENGEL, JOHN CIMBALA, la figura de barra rectangular con esquinas agudas es
equivalente en forma y orientacin a la pila rectangular, por eso eligi el
coeficiente de arrastre de 1.3 [-] para una relacin L/D de 3.75. En la tabla 17, se
observa la forma circular con coeficiente de arrastre de 0.3 [-] para flujo turbulento
y la forma triangular con coeficiente de 0.48 [-].
124
Debido a que, las pilas simuladas poseen la misma base en planta el rea
proyectada o normal al flujo ser igual, la velocidad de aproximacin de corriente
libre es de 9 [m/seg], la densidad del agua para temperatura estndar de 20 C,
es de 998.2 [Kg-m/m3], teniendo todas las variables definidas, (Vase en la
tabla 17) el valor de arrastre de forma para cada seccin de pila.
Tabla 17. Clculo de arrastre de forma
PILA
Rectangular
Circular
Punta de Diamante
Base
[m]
0,8
0,8
0,8
Tirante
[m]
6
6
6
Densidad
[Kg/m^3]
998,2
998,2
998,2
Velocidad
[m/ seg]
9
9
9
CD
[-]
1,6
0,3
0,48
Area
[m^2]
4,8
4,8
4,8
FD
[KN]
310,48
58,22
93,14
Con los resultados obtenidos en la tabla 17; la pila con menor empuje dinmico es
la circular y la pila con ngulos recto presenta el mayor empuje dinmico de agua.
125
8. GEOMETRIA PROPUESTA PARA OPTIMIZACIN HIDRULICA DE

PILARES DE PUENTE
El estudio del diseo de pilas que afecten en menor magnitud las condiciones
hidrulicas del cauce son recientes, los principales mecanismos empleados para
lograr este fin son; el diseo de pilas con formas geomtricas hidrodinmicas, la
reduccin de vorticidad y corrientes secundarias. Esta metodologa consiste en
implementar elementos no estructurales a la pila para la disipacin de energa.
Las pilas simuladas en el proyecto, presentan un importante obstculo al flujo que

se reflejo en el campo de velocidades, presentando zonas donde sucedieron
variaciones, en la velocidad del flujo como se ha visto en el apartado anterior. De
acuerdo con los resultados provenientes de las tres simulaciones se propone
carenar las geometras estudiadas, cumpliendo con las condiciones que se
presentan a continuacin.
La forma de la pila deber ser currentilneo o hidrodinmica. Para alinear su

forma con las lneas de corriente anticipadas al flujo.
El ngulo de incidencia con respecto a la corriente del cauce debe ser cero.
La geometra hidrodinmica de la pila debe iniciar a una profundidad de

desplante de 1 [m] hasta el nivel mximo del agua esperado en la avenida del
rio.
Reduccin de vorticidad y corrientes secundarias.
Trasladar el punto de separacin hacia aguas debajo de la pila, la mayor

distancia posible. Por consiguiente disminucin de la regin separada.
Velocidad aguas arriba de la pila no debe ser menor que el 65[%] del valor
superficial de la corriente libre.
126
Reduccin de socavacin local terica en 30 [%].
Se simularon cuatro nuevas geometras para el elemento estructural (pila),

variando su espesor y longitud, denominadas P1 (Vase la figura 77A), P2
(Vase la figura 77B), P3 (Vase la figura 77C) y P4(Vase la figura 77D). La
forma de estas geometras, fue escogida de acuerdo al perfil ideal supersnico Biconvex Airfoil y se semeja a la estudiada por Roca, que planteo la geometra con
la forma del cuenco de socavacin en sus estudios. Las cuatro formas simuladas
cumplieron con las condiciones provenientes de los resultados de anlisis.
Figura 77. Propuesta de las nuevas geometras para optimizacin hidrulica

Figura 77 A. P1
Figura 77 B. P2
Figura 77 C. P3
Figura 77 D. P4
127
Con los datos de salida de las componentes de velocidad, aguas arriba y aguas
abajo, isocontornos de presin, se realiza un proceso exportacin de resultados a
un procesador de texto y desde all a una hoja de clculo para realizar sus graficas
y anlisis pertinentes. Los resultados del campo de velocidad y perfil de presiones,
provenientes de las simulaciones a los perfiles P1, P2, P3 y P4, se exteriorizan en
las siguientes graficas, para visualizar cual geometra es la que perturba en menor
grado el flujo del agua.
Tabla 18. Componentes de velocidad aguas arriba en lecho de rio para las
nuevas geometras propuestas
P1
S
0
0,24
0,48
0,72
0,96
1,2
1,44
1,68
1,92
2,16
2,4
2,64
2,88
3,12
3,36
3,6
3,84
4,08
4,32
4,56
4,8
V. A ARRIBA
9
9,55
9,68
9,59
9,57
9,53
9,43
9,16
8,76
7,77
4,1
8,09
9,06
9,4
9,55
9,62
9,64
9,58
9,67
9,52
9
P2
S
0
0,24
0,48
0,72
0,96
1,2
1,44
1,68
1,92
2,16
2,4
2,64
2,88
3,12
3,36
3,6
3,84
4,08
4,32
4,56
4,8
P3
V. A ARRIBA
9
9,33
9,35
9,37
9,35
9,28
9,15
8,92
8,48
7,62
7,94
8,56
8,9
9,1
9,25
9,32
9,34
9,31
9,32
9,27
9
128
S
0
0,24
0,48
0,72
0,96
1,2
1,44
1,68
1,92
2,16
2,4
2,64
2,88
3,12
3,36
3,6
3,84
4,08
4,32
4,56
4,8
V. A ARRIBA
9
9,54
9,55
9,5
9,54
9,49
9,41
9,24
8,85
8,15
4,41
8,02
8,81
9,18
9,35
9,44
9,49
9,49
9,5
9,5
9
P4
S
0
0,24
0,48
0,72
0,96
1,2
1,44
1,68
1,92
2,16
2,4
2,64
2,88
3,12
3,36
3,6
3,84
4,08
4,32
4,56
4,8
V. A ARRIBA
9
9,7
9,83
9,68
9,65
9,45
9,34
9,05
8,54
7,3
3,4
7,69
8,64
8,92
9,1
9,25
9,43
9,53
9,65
9,3
9
Grafica 3. Variaciones de las velocidades aguas arriba de las pilas
VELOCIDAD AGUAS ARRIBA [m/seg]
10
9
8
7
6
5
4
3
2
1
0
1
P1
P2
P3
8
P4
10
11
12
13
14
15
16
17
18
P4
19
P3
20
P2
21
P1
En la grafica 3, la geometra que menos variaciones presenta en la componente

de velocidad respecto de la corriente libre, es el P2, debido a que, las
fluctuaciones de velocidades aguas arriba de la pila, se encuentra en una rango
[9 7.6 m/seg], cumpliendo con la condicin donde se indica que, la velocidad
aguas arriba de la pila, no debe ser menor que el 65[%] del valor superficial de
velocidad en la corriente libre.
La distribucin de presin observada fue; al entrar el flujo en contacto con las

diferentes formas geomtricas, se elev hasta un valor mximo en la cara frontal
de la pila y se hizo menor alrededor del obstculo, pero conservando una simetra
en el contorno de la pila, finalmente las presiones inferiores se presentaron en la
parte posterior de cada geometra. La geometra con menor presin dinmica y
extensin radial fue P2, como virtud de de su forma, que presenta alineacin con
las lneas de corrientes anticipada y su velocidad que varia solo en un 20[%] de la
corriente libre. (Vase la figura 78).
129
Figura 78. Isocontorno de presin en geometras propuestas

Figura 78 A. P1
Figura 78 B. P2
Figura 78 C. P3
Figura 78 D. P4
De los resultado anteriores, se concluye que la geometra, que en menor grado

afecta la dinmica del flujo es el P2, a continuacin se dan las respectivas
medidas del espesor y longitud del perfil, el perfil debe poseer su mximo espesor
en el centroide de la pila x , que a su vez, es la distancia de 0.25 veces la Base
de cada lado paralelo al flujo. El punto inicial y final del perfil se ubican en la mitad
de la pila y , se encuentra a una distancia efectiva de 0.6 veces la altura de cada
lado normal al flujo, como se puede ver en la figura 79. En el anexo A se deja la
gua de comandos para simular en ANSYS, el perfil para cualquier forma
geomtrica de pila.
130
Figura 79. Esquema general P2
H/2
0 ,25B
0,60H
B
B/2
0 ,25B
El esquema anterior presenta a P2, currentilneo o hidrodinmico la finalidad de

esta forma es retrasar el punto de separacin por consiguiente disminuir el
arrastre de presin, consecuente la regin separada y desprendimiento de partes
de fluido circundante (Vrtices de estela) se minimizan reduciendo la turbulencia
en el contorno de la pila. Figura 80
Figura 80. Campo de velocidad y lneas de corriente en P2
La presencia de la forma hidrodinmica debilita la accin de los vrtices de

herradura, consecuente al flujo que es capaz de seguir la trayectoria del cuerpo
disminuyendo la divergencia de las lneas de corrientes. Las corrientes
131
secundarias generadas en el flujo separado (regin azul) figura 88, solo pueden
actuar sobre el lecho mediante vrtices laterales ya bien distanciados de la cara
posterior de la pila rectangular original.
A continuacin se presenta el cuadro de clculo para la socavacin local terica

en cada geometra.
Tabla 19. Estimacin de la socavacin local
PILA PUNTA DE DIAMANTE

FR
ds
1,131
1,71
1,008
1,69
0,813
1,64
0,998
1,69
1,137
1,72
1,228
1,73
PILA CIRCULAR
FR
1,263
0,865
0,192
0,885
1,279
1,401
PILA RECTANGULAR
FR
ds
1,679
1,80
0,659
1,60
0,072
1,21
0,670
1,60
1,693
1,80
1,551
1,78
ds
1,74
1,66
1,37
1,66
1,74
1,76
PILA P2
FR
1,145
1,028
1,071
1,156
1,202
1,230
ds
1,23
1,21
1,22
1,23
1,24
1,24
Basados en los resultados de socavacin local podemos concluir que el perfil

geomtrico P2 reduce tericamente 30[%] de la socavacin presentadas para las
pilas rectangular, circular y punta de diamante.
132
9. CONCLUSIONES
Utilizando ANSYS en la fase diseo, se logra experimentar diferentes hiptesis,

con fin de evaluar, cual planteamiento obtiene la adaptacin adecuada para el
problema especifico presentado.
La geometra de pila con aristas que forman ngulos rectos (rectangular)

obstaculiza, en mayor grado al flujo de agua, consecuente con esto, la generacin
de corrientes descendentes, estableciendo presiones dinmicas altas en el lecho
del rio aguas arriba de la pila.
Se verifico la transicin de rgimen rpido a lento en el entorno del obstculo

(pilas), debido a que, las pilas funcionan como seccin de control.
Se confrontaron los resultados provenientes de la simulacin en cada geometra

(rectangular, circular y punta de diamante), concluyendo que la pila en punta de
diamante es la seccin que altero en menor magnitud el flujo del agua aguas
arriba manteniendo un rango de velocidades de [6 9 m/seg], la pila circular
registro la velocidad local menor en el entorno de la pila.
La presencia de la forma hidrodinmica debilita la accin de los vrtices de

herradura, consecuente al flujo que es capaz de seguir la trayectoria del cuerpo
disminuyendo la divergencia de las lneas de corrientes. Las corrientes
secundarias generadas en el flujo separado, solo pueden actuar sobre el lecho
mediante vrtices laterales ya bien distanciados de la cara posterior de la pila.
El porcentaje de reduccin en la profundidad de socavacin local utilizando la

seccin hidrodinmica P2 fue de 30[%] en comparacin con los desarrollados por
las pilas comunes. Adems aumentando las dimensiones para la seccin P2, la
funcionalidad contra la socavacin local se hace ms efectiva.
133
10. RECOMENDACIONES
Realizar un proyecto de investigacin con un modelo a escala de tipo experimental

para corroborar los resultados de la modelacin numrica entregada en el
proyecto. Estimando la socavacin para las geometras comunes
y
confrontndola con la provocada en la geometra P2.
Implementacin de la versatilidad de ANSYS para la comprensin de teoras

estudiadas en las diversas ingenieras que son poco factibles de visualizar en las
aulas de clase y en los laboratorios.
Modelar las pilas de los puentes de San Jos de Ccuta, enfatizando la

interaccin del suelo con el fluido, aplicando el paquete computacional de ansys
FSI, a fin de visualizar la tasa de acarreo, arrastre de material del lecho de rio y la
formacin del cuenco de socavacin. De la misma manera modelar las pilas
comunes despus de ser carenadas con la geometra P2 propuesta en este
proyecto.
Contactar expertos en el manejo del software para profundizar en esta

herramienta en la obtencin y exportacin de datos para su respectivo anlisis.
134
BIBLIOGRAFA
CASTELLET, Ernest Blad. Modelacin del flujo en lmina libre sobre cauces
naturales. Barcelona: Universidad Politcnica de Catalunya, 2005.
CENGEL, Yunus. CIMBALA, John. Mecnica de fluidos. Bogot: Mc Graw Hill
Interamericana, 2006.
CHOW, Ven Te. Hidrulica de canales abiertos. Bogot: McGraw Hill
interamericana, 2004.
KAFI, Mark. ALAM, John. Modification of local scour equations, journal of the
institution of engineers. India: ASCE, 1995.
MOTT, Robert. Mecnica de fluidos. Bogot: Pearson, 2006.
PLATA, Francisco. Revisin de ecuaciones que predicen la socavacin local

alrededor de pilas de puente. Bogot: Universidad de los Andes, 1997.
POTTER, Merle. Mecnica de fluidos. Bogot: Cengage Learning, 2002.
ROCA, Alex. Socavacin en pilas de puente, diseo de estructuras protectoras

con el fin de disminuir la socavacin local y el riesgo de colapso. Bogot:
Universidad de los Andes, 1997.
RODRGUEZ CALDERN, Wilson. PALLARES MUOZ, Myriam Roco.

Modelacin del flujo de aire alrededor de automviles para la determinacin de las
fuerzas lift y drag por medio de elementos finitos. Bogot: Pontificia Universidad
Javeriana, 2005.
U.S DEPARTMENT OF TRANSPORTATION. Evaluating scour at bridges.

Washington: National Highway Institute, 2002.
135
ZAMBRANO, Edwar. Anlisis geomecnico de un pozo de hidrocarburos

empleando elementos finitos. Ccuta: Universidad Francisco de Paula Santander,
2008.
136
ANEXOS
137
Anexo A. Lista de comandos para modelar la seccin p2 utilizando

ANSYS 10.
Los comandos presentados a continuacin permiten modelar la seccin P2 en 3D

con el modulo FLOTRAN en ANSYS. Empleando esta cadena se podr dar forma
hidrodinmica a cualquier pila comn (Rectangular, circular y Punta de diamante).
La lista de comandos simula el flujo en el entorno de la pila a temperatura

estndar de 20 C, si se quiere cambiar la temperatura, se tendr que modificar
las propiedades del agua (Densidad, viscosidad), se discretizar el dominio con
FLUID142. El resultado que se observar con la lista de comandos ser la
solucin nodal de la velocidad en la pila P2.
Las variables a emplear son:
Ba = Base (vista en planta)
He= Altura (vista en planta)
Ti= Tirante de Agua (Profundidad del flujo)
V= Velocidad superficial del cauce
Lista de comandos
FILNAME, P2
/TITLE, Modelacin tridimensional del flujo de agua en la pila
/UNIT,SI
/ PREP7
ET,1, FLUID 142
MP,DENS,1,998.2
MP,VISC,1,1.005E-03
Ba=0.8
He=3
ti=6
v=9
/VIEW,1,1,1,1
/VUP,1,Y
/REPLOT
! Se define los keypoints, Open Keypoints.xls
K,
1,
0,
0,
0
138
K,
2,
0,
0,
2+Ba+2
K,
3,
3+He+5, -0.077, 0
K,
4,
3+He+5, -0.077, 2+Ba+2
K,
5,
3,
-0.0546, 2
K,
6,
3,
-0.0546, 2+Ba
K,
7,
3+He, -0.0546, 2
K,
8,
3+He, -0.0546, 2+Ba
K,
9,
1.2,
-0.0546, 2+Ba/2
K,
10,
3+He/2, -0.0546, 2+Ba+0.2
K,
11,
3+He/2, -0.0546, 1.8
K,
12,
3+He+1.80, -0.0546, 2+Ba/2
BSPLIN,9,10,12,,,,,,,,, ' genera la primer lnea que forma el permetro del perfil 001
BSPLIN,9,11,12,,,,,,,,, ' genera la segunda lnea que forma el permetro del perfil 001
L,1,2
L,2,4
L,3,4
L,1,3
LPLTO
/PNUM,line,1
LSEL,all
LPLOT,all
/PNUM,line,0
AL,1,2
AL,3,4,5,6
/PNUM,area,1 ' activar numeracin reas
ASEl,all
APLOT,all
/PNUM,area,0 'desactivar numeracin reas
VOFFST,1,7, ' generacin de volumen por extruccin
VOFFST,2,Ti,
VSBV, 2,1,,, 'extraer del volumen 2 el 1
/PNUM,line,1
LSEL,all
LPLOT,all
*DO, T, 19,24, 1 ' bucke para discretizar lneas frontera solida
LESIZE,T,0.08,,,,,,,
*ENDDO
LESIZE,19,,76,,,,,,
LESIZE,21,,76,,,,,,
*DO, T, 3,6, 1 ' bucke para discretizar lneas del equilibrio
LESIZE,T,0.3,,,,,,,
*ENDDO
*DO, T, 11,18, 1 ' bucke para discretizar lneas del equilibrio
LESIZE,T,0.3,,,,,,,
*ENDDO
/PNUM,line,0
/PNUM,volu,1
139
VPLOT
VSWEEP, 3
/PNUM,area,1 ' activar numeracin areas
ASEL,all
APLOT,all
*DO, T, 7,8, 1
DA, T, VX, v, , , , , , , ,
*ENDDO
DA, 10, VX, v, , , , , , , ,
DA, 9, PRES, 0
*DO, Q, 11, 12, 1
DA, Q, VX, 0, , , , , , , ,
*ENDDO
*DO, t, 11, 12, 1
DA, t, VY, 0, , , , , , , ,
*ENDDO
FINISH
/SOLU
ACEL, ,-9.81,
FLDATA1, SOLU, TURB, TRUE
SOLVE
/POST1
PRNSOL, VSUM, COMP Velocidad del flujo
140
Anexo B. Sitios de inters para el aprendizaje de ANSYS 10.
A continuacin se presentan sitios de inters en internet donde se puede encontrar

tutorial para la enseanza de ANSYS
http://courses.cit.cornell.edu/ansys/aivan1
http://www.mece.ualberta.ca/tutorials/ansys/IT/Distributed/Distributed.html
http://www.usouthal.edu/engineering/mechanical/faculty/phan/ANSYS_SBeam.pdf
http://www.andrew.cmu.edu/course/24-ansys/problems.htm
http://www.aq.upm.es/Departamentos/Estructuras/e96-622/HTML0102/HTML.html
http://www.engr.uky.edu/~jrbake01/ansys-statics-1-v8p1.pdf
http://ocw.kfupm.edu.sa/user062/AE42801/ANSYS_Tutorials_mk3_Very_Good.pdf
http://www.aaue.dk/bm/ansys/ANSYS-Examples.html
http://www.truveo.com/%E8%BD%B4%E5%AF%B9%E7%A7%B0%E7%BB%93%
E6%9E%84%E9%9D%99%E5%8A%9B%E5%88%86%E6%9E%90/id/356549800
7
141

Modelo Tridimensional Del Flujo en Pilas de Puente 01 PDF

Cargado por

Copyright:

Formatos disponibles

Modelo Tridimensional Del Flujo en Pilas de Puente 01 PDF

Cargado por

Información del documento

Título 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

Modelo Tridimensional Del Flujo en Pilas de Puente 01 PDF

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD FRANCISCO DE PAULA SANTANDER

BIBLIOTECA EDUARDO COTE LAMUS

AUTORES: IVAN HERNANDO RAMIREZ MENDOZA

PLAN DE ESTUDIOS: INGENIERA CIVIL

DIRECTOR: MEIMER PEARANDA CARRILLO

TITULO DE LA TESIS: MODELO TRIDIMENSIONAL DEL FLUJO DE AGUA EN PILAS

El propsito del siguiente trabajo es modelar el flujo de agua en los tipos

MODELO TRIDIMENSIONAL DEL FLUJO DE AGUA EN PILAS DE PUENTES

IVAN HERNANDO RAMIREZ MENDOZA

UNIVERSIDAD FRANCISCO DE PAULA SANTANDER

MODELO TRIDIMENSIONAL DEL FLUJO DE AGUA EN PILAS DE PUENTES

IVAN HERNANDO RAMIREZ MENDOZA

Trabajo de grado presentado como requisito para optar al ttulo de

UNIVERSIDAD FRANCISCO DE PAULA SANTANDER

A mi madre Nhora Esmeralda Mendoza Sanchez por sus oraciones y

A Mariluz Rosado Mendoza porque nunca me ha fallado como madre y padre a la

Los autores del trabajo expresan sus agradecimientos a:

Al semillero de investigacin ACUA, en el cual surgi la idea de centrar nuestra

Al M.Sc. Gustavo Adolfo Carrillo por enfocarnos en el tema especifico de nuestro

Al Ingeniero Mecnico Meimer Pearanda Carrillo por sus sugerencias brindadas

Al PhD Carlos Acevedo por asignarnos el espacio en la sala del grupo

2.1 FLUJO TRIDIMENSIONAL

2.2 PROPIEDADES DEL FLUJO

2.3 PROPIEDADES DEL FLUIDO

2.3.2 Viscosidad dinmica

2.4 TIPOS DE FLUJO

2.4.1 Flujo libre

2.4.2 Flujo variado permanente

2.5 RGIMEN DEL FLUJO

2.5.1 Flujo con rgimen turbulento

2.5.2 Flujo incompresible

2.6 FUNDAMENTOS DE VISUALIZACIN DEL FLUJO

2.6.1 Punto de estancamiento

2.6.3 Lneas de corriente

2.6.4 Graficas de contornos

2.7 VORTICIDAD Y ROTACIONALIDAD

2.8 FUERZA DE ARRASTRE (DRAG)

3. PROCEDIMIENTO GENERAL PARA SOLUCIN DE

3.1 GENERACIN DE LA MALLA ESTRUCTURADA

3.2 CONDICIONES DE FRONTERA

3.3 CONDICIONES DE FRONTERA EN LA SUPERFICIE SLIDA (PILA)

3.4 CONDICIONES DE FRONTERA DE FLUJO DE

4.2 SOCAVACIN LOCAL EN PILAS

5. PLANTEAMIENTO DEL PROYECTO

5.1 DOMINIO DEL PROBLEMA

6. METODOLOGA PARA LA APLICACIN DEL SOFTWARE

6.1 PILA RECTANGULAR MACIZA

6.1.1 Construccin de la geometra de la pila

6.1.3 Condiciones de frontera del flujo de entrada

6.1.4 Condicin de frontera del flujo de salida

6.1.5 Condicin de frontera en la superficie solida (pila)

6.2 PILA CIRCULAR APORTICADA

6.2.1 Construccin de la geometra de la pila

6.2.3 Condiciones de frontera del flujo de entrada

6.2.4 Condicin de frontera del flujo de salida

6.2.5 Condicin de frontera en la superficie solida (pila)

6.3 PILA EN PUNTA DE DIAMANTE

6.3.1 Construccin de la geometra de la pila

6.3.3 Condiciones de frontera del flujo de entrada