El ensayo CBR en detalle

1.- Introduccin
Que es un ensayo CBR? Intentar explicarlo sin meterme en detalles y con un
lenguaje que todos entiendan.
El ensayo CBR sirve para medir la resistencia de un terreno de cara a utilizarlo
en una carretera, es decir, si ponemos ese terreno debajo del asfalto, y lo
apisonamos bien, queremos saber si tras pasar muchos camiones terminarn
salindole baches a la carretera o no.
El procedimiento comienza por medir que significa "lo apisonamos bien", y eso
se hace con el ensayo PROCTOR, que consiste en medir la densidad del
material tras humedecerlo y aplastarlo: se humedece con 3 o 4 cantidades
diferentes de agua, se compacta, y se mide su densidad. Ahora se dibuja una
grfica humedad-densidad uniendo los puntos obtenidos y se busca
"visualmente" el mximo de la curva. Ese mximo corresponde a una cierta
humedad (la ptima) y corresponde a un cierta densidad (la mxima).
Ahora se trata de saber, para ese densidad "mxima" (o un 98% de ella si no
somos muy optimistas respecto de cuanto conseguiremos compactar -apisonar-
el terreno en la prctica), cuanto aguantar el terreno.
Para ello se toman tres cilindros rellenos de este material, y se compactan con
un martillo especial, dndoles martillazos, pero a unos moldes ms y otros
menos, de forma que el que ms se compacte consiga, aproximadamente la
densidad mxima proctor.
Ahora se "mide" el indice CBR de cada uno: Se va poniendo carga sobre el
cilindro hasta que comprimimos el terreno 0.25 mm (0.1" en la norma ASTM), y
comparamos la carga obtenida con la que aguantara un terreno "ideal", que se
fij en su da en 1000 PSI (libras por pulgadas) usando un "terreno ideal" para
carreteras de la zona de California (de aqu lo de C.B.R. o "California Bearing
Ratio", algo as como "porcentaje de aguante Californiano").
Esto nos da un porcentaje: si nuestro terreno compactado con X golpes nos
aguanta 600 PSI (hemos necesitado aplicarle esa presin para que se
hundiese 0.25 mm) y el terreno ideal aguanta 1000 PSI, nuestra muestra de
terreno tiene un ndice CBR del 60%.
Finalmente, tomamos los 3 ndices CBR obtenidos para los diferentes moldes
(con diferentes compactaciones), medimos sus densidades, y trazamos una
grfica densidad-ndice CBR. Como sabemos por el PROCTOR cual es la
densidad que realmente alcanzaremos en obra, solo resta mirar, para esa
"densidad objetivo" qu indice CBR nos da la grfica.
Ese es el indice CBR que tendr nuestro material una vez colocado en la
carretera, humedecido con la humedad ptima, y compactado (apisonado)
correctamente.
Finalmente, en obra solo nos quedara tener cuidado con una cosa: Tras
humedecer y apisonar una parte de terreno y antes de asfaltar, debemos medir
la humedad y densidad real para asegurarnos de que hemos alcanzado el
punto ptimo "PROCTOR" o necesitamos aadir ms agua y volver a apisonar.
Este es el ensayo de "densidad in-situ" que suele hacerse por mtodos
radioactivos para agilizar.
2.- Variaciones segn la norma del CBR
El ensayo CBR incluido en HCLab permite aplicar varias normativas (ASTM,
UNE, NLT y las normas Mexicana, Chilena, Colombiana y Peruana). Aunque
los clculos que se realizan son bsicamente los mismos (la norma ASTM es la
base de todas las dems normas), existen pequeas diferencias que pasamos
a comentar.
Las normas actualmente disponibles para el ensayo CBR son estas:
NLT 111/87 (Norma Espaola antigua).
UNE 103.502:1995 (Norma Europea actual).
ASTM D1883-07 (Norma U.S.A. "original").
M-MMP-1-11-08 (Norma Mexicana).
INV E-148-07 (Norma Colombiana).
NCh1852.Of81 (Norma Chilena, para una sola probeta).
MTC E 132 - 2000N (Norma Peruana, para una sola probeta).
Veamos las pequeas diferencias entre ellas punto por punto.
2.1) Nmero de golpes por capa:
Es una gua de cuantos golpes por capa (los moldes se van rellenando en 2 o
ms capas, segn nos indique la norma para cada tipo de uso previsto) habra
que darle a un molde, con una maza estandarizada y desde una altura fija, para
conseguir una "energia de compactacin" equivalente a la utilizada en el
ensayo proctor con el que obtubimos la densidad mxima, es decir, si le damos
60 golpes "estandar" el molde debera conseguir, aproximadamente, la
densidad mxima indicada en el proctor.
UNE y NLT: 60 golpes.
ASTM, Mexico, Colombia, Per y Chile: 56 golpes.
Nota: Chile no indica golpes pero te remite a normas especficas de
compactacin (NCh1534/1 y 2) donde, para el molde en cuestin, indica 56
golpes por capa.
En la norma ASTM y, en general, en todas las normas americanas, el nmero
de golpes no es totalmente fijo y se permite usar ms o menos golpes si vemos
que no conseguimos alcanzar la densidad mxima proctor, que al fin y al cabo
es el objetivo final de darle todos estos golpes (en los clculos siguientes solo
se necesita saber la densidad alcanzada, no cuantos golpes le tuvimos que
dar).
2.2) Lecturas de Cargas:
Una vez tenemos los 3 moldes correctamente "golpeados", hemos de
presionarlos hasta que cedan 0.25 mm (0.1") y luego seguir hasta que cedan
0.5 mm (0.2"), ya que necesitamos las cargas aplicadas en ambos casos (ya
dije que la introduccin era algo "simplificada").
Respecto de qu unidades se han de usar para anotar estas cargas, cada
norma usa las que ve bien, algunas usan cargas y otras presiones, y la mayora
mezcla 2 o 3 posibles unidades. La verdad es que ninguna hace mencin sobre
que solo valga usar unas unidades, al final se van a convertir a un % por lo que
ninguna norma se preocupa mucho.
Pero por tener una gua, estas son las unidades en que se tratan estos temas:
NLT: Usa kN o kp para las lecturas (unidades de fuerza) pero luego usan PSI,
MPa o kp/cm
2
(unidades de presin) para la comparacin con los valores
ideales.
UNE: Se usa N, kN o kp para todo (no convierte las cargas a presiones
nunca).
ASTM: No indican unidades para las lecturas (usan "loads", cargas), y usan
PSI o MPa para las presiones ideales.
Mexico: usa kN para todo (no convierten las cargas a presiones nunca).
Colombia: Solo hablan de kPa y PSI (es extrao ver los kPa, se suelen usar
MPa).
Chile: No indica unidades de las lecturas, pero pide pasarlas a presiones en
MPa.
Per: Habln de kg e incluyen la conversin a libras.
Nota: El ensayo CBR de HCLab pregunta siempre cargas en kp, pero se puede
eligir una "constante de anillo" adecuada, de forma que si usamos una
constante de 0.102 podemos introducir las cargas en kN.
2.3) Convertir cargas en presiones:
Solo tenemos que dividir la carga por el dimetro del pistn usado para aplicar
la carga, que aunque en la prctica son los mismos, segn las normas su
dimetro varia ligeramente (en teora, la superficie del pistn ha de ser de 3
pulgadas cuadradas):
ASTM, Mexico, Colombia, Per: 49.63 mm.
UNE y NLT: 49.60 mm.
Chile: 49.5 mm. (indica que se utilice un rea de 19.35 cm
2
, lo que
corresponde a 49.635 mm)
2.4) Penetraciones "objetivo":
Ya tenemos unas cargas aplicadas con un pistn de un cierto dimetro, es
decir, presiones, pero necesitamos anotar la presin que corresponde con unas
ciertas deformaciones (o penetraciones) del molde especiales.
ASTM, Mexico, Colombia, Per y Chile: Usan 0.1" (2.54 mm) y 0.2" (5.08
mm).
UNE y NLT: Usan valores redondeados de 2.5 y 5 mm.
Nota 1: En la introduccin solo habl de 2.5 mm. por simplificar la explicacin.
Nota 2: Esta pequea diferencia hace que el valor "ideal" sea ligeramente
distinto aunque se mida a partir del mismo suelo.
2.5) Presiones ideales:
Ya tenemos dos cargas (o presiones) por cada molde, correspondientes a las
penetraciones "objetivo": ahora podemos comparar nuestros moldes con el
"caso ideal", ya que este "caso ideal" viene dado por las dos cargas -o
presiones- necesarias para alcanzar las dos penetraciones "objetivo" en un
molde con ese material "ideal".
En cada norma nos dan el caso ideal usando parejas de cargas y/o presiones,
una para 2.5 mm (0.1") y la otra para 5 mm (0.2").
NLT 2.5 mm: 1000 PSI, 6.90 Mpa. 70.31 kp/cm
2
.
NLT 5.0 mm: 1500 PSI, 10.35 Mpa 105.46 kp/cm
2
.
UNE 2.5 mm: 13.2 kN 1346.9 kp.
UNE 5.0 mm: 20 kN 2040.8 kp (1 kN = 1000/9.8 kp = 102 kp).
ASTM 0.1": 1000 PSI 6.90 MPa.
ASTM 0.2": 1500 PSI 10.30 Mpa.
Mxico 0.1": 13.34 kN.
Mxico 0.2": 20.01 kN. (iguales que ASTM si se usa g=10).
Colombia 0.1": 1000 PSI 6.90 MPa.
Colombia 0.2": 1500 PSI 10.35 Mpa (aparece tambien como 10.3 en otro
punto).
Chile 0.1": 6.90 MPa aprox. 70 kp/cm
2
.
Chile 0.2": 10.30 Mpa aprox. 105 kp/cm
2
.
Per 0.1": 6.90 MN/m
2
( MPa) 70.31 kgf/cm
2
1000 lb/pulgada
2
(PSI).
Per 0.2": 10.35 MN/m
2
( MPa) 105.46 kgf/cm
2
1500 lb/pulgada
2
(PSI).
Nota: La UNE utiliza un valor de g = 9.8 m/s
2
en sus conversin kN a kp,
mientras que la norma Mexicana parece traducir los valores de la ASTM a kN
usando g = 10 m/s
2
.
2.6) ndice CBR de cada molde:
De las dos penetraciones "objetivo" hemos obtenido dos cargas por molde,
pero necesitamos un nico ndice CBR para el molde. En teora habra que
convertir cargas a presiones, compararlas con la dos presiones ideales,
obteniendo dos porcentajes, y finalmente quedarnos con uno de ellos. En la
pctica esto se resume en multiplicar las dos cargas por dos coeficientes que
nos den directamente esos porcentajes, ms un criterio para elegir uno de los
dos porcentajes que hemos obtenido.
Nota: Los coeficientes por los que multiplicar las cargas salen de comparar con
los valores ideales del punto 2.5) aunque hay que convertirlos antes a kp. En
las normas similares a la ASTM (Mxico, Colombia y Chile) he utilizado g = 10
m/s
2
tal como veo que ellos usan, y de esta forma los coeficientes salen iguales
a los de la norma ASTM. Ms abajo se incluyen ejemplos de estos clculos.
UNE : Max(CBR
2.50
* 0.0742, CBR
5.00
* 0.049)
NLT : Max(CBR
2.54
* 0.0735, CBR
5.08
* 0.049)
ASTM : Max(CBR
2.54
* 0.0749, CBR
5.08
* 0.050)
Mexico: Min(CBR
2.54
* 0.0749, CBR
5.08
* 0.050)
Colombia: Max(CBR
2.54
* 0.0749, CBR
5.08
* 0.050)
Chile: Max(CBR
2.54
* 0.0749, CBR
5.08
* 0.050)
Per: Max(CBR
2.54
* 0.0749, CBR
5.08
* 0.050)
Notas:
La norma UNE nos da ya los coeficiente a utilizar precalculados.
La norma NLT usa coeficientes distintos a todas las otras normas porque el
radio de la placa es 0.3 mm inferior a la ASTM y las penetraciones objetivo
son diferentes a las de la UNE, aparte no indica si se debe usar g=10 o 9.8
por lo que utilic el valor "correcto" de 9.8, mientras que ASTM parece preferir
g = 10.
Mxico da las cargas ideales en kN, para convertirlas he usado g = 10
m/s
2
de forma que los coeficientes coincidan con los de la ASTM.
Solo la normativa Mexicana utiliza el menor de los dos valores.
La norma ASTM y la peruana "MTC E 132" comentan que, si el mayor ndice
es el correspondiente a 0.2" no nos lo creamos a la primera y repitamos el
molde antes de dar por bueno el valor para 0.2", pero al final hemos de tomar
el mayor de los dos.
Ejemplo de conversin a kp para la norma ASTM:
Presin mxima a 2.54 mm. Segn ASTM es 6.9 MPa = 6.9 N/mm
2
.
Como dimetro es 49.63 mm, obtenemos superficie 1934.54 mm
2
.
Presin mxima * superficie = carga mxima = 13348.35 N.
Si g = 9.8 carga mxima = 13348.35/9.8 = 1362.076 kp.
Si g = 9.8 coeficiente = 100/1362.076 = 0.0734
Si g = 10 carga mxima = 13348.35/10 = 1334.835 kp.
Si g = 10 coeficiente = 100/1334.835 = 0.0749
Ejemplo de conversin a kp para la norma Mexicana:
Carga mxima 13.34 kN (similar al coef. ASTM para g=10).
13340/10 = 1334.000 kp (1334.835 en ASTM).
100/1334 = 0.0749 (con ms precisin 0.07496251874).
Nota: En la norma Mexicana hubiese sido algo ms exacto usar un coeficiente
redondeado a 0.0750, pero no mereca la pena sobre todo sabiendo que se
basa en la ASTM y que la diferencia se debe a redondeos en la conversin
entre unidades.
2.7) Grfica densidades / ndices CBR:
Una vez se tienen los 3 moldes, disponemos de sus 3 densidades y sus tres
nidices CBR, y tenemos que representarlos en una grfica, que se utiliza para
determinar el nidice CBR que correspondera con la densidad "objetivo"
(densidad proctor * compactacin objectivo).
Pues resulta que algunas normas piden explicitamente situar las densidades en
el eje X de esta grfica (ASTM, Colombia y Mexico), mientras la UNE piede
hacerlo en el eje Y.
La NLT no lo menciona, as que por "proximidad" a la UNE se usa el eje Y para
las densidades, mientras que la norma Chilena tampoco lo menciona, y por
proximidad a la ASTM se utiliza el eje X.
Nota: En el resto de grficas (densidad/humedad, densidad/absorcin y
densidad/hinchamiento) la densidad ocupa el mismo eje que en la grfica
anterior de densidad/CBR.
2.8) Tipo de CBR "por edades":
Existe una variante del CBR algo especial donde, a diferencia de los casos
"normales" comentados, no se hacen tres moldes con diferentes densidades,
sino que se hacen tres moldes para ser ensayados a tres edades distintas, y no
se busca el nidice CBR para una cierta "densidad objetivo", sino para una
cierta "edad objetivo" de 7 das.
La idea es que se ensaya un suelo con un cierto aadido de Cal, y esta mezcla,
una vez humedecida, tarda un tiempo en "fraguar", con lo que se necesita
ensayar el suelo a los 7 das de ser mezclado.
Este sistema aparece comentado, al menos, el el PG-3 espaol, pero
desconozco si existe norma alguna al respecto (el PG-3 es un "pliego de
condiciones" espaol para carreteras, no una normativa).
En estos casos (admitidos por el ensayo CBR de HCLab mediante un tercer
tipo de proctor: "Normal", "Modificado" y "Por edades") lo que se pretende es
ensayar a 1, 3 y 7 das de edad para tener una "curva de fraguado" y con ella
poder, posteriormente, estimar el ndice CBR que se obtendr a 7 das
conociendo solo el de 1 y 3 das (es similar a romper probetas de hormign a 7
y 28 das para obtener un "coeficiente de recuperacin" que nos prediga la
resistencia a 28 a partir de la de 7).
Otra posibilidad es ensayar a ms de 7 das para saber si la mezcla sigue
mejorando con el tiempo o se estanca en el ndice CBR obtenido a 7 das.
En cualquier caso, todo lo comentado aqu servira para el CBR por edades con
solo cambiar "densidades" por "das" al calcular el ndice CBR final y utilizar
como "densidad objetivo" el 7.
3.- Clculos realizados en el CBR
Aparte de las pequeas diferencias comentadas entre normas, el procedimiento
es siempre el mismo, as que vamos a detallar, paso a paso, como se haran
los clculos internos de un ensayo CBR, con especial nfasis en los "puntos
grises" de la norma y como automatizar su resolucin. Para el ejemplo
aplicaremos los coeficientes de la norma UNE.


3.1) Clculo de densidades, humedades y absorciones:
Es lo primero a calcular para cada uno de los tres moldes, y lo ms importante
de cara al ndice CBR final ser la densidad de cada molde.


