9.

2 MODELACIÓN H I D R Á U LI C A DE ACCESORIOS ESPECIALES EN RDAPS

3
3 4

2 7 6

2 6

Figura 9.10 E s q u e m a de la red Ejemplo.

Siguiendo el método del gradiente se obtienen estas matrices:

'
[A12] = matriz de conectividad con dimensión ( 7 x 5 ) .

1 o o o o

-1 1 o o o

o o o o

o -1 1 o o

o o -1 1 o

o o o 1 -1

o o -1 o 1

[A21] = matriz transpuesta de [A12].

1 - 1 o o o o o

o 1 1 - 1 o o o

o o o 1 - 1 o - 1

o o o o 1 1 o

o o o o - 1 o 1

Las otras matrices del método son:

( 709 )
 iEAA ■iiiM

Modelación avanzada de sistemas de tuberías

01
..J
[A10] = matriz topológica tramo a nudo, cuya dimensión es ( 7 x 1 ) .
c. i

�I [Q] = vector de caudales, cuya dimensión es ( 7 x 1 ) . Sus unidades son

.;;-1
[H] = vector de alturas piezométricas desconocidas, cuya dimension es
1

[H] =vector de alturas piezométricas fijas, cuya dimensión es ( 1 x 1

o

[q] = vector de consumo, cuya dimensión es ( 5 x 1 ) . Sus unidades soz

[A10] [Q] [H] [Hº] [q]

-1 0.01

o 0.01 H, 0.03

-1 0.01 0.02
H3

o 0.01 H4 0.03

o 0.01 Hs 0.04

o 0.01 H6 0.03

o 0.01

[N] = matriz diagonal, cuya dimensión es ( 7 x 7 ) . Presenta el valor 2 en la

diagonal, ya que se utiliza la ecuación de Darcy-Weissbach como la

ecuación de fricción.

2 o o o o o o

o 2 o o o o o

o o 2 o o o o

o o o 2 o o o

o o o o 2 o o

o o o o o 2 o

o o o o o o 2

[I] = matriz identidad, cuya dimensión es ( 7 x 7 ) .

1 o o o o o o

o 1 o o o o o

o o 1 o o o o

o o o 1 o o o

o o o o 1 o o

o o o o o 1 o

o o o o o o 1

: Primera iteración

1
Las matrices anteriormente definidas son válidas para todas las iteraciones
1

1 necesarias en el cálculo de la red. Las que se presentan a continuación, varían de

( 7 1 0 )
 1
S una iteración a otra. La forma en que se desarrolla el ejemplo es similar a lo hecho
a. l

�1 en los ejemplos anteriores que hacen uso del método del gradiente.
wl
1

[ A 11 ] = matriz diagonal cuya dimensión es ( 7 x 7 ) . Presenta el valor

a.Q.(n,-1) + A _ + ! ;
' f .., , Q;

en la diagonal; el coeficiente /3 es O para todas las tuberías y el

coeficiente y toma valor diferente de O solo para la tubería 4, que tiene

la VRP.

La siguiente tabla resume el cálculo de los coeficientes a.

Tubería Q (m3/s) f (-) v (m/s) h (m) h +hm (m) a

1 1

1-2 0.1 0.01305 3.084 6.22 6.22 622.46

2-3 0.1 0.01247 5.482 37.61 37.61 3760.52

1-3 0.1 0.01305 3.084 6.22 6.22 622.46

3-4 0.1 0.01305 3.084 6.22 6.22 622.46

4-5 0.1 0.01247 5.482 25.07 25.07 2507.01

6-5 0.1 0.01247 5.482 37.61 37.61 3760.52

4-6 0.1 0.01247 5.482 25.07 25.07 2507.01

Cilculo del coeficiente y de la tubería 4

1
La presión objetivo de la válvula es 30 m, por lo que las pérdidas de altura causadas
1

1 por la válvula tomarán el valor necesario para que esto se cumpla. Como en la
1
primera iteración no se conoce la presión del nudo anterior a la válvula, no se puede

calcular este valor, por lo que inicialmente se supone que y = O. La Matriz [ A 11 ] es

1 la siguiente:

62.246 o o o o o o

o 376.052 o o o o o

o o 62.246 o o o o

[Au]= o o o 62.246 o o o

o o o o 250.701 O o

o o o o o 376.052 o

o o o o o o

( 711 )
 ti·iiiiil■iU

Modelación avanzada de sistemas de tuberías

01
...J
[A11]' = matriz diagonal cuya dimensión es ( 7 x 7 ) . Presenta el valor
o.. 1

�I
la diagonal. Esta matriz se presenta a continuación:
¡;;I

62.246 o o o o o o

o 376.052 o o o o o

o o 62.246 o o o o

[ A 11 ]' = o o o 62.246 o o o

o o o o 250.701 o o

o o o o o 376.052 o

o o o o o o 250.401

En este caso [A11] es igual a [A11]', ya que el término de pérdidas de =rrn.a.a

1

1
la primera iteración es O.
1

alculo de [H,.J
1

1
Siguiendo el método del gradiente se obtienen las siguientes matrices con el ñ=
1

1 calcular las alturas piezornétricas en cada uno de los nudos de la red para la pi

: iteración.

[N)[A11)'

124.492 o o o o o o

o 752.103 o o o o o

o o 124.492 o o o o

o o o 124.492 o o o

o o o o 501.402 o o

o o o o o 752.103 o

o o o o o o 501.402

([N] [ A 1 1 ]')" •

0.008 o o o o o o

o 0.0013 o o o o o

o o 0.008 o o o o

o o o 0.008 o o o

o o o o 0.002 o o

o o o o o 0.0013 o

o o o o o o 0.002

( 712 )
 9.2 MODELACIÓN H I D R Á U LI C A DE ACCESORIOS ESPECIALES EN RDAPS

01 [A21) ([N][A11)')-•

ci'.1

�¡
0.008 -0.0013 o o o o o
;;;-1

1
o 0.0013 o o o o o

o o 0.008 - 0.008 o o o

o o o 0.008 o o o

o o o o - 0.002 o - 0.002

o o o o 0.002 0.0013 o

o o o o o - 0.0013 0.002

[A21) ([N)[A11)')-•[A12)

0.0094 -0.0013 o o o

- 0.0013 0.0174 - 0.008 o o

o - 0.008 0.012 - 0.002 -0.002

o o - 0.002 0.0033 -0.0013

o o -0.002 -0.0013 0.0033

- ([A21) ([N)[A11)')-'[A12))·•

- 109.01 - 15.48 -15.48 - 15.48 - 1 5 -4 8

- 15.48 -109.01 -109.01 - 109.01 -109.01

- 15.48 - 109.01 -203.50 - 203.50 - 203.50

- 15.48 - 109.01 -203.50 - 5 9 1. 6 5 - 376.76

- 15.48 - 109.01 -203.50 - 376.76 - 5 9 1. 6 5

[A11)[Q) [A1o)[H 0) ([A11][Q)) ([A1o][H

0))

6.225 -50 -43.775

37.605 o 37.605

6.225 -50 -43,775

6.225 o 6.225

25.070 o 25.070

37.605 o 37,605

25.070 o 25.070

[A21) ([N)[A11)')''(([A11][Q)) + ([A1o)[H 0))) [A21)[Q) - [q)

-0.402 -0.03

-0.352 [Link]

-0.05 -0.13

0.1 0.16

o -0.03

( 713 )
 li·IAIIHI ii

Modelación avanzada de sistemas de tuberías

01
.., [A21) ([N][A11]')-'(([A11][Q)) + ([A1o][H 0 ))) - ([A21][Q) - [q))
a. 1

-0.372
�·
w'
-0-432

[Link]

-0.060

0.030

Finalmente, al aplicar la Ecuación ( 7 . 6 7 ) :

H, .. = - ([A21) ([N][A11)')"•[A12))·• ([A21) ([N)[A11)')"•(([A11][Q)) +

([A1o][H))) - ([A21)[Q) - [q)))

46.42

47,355

4 1. 1 3 1

58.322

38.982
1

: Con estas alturas ( dadas en m) se procede a calcular los caudales de la -:- -­

, iteración, tal como se muestra a continuación:

1

: Cálculo de [Q,.,)

1 Siguiendo paso a paso el método del gradiente, se obtienen las matrices

1 • • '

1
muestran a continuación:

[A12)[H¡+ 1 J + [A1o][H 0)

46.420 - 3.580

0.935 0.935

47.355 - 2.645

- 6.225 - 6.225

17 . 1 9 1 17 . 1 9 1

19.340 19.340

- 2.149 - 2.149

([N)[A11)')"') x ([A12l[H¡u) + [A1o][H0))

-0.0288

0.0012

-0.0212

-0.05

0.0343

0.0257

-0.0043

( 7 1 4 )
 9.2 MODELACIÓN H I D R Á U LI C A DE ACCESORIOS ESPECIALES EN RDAPS

01 ([N)[Au]')-') [Au]
..J

"-' 0.5 o o o o o o
�' o o o o o o
¡;j'I 0.5

1
o o 0.5 o o o o
1

o o o 0.5 o o o

o o o o 0.5 o o

o o o o o 0.5 o

o o o o o o 0.5

I - (([N)[Au)')·•) [Au))

0.5 o o o o o o

o 0.5 o o o o o

o o 0.5 o o o o

o o o 0.5 o o o

o o o o 0.5 o o

o o o o o 0.5 o

o o o o o o 0.5

(1 - (([N)[Au]')-') [Au))) x [Q)

0.5

0.5

0.5

0.5

0.5

0.5

0.5

3
, Finalmente se obtienen los siguientes caudales (en m /s), utilizando la Ecuación

1 .

' z.68:

Q,., = ((I - (([Nl[Au]')-') [Au))) x [Q)) - (([N)[Au]')-') x ([A12l[H, • .J + [A1o)[H )))

0

0.079

0.049

0.071

0.100

0.016

0.024

0.054

3
Estos caudales están en m /s.

( 715 )