Profundidad Normal Canales Trapeciales

1
Solucin explicita para la profundidad normal en

canales trapezoidales y su implementacin.
Por: Jess A. Monroy y Martin Sortillon V.
Resumen: Sobre la base de la ecuacin adimensional para la
profundidad normal en canales trapezoidales propuesta por
Vatankhah [1] se obtiene una solucin explicita con un error menor
de 0.25%.
Introduccin
En [1] se presenta una ecuacin explicita para el calculo de la
profundidad normal con error menor al 0.25% en canales con
pendientes de talud menor a 760, esta pendiente no cubre en forma
amplia la seccin del trapecio, por esto, el motivo del estudio es
desarrollar una ecuacin que tienda a los 90 0 en el talud con el
objetivo de incluir la solucin de la seccin rectangular.
1) Ecuacin adimensional
La profundidad normal (yn) para un canal de seccin trapezoidal
usando la formula de Manning se expresa de la siguiente forma;
nQ (Byn zyn2 )5/3
S1/2
(B pyn )2/3
o
(1)
donde, Q = gasto, n el nmero
de Ma-nning, So = la pendiente
del fondo del cana, z1 y z2 las
pendientes del talud, z = la
pendiente
media,
p
=
permetro adimensional de los
taludes y = 1 para el sistema
SI y 1.49 para el sistema
ingles.
Figura 1) Variables de la seccin

transversal de un canal trapezoidal.
Segn Vatankhah/1, (1) se puede expresar de forma adimensional de

la siguiente forma;
t2 z
t*
2/5
1 0
(2)
donde, t = (2zyn/B + 1), z* =(2z/p 1) y
z(p/2) 2/3nQ
4
1/2 8/3
So B
3/5
(3)
Sobre la base de estas variables la ecuacin explicita que se propone

en [1] para la aproximacin de t en la 1 iteracin (t1) se basa en
aplicar el mtodo del punto fijo a la ec. (2) en dos iteraciones
anidadas de la siguiente forma;
t1 z* z
t *
0.4
0.4
1
(4)
donde el valor inicial de t se obtiene por el mtodo de regresin;
t z* 0.885z* 0.98 1 z* 1.05
0.033
0.465 11/8
(5)
y finalmente
yn t1 1
B
2z
(6)
El limite de aplicacin de (6) es para valores de z 0.25 y su error es
de 0.25%.
Para obtener la solucin de la seccin rectangular se (z = 0) se
propone el uso de la formula de Srivastava [4];
yn
0.6 1 2.4040.6231 1 2.030.9363

B
0.3929
2/5
(7)
1/ En la referencia [1] solo incluye la solucin de trapecios simtricos ( z1 = z2), en

esta investigacin se generaliza el resultado para trapecios no simtricos, por esto,
la frmula para las constantes z* y difieren entre el presente artculo y la
referencia [1] pero la ec. (2) es la misma.
donde = nQ/(So1/2B 8/3). El limite de aplicacin es para una relacin de

profundidad yn/B 100 y el error de 0.06%.
2) Solucin propuesta
Despejando t de la ec. (2) se obtiene;
z*
2/5
5/8
2/5
(8)
Al aplicar los limites: [t1, z*0], [t, z*0], [t1, z*-1], [t,z*-1]
a (8) se obtiene en forma prctica para estas 4 condiciones la
siguiente ecuacin:
t 1
5/8
(9)
En (7) se observa que para valores de -1< z* <0, en la variable se
debe multiplicar por una funcin, f(z*/t+1), por lo tanto: t = [f(z*/t + 1)
+ 1]5/8, considerando la siguiente funcin para f se obtiene una mejor
estimacin de (9):
az
c * b 1 1
( 1)

5/8
(10)
donde a, b y c son coeficientes empricos, siendo los valores originales;
b = 5/8, c = 1 y a = 1. Para valores grandes de el producto bc debe
ser aproximado al valor original de 5/8, por lo tanto, el valor inicial
que se asigna es: b = c = (5/8)1/2 = 0.791 = y dado que z * es negativo el
valor de a debe ser menor/2 a 1, sobre la base de estos acotamientos,
por ensayo y error estos coeficientes toman los siguientes valores; b
= 0.83, c = 0.77 y a = 0.9994, al sustituir estas cifras en (10) se
obtiene;
0.9994z*
(0.77 1)0.83

5/8
(11)
2/ La variable z* toma valores en el intervalo de: -1 z* 0, por lo tanto, t > 1 para
que el termino (z* + t)2/5 en la ec. (2) tenga una solucin real, por esto, el
coeficiente
a debe ser menor a 1.
Aplicando el mtodo de Newton-Raphson/3 a (2) se obtiene la 1

iteracin cuyo valor inicial t es el resultado obtenido en (11);
2
t(z
t)
2/5
1
*
t1 F(t) t
2t 0.4(z* t)3/5
(12)
y finalmente;
t 1
yn
1
B 2.0035z
(13)
El error que se obtiene en (13) al calcular yn/B es el siguiente:
Figura 2) Porcentaje de error en la
aproxi-macin de la profundidad
normal yn en un canal trapezoidal
usando la ec. (13) con limites de
aplicacin en la pendiente de talud
de: z = 10-8 para
la seccin
rectan-gular y z 0.01 para la

