POTOSÍ

MODELOS MATEMÁTICOS GENERALES PARA LA PREDICCIÓN DEL CBR

(CALIFORNIA BEARING RATIO) EN LOS SUELOS BOLIVIANOS

Tarija, Bolivia

Alberto Benítez Reynoso i

Resumen

Muchas variables, correspondientes a fenómenos geotécnicos, se determinan mediante ensayos de

laboratorio, desde los más simples y económicos, hasta los más complejos y costosos, algunos de los
cuales pueden durar, relativamente, bastante tiempo. Para citar un ejemplo, la determinación CBR, en
laboratorio, toma, como mínimo 4 días, otros ensayos pueden durar más tiempo.
Los ingenieros, a menudo, debemos predecir o estimar algunas propiedades de los suelos cuando hay
poca o ninguna información proveniente de ensayos de laboratorio. Entonces, se plantean las siguientes
preguntas: ¿Cuáles son las leyes de variación de la capacidad portante de los suelos bolivianos,
expresadas en términos del CBR, en función de otra u otras variables, cuya determinación es más
simple?; ¿Cómo se expresan esas leyes en términos de modelos matemáticos y cuál es el grado de
correlación entre las diferentes variables?; ¿Cómo estimar la capacidad portante de los suelos
bolivianos (expresada en términos del CBR) sin recurrir a ensayos de laboratorio, conociendo otra u
otras propiedades, cuya determinación es más simple y toma mucho menos tiempo?; ¿Qué
implicaciones tienen los modelos citados en los proyectos y obras pertenecientes al campo de la
ingeniería de vial?
En consecuencia, en este trabajo de investigación, se determinan las leyes de variación de la capacidad
portante de los suelos, expresada en términos del CBR y representadas por relaciones matemáticas
con otras variables geotécnicas, con referencia particular a los suelos de las diferentes “provincias
fisiográficas” de Bolivia, se formulan modelos matemáticos, en términos de ecuaciones, que permitan
estimar el CBR sin recurrir a ensayos de laboratorio, conociendo otra u otras variables, cuya
determinación es más simple y/o más rápida y se establecen las implicaciones y/o utilidad que tienen
los modelos citados en el contexto de la práctica, cuando se diseña estructuras de pavimentos con
valores del CBR predichos con los modelos, en comparación con los diseños usando valores del CBR
de laboratorio.
El trabajo de investigación (que está sustentado en más de 20000 ensayos de laboratorio y de campo,
perteneciente a los proyectos más importantes de Bolivia, en varios de los cuales el autor ha
desempeñado diferentes funciones), consta de dos fases: una empírica (experimental) y otra teórica
(racional). La fase experimental está referida a los ensayos mencionados, en tanto que, la fase teórica
consiste fundamentalmente en la formulación y validación de los modelos matemáticos simples y
multivariados, mediante los cuales se expresa una variable dependiente o propiedad del suelo que se
desea predecir (CBR), en función de otra u otras, consideradas como variables independientes.
Como metodologías científicas generales, se usan los métodos inductivo y deductivo y, como métodos
específicos, el análisis multivariado, expresado como regresión múltiple. Al final, se formulan las
conclusiones y recomendaciones respectivas.

Palabras Clave: Modelos matemáticos, CBR (California Bearing Ratio), cuadrados mínimos,
fenómenos geotécnicos, variables geotécnicas e idoneidad de modelos.

1
 POTOSÍ

1. INTRODUCCIÓN Y OBJETIVOS: no es posible realizar tales ensayos para un

determinado proyecto o cálculo, porque no se
En el afán de comprender y explicar la tiene el equipo apropiado, porque toma mucho
naturaleza y sus fenómenos, la humanidad ha tiempo (considerando la urgencia del
creado un cuerpo de ideas llamado ciencia. La proyecto), en fin, porque es muy costoso. Estas
primera (naturaleza) existe razones, entre otras, conducen a formular el
independientemente de la humanidad, en tanto PROBLEMA a través de las siguientes
que la ciencia es una creación del hombre, de preguntas, en relación a la predicción o
la humanidad. estimación del CBR, a saber:
Una pequeña parte de esa naturaleza, la ¿Cuáles son las leyes de variación de la
constituye el suelo, que es objeto de estudio de capacidad portante de los suelos bolivianos,
la ciencia denominada Mecánica de Suelos y expresadas en términos del CBR, en función de
cuyas aplicaciones a la Ingeniería son objeto de otra u otras variables, cuya determinación es
la denominada Ingeniería Geotécnica. más simple?
En ese contexto, los ingenieros geotécnicos y ¿Cómo se expresan esas leyes en términos de
los científicos del suelo, han desarrollado una modelos matemáticos y cuál es el grado de
variedad de métodos, procedimientos, equipos correlación entre las diferentes variables?
(de laboratorio y campo) y modelos ¿Cómo estimar la capacidad portante de los
(conceptuales, físicos, matemáticos, etc.), que suelos bolivianos (expresada en términos del
permiten el estudio de los diferentes CBR) sin recurrir a ensayos de laboratorio,
fenómenos que ocurren en el suelo y sus conociendo otra u otras propiedades, cuya
variables inherentes. determinación es más simple y toma mucho
Un fenómeno geotécnico, suele expresarse menos tiempo?
mediante variables que lo representan. Así, una ¿Qué implicaciones tienen los modelos citados
de las formas de expresar la resistencia de un en los proyectos y obras pertenecientes al
suelo (fenómeno mecánico), es mediante el campo de la ingeniería de pavimentos?
esfuerzo de compresión no confinada Consecuentemente, los OBJETIVOS del
(variable). Análogamente, la capacidad trabajo de investigación se pueden precisar en
portante de un suelo (fenómeno mecánico), se los siguientes términos:
expresa en términos de una variable (1) Determinar las leyes de variación de la
denominada CBR (California Bearing Ratio o capacidad portante de los suelos,
Valor Soporte California), de mucha utilidad expresada en términos del CBR y
en el cálculo de la estructura de los pavimentos. representadas por relaciones matemáticas
En síntesis, los fenómenos geotécnicos quedan con otras variables geotécnicas, con
representados por sus variables inherentes que referencia particular a los suelos de las
los caracterizan. diferentes “provincias fisiográficas” de
La mecánica de suelos es una ciencia Bolivia.
sustentada en la experimentación y muchas de (2) Formular modelos matemáticos, en
sus variables se determinan mediante ensayos términos de ecuaciones, que permitan
de laboratorio, desde los más simples y estimar el CBR sin recurrir a ensayos de
económicos, hasta los más complejos y laboratorio, conociendo otra u otras
costosos, algunos de los cuales pueden durar variables, cuya determinación es más
bastante tiempo. Solo la determinación del simple y más rápida.
CBR, en laboratorio, toma, como mínimo 4 (3) Determinar las implicaciones y utilidad
días, otros ensayos pueden durar más tiempo. que tienen los modelos citados en el
En ingeniería geotécnica es frecuente la contexto de los diseños de estructurales de
necesidad de predecir variables de ciertos pavimentos, cuando se los diseña con
fenómenos, cuando hay poca o ninguna valores del CBR predichos con los
información proveniente de ensayos de modelos, en comparación con los diseños
laboratorio o de campo. En muchas ocasiones, usando valores del CBR de laboratorio.

