UNIVERSIDAD PERUANA DE CIENCIAS APLICADAS

FACULTAD DE INGENIERÍA

PROGRAMA ACADÉMICO DE INGENIERÍA CIVIL

<< Análisis tridimensional de estabilidad de taludes en la Costa Verde en

el distrito de Miraflores considerando la influencia de los aspectos
geomorfológicos mediante el método de Elementos Finitos >>

CAPSTONE PROJECT

AUTOR(ES)

Esli Vega, Roni Edison U201614812

Villanueva Bardalez, José Emanuel U201516498

ASESOR

Pérez Zúñiga, Álvaro Javier

Lima, 30 MAYO 2022
I
 CAPÍTULO 1. IDENTIFICACIÓN DEL PROBLEMA

1.1 DESCRIPCIÓN DEL PROBLEMA

En la actualidad, los ingenieros geotécnicos todavía no realizan de manera habitual análisis

tridimensional de estabilidad de taludes. Esto se debe a que los análisis 3D presentan una
complejidad mayor. Los análisis 3D se caracterizan por trabajar con códigos numéricos;
los cuales toman un cierto tiempo mayor para realizar la determinación de parámetros de
material y el cálculo del código (Bovolenta y Bianchi, 2020). Desde la década de 1970,
diversos investigadores propusieron diversos análisis tridimensionales para taludes. Este
intéres está impulsado; ya que en gran parte de fallas en los taludes es de cáracter
tridimensional.

Se han propuesto varios métodos 3D, basados en la extensión simple de los métodos 2D
convencionales, para estimar los efectos 3D en la seguridad de taludes. No obstante, estos
métodos asumieron que el cuerpo deslizante se movía a lo largo del plano de simetría. Esta
suposición podría conllevar a un error de cálculo del factor de seguridad para taludes bajo
condiciones asimétricas. Las variaciones de condiciones geomorfológicas no son
analizadas adecuadamente en los métodos bidimensionales de estabilidad del talud. Por lo
que, el factor de seguridad del talud de la Costa Verde presenta mayor margen de error; y,
en consecuencia, ocurrirán deslizamientos por el análisis no requerido.

Por ello, en el siguiente proyecto, se realizará un análisis tridimensional que considere cada
uno de los factores en la geomorfología del talud; puesto que, modifican las condiciones
mecánicas y de resistencia al corte de la tierra.

1
Figura 1: Desplazamientos permanentes durante el ensayo de vibración.

Fuente: Huang et al (2018)

Figura 2: Mapa causal del análisis del problema

Fuente: Elaboración propia

2
1.2 DELIMITACIÓN DEL PROBLEMA

Nuestro trabajo de investigación será realizado en la Costa Verde ubicada en el distrito

de Miraflores, Lima, Perú. Fue construida en los años 1960 por el señor Luis Bedoya
Reyes. Se encuentra al pie del acantilado frente al océano Pacifico. Sufre constantes
desprendimientos de tierra y rocas, donde la vía más afectada es la de sureste a noroeste. El
talud que se encuentra en el distrito de Miraflores en la Costa Verde presenta una
topografía irregular debido a las salientes (cárcavas) o protuberancias que existen. Esta
zona del talud tiene en promedio una altura de 66 metros.

El período de tiempo de la procedencia de datos para la investigación será no mayor de 10

años. Consideraremos estudios previos geotécnicos del talud de la Costa Verde.

Para la presente investigación, se considerará los ensayos de mecánica de suelos del talud
de la Costa Verde de investigaciones anteriores. Puesto que será necesaria información
como la clasificación del suelo, la cohesión, el ángulo de fricción interna; los cuales
pueden obtenerse en ensayos de corte directo anteriores. Asimismo, será necesario
información de las propiedades físicas como la gravedad específica, contenido de agua,
relación de vacíos, densidad y el índice de plasticidad. Además, consideraremos datos del
SENAMHI relacionado a climatología del lugar (precipitación, temperatura, humedad,
etc).

3
 Figura 3: Modelo mallado de pendiente de talud.

Fuente: Ma et al (2021)

1.3 FORMULACIÓN DEL PROBLEMA

¿De qué manera la variación geomorfológica del talud ubicado en el distrito de

Miraflores, en un análisis tridimensional influye en el factor de seguridad?

INDICADORES DE MEDICIÓN DEL PROBLEMA

Indicador Unidad de Descripción

Medida

Factor de Adimensional Se medirá el factor de seguridad del talud.

seguridad Se investigó diversos autores que propusieron
métodos que involucraran las condiciones
tridimensionales de los taludes. Según
Bovolenta & Bianchi (2020), para las
condiciones geomorfológicas del talud, se
propone el método de elementos finitos. El
cual, simulará comportamientos físicos
precisos. Estos influirán en el cálculo del
factor de seguridad y en el comportamiento de
desplazamientos y tensiones de las cargas
aplicadas en los suelos. En primera instancia,
establecieron caracterizar ambas capas por el
modelo constitutivo de Mohr-Coulomb; el
cual simplifica el comportamiento real del
suelo.

4
5
 CAPÍTULO 2. ESTADO DEL ARTE

Categoría 1: Características de la geomorfología

