1 de 35

Estabilidad de terraplenes sobre columnas de suelo-cemento

Stability of embankments over cement deep soil mixing columns
Pablo Morilla Moar1* y Svetlana Melentijevic2

Palabras clave Sumario

suelo-cemento; Existen diversos métodos de mejora del terreno con características geotécnicas desfavorables. Uno de es-
soil mixing; tos métodos es el de estabilización profunda mediante columnas de suelo-cemento, que se suele denominar
columnas; por sus siglas en inglés DSM (deep soil mixing). En obras lineales proyectados sobre suelos blandos propor-
terraplenes; cionan principalmente una aceleración de la construcción, mejora de la estabilidad, aumento de la capacidad
estabilidad; portante y reducción de los asientos totales y diferenciales.
Existen dos métodos tradicionales de diseño, el método escandinavo (columnas blandas y semi-rígidas) y
el método japonés (columnas rígidas). Los métodos tradicionales de cálculo han sido cuestionados por varios
autores que han demostrado mediante ensayos de laboratorio y métodos numéricos que sobrestiman los fac-
tores de seguridad debido a que no asumen los tipos de rotura más habituales en este tipo de mejora.
En este artículo se realiza un breve repaso de las metodologías de diseño tradicionales para mejora de la
estabilidad de terraplenes sobre columnas de suelo-cemento construidas en suelos blandos, de los estudios
realizados para averiguar los tipos de rotura más probables de dicho tipo de columnas y de metodología que
permite reducir la sobrestimación de los métodos de equilibrio límite y de la que permite detectar los tipos de
rotura apropiados en columnas suelo-cemento (métodos de análisis numérico).
Finalmente se analiza un caso práctico aplicando tanto el método de equilibrio límite como el método de
elementos finitos, cuyos resultados confirman la sobrestimación de los métodos tradicionales en la evaluación
de laestabilidad de terraplenes sobre columnas de suelo-cemento.

Keywords Abstract

deep soil-mixing; The deep soil mixing (DSM) is one of the ground improvement methods used for the construction of embank-
soil-cement columns; ments over soft soils.
embankments; DSM column-supported embankments are constructed over soft soils to accelerate its construction, improve
stability; embankment stability, increase bearing capacity and control of total and differential settlements.
There are two traditional design methods, the Japanese (rigid columns) and the Scandinavian (soft and se-
mi-rigid columns). Based on laboratory analysis and numerical analysis these traditional approaches have been
questioned by several authors due to its overestimation of the embankment stability considering that the most
common failure types are not assumed.
This paper presents a brief review of traditional design methods for embankments on DSM columns construc-
ted in soft soils, studies carried out determine the most likely failure types of DSM columns, methods to decrea-
se the overestimation when using limit equilibrium methods and numerical analysis methods that permit detect
appropriate failure modes in DSM columns.
Finally a case study was assessed using both limit equilibrium and finite element method which confirmed the
overestimation in the factors of safety on embankment stability over DSM columns.

1. INTRODUCCIÓN Japón, mientras que el método seco tiene su origen en los

países escandinavos.
El método de mejora del terreno mediante estabiliza- Los métodos principales de diseño son entonces el mé-
ción profunda o deep soil mixing (DSM) consiste en la todo japonés y el método denominado escandinavo. Sus
mezcla mecánica de un ligante con el terreno natural. Los respectivas guías de referencia son la editada por el Insti-
ligantes pueden introducirse tanto por el “método húme- tuto Tecnológico de Desarrollo de Costas de Japón en 2002
do”, donde se utiliza una lechada de ligante (generalmente (CDIT en sus siglas en Inglés) y el EuroSoilStab (2002).
cemento) y agua, o por el “método seco” donde el ligante Como ligantes en el DSM se utilizan principalmen-
en polvo reacciona químicamente con el agua intersticial te agentes estabilizadores como cal viva, cal apagada, ce-
durante el curado. El método húmedo tiene su origen en mento o una combinación de estos. Recientemente se han
incorporado otros ligantes como yeso, cenizas volantes, es-
corias de alto horno, etc.
*
Corresponding author: [email protected] Para la estabilización con grupos de columnas aisladas
1
Eptisa, Madrid, España. se puede utilizar maquinaria con herramienta de mezcla
2
Grupo Rodio-Kronsa, Madrid, España. de un eje (figura 1), o se pueden usar herramientas de ejes

Ingeniería Civil 175/2014 | 1

23 de 35

Estabilidad terraplenes.indd 1 09/09/14 11:43

 Estabilidad de terraplenes...

