CURSO DE MODELIZACIÓN

DE MUROS PANTALLA CON

PLAXIS 2D

JAVIER MORENO ROBLES

 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

MODULO 2:

COMPORTAMIENTO DEL TERRENO.

MODELOS CONSTITUTIVOS (I)

Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 1

 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

MÓDULO 2. COMPORTAMIENTO DEL TERRENO. MODELOS CONSTITUTIVOS (I)

2.1. INTRODUCCIÓN A LOS MODELOS CONSTITUTIVOS 3

2.2. MODELO MOHR-COULOMB (MC) 5
2.3. MODELO HARDENING SOIL (HSM) 7
2.4. MODELO HARDENING SOIL EN PEQUEÑAS DEFORMACIONES (HS SMALL) 17
2.5. COMPARATIVA RESULTADOS MC-HSM 20

Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 2

 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

MÓDULO 2. COMPORTAMIENTO DEL TERRENO. MODELOS

CONSTITUTIVOS (I)
2.1. INTRODUCCIÓN A LOS MODELOS CONSTITUTIVOS
En la resolución de un problema de elementos finitos es necesario resolver un sistema de
ecuaciones del tipo F = Kd, en el que a partir de las fuerzas impuestas obtenemos los
movimientos asociados o viceversa.

Para poder realizar el proceso completo de cálculo es necesario contemplar diferentes

subprocesos:

 La obtención de las tensiones a partir de las fuerzas exteriores se obtienen con las
ecuaciones de equilibrio, que son independientes del tipo de material. Así, por ejemplo, si
se conoce la fuerza de tracción (T) que actúa sobre una barra de sección conocida (A), se
obtiene la tensión que solicita la barra con una simple consideración de
equilibrio.

 La obtención de las deformaciones a partir de los movimientos se obtienen con las

ecuaciones de compatibilidad de deformaciones que, como las anteriores, son
independientes del material que se está considerando. Por ejemplo, si la barra
anteriormente considerada tiene una longitud de 1 m y se alarga 1 cm = 0,01 m, la
deformación que ha sufrido ha sido de 0,01/1 = 1%.

Es importante considerar que de la barra anterior no se ha indicado si es de acero, de goma o de

roca, y en todos estos casos la tensión que actúa y la deformación indicada estaría correctamente
calculada.

Para poder cerrar la “cadena” de cálculo entre las fuerzas y los movimientos es necesario
introducir los modelos constitutivos que son los que relacionan las fuerzas y las deformaciones y
estos sí que dependen del material.

Figura 1. Diversos tipos de comportamiento en función del tipo de ecuación constitutiva de los materiales

Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 3

 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

Según se puede ver en la Figura 1 existen numerosos tipos de leyes de comportamiento de los
materiales (elástico, elastoplástico, rigido plástico, con endurecimiento, con reblandecimiento,
viscoelástico, etc) y la calidad de los resultados estará en función de la decisión realizada del tipo
de comportamiento a adoptar.

Los elementos más importantes de cualquier modelo constitutivo se definen a continuación:

 Superficie de rotura: es la que delimita los estados tensionales posibles y no posibles.

Habitualmente se emplea la línea de resistencia intrínseca de Mohr-Coulomb (línea roja
de la Figura 2).

 Superficie (o función) de fluencia: es la que define la separación entre el régimen elástico

y plástico (representada en amarillo en la Figura 2). En la misma se ve que la trayectoria
de tensiones en el punto A estará en régimen elástico en carga, la del punto B en régimen
plástico y la del punto C en régimen elástico en descarga.

Figura 2. Representación esquemática de la superficie de rotura y la superficie de fluencia

 Superficie (o función) de potencial plástico: es la que define cómo varían los incrementos
de deformación plástica (volumétrico y de corte) cuando se alcanza la superficie de
fluencia. Está representada en azul en la Figura 3.

 Regla de flujo plástico: es la expresión matemática que define cuánto varía los
incrementos de deformación plástica (módulo) cuando se alcanza la superficie de
fluencia.

Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 4

 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

Figura 3. Representación esquemática de la superficie de fluencia y de potencial plástico

2.2. MODELO MOHR-COULOMB (MC)

El modelo de Mohr-Coulomb es un modelo de plasticidad perfecta muy habitualmente empleado
para materiales friccionales basado en el criterio de rotura de MC.

Figura 4. Superficie de rotura en el modelo MC en el espacio de tensiones

