SciELO - Scientific Electronic Library Online

 
vol.21 issue37Accessibility: fundamental for implementation of transportation infrastructure elementAssessment of the area of influence of San José-Sabanilla-La Campiña: Public transportation trunk corridor Applying the transit oriented development methodology author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • Have no similar articlesSimilars in SciELO

Share


Infraestructura Vial

On-line version ISSN 2215-3705Print version ISSN 2215-3705

Infraestructura Vial vol.21 n.37 San Pedro de Montes de Oca Jan./Jun. 2019

http://dx.doi.org/10.15517/iv.v21i37.42614 

Artículo

Estimación del módulo resiliente para materiales granulares de Costa Rica

Estimation of the resilient modulus for granular materials of Costa Rica

Juan Carlos Munera Miranda1 

Jose Pablo Aguiar Moya PhD.2 

1Universidad de Costa Rica, Costa Rica.juan.muneramiranda@ucr.ac.cr

2Universidad de Costa Rica, Costa Rica.jose.aguiar@ucr.ac.cr

Resumen

El módulo resiliente es una medida de la capacidad de un material de almacenar o absorber energía sin que este experimente deformación permanente. Este se determina mediante ensayos triaxiales de carga repetida y es un parámetro importante en el análisis y diseño de pavimentos flexibles. Con su valor es posible predecir la respuesta mecánica del pavimento como esfuerzos, deformaciones y deflexiones. Sin embargo, su obtención es onerosa por los equipos que involucra, razón por la que no todos los laboratorios de materiales cuentan con el ensayo. Debido a esto, los diseñadores de pavimentos recurren a correlaciones para aproximar un valor de módulo resiliente, que generalmente difiere de forma considerable del obtenido con el ensayo. La siguiente investigación tiene como objetivo la estimación de un modelo universal de módulo resiliente para materiales granulares que considere variables relacionadas con los estados de esfuerzo y las características básicas de los materiales. Para ello, se generó un modelo de regresión lineal mediante mínimos cuadrados ordinarios con datos de al menos 13 diferentes materiales granulares. El modelo se validó mediante indicadores estadísticos obteniendo una bondad de ajuste aceptable.

Palabras-clave: >Módulo resiliente; Materiales granulares; Regresión

Abstract

The resilient modulus is a measure of the ability of a material to store or absorb energy without experiencing permanent deformation. This is determined by triaxial tests of repeated load and is an important parameter in the analysis and design of flexible pavements. With its value it is possible to predict the mechanical response of the pavement as stresses, deformations and deflections. However, obtaining it is expensive for the equipment involved and not all material laboratories have the test. Because of this, pavement designers use correlations to approximate a resilient modulus value, which generally differs considerably from that obtained with the test. The following research aims to estimate a universal model of a resilient modulus for granular materials that considers variables related to stress states and the basic characteristics of the materials. For this, a linear regression model was generated using ordinary least squares with data from at least 13 different granular materials. The model was validated by statistical indicators obtaining an acceptable goodness of fit.

Key words: Resilient modulus; Granular materials; Regression

Introducción

El módulo resiliente es una medida de la capacidad los materiales para almacenar o absorber energía sin deformación permanente, es decir, es una propiedad dentro del límite proporcional del material (Popov, 1990)[1]. Se determina mediante un ensayo de respuesta dinámica basado en la deformación recuperable bajo cargas repetidas y se define como (Ecuación (1)) (Huang, 2004)[2]:

M_R=\dfrac{\sigma_d}{\epsilon_r} (1)

donde:

M_R= módulo resiliente

\sigma_d = esfuerzo desviador

\epsilon_r = deformación unitario recuparable (resiliente)

Cuando un material se somete a un ensayo de carga repetida, al inicio de la aplicación de carga, se da una deformación permanente y al mismo tiempo una deformación recuperable, conforme aumenta el número de repeticiones carga, la deformación plástica va disminuyendo y la deformación pasa a ser prácticamente toda recuperable (Huang, 2004)[2] Figure 1.

