Capítulo 4

Análisis de estabilidad

48

42
Factor de Seguridad = F = 2.44 Centro de giro
36

30
Elevación (m)

24
Grieta de Tensión
18 2.1 m
12 1
4 4
1 2 1
6
3
4 5 6
0
Arena
-6 Roca Fundación de Arcilla
-12
60 48 36 24 12 0 12 24 36 48 60
Distancia en metros desde eje , X


Figura 4.1 Ejemplo de un análisis de estabilidad de taludes (U. S. Corps of Engineeers, 2003).


              
La modelación matemática de los taludes es parte • Investigar los mecanismos potenciales de

de la práctica de la ingeniería geotécnica, con el falla (analizar cómo ocurre la falla).
objeto de analizar las condiciones de estabilidad de

los taludes naturales y la seguridad y funcionalidad • Determinar la sensitividad o susceptibilidad
del diseño en los taludes artificiales (Figura 4.1). de los taludes a diferentes mecanismos de

activación (Efecto de las lluvias, sismos,
Existe una gran cantidad de metodologías
 etc.).
para la modelación matemática, la cual depende
del objetivo del análisis y de los resultados que
 • Comparar la efectividad de diferentes
se deseen obtener. Los objetivos principales
del análisis matemático de los taludes son los opciones de remediación o estabilización, y su
 efecto sobre la estabilidad del talud.
siguientes:
 • Diseñar los taludes óptimos en término de
• Determinar las condiciones de estabilidad del
talud (si es estable o inestable, y el margen de seguridad, confiabilidad y economía.

estabilidad).



128 DESLIZAMIENTOS - ANALISIS GEOTECNICO

Herramientas disponibles La mayoría de los análisis de estabilidad se realizan

Para el análisis de estabilidad de taludes se utilizando programas comerciales de “software”,
dispone de varias herramientas así: los cuales permiten analizar taludes complejos
o con cantidad significativa de información, en
Tablas o ábacos forma eficiente.
Se han elaborado tablas y ábacos para calcular en
forma rápida y sencilla los factores de seguridad Se recomienda en lo posible utilizar siempre
para una variedad de condiciones. programas de computador.

Análisis gráficos Metodologías para análisis de

Históricamente se han utilizado procedimientos estabilidad
gráficos o de polígonos de fuerzas para calcular las Dentro de las metodologías disponibles se
condiciones de estabilidad de los taludes. Estos
encuentran los métodos de límite de equilibrio, los
sistemas gráficos son poco usados actualmente.
métodos numéricos y los métodos dinámicos para
análisis de caídos de roca y flujos, entre otros.
Cálculos manuales
La mayoría de métodos de análisis se desarrollaron
para cálculos matemáticos manuales o con Los métodos numéricos son la técnica que
calculadora, de acuerdo a fórmulas simplificadas. muestra la mejor aproximación al detalle de las
condiciones de estabilidad en la mayoría de casos
Hojas de cálculo de evaluación de estabilidad de taludes.
Algunos autores han desarrollado hojas de cálculo,
las cuales pueden utilizarse para análisis de Sin embargo, los métodos de límite de equilibrio,
taludes sencillos o con bajo nivel de complejidad. son más sencillos de utilizar y permiten analizar
los casos de falla traslacional y de falla rotacional,
Uso de “Software” así como fallas de inclinación (“Toppling”) y fallas
La técnica de análisis que se escoja depende en cuña.
de las características de los sitios y del modo
potencial de falla, dando especial consideración a Igualmente, los métodos de límite de equilibrio
las fortalezas, las debilidades y las limitaciones permiten análisis combinado con técnicas
de cada metodología de análisis. Hasta el año probabilísticas (Stead y otros, 2000).
1975 la mayoría de los análisis de estabilidad
se realizaban en forma gráfica o utilizando En el caso de sistemas de falla complejos, es
calculadoras manuales. conveniente utilizar metodologías de modelación
que tengan en cuenta los factores que producen
Con la llegada del computador los análisis los movimientos.
se pudieron realizar en forma más detallada
inicialmente utilizando tarjetas FORTRAN y
Los factores que generan el deslizamiento
recientemente con programas de software, los
pueden ser complejos y muy difíciles de
cuales cada día son más poderosos.
modelar; Sin embargo, con el objeto de analizar
esas situaciones complejas, existen algunas
Teniendo en cuenta la gran cantidad de
aplicaciones numéricas disponibles en la herramientas utilizando elementos finitos,
actualidad, es esencial que el ingeniero entienda diferencias finitas, elementos discretos y modelos
las fortalezas y limitaciones inherentes a cada dinámicos.
metodología. Existen una gran cantidad de
herramientas informáticas para el análisis Igualmente, se pueden integrar al análisis
de estabilidad de taludes. Dentro de estas modelaciones de hidrogeología y las solicitaciones
herramientas los métodos de equilibrio límite sísmicas.
son los más utilizados, sin embargo, los métodos
esfuerzo - deformación utilizando elementos En la tabla 4.1 se presenta un resumen de las
finitos han adquirido gran importancia y uso en metodologías utilizadas en análisis convencionales
los últimos años. de estabilidad de taludes.
 ANALISIS DE ESTABILIDAD 129

Tabla 4.1 Metodologías utilizadas en la modelación de taludes

Parámetros
Método Ventajas Limitaciones
utilizados

Existe una gran cantidad de

Genera un número único de
paquetes de software.
factor de seguridad sin tener
Topografía del talud, Se obtiene un número de
en cuenta el mecanismo de
estratigrafía, ángulo factor de seguridad. Analiza
Límite de inestabilidad. El resultado
de fricción, cohesión, superficies curvas, rectas, cuñas,
equilibrio difiere de acuerdo al
peso unitario, inclinaciones, etc. Análisis en
método que se utilice. No
niveles freáticos y dos y tres dimensiones con
incluye análisis de las
cargas externas. muchos materiales, refuerzos y
deformaciones.
condiciones de nivel de agua.

Geometría del Permite simular procesos Es complejo y no lineal.

talud, propiedades de deformación. Permite Comúnmente no se tiene
de los materiales, determinar la deformación del conocimiento de los valores
Esfuerzo-
propiedades talud y el proceso de falla. reales a utilizar en la
deformación
elásticas, elasto- Existen programas para trabajar modelación. Se presentan
continuos
plásticas y de “creep”. en dos y tres dimensiones. varios grados de libertad.
Niveles freáticos, Se puede incluir análisis No permite modelar roca
resistencia. dinámico y análisis de “creep”. muy fracturada.

Geometría del talud, Existe poca información

Discontinuos
propiedades del disponible sobre las
Esfuerzo- Permite analizar la deformación
material, rigidez, propiedades de las juntas.
deformación y movimiento relativo de
discontinuidades Se presentan problemas
elementos bloques.
resistencia y niveles de escala especialmente en
discretos
freáticos. taludes en roca.

Es relativamente fácil de Útiles para diseño

Geometría y
Cinemáticos utilizar. Permite la identificación preliminar. Se requiere
características de las
estereográficos y análisis de bloques críticos, criterio de ingeniería para
discontinuidades.
para taludes utilizando teoría de bloques. determinar cuáles son las
Resistencia a las
en roca Pueden combinarse con técnicas discontinuidades críticas.
discontinuidades.
estadísticas. Evalúa las juntas.

Geometría del talud,

Permite analizar la dinámica de Existe muy poca experiencia
Dinámica de tamaño y forma de los
los bloques y existen programas de su uso en los países
caídos de roca bloques y coeficiente
en dos y tres dimensiones. tropicales.
de restitución.

Relieve del terreno. Se puede predecir el Se requiere calibrar

Dinámica de Concentración comportamiento, velocidades, los modelos para los
flujos de sedimentos, distancia de recorrido y materiales de cada región.
viscosidad y sedimentación de los flujos. Los resultados varían
propiedades de la de acuerdo al modelo
mezcla suelo-agua. utilizado.
130 DESLIZAMIENTOS - ANALISIS GEOTECNICO

CARACTERISTICAS DEL ANALISIS La mayoría de los sistemas de análisis asumen un

DE LIMITE DE EQUILIBRIO criterio de “límite de equilibrio” donde el criterio
de falla de Coulomb es satisfecho a lo largo de una
determinada superficie. Se estudia un cuerpo libre
Un análisis de límite de equilibrio permite obtener
en equilibrio, partiendo de las fuerzas actuantes y
un factor de seguridad, o a través de un análisis
de las fuerzas resistentes que se requieren para
regresivo obtener los valores de la resistencia al
producir el equilibrio. Calculada esta fuerza
cortante en el momento de la falla. Una vez se
resistente, se compara con la disponible del suelo
han determinado las propiedades de resistencia
o roca y se obtiene una indicación del factor de
al cortante de los suelos, las presiones de poros y
seguridad.
otras propiedades del suelo y del talud, se puede
proceder a calcular el factor de seguridad del talud. Otro criterio es el de dividir la masa a estudiar
Este análisis de estabilidad consiste en determinar en una serie de tajadas, dovelas o bloques y
si existe suficiente resistencia en los suelos del considerar el equilibrio de cada tajada por
talud para resistir los esfuerzos de cortante que separado. Una vez realizado el análisis de cada
tienden a causar la falla o deslizamiento. tajada se analizan las condiciones de equilibrio de
la sumatoria de fuerzas o de momentos.
La mayoría de los métodos de límite de equilibrio
tienen en común la comparación de las fuerzas
o momentos resistentes y actuantes sobre una

∑ 
determinada superficie de falla. Las variaciones
principales de los diversos métodos son el tipo
∑ 
de superficie de falla y la forma cómo actúan las
fuerzas internamente sobre la superficie de falla.

∑ 
Concepto de factor de seguridad (F. S.) ∑ 
El factor de seguridad es empleado por los
Ingenieros para conocer cuál es el factor de
amenaza de que el talud falle en las peores Concepto de superficie de falla
condiciones de comportamiento para el cual se El término superficie de falla se utiliza para
diseña. Fellenius (1922) presentó el factor de referirse a una superficie asumida a lo largo de
seguridad como la relación entre la resistencia al la cual puede ocurrir el deslizamiento o rotura del
corte real, calculada del material en el talud y los talud (Figura 4.2); Sin embargo, este deslizamiento
esfuerzos de corte críticos que tratan de producir o rotura no ocurre a lo largo de esas superficies si el
la falla, a lo largo de una superficie supuesta de talud es diseñado adecuadamente. En los métodos
posible falla: de límite de equilibrio el factor de seguridad se
asume que es igual para todos los puntos a lo largo
 de la superficie de falla, por lo tanto este valor
 representa un promedio del valor total en toda la

superficie de falla. Si la falla ocurre, los esfuerzos
En superficies circulares donde existe un centro de de cortante serían iguales en todos los puntos a
giro y momentos resistentes y actuantes: todo lo largo de la superficie de falla.

 Generalmente se asume un gran número de

 superficies de falla para encontrar la superficie de
 falla con el valor mínimo de factor de seguridad,
la cual se denomina “superficie crítica de falla”.
Esta superficie crítica de falla es la superficie más
Existen además, otros sistemas de plantear el factor
probable para que se produzca el deslizamiento;
de seguridad, tales como la relación de altura crítica
Sin embargo, pueden existir otras superficies
y altura real del talud y métodos probabilísticos,
de falla con factores de seguridad ligeramente
así como tablas empíricas locales basadas en el
mayores, los cuales también se requiere tener en
comportamiento típico de los taludes. cuenta para el análisis.



ANALISIS DE ESTABILIDAD 131

La profundidad de las grietas de tensión puede
 determinarse de acuerdo a la siguiente expresión:
Superficie de falla


2c  1 
Zc =  2  45 + φ 
 γ  2 

 Donde:
zc = Profundidad de la grieta de tensión.
 c = cohesión.
γ = Peso unitario del suelo.

Figura 4.2 Superficie de falla y dirección de la resistencia f = Angulo de fricción.
al cortante (U. S. Corps of Engineeers , 2003).

 
Formas de la superficie de falla
  R
Las técnicas de límite de equilibrio se utilizan
cuando las fallas corresponden a deslizamientos de
 
traslación o de rotación sobre superficies de falla
determinadas (Figura 4.3). Se pueden estudiar
 
superficies planas, circulares, logarítmicas,
parabólicas y combinaciones de ellas. En los
 
últimos años se han desarrollado algunos modelos
de superficies de falla con forma no geométrica.
  a. Circular

Análisis de superficies planas

 
Cuando existen discontinuidades planas en la Cuña
roca o en suelo del talud, se acostumbra realizar Activa
  Bloque
análisis de falla a traslación. Esta técnica asume Central
el deslizamiento traslacional de un cuerpo rígido a
  Cuña
lo largo de un plano o a lo largo de la intersección Pasiva
de dos planos como el caso de la falla en cuña.
 
Análisis de superficies curvas
 
En suelos o rocas blandas las superficies de falla a
b. Cuña
deslizamiento tienden a tener una superficie curva.
 
Estas superficies se les conoce como “círculos de
falla o superficies de falla rotacionales”. En los
 
análisis de estabilidad se debe determinar la
localización de la superficie crítica de falla y el
 
factor de seguridad a lo largo de esta superficie.
 
Las grietas de tensión
La existencia de grietas de tensión aumenta la
 
tendencia de un suelo a fallar (Figura 4.4), la c. General - No circular
longitud de la superficie de falla a lo largo de la cual
 
se genera resistencia es reducida y adicionalmente
la grieta puede llenarse con agua, en el caso de
  Figura 4.3 Formas de la superficie de falla (U. S. Corps
lluvias, pueden generarse presiones de poros
transitorias que afectan la estabilidad del talud. of Engineeers, 2003).



 132 DESLIZAMIENTOS - ANALISIS GEOTECNICO

La presencia de grietas de tensión dificulta Resistencia al cortante