El comportamiento bajo la superficie de rotura es completamente elástico, con un único valor del
módulo de deformación E’, el coeficiente de Poisson ν o los valores asociados del módulo de
deformación volumétrica K y módulo de deformación de corte G que se relacionan con los
anteriores según las expresiones:

[1]

Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 5

 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

[2]

Figura 5. Esquemas de comportamiento del modelo MC

La función de fluencia se expresa mediante la siguiente formulación:

[3]

y la función de flujo con la siguiente:

[4]

Se puede ver en las expresiones anteriores que la función de fluencia y de flujo plástico son
similares a excepción de la variable empleada. Al ser el ángulo de dilatancia Ψ diferente del
rozamiento interno ϕ se puede decir que es un caso de plasticidad no asociativa, tal y como se ha
comprobado que se comportan lo suelos.

Este modelo constitutivo presenta como principales ventajas su claridad y simplicidad, con una
buena representación de la rotura en condiciones drenadas y es comparable con los resultados
analíticos.

Sin embargo tiene unas muy importantes desventajas como son su naturaleza elástica y lineal con
rigidez constante hasta la rotura, que no distingue entre carga primaria y descarga-recarga y que
en compresión edométrica el comportamiento es elástico, es decir, no “sobreconsolida”.

Este último aspecto se puede ver cuando se simular un ensayo edométrico con el modelo MC
(Figura 6) en el que tras finalizar la carga y descarga, el estado tensional queda igual que antes de
la carga inicial, aspecto este que el suelo real no comparte.

Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 6

 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

Figura 6. Reproducción incorrecta de un ensayo edométrico con un modelo MC.

En definitiva, el modelo MC debe emplearse únicamente en las modelizaciones en las que sólo
nos interese el fenómeno de rotura (estabilidad de taludes, estabilidad de frentes de túneles,
cargas de hundimiento, etc.), pero no es adecuado cuando el aspecto deformacional tenga
importancia en los resultados (diseño de una contención, estudio del sostenimiento de un túnel,
etc.).

Para solventar gran parte de estas desventajas, se muestra a continuación un modelo constitutivo
más avanzado denominado Hardening Soil.

2.3. MODELO HARDENING SOIL (HSM)

2.3.1. Introducción
El modelo Hardening Soil (HSM) fue diseñado por Schanz1 (1998) y Schanz et al2. (1999) para
intentar reproducir algunos de los comportamientos macroscópicos que se había comprobado
que los terrenos exhibían, tales como:

 Variación del índice de huecos e debido a las deformaciones plásticas

 Dependencia de la deformabilidad del suelo con el estado tensional del mismo

 Influencia de la historia tensional de los suelos en fenómenos de preconsolidación

 Aparición de deformaciones volumétricas plásticas negativas durante la aplicación de

tensiones de corte

 Aparición de deformaciones irreversibles al alcanzar el criterio de rotura

1
Schanz, T. (1998). Zur modellierung des mechanischen verhaltens von reinbungsmaterialien. Mitt Inst fur
Geotechnik, Universitat Stuttgart, 1998, 45
2
Schanz, T., Vermeer, P., and Bonier, P., editors (1999). Formulation and verication of the Hardening Soil
model. Beyond 2000 in Computational Geotechnics. Balkema, Rotterdam.
Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 7
 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

Figura 7. Superficies de fluencia en el modelo HS en el espacio de tensiones

Este modelo avanzado, sin embargo, no puede reproducir fenómenos de pérdida de resistencia
debido a fenómenos de remoldeo (típico en arcillas sensibles) o rotura de enlaces cementados
entre partículas (típico de limos yesíferos, por ejemplo).

Como complemento del modelo HSM, según se indicará en el apartado 2.4., se ha desarrollado
un “parche” adicional para mejorar el comportamiento a pequeñas deformaciones observado en
los suelos e introducir un comportamiento histerético necesario en los cálculos dinámicos. Este
“modelo” se denomina HS en pequeñas deformaciones.

Figura 8. Esquemas de comportamiento del modelo HS

Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 8

 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

2.3.2. Reproducción del fenómeno de tensiones de cortante

[Link]. Superficie de fluencia a cortante

La superficie de fluencia a cortante se denomina f1 y está implementada bajo el concepto de la

relación hiperbólica tensión de corte q y deformación total vertical ε1, típica en un ensayo triaxial
(Figura 9).

Matemáticamente la expresión que tiene f1 es la siguiente

[5]