Figure 1 Deformaciones bajo un ensayo de carga repetida (Huang, 2004)[2] 

El módulo resiliente se realiza en materiales cohesivos y granulares mediante un ensayo triaxial de carga repetida, que simula los estados de esfuerzos que experimenta el material por el tránsito de los vehículos. El ensayo consiste en la aplicación de un pulso de carga haversiana de magnitud fija, con una duración de carga de 100 ms y un periodo de reposo de 900 ms, sobre un espécimen de ensayo cilíndrico (102 mm de diámetro y 203 mm de altura, generalmente). Durante el ensayo, la muestra se somete a un esfuerzo cíclico de manera dinámica y esfuerzo de confinamiento de forma estática proporcionado por medio de una cámara de presión triaxial. La respuesta total de la deformación axial resiliente (recuperable) de la muestra se mide y se utiliza para calcular el módulo resiliente, según la norma de ensayo AASHTO T 307 (AASHTO, 2019)[3]. El ensayo se ejecuta para varios estados de esfuerzo, de manera que se pueda calcular un modelo de módulo resiliente para el material que se está analizando, Ecuaciones (2, 3 y 4) (NCHRP, 2004)[4].

M_r=k_1 \cdot P_a \cdot \left( \dfrac{\theta}{P_a}\right)^k_2 \cdot \left(\dfrac{\tau_{oct}}{P_a}\right)^k_3 (2)
\theta = \sigma_1 + \sigma_2 + \sigma_3 (3)
\tau_{oct}=\dfrac{1}{3} \sqrt{(\sigma_1-\sigma_2)^2+(\sigma_1-\sigma_3)^2+(\sigma_2-\sigma_3)^2} (4)

donde:

M_R= módulo resiliente

\theta = invariante de esfuerzos (sumatoria de esfuerzos principales)

\sigma_1, \sigma_2, \sigma_3= esfuerzos principales

P_a= presión atmosférica

\tau_{oct}= esfuerzo cortante octaédrico

k_1, k_2, k_3= parámetros de regresión normalizados log-log

El módulo resiliente se utiliza en el análisis elástico multicapa, y junto con el módulo de Poisson, es una variable indispensable para predecir la respuesta mecánica de los materiales como los esfuerzos, deformaciones y deflexiones dentro de un pavimento flexible (Jiménez, 2009)[5]. Su valor es un parámetro de entrada para el diseño empírico y mecánico-empírico de pavimentos.

Al ser un ensayo costoso por los equipos que implica, no todos laboratorios de materiales convencionales disponen de este ensayo. Muchos diseñadores de pavimentos optan por el uso de correlaciones o nomogramas para determinar el valor aproximado del módulo resiliente. Con los nomogramas se estima el módulo resiliente según los resultados de otros ensayos básicos de caracterización (AASHTO, 1993)[6]. Otra forma, es mediante ecuaciones empíricas como la que propusieron Heukelom y Klomp (Heukelom y Klomp, 1962)[7] (Ecuación (5)).

M_R=1500 \cdot CBR (5)

donde:

M_R= módulo resiliente, psi

CBR = índice de soporte California, %

>El valor de 1500 de la Ecuación (5) puede variar de 750 a 3000 con un factor de 2. Provee mejores resultados con CBR menores a 20, es decir, con suelos de grano fino o subrasantes (Huang, 2004)[2]. Ninguna de estas formas de estimar el módulo resiliente da un resultado cercano al valor real que podría experimentar el material dentro de la estructura del pavimento, ya que, por ejemplo, no toma en consideración el comportamiento esfuerzo-dependiente, es decir, en condiciones reales de carga el material no exhibe un único valor de módulo resiliente, sino que esto va a depender del estado de esfuerzos bajo el cual este sometido en el pavimento.

