EVALUACIN DE CARGAS Y VELOCIDADES EN UN

SISTEMA ACOPLADO RESERVORIO CANAL DE

ADUCCIN DISIPADOR CILNDRICO PARABLICO
TUBERA FORZADA

CURSO: FLUNO NO PERMANENTE

 Ecuacin de Continuidad:

Ecuacin de Cantidad de Movimiento:

 Solucin por el Mtodo de las Caractersticas:
Resolviendo por el mtodo de las caractersticas en la tubera,
tenemos las siguientes ecuaciones:

R S

D
 =-Hp-ZVp

=-Hp+ZVp

Dnde:

Encontramos las siguientes relaciones para la entrada, la tubera y a la salida de la
vlvula

En la entrada de la tubera:

En la tubera y en el canal:

En la vlvula:


 ,

Cierreparcial:

Aperturaparcial:

TR:tiempoderetorno,Tc:tiempodecierre
 Eneldisipador

Se tiene el siguiente sistema acoplado (ver figura) con los datos abajo
mencionados. Se nos pide hallar el valor de la carga (H) y la velocidad (V) en
cada punto sealado (1, 2, 3, 4ED, 4D, 4SD, 5, 6, 7, 8, 9). Para ello se va usar el
mtodo de las caractersticas. Luego se va a proceder a realizar grficos en
funcin del tiempo. Con ello se espera entender cuantitativamente cul es el
comportamiento que tiene el agua dentro de la tubera cuando un golpe de
ariete pasa a travs del sistema acoplado.
Canal Aduccin
Clculos Iniciales

En el Re servorio :
HT 10 m Bn 2.5
En el canal de aduccin
S1 100 m
L1 600 m
D1 3.5 m
f1 0.025
c1 1200 m / s
0
20 o
En el Disipador
HP 50 m HC ?
DP 6 m DC 2.5 m
En la tubera forzada
S 80 m2
L 300 m
2

D 2 m
2
f 0.02
2

c 1000 m / s
2
80 10
1
o
2
o

30 3
o
5 60
4
o
5
o
Variacin del tiempo
Adems :
t 0 s
i

t 400 s (6.67 min .) nf 0 (cierre total)

f

Clculo de t

En el tnel de conduccin

S 100
t 1
1
0.08333 s
c 1200
1
 En la tubera forzada

S 80
t 2
2
0.08 s
c 2
1000
Seescogeelmenorvalor

t 0.08 s
Calculo de la celeridad de la onda promedio en el sistema

L L 600 300
C
1 2
1125 s
m
L L 600 300
1 2

c c 1
1200 1000
2

2 L 2 * 900
TR 1.6 s
C 1125
m

La altura del cilindro que forma parte del disipador

HC HT L1 sen 0 HP 10 600 * Seno(20) 50 165.21 m

 En el tnel/canal de conduccin
G G G G
0 1 1 2 2 3 3 4 ED
t sen 0.08 * Seno(20) 0.02736161 s
0

c1 1200
Z1 122.324159 s
g 9.81

t f 0.08 * 0.025
R c 1200 * 1
0.03494976
2 g D
1 1
2 * 9.81* 3.5
1

3.1416
A D
1
* 3.5 9.6211275 m
2
1
2 2

4 4
 En la tubera forzada
G4 ED 5
t sen 0.08 * Seno(80) 0.07878462 s
1

G t sen 0.08 * Seno(10) 0.0138919 s

5 6 2

G t sen 0.08 * Seno(30) 0.04 s

67 3

G t sen 0.08 * Seno(5) 0.069725 s

7 8 4

G t sen 0.08* Seno(60) 0.06928203 s

8 9 5

c 1000 t f 0.08 * 0.02

Z 101.936799 s
2
R c 2
1000 * 0.04077472
2 g D
2 2
2
g 9.81 2
2 * 9.81* 2

3.1416
A D 2
* 2 3.14159265 m
2 2

4 4
2 2
Clculos de reas transversales
Para el tronco paraboloide (Z < HP)
Ecuacin de la Parbola
P : Z z c X . . . (1)
0
2

De la figura
DC DP
( ,H ) P y (
P
,0) P
2 2
En (1) :
2
DC
HP z 0 c . . . (a)
2
2
DP
0 z0 c . . . (b)
2
De (a ) y (b) :
DP 2
6 2