3.1.1) Empezamos con la Humedad antes de la inmersin:
Agua_antes = Tara+Suelo+Agua - Tara+Suelo
Suelo_antes = Tara+Suelo - Tara
Humedad_antes = 100 * Agua_antes/Suelo_antes


3.1.2) Calculamos la densidad seca (antes de la inmersin):
Suelo+Agua = Molde+Suelo+Agua - Molde
Suelo_Despues = 100 * Suelo+Agua / (100 + Humedad_antes)
Agua = Suelo+Agua - Suelo_Despues
Densidad_Seca = Suelo_Despues / Volumen_molde


3.1.3) Calculamos la energa de compactacin:
Se calcula como un porcentaje sabiendo que 60 golpes corresponden al 100%
de la energa de compactacin (o 56 golpes segn la norma utilizada).
Compactacin = Numero_Golpes / 60 * 100 (redondeado al entero ms prximo).


3.1.4) Volvemos con la humedad despues de la inmersin:
Nota 1: Tara+Suelo despus de la inmersin es un clculo solo si se puso el
100% del suelo en la estufa. Si solo se seca en estufa una porcin, entonces
ser un dato que el operario debe introducir. Aqu suponemos que se introduce
el 100% en la estufa, por lo que se puede calcular. Ademas, la tara despus de
la inmersin coincide con el peso del molde.
Nota 2: Se obtiene sacando el suelo del molde tras la inmersin, y usando los
25 mm. de la parte superior. La normativa Chilena comenta que se puede medir
tambin la humedad promedio de toda la probeta de forma opcional (el ensayo
utilizado en HCLab no contempla una segunda humedad tras la inmersin,
aunque se puede incluir como anotacin en el ensayo).
Tara+Suelo_Despues = Molde + Suelo_Despues
Agua_Despues = Tara+Suelo+Agua_Despues - Tara+Suelo_Despues
Suelo_Despues = Tara+Suelo_Despues - Tara_Despues
Humedad_Despues = 100 * Agua_Despues / Suelo_Despues
Nota: Las normativas ASTM, Mexicana y Chilena piden incluir en los resultados
la "humedad tras inmersin", que es igual a la Humedad_Antes +
Agua_Absorbida, por lo que algunas normas piden el "agua absorbida" y otras
la "densidad tras inmersin".


