U.N.S.C.

H
EscuelaProfesionaldeIngenieraCivil
DiseodeunaReddeAbastecimientodeaguaMtododelGradienteHidrulico Curso :AbastecimientodeAguaPotableyAlcantarillado Profesor :Ing.JoelOrIwanaga Estudiante :CANCHARIGUTIRREZ,Edmundo. Cod.Est. :16005011

Ejemplo #03:
Obtener caudales y presiones en cada tubera y nudos de la red respectivamente, elegir dimetro de las tuberas adecuadas. El sistema de agua potable mostrado est abastecido por dos reservorios como se muestra.

Solucin:
El primer paso es dividir el sistema en una serie de elementos finitos identificando sus puntos extremos como nudos, una tubera debe estar plenamente identificada en la red por su nudo inicial y final estableciendo implcitamente la direccin del flujo del caudal en la tubera. Se debe enumerar nudos y tubera como se muestra.

Donde: Nmero de tuberas Numeracin de nudos Direccin flujo de caudal.

Comentaros: [email protected] Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil 1 de 27

Abastecimiento de Agua Potable Y Alcantarillado Anlisis y Diseo de Redes de Agua Potable Mtodo del Gradiente Hidrulico
Mtodo del Gradiente Hidrulico - Argumentos

2.0 Argumentos 2.1 Definiendo la Red (RED) Cada fila representa la conectividad de la tubera en la red. Donde: Columna #1: Nmero del nudo inicial Columna #2: Nmero del nudo final Columna #3: Longitud de la tubera en metros [m] Columna #4: Dimetro de la tubera en milmetros [mm] Columna #5: Sumatoria de los coeficientes de prdidas locales

RED := 1 2 3 4 5 6 7

1 1 3 3 5 4 5 2

2 3 5 4 4 6 6 7

3 300 300 400 500 400 300 350

4 254 203.2 203.2 203.2 203.2 203.2 254

5 0 0 0.5 0 0.5 0 0

2.2 Cota Topogrfica del terreno (CT) [msnm]

CT :=
1 2 3 4 5 6 1 500 530 470 470 470 470

2.3 Demanda en nudos(Qd) [lt/s]

Qd :=
1 2 3 4 5 6 1 0 0 40 50 40 60

2.4 Rugosidad absoluta de la tubera [m]

k s := 0.06 10
3

2.5 Viscocidad cinemtica [m2/s]

:= 1.14 10
6

Comentaros: [email protected] Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil 2 de 27

2.6 Reservorios que abastecen a la red (RSV) Son los nudos de cota piezomtrica conocida y los argumentos son: Donde: Columna #1: Nmero de nudo de cota piezomtrica conocida Columna #2: Cota piezomtrica [m]
RSV :=

1 1 2 1 2

2 500 530

2.7 Definiendo bombas en la red (BMB) Se debe definir el nmero de la tubera y la altura de agua(presin de agua) adicional con la cual colabora la bomba a la red Donde: Columna #1: Nmero de tubera La ecuacin de la bomba es de la forma: = a(Qac^2) + b(Qac) + c, se debe ingresar: Columna #2: Coeficiente "a" de la ecuacin siempre negativo Columna #3: Coeficneinte "b" de la ecuacin de la bomba Columna #4: Coeficiente "c" de la ecuacin de la bomba
BMB :=

1 1 2 1 5

2 0 0

3 0 0

4 0 0

No existe Bombas en la RED!!!!

Mtodo del Gradiente Hidrulico - Argumentos

Mtodo del Gradiente Hidrulico - Resultados Generales Mtodo del Gradiente Hidrulico - Iteraciones Mtodo del Gradiente Hidrulico - Resultados

Comentaros: [email protected] Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil 3 de 27

Abastecimiento de Agua Potable Y Alcantarillado Anlisis y Diseo de Redes de Agua Potable Mtodo del Gradiente Hidrulico
Mtodo del Gradiente Hidrulico - Argumentos Mtodo del Gradiente Hidrulico - Resultados Generales

