8.

1 FUNDAMENTOS TEÓRICOS

FIGURA 8.10. Cuadro

resumen y ejemplos gráficos
de los posibles perfiles de flujo
gradualmente variado en
canales prismáticos.
 8.1 FUNDAMENTOS TEÓRICOS
8.1.10 Métodos de Cálculo de Perfiles de Flujo Gradualmente
Variado.
8.1.10.1 El Método del Paso Directo. Este método es aplicable a
canales prismáticos.

FIGURA 8.11. Esquema

ilustrativo del Método del
Paso Directo.
 8.1 FUNDAMENTOS TEÓRICOS

Aplicando la ecuación de Bernoulli entre (1) y (2), se tiene:

H1  H1 2  H2 (8.36)

Donde H1-2 = hf1-2 representa la pérdida de carga por fricción

en el tramo comprendido entre las secciones (1) y (2).
2 2
 v1  v2
z1  y 1 cos θ 
2
 hf1  2  z 2  y 2 cos θ 
2

2g 2g
2 2
 v1  v2
z1  z 2  y 1 cos θ 
2
 hf1  2  y 2 cos θ 
2

2g 2g

HIDRÁULICA PROFESOR:
Ramiro Marbello Pérez
 8.1 FUNDAMENTOS TEÓRICOS
Haciendo z = z1 – z2, y sabiendo que
z  E1  hf1  2  E2
2
v
E1  y 1 cos2 θ  1 , se tiene :
2g

E1  E2  hf1  2  z

Además, hf1  2  S f x Entonces,

E1  E2  Sf x  z

HIDRÁULICA PROFESOR:
Ramiro Marbello Pérez
 8.1 FUNDAMENTOS TEÓRICOS
z
Por otra parte, de la Figura 8.11, es claro que senθ 
x
Luego, E1  E2  S f x  x sen
E1  E2  xS f  sen 
E1  E 2
 x  (8.37)
S f  sen θ

El tramo x es la diferencia entre la abscisa de la sección (2) y

la correspondiente a la sección (1). Entonces,
E1  E 2
x  x 2  x1 
S f  sen θ

HIDRÁULICA PROFESOR:
Ramiro Marbello Pérez
 8.1 FUNDAMENTOS TEÓRICOS
E1  E2
De donde, x 2  x1  (8.38)
S f  sen θ

una de las dos formas siguientes:

1ª Forma:
Sf1 2 
1
2

S f1  S f2 
Donde Sf1 y Sf2 son las pérdidas unitarias de carga por fricción,
calculadas en las secciones (1) y (2), respectivamente, con la
siguiente ecuación.
2
 nQ 
Sf   
23 
(8.39)
  A RH 
HIDRÁULICA PROFESOR:
Ramiro Marbello Pérez
 8.1 FUNDAMENTOS TEÓRICOS

2ª Forma: Haciendo

1
y  y 1  y 2  (8.40)
2
A
y calculando con ésta, los elementos A, P, RH 
P

Luego, se calcula Sf con la siguiente ecuación:

2
 nQ 
Sf   
23 
(8.41)
  A RH 

HIDRÁULICA PROFESOR:
Ramiro Marbello Pérez
 8.1 FUNDAMENTOS TEÓRICOS
Finalmente, las ecuaciones (8.37) y (8.38) quedan
transformadas de la siguiente manera:
E1  E 2
x  (8.42)
Sf  sen θ
x 2  x1  x (8.43)

E1  E 2
x2  x1  (8.44)
Sf  senθ

Para ver el procedimiento de este método se remite el lector al

Anexo A4.

HIDRÁULICA PROFESOR:
Ramiro Marbello Pérez
 8.1 FUNDAMENTOS TEÓRICOS
8.1.10.2 Método del Paso Estándar. Este método se aplica a
canales prismáticos y no - prismáticos.

FIGURA 8.12. Esquema ilustrativo del Método del Paso Estándar.

HIDRÁULICA PROFESOR:
Ramiro Marbello Pérez
 8.1 FUNDAMENTOS TEÓRICOS
A igual que con el método anterior, se parte de la ecuación de
Bernoulli, aplicada entre las secciones (1) y (2):
 v 12  v 2 2 
he  K   (8.46)
 2 g 

Donde he representa las pérdidas de carga localizadas, y se

calculan por medio de la siguiente ecuación:
H1  H2  hf1  2  he (8.45)

Para canales convergentes: 0  K  0.1

Para canales divergentes: 0  K  0.2
Para contracciones y expansiones bruscas: K  0.5
Para canales prismáticos y regulares: K  0
HIDRÁULICA PROFESOR:
Ramiro Marbello Pérez
 8.1 FUNDAMENTOS TEÓRICOS
