Columnas de Grava
Columnas de Grava
Columnas de Grava
RESUMEN Se presenta una revisión de los métodos de cálculo de la capacidad de carga y asientos de columnas de grava, y
sus puntos más inciertos o deficientes, así como recientes avances que permiten una solución más ajustada. El análisis se refie-
re al caso de columnas bajo terraplenes o áreas extensas cargadas. En primer lugar, se presenta el problema de la predicción
de la estabilidad de las columnas en suelos blandos, para lo que la aplicación de los distintos métodos conduce a resultados
dispares. El tema central de la conferencia es el análisis del reparto de cargas entre terreno y columnas. Este reparto, varia-
ble a lo largo del proceso de consolidación, está gobernado por la rigidez relativa de columnas y terreno y tiene múltiples con-
secuencias: la evaluación de la estabilidad de las columnas en suelos blandos, la magnitud de los asientos producidos, la velo-
cidad de consolidación radial hacia las columnas como elementos drenantes y la capacidad portante media del terreno
reforzado. En el texto se presentan avances recientes de carácter teórico, contrastados con experimentación en laboratorio y
observación de casos reales, que suponen una importante mejora de los métodos de cálculo tradicionales.
Palabras clave: Mejora del terreno, Columnas de grava, Solución analítica, Experimentación en laboratorio, Consolidación.
Keywords: Ground improvement, Stone columns, Analytical solution, Consolidation, Laboratory experiments.
1 Texto de la conferencia pronunciada por el autor en la Sociedad Española de Mecánica del Suelo e Ingeniería Geotécnica, el 25 de marzo de 2010.
pa
sz sr
σrc σrs d rl
rc
σzc u
σzs
L
COLUMNA
SUELO kz
rl
ES υS kr rc
Ec υc φc ψc
d
rc
FIGURA 1. Esquema de
rl
reparto de tensiones.
Definición de la celda unidad.
2. REPARTO DE TENSIONES ENTRE SUELO Y COLUMNAS • Estabilidad de las columnas: Debido al fuerte aumento
de la carga vertical sobre las columnas, éstas pueden
El fenómeno básico que gobierna el comportamiento de un te- llegar al estado límite último activo. En casos bajo te-
rreno reforzado por columnas de grava es el reparto de tensio- rraplenes el suelo resulta cargado en dirección radial
nes verticales entre ambos. Si se aplica una presión vertical por la columna en expansión, y también verticalmente
media de valor pa sobre la superficie (Figura 1), se reparte de por la carga del terraplén, con lo que se encuentra muy
forma desigual entre las columnas (σzc) y el suelo (σzs). Ello confinado y por tanto sujeto a esfuerzos desviadores pe-
suele cuantificarse mediante el llamado “factor de concentra- queños. En cambio, el problema es crítico cuando las co-
ción de tensiones” (SCF): lumnas se emplean como inclusiones bajo cimentacio-
σ zc nes aisladas, pues al no estar cargada la superficie, la
SCF = (1)
σ zs expansión radial de la columna puede llegar a provocar
también la plastificación del terreno. Los métodos
En general, las columnas son más rígidas que el suelo, por usuales para cuantificar esta posibilidad son muy sim-
lo que la tensión sobre las columnas es mayor que la presión ples (Mitchell, 1981; Van Impe y DeBeer, 1983), y sue-
media aplicada sobre la superficie, y, en paralelo, la tensión len ignorar el reparto de tensiones suelo-columna, o
sobre el suelo es menor. Ello afecta a: mezclar las tensiones iniciales y finales.
• Capacidad portante: Al ser las columnas un material Existe ya una experiencia notable sobre el comporta-
altamente friccional, su resistencia al corte es función miento de suelos reforzados con columnas de grava. Las
directa de la tensión vertical que actúa sobre ellas (σzc), tensiones en suelo y columnas no son fáciles de medir, so-
por lo que la capacidad portante del conjunto depende bre todo en estas últimas, por lo que en la mayor parte de
del reparto de tensiones. casos publicados, sólo se dispone de medidas de asientos.
• Asientos del terreno: Los asientos que sufre el suelo de- La Figura 2 muestra una recopilación de este tipo de me-
penden obviamente de la tensión actuante sobre él (σzs), didas. La dispersión es menor de lo que cabría esperar, ya
por lo que se produce una reducción del asiento respeto al que se trata de casos sin discriminar por el tipo de suelo o
caso de no refuerzo. Suele denominarse “factor de mejora” de columnas, sólo por su área de sustitución. Sobre ellos,
(n) al cociente entre el asiento que se produciría si no se se presentan dos curvas: la contenida en la Guía de Ci-
instalaran columnas y el asiento con columnas. Su in- mentaciones (Ministerio de Fomento, 2002), que está defi-
verso se denomina “factor de reducción de asientos” (α). nida precisamente como una aproximación conservadora a
• Velocidad de asientos: La velocidad del proceso de con- las observaciones reales, y la de Priebe (1976, 1995), de
solidación depende del diámetro y espaciamiento de co- carácter semiempírico, que reproduce una tendencia me-
lumnas, al actuar éstas como drenes radiales. Pero, por dia.
otra parte, la tensión sobre el suelo varía con el tiempo, Las tensiones en suelo y columnas deben cumplir en todo
al ir variando su rigidez aparente, desde la correspon- caso la condición de equilibrio de fuerzas verticales que, si ar
diente a material incompresible en el momento de apli- es la fracción del área total que ocupan las columnas, se ex-
car la carga, hasta su rigidez efectiva drenada, al final presa como:
del proceso. Ello hace que la consolidación sea más rá-
arσ zc + (1 − ar )σ zs = pa (2)
pida que si la carga sobre el suelo fuese constante.
Factor de reducción
de asientos, α=1/n
1.0
Medidas reales
Priebe (φc = 40°)
Guía de cimentaciones (MF, 2002)
0.8 α=(1-ar)2
0.6
0.4
0.2
0
FIGURA 2. Asientos finales. 0 0.2 0.4 0.6 0.8 1.0
Resumen de medidas en casos
reales. Soluciones empíricas o Razón de sustitución, ar
semiempíricas.
Esta condición de equilibrio no basta para determinar el lumna, las condiciones de contorno y de simetría hacen
reparto de tensiones, sino que se requiere hacer intervenir la que la tensión vertical sea siempre uniforme. En el
deformabilidad relativa de ambos elementos. A continuación suelo, en cambio, en algunas situaciones puede haber
se presentan las soluciones existentes, en orden creciente de una variación de la tensión vertical con la distancia a la
complejidad. Los resultados se van analizando en función de columna. Cuando esto ocurre, las expresiones como (1)
su grado de aproximación a la experiencia existente. y (2) tienen sentido entre valores medios de las tensio-
En este estudio se supone terreno arcilloso, en el que nes. En esos casos, cuando es precisa una distinción, se
puede considerarse que la carga se aplica instantáneamente, denotan en lo que sigue dichos valores medios con una
con el suelo en condiciones no drenadas (las columnas sí dre- barra superior (σ–).
nan instantáneamente), analizándose los estados inicial (no • El flujo de agua durante la consolidación es puramente
drenado) y final (drenado), y el proceso de consolidación que radial, desacoplado del vertical mediante la ecuación de
transcurre entre ambos. Dada la rapidez de la consolidación Carrillo (1942).
provocada por la presencia de las columnas, esta situación
debe considerarse representativa sólo en casos de suelos cla- En las soluciones de tipo elástico, las situaciones inicial y
ramente arcillosos. El estudio no pretende abarcar toda la final pueden estudiarse independientemente, pues la solución
gama de condiciones posibles, sino que se centra en las si- final no depende del camino seguido. Sin embargo, con mate-
guientes: riales elastoplásticos es preciso el análisis completo del pro-
ceso de consolidación para obtener la solución para el estado
• Refuerzo de suelos blandos bajo terraplenes o áreas ex- final, como en su momento se comenta.
tensas cargadas, excluyendo el uso como inclusiones
bajo cimentaciones aisladas. 2.1. COLUMNA CONFINADA LATERALMENTE
• Columnas totalmente penetrantes hasta un sustrato rí- La hipótesis más sencilla es la de confinamiento lateral, tanto
gido, excluyendo las columnas flotantes. para la columna como para el suelo. Debe señalarse que la
• Se estudia el comportamiento de una celda cilíndrica celda cilíndrica unidad impone la condición de deformación
unidad (Figura 1), formada por una columna y el suelo lateral nula sólo en el contorno exterior de la célula, pero aquí
circundante, con simetría axial. Se trata de una confi- se supone que esta condición la cumplen además tanto la co-
guración representativa, en las condiciones implicadas lumna como el suelo circundante, deformándose ambos en
por las hipótesis anteriores. condiciones edométricas. Sus deformaciones verticales vienen
• La carga aplicada se supone que produce una distribu- dadas por:
ción uniforme de asientos de columna y terreno (igual σ zc σ zs
deformación vertical), suponiendo que no existe pene- ε zc = ; ε zs = (3)
Emc Ems
tración de las columnas en la base del terraplén.
Igualando ambas deformaciones se obtiene la condición de
• La hipótesis anterior supone por contra que las tensio- compatibilidad, que lleva a un factor de concentración de ten-
nes no son uniformes. En primer lugar, difieren entre la siones igual a la relación modular:
columna y el suelo, como se ha comentado (ecuación 1).
Por otra parte, tanto en uno como en otro elemento, σ zc Emc
= (4)
también pueden ser variar a lo largo del radio. En la co- σ zs Ems
Factor de reducción
de asientos, α=1/n
1.0
Medidas reales
Confinado
Elástico libre
0.8
0.6
0.4
Ec/Es=10
Ec/Es=40
0.2
Esta condición, junto con la de equilibrio (2) permite obte- Los resultados anteriores se refieren a la situación final, dre-
ner las dos tensiones, y el resto de variables se determina di- nada. Para la situación inicial, de carga sin drenaje, la condición
rectamente. La tensión final en el suelo resulta: de que el suelo no cambie de volumen, junto con la de deforma-
ción lateral nula implica que no puede sufrir ningún asiento, por
pa
σ zs, f = lo que la condición de compatibilidad lleva a que la carga sobre
Emc (5)
1 + ar − 1 la columna es, en ese instante, nula. La tensión sobre el suelo es
Ems entonces ligeramente superior a la presión media aplicada, pa, y
y el factor de mejora: el factor de concentración de tensiones inferior a la unidad.
La solución confinada presenta una inconsistencia, pues
1 p E no se cumple el equilibrio radial de tensiones en la interfaz
n= = a = 1 + ar mc − 1 (6)
α σ zs, f Ems columna/suelo: al suponerse que ambos elementos se defor-
man en condiciones edométricas, sus respectivas tensiones
horizontales vienen definidas por la tensión vertical en cada
uno multiplicada por el correspondiente coeficiente de empuje
40
Solución elastoplástica
en reposo, por lo que no tienen por qué ser iguales.
(Castro, 2008)
Ensayos, ar=0.06-0.11 (Cimentada, 2009)
En cualquier caso, el resultado se presenta en la Figura 3,
con líneas discontinuas de trazo largo, para dos valores de la re-
lación modular que cubren el rango habitual. Puede verse que el
Factor de concentración, SCF (σzc/σzs)
20
2.2. COLUMNA ELÁSTICA CON DEFORMACIÓN RADIAL LIBRE
Solución elástica La solución anterior puede mejorarse dejando libre la defor-
(Balaam-Booker, 1981)
mación radial de la columna y del suelo e igualando ambas.
Rango medidas reales
(Barksdale-Bachus, 1983)
De esta forma, se puede cumplir la condición de equilibrio en
10
el contacto entre ambos.
La distribución elástica de tensiones en la columna y en el
terreno se puede obtener directamente utilizando soluciones
conocidas:
0
• La columna es un cilindro sometido a unas tensiones
0 10 20 30 40 50 vertical en sus bases (σzc), y radial en su pared lateral
Relación modular, Ec/Es (σrc), uniformes y desconocidas (situación análoga a una
probeta en un ensayo triaxial convencional, que da
FIGURA 4. Factor de concentración de tensiones. Soluciones teóricas y como resultado tensiones y deformaciones uniformes en
casos reales. toda la columna).
• En cuanto al suelo que la rodea, se trata de un tubo de puede verse que este efecto mejora algo la comparación con
pared gruesa, sometido a unas tensiones vertical en sus las medidas reales, pero sigue sobreestimando la reducción de
bases (σzs), y radial en su pared lateral (σrs), también asientos. En la Figura 4 se puede ver el efecto en el factor de
desconocidas. La solución elástica de este caso se en- concentración de tensiones.
cuentra en cualquier texto básico de Elasticidad. Debe también observarse que en la situación inicial, sin
Para obtener las cuatro tensiones desconocidas anteriores drenaje, el asiento no es nulo (como ocurría en la solución
se dispone (en direcciones vertical y radial) de las dos ecuacio- confinada), sino que el suelo sufre un pequeño asiento, ex-
nes de compatibilidad de deformaciones y las dos de equilibrio pandiéndose al mismo tiempo horizontalmente (para conser-
de tensiones. Con ello se obtiene la solución completa de ten- var el volumen), lo que supone una pequeña contracción ra-
siones y deformaciones en la columna y en el suelo. dial en la columna.
El resultado fue presentado por Balaam y Booker (1981)
de forma completa, tanto para la situación inicial de carga sin 2.3. COLUMNA ELASTOPLÁSTICA
drenaje como para la final, drenada. En la situación final, la La comparación anterior sugiere que la hipótesis elástica, aun
tensión sobre el suelo es: considerando la situación realista de deformación radial libre,
pa conduce a una tensión excesiva sobre la columna. Para mejo-
σ zs, f =
Emc − 2λ c F (7) rar este aspecto, es preciso introducir la existencia de defor-
1 + ar − 1 maciones plásticas.
Ems + 2λ s Far /(1 − ar )
La posibilidad de plastificación afecta en principio sólo a la
siendo: columna, que resulta sometida a una elevada presión vertical y
(λ c − λ s )(1 − ar ) con un confinamiento lateral inferior al de reposo y que podría
F= llegar al límite activo. En cambio, el suelo se encuentra en con-
2[ar (λ s − λ c + Gs − Gc ) + λ c + Gc + Gs ]
diciones de confinamiento lateral importante (la columna su-
y el factor de mejora: pone una inclusión más rígida que el terreno al que sustituye).
E − 2(λ c − λ s ) F Dado que en análisis elastoplástico el estado final depende
1
n=
= 1 + ar mc − 1 (8) del camino seguido, es preciso considerar por separado los ca-
α Ems sos de carga con y sin drenaje, que se hace a continuación de
Comparando estas expresiones con las (5) y (6), respectiva- forma cualitativa, analizando las trayectorias típicas de ten-
mente, se ve que, con respecto a la hipótesis edométrica, au- siones efectivas en la columna en un diagrama de tensiones
menta algo la tensión sobre el suelo y por tanto, los asientos. radial y vertical (Figura 5), a partir de un estado inicial de re-
En la Figura 3 se representan juntas ambas soluciones, y poso, k0 (punto I).
kac kac
Fe Fe
Fp,D
Fp,UC
YUC
k0c k0c
Tensión vertical, σzc
1.0 1.0
YD
U
kpc kpc
I I
FIGURA 5. Esquema de
trayectorias de tensiones en la
columna. Casos con y sin
Tensión radial, σrc Tensión radial, σrc
drenaje.
2.3.1. Carga con drenaje gura 5-b), caracterizado por un aumento relativo mayor de la
En el caso de suelo permeable (Figura 5-a) todo el proceso de tensión radial, por lo que el estado se acerca a la isotropía, y la
carga tiene lugar con drenaje, y las tensiones en la columna columna se mantiene en estado elástico. A continuación, tiene
crecen monótonamente hasta el estado final. Para columna lugar el proceso de consolidación, lo que para columna elástica
elástica, ello lleva al punto Fe, definido en la Figura 5-a pero, llevaría al punto Fe, el mismo que para carga con drenaje. Sin
dependiendo de la magnitud de la carga aplicada y de las ca- embargo, si la resistencia de la columna es limitada, se alcanza
racterísticas del material, puede alcanzarse antes el estado lí- el estado plástico en un cierto punto YUC, durante el proceso de
mite activo en un punto tal como YD. A partir de este momento, consolidación. Al final de éste, se llega a un estado final Fp,UC,
la columna se deforma en régimen plástico, llegando al final en diferente del elástico (Fe) y del estado plástico final para carga
el punto Fp,D, diferente de Fe (la mayor deformabilidad de la co- con drenaje (Fp,D). Comparando ambos casos en la Figura 5
lumna en régimen plástico implica una reducción del factor de puede verse que es diferente el momento de plastificación (Y), y
concentración de tensiones respecto al caso elástico). la proporción relativa de las fases elástica y plástica del proceso.
Este caso ha sido resuelto analíticamente por Pulko y Majes Este caso es más complejo que el anterior, ya que obliga a
(2005). La solución del problema sigue siendo en principio posi- la integración de la ecuación de la consolidación en régimen
ble, pues en el planteamiento sólo cambian las relaciones consti- no lineal, debido a las deformaciones plásticas de la columna.
tutivas en la columna, donde las tensiones y deformaciones son Sin embargo, la solución no es imposible, pues el suelo sigue
uniformes. Aparecen dos funciones incógnitas más (las compo- siempre en régimen elástico, por las razones expuestas más
nentes plásticas de las deformaciones radial y vertical), y dos arriba. En el apartado que sigue se presenta el análisis.
ecuaciones adicionales, las correspondientes al estado plástico:
3. CONSOLIDACIÓN
Condición de plastificación (Rankine) : σ rc = kacσ zc
(9) Al igual que para el estado final, se presentan aquí las solu-
Condición de flujo plástico (dilatancia) : ε zcp = −2kψcε rcp ciones existentes, en orden de complejidad creciente, comen-
zando por la hipótesis de confinamiento lateral.
siendo:
1 − sen φ c 1 − sen ψ c 3.1. COLUMNA CONFINADA LATERALMENTE
kac = ; kψc = ;
1 + sen φ c 1 + sen ψ c La hipótesis de confinamiento lateral implica que la columna
φc= ángulo de rozamiento interno; está lejos del estado límite activo. Por lo tanto, se trata siem-
ψc = ángulo de dilatancia. pre de un proceso elástico, para el suelo y la columna.
Las deformaciones elásticas de la columna siguen regidas 3.1.1. Solución básica para drenes (Barron-Hansbo):
por la ley de Hooke, pero el desarrollo algebraico se simplifica
notablemente si se suponen despreciables frente a las compo- Como referencia, se dispone de la solución clásica de Barron-
nentes plásticas. Hansbo, formulada para la consolidación radial alrededor de
drenes prefabricados. La ecuación (conservación de masa de
2.3.2. Carga sin drenaje fluido y ley de Darcy), es en este caso:
Para terreno de permeabilidad reducida, la carga puede supo- kr ∂2u 1 ∂u ∂εvol,s
+ = − (10)
nerse instantánea, dando lugar a un estado U sin drenaje (Fi- γ w ∂r2 r ∂r ∂t
σz
Tensiones σz, u, σ´z
En todo lo que sigue, se designa como u la sobrepresión in- Esta variación de la carga ha sido tenida en cuenta en el
tersticial (presión intersticial en exceso sobre la hidrostática análisis presentado por Han y Ye (2001). En el desarrollo ante-
inicial). rior, debe modificarse la ecuación (11), expresando la tensión
La deformación volumétrica del suelo, en condiciones con- efectiva en función de la presión intersticial media, en la forma:
finadas, es igual a la vertical, la cual a su vez depende de la σ ' zs (t) σ zs (t) − u(t)
tensión vertical efectiva en cada punto, a través del módulo εvol,s = ε zs == (17)
edométrico. La hipótesis de deformación vertical uniforme im- Ems Ems
plica por tanto que, aunque tanto la tensión total como la in- La tensión total media sobre el suelo (– σzs)para cualquier
tersticial varían a lo largo del radio, su diferencia debe ser tiempo t se puede obtener con las condiciones de equilibrio y
constante, e igual a la existente entre los valores medios de compatibilidad, como se indicó en el Apartado 2.1 para la si-
ambas magnitudes (Figura 6). Por otra parte, los drenes pre- tuación final, resultando:
fabricados tienen una rigidez despreciable, por lo que no so- E
portan ninguna carga, y la tensión media en el suelo es igual pa + mc ar u(t)
Ems
a la presión media aplicada, pa. Todo ello lleva a: σ zs (t) = (18)
E
σ ' zs σ zs − u σ zs − u pa − u 1 + ar mc − 1
εvol,s = ε zs = = = = (11) Ems
Ems Ems Ems Ems Introduciendo este valor en la ecuación (17) se obtiene la
Sustituyendo este valor en la ecuación (10): deformación volumétrica (εvol,s) que, sustituida en la ecuación
_ (10) lleva a la misma ecuación de consolidación de Barron-
∂2u 1 ∂u ∂ u
cvr 2 + = (12) Hansbo, pero con un coeficiente de consolidación equivalente:
∂r r ∂r ∂t
E a
siendo cvr el coeficiente de consolidación radial del suelo: cvrze = cvr 1 + mc r (19)
Ems 1 − ar
(cvr = kr Ems γ w ) . donde el superíndice ze se refiere a la influencia de la defor-
La solución de la ecuación (12) es: mación confinada de la columna elástica. El factor entre pa-
réntesis es siempre mayor que la unidad, lo que implica una
u r 1 r
2
ln − ar − 1 consolidación más rápida que en el caso de carga constante.
u= (13)
f (ar ) rc 2 rc El tratarse de una mera modificación de cvr implica que la
curva de consolidación (U-Tr) es la misma que en el caso de
siendo la sobrepresión intersticial media: carga constante (16), pero trasladada horizontalmente en el
8 Tr eje logarítmico de tiempos.
− (14)
u = pa e f ( ar )
3.2. COLUMNA ELÁSTICA CON DEFORMACIÓN RADIAL LIBRE
Tr es el factor de tiempo (Tr = cvr t 4rl2 ) , y: El método seguido en el caso de confinamiento lateral puede se-
guirse también para deformación radial libre de columna y
ln ar 3 − ar suelo. La ecuación de la consolidación radial sigue siendo la
f (ar ) = − − (15)
2(1 − ar ) 4 misma (10), pero ahora la deformación volumétrica no es igual
a la vertical, como era en las ecuaciones (11) y (17), sino que in-
El grado de consolidación se define en la forma usual: terviene también la deformación horizontal. Para calcularla se
8 Tr debe acudir a la solución elástica general. La solución se ob-
u −
U = 1− = e f ( ar ) (16) tiene como se presentó en el Apartado 2.2 para las situaciones
pa inicial y final, pero considerando que existe una cierta presión
intersticial por disipar. Para estos efectos, se puede sustituir la
presión intersticial en cada punto (u) por su valor medio (u–),
3.1.2. Aplicación directa a columnas de grava (reparto de
variable con el tiempo, aunque en este caso la justificación de
tensiones constante): esta aproximación no es tan directa como en el caso confinado.
La solución de Barron-Hansbo suele aplicarse en la práctica al La solución completa puede verse en Castro (2008). Todas
caso de columnas de grava directamente, sin más que suponer las deformaciones y tensiones resultan necesariamente linea-
que la tensión sobre el suelo no es la presión media aplicada, les con pa y u–, en la forma:
pa, sino la que resulta del reparto de tensiones suelo-columna,
suponiendo que es constante e igual a su valor final, σzs,f, dado
) (σ , ε) = Φ e u + Ψ e pa (20)
por la expresión (5). El resultado, si se expresa mediante el Donde las diferentes funciones Φ e y Ψ e para cada compo-
grado de consolidación adimensional, es obviamente el mismo nente pueden depender del radio r y de los parámetros geo-
anterior, tanto en presiones intersticiales como en asientos. métricos y elásticos, pero no del tiempo, que sólo interviene a
través de u–. Por tanto, la derivada temporal de la deforma-
3.1.3. Reparto de tensiones variable con el tiempo ción volumétrica del suelo que debe introducirse en la ecua-
El reparto de cargas entre columna y suelo varía considera- ción de la consolidación (10), es lineal con la derivada (∂u ∂t) .
blemente durante el proceso de consolidación. Inicialmente, al Ello da como resultado la misma ecuación (12), pero con un
no haber drenaje el suelo no sufre asiento (al haber confina- nuevo coeficiente de consolidación equivalente:
miento lateral), por lo que su módulo edométrico aparente es
infinito. La columna no debe tampoco asentar y, por tanto, no
cvrzre = cvr ⋅
[ar (λc + 2Gc ) + (1 − ar )(λs + 2Gs )] [H − (λc − λs )] − (1 − ar )(λc − λs )2 (21)
recibe ninguna carga, con lo que el suelo está sometido a una (λs + 2Gs )[H − (1 − 3ar )(Gc − Gs )]
carga mayor que la presión media aplicada pa (ecuación 2 con
σzs=0). Con el desarrollo de la consolidación, el módulo edomé- donde el superíndice zre se refiere ahora a la influencia de la de-
trico aparente del suelo va disminuyendo hasta su valor final, formación vertical y radial de la columna elástica no confinada,
drenado. En consecuencia, la carga sobre el suelo va también 1
y siendo: H = (λ c + Gc + Gs ) − (Gc − Gs ) .
disminuyendo, lo que afecta a la velocidad del proceso. ar
7
vr /cvr , cvr/cvr Ec/Es=40
6
ze
5
Ec/Es=20
czre
4
Ec/Es=10
3
Ec/Es=5
2
FIGURA 7. Columna elástica.
Deformación lateral libre. Ec/Es=1
Coeficiente de consolidación 1
equivalente (Castro, 2008). 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
Comparación con el caso
confinado (Han y Ye, 2001). Razón de sustitución (ar)
En la Figura 7 se representa la expresión (21) en función de damente la fase inicial del proceso. Como más adelante se ex-
la relación modular y el área de sustitución. Puede verse que el pone, esto se ha comprobado experimentalmente en laborato-
cociente cvrzre cvr es siempre superior a la unidad, es decir, que rio (Figura 16). En cualquier caso, incluso en esa región, la so-
lución presente supone una mejora sustancial de los
la consideración de la variación del reparto de cargas en el
resultados de las anteriores, en términos de evolución del
tiempo representa una aceleración de los asientos. Sin em-
grado de consolidación.
bargo, este efecto no es tan acusado como en la solución confi-
nada (ec. 19), representada también en la figura. En un caso tí-
3.3. COLUMNA ELASTOPLÁSTICA CON DEFORMACIÓN
pico, con ar=0,2 y Ec/Es=20, el coeficiente de consolidación
aumenta a 6 veces el básico para columna confinada, mientras RADIAL LIBRE
que para columna libre se multiplica sólo por 3,5. Para valores 3.3.1. Tensiones iniciales
muy altos (irreales) de la razón de sustitución, el factor a apli- Toda la formulación anterior, que corresponde al caso de co-
car al coeficiente de consolidación tiende a infinito para confina- lumna elástica, se debe entender como incremental, es decir
miento lateral, y en cambio decrece a la unidad para columna li- que todas las tensiones y deformaciones son incrementos so-
bre lateralmente. Esta aparente ralentización del fenómeno se bre los valores iniciales (antes de aplicar la carga pa). Por
ve compensada por el aumento de la componente instantánea tanto, las soluciones presentadas, tanto en el estado inicial, el
del asiento (que en la solución confinada no existe). final, o durante la consolidación son válidas e idénticas para
La Figura 8 muestra la evolución de la consolidación en un una rebanada de espesor diferencial situada a cualquier pro-
caso típico, analizado por Balaam y Booker (1981) mediante fundidad. En cambio, para columna elastoplástica, la condi-
elementos finitos con consolidación acoplada2. Se muestra la ción de plastificación (9a) debe imponerse con las tensiones
comparación con la solución presente, y también con la solución acumuladas a las iniciales.
confinada lateralmente, así como la que más arriba se ha deno- Por ello, el planteamiento se debe hacer para una profun-
minado “aplicación directa a columnas de grava” de la solución didad genérica z, y luego integrarse para toda la columna.
de Barron-Hansbo. La consolidación bajo carga constante es Las tensiones iniciales dependen del método de instalación de
mucho más lenta que los resultados numéricos. Por contra, la las columnas (vibrosustitución o vibrodesplazamiento), sobre
consideración de carga variable pero manteniendo el confina- todo las horizontales. Para una instalación ideal que no modi-
miento lateral conduce a un proceso excesivamente rápido, so- fique el estado natural se puede tomar:
bre todo para grados de consolidación por encima del 50%.
La solución presente muestra un mejor ajuste con los re- σ zs,i = γ 's z
sultados numéricos. Las diferencias son mayores para grados σ zc,i = γ 'c z
de consolidación inferiores al 30%, deficiencia que es compar- (22)
σ 'rs,i = σ rc,i = k0sσ ' zs,i
tida también por las otras soluciones. Balaam y Booker (1981)
lo atribuyen a que la solución de Barron no reproduce adecua- ui = 0
Ec/Es=10
Asiento as=0.25
20 inicial
υc=υs=0.3
Asiento elástico (%)
40
(cv,eq)
60
elastoplástica, el primer paso es definir el inicio de la plastifi- tica puede hacerse de la forma descrita en el Apartado 2.3, en
cación (punto Y). Este punto se alcanza para un grado de con- función de la sobrepresión intersticial media, lo que resuelve
solidación más pequeño cuanto menor es la profundidad (es el problema en este rango. El estado final (punto Fp,UC) se ob-
decir, cuanto más cerca del origen se sitúe el punto inicial I. tiene para sobrepresión intersticial nula. En lo que sigue, se
Por tanto, la plastificación comenzará en la cabeza de la co- han despreciado las componentes elásticas de la deformación
lumna, progresando hacia abajo al continuar la consolidación. en la zona plástica, como se ha comentado.
La condición de plastificación es: Igual que en el caso elástico, la ecuación de la consolida-
σ rc, y σ rc,i + ∆σ rc, y 1 − senφ c ción (10) se sigue manteniendo. El coeficiente de consolida-
= = = kac (23) ción equivalente es ahora:
σ zc, y σ zc,i + ∆σ zc, y 1 + senφ c
ar λ J
donde el subíndice “y” se refiere al momento de plastificación. (1 − ar ) + s+
Entre los puntos U y Fe, todas las tensiones varían lineal- (λ s + 2Gs ) kψc kac
–) (ec. 20), por lo cvrzrp = cvr ⋅ (26)
mente con la sobrepresión intersticial media (u ar a
1 + 1 − ar + r
(1 − ar ) kψc
que el grado de consolidación puede aplicarse a cualquier kac
componente en la forma:
∆σ y = ∆σ u + U ye (∆σ ef − ∆σ u ) (24) Gs + ar (λ s + Gs )
con : J = λ s + .
siendo ∆σ ef
el valor final de la componente elegida en la hipó- (1 − ar )kψc
tesis de comportamiento puramente elástico (punto Fe), y ∆σu Así como en régimen elástico (confinado o no) el coeficiente
en el estado sin drenaje (punto U). de consolidación dependía de la relación modular columna-
Aplicando esta expresión a las tensiones radial y vertical suelo (ecs. 19 y 21), ahora los únicos parámetros elásticos que
en la columna e introduciéndolas en (23), se puede despejar el intervienen con los coeficientes de Poisson. Esto es consecuen-
grado de consolidación elástico en el momento de plastifica- cia de haber despreciado las deformaciones elásticas frente a
ción, U ey: las plásticas.
U ye =
(σ rc,i − kacσ zc,i) + ∆σ rc,u − kac∆σ zc,u En la Figura 9 se representa la expresión anterior. Puede
(25) verse que, al contrario de lo que ocurre en el tramo elástico,
kac ∆σ zc
e
, f − ∆σ rc, f + ∆σ rc,u − kac ∆σ zc,u
e
en el que el coeficiente de consolidación se multiplicaba por
En esta expresión, sólo el paréntesis del numerador, que es valores de 3 a 6, al pasar al rango plástico, el factor decrece
siempre positivo, depende de la profundidad, por lo que U ey es por debajo de la unidad, lo que supone un brusco retardo del
mínimo en la superficie, como se había anticipado. Alternati- proceso. Esto se muestra claramente en un caso típico (Figura
vamente, puede considerarse como la expresión que define la 10). El punto de plastificación, que como se comentó antes va-
profundidad alcanzada por la zona plástica en función del ría con la profundidad, supone un quiebro brusco en la evolu-
grado de consolidación elástico. ción del proceso.
Este mismo cambio se refleja en el reparto de tensiones
3.3.3. Rango plástico entre suelo y columna. La plastificación de ésta supone una
El punto final real no será Fe, sino Fp,UC, por lo que el grado disminución importante de la rigidez relativa, lo que hace que
de consolidación real será diferente del dado por (25). La de- a partir de este momento la columna va soportando cargas re-
terminación de las tensiones y deformaciones en la zona plás- lativamente menores, como se pone de manifiesto en un caso
1.00
ψc=20°
0.95
vr /cvr
0.90
czrp
ψc=0°
0.85
ψc=10°
φc=50°
φc=40°
φc=30°
0.80 FIGURA 9. Deformación lateral
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 libre. Coeficiente de
consolidación equivalente
Razón de sustitución (ar) en el rango plástico.
(Castro, 2008).
típico representado en la Figura 11. Puede verse que el factor dados por las soluciones elásticas, confinadas o no, y se ajus-
de concentración de tensiones cesa bruscamente de crecer al tan muy bien a los de la solución plástica expuesta.
producirse la plastificación, variando muy poco a partir de en-
tonces.
4. ESTUDIO EXPERIMENTAL EN LABORATORIO
En la Figura 12 se muestra el estado final en el caso gene-
ral, expresado en la forma habitual mediante el factor de re- Como se ha ido viendo a lo largo de los apartados anteriores,
ducción de asientos, como para las soluciones anteriores. existen numerosas evidencias experimentales de medidas de
Puede verse que ahora los asientos finales concuerdan muy asientos finales, y menos de reparto de tensiones entre suelo
bien con la os valores medios de los casos reales. Otro tanto y columna. Aún menos medidas hay del proceso de consolida-
ocurre con el factor de concentración de tensiones, que se re- ción. Por ello, se decidió completar el estudio teórico del pro-
presenta en la Figura 4. Los casos reales (rango definido por ceso expuesto con un estudio en laboratorio (Cimentada,
Barksdale y Bachus, 1983) dan valores muy por debajo de los 2009; Cimentada et al., 2011).
U = 1 - u / uu
Plastificación de la columna
0.2
z=0 m
0.4 pa=50 kPa
ar=0.11
Ec/Es=40
Es=1000 kPa z=5 m
0.6 νc=νs=0.33
φc=40°
ψc=0° z=10 m
γ's=10 kN/m3
0.8 γ'c=10 KN/m3
k0s=0.6
Elástico
σz/pa
3.0
Ec/Es=40
νc=νs=0.33 ψ=20°
2.5 ar=0.25
z=0
φ=40° ψ=0°
2.0 ψ=0–20°
Columna
1.5
punto de plastificación
1.0
Suelo (Total)
ψ=0°
0.5
ψ=20°
Suelo (Efectiva)
1.0
Datos de campo
Ensayos laboratorio (Cimentada, 2009)
Priebe (φc=40°)
0.8 Plástico, Ec/Es=10–40, φc=40°, ψc=10°
Plástico, Ec/Es=10–40, φc=35°, ψc=5°
Factor de reducción
de asientos α=1/n
0.6
0.4
0.2
Para ello, se ha reproducido una rebanada horizontal de la ha utilizado para la columna grava de 4-5 mm de diámetro.
celda unidad en un edómetro de tipo Rowe-Barden, de 254 Para el suelo se ha usado caolín comercial, que supone un ade-
mm de diámetro (Figura 13). Se han ensayado columnas de cuado compromiso entre plasticidad y permeabilidad. El utili-
63,5 mm y 84,67 mm, lo que supone relación de diámetro de zado en este estudio tiene un límite líquido de 73%, índice de
celda a columna N=4 y N=3, respectivamente, y razones de plasticidad de 35%, y deformabilidad edométrica definida por
sustitución, ar, respectivas de 0,0625 y 0,111. unos índices Cc y Cs de 0,44 y 0,14. Para presiones entre 100 y
600 kPa el módulo de elasticidad equivalente era de 600-3000
4.1. MATERIALES Y MÉTODOS kPa y el coeficiente de consolidación medio de 2,5×10-3 cm2/s.
Lo anterior supone una escala aproximada de 1/10 respecto a En cuanto a la grava, se compactó a una densidad de 16,5
las dimensiones reales típicas. Para mantener esta escala, se kN/m3, cercana a la máxima. El ángulo de rozamiento medido
ø258.0
En la Figura 16 se presentan las isocronas radiales de presión
20°
ø274.0
20°
intersticial en un ensayo con cada geometría (N=3 y 4). Puede
ø325.0 verse que, a partir de un cierto tiempo (del orden de 1000-
2000 s) las isocronas tienen sensiblemente el aspecto predicho
20° 20° A por la solución de Barron (ec. 13 y Figura 6), pero no así en
20° 20° las primeras etapas del proceso, como indica la experiencia
pasada.
FIGURA 14. Instrumentación de la base de la célula Rowe. Sensores de La Figura 17 muestra los valores del coeficiente de conso-
presión total e intersticial. lidación equivalente, obtenido mediante ajuste de la solu-
1.0
Valores de t (s):
50.000
0.9 20.000
12.000
0.8 7.000
5.000
4.000
0.7 2.000
1.000
0.6 600
300
100
∆u/∆u0
0.5 10
2
1
0.4
0.5
0.3
0.2
0.1
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
r/re
a) N=4 (ar=0,0625)
1.0
Valores de t (s):
0.9 50.000
20.000
12.000
0.8 7.000
5.000
0.7 4.000
2.000
1.000
0.6 600
300
100
∆u/∆u0
0.5
50
10
0.4 2
1
0.5
0.3
0.2
0.1
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
r/re FIGURA 16. Ensayos en
laboratorio. Isocronas de
b) N=3 (ar=0,111)
presión intersticial.
ción de Barron-Hansbo a los resultados experimentales, ex- a apartarse de los resultados experimentales para razones
presados como cociente respecto al valor básico. Se ha deter- de sustitución altas.
minado para cada geometría (N=3-4) y cada escalón de En cuanto al reparto de tensiones entre suelo y columna, en
carga, utilizando dos métodos: la disipación de presiones in- la Figura 18 se presenta la evolución del factor de concentra-
tersticiales y la evolución de la deformación vertical. Los re- ción, SCF, a lo largo de la consolidación para cada escalón de
sultados muestran una cierta dispersión, pero reflejan la carga. Puede verse que comienza por valores muy bajos, incluso
tendencia a disminuir al crecer la relación de diámetros, N. inferiores a la unidad, como predicen los estudios teóricos para
Como comparación, se presentan los resultados teóricos la situación sin drenaje. Después hay un aumento suave, y
para confinamiento lateral (Han y Ye, 2001) y para deforma- para un tiempo entre 103 y 104 segundos se observa un cambio
ción radial libre (Castro, 2008). La solución confinada tiende de tendencia, manteniéndose sensiblemente constante; al
Estimación a partir de u
Estimación a partir de Ur
Castro (2008)
4 Han y Ye (2001)
3
zre
ceq
(a)
Factor de concentración de tensiones, SCF
30
28
26 Valores obtenidos en el ensayo
24 Castro (2008). Columna elástica
22 Castro (2008). Columna elasto-plástica
20 Han y Ye (2001)
18
16
14
12 N=4
10
8
6
4
2
0
10–1 101 103 105 107
Tiempo (s)
(b)
Factor de concentración de tensiones, SCF
30
28
26 Valores obtenidos en el ensayo
24 Castro (2008). Columna elástica
22
20 Castro (2008). Columna elasto-plástica
18 Han y Ye (2001)
16
14
12
10 N=3
8
6
4
2
0
10–1 101 103 105 107 FIGURA 18. Ensayos en
laboratorio. Evolución del
Tiempo (s) factor de concentración de
tensiones.
mismo tiempo, aparece una notable irregularidad en las cur- Como parte del estudio, se ha realizado también una
vas, debidas a oscilaciones en el registro de las tensiones en la prueba de campo, tratando únicamente de resaltar la impor-
columna, típicas de materiales granulares. Este cambio de ten- tancia de este aspecto. Para ello, se ha contado con la posibili-
dencia podría indicar la plastificación de la columna. En la dad de instrumentar una obra real, en un terraplén de la au-
misma figura se han superpuesto como comparación las curvas tovía Variante de Sueca (Valencia), con una altura de 10 m,
teóricas, observándose que la solución confinada (Han y Ye, más 2 m de precarga. La cimentación se reforzó con columnas
2001) sobrevalora en mucho el factor, y también la elástica li- de grava, de 0,80 m de diámetro, 2,80 m de espaciamiento en
bre. Sólo la consideración de la plastificación de la columna malla triangular, y 9,0 m de longitud.
(Castro, 2008) proporciona ajustes aceptables. El valor final del El terreno consta de una interestratificación de depósitos
factor de concentración toma valores en el rango 4-8, lo que aluviales, predominantemente limosos, de consistencia
concuerda con la experiencia en casos reales (Figura 4). blanda-media, de 8-10 m de potencia. Los ensayos de labora-
En cuanto a las deformaciones finales, definidas por el fac- torio e insitu muestran indicios de una cierta sobreconsolida-
tor de reducción de asientos, se han incluido en la Figura 12 ción, con valores del coeficiente A de Skempton en el rango de
junto con los resultados teóricos y las experiencias reales, ob- +0,20 a -0,17 y una resistencia al corte sin drenaje de 60 a
servándose una buena concordancia, si bien se trata de sólo 160 kPa.
dos puntos, y para razones de sustitución moderadas. La prueba consistió en la instalación de siete columnas, en
celda exagonal (Figura 19). Se instalaron piezómetros en los
baricentros de las subceldas triangulares, a profundidades de
5. EFECTOS DEL PROCESO DE INSTALACIÓN
4 m y 7 m, en las que los sondeos hechos para instalar los pie-
En todo el estudio anterior, no se ha tenido en cuenta el efecto zómetros indicaban niveles más impermeables (el nivel freá-
del proceso real de instalación de las columnas en el terreno, tico se detectó a 3 m de profundidad). El orden de construc-
sino que, tanto en los estudios teóricos como experimentales ción indicado permitía que en la primera columna, instalada
se ha considerado una instalación ideal, que no modifica el es- en un terreno virgen, se dispusiera de medidas a tres distan-
tado tensional del terreno. cias diferentes de la misma (del orden de 1,5, 3,0 y 4,5 m).
0
Profundidad (m)
0
Limos arenosos compactos (cu=80kPa)
1 CPTU–31
Arcilla limosa (cu=50kPa)
–5
2
N.F.
Arcilla blanda (cu=25kPa)
Profundidad (m)
4
–10
Arcillas (cu=60kPa)
6
Arenas (ϕ=30°)
–15
8
10
–20
Arena con niveles arcillosos (ϕ=35°) 0 5 10 15 20
Resistencia por pinta, qc (MPa)
COLUMNAS DE GRAVA
6
PZ01–1–PZ01–2
4m
1 5
6.60m
7m
2 4
PZ03–1–PZ03–2
PZ03–1
PZ01–1 PZ02–1
3
b) Disposición de piezómetros
FIGURA 19. Prueba de c) Orden de instalación
instalación de columnas.
180
PZ 01
140 PZ 02
PZ 03
Presiones intersticiales medidas (kPa)
120
100
80
60
40
20
0
16:00:00 16:30:00 17:00:00 17:30:00 18:00:00 18:30:00
0
Profundidad (m)
2
4
6
8
10
16:00:00 16:30:00 17:00:00 17:30:00 18:00:00 18:30:00
Intensidad (Amperios)
200
150
100
50
0
16:00:00 16:30:00 17:00:00 17:30:00 18:00:00 18:30:00
Tiempo
Pico de presiones
intersticiales (kPa)
100 4m
7m
COLUMNA
50
25
0
0 1.5 3.0 4.5
FIGURA 21. Ensayo de
instalación. Picos de presiones Distancia al eje de la columna (m)
intersticiales.
3•0
Venmans3
1•5
1•0
Castro35
0•5 range
analíticas, considerando las deformaciones vertical y radial factores de corrección en el rango 2-6. Nuevamente, la plas-
de columna y suelo, que llevan a la posibilidad de uso de la tificación de la columna supone un drástico cambio, redu-
solución de Barron-Hansbo, mediante la introducción de coe- ciendo el coeficiente de consolidación incluso ligeramente
ficientes de consolidación equivalentes, que resultan de mul- por debajo del básico.
tiplicar su valor básico por un cierto factor para las distintas Los aspectos anteriores han sido contrastados con los re-
hipótesis. La hipótesis confinada conduce también en este sultados de los ensayos de laboratorio, en los que se ha repro-
aspecto a una sobrevaloración de dicho factor, que se eleva ducido una rebanada horizontal de la celda unidad a escala
hasta en más de un orden de magnitud. La consideración de aproximada 1/10, tanto en cuanto a asientos y concentración
la deformación elástica radial reduce algo este efecto, con de tensiones como al proceso de consolidación.