3. Proceso de clculo Para realizar el clculo de presiones y caudales en la red, es necesario el siguiente planteamiento de matrices y vectores teniendo en ceunta que: Nmero de nudos de cota piezomtrica desconocida:
NN := rows( CT ) rows( RSV) NN = 4

Nmero de tuberas (tramos)

NT := rows( RED) NT = 7

Nmero de nudos de cota piezomtrica conocida

NS := rows( RSV) NS = 2

3.1 Resultados generales Todas las matrices obtenidas en esta seccin se mantienen constante en todo el procedimeinto de diseo. 3.1.1 Obteniendo la matriz de conectividad total (At), su dimensin es NT*(NN+NS) asociada a ca uno de los nudos de la red, con solo dos elementos diferentes de cero en la i-sima fila "-1" en la columna correspondiente al nodo inicial del tramo i "1" en la columna correspondiente al nodo final del tramo i

At :=

for i 1 .. NT ni RED
i, 1 i, 2

nf RED At At At
i , ni i , nf

1 1

1 0 0 At = 0 0 0 0

0 0 0 0 0 0

0 1 0

0 0

1 0 1 1 0 0 0 1 0 0

0 0 0 0 1 0 0 1

1 0 0

1 0 0

1 1 0

1 0

de la matriz At se obtiene las matrices A12 y A10.

3.1.2 Matriz de conectividad A12 asociada a cada uno de los nudos de la red de cota piezomtrica Comentaros: [email protected] desconocida, de dimensin NT*NN
Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil 4 de 27

Los nudos de cota piezomtrica desconocida son(NCPD):

NCPD := submatrix( NODE , rows( RSV) + 1 , rows( NODE) , 1 , 1 ) 3 4 NCPD = 5 6

y la matriz A12 resulta:

A12 := ( NCPD1 , 1) A12 At for n 2 .. rows( NCPD) i NCPD
n, 1

i A12 augment A12 , At A12

1 1 1 A12 = 0 0 0 0
T

0 0 1 1 0 0

0 1 0

0 0 1 0

1 0

1 0

1 1 0

su traspuesta es A21:
1 0 A21 = 0 0 1 1 0 0 1 0 1 0 0 1 0 0 0

A21 := A12

1 0 1 1

1 0

1 0 0

3.1.3 Matriz topolgica tramo a nodo, que asocia a las tuberas con los nodos de cota piezomtrica conocida(Los reservorios) de dimensin NT*NS Los nudos de cota piezomtrica conocida son(NCPC):
1 NCPC := RSV NCPC =

1 2

la matriz A10 resulta:

Comentaros: [email protected] Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil 5 de 27

A10 :=

( NCPC1 , 1) A10 At if rows( NCPC) 2 for n 2 .. rows( NCPC) i NCPC

n, 1

i A10 augment A10 , At A10 A10

1 0 0 A10 = 0 0 0 0

0 0 0 0 0 1

A10 es la matriz topolgica tramo a nodo, para los NS nodos de cota piezomtrica conocida, su dimensin es NT*NS con un valor igual a -1 en las filas correspondientes a los tramos conectados a los reservorios(Nudos de cota piezomtrica conocida) 3.1.4 Vector de Cotas piezomtricas fijas, cuya dimensin es NS*1
2 Ho := RSV Ho =

500 530

3.1.5 Vector de consumo, de dimensin NN*1 En este vector no interviene los nudos de cota piezomtrica conocida.
submatrix( Qd , rows( RSV) + 1 , rows( Qd) , 1 , 1 ) 1000 0.04 0.05 q= 0.04 0.06

q :=

en m3/s

3.1.6 matriz identidad, de dimensin NT*NT

I := identity( NT)

3.1.7 matriz diagonal M, de dimensin NT*NT

Ndw := 2 I