- Título 1.1: 3D Finite Element Seepage and Slope Stability Modeling of a Geomorphic
Landform Reclamation
- Se utilizaron métodos 3d para la filtración y estabilidad de taludes. Se utilizaron
softwares como SVFLux y SVSlope. Se necesitó para modelar la respuesta del suelo a
diversas condiciones climáticas y realizar un análisis de filtración en 3D utilizando la
mecánica del suelo saturado y no saturado. Se utilizó método de equilibrio límite para
la estabilidad del talud. (Lira, 2017).

Categoría 2: Ecuación de reducción al corte

- Título 2.1: Upper Bound Analysis of the Stability of 3D Slopes in the Saturated Soft Clay
Subjected to Seismic Effect.
- El autor estudia los problemas de estabilidad tridimensional del talud de arcilla blanda
saturada bajo cargas sísmicas. Se analizó la influencia del ángulo del talud y de las
fuerzas sísmicas horizontales y verticales sobre el coeficiente de estabilidad del talud
saturado de arcilla blanda. Utilizó una fórmula conformada por la resistencia al corte y
la tensión principal. Este método puede calcular con precisión la dirección de
deslizamiento de la superficie de falla y reducir el tedioso trabajo de los métodos
convencionales. Los resultados muestran que el ángulo de la pendiente tiene una gran
influencia en el coeficiente de estabilidad. (Biao et al, 2021).

Categoría 3: Método de Elementos Finitos

- Título 3.1: Stability Analysis of Ecological Slopes Based on a 3D Finite Element Model.
- Este estudio realizó una simulación por elementos finitos del ensayo de corte directo
del compuesto raíz-suelo. Este artículo utiliza el software de elementos finitos
ABAQUS para llevar a cabo un estudio de simulación numérica del suelo, raíces de
hierbas y raíces de arbustos en pendientes ecológicas. Se comprobó que la estabilidad
del talud bajó debido a las intensas lluvias. Las plantas que sirvieron como efecto de
refuerzo influyeron en un factor de seguridad mayor. (Sui et al, 2021)

6
- Título 3.2: Confiabilidad de taludes de suelo no saturado expuestos a precipitaciones
pluviales
- El autor desea calcular la probabilidad de falla y el índice de confiabilidad de taludes
de suelos no saturados expuestos a precipitaciones pluviales. Los cálculos de
infiltración en el talud se realizan con el software Seep mediante el método de
elementos finitos basados en la ecuación de Richards y el análisis de estabilidad
mediante Slope con el método de Morgenstern-Price. Se concluyó que la topografía
tiene un efecto significativo en el FS calculado en un análisis 3D. Los suelos cohesivos
son más sensibles a la topografía de la superficie del suelo. Las geometrías cóncavas
tienen factores de seguridad más sensibles a los cambios en la resistencia cortante.
(Pérez, 2018).

Categoría 4: Modelamiento 3D

- Título 4.1: Geotechnical Analysis and 3D Fem Modeling of Ville San Pietro (Italy)
- El modelo numérico tridimensional de Ville San Pietro es útil para el estudio del sitio,
permitiendo la identificación de las áreas sujetas a movimiento. El modelo fue
calibrado realizando varios análisis. El modelo de elementos finitos actual reproduce la
estratigrafía del sitio y tiene en cuenta la iglesia, aunque de forma simplificada. La
comparación entre los resultados del análisis numérico y los datos experimentales
muestra una buena concordancia. De hecho, los desplazamientos de suelo calculados
concuerdan con los del inclinómetro. Además, las zonas críticas simuladas
corresponden a las áreas donde se detectaron desplazamientos de taludes y donde se
aprecian daños en edificios e infraestructuras. (Bovolenta y Bianchi, 2020).
- Título 4.2: Stability Analysis of Three-Dimensional Slopes Considering the Earthquake Force
Direction
- Se propone el método 3DSpencer simplificado para evaluar la estabilidad de taludes
bajo la fuerza de un terremoto y se utiliza un algoritmo de búsqueda aleatoria para
determinar la superficie de deslizamiento crítica. Se escribió un código de computadora
en Matlab para determinar la superficie crítica de deslizamiento y calcular su factor de
seguridad y dirección de deslizamiento (DOS). (Wan et al, 2018).
- Título 4.3: 3D Finite Element Seepage and Slope Stability Modeling of a Geomorphic
Landform Reclamation

7
- Se utilizaron métodos 3d para la filtración y estabilidad de taludes. Se utilizaron
softwares como SVFLux y SVSlope. Se necesitó para modelar la respuesta del suelo a
diversas condiciones climáticas y realizar un análisis de filtración en 3D utilizando la
mecánica del suelo saturado y no saturado. Se utilizó método de equilibrio límite para
la estabilidad del talud. (Lira, 2017).

Categoría 5: Comparativa con el método bidimensional

- Título 5.1: 3D stability analyses of Skjeggestad landslide