Siendo γPS el parámetro que gobierna el endurecimiento plástico y qa el valor asintótico (no de
rotura que es qf) de la curva del ensayo que no se obtiene del mismo, sino que es un concepto
matemático que cumple la siguiente expresión

Figura 9. Relación hiperbólica q-ε1. Tomado del Material Manual de Plaxis.

[6]

El valor de la relación de tensiones en rotura (Failure ratio) se adopta por defecto como Rf = 0,9, y
para la mayoría de los suelos se sitúa entre 0,75 y 1.

A medida que el γPS va incrementándose, las superficies de fluencia f1 van aproximándose a la

línea de rotura de Mohr-Coulomb y perdiendo su curvatura.

Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 9

 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

PS
Figura 10. Evolución de la superficie de fluencia f1 a medida que aumenta la deformación plástica de corte γ .
Tomado del Material Manual de Plaxis.

Para los valores más elevados de γPS, es decir, cuando el suelo está próximo a rotura, el valor de qf
se obtiene del criterio de rotura de Mohr-Coulomb según la siguiente expresión

[7]

[Link]. Dependencia de la deformabilidad del suelo y la tensión de confinamiento

El módulo secante correspondiente al 50% de la resistencia de rotura qf se emplea para definir el

valor de la deformabilidad3 del terreno (Figura 9) y se denomina E50. Sin embargo, el HSM permite
definir una dependencia de este módulo con la presión de confinamiento σ3 según la siguiente
expresión:

( ) [8]

Siendo el módulo de deformación secante al 50% de la carga de rotura qf cuando la tensión

de confinamiento es σ3 = pref. EL valor por defecto es pref = 100 kPa.

El exponente m varía entre 0,3 y 1 para los suelos naturales. Para materiales granulares los
valores habituales están en el entorno de m = 0,5, mientras que para los arcillosos y cohesivos m
se aproxima más a m = 1. En el caso de que m = 0, la deformabilidad del material será
independiente de la presión de confinamiento.

En las siguientes gráficas se muestra la variación del módulo de deformabilidad E50 para un suelo
de peso específico γ = 20 kN/m3 y K0 = 0,5, cuando se varía el exponente m, el ángulo de
rozamiento interno ϕ’ y la cohesión c’. Se observa que se produce un pivote en torno al valor del

3
Es importante entender que E 50 no es un módulo de elasticidad, sino de deformabilidad, ya que desde el
inicio de los incrementos de carga se producen tanto deformaciones elásticas como plásticas. Por este
motivo no se puede relacionar E50 con G50 por las formulaciones habituales de la teoría de la elasticidad.
Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 10
 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