1 0 0 I = 0 0 0 0

0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1

2 0 0 Ndw = 0 0 0 0

0 0 0 0 0 0 2 0 0 0 0 0 0 2 0 0 0 0 0 0 2 0 0 0 0 0 0 2 0 0 0 0 0 0 2 0 0 0 0 0 0 2

los elementos de la diagonal principal son iguales al coeficiente "m", que depende de qu ecuacin para la prdida de carga se est utilizando, en este caso utilizar la de Darcy-Weisbach, para lo cual m=2
Comentaros: [email protected] Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil 6 de 27

3.1.6 Ordenando el coeficiente de las ecuaciones para cada tubera

BOMB := f ( x , y) 0 BOMB matrix( NT , 3 , f ) for i 1 .. rows( BMB) t BMB BOMB
t, 1 t, 2 t, 3 i, 1

BMB BMB BMB

i, 2 i, 3 i, 4

BOMB BOMB BOMB

0 0 0 BOMB = 0 0 0 0

0 0

0 0 0 0 0 0 0 0 0 0 0 0

3.2 Valores iniciales para las iteraciones. 3.2.1 Caudales que circulan en cada tubera
f ( x , y ) := 0.2 Q := matrix( rows( RED) , 1 , f ) Q = ( 0.2 0.2 0.2 0.2 0.2 0.2 0.2 )
T

3.2.2 Dimetro de la tuberas [m]

D := 4 RED 1000

D = ( 0.254 0.203 0.203 0.203 0.203 0.203 0.254 )

Mtodo del Gradiente Hidrulico - Resultados Generales Mtodo del Gradiente Hidrulico - Iteraciones Mtodo del Gradiente Hidrulico - Resultados

Comentaros: [email protected] Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil 7 de 27

Abastecimiento de Agua Potable Y Alcantarillado Anlisis y Diseo de Redes de Agua Potable Mtodo del Gradiente Hidrulico
Mtodo del Gradiente Hidrulico - Argumentos Mtodo del Gradiente Hidrulico - Resultados Generales Mtodo del Gradiente Hidrulico - Iteraciones

4. Proceso Iterativo: 4.1 Iteracin #1 El caudal para la iteracin actual es:

El caudal inicial es Q, luego se toma el caudal resultante Qnext para cada nueva iteracin cambiando de signo si alguno resultase negativo

Qac := Q

0.2 0.2 0.2 Qac = 0.2 0.2 0.2 0.2

1. Obteniendo la matriz A11 Esta matriz contiene en su diagonal principal el siguiente valor:
i Qi
mi 1

i + i + Qi

1.1 Obteniendo el coeficiente

:= for i 1 .. NT 4 Qac Re D
i, 1

i, 1

fa 0.01 fa root

1 ks 2.51 , fa + 2 log + 3.7 D Re fa fa i, 1

i, 3

0.08262686 fa RED
i, 1

(Di, 1)

355.534 3 1.113 10 3 1.485 10 = 1.856 103 3 1.485 10 1.113 103 414.79

Comentaros: [email protected]

Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil

8 de 27

1.2 Prdida de carga localizadas

:= for i 1 .. NT 8 Qac
i, 1

i, 1 4

9.807 D

( i, 1)

RED

i, 5

= ( 0 0 4.848 0 4.848 0 0 )
T

1.3 Cuando existe bombas en la red

:= for i 1 .. NT BOMB
i i, 1

Qac

i, 1

+ BOMB

i, 2

Qac

i, 1

+ BOMB

i, 3

= (0 0 0 0 0 0 0 )
T

La matriz A11 resulta:

A11 := for i 1 .. NT A11 A11 0 0 0 0 0 0 71.107 0 0 0 0 0 222.677 0 0 0 0 0 301.75 0 0 0 0 0 371.128 0 0 0 A11 = 0 0 301.75 0 0 0 0 0 0 222.677 0 0 0 0 0 82.958 0 0 0 0 0
i, i

