Cálculo de daño en pavimentos mediante un análisis hidro-mecánico de pavimentos flexibles

<!-- metadata commented in wiki content


<div id="_GoBack" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<big>Cálculo de daño en pavimentos mediante un análisis hidro-mecánico de pavimentos flexibles</big></div>

Carlos Chávez <sup>a</sup>, Manuel A. Hernández <sup>b</sup>, Francisco J. Domínguez <sup>a</sup> y J. Eleazar Arreygue<sup>a</sup>

<sup>a</sup> Facultad de Ingeniería Civil, Universidad Michoacana de San Nicolás de Hidalgo (UMSNH)

<sup>b</sup> Estudiante de la maestría en infraestructura del transporte en la rama de las vías terrestres (UMSNH).

[mailto:cachavez@umich.mx cachavez@umich.mx]; [mailto:man_her_30@hotmail.com man_her_30@hotmail.com]; [mailto:dmota@umich.mx dmota@umich.mx]; <span style="text-align: center; font-size: 75%;">[mailto:arrocha@umich.mx arrocha@umich.mx]</span>;

'''Resumen'''. Se presenta un análisis hidro-mecánico de un pavimento en el que se consideran la variación de las condiciones de contenido de agua para la determinación del módulo de elasticidad de las capas de suelo en el pavimento. Para ello se soluciona la ecuación diferencial de Richards y de equilibrio en condición axi-simétrica de manera simultánea, considerando una elasticidad no lineal que depende de la succión y del estado de esfuerzos en el suelo. Con este tipo de análisis es posible cuantificar la variación del daño con respecto al tiempo y visualizar el deterioro de las estructuras de pavimento. Se considera un caso práctico en el cual se cuantifica el incremento del daño en función de la variación de las condiciones ambientales y se obtienen conclusiones.
-->

'''Abstract'''. A hydro-mechanical analysis of a pavement is presented in which the variation of the water content is considered for the estimation of modulus of elasticity of the layers of soil in the pavement. To do this, the Richards differential equation and equilibrium in the axi-symmetric condition are solved simultaneously, considering a non-linear elasticity that depends on the suction and stress state in the soil. With this type of analysis, it is possible to quantify the variation of the damage with respect to time and visualize the deterioration of the pavement structures. It is considered a practical case in which the increase in damage is quantified in function of the variation of environmental conditions and conclusions are obtained.

'''Resumen'''. Se presenta un análisis hidro-mecánico de un pavimento en el que se consideran la variación de las condiciones de contenido de agua para la determinación del módulo de elasticidad de las capas de suelo en el pavimento. Para ello se soluciona la ecuación diferencial de Richards y de equilibrio en condición axi-simétrica de manera simultánea, considerando una elasticidad no lineal que depende de la succión y del estado de esfuerzos en el suelo. Con este tipo de análisis es posible cuantificar la variación del daño con respecto al tiempo y visualizar el deterioro de las estructuras de pavimento. Se considera un caso práctico en el cual se cuantifica el incremento del daño en función de la variación de las condiciones ambientales y se obtienen conclusiones.

Palabras clave: Análisis de Pavimentos, flujo de agua en suelos, daño en pavimentos

==1 Introducción==

La estructura de pavimentos está sometida a las solicitaciones del tránsito vehicular y clima; en particular la temperatura y el contenido de agua afectan los módulos con los cuales se realiza el análisis de esfuerzos y deformaciones de los pavimentos. Dentro de las metodologías mecanicistas de diseño de pavimentos, el estado de esfuerzos y deformaciones se relaciona con el número de repeticiones para la falla [1]. El comportamiento esfuerzo deformación de la carpeta asfáltica depende de la temperatura y falla por agrietamiento, el de las capas de suelo depende del contenido de agua y fallan por deformación permanente. Si cambian las condiciones anteriores la respuesta del pavimento varía. En el caso en que los módulos de resiliencia disminuyan el daño es mayor, esto sucede cuando la temperatura de la carpeta asfáltica aumenta y cuando existe infiltración de agua hacia las capas de suelo o un incremento en el nivel de aguas freáticas.

