OBRAS GEOTÉCNICAS

CON PLAXIS 2D
Anderson Luis Barreto Goñi
ELEMENTOS FINITOS

• Permiten resolver problemas de mecánica

de suelos y de fluidos, y la interacción entre
estos.
• Análisis de problemas geotécnicos con
elementos estructurales.
• Métodos de elementos finitos, método de
diferencias finitas.
• Permite análisis la distribución de esfuerzos
y deformaciones.
• Adicionalmente puede determinar el factor
de seguridad.
• Su aplicación ha mejorado los estándares
de diseño.
 INTRODUCCIÓN
• El método de elementos finitos
consiste en la discretización de
un cuerpo de estudio en
elementos pequeños que
pueden ser de distintas formas
y tienen nodos definidos, estos
en conjunto forman un
enmallado (conexión entre si)
el cual representa el cuerpo de
estudio con este se calcula
aproximadamente los esfuerzos
y desplazamientos en el
cuerpo.
MÉTODOS DE ELEMENTOS FINITOS APLICADO
A LA GEOTÉCNIA
 FINALIDAD MODELAMIENTO
GEOTÉCNICO

• Comprensión del comportamiento de las estructuras

geotécnicas
• Determinación de los desplazamientos del suelo bajo cargas
de trabajo.
• Cálculo de fuerzas internas para el diseño de estructuras.
• Identificación de posibles mecanismos de falla.
 MODELOS CONSTITUTIVOS
Expresiones matemáticas que relacionan esfuerzos y deformaciones, que
son empleados para modelar el comportamiento de un elemento.
Parámetro Descripción
c Cohesión
• Obedece la ley de Hooke de • Modelo elasto plástico muy
elasticidad lineal isotrópica empleado para primeras φ Ángulo de fricción
aproximaciones
• Incorpora el módulo de Young E y Ángulo de
el coeficiente de Poisson u • Incorpora el módulo de Young Ψ
E, el coeficiente de Poisson u, dilatancia
• Se podría emplear para capas cohesión c, ángulo de fricción E Elasticidad
muy rígidas del suelo y/o para phi, y el ángulo de dilatancia
muy pequeñas deformaciones Modelo elástico Modelo Mohr 𝜈𝜈 Poisson
lineal Coulumb
Parámetro Descripción
E Elasticidad
𝜈𝜈 Poisson

Modelo de suelo
Modelo para con
• Modelo variante elastoplástico del
suelo blando endurecimiento modelo hiperbólico
(Soft Soil Model) (Hardening Soil
• Admite plasticidad de
• Basado en el model Cam Clay Model) endurecimiento por fricción
• Simula el comportamiento de suelos • Incluye endurecimiento por
blandos como arcillas normalmente compresión
consolidadas y turbas • Simula de carga y descarga
• Buenos resultados para compresión • Ideal para arenas, gravas, suelos muy
primaria blandos
IDEALIZACIÓN GEOMÉTRICA
• ESTADO PLANO DE ESFUERZOS
Simplificación de un problema a un plano 2D
(deformación plana)
 Procedimiento

Apuntes del curso Elementos Finitos Aplicados a Geotecnia con el Dr. Ing. Jorge Luis Cárdenas Guillén
Método de elementos finitos
Esta matriz de rigidez es Fuerzas externas {Q}
valida para todo el
𝑒𝑒𝑝𝑝 𝑇𝑇 𝑒𝑒𝑝𝑝 𝑄𝑄 = 𝐾𝐾 𝑞𝑞 Desplazamientos nodales 𝑞𝑞
𝐾𝐾 = � 𝐵𝐵 � 𝐷𝐷 � 𝐵𝐵 � 𝑑𝑑𝑑𝑑 rango elástico y plástico.

𝐾𝐾1 𝑒𝑒
Matriz de rigidez global 𝐾𝐾
Matriz de rigidez de
cada elemento 𝐾𝐾 𝑒𝑒
𝐾𝐾2 𝑒𝑒

𝐾𝐾 𝑒𝑒 𝐾𝐾3 𝑒𝑒 𝐾𝐾6 𝑒𝑒
𝐾𝐾4 𝑒𝑒 𝐾𝐾5 𝑒𝑒

𝐾𝐾7 𝑒𝑒 𝐾𝐾𝑛𝑛 𝑒𝑒
𝐾𝐾9 𝑒𝑒
𝐾𝐾8 𝑒𝑒
𝐾𝐾10 𝑒𝑒

Hecho por Luis Felipe Rojas Gomez

 Método de elementos finitos 4 Derivada de la función de interpolación (B)