módulo a 10 m = 100 kPa de tensión de confinamiento que es el valor introducido por el usuario
( .

Figura 11. Variación de E50 con la profundidad

Es conocido que la rama de descarga-recarga es sensiblemente más vertical que la de la carga

noval, con una relación comprendida entre 3 y 6, por lo que para poder establecer una relación
con el E50 se emplea una relación equivalente a la indicada para el E50 que se muestra en la
siguiente ecuación

( ) [9]

[Link]. Ley de endurecimiento por el cortante

La formulación analítica de f1 se puede descomponer conceptualmente en dos partes. Los dos

primeros términos corresponden a su relación con el estado tensional y el último con la
deformación plástica de corte.

[10]

Si se toma la deformación plástica de corte como

[11]

Si se considera que el suelo no es excesivamente blando, se puede adoptar sin cometer un

excesivo error que por lo que

[12]

Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 11

 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

luego

[13]

[14]

Descomponiendo la deformación axial (ε1) en sus componentes elástica y plástica se obtiene

( ) [15]

Siendo la parte elástica de la deformación en la dirección axial y en

las direcciones perpendiculares.

La relación entre la deformación axial y la carga desviadora de la expresión [X] se puede reescribir
como

[16]

Que representando en un sistema ε1/q y ε1 la hipérbola se transforma en una recta (Figura 12).

Figura 12. Representación de una curva σ-ε de un ensayo triaxial.

[Link]. Ley de flujo plástico y dilatancia

Se conoce como dilatancia a la relación que existe entre los incrementos de deformación
volumétrica plástica ̇ y los incrementos de deformación de corte plástica ( ̇ según la
siguiente expresión:

̇ ̇ [17]

Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 12

 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

La ley de flujo plástico se deriva de la superficie de potencial plástico que se expresa según la
siguiente expresión

[18]

Siendo la dilatancia movilizada, que se calcula según el valor de (ángulo de rozamiento

movilizado) con las siguientes expresiones

[19]

[20]

El concepto del ángulo de fricción movilizado es sencillamente el valor de la pendiente de la

tangente al círculo de Mohr, que lógicamente va evolucionando según el estado tensional de
cada punto, y tiene la siguiente expresión:

[21]

El ángulo de rozamiento en estado crítico es una propiedad del material y no depende de su

estado tensional ni de su índice de huecos, y habitualmente adopta un valor de 30º para los
suelos granulares y valores inferiores para los suelos cohesivos.

La relación entre estos ángulos de rozamiento se puede obtener con la siguiente expresión

[22]

En la Figura 13 se puede ver la evolución de la dilatancia movilizada en función del estado

tensional, con un valor nulo en el dominio contractivo y un valor variable en el dominio dilatante,
hasta alcanzar el valor máximo introducido por el usuario (Ψ).

Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 13

 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

Figura 13. Explicación gráfica de la dilatancia movilizada Ψm

Esto significa que la dilatancia puede darse para valores mayores del ángulo de fricción
movilizado mientras que para menores relaciones tensionales ( ) el material
se comporta como contractivo y el ángulo de dilatancia se limita al valor de 0 (Figura 14).

Figura 14. Tipos de comportamiento de los suelos en función de la dilatancia. Tomado de Obrzud (2018)

La relación entre la deformación axial y la deformación volumétrica total está definida por la
relación

[23]

El valor de d es la pendiente de la curva que se puede obtener en un ensayo triaxial

drenado. Para evitar que el valor de sea muy elevado, se suele introducir un valor máximo del
índice de huecos ( de forma que si se alcanza este valor, se adopta (Figura 15).

Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 14

 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

Figura 15. Gráfico de dilatancia en una simulación de ensayo triaxial. Tomado de Obrzud (2018)

2.3.3. Reproducción del fenómeno de compresión volumétrica

Según se ha indicado anteriormente, el mecanismo de cortante no genera deformaciones

volumétricas plásticas en el dominio contractante de los suelos, por lo que es necesario
incorporarlo para tener presente las deformaciones volumétricas que se observan cuando el
estado tensional supera la tensión de preconsolidación y el suelo se comporta como un material
normalmente consolidado.

Para tener presente este comportamiento el modelo HSM incorpora una segunda superficie de
fluencia f2 con forma de elipsoide cerrado de formulación analítica

[24]

Donde M es el parámetro que controla la forma del elipsoide y que se relaciona con , es la
presión de preconsolidación y define el punto de corte del elipsoide con el eje p’ (Figura 16 y
Figura 17). En esta superficie de fluencia f2 se considera una plasticidad asociada4, siendo el
vector de flujo plástico perpendicular a la superficie de fluencia f2.

Figura 16. Esquema de la trayectoria de tensiones y la superficie de fluencia f2

4
La superficie de fluencia es la misma que la superficie de potencial plástico
Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 15
 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

Figura 17. Representación 3D del elipsoide que incorpora la compresión volumétrica. Tomado de Obrzud (2018)

La evolución de las deformaciones volumétricas plásticas en función de la variación de la

presión de preconsolidación se produce según la siguiente expresión:

( ) [25]

Siendo H el parámetro que controla la variación de deformación plástica volumétrica y está

relacionado con el módulo edométrico a una presión determinada de referencia ( ).

El incremento de deformación volumétrica plástica se obtiene con la siguiente expresión:

( ) [26]

En relación a los parámetros del modelo HSM, y como estimación inicial, se pueden adoptar las
recomendaciones indicadas por Vermeer & Schanz (1997) que proponen para los diferentes
materiales, y a falta de ensayos triaxiales de los mismos, las siguientes aproximaciones:

Materiales cohesivos

E50ref ≈ 2Eoedref

Eurref ≈ 6 E50ref

m ≈ 0,75 a 1

Materiales granulares

E50ref ≈ Eoedref

Eurref ≈ 4 E50ref

m ≈ 0,5

Este modelo constitutivo presenta como principales ventajas una buena capacidad predictiva
para problemas estáticos, incluye plasticidad por corte (como los suelos muestran en los ensayos
triaxiales) y por compresión (según se registra en los ensayos edométricos) y es adecuado para

Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 16

 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

reproducir los procedimientos constructivos. Sin embargo como desventajas tiene que no
permite la anisotropía, no incluye el efecto del tiempo (creep) ni evoluciona el índice de huecos
(e), por lo que no alcanza el estado crítico ni modifica el peso específico de los materiales cuando
se simula una compactación.

Adicionalmente, cuando se trabaja en rangos de deformación muy reducidos, tal y como ocurre
en los rangos iniciales de una carga noval o en descarga/recarga, el valor es controlado por el
módulo Eur, que presenta un valor sensiblemente inferior que el obtenido en los ensayos de
laboratorio. Es decir, el ajuste del valor de deformación en el rango de la relación de tensión
tangencial con la tensión de rotura comprendido entre 0,2 y 1, es adecuado en el modelo HS
(sombreado amarillo en la Figura 19), pero no es válido en el rango 0 a 0,2.

Figura 18. Reproducción correcta de un ensayo edométrico con un modelo HS.

Figura 19. Rango de validez del módulo de descarga en el modelo HS

Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 17

 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

Para corregir este aspecto, se ha implementado un “parche” que se aplica sobre el modelo HS
que se denomina HS en pequeñas deformaciones (HS Small) que se describe en el siguiente
apartado.

2.4. MODELO HARDENING SOIL EN PEQUEÑAS DEFORMACIONES (HS SMALL)

Para mejorar el comportamiento del modelo HS en pequeñas deformaciones se ha incluido una
variación respecto a aquél en el modelo HS Small.

Básicamente consiste en incorporar una nueva formulación para controlar la deformabilidad a

muy baja deformación empleando una variable que discrimina el tipo de formulación a emplear
en cada momento.

Figura 20. Rango de validez del módulo de descarga en el modelo HS Small

( ) [27]

Con esta expresión se controla la rigidez inicial tanto en carga noval como en descarga-recarga
(pendiente de las líneas verdes de la Figura 22), ajustándose una hipérbola de baja deformación
según la siguiente expresión:

[28]
⁄

Siendo γ0,7 el valor de deformación de corte en el que el módulo de corte se reduce al 70% del
valor G0.

Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 18

 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

Figura 21. Ejemplo de la curva de degradación de la deformabilidad del suelo

De esta manera se ajusta el valor de las curvaturas para realizar una buena conexión con la
hipérbola de grandes deformaciones, que sigue la formulación del modelo HS.

Otra ventaja que incorpora el HS Small es la creación de un lazo de histéresis en los procesos de
carga y descarga, que permite incorporar el amortiguamiento histerético en los cálculos
dinámicos a realizar (Figura 22).

En este modelo constitutivo se hereda del modelo HS la función potencial de Vermeer-de Borst,
el endurecimiento por cortante, el criterio de rotura de MC y el comportamiento de plastificación
debido a la compresión.

En relación a la dilatancia a cortante, en lugar de emplear la formulación de Vermeer-de Borst

incluida en el HS, se adopta la de Dafalias que permite un cálculo numéricamente más estable,
pero que apenas tiene influencia en los resultados.

Figura 22. Determinación de la pendiente inicial controlada por G0 y forma de la hipérbola de baja deformación

Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 19

 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

Por lo tanto, el empleo del modelo HS Small únicamente requiere la introducción adicional de dos
valores para controlar el comportamiento a pequeñas deformaciones: G0 y γ0,7. El resto de los
valores requeridos son idénticos a los requeridos por el modelo HS.

2.5. COMPARATIVA ENTRE LOS MODELOS MC Y HSM

Las diferencias entre ambos modelos ya se han explicado de forma detallada en los apartados
anteriores, por lo que en este apartado únicamente se muestran resultados gráficos de dos tipos
de cálculo (excavación al amparo de unos elementos de contención y un túnel).

La simple comparativa de los resultados obtenidos permite establecer que para estos tipos de
cálculo, en los que influye claramente la deformabilidad del terreno, es muy recomendable el
empleo de un modelo más avanzado como es el HS en lugar del MC.

Figura 23. Ejemplo de cálculo de excavación con pantallas

Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 20

 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

Figura 24. Ejemplo de cálculo de excavación de un túnel

Figura 25. Ejemplo de cálculo de excavación de un túnel

Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 21

 CURSO DE MODELIZACIÓN DE MUROS PANTALLA CON PLAXIS 2D

Figura 26. Comparativa de resultados obtenidos con varios modelos y las mediciones realizadas

Módulo 2. Comportamiento del terreno. Modelos constitutivos (I) Página 22