seccin trapezoidal y triangular.
El porcentaje de error se calcula
como;
yn
yn
B aprox / B exacta 1 100

La aplicacin prctica de la pendiente, z = 10-8, es aproximar el valor

de yn a un canal de seccin rectangular asumiendo que, {(z=10-8) =
(z=0)}, esto es, el rea de conduc-cin del trapecio es igual al rea
del rectngulo, la hiptesis, solo es vlida para relaciones yn/B 300
para asegurar que el error sea menor a 0.25%.
Conclusin
Los resultados indican, que usando la ecuacin (2) de Vantankah y el
mtodo de Newton-Raphson en la 1 iteracin
obtiene la
aproximacin de la profundidad normal yn/B con un error de 0.25% sin
limites de aplicacin prcticos/4 que incluyen la solucin de la seccin
rectangular, trapezoidal y triangular, lo cual elimina la necesidad de
3/ En vez de usar el mtodo del punto fijo como lo hace Vatankhah o Wang
[14], se opta por el de Newton-Raphson ya que este converge ms rpido a
la solucin, si se asegura su convergencia.
usar una formula diferente para el clculo de cada seccin (como las
ecuaciones 4 y 7) y esto, simplifica su implementacin como se
presenta en el Anexo 1.
Referencias
[1] Vatankhah AR, Explicit solutions for critical and normal depths in
trapezoidal and parabolic open channels, Ain Shams Eng J (2012),
http://dx.doi.org/10.1016/j.asej.2012.05.002
[2] Swamee PK. Normal-depth equations for irrigation canals. J Irrig Drain
Eng 1994;120(5): 9428.
[3] Swamee PK, Rathie PN. Exact solutions for normal depth problem. J
Hydraulic Res 2004; 42(5):5417.
[4] Srivastava R. Discussion of Exact solutions for normal depth problem
by Prabhata K. Swamee and Pushpa N. Rathie. J Hydraulic Res
2006;44(3):4278.
[5] Chaudhry M. H. Open-channel flow. New Jersey: Prentice Hall; 1993.
[6] Kouchakzadeh S, Vatankhah AR. Discussion of exact solutions for
normal depth problem by Prabhata K. Swamee and Pushpa N. Rathie.
Journal of Hydraulic Research 2007;45(4): 56771.
[7] Guerrero-Angulo, Jos Oscar; Arreguin-Corts, Felipe Ignacio; MorenoTrujillo, Marco Antonio. Ecuaciones explcitas para calcular el tirante normal
en canales rectangulares y trapeciales. Ingeniera Hidrulica en Mxico:
revista interdisciplinaria de ciencia y tecnologa del agua;20(3):117-23, Jul.Set. 2005. Ilus.
[8] Babaeyan-Koopaei, K. (2001). Dimensionless Curves for Normal-Depth
Calculations
in
Canal Sections.. J. Irrig. Drain Eng.,127(6), 386
389.doi:10.1061/(ASCE)07339437(2001)127:6
[9] Attari, J. and Foroghi, A. M. (1995). Explicit Formulae for Estimation of
Normal Depth in Channels. Proceedings, 26th Congr. Int. Assn. for Hydraul.
Res., Thomas Telford, London VoL 5, pp. 7-12.
[10] Shrestha, R. C. and Barkdoll, B. (2005), A Direct Solution for Normal
Depth in Open Channels, Proceedings EWRI of ASCE 'World Water and
Environmental Resources Congress', May 2005, Anchorage, Alaska, doi:
10.1061/40792(173)405.
[11] McLatchy, M. J. ( 1989). Newton-Raphson and Normal Depth in Open
Channel Flow. Channel Flow and Catchment Runoff: S.C. Yen (editor),
Proceedings of the Charlouesville Conference, Virginia, USA, 2, pp. 983-992.
[12] LMNO Engineering, Research, and Software, Ltd.
[13] Heastad Methods, Inc Products, FlowMaster 2005 (open cannel
hydraulic analysis software).
[14] Wang ZZ. Formula for calculating critical depth of trapezoidal open
channel. J Hydraul Eng 1998;124(1):901.
Nomenclatura
A
B
Fr
n
p
= rea de flujo = Byn + zyn2 [m2]

= ancho de fondo [m]
= nmero de Froude
= nmero de Manning [s/m 1/3]
= permetro mojado adimensional = (1 + z 12)0.5 + (1 + z22)0.5
4/ Para pendientes de talud z en el intervalo de [10-6, 10-3] el error aumenta

a 0.7%.
6
P
Q
R
So
t
T
yn
V
z
z1
z2
z*
= permetro mojado = B + pyn [m]