- Se modela el análisis de equilibrio límite de la sección transversal bidimensional (2D)
central. La geometría del deslizamiento de tierra de Skjeggestad no estaba cerca de un
caso 2D, y el modelado tridimensional (3D) es más apropiado para analizar la
pendiente. El documento calcula la estabilidad del talud que falló y compara los
resultados de los análisis de elementos finitos en 3D con el equilibrio límite clásico y
los análisis de elementos finitos en 2D. Los análisis se realizaron en el software
Novapoint GeoSuite Stability. . El trabajo demuestra que para el deslizamiento de
tierra, los efectos 3D fueron significativos (hasta un 40% de aumento en el factor de
seguridad de la sección transversal 2D más crítica). (Jostad et al, 2021).
- Título 5.2: Comparison of limit equilibrium and 2-D, 3-D Finite Element Slope
Stability models: A case study on the slope in Akpinar district, in Ankara
- El autor narra sobre el deslizamiento de tierra ocurrido en el distrito de Akpinar en
Ankara. Debido a ello, realiza una comparación de resultados de métodos
bidimensionales de equilibrio límite y elementos finitos, junto con el método de
elementos finitos en 3D. Analiza los datos de estudios previos y evalúa el rendimiento
de la pendiente de diferentes métodos para comparar los análisis estáticos anteriores.
Lo que da como resultado un FS de 1, 0.89 y 1.11 respectivamente. Se concluyó que
gracias al modelamiento 3D, se demostró la ocurrencia del deslizamiento superficial y
un factor de seguridad mayor a los métodos bidimensionales, contrario a expectativas.
(Can, 2019).
- Título 5.3: 3D-effects in Slope Stability Analysis: Sensitivity of Geometry in 2D and 3D
Stability Programs
- Se utilizaron elementos triangulares con 15 nodos para construir mallas en análisis 2D
suponiendo condiciones de deformación plana, mientras que las cuñas de 15 nodos se

8
 utilizaron para construir mallas en análisis 3D. En el análisis 3D se utilizó en todos los
casos un ancho de 100m. Se usó PLAXIS 3D, ESTABL5M UTEXAS3 entre otros. El
FS mayor fue el de PLAXIS 3D. En la práctica, el análisis bidimensional de pendientes
es más aceptado ya que menor los factores de seguridad resultan de ignorar los efectos
secundarios de las masas deslizantes 3D en el método de equilibrio límite. (Solli, 2020)
- Título 5.4: Contribution to the study of the effect of three-dimensional stress states to
slope stability analysis
- Se utiliza el método de Kulhawy. Las condiciones geométricas de campo en la mayoría
de los problemas de ingeniería se representan mejor mediante modelos 3D. Las
configuraciones geométricas en general suelen estar asociadas con un papel importante
en el análisis de estabilidad de taludes debido a la distribución de tensiones. Desde el
punto de vista de la estabilidad, las pendientes convexas se ven más beneficiadas que
las cóncavas con las mismas condiciones. (Rodrigues, 2018).

9
 CAPÍTULO 3. METODOLOGÍA

3.1 OBJETIVOS

Objetivo General

- Desarrollar un análisis tridimensional de estabilidad de taludes para la Costa Verde

en el distrito de Miraflores mediante un software basado en el método de elementos
finitos para representar de una manera más realista los cambios geomorfológicos
que se dan en un talud del acantilado Miraflores-Costa Verde.

Objetivos Específicos

- Identificar las propiedades mecánicas y físicas del suelo del talud.

- Recopilar los datos del terreno para determinar las condiciones topográficas del
lugar a ser investigado.
- Utilización del software para simular un código numérico de elementos finitos para
el talud.
- Explicar cómo la variación geomorfológica del terreno afecta la estabilidad del
talud en la Costa Verde.
- Evaluar el análisis bidimensional y tridimensional para luego realizar una
comparación de los resultados.

3.2 PROCEDIMIENTO

3.2.1. Nivel y Diseño de la Investigación

Nivel de Investigación: Explicativo

El nivel de investigación es Explicativo; ya que se buscará determinar las causas y efectos

de la variación geomorfológica en el talud y cómo está influyendo en la estabilidad de este.

Diseño de Investigación: Descriptivo

10
El diseño de investigación es Descriptivo; porque realizaremos una evaluación del método
tridimensional para estabilidad de taludes y el efecto de la variación geomorfológica que
modificará la pendiente y geometría general del talud. Además, se realizará la comparación
con el método bidimensional y las diferencias y aporte que trajo consigo el análisis 3D.

3.2.2. Procedimiento de la Investigación

El procedimiento de la investigación empezará con la utilización del modelo de la

topografía de un talud estudiado en investigaciones anteriores de la Costa Verde.

Se utilizará datos de investigaciones anteriores; en la cual, para este estudio la falla

potencial será esférica. Datos como los parámetros de densidad y de resistencia del suelo,
la superficie freática, la resistencia al corte, los pesos unitarios del suelo y el material de
relleno; serán datos de entrada para el modelamiento 3D en el software.

Se procederá a analizar el talud con un software de modelamiento 3D mediante elementos

finitos.

Se hará una evaluación de la ubicación del nivel freático y gradiente hidráulico promedio.

Solución de los sistemas de ecuaciones de reducción de resistencia al corte y cálculo de los

factores de seguridad.

Comparar los resultados con el análisis bidimensional e identificar los valores que influyen
en el factor de seguridad.

11
 Figura 4. Flujograma del procedimiento

Fuente: Propia

3.3 JUSTIFICACIÓN DEL APORTE DE LA INVESTIGACIÓN

Durante años, se han utilizado, convencionalmente, métodos bidimensionales para la