en forma considerable la confiabilidad de los La resistencia al cortante para utilizar en los
análisis cuando no se tiene en cuenta este factor. análisis puede ser medida por alguno de los métodos
Las grietas de tensión son muy importantes y de laboratorio o de campo que se indicaron en el
profundas en cortes de taludes, donde existe un capítulo 3. Se debe tener en cuenta si se trata de
alivio de presiones de confinamiento al ejecutarse condiciones drenadas o no drenadas, o si el análisis
la excavación. es en estado no-saturado. Los parámetros deben
corresponder a los niveles de esfuerzos sobre las
Parámetros utilizados en los análisis de superficies de falla potenciales. En los casos en los
cuales ya ha ocurrido la falla del talud, se recomienda
límite de equilibrio
utilizar las resistencias residuales (Skempton, 1970,
Los modelos tienen en cuenta los factores primarios
1977,1985). Igualmente, debe tenerse en cuenta la
que afectan la estabilidad. Estos factores incluyen
disminución de resistencia con el tiempo.
geometría del talud, parámetros geológicos,
presencia de grietas de tensión, cargas dinámicas Para suelos que son completamente saturados,
por acción de sismos, flujo de agua, propiedades el ángulo de fricción para condiciones no drenadas
de resistencia y peso unitario de los suelos, etc. es igual a cero. La resistencia no drenada para
Sin embargo, no todos los factores que afectan suelos saturados puede ser determinada de
la estabilidad de un talud se pueden cuantificar ensayos no-consolidados no-drenados.
para incluirlos en un modelo matemático de límite
de equilibrio. Por lo tanto, hay situaciones en Para suelos parcialmente saturados tales
las cuales un enfoque de límite de equilibrio no como arcillas compactadas o suelos arcillosos
produce resultados satisfactorios. por encima del nivel freático, las resistencias no
drenadas deben obtenerse utilizando ensayos no-
Pesos unitarios consolidados no-drenados en muestras con el mismo
El peso unitario es tal vez el parámetro más sencillo grado de saturación que el suelo en el campo. La
de medir para el análisis de estabilidad de taludes, envolvente de falla para esos suelos generalmente,
es el que influye menos en el factor de seguridad. es curva y por lo tanto es importante utilizar el
Los pesos unitarios totales son pesos húmedos por mismo rango de presiones de confinamiento en los
encima del nivel freático y saturados por debajo ensayos de laboratorio que en el campo.
del nivel freático. En el caso de que se utilicen
pesos sumergidos, se debe ignorar la presencia Condiciones drenadas o no drenadas
de nivel freático. La densidad saturada se puede Las fallas de los taludes pueden ocurrir en
determinar asumiendo un valor de gravedad condiciones drenadas o no drenadas. Si la
inestabilidad es causada por cambios en la carga,
específica G, el cual se puede suponer igual a 2.68
tal como la remoción de materiales de la parte
para la mayoría de los suelos (Cornforth, 2005).
baja del talud o aumento de las cargas en la parte
superior, en suelos de baja permeabilidad, estos
Grieta de
Tensión pueden no tener tiempo suficiente para drenar
durante el tiempo en el cual ocurre el cambio de
Zc carga. En ese caso se dice que las condiciones son
no drenadas.

Ignore este
Generalmente, los suelos tienen
suelo en los permeabilidades suficientes para disipar las
cálculos de presiones de poros en exceso y se comportan en
estabilidad
condiciones drenadas. Para ratas normales de
carga que equivalen a meses o semanas, suelos con
permeabilidades mayores de 10–4 cm/seg, se pueden
considerar drenadas y suelos con permeabilidades
Figura 4.4 Esquema de una grieta de tensión para menores de 10-7 cm/seg, se consideran no drenadas.
análisis de límite de equilibrio (U. S. Corps of Engineeers, Mientras las permeabilidades intermedias se
consideran parcialmente drenadas.

2003).

 ANALISIS DE ESTABILIDAD 133

Duncan (1996), recomienda que para los taludes en Estabilidad a corto y a largo plazo
los cuales la causa de la falla es el aumento de la En la estabilidad a corto plazo debe tenerse en
presión de poros debida a las lluvias, el problema cuenta que los suelos que no tienen un drenaje
debe analizarse como condición drenada. rápido están sujetos a presiones de poros por
acción de las cargas aplicadas. En la estabilidad
Para determinar las condiciones de drenaje a largo plazo se supone los suelos drenados. La
Duncan (1996) sugiere utilizar la siguiente estabilidad a corto plazo de arcillas normalmente
expresión: consolidadas y limos se recomienda modelar con
Cv t análisis de esfuerzos totales. Aunque se puede
T= realizar el análisis utilizando esfuerzos efectivos,
D2 es muy difícil estimar o medir las presiones de
Donde: poros para la utilización en el análisis.
T = Factor adimensional
Cv = Coeficiente de consolidación Para arcillas sobre-consolidadas el análisis
t = Tiempo de drenaje de estabilidad a corto plazo es prácticamente
D = Longitud del camino de drenaje o distancia imposible de realizar, debido a que la resistencia
de salida del agua al cambio de presiones. del suelo cambia muy rápidamente con el tiempo.
En este caso se recomienda utilizar la experiencia
Si T es mayor de 3 la condición es drenada. local en la formación arcillosa específica analizada
Si T es menor de 0.01 la condición es no y utilizar criterios empíricos (Cornforth, 2005).
drenada.
La estabilidad a largo plazo es más fácil de
Si T está entre 0.01 y 3.0 ocurre drenaje parcial analizar que la estabilidad a corto plazo. Para
durante el tiempo de cambio de cargas. En este todos los casos se recomienda utilizar análisis de
caso deben analizarse ambas condiciones. El caso esfuerzos efectivos.
drenado y el caso no drenado.
Limitaciones de los métodos de límite de
Esfuerzos totales y efectivos equilibrio
Como se estudió en el capitulo anterior los Los análisis de límite de equilibrio tienen algunas
problemas de estabilidad de taludes pueden limitaciones entre las cuales se encuentran las
analizarse suponiendo sistemas de esfuerzos siguientes:
totales o efectivos. En principio, siempre es posible
analizar la estabilidad de un talud utilizando el • Se basan solamente en la estática. Como
método de presión efectiva, porque la resistencia los métodos de límite de equilibrio se basan
del suelo es gobernada por las presiones efectivas solamente en la estática y no tienen en cuenta
tanto en la condición drenada, como en la condición las deformaciones, las distribuciones de
no drenada, pero en la práctica; sin embargo es presiones en muchos casos no son realistas.
virtualmente imposible determinar con precisión Sin embargo, debe tenerse en cuenta que estos
cuáles son los excesos de presión de poros que esfuerzos no realistas generalmente ocurren en
se van a generar por los cambios en las cargas algunas tajadas del análisis y no significa que
(excavaciones, colocación de rellenos o cambios en el factor de seguridad general sea inaceptable.
el nivel de agua).
• Suponen los esfuerzos uniformemente
Debido a esta razón, no es posible desarrollar distribuidos. Debe tenerse cuidado cuando
análisis precisos de estabilidad en estas existan concentraciones de esfuerzos debidos
condiciones, utilizando procedimientos de a la forma de la superficie de falla o a la
esfuerzos efectivos. Sin embargo, se puede trabajar interacción suelo-estructura.
todo el análisis utilizando presiones efectivas,
sin que se requiera especificar los valores de los • Utilizan modelos de falla muy sencillos.
excesos de poros en las condiciones no drenadas. El diseño de taludes utilizando solamente
La mayoría de los modelos de análisis trabajan la modelación con métodos de límite de
con base en presiones efectivas. equilibrio es completamente inadecuado si los
 

134 
DESLIZAMIENTOS - ANALISIS GEOTECNICO


procesos de falla son complejos, especialmente Tajada típica
Superficie
freática
cuando están presentes procesos de “creep”,

deformación progresiva, flujo, rotura por θ

fragilidad, licuación y otras formas de

deterioro de la masa del talud. 
hw Cabeza de Presión
de poros (hwCos2 θ )

• Generalmente se asume el material  como

isotrópico. La mayoría de los trabajos
 el tema
que aparecen en la literatura sobre
Linea Equipotencial
a) Superficie Freática
asumen que el suelo es un material isotrópico
y han desarrollado métodos de análisis de Tajada típica
superficies circulares o aproximadamente Superficie
circulares principalmente. Sin embargo, el piezometrica
mecanismo de falla en materiales residuales
donde aparece el suelo, la roca meteorizada y Cabeza de Presión
hw
la roca sana, así como formaciones aluviales y de poros (hwCos ) de poros (hw)

coluviales no-isotrópicas requieren de nuevos

enfoques y del estudio de superficies de falla
no simétricas.
b) Superficie piezometrica
D
A pesar de las debilidades de un determinado θ
modelo, determinar el factor de seguridad B
asumiendo superficies probables de falla, permite
al ingeniero tener una herramienta muy útil para
la toma de decisiones. Los métodos de límite de h1
C h2
equilibrio son una herramienta muy útil en la
práctica y se recomienda tener cuidado de no
A
abusar en la aplicación del método para casos E
complejos donde la distribución de esfuerzos y las Líneas
deformaciones juegan un papel importante en el de Flujo Líneas Equipotenciales
comportamiento del talud (Krahn, 2004). AB- Superficie freática real
CD- Inclinación asumida del
PRESIONES DE POROS nivel freático dentro de la
tajada

Las condiciones de presión de poros son c) Redes de Flujo

generalmente obtenidas de las características de
las aguas subterráneas y pueden especificarse para Figura 4.5 Representación de la presión de poros.
los análisis utilizando los siguientes métodos:

Superficie freática En el caso de líneas freáticas de gran pendiente,

Esta superficie o línea en dos direcciones se define el cálculo anterior puede resultar sobre estimado
como el nivel libre del agua subterránea. En una y se requiere tener en cuenta que las líneas
superficie freática la presión de poros es calculada equipotenciales tienden a ser curvas.
de acuerdo a las condiciones de estado de régimen

permanente (“Steady-state”). Este concepto Datos piezométricos
se basa en la suposición de que todas las líneas Es la especificación de presiones de poros en
equipotenciales sean ortogonales. Entonces, si la puntos discretos dentro del talud y la utilización
inclinación del segmento de superficie freática es θ de un esquema de interpolación para estimar las
y la distancia vertical entre el punto y la superficie presiones de poros requeridas a cualquier punto.
freática es hw, entonces la presión de poros está Las presiones piezométricas pueden determinarse
dada por la expresión ( Figura 4.5): mediante piezómetros, redes de flujo o soluciones
numéricas, utilizando diferencias finitas o
(
u = γ w hw 2 θ) elementos finitos.
 ANALISIS DE ESTABILIDAD 135

Aunque este sistema está disponible solamente equilibrio, debido a que puede generar valores de
en muy pocos de los programas de computador resistencia no confiables (Abramson y otros, 2002).
existentes, se recomienda por su confiabilidad Sin embargo, con los modelos de computador
para representar las condiciones reales en el actualmente disponibles es relativamente sencillo
campo (Chugh, 1981). incorporar las presiones de poros negativas para
tener en cuenta el escenario de la situación no
Relación de presión de poros saturada.
Este es un método muy simple y popular para
normalizar el valor de la presión de poros en un Efecto de los ductos de agua en la corona
talud de acuerdo a la definición: de los taludes
Siempre que sea posible es imperativo el localizar
u los ductos de agua lejos de la corona de taludes o
ru =
σv laderas donde se requiera su estabilidad. Como
una regla general, la distancia entre la corona
Donde: de los taludes y la localización de todo tipo de
u = Presión de poros tuberías y servicios debe ser igual a la altura total
σv = Esfuerzo total vertical del suelo a una del talud. Aunque este es el estándar mínimo
profundidad z. recomendado (Abramson, 1996), en ocasiones se
requieren aislamientos mayores. En el caso en el
Este factor se implementa fácilmente, pero la cual no es posible mantener estos aislamientos,
mayor dificultad está asociada con la asignación el talud debe diseñarse para tener en cuenta su
de este parámetro a diferentes partes del talud. saturación debida a la muy posible infiltración de
En ocasiones, el talud requiere de una extensiva agua, teniendo en cuenta que en gran cantidad de
subdivisión en regiones con diferentes valores de ru. casos se producen fugas de los ductos.

Superficie piezométrica metodos de LIMITE DE

Esta superficie se define para el análisis de una EQUILIBRIO
determinada superficie de falla. Debe tenerse
claridad en que la superficie piezométrica no es la El análisis de los movimientos de los taludes o
superficie freática y que el método de calcular la laderas durante muchos años se ha realizado
presión de poros es diferente para los dos casos. En utilizando las técnicas del límite de equilibrio.
la superficie piezométrica, la presión de poros es la El sistema de límite de equilibrio supone que
distancia vertical entre la superficie piezométrica en el caso de una falla, las fuerzas actuantes y
indicada y el punto. resistentes son iguales a lo largo de la superficie de
falla equivalentes a un factor de seguridad de 1.0.
Presión de poros constante
Este procedimiento puede utilizarse si el Ingeniero El análisis se puede realizar estudiando
desea especificar una presión de poros constante directamente la totalidad de la longitud de la
a una determinada capa de suelo. Este sistema superficie de falla o dividiendo la masa deslizada
puede utilizarse para analizar la estabilidad de en tajadas o dovelas. Cada día se han mejorado los
rellenos colocados sobre suelos blandos, durante la sistemas de dovelas desarrollados por Petterson y
construcción donde se generan presiones de poros, Fellenius (1936). Algunos métodos son precisos
de acuerdo a la teoría de la consolidación. y otros solamente aproximados (Figura 4.6). Los
métodos de Bishop (1955) y Janbú (1954) han sido
Presiones de poros negativas muy utilizados en los últimos 50 años y se han
En algunos casos el ingeniero desea utilizar en desarrollado métodos de análisis más precisos y
los análisis las presiones de poros negativas para complejos como los de Morgenstern y Price (1965)
aprovechar la resistencia adicional o cohesión y Spencer (1967), ayudados por programas de
aparente debida a la succión en suelos no saturados. software, los cuales permiten realizar análisis
Aunque teóricamente la cohesión aparente es una muy rigurosos. Generalmente los métodos son de
realidad física, algunos autores no recomiendan iteración y cada uno de los métodos posee un cierto
su incorporación a los modelos de límite de grado de precisión.
136 DESLIZAMIENTOS - ANALISIS GEOTECNICO

En la tabla 4.2 se enumeran algunos de los métodos más utilizados.

Tabla 4.2 Métodos de análisis de estabilidad de taludes

Superficies
Método Equilibrio Características
de falla

Bloque delgado con nivel freático, falla paralela a

Talud infinito Rectas Fuerzas
la superficie.

Cuñas con
Cuñas simples, dobles o triples analizando las
Bloques o cuñas tramos Fuerzas
fuerzas que actúan sobre cada cuña.
rectos

Espiral logarítmica Espiral Fuerzas y Superficie de falla en espiral logaritmica. El radio

(Frohlich, 1953) logarítmica momentos de la espiral varía con el ángulo de rotación.

Círculo de falla, el cual se analiza como un solo

Arco circular,
Circulares Momentos bloque. Se requiere que el suelo sea cohesivo
(Fellenius, 1922)
(φ = 0).

Ordinario o
de Fellenius Circulares Fuerzas No tiene en cuenta las fuerzas entre dovelas.
(Fellenius 1927)

Bishop simplificado Asume que todas las fuerzas de cortante entre

Circulares Momentos
(Bishop 1955) dovelas son cero.

Janbú Simplificado Cualquier Asume que no hay fuerza de cortante entre

Fuerzas
(Janbú 1968) forma dovelas.

Sueco Modificado.
Cualquier Las fuerzas entre dovelas tienen la misma dirección
U.S. Army Corps of Fuerzas
forma que la superficie del terreno.
Engineers (1970)

Las fuerzas entre dovelas están inclinadas a un

Lowe y Karafiath Cualquier
Fuerzas ángulo igual al promedio de la superficie del terreno
(1960) forma
y las bases de las dovelas.

Cualquier Momentos La inclinación de las fuerzas laterales son las

Spencer (1967)
forma y fuerzas mismas para cada tajada, pero son desconocidas.

Morgenstern y Cualquier Momentos Las fuerzas entre dovelas sea asume que varían de
Price (1965) forma y fuerzas acuerdo a una función arbitraria.