3.1.5) Lo siguiente es calcular el agua absorbida:
Agua_Absorbida = Humedad_Despues - Humedad_Antes


3.2) Hinchamientos individuales:
El hinchamiento se mide dejando a remojo el molde totalmente relleno de
material, y midiendo tras 96 horas cuanto ha aumentado en altura.
Este incremento en altura se compara con la altura inicial para obtener un % de
hinchamiento. Esta altura inicial corresponde con la altura del molde (177.8
mm) menos la altura del collar de expansin (50.8 mm), es decir, 127 mm de
altura inicial.
Nota: Se pueden llegar a medir hinchamientos negativos.
Hinchamiento = (Lectura_Final - Lectura_Inicial ) * 100 / 127


3.3) Grfica de penetraciones y cargas empleadas:

3.3.1) Cargas aplicadas:
Las penetraciones en milmetros estn fijadas por las distintas normas. Las
cargas introducidas se suponen en kp y se multiplican, si existe, por una
constante de anillo (podemos usar una constante de 0.102 si preferimos
introducir kN, por ejemplo).
Fuerza_Real = Fuerza * Constante_Anillo


3.3.2) Calculo de los desplazamientos del origen "segn norma":
Este es el punto ms conflictivo de todas las normas (todas copian de la
ASTM): matemticamente hablando es un disparate lo que se propone, pero
usando la intuicin se entiende cual es el propsito: Si las primeras cargas solo
sirvieron para "asentar" las piedrecitas de la parte superior del molde, esas
medidas sern "falsas" y las penetraciones que medimos a continuacin todas
"arrastrarn" este error.