El objetivo principal de la investigación, es la estimación de un modelo universal de módulo resiliente para materiales granulares, que considere variables relacionadas con los estados de esfuerzos y las características básicas de los materiales. Para lograr esto, se tomó una base datos con 13 materiales granulares diferentes con resultados de ensayo de módulo resiliente y caracterización básica. Con esta información, se generó un modelo de regresión lineal mediante mínimos cuadrados ordinarios (MCO). El modelo se validó mediante diferentes indicadores estadísticos y se obtuvo una bondad de ajuste aceptable.

Metodología

Para el desarrollo del modelo de predicción del módulo resiliente de materiales granulares, se analizaron 13 muestras provenientes de diferentes zonas de Costa Rica compuesta por 8 bases, 2 subbases y 3 préstamos; de las cuales se tienen datos de ensayo de: módulo resiliente, granulometría, densidad-humedad (Próctor modificado) e índice de soporte California (CBR). Para el análisis se seleccionaron inicialmente las variables mostradas en la Tabla 1. En esta misma tabla se indican los rangos de datos de cada variable para las muestras analizadas.

Table 1 Variables para análisis de regresión 

Variable Abreviatura Unidades Rango de datos
Mínimo Máximo
Módulo resiliente Mr kPa 48100 555000
Esfuerzo desviador \sigma_{d} kPa 10 966
Esfuerzo de confinamiento \sigma_{c} KPa 21 138
Invariante de esfuerzos (suma de esfuerzos principales) \theta kPa 72 1381
Contenido óptimo de humedad \%Opt % 5 15
Densidad máxima seca (Próctor modificado) Dmax kg/m3 1800 2282
Porcentaje pasando el tamiz N°4 (4,75 mm) Pas N°4 % 35 93
Porcentaje pasando el tamiz N°200 (0,075 mm) Pas N°200 % 4 16
Abertura de tamiz donde pasa el 60% D60 mm 1 18
Índice de soporte California CBR \; 95\% % 29 140

Se generó la matriz de coeficientes de correlación Pearson con el objetivo de determinar el grado de asociación lineal entre pares de variables. Este análisis estadístico se utilizó como criterio para descartar variables explicativas que estuvieran correlacionadas (esto se detalla posteriormente en la Tabla 2).

Con las variables definidas se aplicó un modelo clásico de regresión lineal mediante el método de mínimos cuadrados ordinarios (MCO) para determinar el modelo de predicción. Para que los parámetros o predictores de la regresión sean consistentes (no sesgados) y eficientes (varianza mínima) se deben cumplir los siguientes supuestos:

  • Los residuos tienen una media de cero.

  • Todos los predictores no están correlacionados con los residuos.

  • Los residuos no están correlacionados entre sí.

  • Los residuos tienen una varianza constante (homocedasticidad).

  • Ninguna variable predictora está correlacionada perfectamente (r = 1) con otra variable predictora (multicolinealidad).

  • Los residuos están distribuidos normalmente.

Los supuestos se verificaron con el análisis de gráficos de residuos, histogramas e inferencia por medio de pruebas de hipótesis. Para considerar el no cumplimiento del supuesto de varianza constante de las perturbaciones se realizó una corrección de White de los errores estándar por heterocedasticidad con el fin de tener una inferencia robusta. Para garantizar la linealidad de las variables se aplicó una transformación logarítmica ya que algunas variables tienen un comportamiento no lineal.

Se evaluó la bondad de ajuste del modelo, es decir, el grado en que las variables independientes analizadas explicaban la variable dependiente (módulo resiliente) por medio del coeficiente de determinación ajustado.

En la Figure 2 se muestra el esquema que se siguió para validar el modelo de regresión.

Figure 2 Esquema para la estimación de modelo 

Resultados y análisis

Linealidad

Para garantizar la linealidad del modelo se realizó una transformación de variable de tipo logarítmica, debido a que el módulo resiliente con respecto al invariante de esfuerzos o el esfuerzo desviador, puede variar de forma no lineal. Por lo tanto, el modelo no lineal propuesto inicialmente se indica en la Ecuación (6). Los β corresponden a los coeficientes de regresión.