Utiliza el método de las dovelas para calcular la

Cualquier Momentos
Sarma (1973) magnitud de un coeficiente sísmico requerido para
forma y fuerzas
producir la falla.
 ANALISIS DE ESTABILIDAD 137

Métodos
de
Cálculo

Métodos Métodos
de numéricos
Equilibrio Límite

Exactos
Rotura plana Aproximados
Rotura por cuña

Elementos Finitos Diferencias Finitas

Cuña Simple Cuña Doble Tabla de Taylor Tabla de Janbú

Elementos Elementos
Cuña Triple No Exactos
Discretos de Borde

Métodos Métodos de
de Dovelas
estabilidad global

Aproximados Precisos
Espiral Morgenstern-Price,
Arco Circular Janbú, Fellenius,
Logaritmica Spencer,
Bishop simplificado
Bishop riguroso


Figura 4.6 Métodos de análisis de estabilidad de taludes.



TABLAS PARA ANALISIS RAPIDOS
 El uso de tablas no debe reemplazar los análisis
Para taludes simples homogéneos se han
rigurosos, sino que puede servir de base de
desarrollado tablas que permiten un cálculo rápido
 comparación de los resultados, o para la
del factor de seguridad. Existe una gran cantidad
evaluación rápida y general de las condiciones de
de tablas desarrolladas por diferentes a utores.
 estabilidad.
La primera de ellas fue desarrollada por
Las tablas dan una “idea” general del nivel
Taylor en 1966. Desde entonces varias tablas
de estabilidad de un talud. Las tablas de mayor
han sido sucesivamente presentadas por Bishop
utilidad son las que se elaboran para áreas
y Morgenstern (1960), Hunter y Schuster (1968),
homogéneas específicas locales con base en análisis
Janbú (1968), Morgenstern (1963), Spencer
completos de estabilidad y debidamente validadas
(1967), Terzaghi y Peck (1967) y otros, las cuales
en campo.
se resumen en la tabla 4.3.
138 DESLIZAMIENTOS - ANALISIS GEOTECNICO

Tabla 4.3 Listado de tablas para cálculo de estabilidad de taludes disponibles en la literatura.

Método
Inclinación
Autor Parámetros analítico Observaciones
de talud
utilizado
φ=0
cu 0-90o Análisis no drenado.
Taylor (1966) Circulo de
c, φ 0-90 o Taludes secos solamente.
fricción
Bishop y
Primero en incluir efectos del
Morgenstern c, φ,ru 11-26.5 o Bishop
agua.
(1960)

Análisis no drenado con cero

Gibsson y resistencia en la superficie y
cu 0-90 o
φ=0
Morgenstern cu aumenta linealmente con la
profundidad.

Spencer (1967) c, φ, ru 0-34 o Spencer Círculos de pie solamente.

Una serie de tablas para

cu φ=0 diferentes efectos de
Janbú (1968) 0-90 o
c, φ, ru Janbú GPS movimiento de agua y grietas
de tensión.

Análisis no drenado con

una resistencia inicial en
Hunter y Schuster
cu 0-90 o φ=0 la superficie y cu aumenta
(1968)
linealmente con la
profundidad.

Chen y Giger
c, φ 20-90 o Análisis límite
(1971)

Bishop y Morgenstern (1960)

O´Connor y
c, φ,ru 11-26 o Bishop extendido para incluir
Mitchell (1977)
Nc = 0.1
Incluye agua subterránea y
Círculo de
Hoek y Bray c, φ 0-90 o grietas de tensión.
fricción
(1977) c, φ 0-90 o Análisis de bloque en tres
Cuña
dimensiones.

Círculo de Extensión del método de

Cousins (1978) c, φ 0-45 o
fricción Taylor (1966).

Charles y Soares Envolvente de falla no lineal

φ 26-63 o Bishop
(1984) de Mohr-Coulomb.

Extensión de Bishop y
Morgenstern (1960) para un
Barnes (1991) c, φ, ru 11-63 o
Bishop
rango mayor de ángulos del
talud.
 ANALISIS DE ESTABILIDAD 139

TABLA DE TAYLOR Donde:

No = Número de estabilidad que se obtiene de
Una forma rápida para determinar el factor de la tabla
seguridad de un talud es utilizando las tablas Creq = Cohesión requerida para F.S. = 1.0
de Taylor. Es importante tener en cuenta que el γ = Peso unitario del suelo
método de Taylor supone un suelo homogéneo y un H = Altura del talud
manto rígido profundo. Este método solo se utiliza
para suelos cohesivos (φ =0), y se aplica solamente
para el análisis de esfuerzos totales, debido a que Paso 5. Calcular el Factor de seguridad del
no considera presiones de poros. talud

A continuación se presenta el procedimiento de Como paso final se calcula el factor de seguridad

manejo de la tabla de Taylor. con la siguiente fórmula:

Paso 1. Parámetros que se requieren para el Cu

F S  =
análisis. Se requiere conocer: Creq
• Altura del talud H (metros)
• Cohesión del suelo Cu (KN/m2) TABLAS DE JANBÚ
• Pendiente del talud β (grados)
• Peso especifico del suelo γ (KN/m3)
Las tablas desarrolladas por Janbú (1968),
• Profundidad al manto de suelo duro permiten el análisis de diferentes condiciones
impenetrable D (Metros) geotécnicas y factores de sobrecarga en la corona
del talud, incluyendo niveles freáticos y grietas de
Paso 2. Calcular el factor de profundidad d tensión.
El factor de profundidad, d, se calcula por medio
de la fórmula: El método de tablas de Janbú presenta
D dos procedimientos, uno para suelos cohesivos
d= (φ = 0), y otro para suelos friccionantes (φ > 0).
H Para suelos cohesivos el procedimiento es el mismo
Donde: de Taylor. Para suelos friccionantes o mixtos el
D = profundidad del manto de suelo duro procedimiento es un poco más complejo.
impenetrable (Roca).
H = altura del talud. Procedimiento para las Tablas de Janbú
para φ = 0.
Paso 3. Determinar el número de estabilidad (No)
Paso 1. Parámetros que se requieren para el
Del gráfico de Taylor (Figura 4.7) se determina análisis
el valor del número de estabilidad, No, el cual Se requiere conocer:
depende del ángulo del talud, β, y del valor de “d”
que se calculó en el paso anterior. • Altura de cada suelo H (metros).

• Pendiente del talud β (grados).

Paso 4. Calcular Creq para factor de seguridad de
1.0. Se utiliza la siguiente expresión: • Cohesión del suelo Cu (KN/m2).

• Altura del nivel freático HW (m).

γ H • Peso especifico del suelo γ (KN/m3)

NO =
Creq
• Perfil geotécnico incluyendo todos los mantos
de suelo.
 
140 DESLIZAMIENTOS - ANALISIS GEOTECNICO


 11


10 Factor de seguridad
Círculos pie Círculos
 D
Círculos base d = H Talud
Número de estabilidad, No 9 Círculos Talud

β H
8
Base Firme D

0
=
1
0.
d

0.3 .2
γ = Peso unitario total del suelo

0
7 Círculos
base

5
1.0
0.

1.5
2
3
6

5.53 d=α
pie
5 los
í rcu
C

4
3.83 Cotg β

0.25 0.50 0.75 1.0 1.5 2 3 4 6 10 α

0
90 80 70 60 50 40 30 20 10 0
Angulo del Talud - β (grados)

Figura 4.7 Tabla de Taylor (Taylor, 1966).

Profundidad al manto de suelo duro

• Paso 4. Calcular C promedio
impenetrable D (Metros) Utilizando como guía el círculo estimado, se
determina el valor promedio de la resistencia, C.
Paso 2. Calcular el factor d Esto se realiza calculando el promedio ponderado
Calcular el factor de profundidad, d, por medio de de las resistencias a lo largo del arco de falla, con
el número de grados interceptado por cada tipo de
Hw suelo como factor de ponderación.
d=
H
la fórmula: Paso 5. Calcular el factor de reducción
Puede encontrarse factor de reducción por carga
Donde:
adicional, factor de reducción por sumergencia
HW= Altura del nivel freático
e infiltración, factor de reducción por grieta de
H = Profundidad del pie del talud al punto más tracción sin presión hidrostática en la grieta y
bajo del círculo de falla. factor de reducción por grieta de tracción con
presión hidrostática en la grieta. En las figuras
Paso 3. Obtener la localización del círculo crítico 4.10 a 4.13, se muestran las tablas a usar según el
(Xo, Yo). (Figura 4.8) caso que se presente.

De las Figuras 4.8 y 4.9 determinar la localización Paso 6. Calcular Pd

del centro del círculo crítico Xo, Yo. Para taludes Pd se calcula con la siguiente fórmula:
más empinados que 53°, el círculo crítico pasa por
el pié. Para taludes más tendidos que 53°, el círculo
Pd =
(γ  H ) + q − (γ w  H w )
crítico pasa tangente a la superficie firme o roca. µ q  µ w  µt
 ANALISIS DE ESTABILIDAD 141

4
β = 0º
1.0
Xo
3 Centro Crítico
Yo 30º
β H 0.9
Abscisa del centro - Xo

Factor µq
2 60º
Xo = xoH d = 0.5
d=0 0.8
1 Círculo por el pie 90º

0 0.1 0.2 0.3 0.4 0.5

Pie
0 r culo Relación q/ γH
Cí ase (a)
yB
Cot β
-1 0.25 0.50 1.0 1.5 2 3 4 6 10 α

90 80 70 60 50 40 30 20 10 0

Angulo del Talud - β (grados)

d=α
1.0
1.0

Figura 4.8 Coordenada Xo para el círculo crítico. 0.5

Factor µq

0.9
(Janbú 1968).
0

0.8
Círculo por la base

5 0 0.1 0.2 0.3 0.4 0.5

(b) Relación q/ γH
4

Yo = yo H
Ordenada del centro - yo

3 d = 3.0
2.5
Círcu dio

2.0 q
0.3 Leyenda
me
lo pu

2 1.5
Círc
ulo 1.0
P
nto

ie
0 H
1
β
Cot β
D=dH
0.25 0.50 1.0 1.5 2 3 4 6 10 α Base Firme
0
90 80 70 60 50 40 30 20 10 0
Angulo del Talud - β (grados)


Figura 4.9 Coordenada Yo para el círculo crítico. Figura 4.10 Factor de reducción por carga adicional
(Janbú 1968). para tablas de Janbú.

142 DESLIZAMIENTOS - ANALISIS GEOTECNICO

β = 0º
1.0
30º

Factor µw y µ'w
60º
0.9
90º

Círculo por pie H

0.8 Hw
β

0 0.5 1.0 D= dH
Base Firme
(a) Relación Hw / H y H'w / H

d=α
1.0
1.0
Factor µw y µ'w

0.5 H
0.9 Hw
0
D= dH
Base Firme
0.8 Círculo por la base

0 0.5 1.0

(b) Relación Hw / H y H'w / H

Figura 4.11 Factor de reducción por sumergencia (μw) e infiltración (μ’w).

β = 0º
1.0
30º
0.9
 
60º
0.8
Factor µt

90º
0.7

0.6 Grietas de
Círculo por pie Tracción
0.5
0 0.1 0.2 0.3 0.4 0.5 Ht

(a) Relación Ht / H H

β
d=α
1.0 D= dH Base Firme
1.0
0.5
0.9
0
0.8
Factor µt

0.7

0.6 Círculo por la base

0.5
0 0.1 0.2 0.3 0.4 0.5

(b) Relación Ht / H

Figura 4.12 Factor de reducción por grieta de tracción sin presión hidrostática en la grieta. (Janbú, 1968).
 ANALISIS DE ESTABILIDAD 143

Donde:
γ = peso unitario promedio del suelo β = 0º
H = altura del talud
q = sobrecarga 1.0
30º
γw = peso unitario del agua 0.9
Hw = altura de agua fuera del talud 60º

Factor µt
μq = factor de reducción por sobrecarga 0.8

μw = factor de reducción por sumergencia 0.7 90º

μt = factor de reducción por grieta de tensión
0.6
Círculo por el pie
Si no hay sobrecarga, μq = 1; si no hay sumergencia, 0.5
μw = 1, y si no hay grieta de tensión, μt= 1. 0 0.1 0.2 0.3 0.4 0.5

(a) Relación Ht / H
En la fórmula de Pd se toma q = 0, μq =1 para
condición no consolidada

Paso 7. Calcular el número de estabilidad NO d=α

De la Figura 4.14, se determina el valor del número 1.0 1.0
de estabilidad, No, que depende del ángulo del 0.5
talud. 0.9

Factor µt 0.8
Paso 8. Calcular la cohesión requerida 0
Se calcula despejando creq de la fórmula del número 0.7

de estabilidad No.
0.6
Círculo por la base

γ H 0.5
0.1
No = 0 0.2 0.3 0.4 0.5
Creq (b) Relación Ht / H

Paso 9. Calcular el factor de seguridad

Se utiliza la expresión:

N o  Creq Grietas de
F S  = Tracción
Pd
Ht

Procedimiento para las Tablas de Janbú H

para φ > 0. β

A continuación se describe los pasos a seguir para D= dH

Base Firme
este caso, que es similar al anterior desde el paso
1 hasta el paso 6.

Paso 1. Parámetros que se requieren para el

análisis
Figura 4.13 Factor de reducción por grieta de tracción
Paso 2. Calcular el factor d.
con presión hidrostática en la grieta.
Paso 3. Obtener la localización del círculo crítico.

Paso 4. Calcular C promedio
Paso 5. Calcular el factor de reducción

Paso 6. Calcular Pd

144 DESLIZAMIENTOS - ANALISIS GEOTECNICO

11

10 Factor de seguridad Fn = No c
Círculos pie Pd Círculos
Círculos base Talud
9 Círculos Talud
D
d=
Número de estabilidad, No
H H β
8
Base Firme D

0
=
1
0.
d

0. 2
γ = Peso unitario total del suelo

0.
3
7 Círculos
base

0
0.

1.
5
1.
2
3
6

5.53 d=α
s pie
5 ulo
c
Cír

4
3.83 Cotg β

0.25 0.50 0.75 1.0 1.5 2 3 4 6 10 α

0
90 80 70 60 50 40 30 20 10 0
Angulo del Talud - β (grados)

Figura 4.14 Número de estabilidad.

100
300
Para c = 0 50
200
F = Pe b tg φ 30
Pd 20
100
15
Número Crítico de Estabilidad, Ncf

10

50 8
6
4

20 2

1
10
0

5
Valores de λ cφ

F=N C cf
2 Pd
λcφ = Pe tg φ
c
1
0 1 2 3 4 5
Relación de Talud b = cot β

q
γ H + q - γ w Hw
b Pd =
l µq µw µt
Ht
H
Hw H' w
β γ H + q - γ w Hw
Pe =
µc µ'w

Figura 4.15 Número de estabilidad Ncf.

 ANALISIS DE ESTABILIDAD 145

Paso 7. Calcular Pe. 3.0

Pe se calcula con la siguiente fórmula:

Pe =
( γ  H ) + q − (γ w  H w ) yo
λ cφ = 100
µq  µ w 2.0 20
10
5

