HIDRULICA DE CONDUCTOS FIC
HIDRULICA DE CANALES
Flujo Uniforme
Walter La Madrid
[email protected]
Flujo Uniforme en canales abiertos
Flujo uniforme es una condicin de equilibrio de
flujo en un canal:
a) pendiente constante
b) seccin constante
c) rugosidad constante
d) tirante, rea mojada, la velocidad y la
descarga en cada seccin del canal son
constantes
e) pendiente del lecho del canal, de la
superficie del agua y la lnea energa son
paralelas, So = Sw = S
f) y1 = y2 , V1 = V2
Prof. Walter La Madrid
Flujo Uniforme en canales abiertos
Prof. Walter La Madrid
Ecuacin de Manning
Gauckler(1867)ManningStrickler (1923)
En 1889 Ingeniero irlands, Robert Manning present la
frmula:
2 1
3
AR S o 2
Q
n V velocidad del flujo (m/s)
2 1 Q caudal (m3/s) = VA
3
A( A / P) S o 2
n coeficiente de rugosidad de
Q Manning
n R Radio hidrulico (A/P)
5
3
1 So es la pendiente de fondo dado
A So 2
Q 2
como fraccin
P 3n
Prof. Walter La Madrid
Flujo Uniforme en canales abiertos
Valores de n para la Frmula de
Manning
Existen en la literatura especializada tablas
que relacionan los valores del coeficiente
de rugosidad n de la frmula de Manning,
con la natureza de las paredes (permetro
mojado) de los canales, tanto para
conductos naturales como artificiales
Prof. Walter La Madrid
Material Coeficiente de Manning Coef. Hazen- Williams
n CH
Asbesto cemento 0.011 140
Latn 0.011 135
Fierro fundido (nuevo) 0.012 130
Concreto (cimbra metlica) 0.011 140
Concreto (cimbra madera) 0.015 120
Concreto simple 0.013 135
Cobre 0.011 135
Acero galvanizado 0.016 120
Acero (esmaltado) 0.010 148
Acero (nuevo, sin recubrim.) 0.011 145
Acero (remachado) 0.019 110
Plomo 0.011 135
Plstico (PVC) 0.009 150
Madera (duelas) 0.012 120
Vidrio (laboratorio) 0.011 140
(Fuente: Computer Applications in Hydraulic Engineering, 5th Edition,
Haestad Methods)
Prof. Walter La Madrid
Flujo sobre una protuberancia
Qu pasa con el nivel de agua?
Pues resulta que
Izquierda = subcritico Derecho = supercritico
Derivando la Ecuacion de Bernoulli :
Prof. Walter La Madrid
Aplicando la Ecuacion de Bernoulli (p=0):
Para un canal de seccin rectangular:
Prof. Walter La Madrid
Sustituyendo Q = V z b en la Ecuacion de Bernoulli :
Para determinar la forma de la superficie libre, derivamos con
respecto a x:
Despejando dz / dx:
Prof. Walter La Madrid
(de la diapositiva anterior)
Como subcritico: Fr < 1
supercritico: Fr > 1
Flujo Subcritico con dh / dx > 1 dy / dx < 1
Flujo Supercritico con dh / dx > 1 dy / dx > 1
si el flujo es subcrtico si el flujo es supercrtico
Prof. Walter La Madrid
Ejemplo 1
El ancho de un canal rectangular horizontal se reduce
de 3,5 m a 2,5, y el piso se eleva 0,25 m de altitud
en una seccin determinada. En la seccin de aguas
arriba, la profundidad de flujo es 2,0 m y el factor de
correccin de la energa cintica a es 1,15. Si la
cada en la elevacin de la superficie del agua en la
contraccin es de 0,20 m, calcule la descarga si :
(a) la prdida de energa es despreciable, y
(b) la prdida de energa es una dcima parte de la
altura de velocidad en sentido ascendente.
[El factor de correccin de la energa cintica en la
seccin contrada puede suponerse que es la unidad].
Prof. Walter La Madrid
Prof. Walter La Madrid
Referencia figura anterior
y1 2.0m
y2 2.0 0.25 0.20 1.55m
Por continuidad B1 y1V1 B2 y2V2
2.5 1.55
V1 V2 0.5536 V2
3.5 2.0
(a) Cuando no hay prdida de energa :
Aplicando ecuacin de la energa a las secciones 1 a 2,
V12 V22
Z1 y1 a1 Z1 Z y2 a 2
2g 2g
a1 1.15 y a 2 1.0
Prof. Walter La Madrid
... Cont.
V22 1.15V12 y1 y2 Z
2g
V22
2g
1 1.150.5536 2.00 1.55 0.25
2
V22
0.6476 0.2
2 9.81
V2 2.462 m s
Caudal
Q 2.5 1.55 2.462 9.54 m3 s
Prof. Walter La Madrid
(b) Cuando hay prdida de energa:
V12 V12
H L 0.1a1 0.115
2g 2g
Ecuacin de la energa:
V12 V22
Z1 y1 a1 Z1 Z y2 a 2 HL
2g 2g
V22 V12
a 2 a1 H L y1 y2 Z
2g 2g
V12
a 2 1.0, a1 1.15 y H L 0.115
2g
Sustituyendo
V22 V12 V12
1.15 0.115 2.00 1.55 0.25
2g 2g 2g
Prof. Walter La Madrid
De
V1 0.5536 V2
V22
2g
1 0.91.150.5536 0.2
2
0.6826 V22
0.2
2 9.81
2.397 m s
Caudal
Q 2.5 1.55 2.397 9.289 m3 s
Prof. Walter La Madrid
Cantidad de Movimiento en Canales
El principio de momento es un concepto importante en problemas
de flujo en el que necesitamos determinar fuerzas. El principio de
impulso de la segunda ley de Newton que la suma de las fuerzas F
externa sobre un cuerpo lquido o sistema S es igual a la tasa de
cambio del momento lineal mV de ese sistema. Puesto que los
smbolos F y V representan los vectores cambio en el momento
debe estar en la misma direccin que la fuerza.
Fx=P1A1P2A2=Q(V2-V1)
Prof. Walter La Madrid
Fx=P1A1P2A2=Q(V2-V1)
P es igual a y
A es el area
es la densidad del fluido
Q es el caudal
V es la velocidad
Fx= y1A1 - y2A2 =Q(V2-V1)
Prof. Walter La Madrid
Ecuacin de Momento
Flujo Uniforme
Momentum es una magnitud vectorial. La ecuacin de
momento de uso comn en la mayora de los problemas
de flujo de canales abiertos es la ecuacin lineal-
momento. Esta ecuacin establece que la suma
algebraica de todas las fuerzas externas que actan en
una direccin dada en una masa de fluido es igual a la
tasa de tiempo de cambio de lineal-momento de la masa
de fluido en esa direccin. En un flujo constante la
velocidad de cambio de momentum en una direccin
dada ser igual al flujo neto de momentum en esa
direccin.
Prof. Walter La Madrid
La figura 1.15 muestra un volumen de control (un
volumen fijo en el espacio) delimitado por las
secciones 1 y 2, el fondo y una superficie situada
por encima de la superficie libre.
Prof. Walter La Madrid
Las diversas fuerzas que actan sobre el volumen de
control en la direccin longitudinal son :
(i) Fuerzas de Presin que actuan en el volumen de
control, F1 y F2 .
(ii) Fuerza tangencial en el fondo, F3 ,
(iii) La fuerza de la masa, es decir, la componente
del peso del fluido en la direccin longitudinal, F4 .
Por la ecuacin lineal-momento en la direccin
longitudinal para una descarga de flujo estable
de Q ,
F F F
1 2 F3 F4 M 2 M1
Prof. Walter La Madrid
El M1 1 QV1 = momentum de flujo entrando al
volumen de control M 2 2 QV2 = momentum flujo
saliendo del volumen de control.
En las aplicaciones prcticas de la ecuacin de momento, la
correcta identificacin de la geometra del volumen de
control y las distintas fuerzas que actan sobre l son muy
importantes. La ecuacin de momento es una herramienta
particularmente til en el anlisis de flujo rpidamente
variado (RVF) situaciones en las que las prdidas de energa
son complejos y no puede calcularse fcilmente. Tambin es
muy til en la estimacin de las fuerzas en una masa fluida.
Prof. Walter La Madrid
EJEMPLO 2 Estimar la fuerza en la
compuerta que se muestra en la figura. 1.16
Prof. Walter La Madrid
Solucin Considere un canal de ancho unitario. La fuerza
ejercida sobre la compuerta por el fluido es F ' , como se
muestra en la figura. Este es igual y opuesta a la fuerza
ejercida por el fluido sobre la puerta, F .
Considere el volumen de control, como se muestra por
lneas de puntos en la figura. La Seccin 1 esta lo
suficientemente lejos de la seccin de flujo de salida y la
distribucin de la presin se supone hidrosttica. La fuerza
de friccin en el fondo entre las secciones 1 y 2 se
desprecia. La Seccin 2 se encuentra en la vena contrada
del chorro donde las lneas de corriente son paralelos al
fondo.
1 2 1.0 Coeficientes de Boussinesq
Distribucin uniforme de las
velocidades en las secciones
Prof. Walter La Madrid
Las fuerzas que actan sobre el volumen de control
en la direccin longitudinal son :
F1
= Fuerza de presin en el volumen de control
1 2
seccin 11' y1
F2 2
= Fuerza de presin en el volumen de control
1 2
seccin 22' y2 actua en una direccin
2
opuesta
F3
= Fuerza de reaccin en la compuerta en la
seccin 33'
Por la ecuacin de momento,
y1 y2 F qV2 V1
1 2 1 2
2 2
Prof. Walter La Madrid
Donde q = caudal por unidad de ancho V1 y1 V2 y2
Simplificando la ecuacin,
1 y1 y2 2q 2
F y1 y2 y1 y2
2 y1 y2 g
Si las prdidas de energa entre las secciones 1 y 2
se supone que es insignificante, por la ec. de la energa
a1 a 2 1.0
V12 V22
y1 y2
2g 2g
Sustituyendo
q q
V1 y V2
y1 y2
Prof. Walter La Madrid
q2 2 y12 y22
g y1 y2
y
1 y1 y2
3
F
2 y1 y2
La fuerza en la compuerta F ' sera igual y
opuesta a F .
Prof. Walter La Madrid
Fuerza Especifica (Funcin Momentum)
La fuerza especfica o la funcin
momentum (F) se define como:
F o M = yA + Q2/gA
El diagrama de la fuerza especfica
esta en funcin de la profundidad del
canal (y en funcin de F). La funcin
de momentum mnimo (Fmin)
corresponde al flujo crtico (es decir,
el Fr= 1), y la profundidad del canal
correspondiente se conoce como la
profundidad crtica (yc). Para una
funcin momentum dado, dos
tirantes (ysub y ysup) son posibles,
y se denominan tirantes conjugados.
Las partes superior e inferior de la
curva de fuerza especfica
corresponden a flujos subcrtico y
supercrtico, respectivamente.
Prof. Walter La Madrid
Curva Fuerza Especifica (Funcin
Momentum)
Prof. Walter La Madrid
Curva Fuerza Especifica (Funcin
Momentum)
Prof. Walter La Madrid
ECUACIN DE CHEZY
Por definicin no hay aceleracin en flujo uniforme.
Aplicando la ecuacin de momento para un volumen de
control que abarca las secciones 1 y 2, la distancia L de
diferencia, como se muestra en la figura. 3.1,
P1 W sin Ff P2 M 2 M1
31
Prof. Walter La Madrid
Donde P1 y P2 son las fuerzas de presiones y M 1 y M 2
son los momentos del flujo en las secciones 1 y 2
respectivamente, W = peso del fluido en el volumen de
control y F f = fuerza cortante en el lmite.
Como el flujo es uniforme,
P1 P2 y M1 M 2
Tambien,
W AL y Ff 0 PL
Donde 0 = esfuerzo cortante promedio en el permetro
mojado de longitud P y = peso especfico del agua.
Remplazando sin por S 0 (= pendiente del fondo),
Se puede escribir como ALS0 0 PL
A
0 S0 RS 0
P
Prof. Walter La Madrid
Donde R A P se define como el radio hidrulico.
R es un parmetro de longitud que representa la forma del
canal. Desempea un papel muy importante en el desarrollo
de las ecuaciones de flujo que son comunes a todas las
formas de canales.
Expresando 2 el esfuerzo cortante promedio como,
0 kV donde k = un coeficiente que depende de la
naturaleza de los parmetros de superficie y de flujo, la
ecuacin. anterior se escribe como
kV 2 RS 0
dando lugar a
V C RS 0
Prof. Walter La Madrid
1
donde C = un coeficiente que depende de la
k
la naturaleza de la superficie y el flujo. La ecuacin se conoce
como la frmula de Chezy despus de que el ingeniero francs
Antoine Chzy, que se acredita con el desarrollo de esta relacin
bsica simple en 1769. Las dimensiones C son L
12
T 1
y se
puede hacer adimensional dividindolo por g . El coeficiente
es conocido como el coeficiente de Chezy.
Prof. Walter La Madrid
FORMULA DE MANNING
Una frmula de resistencia propuesto por Robert Manning,
ingeniero irlands, para flujo uniforme en canales abiertos, es
1 23 12
V R S0
n
donde n = coeficiente de rugosidad de Manning. Este
coeficiente es esencialmente una funcin de la naturaleza de la
superficie lmite. Cabe sealar que las dimensiones de n son
L1 3T . La ecuacin es conocida popularmente como la
frmula de Manning. Debido a su simplicidad y grado
aceptable de precisin en una variedad de aplicaciones
prcticas, esta frmula de la es probablemente la frmula de
flujo uniforme ms ampliamente utilizado en el mundo.
Comparando esta ecuacin con la frmula de Chezy, tenemos
1 16
C R
n
Prof. Walter La Madrid
RUGOSIDAD EQUIVALENTE
En algunos canales en el permetro pueden tener diferentes
rugosidades. Canales en el que se alinean solamente los lados,
canaletas de laboratorio con paredes de cristal y fondos speros,
ros con una cama de arena en la parte de aguas profundas y las
llanuras aluviales cubiertas de vegetacin, son algunos ejemplos
tpicos. Esta rugosidad equivalente, tambin llamada la
rugosidad compuesta, representa un valor medio ponderado
para el coeficiente de rugosidad. Existen varias frmulas para
calcular la rugosidad equivalente. Todos se basan en ciertas
suposiciones que son aproximadamente efectiva en la misma
medida. Uno de estos mtodos de clculo de la rugosidad
equivalente es la siguiente.
Prof. Walter La Madrid
Considere la posibilidad de un canal que tiene su permetro
formado por N tipos de rugosidades. P1 , P2 ,..., Pi ,..., PN son las
longitudes de estas N partes y
n1 , n2 ,..., ni ,...nN son los respectivos coeficientes de rugosidad.
Que cada parte Pi esta asociada con un rea parcial Ai de tal
forma que
N
A A A
i 1
i 1 2 ... AN A total area
Prof. Walter La Madrid
Se supone que la velocidad media en cada rea parcial es la
velocidad media para toda el rea de flujo, es decir,
V1 V2 ... Vi ...VN V
Por la frmula de Mannning
V1n1 V2 n2 Vi ni VN nN
S 12
0 2 3 2 3 ... 2 3 ... 2 3
R1 R2 Ri RN
Vn
23
R
Donde n = Rugosidad equivalente
23
Ai ni Pi 2 3
A nP 2 3
Prof. Walter La Madrid
ni3 2 Pi
Ai A 3 2
n P
n 32
P
A A A
i
n P
i
32
i
n
n 32
i P i
32
23
P
Esta ecuacin proporciona un medio para estimar la
rugosidad equivalente de un canal que tiene mltiples
tipos de rugosidad en su permetro.
Prof. Walter La Madrid
EJEMPLO Un canal trapezoidal de tierra (n = 0,025) tiene un
ancho de fondo de 5,0 m, pendientes laterales de 1,5
horizontal: 1 vertical y una profundidad de flujo uniforme de
1,1 m. En un estudio econmico para solucionar la filtracin
excesiva de los canales hay dos propuestas, a saber. (a)
revestir la lnea de los lados solamente y (b) revestir la lnea
de la cama solamente. Si el revestimiento es de concreto
liso (n = 0,012), determinar la rugosidad equivalente en los
dos casos anteriores.
solucin
Caso (a): revestimiento solo en los lados
Aqu para la cama
n1 0.025 y P1 5.0 m
Para los lados
n2 0.012 y P2 2 1.10 1 1.5
2
3.966 m
Prof. Walter La Madrid
P P1 P2 5.0 3.966 8.966 m
Rugosidad equivalente
n
5 0.025
1.5
3.966 0.012
1.5 2 3
8.96623
0.085447
0.020
4.31585
Caso b: Revestimiento solo del fondo
P1 5.0 m n1 0.012
P2 3.966 m n2 0.025
P 8.966
Prof. Walter La Madrid
n
5 0.012
1.5
3.966 0.025
1.5 2 3
8.96623
Rugosidad equivalente
0.079107
0.018
4.31585
Prof. Walter La Madrid
Seccin hidrulica ms eficiente
El transporte en una seccin de canal de un rea
dada aumenta con una disminucin en su
permetro. Por lo tanto una seccin de canal que
tiene el permetro mnimo para un rea dada de
flujo proporciona el valor mximo de conduccin.
Con la pendiente, coeficiente de rugosidad y el rea
de flujo fijo, un mnimo permetro de seccin
representar la seccin hidrulica eficiente, ya que
transporta el caudal mximo. Esta seccin del canal
tambin se llama la mejor seccin.
De todas las diferentes secciones de canal
abierto posibles, la forma semicircular tiene la
menor cantidad de permetro para un rea dada.
Prof. Walter La Madrid
Seccin hidrulica ms eficiente
Para una seccin transversal dada determinar la seccin
hidrulica mas eficiente.
1 Canal rectangular:
Canal Rectangular
hidrulicamente mas eficiente
es la mitad de un cuadrado.
Prof. Walter La Madrid
2 Seccin Trapezoidal
Para un rea dada de flujo
Diferenciar con referencia a
y suponiendo que A y m son
constantes.
Sustituyendo para el rea, la
ecuacin anterior se puede
reescribir como
Prof. Walter La Madrid
La mitad de la anchura de la parte superior = distancia
pendiente lateral (para pendiente lateral dado)
Sustituya este valor de b en la ecuacin A y P y
simplificando
Sustituya el valor de y en P
Prof. Walter La Madrid
Cual es el valor de m que hace que P sea menos?
Derivando con respecto a m e igualando a cero
Esto significa que la seccin es un medio
hexgono. Si un semi crculo se dibuja
con radio igual a profundidad y luego
lados de esta seccin son tangentes al
crculo.
Prof. Walter La Madrid
Seccin triangular: Si es el ngulo medio
R debe ser como mximo.
Prof. Walter La Madrid
Anchura de la superficie libre es
igual a la diagonal
Es un medio cuadrado sobre su vrtice y el ancho mximo es
igual a la diagonal.
Derivando alternativa para la seccin triangular:
Prof. Walter La Madrid
Resolviendo para
Prof. Walter La Madrid
SECCIN COMPUESTA
Algunas secciones de canal pueden estar formados como una
combinacin de secciones elementales. Normalmente los
canales naturales, como ros, tienen las llanuras de
inundacin, que son de ancho y poco profundo en
comparacin con el canal principal de profundidad. La Figura
representa una seccin simplificada de un arroyo con sus
bancos. Canales de este tipo se conocen como secciones
compuestas.
Prof. Walter La Madrid
Prof. Walter La Madrid
Considere la seccin compuesta que se divide en
subsecciones por lneas arbitrarias. Estos pueden ser
extensiones de los lmites de los canales profundos o
verticales de las lneas dibujadas en el borde de los canales
profundos. Es fcil ver que las subsecciones tendrn
diferentes velocidades medias dependiendo de la
profundidad y de la rugosidad de los lmites. Generalmente,
llanuras aluviales tienen mayor rugosidad que el canal
principal ms profundo.
Si la profundidad de flujo se limita slo (i.e. y h) al
canal profundo, el clculo de la descarga mediante el uso de
la frmula de Manning es muy simple.
Prof. Walter La Madrid
Sin embargo, cuando el flujo se desborda en la llanura de inundacin
(,i.e. y h) , el problema del clculo de descarga es complicada ya
que el clculo del radio hidrulico puede inducir a un pequeo error
de la seccin de flujo de conjunto y, por tanto, la descarga puede ser
subestimado. Esta subestimacin de la descarga ocurre en un
pequeo rango de y h y ym ,
, Donde ym = valor mximo de
y, del cual no se produce la subestimacin de la descarga como
antes. Para un valor de,
y ym el clculo de la descarga, considerando toda la seccin como
una unidad sera adecuada. Para valores de yen el rango de h y y m el
canal tiene que ser considerado que se compone de sub-reas y la
descarga en cada subzona determinada por separado.
Prof. Walter La Madrid
La descarga total se obtiene como la suma de las descargas a
travs de todas estas subreas. El valor de
ym dependera de la geometra del canal. Sin embargo, para
los propsitos prcticos, el siguiente mtodo de estimacin de
descarga puede ser adoptado.
(i) La descarga se calcula como la suma de la
descargas parciales en las sub-reas; por ejemplo
unidades 1,2 y 3 en Fig.321.
Qp Qi Vi Ai
(ii) La descarga tambin se calcula considerando
toda la seccin como una unidad, (parte
ABCDEFGH en la figura. 3.21), por ejemplo Qm
Prof. Walter La Madrid
(iii) Lamayor de las dos descargas anteriores, QP y
QW , es adoptada como la descarga a la profundidad y.
Para la determinacin de las descargas parciales Q y, por
i
tanto QP , en el paso (i) anterior, dos mtodos estn
disponibles.
El mtodo de Posey
En este mtodo, mientras que el clculo del permetro mojado
para las sub-reas, las divisiones imaginarias (FJ y CK en la fig.
3.21) se consideran como lmites para slo la parte ms
profunda y descuidaron por completo en el clculo
correspondiente a la parte menos profunda. De esta manera el
esfuerzo cortante que se produce en la interfaz de las partes
ms profundas y poco profundas es empricamente en cuenta.
Prof. Walter La Madrid
Mtodo Zero
Ejemplo Para el canal compuesto mostrado en la figura. 3.22
Determinar la descarga para una profundidad de flujo de (a)
1,20 m y (b) 1,60 m. Utilice el mtodo de Posey para el
clculo de las descargas parciales.
Prof. Walter La Madrid
Solucin
El caso (a):
(i) La zona de descargas parciales por el mtodo de Posey
Subarea 1:
A1 7.0 0.3 2.1 m
2
P1 0.3 7.0 7.3 m
R 2.1 7.3 0.288 m
2.1 0.288 0.0002
1
QP1
23 12
0.02
0.647 m s
3
Similarmente
QP3 0.647 m s
3
Sub-area 2:
A2 3.0 1.2 3.6 m2
Prof. Walter La Madrid
P2 3.0 1.2 1.2 5.4 m
R2 3.6 5.4 0.667 m
3.60.667 0.0002
1
QP 2
23 12
0.02
1.943 m3 s
QP total discharge by partial areas
0.647 0.647 1.943 3.237 m3 s
(ii) Por el mtodo de la Seccin Total
A 2.1 2.1 3.6 7.8 m2
P 0.3 7.0 0.9 3.0 0.9 7.0 0.3 19.4 m
R 7.8 19.4 0.402 m
Prof. Walter La Madrid
7.80.402 0.0002
1
QW
23 12
0.02
3.005 m3 s
Dado que QW QP , la descarga en el canal se toma como
Q QP 3.237 m s 3
Caso (b):
y0 1.6 m
(i) rea de descargas parciales por el mtodo de Posey
Subarea 1:
A1 7 0.7 4.9 m2
P1 0.7 7.0 7.7 m
R1 4.9 7.7 0.636 m
4.90.636 0.0002
1
QP1
23 12
Prof. Walter La Madrid
0.02
2.563 m s
3
Similarmente, Q 2.563 m3 s
P2
Sub-area 2: A 3 1.6 4.8 m2
2
P2 1.6 3.0 4.8 m
R2 4.8 6.2 0.636 m
4.8 0.774 0.0002
1
QP 2
23 12
0.02
2.862 m3 s
QP 2.563 2.563 2.862 7.988 m s
2
(ii) Metodo de la SeccinTotal
A 4.9 4.9 4.8 14.6 m2
Prof. Walter La Madrid
P 0.7 0.7 0.9 3.0 0.9 0.7 0.7 20.2 m
R 14.6 20.2 0.723 m
14.6 0.723 0.0002
1
QW
23 12
0.02
8.315 m3 s
Dado que QW QP , la descarga en el canal se
toma como
Q QW 8.315 m3 s
Prof. Walter La Madrid
BUENA SUERTE
63
64