Mr = \beta_{0} \cdot \sigma_{d}^{\beta_{1}} \cdot \sigma_{c}^{\beta_{2}} \cdot \theta^{\beta_{3}} \cdot \%Opt^{\beta_{4}} \cdot Dmax^{\beta_{5}} \cdot Pas3/4"^{\beta_{6}} \cdot PasN°4^{\beta_{7}} \cdot PasN°^{\beta_{8}} \cdot D60^{\beta_{9}} \cdot CBR95%^{\beta_{10}} (6)

Aplicando la transformación de variable (log-log) el modelo lineal sería (7):

\log{Mr}=\log{\beta_{1}}+\beta_{2} \log{\sigma_{d}} + \beta_{3} \log{\sigma_{c}} + \beta_{4} \log{\theta} + \beta_{5} \log{\%opt} + \beta_{6} \log{Pas3/4"} + \beta_{7} \log{PasN°4} + \beta_{8} \log{PasN°200} + \beta_{9} \log{D60} + \beta_{10} \log{CBR95\%} (7)

Correlación Lineal entre Variables

Para seleccionar de la Table 1 las variables independientes a ser incluidas en el análisis de regresión, se calcularon los coeficientes de correlación entre todas las variables (Table 2). Un coeficiente de correlación mayor a 0,70, se considera como una correlación significativa entre variables (Holicky, 2013)[8].

Table 2 Coeficientes de correlación Pearson entre variables 

Variable Mr \sigma_{d} \sigma_{c} \theta \%Opt Dmax Pas \\ N°200 Pas \\ N°4 Pas \\ 3/4°" D60 CBR \\ 95\%
Mr 1,00
\sigma_{d} 0,79 1,00
\sigma_{c} 0,87 0,60 1,00
\theta 0,93 0,89 0,88 1,00
\%Opt -0,18 0,03 0,00 0,02 1,00
Dmax 0,15 -0,05 0,00 -0,04 -0,91 1,00
Pas \\ N°200 -0,19 0,02 0,00 0,02 0,37 -0,29 1,00
Pas \\ N°4 -0,13 0,04 0,01 0,03 0,37 -0,25 0,88 1,00
Pas \\ 3/4°" -0,13 0,05 0,01 0,04 0,33 -0,26 0,72 0,86 1,00
D60 0,09 -0,03 0,00 -0,02 -0,34 0,25 -0,73 -0,90 -0,81 1,00
CBR \\ 95\% 0,07 0,00 0,00 0,00 -0,01 -0,08 -0,59 -0,72 -0,46 0,74 1,00

Se puede observar que las variables de esfuerzo desviador y esfuerzo de confinamiento están correlacionadas con el invariante de esfuerzos de forma positiva, esto se debe a que el invariante de esfuerzos es la suma de esfuerzos principales y estos incluyen al esfuerzo desviador y el esfuerzo de confinamiento. La densidad máxima seca y el contenido óptimo de humedad están correlacionadas de forma negativa, esto es esperable ya que una es función de la otra. Con respecto a las variables granulométricas todas entre ellas están correlacionadas positivamente. Esto es razonable ya que los porcentajes pasando son sumas acumuladas y el tamaño por el que pasa el 60%, es dependiente de los porcentajes pasando. El índice de soporte California se correlaciona negativamente con el porcentaje pasando el tamiz N°4 y positivamente con el tamaño por el que pasa el 60%, esto es esperable, ya que, a mayor contenido de partículas gruesas, mayor valor de CBR.

Con respecto a la correlación del módulo resiliente con las demás variables, este tiene una correlación alta positiva con respecto el esfuerzo desviador, el esfuerzo de confinamiento y el invariante de esfuerzos, siendo este último el que tiene mayor coeficiente de correlación. Se observa una correlación menor del módulo resiliente con las restantes variables.

Entre las variables asociadas al ensayo de Próctor modificado, el contenido óptimo de humedad tiene una mayor correlación con el módulo resiliente. Entre las variables granulométricas, el porcentaje pasando el tamiz N°200 es el que tiene mayor coeficiente de correlación. El índice de soporte California tiene un bajo coeficiente de correlación con el módulo resiliente.

Del análisis realizado se seleccionaron para generar el modelo de regresión las variables independientes: invariante de esfuerzos, contenido óptimo de humedad y porcentaje pasando el tamiz N°200. Con esta selección, se elimina la multicolinealidad entre las variables explicativas.

Generación de Modelo

Con las variables independientes escogidas y bajo el cumplimiento de los supuestos, se realiza el modelo de regresión lineal múltiple con el método de mínimos cuadrados ordinarios. En la Table 3 se presenta un resumen con los resultados obtenidos del análisis de regresión lineal múltiple realizado con el software Microsoft Excel.

Table 3 Resumen de resultados de regresión lineal múltiple 

Variable Coeficientes Error estándar (S.E.) t sutdent P value (normal) P value (t student)
Intercepto 3,98927 0,04172 95,62702 0,00E+00 7,989E-248
\log{\theta} 0,67350 0,01168 57,64553 0,00E+00 1,432E-177
\log{\% Opt} -0,20049 0,02968 -6,75529 1,43E-11 6,149E-11
\log{PasN°200} -0,25639 0,02458 -10,43108 0,00E+00 2,760E-22
Fcrítico 1c: 5,5318E-179 tcrítico 2c: 1,960
SSR SSE SST F R2 R2 ajustado
17.92171 1,73011 19,65182 1173,98626 0,91196 0,91118

Análisis de Residuales

Se verificaron los supuestos que especifican el comportamiento de los residuos (diferencia entre los valores observados y los predichos): no correlación, normalidad y varianza constante para todas las observaciones (homocedasticidad). En la Figure 3, se muestra el gráfico de residuos vs los valores ajustados, se puede apreciar que inicialmente se cumple el supuesto de homocedasticidad ya que los residuos se ubican de forma aleatoria a ambos lados del cero, con patrones no detectables en los puntos. A su vez, en la Figure 4 se observa la variación de los residuos con el orden de los datos, se puede ver que no hay tendencias o patrones cuando se muestran en orden cronológico los residuales, esto indica que se cumple con el supuesto de independencia entre residuos (Montgomery & Runger, 2003)[9].

Figure 3 Residuales con respecto a los valores ajustados 

Figure 4 Residuales con respecto al orden de datos 

Para verificar que los residuos están distribuidos normalmente, se elaboró un histograma de frecuencias y un gráfico de normalidad. En la Figure 5 se observa que los residuales tienden a seguir una distribución normal por la forma simétrica de campana de Gauss, que muestra la distribución de frecuencias. Sin embargo, en la Figure 6 se advierte que los residuos no están distribuidos normalmente puesto que no todos los puntos caen sobre la línea recta, el gráfico tiene forma de S que implica una distribución colas, esto puede indicar que existen pocas o demasiadas observaciones en las colas de la distribución.

Figure 5 Histograma de frecuencia de los residuales 

Figure 6 Gráfico de probabilidad normal de los residuos 

La no normalidad de los residuos no implica que los estimadores MCO no sean insesgados (consistentes) y de varianza mínima (eficientes) (Gujarati & Porter, 2009)[10]. Con respecto a la inferencia, los estadísticos t y F tienen distribuciones que son aproximadamente distribuciones t y F, cuando las muestras son de tamaño grandes (n - k > 30) (Benavente, Otero, & Vásquez, 2007)[11]; es decir, que si se cumple con esta condición, la inferencia estadística sigue siendo confiable, a pesar del no cumplimiento de la normalidad de los residuos3.

Análisis de Estimadores

El intercepto tiene un coeficiente de regresión mayor que el de las variables explicativas, esto implica que una porción importante del modelo está siendo explicado por intercepto. Esto podría indicar que modelo se podría mejorar incorporando otras variables explicativas. El invariante de esfuerzos tiene un coeficiente positivo, lo cual es esperable, ya que conforme aumenta el invariante de esfuerzos aumenta el módulo resiliente. El coeficiente del contenido óptimo de humedad tiene signo negativo. Esto también es razonable, ya que se espera que el módulo resiliente disminuya conforme mayor sea el contenido de humedad óptimo del material. De igual manera, el coeficiente del porcentaje pasando el tamiz N°200 tiene signo negativo, este comportamiento es aceptable, ya que se espera que un material con mayor porcentaje pasando el tamiz N°200 tenga un menor módulo resiliente.

Los valores de t-student calculados para cada estimador, son mayores al valor crítico de 1,96 para un nivel de confianza del 95%, esto quiere decir, que se rechaza la hipótesis nula H0 de que alguno de los estimadores sea igual a cero. De igual forma, los p-value menores al nivel de significancia α = 0,05 indican que estadísticamente ninguno de los estimadores considerados en el modelo sea igual a cero. El valor de F calculado es mayor que el crítico, lo que indica que se rechaza la hipótesis nula H0 de que todas pendientes del modelo son iguales a cero.

Análisis de Bondad de Ajuste

La suma de cuadrados de la regresión (SSR) es mayor a la suma de cuadrados de errores (SSE), esto implica que la varianza de la variable dependiente se explica en su mayor parte por las desviaciones de las variables independientes y no por las desviaciones de los residuos. Este comportamiento es el deseado, y es por esto, que el coeficiente de determinación R2 es cercano a 1. El R2 ajustado brinda un valor más confiable de la bondad de ajuste del modelo ya que el R2 tiende a sobreestimar el poder predictivo del modelo, aunque las variables no sean significativas.

En la Figure 7 se presenta un gráfico donde se comparan los módulos resilientes estimados vs los medidos, si los puntos pasan por la línea a 45°, implica que el modelo propuesto es efectivo prediciendo la variable dependiente y no presenta un sesgo; el modelo propuesto se ajusta satisfactoriamente a los datos medidos, no así la correlación de la Ecuación (5) comúnmente utilizada para estimar el módulo resiliente (Heukelom & Klomp, 1962)[7], que presenta un sesgo importante y tiende a sobreestimar su valor.

Figure 7 Módulo resiliente estimado vs el medido 

Corrección por Heterocedasticidad

Suponiendo que no se cumpliera el supuesto de homocedasticidad (varianza constante de los residuos), se realizó una corrección de heterocedasticidad de White de los errores estándar para las inferencias estadísticas fueran asintóticamente válidas (es decir, para muestras grandes), sobre los verdaderos valores de los parámetros (Gujarati & Porter, 2009)[10]. En la Table 4 se resumen los resultados.

Table 4 Resumen de resultados de regresión lineal múltiple con corrección de White 

Variable Coeficientes Error estándar (S.E.) t sutdent P value (normal) P value (t student)
Intercepto 3,98927 0,04210 94,76426 0,00E+00 1,558E-246
\log{\theta} 0,67350 0,01273 52,91893 0,00E+00 3,344E-166
\log{\% Opt} -0,20049 0,02552 -7,85624 4,00E-15 5,251E-14
\log{PasN°200} -0,25639 0,02196 -11,67264 0,00E+00 1,057E-26
SSR SSE SST R2 R2 ajustado
17,92171 1,73011 19,65182 0,91196 0,9118

Los valores t student calculados son menores al t crítico igual a 1,960 para un nivel de confianza del 95%, por lo tanto, hay suficiente evidencia estadística para rechazar la hipótesis nula de que uno de los estimadores β sea igual cero. Los t student calculados de 3 comparados con los de la Table 4 (sin corrección de White), no son apreciablemente diferentes, por lo tanto, el supuesto de homocedasticidad se puede confirmar nuevamente.

Conclusiones y recomendaciones

De los datos analizados y de la validación del modelo de regresión propuesto para el módulo resiliente se obtuvo que:

  • Para la base de datos y las variables consideradas inicialmente en el estudio, se determinó que las variables más significativas estadísticamente y que mejor describen el modelo de módulo resiliente son: el invariante de esfuerzos (q), el contenido óptimo de humedad (%Opt) y el porcentaje pasando el tamiz N°200 (PasN°200).

  • Se verificaron los supuestos y cumplimientos para la regresión por mínimos cuadrados ordinarios, excepto la normalidad de los residuos, sin embargo, esto no afecta la consistencia y eficiencia, ni tampoco la inferencia del estimador al tener una muestra grande. Este efecto sobre la normalidad de los residuos se puede mitigar aumentando el tamaño de la muestra.

  • El modelo final propuesto se muestra en la Ecuación (8):

Mr=9755.95975 \cdot \theta^{0.67350} \cdot \%Opt^{-0.20049} \cdot PasN°200^{-0.25639} (8)
  • El modelo que se propone posee una bondad de ajuste satisfactoria, y se puede verificar con base en la muestra analizada que la correlación MR = 1500 × CBR, presenta un sesgo importante y tiende a sobrestimar los valores de módulo resiliente y no toma en cuenta el comportamiento esfuerzo-dependiente del material.

  • Se recomienda aumentar la base de datos con más materiales para ampliar el espectro de materiales que incluye el modelo y así mejorar su confiabilidad.

  • Aunque no se incluyó en el análisis debido a que algunas muestras carecían del dato, sería importante evaluar variables como límites de consistencia, compresión inconfinada y parámetros de resistencia al corte, entre otras variables que puedan describir mejor el modelo.

  • El módulo resiliente varía según el contenido de humedad, por lo tanto, se recomienda evaluar esta variable de tal manera que el modelo contemple cambios de humedad a los que podría estar expuesto el material durante su tiempo de servicio.

Reconocimientos

Los ensayos y resultados de laboratorio fueron generados en el Laboratorio Nacional de Materiales y Modelo Estructurales (LanammeUCR) de la Universidad de Costa Rica.

REFERENCIAS

AASHTO. (1993). AASHTO Guide for design of pavement structures. Washigton, D.C.: American Association of State Highway and transportation Officials. [ Links ]

AASHTO. (2019). Determining the resilient modulus of soils and aggregate aaterials. Washington, D.C.: American Association of State Highway and Transportation Officials. [ Links ]

Benavente, J., Otero, A., y Vásquez, J. (2007). Econometría I. Santiago: Universidad de Chile. [ Links ]

Gujarati, D., y Porter, D. (2009). Econometría. México, D.F.: McGraw-Hill Interamericana Editores, S.A de C.V. [ Links ]

Heukelom , W., y Klomp, J. (1962). Dynamic testing as a means of controlling pavements during and after construction. Proceedings, (1st) International Conference on the Structural Design of Asphalt Pavements, 667-685. [ Links ]

Holicky, M. (2013). Introduction for probability and statistics for engineers. Prague: Springer-Verlag Berlin Heidelberg. [ Links ]

Huang, Y. (2004). Pavement analysis and design. New Yersey: Pearson Education, Inc. [ Links ]

Jiménez, M. (2009). Implementación del ensayo de módulo resiliente en bases granulares para pavimentos. San José: LannammeUCR. [ Links ]

Montgomery, D., y Runger, G. (2003). Applied statistics and probability for engineers. New York: John Wiley y Sons, Inc. [ Links ]

Popov, E. (1990). Engineering mechanics of solids. New Yersey: Prentice-Hall, Inc. [ Links ]

NCHRP. (2004). Laboratory determination of resilient modulus for flexible pavement design. Washigton, D.C.: National Cooperative Highway Research Program. [ Links ]

Recibido: 15 de Enero de 2020; Aprobado: 17 de Marzo de 2020

Creative Commons License Este es un artículo publicado en acceso abierto bajo una licencia Creative Commons