múltiples para crear un grupo de columnas solapadas de

una sola maniobra (figura 2).

Figura 3. Ejemplos de patrones de instalación de columnas (Euro-

SoilStab, 2002).

Para comparar varios patrones de instalación de co-

lumnas mediante el área tratada, y para evaluar las carac-
terísticas conjuntas del suelo tratado y el no tratado que lo
rodea, se definió la relación de área mejorada ( ):
Figura 1. Construcción de una columna de suelo-cemento (Filz, 2007).

= [1]
´

Donde:
= área del suelo tratado con columnas de suelo-ce-
mento (m2).
´ = área total correspondiente (m2).
El método de estabilización profunda se podría ejecu-
tar hasta 100 % de la relación de área mejorada, en el caso
de tratar el bloque del terreno completo. En las obras linea-
les, en los países escandinavos el rango de valores de as está
Figura 2. Herramienta de DSM con tres barrenos continuos para entre 10 y 30%, mientras que en Japón se varía de 30 a 50%
hacer columnas solapadas (Topolnicki, 2004). debido a su aplicación en la zona sísmica de alto riesgo (To-
polnicki, 2004).
Según las citadas guías EuroSoilStab (2002) y CDIT
(2002) las aplicaciones más habituales de la estabilización 1.2. Propiedades geomecánicas de las columnas de
profunda de suelos son: suelo-cemento

• Incrementar la resistencia del suelo blando con el fin Para el diseño de las columnas DSM de forma habitual
de aumentar la estabilidad de terraplenes y la capa- se asume un valor de resistencia al corte sin drenaje de
cidad portante del terreno, reducir los empujes acti- cu=0,5∙qu. Sin embargo, Broms (2004) afirma que este va-
vos en muros de contención, prevenir la licuefacción lor puede ser muy alto cuando el suelo estabilizado no está
y reducir movimientos horizontales. completamente saturado. Además, la resistencia al corte
• Remediación de terrenos contaminados creando ba- dependerá de muchos factores, por ejemplo cu se incre-
rreras medioambientales y, estabilizando el terreno menta de manera proporcional al contenido en cemento y
contaminado o creando barreras geohidrológicas. con el descenso del límite líquido del terreno natural. To-
• Sistemas de contención para taludes, excavaciones, polnicki (2004) presenta una serie de correlaciones de cu,
barreras de impermeabilización, deslizamientos de para suelos tratados con cemento por el método húmedo,
laderas, etc. que dependerán de la resistencia a la compresión simple
• Cimiento de estructuras (edificación, depósitos, alcanzada cu=(0,4-0,5)qu para qu<1MPa; cu=(0,3-0,35)qu
etc.). para qu=(1-4)MPa y cu=0,2qu para qu>4MPa.
Para el módulo de elasticidad (Ecol) de las columnas
1.1. Patrones de instalación de las columnas existen muchas correlaciones empíricas mediante la resis-
tencia a compresión simple. Filz y Navin (2006) aconsejan
La mezcla de suelo puede realizarse hasta una relación una relación para método húmedo de E=300qu, pero puede
de un 100% entre el suelo tratado y sin tratar, como es el variar entre 150-1000qu. Para el método seco esta relación
caso de los tratamientos superficiales, pero en el caso de la varía entre 65qu y 250qu. En la tabla 1 pueden verse dife-
mezcla profunda, esta relación es menor. De este porcenta- rentes valores de qu dependiendo del terreno a tratar y del
je de suelo tratado dependerán las capacidades mecánicas contenido en cemento.
del método seleccionado. La disposición de las columnas En general, las columnas de suelo-cemento poseen ele-
se diseñará dependiendo del objetivo del proyecto, las ca- vadas resistencias, y deformaciones en rotura muy bajas,
racterísticas del terreno natural, las condiciones del empla- del orden de 1 al 2%. En el caso de resistencia confinada
zamiento de la obra, y el coste del tratamiento (figura 3). (in situ) incluso con un confinamiento bajo, la resistencia
Las mallas triangulares o cuadradas de columnas ais- residual de la columna se incrementa hasta casi el 80% de
ladas o combinadas, se utilizan habitualmente en obras li- la resistencia a la compresión simple por lo que es aconse-
neales para la reducción de asientos y para la mejora de la jable utilizar la resistencia residual en el diseño (Tatsuoka
estabilidad. y Kobayashi, 1983).

2 | Ingeniería Civil 175/2014

24 de 35