2
 POTOSÍ

2. NOVEDAD CIENTÍFICA : funciones y su comportamiento cuando las

variables adoptan valores extremos
Bolivia es un país diverso en varios aspectos, (máximos y mínimos).
pero, particularmente lo es en relación a su b) La utilidad o significación práctica radica
geografía, en especial con relación al origen y en la solución de uno de los problemas
composición de sus suelos, mismos que planteados anteriormente, en términos de
presentan propiedades geotécnicas distintas, preguntas. Es decir, los modelos permiten
según la zona geográfica o “provincia predecir (calcular) el CBR para diferentes
fisiográfica o geológica” a la que pertenecen. necesidades prácticas, sin la necesidad de
En general, en Bolivia no existen antecedentes realizar el ensayo de laboratorio, mismo
teóricos y prácticos con relación a la que toma, al menos, cuatro días de tiempo.
formulación de modelos teórico-matemáticos, Asimismo, una vez calculado los valores
expresados en términos de ecuaciones, que del CBR, usando los modelos matemáticos
permitan la predicción de variables que correspondan, según la provincia
pertenecientes a diferentes fenómenos fisiográfica o geológica a la que pertenece
geotécnicos. Específicamente, no se han el suelo, se pueden aplicar para diferentes
registrado, publicado, ni patentado modelos fines, por ejemplo, para el dimensionado
matemáticos que permitan predecir el Valor de estructuras de pavimentos, problema
Soporte California, conocido como CBR, en frecuente en los ingenieros de carreteras.
función de otras variables geotécnicas, cuya c) Además de los mencionados en a) y b), hay
determinación es más simple y/o más rápida, un otro elemento, que sustenta, aún más, la
por ejemplo, los límites de consistencia, utilidad práctica de los modelos, que
propiedades granulométricas y otras. consiste en la calificación de la
Consecuentemente, la novedad científica, de la “idoneidad” de los mismos, en términos de
investigación se resume en los siguientes predicción, la cual se determina al realizar
elementos: la comparación de los valores del CBR
a) Desde el punto de vista teórico, se destaca observados (medidos experimentalmente)
el uso de la matemática avanzada (cálculo con los valores predichos (calculados con
diferencial y análisis matricial), mediante los modelos).
la cual formulan los modelos, expresados
en términos de ecuaciones, que permiten 3. FUNDAMENTOS EMPÍRICOS Y
expresar leyes de variación de un TEÓRICOS; EL ENSAYO DEL
fenómeno geotécnico (capacidad portante, CBR Y LA MODELACIÓN
identificado por su variable CBR), en MATEMÁTICA:
función de otros fenómenos, como la
plasticidad, expresada en términos de los Tal y como refiere el autor (Benítez, 2014), la
límites de consistencia (límite líquido, Ingeniería Civil representa una síntesis entre lo
límite plástico e índice de plasticidad), el empírico y racional o teórico, es decir, se nutre
índice de grupo, etc., de más fácil y simple de ambas corrientes del pensamiento científico.
determinación experimental que el CBR. Muchas teorías se han generado a partir de
Los modelos se formulan para cada una de mediciones empíricas (experimentales) y
las seis provincias fisiográficas o muchas situaciones empíricas han sido
geológicas de Bolivia, a saber: Cordillera consecuencia de teorías. En el caso de la
Oriental, Altiplano, Cordillera Oriental, Mecánica de Suelos e Ingeniería Geotécnica, la
Sub Andino, Llanura Chaco-Beniana y situación no es tan distinta. Por ejemplo, se
Escudo Brasileño. tiene una teoría de la consolidación que se
Pero, el uso de la matemática no solo se sustenta empíricamente, pero, a su vez, la
limita a la formulación de los modelos, teoría formulada genera consecuencias que
sino que, el cálculo diferencial se usa para pueden ser probadas empíricamente.
interpretar esos modelos en términos de Análogamente, se tiene lo que sucede en

3
 POTOSÍ

relación a los modelos matemáticos que haciendo alusión a que, en el ámbito de la

expresan leyes de variación (teóricas) entre ciencia, una hipótesis científica es una
algún fenómeno (capacidad portante) y otros, afirmación provisional que surge como una
por ejemplo, plasticidad de los suelos. Es decir, explicación o fundamentación explicativa del
también, en el caso de la capacidad portante de fenómeno observado y estudiado (Gutierrez,
un suelo, expresada en términos de la variable 1996). Cuando la hipótesis logra ser
CBR, hay una combinación (síntesis) entre lo confirmada, se convierte en una tesis o ley
empírico (experimental) y lo teórico (racional). científica; esta confirmación puede llevarse a
Por esa razón, a continuación, se fundamenta el cabo mediante observaciones, mediciones
trabajo, primero desde el punto de vista experimentales (métodos empíricos), por
experimental, en relación al ensayo del CBR y, métodos racionales (raciocinio) y/o por una
en segundo término, se esboza la teoría combinación de más de uno de ellos. Es decir,
(matemática). una vez planteado el problema, como se lo ha
Sobre el ensayo experimental del CBR, se ha hecho en la introducción, se busca una
escrito bastante, por lo que, simplemente se solución, a menos que se demuestre que se
hace referencia a las ideas esenciales, a saber carece de tal. El camino para el hallazgo de la
(Carter y Bentley, 1991): solución depende de la naturaleza del
- Fue originalmente desarrollado por la problema.
División de Carreteras de California en Al formular una hipótesis que, probablemente,
1930 y años siguientes, como parte de un se convierta en tesis o ley científica, se trata de
estudio de las fallas en los pavimentos. relacionar dos o más fenómenos o variables
Posteriormente su uso fue extendido a las involucradas en estos. Por eso es que, una ley
sub rasantes y, actualmente, es usado científica puede definirse como una relación
ampliamente en el diseño de pavimentos a constante entre dos o más hechos, fenómenos o
lo largo de todo el mundo. sus variables inherentes (Bunge, 2011 y
- Durante el ensayo experimental, se hace Gutierrez, 1996).
penetrar un pistón en el suelo, el cual es Hay varios fenómenos en la Mecánica de
contenido en un molde estándar, a una tasa Suelos e Ingeniería Geotécnica, que pueden ser
o velocidad especificada de penetración. expresados mediante sus variables más
- El ensayo prosigue hasta una penetración importantes, consecuentemente, podemos
de 5 mm y se traza la curva de esfuerzos relacionar, mediante una hipótesis y/o ley,
(cargas) en función de las deformaciones estos fenómenos o sus variables relevantes.
(hundimientos). Las magnitudes se llaman también variables
- La curva carga-deformación es comparada cuantitativas o numéricas, porque su campo de
con la obtenida para un material estándar variabilidad es un conjunto de números. Este es
triturado. el tipo de variables con las que trata
- Se llama Valor Soporte California o CBR habitualmente el ingeniero. Se necesita algunos
(por sus siglas en inglés), al tanto por conceptos más para caracterizar el de ley, a
ciento de presión (esfuerzo) ejercida por el saber: variable independiente, variable
pistón sobre el suelo, para una penetración dependiente y parámetros. Sea la expresión
determinada, en relación con la presión (modelo matemático) de la forma:
correspondiente a la misma penetración en
una muestra patrón (el mencionado Y = o + 1X (1)
material triturado).
- Los detalles del ensayo han sido X e Y suelen llamarse la variable
estandarizados y están dados en las independiente y la variable dependiente,
especificaciones estándar ASTM D1883- respectivamente, mientras que  o y 1 son
07e2. parámetros del modelo. En la ciencia, la
Con relación a la fundamentación teórica, variable independiente es, a menudo, la
científico-matemática, se puede comenzar variable de control, o sea, la variable a la cual