A nivel de ensayo se propone dejar una carga inicial suave durante un rato para
que la superficie se asiente y podamos comenzar el ensayo sin falsas
mediciones, pero al final, estas medidas falsas afloran, y si miramos la grfica
"carga-penetracin" se suele ver que la curva no siempre "va a morir" al origen
de la grfica, as que al final, tenemos que calcular ese "desplazamiento del
origen" y restarlo en todas las medidas de penetraciones posteriores si
queremos un resultado fiable.
Todas las normas existentes proponen hacerlo detectando "el punto de
inflexin" de la grfica y usarlo para calcular el corte de la grfica con el origen
(a imagen de la norma ASTM), pero en la practica no es posible usar ese
mtodo de forma automatizada, ya que se basa en la intuicin del tcnico en
varios aspectos.
Por un lado, lo normal es tener varios puntos de inflexin, pero la norma solo
contempla la existencia de uno, por lo que el tcnico debe elegir uno
intuitivamente. Por otro lado, y aunque solo existiera un punto de inflexin, este
puede no ser apropiado para usar su tangente tal como describe la norma para
calcular el desplazamiento del origen. En general, solo en casos idealizados
este mtodo puede ser aplicado.
Se incluye a continuacin un ejemplo grfico de aplicacin de la norma al pie
de la letra a un molde de un CBR real que da lugar a un clculo totalmente
errneo del desplazamiento.