Estabilidad terraplenes.indd 2 09/09/14 11:43

 Estabilidad de terraplenes...

Tabla 1. Valores típicos de resistencia y permeabilidad para

diferente factor de cemento y diferentes tipos de suelo tratados
por el método húmedo (Adaptado de Topolnicki, 2004).
Factor
RCS a 28 días, quf
Tipo de Suelo cemento,
(MPa)
α (kg/m3)
Fango 250-400 0,1-0,4
Turba, arcillas y limos orgánicos 150-350 0,2-1,2
Arcillas blandas 150-300 0,5-1,7
Arcillas duras/medias 120-300 0,7-2,5
Figura 4. Patrones de rotura asumidos por el método de diseño
Limos y arenas limosas 120-300 1,0-3,0
tradicional (Kitazume y Maruyama, 2006).
Arenas finas a medias 120-300 1,5-5,0
Arenas gruesas y gravas 120-250 3,0-7,0 Mientras que el método de diseño japonés asume que
las resistencias al corte de las columnas y del suelo sin tra-
2. METODOS TRADICIONALES DE DISEÑO Y ANÁLISIS tar entre ellas no se movilizan al mismo tiempo. La defor-
mación axial de las columnas de suelo tratado por el DSM
Tal y como se ha indicado anteriormente, existen dos es muy reducida en comparación con la deformación en
métodos de diseño tradicional, el método de diseño escan- rotura del suelo sin tratar.
dinavo, y el método de diseño japonés, ambos tienen as- Esto se debe a que en el método escandinavo se utilizan
pectos en común pero también diferencias importantes. más las columnas de suelo-cal y cal-cemento, mientras que
Ambos métodos consideran el conjunto de las co- en el japonés es más habitual el uso de columnas de suelo
lumnas de DSM y el suelo sin estabilizar como un mate- cemento. Por lo tanto este último es el más adecuado para
rial compuesto. Así mismo, basan el cálculo de estabilidad el diseño de las columnas de suelo cemento.
del terraplén en análisis de estabilidad de taludes mediante Para los análisis de estabilidad se analiza una superficie
métodos de equilibrio límite. de corte circular y los parámetros para dicho análisis se cal-
La estabilidad de un terraplén mejorado con colum- culan transformando el área mejorada por las columnas, en
nas DSM se analiza en el estado sin drenaje, es decir justo un suelo equivalente con características medias de ambos.
después de la construcción dado que representa el estado
más desfavorable o más crítico, debido al aumento con el 3. INVESTIGACIONES SOBRE MECANISMOS DE
tiempo de la resistencia al corte de las columnas y del sue- ROTURA DE COLUMNAS DSM
lo entre ellas. Se utiliza una resistencia al corte sin drenaje
“equivalente”, es decir, se hace una media ponderada entre Se han realizado una serie de análisis de laboratorio
las resistencias del suelo sin tratar y de las columnas. mediante ensayos de centrifuga con el fin de investigar dis-
tintos tipos de rotura en columnas de suelo-cemento (fi-
cu = cu,col as + cu,suelo (1 as ) [2] gura 5).
Los ensayos de centrifuga se han realizado principal-
Donde: mente por investigadores japoneses, y estudian el compor-
cu,col = es la resistencia al corte sin drenaje de la co- tamiento de grupos de columnas ante una carga vertical y
lumna. el empuje horizontal correspondiente.
cu,suelo = es la resistencia al corte sin drenaje del suelo
sin tratar entre las columnas (kN/m2).
as = relación de área mejorada (ver ecuación 1).
Según Kitazume y Maruyama, 2006, se analizan dos pa-
trones de rotura, la estabilidad externa y la estabilidad in-
terna (figura 4):

– Para la estabilidad externa se calcula el factor de

seguridad al deslizamiento en bloque, en el cual el
conjunto columnas-suelo sin tratar entre ellas se
mueven de forma horizontal sobre una capa rígida,
sin ningún movimiento relativo o reordenación de
las columnas.
– En el análisis de la estabilidad interna, se comprueba
la resistencia a cortante de las columnas evaluándo- Figura 5. Aparato de centrifuga para geotecnia (Kitazume y Ma-
se mediante un análisis por círculo de deslizamiento. ruyama, 2007).

La diferencia principal entre estos métodos es que en el Kitazume et al. (1991 y 2000) realizaron una serie de en-
método escandinavo se supone que existe una interacción sayos de centrifuga sobre modelos de grupos de columnas
total entre las columnas y el suelo que las rodea. La resis- y descubrieron que respecto a la estabilidad externa una
tencia al corte de pico de las columnas se moviliza al mis- rotura tipo vuelco (figura 6) puede darse con más frecuen-
mo tiempo que la resistencia al corte del suelo entre ellas. cia que una rotura por deslizamiento en bloque de la zona