estabilidad de taludes; sin embargo, no precisan correctamente el comportamiento de los
cambios geomorfológicos del suelo del talud y cómo varían en el terreno. Debido a que las
condiciones de falla del talud son de condición tridimensional; y los métodos
bidimensionales son una simplificación de la situación real. Por ello, se propone la
utilización del método tridimensional para simular de manera práctica las áreas críticas
calculadas por el código numérico con las porciones reales afectadas por los
desplazamientos y daños posibles del suelo relacionados con el nivel freático. Además, el
modelo 3D podrá ser de utilidad para futuras investigaciones geotécnicas del talud de la
Costa Verde.

12
 CAPÍTULO 4. MARCO TEÓRICO

El suelo:

El suelo es la porción más superficial de la corteza terrestre, constituida en su mayoría por

residuos de roca provenientes de procesos erosivos y otras alteraciones físicas y químicas,
así como de materia orgánica fruto de la actividad biológica que se desarrolla en la
superficie.

Propiedades físicas del suelo:

• La textura es la que determina la proporción en la que se encuentran las partículas

minerales de diversos tamaños que hay presentes en el suelo.
• La estructura es la forma en la que las partículas del suelo se unen para formar
agregados. La densidad influye en la distribución de la vegetación.
• Suelos más densos son capaces de sustentar más cantidad de vegetación.
• La temperatura también influye en la distribución de la vegetación, sobre todo en
altitud.
• El color depende de sus componentes y varía con la cantidad de humedad presente
en el suelo.

Gravedad específica

Gravedad especifica es un tipo particular de densidad relativa, particularmente en mecánica

de suelos es la densidad del suelo entre la densidad del agua. Este cociente nos da un
valora dimensional que nos sirve para clasificar en un tipo específico nuestra muestra de
suelo. (WordPress, 2022)

Contenido de agua

El contenido de humedad del suelo es un valor que determina la cantidad de agua en una
cierta cantidad conocida de suelo; puede expresarse como porcentaje, agua por peso o
volumen, o pulgadas de agua por pie de suelo. (EOS, 2022)

Relación de vacíos

13
Representa la proporción de vacíos de un suelo. Se define como la relación entre el
volumen de sólidos y el volumen de vacíos de un suelo. (INACAP, 2022)

Densidad

Se refiere al peso por volumen del suelo. Existen dos tipos de densidad, real y aparente. La
densidad real, de las partículas densas del suelo, varía con la proporción de elementos
constituyendo el suelo y en general está alrededor de 2,65. Una densidad aparente alta
indica un suelo compacto o tenor elevado de partículas granulares como la arena. Una
densidad aparente baja no indica necesariamente un ambiente favorecido para el
crecimiento de las plantas. (FAO, 2022)

Indice de plasticidad

El índice de plasticidad se expresa con el porcentaje del peso en seco de la muestra de

suelo, e indica el tamaño del intervalo de variación del contenido de humedad con el cual
el suelo se mantiene plástico. (FAO, 2022)

Propiedades mecánicas:

Cohesión

Existen suelos finos como los arcillosos y limosos donde existen esfuerzos friccionales.
Los suelos arcillosos pueden soportar sobrecargas a través del tiempo de su vida geológica.
Estos esfuerzos que actúan en las partículas del suelo provocan una resistencia adicional a
la que llamaremos cohesión. (Leoni, 2015).

La cohesión es la medida de la adherencia que existe entre las partículas del suelo. Es la
representación de la resistencia cortante debido a cementación. (Quispe, 2021).

Ángulo de fricción interna

El ángulo de fricción es la representación matemática del coeficiente de rozamiento (𝝓), el

cual es concepto básico de la física. El ángulo de fricción interna depende de la
compacidad del material, de la forma de granos, grosor y la granulometría. Según el 2do
Congreso Internacional de Rotterdam, se precisan órdenes de magnitud de la influencia de
estos factores. (Costet, 1975).

14
 Figura 5: Valores del ángulo de fricción interna para determinadas situaciones.

Resistencia al corte

Cuando se trata de analizar la estabilidad de un determinado talud constituido por suelos

naturales o artificiales es necesario valorar la resistencia al esfuerzo cortante de éstos, y el
estado tensional en el interior terreno, ya que inevitablemente se producen esfuerzos
cortantes.(IGME, 1986).

Los procedimientos más habituales para determinar la resistencia al corte de los suelos
consisten en efectuar ensayos en laboratorios con muestras, remoldeadas o inalteradas, que
reproduzcan el estado en que el material se encuentra en el terreno.(IGME, 1986)

Nivel freático:

Una baja permeabilidad de los estratos superiores implicará una alta presión porosa de
agua, produciendo a la vez un decrecimiento de los esfuerzos efectivos entre los granos
sólidos. La poca profundidad del nivel freático. la alta presión porosa del agua y la
granulometría de los estratos superiores. pueden producir un flujo hacia arriba. el cual
reducirá los esfuerzos efectivos de los estratos superiores. en especial si se trata de arena.
(INGEMMET, 1997).

Zonificación sísmica:

Según Cañari, el área que abarcó el informe de EMAPE por HIDROENERGIA consultores
en ingeniería SCRL fue la zona desde el distrito de San Miguel (bajada Bertolotto) hasta la
quebrada de Armendáriz (Miraflores), con una longitud aproximada de 8.5 km y un ancho
aproximado de 55 m. Resalta, entre otros aspectos, el ensayo in-situ de corte directo y la

15
Zonificación Geológico-Geotécnica de los acantilados en términos de riesgo. Es uno de los
pocos trabajos que logra encontrar valores de parámetros de la resistencia al corte del
conglomerado de la Costa Verde, mediante ensayos de campo a gran escala. Utiliza los
valores obtenidos para calcular factores de seguridad en la estabilidad de taludes.

En el trabajo de Cañari (2001) podemos apreciar una zonificación geológicageotécnica de

los acantilados en términos de riesgo. Según ésta, el tramo tiene un gran sector catalogado
como de riesgo muy elevado o crítico (rojo), y en menor proporcion, zonas de riesgo
medio (azul) y riesgo elevado (naranja).

Figura 6. Zonificación geológica-geotécnica de los acantilados en términos de riesgo.

Coeficiente sísmico:

Debido a la condición de borde que representa un talud durante un sismo, a la subjetividad

en que se basan muchos informes acerca de la “intensidad” de daños ocurridos y a las
hipótesis que se plantean para el desarrollo de un análisis pseudo estático, se selecciona el
sismo ocurrido el año 1966, sobre otros sismos de gran envergadura como los de 1970 y
1974, pues este presenta la máxima “aceleración máxima” registrada en el conglomerado
de Lima, que es aprox. igual a 0.275 g.

Datos del registro: LIMA, PERU, INSTITUTO GEOFÍSICO DEL PERÚ, 10/17/66,
COMP N08E 3282 POINTS OF A(T) IN CM/SEC/SEC/1000 AT .02 SEC INTERVAL
PEAK VALUES0 ACCEL=-269.3 CM/SEC/SEC, VELOCITY= -21.6 CM/SEC, DISPL=
-16.6 CM. (Díaz, 2008).

Concepto y falla de un talud

Talud: La definición de talud corresponde a una superficie o plano inclinado, se puede

encontrar de manera natural o como una formación antrópica.También se define talud
como la acumulación de materiales rocosos y suelos acumulados a los pies de una
pendiente. Suelen ser estructuras compuestas del mismo material presente en el suelo, roca,

16
concreto armado u otro que pueda contener la presión ocasionada por el suelo. En nuestra
empresa prestamos el mejor servicio en estudios de estabilidad de taludes.

Partes del talud:

Figura 7. Partes de un talud

Fuente: Curso de ingeniería geotecnía UPC

Circulo de falla

El circulo de falla es un estudio físico matemático que se desarrolla en el área de la

geotecnia, simulando el plano de falla que se origina con los deslizamientos en taludes, el
círculo se proyecta con un radio imaginario en el cual se definen los esfuerzos de
resistencia y de cortante.

Los esfuerzos resistentes dependen de las propiedades mecánicas del suelo (cohesión
existente entre sus partículas y su grado de fricción) y como su nombre lo indica, son los
que resisten los esfuerzos generados por el peso de la masa de suelo que se desplaza y las
cargas adicionales que se trasmiten sobre el talud.

Tipos de suelos en taludes:

• No Cohesivo C=0
• Cohesivo C≠0

17
Falla en un talud: Ocurre como un deslizamiento de la masa de suelo, actuando como un
sólido de cuerpo rígido que se desliza a lo largo de la falla, afectando mucho la estabilidad
del talud en gran medida.

Elementos de un deslizamiento:

Figura 8. Elementos de un deslizamiento

Fuente: Curso de ingeniería geotecnía UPC

Métodos bidimensionales:

Estabilidad de taludes 2D

La estabilidad de taludes se calcula tradicionalmente con métodos bidimensionales. A

menudo, se utiliza el método de equilibrio límite que consiste en dividir la resistencia al
corte ante las fuerzas impulsoras del peso de la masa del suelo y otras cargas externas, las
cuales generan una superficie de falla crítica. (Jostad et al, 2021)

18
Para los análisis bidimensionales, se utlizan los parámetros de resistencia, tipo de material,
condiciones de contorno, etc; y posteriormente se realiza el cálculo de la superficie de
deslizamiento potencial mediante una sección transversal. Luego, se obtienen los rangos de
FS de estabilidad. (Briceño et al, 2021).

Figura 9 : Superficie de falla crítica en un análisis de equilibrio límite 2D.

No obstante, los análisis 2D no simulan la complejidad del comportamiento de los

deslizamientos de tierra como se realizaría en un análisis 3D. Esto se debe a que los
deslizamientos pueden fallar hacia diferentes direcciones y velocidades. (Bovolenta y
Bianchi, 2020).

Método del equilibrio límite (LEM)

El método del equilibrio límite consiste en estudiar el equilibrio de un cuerpo rígido,
constituido por el talud y por una superficie de deslizamiento de cualquier forma (línea
recta, arco circular, espiral logarítmica). Con tal equilibrio se calculan las tensiones de
corte (τ) y se comparan con la resistencia disponible (τf), calculada según el criterio de
rotura de Coulomb: De tal comparación deriva la primera indicación de estabilidad, con el
coeficiente de seguridad:

F=τf /τ