4
 POTOSÍ

pueden atribuirse valores o cambiarlos a No se puede terminar esta breve

voluntad dentro de ciertos límites. Esta fundamentación teórica sin plantear dos
distinción tiene una raíz filosófica, ya que, los cuestiones importantes, a saber:
cambios en los valores de la variable de control  Si se plantea la hipótesis de que dos
suelen llamarse causas, mientras que los variables, que reflejan un fenómeno real de
cambios que se producen en la variable la Mecánica de Suelos, están relacionadas
dependiente se llaman efectos. En fin, se da el mediante la ley (1) u otra más compleja
nombre de parámetro a una variable cuyo valor desde el punto de vista matemático, ¿cómo
no cambia por el cambio de los valores de las determinar los parámetros o y 1
demás variables; en el ejemplo de la ecuación considerando la información empírica
anterior, o y 1 son parámetros porque sus (mediciones experimentales) existente?.
valores se mantienen constantes en un contexto  ¿Hay experiencias teóricas y prácticas que
dado, que se determinan con base en la permitan sustentar la posibilidad de dar las
información experimental. respuestas más apropiadas a nuestras
Adicionalmente, para todo par de valores de o preguntas científicas?
y 1, la ecuación (1) puede considerarse como La primera de las preguntas tiene que ver con
la representación analítica de una recta infinita la teoría estadística, asociada con la
en el plano (X,Y). Si o y 1 toman valores metodología que se plantea en el próximo
diferentes, se obtienen un número infinito de apartado, tanto para modelos unidimensionales
esas rectas. En consecuencia, cada miembro de o simples (una variable dependiente y una
este conjunto puede entenderse como una ley variable independiente), como para modelos
individual. múltiples (una variable dependiente y varias
Hasta aquí, la (1) aún no es una ley científica variables independientes).
porque sus variables tienen una interpretación Con relación a la segunda pregunta, se puede
solamente matemática. Únicamente si se afirmar, categóricamente, que hay varias
interpretan al menos las variables X e Y, no investigaciones teórico-prácticas que dieron
solo como números cualesquiera, sino como resultados muy importantes sobre problemas
variables numéricas de propiedades de algún análogos a esta investigación. Así:
sistema real, puede convertirse en una ley - En Carter y Bentley (1991) y Magnan
científica (Bunge, 2011). Esta es una cuestión (1982), se mencionan una importante
clave en el planteamiento del problema y las cantidad de correlaciones y modelos
respuestas a las preguntas científicas hechas en teóricos de regresión entre diferentes
la introducción, pues, un modelo matemático propiedades geotécnicas de los suelos.
como (1) y otros más complejos, como varios También se hace énfasis en la aplicación e
de los obtenidos como resultados en esta importancia de los métodos estadísticos y
investigación, reflejan relaciones entre probabilistas en el campo de la ingeniería
variables que toman valores numéricos y geotécnica.
representan propiedades de un sistema de la - Son particularmente relevantes las
realidad, que son los fenómenos de la conocidas ecuaciones que permiten
Mecánica de Suelos e Ingeniería Geotécnica, estimar el coeficiente de consolidación en
que se han mencionado en el planteamiento del función del límite líquido (Terzaghi y
problema. Entonces, se puede definir a una ley Peck, 1967) y varias otras (mencionadas en
científica como una hipótesis científica Ameratunga y otros, 2016).
confirmada, que establece una relación - No se puede dejar de mencionar las
constante entre dos o más variables, cada una correlaciones y modelos de regresión
de las cuales representa una propiedad de propuestos por Casagrande, entre el índice
sistemas reales concretos, al menos parcial e de plasticidad y el límite líquido, que han
indirectamente; en este caso, el sistema real y sido usados en su carta de plasticidad, base
concreto, además de natural, es el suelo. del sistema unificado de clasificación de
suelos (Carter y Bentley, 1981)

5
 POTOSÍ