Los análisis de pavimentos se pueden realizar con la teoría de elasticidad multicapa de Burmister [2].Debido a las condiciones de contorno de la carga (neumático) se utiliza un análisis simplificado axi-simétrico alternativo al de tres dimensiones [1, 2]. En la <span id='cite-_Ref531712197'></span>[[#_Ref531712197|Fig. 1]] se presenta el sistema multicapa axi-simétrico simplificado que usualmente se utiliza para los análisis de pavimentos flexibles. En la <span id='cite-_Ref531712197'></span>[[#_Ref531712197|Fig. 1]], <math display="inline">a</math>, es el radio del área cargada que representa la llanta del vehículo, ''q'', el esfuerzo uniforme de la llanta (presión de inflado) y ''z'' las distancias a las capas de interés. Existe un número <math display="inline">n</math> de capas que tienen propiedades <math display="inline">E\, y\, \mu</math> . La teoría de Burmister considera que el módulo de elasticidad es constante, pero en la realidad depende del estado de esfuerzos y de las condiciones climáticas.

Alternativamente, se puede utilizar el método de los elementos finitos para resolver las ecuaciones diferenciales que resultan del problema elástico lineal en condiciones axi-simétricas que se plantea [3][4]. Como se menciona el clima afecta el comportamiento de los materiales con los que se construye un pavimento. En este estudio sólo se toma en cuenta la influencia de la variación de los contenidos de agua en las capas de suelo. Para ello en necesario solucionar la ecuación de Richards [5], para suelos parcialmente saturados.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Chavez_729651003-image1.png|474px]] </div>

<div id="_Ref531712197" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig. 1. Sistema Multicapa simplificado para el análisis de pavimentos flexibles [1].</div>

Con la solución de ambas ecuaciones diferenciales de manera simultánea, es posible estimar como varía el daño en pavimento cuando este está sometido a cambios en su contenido de agua, debido a la infiltración e incremento de nivel de agua freática.

==2 Análisis de pavimentos==

En la <span id='cite-_Ref531716444'></span>[[#_Ref531716444|Fig. 2]] se presenta el estado de esfuerzos en un punto en una condición axi-simétrica. El estado de esfuerzos es tridimensional, pero se puede estudiar con el uso de dos coordenadas (''r'', ''z'') si las cargas exteriores se consideran axi-simétricas y generan desplazamientos en direcciones radial (''u'') y axial (''w'') [6].

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">[[File:Fig 2.png|centre|685x685px]]</div>

<div id="_Ref531716444" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig. 2. Componentes de los esfuerzos de una sección axi-simétrico.</div>

===2.1 Teoría de la elasticidad en condiciones axi-simétricas===

De las ecuaciones de equilibrio en condenadas cilíndricas [7], se puede deducir las ecuaciones de equilibrio para problemas axi-simétricos. Debido a la simetría axial u y w son independientes de la coordenada <math display="inline">\theta</math>  [6], por tanto, las ecuaciones de equilibrio se simplifican a dos ecuaciones que son:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>\frac{\partial \left( r\, {\sigma }_{r}\right) }{\partial r}\frac{1}{r}-\frac{{\sigma }_{\theta }}{r}+</math><math>\frac{\partial {\tau }_{rz}}{\partial z}+{F}_{r}=0</math>

<math>\frac{\partial \left( r{\tau }_{rz}\right) }{\partial r}\frac{1}{r}+\frac{{\partial \sigma }_{z}}{\partial z}+</math><math>Fz=0</math>


|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|<span id='EquilibrimEq'></span>(1)
|}


Donde <math display="inline">{\sigma }_{r}</math> y <math display="inline">{\sigma }_{z}</math> son los esfuerzos en las direcciones <math display="inline">r</math> y <math display="inline">z</math> , <math display="inline">{\sigma }_{\theta }</math> en la dirección circunferencial, <math display="inline">{\tau }_{rz}</math> es el esfuerzo de corte y <math display="inline">{F}_{r}</math> y <math display="inline">{F}_{z}</math> son las fuerzas de peso propio en las direcciones <math display="inline">r</math> y <math display="inline">z</math>. Las correspondientes deformaciones quedan definidas por las siguientes ecuaciones:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{\epsilon }_{r}=\frac{\partial u}{\partial r}</math>

<math>{\epsilon }_{z}=\frac{\partial w}{\partial z}</math>

<math>{\gamma }_{rz}=\frac{\partial u}{\partial r}+\frac{\partial w}{\partial z}</math>

<math>{\epsilon }_{\theta }=\frac{u}{r}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(2)
|}


