Mau Al Redes 2007
Mau Al Redes 2007
Mau Al Redes 2007
FACULTAD DE INGENIERÍA
CENTRO DE INVESTIGACIONES EN
ACUEDUCTOS Y ALCANTARILLADOS
CIACUA
Versión Académica
También calcula la edad del agua en la red y la concentración de un soluto (cloro residual)
en los nodos para condiciones de flujo permanente.
La Interfaz Gráfica del programa ofrece numerosas opciones para visualizar las variables de
los nodos y los tubos. Las propiedades de los nodos se pueden ver en 3D, así como la
topografía de la red. Tiene una gran gama de colores que permite apreciar gráficamente las
diferencias en los valores de las variables en los nodos, tubos y embalses. La visualización
en 3D esta mejorada, de manera que la representación de las variables de los nodos se
encuentra ajustada a la forma de la red, sin interpolaciones a formas rectangulares.
Esta versión mejora la representación gráfica de los tubos con la inclusión de los vértices
como accesorios que permite ver los diferentes quiebres que tienen los tubos debido a la
topografía en que se encuentra la red.
La configuración de los valores predeterminados por el programa pueden ser cambiados por
el usuario, para adaptarse a sus necesidades, las cuales pueden ser el valor de los diámetros
comerciales así mismo como la rugosidad para cierto material, que no necesariamente
deben corresponder con los valores que se encuentran en las listas que trae el programa. Las
propiedades físicas como gravedad y temperatura pueden ser modificadas para una mejor
modelación del problema o situación.
Copie la carpeta del programa redes en la ruta deseada. Esta debe contener el ejecutable de
redes (redes2007academica.exe) y tres librerías dll adicionales. Para ejecutar el programa
haga doble clic en el archivo redes2007academica.exe.
Durante el resto de este capítulo se va a desarrollar este ejemplo con el fin de mostrar la
funcionalidad del programa y su forma de uso en un ejemplo concreto.
Demandas en los
Nodos
ID Qdem(l/s)
2 60
3 40
4 30
5 30
6 40
Características de la tuberías
ID Union1 Union2 Diámetro Longitud Rugosidad CoefMenores
(Pulgadas) (m) (mm)
1 1 2 10 500 0.06 0
2 2 3 6 400 0.06 10
3 3 4 4 200 0.06 0
4 5 4 6 400 0.06 0
5 2 5 4 200 0.06 0
6 6 5 8 600 0.06 0
7 1 6 10 300 0.06 0
Se repite el procedimiento hasta completar la adición de todos los tubos. Figura 2-7
Se hace “clic” sobre el tubo que se desee insertar la válvula. En este caso se realizará sobre
el tubo 2. Figura 2-14
Figura 2-14 Adición de una válvula sobre un tubo existente.
Para la adición de un nuevo tubo con una válvula, se debe hacer la primera parte del
procedimiento anterior, es decir, hacer “clic” sobre el botón y seleccionar el tipo de
válvula en la ventana Insertar Válvula. Luego se procede a insertar el tubo con válvula de
manera igual que se hizo para adicionar un tubo.
Figura 2-15 Adición de una válvula en un tubo nuevo.
La adición de bombas al igual que las válvulas, se hace mediante tubos que contienen una
bomba con características específicas, las cuales se explicarán más adelante. Hay dos
formas de llevar a cabo la inserción de una bomba, sobre un tubo ya existente ó con un tubo
nuevo en la red. En ambos casos se debe hacer “clic” sobre para colocar modo
inserción de bombas.
2.4.5.1 Adición de una bomba sobre un tubo existente
En modo inserción de bombas se hace “clic” sobre el tubo que se desea insertar la bomba.
En este caso se eligió el tubo 1. Figura 2-16
Para cargar el editor de red basta con hacer un “clic” sobre el botón , con lo cual se
despliega el Editor de Red. Figura 2-18
Las unidades también pueden modificarse; siguiendo con el ejemplo, se puede editar la
columna de diámetro, inicialmente aparece el diámetro en milímetros, pero la información
que se dispone de los diámetros está en pulgadas, esto NO representa ningún problema, ya
que el programa REDES 2007 cuenta con un poderoso y flexible sistema de unidades que
permite ingresar los datos en las unidades más comunes según la variable correspondiente.
Figura 2-20
Ahora se puede cerrar el programa o presionar el botón para crear una red nueva
sobrescribiendo la anterior, por eso es recomendable guardar el modelo actual antes de usar
este botón.
Para abrir y cargar un archivo de redes (.red) se hace “clic” sobre el botón después de
lo cual se despliega la ventana de la Figura 2-27.
Esto carga un modelo, con el cual se puede continuar el trabajo donde se dejo la última vez
que se guardo el archivo.
Para realizar el cálculo hidráulico bajo la suposición de flujo permanente se debe presionar
el botón , luego se lo cual se despliega el mensaje de CALCULO EXITOSO. Figura
2-29
Para ver los resultados en el mapa hay que hacer “clic” en el botón , al hacer esto se
despliega el Editor de Visualización. Figura 2-31
En esta ventana se puede seleccionar la variable que se desea visualizar en el mapa para
cada tipo de elemento, por ejemplo demanda en el embalse, presión en los nodos y caudal
en los tubos. Figura 2-32
Figura 2-32 Editor de Visualización, selección de variables
Al hacer clic en “Aceptar” se despliega la ventana principal del programa mostrando los
valores de las variables seleccionadas. Figura 2-33
Barra de Herramientas
La Barra del Menús esta localizada en la parte superior de la ventana principal y es una vía
fácil para introducir los datos y correr el modelo, esta contiene las siguientes opciones:
• Archivo
• Editar
• Ver
• Calcular
• Resultados
• Opciones
• Unidades
• Ayuda
3.2.1 Archivo
3.2.2 Ver
El menú Ver, ofrece 3 opciones diferentes: Nodos, Embalses y Tubos. En cada una de estas
opciones aparece un recuadro dividido en dos. La parte izquierda muestra las opciones de
entrada al modelo y la derecha las opciones de salida. Al seleccionar cualquiera de estas se
ven en el modelo los datos deseados. Se debe tener en cuenta que las opciones de la derecha
solo funciona siempre y cuando el modelo ya esté calculado.
3.2.3 Calcular
Calidad de Agua Estado Estable: Calcula las condiciones de calidad del parámetro
introducido en un instante del tiempo.
3.2.4 Resultados
Tabla
En esta tabla ( Figura 3-2 ) se pueden observar los diferentes resultados, ya sea de Nodos o
Tuberías. La tabla muestra dos opciones dependiendo del tipo de cálculo, si es de Estado
Estable o Periodo Extendido. En el primer caso, el usuario tiene la posibilidad de
seleccionar si es Nodo o Tubo y después escoger la variable que quiere estudiar. La tabla
muestra el ID de cada elemento y su respectivo valor. Las unidades en las que se muestran
los valores pueden variarse sobre la misma tabla.
Tubo:
• Caudal
• Velocidad
• Volumen
• Factor de Fricción
• Perdidas Totales
• Perdidas por Fricción
• Perdidas Menores
• Reynolds
Para el caso de un cálculo en Periodo Extendido, la tabla ofrece la opción de Hora Fija o
Serie de Tiempo. En Hora Fija, se ven los resultados de toda la red (Nodos o Tubos) en una
hora dada, la cual puede variarse para poder abarcar la totalidad del periodo. Si se
selecciona Serie de Tiempo, se observan los resultados de un elemento en especial durante
todo el periodo de cálculo. El elemento del cual se están mostrando se puede cambiar en la
casilla de selección de la derecha.
3.2.5 Opciones
En este menú se pueden seleccionar algunas opciones las cuales se listan a continuación:
El botón “Aceptar” incorpora los cambios al modelo, “Cancelar” se sale de la ventana sin
tener en cuenta los valores modificados y “Reset” sirve para volver a cargar los valores
predeterminados, estos valores son los que se encuentran en la literatura para el cloro
residual.
Los siguientes botones forman un grupo, uno de ellos debe estar presionado, según lo cual
Muestra una ventana con vista de red en 3-D y las propiedades en los nodos.
Este botón sirve para ajustar el mapa conservando la proporción entre las
coordenadas X y Y del modelo de la red.
Este botón sirve para cambiar las coordenadas máximas y mínimas del
modelo, se debe usar más que todo cuando se va a empezar un modelo desde
cero o cuando se va a ampliar un modelo ya construido.
Este botón sirve para disminuir el número de decimales con los que se
presentan los resultados en el mapa y en las tablas.
Este botón sirve para aumentar el número de decimales con los que se
presentan los resultados en el mapa y en las tablas.
Muestra el costo de la red que se calcula con una ecuación de costos con
parámetros establecidos.
El programa REDES 2007 cuenta con un poderoso motor gráfico que permite desplegar el
mapa de la red de tuberías. Figura 3-6
Figura 3-6 Red ejemplo, mapa total
Para desplegar el editor de visualización solo hay que presionar el botón para
mostrarlo. Figura 3-8
En esta ventana se pueden ver tanto las variables de entrada, como las variables calculadas.
Luego de hacer “clic” en el botón “Aceptar” se muestra el mapa de la red mostrando las
variables seleccionadas en sus unidades actuales. Figura 3-9
Figura 3-9 Ventana Principal, visualización de LGH y Demanda
Haciendo clic sobre la grafica y manteniéndolo presionado se puede girar y rotar la gráfica
para ver la topografía de la red de diferentes perspectivas. Para ver u ocultar los ejes de la
gráfica se hace clic sobre el botón
3.9.1 Superficie
Esta es una de las características de la Interfaz Gráfica que le permite al usuario tener una
mayor percepción a través de la visualización de las propiedades en los nodos que se
pueden elegir en el menú desplegable.
Para ver gráficamente las propiedades seleccionamos una del menú anterior y le damos clic
4.1 Generalidades
El programa REDES 2004 hace uso de patrones de consumo para multiplicar las demandas
base en los nodos y así obtener una demanda diferente para cada periodo de análisis, cada
nodo puede tener un patrón diferente.
Un patrón es una serie de valores, sirve para modelar el cambio de la demanda a lo largo
del día, generalmente está centrado en el valor de 1, de tal forma que se conserve el caudal
promedio a lo largo del día.
1.4
1.2
1
Valor del patrón
0.8
0.6
0.4
Variable
0.2
Constante
0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
Hora
Para desplegar la ventana del editor de patrones hay que presionar el botón .
La representación gráfica del patrón se actualiza cada vez que se hace clic sobre el
botón .
Para asociar un patrón a todos los nodos se debe presionar el botón luego sale
el Mensaje de confirmación de la asignación.
Para seleccionar el patrón para cualquier nodo, simplemente hace clic en la celda
correspondiente a la fila del nodo deseado y la columna “Patrón”, por ejemplo para el nodo
2 sería. Figura 4-5
Este procedimiento se repite hasta obtener la configuración deseada para cada nodo.
Capítulo 5 Análisis de la red
Para correr el cálculo hidráulico estático hay que presionar el botón en la ventana principal
Luego sale la ventana de selección de periodos, esta ventana se puede cerrar y hacer visible
con el botón . Figura 3-10
La ventana de selección de Periodos sirve para cargar en la red actual el resultado del
cálculo en ese instante del tiempo.
5.2 Análisis de la calidad del agua
Para calcular la calidad del agua es necesario haber calculado la hidráulica previamente.
En el caso de que sea un cálculo de Calidad del Agua o de la Edad en la tubería, estos
valores serán mostrados sobre el mapa tal como el ejemplo anterior. Sin importar el caso,
cada valor que se muestra esta acompañado de la variable y las unidades en las que se
encuentra.
De igual forma los resultados pueden ser observados desde una tabla, esta puede ser
accedida desde la Barra de Menús, Resultados, Tabla o por medio de la Barra de
En el caso de un cálculo en periodo extendido, es posible seleccionar una Hora Fija y ver
los resultados de toda la red a esta hora, la cual puede cambiarse para recorrer todo el
periodo. También se puede elegir la opción de Serie de Tiempo y escoger el elemento que
se quiere analizar durante todo el periodo de análisis.
Capítulo 7 Algoritmos y Ecuaciones
7.1 Hidráulica
7.1.1 Ecuación de pérdidas por fricción
El programa REDES calcula las pérdidas por fricción por medio de las ecuaciones
físicamente basadas (Darcy-Weissbach en conjunto con la ecuación de Colebrook-White).
Las ecuaciones físicamente basadas se terminaron de desarrollar hace poco más de 50 años.
La historia del desarrollo de las ecuaciones físicamente basadas incluye estudios desde
Leonardo Da Vinci (1500) hasta Prandtl, Von Kármán y Nikuradse (1950). Alrededor de
450 años de investigaciones.
La ecuación hidráulica físicamente basada que gobierna las pérdidas de energía por unidad
de peso del flujo al interior de una tubería es la siguiente:
L v2
Hf = f
D 2G
En donde:
Ésta ecuación fue desarrollada por Weisbach en 1854 y por Darcy, de manera
independiente en 1858. Sin embargo, la siguiente parte del desarrollo consistía en calcular
el valor de f. El proceso de explicar el comportamiento de este factor de fricción de Darcy
probó ser muy complejo. Se necesitaron nuevas teorías acerca del la interacción del flujo
con paredes sólidas y acerca de la turbulencia y su efecto sobre las características del flujo
viscoso.
Tras plantear la nueva hipótesis de la longitud de mezcla, y probarla con éxito en los datos
experimentales, Prandtl y Von Kármán resolvieron el problema del valor de f para el flujo
turbulento hidráulicamente liso. Nikuradse lo hizo para el flujo turbulento hidráulicamente
rugoso. El valor de f para flujo laminar, se podía resolver utilizando otras ecuaciones
físicamente basadas (la ecuación de Hagen-Poiseuille, obtenida mucho antes).
1 ks 2.51
= −2Log 10 +
f 3.7D Re f
en donde:
vD
Re =
en donde:
Re : Número de Reynolds
v : Velocidad
D : Diámetro
ν : viscosidad cinemática
V2
Hm = Km
2g
En donde los términos Hf y Hm son las perdidas por fricción y las perdidas menores,
respectivamente.
I S
Qi = Qs + D
i =1 s =1
en donde:
I S
Qi = Qs + ∆A
i =1 s =1
en la cual:
EF = EI − L
Además, algunos métodos tienen en cuenta que los cambios de energía a través de un
circuito suman cero, por ejemplo:
Con la anterior información, la perdida de cabeza en cada tramo de una tubería que conecte
con dos nodos de la red es:
α 1Q (n −1)
1
0 0 ... 0
0 α 2 Q (n2 −1) 0 ... 0
[A11] = 0 0 α 3 Q (n3 −1) ... 0
. . . . 0
0 0 0 ... α NY Q (nNT −1)
Se puede observar que la ecuación de continuidad para todos los nodos es:
[A21][Q] = [q]
En donde:
[A21] = Matriz transpuesta de [A12]
[q] = Vector de consumo (demanda) o de entrada (oferta) en cada nodo de la red,
con dimensión NN x 1
En la red de tuberías se pueden considerar los embalses (tanques de cabeza fija) y los nodos
como uniones. Los tubos se pueden considerar como arcos. Para calcular ciertas variables
de la red se utilizó un algoritmo recursivo implementado iterativo, que en adelante se
conocerá como ARI. La implementación del modelo de la red de tuberías del programa
REDES 2004 incluye este método genérico de cálculo (en la clase TGrafo).
El algoritmo ARI permite realizar cálculos complejos que requieren del recorrido de todos
los nodos y/o tubos de la red de tuberías. Se destaca por ser un método que ha probado su
efectividad en procesos como el cálculo estático de calidad del agua en redes grandes
sustituyendo metodologías tradicionales que se demoran muchísimo más en la finalización
de los cálculos.
El modelo matemático ARI se caracteriza por ser robusto y genérico, permite su adaptación
al cálculo de varios tipos de variables según el interés particular, sin importar si la sustancia
a analizar es o no conservativa.
7.3.1 Descripción del algoritmo
n
(Ci ⋅ Qi )
Ck = i =1
n
Qi
i =1
En donde:
n es el número de tubos que alimentan al nodo k.
Ck es la concentración en el nodo k.
Ci es la concentración al final del i-ésimo tubo que alimenta al nodo k.
Qi es el caudal del i-ésimo tubo que alimenta al nodo k.
C = Co ⋅ e − K ⋅t
En donde:
C es la concentración aguas abajo del tubo.
C0 es la concentración agua arriba del tubo.
K es un coeficiente general de reacción del químico con el cuerpo de agua y las
paredes del tubo.
t es el tiempo promedio de viaje del agua en el tubo.
Es claro que:
Voli
ti =
abs (Qi )
En donde:
ti es el tiempo de viaje promedio del agua en el i-ésimo tubo.
Voli es el volumen del agua en el i-ésimo tubo.
Qi es el caudal del i-ésimo tubo.
υ
Sch =
Dif
El Algoritmo Genético también hace parte de las opciones de esta ventana, y se puede
La ventana se carga con valores predeterminados, los cuales funcionan muy bien en el
diseño en la mayoría de los casos.
- Tamaño de la Población: es el número de individuos para cada generación.
- Número de Generaciones: es el número de generaciones que va a tener el algoritmo.
- Probabilidad de Mutación:
- Costo Constructivo:
- Constante de Reproducción:
- Generación Inicial: para la generación inicial puede el programa generar un aleatoria en
el caso en que se vaya a optimizar por primera vez la red, o puede utilizar cono
generación inicial la red actual para intentar encontrar una mejor y optimizar el
algoritmo.
Selecciona los diámetros inactivos, y deja inactivos los diámetros que están
seleccionados.
Si se pretende hacer un diseño óptimo de la red y no se desea interactuar con las diferentes
opciones que ofrece el programa, se puede hacer uso de la opción Diseño Rápido haciendo
Se debe volver a la parte de Diseñar para realizar un diseño óptimo de la red basado en la
Superficie Óptima de presiones, después de haber seleccionado la ecuación para la Línea de
Gradiente Hidráulico y la Presión Mínima, y en el caso de la ecuación Cuadrática el
porcentaje para la flecha.
8.4 Costo
Esta opción que se encuentra en la ventana de Diseño Óptimo muestra las ecuaciones para
calcular el costo de una red. Figura 8-10
• x: es el exponente del diámetro debido a que su relación con el costo porque la relación
entre el área y el diámetro de un tubo no es lineal.
La segunda Ecuación es para calcular el costo de la red, la cual consiste en la suma de los
costos individuales de cada tubería.
También se puede ver el costo actual de la red, ya que aparece en la opción de Costo de la
ventana Diseño Óptimo.
Programación por Restricciones es otra de las opciones que contiene el Módulo de Diseño
del Programa Redes, el cual consta de dos partes:
• Aumentar Diámetros: esta opción se utiliza cuando las presiones en algunos nodos de la
red están por debajo de la presión mínima. Su algoritmo consiste en aumentar ciertos
diámetros de algunas tuberías para cumplir con la presión mínima en todos los nodos,
pero teniendo como restricciones que el costo de la red no aumente mucho y que la red
no vaya a quedar muy sobrada.
• Disminuir Diámetros: esta opción se utiliza cuando todas las presiones en los nodos de
la red cumplen con la presión mínima. Su algoritmo consiste en disminuir ciertos
diámetros de manera que se pueda reducir el costo pero cumpliendo con las
restricciones de la presión mínima en todos los nodos y con los caudales demandados
en cada nodo.
Para ver el funcionamiento de estos algoritmos se disminuirán los diámetros de la red con
la que se ha venido trabajando de manera que no cumpla con la presión mínima en todos
los nodos. Figura 8-11
Figura 8-11 Red con presiones por debajo de la Presión Mínima en algunos Nodos.
Miremos los costos de la red para rectificar que no suben mucho aumentando ciertos
diámetros para cumplir con la presión mínima en todos los nodos utilizando el la opción de
diseño Aumentar Diámetro por restricciones. Figura 8-12
El cambio en el costo no es muy grande por lo que ahora se tiene una red un poco más
costosa pero que cumple con la presión mínima en todos los nodos sin presentar presiones
muy superiores a la mínima.
Ya que se tiene una presión mayor en todos los nodos a la presión mínima, miremos si se
pueden disminuir los diámetros para reducir los costos conservando la restricción de la
presión. Para llevar a cabo esta operación se debe presionar el botón
Cerremos la ventana de Diseño Óptimo y volvamos a la ventana Principal para mirar los
resultados obtenidos por el Algoritmo de Disminuir Diámetros. Figura 8-17
Para obtener mejores diseños se puede utilizar el Diseño por Programación por
Restricciones junto con el Diseño Genético.
Capítulo 9 Equipo de Trabajo
Dirección:
Juan Guillermo Saldarriaga V, Ingeniero Civil
Profesor Titular, Facultad de Ingeniería
Departamento de Ingeniería Civil, Universidad de Los Andes
Diseño y desarrollo
German Villalba F., Ingeniero Civil, Msc. Ingeniería de Sistemas.
Daniel Salas U., Ingeniero Civil, Msc. Ingeniería de Sistemas.
Documentación
Aura Rodríguez Ruiz.