i, 1

Qac

i, 1

2 1

+
i, 1

i, 1 i, 1

Qac

Vector de cargas piezomtricas

Hnext := A21 Ndw

A11

A12

1 1 A21 Ndw Qac + A11 A10 Ho + q A21 Qac

Hnext = ( 487.201 515.612 488.503 535.889 )

Comentaros: [email protected] Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil 9 de 27

Vector de caudales en las tuberas

Qnext := I Ndw
1 1 1 Qac Ndw A11 ( A12 Hnext + A10 Ho)

0.19 0.097 0.053 0.063 Qnext = 0.066 6.401 10 3 3.294

Comparando los caudales(en listros): 10 102.923 147.077 136.522 ( Qnext Qac ) 1000 = 133.599 193.599 3 3.094 10

la norma del vector es:

Qnext Qac

Error :=

Error = 3.112

Mtodo del Gradiente Hidrulico - Iteraciones Mtodo del Gradiente Hidrulico - Resultados

Comentaros: [email protected] Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil 10 de 27

Abastecimiento de Agua Potable Y Alcantarillado Anlisis y Diseo de Redes de Agua Potable Mtodo del Gradiente Hidrulico
Mtodo del Gradiente Hidrulico - Argumentos Mtodo del Gradiente Hidrulico - Resultados Generales Mtodo del Gradiente Hidrulico - Iteraciones

4. Proceso Iterativo: 4.1 Iteracin #2 El caudal para la iteracin actual es:

El caudal inicial es Q, luego se toma el caudal resultante Qnext para cada nueva iteracin cambiando de signo si alguno resultase negativo

0.19 0.097 0.053 0.063 Qac := 0.066 6.401 10 3 3.294

1. Obteniendo la matriz A11 Esta matriz contiene en su diagonal principal el siguiente valor:
i Qi
mi 1

i + i + Qi

1.1 Obteniendo el coeficiente

:= for i 1 .. NT 4 Qac Re D
i, 1

i, 1

fa 0.01 fa root

1 ks 2.51 , fa + 2 log + 3.7 D Re fa fa i, 1

i, 3

0.08262686 fa RED
i, 1

(Di, 1)

356.565 3 1.157 10 3 1.621 10 = 1.994 103 3 1.589 10 1.68 103 390.43

Comentaros: [email protected]

Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil

11 de 27

1.2 Prdida de carga localizadas

:= for i 1 .. NT 8 Qac
i, 1

i, 1 4

9.807 D

( i, 1)

RED

i, 5

= ( 0 0 1.285 0 1.6 0 0 )
T

1.3 Cuando existe bombas en la red

:= for i 1 .. NT BOMB
i i, 1

Qac

i, 1

+ BOMB

i, 2

Qac

i, 1

+ BOMB

i, 3

= (0 0 0 0 0 0 0 )
T

La matriz A11 resulta:

A11 := for i 1 .. NT A11 A11 0 0 0 0 0 0 67.747 0 0 0 0 0 0 112.205 0 0 0 0 87.194 0 0 0 0 0 125.611 0 0 0 A11 = 0 0 0 106.447 0 0 0 0 10.752 0 0 0 0 0 3 1.286 10 0 0 0 0 0 0
i, i

i, 1

Qac

i, 1

2 1

+
i, 1

i, 1 i, 1

Qac

Vector de cargas piezomtricas

Hnext := A21 Ndw

A11

A12

1 1 A21 Ndw Qac + A11 A10 Ho + q A21 Qac

Hnext = ( 487.128 480.379 478.982 478.646 )

Comentaros: [email protected] Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil 12 de 27

Vector de caudales en las tuberas

Qnext := I Ndw
1 1 1 Qac Ndw A11 ( A12 Hnext + A10 Ho)

0.19 0.085 0.065 Qnext = 0.026 0.041 0.019 1.853