- Agarwal y Ghanekar (1970); Patel y Desai matemáticos), lo que, sin duda alguna,
(2010); Vinod y Reena (2008), han constituye el llamado método inductivo; pero
propuesto modelos de regresión entre el ¿Cómo se usa el método inductivo en esta
CBR y otras propiedades para suelos de la investigación?; de la manera más sencilla, se
India. siguen las siguientes etapas:
- En fin, el autor (Benítez, 2012) ha - Mediciones experimentales que expresen
formulado varios modelos de regresión la variación de un fenómeno geotécnico
para algunos suelos bolivianos. (capacidad portante de los suelos en
En suma, muchos problemas de la Mecánica de términos de la variable CBR), en función
Suelos e Ingeniería Geotécnica se han resuelto, de otros fenómenos, tales como la
precisamente, formulando correlaciones y plasticidad, expresada en términos de las
modelos de regresión entre las distintas variables conocidas como límites de
variables que representan algunos fenómenos consistencia y otras.
de esta rama de la Ingeniería, aunque, no se ha - Analizadas estas observaciones, se detecta
agotado el análisis y la discusión, pues, nunca una regularidad, por ejemplos, el CBR
se consideran todos los elementos, ni la varia de manera inversa con el índice de
totalidad de las variables que intervienen en un grupo, o, de manera inversa con el índice
fenómeno dado, así como tampoco los suelos de plasticidad.
de todas las regiones del mundo. - Si lo expresado en al párrafo precedente se
Consecuentemente, los aportes, como el repite para una gran cantidad de
presente, van a seguir enriqueciendo el acervo observaciones, en diferentes tipos de
del conocimiento científico de la Mecánica de suelos y bajo una amplia variedad de
Suelos y sus implicaciones en la Ingeniería condiciones, entonces, será posible
Geotécnica. construir uno o varios modelos teóricos
(matemáticos) que expresen tales
4. METODOLOGÍA circunstancias.
Eso es precisamente lo que se hace en este
En el ámbito metodológico, es pertinente trabajo, es decir, a partir de los hechos
diferenciar entre metodologías científicas singulares (mediciones experimentales) se
generales y métodos particulares que dependen genera uno o más modelos teóricos, es decir, se
de la naturaleza del problema. En este trabajo, usa el método inductivo.
se han usado ambos tipos de metodologías, las Según el mismo autor (Chalmers, 1999), la
cuales son descritas, sucintamente en los exigencia para que esté justificada una
párrafos siguientes. inferencia inductiva como la señalada, desde
los hechos observables (mediciones
4.1 Metodologías científicas generales experimentales en este caso), hasta las leyes
(modelos matemáticos), debe satisfacer las
La expresión “la ciencia se deriva de los siguientes condiciones:
hechos” puede ser interpretada en sentido de (a) El número de observaciones, que
que el conocimiento científico se genera constituyen la base de la generalización,
estableciendo primero los hechos y edificando debe ser grande.
luego una teoría, un modelo, etc., que se ajusta (b) Las observaciones se deben repetir en una
a ellos (Chalmers, 1999). amplia variedad de condiciones.
En el caso de la Mecánica de Suelos e (c) Ningún resultado observacional aceptado
Ingeniería Geotécnica, es posible partir de debe entrar en contradicción con la ley
mediciones experimentales de variables general (modelo) derivada.
inherentes a ciertos fenómenos. Estas En general, estas tres condiciones son
mediciones son hechos singulares, que, se cumplidas al formular los modelos teóricos
supone, darán lugar a la construcción de una (matemáticos).
teoría o modelo general (modelos

6
 POTOSÍ

Pero, también, en este trabajo, se usa el método Donde β̂0 y β̂1 son los estimadores de los
deductivo (hipotético-deductivo), pues, se parámetros 0 y 1; Yi son los valores
formula una hipótesis respecto de la forma medidos (experimentales); Y ̂i son los
matemática de los modelos y, una vez valores estimados (con el modelo).
comprobada o contrastada la hipótesis, se Resolviendo las ecuaciones (3), se hallan
acepta el modelo como válido. Pero, no solo los parámetros buscados.
eso, sino que, a partir de los modelos generales, 3. Especificar la distribución de probabilidad
se deducen consecuencias, se hacen del componente aleatorio ε.
predicciones, interpretaciones y se realizan 4. Evaluar la “idoneidad” del modelo. Esto es
explicaciones. calcular el coeficiente de correlación, R, y
En suma, en este trabajo de investigación se el coeficiente de determinación, R2.
usan, como metodologías científicas generales, 5. Realizar un análisis de los residuales o
el método inductivo y el método deductivo. ̂i − Yi ).
errores (Y
6. Validar prácticamente el modelo, es decir,
4.2 Métodos particulares comparar los valores observados o
medidos (Y), con los valores estimados
Así como hay metodologías científicas con el modelo (Ŷ).
generales, que usan todas las ciencias, hay 7. Usar el modelo para hacer predicciones de
otros métodos específicos que conducen a la Y.
solución de problemas particulares. En este El método presentado se generaliza para
caso, vale la pena formular la siguiente modelos de regresión múltiple.
pregunta: ¿Cómo construir los modelos Sea Y una variable (dependiente) que
teóricos (matemáticos) que expresen las leyes corresponde a una propiedad mecánica del
de variación del CBR como función de otras suelo. Sean X1, X2. . . . , Xp, las variables
variables, a partir de las observaciones o (independientes) correspondientes a otras
mediciones experimentales? propiedades mecánicas del suelo. Además,
Una respuesta a esta pregunta tiene su base en supóngase que se tienen n mediciones de cada
el uso de los métodos estadísticos, en particular una de las citadas variables. Se plantea un
el método de los cuadrados mínimos, que, por modelo lineal general de la forma:
su importancia, se sintetiza a continuación.
Supóngase que, Y es una variable p
Y = β0 + ∑i=1 βi X i + ε (4)
(dependiente) que corresponde a una propiedad
mecánica de los suelos, por ejemplo, el CBR.
Del mismo modo, sea X una variable Donde:
(independiente) que corresponde a otra  = componente de error aleatorio;
propiedad del suelo, por ejemplo, el índice de p = número de variables independiente;
plasticidad. p
Si se utiliza el método de los cuadrados Y = β0 + ∑i=1 βi X i (5)
mínimos, el problema puede ser tratado
mediante la siguiente estrategia (Mendenhall, es el componente determinista.
W. y Sincich, T., 1997 y Benítez, 2004): o,  1, . . . . , p son los parámetros del
1. Proponer un modelo hipotético, por modelo. Se formulan n ecuaciones:
ejemplo, Y = βo + β1X + ε; siendo ε el
componente aleatorio. Y1 = o + 1X11 + 2X12 + . . . . . . . + pX1p
2. Encontrar los parámetros 0 y 1. Se hace: Y2 = o + 1X21 + 2X22 + . . . . . . . + pX2p
. . . . . . . . . . . . . . . .
SSE = ∑ni=1(Y
2
̂i − Yi ) = mínimo (2) Yn = o + 1Xn1 + 2Xn2+ . . . . . . . + pXnp
(6)
∂SSE ∂SSE Donde:
̂0 = 0; ̂1 =0 (3) Yi = i-ésima observación (medición) de Y;
∂β ∂β

7
 POTOSÍ

Xi,j = i-ésima observación de la j-ésima (Atterberg), propiedades granulométricas

variable independiente. (porcentajes más finos que las aberturas de
Si se usa notación matricial, el modelo lineal diferentes tamices), densidad máxima
general se puede expresar como: seca, humedad óptima, CBR y otras.
- Pero la observación, en el ámbito
Y=X + (7) científico, no solamente es una medición
mecánica, sino que se refiere, además, a
Siendo: una actividad mental, una reflexión, la cual
Y = matriz de los valores de la variable tiene que ver con los aspectos formulados
dependiente; en la introducción.
X = matriz de los valores de las variables
independientes; (2) Formulación del problema:
β = matriz de los parámetros. La segunda etapa, consiste en la formulación
La solución de cuadrados mínimos de esta del problema, que, en este caso, se sintetiza en
ecuación matricial es: términos de las preguntas ya formuladas en la
introducción.
̂ = (𝐗 ′ 𝐗)−𝟏 𝐗′𝐘
𝛃 (8)
(3) Las hipótesis:
Ecuación matricial, que permite la estimación Una vez formulado el problema, se plantea la
de los parámetros del modelo. hipótesis, como una respuesta anticipada
Un modelo de la forma: (modelo teórico) al problema formulado. Esto
es precisamente lo que se hace a continuación,
Y = o X11 X22 . . . . Xpp (9) de manera muy puntual y específica para el
caso investigado.
puede ser transformado, fácilmente, a la forma Sean:
lineal. - Capacidad portante del suelo (fenómeno
Con las particularidades de la regresión geotécnico), expresada en términos de su
múltiple, en general, la secuencia del método variable (dependiente) CBR = Y, al 100%
es análoga a la presentada para el modelo de de la densidad seca máxima;
regresión simple. - X = variable independiente general (que
puede ser Xo, X1, X2 o X3, que se llamarán
4.3 Síntesis de la estrategia metodológica variables independientes específicas);
- Índice de grupo del suelo, representada por
De manera didáctica, se mencionan a la variable (independiente) Xo;
continuación, las etapas seguidas en la - Plasticidad de un suelo (fenómeno
investigación, usando las metodologías geotécnico), representada por las variables
científicas generales y métodos particulares, independientes X1 (límite líquido) y X2
presentados en los párrafos que anteceden. (índice de plasticidad);
- Tamaño de los granos del suelo (propiedad
(1) Observación: física), expresado en términos de la
Se trata de una etapa fundamental en la variable independiente X3 (porcentaje más
investigación científica, sin la observación la fino que la abertura del tamiz N° 200, es
ciencia no avanza. El este caso, la observación decir, 0.075 mm).
tiene dos facetas, a saber:
- Las mediciones experimentales, sin duda (4) Formulación de las hipótesis o modelos
alguna forman parte de la observación. Se teórico-matemáticos, como una respuesta
han medido (en laboratorio) y determinado (solución) al problema formulado, en base
las variables principales, inherentes a las a las premisas establecidas en los párrafos
propiedades geotécnicas de los suelos precedentes.
bolivianos, a saber: límites de consistencia Primera hipótesis:

8
 POTOSÍ

La ley de variación del CBR = Y, en Si estas pruebas, especialmente las dos

función de la variable X (que puede ser Xo, primeras, son satisfactorias, las hipótesis están
X1, X2, o X3), viene dada por el modelo contrastadas o comprobadas y, en
matemático general: consecuencia, se ratifica su cualidad de ley.

A RESULTADOS Y DISCUSIÓN
Y= (10)
(X+B)
Bolivia tiene seis provincias fisiográficas o
Donde A y B son parámetros del modelo, a geológicas, a saber:
ser determinados (en cada caso) en base a (I) Cordillera Occidental;
la información experimental, por el (II) Altiplano;
método de los cuadrados mínimos (análisis (III) Cordillera Oriental;
de regresión simple o unidimensional). (IV) Sub-andino;
(V) Llanura Chaco-Beniana;
Segunda hipótesis: (VI) Escudo Brasileño.
La ley de variación del CBR = Y, en (Los números romanos identifican a las seis
función de más de una variable provincias fisiográficas o geológicas y serán
independiente, puede ser expresada mantenidos como tales en todo el trabajo).
mediante el modelo matemático general: Considerando estas provincias fisiográficas, se
han seleccionado los proyectos viales más
β1 β2 β3
Y = βo X1 . X 2 . X 3 (11) relevantes, que cuentan con información
experimental representativa, relativas a las
Donde βo, β1, β2 y β3, son parámetros del variables consideradas en la investigación.
modelo, a ser determinados (en cada caso) Entre estos proyectos, se tienen: Patacamaya-
en base a la información experimental, por Tambo Quemado, Oruro-Pisiga San Cristobal-
el método de los cuadrados mínimos Frontera, La Paz-Oruro, Oruro-Potosí, Potosí-
(análisis de regresión múltiple, Tarija, Potosí-Villazón, Cochabamba-Oruro,
multivariado o multidimensional). Padcaya-Bermejo, Tarija-Villamontes, Abapó-
Camiri, Santa Cruz-Puerto Suárez, Naureda-
(5) Contrastación o comprobación de las Extrema, Puerto Rico-Sena, Ixiamas-Chivé,
hipótesis: San Borja-Trinidad, Santa Cruz-Puerto Suárez,
En esta etapa, se ha procedido como sigue: Porvenir-El Chorro y Riberalta-Guayaramerín
- Cálculo del coeficiente de correlación, R, (Ministerio de Obras Públicas-ABC, 1970-
o, determinación, R2, que expresa la 2015).
“fortaleza” de la relación matemática entre
la variable dependiente y las variables CBR (Y) en función de una variable (análisis
independientes. Este es un primer de regresión simple o unidimensional):
indicador de “idoneidad” de los modelos.
- Prueba de valores observados (medidos Una vez analizada, seleccionada y
experimentalmente) contra valores sistematizada la información experimental, se
calculados (usando el modelo). El modelo han constituido las series estadísticas de
se dice “idóneo”, cuando los primeros son trabajo. En lugar de presentar estas series de
iguales a los segundos, o, las pequeñas datos, por razones de extensión (ya que son 18
diferencias, que siempre se presentan, se series), se han calculado los indicadores
atribuyen a errores de muestreo. estadísticos más relevantes.
- Análisis de residuales, que deben tener una Estos indicadores, para los modelos simples Y
fluctuación aleatoria, tomando como = f(Xo), se presentan en el cuadro N° 1.
referencia el cero.
- Distribución de residuales, que debe ser
aproximadamente normal).

