Metodo del Gradiente Hidraulico

Aplicacion

Docente: Ing. Anderson Lincol Condori Paytan

Curso: Abastecimiento de Agua y Alcantarillado

UNIVERSIDAD NACIONAL DE HUANCAVELICA
Escuela Profional de Ingenieria Civil Huancavelica
Computational Hydraulics

2017 - II
Huancavelica - Peru

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 1 / 34

Indice

1 Metodo de Gradiente Hidraulico

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 2 / 34

 Metodo de Gradiente Hidraulico

Esquema de la red de distribucion

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 3 / 34

 Metodo de Gradiente Hidraulico

Datos de la red de agua

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 4 / 34

 Metodo de Gradiente Hidraulico

Cotas de Terreno y Cotas piezometricas conocidas

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 5 / 34

 Metodo de Gradiente Hidraulico

Matriz de Datos Iniciales: NT = 14, NN = 11 y NS = 2

Ni Nf L(m) D(mm) Km Ks(m)
1
5 100 337,9 0 6,00E 05

5
6 100 299,6 0 6,00E 05

6
7 100 190,2 0 6,00E 05

6
8 100 292,2 0 6,00E 05

4
7 100 190,2 0 6,00E 05

4
8 100 231,8 0 6,00E 05

2
3 100 292,2 0 6,00E 05

3
4 100 231,8 0 6,00E 05

8
9 150 292,2 0 6,00E 05

9
10 110 231,8 0 6,00E 05

9
11 110 231,8 0 6,00E 05

10
12 110 185,4 0 6,00E 05

11 12 110 185,4 0 6,00E 05
12 13 100 185,4 0 6,00E 05

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 6 / 34

 Metodo de Gradiente Hidraulico

Matriz de cotas y Qdemanda

cota(msnm) Q(Lt/s)

985 30

960 30

985 30

960 30

940 30

930 30

930 30

930 30

930 30
930 30
920 80

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 7 / 34

 Metodo de Gradiente Hidraulico

Vector de carga Piezometrica conocida NSx1

cota(msnm) Q(Lt/s)
1 1002
2 1002

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 8 / 34

 Metodo de Gradiente Hidraulico

Matriz de Conectividad A12 (NT x NN)

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

0
0 0 1 1 0 0 0 0 0 0

0
0 0 1 0 1 0 0 0 0 0

0
1 0 0 1 0 0 0 0 0 0

0
1 0 0 0 1 0 0 0 0 0

1 0 0 0 0 0 0 0 0 0 0

1 1 0 0 0 0 0 0 0 0 0

0
0 0 0 0 1 1 0 0 0 0

0
0 0 0 0 0 1 1 0 0 0

0
0 0 0 0 0 1 0 1 0 0

0
0 0 0 0 0 0 1 0 1 0

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

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 9 / 34

 Metodo de Gradiente Hidraulico

Matriz A21

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

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

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

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

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

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

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

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

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 10 / 34

 Metodo de Gradiente Hidraulico

Matriz A10 (Topologica)

1

0
0 0

0 0

0 0

0 0

0 0

0
1
0 0

0 0

0 0

0 0

0 0

0 0
0 0

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 11 / 34

 Metodo de Gradiente Hidraulico

Matriz [N] con el exponente m = 2 (D-W) NT x NT

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

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

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

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

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

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

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

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

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

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

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

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

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 12 / 34

 Metodo de Gradiente Hidraulico

Matriz Identidad NT x NT

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

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

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

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

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

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

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

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

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

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

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

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

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 13 / 34

 Metodo de Gradiente Hidraulico

Vector de Caudales Iniciales [Qi ] NT x 1 (m3/s)

0,2
0,2

0,2

0,2

0,2

0,2

0,2

0,2

0,2

0,2

0,2

0,2

0,2
0,2

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 14 / 34

 Metodo de Gradiente Hidraulico

Matriz A11, primera Iteracion

0,0826268f L Ks 2,51 2
= , f = (2log( + )) , = 0, = 0
D5 3,7D Re f

5,5 0 0 0 0 0 0 0 0 0 0 0 0

0 10,1 0 0 0 0 0 0 0 0 0 0 0

0 0 103,8 0 0 0 0 0 0 0 0 0 0

0 0 0 11,5 0 0 0 0 0 0 0 0 0

0 0 0 0 103,8 0 0 0 0 0 0 0 0

0 0 0 0 0 37,5 0 0 0 0 0 0 0

0 0 0 0 0 0 11,5 0 0 0 0 0 0

0 0 0 0 0 0 0 37,5 0 0 0 0 0

0 0 0 0 0 0 0 0 17,2 0 0 0 0

0 0 0 0 0 0 0 0 0 41,3 0 0 0

0 0 0 0 0 0 0 0 0 0 41,3 0 0

0 0 0 0 0 0 0 0 0 0 0 130,2 0
0 0 0 0 0 0 0 0 0 0 0 0 130,2
0 0 0 0 0 0 0 0 0 0 0 0 0

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 15 / 34

 Metodo de Gradiente Hidraulico

Calculo de Alturas de Gradiente y Caudales

[Hi+1 ] = ([A21][N ]1 [A11]1 [A12])1 {[A21][N ]1 ([Qi ] + [A11]1 [A10][Ho ])+

+[q] [A21][Qi ]} (1)

Qi+1 = ([I] [N ]1 )[Qi ] [N ]1 [A11]1 ([A12][Hi+1 ] + [A10][Ho ]) (2)

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 16 / 34

 Metodo de Gradiente Hidraulico

