CURSO : DRENAJE 5TO AO EAP INGCIVIL

UNJBG
CLASE 07: CAPITULO 06B - CLCULO DEL DRENAJE
SUBTERRANEO - ECUACION DE ERNST





Prof: F. Grnica T. Tacna, 06 de noviembre del 2012

Figura 13. Esquema de dren subterrneo

(REPASO CLASE ANTERIOR)
FORMULA DE DONANN
FORMULA DE DONANN
FORMULA DE DONANN
FORMULA DE HOOGHOUDT
FORMULA DE HOOGHOUDT

FORMULA DE HOOGHOUDT
FORMULA DE HOOGHOUDT

7.5.3 FORMULA DE ERNST

La ecuacin de Ernst, ofrece mejoras sobre las formulas
anteriores, ya que el lmite entre los dos estratos puede
estar por encima o por debajo del nivel de los drenes.

El principio fundamental de la solucin de Ernst , es el de
considerar tres componentes en el flujo:
-vertical,
-horizontal y
-radial (grficos 6.10 y 6.11).

Estos componentes dan un esquema bastante completo
del sistema total de flujo que esta basado directamente en
la analoga entre las leyes de Darcy y Ohm.

Grfico 6.10 - Esquematizacin del flujo de drenaje de acuerdo a
Ernst


L/2
superficiedel terreno
capa impermeable
L C
r
0
Flujo Vertical
Flujo
Radial
Flujo
Horizontal
K
1
K
2
3 Componentes de la perdida de carga
Este concepto implica que la prdida de carga hidrulica tambin tiene tres
componentes los cuales pueden ser calculados separadamente, tenindose
por superposicin la carga total Ah:


r h v
h h h h A + A + A = A
(6.15)


En donde:

Ah = prdida total de carga hidrulica (m).
Ah
v
= prdida de carga hidrulica debido a la componente vertical del flujo (m).
Ah
h
= prdida de carga hidrulica debido a la componente horizontal del flujo (m).
Ah
r
= Prdida de carga hidrulica debido a la componente radial del flujo.
Grfico 6.11 - Parmetros de la frmula de Ernst para los
componentes vertical y horizontal


L/2
R (m/da)
superficiedel terreno
capa impermeable
D
1
D
2
D
v
Ah
Ah/2
K
1
K
2
u
Parmetros de la frmula de Ernst para los
componentes vertical y horizontal

Componente vertical del flujo (Ah
v
):
La prdida de carga hidrulica debido al componente vertical del flujo, sigue la ley de Darcy:

1
v
v
K
D
. R h = A (6.16)
en el cual:
R = descarga normativa (m/da)
D
v
= distancia vertical entre el nivel del punto medio de la napa fretica entre los drenes y el
fondo de los drenes (m).
K
1
= conductividad hidrulica del estrato superficial (m/da)

La prdida de carga hidrulica debido al componente vertical del flujo puede ser de importancia,
cuando el suelo situado sobre el nivel de los drenes tiene relativamente una baja conductividad
hidrulica. En estos casos el flujo se dirige verticalmente hacia abajo a travs de esa capa y luego
se dirige hacia los drenes como flujo horizontal y radial. Cuando hay varios estratos con
conductividades hidrulicas diferentes situados sobre el nivel de los drenes, se puede adicionar las
prdidas de carga de cada estrato:

n . 1
Vn
2 . 1
2 V
1 . 1
1 V
1
V
K
D
.........
K
D
K
D
.( R
K
D . R
+ + + = (6.17)
Parmetros de la frmula de Ernst para los
componentes vertical y horizontal

Componente horizontal del flujo (Ah
h
):
La formula usada para el clculo de la prdida de carga hidrulica debido al componente
horizontal es la misma que la formula de Donnan, que Ernst presenta en la siguiente forma:
KD 8
RL
h
2
h
E
= A (6.18)
en donde:
n n 2 2 1 1
D K ......... D K D K KD + + + = E (6.19)
siendo:
R = descarga normativa (m/da)
L = espaciamiento de los drenes (m)
D
1
, D
2
, ..... D
n
= espesor de los estratos con diferente conductividad hidrulica (m)
K
1
, K
2
, ..... K
n
= conductividad hidrulica de los diferentes estratos (m/da)

La ecuacin anterior, indica que cualquier nmero de capas u horizontes del perfil del suelo
pueden ser tomados en cuenta por el valor EKD. Sin embargo, el espesor total del acufero no es
tomado en cuenta cuando es mas de un cuarto (1/4) del distanciamiento de los drenes.
L
4
1
D= E (6.20)
En caso que ED sea mayor, se toma el valor ED = (1/4) L.
Componente radial del flujo (Ah
r
):
La prdida de carga hidrulica debida al componente radial, se puede expresar de la manera
siguiente:
W . L . R h = A (6.21)
en donde:
R = descarga normativa (m/da)
L = espaciamiento de los drenes (m)
W = resistencia radial (das/m)

El clculo de la resistencia radial W presenta algunos problemas porque depende de la ubicacin
de los drenes respecto al lmite entre dos estratos con diferente conductividad hidrulica y a la
relacin entre estas conductividades hidrulicas.

Sin embargo, debe anotarse primero que para el permetro mojado de los drenes, que es un valor
importante en la determinacin de la resistencia radial, Ernst usa un valor diferente al de
Hooghoudt. En vez de utilizar Hr
0
como Hooghoudt, Ernst utiliza el permetro mojado de la zanja
vertical de relleno (P) que es igual al ancho de la zanja b mas 4 veces el radio exterior del tubo
(4r
0
) (P = b + 4 r
0
).
En el caso de zanjas con talud, se utiliza
2
z 1 y 2 b P + + =
, donde b es la anchura del
fondo de la zanja, y es el tirante de agua en dren y z es el talud de la zanja
(horizontal/vertical).

Ernst lleg a este valor asumiendo que la conductividad del material de relleno en la zanja es muy
alta y que el agua est justamente por el tope de los drenes.

Respecto a la relacin entre la conductividad hidrulica de los diferentes estratos, se puede
distinguir tres casos, considerando en todos ellos la nomenclatura respectiva en todas las
ecuaciones de resistencia radial:
W (das/m);
K = conductividad hidrulica del suelo (m/da);
D
r
= espesor del estrato en que ocurre el flujo radial (m);
P = permetro mojado (m).
Grfico 6.12: K
1
= K
2
: Suelo homogneo (Caso 1)














Para la resistencia radial en este caso, Ernst encontr:

P
D
ln
K
1
W
r
H
= (6.22)
L/2
superficiedel terreno
capa impermeable
L C
D
r
Caso 1: K
1
=K
2
=K
K
Parmetros de la frmula de Ernst para los
componentes vertical y horizontal

Caso 2.- K1 << K2: El estrato superficial tiene una conductividad hidrulica mucho menor
que el estrato profundo.

El dren puede tener tres ubicaciones respecto al lmite de los estratos:

Caso 2 A: El dren esta ubicado por encima del lmite de los dos estratos (grfico 6.13)
Caso 2 B: El dren esta ubicado en el lmite entre los dos estratos (grfico 6.15)
Caso 2 C: El dren esta ubicado por debajo del lmite de los dos estratos (grfico 6.16).
Parmetros de la frmula de Ernst para los
componentes vertical y horizontal

Caso 2 A: Dren ubicado por encima del lmite de los dos estratos.

En el caso 2A, la conductividad hidrulica es muy importante y se distinguen dos situaciones:

- Si K
2
> 20 K
1
se tiene:
P
D 4
ln
K
1
W
r
1
H
= (6.23)

- Cuando K
2
< 20K
1
se tiene:
P
aD
ln
K
1
W
r
1
H
= (6.24)

El parmetro "a" se determina del baco a = f(D
2
/D
r
, K
2
/K
1
) que se muestra en el grfico 6.14.
Grfico 6.13

L/2
superficiedel terreno
capa impermeable
L
C
D
r
D
2
K
1
K
2
Lmite entre los
Estratos
Caso 2A: K
1
<< K
2
Grfico 6.14

Caso 2 B: Dren ubicado en el lmite entre los dos estratos.

En el caso 2B (el dren esta en el lmite y K
1
<<K
2
) se tiene:

b
D 4
ln
K
1
W
r
2
H H
= (6.25)

En donde "b" es el ancho de la zanja vertical abierta o rellenada.












Grfico 6.15
L/2
superficiedel terreno
capa impermeable
L C
D
r
D
2
K
1
K
2
Lmite entre los
Estratos
Caso 2B: K
1
<< K
2
b
Caso 2 C: Dren ubicado por debajo del lmite de los dos estratos.

En el caso 2C (dren bajo el lmite y K
1
<<K
2
) se tiene:

P
D
ln
K
1
W
r
2
H
= (6.26)














Grfico 6.16
L/2
superficiedel terreno
capa impermeable
L C
D
r
D
2
K
1
K
2
Lmite entre los
Estratos
Caso 2C: K
1
<< K
2
Caso 3.- K
1
>> K
2
: El estrato superficial tiene una conductividad hidrulica mucho mayor
que el estrato mas profundo.

Aqu tambin se tiene que distinguir tres ubicaciones respecto al lmite entre los estratos, en la
misma forma que en el caso 2:

Caso 3 A: El dren esta ubicado por encima del lmite de los dos estratos (grfico 6.17).
Caso 3 B: El dren esta ubicado en el lmite entre los dos estratos (grfico 6.18)
Caso 3 C: El dren esta ubicado por debajo del lmite de los dos estratos (grafico 6.19).
Caso 3 A y 3B: Dren ubicado por encima del lmite de los dos estratos y en el lmite de los
dos estratos.

En los casos 3A y 3B, la formula de Ernst no es muy adecuada. Es mejor utilizar la formula
de Hooghoudt.














Grfico 6.17
L/2
superficiedel terreno
capa impermeable
L C
K
1
K
2
Lmite entre los
Estratos
Caso 3A: K
1
>> K
2
Grfico 6.18

L/2
superficiedel terreno
capa impermeable
L C
K
1
K
2
Lmite entre los
Estratos
Caso 3B: K
1
>> K
2
Caso 3 C: Dren ubicado por debajo del lmite de
los dos estratos.

En el caso 3C se tiene:
P
D
ln
K
1
W
r
2
H
= (6.27)












Grfico 6.19
L/2
superficiedel terreno
capa impermeable
L C
D
r
D
2
K
1
K
2
Lmite entre los
Estratos
Caso 3C: K
1
>> K
2
FORMULA COMPLETA DE ERNST:
La formula completa a partir de la formula (6.15)
r h v
h h h h A + A + A = A ser:

W . L . R
KD 8
RL
K
RD
h
2
1
V
+
E
+ = A (6.28)

En la formula las normas de drenaje estn dadas por Ah y R. Los valores de K y D son obtenidos
de los estudios de campo y W se puede calcular. As, todos los parmetros son conocidos, excepto
el espaciamiento de los drenes L, que se puede calcular a partir de los otros.

En la ecuacin (6.28), se puede despejar L como en la solucin de la ecuacin bicuadrada,
despus de tener el valor numrico de los otros parmetros.

6.5.4 CONCLUSIONES SOBRE LAS FORMULAS DE REGIMEN
PERMANENTE

A mayor conductividad hidrulica, mayor
espaciamiento entre drenes.
En un suelo homogneo, a mayor profundidad de
drenes, corresponde un mayor espaciamiento.
Los drenes deben colocarse en los estratos de mayor
permeabilidad.
El movimiento del agua no afecta prcticamente a una
profundidad superior a cuarta parte del
espaciamiento (L/4).
Los espaciamientos calculados, deben disminuirse en
un 15 a 20% para una mayor garanta.

6.5.5 CAUDAL A ELIMINAR POR LOS DRENES

Una vez determinados el espaciamiento y la profundidad de los drenes, es preciso conocer el
caudal que han de transportar, con el objeto de dimensionar las caractersticas hidrulicas de
estos.

Los drenes son conducciones abiertas o cerradas cuyo caudal aumenta con la longitud. El caudal
en el extremo de un dren si su longitud es "l" (grfico 6.20) ser:

RLl Q= (6.29)













Grfico 6.20 - Caudal a eli minar
L L L
L/2 L/2
Dren colector
Drenes
l
6.6 REGIMEN VARIABLE: FORMULA DE GLOVER DUMN

Las formulas de rgimen variable consideran el movimiento de la capa
fretica, tanto para carga del acufero como durante la descarga, en
consecuencia existe una variacin en el nivel de esta.

Como anteriormente se indic, trataremos la situacin de carga instantnea
del acufero (por ejemplo, en riego por gravedad). Si la cantidad de agua
aportada expresada en altura es "R", la capa fretica se eleva en una altura
R/ y a continuacin comienza a descender. La formula de Glover Dumn
estudia esta situacin.

Las formulas de Donnan, Hooghoudt y Ernst anteriormente tratadas estn
basadas en una situacin de flujo con recarga permanente. Esto significa
que tanto la altura de la napa fretica como la velocidad del flujo que
ingresa a los drenes no cambia durante el proceso de drenaje. La recarga de
la napa, proveniente de las precipitacin y la descarga de los drenes es igual
en estas condiciones y ambas son indicadas con el mismo smbolo "R".
Aunque, estas condiciones no son exactamente satisfechas en la realidad, la
suposicin de una situacin de flujo permanente ha probado producir muy
satisfactoriamente resultados en climas hmedos, cuando las variaciones en
precipitacin no son extremas.

FORMULA DE GLOVER DUMN

Sin embargo, cuando nos enfrentamos con problemas de drenaje
conectados con irrigacin, la recarga es evidentemente no
permanente, ya que ocurre a intervalos que son determinados por la
frecuencia de riegos. En estas condiciones, puede ocurrir una notoria
elevacin de la napa fretica como resultado de la inevitable
precolacin de parte del agua aplicada en el riego, seguido por una
gradual declinacin de la napa en los perodos comprendidos entre
dos riegos consecutivos. En estos casos la recarga R y la velocidad de
descarga R
t
evidentemente no son iguales.

Una formula de drenaje desarrollado para estas condiciones deber
dar respuesta a la pregunta: Cunto descender la napa fretica y
cual ser la velocidad de descarga de los drenes como funcin del
distanciamiento de los drenes y de la recarga no permanente? Tal
formula ha sido derivado por numerosos autores.

La primera formula que aqu se muestra fue obtenida por los estudios
de Glover Dumn y presentada por Dumn en 1954:

FORMULA DE GLOVER DUMN


)
h
h
27 . 1 ln(
t
.
D . K
L
t
0
2 2

H = (6.30)
en la que:
L = espaciamiento de los drenes (m)
K = conductividad hidrulica (m/da)
D = espesor del estrato en donde hay flujo horizontal (m)
= espacio poroso drenable (fraccin)
h
0
= altura de la napa fretica sobre el fondo del dren antes del descenso (m)
h
t
= altura de la napa fretica sobre el fondo del dren despus del descenso (m)
t = tiempo en que debe ocurrir el descenso de la napa fretica (das)
Grfico 6.21 - Parmetros de la Ecuacin
de Glover Dumn

L/2
superficiedel terreno
capa impermeable
L C
D
0
D
h
0
h
t
(h
0
+h
t
)/4
Parmetros de la Ecuacin de Glover Dumn
D es definida como la altura desde la capa impermeable hasta el punto medio de la vertical que va
desde el nivel de los drenes hasta la altura promedio de la napa en el punto equidistante entre los
drenes, o sea:

)
4
h h
( D
2
2
h h
D D
t 0
0
t 0
0
+
+ =
+
+ = (6.31)

K, D y se miden; la combinacin de h
0
, h
t
y t constituye las normas de drenaje. El
espaciamiento L debe ser tal que despus de un riego la napa fretica descienda desde h
0
hasta h
t

en t das.

La derivacin de la formula de Glover se basa en la solucin de la llamada ecuacin del flujo de
calor:
2
2
x
y
.
KD
t
y
c
c

=
c
c
(6.32)

Glover asume que inicialmente, esto es al tiempo t = 0, existe una napa fretica plana, aunque el
agua sobre los drenes se mantiene al mismo nivel formando casi un rectngulo.
Parmetros de la Ecuacin de Glover Dumn
Las condiciones lmites para la solucin de la ecuacin pueden ahora ser especificadas:

Y = Y
0
para 0 < x < L y t = 0

Y = 0 para x = 0, x = L y t > 0

La solucin de la ecuacin diferencial planteada, sujeta a las condiciones arriba mencionadas, para
la altura de la napa fretica en el punto medio entre los drenes (Yt), aparece como:

=
=

(
(

H
=
n
5 , 3 , 1 n
L
KD
n (
0 t
2
t
2
2
e
n
1 4
Y Y (6.33)

Kraijenhoff van de Leur introdujo el "coeficiente de reservorio" j, en su formula, que es mas
completa y que no trata este captulo. Este coeficiente se puede usar en la formula de Glover de la
siguiente manera:

y por tanto:

=
=

H
=
n
5 , 3 , 1 n
)
j
t
. n (
0 t
2
e
n
1 4
Y Y (6.34)

Para depresiones apreciables de la napa fretica cuando x = L/2, por ejemplo si Y
t
/Y
0
es de 0.8, se
puede usar el primer trmino solamente, con una aproximacin suficiente.
Parmetros de la Ecuacin de Glover Dumn
La ecuacin (6.34) puede ser escrita en estos casos de la manera siguiente:

j
t
0 t
e
4
Y Y

H
= (6.35)
o en la forma:
j
t
)
Y
Y
27 . 1 ln(
t
0
= (6.36)

Entonces cuando:
KD
L
j
2
2
H

= y Y
0
= h
0
y Y
t
= h
t


La formula de Glover obtenido despejando "L", viene a ser:

)
h
h
27 . 1 ln(
t
.
D . K
L
t
0
2
2

H
= (6.37)

La ecuacin para la descarga R
t
como funcin del tiempo sera:

=
=

H
=
n
5 , 3 , 1 n
j
t
. n
2
t
2
e
j
1
.
8
. R R (6.38)
Parmetros de la Ecuacin de Glover Dumn
Para disminuciones significativas de la napa fretica puede usarse solo el primer trmino de la
serie, o sea n = 1, tenindose:

j
t
2
t
e
j
8
R R

H
= (6.39)
O de la forma:
j
t
)
R
R
j
8 . 0
ln(
t
= (6.40)

De las ecuaciones anteriores, puede despejarse R
t
en funcin de Y
t
y por sustitucin de Y
0
= R/
se tiene:

H
=
H
H

=

2
j
e 8
je 4
.
R
R
Y
Y
)
j
t
(
j
t
2
n
t
(6.41)

y sustituyendo j por
KD
L
2
2
H

y Y
t
por h
t
se obtiene:

Parmetros de la Ecuacin de Glover Dumn
2
t
t
L
KDh 2
R
H
= (6.42)

Donde R
t
es la velocidad de descarga de los drenes, normalmente expresados en m/da en el
tiempo t, cuando la napa fretica esta descendiendo.

En 1960, Glover present la nueva formula para el espaciamiento de drenes en situacin de flujo
no permanente.

Esta nueva formula fue obtenida asumiendo que el tiempo t = 0, la napa fretica tiene la forma de
una parbola de cuarto grado. Esta suposicin se acerca mas a la realidad que la situacin asumida
para la primera formula, en que se estableca que la napa era inicialmente plana.

La consecuencia de este cambio en la situacin inicial, es que en la ecuacin (6.30) anteriormente
presentada, el termino: 4/H = 1.27, viene a ser 1.16.

Entonces, la ltima formula de Glover se presenta entonces en la siguiente forma:

t .
h
h
16 . 1 ln .
KD
L
t
0
2
2
|
|
.
|

\
|

H
= (6.43)
Parmetros de la Ecuacin de Glover Dumn
Debe notarse de que la ecuacin para el calculo de R
t
debe de ser mantenida y que la relacin h
t
/R
t

no es afectada por el trmino de 4/H.

Las formulas de Glover, se basan en un flujo solamente horizontal. No toman en cuenta la
resistencia vertical y radial.

Sin embargo, la resistencia radial se puede incorporar en la formula de la misma manera que
en la formula de Hooghoudt. Entonces, se puede usar el espesor del estrato equivalente "d"
en lugar de la distancia entre el fondo del dren y la capa impermeable (D
0
). Esta d, depende
de L, r
0
y D y se puede usar las mismas formulas o tablas que para la formula de Hooghoudt.

En este caso, la D de la ecuacin es igual a:

4
h h
d D
0 t
+
+ =
(6.44)
Los parmetros ya son conocidos.
Parmetros de la Ecuacin de Glover Dumn

El coeficiente de reservorio "j" usado por Kraijenhoff Van
der Leur tiene un sentido fsico. Este coeficiente incorpora
las caractersticas fsicas del suelo y las condiciones de
drenaje. Estos factores juntos, dictan el movimiento del
agua en el suelo. La dimensin de "j" es das.

Cuando "j" es grande, la napa fretica baja lentamente,
porque la transmisibilidad del suelo es pequea o el espacio
poroso drenable y el espaciamiento de drenes son grandes.
Al contrario, cuando "j" es pequea, la napa fretica baja
rpidamente.

6.7 CALCULO DEL DIMETRO DE LOS DRENES


Una vez conocido el caudal a eliminar por los drenes, el paso siguiente es el clculo de dimetro y
pendiente. Las formulas a utilizar (Pizarro, 1978) son las siguientes:

Drenes lisos (cermicos, hormign, PVC liso):

211 . 0 368 . 0
i Q 193 . 0 d

= (6.45)

572 . 0 714 . 2
i d 89 Q = (6.46)

Drenes corrugados:

187 . 0 375 . 0
i Q 2557 . 0 d

= (6.47)

5 . 0 667 . 2
i d 38 Q = (6.48)
donde:
d = dimetro interior del dren
Q = caudal a eliminar en m3/s
i = gradiente hidrulico, sin dimensiones.
CALCULO DEL DIMETRO DE LOS DRENES
Por lo tanto, una vez conocido Q existen infinitas
combinaciones satisfactorias entre "d" e "i". En la
prctica "i" suele venir limitado por la topografa.

El dimetro as calculado debe aumentarse
aproximadamente en un 10% (lo que representa
un aumento del 30% aproximadamente del
caudal transportable), con objeto de tener un
margen de garanta frente a defectos en la
alineacin de los drenes y disminucin del
dimetro efectivo por colmatacin del dren.

EJEMPLO FORMULA DE ERNST