Por un lado, se ven varios puntos de inflexin (marcados con un crculo
amarillo), y en cada uno de ellos se ha trazado la tangente (en azul) para
comprobar si corta donde intuitivamente esta el origen de la grfica.
Si elegimos el ms favorable para la aplicacin del mtodo de la tangente (en
este ejemplo parece ser el segundo de los puntos) y trazamos la tangente,
vemos claramente (lnea roja) que no corta al eje donde intuitivamente debera,
por lo que da lugar a un desplazamiento de la grfica del orden de 0.5 mm,
cuando a simple vista debera de ser del orden de 0.25 mm.
Como puede apreciarse en el grfico, no solo el resultado obtenido usando el
segundo punto es errneo, lo peor es que si eligisemos automticamente el
punto de inflexin ms prximo al origen (el primero de ellos), obtendramos un
desplazamiento negativo, que desplazara a la grfica a la derecha en lugar de
a la izquierda, dando lugar a ndices CBR del orden de la mitad del correcto, o
incluso, como se observa en el tercero de los puntos, no podramos obtener un
ndice CBR tras desplazar la grfica ya que esta quedara desplazada ms alla
de los 2.5 mm e incluso de los 5 mm.


3.3.3) Clculo de desplazamientos del origen "mejorado":
El mtodo utilizado por HCLab est diseado para dar unos resultados
similares al "mtodo intuitivo" indicado en la norma, pero sin ser tan sensible a
pequeas variaciones en puntos individuales. De esta forma, el mtodo puede
ser aplicado automticamente sin peligro de caer en un caso como el
comentado arriba.
El mtodo consiste en calcular la parbola que mejor se ajusta por mnimos
cuadrados a los puntos obtenidos, y luego obtener el punto de corte de dicha
parbola con el eje Y=0 ms cercano a X=0 (siempre existe un segundo punto
de corte "lejano").
A continuacin mostramos el mismo caso del punto anterior pero esta vez
usando "el mtodo de la parbola".