z HP 50 * 60.50. . . (c)
DP DC 6 2.5
0

2 2 2

4 HP 4 * 50

c HP 50 * 336.13. . . (d )
DP DC 6 2.5
2 2 2 2

De (1)
Zz
X . . . (2)
0

c
(c) y (d ) en (2)
DP 2 HP ( DP 2 DC 2 ) Z
X
4 HP
Pero,
DP 2 HP ( DP 2 DC 2 ) Z
A( Z ) X
2

4 HP
( DP DC )
2 2

A( Z ) DP 2
Z . . . (3)
4 HP
 Para el cilindro (Z > Hp)

Ecuacin del Cilindro

3.1416
A( Z ) DC
2
* 2.5 4.90873852 m . . . (4)
2 2

4 4
Clculo de H y V para t0 = 0 s Condicion de flujo permanente

Se determina el nivel de referencia, el cual debe estar ubicado en el

punto ms bajo del sistema. Por conveniencia se escoge el punto
ms bajo de la tubera.
 Hallando la altura inicial (H0):

H HT L sen S sen sen sen sen sen

0 1 0 2 1 3 2 5 4
H 10 600 * Seno(20) 80 * ( Seno(80) Seno(30) Seno(10) Seno(60) Seno(5))
0

H 382.414425 m
0

Hallando la altura final (Hf): Hf 0 m

Hallando la altura neta (H):

H H 0 H f 382.414425 m
Hallando el coeficiente de prdida (KFL)
4

f L D
4
f L 0.025 * 600 2 0.02 * 300
K
FL
1
1
2
*2

2
3.45695246 m
D D
1
D 1
3.5 3.5
2
2

Hallando V para la tubera forzada:

2 g H 2 * 9.81* 382.4144
V 41.02963 m / s
1 K 1 3.4567
2
FL
Hallando V para el canal de aduccin por ecuacin de continuidad

2
D
2
2
V V
1
2
* 41.02963 13.3974312 m / s
2

D 1 3.5
Hallando las alturas (H) para cada punto sealado en el canal de
aduccin

Punto 0 H 0t H 0 382.414425 m
0

Punto 1 H1 H 0
t0 t0

1 f 1
S1 V1 2
382.414425
1* 0.025 *100 *13.3974312 2
375.87987 m

2 g D1 2 * 9.81* 3.5

H H
2 f 1
S V
1 1
2

369.345314 m
Punto 2
t0 t0

2 g D
2 0
1
 Punto 3
H3 H0
t0 t0

3 f 1
S1 V1 2

362.810759 m
2 g D1
Punto 4ED
H t0
H t0

4 f 1
S V
1 1
2

356.276204 m
2 g D
4 ED 0
1

Hallando H en el disipador

Punto 4D H 4t D H 4t ED 356.276204 m
0 0

Hallando H en cada punto sealado en la tubera forzada:

Punto 4SD H t0
4 SD
H 356.276204 m
t0
4D
Punto 5
H H t0 t0
1 f S V
2 2 2
2

287.634785 m
2 g D
5 4D
2

Punto 6
H H
t0 t0

2 f 2
S V 2 2
2

218.993367 m
2 g D
6 4D
2

Punto 7
H 7t H 4t D
0 0

.3 f 2
S 2
V2
2

150.351948 m
2 g D2
Punto 8
H H
t0 t0

.4 f 2
S V 2 2
2

81.7105292 m
2 g D
8 4D
2

Punto 9 H H t0 t0 5 f S V
2 2 2
2

13.0691105 m
2 g D
9 4D
2
Clculo de H y V para t > t0 usando el MTODO DE LAS CARACTERSTICAS

Caso TC = 4 TR

Se procede a calcular el TC

TC 4 * TR 4 *1.6 6.4 s

Ahora se procede a usar el mtodo de las caractersticas en cada punto del

sistema acoplado para el tiempo
t1=t0+t.
Para el punto 0 (Entrada al canal de aduccin )

c t 1200 * 0.08
H H H H
to 1 to to
382.414425 * 382.414425 375.87987 376.141252 m
S
S 0 0 1
1
100