La ecuación constitutiva de acuerdo con la teoría de la elasticidad queda definida por:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>D=\frac{E}{\left( 1+\nu \right) \left( 1-2\nu \right) }\left[ \begin{matrix}1-\nu &\nu &\nu &0\\\nu &1-\nu &\nu &0\\\nu &\nu &1-\nu &0\\0&0&0&\frac{1-2\nu }{2}\end{matrix}\right]</math> 
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(3)
|}


Donde <math display="inline">E</math> es el módulo de elasticidad y <math display="inline">\nu</math>  es el coeficiente de Poisson.

En el análisis de pavimentos se utiliza el módulo de resiliencia ( <math display="inline">{M}_{r}</math>) como módulo de elasticidad y se define como:

{| class="MsoTableGrid"
 
  |<math>M_r= \frac{\sigma_d}{\varepsilon_r}</math>                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
  |
<nowiki>  </nowiki>(4)  
 
|}

Donde <math display="inline">{\sigma }_{d}=</math><math>{\sigma }_{1}-{\sigma }_{3}</math>, es el esfuerzo desviador y <math display="inline">\varepsilon_r</math>, es la deformación recuperable después de cien ciclos de carga y descarga.En los protocolos de ensayo del módulo de resiliencia [8] se utilizan combinaciones de esfuerzo de confinamiento ( <math display="inline">{\sigma }_{3}</math>) y esfuerzo desviador cíclico para obtener los ''Mr'', en distintas condiciones de esfuerzos.  Existen ecuaciones de ajuste clásicas como la siguiente:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"  
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>E={k}_{1}\, {P}_{a}{\left( \frac{\theta }{{P}_{a}}\right) }^{{k}_{2}}{\left( \frac{{\tau }_{oct}}{{P}_{a}}+1\right) }^{{k}_{3}}</math>
|}
| style="text-align: center;width: 5px;text-align: right;white-space: nowrap;" |<span id="ModuloMepdgEq"></span>(5)
|}