Lenguaje de programacion Python para calculo

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 17 / 34

 Metodo de Gradiente Hidraulico

Resultados Por Metodo de Darcy - Weisbach (friccion) - Tuberas

T ub Q f vel Hf
1 236,82 0,014540 2,641 1,5301

2 206,82 0,014801 2,934 2,1678

3 27,95 0,017955 0,984 0,4657

4 148,87 0,015162 2,220 1,3039

5 2,05 0,024753 0,515 0,4701

6 81,13 0,016093 1,922 1,3082

7 143,18 0,015207 2,135 1,2096

8 113,18 0,015701 2,682 2,4839

9 200,00 0,014862 2,982 3,4600

10 85,00 0,016032 2,014 1,5737

11 85,00 0,016032 2,014 1,5737

12 55,00 0,016801 2,037 2,1094

13 55,00 0,016801 2,037 2,1094
14 80,00 0,016359 2,963 3,9505

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 18 / 34

 Metodo de Gradiente Hidraulico

Resultados Por Metodo de Darcy - Weisbach (friccion) - Nodos

N Hi P Qdem
3 1000,791 15,791 30,00

4 998,308 38,308 30,00

5 1000,470 15,470 30,00

6 998,303 38,303 30,00

7 997,838 57,838 30,00

8 997,000 67,000 30,00

9 993,541 63,541 30,00

10 991,968 61,968 30,00

11 991,968 61,968 30,00

12 989,859 59,859 30,00
13 985,910 65,910 80,00

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 19 / 34

 Metodo de Gradiente Hidraulico

Resultados Por Metodo de Hazen - Williams (friccion) - Tuberas

T ub Q f vel Hf
1 236,82 0,014540 2,641 1,530

2 206,82 0,014801 2,934 2,168

3 27,89 0,017959 0,982 0,464

4 148,93 0,015162 2,221 1,305

5 2,11 0,024638 0,530 0,495

6 81,07 0,016094 1,921 1,306

7 143,18 0,015207 2,135 1,210

8 113,18 0,015701 2,682 2,484

9 200,00 0,014862 2,982 3,460

10 85,00 0,016032 2,014 1,574

11 85,00 0,016032 2,014 1,574

12 55,00 0,016801 2,037 2,109

13 55,00 0,016801 2,037 2,109
14 80,00 0,016359 2,963 3,950

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 20 / 34

 Metodo de Gradiente Hidraulico

Resultados Por Metodo de Hazen - Williams (friccion) - Nodos

N Hi P Qdem
3 1000,910 15,910 30,00

4 998,730 38,730 30,00

5 1000,636 15,636 30,00

6 998,728 38,728 30,00

7 998,302 58,302 30,00

8 997,555 67,555 30,00

9 994,518 64,518 30,00

10 993,107 63,107 30,00

11 993,107 63,107 30,00

12 991,238 61,238 30,00
13 987,835 67,835 80,00

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 21 / 34

 Metodo de Gradiente Hidraulico

Grafica de Velocidades

Grafica de velocidad
3.0 2.98 2.96
2.93

2.64 2.68

2.5

2.22
2.14
2.01 2.01 2.04 2.04
2.0
1.92
Velocidad

1.5

1.0 0.98

0.53
0.5 Velocidad
2 4 6 8 10 12 14
Tuberias

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 22 / 34

 Metodo de Gradiente Hidraulico

Grafica de Presiones

Grafica de Presiones
70 Presiones 67.56
64.52
63.11 63.11
61.24
60 58.3

50
Presion

40 38.73 38.73

30

20
15.91 15.64

10
4 6 8 10 12
Nodos

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 23 / 34

 Metodo de Gradiente Hidraulico

Grafica de Perdidas

Grafica de Perdidas
4.0 Perdidas 3.95

3.5 3.46

3.0

2.5 2.48
Perdidas

2.17 2.11 2.11

2.0

1.53 1.57 1.57

1.5
1.3 1.31
1.21
1.0

0.5 0.46 0.5

2 4 6 8 10 12 14
Tuberias
Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 24 / 34
 Metodo de Gradiente Hidraulico

Grafica del Factor de Friccion

Grafica del Factor de Friccion

0.0246 Friccion
0.024

0.022

0.020
Friccion(f)

0.018 0.018

0.0168 0.0168
0.0164
0.016 0.0161 0.016 0.016
0.0157
0.0152 0.0152
0.0148 0.0149
0.0145

2 4 6 8 10 12 14
Tuberias

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 25 / 34

 Metodo de Gradiente Hidraulico

Software Watercad (civil 3D)

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 26 / 34

 Metodo de Gradiente Hidraulico

Resultados Watercad - Darcy Weisbach

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 27 / 34

 Metodo de Gradiente Hidraulico

Resultados Watercad - Darcy Weisbach

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 28 / 34

 Metodo de Gradiente Hidraulico

Resultados Watercad - Darcy Weisbach

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 29 / 34

 Metodo de Gradiente Hidraulico

Resultados Watercad - Darcy Weisbach

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 30 / 34

 Metodo de Gradiente Hidraulico

Resultados Watercad - Hazen Williams

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 31 / 34

 Metodo de Gradiente Hidraulico

Resultados Watercad - Hazen Williams

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 32 / 34

 Metodo de Gradiente Hidraulico

Resultados Watercad - Hazen Williams

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 33 / 34

 Metodo de Gradiente Hidraulico

Resultados Watercad - Hazen Williams

Ing. Anderson Lincol Metodo del Gradiente Hidraulico . . . Noviembre, 2017 34 / 34