c t 1200 * 0.08
V V to 1
V V to to
13.3974312 * 13.3974312 13.3974312 13.3974312 m / s
S 100
S 0 0 1
1

C H Z G R V
S S 1 0 1 1 S V
S

C 376.141252 122.324159 0.02736161 0.03494976 *13.3974312 *13.3974312 1256.78166 m
s

H H 382.414425 m
t1
0
t0
0

C H t1

V
0
t1
13.400428 m / s
S 0

Z 1
 Para el punto 1 c1 t
H R H1 to
H 1t H 0t 376.141252 m
o o

S1
c1 t
VR V1
t o
V1 t V0t 13.3974312 m / s
o o

S1
C H Z G R V
R R 1 0 1 1 R V R
2014.3428 m
c1 t
H S H1 to
H 1 H 2 369.606697 m
t t o o

S1
c1 t
VS V1
to
V1 t V2t 13.3974312 m / s
o o

S1 C C
H t1
375.513294 m S R

V
1

C H Z G R V 1263.31621 m 2
C C
S S 1 1 2 1 S S

V t1
13.3974312 m / s
S R

2 Z
1
1
 Para el punto 2

c t
H H H H 375.618487 m
to 1 to to

S
R 2 2 1

C H Z 1 G R V V 2007.8082 m
c t
R R 1 1 2 1 R R

V V V V 13.3974312 m / s
to 1 to to

S
R 2 2 1
1

c t
H H H H to 1 to to
363.072141 m
S
S 2 2 3
1
S S
C H Z G R V V 1269.850077 m
1 2 3 1 S S

c t
V V V V
to 1 to to
13.3974312 m / s
S
S 2 2 3
1

C C CS CR
H t1
368.978739 m S R
V2
t1
13.3974312 m / s
2 Z1
2
2
 Para el punto 3
c t c t
H H H H
to 1 to to
369.083932 m V V V V
to 1 to to
13.3974312 m / s
S
R 3 3 2
S
R 3 3 2
1
1

C H Z G R V
R R 1 2 3 1 R VR
2001.2737 m

c t c t
H H H H
to 1 to to
356.537586 m V V to 1
V V to to
13.3974312 m / s
S S
S 3 3 4 S 3 3 4

1 1

C H Z G
S S 1 3 4 ED
R V 1 S V
S
1276.38532 m

C C C C
H t1
362.444184 m S R
V t1
13.3974312 m / s
S R

2 Z
3
2
3

1
 Para los punto 4ED, 4D y 4SD
c t
H H H H 362.549377 m
to 1 to to

S
R 4 ED 4 ED 3
1

c t
V V V V 13.3974312 m / s
to 1 to to

S
R 4 ED 4 ED 3
1

C H Z G
R R 1 3 4 ED
R V 1 R V R
1994.7391 m

c t
H H to
H
2 to
H to
287.634785 m
S
S 4 SD 4 SD 5
2

c t
V V to
2
V to
V to
41.0296331 m
S
S 4 SD 4 SD 5
2
C H Z G
S S 2 4 SD 5
R V2 S V
S
3829.38577 m

Sea:
Z 2 S 2 sen1 sen 3 sen 2 sen 5 sen 4 167.202339 m

Z t H 0t H 4t D 26.14 m
0
0 0

Zpt0 Z t0 HP HC 189.07 m
Qp t0 0 m / s 3

Entonces:

Z t1 H H t0
0
t0
4D 26.14 m
Como Z t1 HP AD t1 16.0596 m 2 ( Est dentro de cilindro )
2 AD A A
Zp Qp C Z C Z
t1 1 2

t Z Z
t0 t0 R 2 S 2

Zp 1 2
189.07 m
t1
A A 2 AD
1 2 t1

Z Z t 1 2

2 ADt1
Qp t1
t

Z t1 Z t0 -Qp t0 0.1284 m 3 / s

H t1
4 ED
HD Ht1
4D
t1
4 SD
Zp Z 356.27588 m
t1 2
 CR Zpt Z 2
V4 ED t1 1
13.39444 m / s
2 Z1
Qp
V t1 t1
0.00800 m / s
AD
4D
t1

C Zp Z
V t1
S t1 2
41.06134 m / s
2 Z
4 SD
2
 Para el punto 5