Donde <math display="inline">{P}_{a}</math> es la presión atmosférica, <math display="inline">\theta</math>  es la suma de los esfuerzos principales <math display="inline">=</math><math>{\sigma }_{1}+{\sigma }_{2}+{\sigma }_{3}</math>, <math display="inline">{\tau }_{oct}</math> el esfuerzo octaédrico  <math display="inline">=</math><math>1/3\sqrt{{\left( {\sigma }_{1}-{\sigma }_{2}\right) }^{2}+{\left( {\sigma }_{1}-{\sigma }_{3}\right) }^{2}+{\left( {\sigma }_{2}-{\sigma }_{3}\right) }^{2}}</math> y k<sub>1</sub>, k<sub>2</sub>, k<sub>3 </sub>son parámetros experimentales. La ecuación (<span id="cite-ModuloMepdgEq"></span>[[#ModuloMepdgEq|4]]) es propuesta por la guía empírico-mecanicista de diseño de pavimentos (MEPDG 2008). Si el módulo de resiliencia depende del estado de esfuerzos, se debe de resolver el análisis de esfuerzos de manera iterativa, correspondiendo a un análisis con elasticidad no lineal. Debido a la no linealidad de la ecuación (5), la convergencia de los análisis de esfuerzos se complica. Para los análisis presentados se usa la ecuación lineal utilizada por Huang [1], que se reproduce enseguida:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"  
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>E={E}_{0}\left( 1+\beta \, \theta \right)</math> 
|}
| style="text-align: center;width: 5px;text-align: right;white-space: nowrap;" |<span id="HuangEq"></span>(6)
|}

Donde, <math display="inline">\beta</math>  es una constante de ajuste y <math display="inline">{E}_{0}</math> es el módulo de resiliencia para el caso de que el invariante <math display="inline">\theta</math>  tenga un valor de cero. La ecuación anterior no introduce demasiado error ya que los esfuerzos de compresión que se presentan en la base son mayores a cero debido al confinamiento que ejerce la carpeta.

===2.2 Teoría de flujo de agua===

Las capas de suelo de un pavimento están conformadas de suelos compactados. El proceso de compactación condiciona a los suelos a estar en estado de saturación parcial, con tres fases (aire, agua y sólidos). Los análisis de flujo de agua en pavimentos se resuelven en términos de succión matricial o presión negativa de agua. La succión está definida como

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>s={u}_{a}-{u}_{w}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(7)
|}


Donde

<math display="inline">{u}_{a}</math> es la presión del aire que se encuentra en los poros

<math display="inline">{u}_{w}</math>, es la presión del agua que está en los poros

El flujo de agua se resuelve tomando como base la ley de Darcy, que es válida suelos parcialmente saturados y se escribe como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{v}_{x}=-{k}_{x}\frac{\partial h}{\partial x}</math>


|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(8)
|}


Donde

<math display="inline">{v}_{x}</math> es la velocidad del agua

<math display="inline">{k}_{x}</math> es la conductividad hidráulica

<math display="inline">\frac{\partial h}{\partial x}</math> es el gradiente hidráulico en la dirección <math display="inline">x</math>

<math display="inline">h=z+\frac{{u}_{w}}{{\gamma }_{w}}</math>  es la carga hidráulica total

En suelos parcialmente saturados la ecuación diferencial de flujo transitorio que gobierna el flujo de agua es la ecuación de Richards [5]. La formulación clásica se considera en dos o tres dimensiones. La versión de dos dimensiones se puede reescribir para resolverla de manera simultánea con las ecuaciones de teoría de la elasticidad axi-simétrica, antes mencionadas. Para este caso la ecuación de Richards se escribe como:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>\frac{\partial }{\partial z}\left( {k}_{z}\left( s\right) \left( 1+\frac{1}{{\gamma }_{w}}\frac{\partial {u}_{w}}{{\partial }_{z}}\right) \right) +</math><math>\frac{\partial }{\partial r}\left( \frac{{k}_{r}\left( s\right) }{{\gamma }_{w}}\frac{\partial {u}_{w}}{\partial r}\right) +</math><math>\frac{k}{r}\frac{1}{{\gamma }_{w}}\frac{\partial {u}_{w}}{\partial r}\, ={m}_{w}\frac{\partial {u}_{w}}{\partial t}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|<span id='RichardsEq'></span>(9)
|}


Donde  <math display="inline">{m}_{w}</math> es el módulo de almacenamiento y queda determinado por la ecuación:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{m}_{w}=\frac{\partial {\theta }_{w}}{\partial s}</math>


|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(10)
|}


Donde <math display="inline">{\theta }_{w}</math> se obtiene de la curva de retención del suelo. Se puede utilizar la ecuación de Van Genuchten [9], que queda definida por la siguiente ecuación:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>\theta n\left( s\right) =\frac{1}{{\left( {\left( \alpha \, s\right) }^{n}+1\right) }^{m}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(11)
|}


Donde

