‘’PERDIDA
DE AGUA EN REPRESAS HECHAS SOBRE SUELOS ARENOSOS Y DESCRIPCIÓN DE SU COEFICIENTE DE PERMEABILIDAD Y REDES DE FLUJO DESCRITAS POR EL RECORRIDO DEL AGUA SOBRE EL SUELO ARENOSO’’
CURSO
:
DOCENTE :
MECÁNICA DE SUELOS I ING. MILLA VERGARA, Elio
ALUMNO:
CHINCHAY POMA, Gianina Mariela………………………………...09 Mariela………………………………...0902.0904. 02.0904. 301 GRANADOS PALMADERA, Jhonston…………………………………101.0904.381 DOMINGUEZ FALCON, Pedro Alejandro…………………………072.0709.253 GONZALES JAQUE, Juan Benito…………………………………………111.0904.429
PACCINI SÁNCHEZ, Jean Carlos………………………………………….111.0904.424 RAMIREZ VIERA, Ray Robinson…………………………………………..072.0709.258
FECHA:
03 de abril del 2014.
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
PROYECTO DE INVESTIGACIÓN INVESTIGACIÓN – MECÁNICA DE SUELOS I I. GENERALIDADES 1. TÍTULO: „‟Perdida „‟Perdida de agua en represas hechas sobre suelos arenosos y descripción de su coeficiente de permeabilidad y redes de flujo descritas por el recorrido del agua sobre el suelo arenoso‟‟ arenoso ‟‟
2. PERSONAL INVESTIGADOR: - CHINCHAY POMA, Gianina Mariela, Universidad Nacional de Ancash
-
-
-
-
-
“Santiago Antúnez de Mayolo” Facultad de Ingeniería Civil, 092.0904.399 GRANADOS PALMADERA, Jhonston, Universidad Nacional de Ancash “Santiago Antúnez de Mayolo” Facultad de Ingeniería Civil, 101.0904.381 DOMINGUEZ FALCON, Pedro Alejandro, Universidad Nacional de Ancash “Santiago Antúnez de Mayolo” Facultad de Ingeniería Civil, 072.0709.253 GONZALES JAQUE, Juan Benito, Universidad Nacional de Ancash “Santiago Antúnez de Mayolo” Facultad de Ingeniería Civil, 111.0904.429 PACCINI SANCHEZ Jean Carlos, Universidad Nacional de Ancash “Santiago Antúnez de Mayolo” Facultad de Ingeniería Civil, 111.0904.424 RAMIREZ VIERA, Ray Robinson, Universidad Nacional de Ancash “Santiago Antúnez de Mayolo” Facultad de Ingeniería Civil, 072.0709.258
3. TIPO DE INVESTIGACIÓN: 3.1 DE ACUERDO ACUERDO A LA ORIENTACIÓN ORIENTACIÓN APLICADA 3.2 DE ACUERDO ACUERDO A LA TÉCNICA TÉCNICA CONTRASTACIÓN CONTRASTACIÓN TEÓRICA 4. RÉGIMEN DE INVESTIGACIÓN INVESTIGACIÓN LIBRE
2
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
5. INSTITUCIÓN A LA QUE PERTENECE EL PROYECTO UNIVERSIDAD NACIONAL “SANTIAGO ANTÚNEZ DE MAYOLO” FACULTAD DE INGENIERÍA CIVIL ESCUELA DE INGENIERÍA CIVIL 6. LOCALIDAD E INSTITUCIÓN DONDE SE EJECUTARÁ EL PROYECTO PROVINCIA: Huaraz INSTITUCION: Universidad Nacional “Santiago Antúnez De Mayolo”.
II. PLAN DE INVESTIGACIÓN 1. ANTECEDENTES Y JUSTIFICACION: Hace solo 60 años aproximadamente los proyectos de presas y de estructuras de retención de agua hechas con suelos se basaban casi exclusivamente en reglas empíricas que los constructores se transmitían por tradición oral. Se adoptaban las obras que habían resistido satisfactoriamente los estragos a causa del tiempo y de las aguas, independientemente de la naturaleza de los materiales y de las características del terreno de cimentación. Con el nacimiento de la mecánica de suelos y el conocimiento de los materiales, que con esta se adquirió, ha sido posible analizar bajo un nuevo fulgor el comportamiento de las presas y de las estructuras de retención. Fue el francés Henry Darcy quien estableció las bases para un estudio racional de los problemas prácticos acerca de la infiltración del agua a través de los suelos. Darcy en el siglo XIX estudió en forma experimental el flujo del agua a través de un medio poroso y estableció la ley que se conoce con el nombre de ley de Darcy. Dicha ley se basa en las siguientes hipótesis, que condicionan su validez: · Medio continuo, es decir que los poros vacíos estén intercomunicados. · Medio isótropo. · Medio homogéneo. · Flujo del agua en régimen laminar Darcy demostró que el caudal Q es proporcional a la pérdida de carga e inversamente proporcional a la longitud del lecho de arena y proporcional al área de la sección y a un coeficiente que depende de las características del material. De esta manera estableció que: Q = KA(h1 – KA(h1 – h2)/L h2)/L
3
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
En donde K es un coeficiente que se ha denominado coeficiente de permeabilidad con unidades L/T. Esta ley es solo aplicable en la resolución de problemas en que el flujo del agua sea laminar. Es decir que el flujo presente un número de Reynolds inferior a 2000. El número de Reynolds es una relación adimensional entre fuerzas de inercia y fuerzas viscosas, esta relación establece que: R = V D r/ m Posteriormente a Darcy, el siguiente paso fundamental en el conocimiento fue dado por Ph. Forchheimer, quien demostró que la función carga hidráulica que gobierna un flujo en un medio poroso es una función armónica, es decir, que satisface la ecuación de Laplace. Forchheimer desarrolló a principios del siglo XX, las bases para el método gráfico que hoy se conoce con el nombre de método de las redes de flujo, que sigue siendo el arma más sencilla y poderosa de que el ingeniero dispone para la resolución práctica de los problemas que involucren el flujo de agua en suelos. El método de las redes de flujo, que es una solución gráfica de la ecuación de Laplace, fue popularizado a partir de 1937 y desde entonces se ha transformado en el procedimiento normal de trabajo para todos t odos los ingenieros.
2. LA REALIDAD PROBLEMATICA: Las redes de flujo son la combinación de líneas de flujo y líneas equipotenciales; se construyen para calcular el flujo de aguas subterráneas y tienen reglas de trazo específicas.Las dos familias de curvas son ortogonales sola para suelos isotrópicos. Los suelos anisotropicos necesitan transformarse para ser tratados como isotrópicos. La presión intersticial que existe en un suelo con frecuencia no es la que corresponde a las condiciones hidrostáticas, sino aquella creada por el flujo de agua a través de los poros del suelo, flujo que se lograra observar con demostraciones realizadas con anterioridad, a su vez se estudiara la permeabilidad de los suelos. Se dice que un material es permeable cuando tiene vacíos continuos, es decir la facultad con la que el agua pasa a través de los poros. La perdida del agua en represas, debido a la infiltración que se produce a través del suelo ya sea debido a las presiones hidrostáticas como el tipo de material sobre el cual esta construido, identificando el coeficiente de permeabilidad que posee el suelo y asi mismo identificar el comportamiento de un flujo de agua a través del suelo, observando las líneas de corriente.
4
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
2.1. PROBLEMA: La perdida del agua en represas, debido a la infiltración que se produce a través del suelo.
3. OBJETIVOS
Calcular el coeficiente coeficiente de permeabilidad de un un suelo suelo en específico.
Realizar
el trazado de las redes de flujo sea suelos isotrópicos o anisotropicos sea el caso del suelo a utilizar. Poder observar como el líquido ingresa al suelo y como este logra buscar los vacíos continuos para poder salir.
Como
influye el coeficiente de permeabilidad en el tránsito de los fluidos y como se calcula.
Poder observar de manera experimental las líneas de flujo por donde circula el fluido líquido y así poder analizar con mayor detalle las redes de flujo que se trazara.
Poder profundizar más los temas de propiedades hidráulicas de suelos.
4. MARCO TEÓRICO a) LEY DE DARCY: La Ley de Darcy describe, con base en experimentos de laboratorio, las laboratorio, las características del movimiento del agua del agua a través de un medio poroso. La expresión matemática de la Ley de Darcy es la siguiente:
Dónde: = gasto, descarga o caudal en m3/s. = longitud en metros de la muestra = una constante, actualmente conocida como coeficiente de permeabilidad de Darcy, variable en función del material de la muestra, en m/s. = área de la sección transversal de la muestra, en m2.
5
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
= altura, sobre el plano de referencia que alcanza el agua en un tubo colocado a la entrada de la capa filtrante. = altura, sobre el plano de referencia que alcanza el agua en un tubo colocado a la salida de la capa filtrante.
El agua, por relaciones de energía, circula de mayor a menor altura piezométrica. Tal y como se puede ver, la relación se trata del gradiente de alturas priezométricas o gradiente hidráulico y se observa que:
Por lo que adopta un valor negativo. Ello se puede expresar:
Donde h es la altura piezométrica y z la longitud recorrida. Generalizando a 3 dimensiones se obtiene que:
K es la conductividad hidráulica (permeabilidad) y se trata de un tensor simétrico diagonalizable a 3 direcciones principales:
6
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
Y se obtiene:
El agua se desplazará en la dirección donde haya más permeabilidad y esta a su vez indicará a qué velocidad se mueve el agua en condiciones unitarias de gradiente. En terrenos isótropos, las 3 permeabilidades principales serán idénticas. b) PERMEABILIDAD: La facilidad con que el suelo fluye a través de un material se le denomina con el nombre de permeabilidad y el parámetro que permite cuantificar este fenómeno se le llama coeficiente de permeabilidad (K), el cual se define como la velocidad de flujo producida por un gradiente hidráulico unitario. COEFICIENTE DE PERMEABILIDAD (K) Usándose como una medida de la resistencia al flujo ofrecida por el suelo; los factores que intervienen son:
La densidad del suelo La distribución del tamaño de las partículas La forma y orientación de partículas del suelo El grado de saturación / presencia del aire El tipo de cationes y el espesor de las capas absorbidas asociados con los minerales de arcilla (cuando están presentes) La viscosidad del agua del suelo, que varía con la temperatura.
El intervalo de los valores de K es muy amplio y se extiende desde 1000 m/seg en el caso de gravas a grano muy grueso hasta un valor insignificante en el caso de las arcillas. La tabla adjunta presenta un resumen del orden de magnitud del coeficiente de permeabilidad de varios tipos de suelo
7
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
DETERMINACION DEL COEFICIENTE DE PERMEABILIDAD: Puede hacerse de varias maneras. Puede obtenerse un valor aproximado en el laboratorio cuando un ensayo de permeabilidad a carga o altura constante o uno de carga variable; la confiabilidad de las pruebas de laboratorio son: a. Los relativos a la la obtención obtención de buenas buenas muestras representativas. b. Los relativos relativos con con la reproducción de las mediciones de laboratorio. laboratorio. c. Los que dependen dependen de la reproducción reproducción de las condiciones condiciones de campo. PRUEBA DE CARGA CONSTANTE:
La prueba de carga constante se usa para determinar coeficientes de permeabilidad (k) de suelos de granos gruesos tales como gravas y arenas con valores de k mayores a 10-4 m/s. Un arreglo típico de la prueba de permeabilidad bajo carga constante se muestra en la figura.
8
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
En este tipo de arreglo de laboratorio, el suministro de agua se ajusta de tal manera que la diferencia de carga entre la entrada y la salida permanece constante durante el periodo de prueba.
Después que se estableció una tasa constante de flujo el agua es recolectada en una probeta graduada durante cierto tiempo. El volumen de agua Q recolectada se expresa como:
Dónde: A = área de la sección transversal transversal de la muestra del suelo. t = duración de la recolección del agua. L=longitud del espécimen. i = gradiente hidráulica. Además:
9
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
Reemplazando se tiene:
UNASAM – FIC
( )
PRUEBA DE CARGA VARIABLE: Un arreglo típico de la prueba de permeabilidad bajo carga variable se muestra en la figura.
El agua una bureta fluye atreves del suelo. La diferencia inicial de carga hi, en el tiempo t iempo t=0 es registrada y se permite que el agua fluya a través de la muestra de suelo de manera que la diferencia final de carga en el tiempo t= t2 sea h2. La tasa de flujo q de agua, a través de la muestra en cualquier tiempo t se expresa por:
() Dónde: a=área de la sección transversal de la l a bureta A=área de la sección transversal transversal de la muestra del suelo Reordenando la ecuación anterior:
10
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
( ) Al integrar el lado izquierdo i zquierdo de la ecuación con límites de tiempo entre 0 y t y el lado derecho con límites de diferencia de carga entre h1 y h2, se obtiene:
O
c) REDES DE FLUJOS: Las trayectorias del flujo del agua a través de los suelos reales y las correspondientes presiones de poro son extremadamente complejas, debido a la manera errática en la que es probable que varíe de punto a punto y en diferentes direcciones la permeabilidad. Por lo tanto, los análisis exactos de problemas tan comunes, como el efecto de un sistema de desagüe o el flujo bajo una ataguía dentro de la excavación para la pila de un puente rara vez son posibles. Sin embargo, a pesar de las complejidades de los problemas reales, el ingeniero puede mejorar bastante su criterio con respecto a la filtración y sus efectos, estudiando el flujo en condiciones sencillas esquematizadas.
11
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
ECUACIONES HIDRODINAMICAS QUE RIGEN EL FLUJO DEL AGUA A TRAVES DE LOS SUELOS A continuación se presenta un tratamiento matemático que permitirá llegar en forma sencilla a las ecuaciones básicas que se utilizan para plantear teóricamente el problema del flujo de agua a través de los suelos. Considérese un pequeño paralelepípedo de una región de suelo a través de la que fluye el agua, de dimensiones dx, dy y dz, tal como se muestra en la figura:
Supóngase que la velocidad V con la que el agua pasa por el elemento posee tres componentes Vx, Vy y Vz y que estas son solo función de x, y y z respectivamente pero no del tiempo (suponiendo que se trata de un fluido permanente) ni de ninguna otra variable. Suponiendo también que estas componentes son funciones continuas que admiten cualquier orden de derivación necesario al razonamiento expuesto. Dadas las condiciones, si en las caras I las componentes de la velocidad son Vx, Vy y Vz, en las caras II las componentes de la velocidad serán, respectivamente: Vx + (Vx/x) dx Vy + (Vy/y) dy Vz + (Vz/z) dz Supóngase ahora que la porción de suelo a través de la que fluye el agua tiene sus vacíos saturados y que además las partículas que la conforman son incompresibles. Así, durante el flujo, la cantidad que
12
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
entra al elemento tiene que ser igual a la que sale. Por lo tanto, teniendo en cuenta que el caudal que pasa por una sección puede expresarse como el producto del área de la sección por la velocidad del flujo, podrá escribirse: Vxdydz + Vydxdz + Vzdxdy = (Vx + Vx/x dx)dydz + (Vy + Vy/y dy) dxdz + (Vz + Vz/z dz)dxdy
Donde el término del lado izquierdo representa el caudal que entra y el del lado derecho el que sale. Simplificándola se obtiene: Vx/x
dxdydz + Vy/y dxdydz + Vz/z dxdydz = o
De donde: Vx/x
+ Vy/y + Vz/z = o
Esta ecuación es de gran importancia en la teoría de flujo de agua y se conoce con el nombre de ecuación de la continuidad. Es importante recordar que esta ecuación es solo aplicable cuando se cumplen los supuestos anteriormente mencionados, los cuales son: · · · ·
Flujo permanente Suelo saturado El agua y las partículas son incompresibles El flujo no modifica la estructura del suelo
Ahora teniendo en cuenta la ley de Darcy, tenemos que la velocidad de flujo de agua a través del elemento es: Vx = -Kxh/x Vy = -Kyh/y Vz = -Kzh/z
En estas ecuaciones el elemento de suelo se considera anisótropo en lo referente a su permeabilidad K en la dirección de cada eje.
13
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
Introduciendo estas ecuaciones en la ecuación de continuidad, obtenemos: Kx 2h/x2 + Ky 2h/y2 + Kz 2h/z2 = 0 En algunos casos en los cuales la sección transversal (x-y) es mucho mayor que la altura, el problema de flujo podrá estudiarse bidimensionalmente, escribiéndose en una forma más simplificada la ecuación: Kx 2h/x2 + Ky 2h/y2 = 0 Es esta la ecuación fundamental para el análisis de un flujo bidimensional en una región de flujo dada. De encontrarse en un suelo isótropo en lo referente a la permeabilidad, es decir: Kx = Ky = K, podrá ser simplificada para obtener la ecuación de Laplace: 2h/x2
+ 2h/y2 = 2h = 0
Una función que satisface la ecuación de Laplace se le conoce como armónica. Ecuación que para ser aplicada requiere: · ·
Suelo isótropo, en lo relativo a su permeabilidad Flujo bidimensional
Estos dos limitantes no son un gran obstáculo ya que un flujo bidimensional se ajusta a la mayoría de los casos prácticos, en cuanto a la isotropía habrá que considerar que muchas de las estructuras de tierra a través de las que interesa estudiar el flujo se construyen compactando por capas, procedimiento que conduce a permeabilidades horizontales bastante mayores que las que se obtienen para el flujo en la dirección vertical. La ecuación general de Laplace está constituida por dos grupos de funciones que pueden ser representados dentro de la zona de flujo en estudio como dos familias de curvas ortogonales entre sí. La solución general que satisfaga las condiciones de frontera de una región de flujo específica constituirá la solución particular de la ecuación de Laplace para esta región específica.
14
UNASAM – FIC
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
Es conveniente, a partir de la gráfica, una expresión que proporcione el caudal que pasa a través del elemento en el tiempo dt. dq = Kx h/ h/x dydz + Ky h/ h/y dxdz + Kz h/ h/z dxdy
Si el suelo es isótropo en lo referente a la permeabilidad: dq = K( K(h/ h/x dydz + h/ h/y dxdz + h/ h/z dxdy)
En flujo bidimensional: dq = K( K(h/ h/x dy + h/ h/y dx)
SOLUCIÓN DE LA ECUACIÓN DE LAPLACE: Centrándose en el caso de flujo bidimensional, se observa la ecuación de Laplace: 2
h/ x2 +
2
h/ y2 =
2
h=0
y se define una función: f unción: f = -Kh + c Función, conocida como función de potencial, que satisface la ecuación de Laplace, cumpliendo que: 2
/ x2 +
2
/ y2 = 0
Así la función f(x, y)= constante es una solución de la ecuación de Laplace, solución que representa una infinidad de funciones, según el valor de la constante c que intervenga. Esta expresión puede representar una familia de curvas que se desarrollan en la región plana en la que ocurre el flujo, obteniéndose una curva específica de la familia para cada valor de la constante que se tome. Considérese ahora una función y(x, y)= constante llamada función de corriente y definida de modo que: Vx =
/ y
Vy =-
/ x
La función y así definida satisface la ecuación de Laplace, de manera que se cumple:
15
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
2
/ x2 +
2
UNASAM – FIC
/ y2 = 0
Si el conjunto de funciones y(x, y)= constante, se representaran también por una familia de curvas (y= constante) en la región de flujo. La familia y = constante es ortogonal a la familia f= constante. En un problema específico en el que haya unas condiciones de frontera fijas, la solución de la ecuación de Laplace constituida por las dos familias de curvas, más la exigencia de que estas familias satisfagan las condiciones de frontera existentes, produce en definitiva una solución única del problema considerado.
INTERPRETACIÓN FÍSICA Se tendrá entonces que en la curva f= constante, todos los puntos presentaran la misma carga hidráulica, h. Por esta razón es que estas curvas reciben el nombre de líneas equipotenciales. En cuanto a las curvas y= constante, se considerara la la trayectoria del agua que pasa por P(x, y), en este punto el agua presenta una velocidad V.
A lo largo de la curva se tendrá: Tanq = Vy/Vx = dy/dx De donde: Vydx – Vydx – Vxdy Vxdy = 0 Lo que puede ser expresado como: dx/ dx/ x + dy/ dy/ y = 0
16
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
La anterior expresión es precisamente la diferencia total de la función y de manera que se cumple a lo largo de la trayectoria tr ayectoria del agua que: dy = 0 Y por lo tanto: y = constante La familia de curvas y= constante está constituida por las trayectorias físicas y reales del agua a través de la región de flujo.
REDES DE FLUJO La red de flujo es una representación gráfica de la solución de la ecuación de Laplace para f y y con las condiciones de frontera existentes en el flujo. Está constituida por líneas equipotenciales separadas igualmente en f y por líneas de corriente igualmente separadas en y. Esta separación se conoce como canal de flujo o canal de corriente. Todas las intersecciones de la red son ortogonales. Propiedades de las redes de flujo: o
El caudal caudal que que fluye entre dos líneas consecutivas es el mismo por por unidad de ancho.
Ni las líneas equipotenciales pueden cortarse entre sí, dentro del medio fluido, ni las líneas de corriente pueden cortarse entre sí dentro del medio fluido. Se trata entonces de definir en cada caso las condiciones de frontera específicas del problema y trazar, cumpliendo con estas, las dos familias de curvas ortogonales, obteniendo así una verdadera imagen gráfica del problema, que si ha sido realizada con cuidado podrá ser lo suficientemente buena para los fines ingenieriles. Para el trazo de una red de flujo se tienen los siguientes pasos:
o
Dibujar los límites del dominio Fijar tentativamente 3 o 4 líneas de corriente. o Trazar tentativamente equipotenciales, ortogonales a las líneas de o corriente. o Ajustar. Comprobar la bondad bondad del ajuste si al trazar trazar las líneas diagonales o de los cuadros se obtienen también curvas suaves, formando una nueva red. o
CALCULO DEL CAUDAL Al trazar cualquier red de flujo se dibujan las equipotenciales de tal manera que la Dh sea la misma m isma y que el Dq entre dos líneas de corriente sea el mismo.
17
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
Se tendrá entonces que: Dq = KaDh/b Si nf es el número total de canales de la red y nc el número de caídas de potencial que hay en toda la zona de flujo, entonces podrá escribirse: Dq = q/nf
y Dh = h/nc
Donde q y h son el caudal unitario total y la carga total. A partir de lo anterior se puede llegar a que: que: q/nf = = Ka h/nc q = (nf /n /nc) (a/b) kh Puesto que q, k, h, nf y nc son constantes para una red de flujo dada, la relación a/b debe serlo también. Esta condición implica que se estén cumpliendo las dos condiciones iniciales (que la Dh sea la misma y que el Dq entre dos líneas de corriente sea el mismo). El término nf/nc depende únicamente de la forma de la región de flujo, se le conoce como factor de forma y se representa: Ff = nf/nc El cálculo de las presiones hidrodinámicas en el agua que se infiltra a través de la región de flujo, es una de las aplicaciones más útiles de una red de flujo. Por ejemplo, considérese el flujo de agua a través de un material permeable en el que se ha hincado una tablestaca. Se supone que la permeabilidad del suelo es la misma en todos los puntos y que es igual en todas las direcciones; además, la tablestaca y el manto de roca que está debajo del suelo se consideran completamente impermeables. Se acepta como válida la ley de Darcy, y asimismo que tanto el suelo como el agua son incompresibles.
18
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
El agua que entra al suelo aguas arriba del tablestacado, se mueve hacia la superficie del terreno aguas abajo siguiendo trayectorias curvas regulares, como AB, que se conocen con el nombre de líneas de flujo. La circulación es producida por la carga hidráulica h, que impulsa el agua de A a 13. Al moverse una partícula de agua de A hacia B produce un arrastre por fricción en las partículas del suelo; a su vez, este arrastre produce una presión de filtración en la estructura del suelo; la presión de filtración en cualquier punto actúa en la dirección de la línea flujo en ese punto. Debido a esta viscosidad, la carga hidráulica disminuye continuamente de aguas arriba a aguas abajo a lo largo de cada línea de flujo. En consecuencia, el nivel piezométrico en un punto C tiene un valor intermedio entre los de A y B. Entre los extremos de cualquier otra línea de flujo, como A‟B‟ la carga hidráulica es también h y existe un punto C‟, en el en el cual el nivel piezométrico es el mismo que en C. Una línea, como LM, que una puntos de igual nivel piezométrico se conoce como línea equipotencial. Si la permeabilidad es constante y la misma en todas direcciones, la teoría demuestra que las líneas equipotenciales deben ser perpendiculares a las líneas de flujo. Esta conclusión permite resolver problemas en los que interviene el movimiento del agua a través de medios porosos, utilizando un procedimiento gráfico, en el que las líneas de flujo fl ujo y las equipotenciales se dibujan por aproximaciones sucesivas, hasta que se satisfacen las relaciones geométricas necesarias. El diagrama resultante ejemplificado se conoce como red de flujo. El primer paso para construir una red de f lujo consiste en tomar nota de todas las condiciones de frontera que deben satisfacerse; es decir, determinar si se conocen anticipadamente algunas líneas de flujo o equipotenciales. Por ejemplo, la pared formada por el tablestacado mismo constituye una línea de flujo. El agua que entra en el suelo inmediatamente a la izquierda del tablestacado, se mueve verticalmente hacia abajo en dirección de la punta de la tablestacada; pasa a la derecha abajo de la misma, y sube verticalmente a lo largo del paramento de aguas abajo de la propia tablestaca. La superficie del estrato impermeable es también una línea de flujo. El agua que entra a la información infinitamente lejos a la izquierda, fluye a lo largo de esta superficie hasta que ha pasado infinitamente lejos a la derecha. Es evidente que estas dos líneas de flujo marcan las fronteras de la región de flujo. Todas las líneas restantes deben estar situadas entre ellas. Además, debe notarse que la superficie superficie del terreno aguas arriba es una línea equipotencial, porque el nivel del agua en cualquier piezómetro con su extremo inferior en la superficie del terreno, coincidiría con la superficie libre del agua en esa zona. También, la superficie del terreno aguas abajo es una línea equipotencial; el nivel piezométrico coincide con la superficie del agua aguas abajo. Todas las líneas equipotenciales restantes deben localizarse entre estas dos. Las condiciones de frontera del problema se resumen en la figura.
19
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
Están representadas por las líneas equipotenciales ab y cd, y por las líneas de flujo bec yfg. La construcción del resto de la red de flujo se comienza haciendo el croquis de un pequeño número de líneas de flujo, quizá solamente dos; cada línea de flujo empieza en ab, y termina en cd. Como ab y cd son líneas equipotenciales, las líneas de flujo deben interceptarlas en ángulo recto. Las líneas de flujo bosquejadas deben ser curvas suaves cuya forma vaya marcando una transición gradual de una línea de flujo de frontera (bec) a la otra (fg). El tanteo inicial puede parecerse al de la figura. En seguida, se hace el intento de dibujar líneas equipotenciales que cumplan con los requisitos del problema. Estas líneas son también curvas suaves, y deben cruzar a las líneas de flujo en ángulo recto. Además, para simplificar la interpretación de la red de flujo, la separación entre las líneas equipotenciales debe ser tal, que la caída de nivel piezornétrico sea la misma entre cada par de líneas equipotenciales sucesivas. Es también conveniente separar las líneas de flujo, de manera que el gasto en cada canal limitado por dos líneas sucesivas, sea el mismo. Estos dos requisitos pueden satisfacerse haciendo cada área limitada por dos líneas de flujo adyacentes y dos líneas equipotenciales adyacentes, aproximadamente equidimensional. Es decir, las distancias a y b (fig. d) deben ser iguales. Como recurso para juzgar si un área de lados curvos satisface este criterio, puede inscribirse un círculo en el área. Así, en la fig. d es evidente que el área P es razonablemente equidimensional, pero que el área Q no satisface los requisitos.
20
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
a) Red de flujo por debajo de una ataguía de tablestacas. b) Condiciones de frontera que debe satisfacer la red de flujo. c) Primeras líneas de tanteo para la construcción de la red de flujo. d) Primeras líneas equipotenciales de tanteo para la construcción de la red de flujo. El primer tanteo de un conjunto de líneas equipotenciales debe dibujarse haciendo un esfuerzo para que las intersecciones resulten en ángulo recto con las líneas de flujo tanto como sea posible, y subdividir el espacio en áreas que puedan diferir entre sí en tamaño pero en que cada uno sea equidimensional. Ordinariamente, el primer intento no resulta satisfactorio, pero el estudio del croquis sugerirá las modificaciones apropiadas tanto en las líneas de flujo como en las equipotenciales. Como la forma y posición de cada conjunto de líneas depende de las del otro, comúnmente es necesario hacer una serie de ajustes. Puede adquirirse una gran habilidad para dibujar redes de flujo con la
21
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
práctica y el estudio de redes bien dibujadas para varias condiciones de frontera. En la fig. 2.9, se muestran varias redes de flujo para problemas relativos a cimentaciones para obras de ingeniería. Cuando la red de flujo se ha refinado de manera que satisfaga las condiciones de frontera y los criterios geométricos, proporciona la misma información que daría una solución analítica rigurosa del mismo problema. En realidad, con frecuencia las redes de flujo pueden dibujarse fácilmente en problemas demasiado complicados o difíciles para su tratamiento analítico. Con la red de flujo completa puede determinarse la presión en el agua en cualquier punto de un material permeable. Las condiciones en el punto C, fig.a, servirán de ejemplo. De acuerdo con la red de flujo, una partícula de agua que siga la trayectoria AB (o cualquier otra línea de flujo), cruzara ocho espacios limitados por líneas equipotenciales sucesivas. Cada espacio representa una caída equipotencial Δh. Si Nd representa el número de caídas equipotenciales a lo largo de cualquier línea de flujo,
En la fig.a, Δh = 1/8 h. Cuando el agua llegue al punto C, la carga perdida es 6 Δh, o 6/8h. El nivel piezométrico en C es entonces 6/8 h abajo del nivel del agua aguas arriba, o 2/8 h del nivel aguas abajo. La carga piezométrica en C es, por lo tanto,
Y la presión del agua en C es:
La presión de poro en C con respecto al nivel de aguas abajo, es la presión disponible en C para impulsar el agua el resto de la distancia a B, y es igual a 2/8 h γ w. El gasto que pasa debajo del tablestacado por unidad de longitud puede calcularse fácilmente. Considérese el gasto Δq a través del área sombreada en la fig. 2.8a. De acuerdo con la ley de Darcy, el gasto es:
22
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
En la que A es el área de la sección transversal del canal de flujo. El canal tiene la anchura a y un espesor unitario en la dirección del muro. Por lo tanto,
El gradiente hidráulico a través del área sombreada es Δh/a. Sin embargo,
Si el número de canales de flujo es Nf, el gasto total q por unidad de longitud de muro es:
El examen de la red de flujo (fig. 2.8a) muestra que la filtración brota en puntos como E o B‟ en dirección vertical hacia arriba. Por ejemplo, el gradiente hidráulico hacia arriba en E, puede estimarse como Δh dividida dividida por la distancia DE. Si este gradiente excede del valor crítico, el suelo que está inmediatamente aguas abajo del tablestacado se convertirá en movedizo y puede ocurrir una falla. El estudio de las redes de flujo (fig. 2.9) aclarará varias condiciones bajo las cuales pueden producirse condiciones de arena movediza en la práctica, a menos que se tomen precauciones especiales, como la de aumentar la longitud de las tablestacas o añadir bermas, filtros o drenes.
23
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
Figura. Redes de flujo para varias condiciones. a) Ataguía para la construcción de una pila de puente. b) Excavación para cimentación abajo del nivel del agua freática en arena. c) Desagüe de excavación bombeando el agua valiéndose de coladeras de punta.
24
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
En muchos problemas prácticos, el movimiento del agua no está confinado por una frontera artificial superior, sino que tiene lugar abajo de una superficie libre del agua (figs.b y c). La Línea de flujo más alta es entonces una línea de filtración (línea de corriente superior). Como la línea de filtración es una línea l ínea de flujo, las líneas equipotenciales que la intercepten deben hacerlo en ángulo recto, pero cada una de estas líneas termina en la línea de filtración. Para satisfacer los requisitos hidráulicos especiales de la línea de filtración, la componente vertical de la distancia entre las terminaciones de dos líneas equipotenciales adyacentes debe ser igual a la caída equipotencial Δh, como se indica en las figs.b y c. Para dibujar la red de flujo debe suponerse la posición de la superficie libre del agua, construir una red de flujo tentativa, revisar todos los criterios anteriormente discutidos, así como las condiciones especiales que existan a lo largo de la superficie libre, y revisar el diagrama, hasta que se satisfagan todas las condiciones. La construcción de estas redes de flujo es más difícil que si están fijas las condiciones de frontera superiores, pero los principios son los mismos. Los procedimientos que se acaban de describir pueden modificarse para tomar en cuenta la estratificación o valores diferentes de la permeabilidad en direcciones horizontales y verticales (A. Casagrande, 1935), Por supuesto, estas condiciones son las que se encuentran con más frecuencia en la práctica, y no las sencillas que se consideraron en los párrafos anteriores. En realidad, con frecuencia el patrón de permeabilidad es tan variable, que ninguna red de flujo puede representar satisfactoriamente las condiciones reales. Sin embargo, pueden lograrse conclusiones prácticas extremadamente útiles del estudio de las trayectorias de flujo en condiciones simplificadas.
5. HIPÓTESIS - El caudal de infiltración es directamente proporcional al tipo de suelo,sobre el cual esta construido la estructura hidráulica del mismo modo es un caudal constante respecto al tiempo.
6. VARIABLES Variable ind ependiente
Tiempo.- El tiempo es una magnitud física con la que medimos la duración o separación de acontecimientos, de los sistemas sujetos a observación. Para nuestro el tiempo implica, tiempo que demora en salir el agua del suelo (arena gruesa). Agua: Líquido que se hará circular a través del suelo (arena gruesa).
25
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
Suelo: Se denomina suelo a la parte superficial de la corteza terrestre, biológicamente activa, que proviene de la desintegración o alteración física y química de las rocas y de los residuos de las actividades de seres de seres vivos que se asientan sobre ella. Material que se utilizara para el trasporte del líquido. Variable dependiente
Coeficiente de permeabilidad: Es una característica de los suelos, específicamente está ligado a la Ley la Ley de Darcy que se refiere al flujo de fluidos a través de los suelos. El coeficiente de permeabilidad, generalmente representado por la letra k, es extremadamente variable, según el tipo de suelo (arena gruesa).
7. DISEÑO DE CONTRASTACIÓN CONTRASTACIÓN DE EXPERIENCIAS. Para este trabajo se utilizó el siguiente modelo de contrastación: DISEÑO EXPERIMENTAL (GRUPO EN PALALELO) El proyecto consistirá en la elaboración de una maqueta a pequeña escala en la que se pueda apreciar como el fluido, en nuestro caso agua, fluye por lo vacíos continuos del suelo a utilizar, se contara con arena gruesa bien graduada para así poder obtener valores ya establecidos y que se demostrara mediante la práctica y los resultados obtenidos. También se obtendrá el valor del coeficiente de permeabilidad “K “de la arena bien graduada y este valor se obtendrá mediante la demostración que también se realizara. EN LA GRAFICA SE PRESENTA UN BOSQUEJO DE LO QUE SERIA LA MAQUETA.
POBLACION Y MUESTRA:
El suelo con el que se trabajara será la de una Arena bien graduada. No se trabajara con la gran variedad de suelos como la A rcilla o el Limo ya que no se podría visualizar muy bien como fluye el agua a través de él y se incurriría en error ya que el tiempo de salida del agua a través de estos suelos es mayor. TECNICAS DE ANALISIS: De la información que se encontró se optaron por las que mejor se acomodarían a lo que se busca calcular.
26
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
-
UNASAM – FIC
Para calcular el coeficiente coeficiente de permeabilidad de la arena bien graduada graduada se optó por por realizar la prueba de carga constante, la cual es sencilla de realizar.
-
PERMEAMETRO A CARGA CONSTANTE
-
Para la con con fiabilidad de resultados se hará uso de tablas en las que se establezcan rangos para el coeficiente de permeabilidad.
VALORES DEL COEFICIENTE DE PERMEABILIDAD EN DISTINTOS SUELOS (TABLA I)
-
Para la determinación de las líneas de flujo y líneas equipotenciales se realizara el trazado según los pasos ya establecidos.
27
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
IMÁGENES DE REDES DE FLUJO A MODO DE ILUSTRACIÓN:
-
Para la comprobación del dibujo realizado a mano se podría utilizar un software como el Geo - Slope.
IMÁGENES DEL SOFTWARE COMO EL GEO - SLOPE.
28
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
29
UNASAM – FIC
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
COEFICIENTE DE PERMEABILIDAD En el presente informe se muestra el resultado de la prueba de permeabilidad a carga constante para una muestra de suelo con arena lavada sin presencia de materia orgánico que había sido consumida en agua con una anterioridad de 24 horas, donde se introdujo en el molde con un diámetro de 10 cm, con 5 capas del suelo, se aplicó energía a la arena en el molde en forma de 25 golpes por capa. La idea de este ensayo es, antes de realizar cada medición, tener saturado el suelo, para que a la hora de incorporarle agua por medio de la filtración, haya un volumen de agua removido. Donde se realizaron medidas de volumen desalojado
a distintos H y distintos tiempos, para así formar un gradiente, que se relaciona que a mayor gradiente, gradiente, mayor velocidad del flujo de de agua que pasa por el suelo, suelo, se realiza la curva velocidad (Q/At) versus gradiente hidráulico (h/L), que cumple con la condición de flujo laminar, con la que se determina un coeficiente de permeabilidad. Se utilizó el método de carga constante para hallar el coeficiente de permeabilidad de la arena fina cuya arena fue retenida entre los tamices pasante de la malla Nº30 y retenida en la malla Nº100.
PERMEABILIDAD A CARGA CONSTANTE: El ensayo de permeabilidad a carga constante es de gran importancia, ya que el suelo, en su mayoría, está constantemente en contacto con el agua, y es de gran interés para fines ingenieriles conocer el comportamiento del suelo. Como lo sería en el caso en la construcción de taludes, que uno de los factores de gran importancia para que no falle el talud seria conocer la permeabilidad del suelo en su diseño. Los suelos tienen vacíos interconectados a través de los cuales él puede fluir de puntos de alta energía a puntos de baja energía. El estudio del flujo del agua a través de un suelo como medio poroso es importante en la mecánica de suelo, siendo necesario para estimar la cantidad de infiltración subterránea bajo varias condiciones hidráulicos, para investigar problemas que implican el bombeo de
30
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
agua para construcciones subterráneas y para el análisis de estabilidad de las presas de tierra y de estructuras de retención de tierra sometidas a fuerzas de filtración. Este método de prueba cubre la determinación del coeficiente de la permeabilidad con un método constante de cabeza para el flujo laminar de agua a través de suelos granulares. El procedimiento es para establecer los valores representativos del coeficiente de permeabilidad de los suelos granulares que pueden ocurrir en depósitos naturales como colocado en muros de contención.
31
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
32
UNASAM – FIC
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
33
UNASAM – FIC
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
34
UNASAM – FIC
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
PROCEDIMIENTO EN LABORATORIO:
35
UNASAM – FIC
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
36
UNASAM – FIC
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
37
UNASAM – FIC
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
38
UNASAM – FIC
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
39
UNASAM – FIC
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
40
UNASAM – FIC
UNASAM – FIC
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
MEMORIA DE CALCULOS:
Datos obtenidos en laboratorio.
Tiempo (Segundos) 0 30 60 90 120 150 180 210 240 270 300 330 360 390 420 450 480 510 540 570 600
Altura de ascenso (cm) 0 3 13,5 21 26,5 30 31,8 33 33,5 34 34 34,5 34,9 35 35 35 35,2 35 34,8 34,6 34,6
Longitud del permeámetro (L)
10 cm 0.45 cm
Diferencias de Altitudes Altitudes ∆H
volumen del liquido
192 ml
Área del permeámetro (A)
78.5398163 cm2
Demostracion en laboratorio :
Carga hidraulica final
=
,45= 0,045
41
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
Datos corregidos: Tiempo (Segundos) 0 30 60 90 120 150 180 210 240 270 300 330 360 390
42
Altura de ascenso (cm) 0 3 13.5 21 26.5 30 31.8 33 33.5 34 34 34.5 34.6 34.6
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
43
UNASAM – FIC
UNASAM – FIC
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
Tiempo (Segundos)
H/L
Altura de ascenso (cm)
V(cm3)
Q
v=Q/A
0
0
35.05
0.00
0
0
30
3
32.05
16.65
0.554913295
0.007065376
60
13.5
21.55
74.91
1.248554913
0.015897095
90
21
14.05
116.53
1.294797688
0.016485876
120
26.5
8.55
147.05
1.225433526
0.015602704
150
30
5.05
166.47
1.10982659
0.014130751
180
31.8
3.25
176.46
0.980346821
0.012482163
210
33
2.05
183.12
0.872006606
0.011102733
240
33.5
1.55
185.90
0.774566474
0.009862087
270
34
1.05
188.67
0.698779705
0.00889714
300
34
1.05
188.67
0.628901734
0.008007426
330
34.5
0.55
191.45
0.580136626
0.007386529
360
34.6
0.45
192.00
0.533333333
0.006790611
390
34.6
0.45
192.00
0.492307692
0.006268256
44
GRAFICA CURVA VELOCIDAD (Q/At) vs GRADIENTE HIDRAULICO (h/L)
CURVA VELOCIDAD (Q/At) vs GRADIENTE HIDRAULICO (h/L) 0.018 0.016 0.014 y = -0.0002x + 0.0117 R² = 0.0438
0.012 0.01 0.008 0.006 0.004 0.002 0 35.05 32.05 21.55 14.05
8.55
5.05
3.25
2.05
1.55
1.05
1.05
0.55
0.45
0.45
De la gráfica se calcula el coeficiente de permeabilidad del suelo la cual resulta ser la pendiente positiva de la recta lineal.
GRAFICA CURVA VELOCIDAD (Q/At) vs GRADIENTE HIDRAULICO (h/L)
CURVA VELOCIDAD (Q/At) vs GRADIENTE HIDRAULICO (h/L) 0.018 0.016 0.014 y = -0.0002x + 0.0117 R² = 0.0438
0.012 0.01 0.008 0.006 0.004 0.002 0 35.05 32.05 21.55 14.05
8.55
5.05
3.25
2.05
1.55
1.05
1.05
0.55
0.45
0.45
De la gráfica se calcula el coeficiente de permeabilidad del suelo la cual resulta ser la pendiente positiva de la recta lineal.
Como se tuvo un promedio de temperatura ambiente de 18°C se hará la corrección con el coeficiente de viscosidad del agua de 0.9077 con relación a una temperatura ideal de 24°C.
En este laboratorio de permeabilidad de carga constante se utilizó una pequeña muestra de arena, la cual se encontraba saturada y cumplía con los requisitos ya mencionados.
‘‘PROYECTO DE INVESTIGACION – – MECANICA MECANICA DE SUELOS I’’
UNASAM – FIC
8. CONCLUCIONES GENERALES Y RECOMENDACIONES
El coeficiente de permeabilidad depende depende de la porosidad del del suelo, vacíos continuos existente en el suelo.
El coeficiente de de permeabilidad nos permite calcular el caudal de infiltración que fluye por un determinado suelo.
Al trabajar con arenas gruesas y bien graduadas, se buscó trabajar con un flujo laminar ya que la ley de Darcy está dada para flujos laminares y suelos isotrópicos.
Las redes de flujo están formadas por una cuadricula o combinación de familia de líneas ortogonales es decir las líneas de flujo y las equipotenciales.
La línea de flujo es una línea a lo largo largo de la cual una una partícula de agua viaja del lado de aguas arriba al lado de aguas abajo a través de un suelo permeable.
Las franjas formadas entre dos líneas de flujo se les denomina denomina canales de flujo, las cuales se deben de aproximar a un cuadrado, para facilitar el cálculo del caudal de infiltración, aplicando diferencias finitas.
9. REFERENCIAS BIBLIOGRÁFICAS
VÉLEZ Otálvaro, María Victoria. Hidráulica de aguas subterráneas. Facultad de minas, Universidad Nacional de Colombia. Segunda edición, Medellín 1999.
JUÁREZ Badillo, Eulalio. RICO Rodríguez, Alfonso. Mecánica de suelos, tomo III flujo de agua en suelos. Editorial Limusa. México 1980
FUNDAMENTOS DE INGENIERIA GEOTECNICA. Braja M.Das. Editorial THOMSON, 1era Edición, México D.F. Páginas del 99 al 112
MECANICA DE SUELOS, Juárez Badillo, Tomo T omo I „‟Fundamentos de la Mecánica de Suelos‟‟, Editorial LIMUSA, 2005. Páginas del 191 al 227
http://www.biblioteca.udep.edu.pe/BibVirUDEP/tesis/pdf/1_ 136_147_89_1258.pdf
http://fluidos.eia.edu.co/hidraulica/articuloses/interesantes/fl ujodeaguaenelsuelo/flujodelaguaenelsuelo.html.
46