Coordenadas Unitarias Xo e Yo
Donde: 2
100
H´w = altura del agua dentro del talud. 0 20
10
μ´w = factor de reducción por infiltración. 5
xo 2
1.0
λ cφ = 0
Si la sobrecarga se aplica rápidamente de modo
que no hay suficiente tiempo para que los suelos
se consoliden bajo la sobrecarga, se toma q=0 y
μq = 1 en la fórmula de Pe. Si no existe sobrecarga,
μq = 1, y si no existe infiltración, μ’w =1. 0

Coordenadas
Xo = xo H
Yo = yo H
Paso 8. Calcular el parámetro a dimensional λCφ
. Este parámetro es calculado con la siguiente
-1.0
fórmula: 1 2 3 4 5
0
Pe   φ Relación de Talud b
λCφ =
C
Figura 4.16 Coordenadas del centro del círculo
crítico (suelos con φ >0).
Donde: 
tan φ = valor promedio de tan φ.
C = valor promedio de las cohesiones
METODO DEL TALUD INFINITO

Paso 9. Calcular el número de estabilidad Ncf Con frecuencia en deslizamientos de gran

Para calcular este número de estabilidad se usa la magnitud la mayor parte de la masa deslizada
tabla presentada en la Figura 4.15 se mueve en forma aproximadamente paralela
a la superficie del terreno. La naturaleza del
movimiento está controlada por algún elemento
Paso 10. Calcular el factor de seguridad geológico como una capa de roca o una capa de
El factor de seguridad se calcula con la siguiente materiales poco resistentes. Si la longitud relativa
fórmula: del deslizamiento es muy grande con relación a su
espesor, la contribución de resistencia en la cabeza
C y pie del deslizamiento son menores comparados
F S  = N cf 
Pd con la resistencia en el resto de la superficie de
falla.

Paso 11. Obtener la localización del círculo crítico. En las condiciones indicadas se presenta una falla
Para obtener las coordenadas del círculo crítico se paralela a la superficie del talud, a una profundidad
somera y la longitud de la falla es larga comparada
realiza con la tabla mostrada en la Figura 4.16.
con su espesor. Este tipo de deslizamiento se
Se calcula b = cot β Y
puede analizar suponiendo un talud infinito.
146 DESLIZAMIENTOS - ANALISIS GEOTECNICO

El método del talud infinito es un sistema muy

1.0

rápido y sencillo para determinar el factor de 0.9

Relación de presión de poros h/2

seguridad de un talud, suponiendo un talud 0.8 z
h γ =2γ ω
largo con una capa delgada de suelo, en el 0.7 c' = 0, φ β
cual cualquier tamaño de columna de suelo es 0.6

S
S R 1.9
representativo de todo el talud (Figura 4.17). Las

=
SSR= tan φ

2.
0.5

0
1.8
tan β
suposiciones del método de talud infinito son las

1.7
0.4
siguientes: suelo isotrópico y homogéneo, talud

1. 6
1.5
0.3
infinitamente largo, y superficie de falla paralela

1.4
1.3
al talud. El uso principal del método del talud
0.2

1.2
1.1
infinito es para elaborar planos de amenaza a los 0.1

deslizamientos utilizando SIGs. 0

1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0
Factor de seguridad F
Para un talud uniforme y relativamente largo
en el cual el mecanismo de falla esperado no es muy
profundo, los efectos de borde son despreciables y Figura 4.18 Determinación del factor de seguridad
el factor de seguridad puede calcularse para un FS para diferentes alturas del nivel de agua de una
talud infinito de una unidad de área utilizando el 
determinada relación de resistencia para el talud seco
criterio Mohr - Coulomb. (SSR). (Cornforth, 2005).


Analizando el elemento de la figura 4.17 y El ángulo de fricción para factor de seguridad igual
realizando una igualdad de fuerzas resistentes y a 1.0 se le denomina ángulo de reposo.
actuantes se obtiene la siguiente expresión:
Si para el caso anterior el nivel de agua se
c + ( γ z − γ w h ) 2 β  φ  encuentra en la superficie del terreno y por lo
F S  =
γ zsenβ  β tanto el suelo se encuentra totalmente saturado
y la cohesión es cero, se obtiene la siguiente
expresión:
Simplificando para un talud seco de suelos sin
cohesión (c’ = 0) γ   φ 
F S  =
Tanφ γ  β
F S  =
Tanβ Donde:
γ’ = peso unitario sumergido
γ = peso unitario saturado
b
B De la expresión anterior se obtiene que si el suelo
A z se encuentra saturado totalmente, el factor de
h
P
seguridad es aproximadamente la mitad del factor
x
hs de seguridad del talud seco. El factor de seguridad
C disminuye a medida que sube el nivel del agua
E
D
I (Figura 4.18). El factor de seguridad varía con la
β
posición del nivel freático de acuerdo a la relación
PR ru que se denomina coeficiente de presión de poros
W y que relaciona la presión de poros con la altura
PL de suelo.
S u
N U= U I ru =
β
γz
El método del talud infinito también se puede aplicar
Figura 4.17 Diagrama de análisis, método del talud a taludes de suelos cohesivos siempre y cuando la
infinito. (Cornforth, 2005). falla sea paralela a la superficie del talud.
 ANALISIS DE ESTABILIDAD 147

Paso 2. Calcular el factor de seguridad.

Q El factor de seguridad varía con la posición del
nivel freático. El factor de seguridad se determina
por medio de la siguiente expresión:

h c + ( γ z − γ w h ) 2 β  φ 
F S  =
w
γ zsenβ  β
z

Fuerza Resistente
β
Interface
α


Figura 4.19 Talud infinito. a) Cuña Simple

El método del talud infinito cumple condiciones PA Bloque PP

para equilibrio de fuerzas y equilibrios de momentos Analizado

a pesar de que no se considera explícitamente, Superficie Débil

debido a que las fuerzas son colineales y la fuerza
normal actúa en el centro del bloque (Duncan y
Wright, 2005). b) Bloque Deslizante

Este método es muy preciso para analizar suelos

estratificados con falla paralela a la superficie del Graben
terreno.
Cu ñ a
Procedimiento para el método de talud Princip
al
infinito
Zona D
ébil
Paso 1. Parámetros que se requieren para el
análisis. Se requiere conocer: c) Cuña Doble

• Altura de la masa deslizante z (metros).

Graben
• Altura del agua subterránea medida durante Levantamiento
el movimiento h (metros). Cu ñ a
Princip
al

• Angulo de inclinación con la horizontal β Zona

Débil
(grados).

d) Cuña Triple
• Peso especifico del suelo γ (KN/m3).

• Angulo de fricción φ (grados).

Figura 4.20 Tipos de bloques o cuñas para análisis de
• Cohesión C (KN/m ). 2 estabilidad de taludes.


148 DESLIZAMIENTOS - ANALISIS GEOTECNICO

ANALISIS DE BLOQUES O CUÑAS

El análisis de estabilidad de taludes puede

realizarse suponiendo superficies de falla rectas
pre-determinadas. Pueden analizarse superficies
compuestas por una sola línea o por varias líneas,
formando cuñas simples, dobles o triples (Figura
a. Buscar el bloque central crítico
4.20). Este tipo de análisis es apropiado cuando
hay una superficie potencial de falla relativamente
Variar θP para
recta a lo largo de un material relativamente Variar θA para encontrar la fuerza
duro o relativamente blando; por ejemplo mantos encontrar la fuerza mínima en el centro
del bloque
aluviales débiles. Uno de estos métodos es conocido máxima en el centro
del bloque
como “método del bloque deslizante”. αA=θA-φD/2 αP=θP-φD/2

En el análisis de cuñas dobles o triples se

requiere determinar la localización del bloque
central crítico, las inclinaciones críticas de las
cuñas activa y pasiva, y los factores de seguridad
mínimos o críticos.

Los métodos para la localización del bloque

central crítico se muestra en la figura 4.21 (a)
b. Esquema para buscar la inclinación de la cuña
y se refiere a variar en forma sistemática las
coordenadas de los dos extremos de la base
del bloque central hasta encontrar el factor
de seguridad mínimo. Para cada posición del Figura 4.21 Análisis
   de cuñas.
  Suposiciones de
   
localización de cuñas para calcular factores de seguridad
bloque central se varían las inclinaciones de
las (U. S. Army Corps of Engineers, 2003).
cuñas activa y pasiva para encontrar el factor de
seguridad mínimo para cada posición del bloque.
(Figura 4.21 (b)). Una suposición que se efectúa
con frecuencia es la de establecer la inclinación En el caso de tres bloques, la cuña superior se le
de cada cuña activa a un ángulo de 45º + φ’/2, y llama cuña activa y las otras dos, cuña central y
cada cuña pasiva a 45º - φ’/2. Esta suposición solo pasiva, respectivamente. El factor de seguridad
es válida cuando las superficies superiores de puede calcularse sumando las fuerzas horizontales
las cuñas son horizontales, pero puede utilizarse así:
cuando son pendientes suaves. Otra técnica
utilizada es la suposición de cuñas que aumentan
F S  =
(
Pp + c m L + W − u  θ m  )
de inclinación de abajo hacia arriba. P a

Método del Bloque Deslizante

El análisis del bloque deslizante puede utilizarse Donde:
cuando existe a una determinada profundidad Pp = Fuerza pasiva producida por la cuña
una superficie de debilidad relativamente recta y inferior.
delgada sub-horizontal. La masa que se mueve Pa = Fuerza activa producida por la cuña
puede dividirse en dos o más bloques y el equilibrio superior.
de cada bloque se considera independientemente, c’m = Cohesión efectiva del suelo blando en la
utilizando las fuerzas entre bloques (Figura 4.22). base del bloque central.
L = Longitud del fondo del bloque central.
No considera la deformación de los bloques W = Peso total del bloque central.
y es útil cuando existe un manto débil o cuando u = Fuerza total de poros en el fondo del bloque
aparece un manto muy duro sobre el cual se puede central.
presentar el deslizamiento. θm = Fricción del suelo en el fondo del bloque.
 ANALISIS DE ESTABILIDAD 149

El factor de seguridad se determina por medio de

Cuña Activa Bloque Central Cuña Pasiva
la expresión

F S  =
( C  L ) + (W  α  Tanφ  )
Wsenα
Relleno
W
PA
PP
Arena W α

S Cm
90 - α
Material de Baja L Cm = cm L
resistencia L φm + 90
H W
P
P
Figura 4.22 Esquema del método del bloque α φm α − φm
deslizante. β

Polígono de Fuerza

Los valores de las presiones activas y pasivas

pueden obtenerse utilizando las teorías de presión Figura 4.23 Fuerzas que actúan sobre una cuña
de tierras de Rankine o de Coulomb, teniendo simple.
en cuenta el valor de la cohesión movilizada.
A
Una expresión similar puede obtenerse para el C

caso cuando hay dos bloques interrelacionados.

Método de la Cuña Simple

Este método supone una superficie recta de un W
solo tramo, el cual puede analizarse como una
H
cuña simple con la superficie de falla inclinada S

un determinado ángulo con la horizontal (Figuras 3.83 c

Hmáx =
4.23 y 4.24). Una falla de superficie plana puede γ
ser analizada fácilmente con una solución de forma N

cerrada, la cual depende de la geometría de la α

'
pendiente y de los parámetros de fuerza cortante
B
del suelo a lo largo del plano de falla.
Figura 4.24 Análisis de la altura máxima de un talud
Se requiere calcular las siguientes fuerzas: vertical en un suelo cohesivo analizado con cuña simple
(Cornforth, 2005).
• El peso de la cuña (W), descompuesto en la          
fuerza tangente y la fuerza normal, FN y 
FT.
Método de la cuña doble
• FN = W cosα Se analiza una cuña con dos tramos rectos de
superficie de falla (figura 4.25). La cuña superior
• FT = W senα tiene generalmente una pendiente fuerte y la
cuña inferior una pendiente más suave. La cuña
• La fuerza de cohesión, Fc = C x L superior genera una fuerza de empuje sobre la
cuña inferior y esta cuña inferior debe ser capaz de
• La fuerza de fricción, Fφ = FN x Tan φ'. resistir la fuerza impuesta por la cuña superior.
 150 DESLIZAMIENTOS - ANALISIS GEOTECNICO

"Graben" Escarpe
Escarpe reverso
A' A'
Escarpe B E'
D' D
A Escarpe reverso A (α− β)

D D β
α

α (90 − α) (90 − α)
B C
θ
α >> θ B

Figura 4.25 Sección típica de una falla de doble cuña Escarpe

(Cornforth, 2005).
 Escarpe secundario

Generalmente se utiliza para simular fallas sobre

superficies planas duras tales como roca o sobre
superficies planas blandas (manto de arcilla
Superficie de falla basal
blanda).

Debido a que las dos cuñas son geométricamente Grietas

muy diferentes se produce un hundimiento de la

cuña superior (graben) y la cuña inferior se mueve
horizontalmente.
Superficie de falla basal
En el campo este tipo de fallas se reconocen
por la presencia del “graben” (figura 4.26). La
localización, profundidad y extensión del “graben”
permite determinar la profundidad de la falla en Figura 4.26 Formación de “graben” en una falla de

campo. Para el análisis se estudia la estabilidad doble cuña (Cornforth, 2005).
de cada bloque en forma independiente con las
respectivas fuerzas (Figura 4.27).


Adicionalmente a la formación del “graben” Método de la Cuña Triple 

se puede presentar un escarpe secundario en la La falla de triple cuña es común en grandes
parte inferior del deslizamiento formándose en la deslizamientos. Al igual que la falla de doble cuña
práctica tres cuñas. ésta es controlada por los detalles geológicos como
pueden ser una formación de roca o la presencia
A de mantos blandos. En la figura 4.28 se muestra
E cómo en la parte superior del deslizamiento ocurre
un hundimiento (graben) y en la parte inferior
ocurre un levantamiento formándose la tercera
α β cuña.
B

A θ C En la falla de triple cuña las dos cuñas

A E superiores empujan a la cuña inferior para generar
S1
δ P2 el levantamiento del pie del movimiento. Uno de
N1' los factores más importantes para determinar son
α P1 δ P1
α los ángulos de falla de la cuña superior y la cuña
U1 A B
S2 inferior, los cuales no son controlados por las
N2' U2 θ C
características geológicas del talud. El análisis
se realiza estudiando en forma independiente las
Figura 4.27 Fuerzas que actúan sobre las cuñas en
una falla de doble cuña. (Cornforth, 2005).
fuerzas que actúan sobre cada bloque (Figura 4.29).
 ANALISIS DE ESTABILIDAD 151

A
A METODO DE LA ESPIRAL
D
Cuña media LOGARITMICA
A H Cuña inferior
G En el procedimiento de la espiral logarítmica la
superficie de falla se supone que tiene una forma
C
"Graben"
de espiral como se muestra en la figura 4.30.
A
D'
Inicialmente se supone un punto de centro
B B' Levantamiento y un radio r0 para definir la espiral. El radio
H'
de la espiral varía con el ángulo de rotación θ
alrededor del centro de la espiral de acuerdo con
C' G
C la expresión:

Figura 4.28 Esquema típico de una falla de triple cuña

r = ro eθ  φ d
(Cornforth, 2005).
Donde:
φd = es el ángulo de fricción desarrollado el cual
Angulos de las cuñas

depende del ángulo de fricción y del factor de
Cuando se encuentra un caso para el análisis con seguridad.
cuña triple es importante investigar los posibles
ángulos de las cuñas de cabeza y de pie. Existe Los esfuerzos al cortante pueden expresarse
muy poca información de casos históricos y no en esfuerzos totales de acuerdo a la siguiente
existen reglas simples para suponer estos ángulos expresión:
(Cornforth, 2005). Cuando ocurre una falla se
recomienda excavar "apiques" para determinar
c  φ
τ = +σ
estos ángulos con el objeto de poderlos utilizar F F
para el ánalisis de casos similares en la misma
o en términos de las resistencias desarrolladas.
formación geológica.

Generalmente, la inclinación de la cuña τ = Cd + σ  φd

superior es de pendiente fuerte y la de la cuña
inferior es de baja pendiente, la cual puede ser Las ecuaciones de la espiral logarítmica son
hasta de 10º. relativamente complejas para cálculos manuales,
debido a la forma de la superficie de falla.

A Cuña superior

S
W1 Cuña media
S1= c1' I1
P1 δ
P1 W2 F Cuña inferior
α

U1 B W3 G
P3
S2 = c2'I2 δ3 P3
S3 = c3'I3
U2 θ
C
U3

Figura 4.29 Fuerzas que actúan en una falla de triple cuña, (Cornforth, 2005).


152 DESLIZAMIENTOS - ANALISIS GEOTECNICO

Centro
Se calculan los factores de seguridad para todos y
θ tanφd r0 cada uno de los círculos utilizando alguno o varios
r = r0 e
de los métodos existentes y el factor de seguridad
del talud es el mínimo F. S. obtenido de todos los
círculos analizados.

τ Método del Arco Circular

σ
El método del arco circular se le utiliza para
suelos cohesivos solamente (φ = 0). El método fue
propuesto por Petterson en 1916 (Petterson, 1955)
pero solamente fue formalizado por Fellenius en
φd 1922.

Figura 4.30 Talud y superficie de falla espiral

Centros
logarítmica (Frohlich, 1953). R1 = R 2 = R 3 de círculos

R1
Sin embargo, utilizando computadores el R2

análisis es relativamente sencillo. El método R3


de la espiral logarítmica satisface equilibrios
de fuerzas y de momentos y eso hace que el
procedimiento sea relativamente preciso.

Para algunos autores el método de la espiral

logarítmica es teóricamente el mejor procedimiento a) Grilla de centros y círculos de igual radio
para el análisis de taludes homogéneos. Igualmente Centros
el método de la espiral logarítmica es utilizado en de círculos
varios programas de computador para el diseño de
taludes reforzados utilizando geomallas o “nailing”
(Duncan y Wright, 2005).

METODOS DE CIRCULOS DE FALLA Fijar punto común

Las fallas observadas en materiales relativamente

homogéneos ocurren a lo largo de superficies b) Grilla de centros y círculos que pasan por un mismo punto

curvas. Por facilidad de cálculo las superficies Centros

de círculos
curvas se asimilan a círculos, y la mayoría de
los análisis de estabilidad de taludes se realizan
suponiendo fallas circulares.

La localización de los círculos de falla

generalmente se hace dibujando una grilla de
puntos para centros de giro de los círculos y desde
esos puntos se trazan círculos utilizando alguno
de los siguientes criterios (Figura 4.31):
Línea Tangente

• Círculos de igual diámetro. c) Grilla de centros y círculos que son tangentes

a una línea predeterminada

• Círculos que pasan por un mismo punto.

Figura 4.31 Alternativas de procedimiento de
localización de los círculos de falla para análisis de
• Círculos tangentes a una o varias líneas
estabilidad de taludes ( U. S. Corps of Engineers, 2003).
determinadas.
 ANALISIS DE ESTABILIDAD 153

a Métodos de Dovelas
En la mayoría de los métodos con fallas curvas o
circulares la masa arriba de la superficie de falla
se divide en una serie de tajadas verticales. El
número de tajadas depende de la geometría del
r talud y de la precisión requerida para el análisis.
W
Entre mayor sea el número de tajadas se supone
τ que los resultados son más precisos. En los
procedimientos de análisis con tajadas se considera
generalmente equilibrio de momentos con relación
ι al centro del círculo para todas y cada una de las
tajadas (figura 4.33).

Figura 4.32 Fuerzas en un análisis de arco circular Existen una serie de diferencias entre los diversos
(φ = 0) (Duncan y Wright, 2005). métodos que utilizan dovelas, especialmente en
lo referente a las fuerzas que actúan sobre las
paredes laterales de las tajadas (Figuras 4.34 y
En la práctica el método es un caso de la espiral 4.35). El método ordinario o de Fellenius no tiene
 
logarítmica en el   se
cual la espiral  en
convierte   φ     
en cuenta las fuerzas entre tajadas.
 Sin embargo, los análisis son mucho
círculo.
más sencillos para el caso del arco circular y el El método simplificado de Bishop supone que
desarrollo de este método fue anterior al de la las fuerzas laterales entre tajadas son horizontales
espiral logarítmica. y desprecia las fuerzas de cortante y otros métodos
más precisos como los de Morgenstern y Price
En el método del arco circular se supone un utilizan una función para calcular las fuerzas
círculo de falla y se analizan los momentos con entre dovelas.
relación al centro del círculo (Figura 4.32).
Método Ordinario o de Fellenius
El método de Fellenius es conocido también como
clr método Ordinario, método sueco, método de las
F=
Wa Dovelas o método U.S.B.R. Este método asume
superficies de falla circulares, divide el área de
Donde: falla en tajadas verticales, obtiene las fuerzas
c = cohesión. actuantes y resultantes para cada tajada y con la
l = longitud del arco de círculo. sumatoria de los momentos con respecto al centro
r = radio del círculo. del círculo producidos por estas fuerzas se obtiene
W = peso total de la masa en movimiento. el Factor de Seguridad.
a = brazo de la fuerza W con respecto al centro
ai
del círculo

El método del arco circular satisface tanto equilibrio

αi
de fuerzas como equilibrio de momentos. Aunque r
la ecuación fue desarrollada inicialmente para un
valor único de cohesión puede extenderse para Wi
cohesiones diferentes a lo largo del arco circular
y se puede remplazar el término c x l x r por el
αi
término Σ c x l x r. Si

El procedimiento de análisis es sencillo y la

única dificultad es el cálculo del brazo (“a”) para el Figura 4.33 Esquema de un sistema típico de análisis
momento de la fuerza W. Comúnmente el análisis con tajadas o dovelas (Duncan y Wright, 2005).
se realiza en forma manual elaborando gráficos.
154 DESLIZAMIENTOS - ANALISIS GEOTECNICO