<math display="inline">{\theta }_{n}=\frac{{\theta }_{w}-{\theta }_{r}}{{\theta }_{s}-{\theta }_{r}}</math> el contenido de agua normalizado

<math display="inline">{\theta }_{r}</math> es el contenido volumétrico de agua residual

<math display="inline">{\theta }_{s}</math> es el contenido volumétrico saturado de agua

<math display="inline">\alpha</math>  es un parámetro de ajuste relacionado con el valor de entrada de aire (1/kPa)

<math display="inline">n,\, m</math> son parámetros de ajuste, y <math display="inline">m=1-</math><math>\frac{1}{n}</math>

La conductividad hidráulica en un suelo parcialmente saturado depende de la succión a la que se encuentra el suelo y se adopta el modelo de Van Genuchten - Mualem [9][10] para conductividad hidráulica,

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>k\left( {\theta }_{n}\right) ={k}_{s}\sqrt{{\theta }_{n}}{\left( 1-{\left( 1-{\theta }_{n}^{1/m}\right) }^{m}\right) }^{2}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|<span id='ConductivityEq'></span>(12)
|}


<span id='CacacteristicaEq'></span>Donde <math display="inline">{k}_{s}</math> es la conductividad saturada del suelo y <math display="inline">m</math> se describe en la ecuación (11) de la curva característica.

El módulo de elasticidad se puede hacer dependiente de la succión. Una de las ecuaciones de ajuste es la de Liang [11] en términos de esfuerzos efectivos,

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>E={k}_{1}\, {P}_{a}{\left( \frac{{\theta }_{b}+\chi \left( {u}_{a}-{u}_{w}\right) }{{P}_{a}}\right) }^{{k}_{2}}{\left( \frac{{\tau }_{oct}}{{P}_{a}}+1\right) }^{{k}_{3}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|<span id='LiangEq'></span>(13)
|}


Donde

<math display="inline">{\theta }_{b}={\sigma }_{1}+{\sigma }_{2}+{\sigma }_{3}={\sigma }_{1}+</math><math>2{\sigma }_{3}</math> es el invariante,

<math display="inline">{\tau }_{oct}=\frac{1}{3}\sqrt{{\left( {\sigma }_{1}-{\sigma }_{2}\right) }^{2}+{\left( {\sigma }_{1}-{\sigma }_{3}\right) }^{2}+{\left( {\sigma }_{2}-{\sigma }_{3}\right) }^{2}}</math> es el esfuerzo octaédrico de corte

<math display="inline">{\sigma }_{1},{\sigma }_{2},\, {\sigma }_{3}</math>  son los esfuerzos principales

<math display="inline">{k}_{1},{k}_{2},{k}_{3}</math> son parámetros de ajuste

<math display="inline">\chi</math>  es el parámetro de esfuerzos efectivos de Bishop que depende del grado de saturación del suelo cuyo valor oscila entre 0 y 1.

<math display="inline">{P}_{a}</math> es la presión atmosférica (101.325 kPa)

La ecuación (13<span id='cite-LiangEq'></span>) es una modificación de la ecuación (5<span id='cite-ModuloMepdgEq'></span>), la cual toma en cuenta el parámetro de esfuerzos efectivo de Bishop, <math display="inline">\chi</math> . Está ecuación fue desarrollada originalmente para suelos cohesivos.

Para los análisis que se presentan en este artículo se utiliza una modificación de la ecuación (6<span id='cite-HuangEq'></span>). Agregando el efecto de la succión por medio de la modificación lineal del parámetro <math display="inline">E{}_{0}</math> en función de la succión. La ecuación siguiente es la que se utiliza,

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{E}_{0}={E}_{s0}+{\beta }_{s}s</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|<span id='HuangMEq'></span>(14)
|}


Donde <math display="inline">{E}_{s0}</math> valor es del módulo con succión cero, <math display="inline">{\beta }_{s}</math> es una constante que da la variación del módulo con respecto a la succión.

===2.3 Análisis de pavimento y cálculo de daño===