Comparando los caudales(en listros): 1.055 10 12 12.202 12.202 ( Qnext Qac ) 1000 = 37.059 24.858 12.457 1.441 103

la norma del vector es:

Qnext Qac

Error :=

Error = 1.442

Mtodo del Gradiente Hidrulico - Iteraciones Mtodo del Gradiente Hidrulico - Resultados

Comentaros: [email protected] Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil 13 de 27

Abastecimiento de Agua Potable Y Alcantarillado Anlisis y Diseo de Redes de Agua Potable Mtodo del Gradiente Hidrulico
Mtodo del Gradiente Hidrulico - Argumentos Mtodo del Gradiente Hidrulico - Resultados Generales Mtodo del Gradiente Hidrulico - Iteraciones

4. Proceso Iterativo: 4.1 Iteracin #3 El caudal para la iteracin actual es:

El caudal inicial es Q, luego se toma el caudal resultante Qnext para cada nueva iteracin cambiando de signo si alguno resultase negativo

0.19 0.085 0.065 Qac := 0.026 0.041 0.019 1.853

1. Obteniendo la matriz A11 Esta matriz contiene en su diagonal principal el siguiente valor:
i Qi
mi 1

i + i + Qi

1.1 Obteniendo el coeficiente

:= for i 1 .. NT 4 Qac Re D
i, 1

i, 1

fa 0.01 fa root

1 ks 2.51 , fa + 2 log + 3.7 D Re fa fa i, 1

i, 3

0.08262686 fa RED
i, 1

(Di, 1)

356.565 3 1.168 10 3 1.591 10 = 2.202 103 3 1.665 10 1.382 103 391.816

Comentaros: [email protected]

Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil

14 de 27

1.2 Prdida de carga localizadas

:= for i 1 .. NT 8 Qac
i, 1

i, 1 4

9.807 D

( i, 1)

RED

i, 5

= ( 0 0 1.576 0 0.994 0 0 )
T

1.3 Cuando existe bombas en la red

:= for i 1 .. NT BOMB
i i, 1

Qac

i, 1

+ BOMB

i, 2

Qac

i, 1

+ BOMB

i, 3

= (0 0 0 0 0 0 0 )
T

La matriz A11 resulta:

A11 := for i 1 .. NT A11 A11 0 0 0 0 0 67.747 0 0 0 0 0 0 99.246 0 0 0 0 0 104.97 0 0 0 0 0 57.247 0 0 0 A11 = 0 0 69.257 0 0 0 0 0 0 26.263 0 0 0 0 0 726.035 0 0 0 0 0
i, i

i, 1

Qac

i, 1

2 1

+
i, 1

i, 1 i, 1

Qac

Vector de cargas piezomtricas

Hnext := A21 Ndw

A11

A12

1 1 A21 Ndw Qac + A11 A10 Ho + q A21 Qac

Hnext = ( 487.128 479.576 479.381 478.292 )

Comentaros: [email protected] Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil 15 de 27

Vector de caudales en las tuberas

Qnext := I Ndw
1 1 1 Qac Ndw A11 ( A12 Hnext + A10 Ho)

0.19 0.082 0.068 Qnext = 0.011 0.03 0.03 1.291

Comparando los caudales(en listros): 1.055 10 12 3.471 3.471 ( Qnext Qac ) 1000 = 14.703 11.231 11.231 561.504

la norma del vector es:

Qnext Qac

Error :=

Error = 0.562

Mtodo del Gradiente Hidrulico - Iteraciones Mtodo del Gradiente Hidrulico - Resultados

Comentaros: [email protected] Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil 16 de 27

Abastecimiento de Agua Potable Y Alcantarillado Anlisis y Diseo de Redes de Agua Potable Mtodo del Gradiente Hidrulico
Mtodo del Gradiente Hidrulico - Argumentos Mtodo del Gradiente Hidrulico - Resultados Generales Mtodo del Gradiente Hidrulico - Iteraciones