x O
-1 -1
A ng ulo ψ = ta n (tan (1 /F tan φ
α
S
b c'
F I

R
N'

d io
A ta
F nφ

Ra
XL B
XR W
EL ψ
W

N'
ER

N
D
S xL − XR

uI
N

U=
C EL − ER
α

Figura 4.34 Fuerzas que actúan sobre una dovela en un análisis de estabilidad de arco circular con dovelas.
(Cornforth, 2005).

El método ordinario o de Fellenius solamente

Las fuerzas que actúan sobre una dovela son satisface equilibrios de momentos y no satisface
(Figura 4.36): equilibrio de fuerzas. Para el caso de φ = 0 el

método ordinario da el mismo valor de factor de
seguridad que el método del arco circular.
• El peso 
o fuerza de gravedad, la cual se puede
descomponer en una tangente y una normal a
la superficie de falla. Los análisis del método de Fellenius son
muy sencillos y se pueden realizar con métodos
• Las fuerzas resistentes de cohesión y fricción manuales o con computador. Debe tenerse en
que actúan en forma tangente a la superficie cuenta que el método ordinario es menos preciso
de falla. que otros procedimientos y la precisión disminuye
a medida que la presión de poros se hace mayor.
• Las fuerzas de presión de tierras y cortante en
Algunos autores recomiendan que el método
las paredes entre dovelas, no son consideradas
ordinario no se utilice para diseño sino solamente
por Fellenius.
como una base de referencia. Generalmente, el
método ordinario da factores de seguridad menores
Al realizar la sumatoria de momentos con respecto
que otros métodos.
al centro del círculo se obtiene la siguiente
expresión: 0 (Centro de giro)
Q

∑ (
C  ∆l + W  α − u ∆l 2 α Tanφ  
  )
F S  =
∑ Wsenα b

Donde: α T1
α = Angulo del radio del círculo de falla con la
vertical bajo el centroide en cada tajada. E2 E1

W = Peso total de cada tajada. α

T2
u = Presión de poros = γ w h w F. Resistente
Δl = longitud de arco de círculo en la base de
la tajada Fuerza Normal
C’, φ’ = Parámetros de resistencia del suelo.
La ecuación anterior se conoce como ecuación de Figura 4.35 Fuerzas que actúan sobre una dovela en
Fellenius. los métodos de dovelas.
 ANALISIS DE ESTABILIDAD 155

El método simplificado de Bishop es uno de los

métodos más utilizados actualmente para el cálculo
de factores de seguridad de taludes. Aunque el
Desprecia
las método solo satisface equilibrio de momentos, se
fuerzas considera que los resultados son muy precisos en
Desprecia W entre comparación con el método ordinario.
las dovelas
fuerzas
entre Aunque existen métodos de mayor precisión
dovelas que el método de Bishop, las diferencias de los
S factores de seguridad calculados no son grandes.
La principal restricción del método de Bishop
N simplificado es que solamente considera superficies
circulares.

Figura 4.36. Fuerzas que actúan sobre una dovela en

el método ordinario o de Fellenius (Duncan y Wright,
Método de Janbú
2005). El método simplificado de Janbú se basa en la
          que
suposición 
las   dovelas son
fuerzas entre
 horizontales y no tiene en cuenta las fuerzas de
Método de Bishop cortante. Janbú considera que las superficies de
Bishop (1955) presentó un método utilizando falla no necesariamente son circulares y establece
dovelas y teniendo en cuenta el efecto de las fuerzas un factor de corrección fo. El factor ƒo depende
entre las dovelas. Bishop asume que las fuerzas de la curvatura de la superficie de falla (figura
entre dovelas son horizontales (Figura 4.37); o sea 4.38). Estos factores de corrección son solamente
que no tiene en cuenta las fuerzas de cortante. aproximados y se basan en análisis de 30 a 40
casos.
La solución rigurosa de Bishop es muy
compleja y por esta razón se utiliza una versión En algunos casos el suponer f0 puede ser una
simplificada de su método, de acuerdo a la fuente de inexactitud en el cálculo del factor de
expresión: seguridad. Sin embargo, para algunos taludes
el considerar este factor de curvatura representa
 c  ∆l  α + (W ′u ∆l  α )  φ   una mejora en el análisis.
∑   α + ( senα  φ  )  FS

F S  =  
∑ Wsen α

Donde:
Δl = longitud de arco de la base de la dovela Ei+1
Wi
W = Peso de cada dovela
C’, φ = Parámetros de resistencia del suelo. Ei
u = Presión de poros en la base de cada
dovela = γ w x h w
α = Angulo del radio y la vertical en cada Si
dovela.
N
Como se puede observar en la ecuación el término
factor de seguridad FS se encuentra tanto en la
izquierda como en la derecha de la ecuación, se
requiere un proceso de interacción para calcular el
factor de seguridad. Figura 4.37 Esquema de fuerzas sobre una dovela en el
método de Bishop simplificado (Duncan y Wrigth, 2005).

        


156 DESLIZAMIENTOS - ANALISIS GEOTECNICO

El método de Janbú solamente satisface equilibrio

de esfuerzos y no satisface equilibrio de momentos.
De acuerdo con Janbú (ecuación modificada): Zi+1

θ Q
 1  θ
f o ∑  c  b + (W − ub ) Tanφ  
  α  ma 
F S  = θ
∑ (W  α ) Zi

Método del cuerpo de Ingenieros (Sueco

modificado)
En el método del Cuerpo de Ingenieros (1970)
la inclinación de las fuerzas entre dovelas es 
Figura 4.39 Paralelismo de las fuerzas entre dovelas
seleccionada por el analista y tiene el mismo valor en el método de Spencer.
para todas las dovelas. El Cuerpo de Ingenieros
recomienda que la inclinación debe ser igual al
promedio de la pendiente del talud. Este método Método de Lowe y Karafiath
satisface equilibrio de fuerzas pero no satisface El método de Lowe y Karafiath (1960) es
equilibrio de momentos. prácticamente idéntico al del Cuerpo de Ingenieros
con la excepción que la dirección de las fuerzas
entre partículas varían de borde a borde en cada
dovela. Su resultado es menos preciso que los que
satisfacen equilibrio completo y al igual que el
L método del Cuerpo de Ingenieros es muy sensitivo
a la inclinación supuesta de las fuerzas entre
d
partículas. Si se varía el ángulo de estas fuerzas
se varía substancialmente el factor de seguridad.

Superficie curva no circular

Método de Spencer
El método de Spencer es un método que satisface
totalmente el equilibrio tanto de momentos como de
1.2 esfuerzos. El procedimiento de Spencer (1967) se
basa en la suposición que las fuerzas entre dovelas
son paralelas las unas con las otras o sea que tienen
el mismo ángulo de inclinación (figura 4.39).
ƒ
o φ=0
Suelos Cohesivos
La inclinación específica de estas fuerzas entre
1.1 partículas es desconocida y se calcula como una
Suelos Mixtos
C -φ
de las incógnitas en la solución de las ecuaciones
Suelos Granulares
de equilibrio. Spencer inicialmente propuso su
C=0 método para superficies circulares, pero este
procedimiento se puede extender fácilmente a
superficies no circulares.

1.0
0 0.1 0.2 0.3 0.4
Spencer plantea dos ecuaciones una de equilibrio
d/L de fuerzas y otra de equilibrio de momentos, las
cuales se resuelven para calcular los factores de
seguridad F y los ángulos de inclinación de las
Figura 4.38 Diagrama para determinar el factor ƒo
para el método de Janbú. fuerzas entre dovelas θ (Figura 4.40).

ƒ
 ANALISIS DE ESTABILIDAD 157

Se utiliza un sistema de ensayo y error para 10 m γ = 21.5 kN/m 3

Trazado
resolver las ecuaciones para F y θ. Se asumen c = 30 kN/m 2
valores de estos factores en forma repetitiva hasta φ = 15 o
que se alcanza un nivel aceptable de error. 0 10 m Bishop
Spencer FS = 1.005
FS = 1.012
Una vez se obtienen los valores de F y θ se Cauce
Janbu
calculan las demás fuerzas sobre las dovelas FS = 0.987
individuales. El método de Spencer se considera
muy preciso y aplicable para casi todo tipo de
geometría de talud y perfiles de suelo y es tal vez el
procedimiento de equilibrio completo más sencillo
para el cálculo de factor de seguridad. (Duncan y
Wright, 2005).
γ = 21.0 kN/m 3
c = 25 kN/m 2
φ = 34o
Janbu
Método de Morgenstern y Price FS = 0.756
El método de Morgenstern y Price (1965) asume 1m
Spencer
que existe una función que relaciona las fuerzas FS = 0.990
de cortante y las fuerzas normales entre dovelas. 0 1m
Esta función puede considerarse constante como en
el caso del método de Spencer o puede considerarse
otro tipo de función. Esta posibilidad de suponer Figura 4.41 Diferencias entre los resultados de varios
una determinada función para determinar los métodos. En cuál de los casos importa saber cuál de los
valores de las fuerzas entre dovelas lo hace un métodos es el que da el verdadero valor del Factor de
método más riguroso que el de Spencer. Seguridad? (Dibujo de Payá).

Sin embargo, esta suposición de funciones 
diferentes tiene muy poco efecto sobre el cálculo Método de Chen y Morgenstern

de factor de seguridad cuando se satisface el El método de Chen y Morgenstern (1983) es
equilibrio estático y hay muy poca diferencia una refinación del método de Morgenstern y
entre los resultados del método de Spencer y el de Price e intenta mejorar los estados de esfuerzos
Morgenstern y Price. El método de Morgenstern en las puntas de la superficie de falla. Chen y
y Price al igual que el Spencer es un método Morgenstern recomiendan que en los extremos de
muy preciso aplicable a prácticamente todas las la superficie de falla las fuerzas entre partículas
geometrías y perfiles de suelo. deben ser paralelas al talud.

b Método de Sarma
El método de Sarma (1973) es muy diferente a todos
A
los métodos descritos anteriormente porque este
considera que el coeficiente sísmico es desconocido
B
XL y el factor de seguridad desconocido. Se asume
RL
θ un factor de seguridad y se encuentra cual es el
EL XR coeficiente sísmico requerido para producir este
W
factor de seguridad. Generalmente se asume
θ
ER que el factor de seguridad es uno y se calcula el
RR coeficiente sísmico requerido para que se obtenga
D este factor de seguridad. En el método de Sarma
S
la fuerza cortante entre tajadas es una relación
N C con la resistencia al cortante. El procedimiento
α
de Sarma fue desarrollado para análisis sísmicos
de estabilidad y tiene algunas ventajas sobre otros
Figura 4.40 Fuerzas que actúan sobre las dovelas en el
método de Spencer. métodos para este caso.
158 DESLIZAMIENTOS - ANALISIS GEOTECNICO

Tabla 4.4 Comparación de los resultados de cálculo de factor de seguridad para varios métodos. (Fredlund y Krahn, 1977).

Factor de seguridad calculado

Talud
Bishop Spencer Janbú Morgenstern-Price Ordinario
Talud 2H:1V 2.08 2.07 2.04 2.08 1.93

Talud sobre una

capa de suelo 1.38 1.37 1.45 1.38 1.29
débil

Talud con
una línea 1.83 1.83 1.83 1.83 1.69
piezométrica

Talud con
dos líneas 1.25 1.25 1.33 1.25 1.17
piezométricas

COMPARACION DE LOS DIVERSOS Por las razones anteriores se prefieren métodos

más sencillos pero más fáciles de manejar como
METODOS es el método simplificado de Bishop. Todos los
métodos que satisfacen equilibrio completo dan
La cantidad de métodos que se utilizan, dan valores similares de factor de seguridad (Fredlund
resultados diferentes y en ocasiones contradictorios y Krahn, 1977, Duncan y Wright, 1980). No
son una muestra de la incertidumbre que caracteriza existe un método de equilibrio completo que sea
los análisis de estabilidad. Los métodos más significativamente más preciso que otro. El método
utilizados por los ingenieros geotécnicos en todo el de Spencer es más simple que el de Morgenstern
mundo son el simplificado de Bishop y los métodos y Price o el de Chen y Morgenstern. Sin embargo,
precisos de Morgenstern y Price, y Spencer. Cada los métodos de Morgenstern son más flexibles para
método da valores diferentes de factor de seguridad tener en cuenta diversas situaciones de fuerzas
(Figura 4.41). entre dovelas. Sin embargo, debe tenerse en cuenta
que la dirección de las fuerzas entre partículas en
Aunque una comparación directa entre los estos métodos no afectan en forma importante el
diversos métodos no es siempre posible, los resultado del factor de seguridad. Para análisis
factores de seguridad determinados con el método sísmico el método de Sarma tiene ciertas ventajas
de Bishop difieren por aproximadamente el 5% con relación a los demás métodos.
con respecto a soluciones más precisas, mientras
el método simplificado de Janbú generalmente, Alva Hurtado (1994) presenta las siguientes
conclusiones al comparar los diversos métodos
subestima el factor de seguridad hasta valores del
(Tabla 4.5).
30%, aunque en algunos casos los sobrestima hasta
valores del 5%. Esta aseveración fue documentada • Cualquier método que satisface el
por Freddlund y Krahn (1977) Tabla 4.4. Equilibrio de Momentos, da el mismo factor
de seguridad en el análisis de φ = 0 con
Los métodos que satisfacen en forma más superficies de falla circular.
completa el equilibrio son más complejos y
requieren de un mejor nivel de comprensión del • El Método Ordinario de Dovelas (Fellenius),
sistema de análisis. En los métodos más complejos da error en el lado conservador para el caso de
y precisos se presentan con frecuencia problemas φ > 0. Con presiones de poros pequeñas, para
numéricos que conducen a valores no realísticos los análisis en función de esfuerzos totales y de
de FS, por exceso o defecto. esfuerzos efectivos, el error es menor del 10%.
 ANALISIS DE ESTABILIDAD 159

Para pendientes casi planas con presiones de • Amplificación de las cargas sísmicas por la
poros altas, el error puede ser mayor del 50%. presencia de suelos blandos.

• Para el análisis de φ = 0 ó φ > 0 con presiones Se han propuesto cuatro métodos de análisis para la
de poros bajas o altas, el método simplificado evaluación de la estabilidad de taludes y laderas, en
de Bishop es adecuado para el análisis de falla el caso de eventos sísmicos (Houston y otros, 1987):
circular. El método es muy estable.
• Método seudoestático en el cual las cargas del
• Numéricamente, sólo hay problemas de sismo son simuladas como cargas estáticas
convergencia cuando los extremos de la horizontales y verticales.
superficie de falla es muy parada, casi
vertical. • Método del desplazamiento o de las
deformaciones, el cual se basa en el concepto
•En los métodos que satisfacen solamente el
de que las aceleraciones reales pueden
equilibrio de fuerzas, el factor de seguridad
superar la aceleración límite permitida
es muy sensible a la inclinación asumida de
produciendo desplazamientos permanentes
las fuerzas laterales. El método de Lowe y
(Newmark, 1965).
Karafiath es razonable para el análisis de
φ>0, pero no conservador (10-15%) para φ=0.
• Método de la estabilidad después del
• Si todas las condiciones de equilibrio son sismo, la cual es calculada utilizando las
satisfechas, la magnitud del error en el factor resistencias no drenadas en muestras de
de seguridad es muy pequeña, usualmente suelo representativas que han sido sometidas
± 5% de la respuesta correcta. previamente a fuerzas cíclicas comparables a
las del sismo esperado (Castro y otros, 1985).
ANALISIS SISMICO
• Método de análisis dinámico por elementos
finitos. Por medio de un análisis en dos o tres
Los eventos sísmicos son capaces de inducir
dimensiones, utilizando un modelo específico
fuerzas de gran magnitud de naturaleza cíclica,
las cuales pueden producir la falla rápida de se pueden obtener detalles relacionados
taludes y laderas. Además, la resistencia al corte con esfuerzos, deformaciones cíclicas o
de un suelo puede reducirse a causa de cargas permanentes (Finn 1988, Prevost y otros,
oscilatorias que generan deformaciones cíclicas, 1985).
o debido a la generación de presiones altas de
poros. La combinación entre la acción de las Los dos primeros métodos son los más utilizados
cargas sísmicas y la disminución de la resistencia en la práctica de la geotecnia debido especialmente
pueden producir una disminución general de la a su facilidad de implementación.
estabilidad. El caso más crítico es el de materiales
no plásticos de grano fino como son los limos o las ANALISIS DE ESTABILIDAD DE
arenas finas. TALUDES UTILIZANDO METODOS
NUMERICOS
En el análisis de estabilidad se requiere analizar
los cinco factores que se indican a continuación: Los mecanismos de falla de los deslizamientos
son con frecuencia muy complejos e incluyen
• Magnitud de la fuerza sísmica. factores muy difíciles de investigar con análisis
convencionales de límite de equilibrio. Los análisis
• Disminución de la resistencia a causa de las de límite de equilibrio se limitan a problemas
cargas oscilatorias. relativamente simples incluyendo muy poca
información del mecanismo de falla. Las fallas de
• Disminución de la resistencia por aumento de los taludes en su gran mayoría son progresivos, no
la presión de poros. se inicia la falla al mismo tiempo como lo suponen
• Fenómeno de resonancia. los métodos de límite de equilibrio.
160 DESLIZAMIENTOS - ANALISIS GEOTECNICO

Tabla 4.5. Diferencias básicas entre diversos métodos de análisis de estabilidad de taludes (Alva Hurtado, 1994).

Condición de Equilibrio Satisfecha Aplicable A

Ecuaciones Forma de la
Procedimiento e superficie de
Mom. Mom. Incógnitas. falla. Cálculos Cálculos en
Vert Horiz
total Dovela Manuales Computador

Método
ordinario de si no no no 1 circular si si
dovelas
Método
de Bishop si no no no n+1 circular si si
Modificado
Método
de Janbú
Procedimiento si si si si 3n cualquiera si si
generalizado
de dovelas.
Métodos de
Spencer y
si si si si 3n cualquiera no si
Morgenstern
y Price.
Método
de Lowe y no no si si 2n cualquiera si si
Karafiath
Método
espiral
de Espiral si - si si 3 si si
logarítmica
Logarítmica

La mayoría de problemas de estabilidad de Los modelos numéricos son muy útiles para
taludes incluyen complejidades relacionados con analizar fallas en las cuales no existe una superficie
geometría, anisotropía, comportamiento no lineal, continua de cortante como es el caso de las fallas
esfuerzos “in situ” y la presencia de procesos por “volteo”. La incorporación de los defectos o
concomitantes como son las presiones de poros y discontinuidades dentro del modelo permiten
las cargas sísmicas. estudiar el comportamiento del talud.

La principal limitación de los métodos de límite Los métodos numéricos de análisis pueden
de equilibrio está en su inhabilidad para tener clasificarse en varias categorías como se muestra
en cuenta las deformaciones, las cuales pueden en la tabla 4.6.
determinar el proceso de falla particularmente en
los procesos de falla progresiva y los que dependen Modelos numéricos continuos
del factor tiempo. Para resolver estas limitaciones Los modelos continuos son los mejores para
se utilizan técnicas de modelación numérica que analizar taludes de suelo o de roca masiva intacta
permiten soluciones aproximadas a problemas que o rocas blandas, o materiales tan fracturados que
no es posible resolver utilizando procedimientos de se comportan como suelos. De estos se conocen
límite de equilibrio. En este aspecto los modelos los programas FLAC, UDEC (Benko-Stead-1993),
númericos son más precisos. PLAXIS entre otros.
 ANALISIS DE ESTABILIDAD 161

El análisis con masas continuas utilizado en la La herramienta es muy poderosa, su utilización

estabilidad de taludes incluye los métodos de es relativamente compleja y su uso se ha venido
elementos finitos y de diferencias finitas. En ambos popularizando para la solución de problemas
métodos el área problema se divide o discretiza prácticos. Wong (1984) menciona la dificultad
en un grupo de subdominios o elementos. La de obtener factores de seguridad a la falla, pero
solución del problema se basa en aproximaciones esta limitación ha sido resuelta por métodos más
numéricas a las ecuaciones de equilibrio, esfuerzo- recientes (Ugai, 1989).
deformación y deformación-desplazamiento.
Alternativamente el procedimiento puede incluir Un análisis por elementos finitos debe satisfacer
aproximaciones a la conectividad de los elementos, las siguientes características:
la continuidad de los desplazamientos y los
esfuerzos entre elementos.
• Debe mantenerse el equilibrio de esfuerzos
en cada punto, el cual es realizado empleando
Métodos de elementos finitos
la teoría elástica para describir los esfuerzos
El método de elementos finitos fue introducido
y deformaciones. Para predecir el nivel de
por Clough y Woodward (1967). El método
esencialmente divide la masa de suelo en unidades esfuerzos se requiere conocer la relación
discretas que se llaman elementos finitos. En el esfuerzo - deformación.
método UDEC el talud se divide en bloques de
acuerdo al sistema de juntas o grietas, los cuales • Las condiciones de esfuerzos de frontera
pueden ser rígidos o deformables. deben satisfacerse.

Estos elementos se interconectan en sus Existe dificultad en la mayoría de los casos

nodos y en bordes predefinidos. El método prácticos reales para definir la relación esfuerzo
típicamente utilizado es el de la formulación de - deformación, por lo difícil que es describir los
desplazamientos, el cual presenta los resultados depósitos de suelos naturales en términos de
en forma de esfuerzos y desplazamientos a los esfuerzo - deformación. Otra limitante es el poco
puntos nodales. La condición de falla obtenida es conocimiento de los esfuerzos reales “in situ” que
la de un fenómeno progresivo en donde no todos se requieren para incorporar en el modelo.
los elementos fallan simultáneamente.

Tabla 4.6 Métodos numéricos para la estabilidad de taludes (Modificado de Deangeli y Ferrero, 2000)

MÉTODO CARACTERISTICAS UTILIZACIÓN

Se asume una malla de elementos con Se aplica a taludes que puedan

Elementos Finitos
sus respectivos nodos y las propiedades considerarse como masas
(FEM)
elastoplásticas de los materiales. continuas sin bloques.
Se utiliza para modelar masa
Diferencias Se elabora una malla con una variedad
rocosa con un alto grado de
Finitas(FDM) de relación esfuerzo-deformación.
fracturación.
Se divide el talud en elementos con sus
propiedades internas y de las uniones Se aplica para analizar
Elementos Distintos
entre los elementos que se pueden inclinación de bloques.
o Discretos (DEM)
mover libremente.
Se discretizan las áreas para Se utiliza para estudiar
Elementos de Borde
poder modelar la ocurrencia de problemas de propagación de
(BEM)
agrietamientos en el talud. grietas.
 162 DESLIZAMIENTOS - ANALISIS GEOTECNICO

Generalmente, se usa un análisis en dos direcciones Existe en la literatura una gran cantidad de
por la facilidad de su aplicación, de acuerdo a sistemas de elementos finitos con sus respectivos
la capacidad de los computadores sencillos. Sin programas de computador. Los elementos finitos
embargo, las soluciones en tres dimensiones son pueden emplearse para estudiar las diversas
cada día más populares. posibilidades de falla en un talud (Figura 4.43),
o para encontrar los efectos de varios sistemas de
El análisis planar o en dos direcciones asume estabilización para el estudio en casos generales,
cero esfuerzo o cero deformación en las superficies donde las propiedades de los suelos o rocas y
laterales del modelo, por lo tanto para que se condiciones de frontera se pueden suponer. En la
simulen las condiciones de campo se requiere que estabilidad de taludes los métodos de elementos
existan esas condiciones. finitos en 3-D permiten analizar condiciones que
los métodos de equilibrio límite no permiten. El
El empleo de análisis en dos direcciones se análisis en 3-D es tal el mayor aporte de los
puede ampliar aplicando al modelo una carga elementos finitos a la estabilidad de taludes
hidrostática lateral. (Figura 4.44).
En la figura 4.42 se muestra una malla típica
El método de elementos finitos es hoy el más
para el análisis de un talud por elementos finitos
utilizado y probablemente el modelo numérico más
(Ashford y Sitar, 1994). Generalmente las mallas
versátil para el análisis de estabilidad de taludes.
analizadas contienen elementos de tamaño
Las principales ventajas y desventajas del método
uniforme con anchos (w) y alturas (h) iguales. El
de elementos finitos se resumen en los siguientes
tamaño y forma de los elementos influye en forma
párrafos (Carter y otros, 2001).
importante sobre los resultados obtenidos.

Es común que entre más pequeños sean Ventajas de los métodos de elementos finitos:
los elementos se obtienen mayores niveles de
esfuerzos de tensión en la cresta del talud. La • Se puede considerar el comportamiento no
altura del elemento es tal vez el factor más lineal de los materiales en la totalidad del
importante y se recomiendan por lo menos diez dominio analizado.
niveles de elementos entre el pié y la cabeza
del talud para simular en forma precisa el • Es posible modelar la secuencia de excavación
comportamiento del talud. incluyendo la instalación de refuerzos y
sistemas de estructura de soporte.

• La falla es progresiva.
D
2H
• Los detalles estructurales de juntas o fisuras
H cercanas pueden modelarse utilizando una
técnica de homogenización.

2H W
h • Se puede introducir un comportamiento de
los materiales dependiente del tiempo.
Límite
Límite

• El sistema de ecuaciones es simétrico con

excepción de los problemas elastoplásticos y
de flujo.

• Se puede utilizar una formulación convencional

de deformaciones para la mayoría de las
posibilidades de carga.

Figura 4.42 Malla típica 2D para el análisis de un talud • Se han desarrollado formulaciones especiales
vertical por elementos finitos (Ashford y Sitar, 1994). para incluir análisis del agua subterránea.


 163
PLAXIS V8

0.00 5.00 10.00 15.00 20.00 25.00 30.00 35.00

ANALISIS DE ESTABILIDAD
25.00

• Existe mucha experiencia sobre el uso de

estos modelos y los programas de software
20.00
han sido actualizados teniendo en cuenta esas
experiencias.

15.00
Desventajas de los métodos de elementos
finitos.
10.00
• Debido a que el sistema de ecuaciones es muy
grande, se requieren tiempos prolongados y
capacidades altas de memoria dependiendo
de la estructura general de los taludes y la
5.00

implementación de los algoritmos del código

PLAXIS V8 de elementos finitos.
0.00
0.00 5.00 10.00 15.00 20.00 25.00 30.00 35.00

25.00
Connectivities
• La totalidad del volumen del dominio
analizado tiene que discretizarse.
PLAXIS
Project description

desliz
• Algunos modelos requieren de algoritmos
Project name Date User name

Finite Element Code for Soil and Rock Analyses

20.00 16/09/08 Koxhiyoki Kabuto, Japan
sofisticados de acuerdo al tipo de material
Version [Link]

constitutivo utilizado.
15.00

• El método no es apropiado para rocas muy

fracturadas o suelos altamente fisurados
10.00
cuando las discontinuidades se encuentran
distribuidas en forma no uniforme y controlan
el comportamiento mecánico de los taludes.

Las anteriores desventajas son mucho más

5.00

pronunciadas en el análisis 3D y menos fuertes en

PLAXIS V8 el análisis 2D. Sin embargo, teniendo en cuenta
la tendencia a utilizar modelos 3D el manejo de
0.00
0.00 5.00 10.00 15.00 20.00 25.00 30.00 35.00

25.00
Deformed Mesh
Extreme total displacement 69.90*10-3 m
los modelos de elementos finitos es relativamente
(displacements scaled up 50.00 times) complejo.

PLAXIS
Project description

20.00 desliz
Project name Step Date User name

Finite Element Code for Soil and Rock Analyses desliz 9 16/09/08 Koxhiyoki Kabuto, Japan
Version [Link]

15.00

10.00

Y
X
5.00

0.00

Deformed Mesh
Extreme total displacement 69.90*10-3 m
(displacements scaled up 100.00 times)

PLAXIS
Project description

Figura 4.43 Modelación de falladesliz utilizando modelo de

Project name Step
Figura 4.44 Malla típica 3D para un talud utilizando
Date User name

Finite Element Code for Soil and elementos finitos.

desliz (PLAXIS
9 ).
16/09/08 Koxhiyoki Kabuto, Japan
FLAC-3D.
Rock Analyses
Version [Link]


164 DESLIZAMIENTOS - ANALISIS GEOTECNICO

Evaluación del factor de seguridad Si se supone la envolvente de falla no lineal

utilizando elementos finitos (Criterio de Maksimovic), las superficies críticas
Ugai (1989) desarrolló un método para calcular de falla son menos profundas y los factores de
el factor de seguridad utilizando el criterio de seguridad son significativamente menores.
Mohr-Coulomb por medio de elementos finitos. El
factor de seguridad es evaluado realizando una Métodos de diferencias finitas
reducción gradual de los parámetros de resistencia En el método de diferencias finitas los materiales
al cortante c’ y φ' del suelo induciendo una falla son representados por zonas para formar una malla
del análisis. Inicialmente la fuerza de gravedad se de acuerdo a la geometría y se puede seleccionar
aplica en estado elástico para obtener la primera una variedad de relaciones esfuerzo/deformación
distribución de esfuerzos en todo el talud. Luego (FLAC 1998). El método se basa en el esquema
la reducción gradual de resistencia va a producir de cálculo de “Lagrange”, el cual permite modelar
un esfuerzo residual en los elementos fallados deformaciones de gran escala y el colapso de los
evaluando así la fuerza residual. materiales.

El valor inicial de F se asume lo suficientemente El esquema general de análisis consiste en

pequeño para obtener como resultado un problema el re-equilibrio del sistema y el estudio de las
elástico. Luego el valor de F se va aumentando condiciones de falla (Figura 4.45). El método de
etapa por etapa hasta que se desarrolle una falla diferencias finitas es poco utilizado en estabilidad
global del talud (Popescu y otros, 2000) de taludes con excepción de los análisis de flujo,
consolidación y transporte de contaminantes.
Este método se le conoce como modelo de Sin embargo, el método puede manejarse para
elementos finitos de reducción de resistencia al utilizarse en remplazo o en complemento del
cortante (SSRFEM). En forma similar se han método de elementos finitos. El método de
desarrollado procedimientos para calcular el diferencias finitas tiene la ventaja que no requiere
factor de seguridad para envolventes de falla no la solución de gran cantidad de ecuaciones y es
lineales (Tanaka y Sakai, 1993). Los resultados más fácil introducir modelos especiales de suelo.
del círculo crítico de falla y factor de seguridad son Sin embargo, el modelo de diferencias finitas es
diferentes si se asume que la envolvente de falla muy complejo en 3D y existe muy poca experiencia
es o no lineal. de su uso en estabilidad de taludes.

Figura 4.45 Análisis de un talud utilizando un modelo elasto-plástico utilizando diferencias finitas en el código
FLAC (Stead y otros, 2000).
 ANALISIS DE ESTABILIDAD 165

Método de elementos de borde (BEM)

El método de elementos de borde ha adquirido
gran importancia en el análisis de estabilidad de
taludes en materiales discontinuos o fracturados
y es una alternativa al método de elementos finitos
(Figura 4.46). Igualmente permite trabajarlo en
forma conjunta (Beer y Watson, 1992).

Ventajas del método de elementos de borde

a) inicial • La discretización del área y no del volumen

reduce los esfuerzos de procesamiento.

• La discretización de áreas conduce a sistemas

de menor número de ecuaciones y se requiere
menos tiempo de computador y capacidad de
disco.

• Se puede modelar fracturas e interfaces entre

materiales localizados donde se requiera.

Desventajas del método de elementos de

borde
b) intervalo 1
• Solo se puede considerar comportamientos
en materiales elásticos con excepción de las
interfaces y discontinuidades.