= descarga o gasto [m3/s]
= Radio hidrulico = A/P [m]
= pendiente de fondo de canal [sin unidades]
= profundidad normal adimensional = 2zyn/B + 1
= ancho superficial = B + 2zyn [m]
= profundidad normal [m]
= velocidad media = Q/A = R2/3So1/2/n [m/s]
= pendiente media de los taludes = (z1 + z2 )/2 [sin unidades]
= pendiente del talud izquierdo [sin unidades]
= pendiente del talud derecho [sin unidades]
= variable adimensional = 2z/p - 1
= descarga o gasto adimensional
Anexo 1) Implementacin de las formulas (11 y 12).

Figura 3) Implementacin de las
ecua-ciones (11, 12 y 14) en la
Hoja Excel/5.
Datos en unidades del SI
donde:
t2 = F(t1), o sea, la 2 iteracin, y
yn t 2 1
B
2z
(14)
Notas:
1) Para aproximar el clculo de la
seccin rectangular se debe
proporcionar: = 10-8 en las celdas
de talud/6, z1 y z2.
2) Para aproximar el clculo de la
seccin
triangular
se
debe
-5
proporcionar: = 10 en la celda
del ancho del canal/6 B.
5/ En la Figura 3 se observa, que por ms que se reduzca el tamao de las
ecuaciones explicitas (11 y 12) el clculo de las constantes y los resultado
hacen imprctica la solucin manual, por esto se recurre a la
implementacin de estas formulas en una escritura de cdigo en lenguajes
como Basic o Fortran (Chaudry [5]) o en una Hoja de Calculo.
6/ Si se usa la formula de Wang [14] para el clculo de la profundidad critica
y se desea aproximar la sec-cin rectangular o triangular a travs de esta
formula, las cifras: z = 10-8 y B = 10-5 deben proporcionar-se y el resultado
que se obtiene no tiene errores prcticos; z = 0 y B = 0 son puntos de
discontinuidad.
7
3) Solo es necesario la escritura
de (11) y (12), para el clculo de
t1, la ec. (12) se copia o arrastra

para el clculo de t2.
4) Al comparar el valor del gasto Q
en los Datos con el resultado VA
se mide el error/7 en m3/s.
Al incluir la 2 iteracin el error de yn/B disminuye de la siguiente

forma;
Figura 4) Porcentaje de error en la
aproxi-macin de la profundidad
normal yn/B
en un canal
trapezoidal usando la ecuacin (13)
con limites de aplicacin/8 para
pendientes de talud de: z 10-8,
con un error de 0.0005% y sin

los limites sealados en la nota /
4.
7/ Debido a los limites de aplicacin tan diversos, que por lo comn se
olvidan, en los clculos prcticos o finales de la variable yn, ya sea por
mltiples ecuaciones explicitas o mtodos numricos es convenien-te
verificar que tan exacta es la solucin, para el caso de yn, el valor de Q
versus VA indica la exactitud.
8/ El limite de aplicacin en la seccin rectangular sigue siendo la relacin
de profundidad yn/B 300 si este limite es superado, el valor de yn/B se debe
de multiplicar por: [(T/B 1)0.834 + 1] y el limite aumenta a: yn/B 450,000
con un error menor a 0.0005%. Donde T es el ancho superficial que se
obtiene en la Figura 3.
8
La 2 iteracin solo es necesaria
para reducir el porcentaje de error
a niveles de los programas
profesionales de clculos hidrulicos como el de Lmnoeng [12] o
Flow-Master [13].

Profundidad Normal Canales Trapeciales

Cargado por

Copyright:

Formatos disponibles

Profundidad Normal Canales Trapeciales

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Profundidad Normal Canales Trapeciales

Cargado por

Copyright:

Formatos disponibles

1

Solucin explicita para la profundidad normal en

nQ (Byn zyn2 )5/3

Figura 1) Variables de la seccin

Segn Vatankhah/1, (1) se puede expresar de forma adimensional de

Sobre la base de estas variables la ecuacin explicita que se propone

donde el valor inicial de t se obtiene por el mtodo de regresin;

t z* 0.885z* 0.98 1 z* 1.05

0.6 1 2.4040.6231 1 2.030.9363

1/ En la referencia [1] solo incluye la solucin de trapecios simtricos ( z1 = z2), en

donde = nQ/(So1/2B 8/3). El limite de aplicacin es para una relacin de

a debe ser menor a 1.

Aplicando el mtodo de Newton-Raphson/3 a (2) se obtiene la 1

rectan-gular y z 0.01 para la

B aprox / B exacta 1 100

La aplicacin prctica de la pendiente, z = 10-8, es aproximar el valor

= rea de flujo = Byn + zyn2 [m2]

4/ Para pendientes de talud z en el intervalo de [10-6, 10-3] el error aumenta

= permetro mojado = B + pyn [m]

Anexo 1) Implementacin de las formulas (11 y 12).

Datos en unidades del SI

t2 = F(t1), o sea, la 2 iteracin, y

t1, la ec. (12) se copia o arrastra

Al incluir la 2 iteracin el error de yn/B disminuye de la siguiente

con un error de 0.0005% y sin

También podría gustarte