9
 POTOSÍ

CUADRO N° 1 CUADRO N° 3
INDICADORES ESTADÍSTICOS INDICADORES ESTADÍSTICOS
MODELOS SIMPLES: Y = f(Xo) OTROS MODELOS SIMPLES
Var. M̅ S Min. Max. Var. M̅ S Min. Max.
(I) Cordillera Occidental (n = 18) (I) Cordillera Occidental (n = 30)
Y 12.97 18.00 1.60 74.40 Y 21.15 19.83 3.00 83.00
Xo 8.5 5.34 0.00 17.00 X3 42.95 23.06 9.20 94.30
(II) Altiplano (n = 29) (II) Altiplano (n = 28)
Y 3.76 3.02 1.30 15.00 Y 4.16 5.85 1.3 30.8
Xo 15.39 10.45 0.00 37.00 X2 18.59 7.93 1.7 34.5
(III) Cordillera Oriental (n = 19) (III) Cordillera Oriental (n = 22)
Y 8.57 14.65 2.20 66.80 Y 16.64 28.71 2.2 108.3
Xo 7.16 4.65 0.00 16.00 X3 51.52 21.00 9.5 79.9
(IV) Sub-andino (21) (IV) Sub-andino (n = 23)
Y 9.27 10.80 2.20 50.00 Y 14.62 21.00 2.20 80.00
Xo 9.29 6.80 0.00 27.00 X3 56.67 25.25 9.60 96.00
(V) Llanura Chaco-Beniana (n = 33) (V) Llanura Chaco-Beniana (n = 37)
Y 7.27 12.86 1.20 77.20 Y 12.46 24.82 2.8 129.00
Xo 16.58 9.64 0.00 35.00 X2 16.47 8.32 2.3 33.80
(VI) Escudo Brasileño (n = 42) (VI) Escudo Brasileño (n = 36)
Y 3.36 5.71 1.2 37 Y 6.70 12.18 1.30 57.00
Xo 22.38 13.22 0 52 X2 25.81 11.54 3.90 48.80
(Fuente: Elaboración propia) (Fuente: Elaboración propia)