En este caso K  0, y, volviendo a la ecuación (8.45), se tiene:
2 2
 v1  v2
z1  y 1 cos 2 θ   z 2  y 2 cos 2 θ   h f1  2  he (8.47)
2g 2g

Definiendo la cota o altura piezométrica del flujo como

Cp  z  y cos2 θ (8.48)

La ecuación (8.47) queda así:

2 2
 v1  v2
C p1   Cp2   hf1  2  he (8.49)
2g 2g

Además, h f1  2  Sf1  2 x (8.50)

HIDRÁULICA PROFESOR:
Ramiro Marbello Pérez
 8.1 FUNDAMENTOS TEÓRICOS
Donde Sf1- 2 se calcula por medio de una de las dos formas descritas en el método anterior, y la longitud

x del tramo se fija previamente.

Finalmente, se debe comprobar que se satisface la siguiente igualdad:
H2  hf1  2  he  H1 (8.51)

Para ver el procedimiento de este método, el lector debe remitirse al

anexo A4.

Como se pudo observar, los métodos del Paso Directo y del Paso
Estándar sirven para calcular perfiles de F.G.V., por diferencias finitas, sin
intervenir directamente sobre la Ecuación Diferencial (8.5) u (8.15).

HIDRÁULICA PROFESOR:
Ramiro Marbello Pérez
 8.1 FUNDAMENTOS TEÓRICOS
8.1.10.3 Método Numéricos Aplicados a la Integración de
la Ecuación Diferencial de F.G.V.
Estos métodos, a diferencia de los dos anteriores, sí se aplican
directamente a la solución de la Ecuación Diferencial (8.5) u
(8.15).
dy
Dado el P.V.I :  f t, y  ; a  t  b ; C.I. : y t  a   y a    (8.52)
dx

P.V.I. : problema del valor inicial. C.I. : condición inicial

El objetivo es encontrar la función desconocida y(t) que

satisface la ecuación diferencial (8.52).

HIDRÁULICA PROFESOR:
Ramiro Marbello Pérez
 8.1 FUNDAMENTOS TEÓRICOS
Este problema es similar al de la ecuación diferencial ordinaria,
no – lineal y de primer orden que resulta de predecir la
variación de la profundidad, y, del flujo gradualmente variado,
en la dirección del flujo, es decir, la dirección x, en un canal
abierto.

FIGURA 8.13. Esquema ilustrativo de la aplicación de Métodos Numéricos al

cálculo de perfiles de F.G.V.
HIDRÁULICA PROFESOR:
Ramiro Marbello Pérez
 8.1 FUNDAMENTOS TEÓRICOS
En este problema hidráulico, la ecuación diferencial a resolver
es la siguiente:
2
 nQ 
 
dy   A R 2/3   S 0
  H  (8.53)
dx  Q2 T
3
 cos 2
θ
gA
En donde S0, n, Q, , , θ y g son constantes conocidas y
positivas, excepto S0 cuyo signo depende del tipo del canal,
así:
S0 (-) : Para canales tipo M, C y S
S0 (+) : Para canales tipo A
S0 = 0 : Para canales tipo H
HIDRÁULICA PROFESOR:
Ramiro Marbello Pérez
 8.1 FUNDAMENTOS TEÓRICOS

Los demás parámetros hidráulicos: A, RH y T, dependen de la

profundidad, y, y ésta a su vez, de la posición x. Por lo tanto, la
ecuación (8.53) se puede esquematizar como la ecuación
(8.52), de la siguiente manera:

dy
 f x, y  ; x 0  x  x f ; C.I. : y x  x 0   y x 0   y 0 (8.54)
dx
2
 nQ 
 
  A R 2/3   S 0
Siendo f x, y    H  (8.55)
 Q2 T
3
 cos 2
θ
gA
HIDRÁULICA PROFESOR:
Ramiro Marbello Pérez
 8.1 FUNDAMENTOS TEÓRICOS
Aquí, el problema hidráulico es encontrar la solución de la ecuación
diferencial (8.53), es decir, hallar la función desconocida y(x). Dicha
función no será una expresión algebraica, sino una discretización
numérica, de buena aproximación.

Como quiera que la función desconocida (la solución al problema de la

ecuación diferencial (8.53)), y(x), no puede obtenerse analítica o
algebráicamente, se recurre a un método numérico de solución. Éste
puede ser el método de Euler, el método de Heun, o el método de Runge-
Kutta llamados métodos de un paso o de paso simple. También, puede
recurrirse a métodos de pasos múltiples, o multipasos, tales como el
método de Adams – Bashforth (de 2, 3, 4 ó 5 pasos) y el método de
Adams – Moulton (de 2, 3 ó 4 pasos), o algún método predictor –
corrector, como el método de Adams de 4º orden o el de Milne –
Simpson.