𝑑𝑑𝑑𝑑𝑖𝑖
1 Elemento finito 0 0
𝑑𝑑𝑑𝑑
3 Función de interpolación (N) 𝑑𝑑𝑑𝑑𝑖𝑖
0 0
𝜁𝜁 4𝜁𝜁 − 1 4𝜁𝜁 − 2 4𝜁𝜁 − 3
𝑑𝑑𝑑𝑑
𝑁𝑁1 =
6 𝑑𝑑𝑑𝑑𝑖𝑖
𝜉𝜉 4𝜉𝜉 − 1 4𝜉𝜉 − 2 4𝜉𝜉 − 3 0 0
𝑁𝑁2 = Derivando 𝑑𝑑𝑑𝑑
Desplazamientos 6 𝐵𝐵 = 𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑
𝑖𝑖 𝑖𝑖
nodales 𝑞𝑞
𝑁𝑁3 =
𝜂𝜂 4𝜂𝜂 − 1 4𝜂𝜂 − 2 4𝜂𝜂 − 3 𝜕𝜕x, 𝜕𝜕y, 𝜕𝜕z 0
𝒖𝒖 = 𝑁𝑁 � q 6 𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑
𝑁𝑁4 = 4𝜁𝜁𝜁𝜁 (4𝜁𝜁 − 1)(4𝜉𝜉 − 1) 𝑑𝑑𝑑𝑑𝑖𝑖 𝑑𝑑𝑑𝑑𝑖𝑖
𝑁𝑁5 = 4𝜉𝜉𝜉𝜉(4𝜉𝜉 − 1)(4𝜂𝜂 − 1) 0
Desplazamientos 𝑁𝑁6 = 4𝜂𝜂𝜂𝜂(4𝜂𝜂 − 1)(4𝜁𝜁 − 1) 𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑
global {u}
Función de interpolación 𝑁𝑁
𝑁𝑁7 = 𝜉𝜉𝜉𝜉 4𝜁𝜁 − 1 4𝜁𝜁 − 2 ∗ 8/3 𝑑𝑑𝑑𝑑𝑖𝑖 𝑑𝑑𝑑𝑑𝑖𝑖
𝑁𝑁8 = 𝜉𝜉𝜉𝜉 4𝜉𝜉 − 1 4𝜉𝜉 − 2 ∗ 8/3 0
𝑁𝑁9 = 𝜉𝜉𝜉𝜉 4𝜉𝜉 − 1 4𝜉𝜉 − 2 ∗ 8/3 𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑
𝑁𝑁10 = 𝜉𝜉𝜉𝜉 4𝜂𝜂 − 1 4𝜂𝜂 − 2 ∗ 8/3
2 Elemento isoparamétrico 𝑁𝑁11 = 𝜁𝜁𝜁𝜁 4𝜂𝜂 − 1 4𝜂𝜂 − 2 ∗ 8/3 5 Matriz de rigidez (K)
𝑁𝑁12 = 𝜁𝜁𝜁𝜁 4𝜁𝜁 − 1 4𝜁𝜁 − 2 ∗ 8/3
Matriz de rigidez de
𝑁𝑁13 = 𝜂𝜂𝜂𝜂𝜂𝜂 4𝜁𝜁 − 1 cada elemento 𝐾𝐾 𝑒𝑒 𝐾𝐾 𝑒𝑒 = � 𝐵𝐵 𝑇𝑇
� 𝐷𝐷 𝑒𝑒
� 𝐵𝐵 � 𝑑𝑑𝑑𝑑
𝑁𝑁14 = 𝜂𝜂𝜂𝜂𝜂𝜂(4𝜉𝜉 − 1)
𝑁𝑁15 = 𝜂𝜂𝜂𝜂𝜂𝜂(4𝜂𝜂 − 1) Matriz constitutiva 𝐷𝐷
Derivada de la función de interpolación

6 Ecuación global
Fuerzas externas {Q}
𝑄𝑄 = 𝐾𝐾 𝑞𝑞 Desplazamientos nodales 𝑞𝑞
Matriz de rigidez global 𝐾𝐾
Hecho por Luis Felipe Rojas Gomez
INTRODUCCIÓN AL PLAXIS
 ESTRUCTURA DEL PROGRAMA

• Definición del • Definición del • Definición de fases de cálculo

• Definición de la geometría problema proceso • Cálculo
• Generación de la malla (Representación constructivo,
• Definición de condiciones iniciales física y situación estados
inicial) constructivos

Plaxis Plaxis
Input Calculation

Plaxis Plaxis
• Observación de los Curves Output
resultados en los
puntos de • Observación de los
observación resultados del
durante el proceso modelo durante
constructivo los procesos
constructivos

• Inspección de resultados
 COMANDOS

Malla
Elementos físicos Túneles Condiciones de contorno Drenes
y cargas Materiales
 EJEMPLO APLICATIVO 01

Desarrollar el análisis de estabilidad para el siguiente talud

 CONTENIDO