Los resultados de los modelos matemáticos Y, los resultados de los modelos matemáticos
más idóneos (con un 95% de confiabilidad), más idóneos (con un 95% de confiabilidad),
correspondientes, se reflejan en el cuadro N° 2. que corresponden, están en el cuadro N° 4.

CUADRO N° 2 CUADRO N° 4
RESULTADOS PRIMERA HIPÓTESIS RESULTADOS PRIMERA HIPÓTESIS
Y = f(Xo) OTROS MODELOS SIMPLES
P.F. V.I. A B R P.F. V.I. A B R
(I) Xo 64.27 0.85 0.988 (I) X3 508.82 - 3.19 0.993
(II) Xo 42.7 3.01 0.976 (II) X2 49.63 - 0.10 0.989
(III) Xo 26.25 0.39 0.998 (III) X3 422.29 - 5.79 0.987
(IV) Xo 48.42 0.97 0.989 (IV) X3 448.64 - 3.80 0.983
(V) Xo 61.35 0.80 0.998 (V) X2 99.51 - 1.55 0.987
(VI) Xo 29.96 0.71 0.989 (VI) X2 69.11 - 2.71 0.998
(Fuente: Elaboración propia) (Fuente: Elaboración propia)

Análogamente, otros modelos simples (Y en De los cuadros N° 2 y N° 4, se obtienen los

función de X1, X2 o X3) se han probado, modelos matemáticos (ecuaciones) simples.
siempre en términos de la primera hipótesis. Por ejemplo, considerando los valores del
Los indicadores estadísticos se presentan en el cuadro N° 2 (tercera fila), se tiene para el
cuadro N° 3. Altiplano (II):

42.70
Y= (X0 +3.01)
(12)

10
 POTOSÍ

Análogamente, del cuadro N° 4, por ejemplo, han obtenido los resultados reflejados en el
para la Llanura Chaco-Beniana (V), se tiene: cuadro N° 6.

99.51 CUADRO N° 6
Y= (X2 −1.55)
(13)
RESULTADOS SEGUNDA HIPÓTESIS
βo β1 β2 β3
CBR (Y) en función de más de una variable (I) Cordillera Occidental (R2 = 0.973)
independiente (análisis regresión múltiple o
análisis multivariado): 8078 - 0.60 - 1.18
2
(II) Altiplano (R = 0.970)
Los indicadores estadísticos, para cada caso, se 3103 - 0.89 - 1.01
presentan en el cuadro N° 5. (III) 2
Cordillera Oriental (R = 0.940)
466251 - 1.17 - 1.37
CUADRO N° 5 2
(IV) Sub Andino (R = 0.887)
INDICADORES ESTADÍSTICOS
MODELOS MÚLTIPLES 3429 - 1.53
2
Var. M̅ S Min. Max. (V) Llanura Chaco-Beniana (R = 0.968)
(I) Cordillera Occidental (n = 50) 30887 - 0.48 - 1.66
Y 35.06 25.43 3.00 103.00 (VI) Escudo Brasileño (R2 = 0.937)
X1 29.10 8.64 16.80 58.40 3513 - 1.02 - 0.92
X3 31.18 21.53 8.00 94.30 (Fuente: Elaboración propia)
(II) Altiplano (n = 25)
Y 5.56 6.25 1.6 30.8 En los cuadros N° 1, 3 y 5:
X2 16.60 7.80 1.7 33.3 Var. = Variable;
X3 72.92 19.25 35.2 97.0 ̅ = Media aritmética;
M
(III) Cordillera Oriental (n = 23) S = Desviación estándar;
Y 19.65 30.33 2.20 108.30 Min. = Valor mínimo de la serie;
X1 33.92 11.33 19.60 55.70 Max. = Valor máximo de la serie;
X3 51.21 20.47 9.50 79.90 n = tamaño de la muestra (serie de trabajo).
(IV) Sub-andino (n = 26) En los cuadros N° 2, 4 y 6:
Y 15.92 20.72 2.20 80.00 R = Coeficiente de correlación;
X3 55.15 24.20 9.60 96.00 R2 = Coeficiente de determinación.
(V) Llanura Chaco-Beniana (n = 33)
Y 13.15 28.07 2.80 129.00 Las casillas en blanco del cuadro N° 6,
X2 17.40 8.58 2.30 33.80 significan que la variable en cuestión no
X3 82.92 21.38 20.90 99.40 interviene en el modelo. Por ejemplo, el
(VI) Escudo Brasileño (n = 39) modelo (ecuación) correspondiente la
Y 4.17 7.86 1.20 48.00 Cordillera Occidental (I), es:
X2 28.37 9.20 4.20 48.00 8078
X3 70.56 19.26 25.2 98.00 Y= (X0.60
(14)
1 X1.18
3 )
(Fuente: Elaboración propia)
Análogamente, el modelo (ecuación) del
En este caso, si bien se tienen 3 variables
Escudo Brasileño es:
independientes (X1, X2 y X3), cuyos
significados se han establecido en 4.3, los 3513
modelos matemáticos más idóneos no siempre Y= (X1.02
(15)
2 X0.92
3 )
incluyen todas estas variables. Así, con un 95%
de confiabilidad, para la segunda hipótesis, se De la misma manera, se tienen los modelos
para cada una de las provincias fisiográficas.