HIDRÁULICA PROFESOR:
Ramiro Marbello Pérez
 8.1 FUNDAMENTOS TEÓRICOS
Cualquiera sea el método designado, se trata de hallar una solución
aproximada de y(x), para valores específicos de x, así:

y x 0  , yx1  , y x 2  ,  y x i  ,  yx n  ; i  1, 2, 3,, n

Siendo n el número de puntos o segmentos en los cuales se divide
igualmente el intervalo
a  x  b,
así:
b a xf  x0
x  ó x  (8.56)
n n

Donde:
x: es el tamaño del paso.
N: número entero de segmentos iguales en que se divide el intervalo
x0  x  xf
HIDRÁULICA PROFESOR:
Ramiro Marbello Pérez
 8.1 FUNDAMENTOS TEÓRICOS
8.1.10.3.1 Método Runge – Kutta de Orden 4 (MRK4)
aplicado al Cálculo de Perfiles de F.G.V.
Como cualquier otro método de aproximación, el método MRK4
permite encontrar una aproximación de la función desconocida
y( t i), es decir:

wi  y t i   y i , para i  1, 2, 3,, n (8.57)

Si es una buena aproximación de y(t i) = yi, la suposición de

que el problema está bien planteado implica que:
f t i , w i   y t i   f t i , y t i  (8.58)
HIDRÁULICA PROFESOR:
Ramiro Marbello Pérez
 8.1 FUNDAMENTOS TEÓRICOS
Partiendo de la condición inicial conocida, C.I. : w 0   Ana log ía con el F.G.V.: y(x0 )  y 0

en t 0  a Ana logía con el F.G.V.: X 0  a

El MRK4 calcula la solución específica en cada paso o punto i =

1, 2, 3…,n, a partir de la solución específica anterior, así:
1
w i 1  w i  K1  2K 2  2K 3  K 4  Ana logía : y i 1  y 0  1 K1  2K 2  2K 3  K 4  (8.59)
6 6
Donde:
K 1  h f t i, wi  Ana logía : K1  x x i, y i  (8.60)
 h K   x K 
K2  h f  t i  , w i  1  Ana logía : K 2  x f  x i  , y i  1  (8.61)
 2 2  2 2
 h K   x K 
K3  h f  t i  , w i  2  Ana logía : K 3  x f  x i  , y i  2  (8.62)
 2 2   2 2 
K4  h f t i 1 , wi  K 3  Ana logía : K 4  x f x i 1 , y i  K 3  (8.63)
HIDRÁULICA PROFESOR:
Ramiro Marbello Pérez
 8.1 FUNDAMENTOS TEÓRICOS
A continuación, se presenta un diagrama de flujo del algoritmo
del MRK4.

La codificación del MRK4, en Lenguaje de Programación

BASIC, puede verse en el Anexo A2, (PROGRAMA P8).

Dos aplicaciones concretas del Método MRK4 pueden verse en

el Anexo A4.

HIDRÁULICA PROFESOR:
Ramiro Marbello Pérez
 INICIO

TIPO DE CANAL : TC
Q, n,  , g , S 0
TIPO DE SECCIÓN TRANSV. : ST
X 0, Y 0, X F, N

DX  X F
 X0  S
N

nQ Q 2
C1  C2 
 g
SUBRUTINA
  tan 1 S 0  C 3  cos 2  ELEMENTOS
GEOMÉTRICOS
X0 , Y 0

DESDE I = 1 HASTA N RETORNO

X = X0 + I  DX
Y = Y0
S
2
 C1 
   S0
 A R2/3 
 
F  DYDX  H

C2 T
 C3
A3

K1 = DX  F

Y = Y + 0.5  K1

S
2
 C1 
   S0
 A R2 /3 
 
F  DYDX  H

C2 T
 C3
A3

K2 = DX  F

Y = Y + 0.5  K2
S
2
 C1 
   S 0
 A RH
2/ 3

F  DYDX 
C2 T
 C3
A3
K3 = DX  F

Y = Y + K3

S
2
 C1 
   S 0
 A RH
2/ 3

  
F DYDX
C2  T
3
 C3
A
K4 = DX * F

Y = Y 0 + ( K1 + 2 K2 + 2 K3 + K4 ) / 6

X,Y

SIGUIENTE I

FIN