====2.3.1 Definición del problema.====

En este caso se analiza la estructura axi-simétrica del pavimento flexible que se presenta en la <span id='cite-_Ref534280639'></span>[[#_Ref534280639|Fig. 3]]. La estructura tiene un espesor de carpeta de concreto asfáltico de 10 cm, un espesor de base de 20 cm y un espesor de subrasante de 70 cm. La carga aplicada corresponde a un eje equivalente de 8.2 ton aproximadamente. Las condiciones de frontera están esquematizadas mediante triángulos para restringir los desplazamientos en la dirección ''z'' y ''r'', y círculos para restringir los desplazamientos solo en la dirección ''r ''en el eje de simetría.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">[[File:Fig 3.png|centre|600x600px]]</div>

<div id="_Ref534280639" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig. 3. Sección de pavimento para análisis</div>

Las propiedades de las capas de base y de subrasante dependen de los estados de esfuerzos y de la succión de acuerdo con la ecuación (14<span id='cite-HuangMEq'></span>) y se presentan en la <span id='cite-_Ref534282185'></span>[[#_Ref534282185|Tabla 1]]. La carpeta asfáltica se considera que tiene un comportamiento lineal con un módulo elástico de 3,000,000 kPa y con una relación de Poisson de 0.4.

<div id="_Ref534282185" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Tabla 1. Parámetros para el análisis hidro-mecánico</span></div>

{| style="width: 41%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Parámetro</span>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Base</span>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Subrasante</span>
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math display="inline">{E}_{s0}</math> (kPa)
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|80,843
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|34,200
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\beta</math> 
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0.0025
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0.0035
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>{\beta }_{s}</math>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0.5
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|5
|}


Los parámetros de la curva de retención y conductividad hidráulica se presentan en la <span id='cite-_Ref534394538'></span>[[#_Ref534394538|Tabla 2]],

<div id="_Ref534394538" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Tabla 2. Parámetros hidráulicos usados</span></div>

{| style="width: 60%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Parámetro</span>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Carpeta</span>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Base</span>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Subrasante</span>
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>{\theta }_{s}</math>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0.2
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0.35
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0.45
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>{\theta }_{r}</math>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0.01
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0.05
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0.05
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>n</math>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|1.4
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|7.35
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|1.4
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math display="inline">\alpha</math>  (1/kPa)
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0.1
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|1.53
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0.1
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math display="inline">{k}_{s}</math> (m/hr)
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0.01
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|3
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0.1
|}


La simulación considera cuatro etapas:

:1. De 0 a 50 hrs. Equilibrio inicial con un nivel de aguas freáticas localizado a una profundidad de 5 m. La idea es considerar un equilibrio del nivel freático después de una temporada de secas, la succión que se obtiene es la condición de referencia para el daño cuando se presentan los cambios de contenido de agua subsecuentes. En la modelación, esta variación se introduce mediante una condición de frontera de succión (condición de Dirichlet) en la parte inferior, que varía con el tiempo.

:2. De 50 a 74 hrs. Se introduce una infiltración de agua en la superficie del pavimento con una velocidad de 0.85 mm/h. Se implementa con una condición de Neumann en la parte superior.

:3. De 74 a 100 hrs. Se elimina la infiltración y se equilibra nuevamente con las nuevas condiciones.

:4. De 100 a 300 hrs. Se eleva el nivel de aguas freáticas a 3.5 m de profundidad.

====2.3.2 Análisis de resultados.====

Se utilizó el programa FlexPDE 7.0® para solucionar las ecuaciones diferenciales. Es una herramienta que utiliza los elementos finitos para solucionar cualquier sistema de ecuaciones diferenciales, con sus variables dependientes, en función del tiempo y toma en cuenta la no linealidad de las ecuaciones. El programa FlexPDE utiliza un método de Newton-Rapson Modificado para la solución del sistema no-lineal; el cual toma un paso de Newton por cada paso del tiempo, y las no-linealidades se manejan adaptando el tamaño de paso en el tiempo.  Se utilizaron triángulos, con un orden de interpolación cuadrático, para la malla de elementos finitos; se tienen un total de 2737 elementos. La malla se presenta en la <span id='cite-_Ref10530824'></span>[[#_Ref10530824|Fig. 4]], se puede notar que existe un mayor número de elementos en la zona del eje de axi-simetría y las interfaces (zonas con líneas rojas). Es en la zona donde se requiere estimar las deformaciones de tensión y de compresión, además del muestreo de la presión del agua en la base y subrasante.
[[File:Draft_Chavez_729651003-image2.png|centre|600x600px]]

<div id="_Ref10530824" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Fig. 4. Malla utilizada para la solución del problema</span></div>

El problema se resolvió considerando las cuatro etapas mencionadas. Con el objetivo de ver como varían la succión y los módulos de resiliencia cuando se somete el pavimento a variaciones de las condiciones de frontera, se presentan dos gráficas donde se muestra la historia en el eje de simetría y en el centro de las capas de base y subrasante (ver <span id='cite-_Ref534283874'></span>[[#_Ref534283874|Fig. 4]] y <span id='cite-_Ref534366064'></span>[[#_Ref534366064|Fig. 5]]). Si se analiza la <span id='cite-_Ref534283874'></span>[[#_Ref534283874|Fig. 4]], se puede verificar que en la primera etapa de la simulación la succión se equilibra hasta llegar a 67.5 y 63 kPa para la base y subrasante respectivamente. En el periodo de infiltración (segunda etapa) la succión disminuye en mayor medida en la base (11.6 kPa) que en la subrasante (5.8 kPa). En la tercera etapa el suelo tiende a secarse (aumenta la succión), y por último en la etapa cuatro se eleva el nivel de aguas freáticas y como resultado, la succión de la capa subrasante queda en 43 kPa aproximadamente.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Chavez_729651003-image3.jpeg|600px]] </div>

<div id="_Ref534283874" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Fig. 5. Variación de la succión respecto al tiempo, medido a media altura de las capas de base y subrasante.</span></div>

De manera simultánea se presenta la variación de los módulos de resiliencia en los puntos antes mencionados (Ver <span id='cite-_Ref534366064'></span>[[#_Ref534366064|Fig. 5]]). Se puede notar que cuando la succión aumenta los módulos aumentan y viceversa.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Chavez_729651003-image4.jpeg|600px]] </div>

<div id="_Ref534366064" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Fig. 6. Variación del Módulo resiliente respecto al tiempo, medido a media altura de las capas base y subrasante.</span></div>

El número de repeticiones que resiste el pavimento por agrietamiento de la carpeta ( <math display="inline">{N}_{f}</math>) y por deformación permanente ( <math display="inline">{N}_{d}</math>) se puede calcular con las siguientes ecuaciones [1],

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{N}_{f}={f}_{1}{\left( {\varepsilon }_{r}\right) }^{-{f}_{2}}{\left( {E}_{1}\right) }^{-{f}_{3}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(16)
|}


{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{N}_{d}={f}_{4}{\left( {\varepsilon }_{z}\right) }^{{f}_{5}}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(17)
|}


donde <math display="inline">{\epsilon }_{r}</math> es la deformación unitaria de tensión en la fibra inferior de la carpeta, <math display="inline">{E}_{1}</math> el módulo de elasticidad de la carpeta, <math display="inline">{\epsilon }_{z}</math> es la deformación de compresión en la fibra superior de la subrasante, <math display="inline">{f}_{1},{f}_{2},{f}_{3},{f}_{4}</math> y <math display="inline">{f}_{5}</math> parámetros de ajuste. El instituto del asfalto de Estados Unidos propone que estos parámetros sean 0.0796, 3.291, 0.854, 1.365x10<sup>-9</sup> y 4.477 respectivamente [1].

En el análisis del daño que se presenta se toma como referencia la succión de equilibrio al final de la primera etapa (50 hrs), condición de referencia de temporada de secas. Se calculó el número de repeticiones que soporta la estructura de pavimento por deformación permanente ( <math display="inline">{N}_{d}</math>) y por fatiga de la carpeta ( <math display="inline">{N}_{f}</math>). Se definió el daño que causa una pasada como,

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{D}_{i}=\frac{1}{N}</math>
|}
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(18)
|}


Donde <math display="inline">N</math> toma el valor de <math display="inline">{N}_{d}</math> o <math display="inline">{N}_{f}</math> dependiendo del análisis que se realice, deformación permanente o fatiga. El daño que se causa con las condiciones iniciales se consideró como 1 y se realizaron los cálculos necesarios para ver la variación del daño por deformación permanente y agrietamiento. Los resultados de presentan en la <span id='cite-_Ref534368469'></span>[[#_Ref534368469|Fig. 7]], se observa muy claro como el daño aumenta cuando la succión disminuye (ver <span id='cite-_Ref534283874'></span>[[#_Ref534283874|Fig. 5]]). Para este caso el daño final es mayor en ambos casos: 12.5% y 5.7 % por deformación permanente y agrietamiento respecto a la condición inicial.
[[File:Draft_Chavez_729651003-chart1.svg|centre|600x600px]]

<div id="_Ref534368469" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">Fig. 7. Daño por eje equivalente relativo</span></div>

==3 Conclusiones==

La solución simultánea de las ecuaciones diferenciales de flujo y de equilibrio ayuda a la cuantificación del daño en las estructuras de pavimento de manera más precisa. El tema no es sencillo ya que las ecuaciones involucradas son no-lineales. Pero en la actualidad es posible solucionarlas con herramientas numéricas como los elementos finitos y con ello realizar mejores diseños. Para llevar estos análisis a la práctica es importante establecer las condiciones de frontera hidráulicas a las que la estructura estará expuesta, que dependen del clima y la profundidad del nivel de agua freática de la zona. Se necesita realizar más trabajo para establecer condiciones de frontera simplificadas generales para condiciones de diseño correspondientes a la práctica.

==4 Referencias==

[1] Huang, Y.H. (2004). Pavement Analysis and Design. Pearson Prentice Hall

[2] Burmister, D. M. (1945). The general theory of stresses and displacements in layered soil systems. III. Journal of applied Physics, 16(5), 296-302.

[3] Duncan, J. M., Monismith, C. L., & Wilson, E. L. (1968). Finite element analysis of pavements. Highway Research Record, 228, 18-33.

[4] Liu, P., Wang, D., & Oeser, M. (2017). Application of semi-analytical finite element method to analyze asphalt pavement response under heavy traffic loads. Journal of traffic and transportation engineering (English edition), 4(2), 206-214.

[5] Richards, L. A. (1931). Capillary conduction of liquids through porous mediums. physics, 1(5), 318-333.

[6] Oñate, E., (1995). Cálculo de estructuras por medio de elementos finitos. Análisis estático lineal. Centro Internacional de Métodos Numéricos en Ingeniería.

[7] Hamidzadeh, H. R., & Jazar, R. N. (2010). Vibrations of thick cylindrical structures. New York: Springer.

[8] AASHTO-T 307. (2003). Determining the resilient modulus of soils and aggregate materials. Standard specifications for transportation materials and methods of sampling and testing.

[9] Van Genuchten MT (1980). A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci Soc Am J, 44(5):892–8.

[10] Mualem Y (1976). A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resour Res;12(3):513–22.

<span style="text-align: center; font-size: 75%;">[11] </span>Vanapalli, S. K., & Han, Z. (2014). Application of the unsaturated soil mechanics in the design of pavements. In Proceedings of the 6th International Conference on Unsaturated Soils, Sydney, Australia (pp. 1799-1805).