Ingeniería Civil 175/2014 | 3

25 de 35

Estabilidad terraplenes.indd 3 09/09/14 11:43

 Estabilidad de terraplenes...

tratada con columnas de DSM (figura 4a). En este tipo de para reducir los desplazamientos laterales en la zona de
rotura, las columnas vuelcan como piezas de dominó una a corte y de esa manera aumentar la resistencia al corte y a
una. Por lo tanto, el método tradicional de diseño, que no flexión de las columnas.
tiene en cuenta la rotura por vuelco, podría sobrestimar la
estabilidad externa, ya que, es razonable pensar que el te-
rreno tienda a fallar por un tipo de rotura que, bajo ciertas
condiciones, tenga una resistencia menor.

Figura 8. Tipos de roturas en columnas de DSM dependiendo de la

posición de la superficie de deslizamiento (Broms, 2004).

Las columnas aisladas tenderán a fallar antes por vuel-

co o flexión que por cortante o deslizamiento en bloque.

4. PROPUESTAS ALTERNATIVAS DE ANÁLISIS DE

Figura 6. Rotura por vuelco en ensayo de centrifuga (Kitazume y ESTABILIDAD
Maruyama, 2006).
El método de los perfiles.
En el caso de la estabilidad interna, los ensayos de cen-
trifuga (figura 7) revelaron que las columnas de suelo-ce- Este método desarrollado por Barksdale y Bachus
mento muestran varios modos de rotura: a cortante, a (1983) para el diseño de columnas de grava puede aplicar-
flexión y a tracción (figura 8). Debido a que la resistencia se a las columnas DSM analizadas por métodos de equili-
a tracción y a flexión, de columnas de suelo-cemento, es brio límite.
mucho menor que la resistencia a cortante, el método tra- Para reproducir el comportamiento de una malla real
dicional de diseño, basado únicamente en la resistencia al en un programa en 2D se reemplazan las columnas por
corte (figura 4b), podría sobrestimar de manera importan- franjas verticales equivalentes en perfil, es decir, la zona
te la estabilidad interna al asumir un tipo de rotura menos mejorada se simplifica como una serie de franjas alternan-
propenso a producirse. Mediante estos ensayos también se do con suelo sin tratar. El espesor de las franjas se calcula
verificó que las columnas no fallan de forma simultánea, dependiendo de la disposición y de características geomé-
sino que lo hacen de una en una en una secuencia que va tricas de las columnas (figura 9).
desde la columna situada al pie del talud hacia la situada en
el eje del terraplén.

Figura 7. Rotura por flexión en ensayo de centrifuga (Kitazume y

Maruyama, 2007).

Broms (2004) afirma que debe tenerse en cuenta que

la eficacia de las columnas aisladas es baja cuando la carga
axial de la columna también es baja. Las columnas aisla-
das deberían colocarse bajo el cuerpo del terraplén, donde
están cargadas axialmente y donde la resistencia de las co- Figura 9. Conversión de columnas dispuestas a tresbolillo, en
lumnas se utilizará al completo. Mientras que bajo los de- franjas equivalentes para cálculo de estabilidad de taludes (Filz
rrames de los terraplenes se aconseja solapar las columnas y Navin 2006).

4 | Ingeniería Civil 175/2014

26 de 35

Estabilidad terraplenes.indd 4 09/09/14 11:43

 Estabilidad de terraplenes...

En el modelo debe incluirse una “capa ficticia” (figura Caso 1: Terraplén sobre grupo de columnas aisladas
9a) en el límite entre el terraplén y el terreno mejorado con
columnas de suelo-cemento, con el fin de generar una con- Un grupo de columnas de suelo-cemento del diámetro
centración de tensiones en las columnas. de 0,9m se ha dispuesto en malla cuadrada con espaciado
La capa debe tener un elevado peso específico sobre las de centro a centro de 1,8m bajo toda la superficie del terra-
franjas de las columnas y uno negativo sobre el terreno sin plén, resultando en una relación de área mejorada (as) del
tratar. Estos pesos específicos se calculan basándose en el 20%. Las columnas se empotran 0,6 m en el sustrato are-
espesor de la capa, las relaciones de tensiones y la carga del noso (figura 11).
terraplén para poder aplicar la carga correcta en las colum- Utilizando el método de los perfiles las columnas se re-
nas y el suelo que las rodea. (Navin, 2005). presentaron en los programas de cálculo como franjas ver-
En la figura 10 se puede observar un modelo de terra- ticales de 0,4 m paralelas al eje de terraplén y con espaciado
plén sobre columnas convertidas en franjas y con la capa entre centros de 1,8 metros, para reproducir la relación de
ficticia realizado en el programa de equilibrio límite Slide. área mejorada del 20%. En la tabla 2 se resumen los pará-
metros de diseño utilizados para el análisis.