PARÁMETROS SIMBOLOGÍA ARENA 01 ARENA 02 GRAVA 01 UNIDADES

Modelo de material Modelo M.C. M.C. M.C. --
Tipo de comportamiento de
Tipo Drenado Drenado Drenado --
material
Peso unitario de suelo
Y unsat 18 19 21 kN/m3
encima de p.1
Modulo de Young E ref 25,000 35,000 80,000 kN/m2
Relación de Poisson v 0.30 0.30 0.30
Cohesión C ref 5 5 10 kN/m2
Angulo de fricción ϕ 30 32 38 °
Angulo de dilatancia Ψ 0 0 0 °
RGB (COLORES) -- 243,247,24 228,174,37 120,120,104 --
 EJEMPLO APLICATIVO 02 - 03

Desarrollar el análisis de estabilidad para el siguiente talud

Exportar curvas del AutoCad

Hallar el FS Pseudoestatico
 ESFUERZOS INICIALES

K0-procedure
ESFUERZOS INICIALES • Generación de los esfuerzos iniciales
– Se necesita
• Los esfuerzos iniciales representan el estado de Coeficiente de presión lateral del suelo K0.
equilibrio del suelo no disturbado, depende del: – Desventajas
– Peso del suelo No hay equilibrio con superficies inclinadas
– Historial de cargas – Ventajas:
• En Plaxis existen dos posibilidades: No se generan desplazamientos, sólo esfuerzos
– K0 procedure Se puede tener en cuenta el OCR
– Gravity loading
 ESFUERZOS INICIALES

Gravity loading

• El cálculo de los esfuerzos iniciales por el peso del suelo

– Desventajas:
No se crean desplazamientos físicos
Dificultad para relacionarlo con el valor de Ko
Dificultad para introducir y relacionarlo con el OCR
– Ventajas:
Equilibrio satisfactorio en todos los casos

Casos donde gravity loading debe ser usado en lugar de K0-

procedure:
 ELEMENTOS ESTRUCTURALES

Aplicaciones en las que podemos usar placas, anclajes e

interfases

ELEMENTOS ESTRUCTURALES
 Placas / muros

Placas / muros

•3 o 5 nodos en la línea
•3 grados de libertad por nodo
• Placas soportan:
– Fuerzas axiales
– Fuerzas cortantes
– Momentos
• Comportamiento elástico e inelástico
• Para modelar muros, losas, túneles, etc.
• Los parámetros más importantes para representarlos son la
flexión EI y la resistencia axial EA
• El orden de magnitud para un muro de 30cm es: EA de 107 y
EI de 105
 Placas / muros

Placas / muros - Parámetros

 Placas / muros

Placas / muros - Comportamiento con suelo

 Anclajes

Anclajes con fijación (puntales) Anclajes nodo a nodo

Para modelar soportes, puntales y/o anclajes Para modelar anclajes

helicoidales – Comportamiento elasto plástico
– Comportamiento elasto plástico – Conectado en ambos puntos por líneas
– Uno de los puntos es fijo en la geometría y el geométricas
otro es fijo para – No es necesario colocar una línea geométrica
el desplazamiento en su posición
– Posicionamiento para cualquier ángulo – Opción de pretensado
Propiedades de los anclajes
 Anclajes pretensados Geomallas

• Definido durante las etapas constructivas • Con 3 o 5 elementos de línea

• Para tensión (anclajes de cables) o compresión (anclajes de • Comportamiento elasto plástico
barras) • No soportan compresión pero si tracción
Interfase de materiales

Interfase – Factor de reducción

EJEMPLO APLICATIVO 04
 EJEMPLO APLICATIVO 04
Propiedades de los materiales
Parámetros del suelo en base al modelo Mohr – Coulomb
 EJEMPLO APLICATIVO 04

Parámetros de entrada del muro pantalla

Parámetros de entrada de los anclajes

Parámetros de entrada del grout (geotextil)

EJEMPLO APLICATIVO 04
Referencias bibliográficas

Apuntes del curso Taller de Ingeniería Geotécnica con el

Dr. Ing. Zenón Aguilar Bardales

Apuntes del curso Elementos Finitos Aplicados a Geotecnia con el

Dr. Ing. Jorge Luis Cárdenas Guillén

Apuntes del curso Modelos Constitutivos en Materiales Geotécnicos con el

Dr. Ing. Jorge Luis Cárdenas Guillén

Apuntes del curso Presas de Tierra con el

MSc. Ing. Milton Teran Ulloa

Apuntes del curso Dinámica de Suelos con el

MSc. Ing. Denys Parra Murrugarra

Apuntes del Taller Aplicativo - Estabilidad de Taludes del

ECOE
 MUCHAS
GRACIAS !!!

“EL MUNDO ES NUESTRO CAMPO, EL SUELO NUESTRO ELEMENTO”