Como se v, este mtodo nos da el desplazamiento correcto, del orden de 0.25
mm, y el hecho de que algunos puntos aislados no sigan la curva terica no
distorsiona en absoluto el desplazamiento obtenido, como ocurre en el mtodo
de la tangente.
Nota: Fsica y matemticamente, sera ms adecuado hacer un ajuste
logartmico en lugar de usar parbolas, ya que el asentamiento claramente es
un proceso localmente logartmico (cuanto ms se asienta, ms dificil es
aumentar dicho asientamiento). La realidad es que la diferencia entre utilizar un
tipo u otro de funcin es del orden de la milsima de milmetro, por lo que en la
prctica son mtodos equivalentes.


3.3.4) Estabilidad del mtodo de clculo de los desplazamientos:
Para ilustrar esta estabilidad del mtodo frente a pequeas variaciones en las
lecturas individuales, mostramos el caso anterior con dos de los puntos
ligeramente modificados, de forma que la grfica sea menos abrupta y ms
similar a un caso ideal.

Como se puede ver, un ligero cambio en los datos nos lleva a unos clculo de
desplazamiento totalmente diferente al caso original, de nuevo con tangentes
imposibles que cortan al eje X en valores negativos. Es decir, el mtodo original
es muy inestable frente a perturbaciones (mala cosa para un mtodo
matemtico cualquiera).
Ahora vemos este mismo caso modificado con el mtodo de la parbola.

Como se observa, el resultado obtenido vuelve a ser de 0.25 mm sin importar
esos mnimos cambios introducidos. Si se prueba haciendo pequeos cambios
en algunos puntos, siempre se obtiene un desplazamiento casi idntico: el
mtodo es muy estable frente a perturbaciones en los datos.


3.4) Indices CBR individuales:
Una vez se han obtenido los desplazamientos necesarios para aplicarlos a la
grfica, se calculan dos ndices CBR por cada molde, uno con el valor de carga
correspondiente a 2.5 mm de penetracin (0.1" en las normas americanas) y un
segundo en 0.5 mm (0.2" en normas americanas), quedandonos con el mayor
de los dos (el menor en la norma Mexicana).
Nota: Es equivalente desplazar la grfica -0.25 mm (a la izquierda) que sumar
0.25 a la X en la que queremos calcular la altura de la grfica, es decir,
usaremos X = 2.5+0.25 sobre la grfica original en lugar de usar X = 2.5 sobre
la grfica desplazada, ya que el resultado es el mismo.
Para el primero de los dos candidatos a ndice CBR se toma X = 2.5 mm +
desplazamiento, es decir, 2.75 mm en el ejemplo comentado, y se busca el
valor de Y (altura) en la grfica para esta X, obtenindose un valor en kp.
X = 2.5 mm + Desplazamiento = 2.75 mm.
kp = Altura de la grfica en X = 2.75 mm.
CBR
2.5
= kp * 0,0742
Nota: El coeficiente 0,0742 varia ligeramente entre normas como se vio
anteriormente.
El segundo ndice CBR se calcula de forma similar pero para 5 mm (0.2") y con
otro coeficiente:
X = 5.0 mm + Desplazamiento = 5.25 mm.
kp = Altura de la parbola en X = 5.25 mm.
CBR
5.0
= kp * 0,049
Finalmente, se toma el mayor de los dos ndices obtenidos como ndice CBR
del molde.
CBR = Max(CBR
2.5
, CBR
5.0
)
Nota: La normativa Mexicana, por otro lado igual a la ASTM, indica que se ha
de elegir el MENOR de los dos valores, mientras que la ASTM indica que si el
segundo fuese mayor, se repetira el molde por si fuese un error, y luego se
toma el mayor de los dos.


3.5) Indice CBR final:
Llegados a este punto, tenemos tres moldes con tres ndices CBR, pero solo
existe un suelo y por tanto solo necesitamos un nidice CBR "final" que mezcle
los tres obtenidos de forma correcta.


3.5.1) La densidad "objetivo":
Es la densidad prevista en la obra real y para la cual se quiere obtener el ndice
CBR, se calcula en funcin de la densidad ptima proctor y el porcentaje de
compactacin previsto, que suele ser del 95, 98 o 100% segn lo "optimista"
que sea el ingeniero.
Dens_CBR = Dens_Proctor * Compactacion_prevista / 100


3.5.2) El ndice CBR correspondiente a la densidad objetivo:
Una vez obtenidos los tres ndices CBR y sabiendo las densidades secas de
cada uno de los moldes, se ha de calcular el ndice CBR que correspondera a
la densidad prevista que llamamos Dens_CBR.
Este clculo se realiza mediante interpolacin por splines cbicos, tomando los
puntos (Densidad_Seca, Indice_CBR) correspondientes a los tres moldes, y
hallando a continuacin el valor de Y para una X = Dens_CBR, que sera el
valor final del ndice CBR para la compactacin elegida.