11
 POTOSÍ

En suma, se han obtenido 18 modelos fortaleza de la relación entre ambas

matemáticos, 3 para cada una de las seis variables (99.80%).
provincias fisiográficas de Bolivia (Cuadros - La prueba de valores observados
N° 2, N° 4 y N° 6). (medidos) contra valores calculados (con
Con fines didácticos, se ejemplifica, a el modelo), para ser satisfactoria, al ser
continuación, la secuencia general de las fases graficados, deben aproximarse a una recta
más relevantes, seguidas para la obtención de cuya pendiente sea igual a uno y cuya
cada uno de los modelos: ordenada al origen sea cero. Estos pares de
- Si se considera la serie de datos de la valores (Observados, Calculados), se
Provincia Fisiográfica VI, una vez representan en el gráfico N° 2, lo que
sometida al análisis de correlación y indica que la prueba es satisfactoria.
regresión la función Y = f(X2), es decir, - Prueba de residuales. Idealmente, los
CBR en función del índice de plasticidad, residuales o errores, (Y ̂i − Yi ), al ser
se obtiene, como modelo más idóneo (entre graficados, deben mostrar realmente
muchos), aquel cuyos parámetros se fluctuaciones aleatorias alrededor de un
encuentran en la última fila del cuadro N° valor cero, tal y como se muestra en el
4. Entonces, el modelo (a un 95% de gráfico N° 3.
confiabilidad) es: - La distribución de los residuales debe ser
aproximadamente normal. Para realizar
69.11 esta prueba, hay varias opciones. En este
Y= (X2 −2.71)
(16)
caso se ha usado la prueba de Kolmogorov-
Smirnov, que arrojó la siguiente
Tanto la curva teórica (modelo) como los conclusión: No se puede rechazar la idea
pares de puntos correspondientes a la serie de que la variable “residuales” provenga
experimental se presentan en el Gráfico de una distribución normal, con un 95% de
N° 1. confiabilidad.

Gráfico N° 1: Valores medidos y Gráfico N° 2: Prueba de valores

modelo teórico Y = f(X2) - Escudo medidos vs valores calculados con el
Brasileño modelo Y = f(X2) - Escudo
70
Brasileño
60 70
60
50
50
Y (Medidos)

40
40
Y

30
30
20
20
10
10
0 0
0 20 40 60 0.00 20.00 40.00 60.00 80.00
X2 Y (Calculados)

(Fuente: Elaboración propia) (Fuente: Elaboración propia)

- El coeficiente de correlación (Cuadro N°

4) es alto, lo que implica una medida de la

12
 POTOSÍ

Gráfico N° 3: Residuales versus lim Y = pequeño

variable independiente X2 → grande

2.500
lim Y = grande
2.000 X2 →pequeño

1.500
Lo cual es coherente con la realidad empírica,
1.000
Residual

pues, tanto el modelo como los datos

0.500 experimentales usados, expresan que en esta
0.000 provincia fisiográfica (VI), el CBR (Y) varía de
-0.500 0 10 20 30 40 manera inversa con el índice de plasticidad
-1.000 (X2), tal y como lo muestra, también, el gráfico
-1.500
N° 1.
X2 También, puede verse que:

(Fuente: Elaboración propia) lim Y = ∞

X2 →2.71

En consecuencia, se contrasta o comprueba la Pero, el valor mínimo de X2, considerado en la

hipótesis formulada. serie estadística de trabajo (mediciones
De manera análoga, se ha procedido en cada experimentales) es igual a 3.90 (última fila,
uno de los 18 modelos, que responden a las dos cuarta columna del Cuadro N° 3). Por tanto, el
hipótesis formuladas. valor X2 = 2.71, cae fuera del rango de datos
Consecuentemente, las hipótesis se convierten considerado en la formulación del modelo. El
en leyes. mínimo valor que debe usarse es X2 = 3.90.
Para este valor:
Análisis matemático y empírico
lim Y = 58.08
Un breve análisis y discusión de cualquiera de X2 →3.90
los modelos obtenidos, usando elementos del
cálculo diferencial, permite realizar algunas Resultado que es coherente con la realidad
consideraciones, teóricas y empíricas. Para empírica, casi igual al valor máximo del CBR
esto, como ejemplo, se toma el modelo de la experimental considerado, que es igual a 57
ecuación (16). (penúltima fila, última columna del Cuadro N°
Es fácil establecer (teóricamente) que: 3).
Similarmente:
lim Y = 0
X2 → ∞
lim Y = 1.50
X2 →48.80
lim Y = ∞
X2 →0
Siendo 48.80, el máximo valor experimental de
X2 (índice de plasticidad). El valor de 1.50
Es decir, el eje X2 es una asíntota. Igualmente (para Y), es cercano al empírico (1.30),
la variable Y, se vuelve asintótica en X2 = 2.71 reflejado en el mismo cuadro.
(asíntota vertical).
Albert Einstein dijo: “cuando no se refieren a ¿Qué pasa con las pendientes (derivadas) en la
la realidad, las matemáticas son exactas, en curva que corresponde al modelo matemático
tanto que, si se refieren a la realidad no son (Gráfico N° 1)?
exactas”. Consecuentemente, los límites Las respuestas se presentan en los párrafos
expresados anteriormente, si bien son teóricos, siguientes.
sirven para interpretar una realidad empírica, a La derivada de Y, respecto de X2 es:
saber:

13
 POTOSÍ

dY 69.11 (3) Las pruebas, tanto estadísticas, como