Figura 10. Capa ficticia de 15 cm aplicada sobre el terreno mejo-

rado (Slide).
Figura 11. Sección y planta del caso 1, terraplén sobre columnas
Métodos de análisis numérico (MEF, MDF) aisladas (adaptado de Navin y Filz, 2006).

Los métodos de análisis numérico detectan tanto la ro- Caso 2: Terraplén sobre columnas aisladas en la zona
tura por deslizamiento circular del suelo compuesto, como central de terraplén y con paneles continuos de columnas
los métodos de equilibrio límite, como las roturas por fle- solapadas bajo los derrames
xión y vuelco. Además, ofrecen información acerca de las
tensiones inducidas en las columnas, desplazamientos ver- La separación entre paneles en la zona de los derrames
ticales y horizontales, presiones intersticiales, resistencias mantiene la relación de área mejorada del 20% como en el
movilizadas, etc. caso de columnas aisladas.
La resistencia al corte de la zona mejorada con paneles
5. CASO PRÁCTICO se realizó mediante una media ponderada entre las colum-
nas, el suelo sin tratar entre ellas (ecuación 2). También se
Se ha utilizado un caso presentado por Filz y Navin incluyeron 4 juntas de solape (para reproducir una posible
(2006) en el que se analiza la estabilidad de un terraplén en rotura por cizalla de los paneles) de 20 cm de ancho, cuya
diferentes casos: (0) terraplén sobre suelo sin mejorar, (1) resistencia es la mitad de la asignada a los paneles.
terraplén sobre suelo blando mejorado con columnas ais- Los cálculos de estabilidad por el método de equili-
ladas de suelo-cemento en una malla cuadrada (figura 11) brio límite (MEL) se realizarán con el programa Slide por
y (2) terraplén sobre suelo mejorado con columnas aisladas el método de Bishop simplificado, mientras que los cálcu-
en la zona central de terraplén y paneles de columnas sola- los realizados por el método de elementos finitos (MEF)
padas bajo los derrames (figura 12). se efectuaron con el programa Plaxis 2D. Los resultados se
El modelo estudiado es un terraplén de 5,5 metros de compararan con los obtenidos por Filz yNavin (2006) los
altura, inclinación de taludes 2H:1V, cimentado sobre una cuales usaron el programa de equilibrio límite UTEXAS4
capa de arcilla blanda de 8,5 metros de espesor que esta so- (método de Bishop) y el programa de diferencias finitas
bre un sustrato arenoso denso. El nivel freático se localiza (MDF) FLAC. Todos los materiales se han modelizado
en el límite superior de la capa de arcilla blanda. usando el criterio de rotura de Mohr-Coulomb.

Tabla 2 . Parámetros geomecánicos de los materiales del modelo de estudio.

Material γ (kN/m3) E (kPa) Eref (kPa) Ef (kPa) Eincrement (kPa) ν c (kPa) cu,ref (kPa) cu,increment (kPa) Φ (˚)
Terraplén 19,6 29900 - - - 0,30 0 - - 35

Relleno de arena 18,1 12000 - - - 0,33 0 - - 30

Arcilla Blanda 15,1 200cu 2040 4120 245 0,45 - 10,20 1,22 0

Sustrato arenoso 22,0 47900 - - - 0,26 0 - - 40

Columnas 15,1 207000 - - - 0,45 689 - - 0

Ingeniería Civil 175/2014 | 5

27 de 35

Estabilidad terraplenes.indd 5 09/09/14 11:43

 Estabilidad de terraplenes...

columnas modelizadas como franjas introduciendo una

capa ficticia según el método de los perfiles (1c).

1a) Cálculo con los parámetros equivalentes

columnas-suelo
Para calcular la resistencia al corte sin drenaje equiva-
lente (suelo-columnas) se ha utilizado la resistencia al cor-
te de la arcilla en el centro del estrato, la resistencia al corte
de la columna, la relación de área mejorada y la ecuación 2;
el resultado es cu =150kN/m2.
Figura 12. Sección y planta del caso 2 (adaptado de Navin y Filz,
2006).

