ANALISIS DE REDES CERRADAS

En orden cronolgico se presentarn los siguientes mtodos de anlisis y diseo de redes cerradas, sin que se trate de una lista exhaustiva. Mtodo de Hardy-Cross con correccin de caudales en los circuitos Mtodo de Hardy-Cross con correccin de alturas piezomtricas en los nodos Mtodo de Newton-Raphson Mtodo de la teora lineal Mtodo del gradiente hidrulico Con base en mtodos reportados por la literatura tcnica y la facilidad de uso, la historia de los mtodos de anlisis de redes de distribucin de agua potable ha sido dividida en tres perodos: El primero, aproximadamente desde 1850 hasta 1930. El segundo, caracterizado por el uso de ecuaciones de resistencia fluida y de modelos anlogos elctricos se extiende de 1930 a 1960. El tercero, caracterizado por formulaciones matriciales del problema de las redes de distribucin con el fin de hacer un uso intensivo de computadores digitales, va de 1960 hasta nuestros das.

Perodo Perodo I

Ao 1845 1892 1905

Inventor(es)/autor(es)

Mtodo / aplicaciones Frmula para la prdida de altura en un flujo a travs de una tubera simple. Solucin grfica. Frmula para la prdida de altura en el flujo a travs de una tubera simple y un mtodo de tubera equivalente.

Darcy y Weisbach
Freeman Hazen y Williams

Perodo II

1934
1936 1956 1957

Camp y Hazen Cross MCLLROY

Anlisis de una red elctrica.

Tcnica de relajacin. Anlisis del fluido Mcllroy. Adaptacin del mtodo de Hardy Cross para computadores digitales. Mtodo del nodo simultneo. Expansin del mtodo del nodo simultneo. Mtodo del circuito simultneo. Programa comercial para el anlisis de redes con base en el mtodo del circuito simultneo. Teora lineal. - KYPIPE, Programa comercial para el anlisis de

HOAG y Weinberg
Martin y Peters Shamir y Howard EPP y Fowler Jeppson

Perodo III

1963 1968 1970 1977

1972

Wood y Charles

redes.

1980
1987 1994

Wood
Todini y Pilati Rossman

Mtodo del gradiente.

- EPANET Programa comercial para el anlisis de redes.

Principios fundamentales de anlisis de redes cerradas

Si se considera la red cerrada que aparece en la figura 7.1 y se tiene en cuenta que QD1 QD2 QD3 QD4 QDNU son los caudales consumidos en cada uno de los nodos, algunos de los cuales podran tener un valor nulo en un momento dado, y que Qe1 Qe2 Qe3 Qem son los caudales que alimentan la red de distribucin, se puede establecer la siguiente ecuacin de conservacin de la masa:
m
Nu

Q
i 1

e
i

Q (7.1)
D
1

donde Nu es el nmero de uniones (nodos) existentes en la red.

Figura 7.1 Red cerrada. Caudales consumidos en los nodos y caudales de alimentacin a la red.

La Ecuacin 7.1 es una ecuacin de conservacin de la masa para la red como un todo. Sin embargo, para cada uno de los nodos de la red se puede establecer una ecuacin similar, debido a que localmente tambin se debe cumplir el hecho de que la masa se conserve. Dicha ecuacin es:

NTi

Q
j 1

ij

Di

0 (7.2)

donde NTi es el nmero de tubos que llegan al nodo i, y Q representa el caudal que pasa por la tubera Qij hacia el nodo i desde el nodo j. La convencin adoptada por la prctica de la ingeniera hidrulica es que puede ser positivo (va hacia el nodo i) o negativo (sale de dicho nodo). Para cada uno de los caudales Qij de la Ecuacin 7.2 se puede plantear la siguiente ecuacin de conservacin de la energa entre los nodos i y j, incluyendo las prdidas por friccin y las prdidas menores, en trminos de las alturas piezomtricas en dichos nodos:

En donde se ha utilizado la ecuacin de Darcy-Weisbach para el clculo de las prdidas por friccin. Si se despeja Qij de esta ltima ecuacin se obtiene la siguiente expresin, la cual relaciona el caudal que pasa por la tubera ij con las alturas piezomtricas en los nodos i y j.

Si se reemplaza este ltimo resultado en la Ecuacin 7.2, se obtiene:

donde NTi representa el nmero de tuberas que llegan a la unin (nodo) i. A fin de tener en cuenta en forma automtica el signo del caudal ij, la Ecuacin 7.3 se puede cambiar por la siguiente expresin:

Este tipo de ecuaciones para el diseo y anlisis de redes cerradas de tuberas se conocen como ecuacin de altura piezomtrica. En caso contrario se puede suponer alguna de las alturas piezomtricas, ya que los valores absolutos de stas no afectan la distribucin de caudales y adems debe tenerse en cuenta que las ecuaciones de altura piezomtrica son ecuaciones no lineales.

