INTRODUCCIÓN

En la entrega anterior (ver ref. 04), dimos inicio al abordaje matemático-conceptual del fenómeno de consolidación unidimensional de los suelos, donde destacan como principales ideas, el hecho de considerarlos en una condición saturada y de concebir los asentamientos verticales que ocurren diferidos al tiempo de aplicación de las cargas, como una consecuencia de disipación de la presión en exceso a la hidrostática “u”; lo cual se materializa en agua que es expulsada del suelo en cuestión, con una dirección de flujo en sentido vertical, hecho que provoca en la misma proporción que ocurre, el incremento de los esfuerzos efectivos “σ’” y por consiguiente los asentamientos anteriormente descritos. El fenómeno de consolidación unidimensional de los suelos, está por lo tanto íntimamente vinculado a la variación del exceso de la presión a la hidrostática “u”; una variación que está influenciada por el tiempo que ha pasado posterior a la aplicación de las cargas, y por la profundidad “z” del estrato que se esté analizando. Imagen N°01: una visión global de los contenidos a abordar ![picture 02 corregida.png](https://cdn.steemitimages.com/DQmepWq3cCoBpyskzKKTUB81btpNBFC5ZuZj4wTYStut4Ba/picture%2002%20corregida.png)

Fuente: @eliaschess333, año: 2018. Nota: esquema conceptual elaborado por el autor, con ayuda de las herramientas Microsoft PowerPoint y Paint.

Lo anteriormente expuesto, constituye las principales hipótesis de partida, que permitieron a Terzaghi el desarrollo de la “Ecuación Diferencial del Proceso de Consolidación Unidimensional con Flujo Vertical”: ![picture 01.png](https://cdn.steemitimages.com/DQma3qLoJe1FBVRSoXLDoZ6YQ9A2tHfjKcdhpagRcpxP4Q7/picture%2001.png) Esta ecuación es presentada de una forma alternativa, tal como se aprecia en el esquema conceptual de la imagen N°01, y es el punto de partida para deleitarnos con los aportes de las matemáticas enmarcados en la aplicación de técnicas de resolución de ecuaciones diferenciales homogéneas con coeficientes constantes y así posteriormente aplicar las series de Fourier para determinar los valores de las constantes de integración involucradas; todos estos procedimientos los estaremos desarrollando con la finalidad de encontrar la solución analítica de la ecuación N°01, y así seguir profundizando en este interesante fenómeno de importancia en la ingeniería geotécnica, como lo es la consolidación unidimensional de los suelos.

UNA VISIÓN DIDÁCTICA DEL ESTABLECIMIENTO DE LAS CONDICIONES DE BORDE

Dado que deseamos resolver de forma analítica la ecuación N°01, es necesario establecer unas condiciones de borde, lo cual orientará los procesos de cálculo diferencial que desarrollaremos en líneas posteriores. Para tal fin, tomemos como referencia el esquema conceptual de la imagen N°02: Imagen N°02: ilustración didáctica de la condición frontera “u=0” para “z=0 y z=2H” ![picture 03.png](https://cdn.steemitimages.com/DQmbmi5vmq83K9ZfCBBUcVgWds25YHympTF7G6jzzcZEUkB/picture%2003.png)

Fuente: @eliaschess333, año: 2018. Nota: esquema conceptual elaborado por el autor, con ayuda de las herramientas Microsoft PowerPoint y Paint.

Esta primera condición de borde, establece que tiempo después de aplicada la carga “∆σ”, el exceso de presión a la hidrostática “u” en la parte superior e inferior de un estrato que se supone drenado por ambas caras, es cero; hecho que se puede interpretar fácilmente, asociándolo con una salida de agua, la cual ocurre gradualmente por la baja permeabilidad de los suelos que experimentan este fenómeno, como por ejemplo, las arcillas. En el eje de simetría, resaltado por el color “naranja”, se supone que “u” alcanza sus valores máximos, dado que el mismo se considera impermeable, lo cual es un criterio aceptable, ya que el agua tiene dos posibles trayectorias de salida, debido a que el estrato está drenado por ambas caras. La otra condición frontera de importancia se ilustra en el esquema conceptual de la imagen N°03: Imagen N°03: ilustración didáctica de la condición frontera “u=∆σ” para “t=0 y 0Fuente: @eliaschess333, año: 2018. Nota: esquema conceptual elaborado por el autor, con ayuda de las herramientas Microsoft PowerPoint y Paint. Para comprender esta condición, tomemos como referencia la analogía mecánica de Terzaghi, sobre la que puedes profundizar en la referencia N°04 donde realicé un video explicativo al respecto. Se tiene entonces que al aplicar las cargas “∆σ”, en un instante de tiempo teóricamente igual a “cero”, el agua que satura el estrato de suelo se presuriza, es decir, surge una presión en exceso a la hidrostática “u”, que se distribuye uniformemente en el mismo. Esto tiende a ocurrir en suelos de naturaleza arcillosa producto de su baja permeabilidad, siendo muy difícil una expulsión inmediata de agua, hecho que constituye, el sustento de la hipótesis anteriormente planteada. Una vez comprendida la lógica de las condiciones de borde, adentrémonos en la solución analítica de la ecuación N°01.

APLICACIÓN DE LA TÉCNICA DE RESOLUCIÓN DE ECUACIONES DIFERENCIALES HOMOGÉNEAS CON COEFICIENTES CONSTANTES

Lo anteriormente expuesto, permite establecer que “u” depende de la profundidad “z” y del tiempo “t”, lo cual podemos representar de la siguiente forma: ![picture 05.png](https://cdn.steemitimages.com/DQmeaWQVEdbWfQJ6qg9CqT1Fpxe5shm4D5qEd1sdkTKf4Ud/picture%2005.png) Esta ecuación sirve de base para reescribir la ecuación N°01, procedimiento que se ilustra en el esquema conceptual de la imagen N°04: Imagen N°04: rescribiendo la ecuación diferencial de consolidación unidimensional ![picture 06.png](https://cdn.steemitimages.com/DQmeYQsVU2rkZbM5rDp3ybtFEHwCCsWGUDa6wxD5CvJxAkS/picture%2006.png)

Fuente: @eliaschess333, año: 2018. Nota: esquema conceptual elaborado por el autor, con ayuda de las herramientas Microsoft PowerPoint y Paint.

Basados en la forma alternativa de presentar la ecuación N°01, la misma puede ser rescrita de la siguiente manera: ![picture 07.png](https://cdn.steemitimages.com/DQmQaS6NqVEtrv9bDpfk3WBUBunxp4ffrevAvosgLF17JFM/picture%2007.png) El término resaltado en color “rojo” en la ecuación N°03, es una constante que se puede introducir, fundamentándonos en la igualdad que existe entre los cocientes resaltados por el color “naranja”. Esto contribuye al surgimiento de las siguientes ecuaciones: ![picture 08.png](https://cdn.steemitimages.com/DQmWVwv2UVLxyfTbfmtBsTXc9eZR2rjiSS1tov4Eqkj6SVR/picture%2008.png) Las ecuaciones resaltadas en color “azul”, son ecuaciones diferenciales homogéneas con coeficientes constantes, por lo que procedemos a aplicar las técnicas de resolución de ecuaciones diferenciales ajustadas a dicho caso. • ABORDAJE MATEMÁTICO ECUACIÓN N°04 Para tal fin nos valemos de las bondades de la función exponencial “e^x”, representadas por el hecho de que sus derivadas son múltiplos de la misma función. Por lo que para la ecuación N°04 se puede plantear como posible solución: ![picture 09.png](https://cdn.steemitimages.com/DQmQzBUSTpHug1xHAjk4NpwJExqTCAjmP3tc24gSj12TGpE/picture%2009.png) Obtenemos la primera y segunda derivada con respecto a “z” de la ecuación N°06: ![picture 10.png](https://cdn.steemitimages.com/DQmbLmpciKw4zFcgdbL4xSEXo9PprL3qBj399Aq3AKiYNk8/picture%2010.png) Al sustituir las ecuaciones N°06 y N°08 en la N°04 obtenemos lo siguiente: ![picture 11.png](https://cdn.steemitimages.com/DQmWYaGCq1hqQt5nWsJampHL12krT9oF8Dae48kGmw4jVi3/picture%2011.png) Sacamos factor común: ![picture 12.png](https://cdn.steemitimages.com/DQmV4AEbkyhRnYakeLdPdYhmkDpjfuFC3acyH4AEY1HfNf1/picture%2012.png) Dado que el término resaltado en color “verde” de la ecuación N°10 no se anula, para que esta se cumpla, la función auxiliar resaltada en color “rojo” debe ser cero, por lo que: ![picture 13.png](https://cdn.steemitimages.com/DQmZvs2AbK4PWPr5PJbfDAMXFvSA1CPCfQaUx9g9BkYD2vb/picture%2013.png) Aplicamos la técnica de resolvente, para encontrar los posibles valores de “m”, donde: a=1 b=0 c=-λ Entonces: ![picture 15.png](https://cdn.steemitimages.com/DQmTT57Ua3oy7ZcoQiQMSyzxn5VEDT2eGo1U8BEZWHnKA6r/picture%2015.png) Sustituyendo en la ecuación N°12 tenemos: ![picture 16.png](https://cdn.steemitimages.com/DQmVgGvfmBsRZfGBEUsEdB11Zs4vgvscXc4CA6kyKK9Lk4X/picture%2016.png) Por lo que las posibles soluciones de la ecuación N°04 vienen dadas por: ![picture 17.png](https://cdn.steemitimages.com/DQmPLgYdUvGHpNQKwYBod6g47eGfha7DvWfWN5YA47CjqTT/picture%2017.png) Siendo la solución general la suma de ambas, y en cuyo planteamiento se introducen las constantes arbitrarias de integración: ![picture 18.png](https://cdn.steemitimages.com/DQmWCfqouB27bUEUyw8PUtaSaLSpSTduC8nhd78FcNLqSxg/picture%2018.png) Basados en las condiciones fronteras presentadas en el esquema conceptual de la imagen N°02, y aplicadas a la ecuación N°02, se tiene: ![picture 19.png](https://cdn.steemitimages.com/DQmQnMhWHhMGkbrX4op8cTKEJoKeb71h6b3FbzQpABTwmaU/picture%2019.png) Estos resultados permiten abordar la ecuación N°14 con la finalidad de esclarecer el camino en la búsqueda de las constantes de integración “c1” y “c2”, por lo que al evaluar dicha ecuación para “z=0” se tiene: ![picture 20.png](https://cdn.steemitimages.com/DQmaV1U1thJCLsxxr8jNimnmiDLtSnqjvS3za2uJLFsw7W2/picture%2020.png) Esto permite rescribir la ecuación en cuestión de la siguiente manera: ![picture 21.png](https://cdn.steemitimages.com/DQmcoQcoTFnPxq3e7nGYEXgSirhNAHjqghpiw5ZzV4xdZGE/picture%2021.png) La ecuación N°18 nos permite introducir la identidad trigonométrica asociada al seno hiperbólico: ![picture 22.png](https://cdn.steemitimages.com/DQmRrZqASq9xVdoq7GAcFgnqGC3mrknieX5zg4STrXJBvkd/picture%2022.png) Al sustituir la ecuación N°19 en la N°18, obtenemos: ![picture 23.png](https://cdn.steemitimages.com/DQmUcR4FHwVNpKa5CNWmcrM6tw6hTAMjhojTTHgFSpWV5Vt/picture%2023.png) Basado en la condición de borde desarrollada en la ecuación N°16, sustituimos en la ecuación N°20 (ver términos resaltados en color “azul”): ![picture 24.png](https://cdn.steemitimages.com/DQmYgsm8DHNsAmXXvdMgFY87GkEMqbsamWf3a7PpoHLveYJ/picture%2024.png) Estratégicamente procedemos a rescribir la ecuación N°21 de la siguiente manera: ![picture 25.png](https://cdn.steemitimages.com/DQmVYxSr6awn6ZdjgZHxgewL4yjSezi9ThJz5bdNVmNH888/picture%2025.png) Dado que el seno hiperbólico de un argumento imaginario es el seno trigonométrico con el mismo argumento, se puede plantear como valor de “λ” el siguiente: ![picture 26.png](https://cdn.steemitimages.com/DQmbQpiv7KVYCFocg1LkeFYeS9mxN9opxBHCLkn3jT21ynE/picture%2026.png) Sustituimos este valor en la solución en la ecuación N°22: ![picture 27.png](https://cdn.steemitimages.com/DQmWS9ffxF1tqDDSFFGMRyG6vJzkFpnaE9MYv81iMSovWup/picture%2027.png) A partir de estas deducciones podemos reescribir la ecuación N°20, teniendo la precaución de colocar en el argumento de la función seno, una relación de proporcionalidad entre la profundidad de interés “z” y el espesor del estrato “2H”; esta relación se encuentra resaltada por el color “rojo”: ![picture 28.png](https://cdn.steemitimages.com/DQmegMfSvPFPeiCQbgKMcuCsniSsN1Yn7nb9h5WLbdLxyvq/picture%2028.png) Esta ecuación constituye una forma general de presentar la ecuación N°21. A pesar de que el término resaltado en color “naranja” aparece multiplicado por “2”, puede considerarse como la existencia de una única constante arbitraria, por lo que para hacer más cómodo el manejo de esta ecuación, al momento de sustituir en la ecuación N°02, se tiene que: ![picture 29.png](https://cdn.steemitimages.com/DQmfZLnjA1YUtYzjXMN32u2Heh5Mu6zUbhMoYPPHD7gaE6M/picture%2029.png) • ABORDAJE MATEMÁTICO ECUACIÓN N°05 De igual forma como se trató la ecuación N°04, planteamos como posible solución la siguiente: ![picture 30.png](https://cdn.steemitimages.com/DQmTHaGh2GCAFE4Lq6Y9zgJPBo5Mm8Nq3CYjsQ8rEDqqrk4/picture%2030.png) Obtenemos la primera derivada con respecto al tiempo “t”: ![picture 31.png](https://cdn.steemitimages.com/DQmTQGLP7zzJEZULigFECrUxxMijoMdFkAsERQFWUrprXMG/picture%2031.png) Sustituimos las ecuaciones N°27 y N°28 en la ecuación N°05: ![picture 32.png](https://cdn.steemitimages.com/DQmaUScJ6KZbXASM14xP7dA9hQRkYprqM7NYFpUu6jxHQwJ/picture%2032.png) Dado que el término exponencial no se cancela, la ecuación N°29 se satisface para cuando la función auxiliar resaltada por el color “azul” es igual a cero, por lo tanto: ![picture 33.png](https://cdn.steemitimages.com/DQmU8qBTMH7Z1NNPYh4ahSFR5YZwYXruD9hbCBFRCNu4TAj/picture%2033.png) Por lo que: ![picture 34.png](https://cdn.steemitimages.com/DQmYEQGSmq17Cjrixbvnea8puusDZFcQvkDHoGzVRoqPxTb/picture%2034.png) Introduciendo la constante arbitraria de integración y tomando como referencia el valor para “λ” planteado en el tópico anterior, la solución de la ecuación N°05 viene dada por: ![picture 35.png](https://cdn.steemitimages.com/DQmbpak8vWsKq5DP3G1Dqj6kQRrRzhvxJSambpw6PMAGUqg/picture%2035.png) • UNA VISIÓN DE LOS AVANCES LOGRADOS EN LA BÚSQUEDA DE LA SOLUCIÓN ANALÍTICA DE LA ECUACIÓN N°01 Dado que el proceso de consolidación está gobernado por las variaciones de “u” en el tiempo “t” y profundidad “z”, se procedió al desarrollo de la ecuación N°02, y de las que se presentan en el esquema conceptual de la imagen N°04. El abordaje matemático de tales ecuaciones, ha permitido obtener lo siguiente: ![picture 36.png](https://cdn.steemitimages.com/DQmUA8BPEqNJv9W4QH5KaG4E7oUrfsVB8Zj5rC9ReUmnju6/picture%2036.png) Rescribimos: ![picture 37.png](https://cdn.steemitimages.com/DQmcJPno7bQsuW1FERLprFoKNeRAQZUkm6ApZsR4xLeKMEt/picture%2037.png) El término “An” representa las constantes arbitrarias de integración, que depende de los valores de “n”. Como se puede apreciar, el reto que a continuación sigue es encontrar el valor de “An”. Para tal fin, haremos uso de otro de los grandes aportes de las matemáticas como lo son las series de Fourier, lo cual se desarrolla a continuación:

SERIES DE FOURIER ¡AL RESCATE!

Las Matemáticas nos sorprenden con la infinidad de aplicaciones que contribuyen a una mejor comprensión de los fenómenos físicos que ocurren a nuestro alrededor, centrándonos en este caso en el fenómeno de consolidación unidimensional de los suelos. Haremos uso de las series de Fourier para finalmente dar con la solución analítica de la ecuación N°01. Para una mejor comprensión de este tópico, lo he subdivido en las siguientes partes: • UN ENFOQUE TEÓRICO-CONCEPTUAL DE LAS SERIES DE FOURIER La presencia de derivadas parciales en la ecuación N°01 producto de que la presión en exceso a la hidrostática “u”, depende de dos variables como lo son el tiempo “t” y la profundidad “z”, y la necesidad de encontrar el valor de “An” para dar finalmente con la solución general de dicha ecuación, son situaciones que abren las puertas al uso de las series de Fourier. En primer lugar presentemos la forma general de una serie de Fourier: ![picture 38.png](https://cdn.steemitimages.com/DQmXrxxwyR2iEP72juWnbhfKMuTD4UHAwbBxpq974TxsQEM/picture%2038.png) Se establecen como principales hipótesis en torno a la ecuación N°35, que la suma de sus términos converge al valor de la función, y que toda la función es integrable en el rango de valores que va de “-L” a “L”, siendo el cero “0” el punto medio del intervalo y matemáticamente se expresa así: ![picture 39.png](https://cdn.steemitimages.com/DQmWvoBwgv8VnUX232qySYQDrkhmVBgcwQ7wfnGBUiWQe6h/picture%2039.png) El reto es encontrar los valores de los coeficientes resaltados en color “azul” en la ecuación N°35. Para encontrar por ejemplo el coeficiente “ao” se procede a la integración de cada uno de los términos de dicha ecuación, lo cual se ilustra didácticamente en el siguiente esquema conceptual: Imagen N°05: determinación analítica del coeficiente “ao” ![picture 40.png](https://cdn.steemitimages.com/DQmUsNKkSnNy7KLAbbpWCdZ6vEdY4koqdwDCBcQC9e1ebtn/picture%2040.png)

Fuente: @eliaschess333, año: 2018. Nota: esquema conceptual elaborado por el autor, con ayuda de las herramientas Microsoft PowerPoint y Paint.

Procedemos a resolver cada una de las integrales señaladas en la imagen N°05. Vale destacar que en el proceso de obtención de la integral definida, el término resaltado en color “verde” representa los límites de evaluación de la integral. En lo que respecta a la “integral 1” se tiene: ![picture 41.png](https://cdn.steemitimages.com/DQmXG68opT2qcxJVPrB1req7PtQjnL84VyfVMVtSTeZGgRA/picture%2041.png) Para el caso de la “integral 2”: ![picture 42.png](https://cdn.steemitimages.com/DQmek2R61hHnKEfCpnjAqgySr8ZYMyaDA6qiC7DiemoGHvD/picture%2042.png) Aplicamos la técnica de integral por sustitución, donde: ![picture 43.png](https://cdn.steemitimages.com/DQmSBBwyH2cPbxLDP2fudpveU8tEgToQmARH7VyWsaBhzk3/picture%2043.png) Sustituimos las ecuaciones N°38 y N°39 en la N°37: ![picture 44.png](https://cdn.steemitimages.com/DQmSMzeN5oRNEt8gY2F1op9a8FSLXGr88owsppmdV8Z6K49/picture%2044.png) Devolvemos el cambio, sustituyendo la ecuación N°38 en la N°40: ![picture 45.png](https://cdn.steemitimages.com/DQmNsFnN14Kbn1dqji93hrVCvjJ6ACATbEZmWQWqUNhdzjZ/picture%2045.png) Al ser el seno de cualquier múltiplo de “π” igual a cero, se tiene que: ![picture 46.png](https://cdn.steemitimages.com/DQme8Je4GviYVQzL8kyCqNMe7c4doP5meFefoU1R4GBTgQr/picture%2046.png) Algo similar ocurre en la solución de la “integral 3”: ![picture 47.png](https://cdn.steemitimages.com/DQmTLELNAZa7Mvv5Dh4c4pdag32h8cJHQ6eADnGYmZg4YBB/picture%2047.png) Dado que la función a integrar tiene el mismo argumento que el de la “integral 2”, podemos usar las ecuaciones N°38 y N°39 en el proceso de integración por sustitución: ![picture 48.png](https://cdn.steemitimages.com/DQmTZfK8P6NdFD5HPqfuPRZP75mFKq3g6RFgeKKCiyjGCHo/picture%2048.png) Devolvemos el cambio sustituyendo la ecuación N°38 en la N°44: ![picture 49.png](https://cdn.steemitimages.com/DQmTYcPGH5rK9ZHu8r5q1swC2RZf2MiEz11tT9BFabbaxCC/picture%2049.png) Por ser el coseno una función par: ![picture 50.png](https://cdn.steemitimages.com/DQmehT9nPPx9J2fb6nts3Ndh2DnjDgos1xyj3bUmU2tSVKT/picture%2050.png) El término resaltado en color “rojo” de la ecuación N°45 se cancela, obteniendo que: ![picture 51.png](https://cdn.steemitimages.com/DQmWajffxfUsbFrR87ctCwCUq9wtgAMLzBTy8TXVHV1k8Rx/picture%2051.png) Dado que el resultado de “integrales 2 y 3” es cero, estas se omiten al realizar la suma de los resultados del proceso de integración definida en pro de conseguir el coeficiente “ao”, en este orden de ideas: ![picture 52.png](https://cdn.steemitimages.com/DQmcXEVjHzYV1QtZtvM4hE5khxTyzRakMF1c6vGEfncEDzF/picture%2052.png) Contextualicemos todos estos conocimientos al caso en estudio. • APLICACIÓN DE LAS SERIES DE FOURIER AL CASO EN ESTUDIO Basados en la condición de borde del esquema conceptual de la imagen N°03: ![picture 53.png](https://cdn.steemitimages.com/DQmYm7QRuqWWFRUnSi7wHaoHHdKxLANzzzSuZ5GCrF4xtgh/picture%2053.png) Evaluamos la ecuación N°34 para este tiempo: ![picture 54.png](https://cdn.steemitimages.com/DQmatu7z2xS6sTUWg6aSMh3pzAt2XHU4G7DGCApYm4vYSso/picture%2054.png) Por lo tanto: ![picture 55.png](https://cdn.steemitimages.com/DQmRAeyTgdDcWPRtVqPhdpmWZuptMGjnE5BYMtceWUQhX5A/picture%2055.png) Dado que “∆σ” es constante en todo el proceso de consolidación, “u” también lo es para el instante de tiempo teóricamente considerado como cero “0”, por lo que es necesario concebir “u” como una suma infinita de términos sinusoidales, dando así la bienvenida a la serie de Fourier: ![picture 56.png](https://cdn.steemitimages.com/DQmdogmYT9Bk8a1nEz5SGBPEeytcuLDJF3D1ECaU5g1VQHj/picture%2056.png) Al contrastar la ecuación N°50 con la N°35 (forma general de la serie de Fourier), el problema en estudio sólo tiene el término relativo a la función seno, por lo se reduce a encontrar el coeficiente “bn” (Ec.35), que es el mismo “An” (Ec.50). En el tópico anterior demostramos que al aplicar la técnica de integral definida a cada miembro de la ecuación no nos ayuda a encontrar este coeficiente. Para tal fin, lo que se recomienda es multiplicar ambos lados de la ecuación N°50 por el término: ![picture 57.png](https://cdn.steemitimages.com/DQmcU16hPNnRWDZEqS8poCZeTf7fmQdLXbEUc3XvDiVCCVH/picture%2057.png) Obteniendo lo siguiente: ![picture 58.png](https://cdn.steemitimages.com/DQmaGBRc9ExyjHAVuqEe2JpYZU3T6GSupifwjMfRGZ1QGFB/picture%2058.png) Y de este modo proceder a integrar a ambos lados, acotando que los límites de integración, están definidos por el espesor del estrato “ver imágenes N°02 y N°03” ![picture 59.png](https://cdn.steemitimages.com/DQmacWf4TNueEZ3L2iSgY3pQKXHiN8QA8HteB7scXM5eSiW/picture%2059.png) Procedemos a resolver las integrales resaltadas por los colores “azul” y “verde”: • RESOLVIENDO LA INTEGRAL RESALTADA EN COLOR “AZUL” En el esquema conceptual de la imagen N°06 se ilustra este procedimiento: Imagen N°06: solución integral resaltada con el color "azul" (Ec.52) ![picture 60.png](https://cdn.steemitimages.com/DQmXxtnemE3RkqSGiGmrjomPCmSAVc8FrbocT3v8UX2Xch3/picture%2060.png)

Fuente: @eliaschess333, año: 2018. Nota: esquema conceptual elaborado por el autor, con ayuda de las herramientas Microsoft PowerPoint y Paint.

• RESOLVIENDO LA INTEGRAL RESALTADA EN COLOR “VERDE” PARA “n≠m” En primera instancia es necesario rescribir la integral, lo cual se plasma en la siguiente imagen: Imagen N°07: rescribiendo la integral resaltada en color "verde" (Ec.52), caso “n≠m” ![picture 61.png](https://cdn.steemitimages.com/DQmVdnrkxh9TafhqkoLsDemvHKcYx5NmopfCRHaL2Gbc66v/picture%2061.png)

Fuente: @eliaschess333, año: 2018. Nota: esquema conceptual elaborado por el autor, con ayuda de las herramientas Microsoft PowerPoint y Paint.

Procedemos a resolver la “integral 04”, de la siguiente manera: Imagen N°08: resolución de la “integral 04” ![picture 62.png](https://cdn.steemitimages.com/DQmUc5BYwKF5Wpx8Ns7rv4CKQEEMP6kZCrF6dHidojZ4R9d/picture%2062.png)

Fuente: @eliaschess333, año: 2018. Nota: esquema conceptual elaborado por el autor, con ayuda de las herramientas Microsoft PowerPoint y Paint.

La evaluación de la “integral 04” en los límites resaltados por el color “verde” y que representan el espesor del estrato de suelo, se ilustra a continuación: Imagen N°09: evaluación de la integral 04 ![picture 63.png](https://cdn.steemitimages.com/DQmVj9kPdXSjFocwtDioKrwQxdsKDbzmPen3Mm3iNUzJnVA/picture%2063.png)

Fuente: @eliaschess333, año: 2018. Nota: esquema conceptual elaborado por el autor, con ayuda de las herramientas Microsoft PowerPoint y Paint.

En lo que respecta a la “integral 05” tenemos lo siguiente: Imagen N°10: resolución de la integral 05 ![picture 64.png](https://cdn.steemitimages.com/DQmPDJxKpJzSGWB1K4JrtD1mdM6i3As4miBTm1vrivx739J/picture%2064.png)

Fuente: @eliaschess333, año: 2018. Nota: esquema conceptual elaborado por el autor, con ayuda de las herramientas Microsoft PowerPoint y Paint.

Evaluamos el resultado en los respectivos límites de integración: Imagen N°11: evaluación de la “integral 05” ![picture 65.png](https://cdn.steemitimages.com/DQmRSakjTjSmZikx8ncHQdMjgzQc1Ui4GepvN7A66ox4z3i/picture%2065.png)

Fuente: @eliaschess333, año: 2018. Nota: esquema conceptual elaborado por el autor, con ayuda de las herramientas Microsoft PowerPoint y Paint.

Estos resultados parciales se suman, para dar con el resultado de la integral resaltada por el color “verde” en la ecuación N°52; proceso que se muestra en la siguiente imagen: Imagen N°12: agrupación de resultados ![picture 66.png](https://cdn.steemitimages.com/DQmUnuVeijxphkq73VxVMRXvGGZDUY7XwFwW1yvKqhr9dd1/picture%2066.png)

Fuente: @eliaschess333, año: 2018. Nota: esquema conceptual elaborado por el autor, con ayuda de las herramientas Microsoft PowerPoint y Paint.

Se tiene por lo tanto, que el caso donde “m≠n” no contribuye a encontrar el coeficiente “An” dado que el resultado del proceso de integración es cero. • RESOLVIENDO LA INTEGRAL RESALTADA EN COLOR “VERDE” PARA “n=m” En el esquema conceptual de la imagen N°13 se rescribe esta integral: Imagen N°13: rescribiendo la integral resaltada en color "verde" (Ec.52), caso “n=m” ![picture 67.png](https://cdn.steemitimages.com/DQmSs52dZgmewa3hCz8d6NAr7Q6YsZZ3KwdroPKkoo4jY7t/picture%2067.png)

Fuente: @eliaschess333, año: 2018. Nota: esquema conceptual elaborado por el autor, con ayuda de las herramientas Microsoft PowerPoint y Paint.

La solución de las “integrales 06 y 07” se presenta a continuación: Imagen N°14: resolución de las integrales 06 y 07 ![picture 68.png](https://cdn.steemitimages.com/DQmdhKyZeUH7iJAMpdChv8N9guLJDg5yE61QvRbG2rq3j1h/picture%2068.png)

Fuente: @eliaschess333, año: 2018. Nota: esquema conceptual elaborado por el autor, con ayuda de las herramientas Microsoft PowerPoint y Paint.

Finalmente los resultados de las “integrales 06 y 07” se sustituyen en la integral resaltada en color “verde” de la ecuación N°52, obteniendo lo siguiente: ![picture 69.png](https://cdn.steemitimages.com/DQmbEr9zwVZj6s4Nq56vUi6X1f7yakH9HxBDGmuW7fUmNCF/picture%2069.png) • ENCONTRANDO EL VALOR DE LAS CONSTANTES DE INTEGRACIÓN Sustituimos en la ecuación N°52 los resultados de integración, plasmados en el esquema conceptual de la imagen N°06 y de la ecuación N°53: ![picture 70.png](https://cdn.steemitimages.com/DQmbBWg63B9p5H81Fd1Vuh3p6zs8sDitpPPWe7CC6dNHo4J/picture%2070.png) Analizando el término resaltado en color “verde” de la ecuación N°54: ![picture 71.png](https://cdn.steemitimages.com/DQmZoAU3ndxDfSLraSxnDCbhpiUgMPaivsDhds3fx5p3kuC/picture%2071.png) Se concluye por lo tanto que “An” tiene valor cuando “m” es impar, y el resultado en estos casos es dos. Este hecho permite rescribir la ecuación N°54: ![picture 72.png](https://cdn.steemitimages.com/DQmcUhBi3sdCmP9p8eXaoQL6qyLekHjz6Me6wV2SkCu41Wc/picture%2072.png)

Y ASÍ LLEGAMOS A LA SOLUCIÓN ANALÍTICA DE LA ECUACIÓN DIFERENCIAL

Para tal fin, sustituimos el resultado de “An” en la ecuación N°50 ![picture 73.png](https://cdn.steemitimages.com/DQmQuQUKcLQPVhufuGn3ih4igsy1BDR8HzRFpuNtnyoB6F9/picture%2073.png) Rescribimos, tomando al mismo tiempo en consideración de que “n=m” ![picture 74.png](https://cdn.steemitimages.com/DQmdCQ1FdozRfrCCiLMVBnx2K14FZCtoyy9XSinaPPqQ3iz/picture%2074.png) De acuerdo a lo planteado en líneas anteriores, es necesario que “m” sea un número impar, por lo tanto podemos establecer la siguiente ecuación: ![picture 75.png](https://cdn.steemitimages.com/DQmVRAbuhGdpH32j8u2mxyZSnzKbHsp97udfxcKzsmQtpxz/picture%2075.png) En este sentido la ecuación N°57 queda así: ![picture 76.png](https://cdn.steemitimages.com/DQmP9Bq9Sdv5NVD58LXabv7LVy8gM2yMZmyCrcvQHtKxo6n/picture%2076.png) Se desprende así una nueva forma de presentar el valor de “An”, representada por el término resaltado por el color “rojo” de la ecuación N°59: ![picture 77.png](https://cdn.steemitimages.com/DQmbyarUALcGKSiBAFrd7HA5qNdsFzJzHGrVGRQTxpo4fzy/picture%2077.png) La ecuación N°59 está ajustada a la condición de borde N°02, donde para el instante de tiempo “t=0” el incremento de presiones “∆σ” producto de una carga externa aplicada, es igual al exceso de la presión a la hidrostática “u”. Para cualquier instante de tiempo, tenemos lo siguiente: ![picture 78.png](https://cdn.steemitimages.com/DQmRfgU7trUF6sZ1K5pnnXrvPxUP4PYdjmsUwcPHZoqU3hh/picture%2078.png) Partiendo del hecho de que “n=m” tenemos: ![picture 79.png](https://cdn.steemitimages.com/DQmRuTKJtrETVauh5mEg3foKDBDV4LgH2SrTYoPbPjpMgaa/picture%2079.png) Sustituimos el valor de “An” dado en la ecuación N°60 y la forma de presentar “m” según ecuación N°58, llegando así a la solución analítica de la ecuación diferencial de la consolidación unidimensional de los suelos para flujo vertical: ![picture 80 corregida.png](https://cdn.steemitimages.com/DQmfXY6MUdt46CnEedsqJz6QVqa9iiafCqwM7y9DNjFmiYf/picture%2080%20corregida.png)

CONCLUSIONES

Al analizar detalladamente la solución de la ecuación diferencial en cuestión, se deducen las siguientes ideas: 1.- La variación de la presión en exceso a la hidrostática “u”, está influenciada por el coeficiente de consolidación “Cv”, el cual se encuentra vinculado a propiedades del suelo, como son la permeabilidad “k” y el coeficiente de variación volumétrica “mv”, es decir: ![picture 81.png](https://cdn.steemitimages.com/DQmYaJyWMV4aB3TRPsnbfksxRd4CbSwxhArMVzZu73y1Nj8/picture%2081.png) 2.- La baja permeabilidad “k” de las arcillas, repercute directamente en el proceso de disipación de la presión en exceso a la hidrostática. 3.-El incremento de carga “∆σ” es en definitiva lo que propicia el surgimiento de la presión en exceso a la hidrostática “u”; si este valor es cero, esta última también tomará ese valor. Teniendo la solución general de la ecuación diferencial del fenómeno de consolidación de los suelos, el reto es deducir expresiones de interés que contribuyan al tratamiento de la data obtenida a partir de la realización de este ensayo en un laboratorio de Mecánica de Suelos. En una próxima entrega, profundizaré en este aspecto. Escribió para ustedes: @eliaschess333

FUENTES DE INFORMACIÓN CONSULTADAS

1.- BADILLO J. Y RODRÍGUEZ R. 2007. MECÁNICA DE SUELOS TOMO I. FUNDAMENTOS DE LA MECÁNICA DE SUELOS. EDITORIAL LIMUSA S.A. GRUPO NORIEGA EDITORES 2.- BADILLO J. Y RODRÍGUEZ R. 1995. MECÁNICA DE SUELOS TOMO III. FLUJO DE AGUA EN LOS SUELOS. EDITORIAL LIMUSA S.A. GRUPO NORIEGA EDITORES 3.- SIMMONS G.Y ROBERTSON J.1993. ECUACIONES DIFERENCIALES CON APLICACIONES Y NOTAS HISTÓRICAS. SEUNDA EDICIÓN. McGRAW-HILL

LECTURAS RECOMENDADAS

4.- SANTANA E. 2018. ABORDAJE DIDÁCTICO DEL FENÓMENO DE CONSOLIDACIÓN UNIDIMENSIONAL DE LOS SUELOS DESDE UN ENFOQUE MATEMÁTICO CON BASAMENTO EN EL PRINCIPIO DE ESFUERZOS EFECTIVOS Y LA LEY DE PERMEABILIDAD DE DARCY. PARTE I. DISPONIBLE EN: https://steemit.com/stem-espanol/@eliaschess333/ozrzfpbw