• Los sistemas de ecuaciones son generalmente

no simétricos.

• No es posible modelar secuencias de

excavación ni estructuras de soporte.

• La formulación estándar no permite trabajar

con gran cantidad de juntas en la roca
c) intervalo 2 distribuidas en forma aleatoria.

• Existe poca experiencia comparativamente

con el método de elementos finitos.

En la figura 4.47 se muestra como en un talud rocoso

se inicia el proceso de deslizamiento utilizando
una técnica combinada de elementos de borde
con elementos finitos. Se puede observar cómo se
van presentando y ampliando los agrietamientos
a tensión para formar fracturas semi-verticales
normales a la dirección del movimiento.

d) intervalo 3 A medida que la densidad de estas fracturas

aumenta se va desarrollando una superficie de
Figura 4.46 Modelo de fractura utilizando elementos cortante o superficie de falla semi-curva (Eberhard
de borde con modelo ELFEN (Stead y otros, 2006). y otros, 2004).



166 DESLIZAMIENTOS - ANALISIS GEOTECNICO

W
Modelo combinado de Elementos finitos 2500 E

y Elementos de Borde s Escarpe (Mayo 9, 1991)

ei s
Teniendo en cuenta las desventajas que se indicaron 2250 ran ss
Pa e i
anteriormente de los métodos de elementos finitos Or
to n
s
Escarpe (Abril 18, 1991)

y elementos de borde, se pueden minimizar estas 2000 eis

a ran
P
limitaciones utilizando los dos métodos en forma
combinada. Estos modelos combinados se pueden 1750
ei s
s
obtener discretizando el suelo o la roca dentro de to n Topografía
Or iss)
e Actual
una determinada zona particular de interés por 1500 (N
Topografía antes
ejemplo, alrededor de un túnel (Beer y 200 m del deslizamiento
Watson, 1992); Sin embargo, la modelación de
a) Marco Geológico
discontinuidades importantes es complicada
y se genera un sistema de ecuaciones no
simétricas en el modelo combinado. Como este
sistema es relativamente nuevo, todavía se
requiere resolver algunas dificultades, aunque
ya existen algunas experiencias positivas.
Zona de inestabilidad
Métodos de elementos distintos o
Croquis de la
discretos superficie de
Los métodos numéricos continuos (elementos deslizamiento
finitos y diferencias finitas), no permiten analizar
en forma precisa la influencia de la estructura b) Resultado del modelo
geológica. Aunque los métodos continuos pueden
σ1
modificarse para acomodar las discontinuidades,
este procedimiento es difícil y complicado. Los
métodos numéricos discontinuos por su parte σ1
permiten modelar en forma relativamente sencilla σ1
taludes donde el mecanismo de falla está controlado
por el comportamiento de las discontinuidades.
Cuando un talud en roca tiene más de dos grupos τ
de discontinuidades es conveniente utilizar τ
modelos discontinuos. Sin embargo, en ocasiones
se requiere trabajar con modelos que permitan
tanto elementos continuos como discontinuos.
σ1
σ1
Los modelos discontinuos los hay de diferentes
variaciones así: σ1

c) Propagación de Grietas
• Método de elementos distintos o discretos.

• Métodos de análisis de flujo de partículas. Figura 4.47 Resultados de un modelo combinado

de elementos finitos y elementos de borde. A (Marco
• Métodos de deformaciones discontinuas. geológico). B (resultado del modelo). C (Propagación de
grietas) (Eberhard y otros, 2004).
Un modelo discontinuo trata las masas 
de
roca como un ensamblaje de elementos distintos

permitido un avance importante en la modelación
de bloques o cuerpos interactuantes que
están
de taludes en roca. La base del método de elementos
sometidos a cargas externas y se espera que
discretos es que la ecuación dinámica de equilibrio
tengan movimientos significativos en el tiempo
 para cada bloque en el sistema es formulada y
(Figuras 4.48 y 4.49). Esta metodología se le
resuelta repetitivamente hasta que las condiciones
conoce como de elementos discretos. El desarrollo
y leyes de contacto y de borde se satisfacen.
de los procedimientos de elementos discretos ha
 ANALISIS DE ESTABILIDAD 167

Esto representa una interacción no lineal compleja

3
f ni
entre los diversos bloques. Los factores externos f
2
vi
como las presiones de poros y las fuerzas sísmicas
también pueden simularse sobre los elementos f
4
m Fi
discretos.

Los métodos de elementos distintos o discretos f1

es particularmente útil para análisis de caídos,
Fuerzas actuando
inclinaciones, y deslizamientos diversos en macizos sobre la partícula k
de roca (Stead y otros, 2000).

Los elementos discretos se basan en la mecánica

de medios discontinuos donde el comportamiento
del talud está gobernado principalmente por
el efecto de las juntas y grietas. En estos casos k
el método de elementos finitos no es aplicable y
se requiere trabajar con elementos discretos o
independientes.

Este método está caracterizado por lo siguiente:

• Se calculan deformaciones finitas y rotaciones

de cada uno de los bloques suponiendo los
bloques rígidos o deformables. Sistema de partículas discretas

• Los bloques que originalmente se encuentran

conectados pueden separarse en el proceso de Figura 4.48 Fuerzas que actúan sobre un sistema de
análisis. partículas discretas.


Etapa 1 Etapa 2 Etapa 3

 Etapa 4 Etapa 5


Figura 4.49 Esquema de un análisis de falla de un talud con inclinación reversa con elementos discretos.


 168 DESLIZAMIENTOS - ANALISIS GEOTECNICO

• Se pueden desarrollar automáticamente

contactos nuevos entre los bloques que se
desplazan o rotan.

Los códigos UDEC y 3-DEC son los más utilizados,

y ambos emplean esquemas de diferencias finitas
como en el programa FLAC.

El modelo de elementos discretos no es