Entre los métodos del equilibrio último hay algunos que consideran el equilibrio global del
cuerpo rígido (Culman) mientras que otros, por falta de homogeneidad, dividen el cuerpo
en rebanadas y consideran el equilibrio de cada una (Fellenius, Bishop, Janbu, etc.).

19
 Figura 10. Representación del método LEM

Método de las rebanadas

La masa susceptible al deslizamiento se subdivide en un número conveniente de rebanadas.
Si el número de rebanadas es igual a n, el problema presenta las siguientes incógnitas:

- n valores de las fuerzas normales Ni en la base de cada rebanada;

- n valores de las fuerzas de corte en la base de la rebanada Ti;
- (n-1) fuerzas normales Ei en la conexión de las rebanadas;
- (n-1) fuerzas tangenciales Xi en la conexión de las rebanadas;
- n valores de la coordenada del punto de aplicación de las Ei;
- (n-1) valores de la coordenada del punto de aplicación de las Xi;
- una incógnita constituida por el factor de seguridad F.

En total las incógnitas son (6n-2).

Mientras las ecuaciones a disposición son:

- ecuaciones de equilibrio de momentos n;

- ecuaciones de equilibrio en la traslación vertical n;
- ecuaciones de equilibrio en la traslación horizontal n;
- ecuaciones del criterio de rotura n.

Total número de ecuaciones 4n.

20
El problema es estáticamente indeterminado y el grado de indeterminación es igual a:

i=(6n-2)-4n=2n-2

El grado de indeterminación se reduce a (n-2). Al asumir que Ni se aplica en el punto

medio de la franja, esto equivale a crear la hipótesis de que las tensiones normales totales
están distribuidas uniformemente. Los diferentes métodos que se basan en la teoría del
equilibrio límite se diferencian por el modo en que se eliminan las (n-2) indeterminaciones.

Figura 11. Representación del método de las rebanadas

Método de Fellenius (1927)

Con este método (válido solo para superficies de deslizamiento circulares) se pasan por
alto las fuerzas entre las franjas, por lo tanto las incógnitas se reducen a:

- n valores de las fuerzas normales Ni;

- n valores de las fuerzas de corte Ti;
- 1 factor de seguridad.

Incógnitas (2n+1). Las ecuaciones disponibles son:

- n ecuaciones de equilibrio traslación vertical;

21
- n ecuaciones del criterio de rotura;
- ecuación de equilibrio de momentos global.