c t
H H to 2
H H to to
356.276204 m
S
R 5 5 4 SD
2

c t
V V V V
to 2 to to
41.0296331.m / s
S
R 5 5 4 SD
2

C H Z G
R R 2 4 SD 5
R V 2 R V
R
4466.8318 m

c t
H H H H to 2 to to
218.993367 m
S
S 5 5 6
2

C C
c t H 286.303522 m
V V 41.0296331 m / s
t1 S R

V V to 2 to to 5
2
S
S 5 5 6

2
C C
V V 41.0109819 m / s
t1 S R

C H Z G R V
S S 2 56 2 S S
3894.22471 m 5
2 Z 2
Para el punto 6

c t
H H H H
to 2 to to
287.634785 m
S
R 6 6 5
2

c2 t
VR V6 V6t V5t 41.0296331 m / s
to o o

S 2
C H Z G R V
R R 2 5 6 2 R V R
4401.9928 m

c t
H H to 2
H H to to
150.351948 m
S
S 6 6 7
2
C C
H
t1
6
S
218.457763 m
R

c t 2
V V V V
to 2 to to
41.0296331 m / s
S
S 6 6 7
2
C C
V 41.0404789 m / s
V
t1 S R

C H Z G R V
S S 2 67 2 S S
3965.07729 m 6
2 Z 2
 Para el punto 7
c t
H H to 2
H H to to
218.993367 m
S
R 7 7 6
2

c t
V V to 2
V V to to
41.0296331 m / s
S
R 7 7 6
2

C H Z G R V
R R 2 6 7 2 R V R
4331.1402 m

c t
H H H H 81.7105292 m
to 2 to to

S
S 7 7 8
2

c t
V V V V 41.0296331 m / s
to 2 to to
H
t1 C C
S
149.674394 m
R

S
S 7 7 8 7
2 2

C H Z G R V
S S 2 7 8 2 S V S
4031.79145 m
V7t 1
CS CR
41.0201799 m / s
2 Z2
 Para el punto 8
c t
H H H H 150.351948 m
to 2 to to

S
R 8 8 7
2

c t
V V V V 41.0296331 m / s
to 2 to to

S
R 8 8 7
2

C H Z G R V
R R 2 7 8 2 R V R
4264.4261 m

c2 t
H S H8 H 8t H 9t 13.0691105 m
to o o

S 2
C C
c2 t H
t1 S
80.4322597 m
R

VS V8 V8t V9t 41.0296331 m / s 2

t 8
o o o

S 2
CS CR
V8t 41.0449794 m / s

1

2 Z2
C H Z G R V V 4103.56156 m
S S 2 89 2 S S
Para el punto 9 (Localizacin de la vlvula)
c t
H H to
H H
2 to to
81.7105292 m
S
R 9 8 8
2

c2 t
VR V9 to
V9t V8t 41.0296331 m / s
o o

S 2
R R
C H Z G R V2 8 9 2 R V R
4192.656 m
Bn
t1
n 1 1 nf 0.96904236
TC

k
n V2 2
4.1377203 m / s
H
 Bk k Z 2 421.38349
Ck k C Hf 17331.484
R

H 9 CR Z 2 V8 344.706997 m
t1 t1

1/2
Bk2
Bk
V9t
1
Ck 37.7483794 m / s
4 2

Para los dems tiempos t > t1 se procede de manera anloga. Los resultados se
muestran en la tabla I. los grficos H vs. t y V vs. t se muestran
Determinar para el sistema acoplado las H y V utilice los datos del ejemplo, excepto para la altura del
reservorio, los dimetros de la parbola y del cilindro que forman el disipador

La altura del reservorio sern los dos ltimos dgitos del cdigo de su matricula, para aquellos
alumnos que tengan sus dgitos 00,selecciones el sexto y stimo digito que ser la altura

El diametro de la parbola ser la suma del segundo digito con el cuarto digito del cdigo
de su matricula

El diametro del cilindro es el ultimo digito, si fuera cero o 1, tomar el valor de 4 metros.

Dos procesos iterativos para flujo no permanente, cuadro de resultados, graficas y

programa si los hubiera realizado

Fecha de entrega: Mircoles 07 de diciembre horas 9.00 pm- williamchn_59@yahoo.es