comparable con el modelo de elementos finitos,
debido a que en cada uno de estos modelos los
materiales se comportan en forma diferente. La
principal desventaja del método de elementos
discretos es la dificultad para establecer etapas de
construcción. Además, el sistema 3-DEC consume
mucho tiempo de computador. Figura 4.51 Análisis de deformación discontinua (Chen
y Ohnishi, 1999).
Para elaborar un modelo de elementos discretos
se requiere experiencia para determinar los 
valores más apropiados para los parámetros de
entrada tales como la rigidez de las juntas. Estos Métodos de deformación discontinua
parámetros generalmente no se pueden obtener de Los métodos de deformación discontinua permiten
ensayos de laboratorio y al suponerlos conduce a simular deslizamientos en roca, inclinaciones
problemas de cálculo. y caídos (Chen y Omishi, 1999). La figura
4.51 muestra un análisis de falla utilizando
Métodos de flujo de partículas deformaciones discontinuas.
Una variante de los métodos de elementos distintos
es el de la modelación de flujo de partículas (Itasca, Cuál modelo utilizar para cada
1996). Esta metodología permite simular el flujo problema
de partículas granulares debido a la fricción entre Cada problema es diferente y es difícil establecer
partículas (Figura 4.50). También es posible criterios generales sobre qué modelo se debe
simular materiales intactos o bloques dentro utilizar en cada caso. En algunas ocasiones se
del flujo utilizando uniones entre partículas. Se pueden utilizar varios tipos de modelo y se debe
pueden además formar “clusters” de partículas escoger aquel del cual se tenga mayor experiencia
para simular bloques intactos. y familiaridad. En la figura 4.52 se muestra en
forma esquemática que los métodos de límite de
equilibrio son muy útiles para el análisis sencillo
Si los esfuerzos exceden la resistencia de las
de estabilidad de taludes.
uniones se produce la rotura de los elementos
internos. Los métodos de flujo de partículas
Si los patrones de comportamiento del suelo
permiten analizar casos de licuación de suelos. son complejos, se requiere un modelo de elementos
finitos o diferencias finitas y si los materiales se
encuentran fracturados se recomienda utilizar un
modelo de elementos discretos o de elementos de
fuerza aceleración velocidad desplazamiento borde.

Análisis en tres dimensiones

La mayoría de los deslizamientos posee una
geometría en tres dimensiones; varios autores
han presentado métodos de análisis, de los
Figura 4.50 Las fuerzas entre partículas se convierten cuales merece especial interés el de Yamagami
en velocidades y deformaciones en un continuo de flujo. y Jiang (1996).
              

 ANALISIS DE ESTABILIDAD 169

I : Analisis cinemático
y de equilibrio limite

II : Métodos numéricos
continuos y discontinuos
III : Elementos Hibridos finitos
Discretos con fractura

Traslación simple Traslación Rotación y Traslación

o Rotación compleja complejas

Fallas en gradas incluyendo Ablandamiento interno y

Corte sobre superficies
rotura de materiales a lo corte con ensanche de
basales, laterales
largo de superficies con fracturas y degradación de
y traseras suaves.
puentes de materiales resistencia.
intactos.

Mecanismos de daño: Mecanismos de daño:

Mecanismos de daño: Falla plástica y falla fragil.
Rotura de materiales,rotura
Rotura de materiales Degradación progresiva
de asperitas, falla fragil,
y asperitas, falla de la resistencia,
progresiva. fractura de rocas, falla
progresiva. flujo y corte.

Incrementode complejidad

Mecanismo de falla Mecanismo de falla Mecanismo de falla

Falla en gradas Falla profunda de

Falla multiples bloques
Planar multiples con
Puentes de corte
Discontinuidad roca interno
intacta Transición de
frágil a ductil

Figura 4.52 Diagrama que muestra el tipo de modelo que se recomienda utilizar de acuerdo a la complejidad de los
movimientos (Stead y otros, 2006).

Este método utiliza las ecuaciones de factor de Análisis Numérico 3-D

seguridad de Janbú junto con un esquema de Los métodos de elementos finitos utilizan con

minimización basado en programación dinámica. frecuencia análisis 3-D. Estos modelos son muy

Con este programa se obtiene la superficie de falla útiles para la evaluación de estabilidad en macizos
crítica en tres dimensiones, sin restricción a la rocosos donde el efecto de las discontinuidades
forma de la falla, su respectivo factor de seguridad actúa en tres dimensiones, situación que es muy
y la dirección del movimiento (Figura 4.53). díficil modelar utilizando modelos 2-D.

Análisis de equilibrio límite 3-D El análisis de elementos finitos 3-D tiene las
Ocasionalmente se realizan análisis de estabilidad siguientes desventajas:
de equilibrio límite en tres dimensiones. Al igual
que con los métodos 2-D se requiere realizar una • Es muy complejo discretizar el volumen total en 3-D
serie de suposiciones para que el problema sea
estáticamente determinado. La mayoría de métodos • El tiempo de corrida del computador y el
3-D tiene limitaciones importantes y son útiles espacio requerido son muy grandes.
solamente para conocer el efecto de la situación 3-
D sobre una determinada superficie de falla. Los • No son viables para rocas o suelos muy fisurados
métodos de equilibrio límite 3-D se utilizan muy
con fracturas en muchas direcciones.
poco en diseño (U. S. Corps of Engineers, 2003).
170 DESLIZAMIENTOS - ANALISIS GEOTECNICO

• Se necesitan algoritmos muy sofisticados. REFERENCIAS CAPITULO 4

El uso de técnicas 3D con diferencias finitas o con Abramson, L., Lee, T.S., Sharma, S., Boyce, G.M.(2002).
elementos discretos tiene actualmente muchas “Slope Stability and Stabilization Methods”. John
limitaciones. No se han desarrollado hasta el Wiley and Sons, INC. pp. 1- 712.
momento (2008) herramientas eficientes para Abramson, L.W., Lee T.S., Sharma S., Boyce G.M.,(1996),
el análisis 3-D, comparadas con los procesos “Slope stability and stabilization methods”. Wiley-
elaborados para elementos finitos. Interscience, 629 p.
Alva Hurtado, J. E. (1994). “Notas sobre Análisis de
estabilidad de taludes”. Facultad de Ingeniería Civil,
Universidad Nacional de Ingeniería.
Ashford S.A., Sitar N. (1994). “Seismic response of steep
natural slopes” Report UCB/EERC-94/05 University
o of California at Berkeley , 207 p.
θm=0
Barnes, G.E. (1991), “A simplified version of the Bishop
and Morgenstern Slope-Stability Carts”, Canadian
a) Planta O' x'(m)
Geotechnical Journal, Vol. 28, No.4, pp. 630-632.
θ=154 o Beer, G. and Watson, J.O. (1992) “Introduction to Finite
o 10.0
θ n=180 and Boundary Element Methods for Engineers”,
John Wiley & Sons, New York.
Benko B., Stead D. (1993). “The characterization and
5.0
prediction of landslide movements using numerical
9. .0

modelling techniques”. Seventh international

8.317
7. 3
10

6. 50

Dirección de Deslizamiento
5.867

conference and field workshop on landslides in Czech

5.0 3
m

y'(m) 10.0 5.0

and Slovak Republics , pp. 131-136.
Bishop A.W., Morgenstern N. (1960). “Stability
Coefficients for Earth Slopes”. Geotechnique Vol
10, No. 4,.
Bishop A.W. (1955). “The Use of the Slip Circle in the
Stability Analysis of Slopes”. Geotechnique. Vol 5.
Carter, J.P., Desai, C.S., Potts, D.M., Schweiger, H.F.
y Sloan, S.W.(2001). “Computing and Computer
Modelling in Geotechnical Engineering”. pp.1 – 96.
5.0m

Castro, G., Poulos S.J., Leathers F.D., (1985). “Re-

examination of slide of lower San Fernando Dam”.
Journal of the geotechnical engineering division,
z' ASCE, Vol. 111, No. GT-9, pp. 1.093-1.107.
5.0m

Fs ,min=1.11
b) Vista 3-D y' x' Charles, J.A. , Soares M.M., (1984). “Stability of
5.0m .0m
5 compacted rockfill slopes”, Geotechnique, Vol. 34,
No.1, pp. 61-70.
Chen, WF., Giger M.W., (1971). “Limit analysis of slopes”,
Journal of the soil mechanics and foundation division,
Figura 4.53 Dirección del deslizamiento y superficie

ASCE, Vol. 97, No. SM-1, January, pp.19-26.
de falla crítica en un análisis en tres dimensiones
 Chen, Z.Y. y Morgenstern, N.R. (1983). “Extensions to
(Yamagami y Jiang, 1996).
generalized method of slices for stability analysis”,
Canadian Geotechnical Journal, 20(1), 104-119.
Chen, G., Ohnish, Y. (1999). “Slope stability using
discontinuous deformation analysis.” In Amadei,
Proc. 37th. U.S Rock Mech. Sym., Rock Mechanics for
Industry, Vail, Colorado. A.A. Balkema, Rotterdam,
pp. 535-541
Chugh, A.K., (1981). “Pore water pressure in natural
slopes”, International Journal for numerical and
analytical methods in geomechanics, Vol. 5, pp. 449-
454.
 ANALISIS DE ESTABILIDAD 171

Clough, R.W., Woodward R.J., (1967). “Analysis of Janbú N. (1954). “Stability Analysis of Slopes with
Embankment Stresses and deformations”. Journal Dimensionless parameters”. Harvard Soil
of geotechnical division, ASCE, July, pp 529-549. Mechanics Series No. 46.
Cornforth, D.H. (2005). “Landslides in practice”. Janbú, N. (1968). “Slope stability computations”. Soil
Investigation, analysis, and remedial preventative mechanics and foundation engineering report.
options in soils. John Wiley and Sons Inc. Hoboken, Technical University of Norway, Trondheim.
New Jersey. Krahn, J. (2004). “Stability Modeling with Slope/W an
Cousin, B.F., (1978), “Stability charts for simple earth Engineering Methodology”. Pp 408.
slopes”, Journal of Geotechnical division, ASCE, Vol. Lowe, J., Karafiath, L. (1960). “Effect of anisotropic
104, No. GT-2, February, pp. 267-279. consolidation on the undrained shear strength of
Duncan J.M. (1996). “Soil slope stability analysis”. compacted clays”, ASCE Research Conference on
Landslides investigation and mitigation. Special Shear Strength of Cohesive Soils, pp 837-858.
report 247. Transportation research Board. National Morgenstern N. , Price V. (1965). “The Analysis of the
research council. Pp. 337-371. Stability of General Slip Surfaces”, Geotechnique
Duncan J.M., Wright S. G. (2005). “Soil Strength and Vol 15, March.
Slope Stability”. John Wiley and [Link]. Hoboken. Newmark, NM. (1965). “Effects of Earthquakes on
New Jersey 297p. Dams and Embankments”. Geotechnique, Vol. 15,
Duncan, J.M., Wright S.G. (1980). “The accuary of No.2, pp. 139-160.
equilibrium methods of slope stability analysis”. O´Connor, M.J., Mitchell, (1977). “An extension of the
Engineering Geology. 16(1), 5-17. Bishop and Morgenstern slope stability charts”
Fellenius W. (1922). “Statens Jarnjvagars Geoteknniska Canadian Geotechnical Journal. Vol. 14 No.1 pp.
Commission, Stockholm, Sweden” 144-151.
Fellenius, W. (1927), “Erdstatiche Bereshnungen mit Petterson, K.E. (1955). “The early history of circular
Reibung und Kohasion, Ernst, Berlín (En aléman). sliding surfaces”, Geotechnique 5(4), 275.296.
Fellenius W. (1936). “Calculation of the Stability of Popescu M., Ugai K., Trandafir A. (2000). “Linear
Earth Dams”. Second Congress on Large Dams, versus Non-linear Failure envelopes in LEM and
International Commission on Large Dams of FEM slope stability Analysis”. 8th international
the world power conference. Vol. 4, pp 445-462. symposium on landslides Cardiff. pp 1227-1234.
Washington D.C. Prevost, J.H., Abdel-Ghaffar A.M., Lacy S.J., (1985),
Finn, W.D. (1988). “Dynamic analysis in geotechnical “Nonlinear dynamic analysis of an earth Dam”,
engineering”. Proceedings of earthquake engineering Journal of the geotechnical engineering division,
and soil dynamics II-Recent Advanced in ground ASCE, Vol. 111, No. GR-7, pp. 882-897.
motion evaluation, J.L. Von Thun, Edd. Park City, Sarma, S.K. (1973). “Stability analysis of embankments
Utah: ASCE, Geotechnical Special Publication No. and slopes”, Geotechnique, 23(3), 423-433.
20. Skempton, A.W.(1970). “First-time slides in over
Freddlund D.G. , Krahn J.(1977). “Comparison of slope consolidated clays”. Geotechnique, Vol. 20, No. 4,
stability methods of analysis”, Canadian geotechnical pp. 320-324.
Journal, Vol. 14, pp. 429-439. Skempton A.W. (1977). “Slope stability of cuttings
Frohlich, O.K. (1953). “The factor of safety with in London clay”. In Proc., Ninth international
respect to sliding of a mass of soil along the arc conference on soil mechanics and foundation
of a logarithmic spiral”. Proceedings of the 3rd engineering, Tokyo, pp. 261-270.
International Conference on Soil Mechanics and Skempton A.W. (1985). “Residual strength of clays
Foundation Egineering, Switzerland, vol. 2. pp. in landslides, folded strata, and the laboratory”.
230-233. Geotechnique, Vol. 35, No. 1, pp. 3-18.
Hoek, E. , Bray J.W., (1977), “Rock slope engineering”. Spencer, E. (1967), “A method of analysis of the stability
London: Institute of mining and metallurgical of embankments assuming parallel inter-slice
engineering. forces”. Geotechnique, Vol 17 pp 11-26.
Houston, S.L., Houston W.N., Padilla J.M., (1987). Stead D., Benko B., Eberhardt E., Coggan J. S. (2000).
“Microcompute. aided evaluation of earthquake- “Mechanism of failure of complex landslides: A
induced permanent slope displacements”, Numerical Modelling Perspective”. 8th International
Microcomputers in civil engineering, Vol.2, symposium on landslides, Cardiff. pp 1401-1406.
Amsterdam: Elsevier, pp. 207-222. Stead, D., Eberhardt, E., Coggan, J.S. (2006).
Hunter, J.H., Schuster R.L. (1968), “Stability of “Developments in the characterization of complex
simple cuttings in normally consolidated clays”, rock slope deformation and failure using numerical
Geotechnique, Vol. 18, No. 3, pp. 373-378. modelling techniques”. Engineering Geology 83.
Pp.217– 235.
172 DESLIZAMIENTOS - ANALISIS GEOTECNICO

Tanaka, T., Sakai, T. (1993). “Progressive failure and U. S. Corps of Engineers (2003). “Slope Stability” .
scale effect of trap-door problems with granular Engineering and Design. Engineer manual. EM
materials”. Soils and Foundations, 33, No 1. pp. 1110-2-1902.
11-22. Wong, F.S. (1984), “Uncertainties in FE modelling of
Taylor, D. W. (1966). “Fundamentáis of Soil Mechanics”. slope stability”, Computers & structures, New York;
2 edition, John Wiley y Sons, New York, 700 pp. Pergamon, Vol. 19, No.5/6, pp. 777-791
Terzaghi, K. and Peck R.B. (1967). “Soil Mechanics in Yamagami T, Jiang, J.C. (1996). “A generalized
Engineering Practice”. New York: Wiley. method for determining the critical slip surface
Ugai, K. (1989). “A Method of Calculation of total factor in three-dimensional slope stability analysis”.
of safety of slopes by Elasto-Plastic FEM”. Soils Seventh International Symposium on Landslides.
and Foundations, Vol. 29, No 2. pp 190-195. (En Trondheim, pp 1433-1439.
japones)