= − prácticas, han conducido a la contrastación
dX 2 (X 2 − 2.71)2
o comprobación de las dos hipótesis
Es claro que, para X2 = 2.71, la pendiente formuladas, mismas que, en consecuencia,
(derivada) se hace infinita (negativa), lo cual se constituyen en leyes de variación de la
no es práctico, pero basta reemplazar X2 = 3.90 capacidad portante de los suelos (CBR = y)
(mínimo valor experimental de X2) para en función de las otras variables
establecer que, para valores pequeños de X2 consideradas (Xo, X1, X2 y X3), tanto para
(índice de plasticidad), la pendiente se hace los modelos simples (primera hipótesis),
grande (en valor absoluto), igual a – 48.80, en como para los modelos múltiples (segunda
este caso. hipótesis), para cada una de las seis
Similarmente, cuando X2 se hace tan grande provincias fisiográficas o geológicas de
como se quiera, la pendiente (derivada) tiende Bolivia.
a cero. Si se reemplaza X2 = 48.80 (máximo (4) Los modelos matemáticos formulados
valor experimental), se aprecia que la derivada expresan una realidad empírica que puede
(pendiente) se hace muy pequeña, igual a – resumirse en los siguientes términos:
0.033, es decir, prácticamente horizontal, lo - Los resultados de los modelos,
que ratifica la tendencia asintótica con el eje reflejados en el cuadro N° 4 indican
horizontal. que en las provincias fisiográficas
Similares análisis matemáticos pueden hacerse donde se tiene una presencia
para todos los otros modelos. mayoritaria de suelos granulares (I, III
y IV), el CBR depende de una
5. CONCLUSIONES Y propiedad granulométrica (X3).
RECOMENDACIONES : - En tanto que, cuando la realidad
empírica (datos experimentales)
5.1 Conclusiones reflejan una mayor presencia de suelos
finos (zonas II, V y VI), los resultados
(1) Al haber determinado las leyes de del cuadro N° 4 indican que el CBR
variación de la capacidad portante de los depende de una propiedad plástica
suelos, en función de otras variables de (X2).
más simple medición experimental, - Sin embargo, en el caso de los modelos
expresadas en términos de modelos múltiples, solo en el sub-andino el
matemáticos (ecuaciones), mismos que CBR es una función de X3 (propiedad
permiten predecir (estimar) el CBR sin granulométrica), ya que, en las otras
recurrir a ensayos de laboratorio, se han provincias, es una función de una
logrado los dos primeros objetivos propiedad plástica (X1 o X2) y de una
diseñados para la investigación. propiedad granulométrica (X3).
(2) Habiendo comparado los valores del CBR (5) La matemática permite realizar análisis
(Y), observados (laboratorio) con los interesantes. En este caso, conduce a
valores estimados (modelos), siendo las visualizar los valores extremos (máximos
diferencias pequeñas (atribuibles a errores y mínimos), tanto de la variable
de muestreo), se concluye que, a tiempo de dependiente, como de las independientes.
dimensionar estructuras de pavimentos, las
diferencias se pierden en el redondeo, 5.2 Recomendaciones
habitual en el ingeniero, cuando se definen
los espesores de las diferentes capas del (a) Usar los modelos formulados solo para los
paquete estructural. Consecuentemente, se valores (tanto de la variable dependiente
ha alcanzado el tercer objetivo propuesto, como de las variables independientes),
estableciendo, además, la utilidad práctica comprendidos en los rangos de las series
en la aplicación de los modelos. estadísticas utilizadas, cuyos valores

14
 POTOSÍ

máximos y mínimos se precisan en los Benítez, A. (2012). Propiedades mecánicas de

cuadros N° 1, N° 3 y N° 5, pues, la los suelos, correlaciones, modelos
extrapolación fuera de estos rangos, podría matemáticos e implicaciones en la
arrojar resultados alejados a la realidad ingeniería vial. XVI Congreso Argentino
empírica. de Vialidad, 22-26 de octubre, Córdoba.
(b) En todos los modelos formulados y en las Benítez, A (2004). Generalización de las
series estadísticas usadas, se han relaciones entre los procesos del ciclo
considerado suelos que tienen un índice de hidrológico para la cuenca internacional
plasticidad diferente de cero. Se del río de La Plata: “Enfoque Sistémico y
recomienda, en futuras investigaciones, Análisis Multivariado”. Tesis de doctorado,
realizar un trabajo similar, con suelos con Universidad de Sevilla, España.
índice de plasticidad igual a cero. Bunge, M. (2011). La investigación científica.
Siglo Veintiuno Editores, México.
6. AGRADECIMIENTOS : Carter, M. y Bentley, S.P. (1991). Correlations
of soil properties. Pentech Press, London.
El autor agradece al Dr. Adrián Rodríguez- Chalmers, A.F. (1999). ¿Qué es esa cosa
Marek, (Departamento de Ingeniería Civil y llamada ciencia? Siglo Veintiuno de España
Ambiental del Virginia Tech, USA), quien ha Editores, Madrid.
revisado una primera versión de este trabajo de Gutierrez, R. (1996). Introducción al método
investigación, realizando sugerencias y científico. Esfinge, México.
observaciones importantes. Magnan, J.P. (1982). Les méthodes statistiques
et probabilistas en mécanique des sols.
7. REFERENCIAS : Presses de l’ecole nationale des Ponts et
Chaussées, París.
Agarwal, K.B. y Ghanekar, K.D. (1970). Mendenhall, W. y Sincich, T. (1997).
Prediction of CBR from plasticity Probabilidad y Estadística para Ingeniería y
characteristics of soils. Proceedings of 2nd Ciencias. Prentice Hall. México.
South-East Asian Conference on Soil Ministerio de Obras Públicas-ABC (1970-
Engineering, pp. 571-576, Singapore. 2015). Varios estudios de carreteras a nivel
Ameratunga, J., Sivakugan, N. y Das, B. de diseño final o TESA. La Paz.
(2016). Correlation of soil and rock Patel, S.R. and Desai, M.D. (2010). CBR
properties in geotechnical engineering. Predicted by Index Properties for Alluvial
Sringer, India. Soils of South Gujarat. Dec. 16-18, 79-82,
ASTM D1883-07e2 (2007). Standard Test Proceedings of Indian Geotechnical
Method for CBR (California Bearing Ratio) Conference–2010, India.
of Laboratory-Compacted Soils. American Terzaghi, K. y R. Peck (1967). Soil mechanics
Society for Testing of Materials, in engineering practice. John Wiley,
Pennsylvania, USA. Lodon.
Benítez, A. (2014). Dialéctica de la Ingeniería Vinod, P. and Reena, C. (2008). Prediction of
Estructural. Tesis de Doctorado en CBR Value of Lateritic Soils Using Liquid
Ingeniería Civil UNSXX-Convenio Andrés Limit and Gradation Characteristics Data,
Bello, La Paz. Highway Research Journal, IRC, 1 (1): 89-
98.

i [Ingeniero Civil, M.Sc., Ph.D., Dr.Sc.],

15