Por otro lado, a partir de los circuitos de tubos que conforman la red, los cuales pueden ser adyacente o superpuestos, se pueden plantear las siguientes ecuaciones de conservacin de la energa, una para cada uno de los circuitos que conforman la red de distribucin: 1. Ecuacin de continuidad en las uniones que conforman el circuito:

2. Ecuacin de conservacin de la energa alrededor del circuito:

Donde NTi, es el nmero de tubos del circuito I, luego, si se utiliza la ecuacin de Darcy Weisbach en esta ultima ecuacin en conjunto con la expresin de las prdidas menores como funcin de altura de velocidad, se obtiene la siguiente expresin:

Las ecuaciones de tipo 7.7 se conocen como ecuaciones de caudal de la red. Las ecuaciones de tipo 7.7 se conocen como ecuaciones de caudal de la red. En total se tienen NC ecuaciones de caudal, donde NC es el nmero de circuitos que conforman la red. Nuevamente se puede observar que son ecuaciones no lineales. Lo anterior implica que para el clculo de una red cerrada se tiene un nmero total de ecuaciones igual a:

Una vez ms debe establecerse una convencin de signos para las ecuaciones de caudal 7.7. Los caudales en el circuito se consideran positivos si giran en el sentido de las agujas del reloj y negativos si lo hacen en sentido contrario. Para asegurar una correcta asignacin del signo, estas ecuaciones se pueden transformar como se indica a continuacin:

Los mtodos de anlisis de redes que se describen en este captulo estn diseados para llevar a cabo esos clculos de caudales en cada tubera piezomtrica en cada nodo, esto implica que se deben conocer todas las dems variables relacionadas con las tuberas (dimetricos, rugosidades, coeficientes de prdidas menores, accesorios especiales y bombas) y con los nodos (caudales de consumo, altura topogrfica, tanques, diseo y no un diseo en s

MTODO DE HARDY-CROSS CON CORRECCIN DE CAUDALES

Este mtodo para resolver las Ecuaciones 7,4 y 7.7 fue desarrollado en 1936 por el ingeniero norteamericano Hardy Cross, quien era profesor de ingeniera estructural en la Universidad de Illinois, Estados Unidos.

El mtodo en su forma original fue desarrollado para el clculo de estructuras aporticadas de concreto y acero, mediante un mtodo matemtico para llevar a cabo anlisis de distribucin de momentos para estructuras estticamente determinadas. Sin embargo, Cross lo extendi rpidamente al caso de redes cerradas de distribucin de agua potable, y public un artculo en el cual describa la aplicacin de su mtodo (ver Cross, 1936). Como ecuacin de resistencia fluida utiliz la ecuacin de Hazen Williams, establecida en 1906 y popularizada en la dcada de 1930. Para esta poca el material predominante en las redes de distribucin de agua potable era el hierro fundido.

Si se define una altura piezomtrica que incluya la altura piezomtrica perdida por friccin y la altura piezomtrica perdida por accesorios, en la siguiente forma:

la anterior ecuacin se convierte en:

Ahora, utilizando la Ecuacin 7.8 se tiene que:

Esta ltima ecuacin tambin puede escribirse en la siguiente forma:

Mtodo de Hardy Cross con correccin de caudales: pasos que se debe seguir en el anlisis.
El anlisis de una red de distribucin de agua segn el mtodo de Hardy-Cross con correccin de caudales en los circuitos propone los siguientes pasos:

1. Se define claramente la geometra de la red, identificando en forma coherente los nodos y los circuitos. 2. Si existe ms de un nodo con altura piezomtrica constante (tanque en la red o embalse), es necesario conectarlos en pares por medio de tuberas hipotticas que pueden representarse mediante lneas punteadas. En estas tuberas hipotticas se deben suponer dimetros, longitudes y rugosidades absolutas, de tal manera que se pueda calcular el caudal correspondiente a las diferencias de nivel entre los diferentes pares de embalses o tanques. En las correcciones de caudales no deben incluirse los tubos hipotticos, lo cual s debe hacerse en el clculo de las prdidas de altura piezomtrica (por friccin y por accesorios). 3. Se suponen todos los dimetros de la tubera que conforman la red. Tal paso convierte este mtodo en un proceso de comprobacin de diseo. 4. Se supone que la red est compuesta por circuitos cerrados en cualquier orden. Con el fin de acelerar la convergencia se puede suponer que los tubos de dimetros grandes forman circuitos independientes. Se deben utilizar tantos circuitos como sea necesario para asegurar que todos los tubos queden incluidos en por lo menos un circuito.

5. Se supone el caudal a partir de cualquiera de las tuberas de la red. Luego se procede alrededor del circuito que contiene esta tubera para calcular los caudales en las dems tuberas que conforman el circuito teniendo en cuenta los caudales que salen de las uniones (caudales negativos) y los que entran a ellas (caudales positivos). Si los flujos hacia o desde otro circuito son desconocidos, se deben suponer los caudales correspondientes. Esto significa que se deben hacer tantas suposiciones de caudales como circuitos existan en la red que se est analizando. Cuanto mejores sean estas suposiciones ms rpidamente converger el mtodo. La experiencia ayuda mucho en este aspecto. 6. Se calcula la prdida de altura piezomtrica en cada tubera de la red utilizando la siguiente ecuacin (de Darcy-Weisbach). si bien podra emplearse cualquier ecuacin de resistencia fluida, como la de Hazen-Williams.

El factor de friccin f se calcula utilizando la ecuacin de Colebrook-White presentada en el Captulo 1:

1 f

2 log

10

s 3.7d

2.51 f
e

Junto con el Diagrama de Flujo 2. 7. Se calcula la prdida neta de altura piezomtrica alrededor del circuito, es decir, se suman las prdidas altura piezomtrica y se restan las adiciones de altura piezomtrica siempre medidas en el sentido del avance de las agujas del reloj. Si la prdida neta de altura piezomtrica no es cero, se procede a corregir los caudales de cada una de las tuberas del circuito mediante la Ecuacin 7.10:

9. Los pasos 5 a 8 se repiten para todos los circuitos teniendo en cuenta los caudales corregidos en los circuitos calculados previamente. 10. Los pasos 5 a 9 se repiten hasta que el balance de alturas piezomtricas alrededor de todos los circuitos (ecuacin de conservacin de la energa) llegue a valores razonablemente cercanos a cero. Este criterio de convergencia es fijado por el diseador de acuerdo con las caractersticas de la red que est analizando.

EJEMPLO N 01
La red mostrada en la siguiente figura tiene una vlvula en la tubera 2-3, a cual se encuentra parcialmente cerrada y produce una prdida menor local de 10.0 v232 / 2g, la presin en el punto 1 es 100 mca. Analizar los causales y presiones en la red. Los dimetros (en milmetros) y las longitudes (en metros) para cada una de las tuberas son los indicados. Los caudales estn dados en l/s (v=1.141x10-6 m2/s), suponer que las perdidas menores son despreciables, salvo en la tuberia 2-3

Ejemplo del mtodo de Hardy - Cross con correccin de caudales

SOLUCIN
La primera suposicin de caudales en las tuberas puede ser:

Estos datos se pueden representar en forma grfica, tal como se muestra en la Figura, en donde queda claro el por qu de los signos de los caudales de la tabla anterior:

Primera suposicin de caudales en las tuberas

Mediante las Ecuaciones 7.11 y Colebrook-White, se calculan las siguientes tablas, en las cuales se ha supuesto que la viscosidad cinemtica del agua n es 1.14x106 m2/s.

Primer ciclo

Mediante la Ecuacin 7.10 se calcula la correccin de caudal:

Al utilizar la Ecuacin 7.10 se calcula la correccin de caudal:

Segundo ciclo

Mediante la ecuacin 7,10 se calcula la correccin de caudal:

La Ecuacin 7.10 lleva al clculo de la correccin de caudal:

Tercer ciclo

Al utilizar la Ecuacin 7.10 se calcula la correccin de caudal:

Mediante la Ecuacin 7.10 se calcula la correccin de caudal

Cuarto ciclo

La Ecuacin 7.10, lleva al clculo de la correccin de caudal:

Al utiliza la Ecuacin 7.10 se calcula la correccin de caudal:

Como las sumatorias de las perdidas por friccin y las perdidas menores son pequeas en ambos circuitos, el proceso debe parar. Esto significa que se esta cumpliendo la ecuacin de conservacin de la energa en cada circuito. En el momento de parar el proceso las correcciones de caudales (Q) resultantes eran muy pequeas. Los resultados que se obtuvieron son:

Resultados de caudales en las tuberas

Estas alturas piezomtricas de presin pueden variar su valor, dependiendo de la direccin en que se calculen. Tal variacin se debe a que en el momento de parar el proceso iterativo no se tiene que precisin absoluta en los caudales. Sin embargo, las diferencias son pequeas. Este hecho es vlido para los ejemplos de todo este captulo

COMENTARIO
Como se haba anteriormente, el mtodo de Hardy-Cross utilizado en este ejemplo requiri una suposicin inicial de caudales con el fin de iniciar el proceso de clculo. Al ser una red tan pequea las suposiciones de caudal son fciles y rpidas de hacer, con lo cual no se pierde mucho tiempo en la preparacin de los datos. Una vez se tienen estas suposiciones, el proceso necesit nicamente cuatro iteraciones para converger.