Nota: Este mismo sistema de interpolacin es utilizado para calcular el
hinchamiento, la humedad y la absorcin final, tomando como puntos los
formados por (Densidad_Seca, Hinchamiento) , (Densidad_Seca, Humedad) y
(Densidad_Seca, Absorcin) respectivamente, tal como se ve en las tres
grficas pequeas.
Nota: Suele darse el ndice CBR para diferentes % de compactacin, de forma
que si se nos pide alcanzar un cierto ndice CBR en obra, podamos elegir la
mejor compactacin a posteriori: igual pensbamos llegar el 98% pero con el
95% vemos que se obtiene un CBR suficientemente alto y nos ahorramos un
par de das compactando el suelo.


3.6) Los splines cbicos utilizados al interpolar:
Las curvas utilizadas para interpolar un conjunto de n puntos (P
1
, P
2
, ... P
n
),
donde cada punto P
i
= (P
i
x, P
i
y) con P
i
y = ndice CBR del molde i-esimo y P
i
x =
su densidad seca en nuestro caso, se conocen como "Splines cbicos de
Hermite" y, en especial, la variante por "diferencias finitas"
(http://en.wikipedia.org/wiki/Cubic_Hermite_spline).
El nombre es largo e impone: "Spline cbico de Hermite por diferencias finitas",
pero realmente es el sistema matemticamente ms simple para unir los puntos
"sin que se note", y lo de "sin que se note" se resumen en los puntos A, B y C
de ms abajo.
Estas curvas se calculan como la unin de n-1 tramos de curva, donde cada
tramo une dos puntos consecutivos de la lista.
Cada uno de estos n-1 tramos de curva se define en forma paramtrica como
F
i
(t) = (F
i
x(t), F
i
y(t)), con t entre 0 y 1, donde ambas componentes son
polinomios de tercer grado y han de cumplir estas propiedades (definicin de
spline de Hermite):
A) Cada tramo de curva F
i
(t) une P
i
con P
i+1
, es decir:
F
i
(0) = P
i

F
i
(1) = P
i+1

B) Las derivadas coinciden en los puntos de unin, es decir:
F
i
'(1) = F'
i+1
(0)
C) Las derivadas segundas son nulas en los puntos de unin, es decir:
F
i
''(0) = F
i
''(1) = 0
Nota: En el punto B no se concreta cuanto ha de ser la derivada en cada punto
de unin, por eso se utiliza la variante por "diferencias finitas" comentada, que
para entendernos nos indica qu esta derivada (es un vector) ha de ser la mitad
de la lnea que une el punto anterior con el posterior (ver ejemplo ms
adelante).
Por tanto, para el clculo de cada tramo se necesitarn los puntos anterior y
posterior, es decir, necesitamos en total cuatro puntos por tramo: P
i-1
, P
i
, P
i+1
y
P
i+2
.
Es necesario pues duplicar los puntos iniciales y finales de forma que se aada
a la lista de puntos un P
0
y un P
n+1
, para poder calcular los tramos primero y
ltimo. Esta duplicacin -que es un mero artifcio matemtico- se consigue
definiendo P
0
=P
1
y P
n+1
=P
n
, es decir, duplicamos el primer y ltimo punto,
obteniendo as la lista de puntos (P
0
, P
1
, ... P
n
, P
n+1
).
Nota: Normalmente -en la propia wikipedia- se usa otro artifcio matemtico en
lugar de duplicar los puntos inicial y final, que es usar una frmula distinta en
estos dos casos. Personalmente me parece mucho ms sencillo duplicar
puntos que jugar con las frmulas (pero es lo mismo).
Una vez definidos los trminos, pasamos a desglosar el clculo de estas
funciones F
i
(t) = (F
i
x(t), F
i
y(t)) con la ayuda de este esquema.