Para los cálculos de estabilidad por el método de ele-

mentos finitos (MEF) se han modelizado las situaciones
estudiadas según un modelo de deformación plana con 6
etapas constructivas, una de instalación del tratamiento de
DSM y 5 de construcción del relleno (la primera de 1,5m
de espesor y las siguientes de 1,0m).

Caso 0. Terreno sin mejorar

Este cálculo es el primer paso para un diseño del tipo

de actuación necesaria en un terraplén sobre suelos blan-
dos: refuerzo del terreno por estabilidad, por mejora de
capacidad portante, disminución de asientos o una combi-
nación de estas (figura 13).
Figura 14. Caso 1a terreno equivalente, FS=4,3.

1b) Cálculo con las columnas como franjas

Se introducen las columnas como franjas verticales de
0,4 metros de ancho y separadas 1,8m y se realiza el estudio
por el método de equilibrio límite (figura 15a) y de elemen-
tos finitos (figura 15c).

1c) Cálculo con las columnas como franjas y capa ficti-

cia de reparto de tensiones
En el este caso se analiza el modelo por el método de
equilibrio límite pero incluyendo la capa ficticia para re-
partir la tensión del terraplén entre las columnas y el suelo
sin tratar entre ellas (figura 15b). Los pesos específicos ob-
tenidos para la capa ficticia son:
col
f = 2.673 kN/m3
suelo
f = -0,19 kN/m3

Figura 13. Cálculo de estabilidad para talud sin mejorar (a) MEL
(Slide) con FS=0,875 (b) MEF (Plaxis), el talud falla en la 4ª etapa de
construcción.

Caso 1. Grupos de Columnas

En el caso 1 se aborda la solución de estabilización con

grupos de columnas bajo todo el terraplén. Se analiza la
estabilidad por el método tradicional de los parámetros
equivalentes del conjunto terreno-columnas (1a), con las
columnas DSM como modelizadas franjas (1b) y con las

6 | Ingeniería Civil 175/2014

28 de 35

Estabilidad terraplenes.indd 6 09/09/14 11:43

 Estabilidad de terraplenes...

carga axial) hacia la situada más cerca del eje del mismo
(figura 16a).
La columna más cargada axialmente es la situada en el
centro del talud, en este caso la carga axial en dicha colum-
na sería de 468kPa, mientras que las columnas con más es-
fuerzo cortante son las situadas más cerca del pie del talud,
con un cortante de 43kPa.

Caso 2. Paneles de columnas solapadas

2a) Cálculo de estabilidad con paneles

En este caso se introducen pantallas bajo los derrames,
incluyendo cuatro juntas de solape, se realiza el cálculo con
Slide y Plaxis (figuras 17a y 17c respectivamente).

2b) Cálculo de estabilidad con paneles de columnas so-

lapadas y capa ficticia.
Como el caso anterior, pero aplicando la capa ficticia a
la zona donde hay columnas aisladas únicamente para el
análisis por equilibrio límite (figura 17b).

Figura 15. Caso 1b, columnas en todo el terraplén. (a) caso 1a, MEL
(Slide) FS=4,57, (b) caso 1b, MEL (Slide) con capa ficticia FS=2,83 y
(c) caso 1 por MEF (Plaxis) FS=1,46.

El análisis con Plaxis revela un potencial desarrollo de

rotura por vuelco de las columnas (figura 16a). Se observa
que la distribución de las tensiones efectivas medias (figura
16b) se concentra en las columnas situadas bajo el cuerpo
central del terraplén, estando prácticamente sin tensión las
columnas cercanas al pie.

Figura 16. Caso 1, (a) vuelco de columnas (factor de escala 20), (b)
concentración de tensiones efectivas en las columnas (MEF, Plaxis).

Esto coincide con la observación de Kitazume y Maru-