4. Proceso Iterativo: 4.1 Iteracin #4 El caudal para la iteracin actual es:

El caudal inicial es Q, luego se toma el caudal resultante Qnext para cada nueva iteracin cambiando de signo si alguno resultase negativo

0.19 0.082 0.068 Qac := 0.011 0.03 0.03 1.291

1. Obteniendo la matriz A11 Esta matriz contiene en su diagonal principal el siguiente valor:
i Qi
mi 1

i + i + Qi

1.1 Obteniendo el coeficiente

:= for i 1 .. NT 4 Qac Re D
i, 1

i, 1

fa 0.01 fa root

1 ks 2.51 , fa + 2 log + 3.7 D Re fa fa i, 1

i, 3

0.08262686 fa RED
i, 1

(Di, 1)

356.565 3 1.171 10 3 1.585 10 = 2.523 103 3 1.728 10 1.296 103 393.175

Comentaros: [email protected]

Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil

17 de 27

1.2 Prdida de carga localizadas

:= for i 1 .. NT 8 Qac
i, 1

i, 1 4

9.807 D

( i, 1)

RED

i, 5

= ( 0 0 1.648 0 0.727 0 0 )
T

1.3 Cuando existe bombas en la red

:= for i 1 .. NT BOMB
i i, 1

Qac

i, 1

+ BOMB

i, 2

Qac

i, 1

+ BOMB

i, 3

= (0 0 0 0 0 0 0 )
T

La matriz A11 resulta:

A11 := for i 1 .. NT A11 A11
i, i

i, 1

Qac

i, 1

2 1

+
i, 1

i, 1 i, 1

Qac

67.747 0 0 A11 = 0 0 0 0

0 0 0 0 0 96 0 0 0 0 0 109.397 0 0 0 27.754 0 0 0 0 52.581 0 0 0 0 38.89 0 0 0 0 507.589 0 0 0 0 0

0 0 0 0 0 0

Vector de cargas piezomtricas

Hnext := A21 Ndw

A11

A12

1 1 A21 Ndw Qac + A11 A10 Ho + q A21 Qac

Hnext = ( 487.128 479.393 479.516 478.122 )

Comentaros: [email protected] Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil 18 de 27

Vector de caudales en las tuberas

Qnext := I Ndw
1 1 1 Qac Ndw A11 ( A12 Hnext + A10 Ho)

0.19 0.081 0.069 3 Qnext = 7.725 10 0.027 0.033 1.168

Comparando los caudales(en listros): 1.499 10 12 1.355 1.355 ( Qnext Qac ) 1000 = 3.275 2.92 2.92 123.424

la norma del vector es:

Qnext Qac

Error :=

Error = 0.124

Mtodo del Gradiente Hidrulico - Iteraciones Mtodo del Gradiente Hidrulico - Resultados

Comentaros: [email protected] Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil 19 de 27

Abastecimiento de Agua Potable Y Alcantarillado Anlisis y Diseo de Redes de Agua Potable Mtodo del Gradiente Hidrulico
Mtodo del Gradiente Hidrulico - Argumentos Mtodo del Gradiente Hidrulico - Resultados Generales Mtodo del Gradiente Hidrulico - Iteraciones

4. Proceso Iterativo: 4.1 Iteracin #5 El caudal para la iteracin actual es:

El caudal inicial es Q, luego se toma el caudal resultante Qnext para cada nueva iteracin cambiando de signo si alguno resultase negativo

0.19 0.081 0.069 Qac := 7.725 10 3 0.027 0.033 1.168

1. Obteniendo la matriz A11 Esta matriz contiene en su diagonal principal el siguiente valor:
i Qi
mi 1

i + i + Qi

1.1 Obteniendo el coeficiente

:= for i 1 .. NT 4 Qac Re D
i, 1

i, 1

fa 0.01 fa root