3.6.1) Las derivadas V
1
y V
2
:
Primero calculamos los vectores V
1
y V
2
que unen el punto anterior y posterior
en cada uno de los dos puntos a unir (P
i
y P
i+1
) y dividimos las longitudes a la
mitad. Obtenemos as dos vectores V
1
y V
2
, que sern usados como las
derivadas en los puntos inicial y final del tramo.
Si realizamos los clculos coordenada a coordenada, obtendramos estas
formulas:
V
1
x = 0.5*(P
i+1
x - P
i-1
x)
V
1
y = 0.5*(P
i+1
y - P
i-1
y)
V
2
x = 0.5*(P
i+2
x - P
i
x)
V
2
y = 0.5*(P
i+2
y - P
i
y)


3.6.2) Funciones auxiliares:
Para simplificar la formula final de F
i
(t), se definen cuatro funciones polinmicas
cubicas auxiliares:
C
1
(t) = 2*t
3
-3*t
2
+1
C
2
(t) = t
3
-2*t
2
+t
C
3
(t) = -2*t
3
+3*t
2

C
4
(t) = t
3
-t
2



3.6.3) Finalmente, las formulas de F
x
(t) y F
y
(t) seran:
F
x
(t) = C
1
(t)*P
i
x + C
2
(t)*V
1
x + C
3
(t)*P
i+1
x + C
4
(t)*V
2
x
F
y
(t) = C
1
(t)*P
i
y + C
2
(t)*V
1
y + C
3
(t)*P
i+1
y + C
4
(t)*V
2
y
Nota 1: Al ser estas funciones una combinacin lineal de polinomios de grado
tres, sern ellas mismas tambin polinomios de grado tres "a tramos", de aqu
el nombre de cubicas, porque tanto la coordenada X como la Y ambas son
polinomios cbicos.
Nota 2: Las ecuaciones del punto 3.6.2 son consecuencia directa de las
propiedades que exigimos a la curva, pero su deduccin no entra en los
objetivos de este artculo.


4.- Incertidumbre y margen de error
Segn las normativas genricas, como la UNE-EN ISO IEC 17.025=2005, se
debe incluir en las actas de resultados una estimacin del margen de error -la
incertidumbre- del resultado, siempre que esto sea posible. Por ir al grano: en
este ensayo eso no es posible.
La propia ASTM dice al respecto dos cosas: la dispersin de los resultados
ante repeticiones no se ha determinado por dificil y costosa, aunque comenta
que unos test particulares lo sitan en el orden del 9.2 y el 6.9% segn el
mtodo de compactacin utilizado. Estos test fueron sobre siete repeticiones,
por lo que son solo estimaciones.
Pero el problema real viene con el error inherente al mtodo de clculo, no a
imprecisiones de ciertas medidas: no se puede obtener una cota mxima del
error que se produce al seguir este mtodo. La norma ASTM dice al respecto
"There is no accepted reference value for this test method, therefore, bias
cannot be determined".
La explicacin es sencilla: Toda norma que en algn punto pida que "unamos
los puntos obtenidos con una curva ms o menos suave" y la usemos como
parte del proceso de obtencin del resultado, como ocurre tambin con en
PROCTOR, llevan asociado un error en ese paso que es terica y
prcticamente ilimitado.
Dependiendo de como se interprete el prrafo, de si usamos interpolacin o lo
hacemos a mano, segn que tipo de curva usemos al interpolar, o segn el
pulso y la gracia que se tenga si se unen a mano, obtendremos una curva u
otra, y por tanto un resultado diferente. Se pueden encontrar casos donde esta
eleccin produce resultados tan dispares que el error que se produce al elegir
una y no otra es tan grande como queramos, por lo que no hay un "error
mximo" que podamos utilizar.
Si tiene dudas a este respecto, revise las conclusiones del punto 3.3.2 y
siguiente, aqu comprobbamos con un caso real como la forma de "unir los
puntos con una curva suave" puede suponer obtener un CBR del valor que ms
nos interese.
En el caso del CBR, comparado con otros donde la curva se utiliza para
encontrar un valor mximo (el proctor por ejemplo), es an ms problemtico,
ya que se dan como ciertas dos premisas falsas: que la curva que obtengamos
tendr un solo punto de inflexin, y que la tangente a ese punto de inflexin
apuntar directamente al verdadero origen de coordenadas de la grfica.
Ambas suposiciones son falsas y producen resultados disparatados.