yama (2007), la rotura no se produce de forma simultánea Figura 17. Caso 2, (a) paneles bajo los derrames MEL (Slide),
en todas las columnas, sino que fallan una a una empezan- FS=4,07, (b) paneles y capa ficticia, MEL (Slide) FS=2,85, (c) análisis
do por la situada más cerca del pie del talud (la de menor por MEF (Plaxis) con paneles FS=2,73.

Ingeniería Civil 175/2014 | 7

29 de 35

Estabilidad terraplenes.indd 7 09/09/14 11:43

 Estabilidad de terraplenes...

Tabla 3. Comparativa de los factores de seguridad obtenidos con los resultados de Filz y Navin (2006).
Factor de seguridad
Filz y Navin (2006)
Caso Descripción MEL MEF
(SLIDE) (PLAXIS 2D) MEL MDF
(UTEXAS4) (FLAC)

(a) Zona mejorada como terreno equivalente 4,31 - -

4,4
Caso 1 (b) Zona mejorada como columnas representadas por franjas 4,57
1,46 1,4
(c) Terreno mejorado con franjas y capa ficticia 2,83 -
(a) Paneles con juntas de solapes 4,07 4,4
Caso 2 2,73 3,1
(b) Paneles con juntas de solapes y capa ficticia 2,85 -

flexión o vuelco de las columnas. Concluyeron que los méto-

dos tradicionales podrían sobrestimar significativamente la
estabilidad al analizar tipos de rotura menos probables.
Se ha podido demostrar que el uso de cálculos de esta-
bilidad por equilibrio límite usados por los métodos tra-
dicionales (con parámetros equivalentes suelo-columna)
sobreestiman los factores de seguridad, por lo tanto no pa-
recen adecuados para evaluar la estabilidad de terraplenes
sobre suelos mejorados con columnas de DSM.
Figura 18. Deformación de los paneles (factor escala de 60, MEF, En el caso de evaluarse la estabilidad por el método de
Plaxis). equilibrio límite, es preferible utilizar la técnica de la capa
ficticia ya que reduce de manera apreciable la sobrestima-
En la tabla 3 se presenta un resumen de los resultados ción y en algún caso arroja valores similares al análisis por
obtenidos en los análisis de estabilidad mediante equilibrio elementos finitos. Pero con el método de elementos finitos
límite y elementos finitos en los casos 1 y 2, comparados es posible detectar otros tipos de rotura, además del desli-
con los resultados obtenidos por Filz y Navin (2006). zamiento y el corte, como la rotura por vuelco ocurrida en
Como puede apreciarse, los valores obtenidos por el el ejemplo estudiado. Como apuntaron Kitazume y Maru-
método de equilibrio límite sobreestiman notablemente, yama (2007), la instalación de columnas de DSM tiene el
excepto en el caso 2b, los factores de seguridad, respecto a efecto de cambiar el tipo de rotura, de la rotura por desliza-
los obtenidos por métodos numéricos. miento circular a una rotura por vuelco de las columnas o
En el caso 1, tanto el modelo con área equivalente de una rotura por corte a una rotura por flexión. Además,
(1a) como el modelo con columnas aisladas (1b), obteni- se ha comprobado que la rotura no se produce de forma si-
dos mediante MEL, tienen factores de seguridad similares, multánea en todas las columnas, sino de una en una como
que son 3 veces más altos que los obtenidos por métodos fichas de dominó. Todos estos tipos de roturas solo pueden
numéricos (MEF y MDF). En el caso de las columnas con estudiarse mediante métodos numéricos.
capa ficticia (1c), modelizada en MEL, la sobrestimación Por lo tanto, los métodos numéricos resultan más idó-
disminuye hasta casi la mitad. neos para el análisis de estabilidad de terrenos mejorados
En el caso 2, el resultado obtenido por métodos numé- con soil mixing considerando que permite detectar roturas
ricos (MEF y MDF) es 1,5 veces menor que el obtenido por tracción, roturas por flexión, por cortante, deforma-
mediante el método de equilibrio límite sin capa ficticia ciones, desplazamientos, resistencias movilizadas y roturas
(2a). Mientras que usando capa ficticia en las columnas de por cizalla en muros de columnas solapadas.
este mismo modelo (2b) los valores obtenidos son simila- Es importante tener en cuenta que de utilizarse el equi-
res. librio límite para el análisis de estabilidad se aplique el mé-
todo de los perfiles junto con la técnica de la capa ficticia
6. CONCLUSIONES para reducir la sobrestimación de factores de seguridad.
Se ha demostrado que, para la misma relación de área
Existen dos métodos de diseño para terraplenes cimen- mejorada, el uso de paneles de columnas solapadas bajo los
tados sobre terreno mejorado con columnas de suelo-ce- derrames de los terraplenes es mucho más recomendable
mento, ejecutadas por el método de Deep Soil Mixing, el que el de columnas aisladas, ya que, disminuye drástica-
método japonés y el método escandinavo. Ambos métodos mente la posibilidad de rotura de las columnas por vuelco
consideran el cuerpo de terreno mejorado como un mate- o flexión y, por tanto, aumentan considerablemente el fac-
rial equivalente suelo-columnas. Para la estabilidad exter- tor de seguridad.
na se estudian la rotura por deslizamiento horizontal del
bloque de suelo tratado y para la estabilidad interna se ana- 7. BIBLIOGRAFÍA
liza la rotura por una superficie circular de rotura a través
de la zona tratada mediante el método de equilibrio límite. Barksdale, R.D. y Bachus, R.C. (1983). “Design and construc-
Diversos autores han demostrado mediante los ensayos tion of stone columns”. Volume I. US Department of Transporta-
de centrifuga que las roturas más propensas a ocurrir son por tion. Publicación: FHWA-RD-83-026. 194 p.

8 | Ingeniería Civil 175/2014

30 de 35

Estabilidad terraplenes.indd 8 09/09/14 11:43

 Estabilidad de terraplenes...

Broms, B.B. (2004). “Lime and lime/cement columns”. En: Kitazume, M. y Maruyama, K. (2006). “External stability of
Ground Improvement. Editores M.P. Moseley y K. Kirsch 2ª Edi- group column type deep mixing improved ground under em-
ción. Abigdon (Reino Unido): Spon Press, 431 p. bankment loading”. En: Soils and foundations: Japanese Geotech-
CDIT (Coastal Development Institute of Technology), Japan. nical Society.Vol. 46, No 3. pp. 323-340.
(2000). “The deep mixing method. Principle, design and cons- Kitazume, M. y Maruyama, K. (2007). “Internal stability of
truction”. Lisse (Holanda): A.A. Balkema Publishers, 123 p. group column type deep mixing improved ground under em-
EUROSOILSTAB. (2002). “Design guide soft soil stabilisa- bankment loading”. En: Soils and foundations: Japanese Geotech-
tion. Development of design and construction methods to sta- nical Society. Vol. 47, No 3. pp. 437-455.
bilise soft organic soils”. Ministry of Transport Public Works and Kitazume, M. (2011). “Effect of surface improvement layer
Water Management (Holanda). CT97-0351. Project No.: BE 96- on internal stability of group column type deep mixing improved
3177: 94 p. ground under embankment loading”. En: Report of the Port and
Filz, G.M. y Navin, M.P. (2006). “Stability of column-suppor- Airport Research Institute (PARI). Vol. 50. No 1. 20 p.
ted embankments”. Virginia Transportation Research Council. Navin, M.P. (2002). “Stability of embankments founded on
VTRC 06-CR13. 76 p. soft soil improved with deep-mixing-method columns”. Tesis
Filz, G.M. (2007). “Load transfer, settlement, and stability of doctoral (Doctor of Philosophy in Civil and Environmental Engi-
embankments founded on columns installed by deep mixing me- neering). Blacksburg, VA, EE.UU. Virginia Polytechnic Institute
thods”. Geotechnical Engineering Seminar Presentation. Natio- and State University. 247 p.
nal Technical University of Athens. School of Civil Engineering. Navin, M.P. y Filz, G.M. (2006). “Reliability of deep mixing
Geotechnical Department - Foundation Engineering Laboratory. method columns for embankment support”. GeoCongress 2006:
45 p. Geotechnical Engineering in the Information Technology Age.
Kitazume, M., Yamamoto, M. y Udaka, Y. (1999). “Vertical Atlanta, GA (EE.UU). ASCE Proceedings 187, 251. 6 p.
bearing capacity of column type DSM group with low improve- Porbaha, A. (2000). “State of the art in deep mixing tech-
ment ratio”. En: Dry mix methods for deep stabilization. Rotter- nology. Part IV: design considerations”. Proceedings of the ICE
dam (Holanda): Balkema. pp. 115-121. - Ground Improvement, Vol. 4, Issue 3: Enero 2000. pp. 111
Kitazume, M., Okano, K. y Miyagima, S. (2000). “Centrifuge –125.
model tests on failure envelope of column type deep mixing me- Topolnicki, M. (2004). “In situ soil mixing”. En: Ground Im-
thod improved ground”. En: Soils and foundations: Japanese Geo- provement. Editores M.P. Moseley y K. Kirsch 2ª Edición. Abig-
technical Society. Vol. 40, No 4. pp. 43-55. don (Reino Unido): Spon Press. 431 p.

Ingeniería Civil 175/2014 | 9

31 de 35

Estabilidad terraplenes.indd 9 09/09/14 11:43