1 ks 2.51 , fa + 2 log + 3.7 D Re fa fa i, 1

i, 3

0.08262686 fa RED
i, 1

(Di, 1)

356.565 3 1.172 10 3 1.583 10 = 2.696 103 3 1.753 10 1.281 103 393.642

Comentaros: [email protected]

Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil

20 de 27

1.2 Prdida de carga localizadas

:= for i 1 .. NT 8 Qac
i, 1

i, 1 4

9.807 D

( i, 1)

RED

i, 5

= ( 0 0 1.673 0 0.654 0 0 )
T

1.3 Cuando existe bombas en la red

:= for i 1 .. NT BOMB
i i, 1

Qac

i, 1

+ BOMB

i, 2

Qac

i, 1

+ BOMB

i, 3

= (0 0 0 0 0 0 0 )
T

La matriz A11 resulta:

A11 := for i 1 .. NT A11 A11 0 0 0 0 0 67.747 0 0 0 0 0 0 94.917 0 0 0 0 0 110.871 0 0 0 0 0 20.829 0 0 0 A11 = 0 0 47.972 0 0 0 0 0 0 42.27 0 0 0 0 0 459.774 0 0 0 0 0
i, i

i, 1

Qac

i, 1

2 1

+
i, 1

i, 1 i, 1

Qac

Vector de cargas piezomtricas

Hnext := A21 Ndw

A11

A12

1 1 A21 Ndw Qac + A11 A10 Ho + q A21 Qac

Hnext = ( 487.128 479.38 479.524 478.108 )

Comentaros: [email protected] Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil 21 de 27

Vector de caudales en las tuberas

Qnext := I Ndw
1 1 1 Qac Ndw A11 ( A12 Hnext + A10 Ho)

0.19 0.081 0.069 3 Qnext = 7.314 10 0.027 0.033 1.16

Comparando los caudales(en listros): 2.22 10 13 0.442 0.442 ( Qnext Qac ) 1000 = 0.411 0.244 0.244 7.629

la norma del vector es:

Qnext Qac

Error :=

Error = 7.674 10

Mtodo del Gradiente Hidrulico - Iteraciones Mtodo del Gradiente Hidrulico - Resultados

Comentaros: [email protected] Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil 22 de 27

Abastecimiento de Agua Potable Y Alcantarillado Anlisis y Diseo de Redes de Agua Potable Mtodo del Gradiente Hidrulico
Mtodo del Gradiente Hidrulico - Argumentos Mtodo del Gradiente Hidrulico - Resultados Generales Mtodo del Gradiente Hidrulico - Iteraciones

4. Proceso Iterativo: 4.1 Iteracin #6 El caudal para la iteracin actual es:

El caudal inicial es Q, luego se toma el caudal resultante Qnext para cada nueva iteracin cambiando de signo si alguno resultase negativo

0.19 0.081 0.069 Qac := 7.314 10 3 0.027 0.033 1.16

1. Obteniendo la matriz A11 Esta matriz contiene en su diagonal principal el siguiente valor:
i Qi
mi 1

i + i + Qi

1.1 Obteniendo el coeficiente

:= for i 1 .. NT 4 Qac Re D
i, 1

i, 1

fa 0.01 fa root

1 ks 2.51 , fa + 2 log + 3.7 D Re fa fa i, 1

i, 3

0.08262686 fa RED
i, 1

(Di, 1)

356.565 3 1.172 10 3 1.583 10 = 2.726 103 3 1.753 10 1.281 103 393.676

Comentaros: [email protected]

Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil

23 de 27

1.2 Prdida de carga localizadas

:= for i 1 .. NT 8 Qac
i, 1

i, 1 4

9.807 D

( i, 1)

RED

i, 5

= ( 0 0 1.673 0 0.654 0 0 )
T

1.3 Cuando existe bombas en la red

:= for i 1 .. NT BOMB
i i, 1

Qac

i, 1