F={Σì[ci·li+(Wi·cosαi-ui·li)·tanφi}/(Σì·sinαi)

Esta ecuación es fácil de resolver pero se ha visto que da resultados conservadores

(factores de seguridad bajos) especialmente para superficies profundas.

Figura 12. Representación del método de Fellenius

Método de Janbu (1967)

Janbu extendió el método de Bishop a superficies de deslizamiento de cualquier forma.
Cuando se tratan superficies de deslizamiento de cualquier forma el brazo de las fuerzas
cambia (en el caso de las superficies circulares queda constante e igual al radio), por este
motivo es mejor valorar la ecuación del momento respecto al ángulo de cada bloque.

F={Σì[ci·bi+(Wi-ui·bi+ΔXi)·tanφi]·[sec^(2)αi/(1+tanαi·tanφi/F)]}/(ΣìWi·tanαi)

22
 Figura 13. Representación del método de Janbu

Acciones en la i-ésima rebanada según las hipótesis de Janbu y representación de la

totalidad de la masa

Asumiendo ΔXi = 0 se obtiene el método ordinario. Janbu propuso además un método para
la corrección del factor de seguridad obtenido con el método ordinario según lo siguiente:

Fcorretto=f0·F

23
Donde: f0 depende de la geometría y de los parámetros geotécnicos y esto se puede
encontrar en tablas y gráficos. Esta corrección es muy confiable para taludes poco
inclinados.

Figura 14. Gráfica coeficiente vs profundidad

Método de Bishop (1955)

Con este método se toman en cuenta todas las fuerzas actuantes en los bloques. Fue el
primero en describir los problemas relacionados con los métodos convencionales.
Las ecuaciones usadas para resolver el problema son:

Σ Fy=0 Σ M0=0 Criterio de rotura

F={Σì[ci·bi+(Wi-ui·bi+ΔXi)·tanφi]·[secαi/(1+tanαi·tanφi/F)]}/(ΣìWi·sinαi)

Los valores de F y de ΔX que satisfacen esta ecuación dan una solución rigurosa al
problema. Como primer aproximación conviene plantear ΔX = 0 e iterar para el cálculo del

24
factor de seguridad. Este procedimiento se conoce como método de Bishop ordinario y los
errores con respecto al método completo son de alrededor de un 1 %.

Figura 15. Representación del método de Bishop

Estabilidad de taludes 3D

Mediante la realización de la estabilidad tridimensional del talud, podremos evaluar la

probabilidad de falla del talud en función de los parámetros geomecánicos del suelo. En el
método 3D se emplean diferentes estratigrafías y diferentes direcciones de movimiento y
se comparan los resultados con análisis bidimensional. El método tridimensional puede no
ser estrictamente menos conservador que el método de comparación bidimensional, sino
que respondería a una mayor precisión conforme a las distintas configuraciones
geométricas y geológicas de los taludes evaluados. Entre los métdos más usados, tenemos
el método del equilibrio límite y los métodos numéricos que permiten solucionar las
ecuaciones de la mecánica de sólidos. Estos últimos tienen tres versiones principales: los
elementos finitos, los elementos de frontera y el método del punto material. (Ocando,
2020).

25
En la práctica, el análisis 3D de la estabilidad de taludes no se realiza a menos que la
geometría de la pendiente es muy complicado o el mecanismo de falla es complejo.
(Albataineh, 2006).

Además, se requiere de la utilización de otros programas, los cuales están, en primer lugar,
relacionados a la cartografía estructural, toma de datos geotécnicos para el mapeo. (Ramos
Vásquez, 2017).

Uso de software

El software Plaxis utiliza el MEF para el cálculo por simulación de deformaciones,

esfuerzos y estabilidad, utilizando una interfaz gráfica (CAD) para generar los modelos
geométricos. Tiene como limitación la modelación en secciones transversales uniformes
con deformaciones planas donde la carga y los estados tensionales se supongan uniformes,
realizando un análisis paramétrico. (Briceño et al, 2021).

Plaxis 3D es un código de elementos finitos tridimensional utilizado para diseños

geotécnicos.; con el fin de estudiar el comportamiento de taludes e identificar valores de
nodos para el monitoreo de ubicación. (Bovolenta y Bianchi, 2020).

Análisis pseudoestático en taludes

El método utiliza el mismo procedimiento general de cualesquiera de los métodos de

equilibrio límite, con la diferencia de que se incluyen fuerzas pseudoestáticas horizontales
y verticales, debidas al evento sísmico. Estas fuerzas sísmicas se asumen, que son
proporcionales al peso de la masa de deslizamiento potencial y de los coeficientes sísmicos
kh y kv expresada en términos de veces la aceleración g, producida por el sismo, (Aguilar,
2005).

26
 Figura 16. Representación del análisis pseudoestático en taludes

Método de elementos finitos

El método de los elementos finitos es una herramienta computacional muy potente en

ingeniería. Adquiere su poder de la capacidad de simular comportamientos físicos usando
herramientas computacionales sin la necesidad de simplificar el problema, obteniéndose
resultados más precisos y confiables. El método de los elementos finitos puede usarse para
estudiar la estabilidad de taludes usando una definición de falla similar a la de los métodos
de equilibrio límite, éstos proponen en principio una superficie de deslizamiento para
lugeon examinar el valor del coeficiente de seguridad de la misma, el cual se define como
la relación entre la resistencia al corte disponible y la resistencia al corte movilizadora a lo
largo de la superficie. ( Beneyto et al, 2022)

Factor de seguridad

El Factor de Seguridad (F) expresa la relación entre las fuerzas resistentes del terreno y las
inestabilizadoras. Valores de F superiores a 1 indican condiciones estables, y valores
menores de 1 indican condiciones inestables. Cuando F = 1 la ladera o el talud se encuentra
en el punto de rotura.

El factor de seguridad de un talud puede calcularse con un modelo numérico, reduciendo

progresivamente la resistencia al corte del suelo hasta la rotura. Así que, el FS estará

27
definido por el ratio entre la resistencia al corte original y la resistencia al corte para la
rotura.(Galera y Velasco, 2002)

El FS que se usa como parámetro de estabilidad, Plaxis lo evalúa utilizando el método que
se conoce como modelo de elementos finitos, utilizando el criterio de Mohr-Coulomb se
reducen los parámetros de resistencia al cortante de forma gradual, cohesión (C) y ángulo
de fricción interna del suelo (ϕ), provocando la rotura inicial y finalmente la rotura global.
(Briceño et al, 2021)

Figura: Coeficientes de seguridad a emplear en el análisis de estabilidad de taludes

Figura 17. Esquema de estabilidad del talud

REFERENCIAS BIBLIOGRÁFICAS

Albataineh, N. (2006). Slope Stability Analysis Using 2D and 3D methods. University of Akron.

Angelino Calderón, W. L. (2020). Reducción de parámetros de resistencia al corte para la

evaluación de estabilidad de talud de una presa de relave. Universidad San Ignacio de
Loyola.

28
Biao, Z., Yi, J., Hao, C. and Ze, L. (2021). Upper Bound Analysis of the Stability of 3D Slopes in
the Saturated Soft Clay Subjected to Seismic Effect. Frontiers in Earth Science.

Bovolenta, R. and Bianchi, D. (2020). Geotechnical Analysis and 3D Fem Modeling of Ville San
Pietro (Italy). Geosciences.

Bowen Xu (2022) An Analysis of Slope Stability Based on Finite Element Method and Distinct
Element Method, Journal of Physics: Conference Series

Burgess Jesse (2016) Probabilistic seismic slope stability design. Universidad Dalhousie

Can Etiz, M. (2019). Comparison of limit equilibrium and 2-D, 3-D Finite Element Slope Stability
models: A case study on the slope in Akpinar district, in Ankara. Middle East Technical
University.

Cooper Stephen (2017) Estudio de la estabilidad de taludes en una mina operativa excavada en
suelo duro/roca blanda, Universidad Politécnica de Madrid

Fredlund Murray D., Lu Haihua, Fredlund Delwyn D. (2013). Comparison of 3D Finite Element
Slope Stability WIith 3D Limit Equilibrium Analysis. International Conference on Soil
Mechanics and Geotechnical Engineering.

Galera, J. M. & Velasco, P. (2002). Análisis de estabilidad de taludes. Comparativa entre los
métodos de equilibrio límite y los análisis tensodeformacionales. Ingeopres: Actualidad
técnica de ingeniería civil, minería, geología y medio ambiente.

Haihua Lu, Delwyn F., Murray F. (2013) Comparison of 3D Finite Element Slope Stability With
3D Limit Equilibrium Analysis, Proceedings of 18 International Conference on Soil
th

Mechanics and Geotechnical Engineering, Paris 2013

Jostad, H. P., Sivasithamparam N., Lacasse, S., Degago, S. A., Le, T. M. H., Giese, S., Rosenquist
af Åkershult, A., Johansen, T. and Aabøe5, R. (2021). 3D stability analyses of Skjeggestad
landslide. IOP Conference Series: Earth and Evironmental Science.

Lira Santos, I. (2017). 3D Finite Element Seepage and Slope Stability Modeling of a Geomorphic
Landform Reclamation. West Virginia University.

Nusantara D., Putri Y., Sari P. y otros (2020) The Comparison Between 2D and 3D Slope Stability
Analysis Based on Reinforcement Requirements International, Journal on Advanced Science
Engineering Information Technology

Nusantara D.A.D, Sari P.T.K., Putri Y.E., Savitri Y.R., Amalia A.R.., Margini Nastasia F. (2020).
The Comparison Between 2-D and 3-D Slope Stability Analysis Based on Reinforcement

29
 Requirements. International Journal on Advanced Science Engineering Information
Technology .

Ocando Sánchez Darío Alfonso (2020) Universidad Nacional de Colombia, Análisis tridimensional
de estabilidad de taludes y comparación con el análisis bidimensional

Pérez Fajardo, J. S. (2018). Confiabilidad de taludes de suelo no saturado expuestos a

precipitaciones pluviales. Universidad Autónoma del Estado de México.

Ramos Vásquez, A. A. (2017). Análisis de estabilidad de taludes en rocas. Simulación con LS-
DYNA y comparación con Slide.. Universidad Politécnica de Madrid.

Ricci M. (2018) Caracterización Geotécnica y análisis de estabilidad de “Proyecto Continuidad

Mina Gabriela”, Facultad de Ingeniería Escuela de Ingeniería Civil de Minas

Rodrigues Moura, N. (2018). Contribution to the study of the effect of three-dimensional stress
states to slope stability analysis. Universidade Federal de Goiás.

Roohollah K., Mostafa J., Ahmad S. (2018). The Necessity of Employing Three Dimensional
Concepts in Slope Stability Analysis, J8th Int.Conf. on Geotechnique, Construction Materials
and Environment Kuala Lumpur, Malaysia, Nov. 20-22, 2018, ISBN: 978-4-909106991
C3051

Solli, V. A. (2020). 3D-effects in Slope Stability Analysis: Sensitivity of Geometry in 2D and 3D

Stability Programs. Norwegian University of Science and Technology.

Sui, Z., Yuan, W., Yi, W. and Yang, W. (2021). Stability Analysis of Ecological Slopes Based on a
3D Finite Element Model. Advances in Materials Science and Engineering.

Sundari W., Krisnasiwi I.F. (2017). Analysis of slope stability based in the spencer method on the
ring road section, Sikumana.Kupang City. Journal of Physics: Conference Series.

Ttito, D. (2018). Análisis de estabilidad de taludes de la zona “Los Baluartes” en el Parque

Arqueológico de Saqsayhuamán. Rev Yachay, 7(1), 322-333.

Valeriano N., Tica M., Loaiza C. , Holguin C., Japura P., Llanos C. and Obregon M. (2014).
Análisis probabilístico de la estabilidad de talud de cantera Espinal – Juliaca. Universidad
Nacional del Altiplano – Puno, Perú.

Valiente Sanz, R., Sobrecases, S. & Díaz Orrego, A. (2015). Estabilidad de Taludes: Conceptos
Básicos, Parámetros de Diseño y Métodos de Cálculo. Revista Civilizate PUCP.

Wan, Y., Gao, Y. and Zhang, F. (2018). Stability Analysis of Three-Dimensional Slopes
Considering the Earthquake Force Direction. Mathematical Problems in Engineering.

30
Wang, J., Lui S., Xu B. (2022). An Analysis of slope Stability Based on Finite Element Method and
Distinct Elemente Method. Journal of Physics: Conference Series.

Wang, L., Chen L., Han L., Li H., Zhang W. (2021). Slope stability prediction using ensemble
learning techniques:A case study in Yunyang Countym Chonqguing,China. Journal of Rock
Mechanics and Geotechnical Engineering.

Yingfa L. (2015). Deformation and failure mechanism of slope in three dimensions. Journal of
Rock Mechanics and Geotechnical Engineering.

Zabolotnii (2020) Three-Dimensional Slope Stability Effects in the failure at the Mount Polley
Tailings Storage Facility, Department of Civil and Environmental Engineering University of
Alberta

31