+ BOMB

i, 2

Qac

i, 1

+ BOMB

i, 3

= (0 0 0 0 0 0 0 )
T

La matriz A11 resulta:

A11 := for i 1 .. NT A11 A11 0 0 0 0 0 67.747 0 0 0 0 0 0 94.917 0 0 0 0 0 110.871 0 0 0 0 0 19.935 0 0 0 A11 = 0 0 47.972 0 0 0 0 0 0 42.27 0 0 0 0 0 456.664 0 0 0 0 0
i, i

i, 1

Qac

i, 1

2 1

+
i, 1

i, 1 i, 1

Qac

Vector de cargas piezomtricas

Hnext := A21 Ndw

A11

A12

1 1 A21 Ndw Qac + A11 A10 Ho + q A21 Qac

Hnext = ( 487.128 479.379 479.524 478.108 )

Comentaros: [email protected] Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil 24 de 27

Vector de caudales en las tuberas

Qnext := I Ndw
1 1 1 Qac Ndw A11 ( A12 Hnext + A10 Ho)

0.19 0.081 0.069 3 Qnext = 7.301 10 0.027 0.033 1.16

Comparando los caudales(en listros): 1.499 10 12 0.446 0.446 ( Qnext Qac ) 1000 = 0.013 0.253 0.253 0.296

la norma del vector es:

Qnext Qac

Error :=

Error = 7.829 10

Mtodo del Gradiente Hidrulico - Iteraciones Mtodo del Gradiente Hidrulico - Resultados

Comentaros: [email protected] Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil 25 de 27

Abastecimiento de Agua Potable Y Alcantarillado Anlisis y Diseo de Redes de Agua Potable Mtodo del Gradiente Hidrulico
Mtodo del Gradiente Hidrulico - Argumentos Mtodo del Gradiente Hidrulico - Resultados Generales Mtodo del Gradiente Hidrulico - Iteraciones Mtodo del Gradiente Hidrulico - Resultados

5. Ordenado Resultados Programa que Corrige H y Q con los argumentos establecidos en el captulo 3, culmina cuando la norma del vector es menor a 0.0001 H := f ( x , y) 0.2 Q
DQ Qan H Qan Q Qan while DQ > 0.0001 for i 1 .. NT Re 4 Qan D
i, 1

Qan matrix( NT , 1 , f )

i, 1

fa 0.01 fa root

1 ks 2.51 , fa + 2 log + 3.7 D Re fa fa i, 1

i, 3

0.08262686 fa RED

(Di, 1)
8 Qan
2 i, 1

( i, 1) 2 BOMB ( Qan ) + BOMB Qan + BOMB i, 1 i, 1 i, 2 i, 1 i, 3

9.807 D
4

RED

i, 5

A11

i, i

Qan

i, 1

2 1

+ + A12
1

Qan
1 i, 1

H A21 Ndw
1

A11

1 1 A21 Ndw Qan + A11 A10 Ho + q A21 Qan

Q I Ndw Qan Ndw DQ ( Q Qan ) Qan Q

A11

( A12 H + A10 Ho)

H Q
Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil

Comentaros: [email protected] 26 de 27

Los caudales resultantes(que circulan) en cada tubera son(en litros/s):

El signo negativo indica el flujo del caudal en sentido contrario al supuesto inicialmente.

1000 Q = 170 69.744 60.256 6.982 17.238 22.762 1.16 10

Las cotas piezomtricas en cada nudo son(en metros):

Hf := augment Ho , H

Hf = ( 500 530 489.627 483.724 483.858 483.161 )

Las presiones en los puntos son(en metros):

P := Hf CT
T

P = ( 0 0 19.627 13.724 13.858 13.161 )

Mtodo del Gradiente Hidrulico - Resultados

Comentaros: [email protected] Universidad Nacional San Cristbal de Huamanga Escuela Profesional de Ingeniera civil 27 de 27